SIMULATION À L'ÉCHELLE MÉSOSCOPIQUE DES GELS D'ALTÉRATION DES VERRES NUCLÉAIRES
Auteur / Autrice : | Téo Boutin |
Direction : | Alain Cartalade |
Type : | Projet de thèse |
Discipline(s) : | Sciences des Matériaux |
Date : | Inscription en doctorat le 25/10/2021 |
Etablissement(s) : | université Paris-Saclay |
Ecole(s) doctorale(s) : | École doctorale Sciences Mécaniques et Energétiques, Matériaux et Géosciences |
Partenaire(s) de recherche : | Laboratoire : CEA/STMF - Service de Thermohydraulique et de Mécanique des Fluides |
Référent : CentraleSupélec |
Mots clés
Résumé
Ce sujet s'inscrit dans le cadre des études relatives au comportement à long terme des verres de conditionnement des déchets radioactifs. L'un des modèles actuellement développé est basé sur une approche chimie-transport [1] (voir les références dans le sujet détaillé) qui permet de rendre compte des phénoménologies d'altération du verre (inter-diffusion, hydrolyse, formation d'un gel, précipitation). Elle permet de déterminer les vitesses d'altération en considérant une description simplifiée du transport réactif des espèces en solution au sein du gel d'altération. La simulation atomistique, complémentaire à cette approche, permet d'affiner la description des mécanismes et valider des hypothèses de ce modèle [2]. L'approche de simulation proposée dans cette thèse considère une échelle intermédiaire, l'échelle mésoscopique, qui permet de rendre compte du transport réactif et de la maturation au sein du gel. Celui-ci est caractérisé par une couche amorphe, nano-poreuse et hydratée, pouvant ralentir de plusieurs ordres de grandeur la vitesse d'altération [3]. Son évolution est essentielle à considérer sur le long terme. Pour cela, on propose de développer un modèle à champ de phase pour lequel les équations aux dérivées partielles seront établies ou adaptées de la littérature [4, Eqs (4a)-(4d)] pour être les plus représentatives des phénoménologies d'altération. Celles-ci seront ensuite programmées dans le code de calcul LBM_saclay [5] d'ores et déjà dédié aux simulations HPC d'un modèle proche (éqs couplées Allen-Cahn/transport/Navier-Stokes) et peut s'exécuter aussi bien sur les architectures muti-CPUs et -GPUs. L'approche, actuellement mise en oeuvre pour la séparation de phase d'un mélange ternaire, sera adaptée pour développer un modèle d'altération des gels en s'appuyant sur la complémentarité du consortium. Des validations expérimentales seront envisagées avec le LCLT.