Méthode de Monte-Carlo et non-linéarités : de la physique du transfert radiatif à la cinétique des gaz

En physique du transport, en particulier en physique du transfert radiatif, la méthode de Monte-Carlo a été développée à l'origine comme la simulation de l'histoire d'un grand nombre de particules, dont on déduit des observables moyennes. Cette méthode numérique doit son succès à plusieurs qualités : une gestion naturelle des espaces des phases aux nombreuses dimensions, une erreur systématique nulle par rapport au modèle physico-mathématique, les intervalles de confiance donnés avec les résultats, une capacité à prendre en compte simultanément de nombreux phénomènes physiques, la possibilité de calcul de sensibilités simultané, et une parallélisation aisée. En cinétique des gaz, les particules collisionnent entre elles et non pas avec un milieu extérieur ; on dit que leur transport est non-linéaire. Ces collisions mutuelles mettent en défaut l'approche évoquée ci-dessus de la méthode de Monte-Carlo ; car pour simuler des trajectoires indépendantes de multiples particules et ainsi estimer leur distribution, il faut connaître au préalable exactement cette même distribution...Cette thèse fait suite à celles de Jérémi DAUCHET (2012) et de Mathieu GALTIER (2014), consacrées au transfert radiatif. Entre autres travaux, ces auteurs montraient comment la méthode de Monte-Carlo peut s'accommoder de non-linéarités, en gardant son formalisme et ses spécificités habituelles. Les non-linéarités alors franchies étaient respectivement une loi de couplage chimie/luminance, et la dépendance de la luminance envers le coefficient d'absorption. On essaie dans ce manuscrit d'outrepasser la non-linéarité du transport. Pour cela, nos principaux outils sont un suivi des particules en remontant le temps, basé sur des formulations intégrales des équations de transport, formulations largement inspirées des algorithmes dits à collisions nulles. Nous montrons, sur plusieurs exemples académiques, que nous avons en effet étendu la méthode de Monte-Carlo à la résolution de l'équation de Boltzmann. Ces exemples sont aussi l'occasion de tester les limites de ce que nous avons mis en place. Les résultats les plus marquants sont certainement l'absence totale de maillage dans la méthode numérique, ainsi que sa capacité à calculer correctement les quantités de particules de haute énergie cinétique (toujours peu nombreuses par rapport au total, en cinétique des gaz). Au-delà des exemples fournis, ce manuscrit est voulu comme un essai de formalisme et une exploration des bases de la méthode dév