轴对称双喉道气动矢量喷管内特性数值模拟

轴对称双喉道气动矢量喷管内特性数值模拟卿太木;王恒;廖华琳【摘要】Numerical studies on the effect of geometrical parameters on internal performance of axisymmet-ric dual-throat fluidic thrust-vectoring nozzles were performed, and the optimized geometrical parameters of nozzle were obtained. The results indicated that vector performance was influenced by secondary inject angle, cavity divergence angle and cavity length, and thrust performance was greatly influenced by second-ary inject angle and cavity divergence angle. The optimized results indicated that vector angle could reach 15.5°, vector efficiency 5.70°/1%, and thrust ratio 0.92.%采用数值模拟方法，就轴对称双喉道气动矢量喷管主要几何参数对喷管内特性(气动矢量角、气动矢量效率、推力系数、流量系数)的影响进行了分析，并得到了一组内特性较优的喷管几何参数。

结果表明，次流注入角、空腔收敛角、空腔长度等几何因素对喷管矢量特性影响较大，次流注入角、空腔扩张角等对喷管推力特性有较大影响。

优选后的喷管矢量角可达15.5°，矢量效率达5.70°/1%，喷管推力系数为0.92。

【期刊名称】《燃气涡轮试验与研究》【年(卷),期】2014(000)002【总页数】7页(P14-20)【关键词】航空发动机;气动矢量喷管;轴对称;双喉道;几何参数;内特性;第四代战斗机【作者】卿太木;王恒;廖华琳【作者单位】中国燃气涡轮研究院，四川成都610500;中国燃气涡轮研究院，四川成都610500;中国燃气涡轮研究院，四川成都610500【正文语种】中文【中图分类】V235.13+3;V231.31 引言推力矢量技术被认为是第四代战斗机和先进无人作战飞机的必备技术，有着巨大的军事和经济意义。

气动矢量喷管作为推力矢量技术的代表，已成为替代现役和在研机械调节式矢量喷管的革命性技术方案，引起了国内外广泛重视。

双喉道气动矢量喷管，是美国NASA兰利研究中心于2003年提出的一种新型气动矢量喷管概念[1]。

兰利研究中心对二元双喉道气动矢量喷管的几何结构参数、喷管主/次流进口参数变化对喷管内特性的影响，进行了数值模拟和模型试验研究[2]，找出了二元双喉道气动矢量喷管内特性变化规律。

国内对气动矢量喷管也进行了一定研究，文献[3]对二元双喉道气动矢量喷管的几何结构参数变化进行了数值仿真分析，文献[4]对轴对称双喉道气动矢量喷管的进口气动参数进行了数值仿真分析。

本文利用数值模拟方法，对轴对称双喉道气动矢量喷管几何参数变化对喷管内特性的影响进行了分析，并优选了一组喷管几何结构参数。

2 计算模型及内特性参数定义2.1 计算方法轴对称双喉道气动矢量喷管计算几何模型见图1。

采用经文献[4]验算过的FLUENT三维解算器基于有限体积法进行耦合求解，计算模型选用基于湍流能量方程和扩散速率方程的标准两方程k-ω模型，对流项选用二阶迎风差分格式，时间项采用全隐式。

图中，D7为喷管进口直径，D8为喷管喉道直径，D9为喷管出口直径，L为喷管空腔长度，α为次流注入角度，θ1为空腔扩张角，θ2为空腔收敛角，β为次流周向角，d为次流小孔直径，N为次流孔数量。

图1 喷管模型几何参数Fig.1 The geometrical parameters of nozzle2.2 内特性参数定义(1)流量系数式中：mt为喷管主流实际流量，mc为喷管次流实际流量，mi为喷管理论流量。

(2)推力系数式中：Fx为喷管轴向实际推力，Fy为喷管法向实际推力，Fi为喷管理论推力。

(3)气动矢量角(4)气动矢量效率η=δ/[100mc/(mt+mc)]3 几何参数对喷管内特性的影响3.1 次流注入角对次流注入角α分别为30°、60°、90°、120°时的双喉道喷管流场进行模拟，其余几何结构参数中次流周角固定为60°、空腔扩张角固定为10°、空腔收敛角固定为30°、空腔长度固定为132.68 mm。

计算中，喷管主流落压比取1.89和3.00，前者为双喉道喷管设计落压比，后者根据发动机喷管常用落压比值确定。

根据主、次流流量比值范围，将次流落压比定为3.10和5.30。

矢量特性计算结果见图2，推力特性和流量特性见图3。

图2 气动矢量角、气动矢量效率随次流注入角的变化Fig.2 Vector angle and thrust vector efficiency vs.secondary inject angle图3 推力系数、流量系数随次流注入角的变化Fig.3 Thrust ratio and discharge coefficient vs.secondary inject angle可见，随着次流注入角度的增加，两种主流落压比(NPRt)下气动矢量角均呈单调下降趋势。

落压比为3.00时降幅约5.1°，落压比为1.89时降幅约5.6°；落压比为1.89的矢量角总体比落压比为3.00的高1.0°～1.7°。

落压比为 3.00 时，气动矢量角最大为11.2°；落压比为1.89时，气动矢量角最大为12.9°。

这表明次流注入角对喷管气动矢量角的影响非常明显，较小的次流注入角有利于次流流量的注入，提升喷管气动矢量角。

随着次流注入角度的增加，两种主流落压比下的气动矢量效率变化趋势一致，在前三种次流入射角下，喷管矢量效率变化不大；入射角等于90°后，矢量效率急剧下降，NPRt=3.00/1.89时，矢量效率分别为2.65°/1%和3.33°/1%。

这表明入射角90°附近是矢量效率的一个分界点，当次流注入流向与主流流向一致时矢量效率较低，次流注入流向与主流流向相反时能获得较高的矢量效率。

随着次流注入角度的增加，两种主流落压比下Cf和Cm均呈单边上升的趋势，次流的注入改变了喷管流场，进而影响了喷管的Cf和Cm。

但不同的次流注入角影响程度不一样，总体上次流注入角度越大对喷管流场的影响越小，推力损失变越小，Cf和Cm越高。

3.2 次流周向角对次流周向角β分别为50°、60°、72°、85°时的双喉道喷管流场进行模拟，结果如图4～图7所示。

可见，周向角度增加会降低喷管气动矢量角，较小的次流周向注入角有利于次流流量的注入，提升了喷管的气动矢量角。

落压比为1.89时周向角为50°的矢量效率最高，达到了4.84°/1%。

随着次流周向角的增大，气动矢量效率逐渐下降且趋势变缓，较小的次流周向角有利于提升喷管矢量效率。

第一喉道前典型截面的速度云图基本一致，受次流注入的影响微弱，流场呈现较好的轴对称性。

当次流以不同周向角注入后，在接近次流注入的下游区域，可明显看出流场的不对称性。

同时，次流周向角越小，次流注入区域越集中，对流场的影响越大，使得喷管矢量角更大、矢量效率更高，流场结果较好地反应了不同周向角下喷管矢量特性的变化结果。

图4 气动矢量角、气动矢量效率随次流周向角的变化Fig.4 Vector angle and thrust vector efficiency vs.secondary circumferential angle图5 周向角50°下喷管典型截面的速度云图Fig.5 The velocity section contours of 50°ci rcumferential angle图6 不同周向角下喷管次流区域下游截面的速度云图Fig.6 The velocity section contours of different circumferential angles图7 推力系数、流量系数随次流周向角的变化Fig.7 Thrust ratio and discharge coefficient vs.secondary circumferential angle随着次流周向角度的增加，两种主流落压比下喷管推力系数均呈单边微弱下降趋势，且降幅较小，约在0.01以内。

次流周向注入角对喷管推力系数的影响微弱。

两种主流落压比下，除周向角为60°的流量系数较低外，其余流量系数基本一致。

3.3 空腔扩张角对空腔扩张角θ1分别为5°、10°、20°时的双喉道喷管流场进行模拟结果见图8～图10。

可见，10°空腔扩张角有相对较优的矢量角，空腔扩张角增大到20°时，矢量角及矢量效率剧烈下降，5°和10°的矢量效率基本相当。

落压比为1.89的矢量角和矢量效率均比落压比3.00的高。

随着空腔扩张角的增加，喷管推力系数逐渐下降；当落压比为3.00、空腔扩张角由5°升至20°时，推力系数由0.939降至0.907。

由前文可知，扩张角过大不利于喷管矢量偏转，因此双喉道气动矢量喷管的空腔扩张角不宜大于10°。

当落压比为3.00、空腔扩张角由5°升至20°时，流量系数由0.834降至0.780，过大的喷管空腔扩张角不利于喷管流通。

图8 气动矢量角、气动矢量效率随空腔扩张角的变化Fig.8 Vector angle and thrust vector efficiency vs.cavity divergence angle图9 推力系数、流量系数随空腔扩张角的变化Fig.9 Thrust ratio and discharge coefficient vs.cavity divergence angle图10 不同扩张角下喷管的三维流线图Fig.10 The 3D streamline pattern of different cavity divergence angles从不同扩张角下喷管的三维流线图中可看出，在空腔顶部形成了回流区，5°的回流区比10°的明显偏小，因而10°的矢量偏转角较5°的大。