面向个人剂量计校准的小尺度伽玛参考辐射模拟研究

针对伽玛射线个人剂量计基于标准伽玛参考辐射进行校准时检定效率低、校准工作复杂和需要远程送检的关键技术问题,建构了1 Ci ¹³⁷Cs放射源小尺度参考辐射场物理模型,采用蒙特卡罗方法,研究了小尺度参考辐射场内的剂量分布、装置结构和待检剂量计变化导致散射射线对剂量场的影响,获得了待检剂量计形状、数量、类型和装置结构产生的散射伽玛射线对小尺度参考辐射量值定度的影响结果。研究结果表明,1 Ci ¹³⁷Cs可以为小尺度参考辐射辐射场检验点提供1.5 mSv/h的伽玛遂行剂量率,辐照个人剂量计载台直径30 cm束斑上的剂量率相对标准偏差约为0.48%。当载台厚度为20 mm时,散射射线对小尺度参考辐射检验点处剂量率值的影响率为3.27%,高于剂量计尺寸（1.62%）和剂量计数量（0.56%）的影响。     

锶热原子束二维准直的动力学过程的蒙特卡罗模拟及实验研究

在考虑随机因素的情况下, 应用蒙特卡罗方法在理论上详细研究了锶原子束二维激光准直的动力学过程. 综合考虑原子横向发散角、初始原子位置、纵向速度分布、同位素等因素, 获得了激光二维准直后的原子横向空间分布的模拟结果以及随准直光失谐、光功率等参量因素的变化. 通过与实验数据比较, 理论值和实验值很好相符, 显示蒙特卡罗方法可以精确地描述锶原子束二维准直的动力学过程. 为原子束激光二维准直的精确控制, 高精度原子钟系统的优化, 提供了一种理论分析方法. 关键词： 二维准直 蒙特卡罗方法 横向空间分布     

基于法国实验客体的X射线准直技术研究

通过研究法国实验客体(FTO)的闪光照相过程,设计了静态照相过程中采用的准直措施,除传统的一二级准直体外,还通过添加三级准直体来进一步控制入射到客体边缘部分的X射线,从而降低散射的影响。蒙特卡罗方法模拟结果表明,采用三级准直技术可以有效地控制散射的影响,图像接收平面各边界的直散比有了明显的提升。     

基于法国实验客体的X射线准直技术研究

通过研究法国实验客体（FTO）的闪光照相过程, 设计了静态照相过程中采用的准直措施, 除传统的一二级准直体外, 还通过添加三级准直体来进一步控制入射到客体边缘部分的X射线, 从而降低散射的影响。蒙特卡罗方法模拟结果表明, 采用三级准直技术可以有效地控制散射的影响, 图像接收平面各边界的直散比有了明显的提升。     

基于蒙特卡罗的小尺度参考辐射装置屏蔽研究

针对移动式小尺度参考辐射（MRR）装置（移动式校准装置）,在进行射线辐射剂量测量的仪器仪表标定或刻度时,应满足其辐射屏蔽安全限值5 Sv/h的屏蔽技术要求,采用蒙特卡罗输运程序MCNP,开展了移动式小尺度参考辐射装置表面剂量场屏蔽的模拟计算和研究分析工作。研究结果表明,通过MCNP模拟的屏蔽设计方法可以详尽反映MRR装置各个表面的剂量分布特征和规律,实现移动式小尺度参考辐射装置屏蔽设计,采用的铅钢材料复合屏蔽方案能够保证装置硬度且显著地减轻屏蔽体的重量,最终获取的优化MRR屏蔽箱体重量约为271.9 kg。     

基于蒙特卡罗的小尺度参考辐射装置屏蔽研究

针对移动式小尺度参考辐射（MRR）装置（移动式校准装置）,在进行射线辐射剂量测量的仪器仪表标定或刻度时,应满足其辐射屏蔽安全限值5 Sv/h的屏蔽技术要求,采用蒙特卡罗输运程序MCNP,开展了移动式小尺度参考辐射装置表面剂量场屏蔽的模拟计算和研究分析工作。研究结果表明,通过MCNP模拟的屏蔽设计方法可以详尽反映MRR装置各个表面的剂量分布特征和规律,实现移动式小尺度参考辐射装置屏蔽设计,采用的铅钢材料复合屏蔽方案能够保证装置硬度且显著地减轻屏蔽体的重量,最终获取的优化MRR屏蔽箱体重量约为271.9 kg。     

毛细管准直微束过程模拟与微束品质分析

利用GEANT4程序对复旦大学单粒子微束的毛细管准直过程进行了蒙特卡罗模拟。计算结果表明：当前采用的1 m孔径、1 mm长度的毛细玻璃管准直器能够引出峰值能量2.2 MeV、能量分辨130 keV、束径2.4 m的质子微束,可达到装置对微束能量与束径分辨的设计指标,从理论上明确了毛细管准直微束引出的可行性。针对影响入射目标细胞微束品质的主要因素,分析了微束品质与准直系统和入射束流各状态参数的依赖关系,为毛细管准直微束的引出与优化提供了必要的理论依据。     

毛细管准直微束过程模拟与微束品质分析

利用GEANT4程序对复旦大学单粒子微束的毛细管准直过程进行了蒙特卡罗模拟。计算结果表明:当前采用的1 m孔径、1 mm长度的毛细玻璃管准直器能够引出峰值能量2.2 MeV、能量分辨130 keV、束径2.4 m的质子微束,可达到装置对微束能量与束径分辨的设计指标,从理论上明确了毛细管准直微束引出的可行性。针对影响入射目标细胞微束品质的主要因素,分析了微束品质与准直系统和入射束流各状态参数的依赖关系,为毛细管准直微束的引出与优化提供了必要的理论依据。     

一个模拟光子照射人体的蒙特卡罗程序

本文分绍了一个用于辐射防护和医疗目的的蒙特卡罗程序。它能计算出受到光子照射时人体20种器官(组织)的吸收剂量和人体有效剂量当量。程序中使用了一种特殊的方法以加快计算结果的收敛。本文所介绍的方法也可用于其他复杂几何系统的剂量计算。     

蒙特卡罗模拟分析电子束发射度对照射量空间分布影响

加速器靶前正前方1 m处的照射量是衡量加速器光源辐射能力的重要物理参量。电子束轰击高原子序数靶产生的X射线空间分布具有很强的前冲性,照射量空间分布与电子束的发射度密切相关。采用高斯函数对不同电子束发射度条件下入射电子的空间分布和角度分布进行建模。应用蒙特卡罗方法对电子束打靶的轫致辐射过程进行模拟,分析电子束发射度对照射量的影响。同时还对整靶结构和叠靶结构下的轫致辐射光源照射量进行计算比较。结果表明,电子束的发散角是影响轫致辐射光源照射量的主要因素。与采用整靶结构相比,采用叠靶结构所获得的照射量空间分布基本一致,在正前方小角度范围内（0~4）的照射量有2%~3%的降低。     

非定常辐射输运问题的蒙特卡罗自适应偏倚抽样

用蒙特卡罗(MC)方法模拟高温、高压、多介质、大变形辐射输运问题时,由于网格体积悬殊,导致各网格通量的统计误差涨落很大,随着时间步的增加,误差积累甚至会导致计算结果失真.为此,发展了针对全局网格计算的源偏倚抽样技巧.用于源偏倚抽样的价值函数基于上个时间步各网格通量及误差,通过加权构造产生,它比传统MC通过解伴随方程获取价值的性价比要高得多.数值试验表明,全局源偏倚抽样通过自适应分配当前时间步各网格的粒子数,有效地降低了当前步重要网格通量误差. 关键词： 非定常 辐射输运 蒙特卡罗 源偏倚抽样     

Study on CYCIAE-100 radiation field and residual radioactivity

The accelerators should be properly designed to make the radiation field produced by beam loss satisfy the dose limits. The radiation field for high intensity H- cyclotron includes prompt radiation and residual radiation field. The induced radioactivity in accelerator components is the dominant source of occupational radiation exposure if the accelerator is well shielded. The source of radiation is the beam loss when cyclotron is operating. In this paper, the radiation field for CYCIAE-100 is calculated using Monte Carlo method and the radioactive contamination near stripping foil is studied. A method to reduce the dose equivalent rate of maintenance staff is also given.     

Study on CYCIAE-100 radiation field and residual radioactivity

The accelerators should be properly designed to make the radiation field produced by beam loss satisfy the dose limits. The radiation field for high intensity H^－ cyclotron includes prompt radiation and residual radiation field. The induced radioactivity in accelerator components is the dominant source of occupational radiation exposure if the accelerator is well shielded. The source of radiation is the beam loss when cyclotron is operating. In this paper, the radiation field for CYCIAE-100 is calculated using Monte Carlo method and the radioactive contamination near stripping foil is studied. A method to reduce the dose equivalent rate of maintenance staff is also given.     

高空核爆炸瞬发辐射电离效应的数值模拟

采用蒙特卡罗方法模拟了高空核爆炸瞬发辐射中子、γ射线、X射线电离大气的过程,给出了几种爆炸场景下瞬发辐射产生的附加电离电子密度空间分布.针对大气密度随高度非均匀连续变化的特性,采用质量距离抽样方法取代常用的步长抽样方法,无需根据大气密度随高度的变化进行分层处理,提高了计算效率.结果表明:对于不同的爆高,瞬发辐射电离分布存在显著的差异;随着爆高的增加,瞬发辐射附加电离区范围增大,但电子密度的峰值减小. 关键词： 高空核爆炸 瞬发辐射 大气电离 蒙特卡罗方法     

基于小尺度参考辐射的便携式伽玛射线剂量(率)仪校准方法及应用

针对现有伽玛射线剂量(率)测量仪器检定方法送检距离远、周期长、效率低的不足,中国工程物理研究院核物理与化学研究所采用基于样本仪表的机器预测方法,开展了基于小尺度参考辐射的便携式伽玛剂量(率)仪遂行校准技术研究工作。最终成功地研制了一套辐射防护用伽玛射线剂量(率)仪遂行校准装置,实现了测量标准不确定度不大于5%的良好结果。为伽玛空气比释动能量值的遂行定度提出了一条全新的技术思路,为便携式伽玛剂量(率)仪的现场和野外校准提供了一种全新的示范装备。     

紧凑型THz源的史密斯-帕塞尔辐射模拟

 采用模拟和数值计算的方法,研究了THz波段的受激史密斯-帕塞尔辐射特性。实验装置以“上海电子束离子阱”为原型,采用紧凑型设计以便最终实现其可移动性。束流动力学模拟表明,此装置采用强磁场,可以得到平均流强为0.2 A、束流半径为75 μm的高品质电子束,为电子束工作在自由电子激光模式下创造了条件。基于Andrews和 Brau的理论,优化了光栅参数,保证了辐射角度在60°。其中消散场的计算频率为0.365 9 THz。采用particle-in-cell(PIC)程序模拟了光栅表面的辐射场以及电子的动力学特性。模拟结果表明电子有群聚效应,且二次谐波（0.723 THz,约为消散频率的2倍）得到增强。采用后处理方法计算了史密斯-帕塞尔辐射的功率空间分布。计算显示辐射角度与理论角度相一致,表明了方法的有效性。输出的功率约为2 mW。     

基于蒙卡模拟的分段γ扫描无源效率刻度方法

针对200 L核废物桶分段γ扫描（SGS）过程中的效率刻度问题,提出了一种效率刻度函数模型,采用MCNP程序计算不同基质密度和γ射线能量条件下的离散断层效率,经过多元非线性回归获取函数参数,从而建立效率刻度函数,实现核废物桶SGS断层效率刻度。对核废物桶样品进行实验分析,结果表明：对于桶内基质分别为密度0.310 g·cm⁻³的硅酸铝、密度0.595 g·cm⁻³的木质纤维,桶内核素分别为活度3.110×10⁵ Bq的点源¹³⁷Cs、活度1.371×10⁵ Bq的点源⁶⁰Co,在桶内仅有单个点源存在的核素分布极端不均匀情况下,桶内核素活度重建误差在−37.68%～31.52%范围内。本文的方法能够准确有效实现核废物桶SGS断层效率矩阵计算,并确定核废物桶内放射性核素活度,满足实际检测要求。     

Simulation of the infrared radiation characteristics of high-temperature exhaust plume including particles using the backward Monte Carlo method

The analysis of the infrared radiation characteristics of high-temperature free-stream flow including particles is very significant for the field of target detection, combustion diagnosis and temperature measurement of flame. In this paper, the infrared radiation characteristics of high-temperature free-stream flow are calculated and analyzed using the backward Monte Carlo method, considering the effect of the directional radiation heat flux due to the particle scattering and the different boundary conditions. The calculation results of emitting, absorbing and anisotropically scattering media are compared with the forward Monte Carlo and finite-volume methods results, which shows the superiority on computational efficiency with the backward Monte Carlo method.