🙋
「1群拡散理論」って何ですか？「群」って何の群なんですか？
🎓
大まかに言うと、原子炉内を飛び回る中性子を、全部まとめて「一つの平均的なエネルギーを持つ集団」として扱う超シンプルなモデルだよ。実際の中性子は高速から熱中性子までエネルギーがバラバラだけど、それを「一群」にまとめて計算を簡単にしているんだ。このシミュレーターでは、上のスライダーで「拡散長 L」や「移行面積 M²」を変えると、この平均的な中性子の動きやすさが変わるんだ。
🙋
え、そうなんですか！で、画面に出てくる keff が1を超えたり下回ったりするけど、これが臨界ということ？
🎓
その通り！実務ではこの keff（実効増倍率）が設計の生命線だ。keff=1が丁度臨界で、出力が一定に保たれる状態。試しに「炉心半径 a」のスライダーを小さくしてみて。炉が小さすぎると中性子が外に漏れちゃって、keff が1を下回る（未臨界）になるのがわかるよ。逆に「無限増倍率 k∞」を大きくすると、keff>1の超臨界になって中性子束が増えていく。
🙋
なるほど！グラフで山形になる中性子束の形は、どうしてああなるんですか？真ん中で一番高くて、端っこでゼロに。
🎓
良いところに気が付いたね。これは拡散方程式の解で、球炉心だと $\Phi(r) \propto \sin(B_g r)/(B_g r)$ って形になるんだ。中心で核分裂が盛んに起きて中性子が生まれ、端に行くに従って炉外に漏れるから、自然とあの釣り鐘型になる。パラメータを変えると山の広がり方（$B_g$）が変わるから、グラフの形も変わってくるんだ。操作して確かめてみて！
k∞（核分裂生成断面積）を増やすとkeffは上昇し、炉心半径を大きくすると漏れが減りkeffが上がります。逆に反射体厚さを増やすと反射体節約効果でkeffが上昇します。まずk∞または炉心半径を微調整してください。
反射体がない場合と比べ、反射体を設置することで臨界に必要な炉心半径がどれだけ小さくなるかを示す指標です。反射体が中性子を炉心に戻すため、同じkeffを得るのに炉心を小さくできます。シミュレーターで反射体厚さを変えて比較すると直感的に理解できます。
炉心中心は周囲より中性子の漏れが少なく、核分裂で生成された中性子が蓄積しやすいためです。1群拡散理論では、球形炉心の場合、中性子束はcos関数またはBessel関数に従い、中心が最大になります。反射体があると炉心境界付近の勾配が緩やかになります。
拡散長Lは中性子が吸収されるまでに拡散する平均距離の目安です。Lが大きいほど中性子が遠くまで漏れやすく、臨界にはより大きな炉心または高いk∞が必要です。逆にLが小さいと中性子は吸収されやすく、炉心サイズを小さくできます。
原子炉の初期概念設計：新型炉の基本的なサイズや形状を決める第一歩として、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∞に漸近していく。この「漏れ」と「材料」の綱引きを体感できるのが、このツールの強みなんだ。