激波诱导推力矢量喷管不同气体喷注时的性能分析

2018-11-28王晓明刘辉韩龙柱袁修干

北京航空航天大学学报 2018年11期

王晓明，刘辉，韩龙柱，袁修干

(1. 北京航空航天大学航空科学与工程学院, 北京 100083； 2. 中国航发商用航空发动机有限责任公司, 上海 200241； 3. 北京航空航天大学生物与医学工程学院，北京 100083； 4. 北京航空航天大学生物医学工程高精尖创新中心，北京 100083)

现代战争对飞行器的机动性要求越来越高。而推力矢量控制技术的应用可有效提高飞行器的机动性能，减小飞行阻力，提高隐身及生存能力，因此越发受到人们的重视。而传统的机械式推力矢量控制技术系统复杂，质量大，响应时间长，基于减重增效的考虑，研究更高效的推力矢量发动机已成为未来飞行器性能提高的关键技术之一[1-3]。激波诱导推力矢量控制技术是在喷管扩张段注入二次流体形成激波，改变喷管出口气流的方向，来实现推力矢量控制。该技术系统结构简单，响应时间极短，是一种理想的技术方案[4-7]。

国外的激波诱导推力矢量控制技术研究最初始于20世纪60年代，主要针对火箭发动机开展了一系列燃气和液体二次喷射矢量控制的技术研究，其中液体二次喷射激波诱导推力矢量控制方案在一些战术导弹上已成功应用[8]。从20世纪90年代开始，美国的NASA Langley实验室开展了针对航空发动机的流体控制推力矢量控制技术的研究工作，取得了较大进展[5]。而国内相关领域研究开展于20世纪90年代末，理论和试验研究资料相对匮乏，只局限于稳态条件下次流喷射角度、喷射位置等几何因素对推力矢量性能影响规律的分析研究。

激波诱导推力矢量喷管的次流气源可以从发动机燃烧室引气，也可以由单独的气源供气。因此，根据气源的不同，次流的种类也不尽相同。而现有研究主要针对主流和次流为同种气体的情况，本文则主要研究了不同工质次流气体条件下喷管的推力矢量性能。通过对喷管内流场进行数值模拟，分析比较次流的渗透高度、喷管的矢量偏转角和推力系数等参数，研究次流的分子质量对推力矢量性能的影响规律，从而选择合适的次流气体，提高矢量喷管的性能。

1 数值方法

喷管内流场为超声速、可压缩气体，密度脉动不能忽略，因此Navier-Stokes方程中的湍流脉动项采用Favre提出的质量加权平均[9]。为使Favre平均Navier-Stokes方程封闭，Favre平均得出的雷诺扩散项和雷诺导热项由梯度扩散假设封闭，其中湍流施密特数(Sc)和普朗特数(Pr)均设为0.9。雷诺应力项由Boussinesq假设封闭，湍流黏性由k-ωSST(Shear Stress Transport)湍流模型方程获得[10]。

依据上述说明，笛卡儿坐标系下多组分守恒型控制方程的形式如下：

(1)

气体守恒方程为

(2)

式中：p、ρ和T分别为气体压力、密度和温度；R为理想气体常数；Ms为各组分分子质量。

控制方程的无黏通量向量采用AUSM (Advection Upstream Splitting Method)格式的改进型ASUM+格式进行离散[11]：

(3)

式中：各变量符号说明参见文献[11]。

通过对原始变量进行MUSCL (Monotone Upstream-centred Schemes for Conservation Laws)差值，可提高到二阶精度[12]。黏性通量向量采用中心差分离散。时间推进采用LU-SGS (Low-Upper Symmetric Gauss-Seidel)隐式格式[13]。

2 程序验证及算例

2.1 物理模型

将文献[14]中的氮气垂直喷入超声速空气流场的试验结果来验证程序的准确性。其试验台为500 mm×150 mm×150 mm的三维槽道，其结构平面图如图1所示。计算网格为318×69×69，在喷射口和壁面处局部加密。表1给出了主流场中的空气参数和喷射入口处的氮气参数[14]。

图1 氮气喷射试验台几何构型Fig.1 Geometry configuration of nitrogen injection testbed

参数空气流氮气流Ma3.711.0p/MPa1.2360.606T/ K301301YN20.76641YO20.23360

注：Ma—马赫数；YN2—氮气组分质量分数；YO2—氧气组分质量分数。

2.2 计算结果

图2为马赫数等值线，图3为喷射口附近区域的流线，x、y分别为空气流入方向和氮气注入方向。氮气从喷射口欠膨胀喷出，超声速气流通过Prandtl-Meyer膨胀在喷射口外形成马赫盘。由于喷射气流与主流的交叉干扰作用，在喷射口上游区域形成弓形激波。弓形激波前流场压力升高，致使壁面边界层分离，进而形成了一个与弓形激波相交的分离激波。另外，在喷射口下游区域形成再压缩激波。如图2和图3所示，干扰流场复杂的激波结构，马赫盘和分离涡系结构都被数值方法很好的捕捉。

图4为计算所得壁面压力分布与试验数据结果对比，p∞为空气来流总压，L为氮气喷射孔中心线到入口的距离(即330 mm)。由图4可见，计算结果与试验数据吻合得很好，只是在喷射口下游计算所得压力的峰值比试验数据略低。

由上述分析可知，AUSM+格式和k-ωSST湍流模型能够很好地模拟横向喷流的流场结构，计算结果也与试验数据吻合得很好，验证了数值方法的可靠性，因此该方法可以应用到激波诱导推力矢量喷管不同气体喷注时的性能分析研究。

图2 马赫数等值线Fig.2 Mach number contour lines

图3 喷射口附近区域流线Fig.3 Streamlines near injection orifice

图4 壁面压力分布Fig.4 Wall pressure distribution

3 不同气体喷注时性能分析

3.1 物理模型

喷管全长为78 mm，喉部直径47 mm，喷管出口面积与喉部面积之比为2.09。由于在相同条件下，矩形注气口的推力矢量性能优于同面积的圆形注气口[15]，因此本模型在喷管扩张段开设周向角度40°，宽度3 mm的矩形注气口。该注气口的中心线距喷管入口的距离为喷管全长的0.75，次流方向与喷管轴线夹角为θ。由于几何和流场结构的对称性，沿对称面选取一半作为研究对象。计算模型的网格如图5所示。

喷管主流气体为空气，其中氮气的质量分数为0.766 4，氧气的质量分数为0.233 6，入口总压为1 MPa和1.5 MPa，出口环境压强为一个标准大气压，则主流落压比(NPR)分别为10和15。次流气体分别为He、N2和CO2气体，其质量分数均为1，次流与主流总压之比(SPR)分别为0.8和1.0，马赫数为1，注气角度θ分别为90°和105°，垂直和逆向主流喷入。主流和次流的总温均为300 K。

图5 激波诱导推力矢量喷管计算网格Fig.5 Computational grid of shock thrust vector nozzle

3.2 评价标准

通过数值计算可得到喷管矢量偏转角和推力系数作为推力矢量性能指标[16]。

矢量偏转角度为

α=arctan(FN/FA)

(4)

式中：FN和FA分别为轴向力和侧向力，可对喷管出口气流积分得到。

推力系数为

CR,i=FR/(Fi,p+Fi,s)

(5)

次流与主流质量流量比为

ω=Ws/Wp

(6)

式中：Ws和Wp分别为次流和主流的质量流量，可分别对喷管入口和次流入口积分得到。

3.3 计算结果

图6～图8所示的状态:次流为He，NPR=10，SPR=1.0，θ=90°。图6为He注入喷管流场后压强分布。如图6所示，激波诱导推力矢量喷管的流场结构与第2节中验证模型的流场结构相似，二次流气体喷入喷管扩张段，通过Prandtl-Meyer膨胀在喷射口外形成马赫盘，由于次流的阻碍作用，在喷射口上游形成弓形激波，主流通过弓形激波流动方向发生偏转，在喷管出口产生侧向推力，实现矢量控制。

图7为He注入喷管流场后在对称面上质量分数等值线图。如图7所示，He注入喷管流场后渗透到一定高度，压强逐渐降低，y方向速度逐渐减小，受喷管主流的推动作用，然后转弯平行于喷管壁面流动。

图6 流场压强分布Fig.6 Pressure distribution of flow field

图7 He质量分数等值线Fig.7 Mass fraction contour lines of helium

图8 在对称面x=0.066 m处He和N2的质量分数沿y方向的分布Fig.8 Distribution of helium and nitrogen mass fraction along y direction at x=0.066 m on symmetry plane

图8为对称面上注气口下游x=0.066 m处，He和喷管主流中的N2沿y方向的质量分数分布曲线。He在靠近壁面处质量分数很低，随着高度的增加质量分数迅速变大，由于受回流涡结构的影响，在y为-0.031～-0.028之间，质量分数的增大趋势变缓，之后随高度增加质量分数又迅速增大，在y为-0.026～-0.023之间达到相对稳定的接近于1的值，然后随着高度的增大又逐渐减小，直至为零。而主流中N2的质量分数变化与He正好相反。

图9为NPR=10，SPR=1.0，θ=90°和105°时，He、N2和CO2气体在对称面x=0.066 m处沿y方向质量分数分布曲线。由图9可知，3种气体在此处的质量分数分布情况近似相同，在相同高度上，分子质量大的气体质量分数略高。即分子质量越大，渗透高度越大。这说明，分子质量小的气体由于渗透高度较低对主流的阻碍作用就会相应降低，从而减少推力损失。

图9 在对称面x=0.066 m处不同次流气体的质量分数沿y方向的分布Fig.9 Distribution of mass fraction along y direction for different secondary flow gas at x=0.066 m on symmetry plane

表2为NPR=10，SPR=1.0时，不同注气角度下矢量性能参数对比。由表2可知，在不同注气角度下，矢量偏转角均随次流气体分子质量的增大逐渐减小。推力系数则近似相等，随分子质量的增大略微减小。分子质量越小，次流与主流的质量流量比ω越小。说明分子质量小的气体在单位质量流量下产生的矢量偏转角更大，而推力损失更小。

表3为SPR=1.0，注气角度为105°时，不同NPR下矢量性能参数对比。由表3可知，在不同NPR下，矢量偏转角均随次流气体分子质量的增大逐渐减小。另外，在NPR发生变化时，分子质量小的气体其矢量偏转角变化的程度更大。推力系数则近似相等，随分子质量的增大略微减小。

表4为NPR=10，注气角度为105°时，不同SPR下矢量性能参数对比。由表4可知，在不同SPR下，矢量偏转角均随次流气体分子质量的增大逐渐减小。推力系数则近似相等。