mcnp005曲面描述卡

曲面描述卡

用方程定义曲面

格式:  j  n  a  list

j:  曲面的题目编号，写在卡片的第l-5列，l≤j≤99999。可不按顺序填写编号。题目中出现的任何曲面都允证件反射界面。对于反射曲面，只需在它的曲面描述卡第1-5列的曲面序号前加上一个“*”，当粒子打在这样的曲面上时，便按镜面反射后继续输运。有反射曲面的问题，其计数需要各不相同地归一，反射曲面号不能多于4位数字。

n:  相应于TRn卡的整数n它表示该卡定义曲面是在辅助坐标系下描述的，而这个辅助坐标系与基本坐标系之间的关系由TRn卡给出。当曲面直接在基本坐标系下描述时，则应略去该项整数。

a:  曲面助记名（见表2-1）。

list：按照规定次序给出的描述该曲面的数据项（1-10项）。

表2-l列出了曲面的助记名、类型、方程及应给出的数据项次序。在曲面卡上要用方程来描述一个曲面，首先要从表2-1中找出所需的曲面，然后根据其方程的特定形式算出所需的系数，按照上面的书写格式将它们依序写在一张卡片上。

如果一个曲面在点（x,y,z）的计算值是正的，则称该点对于这个曲面具有正的坐向。一个曲面的表达式是表2-l中曲面方程的左边。位于球、柱、锥及环的正坐向是这些曲面的外边。对于垂直坐标轴的平面PX、PY、PZ具有正的尘向的点，分别是x，y或z值大于相应平面截距的那些点。对P、SQ及GQ曲面，用户提供这个表达式的全部系数，于是就能随意确定这个曲面的坐向。这与其它情况不同，其它情况曲面的坐向即使是任意的，也是由表达式的形式唯一确定的。因此，通过一个曲面坐标变换，PX、PY或PZ曲面可能必须由P曲面代替以防止这个曲面的坐向相反。

除去作虚设界面用的曲面外，在有探测器记数或DXTRAN球的情况下，不应使用反射界面。与栅元描述卡一样，曲面描述卡上给出的是曲面的题目编号，MCNP将按照栅元卡的同样方式给出曲面的程序编号。

例2-8  在Y=3处垂直于Y轴的平面，其曲面描述卡为：

             j  py  3

例2-9   j  gq  l  .25  .75  0  -.866  0  -12  -2   3.464   39

这是半径为1cm的圆柱面，其对称轴在垂直与X轴的一个平面上，该平面在X轴上的截距为6cm。与X轴相距2cm的地方绕X轴从Y轴向Z轴旋转30o角。在柱外的点具有正的坐向。先要算出这十项数据并不是很简单的事情。我们可以用另一种方式给出。建立一个辅助坐标系，将柱面的对称轴定义为辅助坐标系的X轴。那么该柱面在辅助坐标系下描述是很简单的。然后再相应给出一张TRn卡，用来定义基本坐标系与辅助坐标系的关系，MCNP将实现曲面描述的转换，将其从辅助坐标系转换成基本坐标系下的描述形式。在辅助坐标系下这个曲面描述如下：

                 j       7   cx   1

                *TR7    6   1   -1.732   0  30  60

用几何点定义曲面

除表2-l中列出的26种曲面外，还可以用助记名为X、Y、Z的三种曲面卡，它们是用几何坐标点来描述曲面的。用这些卡描述的曲面必须是分别关于X轴、Y轴或Z轴对称的，且是唯一的、真正的、连续的曲面。

格式： j    n  a   list

j:  曲面号，1≦j≦99999

n:  TRn卡的号，没有坐标变换缺省

a:  字母X、Y或Z

list:  一至三对坐标。类型为 X、Y及 Z的曲面卡是用坐标点描述曲面而不是用方程系数来描述。用这些卡描述曲面必须是分别关于X、Y、Z轴对称，并且如果这个曲面是由一叶以上组成的，则指定的坐标点必须全都在同一个叶上。

每一对坐标点定义该曲面上的一个几何点。例如在一张Y卡上可以给出：

     j   Y  y1  rl  y2  r2

这里

是点i的柱坐标。给出的坐标对数不同，描述的曲面类型也不同：

（1）给出一对坐标，则定义的是平面，即PX、PY或PZ。

（2）给出两对坐标，定义的是线性曲面，如PX，PY，PZ，CX，CY，CZ，KX，KY或 KZ。

（3）给出三对坐标，则定义一个二次曲面，包括 PX，  PY，  PZ， SO，  SX，  SY，SZ， CX， CY， CZ， KX， KY， KZ或 SQ。

当用两个点定义一个锥面时，只生成一个单叶锥面。

用这个规定所确定的这些曲面（SQ除外）的坐向与由方程指定曲面得到的坐向是等同的。而对 SQ来讲，是这样定义坐向的，即离对称轴足够远的点有正的坐向。这与方程定义的 SQ不同；用方程定义的 SQ可自由选取坐向。

例2-10  曲面卡

      j      X   7   5   3   2   4   3

描述的是关于X轴对称的曲面，该曲面将通过三个点:（x,r）=（7,5）、（3,2）及（4,3）。这个曲面是有两个叶的双曲面。MCNP将对此曲面进行转换，给出标准形式：

      j    SQ  -.083333  1  1  0  0  0  68.52083  -26.5  0  0

例2-11 这一例子定义了一个由一叶锥面、双曲面及一个椭球面界定的栅元。这三个曲面定义的是关于y轴对称的环形栅元，该栅元的截面如图2-10所示。第一个曲面通过点（-3,2）及（2,1），第二个曲面通过点（2,3）、（3,3）及（4,2）。由点（2,1）、（3,1）和（4,2）定义最后一个曲面．这些几何点坐标均指的是(y,r)，那么曲面卡便是：

        1   Y  -3  2  2  1

        2   Y   2  3  3  3  4  2

        3   Y   2  1  3  1  4  2

MCNP指出曲面1是一叶锥面，曲面2是一个椭球，曲面3是一叶双曲面。当使用PRINT输入卡或执行行选择时，打印出各种曲面的标准方程的系数。例如，查看输出文件，便可看到曲面3经MCNP转换成标准曲面形式：

        3   SQ  1  -1.5  1  0  0  0  -0.625  0  2.5  0

由三个点定义一般平面

格式： j   n  P  x1  y1  z1  x2  y2  z2  x3  y3  z3

j：  曲面号，1≤j≤99999

n：  TRn卡的号，没有坐标变换缺省；

(xi,yi,zi)： 定义该平面的点坐标。

MCNP对用户指定的P型曲面，将检查所给的数据个数，若是4项，则作为一般斜置平面方程的系数理解；若多于4项时，便作为三维空间点的坐标值理解，每三个数定义一个空间点。由MCNP把它们转换成平面方程系数以产生平面：

                  AX＋BY＋CZ-D = 0                                

所产生的平面方程系数遵循如下原则：

 （1）应使坐标原点关于该平面有负坐向；

 （2）当该平面通过坐标原点时（D=0），则应使得点（0,0,∞）对该平面有正坐向；

 （3）若上两项都无法做到（即D=C=0），则应使点（0,∞,0）对该平面有正坐向；

   （4）若上三项都无法做到（即D=C=B=0），则应使点（∞,0,0）对该平面有正坐向；       倘若第四项也失败，那么三个点给在一条直线上，则将印出警告信息。