离子输运蒙特卡罗模拟的步长选取

计算物理 ›› 2014, Vol. 31 ›› Issue (4): 417-423.

离子输运蒙特卡罗模拟的步长选取

杨鑫^1,2, 王冠博^1,2, 李润东¹, 刘汉刚¹, 王侃²

1. 核物理与化学研究所, 四川绵阳 621900;
2. 清华大学工程物理系, 北京 100084

收稿日期:2012-11-27 修回日期:2013-05-17 出版日期:2014-07-25 发布日期:2014-07-25
作者简介:杨鑫(1988-),男,贵州,博士生,主要从事反应堆物理及应用研究,E-mail:yangx05@126.com
基金资助:
中国工程物理研究院科学技术发展基金(2010B0103009);国家自然科学基金(11175162)资助项目

Flight Path Selection in Monte Carlo Code ITR for Ion Transport

YANG Xin^1,2, WANG Guanbo^1,2, LI Rundong¹, LIU Hangang¹, WANG Kan²

1. Institution of Nuclear Physics and Chemistry, Mianyang, Sichuan 621900, China;
2. Department of Engineering Physics, Tsinghua University, Beijing 100084, China

Received:2012-11-27 Revised:2013-05-17 Online:2014-07-25 Published:2014-07-25

摘要/Abstract

摘要： 介绍带电离子在物质中输运模拟的详细历史法,开发模拟程序ITR(ions transport and reaction),同时实现抽取离子飞行步长的基本方法和脉冲近似法,计算结果与CORTEO及TRIM程序的结果比较符合很好.讨论两种步长选取方法,对脉冲近似法的最大步长进行限制,分析CORTEO程序与TRIM程序计算结果差别的主要原因.ITR同时具有TRIM和CORTEO的优点,采用脉冲近似法减少需要计算的核碰撞次数,通过索引(Index)方法对散射角的计算进行加速,计算效率大幅提高.

关键词: ITR, 离子输运, 详细历史法, 蒙特卡罗

Abstract: Strategy of detailed history method is summarized, depending on which code ITR(ions transport and reaction) was developed. Both basic method and impulse approximation method to select flight path of ions are complemented respectively. Numerical results agree with CORTEO and TRIM very well. Relation between two methods is discussed. We found reason of CORTEO results differ from TRIM's by limiting maximum of flight length derived by impluse approximation method. ITR has advantages of TRIM and CORTEO. Impluse approximation method reduces nuclear collisions simulated. Index technique is used to calculate scattering angle, which significantly improves efficiency of ITR.

Key words: ITR, ions transport, detailed history method, monte carlo

中图分类号:

TL99

杨鑫, 王冠博, 李润东, 刘汉刚, 王侃. 离子输运蒙特卡罗模拟的步长选取[J]. 计算物理, 2014, 31(4): 417-423.

YANG Xin, WANG Guanbo, LI Rundong, LIU Hangang, WANG Kan. Flight Path Selection in Monte Carlo Code ITR for Ion Transport[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2014, 31(4): 417-423.

[1]	王敏, 申玉清, 陈震宇, 徐鹏. 随机多孔介质的蒙特卡罗重构与渗流特性模拟[J]. 计算物理, 2021, 38(5): 623-630.
[2]	杨波, 衷斌, 徐琪, 成立, 潘流俊, 沈华韵. (α, n)反应中子产额与能谱计算方法[J]. 计算物理, 2021, 38(4): 393-400.
[3]	邓楚楚, 张鑫源, 陈水源, 吴宇晗, 叶晴莹, 黄志高. 小缺陷铁纳米环的磁特性[J]. 计算物理, 2020, 37(4): 473-478.
[4]	任娟, 张宁超, 刘萍萍. 碳气凝胶储氢性能的理论研究[J]. 计算物理, 2019, 36(6): 749-756.
[5]	吴言, 刘思琪, 李胜强. 芯片表面的冷极性分子静电阱与微阱阵列[J]. 计算物理, 2019, 36(4): 483-490.
[6]	王一, 潘流俊, 王瑞宏. 求解α本征值的k-α迭代方法误差估计[J]. 计算物理, 2019, 36(1): 39-46.
[7]	刘聪, 张斌, 张亮, 郑君萧, 陈义学. S_N共轭函数用于蒙特卡罗粒子输运自动减方差的研究[J]. 计算物理, 2018, 35(5): 535-544.
[8]	马先林, 周德胜. 基于水平集的两步MCMC方法对河道形态的识别[J]. 计算物理, 2018, 35(3): 321-329.
[9]	彭钢. 连续能量共轭加权蒙特卡罗动态参数计算[J]. 计算物理, 2018, 35(1): 87-94.
[10]	徐小东, 周彬, 郭金川, 易明皓, 肖飞虎. 微结构X射线源中电子与阳极作用的蒙特卡罗模拟[J]. 计算物理, 2018, 35(1): 95-102.
[11]	李胜强. 一种针对弱场搜寻态冷极性分子的可控静电表面囚禁方案[J]. 计算物理, 2017, 34(6): 731-739.
[12]	李刚, 张宝印, 邓力, 上官丹骅, 李瑞, 马彦, 付元光, 胡小利. 蒙特卡罗区域分解并行计算中确保串并行结果一致的伪随机数应用[J]. 计算物理, 2017, 34(1): 67-72.
[13]	杨师俨, 何曼丽, 蒋丹枫, 王国辉, 苑超, 戴耀东. 反应堆²³⁵U裂变源辐射防护材料的优化设计[J]. 计算物理, 2017, 34(1): 73-81.
[14]	上官丹骅, 许海燕. 非定常输运模拟中基于粒子标识分类的源偏倚算法[J]. 计算物理, 2016, 33(6): 639-644.
[15]	刘雄国, 邓力, 胡泽华, 李瑞, 付元光, 李刚. 基于VENUS-III国际基准模型的JMCT程序验证[J]. 计算物理, 2016, 33(5): 570-580.