目录
MCNP几何体基本概念
详细介绍
我们在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 就可以看到透明的了 我们可以看到就是大球包裹着一个圆柱
与我们编写的程序想要达到的实验效果一致。