原子炉シミュレーションPC｜Serpent+MCNP+OpenMC+中性子輸送

原子炉シミュレーションにおける計算科学の重要性

原子力エネルギーの安全性と効率性を決定づけるのは、炉心内部で発生する複雑な核反応の正確な予測です。原子炉内では、中性子がウラン235などの核分裂性物質に衝突し、連鎖反応を引き起こします。この過程で発生する中性子の挙動、エネルギー分布、核分裂生成物の分布を計算することを「中性子輸送計算」と呼びます。この計算には、膨大な数の粒子の軌跡を確率的に追跡する「モンテカルロ法」が用いられます。

モンテカルロ法は、統計的な手法を用いて物理現象をシミュレートする手法です。中性子が物質（燃料、冷却材、構造材）を通過する際に、どの程度の確率で吸収されるか、あるいは散乱されるかを、膨大な回数のサンプリングによって算出します。この計算は極めて高い精度を誇りますが、一方で計算負荷は天文学的です。中性子一つひとつの挙動を追跡するため、計算量は中性子の数や、扱う幾何学的形状の複雑さに比例して増大します。

202決算後の2026年現在、次世代の小型モジュール炉（SMR）や高速炉の開発が加速しており、これまでの大規模炉とは異なる複雑な冷却材（液体金属や溶融塩など）を用いたシミュレーションが求められています。これに伴い、従来のワークステーションでは対応不可能な、超高性能なHPC（ハイパフォーマンス・コンピューティング）環境としての「原子炉シミュレーションPC」の重要性がかつてないほど高まっています。

主要な中性子輸送計算ソフトウェアの詳細比較

原子炉シミュレーションの世界には、用途や目的に応じていくつかの主要なソフトウェアが存在します。これらは大きく分けて、決定論的手法（Deterministic method）とモンテカルロ法（Monte Carlo method）に分類されます。近年では、計算の正確性を担保するために、モンテカルロ法を用いた高精度な検証が標準となっています。

代表的なソフトウェアとして、業界標準である「MCNP 6.3」、格子計算に特化した「Serpent 2」、そしてオープンソースとして急速に普及している「OpenMC」が挙げられます。MCNPは、非常に複雑な形状の定義が可能であり、放射線防護や遮蔽計算にも広く利用されています。一方、Serpent 2は、核燃料の格子計算や燃焼計算（燃料が時間の経過とともにどのように変化するかを計算すること）に特化しており、非常に効率的な計算が可能です。OpenMCは、Pythonインターフェースを備え、現代的な並列計算環境やGPU活用に適した設計となっており、次世代の計算環境において主役となりつつあります。

以下の表に、これら主要ソフトウェアの機能的な違いをまとめます。

ソフトウェアの選択は、計算の「精度」と「計算時間」のトレードトレード（妥協点）をどのように管理するかによって決まります。例えば、初期の炉心設計段階では、OpenMCやSerpent 2を用いて、高速にパラメータを振った計算を行い、最終的な安全性検証としてMCNPによる厳密な計算を行うといった、多段階のワークフローが一般的です。

計算機構成の核心：AMD EPYC 9654と超多コア環境

原子炉シミュレーションPCの心臓部は、CPU（中央演算処理装置）のコア数と、それらが持つキャッシュメモリの容量に集約されます。モンテカルロ法による中性子輸送計算は、個々の粒子の追跡が独立して行えるため、「並列計算（複数の計算を同時に進めること）」との相性が極めて良いという特徴があります。つまり、コア数が増えれば増えるほど、計算時間は理論的に短縮されます。

ここで推奨されるのが、AMDのサーバー向けプロセッサ「EPYC 9654」です。このCPUは、1つのソケットに96個の物理コア、192個のスレッドを搭載しています。Zen 4アーキテクチャを採用しており、1コアあたりの計算能力が高いだけでなく、膨大なL3キャッシュ（CPU内部の高速メモリ）を搭載している点が重要です。中性子計算では、核データライブラリ（後述）から頻繁にデータを参照するため、キャッシュミス（必要なデータがキャッシュにない状態）が発生すると、計算速度が劇的に低下します。EPYC 9654の巨大なキャッシュ容量は、このボトルネックを解消する鍵となります。

また、AVX-512（Advanced Vector Extensions 512）などの命令セットへの対応も、SIMD（単一命令複数データ）演算を多用する物理計算において、スループットを向上させる重要な要素ですな。2026年現在の最新のシミュレーションコードは、これらの命令セットを最適に利用するように設計されており、コア数と命令セットの恩恵を最大限に受けることができます。

メモリ・ストレージ戦略：512GB RAMと高速NVMeの役割

原子炉シミュレーションにおけるメモリ（RAM）の役割は、単なる作業領域に留まりません。計算に必要となる「核データライブラリ（ENDF/B-VIIIなど）」や、複雑な幾何形状の定義、そして計算結果（タリー）の蓄積のために、膨大なメモリ容量が要求されます。特に、燃料の燃焼計算（Burnup calculation）を行う場合、時間ステップごとの組成変化を保持するために、メモリ消費量は指数関数的に増大します。

本構成では、512GBのDDR5 ECCメモリを採用します。ここで「ECC（Error Correcting Code）」が極めて重要です。ECCメモリは、メモリ内で発生したデータの誤りを自動的に検出し、修正する機能を持っています。数日間、あるいは数週間にわたって数千個のコアで計算を継続する場合、宇宙線などの影響によるビット反転（1が0に変わる現象）が致命的な計算エラーや、誤った物理結果を招く恐える可能性があります。そのため、科学計算用PCにおいてECCメモリは必須のコンポーネントです。

また、ストレージ（SSD）の性能も無視できません。計算の開始時に、数GBから数十GBに及ぶ核データライブラリをメモリへロードする際、ストレージの読み込み速度が全体の待ち時間に直結します。PCIe Gen5対応のNVMe SSDを使用することで、データ転送の遅延を最小限に抑えることが可能です。

GPUアクセラレーション：NVIDIA RTX 6000 Adaによる次世代計算

近年の計算科学における最大のパラダイムシフトは、GPU（Graphics Processing Unit）の活用です。従来のCPU中心の計算に対し、GPUを用いた「GPUアクセラレーション」は、中性子輸送計算の速度を数十倍に引き上げる可能性を秘めています。特にOpenMCのような、CUDA（NVIDIAの並列計算プラットフォーム）をネイティブにサポートするソフトウェアにおいて、その真価が発揮されます。

本構成で採用する「NVIDIA RTX 6000 Ada Generation」は、プロフェッショナル向けの最高峰GPUです。48GBという広大なビデオメモリ（VRAM）を搭載しており、これにより、CPUではメモリ不足となるような巨大な核データセットや、高解像度なメッシュ（空間を細かく分割した格子）を用いた計算も、GPU内部の高速なメモリ空間で完結させることが可能になります。

GPU計算のメリットは、数千個の小さなCUDAコアが、中性子粒子の移動計算を並列に処理することにあります。一方で、課題もあります。GPUへのデータ転送（PCIeバス経由）のオーバーヘッドや、複雑な分岐計算（中性子の衝突判定など）におけるスループットの低下です。そのため、計算内容に応じて「CPUで行うべき計算」と「GPUにオフロード（肩代わり）させるべき計算」を適切に使い分ける設計が求められます。

核データライブラリ（ENDF/B-VIII）の管理と計算負荷

中性子輸送計算の「精度」を決定づけるのは、ソフトウェアのアルゴリズムだけではありません。計算の基盤となる「核データライブラリ」の質が、物理的な正しさを左右します。代表的なものに、アメリカのLANL（ロスアラモス国立研究所）などが提供する「ENDF/B-データを挙げることができます。最新の「ENDF/B-VIII.0」や「ENDF/B-VIII.1」では、中性子のエネルギー（断面積：Cross-section）や、核反応の確率が、より高精度な実験値に基づいて更新されています。

「断面積（Cross-section）」とは、中性子が原子核に衝突して反応を起こす「的な大きさ」を指します。この値は、中性子のエネルギーによって劇的に変化します。シミュレーションでは、このエネルギー依存性を正確に扱うために、エネルギー軸を細かく分割して計算を行います。この分割数（Energy Grids）が増えるほど、精度は上がりますが、参照すべきデータ量が増大し、メモリ消費量と計算負荷が爆発的に増加します。

したがって、シミュレーションPCの設計においては、これらの膨大なデータセットを「いかに高速に、かつメモリ不足を起こさずに演算器へ供給するか」が最大の課題となります。ENDF/B-VIIIのような大規模なライブラリを扱う場合、単なる容量だけでなく、データの並列アクセス性能（メモリ帯域）が、計算の停止（Stall）を防ぐために極めて重要です。

炉心設計・燃料管理シミュレーションのワークフロー

原子炉シミュレーションPCを用いた実務的なワークフローは、単一の計算を実行して終わりではありません。それは、設計の「反復（イテレーション）」のプロセスです。

1. 初期設計（Initial Design）: 燃料の配置、冷却材の流量、制御棒の位置などのパラメータを決定します。ここでは、OpenMCなどの軽量で高速なコードを用い、広範なパラメータ空間を探索しますな。
2. 高精度検証（High-Fidelity Verification）: 設計案が固まった段階で、MCNP 6.3などの高精度コードを用い、安全性のマージン（余裕）を厳密に計算します。ここでは、複雑な幾何形状や、詳細な核データを用いた計算が行われます。
3. 燃焼計算（Burnup Calculation）: 燃料が核分裂を繰り返すにつれ、核分裂生成物（FP）が蓄積し、燃料の組成が変化していきます。この組成変化に伴う反応率（反応度）の低下を、時間ステップごとに計算します。
4. 燃料管理（Fuel Management）: 次回の燃料交換サイクルにおいて、どの燃料集合体を再利用し、どの位置に配置するかを最適化します。このプロセスには、膨大な数の計算パターンを評価する必要があるため、前述したEPYC 9654のような多コア環境が威力を発揮します。

この一連の流れをスムーズに行うためには、計算結果を解析するための「ポストプロセッシング（後処理）」環境も重要です。Pythonを用いたグラフ作成や、VTK（Visualization Toolkit）を用いた3次元的な中性子束（Neutron Flux）の可視化など、計算機と解析ソフトのシームレスな連携が、研究のスピードを左右しますな。

予算別・用途別シミュレーションPC構成案

原子炉シミュレーションPCの構築には、莫大なコストがかかります。研究予算やプロジェクトの規模に応じて、適切な構成を選択する必要があります。

「Entry」構成では、単一のGPUと、比較的コア数の少ないプロセッサを用いますが、基本的なモンテカルロ法の学習や、小規模な幾何形状の計算には十分です。一方、「High」構成は、まさに「計算機としての要塞」であり、企業や大規模研究施設での、実用的な炉心設計・安全性評価を目的とした、極めて高い信頼性とスループットを追求したものです。

よくある質問（FAQ）

Q1: モンテカルロ法と決定論的手法の計算時間の違いはどのくらいですか？ A1: 決定論的手法（Diffusion法など）は、メッシュを細かく分割するほど計算量が増えますが、粒子を追跡するモンテカルロ法に比べれば、圧倒的に高速です。しかし、モンテカルロ法は精度が極めて高く、複雑な幾何形状の扱いに優れています。用途に応じて、決定論的手法で大まかな分布を出し、モンテカルロ法で詳細を詰めるという使い分けが一般的です。

Q2: GPU計算（OpenMCなど）を行う際、CPUとのバランスはどう考えるべきですか？ A2: GPUは並列計算において強力ですが、データの転送（PCIeバス経由）がボトルネックになります。そのため、強力なGPUを搭載する場合でも、それにデータを供給するための高速なCPU（EPYC等）と、広帯域なメモリ（DDR5）を組み合わせることが、トータルの計算時間を短縮するために不可欠です。

Q3: 512GBものメモリが必要なのはどのような時ですか？ A3: 主に、非常に高精細な「エネルギーグリッド」や「空間メッシュ」を使用する場合、および、長期間の「燃焼計算（Burnup）」を行う場合です。核データの詳細度を上げると、メモリ消費量は指数関数的に増大するため、大規模な研究には大容量メモリが必須となります。

Q4: 冷却（クーリング）についての注意点はありますか？ A4: EPYC 9654のような高TDP（熱設計電力）のCPUと、RTX 6000 Adaのような強力なGPUを搭載する場合、発熱量は非常に膨大になります。一般的なデスクトップPCの冷却では不十分なことが多く、水冷システム、あるいはサーバーグレードの強力な空冷・液冷ソリューションが必要です。また、室内のエアコン容量も考慮しなければなりません。

Q5: ソフトウェアのライセンス費用はどの程度見込むべきですか？ A5: OpenMCのようなオープンソースは無料ですが、MCNP 6.3のような商用ソフトウェアは、利用機関の規模や研究内容によって、年間数十万〜数百万円単位のライセンス費用が発生することがあります。ハードウェアの予算だけでなく、ソフトウェアの維持費も予算計画に含める必要があります。

Q6: データのバックアップはどのように行うべきですか？ A6: シミュレーション結果（タリー）は、テラバイト級に達することがあります。RAID構成をとった大容量のHDD/SSD、あるいはネットワーク接続されたNAS（Network Attached Storage）への自動バックアップ体制を構築することが、研究の継続性を守る上で重要です。

まとめ

本記事では、原子炉シミュレーションに特化した、次世代のハイパフォーマンス・コンピューティング環境について解説しました。

• 計算手法の理解: 中性子輸送計算には、精度に優れたモンテカルロ法（MCNP, Serpent 2, OpenMC）が不可欠である。
• ハードウェアの核心: AMD EPYC 9654のような、多コア・大容量キャッシュを持つCPUが、計算効率を劇的に高める。
• メモリの重要性: 512GB以上の大容量かつECC（誤り訂正）機能を持つメモリが、大規模な核データと燃焼計算を支える。
• GPUの活用: NVIDIA RTX 6000 AdaのようなプロフェッショナルGPUの導入が、次世代の並列計算における鍵となる。
• データの精度: ENDF/B-VIIIなどの最新の核データライブラリを、いかに高速に処理するかが、シミュレーションの成否を分ける。
• ワークフローの最適化: 設計、検証、燃焼計算、燃料管理という、一連の反復プロセスを支えるインフラ構築が求められる。

原子炉シミュレーションPCの構築は、単なるPC自作の枠を超え、物理学、計算科学、そして高度なシステムエンジニアリングが融合した、極めて専門性の高い領域です。