Global variance reduction based on forward Monte Carlo calculation

Shi Tao; Ma Jimin; Qiu Youheng; Huang Hongwen; Li Zhenghong; Qian Dazhi

doi:10.11884/HPLPB201830.170163

Volume 30 Issue 1

Jan. 2018

Turn off MathJax

Article Contents

Article Navigation > High Power Laser and Particle Beams > 2018 > 30(1): 016006

Shi Tao, Ma Jimin, Qiu Youheng, et al. Global variance reduction based on forward Monte Carlo calculation[J]. High Power Laser and Particle Beams, 2018, 30: 016006. doi: 10.11884/HPLPB201830.170163

Citation:

Shi Tao, Ma Jimin, Qiu Youheng, et al. Global variance reduction based on forward Monte Carlo calculation[J]. High Power Laser and Particle Beams, 2018, 30: 016006. doi: 10.11884/HPLPB201830.170163

Shi Tao, Ma Jimin, Qiu Youheng, et al. Global variance reduction based on forward Monte Carlo calculation[J]. High Power Laser and Particle Beams, 2018, 30: 016006. doi: 10.11884/HPLPB201830.170163

Citation:

Shi Tao, Ma Jimin, Qiu Youheng, et al. Global variance reduction based on forward Monte Carlo calculation[J]. High Power Laser and Particle Beams, 2018, 30: 016006. doi: 10.11884/HPLPB201830.170163

PDF( 4258 KB)

Global variance reduction based on forward Monte Carlo calculation

doi: 10.11884/HPLPB201830.170163

1.
Institute of Nuclear Physics and Chemistry, CAEP, Mianyang 621999, China
2.
Institute of Applied Physics and Computational Mathematics, Beijing 100094, China

Received Date: 2017-05-08
Rev Recd Date: 2017-09-14
Publish Date: 2018-01-15

Abstract

Abstract

Monte Carlo method is extensively employed in radiation shielding calculation with the advantages of high fidelity geometry modeling, complex radiation source description and continuous-energy cross sections. However, the simulation is impractical for Monte Carlo particles that have little chance of being transported to the far-source area and the statistical error may be unacceptable. Accordingly, effective global variance reduction method (GVR) is significant for Monte Carlo deep-penetration radiation shielding calculation. This work constructs weight window parameters based on forward Monte Carlo calculation to bias the particles position and energy during transport. According to the validation of shielding benchmark, particles have been guided to far-source area by using the weight window. The relative error has been decreased to less than 10%, and almost all the mesh cells have been tallied. The results show that the transport parameters obtained from forward Monte Carlo calculation could realize global variance reduction.
- Monte Carlo,
- forward transport,
- shielding calculation,
- global variance reduction

FullText(HTML)

References(15)

References

[1]	Cooper M A, Larsen E W. Automated weight windows for global Monte Carlo particle transport calculations[J]. Nuclear Science and Engineering, 2001, 137: 1-13. doi: 10.13182/NSE00-34
[2]	Becker T L, Larsen E W. The application of weight windows to global Monte Carlo problems[C]//Proc of International Conference on Mathematics, Computational Methods & Reactor Physics. 2009.
[3]	Becker T L. Hybrid Monte Carlo/deterministic methods for deep-penetration problems[D]. Michigan: University of Michigan, 2009.
[4]	Peplow D E, Blakeman E D, Wagner J C. Advanced variance reduction strategies for optimizing mesh tallies in MAVRIC[J]. Transactions of the ANS, 2007, 97: 595-597.
[5]	Wagner J C, Peplow D E, Mosher S W. FW-CADIS method for global and semi-global variance reduction of Monte Carlo radiation transport calculations[J]. Nuclear Science and Engineering, 2014, 177: 37-57.
[6]	Wijk Van A J, den Eynde Van G, Hoogenboom J E. An easy to implement global variance reduction procedure for MCNP[J]. Annals of Nuclear Energy, 2011, 38: 2496-2503. doi: 10.1016/j.anucene.2011.07.037
[7]	Fan X, Zhang G, Wang K. Development of new variance reduction methods based on weight window technique in RMC code[J]. Progress in Nuclear Energy, 2016, 90: 197-203. doi: 10.1016/j.pnucene.2016.03.023
[8]	邱有恒, 邓力, 李百文, 等. 基于MESH技术的伴随权重窗自动生成方法及其应用[J]. 原子能科学技术, 2014, 48(11): 2148-2152. https://www.cnki.com.cn/Article/CJFDTOTAL-YZJS201411039.htm Qiu Youheng, Deng Li, Li Baiwen, et al. Automatic generation method mesh-based adjoint weight window and its application. Atomic Energy Science and Technology, 2014, 48(11): 2148-2152 https://www.cnki.com.cn/Article/CJFDTOTAL-YZJS201411039.htm
[9]	郑征, 吴宏春, 曹良志, 等. 蒙特卡罗与离散纵标耦合方法在压水堆堆腔漏束计算中的应用[J]. 强激光与粒子束, 2013, 24(12): 2946-2950. doi: 10.3788/HPLPB20122412.2946 Zheng Zheng, Wu Hongchun, Cao Liangzhi, et al. Application of Monte Carlo and discrete ordinate coupling method to pressurized water reactor cavity radiation streaming calculation. High Power Laser and Particle Beams, 2013, 24(12): 2946-2950 doi: 10.3788/HPLPB20122412.2946
[10]	韩静茹, 陈义学, 袁龙军, 等. 基于蒙特卡罗-离散纵标双向耦合方法的堆坑粒子注量率计算[J]. 原子能科学技术, 2014, 48(9): 1621-1626. https://www.cnki.com.cn/Article/CJFDTOTAL-YZJS201409016.htm Han Jingru, Chen Yixue, Yuan Longjun, et al. Particle fluence rate calculation of reactor pit based on MC-S_N bidirectional coupled method. Atomic Energy Science and Technology, 2014, 48(9): 1621-1626 https://www.cnki.com.cn/Article/CJFDTOTAL-YZJS201409016.htm
[11]	肖锋, 应栋川, 章春伟, 等. 离散纵标与蒙特卡罗耦合方法在反应堆屏蔽计算中的应用[J]. 核动力工程, 2014, 35(5): 9-12. https://www.cnki.com.cn/Article/CJFDTOTAL-HDLG201405003.htm Xiao Feng, Ying Dongchuan, Zhang Chunwei, et al. Application of discrete ordinate and Monte-Carlo coupling method for reactor shielding calculation. Nuclear Power Engineering, 2014, 35(5): 9-12 https://www.cnki.com.cn/Article/CJFDTOTAL-HDLG201405003.htm
[12]	X-5 Monte Carlo Team. MCNP—a general Monte Carlo N-particle transport code, version 5[R]. LA-UR-03-1987.
[13]	邱有恒, 邓力, 上官丹骅, 等. JMCT软件的自动源偏倚抽样功能及其应用[J]. 计算物理, 2016, 33(5): 593-598. https://www.cnki.com.cn/Article/CJFDTOTAL-JSWL201605012.htm Qiu Youheng, Deng Li, Shangguan Danhua, et al. Auto source particle biasing sampling function of JMCT and its application. Chinese Jounal of Computational Physis, 2016, 33(5): 593-598 https://www.cnki.com.cn/Article/CJFDTOTAL-JSWL201605012.htm
[14]	Peplow D E. Comparison of hybrid methods for global variance reduction in shielding calculations[J]. Transactions of the American Nuclear Society, 2012, 107: 512-515.
[15]	OECD-NEA. Ispra iron benchmark experiment (EURACOS)[EB/OL]. https://www.oecdnea.org/science/wprs/shielding/sinbad/eurac_fe/eufe-abs.htm.