Vacancy Diffusion Across FeCrAl Alloy Composition Space for Accident-Tolerant Fuel Cladding

Résumé Technique : Diffusion des lacunes à travers l'espace de composition des alliages FeCrAl pour le gainage de combustibles tolérants aux accidents

Énoncé du problème
Les alliages fer-chrome-aluminium (FeCrAl) sont des candidats de premier plan pour le gainage de combustibles tolérants aux accidents (ATF) dans les réacteurs à eau légère (LWR) en raison de leur résistance supérieure à l'oxydation à haute température par rapport au Zircaloy. Cependant, la durée de vie en réacteur de ces gainages est régie par l'évolution microstructurale induite par le rayonnement, où le transport des défauts ponctuels — spécifiquement la diffusion des lacunes — est le mécanisme dominant. Bien que la formation de couches d'oxyde protectrices soit bien caractérisée, la physique fondamentale des défauts dans le FeCrAl sous irradiation reste mal comprise. Les traitements théoriques existants reposent souvent sur des environnements idéalisés ou compositionnellement simplifiés, ne parvenant pas à rendre compte des fluctuations chimiques locales complexes inhérentes aux alliages de solution solide concentrés. Ces environnements locaux modifient considérablement la liaison des lacunes, les voies de migration et les comportements collectifs tels que le regroupement (clustering) et le piégeage. Il existe un besoin critique de démêler les contributions individuelles des éléments constituants (Fe, Cr, Al) au transport global des lacunes afin de prédire l'évolution matérielle à long terme, la stabilité des phases et le fluage sous des charges thermiques et d'irradiation combinées.

Méthodologie
Les auteurs utilisent un cadre de Monte Carlo cinétique (KMC) résolu par espèce pour modéliser la diffusion des lacunes à travers une gamme systématique de compositions Fe/Cr (de riche en Fe à riche en Cr) avec une concentration fixe de 5 at.% d'Al. La méthodologie procède selon les étapes suivantes :

1. Potentiel et génération de données : Les interactions interatomiques sont évaluées à l'aide d'un potentiel de la méthode des atomes enchâssés (EAM) ajusté à une base de données de la théorie de la fonctionnelle de la densité (DFT). Les données de migration des lacunes sont générées via le flux de travail « Hop-Decorate », où une structure de référence de Fe BCC est décorée aléatoirement avec des atomes de Fe, Cr et Al selon les compositions globales cibles. Les chemins de plus basse énergie sont calculés à l'aide de la méthode Climbing-Image Nudged Elastic Band (CI-NEB) pour les compositions globales de Fe$_{80}$Cr$_{15}$Al$_5$, Fe$_{47,5}$Cr$_{47,5}$Al$_5$ et Fe$_{15}$Cr$_{80}$Al_5}.
2. Modélisation de substitution (Surrogate Modeling) : Pour éviter le coût computationnel des calculs NEB en temps réel pendant la KMC, des modèles de substitution linéaires sont entraînés sur la base de données de barrières générée. L'environnement chimique local pour chaque saut est codé à l'aide de descripteurs comptant les atomes de Fe, Cr et Al dans des couches radiales fixes autour du point de selle et de l'atome migrant. Des modèles de régression linéaire distincts sont ajustés pour chaque espèce migrante (Fe, Cr, Al) afin de capturer les dépendances de barrière spécifiques à l'espèce.
3. Simulations KMC : Les modèles de substitution entraînés pilotent des simulations KMC à temps de résidence à des températures allant de 500 K à 1500 K. Les simulations suivent les trajectoires de lacunes uniques sur $10^5$ étapes pour calculer le déplacement quadratique moyen (MSD) et extraire les coefficients de diffusion ($D$), les énergies d'activation ($E_a$) et les facteurs pré-exponentiels ($D_0$).

Résultats clés

• Barrières de migration et environnement local : L'étude confirme que les barrières de migration sont hautement sensibles à l'environnement chimique local. Alors que le Fe pur présente une barrière d'environ 0,63 eV, le Cr présente une barrière de ~0,95 eV et l'Al de ~0,69 eV. Dans les alliages, les distributions de barrières s'élargissent considérablement, particulièrement pour le Cr et l'Al, indiquant que l'asymétrie chimique locale module les voies de migration.
• Rôle de l'aluminium : Malgré sa faible concentration (5 at.%), l'Al exerce une influence disproportionnée sur les barrières de migration. Dans les trois modèles de substitution, la présence d'un atome d'Al au premier voisin (1NN) de la coquille de selle augmente considérablement la barrière énergétique (+0,048 eV pour le Fe, +0,073 eV pour le Cr et +0,151 eV pour l'Al).
• Dépendance de la diffusivité par la composition :
  • De riche en Fe à équiatomique : Dans les compositions allant de Fe$_{80}$Cr$_{15}$Al$_5$ à Fe$_{32,5}$Cr$_{62,5}$Al$_5$, la diffusivité des lacunes est relativement insensible à la composition globale. Les énergies d'activation convergent vers une plage de 0,70–0,84 eV, suggérant que la migration du Fe domine le transport dans ces régimes.
  • Régime riche en Cr : La composition riche en Cr (Fe$_{15}$Cr$_{80}$Al$_5$) présente une énergie d'activation nettement élevée de 1,135 eV. Par conséquent, sa diffusivité macroscopique est de plusieurs ordres de grandeur inférieure à celle des alliages riches en Fe à basses températures (par exemple, ~$10^{-17}$ m$^2$ s$^{-1}$ à ~570 K contre ~$10^{-13}$ m$^2$ s$^{-1}$ pour les alliages riches en Fe).
• Percolation et diffusion « lente » : Les résultats suggèrent un seuil de percolation où le réseau connecté de chemins à faible barrière (dominé par les sauts de Fe) se rompt à mesure que la teneur en Cr augmente. Dans le régime riche en Cr, la migration des lacunes est dominée par les sauts de Cr à haute barrière, conduisant à une « diffusion lente ». Les auteurs observent que les trajectoires des lacunes dans des paysages chimiquement rugueux impliquent une échappée de poches de piégeage, entraînant un comportement non linéaire du déplacement quadratique moyen avant d'atteindre un régime diffusif.

Signification et affirmations
L'article affirme que la diffusivité macroscopique des alliages FeCrAl ne peut être prédite par une simple moyenne pondérée par la composition des barrières des composants purs ; elle est plutôt déterminée par l'hétérogénéité du paysage énergétique local et la connectivité spécifique des chemins de migration à faible barrière.

• Implications de conception : L'étude identifie un régime compositionnel distinct (riche en Cr) où la mobilité des lacunes est supprimée. Cette suppression implique que dans les phases $\alpha'$ riches en Cr (connues pour se former sous irradiation), les lacunes sont moins mobiles, ce qui pourrait réduire le taux de recombinaison lacune-interstitiel et altérer la cinétique d'accumulation des défauts.
• Tolérance aux accidents : Sous des conditions de perte de refroidissement (LOCA) où les températures dépassent 1200 K, la dépendance de la diffusivité par la composition devient prononcée. Les auteurs suggèrent que la composition riche en Cr, avec son énergie d'activation élevée et sa diffusivité supprimée, peut offrir une cinétique de lacunes favorable pour le gainage tolérant aux accidents, complétant potentiellement sa résistance connue à l'oxydation.
• Contribution méthodologique : Ce travail démontre que les cadres de Monte Carlo cinétique à ordre réduit, paramétrés par des modèles de substitution linéaires résolus par espèce et entraînés sur des bases de données NEB, peuvent efficacement cartographier le transport des défauts à travers des espaces de composition complexes. Cette approche fournit un lien mécaniste entre les processus atomistiques et les modèles mésoscopiques utilisés dans les codes de performance des combustibles, allant au-delà de l'approche empirique par essais et erreurs de la conception d'alliages.

Limites notées par les auteurs
Les auteurs notent explicitement plusieurs limites : le modèle traite le réseau comme statique et gelé, ignorant la relaxation locale et l'évolution de la composition induite par le rayonnement (ségrégation) au fil du temps. Le modèle de substitution linéaire ne peut capturer les effets directionnels ou de rupture de symétrie, ni les contributions magnétiques au paysage énergétique. De plus, l'étude se concentre uniquement sur le transport des lacunes et n'aborde pas le transport des amas d'interstitiels ni la transformation des boucles de dislocation induite par le rayonnement, qui peuvent également influencer les densités de puits de défauts.