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

이 논문은 고온 용융염 원자로 (MSFR) 의 중성자물리 및 열수력학적 결합을 분석하기 위해 NekRS 와 MOOSE 프레임워크 기반의 Cardinal 코드를 활용하되, OpenMC 대신 빠른 계산 속도와 정확성을 갖춘 핵분열 행렬 (Fission Matrix) 방법을 적용한 다물리 모델링 접근법을 제안합니다.

핵심

🌟 핵심 아이디어: "원자로의 심장을 두 가지 눈으로 보기"

이 연구의 주인공은 **MSFR(용융염 고속로)**이라는 특별한 원자로입니다. 일반적인 원자로는 핵분열을 일으키는 '연료'와 열을 식혀주는 '냉각수'가 따로 있지만, 이 원자로는 연료 자체가 액체 상태의 소금 (용융염) 이 되어 흐릅니다.

마치 뜨거운 국물이 그 자체로 연료가 되어 순환하는 거대한 냄비라고 생각하시면 됩니다. 이 때문에 '핵분열 (중성자)'과 '열/유체 흐름'이 서로 너무 밀접하게 얽혀 있어, 한 가지 방법으로만 분석하면 안 됩니다. 그래서 연구진은 이 복잡한 관계를 풀기 위해 두 가지 다른 시뮬레이션 도구를 함께 사용했습니다.

🛠️ 사용된 도구들 (비유 버전)

1. NekRS (네커스): "유체 흐름의 마라톤 코치"

   • 이 프로그램은 액체 소금 (연료) 이 어떻게 흐르고, 열을 어떻게 전달하는지 아주 정밀하게 계산합니다.
   • 마치 강물 흐름을 분석하는 전문가처럼, 어디에 물이 멈추는지 (저류 구역), 어디가 가장 뜨겁게 달아오르는지 찾아냅니다. 이 연구에서는 최신 그래픽 카드 (GPU) 를 이용해 매우 빠르게 계산했습니다.

2. Fission Matrix Method (분열 행렬법): "핵분열의 예측 카드 뭉치"

   • 보통 원자로의 핵분열을 계산하려면 'OpenMC'라는 무거운 프로그램을 쓰는데, 이번 연구에서는 대신 **'분열 행렬 (FM)'**이라는 더 가볍고 빠른 방법을 썼습니다.
   • 비유하자면: 원자로를 작은 방 (셀) 들로 나누고, "A 방에서 중성자가 튀어 나오면 B 방에서 몇 번의 핵분열이 일어날까?"를 미리 계산해 놓은 거대한 데이터 카드 뭉치를 가지고 있습니다.
   • 연구진은 이 카드 뭉치를 온도 변화에 따라 실시간으로 섞어서 (보간) 사용했습니다. 마치 날씨에 따라 옷을 갈아입는 것처럼, 소금의 온도가 변하면 그에 맞는 '핵분열 카드'를 꺼내 쓰는 방식입니다.

3. Cardinal (카디널): "두 전문가를 연결하는 중재자"

   • NekRS(유체) 와 FM(핵분열) 이 서로 대화할 수 있게 해주는 다리 역할을 합니다.
   • 작동 원리:
     1. **핵분열 전문가 (FM)**가 "여기서 열이 많이 나옵니다!"라고 알려줍니다.
     2. **유체 전문가 (NekRS)**는 그 열을 받아 "아, 그럼 이쪽 물이 더 뜨거워지고 흐르는 속도가 변하겠네요"라고 계산합니다.
     3. 다시 핵분열 전문가는 "오, 온도가 변했으니 핵분열 패턴도 조금 달라지겠네요"라고 수정합니다.
   • 이 과정을 수천 번 반복하며 (피카르 반복), 시스템이 안정된 상태에 도달할 때까지 계속 맞춥니다.

📊 연구 결과: "우리가 맞췄나요?"

연구진은 슈퍼컴퓨터 (Frontier) 를 이용해 이 모델을 실행했습니다. 결과는 매우 만족스러웠습니다.

• 온도와 흐름: 액체 소금이 흐르면서 가장 뜨거운 곳은 물이 멈추는 '고인 물' 같은 곳이었습니다. 유체가 잘 순환하지 않는 구석구석에서 열이 쌓여 온도가 최고조에 달하는 것을 발견했습니다.
• 정확도 검증: 이 새로운 방법 (FM 사용) 으로 계산한 결과와, 기존에 쓰던 무거운 방법 (OpenMC 사용) 으로 계산한 결과를 비교했습니다.
  • 결과: 두 방법의 온도 차이가 0.1 도밖에 나지 않았습니다. (예를 들어, 997.3 도 vs 997.4 도)
  • 의미: 무거운 프로그램을 쓰지 않아도, 훨씬 가볍고 빠르게 거의 똑같은 정확한 결과를 얻을 수 있다는 뜻입니다.

💡 결론 및 향후 계획

이 논문은 **"복잡한 원자로 시스템을 분석할 때, 무거운 총 대신 정교한 사냥개 (FM 방법) 를 쓰면 훨씬 효율적이다"**라는 것을 증명했습니다.

• 현재: 유체 흐름과 핵분열이 서로 영향을 주고받는 과정을 성공적으로 모사했습니다.
• 앞으로: 계산에 사용한 '방 (격자)'의 크기와 모양을 더 세밀하게 다듬어, 어떤 모양이 가장 정확한지 실험해 볼 계획입니다.

🎯 한 줄 요약

"액체 연료 원자로의 뜨거운 흐름과 핵분열을, 무거운 계산 없이 가볍고 빠르게, 하지만 정확하게 예측하는 새로운 시뮬레이션 기술을 개발했습니다."

기술 요약

논문 요약: NekRS 와 핵분열 행렬 (Fission Matrix) 방법을 활용한 용융염 고속로 (MSFR) 의 다물리 모델링

1. 연구 배경 및 문제 제기 (Problem)

• 배경: 용융염 원자로 (MSR) 는 제 4 세대 국제 포럼 (Gen-IV) 에서 안전성, 신뢰성, 폐기물 감소, 경제성 등의 이유로 선정된 차세대 원자로 개념입니다. 특히 흑연 감속재를 제거한 **용융염 고속로 (MSFR)**는 냉각재가 동시에 연료 역할을 하므로, 연료가 1 차 계통을 순환하며 이동하는 특성상 중성자 물리 (Neutronics) 와 열유체 (Thermal-hydraulics) 간의 결합이 매우 강력합니다.
• 문제: 이러한 강한 결합을 분석하기 위해서는 정밀한 다물리 (Multiphysics) 모델링이 필수적입니다. 기존에 사용되던 고충실도 코드인 OpenMC(몬테카를로 입자 수송 코드) 를 사용하는 방식은 계산 비용이 매우 높습니다. 따라서, 중성자 물리 방정식을 풀기 위해 차원 축소 모델 (Reduced-order model) 인 핵분열 행렬 (Fission Matrix, FM) 방법을 도입하여 계산 효율성을 높이면서도 정밀한 다물리 해석을 수행하는 새로운 접근법이 필요합니다.

2. 방법론 (Methodology)

이 연구는 MOOSE 프레임워크 기반의 오픈소스 다물리 코드인 Cardinal을 활용하여 MSFR 의 중성자 물리 - 열유체 결합 모델을 개발했습니다.

• 주요 도구 및 역할:

  • NekRS: GPU 가속 기반의 스펙트럴 요소 방법 (SEM) 을 사용하는 CFD 코드입니다. 비압축성 Navier-Stokes 방정식과 에너지 방정식을 풀며, RANS $k-\tau$ 난류 모델을 사용하여 유동 및 열전달을 시뮬레이션합니다.
  • Cardinal: MOOSE 프레임워크 내에서 NekRS 와 중성자 물리 코드를 결합하는 래퍼 (Wrapper) 입니다.
  • 핵분열 행렬 (FM) 방법: 기존 OpenMC 대신 FM 방법을 사용하여 중성자 물리 문제를 해결합니다.
    • 핵분열 행렬 (FMDB): 몬테카를로 코드인 Serpent 를 사용하여 다양한 연료 온도 (800K ~ 1200K) 조건에서 3D FM 데이터베이스 (FMDB) 를 생성했습니다.
    • 선형 보간: 실제 시뮬레이션 중의 국소 온도에 맞춰 인접한 두 개의 FMDB 를 선형 보간하여 핵분열 행렬 $A$를 실시간으로 구성합니다.
    • 해석: 생성된 행렬을 MOOSE 내에서 C++ 로 구현된 커스텀 모듈로 풀어 유효증배인자 ($k_{eff}$) 와 핵분열 소스 분포를 구합니다.

• 모델 구성:

  • 기하학적 모델: MSFR 1 차 루프의 3D 기하학 (Gmsh 생성) 을 사용하며, 펌프나 열교환기는 생략하고 연속적인 입/출구로 단순화했습니다.
  • 메쉬: NekRS 는 약 8,600 만 개의 가우스 - 로바토 - 르장드르 (GLL) 점으로 구성된 135 만 개의 육면체 요소 메쉬를 사용했습니다. FM 모델은 13x13x13 의 직교 격자를 사용하며, NekRS 메쉬 내부에 매핑됩니다.
  • 결합 프로세스 (Cardinal Coupling):
    1. FM 모델이 초기 온도로 핵분열 소스를 계산하고 열원 ($q'''_{fis}$) 으로 변환합니다.
    2. Cardinal 이 이 열원을 NekRS 로 전달합니다.
    3. NekRS 가 유동 및 온도장을 계산하여 결과를 Cardinal 로 되돌려줍니다.
    4. FM 모델이 업데이트된 온도를 받아 행렬을 재보간하고 과정을 반복 (Picard iteration) 하여 수렴합니다.

3. 주요 기여 (Key Contributions)

1. 고성능 다물리 결합 프레임워크 구축: MOOSE/Cardinal 환경에서 NekRS(CFD) 와 FM 방법 (중성자 물리) 을 성공적으로 결합한 새로운 워크플로우를 제시했습니다.
2. 계산 효율성 향상: 고비용의 몬테카를로 입자 수송 (OpenMC) 대신 FM 방법을 사용하여 중성자 물리 계산을 수행함으로써, 다물리 시뮬레이션의 계산 부하를 줄이면서도 물리적 정확도를 유지하는 방안을 검증했습니다.
3. 온도 의존적 FM 데이터베이스 (FMDB) 활용: 다양한 온도 조건에서 생성된 FMDB 를 실시간으로 선형 보간하여 국소 온도 변화에 따른 중성자 물리 피드백을 정밀하게 반영하는 기법을 적용했습니다.

4. 결과 (Results)

• 시뮬레이션 환경: Frontier 슈퍼컴퓨터 (GPU: NekRS, CPU: FM 모델) 에서 수행되었습니다.
• 성능 지표:
  • 계산된 유효증배인자 ($k_{eff}$): 1.04632
  • NekRS 의 outlet 온도: 997.4 K, 평균 온도: 973.0 K
• 검증 및 비교:
  • 기존 연구 [6] 의 NekRS-OpenMC 모델 결과와 비교했습니다.
  • Outlet 온도 차이: 0.1 K (997.3 K vs 997.4 K)
  • 평균 온도 차이: 10.1 K (983.1 K vs 973.0 K)
  • 전체 에너지가 동일하므로 온도 분포의 미세한 차이는 예상된 범위 내이며, 열원 분포, 속도장, 온도장의 분포 패턴은 기존 모델과 우수한 일치를 보였습니다.
• 물리적 현상 관찰:
  • 열원 분포는 사인파 (sinusoidal) 형태를 띠었습니다.
  • 중앙 실린더 주변의 저속 영역 (Stagnation regions) 에서 유체가 갇혀 계속 가열됨에 따라 최고 온도가 발생하는 것을 확인했습니다.

5. 의의 및 결론 (Significance & Conclusion)

• 의의: 본 연구는 MSFR 과 같은 복잡한 다물리 시스템에 대해, 고비용의 몬테카를로 시뮬레이션 없이도 FM 방법을 통해 효율적이고 정확한 다물리 해석이 가능함을 입증했습니다. 이는 차세대 원자로 설계 및 안전성 평가에 있어 계산 자원을 절감하면서도 신뢰할 수 있는 도구를 제공한다는 점에서 중요합니다.
• 결론: NekRS 와 FM 방법을 결합한 Cardinal 기반 모델은 열유체 및 중성자 물리 피드백을 성공적으로 구현하여 물리적 거동을 잘 재현했습니다.
• 향후 과제: FM 모델의 메쉬 민감도 분석 (더 굵거나 미세한 메쉬 테스트) 과 원통형 메쉬를 사용한 격자 형태가 결과에 미치는 영향 분석이 추가로 수행될 예정입니다.

이 논문은 차세대 원자로 개발을 위한 고성능 컴퓨팅 (HPC) 기반의 다물리 시뮬레이션 기술의 발전 가능성을 보여주는 중요한 사례입니다.