一種彎曲激波壓縮流場氣流參數(shù)的快速確定方法
【技術(shù)領(lǐng)域】
[0001 ]本發(fā)明涉及一種超聲速或高超聲速彎曲激波壓縮流場氣流參數(shù)的快速確定方法， 屬于流場氣流參數(shù)的確定方法。
【背景技術(shù)】
[0002] 彎曲激波壓縮流場中包含前緣激波和等熵壓縮波，兩者相互作用還會形成反射波 系、滑流間斷，反射波系和滑流間斷又將繼續(xù)與下游的壓縮波產(chǎn)生相互作用，對壁面附近流 動產(chǎn)生影響，因此壓縮面發(fā)出的左行特征線一般不是直線，特征線上的馬赫數(shù)、壓力等參數(shù) 并不均勻，給流場分析帶來了困難，需要借助特征線法、有限體積法等流場計算方法進行計 算，其復雜程度較高、計算耗時較多，不利于流場分析和初步設(shè)計中的快速估算。
[0003] 為了進行快速計算，"Numer i ca 1 me thod and r e su 11 s for inv i sc i d supersonic flow over a compressive ramp"（Emanuel G)、"彎曲激波壓縮面設(shè)計及試驗 研究"（居燕）等文獻中仍假設(shè)彎曲壓縮面附近為PrandtΙ-Meyer流動，直接根據(jù)Prandt 1-Meyer方程計算壓縮面附近參數(shù)，在激波附近則按照激波與離散壓縮波的相交來計算激波 后參數(shù)。而事實上彎曲激波上兩者相互作用產(chǎn)生反射波系，下游流場不再是簡單的 Prandtl-Meyer流動，忽略反射波系的計算將有可能產(chǎn)生較大的誤差。文獻"超聲/高超聲速 非均勻來流下曲面壓縮系統(tǒng)研究"（潘瑾）中嘗試了基于設(shè)計樣本的計算數(shù)據(jù)進行擬合的方 法，建立了彎曲激波坐標的估算公式。但是其結(jié)果依賴于所采用的型面設(shè)計方法，并且建立 過程需要大量的計算樣本，因此應(yīng)用不易擴展，而且擬合過程也會帶來誤差。

【發(fā)明內(nèi)容】

[0004] 發(fā)明目的：為了克服現(xiàn)有技術(shù)中存在的不足，本發(fā)明提供一種超聲速或高超聲速 彎曲激波壓縮流場氣流參數(shù)的快速確定方法，可以近似計算彎曲激波壓縮流場中壁面參 數(shù)、激波形狀、波后參數(shù)、流場內(nèi)流線形狀、流線上參數(shù)以及出口參數(shù)，計算方法簡單、快速。
[0005] 技術(shù)方案:為實現(xiàn)上述目的，本發(fā)明采用的技術(shù)方案為：
[0006] -種彎曲激波壓縮流場氣流參數(shù)的快速確定方法，分別獨立計算壓縮壁面發(fā)出的 壓縮波系與彎曲激波反射的膨脹波系造成的流動參數(shù)的變化，在計算過程中，將原流場中 的壓縮波系簡化為匯聚于點S的Prandtl-Meyer壓縮波，將原流場中的膨脹波系簡化為經(jīng)過 點S的一道膨脹波，根據(jù)簡化后的Prandtl-Meyer壓縮波和膨脹波快速確定流場中的壓縮壁 面參數(shù)、彎曲激波形狀、波后參數(shù)、流場內(nèi)流線形狀、流線上參數(shù)以及出口參數(shù)。
[0007] 設(shè)自由來流的馬赫數(shù)、壓力和流動方向角分別為M〇〇、p〇〇和θ〇〇,壓縮壁面的前緣點與 彎曲激波的前緣點重合，記為點〇st art(〇，〇)，壓縮壁面的末端點記為Wend，彎曲激波的末端點 記為Send;壓縮壁面的形狀為y w=f (Xw)，0 < Xw < L;根據(jù)壓縮壁面形狀計算點0startfe置的流 動方向角9start，然后根據(jù)氣體動力學中的斜激波關(guān)系計算點0 start位置的馬赫數(shù)Mstart和壓 力P start ；
[0008] 對原流場中其他參數(shù)采用如下方法確定：
[0009] (1)在原流場中，經(jīng)過壓縮壁面上點Wdxw^ym)的壓縮波與彎曲激波的交點記為 Sj;從原流場中分離出點Wi之前的壓縮波系并簡化為匯聚于點Sj的Prandtl-Meyer壓縮波， 根據(jù)點Wi與點0 start的位置關(guān)系計算Wi位置的流動方向角0wi，然后采用Prandtl-Meyer公式 計算點Wi位置的馬赫數(shù)M wi和壓力pwi;同時計算自由來流經(jīng)過轉(zhuǎn)折角等于點Wi位置流動方向 角0wi時產(chǎn)生的斜激波之后的馬赫數(shù)、壓力和流動方向角沒wi;
[0010] 對壓縮壁面上從點〇start開始的每一點Wi重復該過程，直至點Wi到達點W end，就確定 了原流場中的壓縮波造成的壓縮壁面的參數(shù)變化；
[0011] (2)從原流場中分離出點&之前的膨脹波系并簡化為經(jīng)過點&的一道膨脹波，該膨 脹波與壓縮壁面的交點記為該膨脹波導致的轉(zhuǎn)折角心為：
[0013]式中:k為氣體比熱容;根據(jù)下式計算點Wj的坐標(Xwj，ywj):
[0018] 對壓縮壁面上從點0start開始的每一點Wj重復該過程，直至點Wj到達點W end，就確定 了原流場中的膨脹波造成的壓縮壁面的參數(shù)變化；
[0019] (3)計算點I位置經(jīng)壓縮波和膨脹波合成后流動相對于點0start的轉(zhuǎn)折角 correction，
[0020] 5wj-correction- ( ^wj-9start) -2δ j
[0021 ] 采用Prandtl-Meyer公式，根據(jù)Ostart位置的馬赫數(shù)Mstart和壓力Pstart計算點Wj位置 經(jīng)壓縮波和膨脹波合成后的馬赫數(shù)Mwj-cOT：rec;tii3r^PStpwj-c;c)rrec;tic)n;
[0022] 對壓縮壁面上從點0start開始的每一點Wj重復該過程，直至點Wj到達點W end，就確定 了原流場中的壓縮波和膨脹波合成后造成的壓縮壁面的參數(shù)變化；
[0023] (4)計算點Sj位置相對于自由來流的轉(zhuǎn)折角3^-。。1^心。11:
[0024] 3sj-correcti〇n- 0ffi+6 j-θ〇〇
[0025] 根據(jù)轉(zhuǎn)折角Ssj-ccirrectiQn和氣體動力學理論計算點Sj位置的馬赫數(shù)M sj-ccirrecticin、S 力Psj-correction和流動萬向角 〇sj-correction ;
[0026]根據(jù)下式計算點Sj的坐標(Xsj，ysj):
[0031] 對壓縮壁面上從點0start開始的每一點Wj重復該過程，直至點Wj到達點W end，就確定 了彎曲激波的形狀和波后氣動參數(shù)；
[0032] (5)在原流場中，經(jīng)過點Wk的壓縮波與以點Sj為起點的流線的交點為R，點R的流動 力-(? @ Θγ-correction 以及點R相對于點Sj的轉(zhuǎn)折角S?rrec;ticiA:
[0033] 0r-correction-^wk
[0034] 5r-correcti〇n= (9wk-25j)-( 0wi-25i)
[0035] 采用Prandtl-Meyer公式計算點1?位置的馬赫數(shù)Mr-correction和壓力Pr-correction;根據(jù) 下式計算點R的坐標(Xr，y r):
[0040] 對壓縮壁面上從點I開始的每一點％重復該過程，直至點％到達點Wend，就確定了 以點Sj為起點的流線的形狀和流線上氣動參數(shù)；
[0041] (6)出口截面的馬赫數(shù)、壓力和總壓取為壓縮壁面、彎曲激波和流線末端參數(shù)的平 均值：
[0043] 式中：η為總的流線數(shù)量，Mw-end、pw-end和為點W end位置經(jīng)壓縮波和膨脹波合成 后的馬赫數(shù)、壓力和總壓，M?d、p?4P^nd為第r條流線末端位置的馬赫數(shù)、壓力和總壓， Msid、Ps-end和為點S-位置的馬赫數(shù)、壓力和總壓。
[0044] 有益效果:本發(fā)明提供的彎曲激波壓縮流場氣流參數(shù)的快速確定方法，計算過程 較為簡單，計算結(jié)果可達到較高的精度，能夠用于對流場的計算和分析。
【附圖說明】
[0045] 圖1為原流場示意圖，包括壓縮壁面、彎曲激波、壓縮壁面發(fā)出的壓縮波（即左行特 征線)和彎曲激波上反射的膨脹波(即右行特征線）；
[0046]圖2為簡化后的Prandtl-Meyer壓縮波和膨脹波示意圖；
[0047]圖3為彎曲激波上的點、壓縮壁面上的點和流線上的點的計算示意圖；
[0048]圖4為出口截面的計算示意圖。
【具體實施方式】
[0049]下面結(jié)合附圖對本發(fā)明作更進一步的說明。
[0050] 如圖1所示為一種圖1是原流場示意圖，其中不僅包含前緣激波和壓縮波，兩者相 互作用還會形成反射波系、滑流間斷，反射波系和滑流間斷又將繼續(xù)與下游的壓縮波產(chǎn)生 相互作用，對壁面附近流動產(chǎn)生影響，導致壓縮面發(fā)出的左行特征線不再是直線，特征線上 的馬赫數(shù)、壓力等參數(shù)不均勻，這是精確的計算較為繁瑣和困難的原因。本發(fā)明提出的方法 對計算過程進行如下簡化:分別獨立計算壓縮壁面發(fā)出的壓縮波系與彎曲激波反射的膨脹 波系造成的流動參數(shù)的變化，在計算過程中，將原流場中的壓縮波系簡化為匯聚于點S的 Prandtl-Meyer壓縮波，將原流場中的膨脹波系簡化為經(jīng)過點S的一道膨脹波，根據(jù)簡化后 的Prandtl-Meyer壓縮波和膨脹波快速確定流場中的壓縮壁面參數(shù)、彎曲激波形狀、波后參 數(shù)、流場內(nèi)流線形狀、流線上參數(shù)以及出口參數(shù)。
[0051] 設(shè)自由來流的馬赫數(shù)、壓力和流動方向角分別為M〇〇、p〇〇和θ〇〇,壓縮壁面的前緣點與 彎曲激波的前緣點重合，記為點〇st art(0，0)，壓縮壁面的末端點記為Wend，彎曲激波的末端點 記為Send;壓縮壁面的形狀為y w=f (Xw)，0 < Xw < L;根據(jù)壓縮壁面形狀計算點0start位置的流 動方向角9start，然后根據(jù)氣體動力學中的斜激波關(guān)系計算點0 start位置的馬赫數(shù)Mstart和壓 文JPstart;
[0052] 對原流場中其他參數(shù)采用如下方法確定：
[0053] (1)在原流場中，經(jīng)過壓縮壁面上點的壓縮波與彎曲激波的交點記為 Sj;從原流場中分離出點Wi之前的壓縮波系并簡化為匯聚于點Sj的Prandtl-Meyer壓縮波， 根據(jù)點Wi與點Ostart的位置關(guān)系計算Wi位置的流動方向角0wi，然后采用Prandtl-Meyer公式 計算點Wi位置的馬赫數(shù)M wi和壓力pwi;同時計算自由來流經(jīng)過轉(zhuǎn)折角等于點Wi位置流動方向 角9wi時產(chǎn)生的斜激波之后的馬赫數(shù)、壓力和流動方向角;
[0054] 對壓縮壁面上從點0start開始的每一點Wi重復該過程，直至點Wi到達點W end，就確定 了原流場中的壓縮波造成的壓縮壁面的參數(shù)變化；
[0055] (2)從原流場中分離出點&之前的膨脹波系并簡化為經(jīng)過點&的一道膨脹波，該膨 脹波與壓縮壁面的交點記為該膨脹波導致的轉(zhuǎn)折角心為：
[0057]式中:k為氣體比熱容;根據(jù)下式計算點Wj的坐標(xwj，ywj):
[0062] 對壓縮壁面上從點0start開始的每一點Wj重