Modélisation à l'échelle atomique de la ségrégation induite par l'irradiation dans les alliages austénitiques

Modelling of radiation induced segregation in austenitic Fe alloys at the atomistic level

• Auteur / Autrice : Jean-Baptiste Piochaud
• Direction : Charlotte Becquart, Christophe Domain
• Type : Thèse de doctorat
• Discipline(s) : Molécules et Matière Condensée
• Date : Soutenance le 09/01/2013
• Etablissement(s) : Lille 1
• Ecole(s) doctorale(s) : École doctorale Sciences de la matière, du rayonnement et de l'environnement (Villeneuve d'Ascq, Nord)
• Partenaire(s) de recherche : Laboratoire : UMET - Unité Matériaux et Transformations

Mots clés

Mots clés contrôlés

Défauts ponctuels

Ségrégation (métallurgie)

Monte-Carlo, Méthode de

Théorie de la fonctionnelle de densité

Propriétés magnétiques

Acier au chrome-nickel

Simulation par ordinateur

Acier inoxydable austénitique

Mots clés libres

Ségrégation induite par l'irradiation

Modèle d'interactions de paires

Résumé

Dans les réacteurs à eaux pressurisée, les internes de cuve sont soumis à d’intenses irradiations induisant le mécanisme de corrosion sous contrainte assistée par l’irradiation initié par le phénomène de ségrégation induite par l’irradiation (SIR). La SIR observée dans les aciers austénitiques 316 est modélisée à l’échelle atomique en considérant un alliage ternaire modèle Fe–10Ni–20Cr. Pour atteindre cet objectif, nous avons construit un modèle d’interactions de paires pour le système Fe-Ni-Cr afin de modéliser la SIR par méthodes Monte Carlo. Le modèle d’interactions de paires fut déduit à la fois des calculs DFT (Density Functional Theory) disponibles dans les systèmes pures cubiques faces centrées, mais aussi en effectuant des calculs DFT dans l’alliage Fe–10Ni–20Cr. Les énergies de formation calculées dans cet alliage modèle dépendent fortement de l’environnement local des défauts ponctuels. Nous avons pu ainsi établir une relation liant ces énergies au nombre et à la position des atomes de Ni et de Cr se trouvant à proximité des défauts ponctuels. Dans cette étude, nous montrons qu’un modèle d’interactions de paires uniquement basé sur les interactions entre éléments positionnés sur les sites du réseau (éléments d’alliage et lacune) est incapable de prendre en compte seul à la fois les aspects thermodynamiques et cinétiques de la SIR. Une estimation plus précise des barrières de migration est requise pour les espèces migrantes que celle proposée qui ne prend pas en compte l’environnement du point col. Ce travail montre ainsi qu’il est nécessaire de réaliser des calculs DFT des énergies de migration afin d’être en mesure de calibrer un modèle utilisable dans le cadre de simulations de Monte Carlo cinétique atomique. Nous montrons également que l’ajustement de notre modèle sur les données DFT obtenues pour les interstitiels est incompatible avec la modélisation de la SIR aux électrons.

In pressurized water reactors, under irradiation internal structures are subject of irradiation assisted stress corrosion cracking which is influenced by radiation induced segregation (RIS). In this work RIS of 316 stainless steels is modelled considering a model ternary Fe–10Ni–20Cr alloy. For this purpose we have built an Fe-Ni-Cr pair interaction model to simulate RIS at the atomistic level using an atomistic kinetic Monte Carlo approach. The pair interactions have been deduced from density functional theory (DFT) data available in the pure fcc systems but also from DFT calculations we have performed in the Fe–10Ni–20Cr target alloy. Point defect formation energies were calculated and found to depend strongly on the local environment of the defect. As a consequence, a rather good estimation of these energies can be obtained from the knowledge of the number and respective positions of the Ni and Cr atoms in the vicinity of the defect. This work shows that a model based only on interaction parameters between elements positioned in perfect lattice sites (solute atoms and vacancy) cannot capture alone both the thermodynamic and the kinetic aspect of RIS. A more accurate of estimating the barriers encountered by the diffusing species is required than the one used in our model, which has to depend on the saddle point environment. This study therefore shows thus the need to estimate point defect migration energies using the DFT approach to calibrate a model that can be used in the framework of atomic kinetic Monte Carlo simulations. We also found that the reproduction by our pair interaction model of DFT data for the self-interstitial atoms was found to be incompatible with the modelling of RIS under electron irradiation.