用于厚屏蔽小探测器的蒙特卡罗模拟减方差方法研究

发布时间：2021-04-02 14:18

　　反应堆屏蔽计算中经常出现厚屏蔽、小探测器问题,常规蒙特卡罗方法难以有效解决。基于自动重要抽样（AIS）方法,本文提出了小探测器自动重要抽样（SDAIS）方法,并针对小探测器问题,优化了AIS方法的虚粒子赌分裂算法。该方法在自主开发的蒙特卡罗屏蔽程序MCShield上进行了实现。使用NUREG/CR-6115 PWR基准题验证该方法的正确性和计算效率。结果表明,SDAIS方法可有效地解决厚屏蔽小探测器问题,相比AIS方法及传统的重要性方法,计算效率提升1～2个量级。 

【文章来源】：原子能科学技术. 2020,54(07)北大核心EICSCD

【文章页数】：7 页

【部分图文】：

AIS方法流程图

DXTRAN球上的虚粒子通过抽样产生,并被赋予相应的权重,以保证无偏性。图3为DXTRAN球上的虚粒子抽样方法。设粒子在P1发生碰撞,DXTRAN球心为P0,DXTRAN外球半径为RO,内球半径为RI,P1P0长度为L。DXTRAN外球、内球在P1P0方向上极角最大值分别为θO、θI,其余弦值分别为ηO、ηI。虚粒子位置PS在DXTRAN外球上,具体坐标通过抽样以P1P0方向为基准的极角和方位角确定。极角η余弦值按式(2)的概率密度函数p(η)抽样。图3 DXTRAN球上虚粒子的抽样方法

DXTRAN球上虚粒子的抽样方法

【参考文献】：
期刊论文
[1]Improved algorithms and coupled neutron-photon transport for auto-importance sampling method[J]. 王鑫,李君利,武祯,邱睿,李春艳,梁漫春,张辉,刚直,徐红.  Chinese Physics C. 2017(01)
[2]粒子输运蒙特卡罗模拟现状概述[J]. 邓力,李刚.  计算物理. 2010(06)
[3]内照射小器官剂量计算中减方差技巧的比较和应用[J]. 张崭,李君利,武祯,范佳锦,曾志.  清华大学学报(自然科学版). 2007(S1)
[4]一种改进的计算探测器校正因子的相关抽样方法[J]. 武祯,李君利,程建平.  高能物理与核物理. 2006(08)



本文编号：3115421