回流区稳定燃烧的近熄火特性理论分析1)

汪洪波 连城阅 张锦成 曾 宇 杨揖心 王亚男

(国防科技大学空天科学学院高超声速技术实验室,长沙 410073)

引言

在高速气流中制造低速回流区是一种最常用的火焰稳定方式,而回流区的产生既可以采用侵入式的钝体[1-2]或支板[3],也可采用非侵入式的台阶[4]或凹腔[5-7].回流区与高速自由流之间形成一个剪切层,剪切层通过同时从回流区和自由流卷吸流体而逐渐增长,由于低速回流区不停地从剪切层中卷吸燃烧产物并将其输运到上游并促进新鲜反应物燃烧,火焰可以被稳定在剪切层中[8-9].为了降低燃气涡轮等燃烧装置的排放(特别是NOx的排放),可以通过使燃料与大量空气掺混来实现贫燃设计,此时燃烧装置通常工作在低当量比条件,有时甚至非常接近贫燃吹熄极限[10].因此,回流区稳定燃烧的近熄火特性是一个非常重要的问题.

一般通过分析流动驻留时间和化学特征时间能够判断熄火是否发生,而当量比是决定化学特征时间的关键参数,因此当量比是影响贫燃预混火焰熄火极限的主导因素.采用凹腔主动喷注的燃烧室中,燃料在凹腔里预混更充分,这种主导性更强.Barnes等[11]采用平面激光诱导荧光(planar laser-induced fluorescence,PLIF)技术研究了凹腔稳焰燃烧室中的燃料分布,发现尽管全局当量比是贫燃的,但回流区内基本上是富燃的,凹腔局部当量比决定着火焰的稳定性.Rasmussen 等[12-15]对碳氢燃料的火焰稳定性进行了实验研究,发现贫燃熄火极限取决于空气质量流率、凹腔几何构型、燃料喷注方案、马赫数和燃料类型.其中喷注位置对贫/富燃熄火极限都有很大的影响,这表明将凹腔建模为一个良好的搅拌反应器具有一定局限性.基于这一结论,Driscoll 等[16]通过Da数引入高速气流中凹腔稳定的火焰非预混特性和熄火极限的关联,建立了预测富燃和贫燃吹熄极限边界的经验公式.

对于凹腔被动燃料喷注,由于射流和凹腔之间的复杂相互作用,当量比对火焰熄火机制的主导作用将受到其他机制的影响.Owens 等[17]采用凹腔上游壁面喷注方案,研究了凹腔稳焰模式下火焰的稳定性.研究发现,局部的当量比是影响火焰稳定性的关键因素,但来流滞止温度可改变当量比对熄火极限的影响规律.Zhang 等[18]对超声速燃烧室中超临界煤油的熄火极限进行了实验研究,发现来流空气滞止温度和燃料喷注位置对熄火极限有很大影响:贫燃熄火极限随着滞止温度的增加而增加.Wang 等[19]研究了超声速气流中凹腔上游喷注煤油诱导的火焰稳定性,发现贫燃极限随着总压的增加或者总温的降低而增加,且总温的影响更为显著.此外,回流区内热产物对混合气体的预热也会使火焰熄火极限改变,但其作用机制复杂且相关研究较少,目前当量比仍然是受主要关注的熄火极限控制变量.

对于回流区稳定的贫燃火焰,其吹熄或近吹熄行为有很多相关的研究,包括实验研究[20-21]、数值模拟[22-24]及基于反应器网络模型的分析[25].Shanbhogue等[1]对钝体稳定火焰的吹熄特性做了系统的分析和总结,认为吹熄过程包括几个阶段: 火焰面的局部熄火、尾迹区的大尺度分裂、尾迹冷却及收缩触发的最终吹熄.近来,Fugger 等[26]和Kumar 等[27]还分析了展向非均匀性及三维效应对吹熄特性的影响.然而贫燃火焰燃烧产物中过剩的氧化剂对剪切层火焰的影响却鲜有报道.一旦火焰稳定在剪切层中,由于剪切层同时从两侧卷吸流体,所以进入火焰的流体将是自由流与回流区流体的混合物.如果是贫燃火焰,回流区流体由燃烧产物和多余的氧化剂组成而没有燃料,因此进入火焰的混合物当量比将低于自由流的当量比.如此一来,即使自由来流是可燃的,剪切层中混合物的当量比也可能超出可燃极限.本文将重新审视回流区稳定燃烧,力图通过理论分析揭示这些近吹熄火焰的一些基本特性.

1 理论模型

不失一般性,考虑一个凹腔稳定的近贫燃预混火焰,如图1 所示.由于在稳定火焰和减少总压损失方面表现俱佳,凹腔火焰稳定器广泛应用于驻涡燃烧室 (TVC)[28]和亚燃/超燃冲压发动机燃烧室[29-32].由于自由流是贫燃预混反应物,其当量比ϕ=＜1,其中Yf和Ya分别为燃料和空气的质量分数,rs为当量燃料−空气比.为了推导分析回流区稳定贫燃火焰的模型,做如下假设.

图1 凹腔稳定贫燃预混火焰示意图Fig.1 Schematic of a cavity-stabilized lean premixed flame

(1) 火焰稳定于回流区与自由流形成的剪切层中[16].对后向台阶形成的回流区,Coats 等[33-34]发现,在很宽的条件范围内,主要的燃烧均以一个抬举火焰的形式发生于距台阶一定距离的回流区剪切层内.此外,随着距凹腔前缘距离的增加,剪切层中的应变率减小而反应物和热回流流体的接触距离增加,因此在近吹熄时火焰很可能稳定在剪切层的下游部分.Stöhr 等[35]的测量也表明贫燃火焰可能稳定在应变率较低的区域.

(2) 剪切层对自由流和回流区流体的卷吸发生在火焰稳定点之前,而剪切层流体进入回流区发生在火焰稳定点之后.Coats 等[33-34]的结果表明剪切层对两侧流体的卷吸主要发生在点火/稳焰点之前.

(3) 进入火焰的剪切层混合物是均匀的.基于假设 (1),可以假定剪切层在火焰稳定点之前已经经历混合转捩,剪切层中的湍流将使来自回流区和自由流的流体快速达到分子尺度混合[33-34].

(4) 回流区内是均匀混合物.亦即,回流区内原有的流体瞬间与来自剪切层的流体达到均匀混合状态,Kalt 等[36]的拉曼测量基本支持该假设.如果火焰稳定在剪切层内,进入回流区的流体仅包含贫燃火焰的产物和剩余空气;如果火焰吹熄,进入回流区的流体包含燃料和空气.

(5) 假设流场为准稳态的.亦即,卷吸进入剪切层和回流区的流体质量流率为常数,但是流体组分可随时间变化.需要注意的是,采用准稳态假设是为了简化分析并突出新的机制,即使在非稳态条件下本文所提出的机制仍然有效.

根据上述的准稳态假设,凹腔回流区内的总质量M为常数.若假定从自由流进入剪切层的质量流率为fr,从回流区进入剪切层的质量流率为rz,则从剪切层进入回流区的质量流率也为rz,从而回流区内空气质量分数的变化率可写为

其中Ya,rz表示回流区内空气质量分数,Ya,sh,af表示剪切层内火焰后的空气质量分数.当贫燃混合物经过火焰后,其中所有燃料和部分空气转化为产物,这里“产物”既包括燃料/氧气的化合物,也包括原来空气中与所消耗氧气对应的那部分惰性组分(如氮气),因此“空气”的成分和燃料−空气恰当当量比rs保持不变.因此火焰后方的空气质量分数应为

将式(2)～式(4) 代入式(1) 得

其中 α=/为剪切层卷吸比,τ=M/rz为凹腔回流区的驻留时间.方程 (5)的解为

其中a=−α/(α+1)/τ;b=(1−ϕ)/(1+ϕrs);c为由初始条件决定的一个常数.

很容易得到剪切层内火焰前方的当量比

值得注意的是,当t→∞ 时,Ya,rz→b且ϕsh,bf→αϕ/(α+1−ϕ).由于 α＞0 且 ϕ＜1,所以αϕ/(α+1−ϕ)＜ϕ,这表明进入火焰的混合物的有效当量比低于自由来流的当量比,且偏离的程度由卷吸比 α 和自由来流当量比 ϕ 共同决定.这意味着此时的贫燃吹熄极限当量比将高于预期值.对于接近富燃吹熄极限的火焰,存在类似情况.

一旦火焰被吹熄,可以发现当t→∞ 时,Ya,rz→1/(1+ϕrs) 且 ϕsh,bf→ϕ.

2 近贫燃熄火特性分析

为了简化分析,这里仅仅关注当量比对近吹熄特性的影响而暂时忽略其他因素,即认为当量比对可燃极限的影响起主导作用,而回流区产物对剪切层的预热作用相对较小,实际中可能对应于温升较小的火焰或回流区的热损失比较大的情况.因此,假定火焰的稳定或吹熄仅仅由剪切层内火焰之前的有效当量比决定.需要说明的是,若考虑产物对反应物的预热,则会出现有效当量比与预燃效应的竞争,问题会变得更复杂一些,但是不会从根本上影响本文的分析结论.

由于点火条件通常比火焰稳定条件更加苛刻[14,37-38],参数空间内将在火焰吹熄和再点火/再稳焰之间出现一个迟滞,形成我们熟悉的具有上点火分支、下熄火分支和不稳定中间分支的“S−曲线”[39],如图2所示.对于贫燃火焰,化学反应特征时间随当量比的增加而减小,若假定流动特征时间基本保持不变,则剪切层内的局部Da数随当量比的增加而增大.因此吹熄当量比 ϕblowoff应该比再点火当量比 ϕreignition低,即 ϕblowoff＜ϕreignition,Ahmed 等[37]的实验支持这一结论.根据前面的分析,剪切层内的有效当量比ϕsh,bf在稳焰期间将趋近于 α ϕ/(α+1−ϕ) 而在吹熄期间趋近于 ϕ,从而可得到回流区稳定燃烧在近贫燃条件下的4 种可能模式.

图2 具有多解的典型折叠S−曲线Fig.2 Representative folded S-curve with multiple solutions

(1) 若 αϕ/(α+1−ϕ)＜ϕblowoff,ϕ＜ϕreignition,剪切层中将不会出现火焰或者火焰最终被吹熄,剪切层中的当量比 ϕsh,bf最终趋近于自由来流的值 ϕ.这是因为 αϕ/(α+1−ϕ)＜ϕblowoff使得火焰最终将被吹熄而 ϕ＜ϕreignition意味着无法出现再点火/再稳焰.

(2) 若 αϕ/(α+1−ϕ)＜ϕblowoff,ϕ ≥ϕreignition,火焰稳定和吹熄交替出现,导致燃烧振荡.因为αϕ/(α+1−ϕ)＜ϕblowoff意味着火焰将被吹熄,而火焰吹熄后剪切层中的当量比 ϕsh,bf将趋近于自由来流的值 ϕ,所以 ϕ ≥ϕreignition又确保一段时间后可实现再点火/再稳焰.

(3) 若 αϕ/(α+1−ϕ)≥ϕblowoff,ϕ＜ϕreignition,燃烧是亚稳定的.此时,剪切层中可能存在稳定的火焰,因为 αϕ/(α+1−ϕ)≥ϕblowoff,但是一旦火焰被吹熄,无法实现再点火/再稳焰.

(4) 若 αϕ/(α+1−ϕ)≥ϕblowoff,ϕ ≥ϕreignition,燃烧是超稳定的.此时,剪切层中可以保持稳定的火焰,即使火焰被吹熄仍可实现再点火/再稳焰,且剪切层中的当量比 ϕsh,bf将趋于 α ϕ/(α+1−ϕ).

回流区稳定燃烧在近贫燃吹熄时的模式如图3所示.其中,仅有直线 αϕ/(α+1−ϕ)=ϕ 以上的部分才具有物理意义,因为对于 ϕ＜1 有 αϕ/(α+1−ϕ)＜ϕ.可以发现,这些模式由4 个参数决定: 自由来流当量比 ϕ、剪切层卷吸比 α、吹熄极限 ϕblowoff和再点火极限 ϕreignition.若 α,ϕblowoff和 ϕreignition保持不变,当 ϕ 从一个中等当量比逐渐减小时,燃烧系统将从一个超稳定火焰分岔为一个亚稳定火焰或振荡火焰,直至完全吹熄.对于近富燃吹熄的回流区稳定燃烧,可以得到一个类似的模式图.

图3 回流区稳定燃烧在近贫燃吹熄时的模式图Fig.3 Regimes of recirculation-zone stabilized combustion near lean blowoff

4 个参数(来流当量比、剪切层卷吸比、吹熄极限和再点火极限) 的典型范围,跟自由来流状态(温度、压力和马赫数等)、燃料种类、火焰稳定器构型(影响回流区形态或剪切层特性)等均有密切关系,因此很难给出一个通用的范围,只能给出一个近似估计,具体情况需要具体分析.根据文献调研和作者已有的研究认识,贫燃预混气来流当量比范围为0＜ϕ＜1;剪切层卷吸比与来流条件及具体稳焰器构型有关,一般在1 附近,对于超声速气流中的凹腔剪切层,通常范围为 0.5 ≤α ≤2.0[40];吹熄极限和再点火极限与燃料种类、稀释成分及具体流动条件有关,对于特定燃料和稀释成分,预混气存在一个可燃极限[39],由于气流的剪切耗散作用,回流区稳定燃烧的吹熄极限和再点火极限比该可燃极限更窄,且由于点火条件通常比火焰稳定条件更加苛刻[14,37-38],再点火极限比吹熄极限更窄,因此其范围为ϕl,b＜ϕblowoff＜ϕreignition＜1,其中ϕl,b为该燃料在相同温度、压力静止气中的可燃极限当量比下限.

值得一提的是,上述分析在参数空间中发现了一个振荡区域,这为近吹熄火焰引入了一种已有研究中不曾关注到的和新的不稳定机制.此外,这是一种不需要外界触发的、自发形成的本质不稳定性.这个发现与这样一个事实相符: 很多在远离吹熄条件时表现稳定的火焰在接近吹熄极限时往往出现振荡.Lieuwen 等[41]指出,若当量比振荡是由燃烧过程和压力振荡驱动的,那么这种当量比振荡可引起燃烧不稳定.在当前机制中,剪切层中的当量比振荡确实是由燃烧过程本身驱动的,因为火焰的位置/振荡影响着进入回流区的流体成分进而影响剪切层中的当量比.因此,当量比振荡与燃烧过程之间的反馈循环由于周期性的火焰吹熄和闪回/再稳定而得以封闭.唯一的区别在于,此时的反馈循环是由回流区的对流过程而不是传统的声波/压力波传播过程所封闭.所以,这里的不稳定可被认为是一种“对流−对流”不稳定性,而不是传统的“对流−声波”不稳定性.本文分析发现的这种不稳定性在某种程度上与Dhanuka 等[21]及Huang 等[24]观察到的“贫燃极限”不稳定类似,其中的振荡伴随着周期性的火焰抬举和闪回.

需要指出的是,上述的近贫燃吹熄振荡场景是一种理想的准稳态情况.实际上,如Shanbhogue 等[1]指出的那样,吹熄过程通常由多个步骤组成,且以火焰面上的局部熄火为第一步.在这种情况下,火焰完全吹熄之前将有未燃燃料被卷吸进入回流区,并在一段时间后导致剪切层中当量比的增加.若剪切层当量比的响应时间比剪切层中部分熄火火焰的存活时间要短,那么整个火焰可能被重新点燃/稳定.此后,没有燃料进入回流区,剪切层当量比降低,又导致局部熄火.这些过程重复发生,引起燃烧振荡.这种不稳定机制可能是钝体燃烧器在贫燃吹熄极限时出现局部熄火和再点火过程[42-43]的一个原因,不过还需要仔细确认和进一步分析.

文献中贫燃吹熄一般是通过逐渐减少燃料流率(即当量比)而实现的.根据上述分析,在此期间火焰应该逐步经历“超稳定”、“振荡”/“亚稳定”和吹熄几种模式,这与实验观察到的贫燃火焰行为[42-43]是基本一致的.在“超稳定”模式下,即使出现短暂的吹熄,仍能出现再点火/再稳焰,并表现为火焰振荡.因此,实际上一个“超稳定”的火焰可以是稳定的或者是振荡的(遇到足够大的扰动时).然而在“亚稳定”模式下,一旦吹熄将无法实现再点火/再稳焰,因此“亚稳定”模式可能很难从“吹熄”模式中分离.此外,虽然上述分析表明回流区对流可改变局部当量比,但Barlow 等[44-45]观察到的优先组分扩散也可能对有效当量比产生附加影响.因此,需要在以后的工作中采用仔细的设计和测量来进一步确认本文的理论分析结论及所提出的模式图.

3 结论

基于火焰稳定在剪切层中的假设,建立了近贫燃熄火极限的回流区稳定燃烧理论分析模型.分析表明,进入剪切层火焰的混合物有效当量比低于自由来流的当量比.根据理论分析建立了近贫燃熄火极限的回流区稳定燃烧的模式图,其中涉及4 个参数: 自由来流当量比、剪切层卷吸比、吹熄极限和再点火极限.理论分析揭示出4 种燃烧模式: 超稳定火焰、亚稳定火焰、振荡火焰和熄火.研究发现了参数空间存在一个振荡区,它可能为近吹熄火焰引入一种新的固有不稳定性,即一种由回流区剪切层内有效当量比振荡引起的“对流−对流”燃烧不稳定性,当量比振荡与燃烧过程之间的反馈循环由于周期性的火焰吹熄和闪回/再稳定而得以封闭.

由于近吹熄行为可能受多种机制影响(如声学、火焰−涡干扰和火焰闪回等等),未来工作中还需要精细的实验或数值模拟来确认本文的理论分析.