0 引言
近年来,诸多文献对含氢燃料混合物的LBV进行了广泛的报道。TUTAK等[13]研究了柴油-天然气-氢气混合替代燃料的燃烧过程,发现氢气的添加虽然使发动机出现了爆震,但是燃烧持续的时间以及燃烧末端的扩散等燃烧过程都得到了明显的改善。JAMROZIK等[14]在双燃料柴油机中对氢气与柴油混合燃料进行了排放与燃烧过程的研究,发现氢气的添加使得最大燃烧压力、热释放率和压力上升率都得到了提升。LI等[15]利用球形火焰传播实验确定了在大气压下不同乙醇比例对乙醇-空气混合物LBV的影响,并在现有实验的基础上建立了经验方程来估算乙醇-氢气-空气预混火焰的LBV。XU等[16]在乙酸乙酯中加入氢气,发现氢气延长了乙酸乙酯的可燃下限,大大提高了其LBV。LI等[17]采用球形火焰法对不同燃料比影响甲烷-氢气预混火焰LBV的因素进行综合分析,发现LBV的变化与热效应变化密切相关。SUN等[18]分析了氢气-空气预混火焰爆炸特性。结果表明,对于氢气-空气混合物,化学计量比的混合物具有最高的绝热火焰温度,即在相同的初始环境下,化学计量爆炸可能释放的热量最大。
综上所述,氢气与其他燃料混合的研究从未停止,其原因在于氢气的加入可以极大地提升混合燃料的LBV,还可以改善其排放特性,拓展燃烧极限。但是由于氢气存在点火能低、易爆炸等危险性,氢气的加入也会不可避免地使得整体燃料的安全性降低。因此,氢气的存储与安全性问题至今仍是热门研究方向。本研究利用Chemkin软件,探究初始压力为0.2 ~ 0.4 MPa、温度为400 K时的氢气-乙醇-空气预混火焰的LBV,结合绝热火焰温度和火焰厚度等参数对燃料特性进行综合分析,并对影响燃料LBV变化的因素进行深入讨论。本研究旨在探寻多元混合燃气的燃烧机理,为氢气混合物在存储和运输上的安全问题及氢混乙醇的实际应用提供理论指导。
1 计算模型和相关理论
1.1 计算模型
一维稳态火焰的控制方程可以用下面四个方程描述[19]。
（1）连续性方程：
$\dot{M}=\rho uA$ (1)
式中：$\dot{M}$、$\rho $、u分别为混合燃气的质量流率、密度、流速;A为预混燃气的火焰在平面的横截面积。
（2）能量守恒方程：
$\begin{align} & \dot{M}\frac{\text{d}T}{\text{d}x}-\frac{1}{{{c}_{p}}}\frac{\text{d}\left( \lambda A \right)}{\text{d}x}\frac{\text{d}T}{\text{d}x}+ \\ & \frac{A}{{{c}_{p}}}\left( \sum\limits_{k=1}^{k}{\rho }{{Y}_{k}}{{v}_{k}}{{c}_{p,k}}\frac{\text{d}T}{\text{d}x}+\sum\limits_{k}^{k}{{{{\dot{\omega }}}_{k}}}k{{h}_{k}}{{W}_{k}} \right)=0 \\ \end{align}$ (2)
式中：x为火焰在空间中的坐标;T为变化的温度;λ为混合燃气的导热率;cp为混合燃气的定压比热容;Wk、Yk、cp,k、hk、vk、${{\dot{\omega }}_{k}}$分别为混合燃气中第k种物质的分子质量、质量分数、定压比热容、比焓、扩散速度和化学反应速率。
（3）质量守恒方程：
$\dot{M}\frac{\text{d}{{Y}_{k}}}{\text{d}x}+\frac{\text{d}}{\text{d}x}\left( \lambda \rho {{Y}_{k}}{{v}_{k}}A \right)-A{{\dot{\omega }}_{k}}{{W}_{k}}=0$ (3)
（4）状态方程：
$\rho =\frac{p\bar{w}}{RT}$ (4)
式中：p为变化的压力;$\bar{w}$为混合燃气的平均分子量;R为理想气体常数。
1.2 火焰化学反应动力学机理
本研究使用Chemkin中的Premix模块进行模型搭建,如图1所示,通过添加气相动力学文件、表面动力学文件和传递过程文件对氢气-乙醇-空气预混燃料的燃烧过程和特性进行仿真研究。
Fig. 1 Model building of Chemkin图1 Chemkin模型搭建 |
2 结果与讨论
2.1 实验与仿真
2.2 层流燃烧速度
H2 比例分别为0%、50%、70%和80%,初始温度为450 K,初始压力为0.2 ~ 0.4 MPa时对氢气-乙醇-空气预混火焰LBV影响的变化规律如图3。当H2比例提高后,LBV峰值的出现逐渐向当量比高的一侧偏移。H2比例的提高会使得LBV明显上升,特别是在H2比例大于50%后。但随着初始压力的升高,H2比例的提高对LBV的影响开始减小。随着初始压力的升高,LBV逐渐受到抑制,在相同Φ下,LBV会随着初始压力的增大而减小。这是由于H2比例的上升使得基元反应中的H自由基浓度得到提升,加快了OH自由基的生成,同时也会使得链式反应放出更多的热量,从而加快了预混火焰的LBV。然而,增加初始压力会减小H自由基的浓度影响,从而降低氢气添加对燃料LBV增加程度的影响。
Fig. 3 Curve of laminar burning velocity changing with hydrogen proportion图3 层流燃烧速度随氢气比例变化的曲线 |
2.3 密度比
密度比 𝜎 定义为：
$\sigma =\frac{{{\rho }_{\text{u}}}}{{{\rho }_{\text{b}}}}$ (5)
式中：${{\rho }_{\text{b}}}$和${{\rho }_{\text{u}}}$分别为燃烧和未燃烧混合物的密度。
燃料在点燃之后,已燃区域的密度会发生改变,这将导致火焰在其前沿和后沿的密度不均匀的情况下发生拉伸,从而在火焰前锋面的气流形成流线偏折[19],使火焰单位体积燃烧速率发生变化,这一变化最终会使得火焰出现流体动力学不稳定。因此与流体动力学不稳定影响正相关的 𝜎 可以很好地反映火焰出现不稳定的程度。
图4展示了Φ、H2比例、初始温度和初始压力对于 𝜎 的影响。由图可知初始温度的变化可以显著影响 𝜎。升高初始温度,氢气-乙醇-空气预混合火焰的 𝜎 呈现明显的下降趋势。当Φ升高时,可以看出 𝜎 在Φ = 1.1左右达到峰值。当H2比例升高时,𝜎 会上升,但上升并不明显。而当初始压力升高时,预混火焰的 𝜎 几乎没有变化。
Fig. 4 Influence of different working conditions on density ratio图4 不同工况对密度比的影响 |
2.4 绝热火焰温度
图5是绝热火焰温度Tad 随着Φ、H2比例和初始压力变化的曲线。在Φ ˂ 1的阶段,氢气-乙醇-空气预混火焰的Tad几乎不变,在Φ = 1.0之后则开始出现明显的差异,并在Φ = 1.1时达到最高值。出现这种现象的原因是H2的摩尔热值相对于乙醇来说非常小,约为其五分之一,因此当H2比例上升时,总混合物的摩尔热值会下降。空气中的主要成分为N2,总混合物的比热主要由N2决定。因此当N2含量减少即Φ上升时,氢气-乙醇-空气预混火焰的Tad呈现上升的趋势。当Φ ≥ 1.0时,混合物可以进行充分的反应,使得Tad对H2比例变化的响应开始明显。初始压力的变化对氢气-乙醇-空气预混火焰的Tad的影响非常小,仅仅在化学计量比时有些许差距。而当预混火焰中H2的比例上升时,初始压力对其影响相比于低H2比例时会稍有提高。
Fig. 5 Influence of different working conditions on the temperature of adiabatic flames图5 不同工况对绝热火焰温度的影响 |
2.5 火焰厚度
$\delta =\frac{{{T}_{\text{ad}}}-{{T}_{\text{u}}}}{{{\left( \text{d}T/\text{d}R \right)}_{\max }}}$ (6)
式中：Tu为初始火焰温度。
图6为火焰厚度随着Φ、H2比例以及初始压力变化的曲线。由于混合物的运动黏度对Φ的变化不敏感,δ表现出来的行为与无拉伸火焰速度正好相反。在H2比例小于50%时,δ最低点出现在Φ = 1.2处,并且H2比例的上升并没有对δ产生明显的影响。而当H2比例大于50%的时候,δ最低点出现的位置向当量较高的一侧偏移,其峰值出现在Φ = 1.3时,H2含量对δ的影响也更为明显。初始压力的升高会使得δ变薄,进而使火焰前锋面两边的密度梯度变化增大,导致火焰表现出不稳定的趋势。
Fig. 6 Influence of different working conditions on flame thickness图6 不同工况对火焰厚度的影响 |
2.6 敏感性分析
火焰的温度与火焰的传播速度和化学反应速率存在线性关系,因此从当量比、初始压力、H2比例变化来对氢气-乙醇-空气预混火焰对温度的敏感性进行分析。
基元小分子反应一般作为碳氢化合物和含氧碳氢化合物的高温氧化反应的支配者,起到促进或抑制的作用,从而影响火焰传播速度。本文主要对前6个敏感反应进行分析。图7显示了氢气-乙醇-空气预混火焰在 Φ = 0.7、1.2和1.4时以及不同H2添加比例下的敏感反应。反应 H + O2 = O + OH (R1)是最敏感的反应,因其可以加速O 和 OH 等活性自由基的产生,并且R1的灵敏度随H2添加量的增加而增加。反应H2O + M = H + OH + M (R8)、HO2 + OH = H2O + O2 (R13)、H + O2 (+M) = H2O (+M) (R9) 和 CH3 + H (+M) = CH4 (+M) (R96) 的敏感性系数为负数,表明这些基元反应会抑制火焰的发展。
Fig. 7 Sensitivity coefficient at different equivalence ratios图7 不同当量比下的敏感性系数 |
图8为氢气-乙醇-空气预混火焰在Φ = 1.2时,分别在6种不同H2比例、3种不同初始压力下的前6个敏感反应。由图可知,R1、HCO + M = H + CO + M (R26)、CO + OH = CO2 + H (R24)和CH3 + O = CH2O + H (R108)为起到正向促进作用的基元反应。其中R1为最主要的反应,该反应消耗了H自由基并产生了O和OH自由基,是消耗H2的主要反应链,有助于提升LBV。在当量比和初始温度相同的情况下,R8和R96这两种主要起到抑制作用的基元反应的敏感性系数随着初始压力的上升而上升,而这两种主要的基元反应消耗了高活性基,导致LBV随着初始压力的上升而下降。
Fig. 8 Sensitivity coefficient at different initial pressure图8 不同初始压力下的敏感性系数 |
3 结论
基于Chemkin中的Premix模块搭建一维火焰模型,从层流燃烧速度、密度比、绝热火焰温度、火焰厚度、敏感性五个方面,探讨了当量比、初始压力和H2比例对氢气-乙醇-空气预混火焰的燃烧特性的影响。结论如下：
（1）氢气-乙醇-空气预混火焰的LBV峰值出现在当量比为1.2左右,并且随着氢气含量的上升,LBV峰值对应的当量比会朝着富燃一侧移动。
（2）随着当量比的上升,密度比和绝热火焰温度都呈现出先上升再下降的规律,而火焰厚度则呈现相反的变化,其峰值都出现在当量比为1.2或者1.3。这表明预混燃气的流体动力学不稳定性随着当量比的变化呈现出先减小后增加的趋势,并且在1.2或1.3时达到最小。
（3）链分支反应H + O2 = O + OH (R1) 的敏感系数在任何工况下都是最高的,该反应可以生成O和OH这两种高活性的自由基,使得链式反应快速进行,释放大量的热从而提高火焰速度,因此该链分支反应的变化规律在较大程度上会影响LBV的变化。而基元反应R8、R13、R9和R96的敏感性系数为负数,对火焰的发展起到抑制作用。
