基于复杂接续面的蒙特卡罗接续计算方法研究

辐射防护 ›› 2023, Vol. 43 ›› Issue (1): 77-82.

基于复杂接续面的蒙特卡罗接续计算方法研究

余鸿, 吕焕文

核反应堆系统设计技术重点实验室,成都 621213

收稿日期:2022-05-06 出版日期:2023-01-20 发布日期:2023-03-20
通讯作者: 吕焕文。E-mail:lvhuanwen@126.com
作者简介:余鸿(1998—),男,2019年毕业于东南大学核工程与核技术专业,2022年毕业于中国核动力研究设计院辐射防护及环境保护专业,获硕士学位。E-mail:2628825629@qq.com

Monte Carlo continuation calculation method based on complex intermediate surfaces

YU Hong, LV Huanwen

Key Laboratory of Nuclear Reactor System Design Technology, Nuclear Power Institute of China,Chengdu 610213

Received:2022-05-06 Online:2023-01-20 Published:2023-03-20

摘要/Abstract

摘要： 为了提高大型复杂问题的计算效率,提出一种适用于具有多个空间区域表面与复杂曲面型表面结构体的蒙特卡罗接续计算方法。该方法以结构体外边界面作为接续面源进行接续计算,将写源过程的计算结果与接续过程的计算结果之和作为计算的最终结果;该方法通过一次写源计算与多次接续计算,避免直接计算过程中对相同结构的重复计算,大大提升了计算效率。以MCNP程序建立典型反应堆结构计算模型,以该接续计算方法进行4种结构方案的计算,并与直接计算方法的计算结果进行对比。结果表明,该接续计算方法的计算结果与直接计算的计算结果一致,且使用该方法比直接计算方法在效率上取得了约3倍的提升效果。因此,本研究提出的蒙特卡罗接续计算方法对于涉及大量重复计算的大型复杂问题可以取得显著的效率提升。

关键词: 蒙特卡罗方法, 接续计算, 效率, 多表面, 复杂曲面

Abstract: In order to improve computational efficiency of large and complex problems, a Monte Carlo Continuation Calculation method is proposed for structures with multiple spatial surfaces and complex curved surfaces. In this method, the external boundary of the structure is used as the source for the continuation calculation, and the sum of the calculation results of the source writing process and the continuation calculation process is taken as the final result. In terms of efficiency improvement, the method can avoid repeated calculation of the same structure in the process of direct calculation by means of one source writing process and several continuation calculation processes, which greatly improves the computational efficiency. A typical reactor structure calculation model is established by MCNP, and the continuation calculation results of four structural schemes are compared with those of direct calculation method. The results show that the calculation of the continuation method is consistent with those of the direct method, and the efficiency of the continuation calculation method is about three times higher than that of the direct calculation method. Therefore, the Monte Carlo continuation Calculation method proposed in this study can achieve significant efficiency improvement for large and complex problems involving a large number of repeated calculations.

Key words: Monte Carlo method, continuation calculation, efficiency, multiple surface, complex surfaces

中图分类号:

TL329

余鸿, 吕焕文. 基于复杂接续面的蒙特卡罗接续计算方法研究[J]. 辐射防护, 2023, 43(1): 77-82.

YU Hong, LV Huanwen. Monte Carlo continuation calculation method based on complex intermediate surfaces[J]. RADIATION PROTECTION, 2023, 43(1): 77-82.

参考文献

[1] 张磊,贾铭椿,龚军军,等.反应堆辐射屏蔽计算方法与程序概述[J].核电子学与探测技术,2018,38(04):516-520.
ZHANG Lei,JIA Mingchun, GONG Junjun, et al. Overview of calculation methods and codes for reactor radiation shielding[J]. Nuclear Electronics & Detection Technology, 2018,38(04):516-520.
[2] Lee D. Hybrid deterministic-monte carlo methods for neutron particle transport[C]// Korean Nuclear Society Spring Meeting. 2012.
[3] 庄思璇. 三维蒙特卡罗—离散纵标能量耦合中子输运计算方法研究[D].中国科学技术大学,2019.
[4] 丁大钊,叶春堂,赵志祥,等.中子物理学-原理、方法与应用[M].北京:原子能出版社,2001.
[5] Briesmeister J F. MCNP-A general monte carlo n particle transport code, version 4C[R]. Report LA-13709-M.Los Alamos National Laboratory, April 2000.
[6] 谢仲生,主编.核反应堆物理分析[M].西安:西安交通大学出版社,2004.
[7] Zheng Q, Shen W, Li Y Z, et al. A deterministic-stochastic energy-hybrid method for neutron-transport calculation[J]. Annals of Nuclear Energy, 2019,128: 292-299.
[8] 韩静茹. 三维蒙特卡罗-离散纵标双向耦合屏蔽计算方法研究[D].华北电力大学,2012.
[9] Zhao J, Li X, Wu B, et al. An automatic adaptive mesh generation method for weight window in Monte Carlo particle transport[J]. Annals of Nuclear Energy,2016,91:105-110.
[10] 杨德锋,程和平.MCNP程序分段-衔接计算方法研究[J].核动力工程,2010,31(S2):48-53.
YANG Defeng, CHENG Heping. Research on MCNP bootstrapping calculation method[J]. Nuclear Power Engineering, 2010,31(S2):48-53.
[11] 黄倩倩,唐邵华,杨寿海.MC-MC耦合计算方法在核电厂中子屏蔽设计中的应用[J].辐射防护,2019,39(04):309-313.
HUANG Qianqian, TANG Shaohua, YANG Shouhai. Application of MC-MC coupled method in neutron shielding analysis of NPP[J]. Radiation Protection, 2019,39(04):309-313.
[12] Wijk A J van, Hoogenboom J E. A priori efficiency calculations for monte carlo applications in neutron transport[J]. Progress in Unclear Science and Technology,2011,2:722-732.DOI:10.15667/pnst.2.722.
[13] Qiu Y S, Aufiero M, Wang K, et al. Development of sensitivity analysis capabilities of generalized responses to nuclear data in Monte Carlo code RMC[J]. Annals of Nuclear Energy, 2016, 97:142-152.