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. Les méthodes LBM par rapport aux méthodes CFD classiques permettent de mieux capturer certains phénomènes complexes (comme les écoulements multiphasiques, turbulents ou en milieux poreux) car elles reposent sur une modélisation mésoscopique du fluide, directement dérivée de la cinétique des particules, plutôt que sur une résolution macroscopique des équations de Navier–Stokes. Cette approche permet une représentation plus fine des interfaces, des effets non linéaires et des interactions locales, souvent difficiles à modéliser correctement avec les méthodes CFD classiques. 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 les dispositifs expérimentaux (DEX) du Réacteur Jules Horowitz, car ils sont refroidis par de l'eau en convection naturelle (dispositifs de type fuel capsule) ou forcée (dispositifs de type boucle). 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. Les études précédentes menées sur un DEX de type fuel-capsule à l’aide du code NEPTUNE-CFD (méthodes CFD classique) ont montré que la modélisation est limitée à une région située loin du flux critique. De façon générale, les écoulements à fort taux de vide (supérieurs à 10%) ne peuvent être résolues aisément par les approches classiques de la CFD.
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 sur le flux critique obtenus sur cette configuration en les comparant aux données disponibles. Enfin, des calculs exploratoires en convection forcée (régime laminaire à turbulent) seront menés.
L’étudiant sera accueilli au sein de l’institut IRESNE.
Conception et Optimisation d'un Procédé Innovant pour la Capture du CO2
Dans une enquête réalisée en 2023 par la BEI, deux tiers des jeunes français ont affirmé que l’impact climatique des émissions de leur potentiel futur employeur est un facteur important au moment de choisir un emploi. Mais pourquoi s’arrêter là quand vous pouvez choisir de travailler activement pour la réduction de ces émissions, tout dans le cadre d’un sujet de recherche riche et passionnant ? Au Laboratoire de Simulation de Procédés et analyse de Systèmes, nous proposons une thèse qui vise à concevoir et ensuite à optimiser un procédé pour la capture du CO2 dans les rejets gazeux des industries. Son principe de fonctionnement dérive du procédé « Benfield » pour la capture du CO2. Nous proposons des conditions opératoires optimisées pour lesquels le procédé Benfield serait plus performant. Le deuxième axe d’innovation consiste dans une étude de couplage thermique avec une installation industrielle disposant de la chaleur à céder.
La recherche sera menée en collaboration avec le CEA de Saclay et le Laboratoire de Génie Chimique (LGC) à Toulouse. Dans un premier temps, le thésard va réaliser des travaux de simulation numérique à l’aide d’un logiciel de simulation de procédé (ProSIM). Ensuite, il pourra explorer et proposer différentes solutions pour minimiser le besoin énergétique du procédé. Les schémas de procédé obtenus pourront être validés expérimentalement au LGC, où le thésard sera encadré par des experts en procédé de transfert liquide-gaz. Il sera responsable de mettre en place un montage expérimental à l’échelle pilote pour acquérir des données sur les étapes d’absorption et désorption en colonne, avec un garnissage de structure innovante conçu par la fabrication additive. Il conduira lui-même les manips et pourrait éventuellement encadrer un stagiaire pour le support aux acquisitions expérimentales.
Si vous êtes passionné du Génie de Procédés et que vous cherchez un sujet de thèse stimulant et de grand impact pour la société, postulez et intégrez nos équipes !
Cinétiques de ségrégation et précipitation dans les alliages ferritiques sous irradiation : couplage des effets magnétiques, chimiques et élastiques
Les aciers ferritiques sont envisagés comme matériaux de structure dans les futurs réacteurs nucléaires à fission et à fusion. Or ces alliages ont des propriétés tout à fait originales, liées aux couplages entre les interactions chimiques, magnétiques et élastiques qui affectent à la fois leurs propriétés thermodynamiques, la diffusion des espèces chimiques et celle des défauts ponctuels du cristal. Le but de la thèse sera de modéliser à l’échelle atomique l’ensemble de ces effets et de les intégrer dans des simulations Monte Carlo pour modéliser les cinétiques de ségrégation et de précipitation sous irradiation, phénomènes qui peuvent dégrader leurs propriétés d’usage. L’approche atomique est indispensable pour ces matériaux soumis à une irradiation permanente, pour lesquelles les lois de la thermodynamique d’équilibre ne s’appliquent plus.
La candidate ou le candidat recherché(e) devra avoir une bonne formation en physique statistique ou en sciences des matériaux, et être attiré(e) par les simulations numériques et la programmation informatique. La thèse se déroulera au laboratoire de métallurgie physique du CEA Saclay (SRMP) dans un environnement de recherche bénéficiant d’une expérience reconnue en modélisation multi-échelles des matériaux, avec une quinzaine de thèses et de contrats post-doctoraux en cours sur ces thématiques.
Un stage de Master 2 sur le même sujet est proposé pour au printemps 2025 et est vivement recommandé.
Compréhension des mécanismes de dissolution oxydante de (U,Pu)O2 en présence d'Ag(II) généré par ozonation
Le recyclage du plutonium contenu dans les combustibles MOx, constitués d’oxydes mixtes d’uranium et de plutonium (U,Pu)O2, repose sur une étape clé : la dissolution complète du dioxyde de plutonium (PuO2). Or, ce dernier se dissout difficilement dans l’acide nitrique concentré utilisé industriellement. L’ajout d’une espèce fortement oxydante, telle que l’argent(II), permet d’accélérer cette dissolution : c’est le principe de la dissolution oxydante. L’ozone (O3) est utilisé pour régénérer en continu l’oxydant Ag(II) dans le milieu. Si ce procédé a démontré son efficacité, les mécanismes mis en jeu demeurent encore mal connus et peu documentés. Leur compréhension constitue un préalable indispensable à toute industrialisation future.
L’objectif de cette thèse est de mieux comprendre les mécanismes d’interaction dans le système HNO3/Ag/O3/(U,Pu)O2. Le travail proposé s’articulera autour d’une étude expérimentale paramétrique de complexité croissante. Dans un premier temps, les mécanismes de génération et de consommation d’Ag(II) seront étudiés dans le système simple HNO3/Ag/O3. Puis dans un second temps, l’influence de divers paramètres sur la dissolution oxydante de (U,Pu)O2 sera examinée. Ces résultats permettront d’élaborer un modèle cinétique de dissolution en fonction des paramètres étudiés.
A l’issue de cette thèse, le(la) candidat(e), de formation initiale en physico-chimie, maîtrisera un large panel de techniques expérimentales ainsi que des méthodes de modélisation pointues. Cette double compétence lui ouvrira de nombreuses perspectives d’emploi en recherche académique ou en R&D industrielle, tant dans le secteur nucléaire que dans d’autres domaines de la chimie et des matériaux.
Développement d'une approche macroscopique pour la dégradation à long terme des structures en béton sous irradiation
Dans les centrales nucléaires, la protection biologique en béton (CBS) est conçue à proximité de la cuve du réacteur. Cet élément, qui joue également le rôle de structure porteuse, absorbe donc les radiations. Il est ainsi exposé pendant la durée de fonctionnement de la centrale à des niveaux élevés de radiations qui peuvent avoir des conséquences à long terme. Ces radiations peuvent notamment entraîner une diminution des propriétés mécaniques des matériaux et de la structure. Etant donné son rôle clé, il est donc nécessaire de développer des outils et des modèles, pour prédire les comportements de telles structures à l'échelle macroscopique.
Sur la base des résultats existants obtenus à une échelle inférieure - simulations mésoscopiques, à partir desquelles une meilleure compréhension de l'effet de l'irradiation peut être obtenue, et des résultats expérimentaux qui viendront alimentés la simulation (propriétés des matériaux en particulier), il est proposé de développer une méthodologie macroscopique pour le comportement de la protection biologique en béton. Cette approche inclura différents phénomènes, parmi lesquels l'expansion volumétrique induite par le rayonnement, le fluage induit, les déformations thermiques et le chargement mécanique.
Les outils seront développés dans le cadre de la mécanique de l'endommagement. Les principaux défis numériques concernent la proposition et l'implémentation de lois d'évolution adaptées, et en particulier le couplage entre l'endommagement microstructural et l'endommagement au niveau structurel dû aux contraintes appliquées sur la structure.
Ce travail numérique pourra être réalisé dans un contexte de collaboration internationale. Il permettra au candidat retenu de développer un ensemble de compétences autour de la simulation de structures en béton armé en environnement complexe.
Impact des paramètres d’irradiation sur la formation de la phase alpha’ dans les aciers renforcés par dispersion d’oxydes (ODS)
Les aciers ferritiques-martensitiques renforcés par dispersion d'oxydes (aciers ODS) sont des matériaux d’intérêt pour la filière nucléaire. Composés majoritairement de fer et de chrome, ces matériaux peuvent être fragilisés par la précipitation sous irradiation d’une phase riche en chrome, la phase alpha prime. Cette phase, réputée sensible aux conditions d’irradiation, en fait un sujet idéal pour mieux comprendre les limites de la transférabilité ions-neutrons. En effet, si les irradiations aux ions sont fréquemment utilisées pour comprendre les phénomènes observés sous irradiation neutronique, la question de leur représentativité est régulièrement soulevée.
Dans cette thèse, nous cherchons donc à comprendre dans quelle mesure les paramètres des irradiations impactent les caractéristiques de la phase alpha’ dans les aciers ODS. Pour cela, des aciers ODS seront irradiés dans différentes conditions (flux, dose, température type de particules (ions, neutrons, électrons)) puis analysés à l’échelle nanométrique. Les caractéristiques des nano-oxydes (taille, densité) et de la phase alpha’ (taille, teneur en Cr), obtenues pour chacune des conditions d’irradiation, seront comparées à celles d’un échantillon de MA957 après irradiation aux neutrons.
Simulation numérique de modèles de turbulence sur des maillages déformés
La turbulence joue un rôle important dans de nombreuses applications industrielles (écoulement, transfert de chaleur, réactions chimiques). Comme la simulation directe (DNS) est souvent d’un coût excessif en temps calcul, les modèles en moyenne de Reynolds (RANS) sont alors utilisés dans les codes de CFD (computational fluid dynamics). Le plus connu, qui a été publié dans les années 70, est le modèle k – e. Il se traduit par deux équations additionnelles non-linéaires couplées aux équations de Navier-Stokes, décrivant le transport, pour l’une, de l’énergie cinétique turbulente (k) et, pour l’autre, de son taux de dissipation (e). Une propriété très importante à vérifier est la positivité des paramètres k et e qui est nécessaire pour que le système d’équations modélisant la turbulence reste stable. Il est donc crucial que la discrétisation de ces modèles préserve la monotonie. Les équations étant de type convection-diffusion, il est bien connu qu’avec des schémas classiques linéaires (Eléments finis, Volumes finis etc…), les solutions numériques sont susceptibles d’osciller sur des mailles déformées. Les valeurs négatives des paramètres k et e sont alors à l’origine de l’arrêt de la simulation. Il s’avère donc nécessaire de rendre monotone les schémas linéaires classiques de la littérature de manière consistante et stable.
Nous nous intéressons aux méthodes non linéaires permettant d’obtenir des stencils compacts. Pour des opérateurs de diffusion, elles reposent sur des combinaisons non linéaires de flux de part et d’autre de chaque arête. Ces approches ont montré leur efficacité, particulièrement pour la suppression d’oscillations sur des maillages très déformés. On pourra également reprendre les idées proposées dans la littérature, où il est par exemple décrit des corrections non linéaires s’appliquant sur des schémas linéaires classiques.
L’idée serait donc d’appliquer ce type de méthode sur les opérateurs diffusifs apparaissant dans les modèles k-e. Dans ce contexte, il sera également intéressant de transformer des schémas classiques de la littérature approchant les gradients en schémas non linéaires à deux points. Des questions fondamentales doivent être examinées dans le cas de maillages généraux à propos de la consistance et de la coercivité des schémas étudiés.
Au cours de la thèse, on prendra le temps de régler les problèmes de fond de ces méthodes (première et seconde année), à la fois sur les aspects théoriques et sur la mise en œuvre informatique. Cette dernière pourra être effectuée dans les environnements de développement Castem, Trust/TrioCFD ou POLYMAC. On s’intéressera alors à des solutions analytiques régulières et aux cas d’application représentatifs de la communauté.
Thermorégulation diphasique pour les composants semiconducteurs Ultra Grand Gap en diamant
Cette thèse porte sur l’étude d’un système de thermorégulation diphasique pour composants à semiconducteur de puissance ultra grand gap en diamant. Les composants en diamant ayant la particularité d’avoir leur résistance à l’état passant qui diminue lorsque la température augmente, cette thermorégulation vise à optimiser les pertes globales du système ainsi que d’assurer un équilibrage des températures ainsi que des contraintes entre plusieurs composants en diamant en parallèle.
Basé sur un cahier des charges qui sera défini en début de thèse (calories à évacuer, plage de température à réguler), le doctorant aura pour objectif de :
- définir une stratégie de contrôle de la température
- définir un couple matériau / fluide le mieux approprié
- Dessiner le système de thermorégulation
- Mettre en œuvre et valider expérimentalement le système proposé
La thèse abordera des aspects de simulations numérique (modélisation de composants et du système de thermorégulation) ainsi que des essais expérimentaux via la réalisation d’un prototype de TRL3-4 intégré à un système de convertisseur intégrant des diodes en diamant.
L’objectif final est de pouvoir mettre en avant un système innovant modélisé et démontré expérimentalement, où la stratégie de contrôle et les éléments dimensionnels et opératoires dimensionnants auront été investigués et optimisés.
Influence de la démontabilité des systèmes batteries sur leurs impacts environnementaux
Avec l'essor de la mobilité électrique et du stockage d'énergie, la demande en batteries explose. Mais cette croissance soulève une question essentielle : comment concevoir des batteries à la fois performantes, durables et plus respectueuses de l'environnement ?
Sans regarder la chimie de la cellule, l'une des pistes prometteuses est la démontabilité : rendre les packs batteries plus faciles à démonter permettrait de faciliter leur réparation, réutilisation ou recyclage. Cependant, un design plus démontable peut aussi augmenter sa masse ou réduire la fiabilité du système, et donc impacter sa durée de vie globale.
Cette thèse propose de relever ce défi en développant une méthode d'analyse capable de relier la conception de batteries démontables à leurs impacts environnementaux réels, tout en intégrant la question de la fiabilité. La ou le doctorant(e) évaluera la démontabilité de différents systèmes batteries, quantifiera les gains et pertes environnementaux associés, et contribuera à la mise au point d'un outil d'aide à la décision pour guider les choix de conception. Le travail proposé fera appel, entre autres, à la modélisation sous logiciel d'Analyse du Cycle de Vie (ACV) couplée à des modèles de performance et de vieillissement de batterie et à des probabilités de défaillance.
Ce travail s'inscrit dans un contexte technologique marqué par la recherche de circularité des ressources, les enjeux d'automatisation du démontage, et les nouvelles réglementations européennes sur les batteries. Il offrira une occasion unique de contribuer à la conception des batteries de demain, plus durables.
Génération assistée de noyaux de calculs complexes en mécanique du solide
Les lois de comportement utilisées dans les simulations numériques décrivent les caractéristiques physiques des matériaux simulés. À mesure que notre compréhension de ces matériaux évolue, la complexité de ces lois augmente.L'intégration de ces lois constitue une étape critique pour la performance et la robustesse des calculs scientifiques. De ce fait, cette étape peut conduire à des développements intrusifs et complexes dans le code.
De nombreuses plateformes numériques telles que FEniCS, FireDrake, FreeFEM, Comsol, proposent des techniques de génération de code à la volée (JIT, pour Just In Time) pour gérer différentes physiques. Cette approche JIT réduit considérablement les temps de mise en oeuvre de nouvelles simulations, offrant ainsi une grande versatilité à l'utilisateur. De plus, elle permet une optimisation spécifique aux cas traités et facilite le portage sur diverses architectures (CPU ou GPU). Enfin, cette approche permet de masquer les détails d'implémentation: une évolution de ces derniers est invisible pour l'utilisateur et est absorbée par la couche de génération de code.
Cependant, ces techniques sont généralement limitées aux étapes d'assemblage des systèmes linéaires à résoudre et n'incluent pas l'étape cruciale d'intégration des lois de comportement.
S'inspirant de l'expérience réussie du projet open-source mgis.fenics [1], cette thèse vise à développer une solution de génération de code à la volée dédiée au code de mécanique des structures de nouvelle génération Manta [2] développé par le CEA. L'objectif est de permettre un couplage fort avec les lois de comportement générées par MFront [3], améliorant ainsi la flexibilité, les performances et la robustesse des simulations numériques.
Le doctorant recherché devra posséder une solide culture numérique et un goût prononcé pour la simulation numérique et la programmation en C++. Il devra faire preuve d’autonomie et être force de proposition. Le doctorant bénéficiera d'un encadrement de la part des développeurs des codes MFront et Manta (CEA), ainsi que des développeurs du code A-Set (collaboration entre Mines-Paris Tech, Onera, et Safran). Cette collaboration au sein d'une équipe multidisciplinaire offrira un environnement stimulant et enrichissant pour le candidat.
De plus, le travail de thèse sera valorisé par la possibilité de participer à des conférences et de publier des articles dans des revues scientifiques à comité de lecture, offrant une visibilité nationale et internationale aux résultats de la thèse.
Le doctorat se déroulera au CEA Cadarache, dans le sud est de la France, au sein du département d'études des combustibles nucléaires de l'Institut REcherche sur les Systèmes Nucléaires pour la production d'Energie bas carbone (IRESNE) [4]. Le laboratoire d'accueil est le LMPC dont le rôle est de contribuer au développement des composants physiques de la plateforme numérique PLEIADES [5], co-développée par le CEA et EDF.
[1] https://thelfer.github.io/mgis/web/mgis_fenics.html
[2] MANTA : un code HPC généraliste pour la simulation de problèmes complexes en mécanique. https://hal.science/hal-03688160
[3] https://thelfer.github.io/tfel/web/index.html
[4] https://www.cea.fr/energies/iresne/Pages/Accueil.aspx
[5] PLEIADES: A numerical framework dedicated to the multiphysics and multiscale nuclear fuel behavior simulation https://www.sciencedirect.com/science/article/pii/S0306454924002408