MCNP6中Fmesh卡使用简单感想

Fmesh 重叠网格记数

 

      Femsh是选择性使用卡，不是一定要使用卡。FMESH卡允许用户定义几何重叠的网格记数。结果写在不同的输出文件，文件的缺省名为MESHTAL。缺省的话，网格记数算出粒子通量在网格栅元的平均轨迹长度估计，单位SI 。如果FMESH前加“*”，为能量时间粒子权重记数，单位Mev/cm2。

   使用格式如下所示：   

   格式：FMESHn:  pl  fmesh 变量=指定

   其中  n =记数号(只有记数类型4可省)   pl = N 或P或 E

 等号可省略。关键词顺序任意。特殊输入I,M,R不能和GEOM、FACTOR、OUT、TR一起使用。表4-3总结了网格记数的关键词，列出了它们的缺省值。默认网格是直角坐标系，原点（0,0,0），OUT为列格式。对于圆柱网格，默认圆柱轴平行于Z轴，定义半平面θ=0为x轴正方向。网格记数可用于DE、DE和FM的组合。

输出格式由关键词OUT指定。默认格式为”col”，按列格式列出箱坐标、记数结果和相对误差。如果OUT设为”cf”，除了输出坐标和记数结果，还有网格体积和相应的记数时间。第二种输出格式由”ij”,”ik”,”jk”指定，记数结果按一系列二维矩阵输出，其中在不同坐标系i=x或r，j=y或z，k=z或θ。第一个矩阵是记数结果，第二个矩阵是相对误差。行和列用相应网格的箱的中点标记。这些对矩阵在第三维坐标打印出来。

例：       FMESH4:n           GEOM=cyl           ORIGIN= –100 0 0

IMESH=5 10 IINTS=5 2

JMESH= 100 200 JINTS 10 5

KMESH .5 1 KINTS=1 2

AXIS= 1 0 0 VEC=0 1 0  OUT=ij

这是一个沿x轴的圆柱网格记数，底面在y轴正方向，且x = −100，θ=0。记数从r=0到r=5分成5箱，从r=5到r=10分成2箱；从x = -100 到x = 0分成10箱，从x = 0 到x = 100分成5箱，θ=0°到θ=180°为1箱，从θ=180°到θ=360°分成2箱。

例如：fmesh104:n   origin=-260 -260 -260

             imesh=260  iints=52

             jmesh=260  jints=52

             kmesh=260  kints=52 

fmesh114:p   origin=-260 -260 -260

             imesh=260  iints=52

             jmesh=260  jints=52

             kmesh=260  kints=52 

fmesh124:e   origin=-260 -260 -260

             imesh=260  iints=52

             jmesh=260  jints=52

             kmesh=260  kints=52 

fmesh134:h   origin=-260 -260 -260

             imesh=260  iints=52

             jmesh=260  jints=52

             kmesh=260  kints=52 

此程序的origin就是它的起始坐标，本次从（-260，-260，-260）开始，然后x轴是到末尾是260，iints=52就是x轴分成52份；y轴的末尾坐标是260，jints=52 z轴的末尾坐标是260。如果程序写了out=ij ik jk 就会输出例如x,y平面，z轴高度是定值。

如下图所示，x,y,z坐标就是这个立方块的中心坐标。

处理数据：可以通过matlab进行处理，或者通过origin进行处理，或者把上图的数据导入到Excel中，然后通过筛选把z轴值固定住，取x,y面的数据。