计算图 1 的模型时,MCNP 会首先根据源的属性描述,抽样出一个起始粒子。
图 1 中的源为单能且单向的点源,所以每次抽样出的粒子都是能量、方向、种类
相同的粒子。这个粒子会沿着它的出射方向(垂直于铁块左表面)飞行,当它入
射到铁块里时,会有一定的概率发生康普顿散射、电子对效应和光电效应。发生
三种反应的概率由 MCNP 的截面库中的微观截面数据、输入文件中铁的密度以
及抽样得到的随机数共同决定。若 X 射线发生了康普顿散射,原来的 X 射线将
被具有新属性的 X 射线取代,它将有不同的出射方向、能量。MCNP 会继续输
运这个新产生的 X 射线直到它发生下一次反应或者飞出感兴区;X 射线还会有
一定的概率不发生任何反应,直接透过铁块。
当基于一个起始粒子的输运过程结束后,MCNP 会重复上述的过程。由于每
次抽样的得到的随机数会不同,所以每个起始粒子的输运过程都是独立而且不同
的。多次的输运过程使得 MCNP 的记录结果越来越接近它的期望值,重复的次
数越多,统计涨落就会越小。
尽管 MCNP 会记录所有通过记录面的光子。但由于光子一旦发生反应,它
的能量就会变化,所以若只观察能量为 1MeV 的光子,我们就可以排除那些经过
多次散射而打到记录面的光子了。
由于每次输运均是基于一个起始源粒子,所以最后记录得到的结果也都是相
对于一个起始粒子的概率。对于图 1 来说,记录结果就是 1 个 X 射线透过 2cm
的铁后打到记录面的概率。
1.2 MCNP 文件结构
MCNP 的输入文件名不可以多于七个字符,所以不建议对输入文件增加扩展
名;MCNP 的输出文件的扩展名为.o;MCNP 内部计算过程的记录文件的扩展名
为.r。若一次计算结果就已经达到满意的统计结果,那么.r 文件就没有什么价值
保存下来;但是若一次计算结果达不到满意的统计可信度,那么我们可以在已得
到的计算结果的基础上继续模拟输运过程,来减小统计涨落。而为了这么做,
MCNP 需要上一次计算结果中使用的.r 文件。为了这个目的,有时候.r 文件还是
值得保留的。在基于上一次的计算结果继续模拟计算时,MCNP 将不再读入输入
文件,而单单读入.r 文件的信息。所以我们想继续计算时,需要保证.r 文件就是
我们所设想的模型的前一次模拟计算中得到的.r 文件。除了上面三种常用的文件
类型,MCNP 还有两种不长用的输出文件。这两种文件只有用户通过输入文件要
求 MCNP 输出这些文件时才会得到。其中一个为.w 文件,这个文件是我们想通
过 MCNP 输出 surface source 的时候才会产生;另外一个为.p 文件,这个文件时
我们想通过 MCNP 追踪粒子的输运过程的时候才会得到。
上面的几种文件是用户常用到的文件类型,下面将介绍 MCNP 程序内部使
用的若干的文件。
1. MCNP5.exe:MCNP 的应用程序,这个程序是在 dos 下运行的程序。所
以用户不要试图双击这个程序来运行它。
2. vised.exe:这个程序可以让用户观察 MCNP 读入的输入文件中的计算模
型。我们可以利用这个程序观察我们所设想的计算模型与实际书写出来
的计算模型是否有偏差。此外,这个程序还可以方便我们找出 MCNP
输入文件中潜在的一些错误,例如:两个物体有相交的部分等。
3. 目录 Documents: 这个文件夹内包含 LANL(Los Alamos National
评论0