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

Este artículo propone un modelo computacional multiphysics para el Reactor de Sales Fundidas de Espectro Rápido (MSFR) que integra el código de alta fidelidad Cardinal (basado en MOOSE y NekRS) con el Método de la Matriz de Fisión para resolver las ecuaciones de neutrones de manera rápida y precisa, aprovechando la precalibración mediante simulaciones Monte Carlo.

Lo esencial

¡Claro que sí! Imagina que este artículo científico es como la receta y el plan de construcción para un motor de coche futurista que no usa gasolina, sino sal fundida (sal derretida) que también actúa como combustible.

Aquí tienes la explicación de cómo funciona este proyecto, usando analogías sencillas:

1. El Gran Reto: Un Motor que se Mueve y Piensa a la vez

Normalmente, en un coche, el combustible se queda quieto en el tanque y el motor lo quema. Pero en este reactor de Sal Fundida, el combustible es líquido y circula por todo el sistema como la sangre en nuestro cuerpo.

• El problema: Como el combustible se mueve, la temperatura cambia constantemente, y eso afecta a cuánta energía nuclear se produce. Es como intentar conducir un coche mientras el motor se calienta y se enfría al mismo tiempo que cambia de forma.
• La solución: Los científicos necesitan un "cerebro" digital que pueda calcular dos cosas a la vez:
  1. Física de fluidos: ¿Cómo fluye la sal caliente? (Como el agua en una tubería).
  2. Neutrónica: ¿Cuánta energía nuclear se está creando? (Como contar chispas).

2. Las Herramientas del Taller: NekRS y la "Matriz de Fisión"

Para resolver este rompecabezas, usaron dos herramientas principales:

• NekRS (El Simulador de Fluidos): Imagina que es un videojuego de física súper avanzado que corre en superordenadores potentes (con tarjetas gráficas de video). Su trabajo es simular cómo se mueve la sal, dónde se calienta y dónde se enfría. Es muy preciso, como una cámara de alta velocidad viendo el agua correr.
• El Método de la Matriz de Fisión (El "Atajo Inteligente"): Aquí está la parte genial. Normalmente, para calcular la energía nuclear, los ordenadores tienen que simular billones de partículas individuales (como contar cada grano de arena de una playa). Eso tarda muchísimo.
  • La analogía: En lugar de contar cada grano de arena, los científicos crearon un "mapa de calor" precalculado (llamado Base de Datos de Matrices). Es como tener un manual que dice: "Si la sal está a 800 grados, la energía se comporta así; si está a 900, así".
  • Cuando el simulador de fluidos (NekRS) dice "¡La temperatura subió!", el sistema nuclear no vuelve a contar todo desde cero. Solo mira en su manual, hace una pequeña interpolación (una estimación rápida entre dos valores) y dice: "Ah, entonces la energía será X". Esto hace que el cálculo sea muchísimo más rápido.

3. El Baile de la Información (Acoplamiento)

El sistema funciona como un bailarín que sigue el ritmo de otro:

1. Paso 1: El sistema nuclear (la Matriz) dice: "Aquí hay mucha energía, calienta la sal".
2. Paso 2: NekRS (el fluido) recibe esa caloría, mueve la sal y dice: "¡Oye, la sal se movió y ahora está más caliente aquí y más fría allá!".
3. Paso 3: NekRS le pasa esa nueva temperatura al sistema nuclear.
4. Paso 4: El sistema nuclear actualiza su cálculo de energía basado en esa nueva temperatura.

Repiten este paso una y otra vez (miles de veces) hasta que todo se estabiliza, como cuando mezclas leche y café hasta que el color es uniforme.

4. ¿Qué Descubrieron? (Los Resultados)

Al final, probaron este nuevo método rápido en un superordenador llamado "Frontier" y compararon los resultados con el método antiguo (el que cuenta cada partícula, que es lento pero preciso).

• El hallazgo: ¡Funcionó casi igual de bien!
  • La temperatura de salida fue casi idéntica (una diferencia de apenas 0.1 grados, ¡como la diferencia entre un día de verano y uno de verano con una brisa!).
  • Encontraron "zonas de estancamiento": Hay partes del reactor donde la sal se mueve muy lento (como un remanso en un río). En esas zonas, la sal se queda "atrapada" y se calienta mucho más que en el resto. El sistema detectó esto perfectamente.

5. ¿Por qué es importante?

Este trabajo es como crear un prototipo de coche eléctrico que se puede probar en una computadora en minutos en lugar de meses.

• Antes: Simular este reactor tomaba mucho tiempo y recursos.
• Ahora: Con este "atajo inteligente" (la Matriz de Fisión), pueden simular el reactor mucho más rápido.
• El futuro: Esto ayuda a diseñar reactores más seguros y eficientes, que podrían darnos energía limpia sin producir tanto residuo nuclear como los reactores actuales.

En resumen: Crearon un sistema digital que combina un simulador de fluidos de alta velocidad con un "manual de instrucciones" inteligente para la energía nuclear. Esto les permite entender cómo se comporta un reactor de sal fundida en tiempo real, encontrando puntos calientes y asegurando que todo funcione de manera segura y eficiente.

Resumen técnico

Resumen Técnico: Modelado Multiphysics del Reactor de Sal Fundida de Rápido (MSFR) usando NekRS y el Método de la Matriz de Fisión

1. Planteamiento del Problema

El Reactor de Sal Fundida (MSR) fue seleccionado por el Foro Internacional de Generación IV como un concepto de referencia debido a su seguridad, reducción de residuos y eficiencia económica. Una variante innovadora, el Reactores de Sal Fundida de Rápido (MSFR), elimina el moderador de grafito, funcionando como un reactor reproductor rápido donde el refrigerante es también el combustible.

Esta característica crea un acoplamiento estrecho entre la neutrocinética y la hidráulica térmica, ya que el combustible circula a través del sistema primario. Desarrollar modelos computacionales precisos para el MSFR es complejo y requiere un enfoque multiphysics. El desafío principal abordado en este trabajo es resolver las ecuaciones de neutrocinética de manera eficiente dentro de un entorno de alta fidelidad, evitando el alto costo computacional de los métodos de transporte de Monte Carlo tradicionales en cada paso de tiempo.

2. Metodología

El estudio propone un modelo computacional acoplado neutrocinético-hidráulico térmico utilizando el marco de software Cardinal (dentro del entorno MOOSE), pero con una modificación clave en la resolución de la neutrocinética:

• Hidráulica Térmica (NekRS): Se utiliza el código NekRS, un código de Dinámica de Fluidos Computacional (CFD) acelerado por GPU que emplea el Método de Elementos Espectrales (SEM). Resuelve las ecuaciones de Navier-Stokes incompresibles y de energía no estacionarias utilizando el modelo de turbulencia RANS $k-\tau$.
  • Geometría: Se modela un bucle primario simplificado (sin bombas ni intercambiadores) con una separación de 50 cm entre entrada y salida.
  • Condiciones: Temperatura de entrada de 898 K, paredes adiabáticas y perfiles de velocidad parabólicos.
• Neutrocinética (Método de la Matriz de Fisión - FM): En lugar de usar el código Monte Carlo OpenMC (como hace Cardinal por defecto), se implementa el Método de la Matriz de Fisión (FM).
  • Base de Datos (FMDB): Se generaron cinco bases de datos de matrices de fisión utilizando el código Monte Carlo Serpent para diferentes temperaturas uniformes de combustible (800 K a 1200 K). La geometría del núcleo es cilíndrica con un manto fértil y un reflector.
  • Interpolación: Para representar el estado real del reactor con una distribución de temperatura no uniforme, se utiliza una interpolación lineal entre las matrices de las bases de datos más cercanas a la temperatura local de cada celda.
  • Implementación: Se desarrollaron funciones en C++ dentro de MOOSE para cargar las bases de datos, realizar la interpolación, resolver el problema de autovalores ($k_{eff}$ y distribución de fuente) y calcular la fuente de calor volumétrica.
• Acoplamiento (Cardinal): Se utiliza un esquema de iteración de Picard "en el tiempo".
  1. El modelo FM calcula la distribución de fuente de fisión y la escala para obtener la fuente de calor ($q'''_{fis}$).
  2. Cardinal envía esta fuente a NekRS.
  3. NekRS realiza múltiples pasos de tiempo ($N=10000$) para calcular la temperatura y velocidad del fluido.
  4. La temperatura del fluido se devuelve al modelo FM para actualizar la matriz de fisiones.
  5. El proceso se repite hasta alcanzar el estado estacionario (convergencia de temperaturas).

3. Contribuciones Clave

• Desarrollo de un Modelo Reducido de Orden (ROM): La principal contribución es la integración exitosa del Método de la Matriz de Fisión dentro del marco MOOSE/Cardinal para acoplarlo con CFD de alta fidelidad (NekRS). Esto permite resolver la neutrocinética de manera mucho más rápida que con Monte Carlo directo en cada iteración de tiempo.
• Validación de Interpolación de Bases de Datos: Se demuestra la viabilidad de interpolar matrices de fisión basadas en perfiles de temperatura uniformes para modelar distribuciones de temperatura complejas y no uniformes en un flujo de sal fundida.
• Implementación en GPU: Se logra ejecutar simulaciones en el superordenador Frontier, utilizando GPUs para la parte de CFD (NekRS) y CPUs para el modelo de neutrocinética (FM), demostrando la escalabilidad del enfoque.

4. Resultados

Las simulaciones preliminares se ejecutaron durante aproximadamente 60 unidades convectivas, seguidas de un promediado temporal de 80 unidades más.

• Parámetros Neutrocinéticos: El factor de multiplicación ($k_{eff}$) obtenido fue 1.04632.
• Distribución de Calor: La fuente de calor instantánea mostró una forma sinusoidal, consistente con la física del reactor.
• Hidráulica y Temperatura:
  • Se identificaron regiones de estancamiento en la periferia del cilindro central donde la velocidad es baja y la temperatura alcanza su máximo debido al "atrapamiento" del flujo.
  • La temperatura de salida calculada fue de 997.4 K y la temperatura promedio de 973.0 K.
• Comparación con Referencia: Al comparar con un modelo previo que usaba NekRS-OpenMC:
  • La temperatura de salida difiere solo en 0.1 K (997.3 K vs 997.4 K).
  • La temperatura promedio difiere en 10.1 K (983.1 K vs 973.0 K).
  • Se concluye que la diferencia es pequeña y esperable, dado que la energía total es la misma en ambos modelos, y las distribuciones de velocidad, temperatura y fuente de calor muestran un buen acuerdo cualitativo.

5. Significado y Conclusiones

El trabajo demuestra que es viable desarrollar un modelo multiphysics de alta fidelidad para el MSFR utilizando el Método de la Matriz de Fisión como aproximación de neutrocinética. Este enfoque ofrece una alternativa computacionalmente eficiente a los métodos de Monte Carlo completos para simulaciones transitorias o de estado estacionario acopladas con CFD.

Los resultados preliminares confirman que el modelo captura correctamente los fenómenos físicos clave, como las zonas de estancamiento térmico. El estudio sienta las bases para análisis futuros que incluirán:

• Análisis de sensibilidad sobre la malla del modelo FM (prueba de mallas más finas y gruesas).
• Evaluación de la influencia de la forma de la malla (comparando mallas cartesianas con mallas cilíndricas) para optimizar la discretización y la precisión del modelo.

Este enfoque podría acelerar significativamente el diseño y la optimización de reactores de sal fundida, permitiendo simulaciones más rápidas sin sacrificar la precisión física crítica.