進氣道流向伸縮對燃燒室入口總壓恢復的影響

丁浩 田立豐 郭美琦 桂斌

摘要：采用進氣道流向伸縮的方法，對某高超聲速飛行器進氣道彎道及隔離段流場開展了數值模擬計算，獲取了不同伸縮量對應的總壓云圖、馬赫數云圖和燃燒室入口總壓恢復。通過對比發現，在Ma=5的巡航速度下，隨著進氣道向下游的伸展，燃燒室入口總壓恢復系數呈現先升高后降低的趨勢，當伸縮量為30mm時，燃燒室入口總壓恢復系數可達到0.2484，比零偏移時提升了0.4%。通過對特定高超飛行器進氣道采用流向伸縮的方法得到了最優總壓恢復系數，為其他研究者在飛行器進氣道優化設計方向提供了一種新型且有效的方法。

關鍵詞：進氣道；高超聲速；數值模擬；燃燒室；總壓恢復

中圖分類號：V235.21文獻標識碼：ADOI：10.19452/j.issn1007-5453.2021.08.004

基金項目：裝備預研重點實驗室基金（6142703180211）

近年來，伴隨著航空航天技術的快速發展，世界各國加大了對高超聲速飛行器研制的投入力度。進氣道是高超聲速飛行器推進系統的重要部件，對超燃沖壓發動機的性能起著關鍵作用，內收縮進氣道憑借其較高的壓縮效率、較強的流量捕獲能力等優點，受到各研究者的重點關注[1]。

現代高性能飛行器大多采用S形進氣道來提高飛行器的隱身性能，但由于進氣道內部曲率較大，使得邊界層容易分離[2]。同時超聲速氣流與進氣道壁面相互作用會產生復雜的激波結構，氣流經過激波后總壓損失較大，減弱了發動機的推力。飛行器在真實高超聲速飛行過程中，實際來流往往不是均勻來流，因此采用試驗方法開展研究難度較大、成本較高。隨著計算流體力學（CFD）技術的快速發展，采用CFD方法研究進氣道流場成為一種經濟且有效的途徑[3-5]。

為了提高飛行器氣動性能，研究者針對進氣道以及隔離段做了大量優化方案研究。其中，王驥飛[6]通過改變壁面型線、優化進氣道唇口平面形狀、引入型面漸變技術和邊界層修正技術改善了進氣道的氣動性能。李蔚霆[7]、劉蕾等[8]運用編制程序和數值仿真方法，分析了進氣道的長寬比和中心角等幾何參數對其氣動性能的影響。甘文彪等[9]采用改進SST湍流模型和尺度自適應模擬方法對原始設計的S形進氣道進行了優化設計和分析，結果顯示與原始設計相比總壓恢復提高1.1%。劉雷等[10]采用Isight優化軟件對S形進氣道擴壓段進行優化設計，結果顯示優化后的進氣道出口截面總壓恢復、周向總壓畸變、旋流畸變都有所改善。王昌盛等[11]對高超聲速軸對稱進氣道型面經過多目標優化設計后，在設計點下進氣道總壓恢復提升2.63%。

本文在上述研究的基礎上提出一種進氣道流向伸縮的優化方法，旨在對已成形的高超聲速飛行器進氣道及隔離段進一步優化。由于S形進氣道內部曲率大，高超聲速氣流通過進氣道和隔離段時形成的激波結構更加復雜，進氣道和隔離段的相對長度會對激波結構的分布有微弱影響，因此通過調整進氣道與隔離段交界面的位置來改變激波結構的分布，從而提高燃燒室入口總壓恢復成為一種切實可行的辦法。本文采用Fluent軟件對某高超聲速飛行器整體內外流進行數值模擬，重點圍繞進氣道彎道及隔離段總壓恢復進行研究與分析。保持進氣道S彎道與隔離段總長度不變，調整兩者交界面的位置來探究這種偏移對燃燒室入口處總壓恢復的影響，探尋燃燒室入口處最高總壓恢復系數所對應的交界面位置，從而實現在原型基礎上對該飛行器進一步優化的目的。

1數值模擬

1.1氣動模型

本文的研究對象為某高超聲速飛行器的進氣道彎道和隔離段整體。圖1為原型進氣道彎道及隔離段對稱面外形，整體長度為1.45m，其中左側進氣道彎道0.85m，右側隔離段0.6m，進氣道彎道與隔離段連接處相切，以保證交界面處平滑過渡。

圖2為進氣道流向伸縮示意圖，本文通過流向伸縮來調整進氣道與隔離段交界面的流向位置，采用數值模擬方法研究這種調整給燃燒室入口的總壓恢復系數帶來的影響，并探尋最佳的調整狀態，以使燃燒室入口總壓恢復系數達到最高。

1.2網格劃分

本文采用結構網格劃分高超聲速飛行器內外流場，其中進氣道和隔離段總網格數為37.5萬個，同時為了使計算結果更為精準，在內流場的縱向上使用Power Law上下對稱加密，壁面第一層網格高度為0.5mm，增長因子為1.5。圖3給出了飛行器內外流整體網格、進氣道與隔離段網格以及局部放大網格示意圖。

1.3計算方法

本文運用Fluent軟件求解三維可壓縮Navier-Stokes方程；耦合能量方程并結合基于密度的求解器對高超聲速飛行器內外流場進行定常數值模擬；選用Roe-FDS迎風格式；湍流模型采用SST k-omega兩方程模型，該模型對自由剪切湍流、邊界層湍流和適度的分離湍流都有較高的計算精度。

選用理想氣體模型，導熱系數取0.0242W/ （m·K），黏性系數取1.7894×10-5kg/（m·s）；計算域進口邊界條件為壓力遠場來流，來流馬赫數Ma=5，來流靜壓為2971.75Pa；進氣道及隔離段計算域邊界采用無滑移絕熱壁面；計算以來流參數初始化，計算收斂準則為：連續方程、動量方程、能量方程、k-omega方程殘差，以及進出口流量相對誤差至少下降三個數量級，且進氣道和隔離段沿程截面流量穩定。

2結果與分析

2.1流場結構分析

通過對該高超聲速飛行器原型及進氣道流向伸縮后的模型內外流場數值模擬，并在模擬結果中選取了0mm、28mm、32mm和36mm這4個偏移量下進氣道及隔離段對稱面流場的馬赫數云圖和總壓云圖進行對比。圖4、圖5分別為4個偏移量時所對應的馬赫數云圖與總壓云圖。

從圖4中4個不同偏移量下的對稱面馬赫數云圖可以看出，超聲速氣流經進氣道擴張段進一步提速，在S彎道中與壁面相互作用生成結構復雜的激波結構，氣流經過這些激波時總壓也會發生相應的變化，這在圖5的總壓云圖中有相應的體現。

對比圖5中4個偏移量下的總壓云圖，結果發現不同偏移量下的總壓云圖變化并不明顯，因此需要更精確的折線圖才能得到偏移量對總壓的影響。

2.2偏移量對燃燒室入口總壓恢復的影響

圖6、圖7分別給出了進氣道S彎道入口處、燃燒室入口處總壓恢復隨交界面偏移量的變化曲線。對比兩圖發現，高超聲速氣流在進氣道S彎道入口處的總壓恢復系數在0.43～ 0.44之間，當氣流通過進氣道S彎道及隔離段后在燃燒室入口處總壓恢復系數降低到0.247左右，這是由于超聲速氣流經過S彎道及隔離段時激波阻力和壁面摩擦阻力的作用會使氣流減壓增速，從而導致總壓恢復系數大幅降低。

由圖7可以看出，燃燒室入口的總壓恢復系數隨交界面偏移總體呈現先增后降的趨勢，在偏移量為30mm左右時總壓恢復系數達到峰值0.2484，相比零偏移時提升了0.4%，說明調整進氣道與隔離段的交界面是一種優化飛行器氣動性能的有效方法。

3結束語

本文通過對某高超聲速飛行器氣動系統進氣道S彎道及隔離段進行數值模擬，在Ma=5的條件下，分別得到了不同伸縮量對應的總壓云圖、馬赫數云圖以及進氣道S彎道入口和燃燒室入口處總壓恢復系數隨進氣道流向伸縮量的變化曲線，得出如下結論：該型高超聲速飛行器燃燒室入口處總壓恢復系數隨進氣道流向伸縮量的變化總體呈現出先增后減的趨勢，在偏移量為30mm左右達到峰值0.2484，比零偏移時提升了0.4%。說明進氣道流向伸縮是一種提高燃燒室入口總壓恢復系數的有效手段。為其他研究者在飛行器進氣道優化設計方向提供了一種新型且有效的參考方法。

參考文獻

[1]代春良，孫波，梁曉揚，等.真實氣體效應下高馬赫數內轉進氣道特性研究[J].推進技術， 2020，41（7）：1473-1483. Dai Chunliang， Sun Bo， Liang Xiaoyang， et al. Study on characteristics of high Mach number inward turning inlet under real gas effect[J]. Journal of Propulsion Technology， 2020， 41（7）：1473-1483.（in Chinese）

[2]田曉平，潘鵬飛，李密. S形進氣道流動控制數值模擬研究[J].航空科學技術，2017，28（5）：42-46. Tian Xiaoping， Pan Pengfei， Li Mi. Number simulation research on S-inlet flow control[J]. Aeronautical Science & Technology， 2017， 28（5）： 42-46.（in Chinese）

[3]Nishizawa U，Kameda M，Watanabe Y，et al. Computational simulation of shock oscillation around a supersonic air-intake[C]//36thAIAAFluid Dynamics Conference and Exhibit，2006.

[4]Bur R，Benay R，Galli A，et al. Experimental and numerical study of forced shock-wave oscillations in a transonic channel[J].Aerospace Science and Technology，2006，10（4）：265-278.

[5]黃輝先，李斌，莊選，等.不同來流條件下高超聲速進氣道/隔離段的數值模擬[J].湘潭大學自然科學學報，2013，35（1）：113-117. Huang Huixian， Li Bin， Zhuang Xuan， et al. Numerical simulation of ramjet inlet/isolator under different airflow condition[J]. Natural Science Journal of Xiangtan University， 2013， 35（1）： 113-117.（in Chinese）

[6]王驥飛.高超聲速飛行器氣動外形一體化設計方法研究[D].西安：西北工業大學， 2018. Wang Jifei. Research on integration design methodology of aerodynamicshapeforhypersonicaircrafts[D].Xian： Northwestern Polytechnical University， 2018.（in Chinese）

[7]李蔚霆.寬速域高效內收縮式進氣道設計技術研究[D].南京：南京航空航天大學，2018. Li Weiting. Design methodology of performance inlet with wide operation range[D]. Nanjing： Nanjing University of Aeronautics andAstronautics， 2018.（in Chinese）

[8]劉蕾，郭雪巖.喉道肩點圓弧半徑對高超聲速進氣道性能的影響[J].能源工程， 2017（6）：20-27. Liu Lei，Guo Xueyan. Impact of radius of isentropic expansion arc surface on the hypersonic inlet[J]. Energy Engineering， 2017（6）： 20-27.（in Chinese）

[9]甘文彪，周洲，許曉平，等. S彎進氣道優化設計及分析[J].推進技術，2014，35（10）：1317-1324. Gan Wenbiao， Zhou Zhou， Xu Xiaoping， et al. S-duct inlet optimal design and analysis[J]. Journal of Propulsion Technology， 2014，35（10）：1317-1324.（in Chinese）

[10]劉雷，宋彥萍，陳煥龍，等.半埋入式S彎進氣道設計及其優化[J].推進技術， 2014，35（10）：1303-1309. Liu Lei， Song Yanping， Chen Huanlong， et al. Design and optimization of half flush-mounted S-shaped inlet[J]. Journal of Propulsion Technology，2014，35（10）：1303-1309.（in Chinese）

[11]王昌盛，額日其太，丁文豪.高超聲速軸對稱進氣道多目標優化設計[J].航空動力學報， 2020，35（7）：1392-1401. Wang Changsheng， Eriqitai， Ding Wenhao. Multi-objective optimization design of hypersonic axisymmetric inlet[J]. Journal of Aerospace Power， 2020， 35（7）： 1392-1401.（in Chinese）

Impact of Air Inlet Flow Expansion on Total Pressure Recovery of Combustor Inlet

Ding Hao，Tian Lifeng，Guo Meiqi，Gui Bin

School of Aeronautics and Astronautics，Sun Yat-sen University，Guangzhou 510275，China

Abstract： The flow field at the inlet elbow and isolator of a hypersonic vehicle was numerically simulated by using the inlet flow direction expansion method， and the total pressure contour， Mach number contour and the total pressure recovery at the inlet of the combustor corresponding to different expansion were obtained. By comparison， it is found that at the cruising speed of Ma=5， the total pressure recovery coefficient at the combustor inlet increases firstly and then decreases with the expansion of the inlet to the downstream. When the expansion is 30mm， the total pressure recovery coefficient at the inlet of the combustor can reach 0.2484， which is 0.4% higher than that of the zero offset. The optimal total pressure recovery coefficient is obtained by using the method of flow direction expansion for a particular hypersonic vehicle inlet， which provides a new and effective method for other researchers in the direction of aircraft inlet optimization design.

Key Words： air inlet； hypersonic； numerical simulation； combustor； total pressure recovery