强激光与粒子束, 2018, 30 (1): 016009, 网络出版: 2018-01-30  

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

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

李捷, 李云召, 吴宏春, 郑琪. 求解中子动力学方程的加权蒙特卡罗方法[J]. 强激光与粒子束, 2018, 30(1): 016009. Li Jie, Li Yunzhao, Wu Hongchun, Zheng Qi. Weighted Monte Carlo solution of neutron kinetics equations[J]. High Power Laser and Particle Beams, 2018, 30(1): 016009.

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!