Contexte : EDF R&D développe des outils de calcul scientifique exploités par les centres d’ingénierie d’EDF dans le cadre de la démonstration de sureté des réacteurs à eau pressurisée (REP) du parc nucléaire français. Parmi ces outils, THYC est le logiciel de référence d’EDF pour la simulation 3D des écoulements monophasiques et diphasiques (eau-vapeur) dans le coeur des réacteurs en conditions normales d’exploitation ou en conditions d’accident.
THYC résout un modèle décrivant l’évolution du mélange diphasique ainsi que l’équation de conservation de la masse de vapeur. Pour modéliser le milieu très encombré des coeurs de réacteurs principalement constitués de faisceaux de crayons combustibles regroupés en assemblages, THYC résout les équations de conservation homogénéisées sur un maillage du domaine contenant le fluide et les obstacles solides, assimilé à un milieu poreux. La discrétisation repose sur un schéma numérique volumes finis sur grilles décalées en espace et semi-implicite en temps à l’aide d’une méthode à correction de pression.
Des travaux récents [1], menés à EDF R&D, ont permis le développement d’une nouvelle approche pour modéliser les écoulements en milieu encombré en discrétisant les équations de conservation non homogénéisées et en prenant en compte la nature du milieu (fluide ou solide) lors de l’étape d’intégration sur les volumes finis (mailles), un choix qui permet de prendre en compte localement plus naturellement la géométrie des obstacles (notamment les sections de passage du fluide entre les crayons).

Objectifs du stage :

L’approche dite, formulation intégrale, a été appliquée dans le cadre du schéma numérique à mailles décalées de THYC et les principaux développements et des évaluations dans THYC ont été réalisés en 2023 et 2024. On souhaite dans ce stage travailler sur les aspects suivants :
- Vérification du code de calcul sur des cas test instationnaires avec solutions analytiques
- Analyse et amélioration des conditions aux limites
- Sensibilité de la formulation intégrale au raffinement du maillage du domaine de calcul
- Amélioration du schéma numérique de l’équation de conservation de l’énergie

Unité d’accueil :

EDF R&D, Département Mécanique des Fluides, Energies et Environnement
6 quai Watier - 78400 CHATOU |

RER Ligne A - Station Rueil Malmaison

Contacts : erwan.lecoupanec@edf.fr - clement.colas@edf.fr

Référence :

[1] - Formulation intégrale implicite pour la modélisation d'écoulements fluides en milieu encombré, thèse de Clément Colas, 2019, http://www.theses.fr/2019AIXM0555

Stage de fin d'études ingénieur ou équivalent Master II de 6 mois à partir de février/mars/avril 2025

Compétences attendues :

Mécanique des fluides, méthodes numériques, Fortran90/Python