Effect of the Nozzle Structure on Combustion Performance of Nose Fuel-inlet Model Solid-fuel Rocket Scramjet

  • LIU Zai ,
  • LI Xuhua ,
  • WANG Liwu ,
  • CHEN Linquan
Expand
  • The Fourth Academy of the China Aerospace Science and Technology Corporation, Xi’an 710025, China

Received date: 2022-04-12

  Online published: 2025-05-29

Abstract

In order to improve the combustion performance of head intake solid rocket scramjet, the effect of rich fuel gas nozzle structure on the engine combustion performance was researched by using numerical simulation. The result showed that compared with the traditional round section nozzle structure, the elliptical section nozzle structure has more sufficient mixing of rich fuel gas and air, higher chemical reaction rate, higher thrust, and specific impulse gain. With the ratio increase of long half axis to short half axis of elliptical section nozzle, the combustion performance of secondary combustion chamber is better. It is recommended to use elliptical section nozzle in engineering.

Cite this article

LIU Zai , LI Xuhua , WANG Liwu , CHEN Linquan . Effect of the Nozzle Structure on Combustion Performance of Nose Fuel-inlet Model Solid-fuel Rocket Scramjet[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2022 , 42(4) : 47 -50 . DOI: 10.15892/j.cnki.djzdxb.2022.04.009

0 引言

固体火箭超燃冲压发动机具有密度比冲高、结构简单、作战反应迅速等优势,因此在高超声速导弹武器中具有较高的应用价值。头部进气固体火箭超燃冲压发动机是将燃气发生器安装在进气道中心锥内,与来流超声速空气流平行进入补燃室发生掺混燃烧。吕仲[1]率先针对该种发动机开展了可行性实验验证及仿真分析,获得了推力及比冲增益。刘仔等[2-3]研究了补燃室长度及扩张比、来流条件等因素对燃烧性能的影响,在采用圆形截面喷孔且无任何增强掺混燃烧措施的条件下,发动机燃烧效率较低,不具备工程化运用条件。因此,如何提高固体火箭超燃冲压发动机的燃烧性能是进入工程化运用前的重点研究工作。
李轩等[4]采用数值模拟方法研究了凹腔、扰流装置对头部进气固体火箭超燃冲压发动机燃烧性能的影响,研究发现扰流装置在提升燃烧性能方面优于凹腔。黄礼铿等[5]采用地面直连实验手段研究了一次燃气喷注位置、喷注结构等参数对掺混燃烧性能的影响,获得了较优的燃气喷注总体方案。朱韶华等[6]针对燃气喷注方式、扰流装置及燃烧室扩张比等因素,基于数值模拟方法研究相关因素对燃烧性能的影响。凌江等[7]研究了燃气进气方向与轴线之间夹角对补燃室燃烧性能的影响,获得了补燃室燃烧性能随夹角的变化规律,优化了补燃室的进气角度。李潮隆[8]基于数值仿真研究了圆形截面与方形截面两种喷孔结构对补燃室性能的影响。
针对固体火箭超燃冲压发动机增强掺混燃烧的研究工作还较少。借鉴固体火箭冲压发动机中非圆形截面喷孔可有效增强补燃室内亚声速空气流与富燃燃气的掺混燃烧[9-12]的成果,研究喷孔结构对掺混燃烧的影响,为后续富燃燃气的喷孔结构设计提供参考。

1 物理模型

补燃室构型采用三段扩张设计方案,富燃燃气采用头部进气模式;进气道采用三维轴对称进气道形式,即超声速空气从环形通道进入补燃室内。其中补燃室的结构如图1所示。3段燃烧室长度分别为0.1 m,0.2 m与0.4 m,扩张半角分别是0°,1°与2°。空气入口外径为0.05 m,内径为0.045 m。
喷孔安装在补燃室轴线上,并采用单喷孔结构设计,喷孔结构如图2所示。为保证计算结果具有可比性,所有模型中喷孔的横截面积相等,且保证富燃燃气与空气的进气参数一致。研究圆形喷孔与椭圆形喷孔(长半轴a与短半轴b)对补燃室掺混燃烧性能的影响,其中,ab的比值取1,2和4三种情况,ab取值如表1所示。
图2 喷孔截面形状
表1 喷孔的结构参数单位:mm
参数 圆形喷孔 椭圆形喷孔
a=2b a=4b
a 5.0 7.1 10.0
b 5.0 3.5 2.5

2 数值计算方法

2.1 计算方法

采用稳态N-S方程描述空气与富燃燃气的掺混燃烧过程。假设燃气满足理想气体状态方程,不考虑气体的辐射作用,且忽略燃气重力的影响效应。
考虑到超声速状态下气体具有明显的可压缩性,在基于密度条件下,利用二阶迎风格式离散对流项,Roe-FDS求解界面通量。采用自适应SST k-ω湍流模型,结合增强壁面函数法可有效预测反压作用下的边界层流动及分离,并保证补燃室壁面的第一层网格高度为2。为考虑湍流对化学反应过程的影响效应,采用广泛使用的涡团耗散模型描述补燃室内燃烧过程,其模型计算时均化学反应速率的公式见文献[13]。
为节省计算量取1/4计算域进行对称计算,采用六面体结构网格,并对壁面区域网格进行加密处理,确保壁面第一层网格满足计算要求,总网格数为75万,如图3所示。
计算边界条件如下:
1) 喷孔入口
采用压力入口边界条件。燃气发生器设计方案与文献[4]保持一致。其中,燃气发生器室温2 200 K,组分为50%的C2H4、25%的CO2与25%的H2O,燃气发生器喷孔出口Ma=1。简化富燃燃气与空气的化学反应过程,采用C2H4与O2的单步不可逆反应描述燃烧过程。
2) 空气入口
采用压力入口边界条件。空气组分简化N2为73.6%、O2为23.2%、H2O为3.2%。入口空气的马赫数1.6,总压为1.61 MPa,总温为1 200 K。
3) 出口、壁面与对称面
由于补燃室内的燃气流动是超声速流动,因此补燃室出口的流场参数采用外推方式得到。壁面选择绝热、无滑移及零压力梯度壁面边界条件。对称面采用对称边界条件,即流场内的各参数在对称面上的梯度均为零。

2.2 算例验证

对Evans等[14]的超声速轴对称平行进气扩散燃烧实验进行仿真,验证仿真分析模型的准确性,仿真计算结果如图4所示。计算得到的组分分布规律与实验吻合,同一监测点组分的计算结果与试验结果最大偏差不超过10%。因此采用的计算模型与计算方法可用于补燃室内的超声速掺混燃烧过程的仿真。
图4 距离燃烧室入口78.6 mm横截面上组分质量分数的径向分布

3 计算结果与分析

图5为补燃室截面的静压分布云图。不同喷孔结构的静压分布差异较大。圆形喷孔的静压分布只在径向存在差异,椭圆形喷孔的静压分布在径向与周向均存在差异,表明椭圆形喷孔对流场的扰动在补燃室径向与周向均存在,具有增强掺混燃烧的潜力。采用椭圆形喷孔时补燃室内的静压明显高于圆形喷孔,且随着喷孔长半轴的增加而增加。分析认为补燃室静压增大是C2H4与空气掺混燃烧更充分所致。
图5 补燃室静压分布云图
图6为补燃室C2H4质量分数分布,采用不同喷孔结构方案补燃室头部均存在低速回流区,回流区的存在能够有效增强掺混燃烧,并起到点火与火焰稳定的作用。采用椭圆形喷孔时,补燃室头部低速回流区环向的非对称性,使得C2H4的分布更分散,且分散程度随着椭圆长半轴的增加而增大,有效提高C2H4与空气的接触面积,对增强掺混燃烧是有利的。同时,随着椭圆形喷孔长半轴的增加,C2H4流线的扭曲程度逐渐增强,有利于C2H4与空气的掺混。采用圆形喷孔时,补燃室头部低速回流区环向的对称性,导致C2H4仅集中在补燃室中心区,且流线没有发生明显扭曲,不利于富燃燃气与空气掺混燃烧。随着椭球比(a/b)的增大补燃室出口C2H4的质量分数越小,表明C2H4与空气的掺混燃烧更加充分。
图6 C2H4的质量分数分布
图7为补燃室截面平均化学反应速率的轴向分布。
图7 平均化学反应速率轴向分布
不同喷孔结构条件下补燃室截面平均化学反应速率沿轴向的分布规律基本一致,最大化学反应速率均出现在补燃室头部区域。采用椭圆形喷孔的平均化学反应速率明显高于圆形喷孔,表明椭圆形喷孔能够有效促进C2H4与空气的掺混燃烧,与图5图6结论一致。随着椭圆形喷孔长半轴增加,平均化学反应速率也增加。
燃烧效率作为衡量燃料化学能转化为热能的程度,是反映发动机燃烧性能的重要指标。其定义为是已燃烧燃料的质量占燃料总质量的百分比:
η = 1 - Y F ρ u d A m · F , i
式中:YF为燃气的质量分数;ρ为燃气密度;u为燃气的轴向速度; m · F , i为补燃室入口燃气的质量流量;A为补燃室横截面面积。
总压损失主要反映燃气的做功能力,总压损失越大,表明燃气的做功能力越弱。燃烧效率计算为:
σ = p t , e ρ u d A / m · e p t , i ρ u d A / m · i
式中:pt,e为补燃室出口燃气总压; m · e为补燃室出口燃气质量流量;pt,i为补燃室入口燃气总压; m · i为补燃室入口燃气质量流量。
使用补燃室推力及比冲增益来衡量燃烧性能。基于动量定理的原理,补燃室推力增益为:
ΔF= m · eue- m · iui+peAe-piAi
式中:ue为燃气出口轴向速度;ui为燃气入口轴向速度;pe为燃气出口静压;pi为燃气入口静压;Ae为补燃室出口横截面面积;Ai为补燃室入口横截面面积。
补燃室比冲增益为:

ΔIpF/ m · F

式中: m · F为补燃室燃气的质量流量。
表2为不同喷孔结构下补燃室总体性能。采用椭圆形喷孔时补燃室的燃烧效率明显高于圆形喷孔,且燃烧效率随着椭圆形喷孔长半轴的增加而提高。由于燃烧效率的大幅增加,使得燃烧加热带来的总压损失明显提高,因此总压恢复系数与燃烧效率的变化规律相反。补燃室推力及比冲增益随着燃烧效率的增大而增大。因此,椭圆形喷孔能有效提高补燃室的性能,后续可采用椭圆形喷孔代替圆形喷孔。
表2 补燃室性能参数
喷孔截
面参数
燃烧效
σ/%
总压恢复
系数η/%
推力增
益ΔF/N
比冲增益
ΔIp/s
a=b 31.64 46.83 167.1 568.6
a=2b 64.05 30.25 186.8 635.7
a=4b 71.13 29.25 198.4 675.1

4 结论

针对头部进气固体火箭超燃冲压发动机,研究了喷孔结构对发动机燃烧性能的影响,结论如下:
1) 补燃室头部均存在低速回流区,且椭圆形喷孔的低速回流区环向不对称,对整个补燃室流场形成扰动,使得流场参数在径向与周向均存在明显差异。
2) 椭圆形喷孔头部回流区环向的不对称效应,促进了C2H4在补燃室内与空气掺混,提高了C2H4与空气的掺混效率。
3) 随着长半轴与短半轴比值(a/b)的增大,补燃室燃烧效率、推力及比冲增益逐渐增大,总压恢复系数逐渐减小,后续可采用椭圆形喷孔取代圆形喷孔。
[1]
吕仲. 固体火箭超燃冲压发动机工作特性研究[D]. 长沙: 国防科技大学, 2012.

[2]
刘仔, 陈林泉, 吴秋, 等. 固体火箭超燃冲压发动机补燃室构型的影响分析[J]. 固体火箭技术, 2017, 40(4): 432-436.

[3]
刘仔, 陈林泉, 吴秋. 空燃比对固体火箭超燃冲压发动机性能的影响[J]. 弹箭与制导学报, 2017, 37(6):93-95.

DOI

[4]
李轩, 马利锋, 赵永涛, 等. 固体火箭超燃冲压发动机性能数值模拟研究[J]. 弹箭与制导学报, 2014, 34(1): 104-108.

DOI

[5]
黄礼铿, 胡广军, 胡豹, 等. 固体火箭超燃冲压发动机燃烧试验研究[J]. 固体火箭技术, 2020, 43(5): 549-553.

[6]
朱韶华, 梁磊, 秦飞, 等. 固体火箭超燃冲压发动机燃烧性能影响因素研究[J]. 推进技术, 2021, 42(3): 638-646.

[7]
凌江, 徐义华, 孙海俊, 等. 燃气喷射角度对含硼固体火箭超燃冲压发动机补燃室燃烧效率的影响[J]. 火箭推进, 2022, 48(1):69-75.

[8]
李潮隆. 固体火箭超燃冲压发动机燃烧组织技术研究[D]. 长沙: 国防科技大学, 2019.

[9]
王同辉, 白涛涛, 莫展, 等. 特型燃气喷口对补燃室掺混燃烧的影响[J]. 弹箭与制导学报, 2015, 35(2): 97-100.

DOI

[10]
陈林泉, 毛根旺, 霍东兴, 等. 燃气喷射方式对冲压发动机补燃室掺混效果的影响[J]. 固体火箭技术, 2005, 28(1): 40-43.

[11]
胡建新. 含硼推进剂固体火箭冲压发动机补燃室工作过程研究[D]. 长沙: 国防科技大学, 2006.

[12]
GUTMARK E, SHADOW K C, WILSON K J. Subsonic and supersonic combustion using noncircular injectors[J]. Journal of Propulsion and Power, 1991, 7(2): 240-249.

[13]
陶欢. 固体燃料超燃冲压发动机燃烧室工作特性研究[D]. 北京: 北京理工大学, 2015.

[14]
EVANS J S, SCHEXNAYDER J C, BEACH H L. Application of a two-dimensional parabolic computer program to prediction of turbulent reacting flows: NASA-TP-1169[R]. Washington: NASA, 1978.

Outlines

/