Ce sujet s’inscrit dans le cadre des études réalisées sur le comportement à long terme des verres nucléaires immobilisant des déchets radioactifs et potentiellement destinés à être placés en stockage géologique. L’enjeu réside en la compréhension des mécanismes d’altération et de formation d'un gel (couche passivante pouvant ralentir la vitesse d’altération du verre) par l’eau et à la prédiction des cinétiques de relâchement des radionucléides sur le long terme.

L’approche de simulation proposée vise à prédire à l’échelle mésoscopique le processus de maturation du gel formé lors de l’altération du verre par l’eau à l'aide d'un 'modèle à champs de phases' ternaire composé du silicium, du bore et de l’eau (lixiviant) auquel il conviendra d'ajouter l'alluminium.

Le modèle mathématique quaternaire sous-jacent est composé d’un ensemble d’Equations aux Dérivées Partielles non-linéaires couplées. Elles sont basées sur les équations de Allen-Cahn et du transport. La résolution numérique des équations associées est réalisée par méthode de Boltzmann sur réseaux (Lattice Boltzmann Method – LBM) programmée en C++ dans le code de calcul massivementparallèle LBM_saclay qui s’exécute sur plusieurs architectures HPC, aussi bien muti-CPUs que multi-GPUs.

Le sujet proposé nécessite de bonnes bases en mathématiques appliquées et en programmation afin de développer les algorithmes nécessaires à la bonne résolution du nouveau système d'équations fortement couplées.