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

Este artigo propõe um modelo computacional multiphysics de alta fidelidade para o Reator de Sal Fundido de Nêutrons Rápidos (MSFR) utilizando o código Cardinal e o Método da Matriz de Fissão para resolver as equações de neutrônica de forma eficiente e acoplada à hidrodinâmica computacional.

O essencial

Imagine que você está tentando prever o clima dentro de uma panela de pressão futurista e supercomplexa, que não apenas cozinha comida, mas gera energia elétrica para uma cidade inteira. Essa "panela" é o Reator de Sal Fundido (MSFR), um tipo de usina nuclear de nova geração que usa sal derretido como combustível e refrigerante ao mesmo tempo.

O desafio? O sal se move, esquenta, esfria e muda de comportamento o tempo todo. Se você tentar calcular tudo de uma vez só, o computador fica tão sobrecarregado que parece uma tartaruga tentando correr uma maratona.

É aqui que entra o trabalho dos pesquisadores da Universidade da Pensilvânia (Maximiliano, Elia e Saya). Eles criaram um "super time" de softwares para simular esse reator, mas com um truque inteligente. Vamos explicar como funciona usando analogias do dia a dia:

1. O Grande Time de Software (Cardinal, NekRS e MOOSE)

Pense no projeto como uma orquestra musical.

• MOOSE é o maestro. Ele garante que todos os músicos toquem juntos e na hora certa.
• NekRS é o músico especialista em fluidos (o "flautista"). Ele é super rápido (usa chips de videogame modernos, chamados GPUs) e simula como o sal derretido flui, gira e transfere calor dentro do reator.
• Cardinal é o maestro que conecta o flautista ao resto da banda.

2. O Truque da "Fotografia Rápida" (O Método da Matriz de Fissão)

Aqui está a parte mais genial e simplificada do trabalho.

Normalmente, para saber onde a energia nuclear está sendo gerada, os cientistas usam um método chamado "Monte Carlo". Imagine que você quer saber onde estão as gotas de chuva em um furacão. O método Monte Carlo simula cada gota individualmente. É extremamente preciso, mas demorado demais para fazer isso em tempo real enquanto o sal está se movendo.

Os pesquisadores trocaram esse método lento por algo chamado Método da Matriz de Fissão.

• A Analogia: Em vez de simular cada gota de chuva, eles tiraram "fotos" (ou mapas) de como a chuva cai em diferentes temperaturas. Eles criaram um banco de dados com 5 fotos diferentes (uma para 800°C, outra para 900°C, etc.).
• O Truque: Quando o reator está em uma temperatura que não está nas fotos (digamos, 950°C), o computador não precisa recalcular tudo. Ele apenas olha para a foto de 900°C e a de 1000°C e interpola (faz uma média matemática simples) para saber o que acontece em 950°C.
• É como se você estivesse no meio de um filme e, em vez de renderizar cada quadro novo do zero, o computador apenas ajustasse levemente o quadro anterior. Isso torna o cálculo muito mais rápido.

3. A Dança entre Calor e Movimento (Acoplamento Multiphysics)

O reator é um sistema de "dança de casais":

1. O Sal (NekRS) flui e carrega calor.
2. O Núcleo (Método da Matriz) gera calor baseado na temperatura do sal.
3. Eles trocam informações constantemente.

O software faz isso em um ciclo:

• O "Maestro" pede ao "Flautista" (NekRS): "Como está o sal?"
• O "Flautista" responde: "Está quente aqui e frio ali".
• O "Maestro" passa essa info para o "Mapa de Fissão": "Ok, ajuste a geração de energia baseada nisso".
• O "Mapa" diz: "Aqui vai gerar mais calor".
• O "Flautista" recalcula o fluxo do sal com esse novo calor.
• E assim por diante, até que tudo se estabilize.

4. O Que Eles Descobriram?

Ao rodar essa simulação em um supercomputador gigante (Frontier), eles viram coisas interessantes:

• Zonas de "Trânsito Parado": Em algumas áreas do reator, o sal fica quase parado (como um engarrafamento no trânsito). Nessas zonas, o calor fica preso e a temperatura sobe muito, atingindo o máximo.
• Precisão: Eles compararam seu novo método "rápido" (Matriz de Fissão) com o método "lento e tradicional" (Monte Carlo). O resultado foi impressionante: as temperaturas e velocidades calculadas foram quase idênticas (diferença de apenas 0,1 grau em alguns pontos).

Resumo Final

Essa pesquisa é como criar um simulador de voo para uma usina nuclear. Em vez de esperar dias para calcular como o reator vai se comportar, eles criaram um método que usa "mapas prévios" e matemática inteligente para obter resultados rápidos e precisos.

Isso é crucial porque, para construir reatores de sal fundido no futuro que sejam seguros e econômicos, os engenheiros precisam testar milhares de cenários em computadores. Fazer isso em horas, em vez de meses, é o que vai permitir que essa tecnologia saia do papel e se torne realidade.

Em suma: Eles ensinaram o computador a "adivinhar" a física nuclear de forma inteligente e rápida, permitindo que o reator seja desenhado com mais segurança e eficiência.

Resumo técnico

Resumo Técnico: Modelagem Multifísica do Reator de Sal Fundido Rápido (MSFR) usando NekRS e o Método da Matriz de Fissão

1. Problema e Contexto

O Reator de Sal Fundido (MSR) foi selecionado pelo Fórum Internacional de Geração-IV como um conceito de referência devido à sua segurança, redução de resíduos e eficiência econômica. Uma variação inovadora, o Reator de Sal Fundido Rápido (MSFR), elimina o moderador de grafite, operando como um reator rápido. Uma característica crítica do MSFR é que o refrigerante também atua como combustível, o que cria um acoplamento extremamente forte entre a neutrônica e a hidrodinâmica térmica, uma vez que o combustível circula continuamente pelo sistema primário.

O desafio central abordado neste trabalho é desenvolver um modelo computacional de alta fidelidade capaz de resolver esse acoplamento multifísica de forma eficiente. Tradicionalmente, códigos de transporte de Monte Carlo (como o OpenMC) são usados para a neutrônica, mas podem ser computacionalmente custosos para simulações acopladas em tempo real. O objetivo deste estudo é substituir o solucionador de neutrônica tradicional por um Modelo de Ordem Reduzida (ROM) baseado no Método da Matriz de Fissão (FM) para acelerar o processo, mantendo a precisão necessária.

2. Metodologia

A pesquisa utiliza o framework Cardinal, que integra o código de Dinâmica dos Fluidos Computacional (CFD) NekRS e o framework de simulação multifísica MOOSE.

• NekRS (Hidrodinâmica): Um código CFD de alta fidelidade acelerado por GPU, utilizando o Método dos Elementos Espectrais (SEM). No estudo, o NekRS resolve as equações de Navier-Stokes incompressíveis e de energia em regime transiente, utilizando o modelo de turbulência RANS $k-\tau$. A geometria do MSFR é simplificada (sem bombas ou trocadores de calor detalhados), com malha de $1,35 \times 10^6$ elementos hexaédricos.
• Método da Matriz de Fissão (Neutrônica): Em vez de usar o OpenMC em tempo de execução, os autores implementaram o Método da Matriz de Fissão (FM) dentro do MOOSE.
  • Banco de Dados de Matrizes de Fissão (FMDBs): Matrizes de fissão foram pré-calculadas usando o código Monte Carlo Serpent para cinco temperaturas uniformes de combustível (800 K a 1200 K) em uma discretização cartesiana de $13 \times 13 \times 13$.
  • Interpolação Linear: Durante a simulação, a matriz de fissão $A$ para a temperatura local do combustível é obtida por interpolação linear entre as FMDBs mais próximas.
  • Resolução: O sistema resolve a equação de autovalor ($F = \frac{1}{k_{eff}} A F$) para obter a distribuição da fonte de fissão e o fator de multiplicação ($k_{eff}$).
• Acoplamento (Cardinal): O esquema de acoplamento utiliza iterações de Picard no tempo. O modelo FM calcula a fonte de calor volumétrica baseada na temperatura atual; o NekRS recebe essa fonte, calcula o novo campo de temperatura e velocidade do fluido, e retorna a temperatura atualizada ao modelo FM. O processo repete-se até a convergência de estados estacionários.

3. Contribuições Chave

• Implementação do Método da Matriz de Fissão no MOOSE: Desenvolvimento de novas funções em C++ dentro do framework MOOSE para carregar FMDBs, realizar interpolação linear de matrizes e resolver o sistema de neutrônica reduzida.
• Integração NekRS-FM: Demonstração bem-sucedida do acoplamento entre um solucionador CFD de alta fidelidade (NekRS) e um solucionador de neutrônica de ordem reduzida (FM) dentro do ecossistema Cardinal, substituindo o OpenMC para fins de eficiência computacional.
• Validação Preliminar: Comparação direta dos resultados do novo modelo (NekRS-FM) com um modelo de referência existente (NekRS-OpenMC), validando a abordagem de ordem reduzida.

4. Resultados

As simulações foram executadas no supercomputador Frontier (utilizando GPUs para o NekRS e CPUs para o modelo FM). Os resultados preliminares indicam:

• Fator de Multiplicação: O modelo FM obteve um $k_{eff}$ de 1.04632.
• Convergência Térmica: A simulação foi executada por aproximadamente 60 unidades convectivas, com média temporal calculada sobre 80 unidades adicionais.
• Temperaturas:
  • Temperatura de saída (NekRS-FM): 997.4 K.
  • Temperatura média (NekRS-FM): 973.0 K.
  • Comparação com NekRS-OpenMC: A diferença na temperatura de saída foi de apenas 0.1 K (997.3 K no modelo de referência), enquanto a temperatura média apresentou uma diferença de 10.1 K (983.1 K no modelo de referência).
• Comportamento do Fluxo: Identificou-se a formação de zonas de estagnação na periferia do cilindro central, onde a velocidade é baixa e a temperatura atinge seus máximos devido ao aquecimento contínuo do fluido "preso".
• Distribuição de Fonte de Calor: A distribuição da fonte de calor apresentou um formato senoidal, consistente com o modelo de referência.

5. Significado e Conclusões

O estudo demonstra que o Método da Matriz de Fissão é uma alternativa viável e eficiente ao Monte Carlo completo para simulações multifísica de reatores de sal fundido em tempo real. A pequena discrepância nos resultados térmicos em comparação com o modelo OpenMC confirma que a abordagem de ordem reduzida captura a física essencial do sistema, especialmente a conservação de energia total.

A principal vantagem é a redução potencial do custo computacional, permitindo simulações mais rápidas e possivelmente a exploração de cenários operacionais complexos que seriam proibitivos com o Monte Carlo acoplado. Trabalhos futuros incluirão análises de sensibilidade da malha do modelo FM (testando malhas mais finas e grossas) e a investigação do impacto da forma da malha (transição de cartesiana para cilíndrica) nos resultados.

Este trabalho representa um passo significativo na direção de ferramentas de simulação de próxima geração para o projeto e análise segura de reatores de geração IV.