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

本論文では、燃料と冷却材が同一である溶融塩高速炉の解析に向け、NekRS と Fission Matrix 法を統合した Cardinal フレームワークを用いて、中性子物理学と熱水力学的な結合を効率的に解く低次モデルを開発することを提案している。

要約

🍳 1. 何を作ろうとしているの？（溶融塩炉とは）

まず、この研究の対象である「溶融塩炉」は、従来の原子炉とは少し違います。
普通の原子炉は、燃料を固形のまま入れて、別の液体（冷却水）で冷やしますが、この炉は**「燃料そのものが液体（溶けた塩）」**になっています。

• 例え話：
  普通の原子炉が「鍋の中で固形のステーキを焼いて、周りに水を流して冷やす」のに対し、溶融塩炉は**「鍋の中身そのものがスープになっていて、そのスープが燃料になっている」**ようなものです。
  このスープ（燃料）がポンプでぐるぐる循環するため、熱の動きと核反応（エネルギーを作る反応）が密接に絡み合っています。そのため、設計には「熱の計算」と「核反応の計算」を同時にやる必要があります。

🧩 2. どうやって計算しているの？（2 つの役割分担）

この研究では、2 つの異なるコンピュータープログラムを組み合わせる「マルチフィジックス（多物理場）」という手法を使っています。

A. 熱と流れの計算：NekRS（ネックス）

• 役割： 燃料がどう流れて、どこが熱くなっているかを計算します。
• 例え話： これは**「流体力学の天才シェフ」**です。スープが鍋の中でどう渦を巻くか、どこで止まって熱くなりすぎるか、壁にどれくらい熱が伝わるかを、GPU（高性能な画像処理チップ）を使ってリアルタイムにシミュレーションします。

B. 核反応の計算：Fission Matrix（ファイション・マトリックス）

• 役割： 燃料の中でどこで核分裂が起きているか、エネルギーがどれくらい出ているかを計算します。
• 従来の方法： 以前は「OpenMC」というプログラムを使って、粒子一つ一つをシミュレーションしていました。これは**「粒子を一つずつ数える meticulous（几帳面な）学者」**のようなもので、正確ですが非常に時間がかかります。
• 今回の新手法（Fission Matrix）： 今回は**「Fission Matrix（核分裂行列）」**という方法を使っています。
  • 例え話： これは**「地図と予測表」**のようなものです。炉を小さなマス目（セル）に分け、「A 地点で核分裂が起きると、B 地点でどれくらい反応が起きるか」を事前に計算した巨大な表（行列）を持っています。
  • シミュレーション中は、この表を参照して「今の温度なら、この表を少し変形させて使えばいい」という**「簡易的な予測」**を瞬時に行います。これにより、計算が劇的に速くなります。

🤝 3. 2 つのプログラムはどう連携するの？（Cardinal という仲介者）

この 2 つのプログラムは、**「Cardinal（カーディナル）」**という仲介役のソフトを通じて会話しています。

• 連携のプロセス（ピカール反復法）：
  1. 核反応チーム（Fission Matrix）： 「今の温度なら、ここからエネルギーが出ますよ」と熱の源（ヒートソース）を伝えます。
  2. 流体チーム（NekRS）： 「なるほど、その熱が加わると、スープはこう流れて、この部分が 997.4 度になります」と温度と流速を計算して返します。
  3. ループ： 核反応チームは「あ、温度が上がったから、反応の強さを変えよう」と調整し、また流体チームに伝えます。
     この「言い合い」を何千回も繰り返して、温度も流れも核反応も安定する（収束する）状態を見つけ出します。

📊 4. 結果はどうだった？

• 成功： この新しい「Fission Matrix」を使った方法は、従来の「粒子を一つずつ数える方法」と比べて、ほぼ同じ精度で、はるかに速く計算できました。
• 発見： シミュレーションの結果、炉の中心部分の周りに「流れが止まっている場所（停滞域）」があることが分かりました。そこは熱が逃げずに**「お湯が溜まって熱くなりすぎやすい場所」**でした。
• 比較： 従来の方法と比べて、出口の温度は 0.1 度しか違いませんでした。これは「同じエネルギーを使っているのだから、当然の結果だ」と言えます。

🎯 まとめ

この論文は、**「溶融塩炉という複雑な料理を、より速く、より安く設計するために、新しい計算方法（Fission Matrix）を開発し、成功させた」**という報告です。

• 従来の方法： 粒子を一つずつ数えて、何日もかかる精密な計算。
• 今回の方法： 事前に作った「予測表」を賢く使いながら、瞬時に答えを出すスマートな計算。

この技術が完成すれば、将来の安全で効率的な原子力発電所の設計が、もっと簡単になることが期待されています。

技術概要

論文要約：NekRS と分岐行列法を用いた溶融塩高速炉（MSFR）のマルチフィジックスモデリング

1. 背景と課題 (Problem)

溶融塩炉（MSR）は、安全性、信頼性、廃棄物削減、経済性などの優位性から、第 4 世代原子炉国際フォーラム（GIF）によって参照炉概念の一つに選定されました。その中でも、黒鉛減速材を排除した「溶融塩高速炉（MSFR）」は、冷却材が同時に燃料としても機能する特徴を持っています。
この設計において、燃料が一次系を循環するため、中性子工学（中性子輸送）と熱水力（流体・熱伝達）の結合が極めて密接になります。従来の高忠実度シミュレーションでは、モンテカルロ粒子輸送コード（OpenMC など）を用いて中性子計算を行っていましたが、計算コストが高く、動的な熱水力フィードバックとの結合において課題がありました。
本研究の目的は、中性子方程式を低次元モデル（Reduced-Order Model）である「分岐行列法（Fission Matrix Method）」で解くことにより、NekRS（CFD）と効率的に結合した MSFR のマルチフィジックス計算モデルを開発することです。

2. 手法 (Methodology)

2.1 使用ソフトウェアとフレームワーク

• Cardinal: MOOSE フレームワーク内で動作するオープンソースアプリケーション。NekRS（CFD）と中性子計算コードを結合し、高忠実度のマルチフィジックスシミュレーションを可能にします。
• NekRS: GPU 加速型のスペクトル要素法（SEM）に基づく CFD コード。非定常非圧縮性 Navier-Stokes 方程式とエネルギー方程式を解きます。本研究では RANS（レイノルズ平均 Navier-Stokes）の $k-\tau$ 乱流モデルを使用しました。
• 分岐行列法（FM 法）: 従来のモンテカルロ法（OpenMC）の代わりに採用された中性子計算手法。
  • 炉心を空間セルに離散化し、$n \times n$ の行列 $A$（分岐行列）を定義します。
  • 行列の固有ベクトルが分岐源分布、固有値が実効増倍率（$k_{eff}$）となります。
  • 燃料温度分布に応じて、事前計算された「分岐行列データベース（FMDBs）」から線形補間を行い、行列 $A$ を構築します。

2.2 計算モデルの構成

1. 熱水力モデル (NekRS):
   • 3 次元幾何学（Gmsh 生成）を使用。ポンプや熱交換器は簡略化され、連続的な入口/出口としてモデル化されています。
   • メッシュ：135 万個の六面体要素、多項式次数 3（約 8600 万の GLL 積分点）。
   • 境界条件：入口温度 898 K、断熱壁、入口で放物線状の速度プロファイル。
2. 中性子モデル (FM 法):
   • 分岐行列データベース（FMDBs）は、Serpent コードを用いて生成されました（温度 800K〜1200K の 5 パターン）。
   • MOOSE 内で C++ 関数として実装され、FMDB の読み込み、温度補間、行列の求解、分岐源分布の取得を行います。
   • メッシュ：13×13×13 の直交格子（FMDB と同じ解像度）。
3. 結合アルゴリズム (Cardinal):
   • ピカルド反復法（時間内反復）を使用。
   • 手順：
     1. FM モデルが初期温度分布で分岐源を計算し、体積発熱源 $q'''_{fis}$ を算出。
     2. 発熱源を NekRS に転送し、流体温度と速度を計算。
     3. NekRS の温度結果を FM モデルにフィードバックし、行列補間を更新。
     4. 収束するまで（温度が定常状態に達するまで）反復。

3. 主要な貢献 (Key Contributions)

• FM 法の MOOSE 実装: 従来のモンテカルロ法に代わり、計算効率の高い分岐行列法を MOOSE/Cardinal フレームワークに統合し、MSFR の動的結合シミュレーションを可能にしました。
• 温度依存性の線形補間手法: 燃料温度分布に応じて FMDBs から分岐行列を線形補間する手法を実装し、リアルタイムな熱水力フィードバックを再現しました。
• GPU 加速と CPU 計算のハイブリッド利用: NekRS を GPU（Frontier スーパーコンピュータ）で実行し、FM 計算を CPU で実行する効率的な並列計算構成を確立しました。

4. 結果 (Results)

• 計算環境: Frontier スーパーコンピュータ上で実行。約 60 対流単位でシミュレーションし、その後の 80 対流単位で時間平均を行いました。
• 物理量:
  • 実効増倍率 $k_{eff}$: 1.04632。
  • 出口温度: 997.4 K。
  • 平均温度: 973.0 K。
• 分布特性:
  • 発熱源: 正弦波状の分布を示しました。
  • 流速と温度: 中央円筒の周辺部に「停滞領域（Stagnation regions）」が確認されました。ここでは流速が低く、燃料が加熱され続けるため、温度が最大値に達しています。
• 既存モデルとの比較:
  • 以前に開発された「NekRS-OpenMC」モデルとの比較を行いました。
  • 出口温度の差は 0.1 K、平均温度の差は 10.1 K でした。
  • 全エネルギーが同等であるため、この差異は許容範囲内と判断され、温度、流速、発熱源分布において良好な一致が確認されました。

5. 意義と結論 (Significance and Conclusions)

本研究は、NekRS と分岐行列法を結合した MSFR のマルチフィジックスモデルの成功した開発を示しています。

• 計算効率: 高コストなモンテカルロ法を FM 法に置き換えることで、熱水力フィードバックを含む動的シミュレーションの計算負荷を削減しつつ、物理的に妥当な結果を得ることができました。
• 現象の解明: 低流速領域が熱的停滞点となり、局所的な高温を引き起こすという MSFR 特有の熱水力・中性子工学の結合現象を明確に捉えました。
• 将来展望: 今後は、FM メッシュの解像度感度解析や、円筒メッシュ形状の影響評価を行う予定です。

このアプローチは、次世代原子炉の設計・安全性評価において、高精度かつ計算コストを抑えたマルチフィジックス解析手法として大きな可能性を秘めています。