2024年5月23日发(作者:完颜思源)
第1章 MCNP概述
1.1 MCNP计算过程
MCNP(Monte Carlo N—Particle Transport code)是计算粒子输运过程的一套
蒙特卡罗模拟计算程序。这个程序需要用户通过输入文件给出计算模型。计算模型中需要
提供源的属性、感兴区内各种物体的属性、记录粒子信息的方法等。
例如,若想计算一个1MeV的X射线透过2cm铁的概率是多少,我们可以通过下面的
模型进行计算,如图1所示。
图 1 计算模型
在上面的计算模型中,感兴区是一个球的内部,其中包含X射线源、铁块和记录面,而
其他位置均为真空。由于当粒子被输运到感兴区外时,它将肯定不会再对记录结果产生贡
献,所以程序会自动停止这个粒子的输运过程,这也正是设定感兴区的原因。
源的属性主要包括位置、能量、出射方向、粒子种类等.图1的计算模型中,源的能量
为单能1MeV,方向为单向垂直于铁块的左表面,粒子种类为光子(Photon)。
感兴区内物体的属性包括几何尺寸、材料成分、密度等。图1中使用了一块铁块,它
的厚度为2cm,其他方向的尺寸对我们的计算结果没有影响,但要保证铁块完整地包含于
感兴区内。
记录方法有多种,其中包括通过某个面的特定种类粒子的个数。在图1中,我们可以
利用MCNP记录通过“记录面"的能量为1MeV的光子个数.
计算图1的模型时,MCNP会首先根据源的属性描述,抽样出一个起始粒子。图1中
的源为单能且单向的点源,所以每次抽样出的粒子都是能量、方向、种类相同的粒子.这个
粒子会沿着它的出射方向(垂直于铁块左表面)飞行,当它入射到铁块里时,会有一定的
概率发生康普顿散射、电子对效应和光电效应。发生三种反应的概率由MCNP的截面库中
的微观截面数据、输入文件中铁的密度以及抽样得到的随机数共同决定。若X射线发生了
康普顿散射,原来的X射线将被具有新属性的X射线取代,它将有不同的出射方向、能量。
MCNP会继续输运这个新产生的X射线直到它发生下一次反应或者飞出感兴区;X射线还
会有一定的概率不发生任何反应,直接透过铁块.
当基于一个起始粒子的输运过程结束后,MCNP会重复上述的过程。由于每次抽样的
得到的随机数会不同,所以每个起始粒子的输运过程都是独立而且不同的。多次的输运过
程会使得MCNP的记录结果越来越接近它的期望值,重复的次数越多,统计涨落就会越小。
尽管MCNP会记录所有通过记录面的光子。但由于光子一旦发生反应,它的能量就会
变化,所以若只观察能量为1MeV的光子,我们就可以排除那些经过多次散射而打到记录
面的光子了。
由于每次输运均是基于一个起始源粒子,所以最后记录得到的结果也都是相对于一个
起始粒子的概率。对于图1来说,记录结果就是1个X射线透过2cm的铁后打到记录面
的概率。
1.2 MCNP文件结构
2024年5月23日发(作者:完颜思源)
第1章 MCNP概述
1.1 MCNP计算过程
MCNP(Monte Carlo N—Particle Transport code)是计算粒子输运过程的一套
蒙特卡罗模拟计算程序。这个程序需要用户通过输入文件给出计算模型。计算模型中需要
提供源的属性、感兴区内各种物体的属性、记录粒子信息的方法等。
例如,若想计算一个1MeV的X射线透过2cm铁的概率是多少,我们可以通过下面的
模型进行计算,如图1所示。
图 1 计算模型
在上面的计算模型中,感兴区是一个球的内部,其中包含X射线源、铁块和记录面,而
其他位置均为真空。由于当粒子被输运到感兴区外时,它将肯定不会再对记录结果产生贡
献,所以程序会自动停止这个粒子的输运过程,这也正是设定感兴区的原因。
源的属性主要包括位置、能量、出射方向、粒子种类等.图1的计算模型中,源的能量
为单能1MeV,方向为单向垂直于铁块的左表面,粒子种类为光子(Photon)。
感兴区内物体的属性包括几何尺寸、材料成分、密度等。图1中使用了一块铁块,它
的厚度为2cm,其他方向的尺寸对我们的计算结果没有影响,但要保证铁块完整地包含于
感兴区内。
记录方法有多种,其中包括通过某个面的特定种类粒子的个数。在图1中,我们可以
利用MCNP记录通过“记录面"的能量为1MeV的光子个数.
计算图1的模型时,MCNP会首先根据源的属性描述,抽样出一个起始粒子。图1中
的源为单能且单向的点源,所以每次抽样出的粒子都是能量、方向、种类相同的粒子.这个
粒子会沿着它的出射方向(垂直于铁块左表面)飞行,当它入射到铁块里时,会有一定的
概率发生康普顿散射、电子对效应和光电效应。发生三种反应的概率由MCNP的截面库中
的微观截面数据、输入文件中铁的密度以及抽样得到的随机数共同决定。若X射线发生了
康普顿散射,原来的X射线将被具有新属性的X射线取代,它将有不同的出射方向、能量。
MCNP会继续输运这个新产生的X射线直到它发生下一次反应或者飞出感兴区;X射线还
会有一定的概率不发生任何反应,直接透过铁块.
当基于一个起始粒子的输运过程结束后,MCNP会重复上述的过程。由于每次抽样的
得到的随机数会不同,所以每个起始粒子的输运过程都是独立而且不同的。多次的输运过
程会使得MCNP的记录结果越来越接近它的期望值,重复的次数越多,统计涨落就会越小。
尽管MCNP会记录所有通过记录面的光子。但由于光子一旦发生反应,它的能量就会
变化,所以若只观察能量为1MeV的光子,我们就可以排除那些经过多次散射而打到记录
面的光子了。
由于每次输运均是基于一个起始源粒子,所以最后记录得到的结果也都是相对于一个
起始粒子的概率。对于图1来说,记录结果就是1个X射线透过2cm的铁后打到记录面
的概率。
1.2 MCNP文件结构