Modélisation du flux critique à l’aide des méthodes de Boltzmann sur réseau : application aux dispositifs expérimentaux du RJH
Les méthodes LBM (Lattice Boltzmann Methods) sont des techniques numériques utilisées pour simuler des phénomènes de transport dans des systèmes complexes. Elles permettent de modéliser le comportement des fluides en termes de particules qui se déplacent sur une grille discrète (un "réseau" ou lattice). Contrairement aux méthodes classiques, qui résolvent directement les équations différentielles des fluides, les méthodes LBM simulent l'évolution des fonctions de distribution des particules de fluide dans un espace discret, en utilisant des règles de propagation et de collision.
Le choix du réseau dans les méthodes LBM est une étape cruciale, car il affecte directement la précision, l'efficacité et la stabilité des simulations. Le réseau détermine la manière dont les particules de fluide interagiront et se déplaceront dans l'espace, ainsi que la façon dont la discrétisation de l'espace et du temps est effectuée.
Les méthodes LBM présentent un parallélisme naturel, car les calculs à chaque point de la grille sont relativement indépendants. Bien que les méthodes classiques de CFD, basées sur la résolution des équations de Navier-Stokes, puissent aussi être parallélisées, les termes non linéaires peuvent rendre le parallélisme plus difficile à gérer, en particulier pour les modèles impliquant des écoulements turbulents ou des maillages irréguliers. Les méthodes LBM permettent donc, à moindre coût, de capturer des phénomènes complexes. Des travaux récents ont notamment montré qu'il était possible, avec les LBM, de retrouver la courbe de refroidissement de Nukiyama (ébullition en vase) et, ainsi, de calculer avec précision le flux critique. Ce flux correspond à une ébullition en masse, appelée crise d’ébullition, qui se traduit par une dégradation soudaine du transfert thermique.
Le flux critique représente un enjeu crucial pour le Réacteur Jules Horowitz, car les dispositifs expérimentaux (DEX) sont refroidis par de l'eau en convection naturelle ou forcée. Ainsi, afin de garantir le bon refroidissement des DEX et la sûreté du réacteur, il convient de s'assurer que, sur la gamme de paramètres étudiés, le flux critique ne soit pas atteint. Il doit donc être déterminé avec précision.
L'étudiant sera amené, dans un premier temps, à définir un réseau pour appliquer les méthodes LBM sur un dispositif du RJH en convection naturelle. Il consolidera les résultats obtenus en les comparant aux données disponibles. Enfin, des calculs exploratoires en convection forcée (régime laminaire à turbulent) seront menés.
Implémentation d’algorithmes parallèles sur GPU pour les simulations du combustible nucléaire sur supercalculateurs exaflopiques
Dans un contexte où le calcul haute performance (HPC) est en perpétuelle évolution, le design des nouveaux supercalculateurs tend à intégrer toujours plus d’accélérateurs ou de cartes graphiques (GPUs), qui fournissent l’essentiel de la puissance de calcul de la plupart des supercalculateurs actuels. En raison de leurs différences architecturales par rapport aux unités centrales de calcul (CPUs) et des environnements logiciels en constante évolution, les GPUs posent de profonds défis de programmation. Une utilisation efficace de leur puissance de calcul demande une refonte des algorithmes et logiciels de simulation existants pour atteindre un parallélisme massif.
Le CEA a développé la plateforme de calcul PLEIADES dédiée à la simulation du comportement des combustibles, depuis la fabrication jusqu’au comportement en réacteur, puis lors du stockage. Elle inclut une parallélisation en mémoire distribuée MPI permettant une parallélisation sur plusieurs centaines de cœurs. Cette plateforme répond aux exigences des partenaires du CEA que sont EDF et Framatome, mais il est nécessaire de l’adapter pour les nouvelles infrastructures de calcul. Proposer une solution flexible, portable et performante pour les simulations sur des supercalculateurs équipés de GPUs est d'un intérêt majeur afin de capturer des physiques toujours plus complexes sur des simulations comportant des domaines de calcul toujours plus grands.
Dans ce cadre, la thèse visera d’élaborer puis évaluer différentes stratégies de portage de noyaux de calculs sur GPU ainsi que l’utilisation de méthodes de répartition dynamique de la charge adaptés aux supercalculateurs GPUs actuels et futurs. Le candidat s’appuiera sur des outils développés au CEA comme les solveurs thermo-mécaniques MFEM-MGIS [1,2] ou MANTA [3]. Les solutions logicielles et algorithmes parallèles proposés avec cette thèse permettront à terme la réalisation de grands calculs de modélisation multi-physique en 3D du comportement des crayons combustibles sur des supercalculateurs comportant des milliers de cœurs de calcul et des GPUs.
Le candidat travaillera au sein du Laboratoire de développement des Outils de Calcul Scientifique (OCS) combustibles PLEIADES (LDOP) au département d'études des combustibles (DEC - Institut IRESNE, CEA Cadarache). Il sera amené à évoluer dans une équipe pluridisciplinaire composée de mathématiciens, physiciens, mécaniciens et informaticiens. A terme, les contributions de la thèse visent à enrichir la plate-forme numérique pour la simulation de combustibles nucléaires PLEIADES.
Références : [1] MFEM-MGIS - https://thelfer.github.io/mfem-mgis/; [2] Th. Helfer, G. Latu. « MFEM-MGIS-MFRONT, a HPC mini-application targeting nonlinear thermo-mechanical simulations of nuclear fuels at mesoscale ». IAEA Technical Meeting on the Development and Application of Open-Source Modelling and Simulation Tools for Nuclear Reactors, June 2022, https://conferences.iaea.org/event/247/contributions/20551/attachments/10969/16119/Abstract_Latu.docx, https://conferences.iaea.org/event/247/contributions/20551/attachments/10969/19938/Latu_G_ONCORE.pdf; [3] O. Jamond et al. «MANTA : un code HPC généraliste pour la simulation de problèmes complexes en mécanique », https://hal.science/hal-03688160
Étude des transitoires de cristallisation des oxalates d’actinides
Les besoins croissants en énergie et l’urgence climatique nécessitent une transition rapide vers une énergie totalement décarbonée, mixant énergies renouvelables et nucléaire durable. Dans ce contexte, la précipitation du plutonium et de l’uranium sous forme d’oxalate constitue une étape clef du procédé industriel de recyclage du combustible usé. La compréhension fine des mécanismes de cristallisation de ces oxalates constitue ainsi un enjeu fort pour un meilleur pilotage de ces opérations.
Or, il est désormais largement admis que les ions en solution s’assemblent en cristaux via une série d’états transitoires non cristallins, ce qui contredit fondamentalement toutes les théories classiques de nucléation utilisées dans les modèles de précipitation. En particulier, nous avons démontré ces dernières années que les cristaux d'oxalate de terres rares (Eu, Nd, Ce, Tb), certains utilisés pour simuler expérimentalement le recyclage de l'uranium et du plutonium, se forment via des nanogouttelettes liquides riches en réactifs qui se séparent du solvant aqueux. Ce comportement modifie la vision jusqu’à présent retenue pour la précipitation de ces oxalates et amène à s’interroger sur le comportement des oxalates d’actinide.
Le but de cette thèse est de confirmer ou infirmer que les gouttelettes minérales transitoires se forment également lors de la formation des oxalates d'uranium et de plutonium, et de déterminer si les transitoires de la cristallisation impactent les modèles de précipitation utilisés pour calibrer le procédé de recyclage du combustible nucléaire. Cette étude aura un impact non seulement sur les processus de précipitation utilisés dans le recyclage, mais fera également avancer une question fondamentale sur la cristallisation « non classique » débattue depuis longtemps.
Essais thermomécaniques jusqu’à de hautes températures sur une céramique nucléaire irradiée
Ces travaux de thèse s’inscrivent dans le cadre des études sur les interactions pastille-gaine dans les crayons combustibles utilisés dans les réacteurs des centrales nucléaires à eau pressurisé. L’exploitant doit assurer et démontrer l’intégrité de l’ensemble des crayons, en toutes situations. Or, les contraintes mécaniques subies par la gaine, première barrière de sûreté, sont reliés aux propriétés viscoplastiques du combustible. Il faut donc connaître ces propriétés et leurs évolutions en fonctionnement.
Le sujet de thèse proposé s’intéresse particulièrement à la caractérisation en cellule de haute activité de combustibles irradiés en réacteur nucléaire. Une des difficultés majeures réside dans le fait que les combustibles issus d’une irradiation en réacteur sont multi-fissurés, ce qui rend leurs caractérisations mécaniques particulièrement compliquées. Toutefois, une thèse en cours (2022-25) a permis de (i) concevoir une machine spécifique d’essais thermomécaniques, (ii) qualifier en partie ce dispositif, (iii) mettre en place des outils et méthode d’extraction d’échantillon fissuré, (iv) ainsi qu’une modélisation intégrale du système (jumeau numérique).
La thèse proposée ici correspond à la suite de ces travaux et sera construite en quatre étapes sur trois plateformes expérimentales disponibles au CEA
1. Prise en main et amélioration des outils numériques et expérimentaux existant,
2. Implémentation du dispositif en cellule de haute activité sur un four existant,
3. Réalisation d’essai thermomécanique sur combustible irradié ; ce qui constituera en soi une première mondiale dans ces conditions.
Les essais demanderont des post-traitements dédiés basés sur une discussion essai-simulation. En premier lieu, cela permettra de décolérer les phénomènes et de se focaliser sur le comportement viscoplastique du matériau. Une fois la base expérimentale suffisamment étoffée et interprétée, il sera alors possible de conforter ou revoir les lois de comportement du combustible. Un lien avec la microstructure des matériaux pourrait être abordé.
Tout au long de ces étapes, le thésard s’appuiera sur les compétences et expertises de différents laboratoires du Département d'études des Combustible (Institut IRESNE, CEA Cadarache) et sur un environnement collaboratif académique. Cette thèse s'inscrit également dans le cadre du projet européen OPERA HPC et en est un enjeu majeur.
Le doctorant devra présenter un goût marqué pour l’approche expérimentale et quelques facilités pour l’utilisation d’outils numériques. Des connaissances en science des matériaux sont le minimum requis. Au long des trois années de doctorat, le doctorant améliorera ses compétences multiphysiques en conception de dispositifs expérimentaux et comportement des matériaux à haute température, ainsi qu’en simulation numérique, ce qui facilitera son insertion professionnelle. Au regard de l’ampleur de la tâche, il serait préférable, mais non obligatoire, que le candidat ait suivi préalablement le stage de master lié à ce sujet (voir par ailleurs).
Étude multi-échelle du transport ionique dans des matériaux nanoporeux hiérarchiques non saturés : application aux matériaux cimentaires
Le transport ionique est crucial pour déterminer la durabilité des matériaux à base de ciment et, par conséquent, l'extension de la durée de vie des (infra)structures en béton. Les phénomènes de transport déterminent la capacité de confinement du béton, essentielle à la conception et à la gestion des infrastructures en béton pour la production d'énergie. Dans la plupart des conditions de service, le béton se trouve dans un état non saturé. Un transport anormal a été observé dans les matériaux à base de ciment, et les raisons de ces écarts par rapport au comportement attendu d'autres matériaux poreux peuvent provenir de processus à l'échelle nanométrique.
A ce jour, la majorité des modélisations prédictives de la durabilité ne tiennent pas explicitement compte des processus à l’échelle nanométrique, pourtant fondamentaux pour déterminer les propriétés de transport. Des progrès récents ont été réalisés dans la quantification du comportement de l'eau confinée dans diverses phases présentes dans les systèmes cimentaires. Les silicates de calcium hydratés (C-S-H) sont la principale phase hydratée dans les matériaux à base de ciment et présentent des nanopores dans les gammes microporeuses et mésoporeuses. Cependant, les effets de la désaturation restent encore à élucider pleinement. Une compréhension fondamentale des processus de transport nécessite un cadre multi-échelle dans lequel l'information de l'échelle moléculaire se répercute à travers d'autres échelles pertinentes (en particulier, l'échelle mésoscopique associée à la porosité du gel C-S-H (~nm), la porosité capillaire et la zone de transition interfaciale (~µm) jusqu'à l'échelle macroscopique des applications industrielles dans les matériaux à base de ciment).
L’objectif de ce travail de doctorat est d’évaluer le transport ionique des chlorures, une espèce critique pour la durabilité du béton, en conditions non saturées en combinant des simulations à petite échelle, une modélisation multi-échelle et des expérimentations dans une approche ascendante. Le travail se concentrera sur le C-S-H. Le projet vise à caractériser les effets de la désaturation sur les processus nanométriques qui gouvernent le transport des chlorures.
Modélisation des équilibres de complexation des actinides en milieu nitrique. Application au procédé PUREX
Le code de calcul PAREX+ est un outil majeur dans le domaine de la chimie séparative. Il permet la modélisation et la simulation des procédés de séparation basés sur l’extraction par solvant. Dans ce code, la distribution des espèces d’intérêt entre les phases aqueuse et organique est calculée en tout point du procédé aussi bien en régime établi que dynamique. L’objectif de la thèse est d’améliorer le modèle de distribution présent au sein de ce code. Pour cela une meilleure compréhension des phénomènes mis en jeu au sein des phases organiques et aqueuses est nécessaire ainsi qu’une nouvelle approche pour les prendre en compte dans le modèle. Cette thèse associe donc l’expérimental et la modélisation. L’étudiant intégrera une équipe d’encadrement composée d’experts en chimie séparative et en modélisation. Son travail sera valorisé par l’émission de publications et de participations à des congrès internationaux. A l’issue de cette thèse, l’étudiant aura de solides connaissances dans le domaine de l’extraction par solvant et de sa modélisation qu’il pourra valoriser auprès des industriels ou des organismes de recherche du nucléaire ou dans les autres domaines de la chimie séparative (séparation des terres rares ou hydrométallurgie).
Modélisation thermo-chemo-mécanique du frittage : effet de l’atmosphère et de la densification différentielle sur le retrait des pastilles
Les combustibles de dioxyde d’uranium (UO2), utilisés dans les centrales nucléaires sont des céramiques, dont le frittage en phase solide est une étape-clé de la fabrication. L’étape de frittage consiste en un traitement thermique sous pression partielle contrôlée de O2 permettant de consolider, densifier le matériau et faire grossir les grains de UO2. La densification induit un retrait macroscopique de la pastille. Si le compact (poudre comprimée par pressage avant le frittage) admet de fortes hétérogénéités de densité, une différence de densification dans la pastille peut avoir lieu entraînant un retrait différentiel et l’apparition de défauts.
Cette thèse se consacre à la mise en place d’un modèle thermo-chemo-mécanique du frittage pour simuler l’impact de la composition et les propriétés physiques de l’atmosphère sur la densification du combustible à l’échelle de la pastille. Cette échelle nous permettra de considérer les gradients de densité issus du pressage, mais également de prendre en compte la cinétique de diffusion d’oxygène impactant localement la vitesse de densification qui elle-même impactera le processus de transport. Une simulation multiphysique est nécessaire pour simuler le couplage de ces phénomènes.
Ce travail de thèse sera mené au sein du Laboratoire commun MISTRAL (Aix-Marseille Université/CNRS/Centrale Marseille et l'institut IRESNE du CEA-Cadarache). Le doctorant valorisera ses résultats au travers de publications et participations à des congrès et aura acquis de solides compétences qui sont recherchées et valorisables dans un grand nombre de domaines académiques et industriels.
Simulations multiphysiques avec estimation d’incertitudes appliquées aux réacteurs rapides refroidis au sodium
La modélisation multiphysique est essentielle pour l'analyse des réacteurs nucléaires, mais la propagation des incertitudes entre différents domaines physiques (comme les comportements thermiques, mécaniques et neutroniques) est souvent négligée en raison de sa complexité. Ce projet de thèse vise à relever ce défi en développant des méthodes innovantes pour intégrer la quantification des incertitudes dans les modèles multiphysiques.
L'objectif principal est de proposer des approches de modélisation optimales, adaptées à différents niveaux de précision. Le projet explorera des techniques avancées, telles que la modélisation d'ordre réduit et l'expansion du chaos polynomial, pour identifier et classer les paramètres d'entrée ayant l'impact le plus significatif sur les sorties du système, indépendamment de leur domaine physique. Une comparaison entre des modèles « haute fidélité », développés à l'aide des outils de simulation de référence du CEA, et des modèles « best-estimate » conçus pour un usage industriel sera réalisée. Cette analyse comparative mettra en lumière comment les erreurs se propagent dans les différentes approches de simulation.
Les modèles seront validés à l'aide de données expérimentales de SEFOR, un réacteur rapide refroidi au sodium. Ces expériences fournissent des repères précieux pour tester les modèles multiphysiques dans des conditions réalistes de réacteur. Ce projet répond directement au besoin croissant de l'industrie nucléaire pour des outils de modélisation fiables et efficaces, visant à améliorer la sécurité et la performance des réacteurs.
Le candidat évoluera dans un environnement dynamique au CEA, avec accès à des ressources de simulation avancées et des opportunités de collaboration avec d'autres chercheurs et doctorants. Le projet offre également la possibilité de présenter les résultats lors de conférences nationales et internationales, avec des perspectives de carrière solides dans la conception de réacteurs nucléaires, l'analyse de la sûreté et la simulation avancée.
mise en ordre de l’oxygène dans le zirconium : mécanismes, cinétiques et propriétés mécaniques associés
Ce travail de thèse à pour but d'étudier les propriétés des alliages binaires zirconium-oxygène (Zr-O), notamment dans le cadre des applications nucléaires. Traditionnellement, l'oxygène est considéré en solution solide dans la matrice de zirconium, sans formation de composés ordonnés comme Zr6O ou Zr3O. Cependant, des études récentes suggèrent qu’à des températures inférieures à 600°C, des composés ordonnés peuvent se former, affectant la limite de solubilité de l'oxygène. Ces composés, observés après des traitements thermiques, pourraient modifier les propriétés mécaniques des alliages Zr-O, en particulier à température ambiante et jusqu’à 350°C. La thèse proposée cherche à comprendre ces mécanismes à travers des expériences de diffraction des rayons X et de microscopie électronique, pour étudier l'arrangement de l'oxygène, la stabilité thermique des composés et leur impact sur la déformation plastique. Cela vise à optimiser l’utilisation de ces alliages dans les réacteurs nucléaires.
PREDICTION PAR LA SIMULATION DES VIBRATIONS DANS LES CENTRIFUGUESES
Les machines tournantes sont des équipements critiques dans de nombreuses installations industrielles et leur exploitation s’accompagne régulièrement de problèmes d’équilibrage qui occasionnent l’apparition de vibrations potentiellement dangereuses pour les opérateurs et les équipements. La décanteuse pendulaire centrifuge est parfois le siège de vibrations qui obligent l’exploitant à ralentir la cadence de production. L’environnement nucléaire dans lequel ces équipements sont exploités ne permet pas de réaliser les mesures et observations nécessaires à une étude purement expérimentale. L’objectif est donc de réaliser une modélisation à partir de données limitées afin d’amener à une compréhension fine des phénomènes en jeu. Le sujet de thèse se propose de coupler des simulations CFD de type Euler-Euler de la répartition de masse dans le bol tournant avec une modélisation masse-ressort des liaisons mécaniques afin de s’approcher des signaux vibratoires mesurés industriellement. Un tel outil numérique serait une aide précieuse pour explorer les diverses pistes, sources potentielle d’apparition d’un déséquilibre de masse, et cela sans avoir à passer par une reproduction expérimentale. Associé à des méthodes de deep learning, ce type de modèle permettrait aussi de bâtir un prédicteur de balourds à partir de courts signaux vibratoires ouvrant la porte à un pilotage actif de la décanteuse.