MCNP使用教程（整理）（一）.pdf

第1 章 MCNP 概述 1.1 MCNP 计算过程 MCNP （Monte Carlo N-Particle Transport code）是计算粒子输运过程的 一套蒙特卡罗模拟计算程序。这个程序需要用户通过输入文件给出计算模型。计 算模型中需要提供源的属性、感兴区内各种物体的属性、记录粒子信息的方法等。 例如，若想计算一个1MeV 的X 射线透过2cm 铁的概率是多少，我们可以通 过下面的模型进行计算，如图1 所示。 图 1 计算模型 在上面的计算模型中，感兴区是一个球的内部，其中包含X 射线源、铁块和 记录面，而其他位置均为真空。由于当粒子被输运到感兴区外时，它将肯定不会 再对记录结果产生贡献，所以程序会自动停止这个粒子的输运过程，这也正是设 定感兴区的原因。 源的属性主要包括位置、能量、出射方向、粒子种类等。图1 的计算模型中， 源的能量为单能1MeV，方向为单向垂直于铁块的左表面，粒子种类为光子（Pho ton）。 感兴区内物体的属性包括几何尺寸、材料成分、密度等。图1 中使用了一块 铁块，它的厚度为2cm，其他方向的尺寸对我们的计算结果没有影响，但要保证 铁块完整地包含于感兴区内。 记录方法有多种，其中包括通过某个面的特定种类粒子的个数。在图1 中， 我们可以利用MCNP 记录通过 “记录面”的能量为1MeV 的光子个数。 1 计算图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 内部计算过程的记录文件的扩展名为. 2 r。若一次计算结果就已经达到满意的统计结果，那么.r 文件就没有什么价值保 存下来；但是若一次计算结果达不到满意的统计可信度，那么我们可以在已得到 的计算结果的基础上继续模拟输运过程，来减小统计涨落。而为了这么做，MCN P 需要上一次计算结果中使用的.r 文件。为了这个目的，有时候.r 文件还是值 得保留的。在基于上一次的计算结果继续模拟计算时，MCNP 将不再读入输入文 件，而单单读入.r 文件的信息。所以我们想继续计算时，需要保证.r 文件就是 我们所设想的模型的前一次模拟计算中得到的.r 文件。除了上面三种常用的文 件类型，MCNP 还有两种不长用的输出文件。这两种文件只有用户通过输入文件 要求MCNP 输出这些文件时才会得到。其中一个为.w 文件，这个文件是我们想通 过MCNP 输出surface source 的时候才会产生；另外一个为.p 文件，这个文件 时我们想通过MCNP 追踪粒子的输运过程的时候才会得到。 上面的几种文件是用户常用到的文件类型，下面将介绍MCNP 程序内部使用 的若干的文件。 1. MCNP5.exe：MCNP 的应用程序，这个程序是在dos 下运行的