張喆 陳志明 張磊揚 郝雯

摘要：本文探討了空氣進氣畸變指數(shù)對雙水平進氣布局固體火箭沖壓發(fā)動機補燃室的性能影響。由于在不同工作姿態(tài)下，固體火箭沖壓發(fā)動機進氣道出口截面的流場結(jié)構(gòu)不同，設(shè)計對應(yīng)不同進氣道流場結(jié)構(gòu)的空氣進氣畸變模擬裝置。根據(jù)進氣畸變圖譜形狀和畸變指數(shù)兩個重要因素進行方案設(shè)計，并采用數(shù)值計算的方法，對不同空氣畸變圖譜、進氣堵塞比的補燃室摻混燃燒性能進行仿真分析。結(jié)果表明，相對于基準工況，隨著堵塞比的增大，補燃室的總壓恢復(fù)系數(shù)呈現(xiàn)出先小幅減小后快速減小的趨勢；在相同畸變圖譜下，隨著堵塞比的增大，補燃室溫升效率基本呈現(xiàn)先小幅減小后增大的趨勢，且堵塞比在42%時溫升效率均低于無畸變情況下的基準工況。

關(guān)鍵詞：固體火箭沖壓發(fā)動機；進氣道；補燃室；進氣畸變指數(shù)

中圖分類號：TP391.9文獻標識碼：ADOI：10.19452/j.issn1007-5453.2022.02.004

固體火箭沖壓發(fā)動機把固體火箭發(fā)動機和吸氣式發(fā)動機結(jié)合在一起，具有比沖高、結(jié)構(gòu)緊湊、機動性好、可靠性高的優(yōu)點[1-3]。進氣道作為固沖發(fā)動機的關(guān)鍵部件之一，其性能和工作狀態(tài)對固沖發(fā)動機補燃室的燃燒性能有很大影響，固沖發(fā)動機進氣道出口流場的畸變效應(yīng)，會傳遞到補燃室，發(fā)動機也更容易發(fā)生不穩(wěn)定燃燒、進氣道回火等故障[4-10]。因此，在進行固沖發(fā)動機設(shè)計時需要關(guān)注空氣進氣畸變對補燃室摻混燃燒性能的影響[11]。

固體火箭沖壓發(fā)動機復(fù)雜的進氣道型面會導(dǎo)致燃燒室入口存在較大的流場畸變。為探討空氣進氣畸變對燃燒性能的影響，采用畸變模擬裝置替代進氣道模擬燃燒室入口流場?，F(xiàn)有空氣進氣畸變裝置多以航空發(fā)動機為研究對象，其中插板式畸變模擬裝置因其具有結(jié)構(gòu)簡單、通用性好等優(yōu)點而被普遍應(yīng)用[12]。除此之外，進氣畸變模擬裝置還包括金屬網(wǎng)、多層板及擾流板等結(jié)構(gòu)。然而這些進氣畸變模擬裝置雖在畸變度大小的模擬上都與實際流場存在較高的一致性，但其對流場分布的模擬度都相對較差[13-14]。

在固沖發(fā)動機中，進氣道與燃燒室緊鄰，進氣道出口流場分布規(guī)律直接影響燃燒室性能，本文設(shè)計了相應(yīng)的固沖發(fā)動機補燃室進氣畸變模擬裝置，對應(yīng)進氣道出口截面的不同流場結(jié)構(gòu)，以此研究空氣進氣畸變對雙水平進氣固沖發(fā)動機補燃室摻混燃燒性能的影響，為固體火箭沖壓發(fā)動機的設(shè)計提供技術(shù)支撐。

1物理模型

1.1固沖發(fā)動機模型

本文所用物理模型以雙水平進氣布局固沖發(fā)動機補燃室為研究對象[15]，如圖1所示。取1/2模型進行計算，該模型僅將固沖發(fā)動機補燃室空氣進氣彎頭前的等直延長段作為空氣進氣通道，為避免進氣道帶來的耦合影響，在模型中不考慮進氣道。空氣進氣畸變模擬裝置安裝于補燃室空氣進氣的等直延長段內(nèi)，距彎頭進口約為3倍等效直徑。

1.2畸變裝置設(shè)計

空氣進氣畸變主要有分布圖譜和畸變指數(shù)兩個具體指標，用以衡量進氣道出口截面上氣流的不均勻程度。本文中，進氣畸變模擬裝置采用多個條板堵塞來流的設(shè)計方式，通過改變條板的數(shù)量、位置和寬度來調(diào)節(jié)畸變圖譜形狀和畸變指數(shù)。

雙水平進氣布局固沖發(fā)動機進氣道出口總壓恢復(fù)系數(shù)分布分別呈單點峰和條形峰兩種狀態(tài)，且即使在無側(cè)滑工況下峰值在截面上的位置也稍有差異。本文設(shè)計了與補燃室徑向呈30o的斜條峰馬赫數(shù)Ma分布圖譜（簡稱30）和與補燃室徑向呈45o的近角區(qū)單點峰馬赫數(shù)Ma分布圖譜（簡稱45）?；诖诉M行不同畸變指數(shù)對固沖發(fā)動機燃燒室性能的影響研究，如圖2所示。

圖2中空氣進氣截面上藍色部分為流通區(qū)域，灰色部分為堵塞區(qū)域。堵塞比是指氣流流動垂直方向上堵塞區(qū)域的截面面積與總截面面積的比值，本文通過設(shè)置不同的堵塞比來調(diào)整空氣進口畸變指數(shù)的變化，針對兩種畸變圖譜分別設(shè)計了5種堵塞比，分別為18%、28%、42%、56%和66%。

1.3計算模型及邊界條件

固體火箭沖壓發(fā)動機的燃燒為復(fù)雜的湍流多相非預(yù)混燃燒，因此在進行數(shù)值模擬時需要對固沖發(fā)動機補燃室的燃燒過程進行一定的假設(shè)和簡化。本文采用PDF湍流燃燒計算模型。

本文采用ICEM CFD16.0軟件對計算區(qū)域進行網(wǎng)格劃分，采用混合網(wǎng)格生成技術(shù)，可大大減少網(wǎng)格數(shù)量，從而提高計算的精度和效率。對摻混裝置之前的部分和畸變模擬裝置的部分采用非結(jié)構(gòu)網(wǎng)格，其余部分均采用結(jié)構(gòu)化網(wǎng)格，網(wǎng)格數(shù)量約為135萬，如圖3所示。

計算模型的入口邊界條件為空氣、燃氣質(zhì)量入口，出口邊界條件為沖壓噴管壓力出口，其余邊界條件為對稱面和固體壁面邊界，具體邊界條件見表1。

1.4截面位置分布

本文將流向方向定義為X軸方向，分別選取固沖發(fā)動機補燃室進氣彎頭內(nèi)沿流向依次作為彎頭進口截面、彎頭出口截面1、彎頭出口截面2。同時將X=0mm作為補燃室頭部，在固沖發(fā)動機補燃室內(nèi)截取6個截面，截面位置分別為250mm、350mm、450mm、550mm、700mm、900mm。其中，X為250mm的截面位于空氣進口下游附近，X為900mm的截面為燃燒室出口附近，截面位置圖如圖4所示。

2計算結(jié)果分析

為分析空氣進氣畸變對固體火箭沖壓發(fā)動機補燃室燃燒性能的影響，首先計算了空氣進氣堵塞比與畸變指數(shù)的對應(yīng)關(guān)系，并對不同畸變指數(shù)和畸變圖譜下的總壓恢復(fù)系數(shù)、補燃室溫升效率及沿程截面靜溫分布進行計算。

2.1畸變指數(shù)和堵塞比的對應(yīng)關(guān)系

通過仿真計算得到空氣進氣堵塞比和畸變指數(shù)之間的對應(yīng)關(guān)系，如圖5所示，可以看出，隨著堵塞比的增大，30畸變圖譜的Ma畸變指數(shù)呈先緩慢增大、再迅速增大、后緩慢增大的趨勢，而其總壓畸變指數(shù)呈緩慢增大的趨勢。45畸變圖譜的Ma畸變指數(shù)變化趨勢與30畸變圖譜相近，但在堵塞比大于56%時存在拐點，而其總壓畸變指數(shù)與30畸變圖譜類似，呈緩慢增大的趨勢，但增速相對較大，且在堵塞比大于56%時也存在拐點。

兩種畸變圖譜時總壓畸變指數(shù)隨堵塞比增大的增速相對于Ma畸變指數(shù)較小，在堵塞比小于等于56%時畸變指數(shù)與堵塞比呈一一對應(yīng)關(guān)系。

2.2不同畸變指數(shù)時總的影響分析

不同畸變指數(shù)（為方便對比，在此采用堵塞比參數(shù)）時補燃室溫升效率和總壓恢復(fù)系數(shù)，如圖6所示，可以看出：相對于基準工況，30畸變圖譜時隨著堵塞比的增大補燃室溫升效率呈現(xiàn)出先小幅減小后增大又減小的趨勢，在堵塞比為56%時高出4.5%左右；45畸變圖譜時隨著堵塞比的增大補燃室溫升效率呈現(xiàn)先減小、后基本持平、再增大、最后又減小的趨勢，在堵塞比為28%和42%時降低3.5%。

兩畸變圖譜時，相對于基準工況，隨著堵塞比的增大，補燃室總壓恢復(fù)系數(shù)呈現(xiàn)出先小幅減小、后快速減小的趨勢，且45畸變圖譜時的減小速率較大，在堵塞比大于56%時存在拐點，相對基準工況最大減小8%左右，可見隨著堵塞比的增大，彎頭內(nèi)局部高速流動增大了總壓損失。

2.3不同畸變圖譜時總的影響分析

30畸變圖譜和45畸變圖譜在不同畸變指數(shù)時補燃室沿程截面溫升效率和總壓恢復(fù)系數(shù)，如圖7和圖8所示，可以看出：畸變指數(shù)不同時，補燃室的沿程溫升效率變化趨勢基本相同，均呈現(xiàn)逐漸升高的趨勢；在堵塞比大于等于42%時，在補燃室上游沿程截面溫升效率增速較高，而在補燃室下游增速較低，在補燃室出口截面溫升效率相對其余工況較高。

畸變指數(shù)不同時，補燃室沿程總壓恢復(fù)系數(shù)變化趨勢基本相同，均呈現(xiàn)先快速減小再緩慢減小的趨勢。

2.4補燃室截面靜溫分布圖

2.4.1 30畸變圖譜

30畸變圖譜在不同堵塞比時補燃室沿程截面靜溫分布圖，如圖9所示，可以看出：在堵塞比小于等于28%時，隨著堵塞比的增大，由于空氣進氣的不均勻性增大，補燃室沿程截面的高溫區(qū)域逐漸向補燃室一側(cè)偏移，但與基準工況差異不大，使得補燃室溫升效率和總壓恢復(fù)系數(shù)均相對于基準工況小幅減小。

在堵塞比大于等于42%時，隨著堵塞比的增大，彎頭內(nèi)空氣高速流動區(qū)域更加集中，在彎頭內(nèi)回流區(qū)域的擠壓下，補燃室沿程截面的高溫區(qū)域逐漸向與堵塞比小于等于28%時相反的補燃室一側(cè)偏移，且由于空氣進氣高Ma區(qū)域相對集中，使得空氣和燃氣相互作用逐漸增強，進而導(dǎo)致補燃室溫升效率逐漸增大，總壓恢復(fù)系數(shù)逐漸減小。

總的來說，對于補燃室沿程截面出現(xiàn)高溫區(qū)域不對稱的情況，是由于在空氣進氣畸變的影響下，空氣入口截面上的Ma分布不均勻，同時在彎頭內(nèi)出現(xiàn)一定的回流區(qū)域，從而導(dǎo)致有一部分空氣先進入補燃室，將燃氣射流壓向一側(cè)。

2.4.2 45畸變圖譜

45畸變圖譜補燃室沿程截面靜溫在不同堵塞比時的分布圖，如圖10所示，可以看出：由于45畸變圖譜相對于30畸變圖譜彎頭內(nèi)高速流動區(qū)域更加偏向于一側(cè)壁面，隨著堵塞比的增大，補燃室沿程截面的高溫區(qū)域逐漸偏向補燃室一側(cè)壁面，空氣和燃氣相互作用逐漸增強，使得補燃室總壓恢復(fù)系數(shù)逐漸降低；此外，高溫區(qū)域不均勻性在一定程度上會減小空氣與燃氣相互作用面積，這兩種作用使得補燃室溫升效率隨堵塞比增大呈現(xiàn)出先小幅減小、后基本持平、再增大，最后再減小的趨勢。

3結(jié)論

通過不同畸變指數(shù)對于補燃室性能的影響分析，可以得到以下結(jié)論：

（1）兩種畸變圖譜時總壓畸變指數(shù)隨堵塞比增大的增速相對于Ma畸變指數(shù)較小，在堵塞比小于等于56%時畸變指數(shù)與堵塞比呈一一對應(yīng)關(guān)系。

（2）兩種畸變圖譜種，隨著堵塞比的增大，補燃室溫升效率呈現(xiàn)出先小幅減小后增大的趨勢，補燃室總壓恢復(fù)系數(shù)呈現(xiàn)出先小幅減小，后快速減小的趨勢。同時，二者均在堵塞比56%時出現(xiàn)拐點，隨著堵塞比的增大，彎頭內(nèi)局部高速流動增強。

（3）兩畸變圖譜時，隨著堵塞比的增大，補燃室的沿程溫升效率和總壓恢復(fù)系數(shù)的變化趨勢基本相同，畸變指數(shù)變化對沿程截面溫升效率和總壓恢復(fù)系數(shù)的影響較小。

（4）兩畸變圖譜時，隨著堵塞比的增大，補燃室沿程截面的高溫區(qū)域逐漸向補燃室一側(cè)偏移，是由于在空氣進氣畸變的影響下，空氣入口截面上的Ma分布不均勻，彎頭內(nèi)出現(xiàn)回流區(qū)域，也會導(dǎo)致部分空氣先進入補燃室將燃氣射流壓向一側(cè)。同時，彎頭內(nèi)高速流動區(qū)域在45畸變圖譜中更加偏向一側(cè)壁面，使得其補燃室沿程截面的高溫區(qū)域?qū)ΨQ性相較30畸變圖譜更差。

參考文獻

[1]曹軍偉，王虎干，蔡選義.整體式固體火箭沖壓發(fā)動機在中遠程空空導(dǎo)彈上的應(yīng)用[J].航空兵器，2002（4）：31-34. Cao Junwei， Wang Hugan， Cai Xuanyi. Application of integral solid rocket ramjet in medium and long range air-to-air missile[J].Aero Weaponry， 2002（4）：31-34.（in Chinese）

[2]陳斌斌，夏智勛，黃利亞，等.含硼固沖補燃室燃燒組織技術(shù)進展[J].航空兵器，2018（4）：3-20. Chen Binbin， Xia Zhixun， Huang Liya， et al. Review on combustion technology of boron-based solid samjet afterburning chamber[J].Aero Weaponry， 2018（4）：3-20.（in Chinese）

[3]楊雄，張志強，郭政波，等.高溫升燃燒室數(shù)值模擬研究[J].航空科學(xué)技術(shù)，2017，28（2）：13-19. Yang Xiong，Zhang Zhiqiang，Guo Zhengbo，et al. Numerical simulation of high temperature rise combustor[J].Aeronautical Science & Technology，2017，28（2）：13-19.（in Chinese）

[4]牛楠，董新剛，霍東興，等.可變流量固沖發(fā)動機一體化流場研究[J].彈箭與制導(dǎo)學(xué)報，2013，33（2）：85-87. Niu Nan， Dong Xingang， Huo Dongxing， et al. Investigation of integrated flow field of variable-ducted rocket[J]. Journal of Projectiles， Rockets， Missiles and Guidance， 2013，33（2）：85-87.（in Chinese）

[5]Kim S，Natan B.Inlet geometry and equivalence ratio effects on combustion in a ducted rocket[J]. Journal of Propulsion and Power，2015，31（2）：619-631.

[6]陳志明，周小旭，張鑫，等.旋流式摻混結(jié)構(gòu)對補燃室二次燃燒性能的影響[J].彈箭與制導(dǎo)學(xué)報， 2017， 37（3）：75-79. Chen Zhiming， Zhou Xiaoxu， Zhang Xin， et al. The effects of swirl mixing structure on the second combustion performance of after-burning chamber[J]. Journal of Projectiles， Rockets， Missiles and Guidance， 2017，37（3）：75-79.（in Chinese）

[7]陳志明，張磊揚，單睿子，等.二維空氣進口結(jié)構(gòu)參數(shù)對固沖補燃室摻混燃燒性能的影響[J].航空兵器，2018（4）：46-51. Chen Zhiming， Zhang Leiyang， Shan Ruizi， et al. Effects of configuration parameters of quadrilateral air entry on the mixing and combustion performance of after-burning chamber[J].Aero Weaponry， 2018（4）： 46-51.（in Chinese）

[8]丁浩，田立豐，郭美琦，等.進氣道流向伸縮對燃燒室入口總壓恢復(fù)的影響[J].航空科學(xué)技術(shù)，2021，32（8）：24-28. Ding Hao， Tian Lifeng， Guo Meiqi， et al. Impact of air inlet flow expansion on total pressure recovery of combustor inlet[J]. Aeronautical Science & Technology， 2021， 32（8）： 24-28.（in Chinese）

[9]王希亮，孫振華，賀永杰，等.頭部兩側(cè)進氣固沖發(fā)動機補燃室內(nèi)流場研究[J].航空兵器，2011（5）：51-55. Wang Xiliang， Sun Zhenhua， He Yongjie， et al. Research on internal flow field in after-burning chamber of ducted rocket with a head bilateral inlet configuration[J]. Aero Weaponry， 2011（5）： 51-55. （in Chinese）

[10]孫海剛，王希亮.針閥流量調(diào)節(jié)裝置對固沖發(fā)動機性能影響的數(shù)值研究[J].航空科學(xué)技術(shù)，2011， 22（3）： 52-54. Sun Haigang， Wang Xiliang. Numerical simulation of solid rocket ramjet performance by using flow regulation vale[J]. Aeronautical Science & Technology，2011， 22（3）： 52-54. （in Chinese）

[11]吳鵬龍.進氣畸變和燃油脈動對燃燒室內(nèi)燃燒流場的影響研究[D].南京：南京航空航天大學(xué)，2012.Wu Penglong. Effects of inlet distortion and fuel pulsation on combustion flow field in combustion chamber[D]. Nanjing： Nanjing University of Aeronautics and Astronautics，2012. （in Chinese）

[12]孫鵬，高海洋，鐘兢軍，等.插板式畸變發(fā)生器后非均勻流場結(jié)構(gòu)數(shù)值模擬[J].推進技術(shù)，2013，34（2）： 173-180. Sun Peng， Gao Haiyang， Zhong Jinjun， et al. Numerical simulation of non-uniform flow field structure behind classic flat baffle[J]. Journal of Propulsion Technology， 2013， 34（2）： 173-180. （in Chinese）

[13]余安遠，王琴芳.一種進氣道流場畸變模擬裝置的研究[J].航空動力學(xué)報，2001，16（2）：123-127. Yu Anyuan， Wang Qinfang. An investigation of distortion simulator for inlet flow field[J]. Journal of Aerospace Power， 2001， 16（2）： 123-127.（in Chinese）

[14]王寅會，張炎，裘云，等.沖壓發(fā)動機燃燒室進氣畸變模擬裝置設(shè)計及試驗[J].實驗流體力學(xué)，2016，30（1）：107-112. Wang Yinhui， Zhang Yan， Qiu Yun， et al. Design and experimental demonstration of distortion simulation board of ramjetcombustor[J].JournalofExperimentsinFluid Mechanics， 2016， 30（1）： 107-112.（in Chinese）

[15]鮑福廷，李進賢，趙飛.固沖發(fā)動機雙下側(cè)二維進氣道設(shè)計研究[J].西北工業(yè)大學(xué)學(xué)報，2002，20（4）：590-593. Bao Fuyan， Li Jinxian， Zhao Fei. A code for designing 2D mixed-compression inlet of solid propellant ramrocket[J]. Journal of Northwestern Polytechnical University， 2002， 20（4）： 590-593. （in Chinese）

Effects of Inlet Distortion Index on the Performance of Post-Combustion in After-Burning Chamber of Solid Ramjet

Zhang Zhe，Chen Zhiming，Zhang Leiyang，Hao Wen China Airborne Missile Academy，Luoyang 471009，China

Abstract： This paper studies the influence of air intake distortion index on the performance of the after-burning chamber of a solid ramjet with double horizontal intake. Due to the flow field structure of outlet section of inlet under different working attitude of ramjet， a corresponding air inlet distortion simulation device is designed. The scheme is designed for the two important factors of the air intake distortion atlas form and the distortion index， and numerical calculation method is used to simulate and analyze the after-burning combustion performance of the combustion chamber with different air intake distortion atlas form and air intake blockage ratio. Results show that relative to the working condition of the benchmark， the total after-burning chamber pressure recovery coefficient slightly decreases first and then rapidly along with the increase of the blockage. Under the same distortion atlas form， the efficiency of temperature rise of after-burning chamber shows a trend of slightly decrease and then increase with the increase of the blockage. When the blockage ratio arrives at 42%， the efficiency of temperature rise is lower than the benchmark condition without distortion.

Received： 2021-10-15；Revised： 2021-11-25；Accepted： 2021-12-25