3. 几何

复杂几何描述是蒙卡程序相对于确定论程序的重要优势之一。与世界上其它绝大多数蒙卡程 序一样，RMC采用基于层级空间的几何描述系统（universe-based geometry system）。

RMC的几何描述系统包括三类基本的几何描述单元，即，曲面（surface），栅元（cell）和 空间（universe）。一般地，物理系统由多个或单个层级空间组成，每个空间由一定数量的 栅元构成，栅元通过曲面方向（sense）的交并运算来定义。RMC几何输入模块的一般具有以下形式：

3.1. 曲面（surface）

曲面是RMC几何描述的最基本单元。考虑到MCNP用户的广泛性，RMC参考了MCNP的曲面定义方 式。曲面输入卡的定义方式为：

Surf <id> <type> {params} [bc = <flag>] [pair = <pair surf>] [value=<v1 v2 ... vn>] [time=<t1 t2 ... tn>]

其中，

  • Surf 是曲面输入卡关键词。

  • id是曲面编号。编号为正整数，且不允许重复。

  • type是曲面类型所对应的关键词，params是曲面方程参数。 表3.1 给出了RMC支持的曲面类型及相应的曲面方程。

  • bc是曲面边界条件（boundary condition）。bc = 0（缺省值）表示真空边界条 件，bc = 1表示全反射边界条件，bc = 2表示白边界条件，bc = 3表示周期边界条件。

  • pair是在周期边界条件下两个周期面的另一个，pair surf 是与该 surf 配对的另 一个周期面。 另外，使用周期边界条件有如下限制：

    • 周期边界条件的面只能是平面，即 P/ PX/ PY/ PZ；

    • 周期边界条件的面只能是最外表面，即周期边界条件的平面的其中一边必须是重要性为零的 Cell

    • 配对的面必须平行，且不能有变形

  • value卡和time卡结合使用，描述曲面方程中最后一个参数随时间的变化规律， 两个卡中输入的值的数目相等，表示当时间超过ti时，对应参数取为vi。

表3.1 曲面类型

类型

关键词

说明

方程

曲面方程参数

平面

P

一般

\(Ax+By+Cz-D=0\)

\(A,B,C,D\)

PX

垂直X轴

\(x-D=0\)

\(D\)

PY

垂直Y轴

\(y-D=0\)

\(D\)

PZ

垂直Z轴

\(z-D=0\)

\(D\)

球面

SO

球心在原点

\(x^2+y^2+z^2-R^2=0\)

\(R\)

S

一般

\((x-x_0)^2+(y-y_0)^2+(z-z_0)^2-R^2=0\)

\(x_0,y_0,z_0,R\)

SX

球心在X轴

\((x-x_0)^2+y^2+z^2-R^2=0\)

\(x_0,R\)

SY

球心在Y轴

\(x^2+(y-y_0)^2+z^2-R^2=0\)

\(y_0,R\)

SZ

球心在Z轴

\(x^2+y^2+(z-z_0)^2-R^2=0\)

\(z_0,R\)

圆柱面

C/X

平行于X轴

\((y-y_0)^2+(z-z_0)^2-R^2=0\)

\(y_0,z_0,R\)

C/Y

平行于Y轴

\((x-x_0)^2+(z-z_0)^2-R^2=0\)

\(x_0,z_0,R\)

C/Z

平行于Z轴

\((x-x_0)^2+(y-y_0)^2-R^2=0\)

\(x_0,y_0,R\)

CX

轴心在X轴

\(y^2+z^2-R^2=0\)

\(R\)

CY

轴心在Y轴

\(x^2+z^2-R^2=0\)

\(R\)

CZ

轴心在Z轴

\(x^2+y^2-R^2=0\)

\(R\)

圆锥面

K/X

平行于X轴

\(\sqrt{(y-y_0)^2+(z-z_0)^2}=\pm t(x-x_0)\)

\(x_0,y_0,z_0,t^2,\pm1\)

K/Y

平行于Y轴

\(\sqrt{(x-x_0)^2+(z-z_0)^2}=\pm t(y-y_0)\)

\(x_0,y_0,z_0,t^2,\pm1\)

K/Z

平行于Z轴

\(\sqrt{(x-x_0)^2+(y-y_0)^2}=\pm t(z-z_0)\)

\(x_0,y_0,z_0,t^2,\pm1\)

KX

轴心在X轴

\(\sqrt{y^2+z^2}=\pm t(x-x_0)\)

\(x_0,t^2,\pm1\)

KY

轴心在Y轴

\(\sqrt{x^2+z^2}=\pm t(y-y_0)\)

\(y_0,t^2,\pm1\)

KZ

轴心在Z轴

\(\sqrt{x^2+y^2}=\pm t(z-z_0)\)

\(z_0,t^2,\pm1\)

椭球面/双曲面/抛物面

SQ

轴平行于X、Y或Z轴

\(A(x-x_0)^2+B(y-y_0)^2+C(z-z_0)^2 +2D(x-x_0)+2E(y-y_0)+2F(z-z_0)+G=0\)

\(A,B,C,D,E,F,G,x_0,y_0 ,z_0\)

圆柱面/圆锥面/椭球面 /双曲面/抛物面

GQ

轴不平行于X、Y 或Z轴

\(Ax^2+By^2+Cz^2+Dxy+Eyz+Fzx+Gx+Hy+Jz+K=0\)

\(A,B,C,D,E,F,G,H,J,K\)

椭圆或圆形的圆环面

TX

平行于X轴

\((x-x_0)^2/B^2+ (\sqrt{(y-y_0)^2+(z-z_0)^2}-A)^2/C^2-1=0\)

\(x_0,y_0,z_0,A,B,C\)

TY

平行于Y轴

\((y-y_0)^2/B^2+ (\sqrt{(x-x_0)^2+(z-z_0)^2}-A)^2/C^2-1=0\)

\(x_0,y_0,z_0,A,B,C\)

TZ

平行于Z轴

\((z-z_0)^2/B^2+ (\sqrt{(x-x_0)^2+(y-y_0)^2}-A)^2/C^2-1=0\)

\(x_0,y_0,z_0,A,B,C\)

3.2. 栅元（cell）

栅元输入卡的定义方式为：

Cell <id> {surf_bool_definition} {cell_info}

其中，

  • Cell 是栅元输入卡关键词。

  • id是栅元编号。编号为正整数，且不允许重复。

  • surf_bool_definition指栅元的曲面布尔定义，由带方向的曲面和布尔运算符组 成，用来定义栅元区域。cell_info定义了该栅元的其它相关信息。下面将分别阐述。

3.2.1. 栅元的曲面布尔定义

栅元的曲面布尔定义由一系列曲面和布尔运算符组成，形如：

<±surf> <boolean> <±surf> <boolean> <±surf> …

曲面方向（sense）的定义为：如果点（x,y,z）在一个曲面方程f (x,y,z)的计算值 为f (x,y,z) > 0，则称该点对于这个曲面是正向的；若计算值为f (x,y,z) < 0，则为负向；若计算值为f (x,y,z) = 0，则表明该点在曲面上。 图3.1 给出了某二次曲面的方向所对应的区域：

../_images/surface_sense.png

图3.1 曲面方向示意图

RMC的布尔运算符包括交集（&）、并集（:）和补集（!）两种，并支持用圆括号调整运算优 先级。补集的优先级高于交集和并集；交集和并集的优先级相同，按照定义的 先后顺序进行逻辑运算；圆括号的优先级最高，并且可以使用多层圆括号嵌套，类似 于算术运算。假设栅元1和栅元2的几何描述分别为：

栅元1： (1 & -2) : 3

栅元2： 4 & -5 : !1

栅元1所表示的几何区域为：（曲面1的正向 ∩ 曲面2的负向）∪ 曲面3的正向

栅元2所表示的区域为：（曲面4的正向 ∩ 曲面5的负向）∪ 非栅元1。栅元2的另一种等价 描述方式为：4 & 5 : !( (1 & -2) : 3)。需要注意的是，若“!”之后紧跟数字，则表示非 栅元；若“!”之后为括号，则表示非曲面。

3.2.2. 栅元信息选项卡

栅元信息选项卡由一系列选项卡组成，主要用于描述栅元的物理和几何参数，包括材料、体 积、温度、层级填充信息、几何变换，等。

Cell … [Mat = <id>] [Vol = <vol>] [Tmp = <tmp>]
[Dens = <dens>] [Void = <flag>] [Fill = <id>] [Inner = <flag>]
[Move = <params>] [Rotate = <params>] [Noburn = <flag>]

其中，

  • Mat选项卡定义该栅元的填充材料，缺省值为Mat=0（真空）。

  • Vol选项卡定义该栅元的体积，单位为cm3，缺省值为Vol=1.0cm3。

  • Tmp选项卡定义该栅元的温度。用户可输入大于0的自然数，或小于0的整数。 当输入的值大于0时，表示输入了该栅元的温度，单位为K；当输入的值小于0时，表示 用户输入了该栅元的温度网格（详见网格章节），单位为K；当用户未输入该选项卡时， 默认用户输入了293.6（K）。 当栅元温度与栅元内填充的材料的温度不匹配时，如果用户未指定任何在线展宽选项， 程序将对核素截面热化区进行简单修正；若用户指定了在线展宽选项，则按照用户指定 的选项进行在线多普勒展宽或在线插值。

  • Dens选项卡定义该栅元的密度。用户可输入大于0的自然数，或小于0的整数。 当输入的值大于0时，表示输入了该栅元的原子密度， 单位为1024原子/cm3；当输入的值小于0时，表示 用户输入了该栅元的温度网格（详见网格章节），单位为g/cm3； 当用户未输入该选项卡时，默认用户输入的密度与所填充材料的密度保持一致。 需要注意的是，当前版本的程序不支持对燃耗栅元使用该选项。

  • Void选项卡指定中子进入该栅元后是否停止跟踪，主要用于描述真空边界以外的 区域。Void=0（缺省值），中子进入该区域继续跟踪；Void=1，中子 进入该区域停止跟踪；

  • Fill 选项卡定义该栅元内部填充的空间，详见后续章节。

  • Inner选项卡指定该栅元是否为内部栅元，即，在填充过程中未被外层边界分割。 Inner= 0（缺省值）表示非内部栅元，Inner =1表示内部栅元。指定 内部 栅元可以加速几何处理，但错误地指定内部栅元会导致几何跟踪出错，因此只建议 高级用户使用。

  • Move选项卡定义该栅元的平移向量，Rotate选项卡定义该栅元几何变换的旋转矩阵的转置， 分别类似于3.3.3 空间几何变换的move和rotate卡。对栅元进行几何变换时，是先旋转，再平移。 对于被填充的栅元，是连同其内部填充结构做整体变换。这里不支持对包含圆环面(TX/TY/TZ)的旋转。

  • Noburn选项卡定义该栅元是否参与燃耗计算。Noburn = 0（缺省值）表示不参与燃耗计算，Noburn = 0表示参与燃耗计算。

3.3. 空间（universe）

3.3.1. 单层空间

空间由一系列栅元组合而成，且这些栅元之间不能存在重叠或未定义区域。单层空间 输入卡的形式为：

UNIVERSE <id> [options]

其中，id是空间编号。options是与空间几何变换、重复结构相关的选项 ，形式如下，后面将具体述及。

[Move = <params>] [Rotate = <params>] [Lat = <params>] [DISP = <params>]
[Pitch = <params>] [Scope = <params>] [Sita = <param>] [Fill = <params>]

对任意的物理系统，至少需要一个空间来描述，这在输入文件中定义为Universe 0。例如， 以下输入文件是一个普通压水堆栅元的几何部分。这里将栅元1顺时针旋转90°。

/////// PWR pin: defined in single universe /////////////
Universe 0
cell 1 -10      mat = 1   move=0 0 0 rotate=0 -1 0 1 0 0 0 0 1// Fuel Pin
cell 2 !1 & -11 mat = 2                // Air
cell 3 11 & -12 mat = 3                // cladding
cell 4 12 & 13 & -14 & 15 & -16 mat= 4 // water
cell 5 -13 : 14 : -15 : 16 void = 1    // outside

Surface
surf 10 cz 0.4096
surf 11 cz 0.4178
surf 12 cz 0.4750
surf 13 px -0.63 bc = 1
surf 14 px 0.63 bc = 1
surf 15 py -0.63 bc = 1
surf 16 py 0.63 bc = 1

若该单层空间用于随机栅格扰动法填充上层重复结合结构，则需要用到DISP关键词对随机颗 粒球心位置进行扰动，具体方法及输入在3.3.3节进行说明。

3.3.2. 多层空间

对于复杂的物理系统，可能需要用到空间填充的描述方式，即，将某个空间填充到另一个空 间的某个栅元当中。注意，填充空间应涵盖被填充的栅元区域，否则该栅元区将存在 未定义的空白区域，造成粒子跟踪错误。

空间填充的选项卡内嵌在栅元输入卡中（参考3.2.2）：

Cell ... [Fill = <universe>]

对上例中的压水堆栅元，我们可以使用空间填充的方式来等价地进行描述，如下所示。首先 ，定义了燃料棒及慢化剂区域（Universe 1），然后将其填充至栅元格（cell 102）。

/////// PWR pin: defined in multilevel universe /////////////
Universe 0
cell 101 13 & -14 & 15 & -16 Fill = 1 // define a cell filled by a universe
cell 102 -13 : 14 : -15 : 16 void = 1 // outside the box

Universe 1
cell 1 -10      mat = 1 // Fuel Pin
cell 2 10 & -11 mat = 2 // Air
cell 3 11 & -12 mat = 3 // cladding
cell 4 12       mat = 4 // water

Surface
surf 10 cz 0.4096
surf 11 cz 0.4178
surf 12 cz 0.4750
surf 13 px -0.63 bc = 1
surf 14 px 0.63  bc = 1
surf 15 py -0.63 bc = 1
surf 16 py 0.63  bc = 1

3.3.3. 几何变换

RMC支持对空间（universe）的平移变换/旋转变换/随机扰动（随机颗粒的球心坐标）。几何变换选 项卡内嵌在空间输入卡当中：

Universe ... [Move = Mx My Mz]
[Rotate =Cx'x Cx'y Cx'z Cy'x Cy'y Cy'z Cz'x Cz'y Cz'z ]

平移变换的表达式为:

\[\mathbf{r'} = \mathbf{r} + \mathbf{m}\]

其中， \(\mathbf{r}=(r_x,r_y,r_z)\) 和 \(\mathbf{r}=(r_x',r_y',r_z')\) 分 别为变换前和变换后的空间任意一点的位置坐标， \(\mathbf{m}=(m_x,m_y,m_z)\) 为 平移变换向量。

旋转变换可以绕任意轴，其表达式为：

\[\mathbf{r'} = \mathbf{R} \cdot \mathbf{r}\]

其中， \(\mathbf{R}\) 为旋转变换矩阵。RMC实际要求用户输入的是旋转变换矩阵的 转置矩阵 \(\mathbf{R}^T\) ，其参数按照以下方式定义：给定某直角坐标 系 \((x,y,z)\) ，它经过该旋转变换后得到新坐标系 \((x',y',z')\) ， 则 \(\mathbf{R}^T\) 可表示为

\[\begin{split}\mathbf{R}^T = \begin{bmatrix} C_{x'x} & C_{x'y} &C_{x'z} \\ C_{y'x} & C_{y'y} &C_{y'z} \\ C_{z'x} & C_{z'y} &C_{z'z} \end{bmatrix}\end{split}\]

其中， \(C_{x'x}\) 表示 \(x\) 和 \(x'\) 两个坐标轴之间的夹角余弦，以 此类推。

注意，如果对某个空间同时做旋转变换和平移变换，应先旋转，再平移。对于多层空间， 对某个空间的几何变换总是连同其内部填充结构做整体变换。此外，Universe 0是基准空间，因此不允许对Universe 0做几何变换。

使用几何变换的方式，我们重新定义上面的压水堆栅元，如下所示。燃料棒和慢化剂区域 （Universe1）定义为与x轴平行，通过平移（move = 0.5 0.5 0）和旋转（rotate = 0 0 -1 0 1 0 1 0 0），填充到栅元格（Cell 102）中。

// PWR pin: defined in multilevel universe with coordinate transformation //
Universe 0
cell 101 13 & -14 & 15 & -16 Fill = 1 // define a cell filled by a universe
cell 102 -13 : 14 : -15 : 16 void = 1 // outside the box

Universe 1 move = 0.5 0.5 0 rotate = 0 0 -1 0 1 0 1 0 0
cell 1 -10      mat = 1 // Fuel Pin
cell 2 10 & -11 mat = 2 // Air
cell 3 11 & -12 mat = 3 // cladding
cell 4 12       mat = 4 // water

Surface
surf 10 c/x -0.5 -0.5 0.4096
surf 11 c/x -0.5 -0.5 0.4178
surf 12 c/x -0.5 -0.5 0.4750
surf 13 px -0.63 bc = 1
surf 14 px 0.63 bc = 1
surf 15 py -0.63 bc = 1
surf 16 py 0.63 bc = 1

对于随机栅格扰动法，DISP参数用于填充重复几何结构的单层空间，球心坐标扰动方程为：

\[\begin{split}x& =x'+(2\xi_{1}-1) \times \delta_{x} \\ y& =y'+(2\xi_{2}-1) \times \delta_{y} \\ z& =z'+(2\xi_{3}-1) \times \delta_{z}\end{split}\]

其中， \(\xi_i\) 为(0,1)之间的随机数， \(\delta_i\) 为对应坐标轴方向的扰动幅度。 注意：坐标转换的随机扰动幅度不能超过栅格的边界。

以下是随机栅格扰动法的一个输入示例：

////////  HTR 5*5*5 lattice, liu-sc, 2014-10-28 ////////
UNIVERSE 2 lat = 1  pitch = 0.1982 0.1982 0.1982    scope = 5  5  1  fill =
  1 1 1 1 1
  1 1 1 1 1
  1 1 1 1 1
  1 1 1 1 1
  1 1 1 1 1

UNIVERSE 1 move = 0.0991 0.0991 0.0991      Disp=0.0536 0.0536 0.0536
cell  3   -1       mat = 3   //fuel
cell  4   1 & -2   mat = 1     //1.1C
cell  5   2 & -3   mat = 4     //1.9C
cell  6   3 & -4   mat = 2     //SiC
cell  7   4 & -5   mat = 4    //1.1C
cell  8   5        mat = 4    //1.1C

3.4. 重复结构（lattice）

重复结构是一类特殊的空间，该空间由规则排列的网格组成。RMC支持常用的四边形重复结 构和六边形重复结构，它们在反应堆堆芯计算分析时最为常见。四边形重复结构可以建立在 1维、2维平面或3维空间，六边形重复结构建立在2维平面。

同时，为了描述弥散燃料介质随机模型，RMC在重复结构选项卡中增加了两个特殊选项，分别 为隐式模拟方法和显式模拟方法。 （注意：RMC的随机介质功能仅企业版本提供） RMC隐式方法采用弦长抽样法（CLS），通过概率分布函数 在线抽样确定随机介质燃料颗粒的位置，其特点是不需要显式构建众多燃料颗粒的位置，因此 在填充率上没有严格限制。其缺点是计算精确性比显式模拟差一些，并且填充率高时，实际 填充率与目标填充率存在差别，需要通过修正改善精度。RMC的显式方法采用随机介质序列添加法 （RSA），在输运计算前会显式确定每一个弥散介质燃料颗粒的空间位置，因此实际填充率即为 目标填充率。但由于方法的固有特性，填充率不宜过高，上限值为38%，越接近上限则产生颗粒的时间越长。

重复结构选项卡内嵌在空间 （universe）输入卡当中：

Universe … [Lat = <type>]

其中，Lat = 1表示四边形重复结构，Lat = 2表示六边形重复结构，Lat = 3表示 显式建模法随机几何，Lat = 4表示弦长抽样法随机几何。下面针对这两种重复结构类型分别阐述。

3.4.1. 四边形重复结构

图3.2 给出了四边形重复结构的示意图。四边形重复网格建立在xyz坐标系，坐标原点O建 立在第一个网格（编号为1）的左下角点。

../_images/lattice_mesh.png

图3.2 四边形重复结构示意图

四边形重复结构的选项卡为：

Universe … [Lat = 1] [Scope = <xNum yNum zNum>]
[Pitch = <xLen yLen zLen>] [Fill = <U1 U2 … UM>]

其中，

  • Lat = 1表示重复结构类型为四边形。

  • Scope选项卡定义重复网格在x，y，z方向的数量。特别地，参数为1表示该方向上 只有一层网格。例如，2维PWR组件的重复网格表示为Scope = 17 17 1。需要指 出的是，尽管程序支持直接定义3维四边形重复结构，但建议用户通过2维重复结构和1维 重复结构的填充方式来生成3维重复结构。

  • Pitch选项卡定义重复网格在x，y，z方向的宽度，参数必须为正。若某方向只有 一层网格，Pitch选项卡中对应的参数没有实际意义。

  • Fill选项卡依次定义网格内填充的空间（universe）的编号，一共 有 \(M = xNum \times yNum \times zNum\) 个编号。Fill选项卡的填充 次序为：先按x方向填充，再按y方向填充，最后按照z方向 填充。图3.2 给出了四边形重复结构的索引下标的编号方式 ，对应Fill选项卡的填充顺序，同时也对应重复结构计数器的编号。

在使用随机栅格扰动法时，需要对四边形重复结构填充的Universe添加disp选项，具体内容参见3.3.3后半部分。

3.4.2. 六边形重复结构

../_images/lattice_hex.png

图3.3 六边形重复结构示意图

六边形重复网格的排列方式见 图3.3 。不难发现，各个六边形的中心 按照平行四边形的方式 排列。平行四边形的两条边所对应的方向向量b1和b2位于xy平面内，b1与x方向重合。

六边形重复结构的选项卡为：

Universe … [Lat = 2] [Scope = <b1Num b2Num>] [Sita = <sita>]
[Pitch = <b1Len b2Len>] [Fill = <U1 U2 … Um]

其中，

  • Lat = 2表示重复结构类型为六边形。

  • Scope选项卡定义重复网格在b1和b2方向的数量。

  • Pitch选项卡定义重复网格在b1和b2方向的宽度。

  • Sita选项卡定义六边形网格其中一对邻边的夹角（如图所示），单位为度°。

  • Fill选项卡依次定义网格内填充的空间（universe）的编号，一共 有 \(M = b_1Num \times b_2Num\) 个编号。 Fill选项卡的填充次序为：先按 \(M = b_1\) 方向（即x方向）填充， 再按 \(M = b_2\) 方向填充，具体次序见 图3.3 给出的编号。

需要指出的是，与四边形重复结构不同的是，六边形重复结构总是建立在xy平面，坐标 原点O建立在第一个重复六边形（编号为1）的中心。通过平移和旋转变换，可以将其转换到 其它平面。

3.4.3. 弦长抽样随机几何（仅限企业版本）

弦长抽样随机几何的选项卡为：

Universe … [Lat = 3] [MATRIC = <UM>] [PARTICLE = <U1 U2 … UP>]
[PF = <pf1 pf2 … pfp>] [RAD = <rad1 rad2 … radp>]
[TYPE = <1/2/3>] [SIZE = <size>] [PFCORRECT]

其中，

  • Lat = 3表示弦长抽样法随机几何。

  • MATRIC选项卡定义基体所在空间（universe）的编号。

  • PARTICLE选项卡定义颗粒所在空间（universe）的编号。通过其后定义的Ui标识 分辨不同类型的燃料颗粒，一种燃料对应一个空间（universe）的编号。

  • PF选项卡定义颗粒占其填充的几何的体积份额。

  • RAD选项卡定义颗粒的半径。

  • PFCORRECT选项卡定义是否使用隐式方法自动修正，若选用该模式， 在临界计算前会对填充率进行自动校正。0表示不用校正方法，为默认选项， 1表示开启自动校正模式。(注： 隐式方法填充率自动修正模式不兼容固定源模式)

3.4.4. 显式建模法随机几何（仅限企业版本）

显式建模法随机几何的选项卡为：

Universe … [Lat = 4] [MATRIC = <UM>] [PARTICLE = <U1 U2 … UP>]
[PF = <pf1 pf2 … pfp>] [RAD = <rad1 rad2 … radp>]
[RSA = <0/1>] [TYPE = <1/2/3>] [SIZE = <size>]
[DEM = <0/1>] [TIME = <time>]

其中，

  • Lat = 4表示显式建模法随机几何。

  • MATRIC选项卡定义基体所在空间（universe）的编号。

  • PARTICLE选项卡定义颗粒所在空间（universe）的编号。通过其后定义的Ui标识 分辨不同类型的燃料颗粒，一种燃料对应一个空间（universe）的编号。

  • PF选项卡定义颗粒占其填充的几何的体积份额。

  • RAD选项卡定义颗粒的半径。

  • RSA选项卡定义随机颗粒的产生方式，RSA = 1 表示颗粒位置由程序内部使用RSA方法产生， 并产生存储颗粒位置的文本,名字为 “ random_geometry_[当前空间号(Lat=4)] ”；RSA = 0表示 球的位置从外部文件读取，外部文件的名字为 “ random_geometry_[当前空间号(Lat=4)] ”。

    注意 颗粒产生文件的格式为：

    x坐标 y坐标 z坐标 颗粒最外层半径 颗粒空间编号（用户编号）

    例如 :

    -1.67308E+00 2.92296E-01 1.07829E+00 4.55000E-02 1

    这里，坐标及半径的单位是cm，位置坐标是颗粒球心坐标，对应的坐标系 原点是随机介质区域中心

  • TYPE选项卡定义被颗粒填充几何的形状，TYPE = 1表示颗粒填充的形状为球体， TYPE = 2表示被颗粒填充的形状为圆柱体，TYPE = 3表示被颗粒填充的形状为长方体。

  • SIZE选项卡定义被颗粒填充几何的尺寸。TYPE = 1时，输入球的半径（一个参数）， TYPE = 2时依次输入圆柱体的半径和高（两个参数），TYPE = 3时，依次输入长方体的长宽高 （x,y,z三个参数）。

  • DEM选项卡选择是否使用DEM方法产生随机颗粒以获得更高的填充率，DEM = 1 表示颗粒位置使用DEM方法产生， 并产生存储颗粒位置的文本,名字同样为 “ random_geometry_[当前空间号(Lat=4)]” ,此时坐标系为随机介质区域中心坐标系， 这个文件可以被RMC直接读取进行计算；DEM= 0表示不使用DEM方法产生颗粒。需要注意的是：选取DEM的 方法获得更高的填充率时必须先进行RSA计算，RSA的输入卡片必须为1，且相关选项卡也必须填写完全。

  • TIME选项卡定义DEM方法产生小球的下落时间(自落体体时间),单位为s，一般建议时间设置为颗粒下落两倍的几何高度所需时间的两倍。 例如几何高度为1m时，建议时间设置为1s。

需要注意的是，显式建模法的Lattice所在的空间的坐标原点不是在球、圆柱、长方体的几何中心 而是在定位小球位置的网格中，坐标原点为网格的左下角角点。如下图所示：

../_images/explicit_coordinate.png

图3.4 显式建模法坐标原点

显式建模法的Lattice所在的空间去填充Cell的时候，要根据被填充Cell 的坐标原点进行必要的坐标平移。Particle空间定义建议球心定位于原点，即选用so类型面。

3.5. 几何模块输入示例

3.5.1. PWR组件

示例3.1 是一个PWR17×17组件（ 图3.5 ） 的输入示例。Universe 1和Universe 3分别为燃料栅 元和管道栅元，中心坐标在（0, 0, 0）。Universe 8为四边形重复结构。因为四边形重复 结构的下角点总是建立在（0, 0, 0），所以Universe 8的第一个网格的中心点位置为 （0.63, 0.63, 0）。将Universe 1和Universe 3按照向量（0.63, 0.63, 0）平移后，填充 到Universe 8的第一个网格中，然后按照四边形重复结构排列展开。

../_images/pwr_assembly_17by17.png

图3.5 PWR17×17组件


示例3.1 17x17组件几何输入示例
// STANDARD WESTINGHOUSE 17*17 ASSEMBLY MODEL. SHE DING : 2012-03-08 //
UNIVERSE 0
CELL 1 6 & -7 & 8 & -9 mat = 0 Fill = 8 // Assembly inside
CELL 2 -6 : 7 : -8 : 9 mat = 0 void = 1 // Assembly outside

UNIVERSE 8 lat = 1 pitch = 1.26 1.26 1 scope = 17 17 1 fill =
    1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
    1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
    1 1 1 1 1 3 1 1 3 1 1 3 1 1 1 1 1
    1 1 1 3 1 1 1 1 1 1 1 1 1 3 1 1 1
    1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
    1 1 3 1 1 3 1 1 3 1 1 3 1 1 3 1 1
    1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
    1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
    1 1 3 1 1 3 1 1 3 1 1 3 1 1 3 1 1
    1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
    1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
    1 1 3 1 1 3 1 1 3 1 1 3 1 1 3 1 1
    1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
    1 1 1 3 1 1 1 1 1 1 1 1 1 3 1 1 1
    1 1 1 1 1 3 1 1 3 1 1 3 1 1 1 1 1
    1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
    1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

UNIVERSE 1 move = 0.63 0.63 0   // Fuel rod
cell 3 -1     mat = 1 inner = 1 // Fuel
cell 4 1 & -2 mat = 3 inner = 1 // Air
cell 5 2 & -3 mat = 4 inner = 1 // Zr
cell 6 3      mat = 5           // water

UNIVERSE 3 move = 0.63 0.63 0 // Guide tube
cell 11 -4     mat = 5 inner = 1 // water
cell 12 4 & -5 mat = 4 inner = 1 // Air
cell 13 5      mat = 5           // water

SURFACE
surf 1 cz 0.4096
surf 2 cz 0.4178
surf 3 cz 0.4750
surf 4 cz 0.5690
surf 5 cz 0.6147
surf 6 px 0     bc = 1
surf 7 px 21.42 bc = 1
surf 8 py 0     bc = 1
surf 9 py 21.42 bc = 1

MATERIAL
mat 1 -10.196
    92235.30c 6.9100E-03
    92238.30c 2.2062E-01
    8016.30c 4.5510E-01
mat 3 -0.001
    8016.30c 3.76622E-5
mat 4 -6.550
    40000.60c -98.2
mat 5 9.9977E-02
    1001.30c 6.6643E-02
    8016.30c 3.3334E-02
sab 5 lwtr.60t

CRITICALITY
PowerIter population = 10000 50 300 // keff0 = 1.0
InitSrc point = 0.63 0.63 0

3.5.2. PWR堆芯

示例3.2 是一个PWR堆芯几何模块输入示例。简单起见，堆芯输入文件仅 包含一种类型的组件 。堆芯（Cell 1）被填充21×21的四边形重复结构（Universe 1），其中包括组件网格和反 射层网格。组件（Universe 3）是17×17的四边形重复结构，填充有燃料栅元（Universe 6） 和管道栅元（Universe 7）。

示例3.2 PWR堆芯几何输入示例
////////// PWR core. SHE Ding 2012-07-01 ////////////
UNIVERSE 0
CELL 1 -10 mat = 0 Fill = 1 // Core inside
CELL 2 10 mat = 0 void = 1 // Core outside

UNIVERSE 1 lat = 1 pitch = 21.42 21.42 1 scope = 21 21 1 Fill = // core lattice zone
    2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
    2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
    2 2 2 2 2 2 2 3 3 3 3 3 3 3 2 2 2 2 2 2 2
    2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2 2
    2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2
    2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2
    2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2
    2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2
    2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2
    2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2
    2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2
    2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2
    2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2
    2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2
    2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2
    2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2
    2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2
    2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 3 2 2 2 2 2
    2 2 2 2 2 2 2 3 3 3 3 3 3 3 2 2 2 2 2 2 2
    2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
    2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2

UNIVERSE 2 // reflector
cell 21 1 mat = 5
cell 22 -1 mat = 5

UNIVERSE 3 lat = 1 pitch = 1.26 1.26 1 scope = 17 17 1 fill = // assembly
    6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
    6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
    6 6 6 6 6 7 6 6 7 6 6 7 6 6 6 6 6
    6 6 6 7 6 6 6 6 6 6 6 6 6 7 6 6 6
    6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
    6 6 7 6 6 7 6 6 7 6 6 7 6 6 7 6 6
    6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
    6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
    6 6 7 6 6 7 6 6 7 6 6 7 6 6 7 6 6
    6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
    6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
    6 6 7 6 6 7 6 6 7 6 6 7 6 6 7 6 6
    6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
    6 6 6 7 6 6 6 6 6 6 6 6 6 7 6 6 6
    6 6 6 6 6 7 6 6 7 6 6 7 6 6 6 6 6
    6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6
    6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6

UNIVERSE 6 move = 0.63 0.63 0   // Fuel rod
cell 3 -1     mat = 1 inner = 1 // Fuel
cell 4 1 & -2 mat = 3 inner = 1 // Air
cell 5 2 & -3 mat = 4 inner = 1 // Zr
cell 6 3      mat = 5           // water

UNIVERSE 7 move = 0.63 0.63 0    // Guide tube
cell 11 -4     mat = 5 inner = 1 // water
cell 12 4 & -5 mat = 4 inner = 1 // Air
cell 13 5      mat = 5           // water

SURFACE
surf 1 cz 0.4096
surf 2 cz 0.4178
surf 3 cz 0.4750
surf 4 cz 0.5690
surf 5 cz 0.6147
surf 10 c/z 224.91 224.91 209 bc = 1 // container

MATERIAL
mat 1 -10.196
    92235.30c 6.9100E-03
    92238.30c 2.2062E-01
    8016.30c 4.5510E-01
mat 3 -0.001
    8016.30c 3.76622E-5
mat 4 -6.550
    40000.60c -98.2
mat 5 9.9977E-02
    1001.30c 6.6643E-02
    8016.30c 3.3334E-02
sab 5 lwtr.60t

CRITICALITY
PowerIter population = 100000 250 500 // keff0 = 1.0
InitSrc point = 224.91 226.17 0

3.5.3. 六边形组件

../_images/lattice_hex_assembly.png

图3.6 六边形组件结构简化示意图

示例3.3 是一个六边形组件的简化输入示例，包含61个六边形燃料栅元。 图3.6 中蓝色为燃 料，黄色为绕线，红色为冷却剂，绿色为反射层。Universe 1为六边形重复结构，其中 Universe 2和Universe 3分别为冷却剂栅元和燃料栅元。由 图3.3 中六边形的排列结构可知， 六边形重复结构的原点在左下角六边形的中心，所以Universe 1需要向x和y方向分别移动 -15和-9.05，才能使得Universe 1的中心与Cell 1的中心重合。


示例3.3 六边形组件输入示例
///////////// MFR ASSEMBLY. FAN Xiao 2012-09-17 /////////////
Universe 0
cell 1 -1 & -2 & -3 & 4 & -5 & -6 & 7 & -8 mat = 0 fill = 1 //Assembly inside
cell 2 16 : -17 : 18 mat = 0 void = 1 //Assembly outside
cell 3 -16 & 17 & -18 & (1 : 2 : 3 : -4 : 5 : 6 : -7 : 8) mat = 5 //reflector

niverse 1 move=-15 -9.05 0 lat=2 pitch=2 2.06787 scope=11 11 sita=63.435 fill=
    2 2 2 2 2 2 2 2 2 2 2
    2 2 2 2 2 3 3 3 3 3 2
    2 2 2 2 3 3 3 3 3 3 2
    2 2 2 3 3 3 3 3 3 3 2
    2 2 3 3 3 3 3 3 3 3 2
    2 3 3 3 3 3 3 3 3 3 2
    2 3 3 3 3 3 3 3 3 2 2
    2 3 3 3 3 3 3 3 2 2 2
    2 3 3 3 3 3 3 2 2 2 2
    2 3 3 3 3 3 2 2 2 2 2
    2 2 2 2 2 2 2 2 2 2 2

Universe 2
cell 21 -15 mat = 1
cell 22 15 mat = 1

Universe 3
cell 31 -15 mat = 2
cell 32 15 mat = 3

Surface
surf 1 py 8.4
surf 2 p 1.732 1 0 16.3
surf 3 p 1.732 -1 0 16.3
surf 4 py -8.4
surf 5 p -1.732 -1 0 16.3
surf 6 p -1.732 1 0 16.3
surf 7 pz -30
surf 8 pz 30
surf 15 cz 0.975
surf 16 cz 30
surf 17 pz -35
surf 18 pz 35

Material
mat 1 -0.8139 // Na
    11023.30c 1.0
mat 2 -10.41 // UO2
    92235.30c -56.5 92238.30c -31.1 8016.30c -12.3 13027.30c -0.02
    20000.60c -0.02 12000.60c -0.02 26000.55c -0.02 14000.60c -0.02
mat 3 -0.8355 // wiry
    11023.30c 2.132E+0 28000.50c 3.223E-3
    24000.50c 4.759E-3 26000.55c 1.634E-2
mat 5 0.1236 // Be9
    4009.30c 1

Criticality
PowerIter keff0 = 1.0 population = 2000 50 300
InitSrc point = 0 0 0

3.5.4. 六边形堆芯

../_images/lattice_hex_core.png

图3.7 六边形堆芯结构简化示意图

示例3.4 是一个六边形堆芯几何模块的简化输入示例。 图3.7 中蓝色为燃料，黄色为绕线，红 色为冷却剂。简单起见，堆芯输入文件仅包含同一类的7个燃料组件，真实堆芯可以自行增 加组件，这里只重点介绍几何变换填充过程。Universe 1是六边形重复结构，其中 Universe 2是冷却剂，Universe 3是六边形燃料组件。Universe 3本身又包含六边形重复 结构排列的燃料栅元（类似于 示例3.3 中的描述），即， 整个系统用到了两层嵌套的六边形 重复结构。在描述该几何结构时，首先把燃料和冷却剂（Universe 4和Universe 5）填充 到六边形栅元重复结构（Universe 3）中；接下来把六边形栅元重复结构（Universe 3）在 x-y平面内逆时针旋转90°（rotate = 0 1 0 -1 0 0 0 0 1）并平移 （move = 9.05, -15, 0），填充到六边形组件重复结构（Universe 1）；接下来再把六边 形重复结构（Universe 1）平移（move = -50.4 -27.942 0）到堆芯（即Cell 1）的中心位 置并填充。需要注意的是，用户在做几何变换时应先做旋转再做平移。


示例3.4 六边形堆芯几何输入示例
//// MFR CORE.FAN Xiao 2012-09-17 ////
Universe 0
cell 1 -21&7&-8 mat=0 fill=1
cell 2 21:-7:8 mat=0 void=1

Universe 1 move=-50.4 -27.942 0 lat=2 pitch=16.8 16.302 scope=5 5 sita=60 fill=
    2 2 2 2 2
    2 2 3 3 2
    2 3 3 3 2
    2 3 3 2 2
    2 2 2 2 2

Universe 2
cell 21 -15 mat=1
cell 22 15 mat=1

Universe 3 move=9.05 -15 0 rotate=0 1 0 -1 0 0 0 0 1
lat=2 pitch=2 2.06787 scope=11 11 sita=63.435 fill=
    4 4 4 4 4 4 4 4 4 4 4
    4 4 4 4 4 5 5 5 5 5 4
    4 4 4 4 5 5 5 5 5 5 4
    4 4 4 5 5 5 5 5 5 5 4
    4 4 5 5 5 5 5 5 5 5 4
    4 5 5 5 5 5 5 5 5 5 4
    4 5 5 5 5 5 5 5 5 4 4
    4 5 5 5 5 5 5 5 4 4 4
    4 5 5 5 5 5 5 4 4 4 4
    4 5 5 5 5 5 4 4 4 4 4
    4 4 4 4 4 4 4 4 4 4 4

Universe 4
cell 41 -15 mat=1
cell 42 15 mat=1

Universe 5
cell 51 -15 mat=2
cell 52 15 mat=3

Surface
surf 5 p 1 1.6632 0 46.474
surf 2 p 1 -1.6632 0 46.474
surf 3 p -1 -1.6632 0 46.474
surf 6 p -1 1.6632 0 46.474
surf 1 px 27.942
surf 4 px -27.942
surf 7 pz -30
surf 8 pz 30
surf 15 cz 0.975
surf 21 cz 25

Material
mat 1 -0.8139 // Na
    11023.30c 1.0
mat 2 -10.41 // UO2
    92235.30c -56.5 92238.30c -31.1 8016.30c -12.3 13027.30c -0.02
    20000.60c -0.02 12000.60c -0.02 26000.55c -0.02 14000.60c -0.02
mat 3 -0.8355 // wiry
    11023.30c 2.132 28000.50c 3.223E-3 24000.50c 4.759E-3 26000.55c 1.634E-2

Criticality
PowerIter keff0=1.0 population = 50000 200 1000
InitSrc point=0 0 0

3.5.5. 随机介质模型

示例3.5 是一个随机介质单棒模型（显式建模法），随机介质内燃料颗粒体积份额为0.32。 对于单棒模型，周围介质为水，单棒包壳外表面半径，内表面半径，燃料棒半径分别为0.7116,0.6546 0.6461cm。

下表给出了所填充TRISO球的几何信息。

表3.2 TRISO颗粒尺寸 :name: triso_sizes

Layer

Radius/cm

Density/ g/cm3

Material

Layer1

0.0450

12.95

UC

Layer2

0.0525

1.05

C

Layer3

0.0555

1.9

C`

Layer4

0.0590

3.18

SiC

Layer5

0.0610

1.9

C

Outside

3.18

SiC

示例3.5 随机介质组件模型（隐式建模法）
  ///////////// SingleRod Explicit Model PF=0.32 /////////////
  Universe 0
  cell 1 -13:14:-15:16:-17:18 mat=0 void=1                           // Outside the Assembly
  cell 3 13&-14&15&-16&17&-18  fill=2                                // Inside the Assembly

  Universe 2                        //Fuel Rods
  cell 21 -30&17&-18    fill=5
  cell 22 30&-31&17&-18  mat=2      // Helium Fill
  cell 27 31&-32&17&-18  mat=3      // Cladding FeCrAl
  cell 23 32&17&-18      mat=4

  // TRISO Particles distribution using explicit model
  Universe 5   lat = 4  MATRIC = 7 move = -0.6461 -0.6461 -176.5
               PARTICLE = 12
               PF = 0.32
               RSA = 1
               RAD = 0.061
               TYPE = 2
               SIZE = 0.6461 353

  Universe 7
  cell 66 -49    mat = 7

  Universe 12
  cell 60 -44 mat=1  vol=0.000381704      tmp=900
  cell 61 44&-45 mat= 5    tmp=900
  cell 62 45&-46 mat= 6    tmp=900
  cell 63 46&-47 mat= 7    tmp=900
  cell 64 47&-48 mat= 8    tmp=900
  cell 65 48 mat=7         tmp=900

  Surface
  surf 30 cz 0.6461
  surf 31 cz 0.6546
  surf 32 cz 0.7116
  surf 13 py  -0.90                   bc=1
  surf 14 py   0.90                   bc=1
  surf 15 px  -0.90                   bc=1
  surf 16 px   0.90                   bc=1
  surf 17 pz  -176.5 bc=1
  surf 18 pz   176.5 bc=1
  surf 44 so  0.0450
  surf 45 so  0.0525
  surf 46 so  0.0555
  surf 47 so  0.0590
  surf 48 so  0.0610
  surf 49 inf
  surf 54 so  0.0130
  surf 55 so  0.0385
  surf 56 so  0.0425
  surf 57 so  0.0460
  surf 58 so  0.0500

  Material
  mat 1 -12.95                      //UC
    92235.90c 16.10097657
    92238.90c 83.89902343
    6000.90c  100
  mat 2 -0.0022                    // Helium
    2004.90c  1.0
  mat 3  1.6701076E-02            //FeCrAl
    26056.60c  1.22593E-02
    24052.60c  3.55342E-03
    13027.60c  8.88356E-04
  mat 4 -0.72    // Water
    8016.60c  1.0
    1001.60c  2.0
  sab 4  lwtr.62t
  mat 5 -1.05
    6000.90c    1.0
  //sab 5 grph.65t
  mat 6 -1.9
    6000.90c    1.0
  //sab 6 grph.65t
  mat 7 -3.18                   //SiC
    6000.90c     1.0
    14028.90c     1.0
  //sab 7 grph.65t
  mat 8 -1.9
    6000.90c    1.0
  //sab 8 grph.65t
../_images/TRISO_Particle.png

图3.8 TRISO球截面几何

../_images/Fuelrod.png

图3.9 燃料棒模型示意图

图3.8 和 图3.9 为 示例3.5 的几何模型示意图。

示例3.6 是随机介质组件模型（隐式建模法），由上述单棒模型9*9棒束阵列构成


示例3.6 随机介质组件模型（隐式建模法）
  ///////////// Array9*9 Implicit Model PF=0.32 /////////////
  Universe 0
  cell 1 -1:2:-3:4    mat=0 void=1                                 // Outside the Assembly
  cell 2 1&-2&3&-4    fill=1                                       // Inside the Assembly

  Universe 1 lat=1 pitch=1.8 1.8 1 scope=9 9 1 fill=
   2*81


  Universe 2  move=0.9 0.9 0       //Fuel Rods
  cell 3  32      mat=4            //water
  cell 21 -30     fill=5
  cell 22 30&-31  mat=2            //Helium Fill
  cell 27 31&-32  mat=3            //Cladding FeCrAl


  // TRISO Particles distribution using implicit model
  Universe 5   lat = 3  MATRIC = 7
               PARTICLE = 12
               PF = 0.32
               RAD = 0.061


  Universe 7
  cell 66 -49    mat = 7

  Universe 12
  cell 60 -44 mat=1  vol=0.000381704      tmp=900
  cell 61 44&-45 mat= 5    tmp=900
  cell 62 45&-46 mat= 6    tmp=900
  cell 63 46&-47 mat= 7    tmp=900
  cell 64 47&-48 mat= 8    tmp=900
  cell 65 48 mat=7         tmp=900

  Surface
  surf 1  px  0                       bc=1
  surf 2  px  16.2                    bc=1
  surf 3  py  0                       bc=1
  surf 4  py  16.2                    bc=1
  surf 30 cz 0.6461
  surf 31 cz 0.6546
  surf 32 cz 0.7116
  surf 13 py  -0.90                   bc=1
  surf 14 py   0.90                   bc=1
  surf 15 px  -0.90                   bc=1
  surf 16 px   0.90                   bc=1
  surf 17 pz  -176.5 bc=1
  surf 18 pz   176.5 bc=1
  surf 44 so  0.0450
  surf 45 so  0.0525
  surf 46 so  0.0555
  surf 47 so  0.0590
  surf 48 so  0.0610
  surf 49 inf

  Material
  mat 1 -12.95                     //UC
    92235.90c 16.10097657
    92238.90c 83.89902343
    6000.90c  100
  mat 2 -0.0022                    //Helium
    2004.90c  1.0
  mat 3  1.6701076E-02             //FeCrAl
    26056.60c  1.22593E-02
    24052.60c  3.55342E-03
    13027.60c  8.88356E-04
  mat 4 -0.72                      //Water
    8016.60c  1.0
    1001.60c  2.0
  sab 4  lwtr.62t
  mat 5 -1.05
    6000.90c    1.0
  //sab 5 grph.65t
  mat 6 -1.9
    6000.90c    1.0
  //sab 6 grph.65t
  mat 7 -3.18                      //SiC
    6000.90c     1.0
    14028.90c     1.0
  //sab 7 grph.65t
  mat 8 -1.9
    6000.90c    1.0
  //sab 8 grph.65t