常用蒙特卡罗程序介绍

(2)
信息块的卡片放在INP文件中标题卡之前。
信息块给出了MCNP的一些运行信息，信息块上各部
分的意思和运行行信息是一样的，当运行行信息与信
息块中所指定的信息相矛盾时，则忽略信息块中相应
的信息，而以运行行信息为准。
(3)
信息块是可选的，信息块的第一张卡片，
必须在第 1～8 列写上“MESSAGE:”，从第一张卡片
定义组成栅元的曲面信息。
其它数据，包括问题类型、源描述、材料描 述、计数描述，问题截断条件等。 选择项
10 可编辑ppt
(3) 接续运行的输入文件
信息块 空行分隔符
选择项
CONTINUE 写在第 1～8 列
数据卡
… 空行分隔符
只允许部分数据卡。（FQ，DD，NPS，CTME，
IDUM，RDUM，PRDMP，LOST，DBCN，PRINT，
7 可编辑ppt
MCNP输入文件中物理量的单位
长度
厘米
能量
MeV
时间
10-8 秒
温度
MeV（kT）
原子密度
1024 个原子 / 厘米3
质量密度
克 / 厘米3
截面
10-24 厘米2
原子量
中子质量的1.008664967倍
阿伏加德罗常数 6.023×1023
8 可编辑ppt
1) 输入文件的基本形式 (1) 信息块
1) MORSE程序 较早开发的通用蒙特卡罗程序，可以解决
中子、光子、中子－光子的联合输运问题。采 用组合几何结构，使用群截面数据，程序中包 括了几种重要抽样技巧，如俄国轮盘赌和分裂 技巧，指数变换技巧，统计估计技巧和能量偏 移抽样等。程序提供用户程序，用户可根据需 要编写源分布以及记录程序。
4 可编辑ppt
(4) 卡片格式 INP输入文件的每一行（称之为一张卡片）
都限于使用第 1～80 列并构成卡片映象。大部 分输入卡片按行填写；然而，对数据卡允许按 列填写。 $ 符号为它所在那行数据的结束符， 在 $ 符号后面的内容作为注释，它可从 $ 符号 后面的任一列开始。
标题卡只占一行，整行都可填入用户需要 的信息，也可以是空行。但要注意在其它地方 使用空行是作为结束符或者分隔符。
2 可编辑ppt
1. 蒙特卡罗方法应用软件的特点
通用蒙特卡罗程序通常具有以下特点： 1) 具有灵活的几何处理能力 2) 参数通用化，使用方便 3) 元素和介质材料数据齐全 4) 能量范围广，功能强，输出量灵活全面 5) 含有简单可靠又能普遍适用的抽样技巧 6) 具有较强的绘图功能
3 可编辑ppt
2. 常用的通用蒙特卡罗程序简介
6 可编辑ppt
3. MCNP程序输入的描述
MCNP 的 输 入 包 括 几 个 文 件 ， 但 主 要 的 一个是由用户编写的INP文件，该文件包括描 述问题所必须的全部输入信息。文件采用卡片 结构，每行代表一张卡片，文件由一系列卡片 组成，对于任一特定的问题，只需用到INP全 部输入卡片的一小部分。
13 可编辑ppt
Fra Baidu bibliotek
为书写方便，可以使用四项书写功能：
A. nR功能，表示将它前面的数据重复n次。
例如：2 4R 等同于 2 2 2 2 2
B. nI功能，表示在与其前后相邻的两个数之间，插入n 个线性插值点。对于 X nI Y 的结构，如果X和Y是整 数，且X－Y刚好是n+1的整倍数，则产生标准的整数 插值，否则产生实数插值，但Y值直接存储。
输入文件中，在标题卡之后及最后的空行 结束卡之前的任何地方都可插入注释卡。注释 卡必须是字母“C”写在 1～5 列中的任意位置， 且至少用一个空格隔开后面的注释内容。
12 可编辑ppt
a) 行输入格式 栅元卡、曲面卡和数据卡的书写格式是相
同的。必须从 1～5 列开始填写这些卡片相应 的名字(或编号)和粒子标识符，后面填写用空 格分隔的数据项。如果 1～5 列为空，则表示 它是前一张卡片的继续卡。如果在一行的末尾 有一个用空格隔开的符号“&”，则表示下一 行是该行的继续卡，数据可填写在 1～80 列。 一个数据项必须在一张卡片上写完，不得跨到 下一张卡片上。完全空白的一行则为两组卡片 的分隔符。 对任何给定的带有粒子标识符的 类型卡只能有一张。需要整数的数据项必须填 写整数，其它数据可填写为整数或浮点数以及 MCNP能读的数据。
KCODE，MPLOT，ZA，ZB，和 ZC）
其它
选择项
接续运行必须在运行行信息或信息块中给出C项选择， 即Cm，表示从RUNTPE文件中读出第m次转储的内容 接着运算，如果m未指定，则读最后一次转储的数据。
如果不需要改变内容，则不需要接续输入文件，仅需
运行RUNTPE以及在运行行加上C选择。
11 可编辑ppt
第八章 蒙特卡罗方法应用程序介绍
1. 蒙特卡罗方法应用软件的特点 2. 常用的通用蒙特卡罗程序简介 3. MCNP程序输入的描述
1 可编辑ppt
第八章 蒙特卡罗方法应用程序介绍
建立完善的通用蒙特卡罗程序可以 避免大量的重复性工作，并且可以在程 序的基础上，开展对于蒙特卡罗方法技 巧的研究以及对于计算结果的改进和修 正的研究，而这些研究成果反过来又可 以进一步完善蒙特卡罗程序。
5 可编辑ppt
3) MCNP程序 MCNP是美国Los Alamos国家实验室开发
的大型多功能通用蒙特卡罗程序，可以计算中 子、光子和电子的联合输运问题以及临界问题， 中子能量范围从10-11MeV至20MeV，光子和电 子的能量范围从1KeV至1000MeV 。程序采用 独特的曲面组合几何结构，使用点截面数据， 程序通用性较强，与其它程序相比，MCNP程 序中的减方差技巧是比较多而全的。
的第 9～80 列到后续卡片的第 1～80 列都可填写运行
信息。在标题卡之前用一个空行分隔符结束信息块。
9 可编辑ppt
(2) 初始运行的输入文件
信息块 空行分隔符
标题卡
栅元卡 … 空行分隔符 曲面卡 … 空行分隔符 数据卡 … 空行分隔符
其它
选择项 仅一行，占用第 1～80 列。作为输出标题。 定义构成整个系统的各个基本介质单元以及 相应的物理信息。
2) EGS程序 EGS是Electron-Gamma Shower 的缩写，
它是一个用蒙特卡罗方法模拟在任意几何中， 能量从几个KeV到几个TeV的电子-光子簇射过 程 的 通 用 程 序 包 。 由 美 国 Stanford Linear Accelerator Center提供。EGS于1979年第一次 公开发表，提供使用。EGS4是1986年发表的 EGS程序的最新版本。