TRIPOLI_4_gamma射线防护模拟要点.docxVIP

TRIPOLI-4?Monte Carlo Transport Code Practical Work. Dimensioning of a shield against gamma rays 2014 Page  PAGE \* MERGEFORMAT 11/ NUMPAGES \* MERGEFORMAT 11 TRIPOLI-4?Monte Carlo Transport Code Practical Work Dimensioning of a shield against gamma rays 周靖丰一．引言 在核反应堆中，很多类型的核反应都会释放出高能光子，比如中子的辐射吸收，一些激发态原子的退激等，这些光子能量高，并且释放频率高，倘若人体接触，会造成难以估计的伤害，所以对于高能伽马射线的屏蔽是非常重要的基础研究。在本次实验中，我们模拟一个伽马点源和一道20cm厚的屏障，二者相距约1m，这与实际生活中实验室对伽马射线防护很接近。我们使用的工具是TRIPOLI，这是由法国原子能协会（CEA）所开发的强大软件，它利用蒙特卡洛方法模拟光子运输过程，得到我们需要的通量密度（flux）和吸收剂量当量率（dose equivalent rate）。 二．不同物质对gamma射线的阻止能力 本次实验的装置如下 本次Tripoli·数值模拟的几何模型如上图所示，所有的实验物件由两个放置在横轴坐标x=10处的放射源，一个厚度为20cm的屏蔽体，以及一个探测器。屏蔽体的中点放置于x=120处，探测器位于x=220处，也就是我们所取的计算点的位置。屏蔽体的材料有以下几个选择： 材料密度(g/cm3)成分空气1.205×10-3714N=80%,816O=20%水111H=11.1%,816O=88.9%混凝土3.511H=0.001,C=0.003,816O=0.513, 27Al=0.02,28Si=0.38,40Ca=0.078,56Fe=0.005钢7.8656Fe铅11.34208Pb 本次模拟选取的gamma放射源有两个，强度均为 1×1010，且均为单能放射源，放出的gamma光子的能量分别为：1.17Mev和1.33Mev。 程序输出的物理量有6个：用三种不同的方法计算的gamma光子的通量与等效剂量率。这三种方法分别为：track、collision以及point kernel。Track方法与collision方法需要用户在其制定的位置定义出一个空间几何体，通过统计在这个定义的空间体中的一些计数来估计该处的gamma光子通量以及等效剂量率。Track方法统计的主要是几何体中所有经过的gamma光子留下的轨迹的长度，而collision方法则是通过统计几何体中粒子的碰撞次数。Point kernel方法不需要定义空间几何体，它是通过计算所有放射出的gamma粒子对该点的通量贡献来估计总的通量的。在接下来的模拟中我们可以清晰的看到在不同情况下各种方法的优劣，这就是本次数值模拟实验的主要目的。 我们选用不同的屏蔽体材料，分别为水，混凝土，钢还有铅，选取size=20000，所得计算结果如下： 表一：不同屏蔽体材料作用下的gamma通量空气 (batch=500)水 (batch=500)混凝土 (batch=500)钢 (batch=1000)铅 (batch=1000)Estimatorsx (cm)平均值方差平均值方差平均值方差平均值方差平均值方差Collision2204.90E+0437.588%3.00E+0437.815%1.27E+0452.312%0.00E+000.000%0.00E+000.000%Track2203.68E+041.499%3.47E+041.388%5.88E+033.526%9.58E+0118.578%0.00E+000.000%Point kernel201.59E+070.096%1.60E+079.481%1.59E+070.100%1.59E+070.066%1.59E+070.066%401.77E+060.099%1.80E+060.097%1.79E+060.095%1.78E+060.067%1.78E+060.083%509.95E+050.095%1.02E+060.089%1.02E+060.095%1.01E+060.071%9.98E+050.070%606.37E+050.103%6.