栅元描述的一般卡片格式
格式: j m d geom params
j:栅元的题目序号,写在第1-5列,1≦j≦99999。“题目序号”是指由用户按任意所希望的次序选定的序号,在一个文件中,可以不连续。然而MCNP将按照读入的顺序对栅元进行编号,从1开始按增量为1单调上升。我们将此称之为栅元的“程序编号”。必须明白题目编号与程序编号的区别,因为后面叙述的一些卡片将会引用到这种或那种编号。不过多数情况下,用户给出的序号往往与程序编号一致,这也是避免混乱,减少错误的最可靠方法。
m:栅元的材料号,它是材料卡(Mm)中相应材料的序号。真空栅元,m=0。
d:栅元材料的密度。填入正值时,表示是原子密度(单位为10 24次方个原子/cm3);填入负值时,则是质量密度(单位为gm/cm3)。对于真空栅元,该项不填,直接列出下一项内容。
geom:栅元的几何描述。它列出界定该栅元的所有曲面号(带有数将,表示坐向),及
描述这些曲面所定义的曲面之间关系的布尔算符,布尔算将包括交(Intersection),用一个空格表示;联(Union),用“:”号表示;余(Complement),用“#”表示。这里所说的曲面号及栅元号都指的是题目编号。如果定义该栅元要用到一个虚设曲面,则它必须作为界定该概元的一个曲面列出(同样带有坐向数符)。
param:任选的栅元参数说明,其形式为:关键词=数值。
在几何描述中,一个曲面的坐向是非常重要的概念。假如有f(x,y,z)=0这样一个S曲面,对几何空间中任一个点(x0,y0,z0),若有f(x0,y0,z0)>0,则称该点对于曲面S有正的坐向,用+S表示,“+”号可不写。反之,若有 f(x0,y0,z0)<0,则称该点对于曲面S有负的坐向,用-S表示。在栅元描述卡上列出的每一曲面,都将整个几何空间分为两个区域;分别对这一曲面具有正的和负的坐向。而所列带有符号的曲面,恰好规定了该栅元选用的是正的或负的坐向区域。
在运用布尔算符描述曲面之间的关系时,与“:”号相邻的空格将无任何意义;但没有“:”号相邻的空格作交处理。栅元描述方式中包括使用括号,每一个括号(‘(’或‘)’)等价于一个空格,所以括号两边的空格可有可无。布尔运算的次序为先#,其次是交,最后是联。如果有括号,先括号内,后括号外。#(- - -)形式,括号内是对某一个栅元的曲面描述—一栅元关系组。这一形式定义的几何区域是由不属于括号内描述区域的点组成的空间。对一个简单栅元(没有联或余运算符)的描述,只用一个空格以分隔界定描述该栅元的曲面即可。
假设块c1中含第16号材料(相应于ml6卡描述的材料),材料密度为 4.2gm/cm3。该栅元由三个曲面s1、s2及 S3界定,那么它的描述卡为:
c1 16 -4.2 s1 -s2 s3
这理,c1及si都是实际数值项的符号表示。如果c1是真空块,则该卡应写成:
c1 0 s1 -s2 s3
栅元参数可以通过在栅元卡上的指定而代替在数据卡上的指定。指定形式为:关键词=数值。可定义的关键词有: imp,vol,pwt,ext,fcl,wwn,dxc,nonu,tmp及与重复结构描述功能有关的 4个栅元参数卡:U卡,TRCL卡,LAT卡和 FILL卡。在栅元参数卡上的等号是可选择的,可用一个空格代替。如果某栅元的参数已在该栅元卡上指定,则不能再在数据卡上指定。允许部分栅元的参数在栅元卡上指定,而另一部分栅元的参数可在数据卡上指定。但与重复结构描述有关的栅元参数最好在该栅元卡上指定。这种方式可减少混乱。
例如:
10 16 -14.2 1 -2 3 imp:n=1 imp:p=8 u=1