基金项目:国家自然科学基金(11602186)
作者简介:李舒欣(1992—),女,硕士,研究领域为液体火箭发动机热过程。
Science and Technology on Liquid Rocket Engines Laboratory,Xi'an Aerospace Propulsion Institute, Xi'an 710100,China
detonation wave; flame acceleration; methane; liquid rocket engine
为了解液氧甲烷火箭发动机爆震波点火器工作特性,采用实验方法研究了低供给压力下甲烷/氧气爆震波的形成过程,分析了点火时序、混合比等参数对爆震波参数的影响、缓燃向爆震转变增强措施对爆震波形成过程的作用规律以及结构尺寸对弯管中爆震波参数的影响。研究结果表明:1)不同当量比下,甲烷氧气爆震管点火器的工作状态可分为爆震、快速火焰、慢速火焰、全局淬熄和未点火等状态。阻塞比在试验涉及的0.12至0.42范围内时,随着阻塞比的增大,能产生爆震的当量比范围逐渐扩大; 当量比在1.5附近,阻塞比的大小对是否能形成爆震影响不大。2)弯管及管中弹簧障碍物对于压力高、速度快,接近于爆震波的火焰表现为抑制作用; 对于压力低、传播速度慢的缓燃火焰,则能起到增强作用。
In order to understand the working characteristics of detonation wave igniter for LOX/CH4 rocket engine,the detonation wave formation of methane/oxygen under low supply pressure was investigated by experiment.The effects of ignition sequence,mixing ratio and other parameters on detonation wave parameters,the effect of enhancement measures of deflagration to detonation transition on detonation wave formation,and the effect of structure size on detonation wave parameters in bent pipe were analyzed.The results show that: 1)Under different equivalent ratios,the working status of the detonation wave igniter can be divided into detonation,fast flame,slow flame,global quenching and unignited states.When the blockage ratio is in the range of 0.12~0.42,the range of equivalent ratio that can produce detonation increases with the increased blockage ratio.When the equivalent ratio is about 1.5,the blockage ratio has little effect on the formation of detonation wave.2)U shape tubes and the spring in it have suppression effect on the flame for the detonation wave with high pressure and fast speed.In addition,for the flame with low pressure and slow propagation speed,it can play an enhanced role.
0 引言
爆震波点火技术是利用低压可燃混气产生的爆震波实现发动机点火启动的一种新型点火技术。该技术应用于液氧甲烷火箭发动机,具有可实现多次点火、降低点火系统供给压力、简化点火系统与保证点火同步性[1]等优势。
爆震波点火的过程如下:甲烷与氧气进入预燃点火室,并在预燃点火室以及爆震波导管中充分混合。随后,电火花塞打火点燃预混点火室中气体,产生缓燃火焰,缓燃火焰在爆震导管中加速,实现缓燃向爆震转变(deflagration to detonation transition,DDT)进而形成爆震波。爆震波及高温燃气传入推力室并将其引燃。
目前,对于液氧/甲烷火箭发动机的点火研究集中在火炬点火方式上[2-8],然而火炬点火对系统供应压力的要求相对较高,使得点火系统相对复杂,这在实现多机并联同步点火时更加明显。为此,有研究者尝试利用爆震波自增压的特性,实现液体火箭发动机低供给压力下的点火。北京航空航天大学郭红杰等人2004—2006年搭建了氢气/氧气爆震波点火试验系统,氢气/氧气供给压力为0.1~0.5 MPa,混合比为1.6~9.4,在常温条件下得到了爆震波在管路中的传播特性以及爆震波点火的单管重复性能试验数据。试验结果显示爆震波点火可在较低的供气压力下获得高温高压的燃烧产物[9-11]。
纵观国内外现有研究成果,面向飞行器发动机的爆震波点火器研究还非常有限。美国NASA在1996—2001年期间组织的X—33飞行器项目采用了塞式发动机[12-13],该塞式发动机中采用了液氧/液甲烷点火剂的爆震波点火方案,全尺寸样机地面低温点火试验验证了方案的可行性。
现阶段针对液氧/甲烷火箭发动机爆震波点火技术的研究较少,对甲烷/氧气缓燃转爆震过程的认识也不够清楚,但国内外在氢气/空气爆震、甲烷/空气爆震、碳氢燃料/空气爆震等方面开展了较为广泛的研究。在爆震理论研究方面,国内外学者们对爆震波的形成给出了不同定性解释:热点理论、梯度理论和相干能量释放的激波放大理论[14-15]。为了适应空间布局的局限性,爆震波在弯管内的传播特性也受到了研究人员的关注。研究者们使用压力测量、烟迹和纹影等方法观察到爆震波在弯管壁面凹侧强度增大,在凸侧由于发散作用产生了膨胀波使得强度降低甚至发生局部淬熄,由于波阵面上横波的作用,爆震波局部淬熄后又得以重新起爆。弯管的曲率半径及管中初始压力是影响爆震波通过弯管时锋面压力、传播速度的重要因素,曲率半径越大,初始压力越高则爆震波越容易通过弯管而不熄灭[16-19]。
从国内外研究现状可知,国内外对液氧/甲烷火箭发动机爆震波点火过程的研究较少,制约了爆震波点火技术在液氧/甲烷火箭发动机中的应用。虽然国内外对爆震的研究成果较多,但其研究更多采用的是氢气/空气、氢气/氧气、甲烷/空气等介质,采用甲烷/氧气的相对较少,对甲烷/氧气燃烧转爆震过程的认识较为缺乏。因此,有必要开展甲烷/氧气爆震波点火技术的研究以推进其在液氧/甲烷火箭发动机中的应用。考虑到工程应用中对点火器低压供给、小包络的要求,本文以试验的形式,对低供给压力下的爆震现象及DDT增强装置和弯管这两种可能的缩减爆震波点火器包络尺寸的结构进行了研究。
形成足够强度的甲烷/氧气爆震波是实现爆震波点火技术的首要条件。液氧/甲烷火箭发动机实际应用中,推进剂贮箱提供的压力相对较低,低压条件下甲烷/氧气爆震波形成过程及其影响因素作用规律尚不十分清楚,本节结合液氧/甲烷火箭发动机中的实际情况,对爆震波形成过程中填充压力、混合比等因素的作用进行考察。
试验所用爆震管装置如图1所示,分为头部、DDT段和加长段3个部分。其中甲烷及氧气由头部接入,通过同轴剪切喷嘴实现掺混,并由火花塞点燃,火焰在DDT段中完成缓燃转爆震的过程。DDT段的圆管内部加装弹簧,起促进DDT过程的作用[20]。DDT段后接加长段,用于测量爆震波通过时的压力数值。整个管道内径为10 mm,尾端未设置喉部。
试验中,利用气源储箱对爆震管进行气体供应,每一路所用孔板固定不变,通过调节供气压力来调节流量。CH4路供气压力在0.58~0.86 MPa范围内,O2路为0.27~2.74 MPa。先后对4种供气时序进行了测试,每种时序详情如表1所示,每种时序均可实现爆震管的满填充。
爆震波波速是试验中获取的重要数据,其计算方法为两脉动压力传感器座间的距离除以爆震波传播通过两个传感器测点所需的时间,即
v=Δx/Δt^-(1)
式中:v为爆震波传播速度; Δx为相邻两传感器距离; Δt^-为多次测量得到的爆震波锋面时间差的平均值。爆震波传播速度v的误差分析方法参见文献[21]。计算可知,当理论爆震波波速为2 389.2 m/s时,试验用装置及测量系统测量结果为(2 389.2±93.64)m/s,可认为测量误差对是否产生爆震的判定不产生影响,在可以接受的范围内。
试验测试了不同点火时序(表1)下当量比在0.47至7.05范围内的一系列工况。记安装在加长段上游的脉动压力传感器为1#,下游为2#,则测得的典型脉动压力曲线如图2所示,该图显示了氧气孔板前压力0.86 MPa、甲烷孔板前压力0.59 MPa、当量比为1.54时加长段中测得的压力曲线。当爆震波经过压力测点时,压力曲线出现一个陡峭的峰值,1#测点测得的峰值压力为2.202 MPa,2#测点测得峰值压力为2.428 MPa。随后曲线缓缓下降,此阶段压力变化符合Taylor波的变化规律。当爆震波从下游管口传出后,膨胀波开始向上游反传,测点压力值不断降低,直至环境压力。记录1#和2#脉动压力传感器捕捉到爆震波峰值压力的时刻,即可由式(1)计算得到爆震波在加长段中的传播速度。计算得到爆震波传播速度为2 625.0 m/s,与CEA计算得到的CJ速度2 591.1 m/s相比略高,这一差别可认为由测量误差造成。
通过上述测量方法,试验共获得38个数据点,各工况下的爆震波速度如图3所示,图中黑实线和虚线分别表示由CEA计算得到的CJ爆震速度和燃气中的声速。从中可以看出按照当量比的大小不同,其对爆震波形成及爆震波参数的影响可分为4个区域。
1)当量比较小(<0.75)时,管中点火成功率不高。
2)当量比在0.75~1.12范围时,火焰传播速度得到一定的提升,但未加速至产生爆震波的情况。火焰被加速至约70%~100%燃气声速。
3)当量比在1.4~3.4范围内时,试验测得的爆震波速度与CJ速度相当。然而并非该当量比范围内的所有工况均能形成爆震,分析认为能否形成爆震与多种因素相关,如管中的湍流强度等,仅有适合的当量比不足以保证爆震的形成。
4)当量比上升至4.2后,出现起爆困难。
试验中涉及到的4种时序皆为氧气和甲烷气同时开启同时关闭,时序2中在氧气、甲烷供气停止后间隔1s再使火花塞放电,在时序1、3、4中均为供气停止后火花塞立即放电,区别在于填充时间长度不同。
图3中展示了不同时序下的点火情况,从图中可见,包含了1 s间隔的时序2无法点燃管中气体,原因在于在停止供气的1 s内,可燃气体向管道下游流动,爆震管头部没有气体补充,使得头部可燃气压力过低,混气浓度低,难以点燃。采用时序4时未形成爆震可能由于工况安排较少及试验中的误差造成,为进一步探明原因需增加试验次数,并综合考虑分析管中的湍流状态。
小型化是爆震波点火器工程化过程中所必须解决的一环。通常会在爆震管中加入各种形式的障碍物,如Shchelkin螺旋、孔板、斜劈等[22-26]来达到加速爆震波形成,缩短DDT距离,实现短距起爆的目的。本节使用弹簧在爆震管内设置Shchelkin螺旋,通过试验,研究了不同几何参数的Shchelkin螺旋对甲烷/氧气DDT过程的影响,分析了DDT增强措施对甲烷/氧气爆震波形成过程的影响规律。
本节研究使用的试验装置为在第1.1节中使用装置的基础上,在其DDT段中替换不同尺寸的弹簧,即变化Shchelkin螺旋的几何尺寸得到。涉及弹簧的规格如表2所示。
试验中向搭载不同阻塞比弹簧的爆震管及不搭载弹簧的光管通入不同当量比的甲烷/氧气气体并使用火花塞点火,通过加长段上的两个测点测得管中气体的压力及燃烧波传播速度。
试验中测得的压力波形可分为如下3类,如图4所示。
图4(a)所示情形加长段2个压力测点测得的峰值压力超过70% CJ速度,且计算得到的压力波传播速度与CJ速度相近,可认为此时管中产生了爆震波; 图4(b)中2个压力测点获得的峰值压力较低,约为CJ压力的20%~40%,该种现象称为快速火焰(fast flame)[27],是一种火焰传播速度被提升,但无法发展成爆震的情况; 图4(c)中,1#压力测点获得的压力数值始终较为平稳,2#压力测点在火花塞放电后的一段时间后出现一次小范围波动,认为是火花塞放电脉冲造成,分析认为该种波形对应的管中现象可以有3种:未点燃、慢速火焰(slow flame)或发生全局淬熄。慢速火焰即管中气体点燃后始终保持较弱的缓燃状态,燃烧波传播速度慢,火焰面前后压差小。全局淬熄指管中气体在点燃后,在与障碍物的作用中发生了全面熄灭。
对搭载不同阻塞比弹簧的爆震管通入不同当量比的甲烷/氧气气体后点火,分析加长段中两个压力测点获得的压力数据,依照上述方法对管中发生的燃烧现象进行判断和分类,获得图5。
为了表述方便,下文中将未点燃、慢速火焰及全局淬熄3种情况统称为“情况三”,分析图中数据点可发现以下特征:
1)阻塞比在试验涉及的0.12~0.42范围内时,无论阻塞比的大小,当量比在1.5左右时管中均能产生爆震。由CEA计算知,该工况下CJ速度的峰值出现在当量比约为1.8时,可认为当混气当量比处于爆震波释放能量最大时,DDT发生的可能性最大,此时阻塞比的大小对是否形成爆震影响较小。
2)阻塞比从0.12向0.42增大的过程中,能产生爆震的当量比范围逐渐扩大。
3)当阻塞比较小(阻塞比小于0.3)时,快速火焰和“情况三”之间的界限不明显,即在当量比偏离化学恰当比较多时,快速火焰和“情况三”均有发生的可能。
图5 不同障碍物阻塞比(BR)及混气当量比下的燃烧状态
Fig.5 Combustion modes under different block ratios(BR)and equivalence ratios
当量比决定了可燃混气的释热能力,混气的释热能力决定了火花塞放电后管中会出现何种燃烧现象,释热能力越强则发生DDT过程的可能性越大,当释热能力过小时,在障碍物的影响下,则有熄灭的可能。障碍物对火焰的作用包含了促进和抑制两个方面,选取合适的阻塞比及当量比将更有益于爆震的形成。本节中得到的结论主要着眼于阻塞比与当量比共同作用于燃烧波时其间的关系,试验中使用了不透明的不锈钢管故不便于测得DDT距离,进一步的研究拟采用光学管并获得障碍物在缩短DDT距离,即缩短爆震波点火器长度方面的影响效果。
点火器产生的热燃气包含能量越多则点火能力越强,为能产生更多的热燃气则需要在爆震管中积累更多的可燃混气,这就使得爆震管的长度不能一味缩短。采用直管形式,相对而言整体尺寸较大,会限制爆震波点火器的使用范围。使用带有弯曲结构的异型爆震管不仅能使点火器在总装布局时的适应性,还有可能增强DDT过程。本节中通过试验,研究了不同曲率下甲烷/氧气爆震波的形成和发展过程,分析爆震管构型对甲烷/氧气爆震波参数的影响规律,为发展异型爆震波点火器形成相应设计准则。
试验所用装置如图6所示。试验装置分为头部、弯管、加长段3个部分,3部分之间由法兰相连。其中头部与1.1节中介绍的相同,加长段由1.1节中的DDT段去掉弹簧充当。整个管道内径为10 mm。弯管段由一组不同弯曲半径的弯管替换件组成,在试验中可变换弯管曲率半径,弯管内侧弯曲半径R(图6所示)分别为40、70、90 mm,弯管截面外径为14 mm,即弯管轴线半径为47、77 mm、97 mm。在弯管的两端设置有两个脉动压力传感器安装座(1#、2#)。对于每种轴线半径的弯管,共制作2根,其中一根管中置入弹簧,另一根管内光滑。置入弹簧规格与1.1节中相同,为丝径2.1 mm、外径10 mm、节距8 mm。加长段用于测量弯管下游爆震波峰值压力及波速,加长段尾部为3#压力传感器安装座。2#与3#间的距离固定不变为236 mm。
考察试验中在3个测点测得的压力曲线,发现所得曲线可归为以下几种典型形式,如图7所示。
图7(a)、(b)、(c)中情况,点火后,燃烧波在管中传播时压力不断降低,然而压力下降的过程有所区别。在图7(a)对应的工况中,1#压力测点处测得的峰值压力高达2.2 MPa,传播314 mm到达2#测点时峰值压力为1.19 MPa,再传播236 mm到达3#测点时峰值压力为0.43 MPa。1#、2#测点间的平均传播速度为1 048.6 m/s,2#、3#测点间的平均传播速度为874.1 m/s。燃烧波压力在整个传播过程中压力的衰减速率较为稳定,在弯管中压力衰减至进入弯管时的一半。在该工况中97 mm弯曲半径的弯管及管中的弹簧对火焰传播速度的增加表现出抑制作用,且火焰在通过弯管段后,压力和传播速度均持续下降。
在图7(b)对应的工况中,1#和2#测点测得的峰值压力均较高,1#处测得峰值压力2.1 MPa,2#处测得的峰值压力为2.0 MPa,即锋面压力可在弯管中维持。1#、2#测点间平均火焰传播速度为1 212.2 m/s。火焰峰值压力接近CJ压力(2.47 MPa),但速度与CJ速度(2 428 m/s)相差较大,可认为此时火焰在弯管中处于一种接近爆震的准爆震状态,在该供气条件下,47 mm弯曲半径弹簧弯管对火焰有一定的加速作用,但还不足以促使DDT过程完成。火焰在传出弯管,进入直管加长段后,压力快速衰减至0.5 MPa,加长段中的平均传播速度为907.7 m/s,这一结果与2.2节中结果一致,即当φ=2.88时,内径10 mm的光滑直管对火焰的加速作用非常有限,较难出现DDT现象。
在图7(c)对应的工况中,1#测点测得的压力较高,为2.0 MPa,然而在弯管中迅速衰减至接近大气压,当火焰传播至2#时,压力仅为0.22 MPa。火焰传至3#时,压力进一步衰减至0.18 MPa。1#、2#测点间的平均传播速度为983.1 m/s,2#、3#测点间的平均传播速度为513.0 m/s,火焰传播速度持续降低。对比图9与图7,图9中燃烧波波后压力的衰减速率远快于图7,鉴于两工况中弯管为相同的97 mm弯曲半径弹簧弯管,可认为造成这一区别的原因仅在于当量比的不同,当量比过高导致管中混气释热能力不足,弯管段造成的能量耗散影响变得更为显著。
图7(d)、(e)对应的工况中爆震管中均形成了爆震波。在图7(d)对应的工况中,火焰在从1#到2#再到3#测点传播的过程中压力不断上升,传播速度越来越快。1#、2#、3#测点测得的峰值压力分别为1.02 MPa、1.78 MPa、2.22 MPa,1#与2#测点间的平均传播速度为1 429.9 m/s,2#与3#测点间的平均传播速度为2 145.46 m/s,可认为在3#处已经形成了爆震波。
在图7(e)对应的工况中,1#、2#、3#测点测得的峰值压力分别为1.95 MPa、1.44 MPa、2.45 MPa,1#与2#测点间的平均传播速度为1 258.9 m/s,2#与3#测点间的平均传播速度为1 966.7 m/s。图7(d)与图7(e)的不同之处在于从1#到2#,图7(d)中火焰压力上升,图7(e)中火焰压力下降,但在这两种情况下最终均能形成爆震波。
对比“火焰衰减”和“形成爆震”两种管中火焰发展的形式,可发现在发生火焰衰减的工况中,1#测点测得的峰值压力均较高,在最终能形成爆震的工况中,1#测点测得的峰值压力有高有低。分析认为弯管及管中弹簧障碍物对火焰的作用随火焰压力高低、传播速度的快慢而有所不同:对于压力高、速度快,接近于爆震波的火焰表现为抑制作用; 对于压力低、传播速度慢的普通缓燃火焰,则能起到增强燃烧,促使DDT过程发生的作用。
图8展示了装有弹簧的不同轴线半径的弯管中火焰传播现象与当量比的关系。总体而言,当量比在1~2范围内时,可通过装有弹簧的弯管形成爆震波,使爆震波点火器的包络尺寸在长度方向大大缩减。如将1.1节中230 mm的DDT段换为R=77 mm的弯管,则长度减小了153 mm。随着当量比的增加,DDT过程逐渐变得困难,弯管对火焰的发展显现出抑制作用。对比不同轴线弯曲半径的弯管,半径较大时(R=97 mm),在高当量比区域弯管对火焰发展的抑制作用十分强烈,这是由于大半径弯管的弯管段长度较大,火焰在弯管中传播的距离更长。
图9就相同弯曲半径下置入弹簧与否进行了对比。在两个97 mm弯管的对比中可以发现,比起置入弹簧,在高当量比范围时光滑管对火焰的抑制作用较弱。对于弯管形式的爆震管,障碍物形式的爆震增强装置易于表现出其抑制火焰传播方面的作用,应根据混气当量比,更为谨慎的选用障碍物形式的爆震增强装置。
1)按照当量比的大小不同,其对爆震波形成及爆震波参数的影响结果可分为4个区域,依次为点火困难、缓燃、爆震、起爆困难,其中缓燃又可细分为快速火焰和慢速火焰两种状态,点火困难可分为全局淬熄和未点火两种状态。
2)甲烷/氧气混气当量比在1.5附近,即混气释热能力最强时,在障碍物的作用下,DDT发生的可能性最大,此时阻塞比在0.12~0.42的范围内时,阻塞比的大小对是否能形成爆震影响不大。即设计爆震波点火器时,若取混合比在1.5附近时,可根据对管路流阻的设计要求在较宽的范围内选择阻塞比。当阻塞比较小(阻塞比小于0.3)时,快速火焰和“未点燃、慢速火焰、全局淬熄”之间的界限不明显。
3)对于压力高、速度快,近似于爆震波的火焰,弯管及管中障碍物表现为抑制作用; 对于压力低、传播速度慢的普通缓燃火焰,则能增强燃烧,促使DDT过程发生; 当量比在1~2范围内时,可通过装有弹簧的弯管形成爆震波,随着当量比的增加,DDT过程逐渐变得困难,弯管对火焰的发展显现出抑制作用。在设计爆震波点火器时应根据发动机总装布局合理使用弯管结构。