一種超音速分離線擺動噴管流動特性數(shù)值研究*

童 悅，鄭 慶，鄒 杰，李修明，陳振陽

(上海航天動力技術(shù)研究所， 上海 201109)

0 引言

與燃?xì)舛媸噶靠刂萍夹g(shù)相比，擺動噴管技術(shù)具有置偏能力強、軸向損失小、結(jié)構(gòu)簡單等優(yōu)點，是高性能固體動力裝置推力矢量控制的重點發(fā)展方向之一。自上世紀(jì)五十年代末起，擺動噴管技術(shù)得到了快速發(fā)展[1]，在現(xiàn)役戰(zhàn)術(shù)導(dǎo)彈固體發(fā)動機中已成熟應(yīng)用：如美國標(biāo)準(zhǔn)-3導(dǎo)彈的MK-72助推發(fā)動機、MK-136雙脈沖發(fā)動機[2]及THAAD導(dǎo)彈一級助推發(fā)動機。

國外對新型擺動噴管的研究較早，美國于1969年開展超音速分離線(supersonic split line，SSSL)推力矢量技術(shù)的研究工作，并通過六分力試驗臺研究了噴管矢量角與擺動角的關(guān)系[3]。挪威Nammo公司使用5 s發(fā)動機先后驗證并獲得了全軸/單軸超音速分離線擺動噴管的密封、燒蝕、擺動力矩及矢量特性[4]。為進一步提高Orbus1發(fā)動機性能，文獻[5]提出了非潛入式超音速分離線柔性擺動噴管概念，并開展了相關(guān)試驗工作，研究結(jié)果表明：優(yōu)化后噴管矢量角放大系數(shù)在1.25以上，擺動比力矩降低29%，噴管重量減輕43%。

國內(nèi)對新型高效擺動噴管的研究起步較晚，目前研究方向仍主要集中在柔性擺動噴管[6-10]、珠承擺動噴管[11-14]及球窩擺動噴管[15-16]等亞音速分離線擺動噴管方面。

小型、輕質(zhì)、高效、低力矩始終是國內(nèi)外擺動噴管技術(shù)的發(fā)展主線。文中對一種超音速分離線擺動噴管開展了全三維數(shù)值仿真計算，分析了其流動結(jié)構(gòu)，并對不同擺角與分離線間隙距離的影響規(guī)律開展研究，為后續(xù)超音速分離線擺動噴管氣動型面及結(jié)構(gòu)設(shè)計提供依據(jù)。

1 研究對象

圖1 擺動噴管三維結(jié)構(gòu)

圖2 擺動噴管內(nèi)型面結(jié)構(gòu)示意圖

2 仿真計算

2.1 計算網(wǎng)格

由于分離線位于喉襯下游超音速區(qū)域，下游流場對喉襯上游亞音速流動無影響，故選取擴張段與喉襯上游的部分流道作為研究對象。同時，待研究的流動關(guān)于擺動噴管鉸接軸中剖面對稱，為此選取模型的一半作為計算域進行網(wǎng)格劃分。整個計算域采用六面體網(wǎng)格單元進行填充，并且在噴管近壁面和分離線等型面變化劇烈的區(qū)域進行網(wǎng)格加密(對稱面及壁面網(wǎng)格見圖3所示)，網(wǎng)格單元總數(shù)約72萬，大部分近壁面網(wǎng)格單元的y+保持在30～100左右。

2.2 計算方法及設(shè)置

采用Fluent軟件進行數(shù)值仿真，其中使用的有限體積法對雷諾平均三維Navier-Stokes(N-S)方程進行離散，無黏對流通量采用基于MUSCL差值的Roe格式進行裂解，而時間推進采用了隱式格式，湍流黏度采用Shear-Stress Transport(SST)k-ω模型進行計算。

圖3 所采用的計算網(wǎng)格(擺角θ=5°)

仿真中各邊界條件給定如下：噴管入口設(shè)置為壓力入口邊界，出口設(shè)置為壓力出口邊界條件。計算過程中對燃?xì)庾魅缦录僭O(shè)：燃?xì)鉃榧儦庀嗖⑶以趪姽苤械牧鲃訛閮鼋Y(jié)流，不考慮熱輻射、內(nèi)壁面燒蝕等影響，計算參數(shù)見表1所示。

表1 計算參數(shù)

2.3 性能監(jiān)控參數(shù)

采用軸向推力系數(shù)Cfx和放大系數(shù)K[3]，對超音速分離線擺動噴管的推力矢量進行定量評價，其定義分別為：

(1)

(2)

式中：Fa為擺動噴管軸向推力，F(xiàn)0為0°時噴管軸向推力，δ為擺動噴管矢量角，θ為擺動噴管幾何擺動角。軸向推力與矢量角按下式計算：

(3)

(4)

3 結(jié)果與分析

3.1 流動分析

對超音速分離線噴管擺動5°，分離線間隙距離Δ=0.03Dt時的內(nèi)流場進行仿真研究。圖4給出了擺動噴管對稱面馬赫數(shù)及壓強分布圖譜，可以看出，燃?xì)庠诤硪r處為跨音速流動，流動較為順暢，燃?xì)鈮簭姴粩嘟档?，速度上升；喉襯后緣分離線處為超音速流動，馬赫數(shù)約為1.5。因噴管擺動后擴張段上下型面不連續(xù)，下端型面外凸，上端型面內(nèi)凹；型面外凸處在喉襯后緣分離線前端形成斜激波，而內(nèi)凹處形成膨脹波；斜激波與膨脹波后燃?xì)夥较蚓驀姽軘[動方向偏轉(zhuǎn)，增強了擺動噴管的矢量效果，但擴張段內(nèi)燃?xì)饬鲃痈鼜?fù)雜，存在明顯的激波/膨脹波/邊界層相互干擾區(qū)，流動損失增大。此外，激波后分離線前燃?xì)獯嬖诘退倩亓鲄^(qū)，其壓強遠(yuǎn)高于膨脹波后燃?xì)鈮簭?，燃?xì)鈩荼卦诟邏簭娞荻茸饔孟拢纬蓹M向流動。圖5較清晰地顯示出超音速分離線擺動噴管內(nèi)流道的流動狀態(tài)：擴張段下端面一部分燃?xì)饨?jīng)過激波后氣流方向發(fā)生偏轉(zhuǎn)，向后流動；一部分燃?xì)庠诜蛛x線處不斷聚積，形成旋渦，并通過分離線間隙通道向擴張段上端面匯集后排出。

圖4 對稱面馬赫數(shù)及壓強分布云圖

3.2 擺角影響規(guī)律研究

為研究超音速分離線擺動噴管擺角對推力矢量特性的影響，對分離線間隙距離Δ=0.03Dt時，擺角θ變化范圍為1°～7°的噴管內(nèi)流場進行數(shù)值研究。

圖6為噴管擺角與軸向推力系數(shù)、矢量角放大系數(shù)的變化關(guān)系曲線。不難看出:研究范圍內(nèi)矢量角放大系數(shù)均大于1.1，即超音速分離線擺動噴管推力矢量角大于噴管擺角，且最大值為噴管擺動角度的2倍。同時，隨著擺動角度增大，超音速分離線擺動噴管軸向推力系數(shù)與矢量角放大系數(shù)均快速下降，擺角θ為7°時，軸向推力系數(shù)僅為0.955，矢量角放大系數(shù)為1.19。對比不同擺角下，超音速分離線擺動噴管對稱面馬赫數(shù)分布云圖(圖7所示)，可以明顯看出：隨著擺角增大，擴張段外凸/內(nèi)凹造成型面不連續(xù)性更加嚴(yán)重，擴張段斜激波激波角增大并向前移動，增大了流動損失；同時，不斷增大的低速高壓旋流區(qū)加強了分離線間隙處的周向流動，使得更多高壓燃?xì)饬鬟M入擴張段上端面，降低了擴張段上下壁面壓差，從而降低了擺動噴管的矢量放大系數(shù)。因此，矢量角放大系數(shù)隨著擺動角度增大而減小。

圖5 超音速分離線擺動噴管內(nèi)流道三維流線

圖6 不同擺角矢量特性曲線

3.3 分離線間隙距離影響規(guī)律研究

擺動噴管為機械式推力矢量調(diào)節(jié)裝置，為保證陰陽半球相對運動順暢，設(shè)置有分離線。分離線間隙過大，擺動噴管動密封可靠性降低，但分離線間隙過小，球面會因高溫高壓結(jié)構(gòu)變形出現(xiàn)“卡死”現(xiàn)象。為研究分離線間隙大小對超音速分離線擺動噴管流動及推力矢量特性的影響，開展相關(guān)數(shù)值分析，研究范圍內(nèi)分離線間隙Δ取值范圍為0～0.06Dt。

圖8給出了不同分離線間隙下擺動噴管矢量特性曲線。各擺角狀態(tài)下，陰陽球體無間隙即Δ=0時，軸向推力系數(shù)與矢量角放大系數(shù)為最大值。隨著分離線間隙的增大，各擺角軸向推力系數(shù)與矢量角放大系數(shù)快速下降，且擺角越大推力系數(shù)下降越明顯，在研究范圍內(nèi)，軸向推力系數(shù)從0.962下降至0.947。不同擺角下分離線間隙引起的矢量角放大系數(shù)變化量基本一致，即間隙從0增至0.03Dt，各擺角下的矢量角放大系數(shù)小幅下降，僅下降2%～7%；間隙距離增至0.06Dt，各擺角下的放大系數(shù)下降幅度較大，下降幅度為11%～20%。

圖7 不同擺角對稱面馬赫數(shù)分布云圖

圖8 不同分離線間隙矢量特性曲線

圖9 壁面靜壓及三維流線分布對比

進一步，對比不同分離線間隙距離下噴管壁面靜壓分布及三維流線圖，明顯看出：分離線無間隙時，喉襯下端低速高壓旋流沿分離線前端周向匯聚后排出，旋流周向跨度大，喉襯壁面高壓影響區(qū)域大(圖9中虛線標(biāo)出)。而間隙增大，大量低速高壓旋流沿分離線間隙處向上排出，喉襯壁面高壓區(qū)域明顯大幅減小，同時，噴管擴張段壁面壓強分布基本不變，造成噴管側(cè)向力降低，矢量角放大系數(shù)下降。因此，小分離線間隙利于擺動噴管矢量性能的保持。

4 結(jié)論

文中對一種超音速分離線擺動噴管的流動特性進行仿真分析，并對比了不同擺角、不同分離線間隙對擺動特性的影響，獲得的結(jié)論主要如下：

1)超音速分離線擺動噴管擺動后，擴張段上下型面不連續(xù)，外凸處形成斜激波，內(nèi)凹處形成膨脹波；斜激波與膨脹波增強了矢量特性，但流動更復(fù)雜，存在明顯的激波/膨脹波/邊界層相互干擾區(qū)，流動損失增大。

2)不同擺角下超音速分離線擺動噴管推力矢量角大于噴管擺角，同時，隨著擺角增大，超音速分離線擺動噴管軸向推力系數(shù)與矢量角放大系數(shù)均快速下降，擺角θ為7°時，軸向推力系數(shù)為0.955，矢量角放大系數(shù)達1.19。

3)分離線間隙增大，喉襯后緣處的低速高壓旋流逐漸沿間隙處向上排出，使得喉襯高壓區(qū)大幅減小，噴管側(cè)向力減小，因此各擺角下，軸向推力系數(shù)與矢量角放大系數(shù)隨著分離間隙增大而快速下降。