内容发布更新时间 : 2025/3/4 19:01:29星期一 下面是文章的全部内容请认真阅读。
(1) MCNP中物理量的默认单位
长度:cm 通量:MeV
时间:刹 shake (10-8秒) 能量:MeV
温度:MeV (kT)
原子密度:1024个原子/cm3 质量密度:g/cm3
截面:巴 barns (10-24/cm2) 加热量:MeV/collision
此外, 原子质量按照中子质量为1.0计算,这种单位下阿佛伽德罗常数是0.59703109; 程序运行时间以分钟为单位。
MCNP的源代码是用FORTRAN语言编写的。
(2) 输入INP文件的基本形式
INP文件由一些被空行分隔的由一些被空行分隔的输入块组成,主要的输入块是信息块、标题和栅元块、曲面块和数据块等。输入块又由一些被称为卡的输入行组成。论坛的朋友可以参照对比一下,上一讲中的greatwall程序内容。 INP文件的格式如下:
信息块 (可选项) 空行分隔符 (可选项) 标题卡 栅元卡 ……
空行分隔符 曲面卡 ……
空行分隔符 数据卡 ……
空行分隔符 其它 (可选项)
说明如下:
? 信息块的第一行,必须在它的1~8列填写MESSAGE:,后面跟着用空格分隔的参数项。其后的9~80列和空行分隔符之前各行的1~80列都看作信息块内容。可用A=B参数项更改输出文件名,如OUTP = MYOUT。信息块是可选的。信息块提供给MCNP一个执行信息的方法。
? 在信息块之后的第一行是问题的标题卡,它仅限于一行,占用1~80列,可以是任何信息,将作为OUTP文件中各个输出表的标题被复制。
? 用户在栅元块和曲面块中描述问题的几何。栅元由栅元卡描述。空间必须由彼此相邻的栅元填满,栅元之间不能重叠,也不能出现无栅元的空区,否则会出现错误。构建栅元的曲面由曲面卡定义,曲面卡在曲面块中给出。曲面卡和栅元卡的填写方法,将在以后的讲座中予以介绍。
? 曲面块之后是数据块,在数据块中用户描述源、记数方式、材料等。数据卡在以后的讲座中予以详细介绍。
? 数据卡后不管有无空行分隔符均可以运行,不同之处是,如果数据卡后面有空行分隔符,则MCNP将不再读后面的附加行(如果附加行存在的话)。这对某些用户喜欢保留附加的内容是很方便的,如改变这个问题的说明或者与这个输入文件本身有关的信息,那么,这个分隔符就能防止读入这些附加信息。用户可以把希望保留的一些附加信息写在数据块的空行之后,MCNP会将它们复制到OUTP文件末尾。
(3)卡片格式
INP文件的每一行(称之为一张卡片)都限于使用1~80列。大多数卡片都可以按行填写,然而数据卡也允许按列填写。
a) 行输入格式
通常卡片的1~5列用于填写栅元号、曲面号或数据卡的助记名,6~72列填写卡片参数,73~80列为注释,$符号之后也为注释。序号或卡片助记名可以写在1~5列的任何地方。带有粒子标识符助记名可能需要5列以上,但冒号必须写在6列以内。如果1~5列空白,表示本行为上一行的接续行。6列之后可以写数据项,多个数据项之间用空格分隔。一个数据项必须在一行上写完,不得接续到下一行。相同编号的卡片只能有一张。
需要注意的地方: 每行最多80个字符?
不含控制字符,比如:Tab键,你是看不见的。它导致程序出错。? 注释行:?
—标题卡之后的任何位置都可插入?C?; —第一列是字母“C”,且随后四个空格; —从输入数据之后的$符号后开始的内容 (以上三种情况可以单独或同时存在)
? $符号为它所