MCNP使用教程（整理）.pptx

第 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 的模型时，MCNP 会首先根据源的属性描述，抽样出一个起始粒子。 图 1 中的源为单能且单向的点源，所以每次抽样出的粒子都是能量、方向、种类 相同的粒子。这个粒子会沿着它的出射方向（垂直于铁块左表面）飞行，当它入 射到铁块里时，会有一定的概率发生康普顿散射、电子对效应和光电效应。发生 三种反应的概率由 MCNP 的截面库中的微观截面数据、输入文件中铁的密度以及 抽样得到的随机数共同决定。若 X 射线发生了康普顿散射，原来的 X 射线将被具 有新属性的 X 射线取代，它将有不同的出射方向、能量。MCNP 会继续输运这个 新产生的 X 射线直到它发生下一次反应或者飞出感兴区；X 射线还会有一定的概 率不发生任何反应，直接透过铁块。 当基于一个起始粒子的输运过程结束后，MCNP 会重复上述的过程。由于每 次抽??的得到的随机数会不同，所以每个起始粒子的输运过程都是独立而且不同 的。多次的输运过程会使得 MCNP 的记录结果越来越接近它的期望值，重复的次 数越多，统计涨落就会越小。 尽管 MCNP 会记录所有通过记录面的光子。但由于光子一旦发生反应，它的 能量就会变化，所以若只观察能量为 1MeV 的光子，我们就可以排除那些经过多 次散射而打到记录面的光子了。 由于每次输运均是基于一个起始源粒子，所以最后记录得到的结果也都是相 对于一个起始粒子的概率。对于图 1 来说，记录结果就是 1 个 X 射线透过 2cm 的铁后打到记录面的概率。 1.2 MCNP 文件结构 MCNP 的输入文件名不可以多于七个字符，所以不建议对输入文件增加扩展 名；MCNP 的输出文件的扩展名为.o；MCNP 内部计算过程的记录文件的扩展名为.;3;4;5;;7;输入文件中有两种注释方法，助记符分别为 c 和$符号。在第一行写 c 就是 把这一行都注释掉；在某一行的某列上写$就是注释该行该列后的内容。但是即 使是注释，请用户也不要超过每一列 80 个字符（包括空格）的要求。 还有三个常用的助记符在这里介绍个用户，分别为 r，i 和 j。r 代表的是 r epeat，例如 1 3r 就等于写了 1 1 1 1；i 代表线性插入点，例如 1 2i 4 就等 于写了 1 2 3 4；j 代表的是 jump，表示该位置使用默认值，例如 1 j 2 就等于 写了 1 默认值 2。 由于 MCNP 对输入文件的规范要求苛刻，所以请用户谨记这样的规范。 2.1 7MeV 轫致辐射 X 射线算例 下面将通过书写计算 7MeV 的轫致辐射 X 射线的角分布、能谱的方法介绍 MC NP 输入文件的书写方法。 在书写 MCNP 输入文件前，请用户首先想好要计算的模型，包括坐标系以及 各个物体在这个坐标系内的相对位置。我的计算模型如下图所示。 图 2 计算模型；在图 2 的计算模型中感兴区是半径为 100cm 的球内；为了书 写上面的计算模型，我们首先写 surface；1，源所在面的书写方法：1PX-2；直 于 X 轴的面；-2 代表的是该面在 X 轴上的轴距；2，同理 Au 的左底面的书写方 法：2PX0；3，同理 Au 的右底面的书写方法：3PX1；4，Au 的侧面的书写方法： 4CX1；到类似地，还会有 CY 和 CZ 的助记符）；5，圈 ;9;10;11;1