[發明專利]節點式非結構網格的邊界修正方法、裝置、設備及介質

有效

• 申請號: 202310525517.1
• 申請日: 2023-05-11
• 公開（公告）號: CN116245049B
• 公開（公告）日: 2023-07-25
• 發明（設計）人: 于廣瀛;李敏;蘇利天;李彬;查浩;王棟志;劉麗麗;何旭
• 申請（專利權）人: 中國空氣動力研究與發展中心計算空氣動力研究所
• 主分類號: G06F30/28
• 分類號: G06F30/28;G06T17/20;G06F113/08;G06F119/14;G06F119/02
• 代理公司: 北京集佳知識產權代理有限公司 11227
• 代理人: 牛亭亭

1.一種節點式非結構網格的邊界修正方法，其特征在于，包括：

針對非結構網格中的每個邊界節點，從所述非結構網格中確定出與所述邊界節點相鄰的各個內節點；所述非結構網格與發動機燃燒室對應；

從各個所述內節點中選擇出與所述邊界節點的邊界面法向最近的目標內節點；

基于所述目標內節點的每個相鄰節點與所述目標內節點的梯度影響因子對所述邊界節點的控制方程的系數進行修正，以便基于修正后的控制方程計算所述邊界節點處流體的物理量；其中：

所述基于所述目標內節點的每個相鄰節點與所述目標內節點的梯度影響因子對所述邊界節點的控制方程的系數進行修正，包括：

針對所述目標內節點的每個相鄰節點，計算所述目標內節點與所述相鄰節點的梯度加權影響系數；

基于每個所述相鄰節點的梯度加權影響系數，得到與所述邊界節點相關的修正系數；

基于所述邊界節點相關的修正系數對所述邊界節點的控制方程的系數進行修正。

2.根據權利要求1所述的節點式非結構網格的邊界修正方法，其特征在于，所述從各個所述內節點中選擇出與所述邊界節點的邊界面法向最近的目標內節點，包括：

針對每個所述內節點，計算所述內節點和所述邊界節點之間連線與所述邊界節點的邊界面法向之間的夾角；

將夾角最小的內節點作為目標內節點。

3.根據權利要求1所述的節點式非結構網格的邊界修正方法，其特征在于，所述基于每個所述相鄰節點的梯度加權影響系數，得到所述邊界節點相關的修正系數的過程為：

基于每個所述相鄰節點的梯度加權影響系數結合所述邊界節點的線性方程，得到所述邊界節點的修正系數，其中，所述線性方程為：

φ_b表示邊界節點b的物理量，a_b表示邊界節點b的修正系數，φ_p表示目標內節點p的物理量，a_p表示目標內節點p的修正系數，φ_i表示目標內節點p的相鄰節點i的物理量，a_i表示目標內節點p的相鄰節點i的修正系數；其中：

其中，表示邊界節點b對于目標內節點p的系數向量，w_b表示邊界節點b的權重系數，表示目標內節點p到邊界節點b的方向向量，w_i表示目標內節點p的相鄰節點i的權重系數，表示目標內節點p的相鄰節點i的系數向量，其中，目標內節點p的相鄰節點i不為邊界節點b，N表示目標內節點p周圍的相鄰節點的總數量。

4.根據權利要求1或3所述的節點式非結構網格的邊界修正方法，其特征在于，所述基于所述邊界節點相關的修正系數對所述邊界節點的控制方程的系數進行修正，包括：

基于所述邊界節點相關的修正系數對隱式求解矩陣中所述邊界節點對應行的各元素系數進行修正。

5.一種節點式非結構網格的邊界修正裝置，其特征在于，包括：

第一確定模塊，用于針對非結構網格中的每個邊界節點，從所述非結構網格中確定出與所述邊界節點相鄰的各個內節點；所述非結構網格與發動機燃燒室對應；

選擇模塊，用于從各個所述內節點中選擇出與所述邊界節點的邊界面法向最近的目標內節點；

修正模塊，用于基于所述目標內節點的每個相鄰節點與所述目標內節點的梯度影響因子對所述邊界節點的控制方程的系數進行修正，以便基于修正后的控制方程計算所述邊界節點處流體的物理量；其中：

所述修正模塊，包括：

第二計算單元，用于針對所述目標內節點的每個相鄰節點，計算所述目標內節點與所述相鄰節點的梯度加權影響系數；

第三計算單元，用于基于每個所述相鄰節點的梯度加權影響系數，得到與所述邊界節點相關的修正系數；

修正單元，用于基于所述邊界節點相關的修正系數對所述邊界節點的控制方程的系數進行修正。