目录
源卡(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
运算部分结果如下
图形
(立方体和圆柱体有厚度,所以为两层)
现象
隐藏最外层
两个外曾立方体都隐藏
我们可以看到圆柱体也是双层的