基于蒙特卡罗方法的几何因子计算程序研究
详细信息 查看官网全文
- 英文论文题名：Research on the geometric factor calculation program based on the Monte Carlo method
- 论文作者：张洁 ; 张莹 ; 陈秀莲 ; 庞蓓蓓 ; 白立新
- 英文论文作者：Zhang Jie ; Zhang Ying ; Chen Xiulian ; Pang Beibei ; Bai Lixin ; Key Laboratory for Neutron Physics of Chinese Academy of Engineering Physics Institute of Nuclear Physics and Chemistry ; College of Physical Science and Technology of Sichuan University
- 年：2014
- 作者机构：中国工程物理研究院中子物理学重点实验室核物理与化学研究所;四川大学物理科学与技术学院;
- 论文关键词：几何因子 ; 蒙特卡罗方法 ; 方差减小技巧 ; NDP ; MCNP5
- 英文论文关键词：Geometric factor ; the Monte Carlo method ; variance reduction method ; NDP ; MCNP5
- 会议召开时间：2014-08-13
- 会议录名称：第十七届全国核电子学与核探测技术学术年会论文集
- 语种：中文
- 分类号：TP312.1
- 学会代码：HDZX
- 会议名称：第十七届全国核电子学与核探测技术学术年会
- 会议地点：中国甘肃兰州
- 主办单位：中国电子学会、中国核学会核电子学与核探测技术分会
- 学会名称：中国电子学会核电子学与核探测技术分会
- 页数：9
- 文件大小：540k
- 原文格式：D
- 会议级别：全国
摘要
几何因子在某些放射测量系统和放射源活度的绝对刻度中发挥着重要作用。对于不能使用解析方法求解的几何因子的计算,往往可以使用蒙特卡罗方法解决。目前罕有文献报道对于圆面探测器对不共轴、不平行的圆面源几何因子的计算。该文研究的基于蒙特卡罗方法的几何因子计算程序使用C++语言编写,可用于各种尺寸的圆面探测器对圆面源各种情况下几何因子的计算。该程序使用了方差减小技巧,并与MCNP5的计算结果进行对比,具有结果准确、计算速度快、使用方便的优点。最终使用该程序成功验证了中子深度分布分析(NDP)能谱测量系统探测器位姿的准确性,并对其移动位置进行了粗略修正。
Geometric factor calculations play an important role in the absolute calibration of radioactivity measurement systems and in the determination of the activity of radioactive sources. The Monte Carlo method can always calculate the geometric factors when there is no analytical algorithm available for geometric factor calculations. There are few papers have reported the calculations of geometric factors when source and detector are both disc, but are un-coaxial and un-parallel. The program studied in this paper is based on the Monte Carlo method and is written in the C++ programming language, and, it can be used to calculate the geometric factor of a disc detector for a disc source, the detector and the source can be any size at any cases. This program integrated a variance reduction method, by comparing with the results calculated using MCNP5, it is can be concluded that the program calculates geometric factors accurately, fast and conveniently. The geometric factor calculation program verified the accuracy of the pose of the detector which applied in the NDP energy spectrum measurement system; furthermore, it corrected the detector position roughly.
Geometric factor calculations play an important role in the absolute calibration of radioactivity measurement systems and in the determination of the activity of radioactive sources. The Monte Carlo method can always calculate the geometric factors when there is no analytical algorithm available for geometric factor calculations. There are few papers have reported the calculations of geometric factors when source and detector are both disc, but are un-coaxial and un-parallel. The program studied in this paper is based on the Monte Carlo method and is written in the C++ programming language, and, it can be used to calculate the geometric factor of a disc detector for a disc source, the detector and the source can be any size at any cases. This program integrated a variance reduction method, by comparing with the results calculated using MCNP5, it is can be concluded that the program calculates geometric factors accurately, fast and conveniently. The geometric factor calculation program verified the accuracy of the pose of the detector which applied in the NDP energy spectrum measurement system; furthermore, it corrected the detector position roughly.
引文
[1]何希,雷震,巨凌军,等.α谱仪探测器和面源不同轴情况下几何因子的计算[J].原子能科学技术,2010,44():58-62.
[2]刘宗良.圆面探测器对圆面源(非同轴)的几何因子的蒙特卡罗计算[D].湖南:南华大学,2009.
[3]A.Dirican,P.E.Erden,M.Seferinoglu,等.The use of solid angle for alpha detector efficiency in226Ra analyses of soil samples[J].Applied Radiation and Isotopes,2012,114():1-3.
[4]S.Pommé.A complete series expansion of Ruby’s solid-angle formula[J].Nucl.Instr.and Meth.A,2004,531():616-620.
[5]Stefaan Pommé,Jan Paepen.A series expansion of Conway’s generalized formulas[J].Nucl.Instr.and Meth.A,2007,579():272-274.
[6]M.J.Prata.Solid angle subtended by a cylindrical detector at a point source in terms of elliptic integrals[].Radiation Physics and Chemistry,2003,67():599-603.
[7]Eduardo Galiano,Matthew Rodrigues.Acomparison of different analytical methods of determining the solid angle of a circular coaxial source-detector system[J].Applied Radiation and Isotopes,2006,64()497-501.
[8]S.Pommé.The solid angle subtended by a circular detector for a linear source[J].Applied Radiation and Isotopes,2007,65():724-727.
[9]牛胜利,彭玉,王建国,等.X光导管传输特性的蒙特卡罗模拟[J].强激光与粒子束,2014,,16(12):1513-1517.
[10]卢玉朝,谢奇林,宋凌莉,等.Rossi-α测量方法的蒙特卡罗直接模拟[J].原子能科学技术,2012,46():370-373.
[11]孙嘉龙,余纲林,佘顶,等.堆用蒙特卡罗程序几何重复结构功能开发[J].强激光与粒子束,2013,25(1):219-222.
[12]周宇璐,李仁顺,张宝玲,等.材料中He深度分布演化的Monte Carlo模拟研究[J].物理学报,2011,60(6):060702.
[13]陈楠,李成刚,戴文华,等.蒙特卡罗方法在X光源尺寸测量中的应用[J].强激光与粒子束,2008,20(6):1031-1033.
[14]刘宗良,周剑良.圆面探测器对圆面源几何因子的蒙特卡罗计算[J].科技资讯,2008,6():251-254.
[15]刘宗良.圆面探测器对圆面源(非同轴)的几何因子的加权蒙特卡罗计算[J].湖南人文科技学院学报,2008,4():1-3.
[16]J.R.Tickner.Algorithm for forcing scattered radiation to arbitrary convex regions in neutral particle Monte Carlo simulation[J].Nuclear Instruments and Methods in Physics Research B,2009,267():2361-2364.
[17]R.P.Gardner,H.K.Choi,M.Mickael,等.Algorithms for Forcing Scattered Radiation to Spherical,Planar Circular,and Right Circular Cylindrical Detectors for Monte Carlo simulation[J].Nuclear Science and Engineering,1987,95():245-256.
[18]Jie Zhang,Xiulian Chen,Changsheng Zhang,等.Development of a software package for solid-angle calculations using the Monte Carlo method[J].Nuclear Instruments and Methods in Physics Research A,2014,736():40-45.
[19]许淑艳.蒙特卡罗方法在实验核物理中的应用[M].北京:原子能出版社,2006.
[20]X-5 Monte Carlo Team.MCNP-A General Monte Carlo N-Particle Transport Code,Version 5.MCNP5用户使用说明,2003.
[21]维基百科.Solid Angle.http://en.wikipedia.org/wiki/Solid_angle,2014/5/26.
[22]陈飞达,汤晓斌,王鹏,等.基于蒙特卡罗方法的中子屏蔽材料设计[J].强激光与粒子束,2012,24(12):3006-3010.
[23]张信一,朱养妮,赵柱民,等.放射源辐射场蒙特卡罗模拟计算[J].强激光与粒子束,2013,25(1):223-226.
[24]马运星,白立新.MCNP程序在HPGe探测器马鞍形样品盒最佳选择中的应用[J].原子能科学技术,2009,43(11):1039-1042.
[25]朱金辉,牛胜利,马继明,等.γ射线在LSO晶体中的能量沉积[J].强激光与粒子束,2010,22(6):1351-1354.
[26]窦海峰,李润东,冷军,等.中子深度分析的能量展宽修正[J].核技术,2011,34(9):689-692.
[2]刘宗良.圆面探测器对圆面源(非同轴)的几何因子的蒙特卡罗计算[D].湖南:南华大学,2009.
[3]A.Dirican,P.E.Erden,M.Seferinoglu,等.The use of solid angle for alpha detector efficiency in226Ra analyses of soil samples[J].Applied Radiation and Isotopes,2012,114():1-3.
[4]S.Pommé.A complete series expansion of Ruby’s solid-angle formula[J].Nucl.Instr.and Meth.A,2004,531():616-620.
[5]Stefaan Pommé,Jan Paepen.A series expansion of Conway’s generalized formulas[J].Nucl.Instr.and Meth.A,2007,579():272-274.
[6]M.J.Prata.Solid angle subtended by a cylindrical detector at a point source in terms of elliptic integrals[].Radiation Physics and Chemistry,2003,67():599-603.
[7]Eduardo Galiano,Matthew Rodrigues.Acomparison of different analytical methods of determining the solid angle of a circular coaxial source-detector system[J].Applied Radiation and Isotopes,2006,64()497-501.
[8]S.Pommé.The solid angle subtended by a circular detector for a linear source[J].Applied Radiation and Isotopes,2007,65():724-727.
[9]牛胜利,彭玉,王建国,等.X光导管传输特性的蒙特卡罗模拟[J].强激光与粒子束,2014,,16(12):1513-1517.
[10]卢玉朝,谢奇林,宋凌莉,等.Rossi-α测量方法的蒙特卡罗直接模拟[J].原子能科学技术,2012,46():370-373.
[11]孙嘉龙,余纲林,佘顶,等.堆用蒙特卡罗程序几何重复结构功能开发[J].强激光与粒子束,2013,25(1):219-222.
[12]周宇璐,李仁顺,张宝玲,等.材料中He深度分布演化的Monte Carlo模拟研究[J].物理学报,2011,60(6):060702.
[13]陈楠,李成刚,戴文华,等.蒙特卡罗方法在X光源尺寸测量中的应用[J].强激光与粒子束,2008,20(6):1031-1033.
[14]刘宗良,周剑良.圆面探测器对圆面源几何因子的蒙特卡罗计算[J].科技资讯,2008,6():251-254.
[15]刘宗良.圆面探测器对圆面源(非同轴)的几何因子的加权蒙特卡罗计算[J].湖南人文科技学院学报,2008,4():1-3.
[16]J.R.Tickner.Algorithm for forcing scattered radiation to arbitrary convex regions in neutral particle Monte Carlo simulation[J].Nuclear Instruments and Methods in Physics Research B,2009,267():2361-2364.
[17]R.P.Gardner,H.K.Choi,M.Mickael,等.Algorithms for Forcing Scattered Radiation to Spherical,Planar Circular,and Right Circular Cylindrical Detectors for Monte Carlo simulation[J].Nuclear Science and Engineering,1987,95():245-256.
[18]Jie Zhang,Xiulian Chen,Changsheng Zhang,等.Development of a software package for solid-angle calculations using the Monte Carlo method[J].Nuclear Instruments and Methods in Physics Research A,2014,736():40-45.
[19]许淑艳.蒙特卡罗方法在实验核物理中的应用[M].北京:原子能出版社,2006.
[20]X-5 Monte Carlo Team.MCNP-A General Monte Carlo N-Particle Transport Code,Version 5.MCNP5用户使用说明,2003.
[21]维基百科.Solid Angle.http://en.wikipedia.org/wiki/Solid_angle,2014/5/26.
[22]陈飞达,汤晓斌,王鹏,等.基于蒙特卡罗方法的中子屏蔽材料设计[J].强激光与粒子束,2012,24(12):3006-3010.
[23]张信一,朱养妮,赵柱民,等.放射源辐射场蒙特卡罗模拟计算[J].强激光与粒子束,2013,25(1):223-226.
[24]马运星,白立新.MCNP程序在HPGe探测器马鞍形样品盒最佳选择中的应用[J].原子能科学技术,2009,43(11):1039-1042.
[25]朱金辉,牛胜利,马继明,等.γ射线在LSO晶体中的能量沉积[J].强激光与粒子束,2010,22(6):1351-1354.
[26]窦海峰,李润东,冷军,等.中子深度分析的能量展宽修正[J].核技术,2011,34(9):689-692.