对每一个MCNP问题都有四个源中的一个。四个源包括:①通用源(SDEF卡);②曲面源(SSR卡);③临界源(KCODE卡);④用户提供的源。除了临界源之外,其它源都能使用源分布函数,分布函数分别由SIn、SPn、SBn和DSn卡指定。
下列卡片用于指定源的信息:
在一些情况下,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”的栅元,则是一个致命性错误。
l) SDEF通用源卡
格式: sdef 源变量=给定量
SDEF卡有三种定义形式:
①直接给值;②分布号前面冠以字母“d”;③源变量前面冠以字母“f”,分布号前面冠以字母“d”。 var=Dn的意思是从分布 n中抽样源变量数值; var fvar’Dn var从分布n中抽样,而n又取决于var’。只允许一级相关,每一个分布仅用于一个源变量。由MODE卡确定粒子的类型。
例如, sdef surf ext dl axs 0 1 0 erg d2
sdef erg d1 cel=d4 rad fcel d2 ext d3 axs=0 1 0 pos=35 80 105
把上述形式解释为三级源描述。当一个源变量有一个显值或者缺省(例如,单能)或者一个缺省分布(例如,各向同性角分布)时,第一级存在。当用一个概率分布给出一个源变量时,第二级出现。这个级要求SI,SP卡或者其一。当一个源变量与另一个源变量有关时,第三级出现,这个级要求DS卡。
MCNP按照特殊问题的需要建立的顺序抽取源变量。在抽取有关变量之后,必须抽取与它相关的每一个变量。如果一个变量的值影响另一个变量的缺省值,或者影晌抽取这个变量的范围;如SUR影响DIR变量,则必须按正确的顺序抽取这些变量。在MCNP中用于建立抽样顺序的方式是复杂的,并且这些方法也不总是有效的,如果失败了,将打印信息。在这种情况下,使用显值或者使用分布代替有关的缺省值。
WGT和EFF与其它源变量不同,它必须是一个数值,不允许是一个分布。大部分源变量是一个标量。VEC、POS及AXS是矢量。一个源变量要求一个值,如在SDEF、SI或者DS卡上那样,一般来讲单个数就够了。但是,VEC、POS及AXS必须是三个数为一组。即矢量的x,y,z分量。
用源变量 SUR、POS、RAD、EXT、AXS、X、Y、Z及CCC的各种组合以确定源粒子初始位置的坐标(x,y,z)。用这些变量可以指定三种不同类型的体分布及三种不同类型的面分布。这些分布的退化形式是线源和点源。使用SIn和DSn卡上的S选择,组合几个简单的分布能够近似比较复杂的一些分布。
三个体分布是长方体、球及圆柱。对每一个体分布SUR的值为0(缺省值)。一个体分布能够被用来配合CEL变量进行整个栅元内部的均匀抽样。指定一个完整包围一个栅元的长方体、球或者圆柱,并在体内均匀抽样。如果抽取的点在这个栅元内,则这个点就被接受。否则舍弃这个点,再重新抽取另一个点。如果使用这个技巧(栅元舍弃),必须确保抽样区域包含这个栅元的每部分,因为 MCNP对这方面没有检查的方法,使用后面描述的Cookie-cutter舍弃代替栅元或者和栅元舍弃联合。
用变量X、Y和Z指定一个长方体的体分布,如果这三个变量都是常数的长方体分布,它描述了一个点源。因此给出SDEF卡上的X、Y和Z三个变量的数值就可以指定单个点源。如果一个问题有几个源点,对每个点使用退化的球分布是比较容易的。长方体分布的其它退化情况是一个线源及一个直角平面源。?当栅元近似于长方体时,长方体分布是栅元舍弃技巧的一个有效形状。当栅元是一个薄板时,圆柱体分布是比较好的。然而,有一个限制,即栅元的面要垂直于坐标轴。
用变量POS和RAD指定一个球的体分布。X、Y、Z及AXS不必指定或者把这个体分布认为是一个长方体分布或者圆柱体分布。矢量POS的抽样值定义为球的中心。RAD的抽样值定义为从该球的中心到粒子抽样位置的距离。然后在半径为RAD的球面上均匀抽取粒子的位置。如果RAD是一个a=2的幂函数分布,则得到在一个球体内的均匀抽样。在这种情况下,a=2是缺省值。一个球的体分布多用于在两个球面之间的体上的均匀抽样。在SIn卡上指定RAD的两个半径,并按缺省得到SPn - 21 2 这样一个卡片。如果没有指定RAD,其缺省值所0。这是很有用的一种方式,因为RAD=0则指定了在POS位置上的一个点源。这样对一个在SIn卡上具有L的POS分布是指定一个问题的一组点源的最容易方法。
用变量POS、AXS、RAD及EXT指定一个圆柱体分布。圆柱的轴沿着AXS方向通过点POS。在一个以RAD的抽样值为半径的圆上具有抽取粒子的位置。圆心在这个圆柱的轴上。这个圆位于在离开POS的一定距离上垂直于AXS的一个平面上,这个距离是EXT的抽样值。圆柱分布常常用于圆柱壳体内的均匀抽样。在SIn卡上填写圆柱两端离开POS的EXT距离,并在SIn上填写RAD的内半径及外半径。按缺省值提供两个EXT之间的均匀抽样及两个RAD之间的幂函数抽样,a=1的RAD给出体内均匀抽样。一个有用的退化情况是EXT=0,这个情况提供一个平面上一个圆对称源。
对一个曲面分布,变量SUR的值不为0。如果指定X、Y及Z,他们的抽样值确定 粒子的位置。但在这种情况下,必须确保这个点在这个曲面上,因为MCNP不作检查。如果没有指定X、Y及Z,则在SUR曲面上抽取这个位置。曲面的形状确定抽取位置的方法。曲面的形状可以是棉球、球或者平面。栅元舍弃抽样在这里还适用。但是使用Cookie-Cutter舍弃能够做栅元舍弃要做的任何事情。
如果SUR的值是椭球曲面的名字,则在这个曲面上按面积均匀抽取粒子的位置。这种意义的椭球轴必须平行于一个坐标轴。目前还没有提供有关椭球曲面的简单的不均匀抽样或者偏倚抽样。可以使用一些Cookie-Cutter栅元的一个分布以产生粗糙的不均匀的位置分布。
如果SUR的值是球的曲面名字,则在这个曲面上抽取粒子的位置。如果没有指定矢量AXS,则在这个曲面上按面积均匀抽取粒子位置。如果指定AXS,则EXT的抽样值用来作AXS方向和球的中心到粒子位置点的矢量之间的夹角余弦。方位角在0o~360o之间的均匀抽取。由于EXT的不均匀分布,那么仅按极角的位置不均匀分布是存在的。使用EXT的偏倚分布使比较多的粒子从球离计数区域的最近的一边出发,通常来讲,EXT用指数分布函数(-31)为好。
如果SUR的值是一个平面曲面的名字,则在这个平面上抽取粒子的位置。抽取POS的值必须是这个平面上的一个点。必须保证POS确实在这个平面上,这主要是为了提高速度,MCNP不做这方面的检查。粒子的抽样位置距POS的距离等于RAD的抽样值。这个位置是在以POS为中心RAD为半径的圆上均匀抽取。如果RAD的分布是一个a=1的幂函数。则按面积均匀抽样,在上述情况下a=1是缺省值。
Cookie-Cutter舍弃对栅元源和曲面源都适用。如果CCC出现,则按上述方法抽取的位置在CCC内则接受,不在CCC内舍弃,重新抽取。这一点很象栅元源的栅元舍弃技巧。在指定Cookie-Cutter栅元,并且把这些栅元放在栅元卡描述的末尾,则不会遇到麻烦,不要使得Cookie-Cutter栅元太复杂。对一个曲面源Cookie-Cutter栅元要与源的曲面相交。用一些平面及柱面界定的均匀截面的无限长栅元是合适作一个平面源的Cookie-Cutter栅元。
注意:栅元舍弃或者CCC舍弃(Cookie-Cutter)与位置偏倚抽样结合几乎是不正确的游戏。如果使用这种组合,必须确保这种组合是一个合理的游戏,因为MCNP没有能力检查这种错误。
源变量SUR、VEC、NRM及DIR是用来确定源粒子的初始飞行方向。根据参考方向VEC抽取飞行方向。当然,也可以从一个分布中抽取。极角是DIR的抽样值。方位角0o~360o 之间均匀抽取。如果对一个体分布(SUR=0)没有指定VEC和 DIR,则按缺省值产生各向同性分布。如果对一个曲面源的分布没有指VEC,按缺省使用垂直于该曲面的矢量作为飞行方向。该曲面具有由NRM所确定的符号。如果在一个曲面源的分布上没有指定DIR,则按缺省值使用余弦分布P(DIR)=2*DIR,0<DIR<l。使用DIR的偏倚分布可使较多的粒子沿着朝向几何计数区域的方向出发。一般来讲指数分布最合适。
允许DIR取离散值。DIR=l给出在VEC方向上的单向源。当实际源与问题几何相距很大一段距离时,可用这种方式。在多数情况下,DIR的离散值将阻止对点探测器的直接记录,只有源是在一个平面曲面上,并在一个圆内(使用由 sp -21 1 抽取的RAD)按面积均匀抽取的源才能记录对点探测器的直接贡献。VEC垂直与这个曲面(这时VEC是缺省的),并且 DIR = l,则允许一个Cookie-Cutter栅元,通常 ARA的值也是必要的。
有效指标EFF适合于CCC和CEL这二者的舍弃。如果在任河源栅元或Cookie-Cutter栅元中接受率太低,这个问题被终止。终止的指标是MAX(成功的数目,.10)<EFF*试验数目。对稍微低效率的问题是假设EFF的缺省值为 0.01,但对低的源效率的问题是不可避免的,因此需要指定较低的EFF值.
例子:sdef(没有项)
这个卡指定一个位置在(0,0,0),时间为0,权重为1的14Mev的各向同性源。
重复结构通用源:当源在几何的重复结构部分时,SDEF卡上的CEL参数这部分说明是不同的。这时的CEL参数可以从0级到n级取值,不一定到底,形式如下:
C0: C1: … Ci … : Cn
Ci是填充栅元Ci-1的Universe里的一个栅元,或者Ci是一些栅元的一个分布Dm,或者Ci是0。如果Ci不是分布Dm时,可以有负号。抽取的位置和方向适用于填充栅元Ck-1的Universe,这儿Ck是第一个小于0的Ci,如果不存在Ck,则K=n。必须全部地指定K级以上的栅元。若Ci=0,则搜索这一级的栅元,如在n级的任一级中的栅元。若Ci是一个栅格,必须显示地表示这个栅格,而且该概格中的特定元素表示为:
… Ci(J1,J2,J3): …
其中(J1,J2,J3)为栅格Ci的特定元素。
在Ci ≠0的所有级的棚元中作舍弃,但是如果Ci在它的栅元卡上有一个负的Universe号,且Ci是在K级或K级以上,则不检查较高级的栅元。
2)SIn 是源分布的变量信息卡
格式为:
SIn x I1 I2 … Ik
n是分布号(n=l,999);
x为操作数,表示后面的参数怎样解释。可取参数如下:
x=(空格),表示矩形分布的边界,仅为标量;
=H,表示矩形分布的边界,仅为标量。
=L,表示离散的源变量数值;
=A,表示定义概率密度分布的点;
=S,已表示分布号。
I1 … … Ik :一些源变量值或一些分布号。
3)SPn 是源分布的概率卡
格式为:
SPn y P1 P2 … Pk or SPn f a b
n 是分布号(n=1,999);
y为操作数,表示后面的参数怎样解释。可取参数如下:
y=(空格),对应SI卡上H或L分布(相同于D的功能),对应SI卡上A分布的概率密度;
=D,表示SI卡上H或L分布的概率;
=C,表示给的是别卡上H或L分布的累积概率;
=V,表示仅对栅元分布,概率与体积成比例。
f = 表示分布函数的标志符(负数);
a,b=分布函数的参数(看表3-2)。
P1,P2,… ,Pk:为源变量概率。
4)SBn 是源变量的偏倚分布卡
格式为:
SBn y B1 B2 … Vk or SBn f a b
n是分布号(n=l,999);
y为操作数,表示后面的参数怎样解释。可取参数如下:
y=(空格),对于H或L分布相同于D的功能,对应SI卡上A分布的概率密度;
=D,表示SI卡上H或L分布的概率;
=C,表示给的是SI卡上H或L分布的累积概率;
=V,表示仅对栅元分布,概率与体积成比例。
f = 表示分布函数的标志符(负数);SB卡的f只允许是-21和-31。
a,b=分布函数的参数(看表3-2),对于SB卡,当f为-21或-31时才允许给值。
B1,B2,… ,Bk:源变量偏倚的概率。
默认值为:
SIn H I1 I2 … Ik
SPn D P1 P2 … Pk
SBn D B1 B2 … Bk
例3-3 sdef erg=dl pos=x y z wgt=w
si1 h E E … E
sp1 d 0 P … P
sb1 d 0 B … B
SP卡的第一种格式的第一项是正的或者非数值的,都表明这个卡及它的SI卡定义一个概率分布函数。SI卡的数据项不是源变量值就是分布号。SP卡上的数据项是SI卡上相应项的概率。
当使用H选择时,SI卡上的数据项的值是箱边界并且这些箱边界是单调上升的。SP卡上的第一个数据项的值必须是0,并且后面的数据项是相应箱的概率或是累积箱概率是取决于使用的选择是D或者是C。这些概率不需归一。抽样的过程是这样的:首先根据箱的概率抽取一个箱,然后在这个箱内均匀抽样。
当使用A选择时,SI卡上的数据项是定义一个源变量的概率密度分布的一些点,这些点的值必须是单调上升的,并且最大值和最高值定义这个源变量的范围。SP卡上的数据项是相应于SI卡上数据的概率密度,这些概率密度分布不需归一。在抽样过程中,概率密度
是在指定值的范围内线性插值。SP卡上的第一项及最后一项为0,但是也允许非0的值。
当使用L选择时,SI卡上的数值项是源变量的离散值,如一些栅元号或者光子能谱的一些能量。SP卡上的数据项是那些离散值的概率或者是累积概率是取决于使用D或C的选择。SI卡上的项不需单调增加。
S选择允许在一些分布中为进一步抽样来抽取其中的一个分布。这个特点不需要把分布合并在一起,并且如果这些分布中某些是离散的,其它是线性插值,这个特点就很重要。在SI卡上带有S选择所列的分布本身也可有S选择。MCNP能够处理大约有20层次这种结构。这种结构对任何特殊问题都是很重要的。在SI卡上的每一个分布号前面加一个D或者可省略D。如果一个分布号为0,则使用这个源变量的缺省值。借助于S选择,一个分布可在多个地方出现,但是一个分布只能用于一个源变量。
仅当源变量是CEL时才使用SP卡上的V选择。这个选择当栅元体积是粒子发射概率的一个因子时是很有用的。如果MCNP不能计算这样一个栅元体积,并且在VOL卡也没有给出这个体积,将出现致命错误。
使用SB卡是为抽样提供一个概率分布,这个概率分布与SP卡上的真的概率分布不同。它的用途是偏倚源变量的抽样以提高统计精度。为了校正这个偏倚调整每个源粒子的权重。适合于SP卡的第一种格式的全部规定也都适合于SB卡。
SP卡的第二种格式的第一项是负数表示使用一个内部分析函数以产生该源变量的一个连续密度函数,仅对标量变量可使用内部函数,表3-2给出了这些函数的描述。
1) SSW写曲面源卡
格式: ssw S1 S2(C1 … Ck)… Sn sym=m pty=n
Si = 带有粒子方向向里或向外的适当坐向的问题曲面号,对其粒子穿过信息写到曲面源文件WSSA;
Ci = 问题的栅元号。正数据项表示进入该栅元,负数据项表示离开该栅元。
m=0 假设不对称;
=l 假设球对称,问题曲面号表必须仅包括一个曲面而且一定是球面;
n=0记录所有轨迹;
=1记录中子轨迹;
=2记录光子轨迹。
利用这个卡写一个曲面源文件,供MCNP后续计算使用。必须注意,为了考虑反射效应,必须包含指定曲面以外的足够大的几何。这个卡允许在任河曲面名后面给出一个括入括号内的一个以上正的或负的栅元名列表。如果没有栅元名列表,将记录按适当方向穿过曲面的任何轨迹。如果给出列表,若轨迹按适当方向穿过曲面,并进入在列表中是正数据项的一个栅元或离开在列表中是负数据项的一个栅元,这个轨迹将被记录。如果使用SYM=l选择,每个粒子只有较少几个信息需要写到曲面源文件上,并且当读曲面源文件的时候,某些偏倚选择成为可利用的。缺省值为sym=0,pty=0,面源文件名为WSSA。
例如: ssw 4 -7 19(45 -46) 16 -83(49)
仅当按适当方向穿过曲面 19的轨迹进入栅元45或离开栅元46时才被记录。仅当按适当方向穿过曲面83的轨迹进人栅元49时才被记录。按适当方向穿过曲面4或7或16的轨迹,不管它进入或离开什么栅元都被记录。
2) SSR面源读卡
格式: ssr keyword=values keyword=values …
式中“=”为任选,下面描述各关键词:
old S1 S2 … Sn
问题曲面号表,它是在文件SSWA(现称为RSSA)中SSW卡上曲面的一个子集,
new Sa1 Sa2 … San Sb1 Sb2 … Sbn
本次运行中粒子开始的面源的曲面号。在a,b,…的变换位置中,对初始曲面源可重复n项。
col m 碰撞选项标志
m=-l从曲面源文件开始,仅那些粒子直接来此于非碰撞源;
m=l从曲面源文件开始,仅那些粒子在穿过记录曲面前已经碰撞;
m=0不考虑碰撞,所有粒子从曲面源文件开始。
wgt x 在输运中被接收的每个粒子权重乘以常数x。
tr n 坐标变换号。
使用当前问题的INP文件中必须存在的TRn卡上的变换,轨迹的位置和速度从辅助坐标系(写曲面源文件的坐标系)变换到当前问题的坐标系:
3.4.3 临界源描述
对于临界计算,只需两张附加卡片KCODE卡与KSRC卡。只要指定了KCODE卡,就不再使用通常的外源,如 MCNP的标准源及用户以 SOURCE子程序描述的源,此时必须用KSRC卡或现成的SRCTP文件给出用于选代的初始猜测源,猜测源只需给出一组源点的空间位置坐标,但任何源点都不可置于几何界面上。
临界计算不能用于MODE P的计算,但可用于MODE N P的情况,最主要的是在MODE N的情况运行。
1)KCODE 临界源卡
凡是计算Keff的问题都需给出此卡,其格式为:
kcode NSRCK RKK IKZ KCT MSRK KNRM
i)NSRCK是对每次迭代标定的源大小,而且常取在300<NSRCK≤3000的范围内,它没有缺省值。”
ii)RKK是对Keff的初始尝试值,它的缺省值为1.0。
iii)IKZ是开始累计记数之前应跳过的迭代次数,其缺省值是5。对给得不好的猜测源,甩掉起初的若干代是很重要的。
iv)KCT用以指定选代总次数,其缺省值为0,这就意味着不限定迭代次数,仅由给定的机时来终止计算。
v)MSRK是要提供存储的最大源点数,它的缺省值是4500或1.5*NSRCK,取较大的一个。若从SRCTP文件中读取初始源,它将被修改。
vi)KNRM为0时,将以权记数;若不为 0。则按粒子数记数。其缺省值为0。这仅适用于以Fn卡指定的记数,不适用于Keff本身的记数。
如果第一次迭代初始的源粒子数超过了本代的源,Keff的尝试值( RKK )可能给的太低,此时程序将打印一个注释运行,但不写新的源粒子,返回去读初始源计算Keff,并且用Keff代替RKK开始运行问题。如果第一次迭代之后,再次出现新产生的源粒子数超过了本代源粒子数,计算失败。
2)KSRC 临界计算的源点卡
格式为:
ksrc x1 y1 z1 x2 y2 z2 …
此卡可用于给出多达NSRCK组(x,y,z),它们是 KCODE临界计算的初始源点的位置。至少有一个点必须在栅元的裂变材料中。这些点不能在栅元边界上。
并非必须输入全部NSRCK个坐标点。例如可在每个裂变区给一个点,MCNP便将在给出的每个点上发出(NSRCK/裂变区数)个粒子。在很多情况下,这样做都是可以满意的,因为MCNP会很快地算出新的裂变源分布供迭代使用。类似地,也可使用以前计算中得到的SRCTP文件。如果该文件中不是NSRCK个粒子,则也将按照需要进行放大或缩小。
KSRC是一个选择卡,若不用此卡,则必须提供一个相应的SRCTP源文件。此卡或SRCTP文件所提供的源点,经过MCNP的处理便成为临界计算初始源的空间分布。
在实用上,当没有SRCTP文件可用时,最容易的办法是用KSRC卡输入一个点,只选代一、二次,就取其产生的SRCTP文件作为源文件,去掉KSRC卡,重新运行MCNP。SRCTP文件是在一个作业完成时或完成每次迭代之后写的;如果至少算完了IKZ次迭代,而且在机器中又至少驻留了一小段时间的话,SRCTP文件必定会写成功的。
在真空或0重要性栅元的点将被去掉。
如果在INP文件中没有SDEF、SSR和KCODE卡作为源的类型卡,就表示要由用户提供一个FORTRAN子程序SOURCE来描述粒子源。
在源子程序SOURCE中必须指定如下参量:
xxx:源粒子位置的x坐标(按cm);
yyy:源粒子位置的y坐标(按cm);
zzz:源粒子位置的z坐标(按cm);
uuu:源粒子方向的x轴方向余弦;
vvv:源粒子方向的y轴方向余弦;
www:源粒子方向的z轴方向余弦;
TME:源粒子的出生时刻(按10-8秒);
ERG:源粒子的能量(按MeV);
IPT:源粒子的类型;
WGT:源粒子的初始权(通常为1.0);
ICL:源粒子所位于的栅元i程序编号(对于面源,则是指将进入的栅元);
JSU:对于面源,它是源粒子所在的球面i程序编号;若不是面源,它必须登“0”。
在调用源子程序之前,MCNP已给出了一个各向同性的方向矢量(uuu,vvv,www),因此源如果是各向同性的,SOURCE中可以不再指定(uuu,vvv,www)。
在SOURCE中可以使用SIn、SPn及SBn卡来定义源粒子的能量、时间、方向及权重。如果需要,SOURCE 子程序中随处都可调用初始以[0,1]区间伪随机数的函数子例程RANG(0)。每一个 IDUM和 RDUN卡最多可输入 50个数据供 SOURCE子程序使用。INUM是整数,RDUM是浮点数。
如果在题目中使用了探测器记数或DXTRAN球,而且粒子源具有各向异性角分布,则用户还需提供一个子程序SRCDX,用来对每个探测器或者DXTRAN球指定PSC值。
在粒子库中有四个备用的变量SPARE1、SPARE2、SPARE3及ISPACE供用户使用。因此,用户可根据需要在SOURCE中对每一新的历史重置它们为“0”; MCNP是不做这种处理的。
格式:mgopt MCAL IGM IPLT ISB ICW FNW RLM
MCAL=F 向前问题
=A 共轭问题
IGM:问题中所有粒子类型的能群总数。负总数表示特殊的电子-光子问题。
IPLT:怎样使用权窗标志。
IPLT=0 表示IMP值是栅元重要性。对栅元重要性分裂和轮盘赌将忽略权窗(即使有)。
IPLT=l 表示必须提供权窗。并转换成与能量有关的栅元重要性。对能群零权窗下界产生重要性,它等于那个能群的最小的非零重要性。
IPLT=2 表示权窗正常执行。
ISB: 控制共轭偏倚,仅对共轭问题(MCAL=A)。
ISB=0 表示由无限平均通量偏倚碰撞。
ISB=l 表示由必须提供的权窗得到的函数偏倚碰撞。
ISB=2 表示没有偏倚碰撞。
ICV:产生权窗的参照栅元的名字。
ICW=0意味着不产生权窗。
ICW=0要求对非零重要性的所有栅元提供或计算体积。
FNW:产生权窗的非归一化值。在栅元ICW的最重要的能群中权窗的下边界值放置到FNW。
RLM:产生权窗的压缩范围。在打印产生权窗之前分别检查每个能群的权窗,看看最高和最低的比值是否小于RLM,如果不小,就压缩它们。
仅在多群计算时要求此卡。MCAL和IGM要参数(没有缺省值),其它是选择的。其缺省值:IPLT = 0, ICW = 0, FNW = l, RLM = 1000。 对任何参数,‘J’是不允许的值。负的IGM值允许一个以上粒子类型的截面数据表。这个特征目前仅适用电子/光子多群计算。在一个表中,跟在30个光子能群后面的50个电子能群的问题,将有IGM =-80。负的IGM值也表示在源和记录中能量变量是群号而不是Mev。根据不同的目的可以按几种不同方法组合 IPLT和 ISB。也要注意,某些选择使组合没有意义。共轭问题的输入文件可以有IMP卡和权窗卡( IPLT = 0, ISB = l )。 在权窗卡上的那些值就成为与能量有关的栅元重要性。直到现在,重要性是与能量无关的。
对于耦合的中子-光子计算的几种类型,所希望的结果是整个光子能谱的一个小子集的强度。两个例子是离散能量(line)光子和连续谱的高能末尾部分。在这种情况下,可利用偏倚中子诱发光子谱,产生仅仅感兴趣所那些部分。这可以通过使用PIKMT卡来完成。
格式: pikmt Z1 IPIK1 MT1,1 PMT1,1 MT1,2 PMT1,2 …
MT1,IPIK1 PMT1,IPIK1 … Zn IPIKn MTn,1
PMTn,1 MTn,2 PMTn,2 … MTn,IPIKn PMn,IPIKn
Zi :第i个输人的ZAID。可以指定ZAID的全部或部分;也就是说29000和29000.50是等效的。
IPIKi:控制ZAIDi的偏倚参数。如果IPIKi=0,ZAIDi就没有偏倚;用常规的抽样技巧产生来自ZAID的光子。如果IPIKi=-1,ZAIDi不产生光子。如果在核ZAID的PIKMT卡上没有数据项,从核 ZAID将没有光子产生,这是缺省。如果 IPIKi> 0,对 ZAIDi有i偏倚。IPIKi的值是所抽取的部分光子产生反应的数目。
MTi,j和PMTi,j:仅仅是具有IPIKi>0的ZAIDs所要求的。在这种情况下,必须输入IPIKi对MTs和PMTs。MTs是所抽样的部分光子产生反应的标识符。PMTs在一定程度上控制着所指定的MTs的抽样次数,这些输入不归一,具有正IPIK值的ZAID,在PIKMT卡上没有用MTs标识的任何反应将不被抽样。
此卡是选择的。如果没有此卡,则没有中子产生光子偏倚。为防止误解,给出以下说明:
l)当使用PIKMT卡的时候,最重要的是要意识到可能完全拒绝接纳零概率事件。换句话说,偏倚游戏未必是相对好的。
2)关于产生离散能量光子反应的MT标识符资料,我们不能使用。
3)这个特征对偏倚中子诱发光子谱产生很高能量的光子(例如,Er≥10Mev)也是有用的。若没有偏倚,则很少产生这些高能光子;因此,推断关于它们的可靠统计信息是困难的。当能量降低到感兴趣能量范围以下时,可以利用能量截断来结束轨迹。
例3-12 pikmt 26000.55 1 1020001 7014 0 29000 2 3001 2
3002 1 8016 -1
这个例子对14N产生所有光子-产生反应的常规抽样。中子与Fe碰撞产生的所有光子是来自具有MT标识符102001的反应。允许与Cu两个光子产生反应。因为PMT参数使具有MT标识符3001的反应被抽样两次,具有MT标识符3002的反应被抽样一次。从16O或在PIKMT卡上没有列表的问题中的其它同位素不产生光子。