0 引言
目前,国内外学者已对燃料的层流燃烧特性进行了诸多研究。HE等[7]研究了在汽油中加入氢气对燃料层流燃烧特性的影响,实验结果显示当量比较低时,火焰传播速度受拉伸率影响显著。在保持当量比不变的前提下,随着氢气比例的增加,火焰传播速度也随着变大,并且LBV的峰值随氢气含量的增加向浓混侧偏移。XU等[8]通过分析发现稻壳油主要成分是乙醇、乙酸乙酯和丙酮,并研究了该燃料的LBV,结果表明该燃料的LBV峰值发生在当量比为1.1的时候。VAN TREEK等[9]考察了初始温度为358 K、含水量为0% ~ 40%、当量比为0.7 ~ 1.4的乙醇-水-空气混合物的LBV,发现加水会影响火焰传播过程。这是由于水在分子间反应或直接参与双分子反应中具有较高的碰撞效率。WANG等[10]测量了不同氢气比例下甲烷/二甲醚的LBV。研究表明,预混气体的LBV并不随氢气的增加而单调增加,甲基基团在反应中的作用比其他反应更重要,这表明氢的加入并不是使LBV增加的唯一途径。COMANDINI等[11]研究了正庚烷和氢混合物的燃烧特性,混合物中含有的氢气比例为0%、25%和50%,初始温度为294 K,初始压力为0.1 MPa,当量比在0.8 ~ 1.35之间。实验结果表明,添加25%和50%氢气的混合物的火焰速度分别增加了3%和10%,而马克斯坦长度略有下降。
为提供新的乙醇-氢气-空气混合物层流燃烧特性数据,本文采用高速纹影摄像系统记录定容燃烧弹内火焰传播过程,通过处理得到层流火焰的燃烧特性参数。本次实验主要研究初始温度为400 K、初始压力为0.1 MPa和0.4 MPa、当量比为0.7 ~ 1.4、氢气比例为0% ~ 90%条件下预混气体的层流燃烧特性。
1 实验装置及数据处理
1.1 实验装置
图1为实验装置图,主要由定容燃烧弹、高速纹影摄像系统、点火控制系统和数据存储系统组成。
Fig. 1 The experimental setup图1 实验装置图 |
1.2 数据处理
图2为使用HALCON软件对火焰前锋半径的处理结果。火焰图像处理步骤如下：（1）强化火焰原始图片的对比度;（2）通过thresshold算子对图像进行分割,使用设置阈值的手段从输入图像中选择满足灰度值条件的像素;（3）通过connection算子连接和计算区域,连接确定区域中给定的输入区域;（4）选择ROI,对已有的区域进行筛选;（5）使用select_shape算子,借助形状特征筛选出实验需要的火焰轮廓;（6）提取火焰轮廓并计算轮廓内的像素点个数。
Fig. 2 Flame front radius图2 火焰前锋面 |
通过式（1）可以将提取的像素点个数转化为等效火焰半径。
${{r}_{\text{u}}}=\sqrt{\frac{{{A}_{\text{sch}}}}{A}}\times r$ (1)
式中：ru为实际火焰半径;Asch为纹影中火焰面积;A为纹影中视窗面积;r为实际视窗半径。
火焰在实际传播过程中会受到拉伸率的作用,通过对火焰半径ru与火焰传播时间进行微分可以计算出火焰在受拉伸作用下的传播速度Sb,具体计算公式如下：
${{S}_{\text{b}}}=\frac{d{{r}_{\text{u}}}}{dt}$ (2)
拉伸率α是火焰前锋面面积变化的百分比速率,定义为火焰前锋面上一个无限小微元面积A的对数对时间t的导数,推导过程如下：
$\alpha =\frac{d\ln A}{dt}=\frac{2}{{{r}_{\text{u}}}}\times \frac{d{{r}_{\text{u}}}}{dt}=2\times \frac{{{S}_{\text{b}}}}{{{r}_{\text{u}}}}$ (3)
采用非线性外推法可以消除火焰α对拉伸火焰速度Sb的影响,从而得到层流状态下的无拉伸火焰传播速度。
${{\left( \frac{{{S}_{\text{b}}}}{S_{\text{b}}^{\text{0}}} \right)}^{2}}\ln {{\left( \frac{{{S}_{\text{b}}}}{S_{\text{b}}^{\text{0}}} \right)}^{2}}=-2\frac{{{L}_{\text{b}}}\alpha }{S_{\text{b}}^{\text{0}}}$ (4)
式中：Lb为马克斯坦长度;$S_{\text{b}}^{0}$为无拉伸火焰速度。
根据上面所求的无拉伸火焰速度和火焰燃烧前后的密度比,可以计算出混合燃料的LBV：
${{S}_{\text{u}}}=S_{\text{b}}^{0}\frac{{{\rho }_{\text{b}}}}{{{\rho }_{\text{u}}}}$ (5)
式中：ρb和ρu分别为已燃气体与未燃气体的密度;Su为混合气的LBV。
考虑到辐射的影响,根据YU等[12]提出的辐射校正公式对Su进行修正。
${{S}_{\text{uF}}}={{S}_{\text{u}}}\text{+}0.82{{S}_{\text{u}}}{{\left( \frac{{{S}_{\text{u}}}}{{{S}_{0}}} \right)}^{-1.14}}\left( \frac{T}{{{T}_{0}}} \right){{\left( \frac{P}{{{P}_{0}}} \right)}^{-0.3}}$ (6)
式中：S0为参考燃烧速度,取10 mm/s;T0为参考温度,298 K;P0为参考压力,0.1 MPa。
火焰厚度可以评估火焰传播过程中的流体动力学不稳定性,用如下公式进行计算：
${{\delta }_{\text{L}}}=\frac{{{T}_{\text{ad}}}-{{T}_{\text{u}}}}{{{\left( dT/dX \right)}_{\max }}}$ (7)
式中：Tad为绝热火焰温度;Tu为未燃混合气体温度;(dT/dX)max为火焰传播过程中温度梯度的最大值。
2 结果与讨论
Fig. 3 Schlieren images of the ethanol-hydrogen-air premixed laminar flames at different equivalence ratios and hydrogen fractions: (a) Xh = 10%; (b) ϕ = 1.2图3 乙醇-氢气-空气预混层流火焰在不同当量比和氢气分数下的纹影图像：（a）Xh = 10%;（b）ϕ = 1.2 |
2.1 层流燃烧速度
图4a和图4b分别展示了初始温度为400 K、初始压力为0.1 MPa和0.4 MPa条件下,6种氢气比例的混合燃料LBV受当量比变化的影响。由图中可以明显地看出,火焰LBV随着当量比的增大先增加后减小,峰值出现在当量比1.2附近。随着氢气含量的增加,火焰的LBV随之变大。特别在富氢情况下,火焰的LBV有着非常大的提高。在图4a中,Xh = 90%比Xh = 70%的最大LBV快1 080 mm/s,而Xh = 30%比Xh = 10%的最大LBV快51 mm/s,这是由于氢气的化学反应速率更快,在富氢状态下LBV受氢气含量变化的影响更加显著。此外,随着初始压力的增加,LBV明显降低。初始压力增加使得H自由基浓度影响变小,减弱了氢气添加对燃料LBV的促进作用。
Fig. 4 The relationship between LBV and equivalent ratio of premixed fuels with different hydrogen ratios: (a) 0.1 MPa; (b) 0.4 MPa图4 不同氢气比例的预混燃料LBV随当量比的变化图：（a）0.1 MPa;（b）0.4 MPa |
图5为当量比1.2时,不同初始压力下预混火焰LBV增长率随氢气比例的变化。从图中可以看出,当Xh ≤ 50%,LBV随氢气含量的变化线性增加,而当Xh ≥ 50%时,LBV随氢气含量的变化指数增长。这说明氢气含量越高,LBV增加效果就越明显。当初始压力上升时,可以观察到增长率也在上升。虽然初始压力的上升会导致LBV下降,但会增加LBV的增长率。
Fig. 5 Plot of the growth rate of the LBV of premixed fuels at different initial pressures versus hydrogen content图5 不同初始压力下的预混燃料LBV增长率随氢气含量的变化图 |
2.2 火焰厚度
火焰厚度是表征火焰流体动力学不稳定的重要参数。图6为不同氢气比例下乙醇-氢气混合燃料火焰厚度随当量比的变化曲线。
Fig. 6 The relationship between flame thickness and equivalent ratio: (a) 0.1 MPa; (b) 0.4 MPa图6 火焰厚度随当量比变化关系：（a）0.1 MPa;（b）0.4 MPa |
由于混合物的运动学黏度对当量比的变化不敏感,火焰厚度表现行为与LBV正好相反。当Xh ≤ 50%,火焰厚度最低点出现在当量比为1.2时。而当Xh ≥ 50%,火焰厚度的最低点向当量比较高的一侧偏移。Xh = 90%时,最低点在当量比1.3时出现。相比于贫氢状态,富氢状态下氢气比例的增大对火焰厚度的影响更加显著。随着氢气比例和初始压力的增加,火焰厚度变得更薄,而较薄的火焰厚度会使火焰前锋面两边的密度梯度变化增大,使得火焰表现出不稳定的趋势。
2.3 马克斯坦长度
热扩散不稳定是由预混燃料中热量和质量的不均等扩散燃烧引起[6]。图7中的Lb反映了拉伸率对火焰传播速度的影响,可以量化火焰的热扩散不稳定性[8]。当Lb > 0时,火焰前锋面结构稳定,而当Lb < 0时则表明火焰前锋面结构不稳定[13]。如图7所示,随着初始压力的增加,Lb呈现降低趋势,火焰不稳定性增强。当氢气比例较小时,Lb随当量比的增加而降低,这说明火焰传播过程中火焰稳定性逐渐降低。当处于富氢状态下,则表现出相反的趋势,火焰稳定性不断上升。由此可以发现氢气比例对改变火焰稳定性有很大的影响。特别是在稀燃侧,贫氢和富氢状态下的Lb有着较大的区别,而在当量比为1.2 ~ 1.4的富燃侧,Lb的区别不大。
Fig. 7 The effect of Lb with equivalence ratio at different hydrogen contents图7 不同氢气含量下Lb随当量比的变化图 |
2.4 数值仿真结果
通过热力学仿真软件Chemkin 18.0 Pro分析计算得到实验中所不能直接测量的绝热火焰温度、已燃气体和未燃气体密度等数据。计算网格长度为3 cm,曲线参数和梯度参数均为0.05,从而保证了计算的准确性。本研究采用的化学机理是通过整合Ó CONAIRE等[14]的氢燃烧反应机理和MARINOV[15]的乙醇氧化反应机理得到。图8a展示了数值仿真与实验所获得的火焰LBV对比,可以看出两者的发展趋势一致。仿真结果与实验结果的平均相对偏差小于12%,最大相对偏差小于25%。图8b为乙醇预混燃气的LBV与文献的实验与仿真数据对比情况。结果表明,本研究采用的仿真机理与文献报道的仿真数据差距较小,趋势规律基本吻合;本实验最终得到的乙醇预混燃气的LBV相较于BRADLEY等[16]的实验结果更大,考虑到本研究初始温度为400 K,高于其研究工况,温度越高LBV越大;另外与DIRRENBERGER等[17]研究数据结果基本吻合,差距在4.8%以内,因此本文实验数据可靠性较高。
Fig. 8 Comparison of experimental and simulation results of LBV for ethanol with different H2 addition (a); the comparison of present data with experimental and computed data published in previous literature (b)图8 不同氢气比例的混合燃料LBV实验数据与仿真结果的对比（a）;现有数据与文献中实验数据和计算数据的对比（b） |
3 结论
研究了含有氢气比例分别为0%、10%、30%、50%、70%、90%的乙醇-氢气燃料的燃烧特性。实验采用定压法在初始温度为400 K、初始压力0.1 MPa和0.4 MPa下进行,考察了燃料当量比在0.7 ~ 1.4区间内的层流燃烧特性,当量比跨度为0.1。主要结论如下：
（1）乙醇-氢气-空气预混火焰的层流燃烧速度峰值出现在当量比1.2附近,并且随着氢气比例的增大,峰值当量比会朝着富燃一侧移动。
（2）在相同工况条件下,氢气的添加可以提高乙醇的层流燃烧速度,特别是在富氢状态下提升明显。贫氢阶段的LBV随氢气比例的增加线性增长,而富氢阶段的LBV随氢气比例的增加呈指数增长。此外,初始压力的增加提高了LBV的增长率。
（3）在贫氢阶段,马克斯坦长度表现出由正转负的下降趋势,在富氢阶段,马克斯坦长度表现出由负转正的趋势。这表明预混火焰在贫氢阶段处于稳定状态但是抗失稳的能力随着当量比的上升而不断下降,在富氢阶段则表现出相反的趋势。
（4）在工况改变的情况下,本文仿真模型较好地预测了层流燃烧速度的变化趋势,计算结果与文献报道的不同氢气含量下乙醇-氢气混合气的实验结果一致。
