原子炉中性子拡散シミュレーター — 1群拡散理論

原子炉の初期概念設計：新型炉の基本的なサイズや形状を決める第一歩として、1群計算でおおよその臨界質量や炉心半径を見積もります。複雑な多群計算を行う前に、パラメータ感度を調べるのに役立ちます。

教育・訓練用シミュレータ：原子炉物理の初学者が、炉心サイズ、材料特性（拡散長）、反射体の効果が臨界にどう影響するかを直感的に学ぶためのツールとして利用されます。机上の数式よりも理解が深まります。

制御棒価値の簡易評価：制御棒を「中性子吸収断面積が大きい領域」としてモデル化し、それを考慮した平均的な $\Sigma_a$ を設定することで、制御棒を挿入した時の $k_{\text{eff}}$ の低下を定性的に評価できます。

反射体設計の基礎検討：シミュレーターの「反射体厚さ T」のパラメータは、炉心を囲む反射体が中性子を炉内に跳ね返す効果を簡易モデルで表現しています。反射体を厚くすると中性子漏れが減り、必要な炉心サイズを小さくできることが確認できます。

まず、「1群」は万能ではないということを肝に銘じておこう。これはあくまで「平均的な中性子」を扱う超簡易モデルだ。例えば、高速中性子と熱中性子では拡散の仕方も吸収のされ方も全く違うのに、それを一つにまとめている。だから、定性的な傾向や感度分析には最適だが、実炉の詳細設計には使えない。次のステップとして多群拡散コードを使う必要があるんだ。

次に、パラメータの設定で陥りがちなのが「拡散長 L」と「移行面積 M²」の混同だ。シミュレーターでは両方スライダーがあるけど、L²は中性子が「減速を終えてから吸収されるまで」の平均移動距離の2乗。一方、M²は「核分裂で生まれてから吸収されるまで（減速過程も含む）」の距離の2乗だ。多くの場合、M² ≒ L² + τ（τはフェルミ年齢）という関係がある。例えば、軽水炉ではLは数cm、Mは数十cmオーダーになる。この違いを意識しないと、パラメータ変更の物理的意味を見誤るぞ。

最後に、「k∞を大きくすれば何でも臨界になる」という思い込み。確かにk∞は材料の性質を表すが、炉のサイズが小さすぎると漏れが支配的で、k∞をどれだけ上げてもkeffを1に持っていけない場合がある。逆に、ある一定以上のサイズ（臨界サイズ）を超えると、keffはk∞に漸近していく。この「漏れ」と「材料」の綱引きを体感できるのが、このツールの強みなんだ。