2012蒙特卡罗模拟02课程.pptVIP

* * 运行的方法： 1）假设输入卡文件名是mcin.inp，那么运行命令可以是： mcnp5 i=mcin.inp o=tesa.out ru=mcrntpe 或 mcnp5 i=mcin.inp o=tesa.out 如果在当前文件夹中存在与mcnp将要产生的文件名相同，则mcnp自动按字母顺序将文件名最后一个字符改变为字母表中下一字符。 2）如果想让在一次运行中所产生的文件，有相似的文件名，使用如下方式： mcnp5 name=job1 运行后产生的OUTP文件将为JOB1O，RUNTPE为JOB1R。如果文件已经存在，MCNP将不会覆盖它们，而是给出一条消息并中断运行。 2.4 MCNP5程序安装、运行（4） * * 看图的方法： 1）如果要看几何绘图，先删除工作目录下刚生成的job1o和job1r两个文件，再输入下面命令并回车： mcnp5 ip n=job1 如果能绘图的话会弹出plot提示符，在此例如你可以输入ex=100，就可以看到自己的几何描述了。 2）用“VisEd”工具软件，装载mcin.inp文件即可。 2.4 MCNP5程序安装、运行（5） * * 输入文件的基本形式： INP文件有初始运行和接续运行两种形式，用户须在INP中描述问题的几何、材料、源、记数和其它要求。INP由一些被空行分隔的输入块组成，主要的输入块是信息块、标题和栅元块、曲面块和数据块等。输入块又由一些被称为卡的输入行组成。初始运行输入文件的格式如右： 2.5 MCNP输入文件（1） 信息卡 } 可选 空行分隔符 标题卡 栅元卡 . . 空行分隔符 曲面卡 . . 空行分隔符 数据卡 . . 空行结束符 其它 ? * * 信息卡的1-8列应填写MESSAGE:，后面跟着用空格分隔的参数项。可用A=B参数项更改输出文件名，如OUTP = MYOUT。信息块是可选的。 在信息块之后的第一行是问题的标题卡，它仅限于一行，占用1-80列，可以是任何信息，将作为OUTP文件中各个输出表的标题被复制。 用户在栅元块和曲面块中描述问题的几何。栅元由栅元卡描述。空间必须由彼此相邻的栅元填满，栅元之间不能重叠，也不能出现无栅元的空区，否则会出现错误。构建栅元的曲面由曲面卡定义，曲面卡在曲面块中给出。 曲面块之后是数据块，在数据块中用户描述源、记数方式、材料等。 无论数据块后是否有空行结束符，MCNP都能运行。用户可以把希望保留的一些附加信息写在数据块的空行之后，MCNP会将它们复制到OUTP文件末尾。 2.5 MCNP输入文件（2） 信息卡 } 可选 空行分隔符 标题卡 栅元卡 . . 空行分隔符 曲面卡 . . 空行分隔符 数据卡 . . 空行结束符 其它 ? 二、MCNP程序入门 2.1 蒙特卡罗方法应用程序介绍 2.2 MCNP简介 2.3 MCNP程序的组成及特点 2.4 MCNP5程序安装、运行 2.5 MCNP输入文件 * * * * 建立完善的通用蒙特卡罗程序可以避免大量的重复性工作，并且可以在程序的基础上，开展对于蒙特卡罗方法技巧的研究以及对于计算结果的改进和修正的研究，而这些研究成果反过来又可以进一步完善蒙特卡罗程序。 2.1 蒙特卡罗方法应用程序介绍(1) * * 通用蒙特卡罗程序通常具有以下特点： 具有灵活的几何处理能力 参数通用化，使用方便 元素和介质材料数据齐全 能量范围广，功能强，输出量灵活全面 含有简单可靠又能普遍适用的抽样技巧 具有较强的绘图功能 2.1 蒙特卡罗方法应用程序介绍(2) * * 1)MORSE程序 较早开发的通用蒙特卡罗程序，可以解决中子、光子、中子－光子的联合输运问题。采用组合几何结构，使用群截面数据，程序中包括了几种重要抽样技巧，如俄国轮盘赌和分裂技巧，指数变换技巧，统计估计技巧和能量偏移抽样等。程序提供用户程序，用户可根据需要编写源分布以及记录程序。 2.1 蒙特卡罗方法应用程序介绍(3) * * 2)EGS程序 EGS是Electron-Gamma Shower 的缩写，它是一个用蒙特卡罗方法模拟在任意几何中，能量从几个KeV到几个TeV的电子-光子簇射过程的通用程序包。由美国Stanford Linear Accelerator Center提供。EGS于1979年第一次公开发表，提供使用。EGS4是1986年发表的EGS程序的最新版本。 2.1 蒙特卡罗方法应用程序介绍(4) * * 3)MCNP程序 MCNP是美国Los Alamos国家实验室开发的大型多功能通用蒙特卡罗程序，可以计算中子、光子和电子的联合输运问题以及临界问题，中子能量范围从10-1