Introduction à l'oxydation dans les matériaux à haute température

Les matériaux à haute température constituent l'épine dorsale des systèmes modernes de propulsion aérospatiale, des turbines à gaz, des réacteurs nucléaires et des fours industriels.Ces environnements exposent les alliages, les céramiques et les composites à des charges thermiques extrêmes tout en les défiant simultanément avec des gaz corrosifs, notamment l'oxygène. L'oxydation à des températures élevées est un mécanisme de dégradation primaire. Lorsqu'une surface métallique ou céramique réagit avec l'oxygène, elle forme une échelle d'oxyde. Idéalement, cette échelle est protectrice, ralentissant davantage les attaques.

La caractérisation expérimentale de l'oxydation à haute température est notoirement difficile. Les vitesses de réaction sont rapides, les distances de diffusion sont courtes et l'observation in situ des réarrangements atomiques reste difficile même avec la microscopie électronique avancée. C'est là que la modélisation computationnelle, en particulier la Théorie fonctionnelle de densité (DFT), devient indispensable. La DFT permet aux chercheurs de simuler les interactions électroniques et atomiques fondamentales qui régissent l'oxydation, fournissant un niveau de détail inaccessible par expérience seule.

Cet article présente un examen détaillé et faisant autorité de la façon dont DFT est appliqué pour simuler l'oxydation dans les matériaux à haute température. Nous aborderons les principes fondamentaux de DFT, les étapes spécifiques de la modélisation des réactions d'oxydation, l'analyse de la formation de couches d'oxyde, les forces et les limites de la méthode, et les directions futures passionnantes qui combinent DFT et apprentissage machine pour accélérer la découverte.

Théorie fonctionnelle de densité: Primer pour l'oxydation des matériaux

La théorie fonctionnelle de la densité est un cadre mécanique quantique qui remplace la fonction d'onde multi-électrons par la densité d'électrons comme variable centrale. Cette approche réduit considérablement la complexité computationnelle tout en conservant une grande précision pour les propriétés du sol.Dans le contexte de l'oxydation, DFT calcule les énergies et les forces électroniques pour un système d'atomes représentant la surface matérielle et adsorbant les espèces d'oxygène.

Les substrats typiques comprennent les superalliages à base de nickel, les alumines de titane, les métaux réfractaires (p. ex. molybdène, tungstène) et les composés céramiques comme le carbure de silicium ou l'alumine. La dalle est généralement de 3 à 8 couches atomiques d'épaisseur avec une région de vide au-dessus pour isoler la surface. Les molécules d'oxygène (O2) ou l'oxygène atomique sont ensuite placés à divers sites d'adsorption. L'optimisation de la géométrie DFT détend les positions atomiques pour trouver des configurations d'énergie minimales. Les énergies résultantes permettent le calcul des énergies d'adsorption : Eads = E(slab+O) – E(slab) – 1⁄2 E(O2). Une énergie d'adsorption fortement négative indique une réaction favorable.

Les recherches sur l'état de transition (par exemple, en utilisant la méthode de bande élastique nudguée) révèlent la barrière énergétique pour la dissociation de l'oxygène, étape limitante dans de nombreuses séquences d'oxydation. Les barrières faibles signifient la formation rapide d'oxyde, qui peut être bénéfique (former une échelle de protection dense) ou préjudiciable (consommer le métal trop rapidement).

Il est important de noter que la DFT n'est pas sans approximations. Le choix de la fonction de corrélation d'échange – comme PBE, PW91 ou les fonctions plus récentes de SCAN et de meta-GGA – affecte la précision. Pour les métaux de transition et leurs oxydes, les fonctions standard sous-estiment souvent les écarts de bande et peuvent ne pas saisir les effets de corrélation d'électrons forts dans certains systèmes d- et f-électron. Néanmoins, pour la majorité des problèmes d'oxydation à haute température, la DFT fournit des tendances fiables et des énergies relatives, ce qui en fait un outil de travail pour la science des matériaux informatiques.

Simulation des réactions d'oxydation sur les surfaces en alliages à haute température

Modélisation de l'adsorption et de la dissociation de l'oxygène

La première étape de l'oxydation est la physisorption et la chimisorption subséquente des molécules d'oxygène sur la surface métallique propre. Les études DFT sur le nickel (Ni), élément de base de nombreux superalliages, montrent que les adsorbages O2 sont préférentiels sur les sites de pont et creux avec des énergies autour de -0,5 à -1,0 eV par atome d'O. La molécule se dissocie ensuite en atomes d'oxygène séparés, une réaction qui est très exothermique et presque sans barrière sur la plupart des surfaces métalliques de transition.

Pour les alliages, la composition de la surface peut différer de celle du vrac en raison de la ségrégation. Le DFT peut expliquer cette situation en modélisant les configurations de dalles avec un ou plusieurs éléments d'alliage (p. ex., Cr, Al, Mo, W) à la surface. Le chrome et l'aluminium sont des anciens oxydes classiques; leur présence réduit l'énergie d'adsorption de l'oxygène et favorise la formation rapide d'échelles Cr2O3 ou Al2O3. DFT prédit que les atomes d'oxygène se lient plus fortement sur ces éléments réactifs que sur le nickel ou le cobalt, ce qui déclenche la nucléation locale de l'oxyde.

Un exemple représentatif est l'oxydation des alliages Ni-Al, la base de nombreux couches de liaison. Les calculs DFT par des chercheurs à [Université de Californie, Santa Barbara](https://www.mrl.ucsb.edu/) ont montré que la barrière d'activation de la diffusion de l'oxygène par une couche d'alumine mince est exceptionnellement élevée (>3 eV), expliquant pourquoi les échelles d'alumine fournissent une protection aussi efficace. Inversement, dans les alliages Ni-Cr, la diffusion de l'oxygène par la chromie est plus rapide (barriers ~2 eV), ce qui conduit à une cinétique d'oxydation parabolique bien prédite par les coefficients dérivés de DFT.

Diffusion et nucléation d'oxydes de surface

Une fois les atomes d'oxygène chémisorbés, ils peuvent pénétrer dans la région subsurface, ce qui est essentiel pour l'oxydation interne, qui peut fragiliser un alliage. Les simulations DFT des sites interstitiaux d'oxygène dans le nickel et le fer révèlent que l'oxygène subsurface n'est stable que lorsque la concentration dépasse un seuil. La barrière de migration pour la diffusion de l'oxygène dans le métal en vrac peut être calculée à l'aide de la méthode de bande élastique à armature en relief.

Les études DFT utilisant des modèles supercellulaires de petits amas d'oxyde (p. ex. (NiO)n, (Al2O3)n pour n=1-10) montrent que l'énergie cohésive de ces amas augmente avec la taille, ce qui correspond éventuellement à l'énergie d'oxyde en vrac. La transition de l'oxygène dispersé à un film d'oxyde continu dépend de l'équilibre entre l'énergie d'adsorption d'oxygène et l'énergie de surface de l'interface métal-oxyde. DFT fournit ces énergies d'interface, permettant de prédire si l'oxyde formera une échelle uniforme ou des précipités isolés.

Rôle des éléments réactifs et des dopants

Les simulations DFT permettent d'élucider le mécanisme : ces éléments se séparent des limites de grain dans l'échelle d'oxyde, où ils renforcent la liaison entre les grains d'oxyde et réduisent le stress de croissance.Les calculs de l'yttrium aux limites de grain Al2O3 montrent une réduction de l'énergie limite de 30 %, inhibant la diffusion des cations d'aluminium vers l'extérieur, ce qui ralentit la croissance de l'oxyde et réduit l'accumulation de vide à l'interface métal-oxyde. Ces observations à l'échelle atomique ont directement influencé la composition des superalliages commerciaux (p. ex., CMSX-4) utilisés dans les pales à turbine fonctionnant au-dessus de 1000 °C.

Formation de couches d'oxyde et propriétés de protection

Structure et stabilité des phases d'oxyde

Les couches d'oxyde formées sur des matériaux à haute température ne sont pas des films simples; elles se composent de plusieurs phases avec des structures distinctes. Par exemple, sur les alliages formant du chrome, la chromie (Cr2O3) avec la structure corundum est thermodynamiquement stable. Les calculs DFT du module en vrac, des constantes élastiques et de l'énergie de surface de Cr2O3 aident à prédire si l'échelle sera dense et adhérente. De même, l'alumine peut exister sous forme de γ-Al2O3, γ-Al2O3, ou α-Al2O3, chacun avec différentes qualités de protection.

Pour les métaux réfractaires et les alliages (p. ex., Mo-Si-B), l'échelle d'oxyde peut être un mélange complexe de verres SiO2, MoO3 et borosilicate. DFT révèle que l'oxyde de molybdène se volatilise au-dessus de 800°C (sous-limité sous le nom de MoO3), ce qui entraîne une oxydation catastrophique connue sous le nom de -pestage. - L'ajout de silicium et de bore favorise la formation d'une couche borosilicate protectrice.

Oxygène et diffusion par cation à l'échelle

La cinétique de croissance de l'oxyde est contrôlée par le transport des ions oxygène vers l'intérieur ou des cations métalliques vers l'extérieur à travers l'échelle. Les calculs DFT des barrières de migration dans les réseaux d'oxydes en vrac sont essentiels. En chromie, les barrières de migration de l'oxygène sont d'environ 1,0 eV, tandis que la migration interstitielle du chrome est d'environ 0,8 eV. En alumine, la diffusion de l'oxygène est significativement plus lente (barrier >3,5 eV) en raison de la structure corundum dense, expliquant pourquoi les échelles d'alumine croissent beaucoup plus lentement que les échelles de chromie.

Les calculs de DFT montrent que le soufre réduit le travail de séparation à l'interface Ni/α-Al2O3 d'environ 8 J/m2 à moins de 3 J/m2, augmentant de façon spectaculaire le risque de spallation. Cette compréhension au niveau atomique a conduit à des recommandations pratiques pour limiter la teneur en soufre dans les superalliages à <10 ppm.

Défis et limites de la DFT dans les études d'oxydation

Malgré sa puissance, DFT est confronté à des défis importants lorsqu'il est appliqué à l'oxydation à haute température. Premièrement, le coût de calcul s'évalue comme O(N3) avec le nombre d'électrons, limitant les tailles du système à quelques centaines d'atomes – bien moins que les millions nécessaires pour modéliser des microstructures réalistes à l'échelle de l'oxyde.

En second lieu, DFT est intrinsèquement une technique de 0 K, bien que les effets de température finis puissent être inclus par des contributions de phonon à l'énergie libre. Pour l'oxydation à 1000°C, l'entropie vibrationnelle et l'anharmonicité deviennent importants, mais ceux-ci nécessitent une dynamique moléculaire ab initio coûteuse ou des approximations quasiharmoniques.

Troisièmement, la précision des prédictions de DFT pour les systèmes fortement corrélés (p. ex. NiO, CoO) est limitée. Les fonctions GGA standard présentent l'écart de bande de NiO comme proche de zéro, alors que la valeur réelle est ~4 eV. Cela affecte la réactivité de surface calculée et la chimie des défauts. DFT+U ou fonctions hybrides améliorent les résultats mais à un coût de calcul plus élevé.

Quatrièmement, l'oxydation implique de multiples phénomènes couplés : transport en phase gazeuse, réactions de surface, diffusion à l'état solide, transformations de phase et contraintes mécaniques. DFT ne traite que de la réaction de surface et de la diffusion à courte portée. Lier les modèles DFT au continuum nécessite une paramétrisation et une validation prudentes.

Orientations futures : DFT avec apprentissage automatique et intégration à plusieurs échelles

Découverte accélérée par le DFT à haut débit

Le flux de travail DFT traditionnel, qui consiste à construire une dalle, à effectuer des calculs pour une poignée de configurations, est trop lent pour permettre de dépister l'immense espace de composition des alliages multicomposants. Le DFT à haut débit automatisé par flux de travail (par exemple AFLOW, Materials Project, Matbench) calcule maintenant les propriétés d'oxydation de milliers de surfaces. Par exemple, une étude de 2023 a analysé les alliages Ni-Co-Cr-Al-Ta pour la ségrégation de surface de l'aluminium et l'énergie d'adsorption de l'oxygène, identifiant Ta comme un puissant amplificateur de la formation d'alumine.

Modèles de substitution pour l'apprentissage automatique

Les potentiels de réseaux neuronaux formés sur les données DFT peuvent simuler la dynamique d'oxydation sur nanosecondes et les distances à l'échelle des microns, impossibles avec DFT seul. Par exemple, des chercheurs de [MIT](https://web.mit.edu/) ont développé un modèle de potentiel profond pour l'oxydation Ni-Cr-Al qui reproduit les énergies DFT à moins de 5 meV/atome.

Correlation de la DFT avec caractérisation expérimentale

L'avenir de la modélisation d'oxydation réside dans une intégration étroite entre la DFT et les expériences avancées. La microscopie électronique à transmission in situ (TEM) combinée à la spectroscopie de perte d'énergie électronique (EELS) peut maintenant mesurer l'état d'oxydation et la liaison locale à l'échelle atomique. Les spectres de DFT-simulés permettent d'identifier directement les phases intermédiaires d'oxyde dans les images TEM. De même, des modèles de diffraction des rayons X de surface peuvent être prédits à partir de modèles de DFT de structure de couche d'oxyde, permettant la validation des voies de réaction calculées.

Conséquences pratiques pour la conception d'alliages à haute température

Les résultats des simulations d'oxydation de DFT ont des conséquences tangibles sur le choix et le traitement des matériaux. Pour les pales de turbines aérospatiales, les lignes directrices de conception de DFT comprennent : a) la fixation d'une concentration critique d'aluminium (>12 à %) à la surface pour former une échelle continue d'α-Al2O3; b) l'ajout de quantités mineures d'éléments réactifs (Y, La) qui se séparent des limites des grains d'oxyde, réduisant ainsi les taux de croissance; c) la réduction de la teneur en soufre inférieure à 5 ppm pour maintenir l'adhérence interfaciale.

Dans le secteur de l'énergie nucléaire, les matériaux de revêtement comme Zircaloy subissent une oxydation dans la vapeur à 1200°C lors d'accidents de perte de refroidissement. Les simulations DFT de diffusion de l'oxygène dans la zircone ont montré que la présence d'étain allié (Sn) augmente l'énergie de formation de l'oxygène, ce qui réduit le taux d'oxydation, ce qui a influencé le développement de formulations avancées de Zircaloy avec une teneur en étain plus faible.

Conclusion

La théorie fonctionnelle de la densité est devenue un outil indispensable pour comprendre et prédire les processus d'oxydation dans les matériaux à haute température.De l'adsorption initiale de l'oxygène à la formation et à la croissance d'échelles d'oxydes de protection, DFT fournit des détails à l'échelle atomique qui sont inaccessibles par seule expérience. Bien que les limites de la taille du système, les effets de la température et la forte corrélation persistent, les développements continus dans les algorithmes, l'informatique à haute performance et l'apprentissage des machines augmentent rapidement l'applicabilité de DFT.