Improving Monte Carlo reactor physics simulations using adaptive sampling of neutron histories
Auteur / Autrice : | Kévin Fröhlicher |
Direction : | Eric Dumonteil, Alain Hébert, Julien Taforeau, Mariya Brovchenko |
Type : | Thèse de doctorat |
Discipline(s) : | Sciences de l'énergie nucléaire |
Date : | Soutenance le 01/02/2023 |
Etablissement(s) : | université Paris-Saclay |
Ecole(s) doctorale(s) : | École doctorale Particules, hadrons, énergie et noyau : instrumentation, imagerie, cosmos et simulation (Orsay, Essonne ; 2015-....) |
Partenaire(s) de recherche : | Laboratoire : Institut de Radioprotection et de Sûreté Nucléaire (Fontenay-aux-Roses ; 2002-....) - Laboratoire de Neutronique (IRSN) (Fontenay-aux-Roses ; 2020-....) |
Référent : Faculté des sciences d'Orsay | |
graduate school : Université Paris-Saclay. Graduate School Physique (2020-....) | |
Jury : | Président / Présidente : Tony Lelièvre |
Examinateurs / Examinatrices : Elsa Merle, Eugene Shwageraus, Margaux Faucher, Cheikh M'Backé Diop | |
Rapporteur / Rapporteuse : Elsa Merle, Eugene Shwageraus |
Résumé
L'Adaptive Multilevel Splitting (AMS) est une méthode de réduction de variance précédemment mise en œuvre pour le transport des neutrons dans les calculs de radioprotection afin d'améliorer l'efficacité des calculs pour la modélisation d'événements rares.Cette thèse vise à présenter une extension du champ d'application de cette technique aux problèmes de physique des réacteurs en régime permanent et en cinétique/dynamique.En ce qui concerne les calculs en régime permanent ou de criticité, l'agorithme inspiré de la méthode d'itération sur la puissance et sur lequel ces calculs sont généralement basés peut présenter des problèmes.En effet, la normalisation de la population effectuée entre les générations amplifie les fortes corrélations spatiales et générationnelles déjà induites par le processus de naissance et de mort (processus de Galton-Watson) dans le transport des neutrons en milieu multiplicatif.Dans la plupart des situations pathologiques, ces corrélations peuvent provoquer la formation de clusters de neutrons, dépeuplant ainsi certaines parties du système et entraînant des problèmes de convergence.Dans cette thèse, l'itération sur la puissance classique a été refondue en un problème de réduction de variance pour lequel l'AMS a été utilisé afin de résoudre ce problème de clustering.Dans le but de réaliser une preuve de principe, l'AMS a été implémentée dans un code Monte Carlo maquette entièrement développé pendant ce doctorat, pouvant modéliser des problèmes de criticité/statiques avec une modélisation physique simplifiée.La méthode a été caractérisée sur un problème homogène avec des neutrons mono-énergétiques, et les résultats obtenus se sont avérés encourageants pour ce qui est de la réduction des corrélations indésirables dans les calculs de criticité.Une fois testée avec succès dans la maquette, l'AMS a donc été implémentée dans le code Monte Carlo SERPENT2 afin d'étendre son utilisation aux problèmes de cinétique et de dynamique.En effet, malgré les récents développements en matière de réduction de variance, les coûts de calcul associés à ce type de calculs demeurrent très élevés.A cet égard, il a été envisagé de s'appuyer sur l'utilisation d'une carte d'importance par l'AMS pour réduire la variance associée à l'estimation des scores locaux lors des transitoires, améliorant ainsi l'efficacité des calculs de cinétique.Pour ce faire, différentes combinaisons de méthodes ont été testées lors de l'implémentation de l'AMS dans SERPENT2.La caractérisation de l'AMS en cinétique a été principalement faite sur des cas simples tels que des compositions homogènes représentatives d'assemblages de Réacteur à Eau sous Pression (REP), ou des problèmes hétérogènes avec neutrons mono-énergétiques.La méthode a ensuite été appliquée à un transitoire opéré sur un groupe de 3x3 assemblages.Une carte d'importance dépendante du temps a été extraite des équations de cinétique ponctuelle adjointes à cette fin.L'objectif étant, avec le temps, de s'appuyer sur des cartes d'importance dépendantes de l'espace et du temps obtenues à partir de la cinétique espace-temps adjointe pour les neutrons et les précurseurs via des méthodes déterministes.