基于Geant4的反应堆临界及动态计算方法

Geant4是一款基于C++面向对象技术的蒙特卡罗开发程序包,可以模拟各种已知粒子与物质之间的相互作用。然而该程序包没有提供临界源功能,无法直接用于反应堆物理计算。因此,利用Geant4提供的基础物理模型和粒子跟踪控制等功能,用两种不同方法实现了临界源的设置,实现了基于Geant4的反应堆静态计算程序G4-RSM和反应堆动态计算程序G4-RDM。两个程序均可用于反应堆临界计算,与MCNP计算结果相对误差在5%以内。G4-RDM程序除可用于临界计算外,还可用于模拟堆内事故工况下的中子学瞬态变化。     

基于Geant4的反应堆临界及动态计算方法

Geant4是一款基于C++面向对象技术的蒙特卡罗开发程序包,可以模拟各种已知粒子与物质之间的相互作用。然而该程序包没有提供临界源功能,无法直接用于反应堆物理计算。因此,利用Geant4提供的基础物理模型和粒子跟踪控制等功能,用两种不同方法实现了临界源的设置,实现了基于Geant4的反应堆静态计算程序G4-RSM和反应堆动态计算程序G4-RDM。两个程序均可用于反应堆临界计算,与MCNP计算结果相对误差在5%以内。G4-RDM程序除可用于临界计算外,还可用于模拟堆内事故工况下的中子学瞬态变化。     

相对论高次谐波的蒙特卡罗模拟

根据相对论条件下的经典电磁辐射理论,采用蒙特卡罗计算方法对类氢离子在强激光场下辐射的高次谐波进行了模拟计算,获得了相对论领域内激光参数（频率、脉冲波形、强度等）对高次谐波影响的各种计算结果,并对结果进行了理论分析。     

相对论高次谐波的蒙特卡罗模拟

 根据相对论条件下的经典电磁辐射理论,采用蒙特卡罗计算方法对类氢离子在强激光场下辐射的高次谐波进行了模拟计算,获得了相对论领域内激光参数（频率、脉冲波形、强度等）对高次谐波影响的各种计算结果,并对结果进行了理论分析。     

400~1500 MeV质子轰击铅靶和钨靶的出射中子能谱的FLUKA和Geant4模拟研究

在ADS散裂靶系统的优化设计中,蒙特卡罗方法结合可靠的散裂反应理论模型进行中子学计算具有重要的作用。本工作利用Geant4程序中的INCLXX模型、BIC模型以及BERT模型和FLUKA程序分别模拟了597 MeV和1 500 MeV质子轰击薄铅靶不同出射角度的中子双微分截面,500,1 500 MeV质子轰击厚铅靶不同出射角度的中子双微分产额,以及400,600,800,1 000和1 200 MeV质子轰击厚钨靶在反角方向（175 °）的中子双微分产额,并与实验数据进行比较。研究表明,对于薄铅靶,Geant4程序的INCLXX模型和FLUKA程序模拟结果与实验符合得更好。能量在10～40 MeV范围内,BIC模型模拟结果明显高于实验数据,而BERT模型模拟结果略微低于实验数据。对于厚铅靶,在40 MeV左右所有的模拟结果都低于实验数据。对于厚钨靶,Geant4程序的BIC模型和FLUKA程序与实验数据符合得较好,INCLXX模型在能量高于60 MeV时模拟结果低于实验数据,BERT模型与实验数据差异较大。总体来看,Geant4程序的INCLXX模型和FLUKA程序进行ADS散裂靶相关的中子学的计算是合理和可靠的。The reliable Monte Carlo simulation codes coupled with nuclear reaction models play an important role in the neutronic calculation for the design and optimization of the ADS spallation target. In this work, the double differential cross sections at different angles produced from a thin lead target bombarded with 597 and 1 500 MeV protons, the neutron energy spectra at different angles produced from a thick lead target bombarded with 500 and 1 500 MeV protons, and the neutron energy spectra in the backward direction(175°) produced from a thick tungsten target bombarded with 400, 600, 800, 1 000 and 1 200 MeV protons are calculated with the Geant4 code coupled INCLXX, BIC and BERT models and the FLUKA code. The calculations are compared with the available experimental data. The results show that, for the thin lead target, the calculations with the Geant4 coupled INCLXX model and FLUKA code well reproduce the experimental results. In a energy range from 10 to 40 MeV, BIC model obviously overestimates the experimental results, and BERT model slightly underestimates the experimental results. For the thick lead target, all of the calculations underestimate the experimental results around 40MeV. For the thick tungsten target, the Geant4 coupled BIC model and FLUKA code well reproduce the experimental results. INCLXX model underestimates the experimental results above 60 MeV. BERT model bad reproduces the experimental results. Overall, the neutronic calculations with the Geant4 code coupled INCLXX model and FLUKA code for the ADS spallation target is reasonable and reliable.     

偏振X射线荧光系统测量痕量Cd元素的Geant4模拟研究

通过传统能量色散X射线荧光分析(EDXRF)对Cd元素进行痕量分析时,X光管发出的原级轫致辐射连续谱对Cd元素的分析具有严重的影响.为了减弱原级X射线对测量结果的影响,通过Geant4程序包模拟不同几何尺寸下偏振激发X射线荧光分析(P-EDXRF)中荧光靶的结构,研究其在减弱测量过程中原级X射线轫致辐射连续谱的影响.为...     

分子轨道波函数相对论效应的高分辨电子动量谱学观测

利用新研制的能量分辨为0.5 eV 的高分辨(e, 2e)谱仪测量了CF3I分子在9-15 eV能区的束缚能谱.新谱仪能较好地分辨该分子碘孤对轨道的自旋轨道劈裂组分5e3/2和5e1/2,并得到了它们各自的电子动量分布以及分支比随动量的变化关系.实验结果清楚地显示了5e3/2和5e1/2态波函数的相对论效应.这是第一个分子轨道波函数相对论效应的直接实验观测.     

分子轨道波函数相对论效应的高分辨电子动量谱学观测

利用新研制的能量分辨为0.5 eV 的高分辨(e, 2e)谱仪测量了CF3I分子在9—15 eV能区的束缚能谱.新谱仪能较好地分辨该分子碘孤对轨道的自旋轨道劈裂组分5e3/2和5e1/2,并得到了它们各自的电子动量分布以及分支比随动量的变化关系.实验结果清楚地显示了5e3/2和5e1/2态波函数的相对论效应.这是第一个分子轨道波函数相对论效应的直接实验观测.     

闪烁光纤在γ射线辐照下部分特性的蒙特卡罗模拟

闪烁光纤在射线成像方面的应用越来越广泛。为了进一步了解闪烁光纤在射线辐照下的基本特性,基于蒙特卡罗方法,利用计算机模拟分析了γ射线在闪烁光纤中的好事例率与光纤长度及射线能量的关系,能量沉积效率与光纤长度及射线能量的关系。此项工作对于闪烁光纤阵列在射线成像,剂量场测量等方面的研究很有价值。     

闪烁光纤在γ射线辐照下部分特性的蒙特卡罗模拟

 闪烁光纤在射线成像方面的应用越来越广泛。为了进一步了解闪烁光纤在射线辐照下的基本特性，基于蒙特卡罗方法，利用计算机模拟分析了γ射线在闪烁光纤中的好事例率与光纤长度及射线能量的关系，能量沉积效率与光纤长度及射线能量的关系。此项工作对于闪烁光纤阵列在射线成像，剂量场测量等方面的研究很有价值。     

Mass attenuation coefficients of composite materials by Geant4, XCOM and experimental data: comparative study

The main goal of this present study is focused on testing the applicability of Geant4 electromagnetic models for studying mass attenuations coefficients for different types of composite materials at 59.5, 80, 356, 661.6, 1173.2 and 1332.5 keV photon energies. The simulated results of mass attenuation coefficients were compared with the experimental and theoretical XCOM data for the same samples and a good agreement has been observed. The results indicate that this process can be followed to determine the data on the attenuation of gamma rays with the several energies in different materials. The modeling for photon interaction parameters was standard for any type of composite samples. The Geant4 code can be utilized for gamma ray attenuation coefficients for the sample at different energies, which may sometimes be impractical by experiment investigation.     

Influence of thyroid volume reduction on absorbed dose in 131I therapy studied by using Geant4 Monte Carlo simulation

A simulation study has been performed to quantify the effect of volume reduction on the thyroid absorbed dose per decay and to investigate the variation of energy deposition per decay due to β- and γ-activity of 131I with volume/mass of thyroid, for water, ICRP- and ICRU-soft tissue taken as thyroid material. A Monte Carlo model of the thyroid, in the Geant4 radiation transport simulation toolkit was constructed to compute the β- and γ-absorbed dose in the simulated thyroid phantom for various values of its volume. The effect of the size and shape of the thyroid on energy deposition per decay has also been studied by using spherical, ellipsoidal and cylindrical models for the thyroid and varying its volume in 1-25 cm3 range. The relative differences of Geant4 results for different models with each other and MCNP results lie well below 1.870%. The maximum relative difference among the Geant4 estimated results for water with ICRP and ICRU soft tissues is not more than 0.225%. S-values for ellipsoidal, spherical and cylindrical thyroid models were estimated and the relative difference with published results lies within 3.095%. The absorbed fraction values for beta particles show a good agreement with published values within 2.105% deviation. The Geant4 based simulation results of absorbed fractions for gammas again show a good agreement with the corresponding MCNP and EGS4 results （±6.667%） but have 29.032% higher values than that of MIRD calculated values. Consistent with previous studies, the reduction of the thyroid volume is found to have a substantial effect on the absorbed dose. Geant4 simulations confirm dose dependence on the volume/mass of thyroid in agreement with MCNP and EGS4 computed values but are substantially different from MIRD8 data. Therefore, inclusion of size/mass dependence is indicated for 131I radiotherapy of the thyroid.     

基于μ子动量谱的多群模型测量平板材料的厚度

在基于子多重库仑散射的断层成像技术中,天然宇宙射线中不同动量的子在同一目标对象中的散射特性存在明显的差异,这在一定程度上影响了成像分辨率。为了减少或消除这种不利影响,建立了一套基于天然宇宙射线子动量谱的多群模型,并检验其用于铅板厚度测量及平板材料识别中的性能和特征。结果表明,采用基于子动量谱的多群模型能够实现1~16 cm范围内铅板厚度的准确测量,相对误差小于5%。此外,该多群模型还可以通过计算辐射长度来识别未知平板的材料组成,尤其是对于高原子序数材料组成的平板,其辐射长度相对误差小于4%。     

Influence of thyroid volume reduction on absorbed dose in ~(131)I therapy studied by using Geant4 Monte Carlo simulation

A simulation study has been performed to quantify the effect of volume reduction on the thyroid absorbed dose per decay and to investigate the variation of energy deposition per decay due to β- and γ-activity of ~(131)I with volume/mass of thyroid,for water,ICRP- and ICRU-soft tissue taken as thyroid material.A Monte Carlo model of the thyroid,in the Geant4 radiation transport simulation toolkit was constructed to compute the β- and γ-absorbed dose in the simulated thyroid phantom for various values of its volume.The effect of the size and shape of the thyroid on energy deposition per decay has also been studied by using spherical,ellipsoidal and cylindrical models for the thyroid and varying its volume in 1-25 cm~3 range.The relative differences of Geant4 results for different models with each other and MCNP results he well below 1.870%.The maximum relative difference among the Geant4estimated results for water with ICRP and ICRU soft tissues is not more than 0.225%.S-values for ellipsoidal,spherical and cylindrical thyroid models were estimated and the relative difference with published results lies within3.095%.The absorbed fraction values for beta particles show a good agreement with published values within 2.105%deviation.The Geant4 based simulation results of absorbed fractions for gammas again show a good agreement with the corresponding MCNP and EGS4 results(±6.667%) but have 29.032%higher values than that of MIRD calculated values.Consistent with previous studies,the reduction of the thyroid volume is found to have a substantial effect on the absorbed dose.Geant4 simulations confirm dose dependence on the volume/mass of thyroid in agreement with MCNP and EGS4 computed values but are substantially diiferent from MIRD8 data.Therefore,inclusion of size/mass dependence is indicated for ~(131)I radiotherapy of the thyroid.     

Geant4在CR39叠层探测器中的应用

利用Geant4程序计算了多种不同组合CR39叠层探测器的中子能量响应函数,研究结果表明,中子能量响应函数曲线的半高宽和探测下限分别由聚乙烯的厚度和吸收片的厚度决定,峰值由聚乙烯和吸收片的厚度共同决定。利用加速器准单能氘氘中子源进行辐照实验,获得了4种不同叠层组合探测器的2.45 MeV中子探测效率,并与计算曲线进行了比较,结果吻合较好,验证了Geant4程序计算结果的可靠性。     

辐射成像系统模拟软件NUCRPD的研发与应用

在辐射成像系统的设计研发中,经常需要计算系统的各项物理性能指标。为了提升产品研发的效率,同方威视联合清华大学,基于蒙特卡罗程序包Geant4,研发了国内首套具有自主知识产权的辐射成像系统模拟软件NucRPD(NUCTECH Radiography Performance Design Tools),其能够快速、精确地模拟系统的各项性能指标。该软件对常用的一些辐射成像系统进行了参数化建模,用户通过修改软件界面上的少量参数,就可以快速建立各类辐射成像系统的几何体、源项、物理模型和统计量,然后在服务器上以并行计算方式完成模拟计算,能在较短的时间内模拟出系统的性能指标,并给出直观的图形帮助用户深入理解模拟结果。NucRPD的模拟结果经过了大量的实验验证,其剂量场分布和物理指标等模拟结果和实验结果符合得很好。NucRPD已经应用于同方威视辐射成像系统产品的设计研发,在产品的物理指标优化和辐射防护优化中发挥了重要作用。     

Photon attenuation coefficients of different rock samples using MCNPX,Geant4 simulation codes and experimental results: a comparison study

ABSTRACT

The photon attenuation coefficients for the rocks (Feldspathic basalt, Compact basalt, Volcanic rock, Pink granite, Sandstone and Dolerite) have been investigated using MCNPX and Geant4 codes for photon energies 122, 356, 511, 662, 1170, 1275 and 1330 keV. Comparison of the simulation and experimental results of mass attenuation coefficients is presented. The exposure buildup factors also have been calculated with the help of the G-P fitting method. The comparison showed reasonable agreement between the simulated and experimental data for all rocks. The results show that sandstone has the highest mean free path as well as exposure buildup factor which signifies that sandstone requires a larger thickness for protection from photons.