蒙特卡罗MCNP学习汇总(二) |
您所在的位置:网站首页 › 几何的意思 › 蒙特卡罗MCNP学习汇总(二) |
几何体定义
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 就可以看到透明的了 我们可以看到就是大球包裹着一个圆柱 与我们编写的程序想要达到的实验效果一致。 |
今日新闻 |
推荐新闻 |
CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |