求解中子动力学方程的加权蒙特卡罗方法

Weighted Monte Carlo solution of neutron kinetics equations

为了实现基于蒙特卡罗方法的中子动力学计算,在传统的直接蒙特卡罗动力学方法的基础上,提出了一种加权蒙特卡罗动力学方法。该方法通过引入粒子权重的概念,隐式考虑中子俘获反应和裂变反应过程中中子数目的变化,避免了模拟粒子的数目随时间的变化,降低了统计偏差,消除了程序计算过程中粒子的存库操作,提高了计算精度。基于单能点堆模型,开发了中子动力学计算程序NECP-Dandi,进行了大量数值验证与分析,包括无缓发中子、单组缓发中子、六组缓发中子、正阶跃反应性引入、负阶跃反应性引入、正脉冲反应性、负脉冲反应性和正线性反应性引入等情况。数值结果表明,相比于直接蒙特卡罗动力学方法,加权蒙特卡罗动力学方法在计算结果的精度和计算效率上有较为明显的改进,程序结构更为简洁。

The solution of time dependent neutronics equations still remains a challenging problem.A weighted Monte Carlo kinetics method（wMCk）is proposed based on traditional analog Monte Carlo kinetics method（aMCk）.The＂implicit capture＂is introduced to avoid the problem of low efficient tallies in aMCk;the definition of particle weighting leads to a more compact simulation flow due to the elimination of stack operation to particle bank.Using this method,a code named NECP-Dandi was developed in mono-energetic point-kinetics model for numerical verification and analysis.11 test cases with different reactivity insertions were employed to verify the method.Numerical results demonstrate that wMCk is superior to aMCk in terms of accuracy,efficiency and code structure.