Fission source convergence of Monte Carlo criticality calculations in weakly coupling fissile arrays

弱く結合した核分裂配列系でのモンテカルロ法臨界計算の核分裂中性子源の収束性

山本 俊弘; 中村 剛実*  ; 三好 慶典

Yamamoto, Toshihiro; Nakamura, Takemi*; Miyoshi, Yoshinori

弱く結合した核燃料ユニット配列でモンテカルロ法により臨界計算を行うと核分裂中性子源の収束が悪くなることを示した。複数ユニットの配列系に結合係数を導入し、収束が、とりわけ対称的な配列系に対して悪いのは、正しい分布に戻そうとする復元力が弱いことと、その統計的変動が原因であることを定量的に示した。配列系における核分裂分布比を核分裂確率行列を用いて求める新たな方法を示した。この方法を用いれば単に比だけではなく、その統計誤差も与えられ、得られた結果の精度評価が可能となる。核分裂分布比がある値に強制的に収束するように領域別にウェイトを調整する方法を開発した。これを用いれば、非常に強い収束を得ることができる。臨界計算で得られた実効増倍率と核分裂分布比との相関を定式化した。二つのユニットからなる対称な系には両者には有意な相関がないことを理論的及び経験的に示した。一般的には、わずかばかり非対称になった系においては、実効増倍率は不正確な核分裂源の分布によってバイアスを持つことがわかった。弱く結合した配列系でモンテカルロ法により臨界計算を行う場合には、一世代あたりの中性子数を多くして計算したほうが望ましい。

no abstracts in English