王鍵靈，徐讓書,1b，侯亞軍，王娟娟

(1.沈陽航空航天大學(xué) a.航空航天工程學(xué)部；b.遼寧省航空推進系統(tǒng)先進測試技術(shù)重點實驗室，沈陽 110136；2.中國航空工業(yè)集團公司中國航空工業(yè)研究院 航空發(fā)動機高空模擬航空科技重點實驗室，四川 江油，621703)

自由射流中飛機進氣道前方亞聲速流場數(shù)值仿真研究

王鍵靈1a，徐讓書1a,1b，侯亞軍1a，王娟娟2

(1.沈陽航空航天大學(xué) a.航空航天工程學(xué)部；b.遼寧省航空推進系統(tǒng)先進測試技術(shù)重點實驗室，沈陽 110136；2.中國航空工業(yè)集團公司中國航空工業(yè)研究院 航空發(fā)動機高空模擬航空科技重點實驗室，四川 江油，621703)

為了預(yù)測航空發(fā)動機高空模擬自由射流試驗中飛機進氣道-發(fā)動機組合體前方的亞聲速流場特性，優(yōu)化試驗艙的氣動設(shè)計，采用CFD方法在亞聲速自由射流和真實大氣飛行條件下對某戰(zhàn)斗機進氣道的外流場進行數(shù)值模擬。分析了進氣道對前方氣流壓縮作用與飛行馬赫數(shù)關(guān)系，對比和分析了分別在自由射流與真實大氣中進氣道前和進氣道入口處的馬赫數(shù)分布，確定了3種馬赫數(shù)下進氣道在自由射流中的最佳安裝位置。比較發(fā)現(xiàn)，亞聲速自由射流對真實高空大氣飛行進行模擬，可以得到馬赫數(shù)相似的流場。

自由射流；發(fā)動機高空模擬；進氣道；馬赫數(shù)；數(shù)值計算

自由射流裝置由于其可以模擬整個飛行包線環(huán)境這一特點，可以被用于航空發(fā)動機高空試驗臺中。早在1990年美國ASTF C-2自由射流裝置就已經(jīng)投入使用，我國的自由射流高空臺也在建設(shè)當(dāng)中[1-4]。

由于自由射流的核心區(qū)能夠保證流體的速度幾乎不變，所以可將進氣道放置于自由射流的核心區(qū)內(nèi)來模擬高空大氣飛行環(huán)境。但因為亞聲速進氣道和發(fā)動機工作情況會對遠前方來流存在擾動作用，這種擾動作用可以使來流減速增壓或者加速減壓，而改變進氣道前的流場情況。所以有必要找到合適的進氣道安裝位置來保證射流核心區(qū)馬赫數(shù)達到要求，并保證能夠?qū)⑦M氣道完全包裹在核心區(qū)內(nèi)[5-7]。本文針對上述問題進行了數(shù)值計算，研究了進氣道對前方流體擾動作用的距離與飛行馬赫數(shù)的關(guān)系，并和對應(yīng)的自由射流條件中的擾動情況進行了對比，確定了0.4、0.6和0.8三種飛行馬赫數(shù)下的進氣道的最優(yōu)安放位置，得到了這三種飛行條件下馬赫數(shù)分布相似的流場。

1 控制方程與計算模型

1.1 控制方程

質(zhì)量守恒定律、動量守恒定律和能量守恒定律是流體力學(xué)中最基本的三大定律。流體力學(xué)中的Navie-Stokes(N-S)方程就是根據(jù)以上3定律建立的。

質(zhì)量方程為：

式中，▽為哈密頓算子。

動量方程為：

能量方程為：

式中，δw為運輸功，δwf為摩擦損耗[8-11]。

本文的流體為可壓縮氣體，氣體的物性參數(shù)以多項式形式給出。

1.2 計算模型及邊界條件

本文采用三維模型，進氣道模型采用美國某戰(zhàn)斗機的腹下式S型進氣道外形，并做一定的簡化和修改，進口截面積為0.653 m2，進氣道部分長5.45 m，發(fā)動機及尾噴管部分長3.96 m。整個計算域為一個長38 m，寬20 m，高20 m的長方體。自由射流計算域模型如圖1所示，大氣環(huán)境的計算域則是將噴管去掉，并延長進氣道前方區(qū)域長度，其他條件不變。

圖1 計算域模型

在真實大氣環(huán)境的計算中計算域邊界均使用壓力遠場邊界條件，而在射流條件的計算中噴管進口采用壓力進口邊界條件，進口的總壓、靜壓、總溫值由流體力學(xué)等熵關(guān)系式求得[7]。計算域的其他邊界依然使用壓力遠場邊界條件。由于模擬的是10 000米高空的飛行環(huán)境，所以給定環(huán)境的靜壓值為26 500 Pa，環(huán)境溫度為223.3 K。噴管和進氣道壁面均為無滑移邊界條件。

1.3 自由射流裝置及其流動特性

自由射流氣體是通過一個二元矩形噴管產(chǎn)生的。噴管豎直方向的兩壁間距離面保持不變，另外兩壁間面距離沿Y方向逐漸減小。噴管進、出口截面面積比為2∶1，出口長1.8 m，寬1 m。氣流通過噴管加速后形成一個速度幾乎不變的核心區(qū)，沿著流動方向核心區(qū)寬度不斷減小，邊界層厚度不斷增大，直到轉(zhuǎn)捩面出現(xiàn)，核心區(qū)消失。本文正是利用這一核心區(qū)來模擬大氣環(huán)境，核心區(qū)形狀如圖2所示。

圖2 自由射流核心區(qū)

2 網(wǎng)格劃分

計算網(wǎng)格均采用ANSYS ICEM進行劃分。對于粘性流體的計算，結(jié)構(gòu)網(wǎng)格要比非結(jié)構(gòu)網(wǎng)格具有優(yōu)勢，但網(wǎng)格劃分的難度更高。所以計算域全部采用結(jié)構(gòu)化網(wǎng)格進行劃分。計算域模型網(wǎng)格總數(shù)量約為200萬，在Deteminant 2×2×2方法下檢查網(wǎng)格質(zhì)量，除少量邊界層網(wǎng)格外網(wǎng)格質(zhì)量均在0.5以上。外流場遠離進氣道前和進氣道的地方網(wǎng)格比較稀疏，在噴管出口以及進氣道部分進行加密處理[12-15]。在進氣道壁的內(nèi)外壁面各使用一層O網(wǎng)格來進行加密和提高網(wǎng)格質(zhì)量。近壁處理采用標(biāo)準壁面函數(shù)，y+滿足要求。進氣道前區(qū)域的網(wǎng)格劃分情況如圖3所示。

圖3 進氣道前區(qū)域網(wǎng)格

3 計算方法

本文采用ANSYS FLUENT 15.0作為計算軟件，使用穩(wěn)態(tài)求解器進行求解。因為Realizablek-ε模型能夠更精確地模擬平面和圓形射流的擴散速度，同時在旋轉(zhuǎn)流計算、帶方向壓強梯度的邊界層計算和分離流計算等問題中，計算結(jié)果更符合真實情況，故本文采用Realizablek-ε湍流模型。由于計算馬赫數(shù)大于0.3，氣體被認為是可壓縮的，所以采用基于密度的求解器，氣體的密度、比熱等參數(shù)以多項式的形式給出。由于考慮到計算時間的經(jīng)濟性，故采用一階迎風(fēng)格式作為離散格式。

3 計算結(jié)果分析

本文計算了0.4、0.6和0.8馬赫數(shù)的大氣環(huán)境和自由射流條件下進氣道前的流場情況，并在相同的馬赫數(shù)下對多個進氣道安裝位置進行了計算。模擬的飛行環(huán)境為10 000米高空，發(fā)動機流量均為30 kg/s，飛行攻角和側(cè)滑角均為0。

3.1 馬赫數(shù)分布比較

飛行馬赫數(shù)為0.6時，兩種條件下流場的馬赫數(shù)分布如圖4所示。

圖4 馬赫數(shù)分布

自由射流條件下，進氣道被完全包涵在自由射流核心區(qū)內(nèi)，在噴管出口處馬赫數(shù)達到0.6，之后沿著Y方向馬赫數(shù)逐漸降低，在進氣道進口處，兩種條件下的馬赫數(shù)均降低到0.4左右，進入進氣道后氣流馬赫數(shù)繼續(xù)降低。兩種條件下進氣道前及進氣道內(nèi)馬赫數(shù)的分布情況十分接近，變化趨勢基本一致。飛行馬赫數(shù)為0.4時，兩種條件下進氣道進口截面的馬赫數(shù)分布情況如圖5所示。

觀察上圖可以發(fā)現(xiàn)，兩種條件下進氣道進口截面馬赫數(shù)分形式布十分相近，變化趨勢一致，都是由上邊緣中心去向外逐漸降低，馬赫數(shù)等值面成圓環(huán)狀分布。本文計算的3種飛行馬赫數(shù)的進氣道進口截面馬赫數(shù)的面積加權(quán)平均值情況如表1所示。

圖5 進氣道截面進口馬赫數(shù)分布

表1 進氣道進口馬赫數(shù)面積加權(quán)平均值

3.2 馬赫數(shù)變化趨勢對比

兩種條件下進氣道中心對稱線上的馬赫數(shù)變化情況如圖6所示。

從圖6中可以發(fā)現(xiàn)，當(dāng)未受擾動截面處馬赫數(shù)大于進氣道進口馬赫數(shù)時，進氣道對前方流體有壓縮作用，即減速增壓作用。隨著飛行馬赫數(shù)提高，進氣道進口的平均馬赫數(shù)反而降低，遠前方未受干擾截面處馬赫數(shù)與進氣道進口馬赫數(shù)差值就越大。而這個馬赫數(shù)差值越大，所需要的壓縮距離就越長，對應(yīng)的在自由射流情況下噴管距進氣道進口的距離就越大。經(jīng)過比較相同飛行馬赫數(shù)下不同安放位置的馬赫數(shù)分布情況，最終確定馬赫數(shù)為0.4、0.6和0.8的進氣道按最佳安放距離分別為1.6 m、1.8 m和1.9 m。

圖6 氣道前方馬赫數(shù)變化

比較兩種環(huán)境下壓縮段的馬赫數(shù)變化情況可以發(fā)現(xiàn)，各對應(yīng)情況的馬赫數(shù)的變化趨勢一致，都是飛行馬赫數(shù)沿Y軸逐漸降低的過程。來流在進氣道進口截面的馬赫數(shù)十分接近，但自由射流條件下馬赫數(shù)的變化更為急促，在相對較短的距離內(nèi)就可以完成對來流的壓縮過程。例如0.6馬赫數(shù)下真實大氣環(huán)境的壓縮過程長度大致為3.5 m，而自由射流下僅需要1.6 m便可以完成。

3.3 噴管出口馬赫數(shù)與進氣道安放位置的關(guān)系

上述壓縮作用距離體現(xiàn)了在自由射流條件下進氣道的安放位置的重要性。若進氣道距離噴管太近，噴管出口馬赫數(shù)便達不到要求，發(fā)生堵塞。若距離太遠，射流核心區(qū)寬度會達不到要求。以0.6馬赫數(shù)為例，進氣道距噴管出口距離為1.2 m、1.4 m、1.6 m和1.8 m時，噴管出口水平中心線馬赫數(shù)分布情況如圖7所示。

從圖上可以看出當(dāng)安放距離小于1.8 m時，噴管出口馬赫數(shù)達不到試驗飛行馬赫數(shù)的要求。所以找到合適的進氣道安放位置對模擬是否成功有很重要的影響。進氣道安放距離隨飛行馬赫數(shù)的變化規(guī)律與大氣環(huán)境下的進氣道對來流壓縮作用距離隨飛行馬赫數(shù)變化的規(guī)律一致，即飛行馬赫數(shù)越大，進氣道安放距離就越遠。但想找到合適的進氣道安放位置，只能通過計算不同的安放距離來找出一個合適的值。

圖7 噴管出口馬赫數(shù)

4 結(jié)論

(1)自由射流條件下進氣道對前方流體馬赫數(shù)變化梯度大于真實大氣條件，使得自由射流條件的流管長度更短。

(2)自由射流條件下，進氣道安放位置必須滿足大于最小安放距離這一條件，否則射流核心區(qū)速度達不到要求值。經(jīng)過比較不同安裝距離最終確定0.4、0.6和0.8飛行馬赫數(shù)下的最佳安裝距離分別為1.6 m、1.8 m和1.9 m。

(3)通過真實大氣環(huán)境和自由射流環(huán)境流場馬赫數(shù)的比較，可以認為自由射流對真實大氣飛行環(huán)境下進氣道前流場馬赫數(shù)分布的模擬可行，能夠獲得進氣道進口流場馬赫數(shù)的相似，進口馬赫數(shù)的面積加權(quán)平均值的誤差均可控制在1.51%以下。

[1]杜鶴齡.航空發(fā)動機高空模擬[M].北京：國防工業(yè)出版社，2002:44-54.

[2]馮國泰,黃家驊,王松濤.航空發(fā)動機數(shù)值仿真試驗臺建立中的幾個關(guān)鍵技術(shù)問題的討論[J].航空動力學(xué)報,2002(10):17-40.

[3]Beale D K,Zelenak M.Development and validation of a free-jet technique for inlet-engine compatibility Testing[J].AIAA,1992(92):3921.AIAA 17th Aerospace Ground Testing Conference,Nashville:TN,1992.

[4]Beale D K,Kelly P G,Capt J E P.Lacasse.Subscale validation Of a free-jet inlet-engine test capability[J].AIAA,1993(93):2179.

[5]Dardis W,Mayhew E,Beale D.An acceptance process for the evaluation of inlet distortion[J].AIAA,1992(92):3918.AIAA 17th Aerospace Ground Testing Conference,Nashville:TN,1992.

[6]Mayward P V.C-2 subscale free-jet development and demonstration[J].AIAA,1993(93):2180.AIAA/SAE/ASME/ASEE 29th Joint Propulsion Conference and Exhibit,Monterey:CA,1993.

[7]Lynn C Sebourn,Frederick L.Shope.Research summary on the AEDC ASTF C-2 aeroacoustic resonance phenomenon[J].AIAA,2005(2932):1-12.11th AIAA/CEAS/Aeroacoustics Conference,Monterey:California,2005.

[8]王福軍.計算流體力學(xué)分析[M].北京：清華大學(xué)出版社,2004:7-9.

[9]劉月玲,喬渭陽,許開富.S彎進氣道內(nèi)流分離數(shù)值仿真.計算機仿真[J],2010(10):29-34.

[10]趙承慶,姜毅.氣體射流動力學(xué)[M].北京：北京理工大學(xué)出版社,1998:37-47.

[11]董志勇.射流力學(xué)[M].北京：科學(xué)出版社,2005:35-40.

[12]廉筱純,吳虎.航空發(fā)動機原理[M].西安:西北工業(yè)大學(xué)出版社,2001:29-36.

[13]劉振俠,郭東明,張麗芬,等.S形進氣道流場數(shù)值模擬[J].航空動力學(xué)報,2006(12):21-26.

[14]羅志煌,李軍,劉超,等.某S型進氣道總壓畸變的數(shù)值模擬研究[J].微計算機信息,2010(26):6-11.

[15]夏楊,李博,王海明.一種兩側(cè)布局的無隔道亞聲速進氣道流場特性[J].航空動力學(xué)報,2013(2):28-32.

(責(zé)任編輯：宋麗萍 英文審校：劉敬鈺)

Numericalcalculationofsubsonicflowfieldinfrontofaircraftinletductinfreejet

WANG Jian-ling1a，XU Rang-shu1a,1b，HOU Ya-jun1a，WANG Jun-jun2

(a.Faculty of Aerospace Engineering;b.Liaoning Key Laboratory of Advanced Measurement and Test Technology for Aviation Propulsion System,1.Shenyang Aerospace University,Shenyang 110136,China;2.Aviation Key Laboratory of Science and Technology and Aero-Engine Altitude Simulation,China Gas Turbine Establishment,Aviation Industry Corporation of China,Jiangyou 621703,China)

To predict the characteristics of subsonic flow field in front of aircraft inlet-engine combination in the free jet test of engine altitude simulation and to optimize the aerodynimic configuration of the test cell,the CFD method has been used to simulate the external flow of a fighter inlet in subsonic free-jet and in the real flight environment in the atmosphere.The relation of compression effect of the inlet to the inflow air and flight Mach number is analyzed.Optimized installation location of inlet has been found in three Mach number of subsonic free-jet.And we compared the Mach number distribution of the forward and entrance of inlet between the condition that in free-jet and in real atmospheric environment.As a result,it can obtain the flow field of similar Mach number that the simulating of real atmospheric flight environment with subsonic free-jet.

free jet;engine simulation of altitude;inlet duct;Mach number;numerical study

2014-06-26

王鍵靈(1988-)，男，遼寧撫順人，碩士研究生，主要研究方向：航空發(fā)動機數(shù)值仿真，E-mail:283410145@qq.com； 徐讓書(1962-)，男，浙江樂清人，教授，主要研究方向：航空發(fā)動機數(shù)值仿真，E-mail:xurangshu@yahoo.com。

2095-1248(2014)05-0039-05

V216.8

A

10.3969/j.issn.2095-1248.2014.05.008