Modèles de sensibilité dans le cadre de la méthode de Monté-Carlo : illustrations en transfert radiatif

par Maxime Roger

Thèse de doctorat en Énergétique et transferts

Sous la direction de Richard Fournier.

Soutenue en 2006

à Toulouse, INPT .


  • Résumé

    En physique du transport, la méthode Monté-Carlo est abordée comme une méthode de simulation numérique de phénomènes stochastiques. Cette approche permet une compréhension intuitive de la méthode grâce aux images statistiques de transport corpusculaire qui accompagnent habituellement le développement d'algorithmes. A partir de ce type d'algorithmes, nous proposons une démarche méthodologique permettant de développer une forulation intégrale de la grandeur estimée. Cette écriture intégrale permet notamment d'utiliser diverses techniques d'optimisation, essentiellement en vue d'une diminution des temps de calculs. Une fois l'intégrale explicitée, nous proposons une méthodologie basée sur la méthode de Monté-Carlo permettant d'ajouter un code existant à des estimations de sensibilités pour un coût supplémentaire en temps de calcul relativement faible. Cette méthodologie s'applique à tous le types de sensibilités. Dans le cas de sensibilités à la déformation du domaine (sensibilité à un paramètre dont le domaine d'intégration est dépendant), un champ de vitesse de déformation du domaine est introduit de manière arbitraire. Les incertitudes statistiques associées aux estimations de sensibilités étant dépendantes du choix de ce champs de vitesse, nous proposons une technique systématique de construction d'un tel champ assurant une bonne qualité de convergence. Cette méthodologie d'estimation simultanée de sensibilités a été évaluée en transfert radiatif sur des exemples académiques inspirés de problématiques courantes en imagerie médicale, en énergétique et en synthèse d'image. Chacun de ces exemples ayant une complexité algorithmique particulière, nous avons pu identifier les principales difficultés de mise en oeuvre dans le cadre du transport de photons et proposer une série de premières solutions effectives pour les surmonter

  • Titre traduit

    Sensitivity models for Monte Carlo methods : illustrations in radiative transfer


  • Résumé

    In the frame of transport physics, the Monte Carlo method is considered as a numerical tool for simulation of stochastic phenomena. This approach allows an intuitive use of the Monte Carlo method following the statistical phisics of corpuscular transport, leading to fast developments of suited algorithms. In this work, we propose a methodology that associates an integral formulation of the estimated quantities to any such intuitive Monte Carlo algorithm. The integral formulation makes it possible in particular to benefit from all available Monte Carlo optimization techniques. Starting from the integral formulation, it is then shown that when a Monte Carlo algorithm is used for the estimation of any physical quantity A, a fast additional procedure can be implemented that simultaneously estimates the sensivity of A to any problem parameter. When a change in the parameter changes the itegration domain (domain-deformation sensivities), a domain-deformation velocity field is requiered. The statistical uncertainties of the sensitivity estimations depend on the choice of this velocity field. Practically speaking, a systematic procedure is proposed here to design velocity fields that insure good convergence qualities. The corresponding Monte Carlo methodology for sensitivity estimation has been tested on radiative transfer applications. Academic configurations related to current problems of medical imaging, energetics and computer grafics are presented. Diferent level of complexity arise from these academics examples and allow us to identify the main implementation dificulties leading us to the proposition of a serie of the first efective solution for algorithmic implementation

Consulter en bibliothèque

La version de soutenance existe sous forme papier

Informations

  • Détails : 1 vol. (142 p.)
  • Notes : Publication autorisée par le jury
  • Annexes : Bibliogr. p. 137-142

Où se trouve cette thèse ?

  • Bibliothèque : École nationale supérieure d'ingénieurs en arts chimiques et technologiques. Médiathèque.
  • Disponible pour le PEB
  • Cote : TG 2006 ROG
  • Bibliothèque : IMT Mines Albi. Centre de documentation.
  • Disponible pour le PEB
  • Cote : P130-ROG
Voir dans le Sudoc, catalogue collectif des bibliothèques de l'enseignement supérieur et de la recherche.