蒙特卡罗全局权窗和均匀裂变源方法在全堆临界计算中的应用研究

张显; 刘仕倡; 强胜龙; 张文鑫; 尹强; 崔显涛; 陈义学

doi:10.7538/yzk.2020.youxian.0760

蒙特卡罗全局权窗和均匀裂变源方法在全堆临界计算中的应用研究

华北电力大学核科学与工程学院,北京102206
中国核动力研究设计院核反应堆系统设计技术国家重点实验室,四川成都610213

计量
- 文章访问数: 117
- HTML全文浏览量: 0
- PDF下载量: 825
出版历程
- 刊出日期: 2021-06-19

Application of Monte Carlo Global Weight Window and Uniform Fission Site Methods in Full Core Criticality Calculation

School of Nuclear Science and Engineering, North China Electric Power University, Beijing 102206, China
Science and Technology on Reactor System Design Technology Laboratory, Nuclear Power Institute of China, Chengdu 610213, China

摘要

摘要: 蒙特卡罗方法在大型裂变反应堆模拟时，由于不同区域的功率密度有较大差异，与高功率区域相比，低功率区域的粒子数相对不足，统计误差也较大。针对这类具有全局特征的蒙特卡罗输运计算问题，为获得中子通量在全空间的准确统计，本文针对全局权窗（global weight window, GWW）方法和均匀裂变源（uniform fission site, UFS）方法两类减方差算法展开了研究，在蒙特卡罗程序RMC中开发了基于通量的GWW功能，将GWW算法以及两类UFS算法应用于Hoogenboom-Martin压水堆全堆基准题，实现了较好的全局减方差效果，降低了堆芯计算的不对称程度，计算效率有较大提高。为结合两种减方差算法的优势，本文在蒙特卡罗程序RMC中开发了GWW与UFS的组合算法，并应用于压水堆全堆临界计算。结果表明，该方法在相同的计算条件下，堆芯外围的方差进一步减小，计算效率得到进一步提升，验证了该方法在全堆临界计算中的有效性。
- 蒙特卡罗 ,
- 全局减方差 ,
- 权窗 ,
- 均匀裂变源 ,
- 全堆
Abstract: When the Monte Carlo method is used to simulate a large-scale fission reactor, the power density of different regions varies greatly. Compared with the high-power region, the number of particles in the low-power region is relatively insufficient, and the statistical error is also larger. For this kind of Monte Carlo transport calculation problem with global characteristics, in order to obtain accurate statistics of neutron flux in the whole space, the global weight window (GWW) method and the uniform fission site (UFS) method were studied in this paper, and the flux-based GWW function was developed in RMC. The GWW function and two UFS functions were applied to the Hoogenboom-Martin PWR whole core benchmark, which achieves a better global variance reduction effect, reduces the asymmetry of the reactor core calculation and greatly improves the calculation efficiency. In order to combine the advantages of the two kinds of variance reduction algorithm, the combined algorithm of GWW and UFS was developed in the Monte Carlo code RMC and applied to the critical calculation of the PWR. The results show that under the same calculation conditions, with adoption of the combined algorithm, the variance of the periphery of reactor core is further reduced and the calculation efficiency is greatly improved, which verifies the effectiveness of the method in the critical calculation of the whole core.
- Monte Carlo ,
- global variance reduction ,
- weight window ,
- uniform fission site ,
- full core

HTML全文

参考文献(11)

[1]	CARTER L L, CASHWELL E D. Particle-transport simulation with the Monte Carlo method[R]. US: USERDA Technical Information Center, 1975.
[2]	BROWN F B. A review of Monte Carlo criticality calculations particle-transport simulation with the Monte Carlo method convergence, bias, statistics[C]∥Proc Int Conf Mathematics, Computational Methods & Reactor Physics. US: [s. n.], 2009.
[3]	IVANOV A, SANCHEZ V. Variance reduction in high resolution coupled Monte Carlo thermal-hydraulics calculations[C]∥ANS MC2015-Joint International Conference on Mathematics and Computation (M&C), Supercomputing in Nuclear Applications (SNA) and the Monte Carlo (MC) Method. US: [s. n.], 2015.
[4]	SMITH K S. Monte Carlo for practical LWR analysis: What’s needed to get to the goal[J]. Transactions of the American Nuclear Society, 2011, 104: 313-318.
[5]	DAVIS A, TURNER A. Comparison of global variance reduction techniques for Monte Carlo radiation transport simulations of ITER[J]. Fusion Eng Des, 2011, 86(9): 2698-2700.
[6]	JESSICA L H, THOMAS M S. A method for reducing the largest relative errors in Monte Carlo iterated fission-source calculations[C]∥M&C 2013. US: [s. n.], 2013.
[7]	HOOGENBOOM J E, MARTIN W R, BOJAN P. The Monte Carlo benchmark test: Aims, specifications and first results[C]∥Proceedings of the International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering. [S. l.]: [s. n.], 2011.
[8]	CHADWICK M B, OBLOINSKY P, HODGSON P E, et al. Pauli-blocking in the quasideuteron model of photoabsorption[J]. Physical Review C, 1991, 44(2): 814-819.
[9]	DANIEL J K, BRIAN N A, BRYAN R H. MC21 analysis of the MIT PWR benchmark: hot zero power result[C]∥M&C 2013. US: [s. n.], 2013.
[10]	HUNTER J L, SUTTON T M. A method for reducing the largest relative errors in Monte Carlo iterated-fission-source calculations[C]∥M&C 2013. US: [s. n.], 2013.
[11]	吴屈,范潇,王刚,等. RMC均匀裂变源方法及其在堆用蒙特卡罗计算不对称性问题中的应用[J]. 原子能科学技术,2017,51(7):1232-1238. WU Qu, FAN Xiao, WANG Gang, et al. RMC uniform-fission-site method and its application to asymmetry of reactor core Monte-Carlo computation[J]. Atomic Energy Science and Technology, 2017, 51(7): 1232-1238(in Chinese).