蒙特卡罗MCNP学习汇总(三)----源项基础（数据卡\源卡编写）

目录

源卡（C data card）

MODE

SDEF

代码

讲解

完整例子

 图形

现象

---

源卡（C data card）

我们经过前面两篇文章的学习，开以构建简单的几何体了，那我我们现在就要学习编写数据卡，来让整个程序变得完整，而不仅仅是单纯的构建一个几何体。

MODE

问题类型卡（联合输运的粒子类型） 格式：

MODE X X= N   中子输运

MODE X X= P 光子输运

MODE X X=E     电子输运

以上几种粒子的联合输运 输入文件中定义类型的种类—那种源粒子将被输运 每个输入文件都必须在数据卡部分有MODE卡 如果MODE卡被省略，则MODE N默认值被使用

SDEF

通用源 SDEF （确定源）

临界源 KCODE/KSRC （裂变源）

用户所需输入项  Position:粒子发射的初始位置 Time: 粒子发射的时间 Energy:粒子的能量 Direction:粒子发射的方向 Weight:粒子的权重 Type: 粒子的类型

代码

C data card
MODE P E
imp:p 1 1 0
imp:e 1 1 0
M1 11000 1 53000 1
M2 82000 1
SDEF POS= 0 0 1 ERG=0.661 PAR=2
F28:P,E 101
E28 0 0.001 510I 1
NPS 10000
PRINT

讲解

体卡 面卡 省略 见前文

C data card $数据卡 源卡
MODE P E $ 定义 P E 粒子类型
imp:p 1 1 0 $ 上面cell卡定义的体对e的光子权重 （上面如果写了可以不写，写在上面或者下面）
imp:e 1 1 0 $   上面cell卡定义的体对e的权重
M1 11000 1 53000 1 $1号物质(对应上面体卡提到的物质数) 前面的是原子序数 这个对应的的就是NaI，11是钠 53是I 000代表的自然界中的比重 （也可以写成 11023 后面代表的就是相对分子质量）1 是比重 1:1 Na是1 I是1所以组成了NaI
M2 82000 1 $ 铀
SDEF POS= 0 0 1 ERG=0.661 PAR=2$源 位置（0，0，1） 能量0.661Mev（137铯源对应的能量） PAR是发射粒子类型 类型为2 
F28:P,E 101$探测器F8（后面具体说明）
E28 0 0.001 510I 1$能量刻度 0 0.001 中间510个均等分布的能量间隔 到 1 
NPS 10000$10000个粒子数
PRINT$打印出来

完整例子

建立一个立方体 并且里面包含六个圆柱 并书写源卡得到数据

C Cell card
101 1 -10.96 -101
102 1 -2.7 101 -102
103 1 -10.96 -103
104 1 -2.7 103 -104
105 1 -10.96 -105
106 1 -2.7 105 -106
107 1 -10.96 -107
108 1 -2.7 107 -108
109 1 -10.96 -109
110 1 -2.7 109 -110
111 1 -10.96 -111
112 1 -2.7 111 -112
113 1 -1 102 104 106 108 110  112 -113
114 1 -1.15 113 -114
115 0 114 

C Surface card
101 RCC 6.5 6 5.5 0 0 14 2.5
102 RCC 6.5 6 5 0 0 15 3 
103 RCC 6.5 15 5.5 0 0 14 2.5
104 RCC 6.5 15 5 0 0 15 3 
105 RCC 6.5 24 5.5 0 0 14 2.5
106 RCC 6.5 24 5 0 0 15 3 
107 RCC 15.5 6 5.5 0 0 14 2.5
108 RCC 15.5 6 5 0 0 15 3 
109 RCC 15.5 15 5.5 0 0 14 2.5
110 RCC 15.5 15 5 0 0 15 3 
111 RCC 15.5 24 5.5 0 0 14 2.5
112 RCC 15.5 24 5 0 0 15 3 
113 RPP  0.5 21.5 0.5 29.5 0.5 24.5
114 RPP  0 22 0 30 0 25

C data card
MODE P E
imp:p 1 13r 0$体卡数量较多且权重相同 可以用nr代替 n为数量 
imp:e 1 13r 0
M1 11000 1 53000 1
SDEF POS= 0 0 -1 ERG=0.661 
F28:P,E 101
E28 0 0.001 510I 1
F18:P,E 101
E18 0 0.001 510I 1
NPS 1000000
PRINT

运算部分结果如下

 图形

 （立方体和圆柱体有厚度，所以为两层）

现象

 隐藏最外层

 两个外曾立方体都隐藏

 我们可以看到圆柱体也是双层的