[发明专利]一种飞行器湍流绕流模拟的壁面距有效单元计算方法

权利要求书

1.一种飞行器湍流绕流模拟的壁面距有效单元计算方法，其特征在于，包括以下步骤：

S1：读入并存储计算网格单元格点坐标与边界条件；

S2：壁面单元预处理；

S2-1：收集壁面单元；

S2-2：确定壁面单元外接长方体的个数：根据计算网格壁面单元数，确定循环盒子法中使用的壁面单元外接长方体个数，壁面单元外接长方体是指包含其内壁面单元的最小长方体；

S2-3：根据所有壁面单元建立初始外接长方体，外接长方体所需记录的信息包括：外接长方体的笛卡尔坐标范围、所含壁面单元数量、所含壁面单元编号；

S2-4：二分壁面点最多的外接长方体；

将包含最多壁面单元的外接长方体沿其最长边二等分为两个外接长方体，并更新这两个外接长方体所需记录的信息，重复步骤S2-4，直至外接长方体的总数达到步骤S2-3确定的个数；

S3：计算壁面距；

S3-1：壁面距归零：分配所有计算网格单元存储壁面距的空间，并将其初值赋为0；遍历所有计算网格单元，对任一单元执行步骤S3-2；

S3-2：确定壁面距计算方式；

若计算网格单元紧邻单元的壁面距均小于设定阈值，则精确计算壁面距，继续步骤S3-3；

若计算网格单元紧邻单元的壁面距大于设定的阈值，当计算网格单元位于其紧邻单元的壁面距减小方向时，则精确计算壁面距，继续步骤S3-3；

当计算网格单元位于其紧邻单元的壁面距增大方向时，则粗略估计壁面距，跳转至步骤S3-5；

S3-3：搜索距离最近的2个外接长方体；

令x_i,c为计算网格单元的格心坐标，x_i,min、x_i,max分别表示外接长方体在x_i方向的最小值和最大值，则计算网格单元到外接长方体的距离d₀为：

（1）

遍历所有外接长方体，根据式（1）筛选出距离最小的2个外接长方体，并对这2个外接长方体执行步骤S3-4；

S3-4：精确计算壁面距；

S3-5：根据梯度估计壁面距：利用单侧差分计算已赋值紧邻单元处的壁面距梯度，根据该梯度、当前网格单元及其紧邻单元的坐标信息，估计当前网格单元的壁面距；

S4：输出结果。

2.根据权利要求1所述的一种飞行器湍流绕流模拟的壁面距有效单元计算方法，其特征在于，所述步骤S2-1包括以下步骤：

S2-1-1：根据步骤S1读入的边界条件，确定计算网格单元中的壁面单元；

S2-1-2：为壁面单元和壁面点编号，确保壁面单元编号与其一尚未编号的顶点使用同一编号以便搜索；

S2-1-3：分别按壁面点编号存储壁面点坐标、按壁面单元编号存储壁面单元与其顶点的对应关系、按壁面单元编号存储其相邻单元的编号；

S2-1-4：遍历壁面单元，计算并存储各壁面单元的单位法向。

3.根据权利要求1所述的一种飞行器湍流绕流模拟的壁面距有效单元计算方法，其特征在于，所述步骤S3-4包括以下步骤：

依据计算网格单元与外接长方体的位置信息，分为两种情况：当计算网格单元所有格点均位于外接长方体内时，执行步骤S3-4-1；当计算网格单元存在不位于外接长方体内的格点时，执行步骤S3-4-2；

S3-4-1：遍历外接长方体内所有壁面单元，逐一计算当前网格单元格点到壁面单元的距离，并存储其中的最小值；

S3-4-2：基于矢量法跳跃式访问外接长方体内壁面单元，在外接长方体内任选一壁面单元，计算并存储当前网格单元到其的距离；计算壁面单元面心指向当前网格单元格心的向量在该壁面单元上的投影方向，该投影方向所指紧邻壁面单元即为下一个计算的壁面单元；跳跃式访问当前外接长方体内壁面单元，直至投影方向指向当前外接长方体的边界。