Docking test verification of high temperature gas ejector system
-
摘要: 引射系统是化学激光器排气系统的核心组件。为提高排气系统引射效率并适应小型化、经济性等要求,提出一种采用燃气发生器产生高温燃气的引射气源方案,建立了"基于航空发动机单喷嘴燃烧室结构的燃气发生器+两级超声速引射器"组合的引射系统试验平台,开展了燃气发生器独立热试车,并实现了与引射器的对接试验。结果表明:燃气发生器点火可靠,运行稳定;在高温燃气的引射作用下,引射器启动迅速,工作平稳,各项性能指标优于设计值;引射燃气温度一定范围内的变化基本不影响引射器工作性能。相关研究结果可为引射系统方案设计及选取提供参考。Abstract: The ejector system is the key component of the chemical laser exhaust system. In order to improve the ejection efficiency of the exhaust system and meet requirements of miniaturization and economy, an ejection gas source scheme based on the gas generator to generate high temperature gas is proposed. A test stand of the ejection system composed of a gas generator based on structure of aero-engine mono-injector combustor and the two-stage supersonic ejector is established. An independent thermal test of the gas generator has been carried out and a docking test with the ejector has been successfully realized. Results show that ignition of the gas generator is reliable and the operation is stable; the ejector starts quickly and works smoothly under the action of high temperature gas injection, and performance indexes are better than design points:the ejector-gas temperature variation within a certain range does not affect the ejector performance. Research results expand the application of the gas generator, and at the same time, provide a powerful reference for the design and selection of ejector system solutions in this field.
-
Keywords:
- ejector system /
- chemical laser /
- gas generator /
- supersonic ejector /
- docking test
-
0 引言
超空泡流动现象常见于水中高速鱼雷、潜水艇发射战略战术导弹和利用高速射弹消灭水雷等军事技术中,因此一直是国内外的热门研究课题之一[1-2],本课题组对此也进行了研究[3-5]。超空泡带迎角地穿过自由面问题源于从潜水艇鱼雷发射管发射战斧式巡航导弹,为此美国海军进行了大量的实验[6]。我国在这方面的研究还不够充分,无论是实验还是数值计算。附着在潜射导弹上的空泡一般是自然空化空泡(Natural Cavitation),空泡内部压力约为水的饱和蒸气压(2.35×103Pa)。导弹出水时,空泡与压力为1个标准大气压(105Pa)的空气接触,因此空气必然进入空泡,使其成为充气空泡(Ventilated Cavitation)。然而必须明确的是,与用人工方法给空泡充气的过程不同(诸如高速鱼雷),在出水时空泡获得的充气仍然是在自然状态下进行的。
胡影影等人[7]使用改进的Youngs算法对航行体在自由面附近出水进行了数值模拟,着重研究了韦伯数We和弗汝德数Fr在出水运动过程中对液面的影响。黄海龙等人[8]通过划分结构网格,对圆盘空化器进行了三维数值仿真,模拟中应用了k-ε湍流模型求解N-S方程,通过改变不同的空化数大小,得到相应的计算参数,然后将计算参数的结果与前人的经验公式相对比。王一伟等人[9]通过将结构力学和计算流体力学2种方法相结合,建立了弹体在带迎角出水时的数值模拟方法。Chu等人[10]模拟了圆柱航行体带空泡的出水过程,根据测量得出的压力数据能够得出空泡溃灭的大致时间和位置,发现空泡溃灭时会产生水拍打液面的现象,通过这种方式产生了高压,并且通过实验拍摄的图片与数值模拟产生的尾部空泡相对比,发现数值模拟产生的尾空泡更小一点;据此推测,这可能是由于在实际压缩气体喷射航行体的过程中,有部分气体会沿着航行体喷出,从而使实验中的尾空泡变得更大一些。张军等人[11]和赵蛟龙等人[12]分别对低速全沾湿物体和局部超空泡倾斜出水流场进行了PIV和高速摄影研究。
本课题组进行了细长钉体(横截面尺寸约为1mm)在小迎角α=5.5°~10°范围的出水实验和数值模拟[13]。本文将报道直径为6mm的钝体在较大迎角α=57°下的出水实验和数值模拟。这里使用“钝体”一词,是为了与以前的尖头钉体有所区别,但其长径比为8~12,实际也属于细长体。
1 实验装置与方法
图 1示出了实验装置。一台倾斜放置的一级轻气炮装置(见图 1中的部件4~12)向水箱自下而上地发射试验模型,模型冲击进入水箱1产生超空泡,然后穿过自由面2。圆柱体模型直径为6mm,长度分别为48、60和72mm(质量为3.69~5.50g),材料为铝镁合金。使用Photron Fastcam SA5高速摄影机记录超空泡穿过自由面的过程,然后将照片导入AutoCAD软件;测量位移等数据,每点测量3次,取平均值,最后计算出速度、阻力系数等物理量。根据冲击工程惯例,定义超空泡运动方向与水平液面的垂直方向之间的夹角为出水迎角α。有关实验技术的更多细节,可参见文献[14-15]。
2 实验结果与理论分析
图 2示出了长径比为10的平头圆柱体带迎角地穿过自由面的高速摄影照片。超空泡在图 2(1)中出现在视场中的左下角,其速度为V0=23.02 m/s;流场的标尺也在图 2(1)中给出。超空泡向右上方运动,逐渐靠近自由面(见图 2(2)~(7))。在图 2(8)所示的时刻,超空泡与自由面接触;然后在图 2(9)~(12)的时间段里,圆柱体伸出水面。从图 2(13)所示的时刻开始,圆柱体彻底脱离水面,进入空气中自由飞行。超空泡与自由面倾斜地相互作用后(见图 2(9)~(10)),出现了2个特殊的现象:一是在水面上沿着运动方向、即沿着超空泡的轴线,出现了不对称的水花。二是负压的超空泡从大气中倒吸空气,形成了充气空泡。如图 2(11)~(15)所示,水花先在轴线的上部形成并向右上方运动,水花头部的速度与圆柱体在空中的飞行速度接近(水花头部速度大于物体速度也是可能的)。从图 2(16)开始,轴线下部的水花才明显地形成,然后逐渐发展,与上部水花合并形成喇叭口的水花(见图 2(23)~(24))。水花以如此的方式形成,就使得空泡在水面上难以产生表面闭合(Surface Closure),空气被吸入负压的空泡中,因此空泡能在自由面附近长期存留而不会破灭。从图 2(11)~(24),空泡已存活了20ms的时间,而且空泡的最大直径从32mm增加到了46mm,尽管空泡后面的空化尾迹已经十分明显(见图 2(21)~(24))。如果对比本文第4节的数值计算结果,就可以理解上述的自然空化空泡向充气空泡转化的过程。
图 3(a)和(b)分别给出了测得的长径比为8、10和12的圆柱体的速度V随时间t的变化,以及阻力系数CD随空化数σ的变化。其中空化数σ是物体在每个时刻下计算出的,因此图 3(b)的横轴不需用时间t表示。从图 3(a)中可以看出,3个圆柱体在出水前后,即在t=10.00、12.86和14.27ms时刻,其速度都有一个跃增。这是因为物体从水中突然进入空气中甩负荷的缘故[5]。如果在穿过水面前后,在水中和在空气中的物体速度分别为Vw和Va,那么实验得到的Va/Vw=1.06~1.11。图 3(b)中给出的σ的最大值即是出水前的空化数,3个圆柱体的该空化数分别为σ=1.134、0.969和0.496。
圆柱体在水中和空气中的运动方程分别为:
(1) (2) 式中:m、aw、A、CD分别为物体的质量、加速度、截面积和阻力系数,ρ、g和α分别是流体密度、重力加速度和出水迎角,下标w和a分别表示在水中和空气中。将两式相除,整理后得到
(3) 根据图 3(a)可知,物体在出水前呈单调地减速,计算出的3个工况的加速度非常接近,其绝对值为1.04×103m/s2。运动物体在空气中的雷诺数Re=(5.93~7.01)×103,所以物体在空气中的阻力系数CDa约为1.00[16-17]。参考图 3(b),得知在出水前,物体的阻力系数CDw约为0.90~1.10,因此不妨认为CDw/CDa →1。然后从式(3)中算出Va/Vw=2.03。这个理论值大于实验值,这是因为物体刚出水时,虽然上部已露在空气中,但下部仍与水沾湿,此时阻力较大,会出现CDw/CDa < 1,这样计算出的Va/Vw的理论值会接近于实验值。
根据图 3(a)还可知,物体在出水后基本上保持恒速,因此式(2)中本应该出现的加速度项maa→0。上述有关出水前后物体速度变化的机理,如果从动量守恒定律来分析,就可以理解为:既然空泡引起的附加质量m*在出水后自动消失,而物体在穿过自由面时又没有受到其他外力作用,那么显然(m+m*)Vw=mVa。所以,
(4) 3 网格划分及计算方法设置
本文需要关注水中物体产生的空泡,所以在物体周围的网格划分较为密集,对物体周围加密时采用了尺寸函数,定义开始尺寸为1mm,增长率为1.2,最大尺寸为2mm。用Gambit软件对计算区域进行划分网格,且对物体周围进行局部加密,计算区域为1500mm×400mm,本文采用三角形非结构网格,网格总数为508 499。边界条件设置为:计算区域入口设置为wall边界,出口设置为压力出口,航行体边界设置为wall边界,如图 4所示。
将模型导入到Ansys14.0中,进行计算参数设置,具体步骤为:计算采用VOF多相流模型(水、空气和水蒸气),使用压力基求解器,计算考虑了重力的影响,重力加速度设置为-9.81m/s2。由于超空泡在运动过程中对时间的依赖性很强,所以选择非定常流动,空化模型采用的是Schnerr and Sauer模型,水的饱和蒸气压为3540Pa。求解方法采用SIMPLE算法。网格运动过程中的更新方法为弹簧近似光滑模型和网格重构模型实现,将编好的profile文件导入Ansys14.0中实现航行体的变速运动。
4 计算结果
图 5给出了水下超空泡物体倾斜出水过程相图的数值模拟结果,包括超空泡从生长到消失的全过程。已计算了各种不同的工况[15],这里只给出一组典型的数据。空泡开始出生时的位置主要发生在物体的尾部和肩部(见图 5(1)~(2)),随着时间的增大,空泡慢慢生长,直到将物体完全包裹、形成超空泡(见图 5(3)~(6))。然后物体带超空泡继续前进,直到穿出水面,并且在物体尾部带出水花(见图 5(15)~(16))。数值计算揭示了当水下超空泡物体接近自由面时(见图 5(5)~(8)),水面必然抬升隆起的、由水的不可压缩性所导致的后果。数值计算的另一个有意思的结果是,从图 5(14)~(16)的不到2ms的时间里,超空泡一接触水面就消失了,大大小于在图 2的实验中观察到的、大于20ms的超空泡的寿命。这是因为计算软件中的空化模型是完全自然空化模型,没有考虑到空泡从大气中吸取空气的机理。因此在今后建立充气空泡模型并改进软件的算法是十分必要的。另外,在第2节中讨论式(3)的计算结果时,认为物体刚出水时上半部分露在空气中,而下半部分与水沾湿,这个推断也被数值计算证实了。今后要注意的是:根据空化数的定义,当空泡浸没在水面下时,空泡壁面上的压力是水的饱和蒸气压;当空泡露出水面进行通气后,空泡壁面上的压力是大气压,有限体积算法不一定能处理好这种问题。Young和Savander[18]用边界元方法计算了部分浸没的空泡以及表面划水螺旋桨叶片的流场,值得参考。
图 6给出了3个不同时刻水下物体周围的压力图(压力的单位是Pa),其相图分别是图 5(3)、(4)和(8)。从图 6(a)中可以看出,在t=2.57ms时,超空泡离自由面还较远,此时从相对应的压力图中可以看出,物体的周围压力对称分布。随着超空泡靠近自由面(见图 6(b)和(c)),物体头部下方的压力比上方大。在研究浅水区的超空泡运动特性时,也发现了空泡下部的压力大于上部的压力[19]。这说明当超空泡靠近自由面时,浮力发挥作用,这对空泡后期的发展,包括形状和溃灭特性,产生一定的影响。
5 结论与讨论
(1) 超空泡倾斜地与自由面相互作用时,会形成非轴对称的、开口喇叭形的水花,因此原来处于负压的空泡更容易地吸取大气中的空气,空泡从自然空化向充气空泡转变,空泡寿命增加(大于20ms)。
(2) 物体出水前后,速度有一跃增,速度增加率达到6%~11%。通过分析运动方程,从理论上证明了当物体从水中进入空气中时,其速度的确是增加的。运用动量守恒定律,也得出相同结论。对于长径比为10的物体,如图 3(a)中所示的、在t=21.43ms时刻出现的二次增速,是由测量误差造成的;更多的测量结果,可见文献[15]。
(3) 当水中超空泡物体接近自由面时,自由面将隆起抬高,这是后期表面波浪形成的初始原因;自由面对超空泡的影响之一,是空泡下部的压力大于上部的压力。严格地说,用本文的方法产生的超空泡里,会夹带空气的进入。但是,实验中发射管出口离水箱侧面隔膜较近,这样对于空气进入有抑制作用。而数值计算的空泡形状与实验结果比较接近,说明即使有空气进入,其量也是有限的。
-
表 1 燃气发生器设计参数
Table 1 Design parameters of gas generator
设计参数 数值 余气系数 3.0 总流量 26 kg/s 燃烧室压力(无量纲) 3.0~4.0 燃气总温 1100 K 出口温度偏差 ±50 K 表 2 引射器设计参数
Table 2 Design parameters of ejector
设计参数 数值 被引射气流总压 4.00 kPa 被引射流量 0.32 kg/s 引射马赫数 一级 5.0 二级 3.9 引射流量 一级 0.62 kg/s 二级 2.63 kg/s 引射压力(无量纲) 一级 2.41 二级 2.83 引射温度 一级 800 K 二级 900 K -
[1] 陈振华, 任泽斌, 赵波, 等.化学氧碘激光器尾气主动冷却技术试验[J].强激光与粒子束, 2014, 26(8):081013. DOI: 10.11884/HPLPB201426.081013 CHEN Z H, REN Z B, ZHAO B, et al. Experimental research on active cooling technique of exhaust gas from COIL[J]. High Power Laser and Particle Beams, 2014, 26(8):081013. DOI: 10.11884/HPLPB201426.081013
[2] 李金雪, 王杰, 颜飞雪, 等. DF化学激光器扩压器流场仿真及优化[J].强激光与粒子束, 2018, 30(10):101002. DOI: 10.11884/HPLPB201830.180113 LI J X, WANG J, YAN F X, et al. Simulation and optimization of DF chemical lasers' diffuser[J]. High Power Laser and Particle Beams, 2018, 30(10):101002. DOI: 10.11884/HPLPB201830.180113
[3] 符澄, 彭强, 刘卫红, 等.光腔与扩压器化学反应流场优化数值模拟[J].强激光与粒子束, 2015, 27(11):111009. DOI: 10.11884/HPLPB201527.111009 FU C, PENG Q, LIU W H, et al. Numerical simulation of chemical reaction flow optimization in cavity and diffuser[J]. High Power Laser and Particle Beams, 2015, 27(11):111009. DOI: 10.11884/HPLPB201527.111009
[4] CONNAUGHTON J C. Application of a hydrazine gas generator to vacuum ejector pumping of a chemical laser[R]. AIAA-77-0892, 1977.
[5] 邹建军, 周进, 徐万武, 等.大压缩比两级燃气引射系统参数匹配试验研究[J].强激光与粒子束, 2008, 20(10):1589-1592. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs200810002 ZOU J J, ZHOU J, XU W W, et al. Experimental investigation on parameter matching of high compression-ratio two-stage gas ejector system[J]. High Power Laser and Particle Beams, 2008, 20(10):1589-1592. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs200810002
[6] 廖达雄.气体引射器原理及设计[M].北京:国防工业出版社, 2018:72-73. [7] 徐万武.高性能、大压缩比化学激光器压力恢复系统研究[D].长沙: 国防科学技术大学, 2003. XU W W. Study of high performance, high compression ratio pressure recovery system for chemical laser[D]. Changsha: National University of Defense Technology, 2003.
[8] 徐万武, 邹建军, 王振国, 等超声速环型引射器启动特性试验研究[J].火箭推进, 2005, 31(6):7-11. DOI: 10.3969/j.issn.1672-9374.2005.06.002 XU W W, ZOU J J, WANG Z G, et al. Experimental investigation of the start performances of the supersonic annular ejector[J]. Journal of Rocket Propulsion, 2005, 31(6):7-11. DOI: 10.3969/j.issn.1672-9374.2005.06.002
[9] 徐万武, 谭建国, 王振国.高空模拟试车台超声速引射器数值研究[J].固体火箭技术, 2003, 26(2):71-74. DOI: 10.3969/j.issn.1006-2793.2003.02.020 XU W W, TAN J G, WANG Z G. Numerical study on supersonic ejector of simulated altitude test stand[J]. Journal of Solid Rocket Technology, 2003, 26(2):71-74. DOI: 10.3969/j.issn.1006-2793.2003.02.020
[10] 徐万武, 谭建国, 王振国.二次流对超声速环型空气引射器真空度的影响[J].国防科技大学学报, 2003, 25(3):6-9. DOI: 10.3969/j.issn.1001-2486.2003.03.002 XU W W, TAN J G, WANG Z G. The effects of second-flow on the vacuum-degree of annular supersonic air ejector[J]. Journal of National University of Defense Technology, 2003, 25(3):6-9. DOI: 10.3969/j.issn.1001-2486.2003.03.002
[11] 陈志强, 廖达雄, 刘宗政, 等.过氧化氢加酒精补燃气体发生器实验研究[J].强激光与粒子束, 2007, 19(9):1409-1412. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs200709001 CHEN Z Q, LIAO D X, LIU Z Z, et al. Experimental studies of gas generator based on hydrogen peroxide and hypergolic ethanol[J]. High Power Laser and Particle Beams, 2007, 19(9):1409-1412. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs200709001
[12] 刘盛田, 胡兴伟, 柳琪, 等.一氧化二氮/乙醇燃气发生器试验研究[J].舰船科学与技术, 2008, 30(z6):223-226. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jckxjs2008z2012 LIU S T, HU X W, LIU Q, et al. The experimental studies of gas generator based on nitrous oxide and ethanol[J]. Ship Science and Technology, 2008, 30(z6):223-226. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jckxjs2008z2012
[13] 刘盛田, 胡兴伟, 巴俊州, 等. N2O/C2H5OH单喷嘴燃气发生器研究[J].舰船防化, 2007(2):30-34. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK200702210228 LIU S T, HU X W, BA J Z, et al. The Study of N2O/C2H5OH mono-injector gas generator[J]. Chemical Defence on Ships, 2007(2):30-34. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK200702210228
[14] 冯军红, 沈赤兵, 赵芳.液氧与气氧对空气加热器燃烧流场的影响分析[J].国防科技大学学报, 2012, 34(4):43-48. DOI: 10.3969/j.issn.1001-2486.2012.04.009 FENG J H, SHEN C B, ZHAO F. Analysis of liquid and gaseous oxygen influence on the combustion flow field of air heater[J]. Journal of National University of Defense Technology, 2012, 34(4):43-48. DOI: 10.3969/j.issn.1001-2486.2012.04.009
[15] 罗智锋, 刘重阳, 冯大强, 等.一种以乙醇为燃料的单管燃烧室燃气发生器: 201420794555.3. 2015-08-05. LUO Z F, LIU C Y, FENG D Q, et al. A single tube combustor gas generator fueled by ethanol: 201420794555.3. 2015-08-05.
[16] 廖达雄, 任泽斌, 余永生, 等.等压混合引射器设计与实验研究[J].强激光与粒子束, 2006, 18(5):728-732. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs200605007 LIAO D X, REN Z B, YU Y S, et al. Design and experiment of constant pressure mixing ejector[J]. High Power Laser and Particle Beams, 2006, 18(5):728-732. http://d.old.wanfangdata.com.cn/Periodical/qjgylzs200605007
-
期刊类型引用(3)
1. 田建辉,胡晨明. 高速弹体出水过程数值模拟. 兵器装备工程学报. 2023(12): 73-78+176 . 百度学术
2. 李易松,袁航,张建宇,冯会方. 航行体跌落试验研究. 舰船科学技术. 2023(24): 218-221 . 百度学术
3. 侯健,施红辉,孙亚亚,高聚瑞. 高速射弹出水过程中水弹道问题研究. 弹道学报. 2017(01): 51-56 . 百度学术
其他类型引用(2)