三維雙楔面定常超聲速流動(dòng)研究

楊 旸，滕宏輝，姜宗林

（中國科學(xué)院力學(xué)研究所高溫氣體動(dòng)力學(xué)重點(diǎn)實(shí)驗(yàn)室，北京 100190）

0 引言

三維雙楔面超聲速流動(dòng)是一個(gè)具有廣泛應(yīng)用背景的研究問題。對(duì)于超聲速／高超聲速飛行器，其翼身結(jié)合、垂尾／平尾夾角等處會(huì)形成復(fù)雜的三維激波／激波相互作用和激波／邊界層相互作用，影響飛行器的氣動(dòng)性能。又如對(duì)于超燃沖壓發(fā)動(dòng)機(jī)三面壓縮進(jìn)氣道，雙楔面超聲速流動(dòng)形成的三維激波結(jié)構(gòu)會(huì)導(dǎo)致低總壓區(qū)的形成，降低進(jìn)氣道的總體性能。過去的幾十年里，人們對(duì)于該問題給予了足夠的重視并開展了較多的研究［1－7］。Charwat等［1］對(duì)該問題做出了較早的研究，并依據(jù)實(shí)驗(yàn)結(jié)果闡述了三維激波相互作用后形成的橋激波、透射激波等基本波系結(jié)構(gòu)。Watson等［2］通過實(shí)驗(yàn)發(fā)現(xiàn)透射激波與邊界層相互干擾形成的三維λ波及三維渦結(jié)構(gòu)，并依此解釋了壁面熱流和壓力的變化。Marsilio［5］對(duì)該問題中渦的不穩(wěn)定性做出了數(shù)值研究，結(jié)果顯示隨著來流馬赫數(shù)的增加，截面上的滑移線逐漸卷起形成對(duì)稱結(jié)構(gòu)的渦，并隨著馬赫數(shù)的進(jìn)一步增加而失穩(wěn)并不再對(duì)稱。岳連捷等［6］從進(jìn)氣道角度研究該問題，發(fā)現(xiàn)激波相互作用形成的橋波區(qū)總壓恢復(fù)系數(shù)較低，并研究了楔面角度變化對(duì)橋波區(qū)的影響。該雙楔面超聲速流動(dòng)問題主要包括遠(yuǎn)離邊界層近似無粘的激波／激波相互作用和粘性主導(dǎo)的激波／邊界層相互作用兩個(gè)方面，在具體研究過程中可以將其分別討論。圖1是該三維無粘流動(dòng)的示意圖，其中來流馬赫數(shù)為M0，兩個(gè)楔面的傾角和激波角分別為θ和β，前緣后掠角及雙楔面間夾角分別為χ和ν。θn和βn分別為垂直于楔面前緣方向上的楔面角度和斜激波角，當(dāng)后掠角χ趨于0時(shí)θn和βn分別趨于θ和β。超聲速流動(dòng)使得兩個(gè)楔面上分別形成斜激波面i，其相互作用后在對(duì)稱面反射形成規(guī)則或馬赫反射。對(duì)于馬赫反射，在兩楔面夾角處形成橋激波m。圖2為Goonko等［7］的無粘數(shù)值模擬結(jié)果中垂直于來流方向的二維截面上的激波結(jié)構(gòu)?？梢钥吹?，截面中的定常激波在對(duì)稱面上反射形成了過渡馬赫反射和雙馬赫反射結(jié)構(gòu)，然而這兩種反射結(jié)構(gòu)在二維定常激波反射中是無法形成的［8］。因此，三維雙楔面超聲速流動(dòng)形成的激波／激波相互作用在二維截面上的定常激波反射結(jié)構(gòu)無法用二維定常激波反射理論解釋。之前人們的研究大多是通過實(shí)驗(yàn)和數(shù)值模擬得出結(jié)論，并沒有從理論上分析和解釋三維空間內(nèi)形成這種激波結(jié)構(gòu)的原因，也沒有給出形成各種不同馬赫反射結(jié)構(gòu)的來流馬赫數(shù)和楔面傾角組合范圍。本文著眼于遠(yuǎn)離邊界層的三維激波相互作用，通過理論分析將該三維定常激波相互作用與二維非定常激波反射建立聯(lián)系，解釋了二維截面上形成過渡馬赫反射和雙馬赫反射的原因。同時(shí)，通過二維簡化，利用非定常激波楔面反射理論求解得出截面上形成各種反射結(jié)構(gòu)的馬赫數(shù)和楔面傾角組合范圍以及雙楔面間夾角、前緣后掠角的影響，并通過數(shù)值模擬加以驗(yàn)證。

圖1 三維雙楔面超聲速流動(dòng)示意圖Fig.1 Schematic of supersonic flow over 3－D intersecting wedges

圖2 二維橫截面上過渡馬赫反射和雙馬赫反射Fig.2 The transitional－Mach reflection and double－Mach reflection on 2－D cross section

1 理論分析方法

由于本文所探討的雙楔面是對(duì)稱結(jié)構(gòu)，因此在兩楔面之間三維空間中存在有對(duì)稱面，理論分析時(shí)只需考慮對(duì)稱面的一側(cè)即可。由于目前為止還沒有三維激波相互作用的理論，因此該三維問題需要二維簡化后才能進(jìn)行理論分析。將三維激波簡化到二維截面進(jìn)行分析的前提是截面內(nèi)的激波只具備二維特性，即在截面內(nèi)二維激波前后垂直于截面的氣流速度分量不因激波的存在而改變，二維激波的存在只改變了截面內(nèi)平行于截面的氣流速度分量?；谝陨峡紤]，選取垂直于斜激波面i與對(duì)稱面的交線AB的截面作為二維特征截面簡化該三維問題，如圖3 所示。圖4是該特征截面的示意圖，截面內(nèi)激波i′前后的氣流速度分別被分解為垂直于截面的速度分量和平行于截面的速度分量。若能證明激波前后垂直于截面的速度分量相等，即M0n＝M1n，則此二維簡化是可行的。以下對(duì)此進(jìn)行證明：

對(duì)于來流馬赫數(shù)M0，前緣后掠角和雙楔面間夾角分別為χ和ν的流動(dòng)，來流經(jīng)過斜激波面i后的速度矢量為：

其中：

AB方向的矢量為：

圖3 理論分析方法示意圖Fig.3 The schematic of theoretical analysis

圖4 二維特征截面Fig.4 The 2－D characteristic cross section

又有：

因此，可以求得M1在AB方向的投影為：

同時(shí)，不難得出M0n＝M0cosγ，由此證明了二維特征截面上激波兩側(cè)垂直于截面的氣流速度分量相等，即將三維激波結(jié)構(gòu)簡化到該二維特征截面進(jìn)行分析是正確的。

在二維特征截面內(nèi)，激波i′前平行于截面的氣流速度分量為：

該速度分量方向平行于二維特征截面與對(duì)稱面的交線s。DF為二維特征截面內(nèi)垂直于激波i′的直線，通過幾何關(guān)系不難證明DF 垂直于斜激波面i。由點(diǎn)A、B、C的坐標(biāo)可得斜激波面i及直線AB的解析式，即可通過點(diǎn)到平面和直線的距離公式求出EF 和DE，最終得出二維特征截面內(nèi)激波i與對(duì)稱面s的夾角：

在二維特征截面內(nèi)，定常激波i′靜止，激波前的空氣以馬赫數(shù)M0τ沿二維特征截面與對(duì)稱面的交線s朝著激波方向運(yùn)動(dòng)。通過伽利略變換，截面內(nèi)激波i′前方的空氣變換為靜止?fàn)顟B(tài)，激波i′以馬赫數(shù)M0τ沿著二維截面與對(duì)稱面的交線s方向運(yùn)動(dòng)并形成反射。至此，通過伽利略變換，二維特征截面內(nèi)定常激波反射問題轉(zhuǎn)變?yōu)榱朔嵌ǔ＜げǖ男娣瓷鋯栴}，如圖5所示。由于在一定的激波馬赫數(shù)和楔面傾角組合下非定常激波楔面反射可以形成過渡馬赫反射和雙馬赫反射，因此二維截面內(nèi)形成以上兩種馬赫反射結(jié)構(gòu)的現(xiàn)象也得到了解釋。同時(shí)，激波馬赫數(shù)M0τ以及激波與對(duì)稱面夾角 都已求得，因此可以通過二維非定常激波楔面反射的結(jié)論求解該截面內(nèi)激波反射構(gòu)造，并得出三維空間內(nèi)激波相互作用結(jié)構(gòu)。

圖5 伽利略變換后運(yùn)動(dòng)激波反射示意圖Fig.5 Schematic of the pseudo－steady shock wave reflection after the Galilean transformation

2 數(shù)值模擬方法

本文研究三維空間內(nèi)激波相互作用結(jié)構(gòu)，粘性作用可以忽略，因此采用三維Euler方程作為控制方程。對(duì)于理想氣體，直角坐標(biāo)系下的控制方程可以表達(dá)為：

上述表達(dá)式中，ρ、p和E分別代表氣體密度、壓力和單位質(zhì)量總內(nèi)能，u、v和w為氣體在x、y和z方向上的速度分量。單位質(zhì)量總內(nèi)能表達(dá)式為：

控制方程的離散為如下方程：

其中的對(duì)流項(xiàng)通量項(xiàng)F、G、H利用二階NND 格式［9］求解，采用minmod函數(shù)作為限制器。計(jì)算域網(wǎng)格采用正交化處理后的均勻結(jié)構(gòu)網(wǎng)格。由于采用并行程序計(jì)算，網(wǎng)格數(shù)達(dá)千萬量級(jí)。計(jì)算域入口給定來流條件；出口邊界假定沿來流方向流場(chǎng)參數(shù)變化很小，邊界條件由外推法給出；上邊界和右邊界給定無反射邊界條件；固壁邊界采用無穿透邊界條件。

3 結(jié)果與討論

3.1 無后掠角垂直雙楔面

對(duì)于兩個(gè)楔面間夾角ν＝90°且后掠角χ＝0°的情況，在不同馬赫數(shù)和楔面傾角的組合下，入射斜激波在對(duì)稱面上反射形成各種不同的反射結(jié)構(gòu)。圖6為本文通過數(shù)值模擬得出的在垂直于自由來流速度方向的截面上形成的規(guī)則反射（RR）、單馬赫反射（SMR）、過渡馬赫反射（TMR）以及雙馬赫反射（DMR）結(jié)構(gòu)。由于二維非定常激波楔面反射理論中弱激波反射理論尚不完善，本文中暫不考慮弱激波反射問題。對(duì)于無后掠角垂直雙楔面，M0τ和η的表達(dá)式可簡化為：

圖6 截面上激波反射結(jié)構(gòu)Fig.6 The shock reflection patterns on the cross section

這樣，對(duì)于不同的來流馬赫數(shù)以及楔面傾角所對(duì)應(yīng)的激波角，M0τ和η就可以求得。通過伽利略變換，可以在二維特征截面內(nèi)將定常的激波反射問題轉(zhuǎn)換為以M0τsinη的馬赫數(shù)運(yùn)動(dòng)的激波在傾角為90°－η的楔面上的反射，如圖5所示。由此通過二維非定常激波反射已有的結(jié)論就可以求得二維特征截面內(nèi)激波反射構(gòu)造。圖7為理論分析及數(shù)值驗(yàn)證結(jié)果，圖中橫坐標(biāo)為自由來流馬赫數(shù)，縱坐標(biāo)為楔面傾角，曲線代表通過上述理論分析得出的各種激波反射結(jié)構(gòu)的分界線，各圖標(biāo)分別代表在各馬赫數(shù)和楔面傾角組合下的數(shù)值模擬結(jié)果，其中實(shí)心圓點(diǎn)表示在給定馬赫數(shù)下由于楔面傾角過大而形成了脫體激波?？梢钥吹?，理論分析結(jié)果與數(shù)值結(jié)果吻合良好，由此可以證明通過二維簡化分析該三維問題是正確的，以及通過二維非定常激波反射解釋該三維激波構(gòu)造是可行的。

圖7 不同楔面傾角和馬赫數(shù)組合下截面上的激波反射結(jié)構(gòu)（υ＝90°，χ＝0°）Fig.7 The shock reflection patterns on the cross section for different combinations of wedge angles and Mach numbers（υ＝90°，χ＝0°）

3.2 雙楔面間夾角及前緣后掠角的影響

當(dāng)兩楔面間夾角ν≠90°或前緣后掠角χ≠0°時(shí)，二維截面上各種激波反射結(jié)構(gòu)的形成范圍與無后掠角的垂直雙楔面的結(jié)果有所不同。以下分別討論兩楔面間夾角以及前緣后掠角的影響。對(duì)于ν≠90°，χ＝0°的情況，M0τ和η的表達(dá)式可簡化為：

圖8為通過上述理論分析所得出的兩楔面間夾角對(duì)截面上激波構(gòu)造的影響及數(shù)值模擬對(duì)此結(jié)果的驗(yàn)證。可以看到，隨著兩楔面間夾角的變化，二維截面上各種激波反射結(jié)構(gòu)出現(xiàn)的范圍也有所變化，即隨著兩楔面間夾角的減小，單馬赫－過渡馬赫反射以及過渡馬赫－雙馬赫反射的分界線朝著來流小馬赫數(shù)方向移動(dòng)，而馬赫－規(guī)則反射分界線朝大楔面傾角方向移動(dòng)。注意到在圖8（d）中由理論分析得出的規(guī)則反射區(qū)中出現(xiàn)了馬赫反射的數(shù)值結(jié)果。值得注意的是，圖8中的規(guī)則－馬赫反射的實(shí)線分界線均是由聲速準(zhǔn)則（θS）得出的。von Neumann準(zhǔn)則要求來流馬赫數(shù)大于Mc0＝2.202，當(dāng)兩楔面夾角為120°和105°時(shí)各種來流馬赫數(shù)和楔面傾角組合下二維截面內(nèi)來流馬赫數(shù)均小于Mc0，而當(dāng)ν＝60°及75°時(shí)部分組合下截面內(nèi)來流馬赫數(shù)達(dá)到Mc0，可以得出該情況下的von Neumann準(zhǔn)則（θN），即圖8（c）、（d）中的虛線。通過數(shù)值模擬可以得出，在該虛線上下，激波結(jié)構(gòu)分別呈馬赫反射和規(guī)則反射結(jié)構(gòu)，即此時(shí)二維定常激波反射所遵循的von Neumann準(zhǔn)則代替聲速準(zhǔn)則成為規(guī)則－馬赫反射轉(zhuǎn)變的分界線。

圖8 不同楔面傾角和馬赫數(shù)組合下截面上的激波反射結(jié)構(gòu)（υ≠90°，χ＝0°）Fig.8 The shock reflection patterns on the cross sectionfor different combinations of wedge angles and Mach numbers（υ≠90°，χ＝0°）

對(duì)于ν＝90°，χ≠0°的情況，M0τ和η的表達(dá)式可簡化為：

圖9是通過上述理論分析所得出的前緣后掠角對(duì)截面上激波構(gòu)造的影響以及數(shù)值模擬對(duì)此結(jié)果的驗(yàn)證。隨著后掠角增大，單馬赫－過渡馬赫反射以及過渡馬赫－雙馬赫反射分界線向來流高馬赫數(shù)方向移動(dòng)，而規(guī)則－馬赫反射分界線向小楔面傾角方向移動(dòng)。注意到圖9的幾種情況下二維截面內(nèi)的激波馬赫數(shù)均達(dá)不到Mc0，即不存在由von Neumann準(zhǔn)則得出的規(guī)則－馬赫反射分界線。

圖9 不同楔面傾角和馬赫數(shù)組合下截面上的激波反射結(jié)構(gòu)（υ＝90°，χ≠0°）Fig.9 The shock reflection patterns on the cross section for different combinations of wedge angles and Mach numbers（υ＝90°，χ≠0°）

4 結(jié)論

本文采用理論分析和數(shù)值模擬相結(jié)合的方法，對(duì)三維雙楔面超聲速流動(dòng)無粘激波結(jié)構(gòu)做出了研究，得出以下結(jié)論：

（1）三維雙楔面超聲速流動(dòng)形成的定常激波結(jié)構(gòu)可以通過二維簡化利用二維非定常激波反射理論分析得出。理論分析結(jié)果指出，入射斜激波會(huì)在對(duì)稱面上反射形成規(guī)則反射、單馬赫反射、過渡馬赫反射以及雙馬赫反射等多種激波反射結(jié)構(gòu)，并通過數(shù)值模擬對(duì)此結(jié)果進(jìn)行了驗(yàn)證；

（2）通過理論分析探討了兩楔面間夾角和前緣后掠角對(duì)該激波反射結(jié)構(gòu)的影響，得出了在不同雙楔面夾角和后掠角下各種激波反射構(gòu)造的形成范圍，并分析了各種反射構(gòu)造形成范圍的變化趨勢(shì)；

（3）在二維截面上兼有二維非定常激波和定常激波反射的特性，即形成了類似于二維非定常激波楔面反射的多種馬赫反射結(jié)構(gòu)，同時(shí)其規(guī)則－馬赫反射轉(zhuǎn)變卻遵循二維定常激波反射所遵循的von Neumann準(zhǔn)則。

［1］CHARWAT A，REDEKEOPP L.Supersonic interference flow along the corner of intersecting wedges［J］.AIAAJournal，1967，5（3）：480－488.

［2］WATSON D，WEINSTEIN M.A study of hypersonic corner flow interaction［J］.AIAAJournal，1971，9（7）：1280－1286.

［3］KUTLER P.Supersonic flow in the corner formed by two intersecting wedges［J］.AIAAJournal，1974，12（5）：577－579.

［4］MARCONI F.Supersonic，inviscid，conical corner flowfields［J］.AIAAJournal，1980，18（1）：78－84.

［5］MARSILIO R.Vortical solutions in supersonic corner flows［J］.AIAAJournal，1993，31（9）：1651－1658.

［6］岳連捷，徐顯坤，張新宇，等.三面壓縮進(jìn)氣道三維激波相干研究［A］／／中國力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì)文集［C］.鄭州，2009.

［7］GOONKO Y，KUDRYAVTSEV A，CHPOUN A.3D interaction of shock waves in corner flow ［A］／／Shock Waves，Proceeding of the 24th international symposium on shock waves［C］.Beijing，2004.

［8］BEN－DOR G.Shock wave reflection phenomena［M］.New York：Springer－Verlag Inc.，1992.

［9］張涵信.無波動(dòng)、無自由參數(shù)的耗散差分格式［J］.空氣動(dòng)力學(xué)報(bào)，1988，6（2）：143－165.