基于進氣道總壓畸變試驗的插板數值模擬研究

朱愛迪，鐘易成

(南京航空航天大學 能源與動力學院, 江蘇 南京 210016)

?

基于進氣道總壓畸變試驗的插板數值模擬研究

朱愛迪，鐘易成

(南京航空航天大學 能源與動力學院, 江蘇 南京 210016)

摘要：通過數值計算方法研究了進氣道插板式總壓畸變試驗技術。在已有試驗數據的基礎上，應用數值仿真的手段模擬了插板擾流后的流場特性，分析了插板后不同位置的流場結構和畸變度。通過和試驗數據對比，得出計算結果在AIP截面和試驗結果在趨勢上吻合良好，數值上則存在較大差異。在插板后異于AIP截面，存在能滿足周向總壓畸變指數與試驗值相當的截面，從而可通過插值處理，計算出不同進口馬赫數對應的計算截面，對試驗進行初步預估。

關鍵詞：數值計算；插板；總壓畸變

0引言

隨著飛機飛行速度、高度的不斷提高，機動性的不斷增加，導彈武器的使用，進氣道出口氣流產生畸變，嚴重影響了發動機的穩定工作。目前國內主要通過插板畸變試驗來確定所允許的發動機進口流場畸變程度[1]。插板試驗技術最初是由俄羅斯發展起來的，它是通過在發動機前放置一定高度的擾流器，產生動態+穩態的畸變流場，以得到與進口流場一致的畸變圖譜與指數，來考核發動機的穩定性[2]。

插板畸變試驗屬于風險比較高的試驗科目，如果對插板深度與進氣畸變程度沒有清晰的認識，在試驗中可能會對發動機造成不必要的損害。而且試驗的周期長成本高，不能對設計加工起到預見和指導作用?？紤]到數值模擬技術經濟而且快速的特點，可以運用數值方法對試驗過程進行模擬，預估試驗方案[3]。

在已有總壓畸變試驗數據基礎上，運用Fluent軟件進行數值模擬，獲取畸變流場結構以及相關畸變評定參數的變化規律，并和試驗結果進行了對比。

1數值計算

數值模擬的對象簡圖如圖1所示。

圖1　試驗裝置

試驗裝置為一內徑D=206mm的足夠長的等直段，前部安裝雙扭線型進氣頭部，以保持均勻來流。頭部與試驗段采用螺栓聯接，聯接處即為放置插板位置，便于更換插板。插板后2D處為測量截面，總壓測耙就安裝在這個截面上。插板有效厚度為12mm，其幾何特性由相對堵塞深度L=H/D來描述，如圖2所示。

圖2　擾流插板幾何特性簡圖

由于該幾何模型結構相對簡單，因此在生成網格時在整個計算域內均選擇結構化網格。在插板附近及近壁面區域對網格進行加密，使得Y+符合計算模型要求，其他通流部分則采用較粗的網格。計算網格也進行了網格無關性驗證，對不同網格密度的模型進行了試算，結果表明當網格總數超過90萬后，AIP截面氣動參數無明顯變化，因此文中計算模型網格數都保持在90萬左右。

數值計算采用Fluent軟件，基于雷諾平均的N-S方程組，有限體積法離散控制方程，它全面的考慮了流體的粘性、熱傳導和可壓縮性等的影響。層流粘性系數由Sutherland公式求出，湍流模型選取SSTk-ω模型。進口采用壓力遠場邊界，出口采用壓力出口邊界，其余均為無滑移絕熱壁面邊界。

2計算結果與分析

文中對L=30%和L=40%的插板分別計算了4個試

驗狀態點。通過調節出口背壓，在雙扭線出口獲取所需要的來流馬赫數M1。在AIP截面后氣流方向沿程截取若干個截面，觀察流場變化情況。由于計算狀態比較多，只選取一個狀態點的流場結構進行分析。

圖3和圖4分別為L=30%，M1=0.282時進氣系統的二維流線圖和沿氣流方向不同截面的總壓畸變圖譜。從流線圖可以看出，插板后產生了回流區域，由于氣流的摻混造成了總壓的損失。從圖4可以看出，總壓畸變圖譜基本分為明顯的三個區域，其中2個壓力梯度較小，分別為高壓區和低壓區，一個壓力梯度較大，為摻混區。隨著氣流沿程的不斷摻混，在達到插板后某截面時，總壓畸變圖譜形狀以及等值線的大小基本不再變化。

圖3　L=30%，M1=0.282時進氣系統流線圖

圖4　L=30%，M1=0.282時插板后沿程不同截面總壓畸變圖譜

對數值計算得到的流場數據進行了分析，并同試驗數據進行了對比，以考察數值計算結果同試驗測量值之間的區別。

從圖5可以看出，對于L=30%和L=40%的插板，隨著M1的增大，不論是試驗值還是計算值，總壓恢復系數呈單調下降趨勢；圖6則表明低壓區角度隨M1的增大基本不變，只與插板高度H有關。圖7表明隨著M1的增大，不論是試驗值還是計算值，周向總壓畸變指數呈單調上升趨勢。

從圖5-圖7可以看出，在AIP截面，數值計算結果與試驗數據在趨勢上吻合較好，可以滿足定性分析，但是定量則存在較大差異，尤其當插板高度變大、M1變大時。

圖5　總壓恢復系數隨M1變化情況

圖6　低壓區角度隨M1變化情況

圖7　Δσ0隨M1變化情況

數值模擬對板后氣流高壓區和低壓區的摻混能力不足，導致在AIP界面，數值計算數據普遍大于試驗數據。氣流在板后經過沿程摻混后會逐漸趨向均勻，周向總壓不均勻度也會逐漸較小。文中對沿程不同截面數據進行了統計分析。

圖8-圖9分別為插板高度L=30%和L=40%時沿程不同截面低壓區角度隨M1變化曲線圖。從圖中可以看出，隨著M1的增大，同一截面對應的低壓區角度基本不變；對應于同一M1，沿程不同截面低壓區角度變化不大。

圖8　L=30%低壓區角度沿程曲線圖

圖9　L=40%低壓區角度沿程曲線圖

圖10-圖11為插板高度L=30%和L=40%時Δσ0在不同進口M1時沿程變化曲線圖。從圖中可以看出，Δσ0沿程單調降低，而且趨勢逐漸變緩，這主要是高低壓氣流之間由于不斷摻混而趨于均勻的緣故。這種趨勢與插板高度和氣流進口馬赫數無關穩態計算時，衡量畸變程度主

圖10　L=30%Δσ0沿程曲線圖

圖11　L=40%Δσ0沿程曲線圖

要考慮低壓區范圍和周向總壓畸變指數兩個參數。低壓區范圍沿氣流方向基本無變化，鑒于周向總壓畸變指數沿程單調降低，考慮對沿程不同截面畸變指數進行曲線擬合處理，將試驗值進行插值，找出試驗值所對應的計算截面。結果如圖12所示。

從圖中可以看出，隨著M1的增加，Δσ0試驗值所對應的截面也在不斷后移；隨著插板高度H的增加，Δσ0試驗值所對應的截面也在后移。

當試驗要求不同的進口M1時，可以根據該圖進行插值，求得計算對應截面值，然后進行CFD計算，找出該截面的各種氣動參數，對試驗進行初步預估。

4結語

圖12　Δσ0試驗值所對應的計算截面試驗裝置

運用CFD技術數值模擬了不同高度插板的畸變試驗的進氣流場，獲得了畸變流場的氣動特性、總壓圖譜。在AIP截面，利用數值仿真數據進行精確的定量分析還難以實現。文中考慮利用插板后其它截面數值仿真數據和試驗AIP截面數據進行對比。對比結果表明，在插板后異于AIP截面，存在能滿足Δσ0與試驗值相當的截面，從而可通過插值處理，計算出不同進口M1對應的計算截面，對試驗進行預估。

參考文獻:

[1] 江勇. 航空發動機插板式進氣畸變實驗技術研究[D]. 西安:空軍工程大學, 2005.

[2] 劉大響, 葉培梁. 俄羅斯的發動機進口流場畸變評定指南[J]. 燃氣渦輪試驗與研究, 1994，(3): 1-10.

[3] 薛海波. 插板后畸變流場的數值模擬及參數分析[D]. 南京：南京航空航天大學, 2009.

Research on Flashboard Numerical Simulation Based on Inlet

Total Pressure Distortion Experiment

ZHU Ai-di, ZHONG Yi-cheng

(College of Energy and Power, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China)

Abstract:The numerical calculated method is used to research on the technology of the total pressure distortion with flashboard. Based on experimental data, the method of numerical simulation is used to simulate the flow field. Then flow field structure and distortion degree of different positions behind the flashboard are analyzed. Comparing the numerical result with experimental data, they trend to great agreement, but a great difference is made in the mumerical value. Behind AIP, a section -where Δσ0is the same with experimental data. The interposition method can be used to calculate the section corresponding to different M1, and complete the preliminary estimation of the experiment.

Keywords:numerical simulation; flashboard; total pressure distortion

基金項目：國家863課題(2011AA110506)；國家科技支撐計劃項目(2011BAG01B05)

收稿日期：2014-11-29

中圖分類號：TH123+.1

文獻標志碼：A

文章編號：1671-5276(2015)03-0143-03

作者簡介：朱愛迪(1989-)，男，湖北天門人，碩士研究生，研究方向為流場數值模擬。