在信息块卡、栅元描述卡及曲面描述卡之后输入的卡片统称为数据卡,它包括相当丰
富的内容,大体分为十二类:
1)问题类型卡; 2)栅元和曲面参数;
3)坐标变换; 4)源的描述;
5)多群特征描述; 6)光子产生偏倚描述;
7)记数指定; 8)材料描述;
9)能量及热中子的处理; 10)截断条件;
11)用户数据组; 12)外围卡。
在接续运行输入文件中只允许有10)、11)两类卡,还有12)类中的某些卡。在一个
输入文件中;同名同号码的数据卡至多使用一次,例如 M1及 M2卡是可以同时接受的,
但不允许同时给出两个 M1卡或两个 MODE卡。
3.l 问题类型(MODE)卡
MCNP有三种不同的运行方式:
i)方式n:中子输运。
ii)方式n p:中子-光子耦合输运。
iii)方式p:光子输运。
MODE卡的第 l-5列填入助记名mode,第6-80列上填入n或 p或 n p,这由采用那种运行方式决定。如果不给 MODE卡,则按缺省运行方式,缺省运行方式为 n。如果选用方式n p运行时,所有源粒子都必须是中子。此时只考虑中子核反应诱发的光子,且不考虑光激中子的产生。由于不是所有核素都配备了光产生截面,因此按方式np运行时,如果用到的某种核素没有光子产生截面,在运行XACT时会给出警告信息。下面给出一个MODE卡的格式:
mode p
3.2 栅元参数及曲面参数卡
此类卡片大部分是定义栅元参数的值, AREA 定义曲面参数的值。 这儿也包括了
VECT、WWE、WWP及WWGE等卡。参数卡包括以下卡片:
以上各项只有重要性卡是必须给出的(使用权窗处理时有例外),其他大多数项是对使用降低方差技巧的选择项。
该类卡片上填写数据的次序应与栅元(或曲面)表相对应,卡片上的第n个数据项将是赋给程序编号(描述卡在INP文件出现的顺序)为n的栅元(或曲面)的值。每项栅元参量卡上的数据项数应当等于题目中的栅元数(或曲面数,对AREA卡)。
imp:n x1 x2 … xi … xI
栅元的重要性有两个用途:
a)当重要性为0时,用它来终止粒子的历史;
b)用于几何分裂与轮盘赌技巧,促迸粒子轨迹朝向几何上较为重要的区域运动。
在中子-光子耦合输运问题中,用户可能希望对光子的分裂与中子有所不同,所以可能需要用到两组重要性(中子和光子)。如果在运行方式为n p的问题中,未给出 imp:p卡,那么所有栅元的光子重要性均置为1。
当IMP卡上的数据项数不等于问题中的栅元数时,MCNP将指出为致命性错误。
对IMP卡,nR及xM书写功能是很有用的,例如给出一张IMP卡为:
imp:n l 2 2m 0 1 20r
则使栅元 1的重要性为 l,2为2,3为 4,4为 0,而 5-25,则都是 l。此时若有一条轨迹从块2走进块3,则将被分裂成两条轨迹,每条新轨迹的权是分裂前原轨迹权的一半。相反若轨迹是从块3走进块2,则将以50%的概率终止这条轨迹;如果未被终止,则该轨迹的权将加倍。
应记住,凡进入“0”重要性栅元的轨迹都将被杀死。对探测器进行记数或向DXTRAN球上生成假粒子时,如果轨迹经过一个“0”重要性栅元,则也将被杀死。
当粒子轨迹进入真空栅元时,即使重要性增大,也不进行分裂;但对离开真空块的轨迹,若进人了较低重要性的栅元,则仍将常规实行轮盘赌。
格式: vol x1 x2 … xi … xn
或者: vol no
xi = 相应栅元的体积;
no = 不计算体积或面积。
VOL卡是输入栅元体积的可选择卡。如果没有此卡,MCNP将计算全部栅元的体积。当 VOL卡上填写 no时,表示不计算体积。对于在 VOL卡上没有填写的栅元体积 MCNP
将使用计算的体积。若不使用VOL卡,而栅元的体积不能由MCNP计算,该体积又在计数中要被运用时,程序给出致命性错误信息。可使用nj功能跳过不想填入体积的栅元。栅元的体积按cm3单位填入。
凡是具有任意对称轴(即使轴是斜置的)的旋转栅元,以及多面体,MCNP都将自动算出它们的体积。在计算体积的同时,也可以计算面积和质量。在输出文件 OUTP中输出这些体积、面积和质量。用户可用VOL卡填写的体积代替计算的体积。
如果需要指定某一栅元的片段体积,则应使用SDn卡。VOL卡涉及到几何系统中的所有栅元,而SDn卡只涉及因特殊记数要求而需切割的栅元,给出栅元及分割部分的体积及质量。
格式: area x1 x2 … xi … xn
xi = 曲面i的面积;
它类似与VOL卡,用以补充MCNP计算曲面面积的不足。该卡是输入曲面面积的选择卡。如果**NP能够计算出一个曲面的任何一侧的所有栅元的体积,则该曲面的面积也将被计算,否则该曲面的面积将不被计算。对于无法计算的面积,而且又未在AREA卡或SDn上给出,该面积又在计数中要被运用时,程序给出致命性错误信息。
AREA卡将涉及几何描述中的所有曲面,并非只涉及一个特定的记数界面。如果需要一个曲面片断的面积,则需使用FSn卡来切割整个曲面,SDn卡可提供曲面及分割部分曲面的面积。
格式: pwt w1 w2 … wi … wI
Wi : 中子在第i个栅元碰撞所产生的光子权重的值;
I : 问题中的栅元总数。
如果缺省此卡,则所有栅元的 Wi均被置为-l。PWT卡仅用于MODE n p的问题。
对每个栅元i , 可以指定一个最小的光子权Wi (>0),中子诱发的光子将以高于wi*Ii/Is的权出生,其中Is和Ii分别表示中子从第s个栅元出发在第i个栅元碰撞的两个栅元的重要性。当出生光子的权重低于这个值时,就用轮盘赌确定其是否存活。
对于每一个具有负的wi(<0)项的栅元i,仅当中子诱发的光子权重大于- wi*ws*Ii/Is时,则产生光子,其中ws是导致当前光子产生事件中的源中子的初始权,Is和Ii分别表示中子从第s个栅元出发在第i个栅元碰撞的两个栅元的重要性。如果诱发光子的权重低于这个值时,就用轮盘赌确定其是否存活。
当PWT卡上指定的wi为0,那么在光子产生是可能的情况下,则导致中子在i栅元中的每次碰撞都考虑产生光子。想要关掉i块中的光子产生过程时,应置wi =-1.0e6。
PWT卡通过诱发光子总的权重与在PWT卡上指定的有关的权重阈值的比较来控制中子在碰撞时诱发光子的数目。权重阈值与柳元重要性有关,并且如果wi<0,也与源中子权重有关。如果希望有较多的诱发光子,应降低wi的绝对值以减少该权重,这样就会增加光子数。如果希望有较少的诱发光子,应增加wi的绝对值。
下面是关于如何选取PWT卡上数据项的一些要注意的问题:
对与能量无关的光子权窗问题,在 PWT卡上的数据项应与在 WWN1:P卡的数据项相
同,而且诱发光子应尽可能多的在光子权窗界内出生。对与能量有关的光子权窗问题,在PWT卡上的数据项应与在WWN1:P卡上的任何一个权窗能群的最低非零的权窗数据项一样。
对使用光子栅元重要性(imp:p),而不使用权窗(WWNn:P)的问题,PWT卡上的数据项最好为缺省项,即Wi =-1,或者把每个栅元的wi是设置为平均源的权。
当省略了PWT卡时,则意味着所有栅元的最小光子权置为一1.这个缺省值对很多问题都是适用的。
格式: ext:p A1 A2 … Ai …AI
P: n,中子,或p,光子;
Ai : 第i个栅元的数据项,每个Ai的格式均为A=QVm , 其中Q为拉伸的长度,
Vm定义拉伸的方向;
I:问题中的栅元总数。
如果省略卡,则所有栅元的 Ai=0 ,此即非偏移输运的情况。EXT卡是选择卡 , 当有权窗卡时,使用该卡要小心.
指数变换方法是通过调整总截面在所感兴趣的方向上拉伸碰撞之间的径迹长度。总截面调整如下:
其中: ∑t*:调整后的总截面;
∑t :真的总截面;
P:拉伸参数;
μ:粒子方向和拉伸方向角度余弦。
根据拉伸长度Q的三种方式能够指定拉伸参数P,即
Q=0;P=0,不使用指数变换;
Q=p;0<P<l,固定的拉伸参数;
Q=S;P=∑a/∑t ,这里的∑a是俘获截面。
沿飞行方向的隐式俘获使用P=∑a/∑t。
根据EXT卡上的每个Ai项的Vm的三个选择来定义拉伸的方向。
(1)Ai项没有Vm部分,也就是对一个给定的栅元仅填写拉伸度Ai=Q。这样只在粒子的方向上伸长(μ=l),与粒子的方向无关,我们不推荐使用这种方法,除非在沿飞行的径迹上作隐式俘获,在这种情况下人Ai=Q=S,而且抽取的距离是散射距离而不是抽取碰撞的距离。
(2)指定拉伸方向为Vm Vm是从碰撞点到点(Xm,Ym,Zm)的方向,其中(Xm,Ym,Zm)是在VECT卡上指定的。方向μ是粒子方向和拉伸方向Vm之间的夹角余弦,Am的符号给出拉伸的方向是朝向(Xm,Ym,Zm)还是远离(Xm,Ym,Zm)。
(3)指定拉伸方向为Vm=X,Y或Z,那么方向μ分别为粒子方向与X-,Y-或Z-轴之间的夹角余弦,Ai的符号给出拉伸的方向是朝向X,Y或Z轴,还是远离这三个轴。
例子: ext:n 0 0 .7V2 S -SV2 -.6V9 0 .5V9 SZ -.4X
vect V9 0 0 0 V2 1 1 1
卡上的10项是这个问题的10各个栅元的数据项。这个问题对光子和对栅元1,2,7
都不作指数变换。下面给出其它7个栅元的径迹拉伸一览表。
栅元 Ai Q V 拉申参数 方 向
3 .7V12 .7 V2 P=.7 指向点(1,1,1)
4 S S . P=∑a/∑t 粒子方向
5 -SV2 S -V2 P=∑a/∑t 远离点(1,1,1)
6 -.6V9 .6 -V9 P=.6 远离原点
8 .5v9 .5 V9 P=.5 指向原点
9 SZ S Z P=∑a/∑t 沿着Z轴
10 -.4X .4 -X P=.4 沿着-X轴
在强迫碰撞的栅元中或者在没有权窗控制的棚元中不使用指数变换。仅当粒子通量有
一个指数分布时,如强吸收问题,使用指数变换合适。
格式: vect vm xm ym zm … vn xn yn zn …
m,n:标记向量vm,vn任意号,这个号是任意的;
xm,ym;zm:为对应向量vm的定义坐标。
VECT卡的数据项是四个为一组,这些组定义指数变换中的若干个向量或者用户修补的若干个向量。该卡是选择卡,无缺省值。
格式: fcl:n x1 x2 … xi … xI
其中n为n,对于中子,n为p,对于光子;xi为栅元i的强迫碰撞控制,-l<xi<l,缺省值为xi=0,无强迫碰撞; I为问题中的栅元数。
FCL卡指定中子或光子在栅元内的强迫碰撞。该卡对点探测器或DXTRAN球是非常有用的。应用了强迫碰撞功能的栅元不再进行权窗游戏。
如果xi=0,就关掉了对这一栅元的强迫碰撞功能。用适当的权重控制把进入栅元i的全部粒子分为碰撞部分和不碰撞部分。如果|xi|<1以|xi|的存活概率对碰撞部分施行俄罗斯轮盘赌以防止碰撞的历史数变得太大。如果一些强迫碰撞的栅元彼此邻近,则建议xi选为小数。
当不给FCL卡时,各栅元的强迫碰撞数均为“0”。因强迫碰撞花费机时多,强迫碰撞数一般只给成l。
当 xi<0时,强迫碰撞过程仅适用于进入栅元i的粒子。强迫碰撞之后,不考虑权重截断并且用通常的类似方式处理以后的全部碰撞。不忽略权窗并且在对探测器和 DXTRAN球作贡献之后使用这些权窗。
当 xi>0时,强迫碰撞过程适用于进入栅元i的粒子以及权截断或权窗游戏幸存的碰撞粒子。粒子将继续被分为非碰撞和碰撞两部分(概率为|xi|),直到这些粒子被权截断或权窗杀死为止。
习惯应用:
一般情况下令xi =l或xi =-1。当一些强迫碰撞栅元彼此邻近或强迫碰撞粒子产生的数比所要的粒子数高些,应置xi为小数。
当权窗边界把一个进入栅元i的典型权重分类时,选xi>0。当权窗界把强迫碰撞粒子的典型权重分类时,选xi >0。
当应用重要性时,选xi >0,因为xi <0关掉权截断游戏。
在使用指数变换的栅元中,不可做强迫碰撞。
权窗是对指定的空间和能量重要性函数提供重要性(IMP:n卡)和能量分裂(ESPLT:n卡)的另一种方法。权窗的优点是:
(1)权窗提供一个空间和能量的重要性函数;
(2)权窗控制粒子的权重;
(3)与其他一些降低方差技巧不矛盾;
(4)在曲面跨越、碰撞或者这两者都可以使用权窗;
(5)能够控制分裂或者轮盘赌的程度;
(6)在指定的空间或者能量区域上关掉权窗;
(7)用权窗产生器自动产生权窗。
权窗的缺点是:权窗不如重要性简单;当改变源的权重时,权窗必须重新归一。
a)WWE权窗能量卡
格式: wwe:n E1 E2 … Ei … Ej , j≦99
Ei 为能量间隔, Ei-1<E<Ei ,E0=0;
n:n为中子, p为光子。
该卡与WWN卡配合使用,它被用来定义由WWN卡指定的权窗边界的能量间隔。当无WWE卡而有WWN卡时,将根据问题运行的能量范围建立一个能量间隔。最小能量为0,在该卡上不输入。该卡是一个选择卡。
b)WWN权窗下限卡
格式: wwni:n wi1 wi2 … wij … wiJ
n:n为中子,p为光子;
wij为栅元j的权窗下限和由WWE卡定义的能量间隔Ei-1<E<Ei, E0=0。如果没有 WWE卡,i=1
j为问题的栅元数。
当使用重要性卡时,不要求此卡。除此之外,要求权窗(WWN和WWP)卡。
WWN卡指定与空间和能量有关的权窗下限,它必须和WWP卡一起使用,并且如果权窗与能量有关,应有WWE卡。如果使用 WWN: n卡,将不再使用 IMP: n卡,其中n是同一个拉子类型。
如果Wij<0,任何进入栅元j的粒子被杀死,因为负输入相当于0重要性。如果负项用于一个能量区间,那么它们也将用于同一个栅元的其它能量区间。
如果Wij>0,进入栅元j或在栅元内碰撞的粒子将根据WWP卡的选择进行分裂或轮盘赌。
如果 Wij=0,关闭栅元j的能量箱i的权窗游戏同时打开权截断游戏。应在 CUT卡上把权截断指定为这个问题所要求的最低的所允许的权重。否则,进入Wij=0栅元的大部分粒子将被用权截断杀死。当用单个重要性函数或者用一组权窗边界不能表示各种能量及空间区域的特征时,则关掉这些区域的权窗游戏的能力是非常有用的。
对权窗来讲,粒子权的界限总是绝对的,而不是相对的。因此在使用它时,必须明显地顾及到其它一些降低方差技巧引起的权变化,对每一个栅元每一能量间隔必须指定权的下限。在描述时,WWNi卡必须从i=l开始,依序给出,不得有空缺。也就是,若指定WWNi卡,则必须指定WWNi,1<j<i。
例3-l 定义中子问题的三个能量间隔和四个栅元的权窗边界。
wwe:n E1 E2 E3
wwn1:n w11 w12 w13 w14
wwn2:n w21 w22 w23 w24
wwn3:n w31 w32 w33 w34
例3-2 wwn1:p w11 w12 w13
这个例子没有 WWE卡,只用WWN:P卡对一个有三个栅元的问题定义了与能量无关的光子权窗。
c)WWP权窗参数卡
格式: wwp:n WUPN WSURVN MXSPLN MWHEM SWITCHN
n: n为中子, p为光子。
此卡上给出用于权窗处理的其他五个参量。这些量都有缺省值,如果全部使用缺省值,则可不给此卡。此卡上的参量适用于问题中的所有栅元。这些参量依序是:
i)WUPN:当粒子权超过WUPN与权窗下限值之积时,粒子将被分裂。WUPN的缺省值为 5,当用户指定值小于2时,MCNP将其重置为 2。
ii)WSURVN:对于经轮盘赌存活下来的粒子,它的权将成为WSURVN与权窗下限值之积(当WGT与MXSPLN之积不小于该下限值之时)。若指定的WSURVN小于1,MCNP将重置为WUPN,并给出一个警告下限。WSURVN的缺省值为0.6*WUPN。
iii)MXSPLN:粒子分裂时,分裂数不会超过指定的MXSPLN值;做轮盘赌时,赌胜概率不会小于l/MXSPLN。MXSPLN的缺省值为5,MXSPLN>l。
iv)MWHERE:此参数用于决定对粒子权进行核查的地方。 MWHERE的值为- l时,仅在碰撞点处核查;为0时,在界面和碰撞点处检查,为 1 时,仅在界面核查。MWHERE的缺省值为O。
v)SWITCHN:当它的值大于 0时;权窗的下限值将被指定为 SWITCHN除以粒子所在栅元的重要性。SWITCHN的缺省值为0。
WWP卡上的参数是控制WWN卡上指定的权窗下限的使用, 或者如果SWITCHN>0,并且IMP卡存在,此时要有WWNi卡,则出现致命错误。当SWITCHN=0时,权窗下限应由WWNi卡指定。
d)利用现成的栅元重要性来定义权窗的下限值
当采用与能量无关的权窗处理时,则可以利用重要性卡上的数据给出权窗参量,此时
应将WWP卡上的 SWITCHN置一个正的常数 C。在这种处理方式下,栅元的权窗下限值将是C/I,I是该栅元的重要性。WWP卡上的其他项没有改变。C的值应当给得使源粒子的出生权位于权窗之内。若做不到这一点,很可能是权窗太窄,或者是需要重新给出源的描述.
权窗产生器估算由用户指定的与位置、能量有关区域的重要性。然后按与该重要性成反比计算与位置、能量有关的权窗参数。产生器应用于这样一些问题,这个问题没有大角度相关而有为抽样所说明的适当几何。
产生器估算一个相空间栅元的平均重要性。如果栅元太大,则栅元内部的重要性变化也大。这样平均重要性不能代表这个栅元。不适当几何说明也会引起邻近栅元之间的重要性有大的区别。恰好产生器能提供为抽样使用的几何说明是否适当的信息。
a)WWG权窗户生器卡
格式: wwg It Ic Wg xr yr zr J
It:问题的记数号(Fn卡的n)。由TFn卡定义特殊计数箱,对这个特殊计数箱权窗产生器是合适的;
Ic:参照栅元,典型的参照栅元是源栅元;
Wg:产生栅元Ic的权窗下限值。如果Wg=0,它将是平均源权重的一半;
xr,yr,zr:粒子朝向点移动或者离开点的坐标。
J:调试打印标志,J=0/1=off/on
此卡是选择卡。缺省此卡,不产生权窗值。
WWG卡引起在WWP和WWNi卡上产生计数It的最佳重要性函数。在输出文件OUTP中打印评价及综述这些最佳重要性函数,所产生的权窗重要性函数可用于以后的运行中。对许多问题这种产生的重要性函数要比有经验的用户在IMP卡上猜想的要好。为了产生与能量有关的权窗,使用下面描述的WWGE卡。
b)WWGE权窗产生能重卡
格式: wwge:n E1 E2 … Ei … Ej; j≤99
n: n为中子, p为光子;
Ei:产生权窗组的能量上限,Ei+1>Ei。
该卡是选择卡。如果没有此卡,并且又使用权窗,则相应于要运行问题的能量范围产生一个能量区间。如果有此卡,但无数据项,那么将产生10个能量箱,每个Ei=10i-8Mev,J=10。
格式: Pdn P1 P2 … Pi … PI
n:计数号;
Pi:栅元i对探测器贡献的概率,其缺省值为1。
I:栅元总数量。
PDn卡用于记数类型号为n的探测器记数,即它与Fnα卡配合使用。PDn卡可减少对指定探测器相对不太重要的栅元对探测器计数的贡献,节约了计算时间。此卡上对每一栅元给出一个概率值Pi,它表示对于栅元i中的任一碰撞点或源点,将以概率Pi判定它实行对探测器的记数。为保证无偏性,当判定实行记数时,原有的记录将乘上修正因子1/Pi。用户通过对离测点较远(以平均自由程计)的栅元置小于1的Pi值,便可提高运行速度。也可以将某些栅元的Pi值置为0,有选择地遏制了这些栅元对探测器的贡献。
此卡有一个很好的用处,即对几个完全一样的探测点(Fnα的n值给的不同),给出不同的PDn卡,便可将各个栅元对该探测器的贡献区分出来。
使用这个输入卡时,它上面的数据项可以都不是0或1,但不推荐这种用法。因为Pi值与粒子能量无关,若 Pi<l,对大栅元,将增大实现可靠抽样的负担。如果整个栅元离开探测器很远,倒是没有问题.
缺省此卡,所有栅元的Pi都取缺省值1。但对所有的探测器记数,使用PD0卡建立一组缺省的概率值以改变原来的缺省值。但这些缺省值又可被在DDn卡上填写的值所代替。
格式: dxc:n P1 P2 … Pi … PI
n: n为中子, p为光子;
Pi:栅元 i对DXTRAN球贡献的概率,其缺省值为 1;
I:栅元总数量。
DXC卡是选择卡。该卡的作用类似于PDn卡,但它是用于对DXTRAN球的贡献。也可以考虑使用DDn卡。