MCNP输入的描述资料下载.pdf
《MCNP输入的描述资料下载.pdf》由会员分享,可在线阅读,更多相关《MCNP输入的描述资料下载.pdf(73页珍藏版)》请在冰点文库上搜索。
有些用户喜欢保留附加行,如改变这个问题的不同版本或者与这个输入文件本身有关的文本信息,这个终止符能防止读入这些附加的信息。
C.接续运行接续运行接续运行文件是用于继续计算一个早先被终止计算的问题(例如,对于一个工作先运行2小时,接着再运行一段时间),也可以用于对早先终止计算的问题进行重新编辑输出文件。
一个接续运行文件必须在信息块或者在MCNP的执行行中含有C或者CN以指明是一个接续运行方式,除非用Cm转储文件M处开始接序运行,否则将从最后一个转储文件处开始。
除了在MCNP的执行行中用C或CN选项,还有两个文件对接续运行非常重要:
(1)省缺名为RUNTPE的转储文件;
(2)一个可选的接续运行文件(省缺名为INP)。
转储文件,它是在初始运行时生成的,这个文件内容包括几何参数、截面参数、问题参数及重新开始运行所必需的全部其它信息,此外该运行在各个阶段的结果也记录在一系列的转储文件中。
转储文件存储的时机将在PRDMP卡中详细讨论,用户可从这些转储中的任何一个启动运行。
CN执行信息选项与C选项只有一点不同,即在接续计算中的转储数据直接存放在转储文件RUNTPE文件的固定数据部分之后,而不是放在原接续位置之后,新的转储内容覆盖旧的转储内容,这为用户防止RUNTPE文件不受控制地增长提供了方法。
RUNTPE的增长可在PRDMP卡中NDMP输入上控制。
在可选的接续运行输入文件中,必须把CONTINUE这个词作为第一行(标题卡)的第一个输入,或者也可放在信息块和空行分隔符之后,字符可以是大写、小写或大小写混合。
接续运行的输入文件格式如下:
信息块可选空行分隔符可选CONTINUE数据卡空行分隔符推荐使用其它可选接续运行的输入文件中所允许的数据卡是初始运行输入文件中数据卡的子集,允许使用的数据卡是FQ、DD、NPS、CTME、IDUM、RDUM、PRDMP、LOST、DBCN、PRINT、KCODE、MPLOT、ZA、ZB以及ZC。
如果上述项都没有改变(并且计算环境允许执行行信息),则不需要接续输入文件;
仅需要运行文件RUNTPE及执行行上的C选项。
例如,如果一个作业运行一分钟左右,因时间限制中断,但又想多运行一些粒子,这时可简单地执行命令行上的C或者CN信息,作业将从中断的地方重新开始运行,直到另一个新给的时间限制或者粒子截断或者人为停止。
这个例子假定来自初始运行的一个名为RUNTPE覆盖文件在当前目录下。
完整的接续运行执行行选项为Cm或者CNm,m特指从转储文件RUNTPE接续运行第页2006年6月MCNP4B教程教程的覆盖号,如果m没有特指,最后覆盖号被默认使用。
如果前次运行是因处理完指定的总粒子数(NPS)而终止的,那么在接续运行时,必须给出NPS卡,增大需处理的总粒子数NPS,NPS包括接续运行和初始运行的所有的初始粒子历史,接续运行中的CTME是指这个运行所花的时间,而不是总时间。
为了运行更多的KCODE循环,只有第四个输入KCT起作用,和NPS卡类似,KCT包括以前的循环在内的所有循环。
在接续运行文件中,NPS卡上给出是负值时,将只是对前次转储的中间结果进行编辑输出,并不再多计算历史。
当原有的输出文件已经丢失。
或者想通过PRINT卡或FQ卡改变输出内容时,使用这一方式是十分方便的。
但在初始运行使用FILES卡时要谨慎小心。
见3-119页。
D.卡片格式卡片格式INP输入文件的每一行(称之为一个卡)都限于使用1-80列。
英文字符可以用大写、小写或大小写混合。
大部分输入卡按行填写;
然而,对数据卡允许按列格式填写。
注释可以加入到任意一个输入卡中。
$符号(美元符号)为它所在那行数据的结束符,在$符号后面的内容作为注释。
空行可以被用作分界符和终止符,数据项用一个或多个空格分开。
输入文件中,在标题卡之后及最后的空行结束卡之前,任何地方都可插入注释卡。
注释卡的第15列必须有字母“C”,且后面至少要有1个空格。
注释卡仅在输入内容的原形输出部分打印出,不再在输出文件的其它任何部分出现。
FCn卡是作为用户注释用,它作为记数类型n的表头文字,作为记数的标题。
SCn卡也是作为注释用,它作为源概率分布n的表头文字。
标题卡被限定在180列的范围内使用,可以包括任何用户所期望的信息甚至也可以是空白。
它通常描述的是特别问题的信息,但要注意其它地方使用空格或空行表示分界符和终止符。
1.行输入格式行输入格式栅元卡,曲面卡,及数据卡都必须从第15列开始,数据输入紧跟在卡名或卡号和指定粒子之后,而这些数据项又被一个或多个空格分开。
如果15列为空白,则表示它是前一张卡的继续卡。
而至少一个空格以后的符号&
表示随后的卡是继续卡,继续卡数据限定在1-80列。
完全空白卡为在两个输入块之间的定界符。
一个数据项必须在一张卡上写完,不得跨到下一张卡上。
对任何给定的带有粒子标识符的类型卡只有一张。
需要整数输入的地方必须填写整数。
其它数据可填写为整数或浮点数,MCNP能正确地读入。
实际上,非整数数据能用任何FORTRAN语言接受的E格式描述。
为书写方便,可以使用四项书写功能:
AnR功能功能,表示将它前面的一个数据再重复n遍(R表示重复)。
例如:
24R和22222是一样的。
BnI功能功能,表示在与其前、后相邻的两个数之间插进n个线性插值点(I表示插值)。
1.52I3.0就相当于给出1.52.02.53。
在XnIY这种结构中,如果X和Y是整数,并且Y-X刚好是n+1的倍数,则产生标准的整数插值。
否则产生实数插值。
而y被保存到所有的事例中。
在上述例子中,2.0是实数,但不一定精确,而14I6=123456的全部插值点都是精确的。
CXM功能功能:
它表示的数值是等于前面的数据与X之积。
112M2M2M4M2M2M相当于112481664128256DnJ功能功能:
表示在它出现的卡片上,从它所在位置开始跳过n项数据,使得它们使第页2006年6月MCNP4B教程教程用缺省值。
例如下面两个卡是等同的:
DD.11000DDJ1000JJJ也等于3J。
这个功能使你能够简单地跳到卡片上的特定数据项,而不须给出前面的数据项数据,当想要使用缺省值,而又记不清它是什么值时,使用这一功能是非常方便的。
DBCN7J5082是另外一个例子。
上述四项功能对整数及浮点数的数据项都适用。
如果nR,nI及nJ项中的n缺省,则假设n为1。
如果XM项省略X,则是致命错误。
处理这几个特殊相邻项的规定如下:
(1)nR前面必须有一个数或有由R或M所产生的数据项。
(2)nI前面必须有一个数据或者有由R或M所产生的数据项,而且后面还必须有一个数。
(3)XM前面必须有一个数或者有由R或M所产生的数据项。
(4)nJ前面可以是除了I以外的任何内容,或作为输入行的开始。
例:
13M2R=133313MI4=133.5413M3M=13912R2I2.5=1111.52.02.51R2M=1121RR=11112I43M=12341212I42I10=123468103J4R是不合法的14I3M是不合法的14IJ是不合法的2.列输入格式列输入格式列输入对栅元参数及源描述是非常有用的,按行排列的栅元重要性及体积可读性不好,而且当增加或者删掉栅元时常常引起错误。
用列输入格式,一个栅元的所有栅元参数是放在标有这个栅元名字的那行上。
如果删掉一个栅元,用户只需删掉栅元参数这一行而不需在每一个栅元参数卡上寻找属于这个栅元的数据项。
对源描述,相应SI、SP及SB的数据逐个放在每一行上。
源选项和其他缺省值都在下一行上,显式输入,不用续行符号&
,如果有则被忽略。
用列格式,卡片名字逐个放在一个输入行上并且在这些卡片名字下面按列列出数据项。
#号是放在有卡片名字的这一行的15列的某列上。
卡片名字必须全部都是栅元参数、或全部都是曲面参数、或全部其它。
如果一个卡片名字在一个#卡上出现,则在同一个输入文件内这个卡片决不能用行格式。
如果在#卡上输入数据行的数据比卡名还要多,那第一个数据认为是一个栅元或曲面号。
如果填写了某一个栅元名字,则必须填写全部的栅元名字。
但填写栅元名字时,不必按栅元卡的同一个顺序。
如果没有指定栅元名字,则缺省为按栅元卡描述的顺序的卡片。
这些约定同样适合于曲面参数,但目前仅有一个曲面参数(AREA),因此曲面参数的列输入用途不大。
在一个输入文件中,可有多个列数据块,一般有一个栅元参数列数据块,及对每个源分布描述对应的一个个列数据块。
如果使用大量栅元参数选项,将需要增加列数据块。
在每一列上的数据项不需恰好就在卡片名字下,但是为了易读最好放整齐些。
输入数据的列格式当然适合等长度的列,但是不整齐的数据列不禁止如果一个较长的列是在较短列的第页2006年6月MCNP4B教程教程右边,则较短的列必须用足够的J填齐,以消除那一列数据栏中的不明确性。
在列格式输入中使用特殊功能项(R、M、I及J)不如在行格式输入中使用合适,但是也不禁止使用。
当然它们是被作为列功能而不是作为行功能。
如果有栅元或曲面号,则象9R这样的多重特殊功能都不能使用。
一个列输入块的格式是:
#S1S2-SmK1D11D12-D1mK2D21D22-D2mKnDn1Dn2-Dnm
(1)#占用15列。
(2)每一行仅有80列。
(3)Si的每一列,l不超过n时Dli为一个规则的输入卡。
(4)Si必须是MCNP卡片名字。
它们必须全部是栅元参数、全部是曲面参数或者全部其它。
(5)D1iDni必须是Si卡的有效项,但Dl+1,iDni当中可以是某些J或者J项后面跟有一些空格的除外。
(6)如果Dji不是空格,则Dj,i-1也必须不是空格。
如果有必要可以使用J功能使Dj,i-1不为空格。
(7)Si不能再在输入文件的其它地方出现。
(8)Kj是整数选择项,如果任何一个Kj不是空格,则全部Kj也都必须不是空格。
(9)如果Si是栅元参数卡片名字,那么,若要填写Kj,则Kj必须是栅元名字。
这个规定也适合曲面参数。
(10)如果Kj存在,Dji必不能是多重特殊功能项(如9R)。
E.粒子标识符粒子标识符有几个输入卡需要粒子标识符以区别中子和光子及电子的输入数据,这些卡是IMP、EXT、FCL、WWN、WWE、WWP、WWGE、DXT、DXC、F、F5X、F5Y、F5Z、PHYS、ELPT、ESPLT及CUT。
粒子标识符是由上述卡片名字后面的冒号、字母N及P或E直接组成的。
在粒子识别符后必须有至少一个空格。
例如,填写中子重要性的卡为IMPN;
填写光子重要性的卡为IMPP,对于一种以上的粒子指定相同的项可以用同一个单卡代替多个卡,例如IMPE,P,N110。
对于计数卡,是在含有计数号的卡片名字后面有粒子标识符,如*F5N表明一个中子对点探测器能量计数。
在热计数情况下,二个粒子标标识符可同时出现,F6N,P表明中子和光子的热计数混合。
F.缺省值缺省值MCNP的许多输入参数都有缺省值,因此当缺省值符合用户要求时,用户并不需要每次都必须给出各个输入参量的值。
当省略某张输入卡时,则该卡上的全部参数均使用缺省值(例如,带有5个入射光子的中断卡)。
然而,只想改变一张卡上的某一个特定缺省参量时,在它前面的参量仍需指明,或者用nJ方式跳过前面那些使用缺省值的参量也可。
例如,光子截断卡CUTP3J-10表示前3个参量使用缺省值,只改变第四项的值。
G.输入错误信息输入错误信息MCNP对输入文件出现的错误作广泛地检查(超过400项),如果用户违反了输入说明第页2006年6月MCNP4B教程教程的规定时,将在终端上及输出文件中打印致命错误信息,并MCNP将不再做粒子输运计算,作业中断。
第一个出现的致命错误是真的;
其后的一系列错误有的可能是真的,也可能是假的,这要根据前面出现的致命错误的情况而定。
若在MCNP执行行上指定了FATAL项,则MCNP忽略致命错误,运行照常,但是这时用户应该特别注意。
MCNP还给出很多警告性的而不是致命的错误,对这些警告用户也不应忽视,在进行重要计算之前,应当搞清它们的定义再作决断。
除了致命错和警告信息以外,MCNP在遇到任何严重问题之前(如零作除数)立即出现BADTROUBLE信息,并运行终止。
在INP文件中,用户输入错误是给出严重问题BADTROUBLE信息的最常见原因,这些错误信息提示表明检验校正的必需性。
H.检查几何错误检查几何错误MCNP在处理INP文件的数据时,不能检查出一种非常严重的输入错误,而在粒子丢失之前,MCNP不能查出栅元的重叠及栅元之间的空隙,并且即使出现了这种情况,准确的错误性质可能仍然不清楚,因此在大型的正式计算之前,我们必须检查几何错误。
利用MCNP的几何绘图的功能从几个不同的方向和尺寸检查输入问题几何,确保你看到的和真实情况一样,也许任何几何中的缺口和重叠都将象破折号那样出现,如果曲面的每一边的每个点没有对应的栅元,那曲面与平面的交线为破折线,如果绘图平面与问题某一平面重合或在源分布在许多栅元问题中有DXTRAN球,破折线都会出现。
另一种方法是建立和运行一个简短的问题,对这个问题,用外源的粒子轨迹注满这个系统。
对INP文件需要做如下改动:
(1)增加一个VOID卡,这个卡将废弃这个问题的某些其它说明使全部栅元为真空,把加热的计数转化为通量计数,并关掉所有FM卡。
(2)对这个问题要增加另一个栅元及一个大的球面,这个球面包围这个系统并且这个系统的外界栅元被这个新的曲面分成二个新的栅元:
一个是新的球面和原系统之间空间的一个栅元,另一个是这个新球面的外边界空间,这个空间是现在的外界栅元。
必须使新的栅元重要性不为零,实际上,最好使所有非零重要性相等。
如果这个系统是一维或二维无限,使用一个或一个以上平面代替一个球。
(3)用一个内向曲面源原来代替原来的源,用来注满几何体,源的说明为:
SDEFSUR=mNRM=-1其中m是在上面第二步中新增加的球面号。
如果新的曲面是一个平面,则必须用POS和RAD或者X、Y及Z的源描述方法指定源分布的部分。
由于没有碰撞,一个短时间运行将产生大量的轨迹穿过系统,如果有几何错误,将会引起某些粒子丢失。
当一个粒子首次丢失时,不管是在一个具有VOID卡的特殊运行或是一般运行,都要重新运行这个历史并在OUTP文件上产生某些特殊的输出,在重新运行期间打印出事件记录。
这个记录将显示跨越的所有曲面并将告诉你这个粒子在几何上走向坏点的径迹。
当这个粒子再次丢失时,打印出那个坏点上粒子的状态描述。
这样,你可从输出结果上推出引起粒子丢失的原因。
如果粒子丢失的原因仍然不清楚,一般常用的有效方法是画出几何图形,几何图的坐标原点放在丢失粒子的坐标点上,然后将绘图平面的水平轴选为丢失粒子运动的方向。
这样错误的原因很可能就在画图上出现破折号,或者画图和所想要的图之间有些不一致。
第页2006年6月MCNP4B教程教程.栅元卡栅元卡格式:
jmdgeomparams或jLIKEnBUTlistj:
栅元序号,1j99999。
m:
m=0如果是真空栅元m=材料号,如果栅元不是外真空。
这表明栅元包含Mm卡上指定的材料m。
d:
栅元材料的密度。
填入正值时,表示是原子密度(1024个原子/cm3),填入负值,则是质量密度(g/cm3)。
geom:
栅元的几何说明。
它列出界定该几何块的所有曲面号,以及描述这些曲面所定义的区域之间关系的布尔算符。
params:
任选的栅元参数说明,其形式为关键词=某个值。
n:
其他栅元名List:
栅元n与栅元j不同属性的定义,其形式为关键词=指定值。
在几何说明中,曲面号的正负号表示指定点在曲面旁边的那一个区域内。
对于确定点前所加的正负号是可选的。
布尔算符包括交(Intersection),联(Unin),及余(Complement)。
交运算符是隐式的,用一个空格表示;
联用号表示;
余用#号表示。
使用括号可控制布尔运算的次序,即先括号内的运算。
括号、算符都起到分界符的作用,有它们的地方,不用再加空格。
缺省的运算顺序是先#,其次是交,最后是联。
一个余运算操作符后的不带括号的数字被认为是栅元号,并且认为是栅元号指定几何的缩写。
例子:
30124$定义栅元3#3$等价于下一行#(-12-4)对一个简单栅元(没有联或余算符)的描述用空格来分隔列出的界定曲面和模糊曲面(用来使得几何定义更明确),这些模糊曲面和栅元曲面的正负号又由有栅元相对于这些曲面的坐向决定。
在栅元卡上可定义栅元参数以代替在输入文件中数据卡部分定义的栅元参数。
空格与等号是等价的。
如果一个栅元参数被输入栅元卡,那么该栅元参数卡不能再出现,甚至在在列形式的输入卡上也不能出现。
栅元卡上定义栅元参数的格式是:
关键词=某个值。
这儿所允许的关键词是:
带有粒子标识符的IMP、VOL、PWT、EXT、ECL、WWN、DXC、NONU、TMP,及与重复结构能力有关的4个栅元参数卡:
U卡,TRCL卡,LAT卡,FILL卡。
像其他任何栅元卡一样,这四个卡也能被放入INP文件的数据卡部分,我们推荐把每个栅元的助记符和输入都放在栅元描述之后的栅元卡行上。
特别当TRCL和FILL卡的输入放在栅元卡行上时,他们很长而且复杂但似乎又能被简单的概括起来。
LIKEnBUT用栅元材料号和密度等关键字,它们各自的助记符是MAT和RHO,这些关键字仅能被用在LIKEnBUT结构中,在通常的栅元描述中,材料号和密度数仍是栅元卡的第二和第三个输入。
在栅元卡上用两个方法填写TMP和WWN数据。
一种是使用关键词=一个值的形式,如TMP1:
n=一个值,TMP:
n=一个值等等,或者使用一个特殊形式,即在关键词TMP:
n后面按照THT卡上相应项顺序填写这个栅元的全部温度。
WWN卡的形式和THT卡类似,即WWN1:
n=一个值或者WWN:
n后面跟有这个栅元的能量区间的所有权重下限。
10164.2123IMP:
N=4IMP:
P=8EXT:
N=-.4X就是栅元10被密度为4.2g/cm3的材料来填充。
栅元是有曲面1的材料16填充,栅元是由曲面1的正向,曲面2的负向和曲面3的正向组成。
在栅元10中,中子重要性为4,第页2006年6月MCNP4B教程教程光子的重要性是8。
栅元10中的中子在负x方向受参数为0.4的指数控制。
当你准备栅元卡时,这里有一些预防措施可以借鉴:
(1)避开过于复杂的栅元。
当问题是由大量简单栅元而不是少量复杂栅元构成时,MCNP会处理的快一些,栅元描述中的最大值项(曲面号、算子、括号)包括所有求余运算符的展开是100,绝大部分栅元描述都应该比它短得多。
(2)避免将不必要的曲面通过勉强的求余运算符加到栅元的几何描述中,额外的曲面会使处理问题的速度变慢,而且会破坏体积和面积计算的必要条件。
(3)经常用MCNP的几何绘图能力来检查问题的几何是否正确。
(4)为了发现问题的几何错误,可以用外部源粒子涨满系统。
(5)如果你要增加或删除任一个栅元,就要相应地改变所有的栅元参量。
如果在栅元参量卡中使用列格式,就会使这样的处理容易。
可以在栅元卡上定义栅元参数值,便于完全地删除相应的栅元参数卡。
A.栅元描述的简写形式LIKEnBUT特征对具有大量重复结构的问题是非常有用的,栅元j继承了栅元n中没有被定义的所有属性的值,在输入文件INP中,栅元n的栅元卡必须在栅元j的栅元卡前面。
在BUT后面出现的任何卡的名字均为栅元卡参数,并且它只能在栅元卡上出现,不能是输入文件INP中数据块的任何卡。
23-3.7-1IMP:
N=2IMP:
P=43LIKE2BUTTRCL=1IMP:
N=10这就是说除了位置(TRCL=1)和中子重要性不同之外,栅元3和栅元2都一样,即栅元3的材料、密度、光子的重要性和2都是相同的。
.曲面卡曲面卡A.由方程定义曲面由方程定义曲面格式:
jnalistj:
j为曲面号,1j99999,如果这个曲面是反射面,在这个曲面号前面加一个*号或加一个+号的白边界。
n:
n为缺省或零,表示没有坐标变换。
n0特定TRn卡的卡号。