一种应用于屏蔽计算曲面衔接中的S

技术特征：
1.一种应用于屏蔽计算曲面衔接中的s
n-mc耦合计算方法，其特征在于：包括如下步骤：步骤1：离散纵标方法s
n
输运计算得到反应堆几何结构包含曲面面源区域在内的网格角注量率；步骤2：根据s
n
的z方向分层将曲面划分为环带，计算环带投影与s
n
网格面的交点，根据交点数量将环带等分为相同数量的环面网格，进而获得环面上网格与s
n
网格的映射关系；具体映射关系用环面网格编号与s
n
网格编号来对应：环面网格编号格式，s
n
网格编号格式，其中与一一对应，号通过公式映射得到：（1）式中：——某一环带对应的分层号；——在某一环带上的环面网格编号；——s
n
输运计算划分在方向上划分的网格编号；——s
n
输运计算在x方向划分的下边界、上边界网格号；——s
n
输运计算在y方向划分的下边界、上边界网格号；步骤3：获得网格映射关系后，将对应s
n
网格中心点的角注量率视为从环面网格处出射的角注量率；同时由于仅考虑沿面源正法线方向的出射流，需要计算曲面上不同位置的法向量：根据步骤2中计算的环带投影与s
n
网格面的交点以及曲面的中心坐标获得交点的法向量：球面某点法向量：，椭球面某点法向量：，其中：——曲面与s
n
网格的交点坐标；——球面和椭球面的球心坐标；——椭球面的三个方向的轴长；将法向量和s
n
计算的离散方向m进行向量点积运算，结果为正判定为正向出射；步骤4：根据角注量率依次求得方向m、网格q、分层k、能群g、面s、粒子类型p上的出射流；曲面计算出射流的流程与平面、柱面类似，有以下三点不同点：
（1）通过对同一分层上网格出射流的求和，获得某一分层出射粒子流；（2）对曲面而言，法向量随网格位置变动，需根据粒子飞行方向与不同位置的法向量的点积判断是否为正向出射；（3）计算从某一网格出射的粒子流时，网格面积应为曲面上的网格面积，推导得到的公式如下：球面网格：（2）椭球面网格：（3）（4）
ꢀꢀꢀ
（5）（6） （7）式中：——曲面某一分层上某网格的面积；——椭球壳表面积；——某一分层的高度；——椭球面三个方向的轴长；——某一分层上的曲面网格数目；
——横向截面的椭圆周率，；——与横向截面垂直的纵向椭圆的椭圆周率；——椭球面在不同高度下的缩放系数；——极坐标系下某点所处位置的极角；步骤5：根据粒子出射流，计算得到方向m、网格q、分层k、能群g、面s、粒子类型p的累积分布函数cdf；步骤6：利用累积分布函数cdf进行源抽样，进行mc续算；为确保抽样获得的粒子几何位置能均布在球或椭球的表面，先根据分层k的累积分布函数随机抽样得到高度z，根据高度z确定截面方程，对于球面而言，其在截面圆网格边界的分布与圆心角成正比，对于椭球而言，为简化抽样过程，近似认为其在截面椭圆网格边界的分布与圆心角成正比，具体几何位置的抽样公式如下：球面：（8）（9）（10） （11）
ꢀꢀ
（12）椭球面：（13）（14） （15）（16）
（17）
ꢀꢀꢀ
（18）式中：——几何位置抽样时的粒子位置坐标；——椭球面三个方向的轴长；——球面或椭球面某一分层z处截面的高度；、 ——球面或椭球面某一高度分层的下边界，上边界；、——某一截面椭圆的长轴、短轴；——抽样时产生的随机数，在之间；——圆心角，；——曲面网格极坐标系下对应的极角下边界、上边界；——球面的半径；——某一截面圆的半径。

技术总结
一种应用于屏蔽计算曲面衔接中的S


技术研发人员：王永平 张运升 郑友琦 杜夏楠 吴宏春 曹良志
受保护的技术使用者：西安交通大学
技术研发日：2022.09.28
技术公布日：2022/11/18