🧑🎓
「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$）が変わるから、グラフの形も変わってくるんだ。触って確かめてみて！
1群拡散理論の基礎となる、定常状態の中性子の収支を表す拡散方程式です。生成（核分裂）、吸収、拡散（漏れ）のバランスを記述します。
$$D\nabla^2\Phi - \Sigma_a\Phi + \frac{\nu\Sigma_f}{k_{\text{eff}}}\Phi = 0$$
$\Phi$: 中性子束 (neutron flux)、$D$: 拡散係数、$\Sigma_a$: 吸収断面積、$\nu\Sigma_f$: 核分裂生成断面積、$k_{\text{eff}}$: 実効増倍率。$\nabla^2\Phi$は中性子の拡散（漏れ）を表します。
球形炉心の幾何学的バックリングと、実効増倍率 $k_{\text{eff}}$ を直接計算する実用的な公式です。シミュレーターの核心となる式です。
$$B_g^2 = \left(\frac{\pi}{a}\right)^2, \quad k_{\text{eff}}= \frac{k_\infty}{1 + M^2 B_g^2}$$
$B_g^2$: 幾何学的バックリング（炉の形状とサイズによる中性子漏れのしやすさ）、$a$: 炉心半径、$k_\infty$: 無限媒体増倍率（炉が無限に大きいと仮定した時の増倍率）、$M^2$: 移行面積（中性子が発生してから吸収されるまでの平均移動距離の2乗）。$M^2 B_g^2$ が漏れの影響を表し、分母が1より大きくなるほど $k_{\text{eff}}$ は小さくなります。
よくある誤解と注意点
まず、「1群」は万能ではないということを肝に銘じておこう。これはあくまで「平均的な中性子」を扱う超簡易モデルだ。例えば、高速中性子と熱中性子では拡散の仕方も吸収のされ方も全く違うのに、それを一つにまとめている。だから、定性的な傾向や感度分析には最適だが、実炉の詳細設計には使えない。次のステップとして多群拡散コードを使う必要があるんだ。
次に、パラメータの設定で陥りがちなのが「拡散長 L」と「移行面積 M²」の混同だ。シミュレーターでは両方スライダーがあるけど、L²は中性子が「減速を終えてから吸収されるまで」の平均移動距離の2乗。一方、M²は「核分裂で生まれてから吸収されるまで（減速過程も含む）」の距離の2乗だ。多くの場合、M² ≒ L² + τ（τはフェルミ年齢）という関係がある。例えば、軽水炉ではLは数cm、Mは数十cmオーダーになる。この違いを意識しないと、パラメータ変更の物理的意味を見誤るぞ。
最後に、「k∞を大きくすれば何でも臨界になる」という思い込み。確かにk∞は材料の性質を表すが、炉のサイズが小さすぎると漏れが支配的で、k∞をどれだけ上げてもkeffを1に持っていけない場合がある。逆に、ある一定以上のサイズ（臨界サイズ）を超えると、keffはk∞に漸近していく。この「漏れ」と「材料」の綱引きを体感できるのが、このツールの強みなんだ。
関連する工学分野
このシミュレーターの根幹を成す「拡散方程式」は、実は原子炉物理以外の多くの輸送現象でも顔を出す。例えば、半導体プロセスにおける不純物の拡散だ。シリコンウェハーにドーパントを注入した後、熱処理で拡散させる過程は、まさに拡散方程式で記述される。生成項が無い（$\nabla^2 C = (1/D) \partial C/\partial t$）など違いはあれ、数学的には兄弟のようなものだ。
また、中性子束分布の釣り鐘型プロファイルを求める計算は、構造力学の固有値解析と数学的に酷似している。拡散方程式を変形すると、$\nabla^2\Phi + B_g^2\Phi = 0$というヘルムホルツ方程式になるが、これは固定端の膜の振動モード（固有振動数と振動形）を求める式と同一形だ。炉心の「臨界」が方程式の「固有値」、中性子束分布が「固有関数」に対応する。CAEエンジニアなら、FEMソルバーが解く問題の本質がここで繋がるはず。
さらに応用を広げると、光子拡散理論（生体光学イメージング）や地下水流動・汚染物質移行のモデリングでも同様の式が使われる。つまり、このツールで学ぶのは、「発生・拡散・消失のバランス」という普遍的な物理概念の第一歩なんだ。
発展的な学習のために
まず次の一歩は、「2群拡散理論」を学ぶことだ。高速群と熱群に分けて、その間に「減速」というプロセスを明示的に入れる。これだけで、現実の炉心で重要な「スペクトル効果」や「減速材の役割」が理解できるようになる。NovaSolverで遊んだ後は、「高速群と熱群の拡散長はどっちが大きい？」「なぜ軽水炉では燃料を格子状に配置するのか？」といった問いに挑戦してみよう。
数学的な背景を深めたいなら、「変数分離法」と「境界条件の扱い」をしっかり復習してほしい。このシミュレーターが解いている球幾何の方程式は、変数分離で $\Phi(r) = R(r)$ と置き、$R(r)$ に関する常微分方程式を解く。その解が $\sin(B_g r)/(B_g r)$ という形になり、境界条件 $\Phi(a)=0$ から $B_g = \pi/a$ が求まる。この流れを自分で数式を追いながら理解できれば、円柱炉心など他の幾何への応用も見通せるようになる。
最終的には、実務で使われるコード（例えば、MONJU、SRAC、MVPなど）の入力デッキを読んでみることをお勧めする。そこには、数十～数百ものエネルギー群、詳細な燃料組成、複雑な炉心配置が定義されている。NovaSolverの「一群」「均質炉心」という仮定が、いかに現実を単純化したものかが痛感できると同時に、その「単純化の本質」を学んだからこそ、複雑な入力データの一つ一つの意味が理解できるようになるんだ。