最新MCNP计数心得.docx

《最新MCNP计数心得.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、精品资料MCNP计数心得........................................精品资料MCNP计数心得   一.计数器的选择及其意义F2面通量密度F4体通量密度F7栅元内的裂变沉积能1.TallyF4对于各种计数器相应的SD卡的问题：F4:n12sd1注：cell12的实际体积为1.55345E+07cc。并且该体积是mcnp能够自己计算出来的话，可以不用sd卡。1)       不用sd卡得出的结果为2.76091E-061/cm^2，表示的是在cell12内，1个源中子的径迹长度总和在体积上的平均值，即可以理解为体通量密度；

2、2)       假设用了sd卡，且为sd1，则计数结果为4.28893E+01，表示在cell12的全部体积内，1个源中子的径迹长度总和，即可以理解为体通量=2.76091E-06*1.55345E+07。（注：因此我们可以知道一个一般性的规律，即涉及到“密度”字样的一般都是指单位体积；涉及到“率”字样的一般都指单位时间。）3)       假设cell12的体积不能被mcnp自己计算出（主要原因是在使用U卡和fill卡，一般都是使用的无限长度的空间与某一曲面截断从而封闭的，因此mcnp不能计算体积；还有就是某些非旋转体，mcnp也不能计算体积），此时

3、需要使用sd卡。如果不使用sd卡，运行程序时会报错，“所计数栅元的体积不能计算或没有给出”。这时，一般来说，要计算通量密度的话，都需要在sd卡上给出栅元的实际体积。若给出的是sd1，意义与上2）中分析的相同，具体实例如下所示：f14:n(120<13[000])(20<13[100]) (110<13[200])(70<13[300])    (80<13[400]) (150<13[500])(100<13[600])    (20<13[010]) (60<13[110]) (70<13[210]) (80<13[310])    (70<13[41

4、0]) (90<13[510]) (100<13[610])    (110<13[020])(70<13[120]) (120<13[220])(70<13[320])    (80<13[420]) (130<13[520])    (70<13[030]) (80<13[130]) (70<13[230]) (60<13[330])    (140<13[430])(100<13[530])........................................精品资料    (80<13[040]) (70<13[140]) (80<13

5、[240]) (140<13[340])    (100<13[440])    (150<13[050])(90<13[150]) (130<13[250])(100<13[350])    (100<13[060])(100<13[160])$组件通量计数fm14(8.78554E+19)fc14fluxtallyineacharraiesofsoutheastcore,totally37.sd141.27106E+0536r$thevolumeofeachfuelarray,andtheyareidentical.注：各组件体积不能由mcnp自己计

6、算出，因此需要在sd卡中给出。此处使用的fm乘子卡，只是乘以一个经过计算出来的归一化源强常数：在某一稳态功率水平下（此处为1035MWt），考虑一次裂变释放的热量1fission放出的热量为180.912MeV（该值也是经过计算得出的，mcnp中提供该值的计算方法，只需要应用相应的反应截面就行），1MeV=1.602e-13joules，1Watt=1joules/second因此，要产生1joule的热量，需要的裂变次数为： 2.TallyF7f7:n12........................................精品资料sd1fm

7、7(1.40744e1)此计数得出的是堆芯的平均裂变沉积能量。其中的fm乘子卡上的数据为，单位转换与归一化源强常数的乘积=1.602e-13joules/MeV*8.78554e19nutrons/sec*1e-6MW/W=1.40744e1。这是由于f7得出的计数值单位为MeV/gram，当使用sd1时，计算cell12中所有质量的沉积能总和，单位变为MeV。并且，计算能量时，一般我们都使用sd1的格式，而不关心单位质量的能量值。而如果要在sd卡上输入值的话，应该写该cell的质量而不是体积了。二、En计数能量卡1.若不使用该卡，则被计数的粒子的整个能

8、量范围将视为一个箱处理；2.若使用该卡，但是为e0卡的话，则建立一个新的能量箱结