魏祥庚，曹東剛，秦 飛，吳繼平

(1.高超聲速?zèng)_壓發(fā)動(dòng)機(jī)技術(shù)重點(diǎn)實(shí)驗(yàn)室(國防科學(xué)技術(shù)大學(xué))，長沙 410073；2.燃燒、熱結(jié)構(gòu)與內(nèi)流場(chǎng)重點(diǎn)實(shí)驗(yàn)室(西北工業(yè)大學(xué))，西安 710072)

受限空間內(nèi)超聲速混合層生長特性

魏祥庚1,2，曹東剛2，秦 飛2，吳繼平1

(1.高超聲速?zèng)_壓發(fā)動(dòng)機(jī)技術(shù)重點(diǎn)實(shí)驗(yàn)室(國防科學(xué)技術(shù)大學(xué))，長沙 410073；2.燃燒、熱結(jié)構(gòu)與內(nèi)流場(chǎng)重點(diǎn)實(shí)驗(yàn)室(西北工業(yè)大學(xué))，西安 710072)

為獲得受限空間內(nèi)激波作用下的超聲速混合層生長規(guī)律，以支板噴射超燃沖壓發(fā)動(dòng)機(jī)典型流道為研究對(duì)象，開展了2.3Ma氫氣射流與2.0Ma空氣來流所形成的超聲速混合層的生長特性研究.基于OpenFOAM計(jì)算平臺(tái)，采用大渦模擬方法，數(shù)值研究了超聲速混合層的流場(chǎng)結(jié)構(gòu)和特征，流場(chǎng)結(jié)構(gòu)和組分分布與實(shí)驗(yàn)結(jié)果吻合較好.通過超聲速混合層組分濃度、厚度、可壓縮效應(yīng)及總壓損失的分析，獲得了超聲速混合層的生長特性.研究結(jié)果表明：受限空間內(nèi)超聲速混合層的生長過程具有4個(gè)典型階段，支板末端的膨脹波/激波結(jié)構(gòu)會(huì)顯著減低對(duì)流馬赫數(shù)，從而降低混合層的可壓縮性，促進(jìn)混合層的生長；激波與混合層的相互作用能夠增強(qiáng)局部湍流強(qiáng)度，獲得渦量增益，加快混合層的生長速率，促進(jìn)混合效率，但同時(shí)會(huì)引起較大的總壓損失，降低發(fā)動(dòng)機(jī)性能.發(fā)動(dòng)機(jī)設(shè)計(jì)時(shí)要綜合考慮波系結(jié)構(gòu)與混合層相互作用帶來的混合增強(qiáng)和總壓損失，實(shí)現(xiàn)性能優(yōu)化.

超聲速混合層；激波；波系結(jié)構(gòu)；受限空間；可壓縮性；大渦模擬

湍流摻混是動(dòng)力系統(tǒng)中最為常見且極為重要的流動(dòng)現(xiàn)象，摻混過程的優(yōu)化組織可以強(qiáng)化燃料與氧氣的混合(非預(yù)混燃燒)或是促進(jìn)燃?xì)馀c周圍氣流的熱量交換(預(yù)混燃燒)，對(duì)于縮短燃燒室長度、提高燃燒效率、優(yōu)化發(fā)動(dòng)機(jī)能量管理有著重要的指導(dǎo)意義[1-2].隨著高超聲速飛行技術(shù)的發(fā)展，高速可壓流動(dòng)中的湍流摻混成了該領(lǐng)域的研究熱點(diǎn)，也是高超聲速推進(jìn)系統(tǒng)中的關(guān)鍵技術(shù)之一.德國宇航院的Papamoschou等[3]以支板噴射超燃沖壓發(fā)動(dòng)機(jī)為對(duì)象，開展了激波作用下混合層結(jié)構(gòu)的實(shí)驗(yàn)研究，研究指出由于高速條件下氣流的可壓縮性增強(qiáng)，使得混合層的生長速率減弱，湍流強(qiáng)度與雷諾應(yīng)力也相應(yīng)減小，混合層的穩(wěn)定性增強(qiáng)，并建議采用對(duì)流馬赫數(shù)Mc作為研究參數(shù)來表征混合層的可壓縮性.Furby等[4]以HyShot為研究對(duì)象，分析了混合層發(fā)展過程中的旋渦類型及演化過程.Gutmark等[5]分析了自由空間內(nèi)可壓縮混合層生長速率與對(duì)流馬赫數(shù)之間的關(guān)系，指出當(dāng)對(duì)流馬赫數(shù)增大時(shí)，可壓縮混合層的生長速率趨于不可壓縮混合層生長速率的20%.文獻(xiàn)[6-8]對(duì)比了自由空間中不同對(duì)流馬赫數(shù)條件下混合層的生長過程，研究結(jié)果均表明隨著對(duì)流馬赫數(shù)的增加，當(dāng)混合層中的大尺度結(jié)構(gòu)相對(duì)于當(dāng)?shù)亓鲃?dòng)的馬赫數(shù)大于1時(shí)，當(dāng)?shù)亓鲃?dòng)受到大尺度結(jié)構(gòu)的阻礙而產(chǎn)生激波，混合層中會(huì)出現(xiàn)小激波結(jié)構(gòu).文獻(xiàn)[9-11]研究結(jié)果表明當(dāng)激波打到混合層上之后能夠強(qiáng)化流體微團(tuán)的旋轉(zhuǎn)運(yùn)動(dòng)，增強(qiáng)局部的組分輸運(yùn)，進(jìn)而促進(jìn)混合過程.Zhang等[12]以自由空間內(nèi)的可壓混合層為研究對(duì)象，討論了激波對(duì)混合層發(fā)展過程的影響，研究表明在激波/混合層的作用點(diǎn)附近，由于激波的壓縮作用會(huì)使得混合層厚度降低，但在作用點(diǎn)下游，由于渦量的增加使得混合層厚度顯著增大.

對(duì)于高超聲速飛行器而言，由于燃燒室進(jìn)口氣流為超聲速，燃料與空氣需要在有限的空間和極短的時(shí)間內(nèi)完成高效摻混，進(jìn)而組織燃燒.與自由空間內(nèi)高速可壓流動(dòng)中的湍流摻混相比，高超聲速推進(jìn)系統(tǒng)中的混合層發(fā)展不僅受其自身的可壓縮性影響，同時(shí)也受到流道幾何結(jié)構(gòu)的約束，其生長特性更為復(fù)雜.此外，由于沖壓流道內(nèi)往往具有復(fù)雜的波系結(jié)構(gòu)，波系/混合層的相互作用也是混合層生長的重要影響因素.為了研究受限空間內(nèi)超聲速混合層的生長特性，本文以支板噴射超燃沖壓發(fā)動(dòng)機(jī)典型流道為研究對(duì)象，基于計(jì)算軟件(open field operation and manipulation, OpenFOAM)，采用大渦模擬(large eddy simulation, LES)對(duì)2.3Ma氫氣射流與2.0Ma空氣來流所形成的超聲速混合層進(jìn)行了研究，獲得了混合層的流動(dòng)結(jié)構(gòu)和流場(chǎng)特征，分析了受限空間內(nèi)超聲速混合層的生長特性，有助于加深對(duì)受限空間內(nèi)超聲速混合層特征的認(rèn)識(shí)和對(duì)復(fù)雜波系作用下湍流摻混過程的理解.

1 物理模型和計(jì)算方法

1.1 物理模型

本文以支板噴射超燃沖壓發(fā)動(dòng)機(jī)為研究對(duì)象，其流道構(gòu)型如圖1所示[13]，定義空氣設(shè)備噴管喉部為x=0，支板前端位于x=95 mm處，支板末端位于x=175 mm處.冷態(tài)空氣經(jīng)噴管加速后以2.0Ma的速度進(jìn)入燃燒室內(nèi)，中心支板對(duì)來流空氣進(jìn)行適當(dāng)?shù)膲嚎s，并為燃料噴注提供載體，氫氣通過支板內(nèi)置的噴管加速后以2.3Ma的速度進(jìn)入燃燒室，與空氣來流形成超聲速混合層.設(shè)備噴管的喉部高度為20.12 mm，燃燒室入口截面為35.4 mm×40.0 mm，中心支板全長80 mm，支板末端高度為2 mm，氫氣噴嘴出口高度為0.6 mm.氫氣流量和空氣流量分別為2.4、740 g/s.

圖1 流道構(gòu)型示意

1.2 計(jì)算方法

大渦模擬結(jié)合了直接數(shù)值模擬的準(zhǔn)確性和雷諾時(shí)均模式的快速性，計(jì)算精度較高，是湍流摻混燃燒數(shù)值模擬的有效方法，LES作為研究燃燒流動(dòng)的有效手段已經(jīng)得到了廣泛的認(rèn)可，并且逐漸被用于研究各種類型發(fā)動(dòng)機(jī)的燃燒流動(dòng)細(xì)節(jié).描述湍流混合過程的基本控制方程為可壓Navier-Stokes方程組，對(duì)連續(xù)方程、動(dòng)量方程、能量方程、組分方程采用Favre平均并進(jìn)行濾波運(yùn)算后可以得到描述湍流摻混的LES控制方程組[14-15]：

式中：“-”為空間濾波；“～”為采用Favre平均后大于濾波尺度的求解參數(shù)項(xiàng)；上標(biāo)“sgs”為小于濾波尺度的亞格子模型求解參數(shù)項(xiàng)；ρ為密度；uj為速度；p為壓力；δij為單位張量函數(shù)；τij為應(yīng)力項(xiàng)；E為能量；Ym為組分m的質(zhì)量分?jǐn)?shù)；q為熱通量；θ為擴(kuò)散通量.在LES的求解過程中，對(duì)于大尺度旋渦進(jìn)行直接求解，對(duì)于比濾波尺度小的流動(dòng)結(jié)構(gòu)通過亞格子模型進(jìn)行模擬.本文以O(shè)penFOAM為計(jì)算平臺(tái)，使用對(duì)于可壓縮流動(dòng)具有較好適用性的rhoReactingFoam求解器進(jìn)行計(jì)算，采用Menon等[16]使用的oneEqEddy亞格子模型對(duì)亞格子項(xiàng)進(jìn)行封閉：

空間離散采用二階精度的TVD格式，考慮背壓對(duì)于流動(dòng)的影響，空氣來流與燃料射流均從噴管喉部開始進(jìn)行計(jì)算，計(jì)算入口的邊界參數(shù)(壓力、速度、溫度、組分)全部給定，見表1.流道出口仍以超聲速流動(dòng)為主，根據(jù)特征線理論，超聲速出口參數(shù)可以根據(jù)場(chǎng)內(nèi)信息進(jìn)行外插計(jì)算得到.時(shí)間離散采用二階隱式Crank-Nicholson格式，時(shí)間步長選取為1×10-8s，計(jì)算的最大CFL數(shù)不超過0.3.網(wǎng)格劃分采用結(jié)構(gòu)化網(wǎng)格，對(duì)氫氣射流出口和燃燒室中心區(qū)域的網(wǎng)格進(jìn)行了加密以便更好地捕捉混合層中的流動(dòng)結(jié)構(gòu).為了更為真實(shí)地模擬受限空間內(nèi)的流動(dòng)狀況，采用無滑移絕熱壁面，結(jié)合大渦模擬對(duì)壁面處網(wǎng)格的要求，通過附面層對(duì)壁面網(wǎng)格質(zhì)量進(jìn)行了提升，壁面處第1層網(wǎng)格的厚度約為0.01 mm，從而使得y+≤1，并使用文獻(xiàn)[13]中通過實(shí)驗(yàn)和計(jì)算所提取出的二維平面特征進(jìn)行計(jì)算，網(wǎng)格總數(shù)為74萬，燃燒室計(jì)算區(qū)域和網(wǎng)格劃分如圖2所示.

表1 空氣來流和氫氣射流的入口邊界條件

圖2 計(jì)算網(wǎng)格

2 結(jié)果與分析

在支板噴射超燃沖壓發(fā)動(dòng)機(jī)中，空氣來流經(jīng)中心支板適當(dāng)壓縮后仍以超聲速進(jìn)入燃燒室，燃料經(jīng)過支板內(nèi)置的噴管膨脹加速后射入燃燒室，空氣來流與氫氣射流形成超聲速混合層，進(jìn)而主導(dǎo)燃燒室內(nèi)的湍流摻混過程.圖3給出了實(shí)驗(yàn)紋影圖[13]，計(jì)算所得的數(shù)值紋影結(jié)果以及5 000個(gè)時(shí)間步統(tǒng)計(jì)平均紋影結(jié)果，可以看出計(jì)算捕捉到了與實(shí)驗(yàn)觀察結(jié)果相一致的主要流場(chǎng)特征：由于中心支板的擠壓作用，超聲速空氣來流進(jìn)入流道后會(huì)在支板前緣形成斜激波，由于受到?jīng)_壓流道自身幾何構(gòu)型的約束，激波在燃燒室壁面進(jìn)行反射并向下游傳播；中心支板的菱形結(jié)構(gòu)使得支板前部形成壓縮型面，后部形成擴(kuò)張型面，因此會(huì)在轉(zhuǎn)折點(diǎn)(x=135 mm)處形成膨脹波；氣流進(jìn)入燃燒室后，由于流道的突然擴(kuò)展，會(huì)在支板后緣形成一系列的膨脹波系，氣流經(jīng)過支板尾部的膨脹波系后向噴流中心偏轉(zhuǎn)，進(jìn)而在噴嘴附近又形成兩道斜激波；激波在向下游的傳播過程中，不僅受流道約束會(huì)在固壁面進(jìn)行反射，而且會(huì)與混合層相互作用，整個(gè)流場(chǎng)表現(xiàn)出以混合層為主導(dǎo)且伴隨多種波系結(jié)構(gòu)的復(fù)雜特征.

圖3 實(shí)驗(yàn)紋影圖、數(shù)值紋影圖與平均紋影圖對(duì)比

Fig.3 Comparisons of experimental schlieren photograph, averaged numerical schlieren image, and instantaneous numerical schlieren image

圖4給出了沿流向4個(gè)不同位置處H2的摩爾分?jǐn)?shù)分布曲線.

圖4 不同位置處H2的摩爾分?jǐn)?shù)分布實(shí)驗(yàn)與計(jì)算對(duì)比

Fig.4 Calculated and measured H2profiles at different locations

從圖4中可以看出計(jì)算曲線與實(shí)驗(yàn)數(shù)據(jù)吻合良好，驗(yàn)證了計(jì)算方法的合理性與可行性.流道內(nèi)混合層的發(fā)展主導(dǎo)著混合過程的進(jìn)行，湍流摻混首先發(fā)生在氫氣射流與空氣來流相接觸的剪切層中，伴隨著剪切層的生長過程，氫氣射流與空氣來流實(shí)現(xiàn)質(zhì)量、動(dòng)量和能量的交換.在初始階段(x=178 mm)各組分分布區(qū)域較為固定且燃料區(qū)與氧化劑區(qū)域之間存在著巨大的濃度梯度，燃料經(jīng)支板噴射進(jìn)入燃燒室后與空氣來流之間形成速度不連續(xù)的間斷面，由于間斷面的不穩(wěn)定性，剪切層產(chǎn)生波動(dòng)進(jìn)而發(fā)展成為旋渦，強(qiáng)化湍流.伴隨著流動(dòng)的發(fā)展，氫氣射流與空氣來流相互卷吸，射流與周圍流體的摻混自邊緣逐漸向中心發(fā)展，混合層逐漸向兩側(cè)擴(kuò)展，混合層厚度逐漸增大，當(dāng)混合層充分發(fā)展以后(x=325 mm)，中心射流與周圍流體之間大的參數(shù)梯度逐漸被抹平，主要組分的濃度曲線逐漸變得平穩(wěn).

為了進(jìn)一步描述復(fù)雜波系作用下受限空間內(nèi)超聲速混合層的生長特性，定義混合層厚度δ作為表征參數(shù)[17]為

δ(x)=R1(x)φ=0.5-R2(x)φ=1.5,

其含義為局部當(dāng)量比φ=0.5(氫氣質(zhì)量分?jǐn)?shù)YH2=1.431%)的等值線到局部當(dāng)量比φ=1.5(氫氣質(zhì)量分?jǐn)?shù)YH2=4.25%)的等值線的距離.圖5給出了混合層厚度沿流道方向的分布曲線，從圖5中可以看出，混合層的生長過程可以分為4個(gè)階段：回流生長區(qū)(175 mm

圖5 混合層厚度沿流道方向的分布曲線

圖6 x方向速度云圖

圖7給出了流場(chǎng)中通過壓力梯度表征的波系結(jié)構(gòu)，圖8給出了流場(chǎng)中相應(yīng)的渦量分布，當(dāng)激波打到混合層上以后(如x=223 mm)，由于壓縮作用會(huì)使得混合層的厚度有所減低，與此同時(shí)入射激波會(huì)迫使外側(cè)的空氣來流進(jìn)入內(nèi)側(cè)的燃料射流中，而且激波引起的局部逆壓梯度會(huì)強(qiáng)化剪切層上、下兩側(cè)流體微團(tuán)的旋轉(zhuǎn)運(yùn)動(dòng)，誘導(dǎo)形成渦量增益，增強(qiáng)兩股氣流在垂直方向的輸運(yùn)過程.入射激波引起的劇烈擾動(dòng)能夠激勵(lì)剪切層K-H不穩(wěn)定性，會(huì)增大局部的湍流強(qiáng)度，獲得渦量增益，增大混合層中旋渦結(jié)構(gòu)的卷起與合并速度，促進(jìn)混合層上、下兩側(cè)燃料與空氣的摻混，使得臨近作用點(diǎn)下游的混合層生長速率增大.需要指出的是激波/混合層的相互作用受到激波強(qiáng)度與混合層內(nèi)外兩側(cè)氣體參數(shù)梯度的共同影響，隨著流動(dòng)向下游的發(fā)展，激波強(qiáng)度逐漸減弱，混合層內(nèi)外兩側(cè)氣體的參數(shù)梯度也逐漸減小，當(dāng)激波強(qiáng)度較大且氣體參數(shù)梯度較小時(shí)，激波會(huì)穿透混合層，表現(xiàn)出折射特征，對(duì)混合層的壓縮作用也較為顯著，但整體上還是可以促進(jìn)混合過程的進(jìn)行，這與文獻(xiàn)[15]中觀測(cè)的趨勢(shì)一致.在快速生長區(qū)之后，混合層經(jīng)過充分發(fā)展，進(jìn)入飽和生長區(qū)，混合層的生長速率減慢，其厚度趨于穩(wěn)定，標(biāo)志著混合過程的基本結(jié)束.

圖7 壓強(qiáng)梯度云圖

圖8 渦量分布云圖

根據(jù)文獻(xiàn)[15]的分析，該工況中氫氣射流與空氣來流的對(duì)流Ma>1[13]，這就意味著可壓縮效應(yīng)比較顯著，混合層的生長速度應(yīng)該比較緩慢，實(shí)現(xiàn)充分摻混需要更長的距離.然而，從圖5中可以看出，初始生長區(qū)中混合層的生長速率的確比較緩慢，但是經(jīng)過很短的距離后混合層就開始迅速發(fā)展，進(jìn)入快速生長區(qū)，混合層的生長速率顯著增大.為了分析這一現(xiàn)象，圖9、10分別給出了流場(chǎng)中時(shí)均化的馬赫數(shù)和聲速分布云圖.氫氣射流通過燃料噴嘴以2.3Ma的速度進(jìn)入燃燒室，與2.0Ma的空氣來流進(jìn)行摻混，從圖9可以看出，在支板末端形成了亞聲速的回流區(qū)，有利于燃料與空氣的摻混.隨著流動(dòng)向下游的發(fā)展，亞聲速的回流區(qū)逐漸減小，兩股氣流通過界面處的剪切作用進(jìn)行動(dòng)量交換，混合區(qū)域內(nèi)呈現(xiàn)出以跨聲速流動(dòng)為主的流場(chǎng)特征.

圖9 時(shí)均化馬赫數(shù)云圖

圖10 時(shí)均化聲速云圖

根據(jù)聲速的定義為

可知?dú)怏w中聲速c的大小直接代表了氣體的可壓縮性的大小，氣體中的聲速越大，則氣體的可壓縮性就越小.從圖10中可以看出，氫氣射流與空氣來流通過混合層實(shí)現(xiàn)摻混，混合核心區(qū)的聲速相對(duì)于氫氣射流的聲速有所降低，相對(duì)于空氣來流的聲速卻有所提高，這表明在剪切層內(nèi)混合氣體的可壓縮性發(fā)生了顯著地變化.為了量化可壓縮效應(yīng)，圖11給出了文獻(xiàn)[3]提出的對(duì)流馬赫數(shù)Mc和Slessor等[18]提出衡量混合層可壓縮效應(yīng)的參數(shù)Πc沿流道方向的分布.Mc和Πc的定義式分別為：

式中：U為氣體速度；c為氣體中的聲速；γ為氣體比熱比.從圖11中可以看出，在噴射截面x=175 mm處，氫氣射流(Ma=2.3)與空氣來流(Ma=2.0)的對(duì)流馬赫數(shù)約為1.3，Slessor數(shù)約為3.75，說明兩股超聲速氣流所形成的混合層在開始階段具有很強(qiáng)的可壓縮性，經(jīng)過很短的距離后，對(duì)流馬赫數(shù)與Slessor數(shù)迅速降低，當(dāng)初始生長區(qū)結(jié)束后，Mc降為0.3，降幅約為77%，Πc降為0.4，降幅約為88%，這都表明混合層的可壓縮效應(yīng)顯著降低，混合層的生長速率大幅度提高.結(jié)合圖3中的流場(chǎng)結(jié)構(gòu)可以看出，由于流道擴(kuò)張?jiān)谥О迥┒诵纬膳蛎洸ㄏ?，氣流方向偏離流道方向，受到流道幾何構(gòu)型和兩股氣流所形成的氣動(dòng)壓縮面的約束，在膨脹波下游會(huì)形成再附激波使得氣流的整體流動(dòng)方向與流道方向一致，支板末端的這種膨脹波/激波結(jié)構(gòu)會(huì)使得混合層的可壓縮效應(yīng)大幅度減小，從而改變混合層的生長特性.這與通常所說的高對(duì)流馬赫數(shù)混合層的生長速率較小并不矛盾，只是與自由空間相比，受限空間內(nèi)混合層的對(duì)流馬赫數(shù)成為受復(fù)雜波系結(jié)構(gòu)影響的空間函數(shù)，初始對(duì)流馬赫數(shù)并不能全面描述混合層的生長特性.

圖11 Mc和Πc沿流道的分布

需要強(qiáng)調(diào)的是，雖然受限空間內(nèi)的波系結(jié)構(gòu)一方面能夠降低混合層的可壓縮性，整體上增大混合層的生長速率，另一方面激勵(lì)K-H不穩(wěn)定波的發(fā)展，強(qiáng)化局部湍流強(qiáng)度，獲得渦量增益，使得局部混合層的生長速率進(jìn)一步增大，從而促進(jìn)湍流摻混，提高混合效率，但沖壓流道中的復(fù)雜波系結(jié)構(gòu)不可避免地會(huì)引起較大的流動(dòng)損失，降低發(fā)動(dòng)機(jī)的整體性能.以噴射截面x=175 mm處的總壓為參考，圖12給出了總壓損失沿流道方向的分布曲線.

圖12 沿流道方向的總壓損失分布曲線

從圖12中可以看出沿著流道方向總壓損失不斷增大，在流道上游，激波強(qiáng)度較大，引起較大的總壓損失，總壓損失的上升速率較快，隨著流動(dòng)向下游的發(fā)展以及波系向下游的傳播，激波強(qiáng)度逐漸減小，總壓損失的上升速率趨于緩慢，在混合層生長階段，總壓損失接近9%.因此在實(shí)際的發(fā)動(dòng)機(jī)設(shè)計(jì)過程中，需要對(duì)激波提高的混合效率和引起的總壓損失進(jìn)行綜合考慮，實(shí)現(xiàn)優(yōu)化設(shè)計(jì).

3 結(jié) 論

1)受限空間內(nèi)混合層的對(duì)流馬赫數(shù)是受復(fù)雜波系結(jié)構(gòu)影響的空間函數(shù)，支板末端所形成的膨脹波/激波結(jié)構(gòu)會(huì)使對(duì)流馬赫數(shù)迅速降低，顯著減弱可壓縮效應(yīng)，從而提高混合層的生長速率，使得混合層經(jīng)過較短距離后進(jìn)入快速生長區(qū).

2)激波強(qiáng)度與混合層內(nèi)外兩側(cè)氣流的參數(shù)梯度共同影響著兩者之間的相互作用，由于激波的壓縮效應(yīng)，在作用點(diǎn)附近混合層的厚度會(huì)略微減小，但在作用點(diǎn)下游，激波會(huì)強(qiáng)化湍流，獲得渦量增益，促進(jìn)混合層的生長.

3)受限空間內(nèi)的復(fù)雜波系結(jié)構(gòu)在促進(jìn)湍流摻混，提高混合效率方面具有積極的作用，但同時(shí)會(huì)引起較大的總壓損失，在實(shí)際發(fā)動(dòng)機(jī)設(shè)計(jì)過程中需要對(duì)其進(jìn)行綜合考慮，實(shí)現(xiàn)優(yōu)化設(shè)計(jì).

[1] KUMAR A, BUSHNELL D H, HUSSAINI M Y. Mixing augmentation techniques for hypervelocity scramjets[J]. Journal of Propulsion and Power, 1989, 5(5): 514-522.DOI:10.2514/3.23184.

[2] 張漫.RBCC燃燒室湍流噴霧燃燒及其穩(wěn)定機(jī)制研究[D]. 西安: 西北工業(yè)大學(xué), 2010.

[3] PAPAMOSCHOU D, ROSHKO A.The compressible turbulent shear layer: an experimental study[J]. Journal of Fluid Mechanics, 1988, 197: 453-477. DOI:10.1017/S0022112088003325.

[4] FUREBY C, CHAPUIS M, FEDINA E, et al. CFD analysis of the HyShotⅡscramjet combustor[J]. Proceeding of the Combustion Institute, 2011, 33(2): 2399-2405.DOI:10.1016/j.proci.2010.07.055.

[5] GUTMARK E J, SCHADOW K C, YU K H. Mixing enhancement in supersonic free shear flows[J]. Annual Review of Fluid Mechanics, 1995, 27: 375-417.DOI:10.1146/annurev.fl.27.010195.002111.

[6] LI Qibing, FU Song. Numerical simulation of high-speed planar mixing layer[J]. Computers & Fluids, 2003, 32(10):1357-1377.DOI:10.1016/S0045-7930(02)00104-4.

[7] SANDHAM N D, REYNOLDS W C. Three-dimensional simulations of large eddies in the compressible mixing layer[J]. Journal of Fluid Mechanics, 1991, 224: 133-158.DOI:10.1017/S0022112091001684.

[8] VREMAN A W, SANDHAM N D, LUO K H. Compressible mixing layer growth rate and turbulence characteristics[J]. Journal of Fluid Mechanics, 1996, 320(1): 235-258.DOI:10.1017/S0022112096007525.

[9] MARBLE F E, HENDRICKS G J, ZUKOSKI E E. Progress toward shock enhancement of supersonic combustion process[C]//AIAA/SAE/ASME/ASEE 23rd Joint Propulsion Conference. San Diego, California: AIAA, 1987: 932-950.DOI: 10.2514/6.1987-1880.

[10]LU P J, WU K C. On the shock wave enhancement of confined supersonic mixing flows[J].Physics of Fluids A, 1991, 3(12): 3046-3062.DOI:10.1063/1.857849.

[11]HUH H, DRISCOLL J F. Shock-wave-enhancement of the mixing and the stability limits of supersonic hydrogen-air jet flames[J]. PSymposium (International) on Combustion, 1996, 26(2): 2933-2939.DOI:10.1016/S0082-0784(96)80135-4.

[12]ZHANG Yunlong, WANG Bing, ZHANG Huiqiang, et al. Mixing enhancement of compressible planar mixing layer impinged by oblique shock waves[J]. Journal of Propulsion and Power, 2015, 31(1): 156-169.DOI: 10.2514/1.B35423.

[13]GERLINGER P, BRUGGEMANN D. Numerical investigation of hydrogen strut injections into supersonic airflows[J]. Journal of Propulsion and Power, 2000, 16(1): 22-28.DOI: 10.2514/2.5559.

[14]張兆順, 崔桂香, 許春曉. 湍流大渦數(shù)值模擬的理論和應(yīng)用[M]. 北京: 清華大學(xué)出版社, 2008.

[15]王振國, 孫明波. 超聲速湍流流動(dòng)、燃燒的建模與大渦模擬[M]. 北京: 科學(xué)出版社, 2013.

[16]MENON S, YEUNG P K, KIM W W. Effect of subgrid models on the computed interscale energy transfer in isotropic turbulence[J]. Computers & Fluids, 1996, 25(2): 165-180.DOI:10.1016/0045-7930(95)00036-4.

[17]KIM J H, YOON Y, JEUNG I S. Numerical study of mixing enhancement by shock waves in model scramjet engine[J]. AIAA Journal, 2003, 41(6): 1074-1080.DOI: 10.2514/2.2047.

[18]SLESSOR M D, ZHUANG M, DIMOTAKIS P E. Turbulent shear-layer mixing: growth-rate compressibility scaling [J]. Journal of Fluid Mechanics, 2000, 414(1): 35-45.DOI:10.1017/S0022112099006977.

Studyonthesupersonicmixinglayergrowthinconfinedspace

WEI Xianggeng1,2, CAO Donggang2, QIN Fei2, WU Jiping1

(1.Science and Technology on Scramjet Laboratory(National University of Defense Technology), Changsha 410073, China; 2.Science and Technology on Combustion, Internal Flow and Thermal-Structure Laboratory (Northwestern Polytechnical University), Xi’an, 710072, China)

The supersonic mixing layer formed by a planar thin hydrogen jet at 2.3Maand a 2.0Masurrounding airflow in a scramjet engine model is studied in order to investigate the growth characteristics with consideration of shocks. Flow field structure and features of the supersonic mixing layer are achieved by using large eddy simulation method with the OpenFOAM software. Reasonable agreements are obtained between calculation and experiment in terms of flow field structure and component distribution. The component concentration, the thickness, the compressibility effect and the total pressure loss are analyzed. Results show that four developing regions can be observed for the growth of the mixing layer. The expansion-fan/shock-wave pattern at the injector exit makes the convective mach number decrease dramatically, leading to a reduction in compressibility effects and a contribution to the development of the mixing layer. The interaction of shock/mixing layer results in local amplification of turbulence and gain of vorticity, which is beneficial to the supersonic mixing. However, the increasement in total pressure loss is unavoidable in the presence of shocks because they can bring performance losses of the scramjet. Thus a tradeoff between the enhanced mixing efficiency and the decreased total pressure recovery should be considered in the scramjet optimization design.

supersonic mixing layer; shock wave; shockwave series; confined space; compressibility; large eddy simulation

10.11918/j.issn.0367-6234.201607066

V431

A

0367-6234(2017)10-0072-06

2016-07-16

高超聲速?zèng)_壓發(fā)動(dòng)機(jī)技術(shù)重點(diǎn)實(shí)驗(yàn)室開放基金(20110302021)

魏祥庚(1979—)，男，博士,副教授

魏祥庚，realysnow@nwpu.edu.cn

(編輯張 紅)