MCNP计数心得
一.计数器的选择及其意义 F2 面通量密度 F4 体通量密度 F7 栅元内的裂变沉积能 1.Tally F4
对于各种计数器相应的SD卡的问题: F4:n 12 sd 1
注:cell 12 的实际体积为1.55345E+07cc。并且该体积是mcnp能够自己计算出来的话,可以不用sd卡。
1) 不用sd卡得出的结果为2.76091E-06 1/cm^2,表示的是在cell 12内,1个源中子的径迹长度总和在体积上的平均值,即可以理解为体通量密度;
2) 假设用了sd卡,且为sd 1,则计数结果为4.28893E+01,表示在cell 12的全部体积内,1个源中子的径迹长度总和,即可以理解为体通量=2.76091E-06*1.55345E+07。(注:因此我们可以知道一个一般性的规律,即涉及到“密度”字样的一般都是指单位体积;涉及到“率”字样的一般都指单位时间。)
3) 假设cell 12的体积不能被mcnp自己计算出(主要原因是在使用U卡和fill卡,一般都是使用的无限长度的空间与某一曲面截断从而封闭的,因此mcnp不能计算体积;还有就是某些非旋转体,mcnp也不能计算体积),此时需要使用sd卡。如果不使用sd卡,运行程序时会报错,“所计数栅元的体积不能计算或没有给出”。这时,一般来说,要计算通量密度的话,都需要在sd卡上给出栅元的实际体积。若给出的是sd 1,意义与上2)中分析的相同,具体实例如下所示:
f14:n (120<13[0 0 0]) (20<13[1 0 0]) (110<13[2 0 0]) (70&l