用于快速計算飛行器高超聲速粘性力的方法和裝置制造方法
【專利摘要】本發明公開了一種用于快速計算飛行器高超聲速粘性力的方法和裝置。所述方法包括：步驟S1.對無粘Euler方程采用計算流體動力學CFD方法獲得所述飛行器高超聲速的表面無粘流場物理參數，其中，所述表面無粘流場物理參數包括：表面的壓力、表面的絕熱溫度和表面的速度；步驟S2.獲取所述飛行器高超聲速的表面的所有表面單元的雷諾數；步驟S3.根據所述表面無粘流場物理參數和所述表面的所有表面單元的雷諾數得到所述飛行器高超聲速的粘性力。所述裝置包括：獲得裝置、獲取裝置和得到裝置。本發明通過上述技術方案提高了工程初期設計中飛行器高超聲速粘性力的計算效率和計算精度，降低了工程使用成本。
【專利說明】用于快速計算飛行器高超聲速粘性力的方法和裝置
【技術領域】
[0001]本發明屬于空氣動力學【技術領域】，特別涉及一種用于快速計算飛行器高超聲速粘性力的方法和裝置。
【背景技術】
[0002]隨著計算機運算速度的大幅提高和CFD (Computational Fluid Dynamics,計算流體動力學)方法的日趨完善，基于無粘Euler方程的CFD方法已經能夠精確、迅速地完成無粘流場特性的計算。而基于粘性的完全Navier-Stokes方程CFD方法獲得更高精度和粘性底層流動特性仍然需要付出很大的代價:繁雜的高質量網格生成、流場完全收斂的極大耗時等。因此在復雜飛行器外形初期選型設計中，流場特性的預測主要還是采用高效率并具有良好精度的無粘Euler方程CFD方法，而粘性力部分通常采用與無粘流場結果完全獨立的工程經驗修正方法計算，如等效平板層流Blassius近似方法和瑞流Van Driest II方法等。這些方法計算效率高，工程應用于一定馬赫數和雷諾數范圍內具有較好的預測能力和準度，但是對于高超聲速粘性流動的預測精度很差。國外在高超聲速飛行設計中已經發展了參考溫度法來進行粘性力的計算，但仍然以自由來流條件作為粘性層外的無粘流動參數，沒有考慮復雜高超聲速繞流對粘性層外緣流動參數的影響，因此準確性不高。
【發明內容】
[0003]為了解決現有技術的問題，本發明實施例提供了一種用于快速計算飛行器高超聲速粘性力的方法和裝置。所述技術方案如下:
[0004]一方面，提供了 一種用于快速計算飛行器高超聲速粘性力的方法，所述方法包括:
[0005]步驟SI，對無粘Euler方程采用計算流體動力學CFD方法獲得所述飛行器高超聲速的表面無粘流場物理參數，其中，所述表面無粘流場物理參數包括:表面的壓力、表面的絕熱溫度和表面的速度；
[0006]步驟S2，獲取所述飛行器高超聲速的表面的所有表面單元的雷諾數；
[0007]步驟S3，根據所述表面無粘流場物理參數和所述表面的所有表面單元的雷諾數得到所述飛行器高超聲速的粘性力。
[0008]在如上所述的方法中，優選，，所述步驟S2具體包括:
[0009]步驟S20，選取所述飛行器高超聲速的表面單元，根據經過所述表面單元的表面流線得到所述表面單元與所述表面流線的出發點之間的流線長度；
[0010]步驟S21，根據所述表面單元的絕熱溫度、所述表面單元的速度和所述飛行器高超聲速的表面溫度得到粘性邊界層內的參考溫度；
[0011]步驟S22，根據所述粘性邊界層內的參考溫度得到所述粘性邊界層內的粘性系數；
[0012]步驟S23，根據所述表面單元的壓力和所述粘性邊界層內的參考溫度得到所述粘性邊界層內的密度；
[0013]步驟S24，根據所述粘性邊界層內的密度、粘性系數、所述表面單元的速度、所述流線長度得到所述表面單元的雷諾數；
[0014]循環執行步驟S21、步驟S22、步驟S23和步驟S24，得到所述飛行器高超聲速的表面的所有表面單元的雷諾數。
[0015]在如上所述的方法中，優選，所述步驟S20中的經過所述表面單元的表面流線和所述表面流線的出發點均可由所述表面單元在空間直角坐標系中的坐標和所述表面單元的速度根據流線定義采用四階顯示龍格庫塔方法求得。
[0016]在如上所述的方法中，優選，所述步驟S21具體包括:
[0017]根據
【權利要求】
1.一種用于快速計算飛行器高超聲速粘性力的方法，其特征在于，所述方法包括: 步驟SI，對無粘Euler方程采用計算流體動力學CFD方法獲得所述飛行器高超聲速的表面無粘流場物理參數，其中，所述表面無粘流場物理參數包括:表面的壓力、表面的絕熱溫度和表面的速度； 步驟S2，獲取所述飛行器高超聲速的表面的所有表面單元的雷諾數； 步驟S3，根據所述表面無粘流場物理參數和所述表面的所有表面單元的雷諾數得到所述飛行器高超聲速的粘性力。
2.根據權利要求1所述的方法，其特征在于，所述步驟S2具體包括: 步驟S20，選取所述飛行器高超聲速的表面單元，根據經過所述表面單元的表面流線得到所述表面單元與所述表面流線的出發點之間的流線長度； 步驟S21，根據所述表面單元的絕熱溫度、所述表面單元的速度和所述飛行器高超聲速的實際表面溫度得到粘性邊界層內的參考溫度； 步驟S22，根據所述粘性邊界層內的參考溫度得到所述粘性邊界層內的粘性系數； 步驟S23，根據所述表面單元的壓力和所述粘性邊界層內的參考溫度得到所述粘性邊界層內的密度； 步驟S24，根據所述粘性邊界層內的密度、粘性系數、所述表面單元的速度、所述流線長度得到所述表面單元的雷諾數； 循環執行步驟S21、步驟S22、步驟S23和步驟S24，得到所述飛行器高超聲速的表面的所有表面單元的雷諾數。
3.根據權利要求2所述的方法，其特征在于，所述步驟S20中的經過所述表面單元的表面流線和所述表面流線的出發點均可由所述表面單元在空間直角坐標系中的坐標和所述表面單元的速度根據流線定義采用四階顯示龍格庫塔方法求得。
4.根據權利要求2或3所述的方法，其特征在于， 所述步驟S21具體包括: 根據
5.根據權利要求4所述的方法，其特征在于，所述步驟S3具體包括: 步驟S31，根據所述表面單元的雷諾數得到所述表面單元的粘性摩擦力系數； 步驟S32，根據所述表面單元的粘性摩擦力系數、所述表面單元的速度、所述粘性邊界層內的密度和所述表面單元的面積得到所述表面單元的粘性力； 循環執行步驟S31和步驟S32，并將飛行器表面的所有表面單元的粘性力進行求和，得到所述飛行器高超聲速的粘性力。
6.根據權利要求5所述的方法，其特征在于， 所述步驟S31具體包括: 如果所述表面單元的雷諾數小于等于轉捩雷諾數，則所述表面單元的粘性摩擦力系數
7.一種用于快速計算飛行器高超聲速粘性力的裝置，其特征在于，所述裝置包括: 獲得裝置，用于對無粘Euler方程采用計算流體動力學CFD方法獲得所述飛行器高超聲速的表面無粘流場物理參數，其中，所述表面無粘流場物理參數包括:表面的壓力、表面的絕熱溫度和表面的速度； 獲取裝置，用于獲取所述飛行器高超聲速的表面的所有表面單元的雷諾數； 得到裝置，用于根據所述獲得裝置獲得的表面無粘流場物理參數和所述獲取裝置獲取的所述表面的所有表面單元的雷諾數得到所述飛行器高超聲速的粘性力。
8.根據權利要求7所述的裝置，其特征在于，所述獲取裝置包括: 第一獲取單元，用于選取所述飛行器高超聲速的表面單元，根據經過所述表面單元的表面流線得到所述表面單元與所述表面流線的出發點之間的流線長度； 第二獲取單元，用于根據所述表面單元的絕熱溫度、所述表面單元的速度和所述飛行器高超聲速的實際表面溫度得到粘性邊界層內的參考溫度； 第三獲取單元，用于根據所述粘性邊界層內的參考溫度得到所述粘性邊界層內的粘性系數； 第四獲取單元，用于根據所述表面單元的壓力和所述第二獲取單元獲取的參考溫度得到所述粘性邊界層內的密度； 第五獲取單元，用于根據所述第四獲取單元獲取的密度、所述第三獲取單元獲取的粘性系數、所述表面單元的速度、所述第一獲取單元獲取的流線長度得到所述表面單元的雷諾數； 獲取循環單元，用于使所述第一獲取單元、所述第二獲取單元、所述第三獲取單元、所述第四獲取單元和所述第五獲取單元循環運行，得到所述飛行器高超聲速的表面的所有表面單元的雷諾數。
9.根據權利要求8所述的裝置，其特征在于， 所述第二獲取單元具體用于根據
10.根據權利要求9所述的裝置，其特征在于，所述得到裝置包括: 第一得到單元，用于根據所述獲取裝置獲取的所述表面單元的雷諾數得到所述表面單元的粘性摩擦力系數 第二得到單元，用于根據所述第一得到單元得到的粘性摩擦力系數、所述表面單元的速度、所述粘性邊界層內的密度和所述表面單元的面積得到所述表面單元的粘性力； 得到循環單元，用于使所述第一得到單元和所述第二得到單元循環運行，并將飛行器表面的所有表面單元的粘性力進行求和，得到所述飛行器高超聲速的粘性力。
【文檔編號】G06F17/50GK103617338SQ201310718622
【公開日】2014年3月5日 申請日期:2013年12月24日 優先權日:2013年12月24日
【發明者】龔安龍, 紀楚群, 劉周 申請人:中國航天空氣動力技術研究院