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

Diese Arbeit stellt ein Multiphysik-Modell für den Flüssigsalz-Brutreaktor vor, das den hochauflösenden CFD-Code NekRS mit der schnellen und präzisen Spaltmatrix-Methode zur Neutronenphysik innerhalb des MOOSE-Frameworks kombiniert.

Das Wesentliche

Das große Ziel: Ein flüssiger Kernreaktor

Stellt euch einen Atomreaktor nicht als riesigen Steinblock mit festen Brennstäben vor, sondern als einen riesigen, kochenden Topf Suppe. Das ist die Idee hinter dem Molten Salt Fast Reactor (MSFR).

In diesem "Topf" ist das Brennstoffgemisch (eine spezielle Salzflüssigkeit) gleichzeitig auch das Kühlmittel. Es zirkuliert ständig durch das System. Das ist genial, aber auch kompliziert: Wenn sich die Flüssigkeit bewegt, ändert sich die Temperatur. Wenn sich die Temperatur ändert, verändert sich, wie die Atomkerne spalten (Neutronen). Und wenn sich die Spaltung ändert, wird es wieder heißer.

Das ist wie ein Tanz zwischen Hitze und Bewegung. Man kann das eine nicht ohne das andere verstehen. Genau das wollen die Forscher von der Pennsylvania State University simulieren.

Die Werkzeuge: Ein digitales Labor

Um diesen Tanz am Computer nachzustellen, brauchen sie zwei spezielle Werkzeuge, die sie zu einem Team zusammengeschweißt haben:

1. NekRS (Der Strömungs-Experte):
   Stell dir NekRS als einen super-schnellen Wettervorhersage-Computer vor, nur für den Reaktor. Er berechnet, wie die heiße Salzflüssigkeit fließt, wo sie sich staut und wie heiß sie wird. Er nutzt dabei die Kraft von Grafikkarten (GPUs), also den gleichen Chips, die in Gaming-Computern stecken, um die Berechnungen blitzschnell durchzuführen.

2. Die Fissions-Matrix (Der Neutronen-Zähler):
   Normalerweise würde man für die Atomspaltung einen Monte-Carlo-Simulator nutzen (wie OpenMC), der jedes einzelne Neutron einzeln verfolgt. Das ist extrem genau, aber auch wie das Zählen jedes einzelnen Sandkorns am Strand – sehr langsam.
   Die Forscher nutzen hier eine Abkürzung: Die Fissions-Matrix.

   • Die Analogie: Stell dir vor, du hast ein Schachbrett (den Reaktor). Anstatt zu berechnen, wie sich jedes einzelne Neutron bewegt, hast du eine Tabelle (die Matrix), die dir sagt: "Wenn in Feld A ein Neutron geboren wird, wie viele neue Neutronen landen wahrscheinlich in Feld B?"
   • Das ist wie ein Kochrezept, das dir sofort sagt, wie viel Salz du brauchst, ohne dass du jedes Salzgranulat einzeln wiegen musst. Es ist schneller und für ihre Zwecke genau genug.

Wie das Team zusammenarbeitet (Der "Tanz")

Die Software Cardinal ist der Choreograf, der diese beiden Systeme zusammenbringt. So läuft das Spiel ab:

1. Der Start: Die "Neutronen-Matrix" schaut sich die aktuelle Temperatur an und sagt: "Okay, bei dieser Hitze passiert hier und dort Spaltung." Sie berechnet, wo die Hitze entsteht.
2. Die Weitergabe: Diese Hitze-Information wird an den "Strömungs-Experten" (NekRS) gegeben.
3. Die Reaktion: NekRS berechnet: "Oh, hier wird es heiß! Die Flüssigkeit strömt jetzt anders, sie kühlt sich hier ab und staut sich dort."
4. Der Rückkopplung: NekRS schickt die neuen Temperaturdaten zurück an die Matrix.
5. Wiederholung: Dieser Kreislauf läuft tausende Male pro Sekunde, bis sich ein stabiles Gleichgewicht einstellt.

Was haben sie herausgefunden?

Die Forscher haben ihren neuen, schnelleren Weg (Matrix statt Einzel-Neutronen-Simulation) getestet und mit dem alten, langsamen Weg verglichen.

• Das Ergebnis: Es funktioniert fast perfekt! Die Temperaturen und Strömungsmuster waren fast identisch mit dem alten Modell.
• Ein kleiner Unterschied: Die neue Methode war etwas schneller, aber die Temperatur am Ausgang war nur winzig anders (0,1 Kelvin Unterschied). Das ist wie wenn zwei Uhren eine Sekunde unterschiedlich gehen – für den Reaktorbetrieb völlig egal.
• Interessante Entdeckung: Sie sahen "Stagnationszonen". Das sind Stellen im Reaktor, wo die Flüssigkeit fast stehen bleibt (wie ein Stau im Verkehr). Dort wird die Flüssigkeit immer weiter erhitzt, weil sie nicht wegfließt. Das ist wichtig zu wissen, damit der Reaktor nicht an einer Stelle zu heiß wird.

Fazit

Die Forscher haben einen neuen, effizienteren Weg gefunden, um diese komplexen Atomreaktoren am Computer zu simulieren. Anstatt jedes einzelne Neutron zu jagen, nutzen sie eine intelligente Tabelle (die Matrix), die den Tanz zwischen Hitze und Strömung fast genauso gut beschreibt wie der alte, langsame Weg.

Das ist ein großer Schritt, um sicherere und effizientere Atomkraftwerke der Zukunft zu bauen, die weniger Abfall produzieren und sicherer sind als die heutigen.

Technische Zusammenfassung

Technische Zusammenfassung: Multiphysik-Modellierung des Molten Salt Fast Reactors (MSFR) mit NekRS und der Spaltmatrix-Methode

1. Problemstellung und Motivation
Der Molten Salt Reactor (MSR) wurde vom Generation-IV-International-Forum als Referenzkonzept ausgewählt. Eine spezifische Weiterentwicklung, der Molten Salt Fast Reactor (MSFR), verzichtet auf einen Graphitmoderator und fungiert als schneller Brutreaktor. Ein entscheidendes Merkmal des MSFR ist, dass das Kühlmittel gleichzeitig als Brennstoff dient. Dies führt zu einer starken Kopplung zwischen der Neutronenphysik (Neutronik) und der Thermohydraulik, da der Brennstoff durch das primäre System zirkuliert. Die Entwicklung von Rechenmodellen zur Analyse des MSFR erfordert daher einen Multiphysik-Ansatz. Das Ziel dieses Papers ist die Entwicklung eines gekoppelten neotronisch-thermohydraulischen Modells, das eine reduzierte Ordnungsmethode zur Lösung der Neutronikgleichungen verwendet, um die Rechenkosten im Vergleich zu vollständigen Monte-Carlo-Simulationen zu senken.

2. Methodik
Das vorgestellte Modell nutzt den Cardinal-Code, eine Open-Source-Anwendung innerhalb des MOOSE-Frameworks (Multiphysics Object-Oriented Simulation Environment), die die Kopplung verschiedener physikalischer Codes ermöglicht.

• Thermohydraulik (NekRS):
  Für die Simulation der Strömung und des Wärmetransports wird NekRS verwendet, ein GPU-beschleunigter CFD-Code, der auf der Spektral-Elemente-Methode (SEM) basiert. Er löst die instationären, inkompressiblen Navier-Stokes- und Energiegleichungen. In dieser Studie wird das RANS-Modell ($k-\tau$ Turbulenzmodell) angewendet. Das geometrische Modell des MSFR-Hauptkreislaufs ist vereinfacht (kontinuierlicher Einlass/Auslass, keine Pumpen/Wärmetauscher) und besteht aus einem 3D-Gitter mit ca. 1,35 Millionen hexaedrischen Elementen (Polynomordnung 3).

• Neutronik (Spaltmatrix-Methode - FM):
  Anstelle des üblichen Monte-Carlo-Codes OpenMC wird die Spaltmatrix-Methode (Fission Matrix, FM) verwendet.

  • Prinzip: Die Geometrie wird in $n$ räumliche Zellen diskretisiert. Eine Spaltmatrix $A$ beschreibt die Anzahl der Neutronen in Zelle $i$, die aus einer Spaltung in Zelle $j$ stammen. Der Eigenvektor entspricht der Spaltquellenverteilung, der Eigenwert dem effektiven Multiplikationsfaktor ($k_{eff}$).
  • Datenbanken (FMDBs): Um den temperaturabhängigen Zustand des Reaktors abzubilden, wurden Spaltmatrizen mit dem Monte-Carlo-Code Serpent für verschiedene uniforme Brennstofftemperaturen (800 K bis 1200 K) generiert.
  • Interpolation: Da die tatsächliche Temperaturverteilung im Reaktor nicht uniform ist, wird die Spaltmatrix für den aktuellen Zustand durch lineare Interpolation zwischen den nächstgelegenen Einträgen der FMDBs berechnet.
  • Implementierung: Der FM-Löser wurde als C++-Modul in MOOSE integriert. Er arbeitet auf einem kartesischen Gitter (13x13x13 Elemente) und berechnet basierend auf der Temperaturverteilung von NekRS die Spaltquellenverteilung und den $k_{eff}$.

• Kopplung (Cardinal):
  Die Kopplung erfolgt über eine Picard-Iteration "in der Zeit".

  1. Der FM-Modell berechnet die Spaltquellenverteilung basierend auf einer initialen Temperatur.
  2. Die Quelle wird skaliert und als volumetrische Wärmequelle ($q'''_{fis}$) an NekRS übergeben.
  3. NekRS führt $N$ Zeitschritte durch, um neue Strömungs- und Temperaturfelder zu berechnen.
  4. Die neuen Temperaturen werden an das FM-Modell zurückgespeist.
  5. Der Prozess wiederholt sich, bis Konvergenz (stationäre Werte für Einlass-, Durchschnitts- und Auslasstemperatur) erreicht ist.

3. Wichtige Beiträge

• Ersetzung von OpenMC durch FM: Der erste Nachweis einer erfolgreichen Kopplung von NekRS mit einem reduzierten Neutronik-Modell (FM) innerhalb des Cardinal-Frameworks für einen MSFR.
• Entwicklung einer FM-Datenbank: Erstellung und Validierung von Spaltmatrix-Datenbanken für verschiedene Temperaturen mittels Serpent, die eine schnelle lineare Interpolation im gekoppelten Lauf ermöglichen.
• Implementierung in MOOSE: Entwicklung von C++-Funktionen innerhalb von MOOSE zum Laden der FMDBs, zur Interpolation und zur Lösung des Eigenwertproblems im Multiphysik-Kontext.
• Validierung: Vergleich der Ergebnisse mit einem bestehenden NekRS-OpenMC-Modell aus der Literatur.

4. Ergebnisse
Die Simulationen wurden auf dem Supercomputer Frontier (GPU für NekRS, CPU für FM) durchgeführt.

• Kinetische Parameter: Der erhaltene Multiplikationsfaktor ($k_{eff}$) beträgt 1,04632.
• Thermische Ergebnisse:
  • Die Auslasstemperatur beträgt 997,4 K.
  • Die durchschnittliche Temperatur beträgt 973,0 K.
  • Im Vergleich zum Referenzmodell (NekRS-OpenMC) weicht die Auslasstemperatur nur um 0,1 K ab, die Durchschnittstemperatur um 10,1 K. Diese geringen Abweichungen werden auf die gleiche Gesamtenergie in beiden Modellen zurückgeführt.
• Strömungsfeld: Die Ergebnisse zeigen Geschwindigkeitsstagnationszonen am Rand des zentralen Zylinders. In diesen Bereichen, wo die Strömung "eingefangen" ist, erreicht die Temperatur ihr Maximum.
• Wärmequellenverteilung: Die momentane Wärmequellenverteilung zeigt eine sinusförmige Struktur, die gut mit den Erwartungen übereinstimmt.

5. Bedeutung und Ausblick
Das Paper demonstriert die Machbarkeit eines effizienten Multiphysik-Modells für den MSFR, bei dem die rechenintensive Monte-Carlo-Neutronik durch eine reduzierte Spaltmatrix-Methode ersetzt wird, ohne signifikante Genauigkeitsverluste in den thermohydraulischen Ergebnissen. Dies ermöglicht schnellere Iterationen für Designstudien und Sensitivitätsanalysen.

Zukünftige Arbeiten umfassen:

• Sensitivitätsanalysen bezüglich der Gitterauflösung des FM-Modells (feineres vs. gröberes Gitter).
• Untersuchung des Einflusses der Gittergeometrie (z. B. zylindrisches Gitter statt kartesisch) auf die Ergebnisse.
• Weitere Verfeinerung der Kopplung und Validierung.

Das Forschungsvorhaben wurde teilweise durch eine Förderung der US Nuclear Regulatory Commission (NRC) finanziert.