張恩來, 李祝飛, 李一鳴, 楊基明

(中國科學(xué)技術(shù)大學(xué) 近代力學(xué)系, 合肥 230027)

0 引 言

激波干擾現(xiàn)象，因其產(chǎn)生的激波、剪切層、膨脹波、超聲速射流以及伴隨的旋渦和湍流等復(fù)雜流動結(jié)構(gòu)，不僅是高超聲速飛行器研究與發(fā)展過程中需要重點關(guān)注的空氣動力學(xué)問題之一，也是困擾工程應(yīng)用的一大難題[1,2]。以斜激波與圓柱/圓球弓形激波之間的相互作用為典型模型，根據(jù)激波干擾位置和激波強度不同，Edney[3]總結(jié)出6種具有普適性的激波干擾類型(Ⅰ～Ⅵ類)。其中，第Ⅳ類激波干擾因其產(chǎn)生的超聲速射流沖擊壁面會帶來嚴(yán)峻的熱/力載荷而受到設(shè)計和研究人員的廣泛關(guān)注[4-7]。在某些情況下，激波干擾產(chǎn)生的超聲速射流高度彎曲，不會直接沖擊壁面，這種激波干擾被稱為第Ⅳa類[8]。

高超聲速進(jìn)氣道流動存在多種形式的激波干擾現(xiàn)象，尤其以設(shè)計工況下前體激波入射進(jìn)氣道唇口[9]形成的激波干擾最為突出。針對前體激波與二元進(jìn)氣道唇口弓形激波的干擾問題，Wieting等[10]進(jìn)行了系統(tǒng)性研究，獲得了豐富的壓力和熱流數(shù)據(jù)。近年來，前緣高度后掠的內(nèi)轉(zhuǎn)式進(jìn)氣道[11-13]逐漸成為研究熱點之一。然而，其普遍采用的V形溢流口(唇口)與二元進(jìn)氣道唇口具有明顯的區(qū)別。最近，Xiao等[14]在不考慮前體激波的情況下，將內(nèi)轉(zhuǎn)式進(jìn)氣道溢流口簡化為V形鈍前緣，研究發(fā)現(xiàn)：在一定的幾何參數(shù)下，V形鈍前緣自身就會形成類似于Edney第Ⅳ類激波干擾的流動結(jié)構(gòu)。顯然，當(dāng)前體激波入射V形鈍前緣溢流口時，將產(chǎn)生更加復(fù)雜的波系結(jié)構(gòu)，很可能加劇溢流口的氣動熱。

另一方面，根據(jù)來流條件的不同，進(jìn)氣道會出現(xiàn)亞額定、額定和超額定等3種工作狀態(tài)[15]。針對二元進(jìn)氣道的研究表明[16-17]：在超額定狀態(tài)下，前體激波入射在唇口內(nèi)部，當(dāng)前體激波與唇口激波發(fā)生Edney第Ⅱ類激波干擾時，形成的亞聲速區(qū)可能導(dǎo)致進(jìn)氣道局部不起動。相對而言，內(nèi)轉(zhuǎn)式進(jìn)氣道唇口更加復(fù)雜，前體激波與其相互作用的研究尚未見系統(tǒng)報道，亟需針對性的機理揭示，以補充認(rèn)識上的不足。

采用完整的內(nèi)轉(zhuǎn)式進(jìn)氣道對V形溢流口激波干擾特性進(jìn)行研究，不僅會耗費大量的計算和實驗資源，更會帶來多種復(fù)雜因素的影響。抓住內(nèi)轉(zhuǎn)式進(jìn)氣道V形溢流口的主要流動特征，有側(cè)重地進(jìn)行簡化，繼而揭示其流動機理和關(guān)鍵影響因素，當(dāng)屬高效簡潔的研究方式。因此，本文將內(nèi)轉(zhuǎn)式進(jìn)氣道V形溢流口簡化為V形鈍前緣平板。鑒于內(nèi)轉(zhuǎn)式進(jìn)氣道的前體激波為復(fù)雜的三維曲面激波，為便于準(zhǔn)確描述和控制，本文暫不考慮入射激波的展向曲率，而是采用簡化的斜激波，以突出V形鈍前緣溢流口自身帶來的主要特征。繼而采用激波風(fēng)洞實驗和數(shù)值模擬相結(jié)合的方法，揭示入射激波與V形鈍前緣溢流口相對位置變化引起的激波干擾形態(tài)演變，以期為內(nèi)轉(zhuǎn)式進(jìn)氣道溢流口設(shè)計提供有價值的參考。

1 實驗和計算方法

1.1 實驗設(shè)備和模型

V形鈍前緣平板由較長的平直段和交叉位置倒圓部分組成，如圖1所示。交叉位置倒圓半徑R=6.5mm，前緣鈍化半徑r=2mm，前緣擴張角β=48°。前緣平直段的長度L=31r，以使其產(chǎn)生的脫體激波充分發(fā)展。以V形前緣駐點位置為坐標(biāo)原點，定義流向為x方向，平板展向為z方向，垂直于平板方向為y方向。

圖1 V形鈍前緣平板示意圖

為了更加清楚地觀察波系結(jié)構(gòu)，對V形鈍前緣平板實驗?zāi)Ｐ瓦M(jìn)行了2.5倍放大。實驗在中國科學(xué)技術(shù)大學(xué)的KDJB330反射型激波風(fēng)洞[18-20]中進(jìn)行。實驗來流馬赫數(shù)為6，靜壓為891Pa，靜溫為101K。紋影系統(tǒng)采用Z字形光路布置，如圖2所示。通過高速攝影記錄流場紋影圖像,拍攝幀率為15000Hz，曝光時間為2μs。入射斜激波由氣流偏轉(zhuǎn)角為12°的斜劈產(chǎn)生，斜劈流向長度為200mm，展向?qū)挾葹?60mm。為了獲得豐富的流場信息，考慮到模型對光路的遮擋，分別從平行于平板(z方向，見圖2)和垂直于平板(y方向，見圖3)2個方向?qū)α鲌鲞M(jìn)行拍攝。通過前后移動斜劈的方式，改變?nèi)肷湫奔げㄅcV形鈍前緣溢流口的相對位置。為了便于分析，斜激波入射位置采用不考慮激波干擾時無粘斜激波在平板上的入射點坐標(biāo)。

圖2 紋影光路圖

1.2 計算方法

數(shù)值計算采用Fluent軟件求解Navier-Stokes方程，并使用具有二階精度的格式離散。對流通量采用AUSM格式，湍流模型采用k-ωSST模型，氣體分子粘性采用Sutherland公式，氣體定壓比熱使用CHEMKIN溫度多項式進(jìn)行計算。對于V形鈍前緣平板流場的計算，由于來流和幾何的對稱性，為了節(jié)省計算資源，采用1/4計算域進(jìn)行模擬。采用結(jié)構(gòu)化六面體網(wǎng)格，計算域及局部網(wǎng)格如圖4和5所示，近壁面采用等比加密，壁面y+<1，總網(wǎng)格量約3000萬。有斜激波入射時的計算域如圖6所示，此時六面體網(wǎng)格總數(shù)在5000萬左右。為了與風(fēng)洞實驗進(jìn)行對比，來流條件設(shè)置與實驗一致。在數(shù)值模擬中監(jiān)視各方程殘差，監(jiān)測出口位置的流量及質(zhì)量平均馬赫數(shù)、駐點位置的壁面壓力和熱流量。待殘差收斂并且各監(jiān)測參數(shù)穩(wěn)定，認(rèn)為流場收斂。

由于三維流動的復(fù)雜性，數(shù)值模擬獲得的流場結(jié)構(gòu)難以直接與具有光路積分效應(yīng)的實驗紋影進(jìn)行對比。為此，本文按照紋影法的原理對數(shù)值模擬流場進(jìn)行處理。如為了獲得z方向穿過全流場的數(shù)值紋影，將流場沿z方向劃分為k個等間距的區(qū)域，假設(shè)每個小區(qū)域內(nèi)的流場參數(shù)沿z方向均勻分布，則經(jīng)過流場的光線偏轉(zhuǎn)角εx由式(1)[21]得到：

圖5 V形鈍前緣平板計算網(wǎng)格

圖6 斜激波入射V形鈍前緣平板計算域

Fig.6ComputationaldomainoftheV-shapedbluntleadingedgeplatewithshockincidence

ξ2-ξ1)(1)

式中，n為折射率，ξ1和ξ2分別為光線進(jìn)入流場和離開流場位置的z坐標(biāo)。由此獲得的光強變化可以與z方向穿過流場的實驗紋影進(jìn)行直觀對比。

2 結(jié)果與分析

V形鈍前緣自身的三維繞流會導(dǎo)致復(fù)雜的激波干擾，因此在研究斜激波入射V形鈍前緣平板激波干擾前，有必要對V形鈍前緣基本波系結(jié)構(gòu)進(jìn)行研究。

2.1 V形鈍前緣平板基本流場

從y方向拍攝的流場紋影與流場x-z對稱面數(shù)值紋影的對比如圖7(a)所示，可以看到數(shù)值模擬較好地捕捉到了基本的波系結(jié)構(gòu)。圖7(b)為上半部分波系結(jié)構(gòu)的示意圖，其中DS代表V形鈍前緣斜直段產(chǎn)生的脫體激波，BS代表V形前緣交叉位置形成的弓形激波，脫體激波DS入射弓形激波BS，相交于三波點T，在三波點位置產(chǎn)生透射激波TS和滑移線Σ。脫體激波DS波后的超聲速氣流在V形前緣交叉倒圓位置被壓縮，壓縮波匯聚形成彎曲激波CS。透射激波TS與彎曲激波CS發(fā)生異側(cè)激波干擾，形成馬赫桿MS。激波BS以及馬赫桿MS波后氣流均為亞聲速，而透射激波TS波后氣流為超聲速，超聲速氣流在亞聲速區(qū)域之間不斷地壓縮和膨脹，交替出現(xiàn)激波和膨脹波，形成超聲速射流結(jié)構(gòu)。超聲速射流沿壁面運動，在駐點前發(fā)生碰撞，并卷起形成反向旋轉(zhuǎn)的渦對結(jié)構(gòu)。渦對結(jié)構(gòu)與弓形激波BS相互作用，使得弓形激波波面向上游凸起，并出現(xiàn)Kink。將脫體激波DS視為入射激波，則V形鈍前緣激波干擾結(jié)構(gòu)與Edney第V類激波干擾類似。流場中存在射流失穩(wěn)和剪切層脫渦等流動過程，數(shù)值模擬雖然不能捕捉這些非定常過程，但能夠刻畫基本波系結(jié)構(gòu)，流動的非定常性沒有顯著影響本文對波系干擾類型的研究。

圖7 x-z對稱面波系結(jié)構(gòu)

為了探討流場波系的三維結(jié)構(gòu)，從z方向也進(jìn)行了紋影拍攝。使用本文1.2節(jié)所描述的三維流場后處理方法得到的數(shù)值紋影與實驗紋影圖像進(jìn)行了對比，如圖8(a)所示。由于流場的三維性，在z方向紋影圖像上可以看到3道激波，從左到右依次為弓形激波BS上凸起的激波面(Bulge)、弓形激波BS和彎曲激波CS。圖8(b)進(jìn)一步定量比較了數(shù)值模擬與實驗獲得的激波BS及其凸起的激波面的位置。結(jié)果表明：激波BS的位置符合較好，而凸起的激波面(Bulge)略有差別。這是由于激波BS相對穩(wěn)定，而凸起的激波面受渦對結(jié)構(gòu)的影響(見圖7)，具有一定的非定常性。整體而言，數(shù)值模擬與實驗紋影吻合較好，可用于本文后續(xù)的激波干擾研究。

圖8 z方向波系結(jié)構(gòu)比較

圖9比較了相同鈍化半徑下二維圓柱產(chǎn)生的弓形激波和V形鈍前緣平板弓形激波。其中，Δ1和Δ2分別代表V形鈍前緣和圓柱弓形激波的脫體距離，V形鈍前緣弓形激波由x-y對稱面的數(shù)值模擬給出，圓柱產(chǎn)生的弓形激波形狀由Billig[22]公式給出，而弓形激波面上聲速點位置的激波角β*由以下公式[23]給出：

式中，Ma∞為來流馬赫數(shù),γ為氣體比熱比。從圖9可以看出，由于自身的激波干擾，V形鈍前緣交叉位置弓形激波的脫體距離較大，并且駐點前存在大范圍的亞聲速區(qū)。進(jìn)一步采用馬赫數(shù)為5.95的等值面給出V形鈍前緣平板的三維激波結(jié)構(gòu)(圖10)，可以看出：前緣脫體激波DS和弓形激波BS均為三維曲面激波，并且BS波面有較大范圍的凸起。這表明V形鈍前緣弓形激波的近正激波區(qū)域較大，當(dāng)斜激波與弓形激波干擾時，出現(xiàn)第Ⅲ和第Ⅳ類干擾的范圍更大[24]。

圖9 V形鈍前緣脫體激波與圓柱脫體激波對比

Fig.9ComparisonbetweenthebowshockofthecylinderandV-shapedbluntleadingedge

圖10 流場三維激波面

2.2 斜激波入射V形鈍前緣溢流口激波干擾流場

當(dāng)入射斜激波距離V形溢流口上游較遠(yuǎn)時，對應(yīng)于亞額定工況下內(nèi)轉(zhuǎn)式進(jìn)氣道前體激波與唇口激波相交于唇口外的情況。圖11(a)為斜激波無粘入射平板位置坐標(biāo)為x/r=-10.5時的流場紋影圖片，圖11(b)為波系示意圖。斜激波IS從平板下側(cè)入射弓形激波BS超聲速區(qū)域，發(fā)生同側(cè)激波相交，出現(xiàn)Edney第Ⅵ類激波干擾。激波相互作用導(dǎo)致弓形激波的激波角減小，入射激波干擾后形成透射激波TS和剪切層Σ。

圖11 斜激波入射x/r=-10.5流場波系結(jié)構(gòu)

Fig.11Shockstructureoftheflowfieldwiththeshockincidenceatx/r=-10.5

向下游移動斜劈，當(dāng)斜激波無粘入射位置為x/r=-3時，斜激波入射在V形前緣駐點前的亞聲速區(qū)，此時對應(yīng)于內(nèi)轉(zhuǎn)式進(jìn)氣道額定工況下前體激波入射唇口的情況。僅從z方向拍攝紋影時，由于模型遮擋，難以獲得干擾點附近激波干擾形態(tài)。因此，需要結(jié)合y方向的紋影來判斷干擾類型。從圖12(a)流場z方向的紋影可以看出,入射斜激波IS分為了兩道斜激波，分別為IS和IS′。其中，IS′是入射斜激波IS與V形鈍前緣的脫體激波DS干擾產(chǎn)生的透射激波。IS′入射V形前緣平板壁面形成反射激波RS，而入射斜激波IS與弓形激波BS的干擾點被模型遮擋。從圖12(b)流場y方向的紋影可以看到，在駐點前存在兩道弓形激波BS和BS′。結(jié)合圖12(a)和(b)兩個方向的紋影，推測激波干擾類型可能為Edney第Ⅳ類激波干擾，其流場結(jié)構(gòu)如圖13所示。此時，由于斜激波IS入射，原有的弓形激波BS轉(zhuǎn)變?yōu)橥ㄟ^透射激波TS連接的兩段弓形激波BS和BS′，并伴隨著超聲速射流的產(chǎn)生，該流場結(jié)構(gòu)很復(fù)雜，將在2.3節(jié)進(jìn)一步分析。因此，在內(nèi)轉(zhuǎn)式進(jìn)氣道額定工況下，前體激波入射溢流口駐點前亞聲速區(qū)域，可能帶來復(fù)雜的激波干擾流場，對進(jìn)氣道流動帶來負(fù)面影響。

圖12 入射位置x/r=-3流場紋影

Fig.12Schlierenimagesoftheflowfieldwiththeshockincidenceatx/r=-3

圖13 入射位置x/r=-3 流場示意圖

Fig.13Schematicoftheflowfieldwiththeshockincidenceatx/r=-3

進(jìn)一步向下游移動斜劈，使得斜激波無粘入射位置為x/r=0.2，此時對應(yīng)于內(nèi)轉(zhuǎn)式進(jìn)氣道超額定狀態(tài)下前體激波入射溢流口內(nèi)部的情況。流場紋影如圖14(a)所示，圖14(b)為其流場示意圖。入射激波IS與弓形激波BS干擾形成的剪切層Σ并沒有撞擊壁面，而是形成一種有別于經(jīng)典的Edney第Ⅲ類激波干擾的結(jié)構(gòu)。其原因在于：V形鈍前緣脫體激波DS與弓形激波BS干擾形成的脫體距離和亞聲速區(qū)域較大，使得發(fā)生第Ⅲ類激波干擾的區(qū)域增大[23]。

圖14 入射位置x/r=0.2流場紋影

Fig.14Schlierenimagesoftheflowfieldwiththeshockincidenceatx/r=0.2

當(dāng)斜激波無粘入射位置為x/r=3時，斜激波與弓形激波BS的超聲速區(qū)域干擾，流場紋影如圖15(a)所示，流場如圖15(b)所示。入射斜激波IS與弓形激波BS相交，形成長度較短的馬赫桿MS、透射激波TS以及剪切層Σ1和Σ2。這種由Edney第Ⅱ類激波干擾形成的亞聲速區(qū)，可能會帶來進(jìn)氣道局部不起動[16-17]，應(yīng)避免出現(xiàn)。

圖15 入射位置x/r=3流場紋影

Fig.15Schlierenimagesoftheflowfieldwiththeshockincidenceatx/r=3

2.3 斜激波/弓形激波第Ⅳ類激波干擾流場分析

當(dāng)斜激波入射V形鈍前緣時，三維波系干擾的情況如圖16所示(這里采用斜激波在上)。根據(jù)激波位置和特點，可能會出現(xiàn)3種干擾情況：第一種是入射斜激波IS與脫體激波DS相交，第二種是入射斜激波IS與弓形激波BS相交，第三種是入射斜激波與激波BS和DS的交線相交。鑒于三維激波干擾的復(fù)雜性，為了深入理解斜激波入射V形前緣時的激波干擾流場，本節(jié)以斜激波無粘入射壁面位置坐標(biāo)為x/r=-5時作為典型工況，結(jié)合數(shù)值模擬進(jìn)行分析。與入射位置為x/r=-3的工況(圖12)相比，此時激波干擾點位置更接近BS的正激波部分。

圖16 斜激波入射V形鈍前緣三維流場示意圖

Fig.16SchematicoftheshockinteractionflowfieldoftheV-shapedbluntleadingedgewithshockincidence

圖17(a)采用流向截面和展向?qū)ΨQ面(z/r=0)展示了斜激波IS與V形鈍前緣波系干擾三維結(jié)構(gòu)的演化情況。其中，流向截面的馬赫數(shù)云圖僅顯示了一半流場，而展向?qū)ΨQ面通過馬赫數(shù)等值線顯示。在x/r=-5截面上(局部放大見圖17(b))，入射斜激波IS與脫體激波DS的較弱部分發(fā)生異側(cè)相交，形成Edney第Ⅰ類干擾，產(chǎn)生的透射激波TS相對于激波IS向DS激波面偏折。在x/r=-2.5截面上(局部放大見圖17(c))，入射斜激波IS與脫體激波DS的較強部分發(fā)生異側(cè)相交，形成Edney第Ⅱ類激波干擾，并產(chǎn)生馬赫桿MS，此時的透射激波TS1入射鈍前緣壁面。圖17沿流向發(fā)展的波系結(jié)構(gòu)表明：由于入射斜激波IS與前緣脫體激波DS相交處DS激波強度的變化，激波干擾類型由Edney第Ⅰ類轉(zhuǎn)變?yōu)榈冖蝾悺６赩形鈍前緣溢流口處，如圖17(a)中的藍(lán)色線框所示，入射激波IS與弓形激波BS發(fā)生干擾并產(chǎn)生超聲速射流，這一區(qū)域的流場結(jié)構(gòu)將在圖18中進(jìn)行分析說明。

在斜激波入射條件下，多股氣流共同匯集在V形鈍前緣交叉位置，形成更加復(fù)雜的三維激波干擾流場。圖18給出了V形鈍前緣交叉位置z/r=0和y/r=-0.4兩個正交截面上的流場結(jié)構(gòu)。從z/r=0截面可以看出，入射斜激波與弓形激波BS干擾后，形成通過透射激波TS4連接的兩段弓形激波BS和BS′以及超聲速射流。該超聲速射流并未直接入射壁面，而是發(fā)生高度彎曲沿壁面方向溢流，此干擾類型是第Ⅳa類干擾。而y/r=-0.4截面展示了斜激波IS同時與弓形激波BS和脫體激波DS發(fā)生多波系干擾的情況，由于IS與DS干擾已經(jīng)形成了馬赫桿MS(見圖17(c))，MS再與弓形激波BS干擾形成透射激波TS3以及超聲速射流區(qū)域。可見，入射激波IS、弓形激波BS與脫體激波DS三者之間的干擾互相耦合。

(a) 三維流場

(b) x/r=-5截面

(c) x/r=-2.5截面

Fig.17ShockinteractionsbetweenISandtheV-shapedbluntleadingedge

將數(shù)值模擬流場沿y方向進(jìn)行三維紋影處理，如圖19(a)所示，與圖19(b)給出的實驗紋影圖像符合較好。從三維紋影結(jié)果中可以看到，V形鈍前緣溢流口弓形激波BS被分為兩道激波BS和BS′。實驗紋影圖像與數(shù)值計算結(jié)果相互印證了斜激波與V形鈍前緣弓形激波發(fā)生了第Ⅳa類激波干擾。

Fig.18ShockinteractionsnearthecrotchoftheV-shapedbluntleadingedgewithshockincidence

Fig.19Comparisonbetweentheexperimentalschlierenandthethree-dimensionalnumericalschlierenviewingfromtheydirection

3 結(jié) 論

通過簡化構(gòu)型對內(nèi)轉(zhuǎn)式進(jìn)氣道V形鈍前緣溢流口激波干擾特性進(jìn)行了研究，得到以下結(jié)論：

(1) V形鈍前緣自身形成的激波干擾，導(dǎo)致其駐點前弓形激波的脫體距離較大，波后存在大范圍的亞聲速區(qū)，因此在估算內(nèi)轉(zhuǎn)式進(jìn)氣道溢流口激波脫體距離時，需慎用圓柱脫體激波距離計算公式。

(2) 斜激波入射在V形鈍前緣駐點前弓形激波接近正激波的位置時，發(fā)生Edney第Ⅳa類激波干擾。此時，入射激波與弓形激波BS的干擾以及弓形激波BS與脫體激波DS的干擾互相耦合，形成多處超聲速射流區(qū)域。因此，在內(nèi)轉(zhuǎn)式進(jìn)氣道額定工況下，前體激波入射V形溢流口駐點前亞聲速區(qū)域可能帶來復(fù)雜的流場結(jié)構(gòu)，在實際應(yīng)用中應(yīng)當(dāng)予以關(guān)注。

(3) 當(dāng)斜激波入射在V形鈍前緣弓形激波亞聲速區(qū)的聲速點附近時，出現(xiàn)一種有別于Edney第Ⅲ類激波干擾的波系結(jié)構(gòu)，此時干擾產(chǎn)生的剪切層未入射壁面。

(4) 入射斜激波與V形鈍前緣脫體激波DS、弓形激波BS的超聲速區(qū)域干擾時的形態(tài)與Edney第Ⅰ、Ⅱ、Ⅵ類激波干擾類似，其中第Ⅱ類干擾形成的馬赫桿后局部亞聲速區(qū)域，可能造成內(nèi)轉(zhuǎn)式進(jìn)氣道局部不起動。

鑒于V形溢流口三維流動的復(fù)雜性，在激波干擾引起的氣動熱以及三維前體激波入射方面，還有待進(jìn)一步的細(xì)致研究。