EGSWIN计算系统的开发

针对EGS4在用户接口方面存在的缺陷，利用VC^ 、VF开发了基于Windows平台的EGSWIN计算程序包，实现了图形化输入，几何区域处理，并行计算，多记录方式输出，几何实体与粒子径迹显示等，极大地克服了原有EGS4使用上的局限，体现了当前国际上EGS的开发方向，并在一定程序上推动了EGS系列的发展。     

6 MV医用电子直线加速器轫致辐射谱研究

用衰减法对6MV医用电子直线加速器的轫致辐射X射线能谱进行测量，并利用模拟电子—光子耦合输运的Monte-Carlo程序EGS4对已知的治疗头几何结构、靶材料和电子源等特征的医用电子直线加速器轫致辐射谱进行研究。实验测量结果和EGS4模拟计算结果符合较好。     

蒙特卡罗程序EGS4在γ照射量积累因子计算中的应用

使用蒙特卡罗程序EGS4计算了入射光子能量范围从0.015 MeV到15 MeV、屏蔽厚度达到40个平均自由程的普通混凝土的γ照射量积累因子,程序计算中使用的是PHOTX光子截面库、各向同性点源和无限均匀介质模型,针对深穿透问题采用了粒子分裂技巧,并考虑了轫致辐射、荧光效应和相干(瑞利)散射对积累因子的影响。并将本文的计算结果与相关文献进行了比较分析,表明本文使用EGS4计算的积累因子数据是可靠的,可以更广泛的使用蒙特卡罗方法来计算γ照射量积累因子。     

在MARC程序中导入ANSYS程序三维网格

MARC程序和ANSYS程序，都是我国核安全局认可的可以用于核电站结构和设备分析的国际知名的大型通用有限元分析程序。在实际工作的应用中发现MARC程序建立三维实体元模型操作复杂，而使用ANSYS程序建立三维实体元模型的操作就很简单。因此，使用ANSYS的APDL语言编制专用程序，可实现ANSYS程序三维实体元网格转入到MARC程序中的功能。     

用于热核巨共振研究的大面积BGO闪烁探测器

报道了大面积BGO闪烁控制器用于132MeV^16O ^197Au反应中形成的热核巨共振研究，利用飞行时间法鉴别中子、γ。用EGS4程序模拟探讨器对γ连续谱的响应，利用BGO闪烁探测器中Ge俘获中子发射的γ对探测器进行了能量刻度。     

蒙特卡罗方法计算高纯锗γ谱仪对面源的探测效率和符合相加修正

本工作使用蒙特卡罗方法，借助于EGS4软件包，通过自编用户程序，研究γ谱仪对面源的效率刻度和符合相加修正。该程序在擦拭与气溶胶样品的实际测量中，可解决缺少标准效率刻度源问题，并使得符合相加修正计算简单实用。     

核数据理论计算程序处理系统的研制

研制了核数据理论计算程序处理系统，其中包括光学模型自动调参程序系统及ＵＮＦ程序输入数据的自动填充和输出结果的自动绘图程序系统。系统安装在中国核数据中心的ＭＩＣＲＯ－ＶＡＸ－Ⅱ上。简要地描述了系统的主要功能并给出使用实例。     

9 MeV行波电子直线加速器屏蔽设计与评价

采用Monte Carlo模拟计算程序EGS4,对海关大型集装箱检测系统用的9MeV行波电子直线加速器的机头和机房屏蔽进行了计算和辐射模拟分析,并分别与经验公式计算数据和实测值进行了比较。结果表明,在加速器机房内外主要参考点上,用EGS4程序计算得剂量当量率与实测结果符合较好,为EGS4的灵活应用提供了范例。     

铀溶液核临界装置物理试验研究进展

对于使用液体燃料的溶液堆，其结构比较复杂，现有计算程序的计算精度有限。MCNP程序对各种几何结构虽均有较强的适应性，但对于该程序的使用，尤其是数据库的选取需要通过校算基准例题和临界试验进一步确认。通过零功率物理试验，也可以为溶液堆的设计与运行提供一定的参考数据。     

热中子粉末衍射轮郭精修程序RIETAN运行环境的设置及使用方法

本文介绍了中子粉末衍射轮廓精修程序ＲＩＥＴＡＮ在微机上的运行环境和使用方法，并以亚稳态化合物Ｎｄ２Ｆｅ２３Ｂ３的精修过程为例，具体说明了输入参数文件和数据文件的编写方法和ＲＩＥＴＡＮ程序运行方法。     

空腔电离理论发展概述

简述了空腔电离理论的基本原理及其在辐射剂量学中的地位;介绍了Bragg—Gray理论、Spencer—Attix理论、Burch理论、Burlin模型、Kearsley模型、等效电子源理论等空腔电离理论及各种理论的适用范围;同时,介绍了用蒙特卡罗方法(EGS4程序)对不同理论及模型的计算结果。     

用EGS4分析9MeV行波直线加速器屏蔽系统

采用了国际上广泛使用的ＭｏｎｔｅＣａｒｌｏ模拟计算程序ＥＧＳ４对海关集装箱检测系统用加速器屏蔽系统进行了辐射模拟分析，与实际测试数据进行分析、比较，解释了特殊剂量分布的原因，提出了系统改进的参考方案。     

Incorporation of Landau-Pomeranchuk-Migdal effect and dielectric suppression effect in EGS5 code

We incorporated the Landau-Pomeranchuk-Migdal (LPM) effect and the dielectric suppression effect for bremsstrahlung and the LPM effect for pair production in the EGS5 code. To verify the validity of the EGS5 code with the LPM plus dielectric cross section for bremsstrahlung, we compared the bremsstrahlung spectrum calculated using this code with that obtained considering the LPM effect. The values calculated using the EGS5 code reproduced the experimental value of the LPM effect well.     

G（E）函数的蒙特卡罗模拟

G（E）函数法是通过软件来调整仪器能量响应的方法,可以直接从脉冲幅度谱得出剂量率。本文对一箱航测用NaI（Tl）探测器进行了地面实验和EGS4模拟计算,用地面实验测得的^137Cs和^60Co源能谱作为EGS4模拟谱的对照,EGS4计算了225个能谱,用这些能谱作为标准谱计算出了该探测器的G（E）函数。     

Gamma Ray Albedo Data Generated by the Invariant Embedding Method

The energy-and-angle dependent doubly differential γ-ray albedos for homogeneous semi-infinite medium have been calculated for water, ordinary concrete, soil, heavy concrete, iron, tin and lead. The procedure of calculation employed in the present paper is the invariant embedding method which is being developed to solve neutral particle transport problems in homogeneous one dimensional medium. The calculated γ-ray albedos are stored in the data base. One can easily obtain from this data base the various kinds of albedos (number, dose, energy, etc.) with a simple interface program. The accuracy of the calculated γ-ray albedos is ascertained by comparing with the Monte Carlo calculations (MCNP4A and EGS4).     

Monte Carlo simulation of Tabata’s electron backscattering experiments

Electron backscattering coefficients, η, obtained from several targets in the MeV range were calculated by using electron-photon Monte Carlo transport calculation codes, i.e., EGS5 and ITS 3.0. These calculated values were compared with those obtained from the electron backscattering experiment performed by Tabata using an ionization chamber [15]. We found that Tabata’s estimation of the multiplication factor of the ionization chamber, f, had a non-negligible error. Then, we calculated the ionization chamber output, I, which is a product of η and f. The ratios of I between the experimental and the calculated values were within 1.5 and 1.3 for the EGS5 code and the ITS 3.0 code, respectively. The ratios of η between the experimental and the calculated values were within 2.4 and 1.5 for the EGS5 code and the ITS 3.0 code, respectively. The differences between the experimental and the calculated values of I and η are large for low-Z targets (Be and C). Here, the ratios obtained by using the ITS 3.0 code are closer to unity than those obtained by using the EGS5 code. The reason of this is the fact that the calculated value obtained by using the ITS 3.0 code is underestimated for low-Z targets; this underestimation can, in turn, be attributed to the use of the default value of the number of steps in the electron transport algorithm in the ITS 3.0 code.     

利用EGS4计算固体电离室的吸收剂量与光子注量的转换因子

利用蒙特卡罗程序EGS4,针对几种特定的γ光子谱,计算了特定形状的LiF固体电离室中空腔辐射敏感元件LiFTLD的吸收剂量-光子注量转换因子;计算结果表明,该转换因子的不确定度主要取决于谱的不确定度;通过对单能光子源的蒙特卡罗计算结果与解析公式计算结果的比较,证明该方法是可靠的,是目前解决此类问题的有效而可行的途径。     

电子束氨法脱硫反应器内部剂量率的模拟计算

为验证MCNP用于计算电子束氨法脱硫反应器内部剂量率的可行性,本文借用“小孔成像”的原理对能量为0.57 MeV的电子加速器进行模拟,计算电子束氨法脱硫反应器内部的剂量率分布,并与文献中的实验测量值以及EGS4程序模拟值进行比较。结果表明:MCNP程序计算所得到的反应器内部剂量率分布在25.7%误差下的结果与EGS4程序模拟结果一致。