王学栋,郭丁情,佟立丽,曹学武

• 1:上海交通大学机械与动力工程学院
• 2:中广核工程有限公司

摘要(Abstract)：

本文基于核电厂严重事故下辐射环境的蒙特卡罗模拟需求,针对模拟计算中经常存在的深穿透区域计数结果很难或者无法收敛的问题,在相关理论基础上,结合蒙特卡罗正算输运的特点建立了迭代计算产生偏倚参数的具体方案,形成了一套完整的全局减方差方法的应用流程。为了验证方法的应用效果,本文针对核电厂模型开展了相关迭代计算,并与直接模拟的情形进行了对比,对比结果表明,采用通量权窗指导粒子输运,问题几何的全局范围内都实现了减方差的效果,模拟质量和效率(品质因数)都获得了大幅度提升,初步验证了方法在解决深穿透问题方面的有效性。

关键词(KeyWords)： 蒙特卡罗;减方差;深穿透;迭代

Abstract:

Keywords:

基金项目(Foundation): 国家科技重大专项课题(2019ZX06004013,2017ZX06002003-001-002)

作者(Author): 王学栋,郭丁情,佟立丽,曹学武

参考文献(References)：

• [1]国家核安全局.HAF 102-2016核动力厂设计安全规定[S].北京:中华人民共和国生态环境部,2009.
• [2]Nuclear Energy Institute.Diverse and flexible coping strategies(FLEX)implementation guide[R].NEI 12-06,Washington DC,2015.
• [3]Lux I,Koblinger L.Monte Carlo particle transport methods:neutron and photon calculations[M].Florida:CRC Press,1990.
• [4]Briesmeister F.MCNPTM,A General Monte Carlo N-particle transport code,version 4B-Manual[R].New Mexico:Los Alamos National Laboratory,1997.
• [5]Wagner C,Peplow E D,Scott W.FW-CADIS method for global and regional variance reduction of Monte Carlo radiation transport calculations[J].Nuclear Science and Engineering,2014,176:37-57.
• [6]Wagner C,Haghighat A.Automated variance reduction of Monte Carlo shielding calculations using the discrete ordinates adjoint function[J].Nuclear Science and Engineering,1998,128:186-208.
• [7]Tang S,Stevens P N,Hoffman T J.Methods of Monte Carlo biasing using two-dimensional discrete ordinates adjoint Flux[R].New Mexico:Los Alamos National Laboratory,1977.
• [8]Miller C,Wright G A,Boyle G B,et al.The use of an inbuilt importance generator for acceleration of the Monte Carlo code MCBEND[C].Paris:International Conference on the Physics of Reactors,1990.
• [9]Kenneth A,Riper V,Urbatsch J,Soran P D,et al.AVATAR-Automatic variance reduction in Monte Carlo calculations[R].New Mexico:Los Alamos National Laboratory,1997.
• [10]Vanijk J,Eynde V G,Hoogenboom J E.An easy to implement global variance reduction procedure for MCNP[J].Annals of Nuclear Energy,2011,38:2496-2503.
• [11]Brewmeister F.MCNP-A general Monte Carlo N-particle transport code[R].Version A,LA-12625,Los Alamos National Laboratory,1993.
• [12]US NRC.Accident source terms for light-water nuclear power plants[R].NUREG-1465,1995.