陳航

摘要：采用數(shù)值模擬方法計(jì)算分析了射流角度對(duì)氣動(dòng)矢量噴管流動(dòng)結(jié)構(gòu)和性能的影響?；谠囼?yàn)設(shè)計(jì)和近似模型方法，結(jié)合數(shù)值模擬，構(gòu)建了氣動(dòng)矢量噴管模型。將氣動(dòng)矢量噴管模型集成到發(fā)動(dòng)機(jī)總體性能程序中，以某航空發(fā)動(dòng)機(jī)為研究對(duì)象，分析了發(fā)動(dòng)機(jī)不同狀態(tài)下的氣動(dòng)矢量噴管性能情況。研究結(jié)果表明：在一定的幾何與氣動(dòng)參數(shù)條件下，存在使氣動(dòng)矢量角最大的射流角度。隨著發(fā)動(dòng)機(jī)進(jìn)口溫度和引氣量的增加，氣動(dòng)矢量角最大模式下的射流角度均逐漸減小，垂直噴射或順流向?qū)μ嵘龂姽芡屏ο禂?shù)有利。

關(guān)鍵詞：氣動(dòng)矢量噴管;試驗(yàn)設(shè)計(jì);近似模型;射流角度;氣動(dòng)矢量角

中圖分類號(hào)：V233.757? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼：A? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文章編號(hào)：1674-957X（2021）11-0004-02

0? 引言

氣動(dòng)矢量噴管在減輕發(fā)動(dòng)機(jī)重量、減少成本、實(shí)現(xiàn)飛機(jī)后體與發(fā)動(dòng)機(jī)高度一體化方面具有明顯的技術(shù)優(yōu)勢(shì)，是未來高推重比渦扇發(fā)動(dòng)機(jī)的首選技術(shù)方案之一[1]。相對(duì)于機(jī)械式矢量噴管，氣動(dòng)矢量噴管內(nèi)部流動(dòng)復(fù)雜，涉及到二次流噴射對(duì)超聲速主流的摻混、流動(dòng)分離和激波，且影響因素更多，并呈現(xiàn)高度的非線性和多因素交互性。對(duì)于氣動(dòng)矢量噴管，國(guó)內(nèi)外學(xué)者基于數(shù)值模擬和試驗(yàn)技術(shù)開展了大量的研究工作[2-5]，對(duì)其性能和響應(yīng)規(guī)律進(jìn)行不斷的探索，得到了主、次流氣動(dòng)、幾何參數(shù)對(duì)噴管推力矢量的定量影響關(guān)系。

但目前的研究大都集中在氣動(dòng)矢量噴管本身，沒有考慮發(fā)動(dòng)機(jī)和氣動(dòng)矢量噴管的耦合，且不涉及二次流射流角度的控制研究。本文從氣動(dòng)矢量噴管模型入手，開展不同射流角度對(duì)噴管流動(dòng)情況的數(shù)值模擬計(jì)算分析?；谠囼?yàn)設(shè)計(jì)和近似模型構(gòu)建了氣動(dòng)矢量噴管模型，建立了主、次流氣動(dòng)、幾何參數(shù)和噴管性能參數(shù)之間的關(guān)聯(lián)關(guān)系。將氣動(dòng)矢量噴管模型集成到航空發(fā)動(dòng)機(jī)總體性能程序中，使噴管主次流進(jìn)口參數(shù)與發(fā)動(dòng)機(jī)模型相關(guān)聯(lián)，對(duì)發(fā)動(dòng)機(jī)不同狀態(tài)下的噴管推力矢量特性進(jìn)行分析，對(duì)二次流射流角度變化規(guī)律進(jìn)行研究。

1?; 數(shù)值方法

1.1 物理模型與網(wǎng)格劃分

本文研究的氣動(dòng)矢量噴管是結(jié)構(gòu)參數(shù)如表1所示的二元矩形收-擴(kuò)噴管，并通過自編程序控制射流口位置和角度（下壁面為單排射流口）。射流角度的規(guī)定見圖1，以主流方向?yàn)榛鶞?zhǔn)（0°），逆主流方向?yàn)樯淞鹘谴笥?0°方向。

1.2 邊界條件和求解方法

主流和二次流進(jìn)口選用壓力進(jìn)口邊界，給定總溫、總壓，出口選用壓力出口，給定靜壓，噴管外流入口給定黎曼邊界條件，壁面條件設(shè)定為絕熱、無滑移條件。

采用變步長(zhǎng)四階Runge-Kutta方法進(jìn)行時(shí)間推進(jìn)求解，對(duì)流項(xiàng)采用二階迎風(fēng)格式離散求解，粘性項(xiàng)采用中心差分格式。湍流模型選擇SST k-ω湍流模型，同時(shí)由于擴(kuò)張段噴管氣流達(dá)到了超聲速，因此選用基于密度的穩(wěn)態(tài)方程，流體介質(zhì)選用理想氣體。

2? 射流角度對(duì)氣動(dòng)矢量噴管性能影響

射流角θ會(huì)對(duì)次流與主流摻混時(shí)的動(dòng)量及能量交換產(chǎn)生較大影響，從而影響到摻混流場(chǎng)結(jié)構(gòu)，并最終影響噴管矢量性能。選取二次流相對(duì)位置0.7，NPR=8、10、12，SPR=0.8、1.0、1.2時(shí)，對(duì)θ取90～130°時(shí)的噴管矢量性能作了研究。

隨著射流角度θ的增大，氣動(dòng)矢量角增大，可見射流角度θ對(duì)噴管流場(chǎng)結(jié)構(gòu)有較大影響。但當(dāng)射流角度θ增大到一定程度時(shí)，在射流角度θ=110°、120°時(shí)氣動(dòng)矢量角基本相同，噴管內(nèi)部流場(chǎng)結(jié)構(gòu)變化不再明顯，繼續(xù)增大射流角度，氣動(dòng)矢量角反而降低，即在一定的幾何與氣動(dòng)參數(shù)條件下，存在一個(gè)最佳的噴射角度，使得噴管氣動(dòng)矢量角達(dá)到最大。

3? 氣動(dòng)矢量噴管模型研究

3.1 氣動(dòng)矢量噴管試驗(yàn)設(shè)計(jì)方法

影響氣動(dòng)矢量噴管的因素很多，DOE試驗(yàn)設(shè)計(jì)方法以其獨(dú)有的多維數(shù)組均值算法可以有效地應(yīng)對(duì)多變量、大范圍的影響參數(shù)，合理安排試驗(yàn)參數(shù)的組合，使得試驗(yàn)點(diǎn)大大減少，且具有較好的可信度。本文中選用的最優(yōu)拉丁采樣算法，該方法保證每個(gè)因子選取的不可重復(fù)性，并可以使所有的試驗(yàn)點(diǎn)盡量均勻地分布在設(shè)計(jì)空間，具有非常好的空間填充性和均衡性。

為研究二次流射流口位置和角度對(duì)氣動(dòng)矢量噴管性能的影響，通過自編程序用來控制二次流射流口位置和角度，搭建了基于DOE試驗(yàn)設(shè)計(jì)方法的氣動(dòng)矢量噴管特性計(jì)算流程。

影響氣動(dòng)矢量噴管性能的因素主要有：

主流落壓比（NPR）：4.5～18;二次流總壓比（SPR）：0.8～1.2;二次流射流口相對(duì)位置（x/x1，相對(duì)于擴(kuò)張段軸向長(zhǎng)度）：0.5～0.8;射流角度（θ，相對(duì)于主流方向）：90～140°。

對(duì)以上因素對(duì)噴管性能的影響進(jìn)行計(jì)算和分析。

對(duì)以上4個(gè)因素分別取值86個(gè)試驗(yàn)點(diǎn)，生成設(shè)計(jì)矩陣，并進(jìn)行計(jì)算。

3.2 氣動(dòng)矢量噴管近似模型構(gòu)建

為對(duì)比不同近似模型的精度，本文分別采用徑向基RBF神經(jīng)網(wǎng)絡(luò)模型及多項(xiàng)式回歸模型對(duì)上述得到的氣動(dòng)矢量噴管數(shù)值模擬結(jié)果進(jìn)行近似建模，同時(shí)做了二階、三階和四階的對(duì)比分析，對(duì)擬合精度和可信度進(jìn)行比較分析。

對(duì)樣本數(shù)據(jù)擬合度可用均方根誤差（RMSE）相對(duì)值和決定系數(shù)R2兩個(gè)標(biāo)準(zhǔn)檢驗(yàn)，RMSE→0表示響應(yīng)值誤差小，R2→1則表明響應(yīng)值與原模型相似度高。

徑向基神經(jīng)網(wǎng)絡(luò)模型以其強(qiáng)大的逼近復(fù)雜非線性函數(shù)的能力偏差最小，擬合精度最高。也說明氣動(dòng)矢量噴管具有高度的非線性及相互的關(guān)聯(lián)性，利用多項(xiàng)式回歸方法很難完整的表達(dá)出其特性。

4? 氣動(dòng)矢量噴管二次流射流角度的優(yōu)化

對(duì)于氣動(dòng)矢量噴管，二次流射流位置是固定的，可調(diào)節(jié)幾何參數(shù)的只有射流角度，開展氣動(dòng)矢量噴管射流角度控制規(guī)律的優(yōu)化，可以為氣動(dòng)矢量噴管的設(shè)計(jì)和控制提供參考。

將構(gòu)建的氣動(dòng)矢量噴管模型集成到發(fā)動(dòng)機(jī)總體性能程序中，構(gòu)建二次流引氣流路及噴管主流進(jìn)口參數(shù)接口，以某高推重比發(fā)動(dòng)機(jī)為研究對(duì)象，進(jìn)行氣動(dòng)矢量噴管二次流射流角度規(guī)律的研究。

計(jì)算條件：二次流射流口相對(duì)位置x/x1=0.7;優(yōu)化變量：二次流入射角度θ;變量范圍：90～140°;引氣位置：外涵道進(jìn)口。

二次流引氣量為外涵進(jìn)口截面的5～15%，開展不同引氣量對(duì)氣動(dòng)矢量噴管性能影響的計(jì)算分析。隨著引氣量的增加，使得噴管內(nèi)產(chǎn)生的斜激波強(qiáng)度增強(qiáng)，氣動(dòng)矢量角增加，推力系數(shù)減少，氣動(dòng)矢量角最大優(yōu)化模式的二次流射流角度逐漸減小，在引氣量為外涵進(jìn)口截面的5%增大到15%時(shí)，矢量角最大模式下的二次流射流角度從130°減小到103°，推力系數(shù)最大模式下二次流射流角度依然保持優(yōu)化范圍最小角度90°。

在引氣位置固定之后，二次流射流角度是氣動(dòng)矢量噴管的重要影響因素，因此，在氣動(dòng)矢量噴管的設(shè)計(jì)過程中，需要對(duì)發(fā)動(dòng)機(jī)不同工作狀態(tài)下的氣動(dòng)矢量角和推力系數(shù)進(jìn)行綜合考慮，提升發(fā)動(dòng)機(jī)性能。

5? 結(jié)論

①在一定的幾何與氣動(dòng)參數(shù)條件下，存在一個(gè)最佳的噴射角度，使得噴管氣動(dòng)矢量角最大。②氣動(dòng)矢量噴管具有高度的非線性及多影響因素的交互性，通過與多項(xiàng)式回歸模型的對(duì)比，徑向基神經(jīng)網(wǎng)絡(luò)模型更適合用于表達(dá)氣動(dòng)矢量噴管特性。③隨著發(fā)動(dòng)機(jī)進(jìn)口溫度和引氣量的增加，氣動(dòng)矢量角最大對(duì)應(yīng)的二次流射流角度均逐漸減小，垂直主流噴射或順流向?qū)μ嵘龂姽芡屏ο禂?shù)有利，氣動(dòng)矢量噴管的設(shè)計(jì)要在氣動(dòng)矢量角和推力系數(shù)中進(jìn)行綜合選取。

參考文獻(xiàn)：

[1]Deere K A Computational Investigation of the Aerodynamic effects on Fluidic Thrust Vectoring. AIAA2000-3598，2000.

[2]Deere K A. Summary of fluidic thrust vectoring research conducted at NASA Langley Research Center[R]. AIAA 2003-3800.

[3]賈東兵，周吉利，鄧洪偉.固定幾何氣動(dòng)矢量噴管技術(shù)綜述[J].航空發(fā)動(dòng)機(jī)，2012（6）：29-33，42.

[4]王占學(xué)，李志杰.噴管主次流氣動(dòng)參數(shù)對(duì)流體推力矢量影響的數(shù)值模擬[J].推進(jìn)技術(shù)，2008，29（2）.

[5]王占學(xué)，王玉男，李志杰，辛?xí)札?基于激波控制的流體推力矢量噴管試驗(yàn)[J].推進(jìn)技術(shù)，2010，31（6）：751-756.