蒙特卡罗方法简介 联系客服

发布时间 : 星期日 文章蒙特卡罗方法简介更新完毕开始阅读2bd812d728ea81c758f5785a

j:栅元的题目序号,写在第 1-5 列,1≦j≦99999。 “题目序号”是指由用户按任意所希望的次序选定的序号,在一个文件中,可以不连续。然而 MCNP 将按照读入的顺序对栅元进行编号,从 1 开始按增量为 1 单调上升。我们将此称之为栅元的“程序编号” 。必须明白题目编号与程序编号的区别,因为后面叙述的一些卡片将会引用到这种或那种编号。不过多数情况下, 用户给出的序号往往与程序编号一致, 这也是避免混乱, 减少错误的最可靠方法。 m:栅元的材料号,它是材料卡(Mm)中相应材料的序号。真空栅元,m=0。 d:栅元材料的密度。填入正值时,表示是原子密度(单位为 1024个原子/cm3);填入负值时,则是质量密度(单位为 gm/cm3)。对于真空栅元,该项不填,直接列出下一项内容。

geom:栅元的几何描述。它列出界定该栅元的所有曲面号(带有数将,表示坐向) ,及描述这些曲面所定义的曲面之间关系的布尔算符,布尔算将包括交

(Intersection) ,用一个空格表示;联(Union) ,用“: ”号表示;余(Complement) ,用“#”表示。这里所说的曲面号及栅元号都指的是题目编号。如果定义该栅元要用到一个虚设曲面,则它必须作为界定该概元的一个曲面列出(同样带有坐向数符) 。

param:任选的栅元参数说明,其形式为:关键词=数值。

3.3.4 源的描述卡

对每一个 MCNP 问题都有四个源中的一个。四个源包括:①通用源(SDEF 卡);②曲面源(SSR 卡);③临界源(KCODE 卡);④用户提供的源。除了临界源之外,其它源都能使用源分布函数,分布函数分别由 SIn、SPn、SBn 和 DSn 卡指定。

下列卡片用于指定源的信息: 助记名(1-5列)

SDEF SIn SPn SBn

卡片内容 通用源 源信息 源概率 源偏倚

助记名(1-5列)

SCn SSW SSR KCODE

卡片内容 源注释 写曲面源 读曲面源 迭代源

DSn 相关源 KSRC 源的一些特点

在一些情况下,MODE 卡也作为源描述部分来说明从源出发粒子的隐含类型。 对源粒子的描述应包括如下信息: i)几何位置(xxx,yyy,zzz); ii)飞行方向(uuu,vvv,www) iii)能量(ERG); iv)开始时间(TME); v)粒子类型(IPT); vi)粒子统计权重(WGT); vii)粒子开始栅元(ICL);

viii)粒子开始曲面(JSU),如果开始点不在曲面上该值为0。

如果对点探测器或者 DXTRAN球的问题还必须定义另外一些变量。源粒子的类型(中子或光子)则自动由 MODE 卡指定了。上述变量的任何一个都可能有一个概率分布。源粒子若出生在重要性为“0”的栅元,则是一个致命性错误。

3.3.5 FSn分段记数卡(记数类型 1,2,4,6,7)

格式: FSn S1 S2 ? Sk n:记数号;

Si:带数符的分段曲面的题目编号。

此卡不用于探测器。可以要求 Sdn 卡。该卡允许用户把一个栅元或曲面分成若干段来记数,其优点在于不必专为记数设置额外的栅元。在 FSn 卡上指定的分段曲面应与正规的问题曲面一起列出曲面描述卡,但它们不一定是实际几何的组成部分。因此无需按栅元与曲面间的关系把它们考虑得很复杂。

如果在FSn卡上填写k个曲面,就把曲面或栅元分成了k+1段。按照在 FSn卡上列出的分段曲面的次序和坐向把记数n分成k+l段记数箱。采用分段记数时,通常只给出每段的结果,当在FSn卡末尾加上“T”时,则增加一个总的记数箱,同时给出整个曲面或栅元上的记数。分段是按如下方法做的: Fn:n S(或C) FSn S1 S2 ? Sk

在曲面S(或栅元C)上的记数n将被FSn卡分成k+l个箱: (1)相对曲面S1的坐向与数符相同的部分;

(2)相对曲面S2的坐向与数符相同的部分,但在其他段上记录过的除外;

????

(k)相对曲面Sk的坐向与数符相同的部分,但在其他段上记录过的除外; (k+1)S(或C)上所有其余部分;

(k+2)如果在FSn卡末尾加上字符T,将给出整个曲面或栅元上的记数。

3.3.6 SDn分段除数卡(记数类型 1,2,4,6,7)

格式:SDn (D11 D12 ? D1m)(D21 D22 ? D2m)?(Dk1 Dk2 ? Dkm)

n:记数号;

k:Fn卡上的栅元或曲面个数,包括 T(如果存在);

m:FSn 卡上分割的箱数,包括剩余部分,如果FSn卡上有 T,也包括总的部分;

Dij:供记数n的第i个曲面或栅元的第j段的面积、体积或质量。括号是选择的。如果缺省此卡,若对栅元或曲面不分段,既没有给出 VOL 卡或 AREA卡,并且 MCNP也不能自动计算体积或面积;或者对栅元或曲面分成段,但MCNP也不能自动计算体积或面积;均将产生致命错误(记数类型 1 除外)。

[6]