， ，

(中國飛行試驗(yàn)研究院 發(fā)動(dòng)機(jī)所，陜西 西安 710089)

0 引言

航空發(fā)動(dòng)機(jī)的飛行推力確定是飛機(jī)與發(fā)動(dòng)機(jī)性能試飛的重要內(nèi)容。目前所常用的飛行推力確定方法都需事先已知被試發(fā)動(dòng)機(jī)的尾噴管特性[1]。獲取尾噴管特性的方法有：縮比模型吹風(fēng)試驗(yàn)，全尺寸部件試驗(yàn)、CFD數(shù)值計(jì)算等。隨著計(jì)算機(jī)技術(shù)和數(shù)值計(jì)算方法的快速發(fā)展，采用CFD計(jì)算手段獲得發(fā)動(dòng)機(jī)尾噴管在各種飛行條件下的特性已經(jīng)成為非常重要的方法。美國GE公司的CF34-10A渦扇發(fā)動(dòng)機(jī)和CFM公司的LEAP-X1C渦扇發(fā)動(dòng)機(jī)以及某型國產(chǎn)發(fā)動(dòng)機(jī)等在尾噴管特性研究時(shí)均采用了CFD數(shù)值計(jì)算的方法[2-6]。

我國某型機(jī)用小型分開排氣的渦扇發(fā)動(dòng)機(jī)的尾噴管為長外涵道形式(即內(nèi)外涵道出口共面)，為了計(jì)算該發(fā)動(dòng)機(jī)的標(biāo)準(zhǔn)凈推力，需已知其尾噴管的特性。文獻(xiàn)[7]中給出了渦扇發(fā)動(dòng)機(jī)的不同形式的尾噴管的特性曲線，但沒有介紹長外涵道形式的尾噴管特性。國內(nèi)的文獻(xiàn)[5-6]計(jì)算的均為分開排氣的短外涵道形式的尾噴管特性，而由于從長外涵道形式的噴管的結(jié)構(gòu)特點(diǎn)，需計(jì)算分析其內(nèi)外涵出口氣流的互相干涉作用，以確定內(nèi)外涵道特性的影響參數(shù)及變化規(guī)律，為標(biāo)準(zhǔn)凈推力計(jì)算提供數(shù)據(jù)支持。因此，本文以該型發(fā)動(dòng)機(jī)的尾噴管為研究對象，計(jì)算了該尾噴管在不同飛行馬赫數(shù)、不同工況下的三維粘性流場，分析了內(nèi)外涵道的流量特性與推力特性與飛行馬赫數(shù)、內(nèi)外涵進(jìn)口總壓比、噴管落壓比之間的關(guān)系。

1 計(jì)算方法

1.1 CFD計(jì)算方法

圖1 噴管結(jié)構(gòu)示意圖

研究的噴管的結(jié)構(gòu)示意圖如圖1，可以看出外涵道的出口位置與內(nèi)涵出口幾乎共面，明顯不同于分開排氣的民用大涵道比渦扇發(fā)動(dòng)機(jī)。本文利用ICEM軟件對計(jì)算模型進(jìn)行結(jié)構(gòu)化網(wǎng)格剖分，計(jì)算域?yàn)橹芟?/6區(qū)域。計(jì)算采用ANSYS CFX軟件，計(jì)算時(shí)，噴管出口為壓力遠(yuǎn)場，改變遠(yuǎn)場來流速度以模擬不同飛行馬赫數(shù)，改變內(nèi)外涵道進(jìn)口總壓，以模擬不同的工況和不同的內(nèi)外涵進(jìn)口總壓比。湍流模型為k-epsilon模型，對于噴管壁面的流動(dòng)模擬選擇CFX提供的可擴(kuò)展壁面函數(shù)方法(Scalable Wall Functions)。

1.2 尾噴管特性獲取方法

以內(nèi)涵道為例，其流量系數(shù)Cd1和推力系數(shù)Cfg1的定義如下：

(1)

式中，實(shí)際流量Wact1由CFD計(jì)算結(jié)果直接得到，實(shí)際總推力Fgact1：

Fgact1=Wact1Vaxial+(P81-P0)A81

(2)

其中：A81內(nèi)涵出口面積，P81為內(nèi)涵出口靜壓，其值由計(jì)算結(jié)果得到。

由于該尾噴管的內(nèi)、外涵道均為收斂式的，因此式(1)中理想流量Wideal1與理想總推力Fgideal1的定義如下：

(3)

(5)

(6)

式中，Pt1、Tt1為內(nèi)涵道進(jìn)口總壓總溫，P0為外界大氣壓力，γ=1.4，R=287.06(kJ/kg·k)。

根據(jù)相似原理分析，固定收斂尾噴管的流量特性和推力特性表征為外流馬赫數(shù)、尾噴管落壓比的函數(shù)[8]，而對于雙涵道噴管，其影響因素還應(yīng)包括內(nèi)涵與外涵進(jìn)口壓力之比，即：

流量系數(shù)：

(7)

推力系數(shù)：

(8)

因此，獲得噴管內(nèi)涵道的特性參數(shù)的計(jì)算過程如下：

1)根據(jù)內(nèi)涵道A81、Pt1、Tt1、P0，計(jì)算其理想的流量和理想的總推力；

2)根據(jù)CFD計(jì)算結(jié)果，計(jì)算實(shí)際流量和實(shí)際總推力；

外涵道的計(jì)算方法與內(nèi)涵道相同。

2 結(jié)果分析

本文計(jì)算了不同Ma，不同落壓比、不同內(nèi)外涵進(jìn)口壓力比工況下的內(nèi)外涵特性參數(shù)。下面對結(jié)果加以說明。

2.1 內(nèi)外涵壓力比對內(nèi)、外涵噴管特性的影響

當(dāng)內(nèi)外涵工況不同時(shí)，各自出口流場壓力值不同，在噴管出口附近，內(nèi)外涵氣流會(huì)互相干涉，進(jìn)而影響內(nèi)、外涵的特性。圖2為Ma=0、0.5，rpt=0.8、1.0、1.2時(shí)，外涵的流量系數(shù)和推力系數(shù)隨其落壓比的變化關(guān)系。

由圖2(a)、(c)可以看出，地面靜止條件下，外涵的流量系數(shù)與推力系數(shù)均隨其落壓比的增大而增大，達(dá)到阻塞狀態(tài)后，流量系數(shù)不再變化，推力系數(shù)略有增加。同時(shí)，不同的內(nèi)外涵壓力比下，外涵道的流量系數(shù)與推力系數(shù)不變，及內(nèi)涵工況變化不影響外涵的特性參數(shù)。由圖2(b)、(d)，當(dāng)有外流馬赫數(shù)時(shí)，亞臨界時(shí)的外涵道流量系數(shù)與推力系數(shù)隨著落壓比的減小而增大；當(dāng)外涵道處于臨界與超臨界狀態(tài)時(shí)，其流量系數(shù)與推力系數(shù)也不因內(nèi)涵工況的變化而不同，當(dāng)外涵道處于亞臨界狀態(tài)時(shí)，其流量系數(shù)與推力系數(shù)隨內(nèi)涵工況變化的變化量非常小，可以認(rèn)為，外涵道的特性參數(shù)不隨內(nèi)涵工況的變化而變化。

圖2 不同的飛行馬赫數(shù)與內(nèi)涵工況，外涵特性參數(shù)與外涵壓力比的關(guān)系圖

由圖3(a)、(c)可以看出，地面靜止條件下，內(nèi)涵道流量系數(shù)與推力系數(shù)隨內(nèi)涵壓力比的變化趨勢與外涵道相同，但外涵工況對內(nèi)涵的特性參數(shù)有影響。即當(dāng)內(nèi)涵噴管處于亞臨界狀態(tài)時(shí)，rpt越大，則內(nèi)涵的流量系數(shù)與推力系數(shù)越大。當(dāng)內(nèi)涵處于臨界與超臨界狀態(tài)時(shí)，內(nèi)涵的特性參數(shù)不隨rpt的變化而變化。由圖3(b)、(d)，Ma=0.3時(shí)，內(nèi)涵處于亞臨界狀態(tài)時(shí)，rpt越大，內(nèi)涵的流量系數(shù)與推力系數(shù)越大，當(dāng)內(nèi)涵處于臨界與超臨界狀態(tài)時(shí)，內(nèi)涵的流量系數(shù)與推力系數(shù)不隨rpt的變化而變化。

圖3 不同的飛行馬赫數(shù)與外涵工況，內(nèi)涵特性參數(shù)與內(nèi)涵壓力比的關(guān)系

2.2 不同飛行馬赫數(shù)對內(nèi)、外涵噴管特性的影響

飛行馬赫數(shù)不同，對噴管的特性會(huì)產(chǎn)生一定影響。這是因?yàn)閬砹魉俣葧?huì)與外涵出口氣流速度互相干涉， 同時(shí)，由于機(jī)體的存在，在機(jī)體到外涵出口附近會(huì)形成低壓區(qū)，進(jìn)而影響了外涵的落壓比。有關(guān)馬赫數(shù)對內(nèi)外涵道特性參數(shù)的影響見圖4和圖5。圖4分別是內(nèi)涵處于亞臨界和超臨界狀態(tài)、不同馬赫數(shù)時(shí)外涵的流量系數(shù)與推力系數(shù)隨外涵落壓比的變化關(guān)系。可以看出，無論內(nèi)涵阻塞與否，超臨界工況下的外涵道的特性參數(shù)不隨馬赫數(shù)變化而變化，而亞臨界工況下的外涵道的流量系數(shù)與推力系數(shù)則會(huì)隨馬赫數(shù)的增大而增大。

圖5是外涵處于亞臨界和超臨界狀態(tài)、不同馬赫數(shù)時(shí)內(nèi)涵的流量系數(shù)與推力系數(shù)隨內(nèi)涵落壓比的變化關(guān)系?？梢钥闯?，亞臨界工況下的內(nèi)涵道的流量系數(shù)與推力系數(shù)則會(huì)隨馬赫數(shù)的增大而增大，超臨界工況下的內(nèi)涵道的流量系數(shù)與推力系數(shù)不隨馬赫數(shù)的變化而變化，其與外涵的變化趨勢是一致的。

圖4 飛行馬赫數(shù)對外涵特性參數(shù)的影響

圖5 飛行馬赫數(shù)對內(nèi)涵特性參數(shù)的影響

3 結(jié)論

1)本文介紹了獲取尾噴管特性的CFD數(shù)值模擬方法，計(jì)算了某渦扇發(fā)動(dòng)機(jī)的尾噴管特性。分析了長外涵道形式的尾噴管的內(nèi)外涵的耦合影響規(guī)律。

2)研究結(jié)果表明，該形式的尾噴管的特性與外流馬赫數(shù)、噴管落壓比、內(nèi)外涵壓力比有關(guān)。地面靜止條件下，內(nèi)外涵的流量系數(shù)與推力系數(shù)均隨各自落壓比的增大而增大，達(dá)到阻塞狀態(tài)后，流量系數(shù)不再變化，推力系數(shù)略有增加。

3)在靜止與飛行條件下，內(nèi)涵工況對外涵的流量系數(shù)與推力系數(shù)沒有影響，外涵工況對亞臨界狀態(tài)的內(nèi)涵的流量系數(shù)與推力系數(shù)有影響。

4)飛行馬赫數(shù)越大，亞臨界狀態(tài)的內(nèi)涵道與外涵道的流量系數(shù)與推力系數(shù)越大。臨界與超臨界狀態(tài)的內(nèi)涵道與外涵道的流量系數(shù)和推力系數(shù)不受馬赫數(shù)影響。

[1] BURCHAM J F.Use of the gas generator method to calculate the thrust of an afterburning turbofan engine[J].Aiaa Journal，1971.

[2] ASBURY S C，YETTER J A.Static performance of six innovative thrust reverser concepts for subsonic transport applications[R].NASA/TM-2000-210300，2000.

[3] GEATZ AM.A prediction code for the thrust performance of two-dimensional，non-axisynnetric，converging diverging nozzles[C].AFIT/GAE/ENY/06-03，2005.

[4] STRATFORD B S.The calculation of the discharge coefficient of the profiled choked nozzles and the optimum profile for absolute air flow measurement[J].Aeronautical Journal，1964，68 (640) ：237-245.

[5] 朱彥偉，袁長波.大涵道比發(fā)動(dòng)機(jī)噴管流量系數(shù)特性數(shù)值計(jì)算與分析[J].計(jì)算機(jī)仿真，2013，30(1)：159-164.

[6] 齊海帆，高揚(yáng)，郝曉樂，等.某型渦扇發(fā)動(dòng)機(jī)尾噴管流動(dòng)特性研究[J].航空發(fā)動(dòng)機(jī)，2015，41(1)：48-52.

[8] 李密，汪濤，張曼.某型發(fā)動(dòng)機(jī)尾噴管特性確定方法[J].現(xiàn)代機(jī)械，2017(4)：102-106.