第53卷第1期2021年2月
南京航空航天大学学报
Journal of Nanjing University of Aeronautics & Astronautics
Vol. 53 No. 1Feb.2021
DOI : 10. 16356/j. 1005-2615. 2021. 01. 018
三头部燃烧室点火过程中火焰传播特性研究
冯剑寒*
1,郎旭东2,金 武打任勇智打李建中1收稿日期:2019-12-13;修订日期:2020-01-26第一届世界杯冠军
通信作者:李建中,男,教授,博士生导师,E-mail :***************** 。
引用格式:冯剑寒,郎旭东,金武,等.三头部燃烧室点火过程中火焰传播特性研究[J ].南京航空航天大学学报,2021,53(1) : 150-157. FENG Jianhan ,LANG Xudong ,JIN Wu ,et al. Study on flame propagation characteristics of a triple-dome combustor during ignition process] J]. Journal of Nanj
ing University of Aeronautics & Astronautics ,2021,53(1): 150-157.
(1.南京航空航天大学能源与动力学院,南京210016; 2.中国航空发动机集团有限公司湖南动力机械研究所,
株洲 412002)
摘要:为了探讨三头部燃烧室点火过程中的火焰传播特性,建立了矩形模型燃烧试验系统,开展了流动与点火过
程中火焰传播特性、当量比对周向点火过程影响规律的研究。研究表明:初始火核主要在中心回流区边界形成,
其传播过程具有阶段性,先顺流线方向进行传递,待燃烧强度增天到一定值后,火焰才开始向其他区域延伸;在 相同压损下,当量比减小会使得燃烧室周向点火时间增加,主要为初始火焰发展阶段时间增加,当量比由0.89减
阮玲玉电影小到0.6时,周向点火时间由33.5 ms 增加到42 ms,流动加速系数由1.75增加到3.77,即随当量比减小,湍流褶皱 等流动因素对火焰传播的加速效应增天;在初始火焰的发展阶段,气体热膨胀效应对火焰传播加速效应占主导 作用。
关键词:航空发动机;三头部燃烧室;流场结构;点火过程;火焰传播
中图分类号:V231.2 文献标志码:A 文章编号:1005-2615(2()21)01-0150-08
Study on Flame Propagation Characteristics of a Triple-dome Combustor
During Ignition Process
FENG Jianhan 1, LANG Xudong , JIN Wu l , REN Yongzhi 1, LI Jianzhong 1
(1. College of Energy and Power ,Nanjing University of Aeronautics & Astronautics ,Nanjing 210016,China ; 2. Hunan
Aviation Powerplant Rearch Institute ,Aero Engine C-orporation of China, Zhuzhou 412002,China )
Abstract: A ries of optical experiments were conducted to investigate the effect of flow filed on the flame
ignition and propagation process between adjacent domes in triple-dome combustor. Results show that the
initial fire core is formed mainly at the boundary of the recirculation zone ,and the ignition process is staged. Initial flame core develops in the direction of local flow and it starts to propagate into other regions with the increa of combustion intensity. Under the same total pressure loss ,the combustor's circumference ignition
time increas when the equivalence ratio decreas ,especially in the time of the initial flame development.微波炉烤花生米
As the equivalence ratio decreas from 0.89 to 0.6,the circumference ignition time increas from 33.5 ms to
42 ms. During the circumferential ignition flame propagation ,the coefficients of flow-induced acceleration increas with the decrea of equivalent ratio. When the equivalent ratio decreas from 0.89 to 0.6,the
coefficients of flow-induced acceleration increas from 1.75 to 3.77,which means the turbulent wrinkle and other flow factors accelerate the flame propagation ; in the initial flame development stage , the thermal
expansion of gas dominates the flame propagation acceleration effect.
Key words : aero-engine ; triple-dome combustor ; flow field structure ; ignition process ; flame propagation
随航空工业的迅速发展,对航空发动机也提出
了越来越高的要求,而航空发动机燃烧室作为航空 发动机的三大核心部件之一,其性能直接影响整个
发动机的性能。其中,点火启动过程关系到飞行器
的地面起飞与高空再点火等性能,对航空发动机运
行的可靠性和安全性至关重要L1J o 现代航空发动
第1期冯剑寒,等:三头部燃烧室点火过程中火焰传播特性研究151
机多使用航空煤油作为燃料,燃烧室长度也越来越短,燃料驻留时间减少,这对点火性能提出更大挑战,因此,缩短点火延迟时间,探究点火过程中火焰传播的机理对飞行器的安全至关重要。
微雨燕双飞
点火是一个综合了物理和化学反应的复杂过程,影响点火过程及周向点火时间的因素错综复杂,国内外进行了大量的研究。Bourgouin等曰使用丙烷和空气作为预混气,对实验室尺度的燃烧室模型的点
火与火焰传播进行了机理研究,鉴别出周向点火的不同发展阶段特征。Philip等内采用大涡模拟(Large eddy simulation,LES)计算重现了Bourgouin等的实验,并且将周向点火过程的不同阶段进行了进一步的细化分析。Cordier等[4-5」对多头部矩形燃烧室模型进行研究,发现随着喷嘴间距的变化,相邻喷嘴间将出现“展向传播”和“混合传播”两种火焰传播模式。令狐昌鸿等在环形旋流燃烧实验平台上也观察到了类似的“混合传播”周向火焰传播模式。但对于多头部和环形燃烧室点火过程的研究,目前多采用实验室尺度的模型,而实际燃烧室加入旋流器、主燃孔、掺混孔等结构后,流场 结构变得更加复杂,火焰传播过程也更加复杂,因此,使用先进激光技术对燃烧室流场结构进行测量,结合流场对点火过程进行分析具有重要的意义。
Prieur等⑺研究了不同当量比及流速下的周向点火时间,发现周向点火时间随当量比的增大而减小,同时他还比较了3种不同燃料的周向点火时间,发现不同燃料特性对点火时间存在影响。Bourgouin等日通过数值模拟和实验的方法研究了流速对周向点火时间的影响,发现周向点火时间随流速增大而减小。众多研究表明当量比与燃烧室周向火焰传播机理有着密不可分的关系,但对其定量关系的研究依然十分缺乏,需要进行进一步的探究。
本文针对三头部燃烧室,使用二维粒子图像测速仪(Particle image velocimetry,PIV)测量燃烧室内不同截面位置的冷态时均流场,并使用高速摄像机拍摄燃烧室不同当量比下的点火过程。结合冷态时均流场与点火过程测量结果分析流动与点火过程中火焰发展的关联,以及当量比对周向点火过程的影响。
1试验系统及方法
如图1所示为燃烧室的整体试验系统图,建立燃烧室的试验系统主要包括:燃烧室、供气系统、供油系统、点火系统及光学测量系统。气源由最大流量为1.5kg/s,供气压力在0.005〜0.07MPa范围内可调的罗茨风机供给,然后通过5路支管分别对燃烧室进行供气,在每一路支管上都安装有阀门、涡街流量计以及浮子流量计,能够实现对每一支管流量的精准控制;同时在燃烧室各进口处均开有测压管,通过与水排连接来对进口总压进行测量。试验采用RP-3液态航空煤油作为燃料,通过3条油路对供油分别进行控制,每一条油路上都安装有油泵以及压力表,试验选取的压力表量程为0〜1.4MPa。试验采用高能点火器进行点火,点火能量为0.4J,点火频率为2Hz o
试验使用PIV来测量流场结构,采用平均直径为10“m的玻璃微珠作为示踪粒子,喷注位置设置在进气总管上,通过长距离的混合来保证示踪粒子分布的均匀性。试验采用的激光器为Nd:YAG双脉冲激光器,激光波长532nm,最大工作频率为15 Hz,数字相机为Bobcat B2041型,同步控制器为Micro Pul752型,保证激光器、数字相机和图像采集器之间的协调工作。
试验使用高速相机来记录点火和火焰传播发展的过程,相机分辨率为1280像素X1088像素,拍摄频率为1000Hz,试验时通过相机对焦至燃烧室特定截面位置来记录燃烧过程化学自发光成像的发光强度值。
试验研究的三头部燃烧室矩形件如图2(a)所示,模型采用分布式进气的方式,通过头部、上下主燃孔、上下掺混孔单独进气来实现对各路流量的精准控制。如图2(b)所示燃烧室3个旋流器分别安装在旋流器安装板上,各旋流器间中心间距为60mm,燃烧室采用喷雾锥角为78°的气动雾化喷嘴,安装在3个旋流器的中心位置,点火电嘴分别布置在中间和左侧旋流器位置,其中心距离旋流器出口7.5mm,其放电端面距离燃烧室上壁面2mm。在燃烧室主体结构上,布置了主燃孔和掺混孔,它们与旋流器之间的比例分别为1:1和1:3,其分布方式均正对旋流器中心位置。为了实现燃烧室轴向平面和周向平面的光学测量,在燃烧室的侧面布置有观察窗1,其高度与燃烧室平齐,长度为110m m,能够对燃烧室主燃区、中间区及掺混区上游轴向截面的流场与火焰结构进行测量,在燃烧室末端布置有观察窗2,它能对燃烧室整个周向截面的流场与火焰结构进行测量。燃烧室出口设置在模型的侧面。
神笔马良主要内容PIV与点火试验的测量截面及点火位置如图3所示,L1为中心对称截面丄2为点火位置截面,点火位置1布置在中间头部,点火位置2布置在左侧头部,点火模式为SFFL模式,即先启动点火器形成电火花后再通燃气。
152
南京航空航天大学学报第53卷
Computer of P IV
文夫blowe
六月雪叶图1试验系统图
Fig.1 Schematic of experimental system
Primary hole inlet (a) Combustor model
(b) The direction of air flow
图2三头部燃烧室模型
Fig.2 Triple-dome combustor model
aition position
* o o
Dilution hole
Ignition position
High-speed ca mera
High-speed camera
Fig.3图3测量截面及点火位置示意图
Schematic diagram of measuring ctions and ignition positions
2结果与讨论
2.1燃烧室点火过程分析
如 图 4 所示为 燃烧室中 心 对称截面点 火过程 示意图(点火器位置用白色虚线标出),进口为常温 常压,总压压损系数d =1.5%,点火油气比FAR=
0.05(当量比卩=0.744),如图4(a)所示点火过程的
原始图像以灰度图像的形式输出,为了突出点火过
程中火焰发展的细节以及边界轮廓,对图像进行了
归一化伪色彩处理(图4(b )),图中标尺代表图片 任意像素的发光强度与图片最大发光强度之比I o / I ”””,量纲为一。同时将点火前冷态时均流场矢量 图与火焰图像相叠加,可反映出燃烧室点火过程中
流动对火焰发展的影响。从图中可以看到:电火花
第1期冯剑寒,等:三头部燃烧室点火过程中火焰传播特性研究153
放电加热附近混气,形成了初始火核(t=1ms);由于电火花放电能量被周围未燃混气所吸收,主燃区附近燃油蒸发也需要吸收热量,初始火核的发光强度首先会减小,同时,因为中心回流区附近燃油与空气掺混效果最佳,初始火核的位置由点火器附近移动到了中心回流区的上侧,并在这里形成了点火热源(t=3〜5ms);点火热源不断放热,此时虽然主燃区附近火焰燃烧强度很小,但是未燃混气的温度在不断升高(t=5〜9ms);中心回流区附近区域开始发生化学反应,形成了初始火焰(t=9〜11ms);初始火焰在中心回流区高度以内向径向传播,并沿回流区流线向旋流器附近进行轴向传播(t=11〜13ms);初始火焰在主燃区的燃烧强度提高,同时部分火焰越过主燃孔射流传播到了中间区部分,主要分布在中间区的局部回流区附近(t=13〜15ms);火焰燃烧强度迅速提高,此时火焰主要分布在主燃孔射流上下
两侧的区域,因为此时主燃区旋流器附近已经处于富油状态,而主燃孔射流带入的空气形成快速猝熄作用使得这部分区域最接近当量比(t=15〜17ms);火焰持续扩大,但由于上下主燃孔和上下掺混孔交汇处会产生较大的速度梯度,此区域内气流负向速度增加,燃油驻留时间减小,所以火焰扩张趋势呈环形,其结构类似"C"形(t=17〜19ms);随着火焰不断扩大,化学反应速率不断提高,燃油蒸发率及未燃混气温度大幅增加,又将进一步促进化学反应速率的提高,在这一对正反馈作用下,燃烧室点火成功(t= 19〜27m s)。
图4L1截面点火过程示意图Fig.4Ignition process at L1
由于可近似认为点火过程中火焰发光强度反映燃烧反应的瞬时放热速率,绘制出点火过程中的发光强度积分随时间变化的曲线如图5所示。可以发现燃烧室点火过程在发光强度积分曲线上有阶段性的变化,大致分为以下3个阶段,通过积分曲线能够更加准确地判断点火时间:
(1)点火热源形成(t=1〜6ms):初始火核能量被未燃混气及燃油蒸发所吸收后,燃烧强度迅速减小,在中心回流区边界处形成点火热源。
(2)初始火焰形成(t=6〜11ms):在点火热源的不断加热下,中心回流区内开始形成了初始火焰,这段时间燃烧强度基本不变。
(3)燃烧室点火成功(t=11〜27ms):分为两个阶段,第一阶段初始火焰燃烧强度很小,火焰仅在中心回流区范围内发展,曲线增加速率很低;第二阶段燃烧强度增加到一定大小后,部分火焰传递到中间区,并呈“C”形环形扩张,曲线增加速率迅速增大,最终燃烧室点火成功,之后随下游火焰扩散出燃烧室,亮度积分曲线下降至稳定。
如图6所示为点火位置截面周向点火过程示意图,压损为1.5%,点火油气比为0.05(当量比0.744),白色虚线标示为旋流器中心位置,该截面能够更加清晰地反映三头部燃烧室周向点火过程中火焰的发展过程。从图中可以看到:
电火花放电
154南京航空航天大学学报第53卷
1
、
雨果的简介.^S U
B
U
二q
w)=
leJSEUI
图5发光强度积分随时间变化曲线图(L1截面) Fig.5Variation of integrated chemiluminescent intensity of flame with time(L1)
加热附近混气,形成了初始火核(t=1ms);未燃混气与燃油蒸发吸热,初始火核燃烧强度减小,在点火器位置附近形成点火热源(t=3〜5ms);点火热源不断放热,未燃燃气温度与燃油蒸发速率不断升高,达到化学反应边界,在点火热源附近形成了初始火焰(t=5〜9ms);初始火焰一方面沿径向传播,这部分火焰燃烧强度增加,另一方面也沿周向向右传播到两旋流器之间的部分,火焰基本没有向旋流器左侧传播,可以发现,燃烧强度较低的初始火焰更易顺流线方向进行传递(t=9〜11ms);火焰在两旋流器
之间部分迅速扩张,由于 中间旋流器中心形成的局部回流区的旋向为顺时针方向,其上侧流线与火焰边界平行,所以火焰延伸了过去,但其下侧流线与火焰传递方向相反,所以火焰并没有传递过去(t=11〜15ms);火焰燃烧强度继续增加,达到一定大小时,火焰开始向左侧旋流器左侧区域传递,直到左侧旋流器点火成功(t=15〜21ms);中间旋流器上侧的火焰向右传递到了右侧旋流器中心,同时其下侧也被点燃,这时中间旋流器点火成功(t=21〜27ms);此时火焰燃烧强度高,火焰传播受流场结构的影响较小,由已经传递到右侧旋流器中心的火焰在右侧旋流器内呈环状扩散传播,最终燃烧室整体点火成功(t=27〜37ms)o
图6L2截面点火过程示意图(气比=0.05) Fig.6Ignition process at L2(FAR=0.05)
如图7所示燃烧室左侧周向点火过程在发光强度积分曲线上也有阶段性的变化,大致分为以下5个阶段,由于周向点火过程也包含了单个头部点火过程,所以前3个阶段是一致的:
(1)点火热源形成(t=1〜5ms):初始火核能量被未燃混气及燃油蒸发所吸收后,燃烧强度迅速减小,并在点火器附近形成点火热源。
(2)初始火焰形成(t=5〜9ms):随未燃燃气温度与燃油蒸发速率不断升高,达到化学反应边界,形成初始火焰,这段时间火焰燃烧强度变化不大。
(3)左侧旋流器点火(t=9〜21ms):初始火焰顺流线方向传递到两旋流器之间的区域,并迅速扩张,燃烧强度也同时迅速增高,直到左侧头部点火成功。
(4)中间旋流器点火(t=21〜27ms):左侧旋流器点火成功后,随下游火焰扩散出燃烧室,左侧区域亮度下降,而中间旋流器还未完全点燃,此时曲线出现一个驻点,燃烧强度一段时间内不再增加甚至减小,之后再继续增加,直到中间旋流器点火成功。
(5)右侧旋流器点火(t=27〜37ms):中间旋流器点火成功后,与上一过程相同,曲线又出现一个驻点,之后再继续增加,直到右侧旋流器点火成功,燃烧室整体点火成功后,随下游火焰扩散出燃烧室,亮度积分曲线下降至稳定。
©
I
/
应
曽
甘
二
q
曰
二
旨3
U I
图7发光强度积分随时间变化曲线图(L2截面) Fig.7Variation of integrated chemiluminescent intensity of flame with time
(L2)