[发明专利]一种基于非结构网格的流动过程数值离散方法

权利要求书

1.一种基于非结构网格的流动过程数值离散方法，其特点在于，包括：

S1：对待求解域进行非结构网格的划分，得到三角形线性单元和对应的顶点，对所述顶点进行编号；

S2：选取所述编号为i的顶点为待处理顶点，并将与所述待处理顶点直接相连的节点作为待离散体，所述待离散体中包含N个三角形线性子单元，将所述N个三角形线性子单元的外心连接，构建与所述待处理顶点对应的控制体；

S3：所述三角形线性子单元将所述控制体划分为N个四边形子控制体；

S4：关联一个所述子控制体与对应的稳态扩散的控制方程的积分域，并采用中点积分近似法将所述稳态扩散的控制方程近似为离散扩散方程；

S5：关联一个所述子控制体与对应的稳态对流扩散控制方程，并采用中点积分近似法将所述稳态对流扩散控制方程近似为离散对流扩散方程；

S6：关联一个所述子控制体与对应的包含体积贡献的稳态扩散方程，并将所述离散扩散方程和所述离散对流扩散方程带入所述包含体积贡献的稳态扩散方程，离散为第一体积贡献稳态扩散离散方程。

2.根据权利要求1所述的基于非结构网格的流动过程数值离散方法，其特征在于，所述对所述顶点进行编号具体包括：

若所述顶点为求解域的非边界节点，确定与所述非边界节点直接相连接的从顶点，并按照逆时针方向对所述从顶点进行排序；

若所述顶点为求解域的边界节点，确定与所述边界节点直接相连接的从顶点，并按照逆时针方向对所述从顶点进行排序，其中，排序的起点为边界节点，且所述逆时针方向的连线从所述求解域内穿过；

对边界进行分段，所述分段后得到的各个段是连续的，且按逆时针排列，将所述段的所述边界节点进行局部坐标编号，并按逆时针方向进行排列。

3.根据权利要求2所述的基于非结构网格的流动过程数值离散方法，其特征在于，所述步骤S3之后，步骤S4之前还包括：

确定一个所述四边形子控制体中与所述待处理顶点非直接相连的两条边分别为第一求解边和第二求解边；

所述三角形线性子单元中与所述待处理顶点直接相连的垂直平分线将所述三角形线性子单元中的非四边形子控制体部分划分为第一求解体和第二求解体，所述第一求解体和所述第二求解体为四边形，所述第一求解体中包含所述第一求解边，所述第二求解体中包含所述第二求解边。

4.根据权利要求3所述的基于非结构网格的流动过程数值离散方法，其特征在于，所述三角形线性子单元中与所述待处理顶点直接相连的垂直平分线将所述三角形线性子单元中的非四边形子控制体部分划分为第一求解体和第二求解体，所述第一求解体和所述第二求解体为四边形，所述第一求解体中包含所述第一求解边，所述第二求解体中包含所述第二求解边之后还包括：

确定沿着所述第一求解体的非第一求解边的逆时针方向移动时对应横坐标和纵坐标的变化值分别为和确定沿着所述第二求解体的非第二求解边的外边沿逆时针方向移动时对应横坐标和纵坐标的变化值分别为和其中，f1表示第一求解体，f2表示第二求解体。