Cette thèse s'inscrit dans le cadre de la recherche nécessaire pour l'utilisation durable de l'énergie nucléaire dans un mix énergétique décarboné respectueux du climat. Les réacteurs de 4° génération refroidis au sodium sont alors des candidats de grand intérêt pour l'économie de la ressource en uranium et la minimisation du volume des déchets ultimes.

Dans le cadre de la sûreté de tels réacteurs, il est important de pouvoir décrire avec précision les conséquences d'une éventuelle dégradation du cœur. Une collaboration avec son homologue japonais JAEA permet au CEA de développer le code SIMMER-V dédié à la simulation de la dégradation du cœur. Le code calcule la thermohydraulique du sodium, la dégradation des structures et la neutronique du cœur pendant la phase accidentelle. L'objectif est de pouvoir représenter non seulement le cœur mais aussi son environnement direct (circuit primaire) avec précision. La prise en compte de cette topologie requiert de partitionner le domaine et d'utiliser une méthode de couplage aux frontières. La limite de cette approche réside généralement dans la qualité et la robustesse de la méthode de couplage, en particulier lors de transitoires rapides au cours desquels des ondes de pression et de densité sont amenées à traverser les frontières.

Une méthode de couplage a été initiée au Laboratoire de Modélisation des Accidents Graves de l’institut IRESNE de Cadarache (Annals of Nuclear Energy 2022,Implementation of multi-domains in SIMMER-V thermohydraulic code https://doi.org/10.1016/j.anucene.2022.109338), qui consiste à fusionner les différentes décompositions de chacun des domaines, dans le but de constituer une décomposition unique du calcul global. Cette méthode a été développée dans un cadre simplifié où les maillages (cartésiens) se raccordent de manière conforme au niveau des frontières. L'opportunité qui s'ouvre est d'étendre cette méthode au cas des raccords quelconques en utilisant la librairie MEDCoupling. Cette première étape dont la faisabilité a été acquise permettra d'assembler des composants pour constituer un système de type 'loop' (réacteur à boucles). La deuxième étape consistera à étendre la méthode de sorte qu'un domaine de calcul puisse être totalement emboîté dans un autre. Cet emboîtement permettra alors de constituer un domaine par juxtaposition ou par emboîtement avec des maillages et des décompositions de domaine non conformes. Après avoir vérifié les qualités numériques de la méthode, la dernière étape applicative consistera à construire une simulation de la dégradation d'un cœur plongé dans sa cuve primaire (configuration 'pool') permettant de valider la méthode suivie.

Ce travail permettra à l’étudiant de développer des connaissances en techniques numériques et modélisation pour les systèmes physiques complexes avec écoulements. Il mettra en œuvre des techniques allant de la conception à la validation de méthodes, dans une équipe pluridisciplinaire et dynamique au CEA Cadarache.