单击此处编辑母版样式,单击此处编辑幻灯片母版样式,第二层,第三层,第四层,第五层,*,*,*,第八章,蒙特卡罗方法应用程序介绍,蒙特卡罗方法应用软件的特点,常用的通用蒙特卡罗程序简介,MCNP,程序输入的描述,第八章,蒙特卡罗方法应用程序介绍,建立完善的通用蒙特卡罗程序可以避免大量的重复性工作，并且可以在程序的基础上，开展对于蒙特卡罗方法技巧的研究以及对于计算结果的改进和修正的研究，而这些研究成果反过来又可以进一步完善蒙特卡罗程序蒙特卡罗方法应用软件的特点,通用蒙特卡罗程序通常具有以下特点：,具有灵活的几何处理能力,参数通用化，使用方便,元素和介质材料数据齐全,能量范围广，功能强，输出量灵活全面,含有简单可靠又能普遍适用的抽样技巧,具有较强的绘图功能,常用的通用蒙特卡罗程序简介,MORSE,程序,较早开发的通用蒙特卡罗程序，可以解决中子、光子、中子光子的联合输运问题采用组合几何结构，使用群截面数据，程序中包括了几种重要抽样技巧，如俄国轮盘赌和分裂技巧，指数变换技巧，统计估计技巧和能量偏移抽样等程序提供用户程序，用户可根据需要编写源分布以及记录程序EGS,程序,EGS,是,Electron-Gamma Shower,的缩写，它是一个用,蒙特卡罗方法,模拟在任意几何中，能量从几个,KeV,到几个,TeV,的电子-光子簇射过程的,通用程序包,。
由美国,Stanford Linear Accelerator Center,提供EGS,于1979年第一次公开发表，提供使用EGS4,是1986年发表的,EGS,程序的最新版本MCNP,程序,MCNP,是美国,Los Alamos,国家实验室开发的大型多功能通用蒙特卡罗程序，可以计算中子、光子和电子的联合输运问题以及临界问题，中子能量范围从10,-11,MeV,至20,MeV,，,光子和电子的能量范围从1,KeV,至1000,MeV,程序采用独特的曲面组合几何结构，使用点截面数据，程序通用性较强，与其它程序相比，,MCNP,程序中的减方差技巧是比较多而全的MCNP,程序输入的描述,MCNP的输入包括几个文件，但主要的一个是由用户编写的INP文件，该文件包括描述问题所必须的全部输入信息文件采用卡片结构，每行代表一张卡片，文件由一系列卡片组成，对于任一特定的问题，只需用到INP全部输入卡片的一小部分,MCNP输入文件中物理量的单位,长度,厘米,能量,MeV,时间,10,-8,秒,温度,MeV（kT）,原子密度,10,24,个原子/厘米,3,质量密度,克/厘米,3,截面,10,-24,厘米,2,原子量,中子质量的1.008664967倍,阿伏加德罗常数,6.02310,23,输入文件的基本形式,信息块,信息块的卡片放在,INP,文件中标题卡之前。
信息块给出了,MCNP,的一些运行信息，信息块上各部分的意思和运行行信息是一样的，当运行行信息与信息块中所指定的信息相矛盾时，则忽略信息块中相应的信息，而以运行行信息为准信息块是可选的，信息块的第一张卡片，必须在第,18,列写上“,MESSAGE:,”，,从第一张卡片的第,980,列到后续卡片的第,180,列都可填写运行信息在标题卡之前用一个空行分隔符结束信息块初始运行的输入文件,信息块,空行分隔符,选择项,标题卡,仅一行，占用第,180,列作为输出标题栅元卡,空行分隔符,定义构成整个系统的各个基本介质单元以及相应的物理信息曲面卡,空行分隔符,定义组成栅元的曲面信息数据卡,空行分隔符,其它数据，包括问题类型、源描述、材料描述、计数描述，问题截断条件等其它,选择项,接续运行的输入文件,接续运行必须在运行行信息或信息块中给出,C,项选择，即,Cm,，,表示从,RUNTPE,文件中读出第,m,次转储的内容接着运算，如果,m,未指定，则读最后一次转储的数据如果不需要改变内容，则不需要接续输入文件，仅需运行,RUNTPE,以及在运行行加上,C,选择信息块,空行分隔符,选择项,CONTINUE,写在第,18,列,数据卡,空行分隔符,只允许部分数据卡。
FQ，DD，NPS，CTME，IDUM，RDUM，PRDMP，LOST，DBCN，PRINT，KCODE，MPLOT，ZA，ZB，,和,ZC,）,其它,选择项,卡片格式,INP,输入文件的每一行（称之为一张卡片）都限于使用第,180,列并构成卡片映象大部分输入卡片按行填写；然而，对数据卡允许按列填写符号为它所在那行数据的结束符，在$符号后面的内容作为注释，它可从$符号后面的任一列开始标题卡只占一行，整行都可填入用户需要的信息，也可以是空行但要注意在其它地方使用空行是作为结束符或者分隔符输入文件中，在标题卡之后及最后的空行结束卡之前的任何地方都可插入注释卡注释卡必须是字母“,C,”,写在,15,列中的任意位置,，,且至少用一个空格隔开后面的注释内容行输入格式,栅元卡、曲面卡和数据卡的书写格式是相同的必须从,15,列开始填写这些卡片相应的名字(或编号)和粒子标识符，后面填写用空格分隔的数据项如果,15,列为空，则表示它是前一张卡片的继续卡如果在一行的末尾有一个用空格隔开的符号“,&,”，则表示下一行是该行的继续卡，数据可填写在,180,列一个数据项必须在一张卡片上写完，不得跨到下一张卡片上。
完全空白的一行则为两组卡片的分隔符对任何给定的带有粒子标识符的类型卡只能有一张需要整数的数据项必须填写整数，其它数据可填写为整数或浮点数以及,MCNP,能读的数据为书写方便，可以使用四项书写功能：,n,R,功能，表示将它前面的数据重复,n,次,例如：2 4,R,等同于 2 2 2 2 2,n,I,功能，表示在与其前后相邻的两个数之间，插入,n,个线性插值点对于,X,n,I,Y,的结构，如果,X,和,Y,是整数，且,X,Y,刚好是,n,+1,的整倍数，则产生标准的整数插值，否则产生实数插值，但,Y,值直接存储例如：1.5 2,I 3.0 1.5 2.0 2.5 3,2.0,可能不精确,而 1 4,I 6 1 2 3 4 5 6,都是精确定整数,X,M,功能，它表示的数值为前面的数据乘上,X,例如：1 1 2,M 2M 4M 2M 1 1 2 4 16 32,n,J,功能，表示其后,n,个数据项使用缺省值例如：,DD .1（,缺省值,）,1000,DD J 1000,如果,n,R,、,n,I,、,及,n,J,项,中缺省,n,，,则假设,n,1列输入格式,列输入块的格式：,S,i,必须是,MCNP,卡片名字，它们必须全部是栅元参数、或者全部是曲面参数、或者全部是其它参数。
15,列,672,列,S,1,S,2,S,m,K,1,D,11,D,12,D,1m,K,2,D,21,D,22,D,2m,K,n,D,n1,D,n2,D,nm,粒子标识符,几个输入卡片都需要粒子标识符以区别中子、光子和电子的输入数据这些卡片是：,IMP、EXT、FCL、WWN、WWE、WWP、WWGE、DXT、DXC、F、F5X、F5Y、F5Z、PHYS、ELPT、ESPLT、CUT,和,PERT粒子标识符由上述卡片名字后面的冒号、字母,N,、,P,或,E,组成例如：中子重要性卡为,IMP:N,光子重要性卡为,IMP:P,缺省值,MCNP,的许多输入参数都有缺省值，因此用户不需要每次都给出各个输入参量的值当缺省值符合用户要求时，便可不在输入文件中指定当省略某张输入卡时，则该卡上的全部参数均使用缺省值如果只想改变一张卡上的某一个特定参量时，则它前面的参量仍需指明，或者用,nJ,方式跳过前面那些使用缺省值的参量例如：光子截断卡,CUT:P 3J -.10,表示前3个参量使用缺省值，只改变第四项参量的值输入错误信息,MCNP,对输入文件出现的错误作广泛的检查，如果用户违反了输入说明的规定，将在终端上以及输出文件中打印致命错误信息，,MCNP,不再进行粒子输运计算，作业中断。
第一个出现的致命错误是真的，而后面的错误可能不一定是真的，这取决于前面出现的致命错误的情况若在,MCNP,运行行上指定,FATAL,项，则,MCNP,忽略致命错误，照常运行对于,MCNP,的警告信息，用户不应忽视，应搞清楚它们的含义检查几何错误,MCNP,在处理输入文件的数据时，不能检查一种非常重要的输入错误即,MCNP,无法查出各栅元之间的重叠和空隙，只有当粒子丢失时，才会发现几何错误即使如此，可能仍然无法准确判断错误性质栅元描述卡,格式：,j,m d geom params,或：,j,LIKE,n,BUT,list,j,栅元号，1,j,99999，写在第,15,列上m,栅元材料号，与材料卡（M,m,）中的序号对应m,0 为真空栅元d,栅元材料密度正值为原子密度，负值为质量密度对于真空栅元，该项缺省，不填写geom,栅元的几何说明由一系列带符号的曲面号经过布尔运算组成params,任选的栅元参数说明n,另一个栅元的名字（编号）list,描述栅元,j,和栅元,n,之间差别的栅元参数在栅元的几何说明中，关于曲面的指向是一个很重要的概念假定曲面,S,的曲面方程为,f,(,x,y,z,)0，则对于,f,(,x,y,z,)0的区域对于曲面,S,具有正的指向；而对于,f,(,x,y,z,)0的区域对于曲面,S,具有负的指向。
正指向的区域用,+S,表示，“,+,”号可不写；负指向的区域用,-S,表示栅元用各相关曲面的布尔运算表示，布尔算符包括,交,（用空格表示）、,并,（用冒号,:,表示）和,非,（用,#,表示）缺省的运算顺序是先,非,，其次是,交,，最后是,并,，使用括号可控制布尔运算的次序非,运算有两种形式：,#,n,，,n,是某个栅元号，,#,n,表示一个由不在栅元,n,内的点组成的空间区域),，,括号内是对某一个栅元进行描述的曲面栅元关系组，这一形式定义的几何区域由不属于括号内描述区域的点组成的空间例如：3 0 -1 2 -4$,定义栅元3,#3$,与下行相同,#(-1 2 -4),在栅元卡上可定义栅元参数以代替在输入文件中数据卡部分定义的栅元参数格式为：关键词值这儿允许的关键词是：带有粒子标识符的IMP、VOL、PWT、EXT、FCL、WWN、DXC、NONU、PD和TMP，以及关于重复结构的4个栅元参数卡：U卡、TRCL卡、LAT卡和FILL卡在,LIKE,n,BUT,格式中，还有两个关键词MAT和RHO，分别表示栅元的介质号和密度例如：,10 16 -4.2 1 -2 3 IMP:N=4 IMP:P=8,表示栅元10由曲面1的正面、曲面2的负面和 曲面3的正面的交集组成，填充质量密度为 4.2 克/厘米,3,的16号材料。
该栅元的中子重要性为4，光子重要性为8例如：2 3 -3.7 -1 IMP:N=2 IMP:P=4,3 LIKE 2 BUT TRCL=1 IMP:N=10,曲面描述卡,由方程定义曲面,格式：,j,n a list,j,曲面号，1,j,99999，写在第,15,列上如果曲面号前有,*,号，则该曲面为反射面n,对应坐标变换卡TR,n,，表示该曲面是在辅助坐标系下描述的，而该辅助坐标系与基本坐标系之间的关系由TR,n,卡给出如果没有坐标变换，即曲面是在基本坐标系下描述的，则该项缺省a,曲面助记符list,曲面方程参数，110项，取决于曲面类型参见MCNP手册，表3.1用点定义轴对称曲面,类型为,X,、,Y,或,Z,的曲面卡是用坐标点描述曲面而不是用方程系数描述用这些卡描述的曲面必须是分别关于,X,、,Y,或,Z,轴对称的，并且如果该曲面是由多叶组成的，则指定的坐标点必须全都在同一个叶上格式：,j,n a list,j,曲面号，1,j,99999，写在第,15,列上n,TR,n,卡的号，如。