mcnp从入门到精通讲座

第一讲： 如何运行MCNP程序

执行命令语句

\\mcnp i=(输入文件名) o=(输出文件名) …[选项]

选项 [i|p|x|r|z]

i=process input file 过程输入文件 p=plot geometry 几何画图

x=process cross-sections 截面处理 r=run particle transport 运行粒子输运

z=plot tally results,plot cross-sections 画记数结果,画截面。

具体执行过程：

假定mcnp.exe在C:\\\\mcnp\\bin\\子目录下， 在DOS窗口里，在命令行提示符>后敲入，

C:\\\\mcnp\\bin>mcnp i|p|x|r|z i=(输入文件名) o=(输出文件名)╛ 其中╛为回车键。 缺省文件名：

inp→ascii input file ascii码输入文件 outp→ascii output file ascii码输出文件 runtpe→binary restart file 二进制转储文件 举例如下：

C:\\\\mcnp\\bin>mcnp ixr i=greatwall name=greatwall.1 其中greatwall是你自己编写的源程序输入文件。

这个命令行语句，将读入输入文件greatwall，调用XS数据库，并且运行mcnp。

.1是你自己选择的输出文件名，每次输出必须更改输出文件名，或者不更改但删除上一次的输出文件名，否则程序运行出错。你也可以改成.2,.8,等等，任你选择。 这样运行的结果，就会有如下输出：

greatwall.1.o→ascii码输出文件 greatwall.1.r→转储文件

其中，转储文件也可以作为下一次运行的输入文件，进行接续运行。关于接续运行，将在以后的讲座中详细介绍。

你自己定义的文件名，将使你能够跟踪上述的辅助文件。

练习：

众坛友，请运行下面的greatwall文件，暂时不要计较文件里的内容，以后的讲座，会陆续予以讲解。程序的目的是，画出在一个立方体当中有一个球体。 message: datapath=C:\\mcnp\\bin

title --greatwall: U sphere in Al shpere in water cube--- c

c quick intro run, plot, interrupts etc c

c *** cell descriptions *** c

10 1 -18.7 -1 imp:n=1 $ Pu Sphere c

20 2 -7.8 +1 -2 imp:n=1 $ Al shell c

30 3 -1 +2 -3 imp:n=1 $ water c

40 0 +3 imp:n=0 $ outside world c

1 sph 0 0 0 5 $ Sphere at origin radius 2 cm c

2 sph 0 0 0 5.02 $ outter Al spherical shell 0.5 cm thick c

3 box +20 -20 -20 0 40 0 0 0 40 -40 0 0 $ surrounding H2O

c *** Source Card*** c

sdef $ default source c

c *** Material Cards *** c

m1 92238 1 $ U238 m2 26000 1 $ Al

m3 1001 +1 8016 +2 $ water c

c *** Print first 50 histories *** print 110 c dbcn 17j 1 j

prdmp 2j 1 1 c

fc4 flux integrated over the Pu/Al interface f4:n 10 f14:n 20

c *** time and nps cut cards *** nps 10 ctme 1

End of File

你可以用Wordpad或者Notepad打开任何一个你已经有的INP文件，删去其所有内容，拷贝上述程序，另存为greatwall文件。注意，不要敲Tab键，否则程序出错。

尝试下面的命令，看看你都得到了什么输出结果，以及输出的文件都在哪里？ (1)执行

C:\\\\mcnp\\bin>mcnp i=greatwall╛ 其中╛为回车键。 (2)再次执行一遍上述命令 C:\\\\mcnp\\bin>mcnp i=greatwall╛ (3)执行

C:\\\\mcnp\\bin>mcnp i=greatwall name=greatwall.1╛ (4)执行

C:\\\\mcnp\\bin>mcnp i=greatwall name=greatwall.2╛ (5)执行

C:\\\\mcnp\\bin>mcnp ixr i=greatwall name=greatwall.2╛

下面是mcnp程序运行的结构框图，从中你可以看到，用户输入文件被mcnp.exe执行后，要根据执行命令语句的设置，来调用XS数据库，然后计算输运过程，把输出文件放到所执行的子目录里。

我不知道，国内怎样翻译XS数据库，索性就把XS放到这里，不再忽发奇想，闭门造车，以免?毁‘人不倦。如果有哪位朋友知道的话，请跟帖加以说明，非常感谢。

这一讲，就说到这里。欲知后事如何，请听下一讲分解：mcnp输入文件。*^_^*

第二讲，MCNP输入文件

(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‖，且随后是个空格; —从输入数据之后的$符号后开始的内容 (以上三种情况可以单独或同时存在)

? $符号为它所在的那行数据的结束符，在$符号后面的内容作为注释，它可从$符号后面的任一列开始。

填写卡片参数时可以使用以下输入简写功能：

1>. nR功能，表示它前面的数据再重复n遍。例如：2 4R => 2 2 2 2 2 2>. nI功能，表示与其前后相临两个数之间插进n个线性插值点。 例如：1 5I 7 => 1 2 3 4 5 6 7

3>. xM功能，表示数值等于它前面数据的x倍。 例如：5 4M => 4 20

4>. nJ功能，表示从它所在位置跳过n项不指定的数据而使用缺省值。 这四项输入简写功能可以综合运用。

另外，FCn卡是作为注释用的，但它将作为记数类型n的表头文字印出，比如可作为记数的标题。

SCn卡也是作为注释用的，它作为源概率分布n的表头文字印出。具体使用方法请见手册。

b) 列输入格式

列输入只能用于数据块中，对栅元参数和源的描述比较有用。按行输入的栅元重要性、体积、权窗等数据项可读性较差，而且增加或删除栅元时要在行输入卡上仔细寻找相应项。列输入的可读性有很大提高，删除或增加与某一栅元相对应的数据项时也比较方便。

列输入格式的第一行以#开始，#可以放在1～5列的任意位置，卡片助记名逐个放在该行6列以后，在这些助记名之下按列给出数据项。同一个列输入格式块中的卡片必须是同一类卡片，比如都是栅元参数卡、都是曲面参数卡或都是源参数卡等，在#号下面的1～5列放置栅元号、曲面号或源分布号。

c) 粒子标识符

粒子标识符是卡片助记名的一部分，:N表示中子，:P表示光子，有时也能遇到表示中子-光子联合输运的:N P。下面一些数据卡需要粒子标识符：IMP、EXT、FCL、WWN、WWE、WWP、WWCE、DXT、DXC、Fna、PHYS、ESPLT和CUT。

d) 缺省值

MCNP许多输入卡的参数项有缺省值，用户不必每次都给出这些参数，如果卡片输入项有固定顺序，可以使用nJ功能跳过n个输入项。如果卡片上所有数据项都想缺省，只给出卡片助记名即可。有些卡片不给出也有缺省值，如MODE N卡就可以省略。

练习：

请找出上一讲中的greatwall文件里的信息卡MESSAGE:，标题卡Title,栅元卡cell Card(从cell descriptions开始)，曲面卡Surface card(从1 sph 0 0 0 5 $ Sphere at origin radius 2 cm 开始的三行),数据卡Data card，包括源卡Source Card，材料卡Material Cards 。论坛朋友们不必彻底搞懂它们的作用，以后讲座会详细讲解；只需体会书写的格式，看看是否符合上述的规则。

这一讲，就说到这里。欲知后事如何，请听下一讲分解：MCNP几何结构描述。*^_^*

MCNP从入门到精通讲座，第03讲

再举一个例子：有一个大球面S1，它的里面有一个小球面S2。在小球面S2外且在大球面S1里的部分，是这样定义的： 在小球面S2里的部分，是这样定义的：

(2) 交集和并集(intersection & union)

在MCNP里，布尔表达式是由下述符号构成的： ：冒号，表示并集； # 井字号，表示补集；

表示交集。这个符号你们看见了吗？我也没看见。交集的表达式是一个空格键，也是隐式的。

缺省的运算顺序是先#补集，其次是交集，最后是：并集。

交集是两个集合的公共部分，如图所示。

并集是两个集合的合集，见图。

若想仅用交算符定义栅元，栅元内所有点对特定的界面必须有相同的坐向。如图所示，1、2、3、4、5面围成一个实体，用密度为1.0g/cm3的材料1填充，外部是真空。

实体内P1和P2点相对于4面有负坐向，P3点相对于4面有正坐向；P1点相对于3面有正坐向，P2和P3点相对于3面有负坐向。这时最好借助另一平面6把实体分成两部分，从而顺利定义栅元：

1 1 -1.0 1 –2 –3 6 2 1 -1.0 1 –6 –4 5

然而，只用交算符还不便定义实体之外的区域，这时用并算符比较方便，实体外可由3、4面以上区域的交集和其它几个面以外区域的并集构成一个栅元：

3 0 -5 : -1 : 2 : 3 4

当然实体之内的区域也可以借助并算符定义，把1、2栅元合并(这样定义会降低计算速度)。由于交运算优先，所以要用括号把并运算括起来：

1 1 -1.0 1 -2 (-3 : -4) 5

(3) 补集(complement) 补算符#的用法有两种：

<1>. #n 表示当前定义的栅元是栅元n其余的部分。 <2>. #(...) 表示括号内曲面号定义区域之外的部分。 引入坐向概念后，补运算就不是交运算和并运算之外的新概念，而是交和并的另一种表示方法。例如，对上面第一种情况的处理过程是，去掉#号，把n用括号括上，n中的交算符变成并算符，n中的并算符变成交算符而且加上背靠背的括号) (，曲面之前的坐向取相反符号。数学中常用?表示交运算、用+表示并运算、ˉ

表示补运算，有 ， ，此处 和a则用坐向指定。 使用补算符时要格外小心，用户常因对补算符理解不透而弄出几何逻辑错误，因此建议不使用补算符。另外，#(n)是不合规定的。

练习题：请写出下面图形中的各个栅元几何描述。其中栅元1是外部世界栅元；栅元2是大圆柱体中除了那个圆柱空洞的部分；栅元3是那个圆柱空洞；栅元4是侧面贴上去的那个圆柱体。

这一讲，就说到这里。欲知后事如何，请听下一讲分解：栅元卡Cell card。*^_^*

第04讲，栅元卡Cell card

第二节 栅元卡

栅元卡的格式： j m d geom. Param.

分别对应于：栅元标号 材料号 材料密度 几何描述 参数说明 举个例子： 1 2 -.98 1 -2 3 -4 5 -6 imp:p=1

解释如下：

j: 在1～5列上填写的栅元标号，可以不连续。MCNP按照读入的顺序对栅元另行编号，称为栅元的程序编号，从1开始单调上升。为避免程序编号和栅元标号

混淆，建议按顺序给出栅元标号。栅元标号的范围：1≤j≤99999。

m: 该栅元的材料号，是材料卡(Mm)中相应材料的序号。真空栅元，m=0。 d: 栅元材料的密度。填入正值时，表示是原子密度，单位1024个原子/cm3。填入负值，则表示是质量密度，g/cm3。对于真空栅元，该项缺省，不填写。 geom: 栅元的几何描述。它列出界定该栅元的所有曲面的号(有数值符号，表示坐向)，以及描述这些曲面所定义的区域之间关系的布尔算符。布尔算符包括交算符(缺省)、并算符( : )和补算符(#)。缺省运算顺序是，先补，再交，最后并。可用括号改变运算次序。

param: 可选的栅元参数说明，其形式为：关键字=某一值。

在栅元卡上可以定义栅元参数，以替代在输入块中定义栅元参数。这里允许的关键字是：带有粒子标识符的IMP、VOL、PWT、EXT、ECL、WWN、DXC和TMP。栅元卡参数项上的等号可用空格代替。若在栅元卡上指定了某种栅元参数，则在数据块中不能重复指定。

即使是十分有经验且细心的用户，在构建比较复杂的几何时，也经常会犯一些错误。有些错误程序可以检查出来，有些则不能，以致造成错误结果。几何构建还直接影响程序运行效率，不合适的几何构建将浪费机时。

A.定义栅元的原则

原则1：栅元界面要少

尽管可用大量的并算符定义一个几何上比较复杂的栅元，但这样做并不明智。问题在于，每次计算栅元内碰撞之间的径迹时，都计算径迹和栅元所有界面是否相交以判断粒子是否进入其它栅元，太多的栅元界面会浪费机时。如，图a中的几何只是一些平行的柱面，较容易定义，然而两个大柱面之间的栅元有14个界面(包括前后界面)；图b则是比较有效的栅元定义，复杂界面的栅元被分隔成一些

小栅元，每个栅元的界面有所减少。

原则2：栅元大小要合适

构建物理空间时，除按不同材料区域分割栅元外，还必须考虑效率因素。为了节省运算时间，栅元不能定义得太大或太小。在比较重要的区域，建议调整栅元大小，使粒子进入相邻栅元后数量约减少一半，以利于使用几何分裂和轮盘赌技巧。

原则3：尽量用简单曲面定义栅元

如前所述，MCNP频繁计算碰撞点到界面距离。点到二阶以上曲面的距离一般不能用解析法求解，用牛顿下山法求解时，多次叠代很浪费时间。如果问题条件允许，尽量使用二阶以下曲面分段替代二阶以上曲面。

原则4：避免曲面和曲面相切

即使使用双精度变量，也不能计算出点到曲面的确切距离，因而在曲面和曲面相切处有时会因计算精度不够导致栅元重叠或空缺，造成粒子丢失。虽然这种粒子丢失几率很小，一般不会影响记数结果，但计算大量粒子问题时可能会因为丢失粒子数量超过限制，使程序停止运行。

B.检查几何错误

MCNP检查INP文件的数据时，不能检查出几何逻辑错误。而运行中粒子丢失时，才能检查出栅元的重叠或栅元之间的空缺。因此在正式运算前，应采取以下办法检查几何逻辑错误： 1)增加一个VOID卡。这个卡将废弃这个问题的材料说明，将全部栅元设为真空，并把加热记数转化为通量记数。

2)对这个问题增加一个大球面，这个球面包围需计算的系统，球面把系统之外的栅元分成两个栅元。球内面是现在问题的边界，内部所有栅元的重要性都置为1。 3)源描述改为：SDEF SUR=m NRM=-1，其中m是刚才定义大球面的标号。 由于上述改动使粒子在计算系统中没有碰撞，一个短时间运行就会使大量粒子的轨迹穿过系统，可能在出现几何错误的地方引起粒子丢失。 粒子丢失时，OUTP文件上将产生丢失粒子打印表。表中列出丢失粒子跨越的所有曲面，并告诉你粒子在什么位置上向什么方向运行后丢失。你可以从这些信息中推断出粒子丢失的原因。

下面是栅元卡的一个例子，请各位朋友以此体会一下栅元卡的使用。 c *** cell descriptions *** c

c **** cell descriptions ****

10 1 -0.95 +10 -20 +30 -40 +50 -60 imp:n=1 imp:p=1 c

20 0 #10 imp:n=0 imp:p=0

c !!!! Notice the imp:p=1 -ensures that PHOTONS ARE NTRANSPORTED c

10 pz -5 $ bot c

20 pz +5 $ top c

30 px -5 $-x c

40 px +5 $+x c

50 py -5 $-y c

60 py +5 $+y c

注意：其中#10 是说除了栅元10的那些空间，即非10，或者10的补集。

这一讲，就说到这里。

欲知后事如何，请听下一讲分解：曲面卡Surface card。*^_^*

言归正传， 第05讲，曲面卡Surface card

第三节 曲面卡Surface card

A. 由方程定义的曲面

格式： j n a list 解释如下：

j ：曲面的序号，写在1～5列上。它的范围 1≤j≤9999。可以不连续，即可不按顺序填写曲面的序号。如果这个曲面是反射面，在这个曲面序号前边加一个―*‖号，来表示反射面。加―+‖号，则表示空白界面。如果由曲面定义的几何块有坐标转换TRCL，则它的范围变成 1 ≤j≤999。

n ：相应于TRn卡的整数n。请见附表。如果缺省，或为零，表示没有坐标转换。

它如果大于0，就是指定 TRn 卡的序号。它表示该卡定义曲面是在辅助坐标系下描述的，而这个辅助坐标系与基本坐标系之间的关系由 TRn 卡给出。 它如果小于0，就是指定面j 是面n 的循环。 a ：方程助记名，见下表。

List ：按照规定次序给出方程描述的数据项(1～10项)。

下表列出了MCNP识别的曲面的类型、方程、助记名及相应的卡片数据项次序。

用这种方法定义一个曲面首先要从表中找到所需的曲面助记名，然后根据其方程的形式算出所需的系数(参照解析几何)。按照前述的书写格式将他们依次写在一张卡片上。在某些情况下，可以用几何系统中的点定义曲面，见B、C 部分的讨论。曲面还可以用宏观几何体的组合来定义，见D部分。

如果一个曲面在点（x，y，z）的计算值是正的，则称这个点对于这个曲面具有正的方向；一个曲面的表达式是表中曲面方程的左边。位于球、柱，锥及环的正的方向，统一是这些曲面的外边。对位于垂直于x、y、z 轴的平面px 、 py、pz 正的方向分别是 x，y 或z 值大于相应平面截距的那些点。对P，SQ 及GQ 曲面，用户提供这个表达式的全部系数，于是用户就能够随意确定这个曲面的方向。这与其它情况不同。其它情况的方向即使是任意的，也是由表达式唯一确定的。因此，通过一个曲面变换（见TRn 卡的手册说明），PX，PY，或PZ 曲面可能必须由P 曲面代替以防止这个曲面的方向相反。

如果曲面序号在*号之后，则定义的这个曲面为反射面。当粒子打在这样的曲面时，便按镜面反射。如果曲面号前加+号，则这个曲面被定义为空白边界。通常在使用反射面或空白边界时，不应有探测器计数或 DXTRAN 球（下一事件估计）。具体请见手册。对于反射面的计数问题需要不同地对其标准化。

如果第二个输入n 为负值，则定义的曲面j 关于曲面k 周期变化。遵循的规则如下：

? 曲面k 和j 必须是平面。

? 平面转换不允许是周期性的平面。

? 周期性的块可以是无穷的也可以由在它的顶部和底部的平面界定，界定的平面可以是反射面也可以是空白面，但不能是周期的。 ? 周期平面仅可以界定其它的周期平面或上下底面。

? 每个周期平面的某个侧面必须是一个单一的零重要性的块。

? 所有的周期平面必须有一个共同的垂直于几何的顶部和底部的旋转顶角。 ? 下一事件估计，例如探测器计数或DXTRAN 球在此不能使用。

大家来看下面的例题：

例1： 1 PY 3 这是描述：曲面1，在y=3 处，垂直于Y 轴的平面,y>3 的全部点都有正的方向。你能想象出吗？拿笔画一下，就领悟了。

例2： 2 S 4 1 －3 3.62

这是描述：曲面2，球心在(4,1,-3)点，半径为3.62 cm 的球体。就是球心的坐标为 x=4,y=1,z=-3。这个比较容易想象。

例3： 3 SX 10 3

这是描述：曲面3，球心在x轴上，(10,0,0)点，半径为 3 cm 的球体。你画一下试试。

例4: 4 K/Y 0 0 2 .25 1

这个曲面4，指定顶点在(0,0,2)处的一个圆锥面,对称轴平行于Y 轴。锥的斜率t 是0.5 (注意，在卡片上输入的是t的平方，要小心啊)，而且只使用其正方向（右边）的一支。锥外的点具有正方向。

例5： 5 GQ 1 .25 .75 0 -866 0 -12 -2 3.464 39

这个曲面5，是半径为1cm 的圆柱面,其对称轴在垂直于X 轴的一个平面上,该平面在X 轴上的截距为6cm。与X 轴相距2cm 的地方围绕X 轴从Y 轴向Z 轴旋转30度角。在柱外的点具有正的方向。该曲面在辅助坐标系下描述是很简单的，将这个柱面的对称轴定义为辅助坐标系的X轴，然后再给一张TRn 输入卡(请见手册)用来定义基本坐标系与辅助坐标系的关系。

在辅助坐标系下这个曲面5 描述如下： 5 7 CX 1

*TR7 6 1 －1.732 0 30 60

等到我们讲述TRn 的讲座时，将会详细讲解这个例子。

TX，TY，及TZ 输入卡描述的是椭圆形面(4 阶曲面)，其旋转对称轴分别平行与X 轴，Y 轴，Z 轴。下图a 给出一个TY 环面。注意到输入参数x , y , z , a, b, c 给定一个在(r,s)柱坐标系统围绕S 轴旋转的椭圆(下图b)。其原点位于原坐标系的点(x , y , z )处，这椭圆环面方程可写成： 当|A|椭圆环面的坐标变换是局限于辅助坐标系统的每个轴都要平行于基本坐标系的一个轴。

呵呵，这一讲，我基本上是按照MCNP 5.1.40手册 来讲解的。 课后练习：请画出下列曲面卡所代表的曲面。 1 PX 8 2 SO 11.1 3 C/Z 0 1 3 4 KZ

5 C/X 0 1 2

答案：(注意：对于助记名来说，大小写都行，对大小写不敏感。) 1 px 8

–Plane normal to the x-axis (x-D=0) data=D

2 so 11.1

–Sphere at the origin (x2+ y2+z2-R2=0 data R

3 c/z 0 1 3

Cylinder parallel to the z axis (x y R)

4 kz

Cone on z-axis (z, t2+1)

5 c/x 0 1 2

cylinder || to x(y,z,r)

This specifies a cylinder lying parallel to the x-axis at a distance y=0, z=1 cm.

The radius of the cylinder is 2cm.

这一讲，就说到这里。欲知后事如何，

请听下一讲分解：曲面卡Surface card (续)。*^_^*

第06讲，曲面卡Surface card (续)

B.用点定义对称曲面

格式： j n a list

j: 在1～5列上填写的曲面序号，它的范围：1≤j≤9999，可以不连续。如果曲面用来定义几何块，则用TRCL转换，这时它的范围：1≤j≤999。

n：若缺省表示无坐标转换，或相应的TRn卡n 值。 a: 字母X、Y、或Z

list: 一至三对点坐标。类型为X,Y,及Z的曲面卡是用坐标点描述曲面而不是用方程系数来描述曲面。用这些卡描述曲面，必须是分别关于X,Y,Z 轴对称，并且如果这个曲面是由一叶以上组成的，则指定的坐标点必须全都在同一个叶上。 每一对坐标点定义在这个曲面上的一个点。例如，在一张Y 卡上可以给出： j Y y1 r1 y2 r2

这里ri=((xi)2+(yi)2)1/2，yi 是点i的坐标。

给出的点数不同，曲面类型也不同，例如：

?

若给出一对坐标，则定义的是平面，即PX，PY，或PZ。

?

若给出两对坐标，则定义的是线性曲面，如PX,PY,PZ,CX,CY,CZ,KX,KY 或KZ。

?

若给出三对坐标，则定义的是二次曲面，包括PX，PY，PZ，SO，SX，SY，SZ，CX，CY，CZ，KX，KY，KZ 或SQ。

当用两个几何点定义锥面时，则只生成一个单叶锥面。

用这个规定所确定的这些曲面（SQ 除外）的方向与方程指定曲面得到的方向是等同的。

而对于SQ 来讲，是这样定义方向的，即离对称轴足够远的点有正的方向。这与方程定义的SQ 不同，用方程可自由选取方向。

例1：1 x 7 5 3 2 4 3

这是描述关于x 轴对称的曲面，该曲面通过三个(x，r)的点(7，5)，(3，2)，及(4，3)。这个曲面是有两个叶的双曲面，在MCNP 中这个曲面被转换成它的标准方程形式 1 SQ -.083333333 1 1 0 0 0 68.52083 -26.5 0 0

例2： 2 Y 1 2 1 3 3 4

这是描述了在Y=1 及Y=3 的两个平行平面，但这个描述是错误的，FATAL error,因为没有满足所有的点都在同一个叶的要求。

例3： 3 Y 3 0 4 1 5 0

这是描述半径为1 的一个球，球心在(x，y，z)＝(0，4，0)。

例4：4 Z 1 0 2 1 3 4

这个曲面被舍弃，因为这些点在双曲面

x2+y2-7z2+20z-13=0 的两个不同叶上。然而，曲面：

4 Z 2 1 3 4 5 9.380832 和上面曲面有同样曲面方程，而且是可以接受的，因为所有坐标点都在双曲面右叶的一个曲面上。

例5： 1 0 1 －2 3 $几何块1

1 Y -3 2 2 1 2 Y 2 3 3 3 4 2 3 Y 2 1 3 1 4 2

最后这个例5定义了一个由一叶锥面，双曲面及一个椭球面界定的几何块。呵呵，比较复杂一点。这三个曲面定义的是关于y 轴对称的环形几何块。这个几何块的横断面如下图所示。要绘出此图，输入PX=0 EX=5。第一个曲面通过点(-3，2)及(2，1)。第二个曲面通过(2，3),(3，3)及(4，2)。由点(2，1)，(3，1)，(4，2)定义最后一个曲面。这些坐标点都是用(y，r)形式。利用这些卡，MCNP 指出曲面1 是一叶锥面，曲面2 是一个椭球，曲面3 是一叶双曲面。当使用PRINT 输入卡或执行选择时，将打印出各种曲面的标准的方程的系数。例如，SQ 卡定义的曲面3 是：

3 SQ 1 -1.5 1 0 0 0 -.625 0 2.5 0

C．由三点定义一般平面

格式：j n P X1 Y1 Z1 X2 Y2 Z2 X3 Y3 Z3

j ：在1～5列上填写的曲面序号。它的范围：1≤j≤99999，可以不连续。或者在重复结构中是1≤j≤999。

n ：没有坐标变换时则省略，或为0。 n大于0，为指定TRn卡的n 值。

n小于0，是指定曲面j 是曲面n 的循环。 (Xi,Yi,Zi) : 定义该平面的点的坐标。

MCNP 对于用户指定的P 型曲面，将检查所给的数据项个数，若是4 项，则作一般平面方程Ax+By+Cz-D=0的系数理解；若多于4 项时，便作为三维空间点的坐标值理解。每三个数定义空间一个点，由MCNP 把它们转换成所需要的曲面系数以产生平面：AX+BY+CZ-D=0。 所产生的平面方程系数遵循如下原则：

?

应使坐标原点关于该平面有负坐向；

?

当该平面通过坐标原点(D=0)时，则应使得点(0，0，∞)相对于该平面具有正坐向；

?

若以上两项都无法满足做到(即D=C=0)，则应使点(0，∞，0)相对该平面具有正坐向；

?

若以上三项都无法满足做到(即D=C=B=0)，则应使点(∞，0，0)相对该平面具有正坐向；

?

若以上四项都无法满足做到，那么用户指定的三个点一定处于一条直线上，则打印出致命错误FATAL error信息。

D.由宏观物体Macrobody定义的曲面

为什么要用宏观物体来定义曲面？答案是简洁方便。我给大家举个例子：对下面的图形，使用栅元和曲面，对每一个面来定义，就需要12个曲面。但如果我们使用宏观物体，仅用2个宏观物体就够了。

利用组合几何体，像宏观物体定义几何块或者曲面是一种可选的方法。在综合泰格系列(ACCEPT)编码中，这些组合几何体的可用性相似。这些宏观物体能与标准的几何块和曲面相混合。这些宏观物体的面被分解为曲面方程，各个小的曲面根据预先的次序重新赋予序号。序号的选择是用户指定的，其后跟有十进制点和数字。这些小的曲面能用来计数、计数分割、定义其他的几何块、SDEF 卡等。但是不能用在SSW/SSR卡、衰减面卡、PTRAC、或者MCTAL文件。

考虑到宏观体的曲面和所有小的面，其内部的空间具有负的坐向，外部空间则具

有正的坐向。这些小的曲面的坐向由宏观体―主‖几何块确定，在其他几何块描述中用到这些小曲面时，其将保持坐向值，这时要做适当的注释。见后面的例题说明。

以下宏观物体是可用的，其定义如下: BOX 任意方向的立方体 RPP 长方体 SPH 球

RCC 正圆柱

RHP or HEX 正六边棱柱 REC 正椭圆柱 TRC 截断直角锥 ELL 椭球 WED 楔体

ARB 任意多面体

下一讲，我们将分别予以解释。

可能有的朋友要问了，既然宏观物体来定义曲面，既简洁又方便，那我们就都用宏观物体，不用曲面卡行不行？我的答案是不行，我们还离不开曲面卡。大家看看上面的宏观物体，只有10种，是有限的。而现实世界的曲面形状千变万化，不可能仅用这10种就能定义所有曲面。在那些比较复杂的形状里，我们还是要用曲面卡来定义曲面的。

练习题1：请画出下面的几何形状。

c *** cell descriptions *** c

10 1 -18.7 -1 imp:n=1 $ Pu Sphere c

20 2 -7.8 +1 -2 imp:n=1 $ Al shell c

30 3 -1 +2 -3 imp:n=1 $ water c

40 0 +3 imp:n=0 $ outside world c

1 sph 0 0 0 5 $ Sphere at origin radius 2 cm c

2 sph 0 0 0 5.02 $ outter Al spherical shell 0.5 cm thick c

3 box +20 -20 -20 0 40 0 0 0 40 -40 0 0 $ surrounding H2O

这是第一讲中的例子，在greatwall文件里。请大家体会一下其中曲面的定义。栅元10里的-1，是指下面曲面1里sph的负坐向；同理，栅元20里的-2，是指下边曲面2里sph的负坐向；栅元40里的+3，是指下边曲面3里box的正坐

向。大家看看那4个栅元都是什么形状？

练习题2：请画出下面的几何形状。

c *** cell descriptions *** c

1 3 -0.98 -1 -2 +3 imp:P=1 2 0 +1:+2:-3 imp:p=0 c

1 cz 10 c

2 pz +10 c

3 pz -10

大家看看这两个栅元又都是什么形状？

这一讲，就说到这里。欲知后事如何，

请听下一讲分解：宏观物体定义曲面详解。*^_^*

第07讲，宏观物体定义曲面详解

(1) BOX 任意方向的直角盒子(所有角都是90度)

它的规范格式是：

BOX Vx Vy Vz Alx Aly Alz A2x A2y A2z A3x A3y A3z

这里 Vx Vy Vz = x,y,z 角度的相应数据 Alx Aly Alz = 第一个面的方向矢量； A2x A2y A2z = 第二个面的方向矢量； A3x A3y A3z = 第三个面的方向矢量。

例如，BOX －1 －1 －1 2 0 0 0 2 0 0 0 2

就是一个中心在原点，每边长2cm，面的法线是主轴的立方体。

我再给大家举个不是立方体的例子，这样便于大家深刻理解。

图中(3, 9, 8)就是长方体的起始点C1的坐标，即在C1点，x=3,y=9,z=8. 从C1点到C2点，是沿着x轴方向移动，y坐标和z坐标没有变化，只有x坐标增加了1cm。所以，各坐标的变化量是(+1, 0, 0),它们就是Alx Aly Alz。同样的道理，从C2点到C3点，是沿着y轴方向移动，x坐标和z坐标没有变化，只有y坐标增加了4cm。所以，各坐标的变化量是(0, +4, 0),它们就是A2x A2y A2z。从C3点到C4点，是沿着z轴方向移动，x坐标和y坐标没有变化，只有z坐标增加了5cm。所以，各坐标的变化量是(0, 0, +5),它们就是A3x A3y A3z。 所以，最后这个长方体的定义，就是： BOX 3，9，8 +1 0 0 0 +4 0 0 0 +5

回过头来，再看上面立方体的例子，是不是太简单了？ BOX －1 －1 －1 2 0 0 0 2 0 0 0 2

注意BOX的定义，可以不一定是沿着x,y,z轴的方向，换句话说，可以是任意方向，只要你知道那4个点的具体坐标，就可以定义它。实际上，对于曲面垂直于主轴的长方体，有特殊的定义RPP。

(2) RPP 面的法线是主轴的长方体, x,y,z的值源于坐标原点

它的规范格式是：

RPP Xmin Xmax Ymin Ymax Zmin Zmax

例如，RPP －1 1 －1 1 －1 1

就是上述立方体的中心在原点，每边长2cm，面的法线是主轴的立方体。

(3) SPH 相当于一般球的曲面方程

它的规范格式是：

SPH Vx Vy Vz R

其中，x,y,z=Vx Vy Vz 是球的中心坐标值，以及半径是R。

(4) RCC 像罐子一样的正圆形柱面

它的规范格式是：

RCC Vx Vy Vz Hx Hy Hz R

其中，正圆形柱面一个底面中心在(Vx ，Vy ，Vz)点，圆柱轴矢量(底面中心为起点)是(Hx, Hy, Hz)，以及半径R。

例如，RCC 0 -5 0 0 10 0 4

就是一个底座平面在y=-5处，轴矢量是正向Y轴，高是10cm，半径是4cm的圆柱面。

我再给大家举一个比较复杂一点儿的例子。

RCC 0 8 1 0 5 8.66 3

这个正圆形柱面的半径是3cm,这很简单。它沿z轴方向，却绕x轴顺时针旋转30度角。(0, 8, 1)是正圆形柱面的底面中心P点的坐标。而从P点到正圆形柱面的另一个面中心P2点，x坐标没有变化，都是0，y坐标从8变到13，变化量是13-8=5 cm。z坐标从1变到9.66，变化量是9.66-1=8.66 cm。所以，各坐标的变化量是(0, 5, 8.66),它们就是Hx, Hy, Hz。

(5) RHP或者HEX 正六边棱柱

RHP不是严密的符合Integrated Tiger Series (ACCEPT)公式。 它们是与RCC相似但底面为六边形的棱柱。

它的规范格式是：

RHP或HEX v1 v2 v3 h1 h2 h3 r1 r2 r3 s1 s2 s3 t1 t2 t3

这里，v1 v2 v3 就是x,y,z，也就是六角形底面相应的中心坐标； h1 h2 h3 是从底面中心到顶面中心的矢量；

对于高是h的z轴矢量，h1 h2 h3 将是 0 0 h。 对于高是h的y轴矢量，h1 h2 h3 将是 0 h 0。

r1 r2 r3是从轴到第一个侧面中心的矢量。对于斜度为2p而且垂直于y轴的面， r1 r2 r3=0 p 0。

s1 s2 s3 是指从轴到第二个侧面中心的矢量。 t1 t2 t3是指从轴到第三个侧面中心的矢量。

例题：RHP 0 0 －4 0 0 8 0 2 0 它是围绕Z 轴的六棱柱，底面在z＝－4 处，高为8cm，第一个面垂直于Y 轴，过y＝2。

宏观物体(Macrobody)的?面‘

这些几何体底面是有限的，能被用于MCNP 卡中。BOX 和RPP 在维数上是无限的，这样，将跳过这两个面，所保存的面的编号将减少两个。RHP 在轴向维度上是无穷的，这样就不存在面7、8。面的编号顺序紧随每个几何体。面的编号在几何绘图中能用MB0DY＝OFF 命令显示出来。

例如：

BOX

1 平面法线终于A1x A1y A1z 2 平面法线始于A1x A1y A1z

3 平面法线终于A2x A2y A2z 4 平面法线始于A2x A2y A2z 5 平面法线终于A3x A3y A3z 6 平面法线始于A3x A3y A3z

请看一个具体的实例。

又例如： RPP

1 平面Xmax 2 平面Xmin 3 平面Ymax 4 平面Ymin 5 平面Zmax 6 平面Zmin

SPH 作为一个规则面用，因而没有小曲面。

再例如：

RCC

1 半径为R 的圆柱面 2 平面法线终于Hx Hy Hz

3 平面法线始于Hx Hy Hz

还有：

RHP or HEX

1 面法线终于r1 r2 r3 2 面1 的对立面

3 面法线终于s1 s2 s3 4 面3 的对立面 5 面法线终于t1 t2 t3 6 面5 的对立面 7 面法线终于t1 t2 t3 8 面法线始于h1 h2 h3

以上的输入文件描述了五个几何块，图解了不同几何体和几何块/曲面定义的联合定义。曲面序号以斜体列在所描述的曲面旁边。注意的是几何块和曲面的序号可以不是从1 开始的，也可以不是连续的。

这一讲，就说到这里。

欲知后事如何，请听下一讲分解。*^_^*

第08讲，MCNP 几何部分栅元和曲面 练习题

MCNP 几何部分栅元和曲面 练习题

Geometry problems I – Cells and Surfaces.

Construct the MCNP input files (use the MCNP plotter to check the geometry in each case) for the following four problems.

Your input file should explicitly use: ? mode: n what does this mean?

? use a standard default standard source – What is this source?

You should ensure that tables showing source, cell, importance and the first 50 histories are produced. Use a time cut-off of 60 seconds and a limit on the number of particles of 60.

You should use cells and surfaces (not Macrobodies) in each case.

PROBLEMS:

1. A cylinder lying along the z-axis. Radius 5 cm height 10 cm. Make the cylinder from material U235 weight fraction 94.73, U238 weight fraction 5.27, material density 18.7g/cm3。

What changes would you make if you have atom fractions?

2. The cylinder in 1. Surrounded by a 20 cm side length cube of water. Position the cylinder centrally in the water.

3. The geometry of problem 2, but with the cylinder rotated about the x-axis by 45 degrees.

4. An Al box, 2 cm wall thickness and internal side 20 cm. The box has a top and a bottom. Inside of the box place a 5 cm radius cylinder of height 10 cm laying long the z-axis. The cylinder should be central within the box. The space between the cylinder and box is filled with water. The Al walls are of density 2.7 g/cm3. The cylinder has the material properties described in problem 1.

练习题答案： 例一：

message: datapath=D:\\transport\\data\\mcnpxs

title -- Problem 1: cylinder along z axis using surfaces --- c

c ************************************************************ c

c !!! REMEMBER TO CHANGE THE MESSAGE BLOCK FOR YOUR SET-UP !!! c

c ************************************************************ c

c *** cell descriptions *** c

10 1 -18.72 -1 -5 +10 imp:n=1 $ cylinder c density in g/cm^3 (-ve) c

20 0 +1:+5:-10 imp:n=0 $ outside world c

1 cz 5 $ infinite cyl along z-axis rad 5 cm c

5 pz +5 $ top of cylinder at 5cm

10 pz -5 $ bottom of cylinder at -5 cm c

c *** Data Cards after blank line***

c *** Source Card*** c

sdef $ default source: 14 MeV point neutron at (0,0,0) c

c *** Material Cards *** c

m1 92235 -94.73 92238 -5.27 $ U235/238 using weight fractions (-ve) c

c ** Print first 50 histories print table 110 *** print 110 $ add table 110 to minimum print-out c

c *** time and nps cut cards ***

nps 60 $ nps 60 - stop simulation after 60 source particles ctme 1 $ ctme 1 - stop simulation after 1 min (60 sec) c

c *** Blank line terminator ***

End of File

例二：

message: datapath=D:\\transport\\data\\mcnpxs

title -- Ex 1b cylinder along z axis, in water using surfaces --- c c

c ************************************************************ c

c !!! REMEMBER TO CHANGE THE MESSAGE BLOCK FOR YOUR SET-UP !!! c

c ************************************************************ c

c *** cell descriptions *** c

10 1 -18.72 -1 -5 +10 imp:n=1 $ cylinder c

c Use of the not operator - use with a cell ie #10 means not cell 10 c 15 2 -1 -15 +20 -25 +30 -35 +40 #10 imp:n=1 $ box but not cylinder c

c Expanded definition use of : operator, use of parenthesis and continuation

15 2 -1 (+1 -15 20 -25 30 -35 40):(-15 20 -25 30 -35 5) :(-15 20 -25 30 40 -10)

imp:n=1 $box but not cylinder

20 0 +15:-20:+25:-30:+35:-40 imp:n=0 $ outside world c

1 cz 5 $ infinite cyl along z-axis rad 5 cm c

5 pz +5 $ top of cylinder at 5 cm

10 pz -5 $ bottom of cylinder at -5 cm c

15 px +10 $ +ve x extent of warer 20 px -10 $ -ve x extent of water 25 py +10 $ +ve y extent of water 30 py -10 $ -ve y extent of water 35 pz +10 $ +ve z extent of water 40 pz -10 $ -ve z extent of water c

c *** Source Card*** c

sdef $ default source c

c *** Material Cards *** c

m1 92235 -94.73 92238 -5.27 $ U235/238

m2 1001 +1 8016 +2 $ water (using atome fractions (+ve)) c

c *** Print first 50 histories *** print 110 c

c *** time and nps cut cards *** nps 60 ctme 1

End of File 例三

message: datapath=D:\\transport\\data\\mcnpxs

title --Ex 1c cylinder in z-y plane rot. about x in water using surfaces -- c

c ************************************************************ c

c !!! REMEMBER TO CHANGE THE MESSAGE BLOCK FOR YOUR SET-UP !!! c

c ************************************************************ c

c *** cell descriptions *** c

10 1 -18.72 -1 -5 +10 imp:n=1 $ cylinder c

c Use of the not operator - use with a cell ie #10 means not cell 10 c 15 2 -1 -15 +20 -25 +30 -35 +40 #10 imp:n=1 $ box but not cylinder c

c Expanded definition use of : operator, use of parenthesis and continuation

15 2 -1 (+1 -15 20 -25 30 -35 40):(-15 20 -25 30 -35 5) :(-15 20 -25 30 40 -10)

imp:n=1 $box but not cylinder

20 0 +15:-20:+25:-30:+35:-40 imp:n=0 $ outside world c

1 1 cz 5 $ infinite cyl along z-axis rad 5 cm

c

5 1 pz +5 $ top of cylinder at 5 cm

10 1 pz -5 $ bottom of cylinder at -5 cm c

15 px +10 $ +ve x extent of water 20 px -10 $ -ve x extent of water 25 py +10 $ +ve y extent of water 30 py -10 $ -ve y extent of water 35 pz +10 $ +ve z extent of water 40 pz -10 $ -ve z extent of water c

c *** Source Card*** c

sdef $ default source c

c *** Material Cards *** c

m1 92235 -94.73 92238 -5.27 $ U235/238 m2 1001 +1 8016 +2 $ water c

c *** Rotation matrix for cylinder *** c

c o1 o2 o3 b1 b2 b3 b4 b5 b6 b7 b8 b9 M c x y z xx' yx' zx' xy' yy' zy' xz' yz' zz' +1|-1

c Tr1 0 0 0 1 0 0 0 .7071 -.7071 0 .7071 .7071 1 $ cosines of angles *tr1 0 0 0 0 90 90 90 315 135 90 45 45 1 $ angles in degrees c

c *** Print first 50 histories *** print 110 c

c *** time and nps cut cards *** nps 60 ctme 1

End of File

例四：

message: datapath=D:\\transport\\data\\mcnpxs

title -- cylinder along z axis, in box using surfaces --- c c

c *** cell descriptions ***

c

10 1 -18.72 -1 -5 +10 imp:n=1 $ cylinder c

c Use of the not operator - use with a cell ie #10 means not cell 10 c 15 2 -1 -15 +20 -25 +30 -35 +40 #10 imp:n=1 $ H2O but not cylinder c

c Expanded deffinition use of : opperator, use of paranthesis and continuation

15 2 -1 (+1 -15 20 -25 30 -35 40):(-15 20 -25 30 -35 5) :(-15 20 -25 30 40 -10)

imp:n=1 $ H2O but not cylinder c

20 3 -2.7 -45 +50 +60 -55 -65 +35 imp:n=1 $ Al top 25 3 -2.7 -45 +50 +60 -55 -40 +70 imp:n=1 $ Al bottom 30 3 -2.7 -15 +20 +60 -30 -35 +40 imp:n=1 $ -y side 35 3 -2.7 -15 +20 +25 -55 -35 +40 imp:n=1 $ +y side 40 3 -2.7 -45 +15 +60 -55 -35 +40 imp:n=1 $ +x side 45 3 -2.7 -20 +50 +60 -55 -35 +40 imp:n=1 $ -x side 50 0 +45:-50:+55:-60:+65:-70 imp:n=0 $ outside world c

1 cz 5 $ infinite cyl along z-axis rad 5 cm c

5 pz +5 $ top of cylinder at 5 cm

10 pz -5 $ bottom of cylinder at -5 cm c

15 px +10 $ +ve x extent of warer 20 px -10 $ -ve x extent of water 25 py +10 $ +ve y extent of water 30 py -10 $ -ve y extent of water 35 pz +10 $ +ve z extent of water 40 pz -10 $ -ve z extent of water 45 px +12 $ +ve x extent of Al 50 px -12 $ -ve x extent of Al 55 py +12 $ +ve y extent of Al 60 py -12 $ -ve y extent of Al 65 pz +12 $ +ve z extent of Al 70 pz -12 $ -ve z extent of Al c

c *** Source Card*** c

sdef $ default source c

c *** Material Cards *** c

m1 92235 -94.73 92238 -5.27 $ U235/238 m2 1001 +1 8016 +2 $ water m3 13027 +1 $ Al c

c *** Print first 50 histories *** print 110 c

c *** time and nps cut cards *** nps 60 ctme 1

End of File

这一讲，就说到这里。

欲知后事如何，请听下一讲分解。*^_^*

第09讲，MCNP 几何部分 宏观物体定义曲面 练习题

MCNP 几何部分 宏观物体定义曲面 练习题 Geometry problems II – Macrobodies.

Repeat the first set of exercises - You should use Macrobodies in each case.

In each case you should try to predict how MCNP will assign each facet number - then check this using the plotter.

PROBLEMS:

5. A cylinder lying along the z-axis. Radius 5 cm height 10 cm. Make the cylinder from material U235 weight fraction 94.73, U238 weight fraction 5.27, material density 18.7g/cm3

6. The cylinder in 1. Surrounded by a 20 cm side length cube of water. Position the cylinder centrally in the water.

7. The geometry of problem 2, but with the cylinder rotated about the x-axis by 45 degrees.

8. An Al box, 2 cm wall thickness and internal side 20 cm. The box has a top and a bottom. Inside of the box place a 5 cm radius cylinder of height 10 cm laying long the z-axis. The cylinder should be central within the box. The space between the cylinder and box is filled with water. The Al walls are of density 2.7 g/cm3. The cylinder has the material properties described in problem 1.

Extension - Try problem 4 using a mixture of cells, surfaces and macrobodies.

练习题答案： 例五：

message: datapath=D:\\transport\\data\\mcnpxs

title -- P5 cylinder along z axis using MBs --- c c

c *** cell descriptions *** c

10 1 -1 -1 imp:n=1 $ cylinder c

20 0 +1 imp:n=0 $ outside world c

1 rcc 0 0 -5 0 0 10 5 $ Right circ. cyl along z-axis rad 5 cm, height 10 cm c

c *** Source Card*** c

sdef $ default source c

c *** Material Cards *** c

m1 92235 -94.73 92238 -5.27 $ U235/238 c

c *** Print first 50 histories *** print 110 c

c *** time and nps cut cards ***

nps 60 ctme 1

End of File

例六：

message: datapath=D:\\transport\\data\\mcnpxs

title -- P6 cylinder along z axis in H2O using MBs --- c c

c *** cell descriptions *** c

10 1 -1 -1 imp:n=1 $ cylinder c

20 2 -1 1 -2 imp:n=1 $ water c

30 0 +2 imp:n=0 $ outside world c

1 rcc 0 0 -5 0 0 10 5 $ Right circ. cyl along z-axis rad 1 cm, height 10 cm c

2 box +10 -10 -10 0 20 0 0 0 20 -20 0 0

c *** Source Card*** c

sdef $ default source c

c *** Material Cards *** c

m1 92235 -94.73 92238 -5.27 $ U235/238 m2 1001 +1 8016 +2 $ water c

c *** Print first 50 histories *** print 110 c

c *** time and nps cut cards *** nps 60 ctme 1

End of File

例七：

message: datapath=D:\\transport\\data\\mcnpxs

title --P7 cylinder in z-y plane rot. about x in H2O using Mbs ---MBs c c

c *** cell descriptions *** c

c ***************** Cylinder***** c

10 1 -6.45 -1 imp:n=1 $ cylinder c

20 2 -1 +1 -2 imp:n=1 $ water c

30 0 +2 imp:n=0 $ outside world c

1 rcc 0 -3.5355 -3.5355 0 7.0710 7.0710 5 c

2 box 10 -10 -10 0 20 0 0 0 20 -20 0 0

c *** Source Card*** c

sdef $ default source c

c *** Material Cards *** c

m1 92235 -94.73 92238 -5.27 $ U235/238 m2 1001 +1 8016 +2 $ water print 110 c

c *** time and nps cut cards *** nps 60 ctme 1

End of File

例八：

message: datapath=D:\\transport\\data\\mcnpxs

title -- P8 cylinder along z axis in box using MBs --- c c

c *** cell descriptions *** c

10 1 -1 -1 imp:n=1 $ cylinder c

20 2 -1 1 -2 imp:n=1 $ water c

25 3 -2.7 +2 -3 imp:n=1 $ Al box c

30 0 +3 imp:n=0 $ outside world c

1 rcc 0 0 -5 0 0 10 5 $ Right circ. cyl along z-axis rad 1 cm, height 10 cm c

2 box +10 -10 -10 0 20 0 0 0 20 -20 0 0 c

3 box +12 -12 -12 0 24 0 0 0 24 -24 0 0

c *** Source Card*** c

sdef $ default source c

c *** Material Cards *** c

m1 92235 -94.73 92238 -5.27 $ U235/238 m2 1001 +1 8016 +2 $ water m3 13027 1 $ Al c

c *** Print first 50 histories *** print 110 c

c *** time and nps cut cards *** nps 60 ctme 1

End of File

这一讲，就说到这里。

欲知后事如何，请听下一讲分解：MCNP源Source的描述。*^_^*

第10讲，MCNP源Source的描述

本文来源：https://www.bwwdw.com/article/k3k6.html