Multiphysics Modelling of the Molten Salt Fast Reactor using NekRS and the Fission Matrix Method

Cet article propose un modèle multiphysique à couplage serré pour le réacteur à sels fondus à neutrons rapides (MSFR) en intégrant la méthode de la matrice de fission, une approche neutonique rapide et précise, au sein du code Cardinal basé sur le framework MOOSE pour coupler la dynamique des fluides via NekRS et la thermique.

L'essentiel

🌟 Le Projet : Une Centrale Nucléaire qui "Respire"

Imaginez une centrale nucléaire traditionnelle comme un gros moteur à piston : le combustible est solide, il reste immobile dans des barres, et un liquide froid passe autour pour l'évacuer. C'est un peu rigide.

Les chercheurs de l'Université d'État de Pennsylvanie travaillent sur une idée plus moderne : le Réacteur à Sels Fondus (MSFR). Ici, le combustible n'est pas solide, c'est un sel liquide bouillant qui circule partout dans le réacteur. C'est comme si le carburant de la voiture était le liquide qui coule dans les tuyaux !

Le problème : Comme le combustible bouge, il emporte la chaleur avec lui. Cela crée une danse très complexe entre la physique des neutrons (qui font exploser les atomes) et l'hydraulique (comment le liquide coule et chauffe). C'est ce qu'on appelle un problème "multiphysique".

🛠️ La Solution : Une Équipe de Super-Héros Numériques

Pour simuler ce réacteur sans le construire réellement, les chercheurs ont créé un modèle informatique très puissant. Ils ont assemblé trois outils principaux, comme une équipe de super-héros :

1. NekRS (Le Coureur de Vitesse) : C'est un logiciel ultra-rapide (qui utilise des puces graphiques de jeux vidéo) qui simule comment le sel liquide coule, tourbillonne et chauffe. Il est très précis pour voir les petits détails de l'écoulement.
2. MOOSE (Le Chef d'Orchestre) : C'est le cadre qui permet aux différents logiciels de se parler et de travailler ensemble.
3. La Méthode de la "Matrice de Fission" (Le Calculateur Rapide) : C'est ici que réside l'innovation. Habituellement, pour savoir où les atomes explosent, on utilise des simulations de Monte Carlo (comme lancer des millions de dés virtuels), ce qui est très précis mais très lent.
   • L'analogie : Imaginez que vous voulez prédire la météo. La méthode classique consiste à lancer des millions de balles dans le ciel pour voir où elles atterrissent. La méthode utilisée ici, c'est comme avoir une carte météo pré-calculée pour différentes températures. Au lieu de tout recalculer, on regarde la carte la plus proche de la température actuelle et on fait une petite estimation (une interpolation). C'est beaucoup plus rapide !

🔄 Comment ça marche ? (La Danse en Boucle)

Le modèle fonctionne comme une conversation en boucle entre deux amis :

1. Le Calculateur (Matrice) dit : "Avec cette température, voici où l'énergie est produite."
2. Le Coureur (NekRS) reçoit cette information et dit : "D'accord, je vais faire couler le sel liquide avec cette chaleur. Regarde, le liquide chauffe ici et ralentit là !"
3. Le Coureur renvoie les nouvelles températures au Calculateur.
4. Le Calculateur ajuste son estimation et on recommence.

Ils répètent cette boucle des milliers de fois jusqu'à ce que tout se stabilise, comme un thermostat qui cherche la température parfaite.

📊 Les Résultats : Une Danse Parfaite

Les chercheurs ont fait tourner leur simulation sur un supercalculateur géant (Frontier). Voici ce qu'ils ont découvert :

• Les Zones de "Trappe" : Ils ont remarqué que dans certains coins du réacteur, le liquide coule très lentement. C'est comme un embouteillage sur une autoroute. Comme le liquide stagne, il continue d'être chauffé sans être refroidi. C'est là que la température atteint son maximum.
• La Précision : Ils ont comparé leur nouvelle méthode (rapide) avec l'ancienne méthode (lente mais très précise).
  • Résultat : Les deux méthodes donnent des résultats presque identiques pour la température de sortie (à 0,1 degré près !).
  • Conclusion : La nouvelle méthode est aussi bonne mais beaucoup plus rapide à calculer.

🚀 Pourquoi c'est important ?

Ce papier montre qu'on peut utiliser des "modèles réduits" (comme la carte météo pré-calculée) pour simuler des réacteurs nucléaires complexes sans perdre en précision.

C'est comme si, au lieu de reconstruire un avion en bois pour tester son vol, on utilisait un simulateur de vol ultra-réaliste qui s'adapte en temps réel. Cela permet aux ingénieurs de tester des idées plus vite, de concevoir des réacteurs plus sûrs et plus efficaces pour l'avenir de l'énergie propre.

En résumé : Ils ont créé un simulateur numérique où un liquide chaud et des réactions atomiques dansent ensemble, en utilisant un raccourci mathématique intelligent pour aller plus vite, tout en restant aussi précis qu'une horloge suisse.

Résumé technique

Titre

Modélisation multiphysique du Réacteur à Sels Fondus à Spectre Rapide (MSFR) utilisant NekRS et la méthode de la matrice de fission.

1. Problématique et Contexte

Le Réacteur à Sels Fondus (MSR), et plus spécifiquement le concept de Réacteur à Sels Fondus à Spectre Rapide (MSFR), a été sélectionné par le Forum International Génération IV pour sa sécurité, sa réduction des déchets et son efficacité économique. La particularité du MSFR réside dans le fait que le caloporteur est également le combustible, ce qui crée un couplage étroit entre la neutronique et l'hydraulique thermique (la circulation du combustible modifie la distribution de puissance et vice-versa).

Le défi principal est de développer des modèles de calcul capables de gérer ce couplage multiphysique complexe de manière efficace. L'approche standard utilise souvent des codes de transport Monte Carlo (comme OpenMC) couplés à la dynamique des fluides (CFD), mais ces méthodes peuvent être coûteuses en temps de calcul. L'objectif de cet article est de proposer une alternative utilisant un modèle d'ordre réduit pour résoudre les équations de neutronique, tout en conservant une haute fidélité pour l'hydraulique.

2. Méthodologie

L'étude repose sur le développement d'un modèle couplé neutronique-hydraulique thermique intégré dans le cadre logiciel Cardinal, qui s'appuie sur l'environnement MOOSE (Multiphysics Object-Oriented Simulation Environment).

• Hydraulique Thermique (CFD) : Le code NekRS est utilisé. Il s'agit d'un code CFD accéléré par GPU basé sur la méthode des éléments spectraux (SEM). Il résout les équations de Navier-Stokes incompressibles et de l'énergie en régime transitoire. Pour cette étude, le modèle de turbulence RANS $k-\tau$ est employé.
• Neutronique (Méthode de la Matrice de Fission - FM) : Au lieu d'utiliser le code Monte Carlo OpenMC (habituellement couplé dans Cardinal), les auteurs implémentent la Méthode de la Matrice de Fission (FM).
  • Une matrice de fission $A$ est définie sur une discrétisation spatiale de la géométrie du réacteur.
  • Les éléments de la matrice représentent le nombre de neutrons générés dans une cellule $i$ à partir d'un neutron de fission né dans la cellule $j$.
  • La résolution de cette matrice fournit le facteur de multiplication effectif ($k_{eff}$) et la distribution de la source de fission.
• Génération des Bases de Données (FMDB) : Des bases de données de matrices de fission sont pré-calculées à l'aide du code Monte Carlo Serpent pour différentes distributions de température uniforme du combustible (de 800 K à 1200 K).
• Interpolation Linéaire : Pour représenter l'état réel du réacteur (avec des températures non uniformes), la matrice de fission est reconstruite en temps réel par interpolation linéaire entre les deux bases de données les plus proches, en fonction de la température locale du combustible fournie par le code CFD.
• Couplage (Itération de Picard) : Le schéma de couplage suit une boucle itérative :
  1. Le modèle FM calcule la source de fission et la puissance volumique ($q'''$) basée sur une température initiale.
  2. Cette source de chaleur est transférée à NekRS.
  3. NekRS calcule l'écoulement, la vitesse et la température du fluide sur plusieurs pas de temps.
  4. La nouvelle distribution de température est renvoyée au modèle FM pour mettre à jour la matrice de fission.
  5. Le processus se répète jusqu'à convergence des températures (entrée, sortie et moyenne).

3. Contributions Clés

• Implémentation d'un modèle d'ordre réduit : Développement d'une interface C++ dans MOOSE pour intégrer la méthode de la matrice de fission, permettant de remplacer le calcul Monte Carlo lourd par une résolution matricielle rapide.
• Couplage CFD-Neutronique avancé : Réalisation d'un couplage fort entre NekRS (GPU) et le modèle FM (CPU) via Cardinal, démontrant la faisabilité de l'utilisation de la méthode FM pour des simulations multiphysiques transitoires de MSFR.
• Validation préliminaire : Comparaison directe des résultats avec un modèle de référence existant (NekRS couplé à OpenMC) pour valider la précision de l'approche FM.

4. Résultats

Les simulations ont été exécutées sur le supercalculateur Frontier (utilisant des GPU pour NekRS et des CPU pour le modèle FM).

• Paramètres de convergence : Après environ 60 unités convectives, une moyenne temporelle a été effectuée sur 80 unités supplémentaires. Le facteur de multiplication obtenu est de 1,04632.
• Distribution de température et de vitesse :
  • La distribution de vitesse montre des zones de stagnation périphériques autour du cylindre central.
  • Ces zones de stagnation correspondent aux pics de température, car le fluide y est « piégé » et continue d'être chauffé.
  • La température de sortie calculée est de 997,4 K et la température moyenne de 973,0 K.
• Comparaison avec le modèle de référence (NekRS-OpenMC) :
  • La différence de température de sortie est négligeable (0,1 K).
  • La différence de température moyenne est de 10,1 K.
  • Les auteurs attribuent cette légère divergence au fait que l'énergie totale est identique dans les deux modèles, mais la répartition spatiale de la source de chaleur diffère légèrement.
• Accord global : Les distributions de source de chaleur, de température et de vitesse montrent un bon accord qualitatif et quantitatif avec le modèle de référence, validant l'approche FM pour ce type de géométrie.

5. Signification et Perspectives

Cette étude démontre que la méthode de la matrice de fission est une alternative viable et potentiellement plus efficace en termes de temps de calcul pour les simulations multiphysiques des réacteurs à sels fondus, par rapport aux méthodes Monte Carlo traditionnelles couplées.

• Avantages : Réduction de la charge computationnelle tout en conservant une précision suffisante pour l'analyse de la sécurité et de la performance thermique.
• Travaux futurs : Les auteurs prévoient d'effectuer une analyse de sensibilité sur la finesse du maillage du modèle FM (cartésien vs cylindrique) afin d'optimiser la discrétisation spatiale et d'affiner la précision des résultats.

En conclusion, ce travail ouvre la voie à l'utilisation de modèles réduits pour le couplage multiphysique haute fidélité des réacteurs de nouvelle génération, facilitant ainsi les études paramétriques et les analyses de sûreté.