MCNP程序简介—张毅综述.pptVIP

MCNP 程序  MCNP程序介绍 张毅 兰州大学核科学与技术学院 主 要 内 容 蒙特卡罗（Monte Carlo）方法 MCNP运行过程 栅元卡 MCNP输入文件 曲面卡 数据卡 程序演示 蒙特卡罗（Monte Carlo）方法，又称随机抽样或统计试验方法属于计算数学的一个分支 MCNP 是一个通用的Monte Carlo 粒子输运程序。可用于中子、光子、电子，或耦合中子、光子、电子输运。是由美国洛斯阿拉莫斯国家实验室(LosAlamos National Laboratory)在蒙特卡罗方法的基础上开发的。 Monte Carlo 运行过程： MCNP运行 MCNP程序输入的描述 MCNP的输入包括几个文件，但主要的一个是由用户编写的INP文件，该文件包括描述问题所必须的全部输入信息。文件采用卡片结构，每行代表一张卡片，文件由一系列卡片组成，对于任一特定的问题，只需用到INP全部输入卡片的一小部分。 具体描述： (1) 空间几何的描述说明； (2) 几何体的使用材料描述； (3) 光子、中子以及电子这3 种粒子源的位置和特性说明； (4) 必要的标记卡及回答卡的类型； (5) 必需的冗余量消除技术以提高计算效率 MCNP输入文件 输入文件结构 MCNP输入文件结构 信息块 空行分隔符 选择项 标题卡 仅一行，占用第 1～80 列。作为输出标题。 栅元卡 … 空行分隔符 定义构成整个系统的各个基本介质单元以及相应的物理信息。 曲面卡 … 空行分隔符 定义组成栅元的曲面信息。 数据卡 … 空行分隔符 其它数据，包括问题类型、源描述、材料描述、计数描述，问题截断条件等。 其它 选择项 数据卡 一 、栅元卡 格式： j m d geom param 或 j LIKE n BUT list j：栅元序号，可在1到9999之间 m：材料序号， 从材料卡中指定材料m，若m=0表示真空 d: 材料密度，正直时表示原子密度（1024个原子/cm3)，当为负值时，表示质量密度（g/cm3) geom:栅元几何说明，由几何曲面号组成 param：栅元参数说明 n：其他栅元 list：栅元j与n的不同之处 例子：2 3 -3.7 -1 2 IMP：N=2 3 LIKE 2 BUT TRCL=1 IMP：N=1 MCNP栅元卡 二、曲面卡 格式： j n a list j：曲面号，1≤j≤99999， 反射面加*号，白边界加+号 n: 一般为缺省或零，表示没有坐标变换，若n≥0特定的TRn卡卡号， 若n≤0曲面j是曲面n的周期变化 a：方程记助名 list：按照规定次序方程描述的数据项 例：1 py 0 2 py 8 1 -2 -3 3 cy 2 或者： 4 RCC 0 0 0 0 1 0 2 MCNP曲面卡 MCNP数据卡 三、数据卡 在信息卡、栅元描述卡和曲面描述卡之后输入的是数据卡，数据卡可分为10类： 1、问题类型 2、几何卡 3、减方差 4、源描述 5、计数描述 6、材料及截面描述 7、能量及热处理 8、问题截断条件 9、用户数据数组 10、外围卡 MCNP数据卡 问题类型（MODE）卡 格式： MODE x1 … xi xi ＝ N，中子输运。 P，光子输运。 E，电子输运。 Mode N neutron transport only (default) N P neutron and neutron-induced photon transport P