[发明专利]一种根据流场数值计算结果确定激波位置的方法

摘要

本发明公开了一种根据流场数值计算结果确定激波位置的方法，对于待判断的当前网格，查找该当前网格压力梯度方向的上游和下游两个网格，通过对这两个网格内参数的比较，判断当前网格是否应当标记为激波核心；若当前网格标记为激波核心，则将当前网格周围的两圈网格标记为激波边缘。本发明提供的根据流场数值计算结果确定激波位置的方法，不需要针对特定流场调整阈值即可实现激波位置的确定，因此适用于完全自动化的分析；本发明所确定的激波是一片网格区域而不是无厚度的线，因此有利于统计激波造成的参数变化。

主权项

1.一种根据流场数值计算结果确定激波位置的方法，其特征在于：对于待判断的当前网格，查找该当前网格压力梯度方向的上游和下游两个网格，通过对这两个网格内参数的比较，判断当前网格是否应当标记为激波核心；若当前网格标记为激波核心，则将当前网格周围的两圈网格标记为激波边缘；该方法具体包括如下步骤：(1)对于当前网格c0，从周围相邻的8个网格中查找出当前网格c0压力梯度方向的上游和下游两个网格，上游网格记为cup，下游网格记为cdown；将当前网格c0压力梯度向量记为上游网格为8个网格中网格中心指向c0中心的向量与的夹角最小的网格，下游网格为8个网格中网格中心指向c0中心的向量与的夹角最大的网格；上游网格的速度向量记为下游网格的速度向量记为(2)判断上游网格和下游网格内的参数是否满足以下不等式： V → u p · ▿ → p ≥ a u p | ▿ → p | V → d o w n · ▿ → p ≥ a d o w n | ▿ → p | ρ u p ≤ ρ d o w n ]]>式中：aup和adown表示上游网格和下游网格的声速，ρup和ρdown表示上游网格和下游网格的密度若满足上述不等式，则判断当前网格c0为激波核心，并将前网格c0周围的两圈网格标记为激波边缘。