Effet de la gravité sur l’agitation au sein d’un écoulement turbulent à bulles en canal
La compréhension des écoulements diphasiques et du phénomène d’ébullition représente un enjeu majeur pour le Commissariat à l’Énergie Atomique et aux Énergies Alternatives (CEA) à la fois pour la conception et pour la sûreté des centrales nucléaires. Dans un Réacteur à Eau Pressurisée (REP), la chaleur dégagée par le combustible nucléaire est transférée à l’eau du circuit primaire. En situation accidentelle, l’eau du circuit primaire peut passer en régime d’ébullition nucléée, voire évoluer jusqu’à la crise d’ébullition. Si le phénomène d’ébullition fait l’objet de nombreuses études, la dynamique des bulles générées retient également une attention particulière au CEA. Cette thèse s’intéressera au couplage entre la turbulence générée par un écoulement cisaillé et l'agitation induite par les bulles. Son originalité réside dans l’étude de l’effet de la gravité, obtenue par l’inclinaison du canal, un paramètre susceptible de générer des régimes d’écoulement complexes.
Le travail, de nature expérimentale, s'appuiera sur le nouveau dispositif CARIBE du CEA Saclay. La mission du doctorant consistera à caractériser les différents régimes d’écoulement, puis à mener une étude détaillée de l’écoulement en mettant en place une métrologie spécifique (notamment Particule Image Velocymetry (PIV), anémométrie à film chaud, sondes optiques). Mené au sein du laboratoire LE2H, le projet bénéficiera d'une collaboration étroite avec le LDEL (CEA Saclay) et l’IMFT (Toulouse). Le doctorant évoluera dans un environnement dynamique avec d'autres doctorants et présentera ses travaux dans des conférences nationales et internationales.
Nous recherchons un(e) candidat(e) en mécanique des fluides avec un intérêt marqué pour l’expérimentation (stage M2 possible). Cette thèse offre l’opportunité de développer une expertise en instrumentation, analyse de données et écoulements diphasiques turbulents, des compétences très valorisées dans les secteurs de l'énergie, de l'industrie et de la recherche académique.
Étude des solidifications locales dans un Réacteur à Sels Fondus
Dans un Réacteur à Sels Fondus (RSF), le combustible nucléaire se présente sous forme de sel liquide à haute température, qui est son propre caloporteur. Certains transitoires accidentels (sur-refroidissement du sel, fuite) peuvent causer des solidifications locales d'une partie du sel combustible. Ces solidifications ont un impact sur l'écoulement du sel dans le cœur, ainsi que son comportement neutronique, et peuvent mener à des échauffements locaux importants de parois. Ces transitoires sont encore peu étudiés, alors qu'ils ont un impact majeur sur la sûreté et le design d'un RSF.
L'objectif de la thèse est d'étudier différents transitoires accidentels qui peuvent mener à des solidifications locales, et d'étudier l'impact neutronique, thermique et thermo-hydraulique de ces solidifications sur les divers composants du cœur. Ces analyses seront menées en utilisant des outils multi-physiques adaptés aux RSF, tel que le code de CFD TrioCFD, ses extensions neutroniques TRUST-NK et de transport réactif Scorpio, ainsi que le code de neutronique déterministe APOLLO3. Afin de concilier précision et temps de calcul, plusieurs modélisations pourront être envisagées en fonction des transitoires étudiés : modélisation 1D / 3D turbulente (RANS, LES) pour la thermo-hydraulique, diffusion / transport SPn / transport Sn pour la neutronique.
Modélisation numérique de la déchirure ductile sur de longues distances en vue de quantifier les marges des méthodes d’ingénierie
La prédiction des modes de ruine des structures métalliques est une étape essentielle de l’analyse de fonctionnement des composants industriels où des éléments mécaniques sont soumis à des sollicitations importantes (par exemple composants des centrales nucléaires, pipelines, éléments structurels d’aéronefs …). Pour procéder à de telles analyses, il est essentiel de simuler correctement le comportement d’un défaut en régime ductile, c’est-à-dire en présence d’importantes déformations plastiques avant et durant la propagation.
La simulation numérique prédictive de la déchirure ductile est encore une problématique scientifique et technique ouverte malgré des progrès importants réalisés ces dernières années. L’approche dite locale de la rupture, notamment le modèle de Gurson (et sa version modifiée GTN), est largement utilisée pour modéliser la déchirure ductile.Mais son utilisation présente des limites : temps de calcul importants, arrêt de simulation suite à la présence d‘éléments complétement endommagés dans le modèle et non-convergence du résultat lorsqu’on diminue la taille des mailles.
Cette thèse a pour but de faire évoluer le modèle de simulation de déchirure ductile utilisé au LISN, pour l'appliquer aux grandes propagations de fissures sur structures complexes. Et de comparer les résultats obtenus avec les méthodes d'ingénieries qui sont plus simples à mettre en œuvre.
Modélisation de la propagation de fissure en fatigue en présence de contraintes résiduelles – Amélioration de la méthode G-theta
Les contraintes résiduelles sont des champs de contraintes auto-équilibrées que l’on retrouve dans certains composants mécaniques en l’absence de chargement extérieur. Dues au soudage, par exemple, ces contraintes peuvent potentiellement avoir un effet sur le comportement de la structure et sur sa résistance à la rupture. Lorsque l’on doit justifier de l’intégrité d’un composant mécanique, dans le cadre d’une démonstration de sûreté dans le nucléaire, il est impératif de connaître précisément le rôle de ces champs de contrainte sur la résistance du composant. Dans le cas de la propagation de fissure en fatigue, pour modéliser avec précision tous les phénomènes en jeu (redistribution des contraintes, évolution de la plasticité, effet de fermeture), il sera nécessaire d’améliorer les outils numériques, comme les méthodes de maillage et propagation de fissure (AMR, X-FEM…) et l’interpolation de l’intégrale J en cas de fissure débouchantes (méthode Gtheta). La thèse comportera deux volets complémentaires : (a) le développement numérique visant l’amélioration de la méthode Gtheta dans Castem associée à une modélisation de la propagation de fissure en 3D avec AMR et (b) la poursuite des essais applicatifs de propagation de fissure en fatigue dans différentes configurations de contraintes résiduelles.
Effet de la radiolyse de l’eau sur le flux d’absorption d’hydrogène par les aciers inoxydables austénitiques en réacteur nucléaire à eau pressurisée
Dans les réacteurs nucléaires à eau pressurisée, les éléments constitutifs du cœur sont exposés à la fois phénomènes de corrosion en milieu primaire, de l’eau pressurisée sous 150 bar et 300 °C environ, et à un flux neutronique. Les aciers inoxydables du cœur subissent des dommages dus à la combinaison du bombardement neutronique et de la corrosion. De plus, la radiolyse de l’eau peut impacter les mécanismes et cinétiques de corrosion, la réactivité du milieu et a priori les mécanismes et cinétique d’absorption d’hydrogène par ces matériaux. Ce dernier point, encore inexploré, peut s’avérer problématique car l’hydrogène en solution solide dans l’acier peut conduire à la modification (et la dégradation) des propriétés mécaniques de l’acier ou induire une fissuration prématurée de la pièce. Les travaux pionniers développés dans cette thèse très expérimentale seront centrés autour de l’impact des phénomènes de radiolyse sur les mécanismes et cinétiques de corrosion et surtout de prise d’hydrogène d’un acier inoxydable 316L exposé au milieu primaire sous irradiation. L’hydrogène sera tracé par le deutérium, l’irradiation neutronique simulée par irradiation électronique sur accélérateurs de particules. Une cellule de perméation existante sera reconfigurée en un dispositif unique pour permettre de mesurer in operando par spectrométrie de masse le flux de perméation de deutérium à travers un échantillon exposé au milieu primaire simulé en conditions de radiolyse. La distribution de l’hydrogène dans le matériau, ainsi que la nature des couches d’oxydes formées, seront analysées finement à l’aide des techniques de pointe disponibles au CEA et dans les laboratoires partenaires. Le(a) doctorant(e) devra in fine (i) identifier les mécanismes en jeu (corrosion et entrée d’hydrogène), (ii) en estimer les cinétiques et (iii) modéliser l’évolution du flux d’hydrogène dans l’acier fonction de l’activité de la radiolyse.
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é.
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é.
Schéma décalé pour les équations de Navier-Stokes avec des mailles de forme générale
La simulation des équations de Navier-Stokes demande de disposer de méthodes numériques précises et robustes prenant en compte des opérateurs de diffusions, des termes de gradient et de convection. Les approches opérationnelles ont montré leur efficacité sur des simplexes. Cependant, dans certaines modélisations ou certains codes (TrioCFD, Flica5), il peut être utile d’améliorer localement la précision des solutions à l’aide d’un estimateur d’erreur ou bien de prendre en compte des mailles de forme générale.
Rappelons que nous nous intéressons ici à des schémas décalés. Cela signifie que la pression est calculée au centre des mailles et les vitesses sur les arêtes (ou les faces) du maillage. On obtient alors des méthodes naturellement précises à bas nombre de Mach.
De nouveaux schémas ont été présentés récemment dans ce contexte et ont montré leur robustesse et leur précision. Cependant, ces discrétisations peuvent être très coûteuses en place mémoire et en temps calcul en comparaison aux schémas MAC sur des maillages réguliers.
Nous nous intéressons aux méthodes de type « gradient ». Certaines d’entre elles reposent sur une formulation variationnelle avec des inconnues de pression aux centres des mailles et des inconnues de vecteur vitesse sur les arêtes (ou les faces) des cellules. Cette approche a montré son efficacité, particulièrement en termes de robustesse.
Notons également qu’un algorithme avec les mêmes degrés de libertés que les méthodes MAC a été proposé et donne des résultats prometteurs.
L’idée serait donc de combiner ces deux approches, à savoir la méthode « gradient » avec les mêmes degrés de libertés que les méthodes MAC. Dans un premier temps, on s’attachera à retrouver les schémas MAC sur les maillages réguliers. Des questions fondamentales doivent être examinées dans le cas de maillages généraux : stabilité, consistance, conditionnement du système à inverser, verrouillage numérique.
On pourra également essayer de retrouver les gains en précisions obtenus à l’aide des méthodes pour discrétiser les gradients de pression.
Au cours de la thèse, on prendra le temps de régler les problèmes de fond de cette méthode (première et seconde année), à la fois sur les aspects théoriques et sur la mise en oeuvre informatique. Cette dernière pourra être effectuée dans les environnements de développement Castem, Trust/TrioCFD, ou POLYMAC. On s’intéressera alors aux cas d’application représentatifs de la communauté.
Développement d’un modèle de chimie transport, sous radiolyse de l’eau, d’un combustible usé en stockage géologique profond
Le stockage direct des combustibles usés (CU) est une solution alternative à leur retraitement pour la gestion des déchets nucléaires. Ce stockage direct des CU en milieu géologique profond pose des défis scientifiques liés à la compréhension fine des processus de dissolution et de libération des radionucléides. Ce sujet de thèse propose de développer un modèle scientifique détaillé, capable de décrire les mécanismes physico-chimiques complexes, tels que la radiolyse de l'eau et l'interaction entre le combustible irradié et son environnement. L'objectif est de proposer une modélisation du transport réactif précise pour améliorer la prédiction à long terme des performances du stockage. En utilisant des approches mécanistes, cette thèse s’inscrit dans une démarche d’aller-retour entre modélisation et expérimentation, visant à affiner la compréhension des mécanismes d'altération et à valider les hypothèses par des données expérimentales. Finalement, en s'appuyant sur des modèles existants, comme le modèle opérationnel radiolytique, ce travail proposera des améliorations pour réduire les hypothèses simplificatrices actuelles. Le candidat contribuera à des enjeux industriels et sociétaux majeurs liés à la gestion des déchets nucléaires et permettra d'apporter des solutions aux enjeux de sûreté associés.
Transfert de chaleur par rayonnement : résolution numérique efficace de problèmes associés en milieu Beerien ou non pour les besoins de validation de modèles simplifiés
Cette proposition de recherche se place dans le cadre de l’étude, par le biais de la modélisation et de la simulation numérique, des transferts thermiques au sein d’un milieu hétérogène constitué de solides opaques et d’un fluide transparent ou semi-transparent. Les modes d’échange considérés sont le rayonnement et la conduction.
Suivant l’échelle considérée, la luminance (ou radiance) est solution de l’Equation de Transfert Radiatif (ETR - Radiative Transfer Equation). Dans sa forme classique, cette ETR décrit les phénomènes de transferts thermiques à l’échelle dite locale où les solides sont distinctement répartis dans le domaine, tandis que, à l’échelle mésoscopique d’un milieu homogène équivalent, la radiance est solution d’une ETR généralisée (ETR(G)) quand le milieu nobéit plus à la loi de Beer-Lambert. Dans notre contexte, nous nous intéresserons à la résolution numérique de cette ETR dans ces deux configurations avec in fine un couplage à la résolution d’une équation de conservation de l’énergie pour la température.
Dans le cadre de la résolution déterministe de l’ETR, une méthode usuelle de traitement de la variable angulaire de cette équation est la méthodes des ordonnées discrètes (Sn) qui utilise une quadrature pour la sphère unité. Lorsque l’on considère un milieu non-Beerien, la résolution de l’ETR(G) est un sujet de recherche très actuel où l’approche Monte-Carlo semble recevoir plus d’attention. Pour autant, on peut rapprocher cette ETR(G) de l’équation de transport généralisée telle que formulée dans le contexte du transport de particules et appliquer une méthode spectrale pour sa résolution déterministe Sn. C’est la piste poursuivie dans cette thèse.
Le cadre applicatif direct de ces travaux est l’étude par simulation numérique des accidents des Réacteurs nucléaires à neutrons thermiques refroidis à l’Eau Légère (REL). En effet, la modélisation des échanges par rayonnement est primordiale car, en cas de dénoyage du coeur et d’assèchement des gaines de combustible, c’est un mécanisme d’extraction de puissance qui devient rapidement, à mesure que la température augmente, important à prendre en compte, au même titre que la convection par le gaz (vapeur d’eau). Par ailleurs, cette thematique est aussi importante avec le renouveau du nucléaire par le biais de startups proposant des réacteurs calogènes de type High Temperature Reactor (HTR) refroidis par un gaz.
L’objectif de cette thèse est l’analyse et le développement d’une méthode numérique innovante et efficace de résolution de l’ETR(G) (dans un environnement de simulation numérique haute performance) couplée à la résolution de la conduction thermique. Du point de vue applicatif, une telle méthode permettrait de réaliser des calculs de référence pour la validation et la quantification du biais des modèles simplifiés mis en oeuvre dans des simulations d’ingénierie.
Un travail réussi dans le cadre de cette thèse permettra à l’étudiant de prétendre à un poste de recherche en simulation numérique haute performance de problèmes physiques complexes, par-delà la seule physique des réacteurs nucléaires.