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é.
Vers une méthode de caractérisation des propriétés électrocinétiques de particules dans l’eau à haute température
Dans le domaine de l’industrie et notamment de l’énergie, les circuits en eau liquide sont omniprésents. Les fluides, en interagissant avec les tuyauteries à base d’alliages métalliques, conduisent inévitablement à la formation de produits de corrosion. Des particules ainsi formées se déplacent dans les circuits sous l’effet d’un écoulement. En fonction des propriétés surfaciques physiques ou chimiques des parois, du milieu et des particules elles-mêmes, ces dernières peuvent s’agréger, se disperser, s’adsorber ou se déposer dans d’autres parties du circuit et conduire, par exemple, à des phénomènes d’encrassement et in fine à la perte de rendement des procédés industriels.
La prédiction du comportement des particules de petite taille (ordre de grandeur du µm) revêt donc un intérêt particulier. En effet, de par leur dimension, le comportement de ces dernières est régi par des forces d'origine électrique responsables de leur adhésion sur les surfaces. Les propriétés électrocinétiques et notamment le potentiel de surface pilotent ainsi le devenir de la particule et peuvent être définies par le biais du potentiel zêta. Cette grandeur caractérise un couple solide/solution et prend en considération à la fois la particule et ses propriétés chimiques de surface ainsi que la solution dans laquelle se trouve la particule.
Si la caractérisation du potentiel zêta à température ambiante est assez répandue, sa détermination à haute température se cantonne aujourd’hui, à quelques exemples (thèses de C. Cherpin 2022 [1] et de M. Barale 2006 [2], les études de VTT [3] et celle d’EDF avec l’université de Besançon 2002 [4] ainsi que le brevet de l’EPRI 1994 [5]). Le CEA (LC2R) a développé un moyen de mesure innovant en cours de brevetage pour explorer des techniques expérimentales peu développées et basées sur des hypothèses théoriques à approfondir.
A travers des approches multi-physique (écoulement, température, chimie, électrochimie, etc.) et multi-échelle (particules microscopiques influant sur un état macroscopique), l’objectif de la thèse est donc de réaliser les mesures des propriétés de surface de particules dans l’eau à haute température en fonction des conditions physico-chimiques (pH, RedOx et température), d’adapter les modèles existants ou en proposer de nouveaux puis de les valider avec les données expérimentales.
Les données ainsi obtenues ont vocation à alimenter les codes de simulation afin de mieux appréhender et maîtriser le vieillissement des circuits.
[1] C. Cherpin, PhD, 2022, Modelling the behaviour of colloidal corrosion products in the primary circuit of Pressurized Water Reactors
[2] M. Barale, PhD, 2006, Etude du comportement des particules colloïdales dans les conditions physico-chimiques du circuit primaire des réacteurs à eau sous pression
[3] E. Velin, Master’s Thesis, 2013, The effect of Temperature on the Zeta Potential of Magnetite Particles in Ammonia, Morpholine and Ethanolamine Solutions
Apprentissage de Modèles Interprétables pour la Corrosion sous Contrainte des aciers inoxydables exposés en milieu primaire des REP
La corrosion sous contrainte (CSC) des aciers inoxydables est l'un des principaux phénomènes de dégradation des composants du circuit primaire des Réacteurs à Eau Pressurisée(REP). La compréhension de ce mécanisme de fissuration est d’une nécessité absolue pour la prolongation de la durée d’exploitation des réacteurs. Avec un nombre important de paramètres critiques qui influent sur la sensibilité du matériau à la CSC et la présence de forts effets de couplage, une grille d’essais expérimentaux assez conséquente est souvent envisagée pour aider à la compréhension du mécanisme. Il est proposé dans ce projet d’adopter une approche nouvelle basée sur l’utilisation de modèles interprétables, avec pour but d’éviter les longues et couteuses étapes de recherches en ciblant des essais pertinents et des paramètres matériaux pouvant améliorer les performances en environnement. L’enjeu ici sera d’ajouter à l’approche expérimentale les performances d’un outil d’intelligence artificielle avec pour objectifs de définir des domaines de sensibilité à l’amorçage de CSC en fonction des paramètres critiques identifiés dans le modèle, et de fournir des données relatives à l’élaboration de nouveaux matériaux par fabrication additive.
La thèse sera consacrée au développement d’un outil numérique adapté à ce nouveau cas d’usage et à la poursuite des activités expérimentales nécessaires à la validation de cette nouvelle approche. Il s’agira d’explorer les contributions de l’intelligence artificielle dans le domaine de la corrosion sous contrainte sur plusieurs volets : l'identification des paramètres au premier ordre sur la sensibilité du matériau, l'évaluation des domaines de criticité à la CSC et l'aide à la compréhension des mécanismes physiques à l’origine de la fissuration.
Elaboration, caractérisation et modélisation de films minces à base d’oxydes (Mn,Co)3O4 appliquées aux revêtements contre la corrosion et à la spintronique
Les spinelles de métaux de transition apparaissent spontanément lors de la corrosion généralisée des aciers ou alliages en milieu aqueux ou gazeux à haute température. Ces phases spinelles de type AB2O4 forment une couche de corrosion continue et régissent de ce fait les processus de corrosion car elles régulent la conductivité et le transport de matière entre le matériau et le milieu. Ces spinelles sont aussi déposés volontairement comme revêtements de protection contre les phénomènes de dégradation. En particulier, le système spinelle Mn-Co-O est très prometteur en tant que couche conductrice protectrice sur l’acier inoxydable ferritique utilisé pour fabriquer des interconnexions dans les piles à combustibles à oxyde solide pour la production d’hydrogène vert. Le choix de la composition de ces phases spinelle détermine évidement les caractéristiques de protection des revêtements. Ces caractéristiques sont particulièrement délicates pour les matériaux des électrolyseurs à haute température car le transport électronique doit être optimal (électrolyse importante) mais ne doit pas s’accompagner de transport de matière (diffusion des cations faible).
Paradoxalement, les propriétés de transport des spinelles de métaux de transition sont en général mal connues. Les mesures sont faites sur des couches de corrosion ou des revêtements de compositions variables, de faibles cristallinités, de microstructures complexes et de surcroit de faibles épaisseurs. Par ailleurs, les spinelles montrent des propriétés magnétiques et de désordre cationique en fonction de la composition largement ignorés alors qu’elles ont un impact fort sur le transport. Ce sont précisément ces propriétés magnétiques et de transport qui présentent un intérêt majeur dans le domaine de la spintronique. Ainsi, la manipulation de la composition chimique de ces oxydes de structure spinelle (normale, inverse ou mixte) offre une large gamme de propriétés magnétiques (ferrimagnétique, antiferromagnétique) et électroniques (semi-métallique, semi-conducteur, isolant). En particulier CoMn2O4, est prédit avoir une configuration magnétique complexe [1], reliée principalement à l’arrangement des cations Co2+ et Mn3+ dans les sites interstitiels, qu’il convient d’analyser en détail. Ces études physiques requièrent, à l’inverse des couches de corrosion, la synthèse de couches minces de composition et de cristallinité bien maitrisées.
L’objectif de la thèse est d’une part d’apporter les connaissances structurales et physiques des oxydes modèles de composition chimique (Mn,Co)3O4 pour contribuer à l’élaboration de diagramme de phase (Mn-Co-O) et d’autre part de développer un modèle de transport électronique reposant sur la relation entre ordre/désordre – propriétés magnétiques et résistivité pour les spinelles (Mn,Co)3O4 et à terme sur l’ensemble des spinelles (Fe,Cr,Mn,Co)3O4. Les investigations seront conduites sur des couches minces de composition parfaitement maitrisées, de grande cristallinité et seront complétées par des simulations numériques. L’ensemble des travaux expérimentaux et de modélisation s’appuiera sur les résultats des études précédentes sur les couches simples de spinelles de composition (Ni,Fe,Cr)3O4 [2,3].
La thèse comportera plusieurs volets :
- Croissance de couches minces et multicouches par MBE (Molecular Beam Epitaxy) (J.-B. Moussy)
- Caractérisations spectroscopiques par XPS (X-ray photoemission spectroscopy) (F. Miserque)
- Caractérisations structurales fines par DRX et par absorption X (XMCD) (P. Vasconcelos)
- Modélisation des spectres XPS et d’absorption X, et modélisation atomistique (A. Chartier)
- Caractérisations magnétiques par magnétométrie SQUID/VSM et électriques (J.-B. Moussy)
[1] Systematic analysis of structural and magnetic properties of spinel CoB2O4 (B= Cr, Mn and Fe) compounds from their electronic structures, Debashish Das, Rajkumar Biswas and Subhradip Ghosh, Journal of Physics: Condensed Matter 28 (2016) 446001.
[2] Stoichiometry driven tuning of physical properties in epitaxial Fe3-xCrxO4 thin films, Pâmella Vasconcelos Borges Pinho, Alain Chartier, Denis Menut, Antoine Barbier, Myrtille O.J.Y. Hunault, Philippe Ohresser, Cécile Marcelot, Bénédicte Warot-Fonrose, Frédéric Miserque, Jean-Baptiste Moussy, Applied Surface Science 615 (2023) 156354.
[3] Elaboration, caractérisation et modélisation de films minces et multicouches à base d’oxydes (Ni,Fe,Cr)3O4 appliquées à la corrosion et à la spintronique, A. Simonnot, thèse en cours.
Etude de la catalyse de l’acide nitrique sur les aciers inoxydables
Le vieillissement des matériaux (principalement des aciers inoxydables) de l’usine de retraitement des combustibles nucléaires usés fait l’objet d’une importante activité de R & D au CEA. Le contrôle de ce vieillissement sera réalisé par une meilleure compréhension des mécanismes de corrosion des aciers inoxydables en acide nitrique (l'agent oxydant utilisé dans les étapes de retraitement).
L'objectif de la thèse est de développer un modèle de corrosion d’un acier inoxydable en fonction de la température et de la concentration en HNO3 via la quantification des produits de corrosion. Cette thèse représente un réel challenge technologique car actuellement peu d’études existent sur des mesures électrochimiques in situ dans l’acide nitrique chaud et concentré. Le doctorant réalisera également un travail expérimental important en couplant des mesures électrochimiques, des analyses chimiques (spectrométrie UV-visible-IR ...) et des analyses de surfaces (SEM, XPS,…). Sur la base de ces résultats expérimentaux, un modèle sera développé, qui sera incorporé à l'avenir dans un modèle plus global du vieillissement des équipements industriels de l'usine.
Le laboratoire est spécialisé dans l'étude de la corrosion dans des conditions extrêmes. Il est composé d'une équipe scientifique très dynamique et motivée qui a l'habitude de recevoir des étudiants.
Simulations mésoscopiques et développement de modèles simplifiés pour le comportement mécanique des bétons irradiés
Dans les centrales nucléaires, le puits de cuve en béton sert de support pour la cuve du réacteur et d’écran de protection contre les radiations. A long terme, l’exposition à des radiations neutroniques peut causer une expansion des granulats du béton par amorphisation, et provoquer une microfissuration et une dégradation de ses propriétés mécaniques. Cette problématique est importante dans les études visant à prolonger la durée de vie des centrales. À l’échelle mésoscopique, ces phénomènes peuvent être modélisés en dissociant le comportement des granulats, de la matrice cimentaire, et des interfaces entre ces phases. Cependant, il est difficile de décrire l’initiation et la propagation des microfissures dans de tels systèmes multi-fissurés hétérogènes complexes. L'objectif de cette thèse, menée dans le cadre d’un projet ANR franco-tchèque, est de développer un outil de simulation numérique performant pour analyser les effets de l’irradiation neutronique sur le béton à l’échelle mésoscopique. Une approche couplée thermo-hydro-mécanique dans laquelle le comportement de la matrice prendra en compte retrait, fluage et microfissuration sera utilisée. Les simulations seront validées à partir de données expérimentales obtenues sur des échantillons testés, et l’outil numérique permettra par la suite d’estimer l’impact de différents facteurs sur le comportement et les performances du béton soumis à une irradiation neutronique.
Ce projet de recherche s'adresse à un doctorant souhaitant développer ses compétences en science des matériaux, avec une forte composante en modélisation et simulations numériques multiphysiques et multi-échelles.
Etude du comportement en corrosion dans NaCl-MgCl2-CeCl3 d’un alliage base nickel en présence de produits de fission (Te,S) pour les réacteurs à sels fondus
L’accès à une énergie propre et peu coûteuse semble plus que jamais primordial dans le contexte actuel d’urgence climatique. Plusieurs pistes sont envisagées depuis plusieurs années déjà mais de nombreux verrous technologiques restent à lever pour les concrétiser, tant elles représentent des ruptures technologiques. Que ce soit pour le stockage d’énergie ou les réacteurs nucléaires de 4ème génération, le milieu sel fondu utilisé comme caloporteur et/ou comme combustible est fortement corrosif rendant le choix des matériaux de structure très complexe.
L’objectif du sujet de thèse proposé au sein du Service de Corrosion et du Comportement des Matériaux (S2CM) consiste en l’étude intégrale du comportement d’alliages base nickel prometteurs dans le ternaire NaCl-MgCl2-CeCl3, représentatif du sel utilisé dans le concept français de réacteurs à sels fondus, à 600°C. Par intégrale, il est ici entendu depuis la préparation d’éprouvette à la caractérisation multi-échelle et multi-techniques des produits de corrosion. Cette thématique revêt un haut caractère expérimental et de compréhension des mécanismes de corrosion. L’influence des produits de fission, tels que le tellure ou le soufre sur les mécanismes de corrosion sera particulièrement étudiée.
Prédire les propriétés thermodynamiques des défauts dans des alliages métalliques multi-composants à partir de l'échelle atomique par apprentissage statistique
Les propriétés et le comportement des matériaux dans des conditions extrêmes sont essentiels pour les systèmes énergétiques tels que les réacteurs de fission et de fusion. Cependant, prédire avec précision les propriétés des matériaux à haute température reste un défi. Les mesures directes de ces propriétés sont limitées par les instruments expérimentaux, et les simulations à l'échelle atomique basées sur des champs de force empiriques sont souvent peu fiables en raison d'un manque de précision. Ce problème peut être résolu à l'aide de techniques d'apprentissage statistique, qui ont récemment vu leur utilisation exploser en science des matériaux. Les champs de force construits par apprentissage statistique atteignent le degré de précision des calculs {it ab initio} ; cependant, leur mise en œuvre dans les méthodes d'échantillonnage est limitée par des coûts de calcul élevés, généralement supérieurs de plusieurs ordres de grandeur à ceux des champs de force traditionnels. Pour surmonter cette limitation, deux objectifs seront poursuivis dans cette thèse : (i) améliorer les champs de force par apprentissage statistique actif en trouvant un meilleur compromis précision-efficacité et (ii) créer des méthodes accélérées d'échantillonnage de l'énergie libre et des chemins cinétiques afin de faciliter l'utilisation de champs de force d'apprentissage statistique coûteux en termes de calcul. Pour le premier objectif, nous améliorons la construction des champs de force d'apprentissage statistique en nous concentrant sur trois facteurs clés : la base de données, le descripteur de l'environnement atomique local et le modèle de régression. Pour le deuxième objectif, nous mettrons en œuvre un schéma d'échantillonnage bayésien rapide et robuste pour estimer l'énergie libre anharmonique, qui est cruciale pour comprendre les effets de la température sur les solides cristallins, à l'aide d'une méthode de force de biais adaptative qui améliore considérablement la vitesse de convergence et la précision globale. Nous appliquerons les méthodes développées au calcul de l'énergie libre et de ses dérivées, des grandeurs physiques qui donnent accès aux propriétés thermo-élastiques des alliages et aux propriétés thermodynamiques des défauts ponctuels. Pour cela, nous utiliserons des extension algorithmiques qui permettent d'échantillonner un état métastable spécifique et aussi les chemins de transition vers d'autres bassins d'énergie et donc d'estimer les énergies libres de formation et de migration de défauts lacunaires. Les grandeurs thermodynamiques calculées seront ensuite utilisées comme données d'entrée de méthodes de Monte Carlo cinétique, qui permettra de mesurer les coefficients de diffusion dans les alliages complexes en fonction de la température. Un but sera d'essayer de relier les propriétés de transport atomique à la complexité de l'alliage. Notre approche ayant une rapidité considérablement supérieure à celle des méthodes standard, nous pourrons envisager de l'appliquer à des alliages complexes comprenant les éléments W, Ti, V, Mo et Ta à des températures et des compositions qui n'ont pas été étudiées expérimentalement.
Les super-réseaux pour la caractérisation de la diffusion sous irradiation à l’échelle atomique
Les alliages métalliques utilisés dans les applications nucléaires sont soumis à des températures relativement basses (inférieures à 450°C) pendant des temps importants (supérieurs à 10 ans). A ces températures, les cinétiques de transformation des microstructures contrôlées par la diffusion sont lentes. L’apparition de certaines phases indésirables, susceptibles de fragiliser le matériau, peut survenir après plusieurs années de service. Les coefficients de diffusion jouent donc un rôle crucial en tant que données d'entrée pour modéliser l'évolution de ces microstructures à l’aide de modèles phénoménologiques. Or, la détermination expérimentale des coefficients de diffusion à basse température (T < 600°C) est extrêmement délicate, notamment en raison de la nécessité de caractériser des longueurs de diffusion nanométriques, une difficulté accrue en présence d'irradiation.
Avec le développement de l’analyse chimique en microscopie électronique en transmission (MET) et de la sonde atomique tomographique (SAT), il est désormais possible d’accéder expérimentalement à de très faibles longueurs de diffusion et donc de déterminer des coefficients de diffusion à basse température à l’aide de super-réseaux, qui sont des empilements de couches nanométriques de compositions chimiques différentes. On peut même caractériser l’effet de l’irradiation sur la diffusion en réalisant des irradiations aux ions, permettant de simuler les modifications causées par l’irradiation neutronique sans activer les matériaux. L’objectif de la thèse porte sur le développement d’une méthodologie et la caractérisation de la diffusion hors et sous irradiation dans un système ternaire d’intérêt (Ni–Cr–Fe), représentatif des aciers et des alliages à haute entropie envisagés dans l’industrie nucléaire.
Ce sujet de thèse est une opportunité de travailler avec des techniques expérimentales de pointe, en étroite collaboration avec une équipe de théoriciens du même département, ainsi qu’avec des équipes spécialisées dans l’élaboration de super-réseaux de l’UTBM à Belfort et du CINAM à Marseille.
Comportement de nanocavités sous chargement mécanique : de la compréhension des mécanismes physiques à l’homogénéisation de matériaux nanoporeux
Des nanocavités - typiquement de quelques nm à quelques dizaines de nm - sont souvent observées dans les métaux, par exemple dans les applications hautes températures suite à la condensation de lacunes ou dans les alliages métalliques utilisés dans les réacteurs nucléaires du fait de l’irradiation. La présence de ces nanocavités dégrade le comportement mécanique des matériaux et contribue à la rupture. Il est donc nécessaire de déterminer les mécanismes physiques associés au comportement de ces nanocavités sous chargement mécanique et d’obtenir des modèles homogénéisés décrivant le comportement macroscopique des matériaux nanoporeux. Les résultats disponibles dans la littérature restent à ce jour limités, notamment en ce qui concerne la représentativité des simulations réalisées et des modèles proposés vis-à-vis des applications d’intérêt. Cela inclut par exemple la prise en compte des défauts cristallins entourant les cavités, l’effet des chargements cycliques et la localisation des nanocavités au niveau des joints de grains. Les objectifs de cette thèse sont donc de déterminer le comportement de nanocavités sous chargement mécanique et les mécanismes physiques associés en considérant des situations réalistes vis-à-vis des applications, de développer des modèles analytiques à bases physiques permettant de décrire le comportement de nanocavités sous chargement mécanique, et enfin de proposer des modèles homogénéisés adaptés aux nanocavités et utilisables pour simuler la rupture par croissance et coalescence de cavités. Les cas d’applications visés sont ceux liés aux alliages métalliques sous irradiation mais les éléments de compréhension obtenus et les modèles développés pourront être utilisés dans un contexte plus large. Afin d’atteindre ces objectifs, des simulations de Dynamique Moléculaire (DM) seront réalisées, analysées à partir de la théorie élastique des dislocations et utilisées pour proposer des modèles homogénéisés pertinents pour les matériaux nanoporeux.