MCNP一个简单案例的应用

本次计算主要是用Fmesh卡来计算钴-60源的剂量当量的计算过程，要求如下：

一个钴60圆柱体放射源，活度是10Ci，直径是5cm，高度是10cm. 2.外面有5层各2cm厚度的铝、铜、铁、铅和钨，是圆柱桶型状的。从里面往外尺寸分别如下所示： 第一层包裹圆柱是铝的内直径是5cm，内高度是10cm；外直径是9cm，外高度是14cm; 第二层包裹圆柱的铜的内直径是9cm，内高度是14cm；外直径是13cm, 外高度是18cm； 第三层包裹圆柱的铜的内直径是13cm，内高度是18cm；外直径是17cm, 外高度是22cm； 第四层包裹圆柱的铜的内直径是17cm，内高度是22cm；外直径是21cm, 外高度是26cm； 第五层包裹圆柱的铜的内直径是21cm，内高度是26cm；外直径是25cm, 外高度是30cm； 3.计算从外表面0-1cm处的剂量当量率的云图。分别使用Tmesh卡和Fmesh卡来计算统计。

代码如下：

C Co60剂量当量计算
1 1 -8.9 -1 7 -8
2 2 -2.72 -2 9 -10#1
3 3 -8.902 -3 11 -12#1#2
4 4 -7.874 -4 13 -14#1#2#3
5 5 -11.34 -5 15 -16#(-4 13 -14)
6 6 -19.26 -6 17 -18#(-5 15 -16)
7 0 -19#(-6 17 -18)
8 0 19

1  cz  2.5
2  cz  4.5
3  cz  6.5
4  cz  8.5
5  cz  10.5
6  cz  12.5
7  pz  -5
8  pz   5
9  pz  -7
10 pz   7
11 pz  -8
12 pz   8
13 pz  -10
14 pz   10
15 pz  -12
16 pz  12 
17 pz  -14
18 pz  14
19 RPP -150 150 -150 150 -150 150

mode p
imp:p 1 6r 0
sdef erg=d1 pos=0 0 0 RAD=d2 EXT=d3 CEL=1 AXS=0 0 1 WGT=3.7E11
si1 L 1.17 1.33
sp1 1 1
si2  2.5
sp2 -21 1
si3 5
sp3 -21 0
m1 27060 1 $Co60 8.9g/cm3
m2 13027 1 $Al 2.72g/cm3
m3 29000 1 $Cu 8.902g/cm3
m4 26056 1 $Fe 7.874g/cm3
m5 82000 1 $Pb 11.34g/cm3
m6 74000 1 $Wu  19.26g/cm3
FMESH4:p GEOM=CYL ORIGIN=0 0 -15 AXS=0 0 1 vec=0 0 1
         imesh=100 iints=100
         jmesh=30  jints=60
         kmesh=1   kints=1 out=ij
DE4  0.01 0.03 0.05 0.07 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55
     0.6 0.65 0.7 0.8 1.0 1.4 1.8
DF4  3.96E-6 5.82E-7 2.9E-7 2.58E-7 2.83E-7 3.79E-7 5.01E-7 6.31E-7
     7.59E-7 8.78E-7 9.85E-7 1.08E-6 1.17E-6 1.27E-6 1.36E-6 1.44E-6
     1.52E-6 1.68E-6 1.98E-6 2.51E-6 2.99E-6
nps 1E7