SP10 -21 0 $ axial sampling weight for src2(2号源轴向抽样权重) c -- set energies of photons for each source(设置每个源的光子能量) DS20 S 21 22
SI21 L 0.6938 1.1732 1.3325 $ Co-60 spectra for src 1(1号源为Co-60能谱) SP21 D 1.6312E-4 1 1 $ frequencies of gammas(伽玛光子频率)
SI22 L 0.667 $ Cs-137 spectrum for src 2(2号源为Cs-137能谱) SP22 D 1
3.4结果说明 在MCNP5的计算中允许的各种类型的输出结果的技术说明,在使用手册的第二章的Section V中给出。
详细指定的各种输出结果,包括结果卡和结果修正卡的细节,将在第三章的Section IV中给出。MCNP5中有效输出结果的摘要如下所示。
使用频率最多的标签是面流量F1,平均面通量F2,点或环通量F5,体通量F4。类似于体通量标签,还有能量沉积的各种标签(F6和F7)。除非使用FM卡特别说明,标签只用于单一粒子源。除了标签F6和F7外,像*F1:P指定一个标签,例如,一个光子能量,使每个事件的输出结果变得复杂。能量通量和能量流
量在输出结果标签中。标签F6和F7已经在能量单元中。
同样可以使用多重标签Fn:pl卡,每一个都有一个唯一的n值。最后的n决定了输出结果的类型。因此,
例如我们可以使用F2:N,F12:P和F22E分别给出中子、光子和电子的平均面通量。
下面介绍各种标签的物理性质。在描述中,假设时间不是独立的,这在MCNP计算仿真中是正常的。通量
经过一段时间就累积,也许称为注量更为恰当。 3.4.1面流标签(F1)
每当粒子穿过特定的面,它的权重就被记录,然后这些权重之和就被F1记录并在MCNP输出文件中显示。注意这没有被区域面积A区分。也没有被截面的方向区分。当应用于几何真空问题时,这个标签对检验能量守恒和粒子数守恒有很好的应用。理论上,如果J(r,E,?)???(r,E,?)是当做位置函数的流量的能量和角度分布,F1标签可有以下得到:
其中,向量n为面在向量rs 外方向的向量。 3.4.2平均面通量结果(F2)
假设一个权重为W的粒子,以与面法向量成?角的方向穿过此面。此粒子在这个面上的通量贡献为W|sec?|/A。贡献的总和在MCNP输出文件中F2标签中给出。 理论上,如果作为位置函数的?(r,E,?)同时与能量和角度相关,则标签F2可以有下士得到:
3.4.3平均体通量(F4)
假设粒子权重为W,能量为E,在体积为V的体元中的分割为长度为T的特定栅元。这部分对栅元通量的贡献为WT/V。贡献的总和在MCNP的输出文件中有标签F4核查。理论上,作为位置函数的?(r,E,?)也受到能量和角度分布的影响,标签F4可以由下式得到:
3.4.4点通量或面通量(F5) 这种类型的标签使用了所谓的方差减小技术,也就是使用“紧邻事件估计”。对于每一个源粒子或是碰撞事件,确定性的评估是通过在点探测器的注量贡献(或是在轴对称问题中的环)。为简化标签的描述,假设模拟只是在同一种煤质中。假设在距离探测器r的地方各向同性源释放了能量为E、权重为W的粒子。应用与点核方法的射线理论,根据下式控制??在点探测器的影响贡献,
其中,?(E)粒子能量相关的线性相互作用系数。注意14?每单位立体角为各向同性点源的角度分布。现在假设在距离点探测器r的地方发生碰撞,到达探测器的过程中,需要考虑散射角?s。这里的E为粒子碰撞后的能量,W为权重。如果?(E,?s)是在角度?s方向上的单位立体角内的线性相关系数,那么
?(E,?s)?(E)是在角度?s方向上的单位立体角的概率。空间按1r2衰减,??在点探测器的影响贡献为
3.4.5标签详细说明卡
至少有一个输出卡是必需的,Fn:pl,其中n是输出的编号(最后一位决定输出卡的类型),pl代表N(中子输出),P(光子输出),N,P是中子光子联合输出,E代表电子输出。下面的标签类型是为输出标签(类型F1或是F2)或是体元(F4类型)指定的栅面。对于类型F5的探测器标签,它跟指定的探测器的位置有关。能量沉积、脉冲高度和其他专门的计算结果在本使用手册中不加讨论。在下面一部分,给出几个关于Fn:pl的例子对各种参数的用法进行展示。 3.4.6栅元和栅面卡 F1:E 1 2 T $ current through a surface(面通量) 这个卡指定电子通过两个栅面的流量,然后T为通过两个面的总量。注意:流量计算没有被面区域所分割。 F2:P 1 (1 2) (2 3 4) T $ fluence averaged over surfaces(在面上的平均影响)