蒙特卡罗MCNP学习汇总（二）--几何体基础，构建几何体

目录

MCNP几何体基本概念

详细介绍

常用几种面

程序

讲解

现象

---

MCNP几何体基本概念

  空间可以划分为不同材料和性质的区域

几何模型必须包括一个无穷空间的任何部分

基本的几何体单元为 Cell

所有的Cell 由闭合的面 ( Surfaces )组成

所有的面将空间划分为两个区域

— “ left ”VS “right”

— “ inside ” VS “outside

几何体定义

MCNP基于面（ Surface ）的组合几何

 定义面 ( Surface )

由面（ Surface ）和逻辑操作定义栅元 ( Cell )

可以由一组栅元(Cell)定义为一个宇宙（ Universe ），可

以将这个宇宙（ Universe ）装入其它栅元 (Cell) 中 [U卡，

FILL卡]

可以由一组栅元(Cell) 定义为一个宇宙（ Universe），重

复这个宇宙（ Universe ）形成格子排列（ lattice ） ,然后

将这个宇宙(Universe) 装入到一个栅元中(Cell)

将 Material 与栅元 Cell相对应

将其它的 属性 与栅元 Cell 相对应（ e.g., importance weights）

采用栅元Cell 或者面 surface 的编号定义计数器（ Tallies ） 

详细介绍

我们在surface card中定义我们需要用到的面

常用几种面

1.PZ 与Z轴垂直的面  例如 PZ 10 意思是在Z轴坐标为10的与Z轴垂直的面

2.PY 与Z轴垂直的面  例如 PY -10 意思是在Z轴坐标为-10的与Y轴垂直的面

3.PX 与Z轴垂直的面  例如 PX 2 意思是在Z轴坐标为2的与X轴垂直的面

4.SO 圆心在原点的圆面 例如 SO 20 意思是 圆心在原点半径为20的圆（cm）

5.C/X C/Y C/Z 圆柱面 分别与X Y Z 轴平行的圆柱面 例如 C/Z 5 5 5 X Y 坐标为 5 5 的半径为 5 的圆柱面 Z轴上无限长 不封底面

6.RCC 圆柱面 例如 0 0 2 0 0 5 3 底面圆心在（0，0，2）高度为5 Z轴方向的圆柱体 4 5 6 三个数字表达向量 指向圆柱朝向

7.RPP 长方体面 例如 2 4 6 7 1 5  X从2到4 Y 从 6到7 Z 从 1到5

利用面构成体

若在面内 则为- 若在面外则为+

若在面左则为- 若在面外则为+

程序

编写一个大球包裹住小圆柱的程序

This is the input file for the Test1
C cell card
101 1 -17.18 -101 102 -103 IMP:N=1
102 2 -1.0 #101 -104 IMP:N=1
103 0 104 IMP:N=0

C surface card
101 C/Z 4 6 5
102 PZ 0
103 PZ 10
104 S 4 6 5 30

C data card
MODE N P
M1 92235 -94.73 92238 -5.27
M2 1001.60c 2 8016.60c 1
KCODE 1000 1.0 15 115
KSRC 4 6 5
PRINT

注：imp:n=1为对中子的权重 因为本软件最终还是要达到仿真运算的目的 所以要尽量避免不必要的计算 所以在我们需要得到数据的地方设置权重为1 在不需要的地方 可以设置为真空 权重为零

讲解

前面文章讲解过的地方 我就尽量不重复了 不懂的可以复习前面

This is the input file for the Test1
C cell card
101 1 -17.18 -101 102 -103 IMP:N=1 $编号101 1号物质 密度17.18 101面里面交102面外面交103面里面 对中子权重为1 （先看下面面的定义，为Z轴朝向圆柱面里交两个垂直Z轴的面，围成了一个圆柱体）

102 2 -1.0 #101 -104 IMP:N=1$编号102 2号物质 （#非 对体的非 除了这个体以外的地方）101体外交104面里 对中子权重为1 （对圆柱体外与圆面里交的空间进行定义）
103 0 104 IMP:N=0$编号103 104面外面 真空 对中子权重为零（定义圆的外面）

C surface card
101 C/Z 4 6 5$圆心为（4，6）半径为5 与Z周平行的圆柱面 不带底面
102 PZ 0$Z轴上坐标为零 与Z轴垂直的面
103 PZ 10$Z轴上坐标为10 与Z轴垂直的面
104 S 4 6 5 30￥圆心为（4，6，5）半径为三十的圆面

C data card（以后再介绍）
MODE N P
M1 92235 -94.73 92238 -5.27
M2 1001.60c 2 8016.60c 1
KCODE 1000 1.0 15 115
KSRC 4 6 5
PRINT

现象

 右键几何体 make all transport 就可以看到透明的了 我们可以看到就是大球包裹着一个圆柱

与我们编写的程序想要达到的实验效果一致。