王希亮，孫振華

(中國空空導彈研究院，河南 洛陽　471099)

?

導彈技術(shù)

頭部兩側(cè)和單側(cè)進氣對固沖發(fā)動機燃燒影響*

王希亮，孫振華

(中國空空導彈研究院，河南 洛陽471099)

摘要：為了研究進氣道布局形式對固沖發(fā)動機二次燃燒性能的影響，分別對頭部兩側(cè)進氣和單側(cè)進氣2種布局方案的固沖發(fā)動機開展了二次燃燒仿真和試驗對比研究。結(jié)果表明：單側(cè)進氣時補燃室摻混燃燒效果優(yōu)于兩側(cè)進氣；單側(cè)進氣會導致補燃室局部出現(xiàn)高溫高速區(qū)域，對補燃室熱防護不利；在相同的工況下單側(cè)進氣獲得的推力與兩側(cè)進氣相當。

關(guān)鍵詞：固沖發(fā)動機；二次燃燒；進氣道布局；熱防護；數(shù)值模擬；試驗

0引言

采用固沖發(fā)動機的導彈為了滿足規(guī)定的戰(zhàn)術(shù)技術(shù)指標要求，必須進行導彈/發(fā)動機/進氣道的一體化優(yōu)化設(shè)計才能實現(xiàn)，由此出現(xiàn)了不同的進氣道布局形式的沖壓發(fā)動機方案。目前國內(nèi)外對各種進氣布局方案的研究主要集中在進氣道特性、導彈/進氣道氣動特性等方面[1-4]。實際上，對固沖發(fā)動機尤其是采用含硼富燃料推進劑的固沖發(fā)動機而言，進氣方式對發(fā)動機的二次燃燒特性也有很大的影響。

受歐洲“流星”空空導彈的影響，國內(nèi)多家單位對采用雙下側(cè)進氣布局方案固沖發(fā)動機二次燃燒性能進行了理論和試驗研究，獲得了大量研究成果[5-12]，同時也有單位對采用中心進氣、環(huán)向進氣和頭部兩側(cè)進氣的固沖發(fā)動機方案開展了研究[13-15]，但對采用單側(cè)進氣道或顎下進氣道形式的固沖發(fā)動機研究尚未見相關(guān)報道。

本文對頭部兩側(cè)進氣和單側(cè)進氣(模擬顎下進氣)2種不同結(jié)構(gòu)形式的固沖發(fā)動機方案開展了二次燃燒仿真和試驗對比研究，其結(jié)果可為固沖發(fā)動機總體設(shè)計提供參考。

1計算模型

1.1物理模型

為便于與試驗結(jié)果進行比較，物理模型采用地面連管試驗狀態(tài)，計算流場范圍包括進氣道、燃氣發(fā)生器(部分)、燃氣通道、沖壓補燃室等部分。為了便于造型與網(wǎng)格劃分，對真實的發(fā)動機結(jié)構(gòu)進行了合理的簡化，由于研究對象是對稱體，流動也具有一定的對稱性，在不考慮迎角和側(cè)滑角組合的情況下，為了減少計算量，對頭部兩側(cè)進氣取整個結(jié)構(gòu)的1/4作為計算域，對單側(cè)進氣取整個結(jié)構(gòu)的1/2作為計算域，其結(jié)構(gòu)如圖1所示。

圖1　發(fā)動機簡化物理模型Fig.1　Simplified physical model of ducted rocket

1.2數(shù)學模型

沖壓發(fā)動機補燃室內(nèi)流場包含非常復(fù)雜的湍流流動和化學反應(yīng)，其流動是三維非定常多組分氣固兩相流動。為了簡化分析，對流場作如下假設(shè)：

(1) 含硼推進劑各組分含量為B(30%)、C4H6O0.03(30%)、AP(40%)；

(2) 補燃室內(nèi)所有氣相為理想氣體，符合氣體狀態(tài)方程；

(3) 補燃室內(nèi)的流動為定常三維流動；

(4) 絕熱層為固定界面，不參與變化和反應(yīng)。

補燃室入口富燃燃氣組分及摩爾分數(shù)通過熱力計算獲得。

根據(jù)以上假設(shè)和流體力學中的質(zhì)量、動量、能量和組分的守恒方程，采用商業(yè)CFD軟件FLUENT進行計算，使用有限體積法求解雷諾平均后的三維N-S方程。紊流模型為RNGκ-ε模型，采用二階迎風格式離散。壁面附近流體計算采用標準壁面函數(shù)處理。

氣相湍流燃燒采用PDF模型，其中熱化學計算在prePDF中進行，并列成表以便于在FLUENT中查詢。紊流和化學的相互作用考慮為一個概率密度函數(shù)(PDF)。

計算收斂準則為：連續(xù)方程、動量方程、能量方程以及κ-ε方程的殘差至少下降3個數(shù)量級，且沖壓噴管出口截面流量穩(wěn)定。

1.3邊界條件

計算中用到的邊界條件類型有：質(zhì)量入口邊界、壓力出口邊界、對稱邊界和無滑移絕熱固壁等，計算模擬狀態(tài)為15 km，3.2Ma，具體邊界條件設(shè)置如表1所示。

表1　邊界條件

2計算結(jié)果與分析

通過仿真獲得了頭部兩側(cè)進氣和單側(cè)進氣2種進氣形式的流場圖譜(見圖2～7所示)。

圖2為2種進氣方式空氣和一次燃氣的流線對比，從中可以看出在單側(cè)進氣時大量空氣直接撞擊在進氣道對側(cè)，然后沿兩側(cè)壁面回卷，在補燃室內(nèi)形成了明顯的螺旋狀流態(tài)，一次燃氣受空氣流的阻擋，在補燃室頭部形成了比較大的回流區(qū)，在補燃室后部一次燃氣與空氣的摻混較好。而在兩側(cè)進氣時大部分一次燃氣都被空氣流包裹在中心區(qū)域，頭部只有很小的回流區(qū)，在整個補燃室內(nèi)一次燃氣與空氣摻混較少。

圖2　2種進氣方式流線對比Fig.2　Comparison of air inlet streamline of two modes of air inlet

圖3和圖4分別是2種進氣方式對稱面和壁面的靜溫分布云圖，從圖3可以看出單側(cè)進氣頭部存在較大的高溫區(qū)域，而兩側(cè)進氣的高溫區(qū)域主要集中在中間的帶狀區(qū)域內(nèi)。從圖4可以看出單側(cè)進氣在進氣道對側(cè)存在大面積的高溫區(qū)域，最高溫度在2 500 K左右，而兩側(cè)進氣時補燃室壁面高溫區(qū)主要集中在兩進氣道所夾的區(qū)域內(nèi)，且最高溫度只有2 100 K左右，相比而言兩側(cè)進氣有利于補燃室的熱防護設(shè)計。

圖3　2種進氣方式對稱面靜溫云圖對比Fig.3　Comparison of static temperature contours on symmetric for two modes of air inlet

圖4　2種進氣方式壁面靜溫云圖對比Fig.4　Comparison of static temperature contours on wall for two modes of air inlet

圖5為2種進氣方式對稱面馬赫數(shù)云圖。從中可以看出由于假設(shè)2種進氣方式捕獲空氣流量相同，單側(cè)進氣進氣道出口的氣流馬赫數(shù)較高，達到了0.8～0.9，略高于一次燃氣的出口馬赫數(shù)，因此單側(cè)進氣發(fā)動機的空燃動量比相對較高，空氣與燃氣接觸時的沖擊增強，強化了燃氣摻混，有利于二次燃燒，同時由于燃氣流無法直接穿透空氣流，在補燃室頭部形成了回流區(qū)，有利于硼粒子的點火。兩側(cè)進氣進氣道出口的氣流馬赫數(shù)為0.4～0.5，低于一次燃氣的出口馬赫數(shù)，因此一次燃氣能夠穿透空氣流，并在補燃室中心區(qū)域取得了較好的摻混效果。在補燃室的大部分區(qū)域，單側(cè)進氣時進氣道對側(cè)的氣流馬赫數(shù)相對較高，而兩側(cè)進氣馬赫數(shù)分布較為均勻。

圖5　2種進氣方式對稱面馬赫數(shù)云圖對比Fig.5　Comparison of Mach number contours on symmetric for two modes of air inlet

圖6和圖7分別是2種進氣方式對稱面和補燃室不同截面的靜壓分布云圖，從中可以看出單側(cè)進氣補燃室頭部靜壓差異較大，而兩側(cè)進氣補燃室頭部靜壓分布較為均勻，這與前面的流動分析結(jié)果相對應(yīng)。在補燃室后部，單側(cè)進氣的高壓區(qū)域明顯大于兩側(cè)進氣，說明單側(cè)進氣的二次燃燒效率較高。從圖7可以看出，在補燃室出口截面單側(cè)進氣的靜壓分布不均勻性要高于兩側(cè)進氣。

圖6　2種進氣方式對稱面靜壓云圖對比Fig.6　Comparison of static pressure contours on symmetric for two modes of air inlet

圖7　2種進氣方式補燃室不同截面靜壓云圖對比Fig.7　Comparison of static pressure contours on different sections for two modes of air inlet

3試驗研究

3.1試驗方案

針對單側(cè)進氣和兩側(cè)進氣開展了對比試驗研究，兩次試驗使用的發(fā)動機除進氣方式不同外其余結(jié)構(gòu)完全相同。在模擬單側(cè)進氣時將一側(cè)進氣道去除，并在相應(yīng)的進氣道出口安裝堵蓋密封。2次試驗?zāi)M空氣來流狀態(tài)相同，并采用同一批次的推進劑以保證一次燃氣的性質(zhì)和流量一致。圖8為2次試驗采用發(fā)動機結(jié)構(gòu)示意圖。

圖8　2次試驗發(fā)動機結(jié)構(gòu)示意圖Fig.8　Schematic diagram of the two test engine structures

3.2試驗結(jié)果分析

根據(jù)試驗方案，完成了2次點火試驗，試驗過程中發(fā)動機工作正常，試驗數(shù)據(jù)有效。圖9為2次試驗發(fā)動機工作過程中燃氣流量與推力變化曲線，其中推力按最大值進行了歸一化處理。從中可以看出2次試驗燃氣流量大小基本相同，變化規(guī)律一致，推力曲線波動幅度相對較大，總體而言單側(cè)進氣的平均推力要略高于兩側(cè)進氣的平均推力。

圖9　2次試驗燃氣流量與推力變化曲線Fig.9　Curve of gas flow and the thrust in two tests

對試驗結(jié)果進行了處理，其結(jié)果見表2。其中試驗值按工作時間內(nèi)的平均值計算，燃燒效率是指溫升效率，用補燃室尾部靜壓數(shù)據(jù)經(jīng)處理獲得，并以兩側(cè)進氣的結(jié)果為基準進行歸一化處理。

表2　主要性能參數(shù)處理結(jié)果

從表2可以看出進氣方式對補燃室前后壓強分布有較大的影響，單側(cè)進氣時補燃室靜壓頭部低、尾部高，而兩側(cè)進氣時補燃室靜壓頭部高、尾部低。單側(cè)進氣的燃燒效率高于兩側(cè)進氣，試驗工況下按溫升效率計算的燃燒效率高14.7%，與前文的數(shù)值仿真結(jié)論一致。從臺架推力可以看出，2種進氣方式獲得的推力基本一致，說明雖然單側(cè)進氣能獲得較高的燃燒效率，但由于補燃室內(nèi)流場分布不均勻，總壓損失較大，導致最終得到的推力并沒有明顯增加。

2次試驗均在尾部截面2個相同的測點位置進行了總壓測試(見圖10)，測點位置距內(nèi)壁面高度分別為90 mm和60 mm，總壓測量結(jié)果見表3。從表中可以發(fā)現(xiàn)兩者有較大的區(qū)別。雙側(cè)進氣補燃室尾部2個測點的總壓在點火前后均基本相同，單側(cè)進氣時2個測點的總壓在點火前后存在明顯差異，說明該截面流場分布極不均勻。

圖10　補燃室內(nèi)總壓測點布置圖Fig.10　Afterburning chamber total pressure measuring point arrangement

表3　補燃室尾部總壓測量結(jié)果

備注：測點1插入深度90 mm，測點2插入深度60 mm

圖11為仿真得到的補燃室尾部截面的總壓分布云圖(壓強基準值為100 000 Pa)，從中可以看出，雙側(cè)進氣時該截面的壓強分布相對較為均勻，最高壓強與最低壓強相差約10%，相反單側(cè)進氣時該截面的壓強分布差異明顯，最高壓強與最低壓強相差高達37%左右，且單側(cè)進氣該截面的平均總壓要略高于雙側(cè)進氣相同截面的平均總壓。

圖11　仿真得到的補燃室尾部測壓截面總壓云圖Fig.11　Simulation result of the afterburner tail pressure section total pressure contours

4結(jié)論

(1) 在設(shè)計的工作條件下，單側(cè)進氣比雙側(cè)進氣的二次燃燒性能要高，以補燃室尾部靜壓為基準，按溫升效率計算單側(cè)進氣比雙側(cè)進氣的燃燒效率高14.7%；

(2) 單側(cè)進氣使得補燃室絕熱層局部承受高溫、高速燃氣的沖刷，在熱防護設(shè)計時必須加以考慮；

(3) 單側(cè)進氣雖然能獲得較高的燃燒效率，但由于補燃室內(nèi)流場分布不均勻，總壓損失較大，導致最終得到的推力并沒有明顯增加。

參考文獻：

[1]蔣妮，王永衛(wèi)，凌文輝，等. 典型布局的進氣道/導彈一體化性能分析[J]. 推進技術(shù), 2011, 32(4):451-454.

JIANG Ni, WANG Yong-wei, LING Wen-hui，et al. Integration Performance Analysis of Typical Arrangement Inlet/Missile[J]. Journal of Propulsion Technology, 2011, 32(4):451-454.

[2]孟宇鵬，閆曉娜，朱守梅，等. 旁側(cè)四超聲速進氣道彈體內(nèi)外流一體化數(shù)值研究[J]. 戰(zhàn)術(shù)導彈技術(shù), 2008(5):24-30.

MENG Yu-peng, YAN Xiao-na, ZHU Shou-mei, et al. Numerical Investigation into Integrated External Flow Around Missile Body and Internal Flow for Four Supersonic Inlets[J]. Tactical Missile Technology, 2008(5):24-30.

[3]李雪平, 楊茂. 進氣道布局對導彈氣動性能影響的數(shù)值研究[J].科學技術(shù)與工程, 2010, 10(35):8747-8751.

LI Xue-ping, Yang Mao. Numerical Study of Inlet Configuration Effects on Aerodynamic Performance for Air-Breathing Missiles[J]. Science Technology and Engineering, 2010, 10(35): 8747-8751.

[4]謝文忠, 郭榮偉. 4種布局形式下超聲速飛行器進氣道氣動特性實驗對比[J]. 南京航空航天大學學報, 2011, 43(1):13-17.

XIE Wen-zhong, GUO Rong-wei. Mixed-Compression Supersonic Inlets Based on Four Air-Breathing Aircraft Configurations [J]. Journal of Nanjing University of Aeronautics & Astronautics, 2011, 43(1):13-17.

[5]STOWER A. Performance Prediction of a Ducted Rocket Combustor [D].University of Laval, 2001.

[6]張磊，周長省，鞠玉濤，等. 燃氣側(cè)噴固沖發(fā)動機補燃室流場三維數(shù)值研究[J].計算機仿真, 2008, 25(5):61-64.

ZHANG Lei, ZHOU Chang-sheng, JU Yu-tao, et al. 3D Numerical Simulation of Flow Field in Secondary Combustion Chamber of Fuel Gas Dual-Injection Ramjet Engine[J]. Computer Simulation, 2008, 25(5):61-64.

[7]VANKA S P.Analytical Characterization of Flow Fields in Side Inlet Dump Combustors[R].AIAA 83-1399.

[8]RISTORI A，DUFOUR E. Numerical Simulation of Ducted Rocket Motor[R].AIAA 2001-3193.

[9]李澤勇，胡建新，夏智勛，等. 進氣道角度對含硼推進劑固沖發(fā)動機性能的影響[J]. 國防科技大學學報, 2008, 30(2):1-4.

LI Ze-yong, HU Jian-xin, XIA Zhi-xun, et al. Effects of Side-Inlet Angle on the Performance of Boron-Based Propellant Ducted Rocket[J]. Journal of National University of Defense Technology, 2008, 30(2):1-4.

[10]霍東興，何國強，陳林泉，等. 固沖發(fā)動機補燃室冷流摻混效果與燃燒效率對比研究[J]. 固體火箭技術(shù), 2006, 29(5):329-332.

HUO Dong-xing, HE Guo-qiang, CHEN Lin-quan, et al. Contrast Study on Cold-Flow Mixing Effect and Combustion Efficiency of Afterburning Chamber of Solid Rocket Ramjet[J]. Journal of Solid Rocket Technology, 2006, 29(5):329-332.

[11]萬少文，何國強. 摻混裝置對固沖發(fā)動機摻混燃燒性能的影響[J]. 彈箭與制導學報, 2011, 31(4):108-109.

WAN Shao-wen, HE Guo-qiang. The Effect of Mixing-Enhanced Device on the Combustion Performance of Solid Rocket-Ramjet Motor[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2011, 31(4):108-109.

[12]李綱，何國強，孫振華，等. 固沖發(fā)動機補燃室二次燃燒實驗研究[J]. 固體火箭技術(shù), 2007, 30(5):400-403.

LI Gang, HE Guo-qiang, SUN Zhen-hua, et al. Experimental Investigation on Secondary Combustion in Solid Rocket Ramjet[J]. Journal of Solid Rocket Technology, 2007, 30(5):400-403.

[13]鄭凱斌，陳林泉，張勝勇. 中心進氣式固體火箭沖壓發(fā)動機試驗研究[J]. 固體火箭技術(shù), 2007, 30(2):124-127.

ZHENG Kai-bin, CHEN Lin-quan, ZHANG Sheng-yong. Experimental Study on Solid Rocket Ramjet with Central Air-Inlet Mode[J]. Journal of Solid Rocket Technology, 2007, 30(2):124-127.

[14]趙春宇，李斌，鞠玉濤. 環(huán)向進氣固體火箭沖壓發(fā)動機補燃室流場數(shù)值模擬[J]. 彈箭與制導學報, 2008, 28(2):136-138.

ZHAO Chun-yu, LI Bin, JU Yu-tao. Numerical Simulation of Flow Field in Afterburning Chamber for Ringed-propellant Solid Rocket Ramjet Engine[J]. Journal of Projectiles, Rockets, Missiles and Guidance,2008, 28(2):136-138.

[15]王希亮，孫振華，賀永杰，等. 頭部兩側(cè)進氣固沖發(fā)動機補燃室內(nèi)流場研究[J]. 航空兵器, 2011(5):51-55.

WANG Xi-liang, SUN Zhen-hua, HE Yong-jie, et al. Research on Internal Flowfield in Afterburning Chamber of Ducted Rocket with a Head Bilateral Inlet Configuration[J]. Aero Weaponry, 2011(5):51-55.

Effect of Head Bilateral or Lateral Inlet Configuration on Combustion Performance of Ducted Rocket

WANG Xi-liang, SUN Zhen-hua

(China Airborne Missile Academy,Henan Luoyang 471099, China)

Abstract:In order to obtain the effect of inlet configuration on the performance of after burning of ducted rocket, both numerical simulation and experiment are performed on two different inlet configurations, with a head bilateral inlet configuration and a head lateral inlet configuration respectively. The results show that the performance of after burning of ducted rocket with a head lateral inlet configuration is better than that with a head bilateral inlet configuration. However, high temperature and high velocity appeared locally in head lateral inlet configuration, which is unfavorable to after burner insulator. Under the same condition, the thrust for two different configurations are almost the same.

Key words:ducted rocket; after burning; inlet configuration; thermal protection; numerical simulation; experiment

*收稿日期：2015-03-03；修回日期：2015-07-07

作者簡介：王希亮(1975-)，男，湖南益陽人。高工，碩士，研究方向為固體火箭沖壓發(fā)動機設(shè)計。

通信地址：471099河南省洛陽市030信箱4分箱E-mail：wangxl04764@163.com

doi:10.3969/j.issn.1009-086x.2016.02.011

中圖分類號：V435

文獻標志碼：A

文章編號：1009-086X(2016)-02-0068-06