吸氣式旋轉爆震組合循環發動機性能

計自飛，李天琦，張會強

(1. 中國航發沈陽發動機研究所，沈陽 110015； 2. 清華大學 航天航空學院，北京 100084)

0 引言

渦輪基組合循環發動機(TBCC)是最典型的吸氣式組合循環發動機，在過去的幾十年里，國內外學者通過理論分析、數值仿真以及整機和部件實驗等手段開展了大量研究工作，并取得了一系列重要成果。然而，由于常規渦輪的高速性能以及常規沖壓的低速性能均較差，常規TBCC發動機在模態轉換期間存在“推力間隙”問題。為了解決上述問題，往往需要引入復雜系統或熱力循環。西安航天動力研究所提出了預冷空氣渦輪火箭發動機(PATR)方案。受壓氣機出口溫度限制，常規的ATR發動機最高飛行馬赫數僅為3.5左右，PATR發動機的進氣道與壓氣機之間設置有預冷器，可降低高馬赫數條件下壓氣機出口溫度，從而具有更寬廣的工作包線。北京動力機械研究所提出的渦輪輔助火箭增強沖壓組合循環發動機(TRRE)則是通過將并聯式TBCC發動機的沖壓燃燒室替換為火箭沖壓復合燃燒室來實現渦輪模態與沖壓模態的平穩接力。英國反作用發動機公司的吸氣火箭協同式發動機(SABRE)是一類特殊的預冷空氣渦輪發動機，該發動機存在吸氣和火箭兩個工作模塊，涉及液氮、氫、空氣、富燃燃氣等4種工質的循環，工作過程相當復雜。

旋轉爆震是一種超聲速增壓燃燒方式，同時具有熵增小、反應強度高、放熱速率快等優勢，將爆震技術引入推進系統可顯著提升發動機的推力和油耗特性，同時使得發動機結構更加緊湊，因而，旋轉爆震發動機是最具潛力的新型航空航天動力裝置之一。旋轉爆震技術在燃氣渦輪發動機和沖壓發動機等單一循環動力系統中應有的優勢已得到驗證，用基于旋轉爆震的渦輪動力單元和沖壓動力單元分別代替組合循環發動機中的常規渦輪動力單元和沖壓動力單元，將可能在不增加系統和循環復雜度的條件下改善吸氣式組合循環發動機在過渡模態的性能。然而，目前有關旋轉爆震發動機的研究主要是針對單一循環方式，有關基于旋轉爆震的渦輪基組合循環發動機的研究十分匱乏。

本文以吸氣式組合循環發動機為研究對象，提出了一種基于旋轉爆震的渦輪基組合循環發動機系統方案，能夠充分發揮每種組成單元在不同飛行條件下的性能優勢，實現飛行器的全速域、跨空間的快速飛行。隨后，建立了參數化總體性能分析模型，在此基礎上研究了發動機在不同飛行參數下的總體性能特征。

1 數學物理模型

1.1 旋轉爆震燃燒室低階計算模型

旋轉爆震過程具有顯著的時空非定常效應，針對常規燃燒方式的低階模型不再適用，建立連續旋轉爆震過程的低階模化方法，是開展旋轉爆震發動機參數化總體性能分析的關鍵環節。旋轉爆震燃燒室的典型二維流場結構如圖1所示，其中的熱力學過程可看作是旋轉爆震過程緊跟著膨脹過程。在靜止坐標系中，由于膨脹波的影響，爆震波前預混氣具有一定的周向速度分量，因而，爆震波面并不與噴注面垂直。CJ速度

v

與進氣速度

v

以及爆震波的周向傳播速度

U

之間的關系可由圖1(b)所示的速度三角形表征。

圖1 旋轉爆震燃燒室典型流場示意圖Fig.1 Schematic diagram of typical flow field in rotating detonation combustor

對于旋轉爆震過程，若來流參數

p

、

T

已知，利用MATLAB軟件耦合Shock and Detonation Tool Box以及CANTERA軟件基于一維ZND理論可以求解爆震波后參數

p

、

T

以及理論速度

v

。假定新鮮混合氣等速進入燃燒室，則爆震波面與新鮮混氣和產物的分界面垂直。Sichel等基于特征線方法得到了在爆震波傳播方向上距離其

x

處的壓力與無量綱長度

ξ

的關系，其中

ξ

是位置坐標

x

與爆震波頭的高度

h

的比值。本模型選用Sichel等的研究結論，結合一維ZND理論和爆震波后壓力的衰減模型，可以得到噴注面的壓力分布

p

(

θ

)。根據

p

(

θ

)可以得到進氣區的范圍Δ

θ

。燃燒室的壓力從爆震波所處的角向位置

θ

開始沿周向衰減，在某一角向位置

θ

處壓力衰減到噴注總壓

p

，可燃混氣開始進入爆震室，從此處開始燃燒室壓力等于新鮮混氣的壓力。相關研究表明，附體斜激波的強度較弱，由于斜激波造成的熵增相比于爆震過程的熵增非常小，同時，由斜激波造成的總壓損失與加熱造成的總壓損失相比也非常小，因此，爆震產物的膨脹過程可按照等熵過程處理。圖2給出了根據本模型計算的燃燒室壓力與角度坐標的函數關系和數值計算結果的對比情況，可見本模型的預測結果與數值仿真的結果吻合很好。

圖2 旋轉爆震燃燒室低階模型計算結果與數值模擬結果的比較Fig.2 Comparison of low-order model calculation results with numerical simulation results for rotating detonation combustor

對于進氣道、壓氣機、渦輪、尾噴管中的氣動熱力過程采用常規航空燃氣渦輪發動機的通用方法計算和分析，具體模型此處不再贅述。

1.2 系統布局與工作原理

渦輪基組合循環發動機的布局有串聯式和并聯式兩種，由于串聯式布局在單元匹配工作方面存在較大困難，本文選用并聯布局。考慮到旋轉爆震燃燒室的工作特性及其與進排氣部件的匹配工作要求，本文提出一種基于旋轉爆震的渦輪基組合循環發動機系統方案，如圖3所示。整體采用雙通道布局，旋轉爆震沖壓通道位于旋轉爆震燃氣渦輪通道的外側，兩個通道共用進氣道和尾噴管膨脹段。在進氣道部分設有模式選擇閥，通過調節其位置可以實現組合發動機在不同工作模式間的轉換。在尾噴管部分設有尾噴管調節閥，通過調節該閥門的位置實現單元間的匹配工作。

圖3 基于旋轉爆震的渦輪基組合循環發動機系統方案Fig.3 Turbine-based combined cycle engine based on rotating detonation

根據高超聲速飛行器在不同飛行條件下對動力系統的需求，該組合循環發動機可實現3種不同的工作模式：①渦輪模式，當飛行馬赫數較低時，組合循環發動機以渦輪模式工作，如圖4(a)所示，調節模式選擇閥使得雙通道均處于打開狀態，然而僅有旋轉爆震渦輪單元起動工作；②過渡模式，隨著飛行馬赫數的增大，渦輪單元的性能顯著下降，當其推力難以滿足飛行器需求時，沖壓通道點火，此后渦輪單元的推力占比逐漸減小，直至完全關機,如圖4(b)所示；③沖壓模式，當飛行馬赫數較高時，組合循環發動機以沖壓模式工作，如圖4(c)所示，調節模式選擇閥使得渦輪單元通道處于關閉狀態，僅沖壓單元啟動工作。

圖4 組合循環發動機的工作模式Fig.4 Operating mode of combined cycle engine

2 結果與討論

2.1 旋轉爆震渦輪發動機總體性能特征

圖5 旋轉爆震渦輪發動機與常規渦輪發動機的性能比較Fig.5 Performance comparison between rotating detonation turbine engine and traditional turbine engine

2.2 旋轉爆震沖壓發動機總體性能特征

圖6 旋轉爆震沖壓發動機與常規沖壓發動機的性能比較Fig.6 Performance comparison between rotating detonation ramjet engine and traditional ramjet engine

2.3 旋轉爆震組合循環發動機總體性能特征

由2.1節和2.2節的結論可知，在

Ma

較低時，爆震渦輪模式的比推力和比沖均具有顯著優勢，隨著

Ma

的增大，爆震渦輪模式的比推力優勢逐漸被爆震沖壓模式取代，然而，前者的比沖仍高于后者。若在某一固定馬赫數下完成模態轉換，模態轉換馬赫數過小，則比沖偏低，模態轉換馬赫數過大，則比推力偏小。考慮到上述問題，本文提出一種等推力等流量模態轉換策略，即模態轉換在一個速度區間內完成，在該區間內，組合發動機處于過渡工作模態，且推力和進氣流量恒定。這種模態轉換策略可兼顧渦輪單元比沖高和沖壓單元比推力高的優勢。按照等推力模態轉換策略，在模態轉換期間，通過調節渦輪單元和沖壓單元的供油規律和噴管調節閥的位置，實現在模態轉換期間推力和流量保持不變。定義渦輪通道的空氣流量占組合循環發動機總進氣流量的比例為流量分配系數,用符號

Λ

表示，其表達式為

(1)

式中

F

表示過渡模態比推力的設計值。根據吸氣式發動機比沖的定義，組合循環發動機在過渡模態的等效比沖可由

Λ

、渦輪單元比沖

I

和沖壓單元比沖

I

表示為

(2)

不同飛行高度對應的旋轉爆震組合循環發動機過渡模態性能特征如圖7所示。可見，不同飛行高度條件下，相同的模式轉換起始馬赫數(

Ma

)對應不同的終止馬赫數(

Ma

)，

H

越大，過渡模式對應的速度區間越小。由于沖壓模式的比推力隨飛行高度

H

的變化不如渦輪模式顯著，因而，相同

Ma

對應的

Ma

隨

H

的增大而減小。從圖7中可以看出，

H

越大，相同

Ma

對應的過渡模式的比推力越大。渦輪單元和沖壓單元的比沖均隨

H

的增大而增大，但是，渦輪單元流量占比

Λ

隨

H

的增大而減小，在兩者的共同作用下，在過渡模式的起始階段，等效比沖幾乎不隨

H

而變化，如圖7(b)所示。此外，由于11～20 km為標準大氣的等溫層，因此，發動機在15 km和20 km的性能非常接近。

圖7 定高巡航工況下H取值對組合循環發動機總體性能和進氣道流量分配的影響Fig.7 Influence of H value on the overall performance of combined cycle engine and the inlet flow distribution under the constant height cruise condition

3 結論

通過熱力循環分析方法研究了旋轉爆震渦輪發動機和旋轉爆震沖壓發動機的總體性能特征，與相同使用條件的基于常規燃燒方式的發動機性能對比表明，在寬廣的飛行馬赫數范圍內，旋轉爆震渦輪發動機的比推力和比沖均有明顯優勢，且這些優勢隨著

Ma

的增大而減小；在低

Ma

工況下，旋轉爆震沖壓發動機的比推力和比沖具有優勢，隨著

Ma

的增大，比推力優勢逐漸減小直至消失。建立了基于旋轉爆震的渦輪基組合循環發動機參數化總體性能分析模型，針對爆震渦輪單元和爆震沖壓單元的性能特征，提出一種等推力等流量模態轉換策略，并研究了在定高巡航飛行方式下，模態轉換對應的飛行馬赫數范圍以及過渡模態性能特征，結果表明，當模態轉換起始馬赫數相同時，模態轉換終止馬赫數隨飛行高度的增大而減小；隨著飛行高度的增大，相同

Ma

對應的比推力單調增大，渦輪單元流量占比單調減小，在過渡模式的起始階段，等效比沖幾乎不隨高度而變化。