Thèse soutenue

Simulation atomistique Monte Carlo Cinétique des processus de croissance de couches passives sur alliage métalliques : cas des alliages Fer-Chrome

FR  |  
EN
Auteur / Autrice : Yves Alain Beh Ongueng
Direction : Philippe Marcus
Type : Thèse de doctorat
Discipline(s) : Simulation informatique et matériaux
Date : Soutenance en 2008
Etablissement(s) : Paris 6

Résumé

FR  |  
EN

La croissance de couches minces d’oxydes sur les alliages métalliques et les métaux purs est un phénomène ayant fait l’objet d’un grand nombre d’études expérimentales. Les structures des couches d’oxydes sont bien connues, ainsi que les modèles mathématiques servant à modéliser les aspects macroscopiques de la croissance. Cependant, les détails des mécanismes de germination de la couche d’oxyde dans les premiers stades de la corrosion ainsi que sa croissance ultérieure demeurent peu ou mal connus. La simulation atomistique apparaît comme une alternative pour évaluer les différents mécanismes proposés et appréhender l’influence des différents paramètres physico-chimiques. Le développement d’un tel outil de simulation a démarré au LPCS avec la thèse de M. Legrand. En se basant sur l’exemple de l’alliage FeCr, un modèle informatique tridimensionnel dit «modèle de Legrand», permettant de simuler la dissolution sélective et la passivation des alliages binaires a été réalisé. L’évolution dynamique est basée sur une technique de type Monte Carlo classique. Le logiciel permet de simuler l’évolution d’un alliage quelconque, d’une composition et d’une structure cristallographique donnée. Il prend en compte la diffusion des atomes sur la surface, leur dissolution et le blocage de la dissolution par formation d’une couche de passivation. Cet outil était adapté pour la simulation des premiers stades de la corrosion. L’objectif de ce travail est d’améliorer ce modèle existant, afin de simuler l’évolution de la couche passive sur une échelle de temps plus longue. A l’issue de ce travail, de nombreux apports ont été effectués. Ainsi, l’introduction d’un champ de force MEAM (Modified Embedded Atom Method) pour le calcul des barrières de diffusion et de dissolution, a permis de remplacer les probabilités de diffusion empiriques par des probabilités calculées, et de mettre en évidence la diffusion préférentielle des Cr vers leurs semblables. L’introduction d’une dynamique de simulation Monte Carlo Cinétique (KMC) a permis une prise en compte réaliste de l'évolution cinétique du modèle avec des temps de simulation reliés au temps réel. Enfin Un second réseau cristallographique RVO (réseau virtuel d’oxyde) tridimensionnel, correspondant à celui de la couche passive (Cr2O3) a été implémenté, ainsi qu’une interface graphique pour un meilleur suivi de la simulation. Les résultats obtenus lors des simulations sont en accord avec les observations expérimentales: passivation totale à partir du Fe-16Cr, enrichissement en Cr de la couche passive, allures des courbes cinétiques, influence du champ électrique, mise en évidence l’apparition de cavités sous la couche passive.