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

Questo articolo presenta un modello computazionale multiphysics per il Reattore a Sali Fusi (MSFR) che integra il codice ad alta fedeltà Cardinal con il metodo della Matrice di Fissione per risolvere le equazioni di neutronica in modo rapido e accurato, sfruttando database precalcolati tramite simulazioni Monte Carlo.

L'essenziale

Immagina di dover progettare il motore perfetto per un'auto futuristica, ma invece di benzina, questo motore usa sale fuso (come quello che usi per le patatine, ma a temperature altissime) che circola continuamente. Questo è il concetto alla base del Reattore a Sale Fuso (MSFR).

Il problema è che questo "motore" è molto complicato: il sale fa due cose contemporaneamente.

1. È il carburante: contiene l'uranio che produce energia (come il motore che scoppia).
2. È il refrigerante: circola per portare via il calore (come l'acqua che raffredda il motore dell'auto).

Poiché il carburante si muove, il calore e l'energia sono strettamente legati. Se il sale si muove veloce, il calore si disperde; se si muove lento, si surriscalda. Per progettare questo reattore in sicurezza, gli scienziati devono simulare al computer come si comportano queste due cose insieme.

La Sfida: Due Mondi che Devono Parlare

Per fare questa simulazione, gli autori (ricercatori della Penn State) hanno dovuto unire due mondi che normalmente non si parlano:

• Il mondo della Fisica Nucleare (Neutroni): Chi produce l'energia? Dove?
• Il mondo della Fluidodinamica (Sale liquido): Come scorre il sale? Quanto è caldo?

Fino a poco tempo fa, per unire questi due mondi si usava un metodo molto preciso ma lentissimo, come cercare di contare ogni singolo granello di sabbia su una spiaggia (il metodo Monte Carlo/OpenMC). È preciso, ma richiede anni di tempo di calcolo.

La Soluzione Creativa: La "Mappa dei Ricordi" (Il Metodo della Matrice di Fissione)

Invece di contare ogni granello di sabbia ogni volta, gli scienziati hanno usato un trucco intelligente: la Metodo della Matrice di Fissione.

Immagina di voler sapere quanto è calda una stanza.

• Il metodo vecchio: Misurare la temperatura di ogni singolo centimetro cubo d'aria in tempo reale. Lento e faticoso.
• Il metodo nuovo (quello usato in questo paper): Hai già fatto degli esperimenti in passato. Sai che se la stanza è a 80 gradi, la mappa del calore è "A". Se è a 90 gradi, è "B". Se è a 100 gradi, è "C".
  Quando il computer deve simulare la stanza a 85 gradi, non ricomincia da zero. Prende la mappa "A" e la mappa "B" e le mescola (interpolazione) per creare una mappa "A+B" perfetta per 85 gradi.

Questo è esattamente quello che fanno loro: creano una biblioteca di mappe (chiamata Fission Matrix Databases) a diverse temperature. Durante la simulazione, il computer guarda la temperatura attuale del sale e "mescola" le mappe giuste per sapere istantaneamente dove sta producendo energia. È come avere una ricetta culinaria che si adatta automaticamente in base alla temperatura del forno, senza dover riscrivere tutto il libro di cucina ogni volta.

Gli Strumenti del Mestiere

Per far funzionare questo gioco di squadra, hanno usato tre strumenti digitali:

1. NekRS: È il "pilota" che simula come scorre il sale liquido (come l'acqua in una doccia). È velocissimo perché usa le schede grafiche dei computer (GPU), proprio come fanno i videogiochi moderni.
2. MOOSE/Cardinal: È il "direttore d'orchestra". Fa in modo che il pilota (NekRS) e il cuoco (il metodo della Matrice) si parlino. Dice al pilota: "Ehi, qui il sale è caldo, cambia la mappa del calore!" e al cuoco: "Ehi, il sale scorre veloce, aggiorna la produzione di energia!".
3. Il Metodo della Matrice: Il "cuoco" veloce che usa le ricette pre-calcolate.

Cosa Hanno Scoperto?

Hanno fatto una simulazione su un supercomputer potentissimo (Frontier) e i risultati sono stati promettenti:

• Il sale circola e crea delle zone di "stallo" (come un fiume che rallenta in una curva). In queste zone, il sale rimane fermo più a lungo e si scalda di più, raggiungendo le temperature massime.
• La temperatura finale che hanno ottenuto (circa 997 Kelvin, o 724 gradi Celsius) è quasi identica a quella ottenuta con i metodi vecchi e lenti.
• La differenza è minima (pochi gradi), ma il vantaggio è enorme: il loro metodo è molto più veloce e efficiente.

In Sintesi

Questo paper racconta come un gruppo di ricercatori abbia inventato un modo più intelligente per simulare un reattore nucleare del futuro. Invece di fare calcoli enormi e lenti ogni secondo, usano una "biblioteca di ricordi" (le mappe pre-calcolate) per adattarsi istantaneamente alle condizioni del sale fuso.

È come se invece di calcolare la traiettoria di ogni goccia d'acqua in una cascata ogni volta che piove, avessimo imparato a prevedere il flusso dell'acqua basandoci su come si comporta quando piove leggermente, molto, o a tempesta. Il risultato è un modello più veloce, sicuro e pronto per il futuro dell'energia nucleare.

Sommario tecnico

Titolo

Modellazione Multifisica del Reattore a Sali Fusi (MSFR) utilizzando NekRS e il Metodo della Matrice di Fissione

1. Il Problema

Il Reattore a Sali Fusi (MSR) è stato selezionato dal Forum Internazionale Gen-IV come concetto di riferimento per la sua sicurezza, affidabilità e sostenibilità economica. In particolare, il Reattore a Sali Fusi ad Alta Velocità (MSFR) rappresenta un'evoluzione che elimina il moderatore di grafite, operando come reattore veloce autofertilizzante.
La sfida principale nella modellazione dell'MSFR risiede nella sua architettura unica: il combustibile è anche il refrigerante e circola attraverso il sistema primario. Questa caratteristica crea un accoppiamento stretto e complesso tra la neutronica (distribuzione della potenza e della fissione) e l'idrodinamica termica (trasferimento di calore e flusso del fluido). Sviluppare modelli computazionali capaci di gestire questa interazione multifisica in tempo reale è essenziale per l'analisi di sicurezza e le prestazioni del reattore.

2. Metodologia

Il lavoro propone un modello computazionale accoppiato che integra la fluidodinamica ad alta fedeltà con una soluzione neutronica a ordine ridotto.

• Strumento Computazionale: L'implementazione avviene all'interno del framework Cardinal, che opera sopra l'ambiente MOOSE (Multiphysics Object-Oriented Simulation Environment). Cardinal è progettato per accoppiare codici di diverse fisica.
• Idrodinamica (CFD): Viene utilizzato NekRS, un codice CFD open-source accelerato da GPU basato sul metodo degli elementi spettrali (SEM). NekRS risolve le equazioni di Navier-Stokes non stazionarie e l'equazione dell'energia per il flusso incompressibile, utilizzando il modello di turbolenza RANS $k-\tau$.
• Neutronica (Metodo della Matrice di Fissione - FM): Invece di utilizzare il codice Monte Carlo standard (OpenMC) per la neutronica, gli autori implementano il Metodo della Matrice di Fissione (FM) come modello a ordine ridotto.
  • Una matrice di fissione $A$ ($n \times n$) descrive la generazione di neutroni tra celle spaziali.
  • Per gestire le variazioni di temperatura del combustibile, vengono generati Database di Matrici di Fissione (FMDB) utilizzando il codice Monte Carlo Serpent per diverse temperature uniformi (da 800 K a 1200 K).
  • Durante la simulazione, la matrice $A$ viene ricostruita in tempo reale tramite interpolazione lineare tra i database più vicini in base alla distribuzione di temperatura locale calcolata da NekRS.
• Schema di Accoppiamento: Viene utilizzato uno schema iterativo di Picard ("in time").
  1. Il modello FM calcola la distribuzione della fonte di fissione e la potenza volumetrica ($q'''_{fis}$) basandosi sulla temperatura attuale.
  2. Questa sorgente termica viene passata a NekRS.
  3. NekRS esegue passi temporali per aggiornare velocità e temperatura del fluido.
  4. Le nuove temperature vengono restituite al modello FM per aggiornare la matrice di fissione.
  5. Il ciclo si ripete fino al raggiungimento dello stato stazionario.

3. Contributi Chiave

• Sostituzione di OpenMC con il Metodo FM: L'innovazione principale è l'uso del metodo della Matrice di Fissione all'interno del framework Cardinal per risolvere le equazioni neutroniche. Questo approccio riduce drasticamente il costo computazionale rispetto ai metodi Monte Carlo tradizionali, rendendo possibile l'accoppiamento diretto con simulazioni CFD ad alta risoluzione.
• Integrazione GPU-CPU: La simulazione è stata eseguita sul supercomputer Frontier, sfruttando le GPU per il calcolo CFD (NekRS) e le CPU per il modello FM, dimostrando l'efficienza di un'architettura ibrida.
• Validazione Preliminare: Il modello è stato validato confrontando i risultati con uno studio precedente che utilizzava l'accoppiamento NekRS-OpenMC, confermando la fattibilità dell'approccio a ordine ridotto.

4. Risultati

Le simulazioni preliminari sono state condotte su una geometria semplificata del loop primario dell'MSFR (1,35 milioni di elementi esaedrici per NekRS e una griglia cartesiana 13x13x13 per il FM).

• Parametri di Funzionamento: Il fattore di moltiplicazione ($k_{eff}$) calcolato è 1.04632.
• Distribuzione Termica e di Flusso:
  • La sorgente di calore mostra una forma sinusoidale.
  • Sono state identificate zone di ristagno (stagnation regions) nella periferia del cilindro centrale, dove la velocità è bassa e la temperatura raggiunge i massimi valori a causa del continuo riscaldamento del fluido "intrappolato".
  • La temperatura di uscita calcolata è 997.4 K e la temperatura media è 973.0 K.
• Confronto con il Modello di Riferimento (NekRS-OpenMC):
  • La differenza nella temperatura di uscita è minima (0.1 K).
  • La differenza nella temperatura media è di 10.1 K.
  • Le distribuzioni di velocità e sorgente di calore mostrano un buon accordo qualitativo e quantitativo.
  • Le piccole discrepanze sono attribuite al fatto che l'energia totale è mantenuta costante in entrambi i modelli, ma le distribuzioni spaziali locali variano leggermente a causa della discretizzazione della griglia FM.

5. Significato e Conclusioni

Questo studio dimostra che l'uso del Metodo della Matrice di Fissione come sostituto a ordine ridotto per i codici Monte Carlo in un contesto di accoppiamento multifisico è una strategia valida ed efficiente per la modellazione dei reattori a sali fusi.

• Efficienza: Permette di ottenere risultati fisicamente coerenti con costi computazionali ridotti, facilitando l'analisi di scenari dinamici complessi.
• Accuratezza: I risultati preliminari confermano che il modello cattura correttamente i fenomeni chiave, come le zone di ristagno termico e l'accoppiamento feedback temperatura-potenza.
• Prospettive Future: I lavori futuri includeranno un'analisi di sensibilità sulla griglia del modello FM (testando mesh più fini o più grezze) e l'analisi dell'influenza della forma della mesh (ad esempio, l'uso di mesh cilindriche) sui risultati finali.

In sintesi, il paper fornisce una base solida per lo sviluppo di strumenti di simulazione avanzati per la prossima generazione di reattori nucleari, bilanciando fedeltà fisica e costi di calcolo.