x

## 留言板

 引用本文:
 Citation:

## Efficient Monte Carlo algorithm of time-dependent particle transport problem in multi-physics coupling calculation

Shangguan Dan-Hua, Yan Wei-Hua, Wei Jun-Xia, Gao Zhi-Ming, Chen Yi-Bing, Ji Zhi-Cheng
PDF
HTML
• #### 摘要

多物理耦合计算在众多领域都有重要应用. 如果其包含粒子输运过程, 用蒙特卡罗方法模拟粒子输运常占据大部分的计算时间, 因此多物理耦合计算中动态输运问题的高效蒙特卡罗模拟方法意义重大, 其不可避免地依赖于大规模并行. 基于动态输运问题的特点, 本文提出了两种新方法: 一是针对输运燃耗耦合计算的新型计数规约算法; 二是动态输运计算样本数自适应算法. 两种算法都能在保持计算结果基本不变的前提下使计算时间大幅减少, 从而提高了效率.

#### Abstract

Multi-physics coupling calculation has applications in many important research fields. If particle transport process is included in this calculation, Monte Carlo method is often used to simulate this process and usually a large amount of calculation time is needed. So, efficient Monte Carlo algorithm for time-dependent particle transport problem is important for an efficiently coupling calculation, which inevitably relies on large-scale parallel calculation. Based on the characteristic of time-dependent particle transport problem, two methods are proposed in this paper to achieve high- efficiency calculation. One is a tally-reducing algorithm which is used in the coupling of transport simulation and burnup calculation. By reducing the quantity of data which should be reduced necessarily, this method can reduce the calculation time largely. It can be seen that a new coupling mode for these two processes in MPI environment has a larger value when model scale is larger than the sample size. The other method is an adaptive method of setting the sample size of Monte Carlo simulation. The law of large number assures that the Monte Carlo method will obtain an exact solution when the sample scale tends to infinity. But generally, no one knows which sample scale is big enough for obtaining a solution with target precision in advance. So, the common strategy is to set a huge-enough sample scale by experience and conduct the posterior check for all results. Apparently, this way cannot be efficient because the calculation will go on after the precision of solution has reached an object value. Another popular method is to set the sample size to rely on the relative error of some single calculation. The sample size is enlarged without a break until the relative error is less than some presetting value. This method is not suitable either, because Monte Carlo particle transport simulation will gives feedbacks to other process which is composed of many tallies. It is inappropriate to adjust the sample size according to the relative error of any calculation. Relying on the generalization of the Shannon entropy concept and an on-the-fly diagnosis rule for a entropy value sequence, the adaptive method proposed in this paper can reduce the original huge sample scale to a reasonable level. By numerically testing some non-trivial examples, both algorithms can reduce the calculation time largely, with the results kept almost unchanged, so the efficiency is high in these cases.

#### 作者及机构信息

###### 通信作者: 姬志成, ji_zc@139.com
• 基金项目: 中国工程物理研究院科学基金(批准号: CX20200028)和国家自然科学基金青年科学基金(批准号: 11705011)资助的课题

#### Authors and contacts

###### Corresponding author: Ji Zhi-Cheng, ji_zc@139.com
• Funds: Project supported by the Science Foundation of China Academy of Engineering Physics, China (Grant No. CX20200028) and the Young Scientists Fund of the National Natural Science Foundation of China (Grant No. 11705011)

#### 施引文献

• 图 1  新旧方法计数与燃耗计算时间之和的对比

Fig. 1.  Comparison of the total time of tally reduce and burnup calculation by using old and new methods.

•  [1] 张鸿, 郭红霞, 潘霄宇, 雷志峰, 张凤祁, 顾朝桥, 柳奕天, 琚安安, 欧阳晓平. 重离子在碳化硅中的输运过程及能量损失. 物理学报, 2021, 70(16): 162401. doi: 10.7498/aps.70.20210503 [2] 邓力, 李瑞, 王鑫, 付元光. 特征γ射线谱分析的蒙特卡罗模拟技术. 物理学报, 2020, 69(11): 112801. doi: 10.7498/aps.69.20200279 [3] 陈忠, 赵子甲, 吕中良, 李俊汉, 潘冬梅. 基于蒙特卡罗-离散纵标方法的氘氚激光等离子体聚变反应率数值模拟. 物理学报, 2019, 68(21): 215201. doi: 10.7498/aps.68.20190440 [4] 上官丹骅, 姬志成, 邓力, 李瑞, 李刚, 付元光. 蒙特卡罗临界计算全局计数问题新策略研究. 物理学报, 2019, 68(12): 122801. doi: 10.7498/aps.68.20182276 [5] 李洪, 艾倩雯, 汪鹏君, 高和蓓, 崔毅, 罗孟波. 外力驱动作用下高分子链在表面吸附性质的计算机模拟. 物理学报, 2018, 67(16): 168201. doi: 10.7498/aps.67.20180468 [6] 李树. 光子与相对论麦克斯韦分布电子散射截面的蒙特卡罗计算方法. 物理学报, 2018, 67(21): 215201. doi: 10.7498/aps.67.20180932 [7] 上官丹骅, 邓力, 李刚, 张宝印, 马彦, 付元光, 李瑞, 胡小利. 蒙特卡罗临界计算全局计数效率新算法研究. 物理学报, 2016, 65(6): 062801. doi: 10.7498/aps.65.062801 [8] 上官丹骅, 邓力, 张宝印, 姬志成, 李刚. 非定常输运问题适应于消息传递并行编程环境的香农熵计算方法. 物理学报, 2016, 65(14): 142801. doi: 10.7498/aps.65.142801 [9] 林舒, 闫杨娇, 李永东, 刘纯亮. 微波器件微放电阈值计算的蒙特卡罗方法研究. 物理学报, 2014, 63(14): 147902. doi: 10.7498/aps.63.147902 [10] 杨亮, 魏承炀, 雷力明, 李臻熙, 李赛毅. 两相钛合金再结晶退火组织与织构演变的蒙特卡罗模拟. 物理学报, 2013, 62(18): 186103. doi: 10.7498/aps.62.186103 [11] 文德智, 卓仁鸿, 丁大杰, 郑慧, 成晶, 李正宏. 蒙特卡罗模拟中相关变量随机数序列的产生方法. 物理学报, 2012, 61(22): 220204. doi: 10.7498/aps.61.220204 [12] 李鹏, 许州, 黎明, 杨兴繁. 金刚石薄膜中二次电子输运的蒙特卡罗模拟. 物理学报, 2012, 61(7): 078503. doi: 10.7498/aps.61.078503 [13] 赵学峰, 李三伟, 蒋刚, 王传珂, 李志超, 胡峰, 李朝光. 超热电子与金黑腔靶作用产生硬X射线的蒙特卡罗模拟. 物理学报, 2011, 60(7): 075203. doi: 10.7498/aps.60.075203 [14] 张鹏飞, 苏兆锋, 孙剑锋, 杨海亮, 李永东, 高屹, 孙江, 王洪广, 尹佳辉, 梁天学, 孙凤举, 王志国. 阳极杆箍缩二极管产生X射线能谱的模拟计算. 物理学报, 2011, 60(10): 100204. doi: 10.7498/aps.60.100204 [15] 张宝武, 张萍萍, 马艳, 李同保. 铬原子束横向一维激光冷却的蒙特卡罗方法仿真. 物理学报, 2011, 60(11): 113701. doi: 10.7498/aps.60.113701 [16] 代秋声, 漆玉金. 针孔单光子发射计算机断层成像的空间分辨率研究. 物理学报, 2010, 59(2): 1357-1365. doi: 10.7498/aps.59.1357 [17] 金晓林, 黄桃, 廖平, 杨中海. 电子回旋共振放电中电子与微波互作用特性的粒子模拟和蒙特卡罗碰撞模拟. 物理学报, 2009, 58(8): 5526-5531. doi: 10.7498/aps.58.5526 [18] 赵宗清, 丁永坤, 谷渝秋, 王向贤, 洪 伟, 王 剑, 郝轶聃, 袁永腾, 蒲以康. 超短超强激光与铜靶相互作用产生Kα源的蒙特卡罗模拟. 物理学报, 2007, 56(12): 7127-7131. doi: 10.7498/aps.56.7127 [19] 孙贤明, 韩一平, 史小卫. 降雨融化层后向散射的蒙特卡罗仿真. 物理学报, 2007, 56(4): 2098-2105. doi: 10.7498/aps.56.2098 [20] 郝樊华, 胡广春, 刘素萍, 龚 建, 向永春, 黄瑞良, 师学明, 伍 钧. 钚体源样品γ能谱计算的蒙特卡罗方法. 物理学报, 2005, 54(8): 3523-3529. doi: 10.7498/aps.54.3523
• 文章访问数:  1424
• PDF下载量:  58
• 被引次数: 0
##### 出版历程
• 收稿日期:  2021-08-10
• 修回日期:  2022-01-24
• 上网日期:  2022-01-26
• 刊出日期:  2022-05-05

## 多物理耦合计算中动态输运问题高效蒙特卡罗模拟方法

• 北京应用物理与计算数学研究所, 北京　100094
• ###### 通信作者: 姬志成, ji_zc@139.com
基金项目: 中国工程物理研究院科学基金(批准号: CX20200028)和国家自然科学基金青年科学基金(批准号: 11705011)资助的课题

/