非定常输运模拟中基于粒子标识分类的源偏倚算法

计算物理 ›› 2016, Vol. 33 ›› Issue (6): 639-644.

非定常输运模拟中基于粒子标识分类的源偏倚算法

上官丹骅, 许海燕

北京应用物理与计算数学研究所, 北京 100094

收稿日期:2015-09-07 修回日期:2016-03-03 出版日期:2016-11-25 发布日期:2016-11-25
作者简介:上官丹骅(1977-),男,博士,副研究员,从事蒙特卡罗算法及应用研究,E-mail:sgdh@iapcm.ac.cn
基金资助:
中国工程物理研究院科学基金（2014B0202029）资助项目

Particle-Flag Based Source Bias Algorithm for Simulating Time-Dependent Particle Transport

SHANGGUAN Danhua, XU Haiyan

Institute of Applied Physics and Computational Mathematics, Beijing 100094, China

Received:2015-09-07 Revised:2016-03-03 Online:2016-11-25 Published:2016-11-25

摘要/Abstract

摘要： 在非定常输运问题的多步蒙特卡罗模拟中，根据粒子的不同属性进行标识分类可以得到非常细致的系统相关标识物理量.对于某些目标标识物理量，模拟的样本中造成有效贡献的样本相对较少，导致这些物理量模拟结果的涨落较大，而单靠增加总样本数不能高效地使有效样本达到一个合理的水平.本文提出一种基于标识分类的源偏倚算法，将增加的所有样本定向赋予目标类粒子，从而高效地降低目标标识物理量的统计误差且不影响非目标标识物理量的计算.以一维多层介质非定常输运模型验证上述结论.

关键词: 非定常输运问题, 蒙特卡罗方法, 标识, 源偏倚

Abstract: In multi-step Monte Carlo simulation of time-dependent particle transport problems, particle-flag based physical quantity can be calculated by appropriate classification of diverse particle's attributes. Some particle-flag based physical quantities' fluctuation are strong since only very small fraction of total histories can make non-zero contribution and it is inefficient to deal with this problem by increasing purely total history number. A source bias algorithm is proposed to decrease stochastic error of target quantity by increasing number of source particle with a specific type only. Meanwhile, precision of non-target quantities are hardly decreased. A one-dimensional multi-layer model is utilized to display effect of the method.

Key words: time-dependent particle transport, Monte Carlo method, particle-flag, source bias

中图分类号:

O571.51

上官丹骅, 许海燕. 非定常输运模拟中基于粒子标识分类的源偏倚算法[J]. 计算物理, 2016, 33(6): 639-644.

SHANGGUAN Danhua, XU Haiyan. Particle-Flag Based Source Bias Algorithm for Simulating Time-Dependent Particle Transport[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2016, 33(6): 639-644.

[1]	邓楚楚, 张鑫源, 陈水源, 吴宇晗, 叶晴莹, 黄志高. 小缺陷铁纳米环的磁特性[J]. 计算物理, 2020, 37(4): 473-478.
[2]	任娟, 张宁超, 刘萍萍. 碳气凝胶储氢性能的理论研究[J]. 计算物理, 2019, 36(6): 749-756.
[3]	刘聪, 张斌, 张亮, 郑君萧, 陈义学. S_N共轭函数用于蒙特卡罗粒子输运自动减方差的研究[J]. 计算物理, 2018, 35(5): 535-544.
[4]	徐小东, 周彬, 郭金川, 易明皓, 肖飞虎. 微结构X射线源中电子与阳极作用的蒙特卡罗模拟[J]. 计算物理, 2018, 35(1): 95-102.
[5]	李刚, 张宝印, 邓力, 上官丹骅, 李瑞, 马彦, 付元光, 胡小利. 蒙特卡罗区域分解并行计算中确保串并行结果一致的伪随机数应用[J]. 计算物理, 2017, 34(1): 67-72.
[6]	杨师俨, 何曼丽, 蒋丹枫, 王国辉, 苑超, 戴耀东. 反应堆²³⁵U裂变源辐射防护材料的优化设计[J]. 计算物理, 2017, 34(1): 73-81.
[7]	郑征, 王梦琪, 黎辉, 梅其良. 三维离散纵标-蒙特卡罗耦合方法在核电厂堆腔漏束计算中的应用[J]. 计算物理, 2016, 33(5): 599-605.
[8]	单卿, 蔡平坤, 褚胜男. PGNAA煤质在线测量中同时提高C、O和其它元素测量精度的慢化体模拟[J]. 计算物理, 2016, 33(5): 625-630.
[9]	王红霞, 刘代志, 宋仔标, 杨倪晨, 杨成莱. 激光在烟幕中传输的蒙特卡罗模拟[J]. 计算物理, 2013, 30(3): 415-421.
[10]	李满仓, 王侃, 姚栋. 连续能量蒙特卡罗方法组件均匀化中的临界效应[J]. 计算物理, 2012, 29(5): 727-732.
[11]	李树, 田东风, 邓力. 提高超高能中子计算效率的蒙特卡罗抽样技巧[J]. 计算物理, 2011, 28(3): 323-328.
[12]	尹延朋, 郑春, 黄坡. CFBR-Ⅱ堆增殖反应性转换系数计算[J]. 计算物理, 2010, 27(6): 799-804.
[13]	李树, 田东风, 邓力. 中子中子碰撞产生超高能中子问题模拟[J]. 计算物理, 2010, 27(5): 717-721.
[14]	许海燕, 黄正丰, 蔡少辉. 蒙特卡罗粒子输运问题的全局降方差方法[J]. 计算物理, 2010, 27(5): 722-732.
[15]	陶应龙, 朱金辉, 王建国, 牛胜利, 范如玉. 非均匀大气中γ射线输运的蒙特卡罗模拟方法[J]. 计算物理, 2010, 27(5): 740-744.