应用数学计算机技术.docVIP

应用数学与计算机技术 单色器屏蔽的蒙特卡罗模拟计算 刘保杰，许淑艳 本工作使用蒙特卡罗方法，借助于MCNP软件包计算了在旋转门关闭打开时单色器屏蔽体外侧各面的中子、光子泄率及剂量 为提高计算效率，在计算中利用几何系统的对称性，使用镜面反射使模型简化使用分裂及轮盘赌技巧，减小方差，提高计算的精确度。 光子电子的输运问题在粒子输运问题中占有相当比重在实验核物理和反应堆物理中大量存在。如探测器（NaIGe等）的探测效率及能量沉积谱等参数，通常采用实验的方法测量。测量时，需各种不同能量的光子源。由于低能光子源E0.1 MeV）和高能光子源（E2.0 MeV）大多寿命较短，获取困难，因而在大多情况下，使用外推方法获得，准确性较差且实验工作量大。因此需借助理论计算。 由于光子（或电子）入射到系统后，其输运过程是光子电子的输运过程，对于如此复杂的问题使用数值方法无法解决，使用蒙特卡罗方法是最行之有效的 对于两种粒子的输运问题，使用蒙特卡罗方法模拟时，采用“字典编辑式多分支”方法处理。由于该模拟过程非常复杂国内外相关软件如MCNPEGS4等均使用直接模拟方法模拟尚未使用减小方差技巧对于有些问题如小系统或吸收性较强的系统计算结果误差大效率低 本文从光子电子耦合输运问题的积分方程理论模型出发，将方差减小技巧（如限制碰撞法、加权法）应用到光子电子耦合输运的多分枝模拟方法中，并将方差减小技巧引入到EGS4软件包中。通过实例计算，大幅减小了方差，提高了效率。 许淑艳，丰树强，赵桂植 为生产123I建立124Xe气体靶靶室。30 MeV质子由回旋加速器射出，入射到12Xe气体靶上，通过发生（p, 2n）和（p, p, n）反应，产生快中子。本文使用蒙特卡罗方法，借助于MCNP及GEANT4软件包计算了气体靶靶室周围水泥墙的中子泄率为建立124Xe气体靶靶室的防护问题提供了理论依据。 气体靶靶室的结构及靶体的结构图1、2。 图1 气体靶靶室结构平面示意图 图2 气体靶靶体结构平面示意图 所求物理量为气体靶靶室周围三面水泥墙墙体的中子泄率，进而计算剂量。计算方法GEANT4计算 使用GEANT4软件包计算30 MeV质子与靶体中的介质发生（p, 2n）和（p, p, n）反应产生的快中子的能谱、角分布。 MCNP源抽样程序 编写MCNP软件包所需的用户源程序。包括质子由源出发，经过输运，在靶体内与Xe碰撞后发生（p, 2n）或（p, p, n）反应，抽样确定快中子的能量及运动方向。 填写MCNP输入卡片 填写MCNP软件包所需的输入文件 INP（其中所需的源程序由生成。） 计算 运行MCNP执行程序进行计算。 本文计算了所设计靶室不同厚度水泥墙的中子泄漏率。水泥墙厚度1.0 m下的计算结果略，不同能级下的中子穿透率图3由以上计算结果可知，将蒙特卡罗方法用于计算124Xe气体靶靶室周围水泥墙中子泄露率是可行的，它为建立124Xe气体靶靶室的防护问题提供了理论依据。 图3 不同能级下的中子穿透率 丰树强，刘保杰，程金生1 （1 中国疾病预防控制中心 辐射防护与核安全医学所） 基于60Co放射源的远程治疗机是常用的放疗设备。对治疗机的靶区做精确的剂量计算有助于确保放疗质量。现今放疗设备剂量分布的计算基本是以半经验解析方法为主，而半经验解析方法的适用范围是有限的，用蒙特卡罗方法可很好地克服以上缺陷，因此有必要对蒙特卡罗方法在放射医学领域中的应用进一步研究，建立一种更为准确和实用的剂量场分布的计算方法。此外，在放疗中剂量数据通常是通过扫描水箱直接测得，代表了射野内某些平面上有限点的数据，对于任意平面上任意点的剂量数据，受测量条件限制，费时费力，且难以保证测量结果的准确性。本工作使用蒙特卡罗方法，借助于EGS4软件包计算了不同射野下60Co远程治疗机在水膜体中的三维剂量场分布，获得了比实际测量更多、更全面的三维剂量数据。同时，为提高抽样效率，对源粒子的运动方向采用了解析方法确定。最后，将得到的计算结果与用RTS-200S自动扫描水箱测量的离轴剂量和中心轴深度剂量比较，二者符合很好。 60Co远程治疗机机头主要包括源和准直器两部分。其中60Co源为圆柱形结构，分为1.17 MeV、1.33 MeV两个能级。治疗机的准直器部分如图1所示：包括一级和二级准直器，控制射线在水膜体中的射野。本工作中采用立方体水膜，几何区域划分如图2所示。其中，准直器主要为钨层和空气两个区域。模型中涉及到的介质有水、空气、有机玻璃、钨，计算前需填写输入卡片USERINP.DAT，计算所需的截面数据由PEGS程序给出。 对于源抽样，为提高抽样效率，将采用确定源粒子入射方向的解析方法。具体做法是：根据给定的射野，确定准直器的位置以及源粒子直