选用Fluent进行数值模拟计算, 燃烧模型选用有限速率/涡耗散模型, 化学反应及燃烧模型中守恒方程为

∂∂t(ρYi)+∇(ρv→Yi)=-∇J→i+Ri+Si(1)

式(1)中, ρ 为密度, Y_i是组分i的质量分数, v为速度, J_i为物质i的扩散通量, R_i是化学反应中组分i的净产生速率, S_i为源项。

J_i由浓度梯度产生, 可由式(2)计算得到

Ji=-ρDi, m∇Yi(2)

式(2)中, D_{i, m}是混合物中物质i的扩散系数。 多个化学反应的组分i的净生成速率由式(3)得到

Ri=Mw, i∑r=1NRR^i, r(3)

式(3)中, M_{w, i}是组分i的分子量, R^i, r是组分i在反应r中产生/分解Arrhenius摩尔速度。

以甲烷/空气同轴射流燃烧为研究对象, 为准确模拟甲烷的燃烧过程, 考察了包含激发态粒子反应在内的多组分甲烷燃烧机理。 基态粒子化学反应模型采用GRI-Mech3.0机理, 共包含53个组分, 325个基元反应, 完整地描述了基态粒子的产生与消耗。

探究OH^* 和CH^* 在燃烧中的分布特性以及与反应热的关系, 因此燃烧机理中还需考虑OH^* 和CH^* 的生成与淬灭, 详细机理如表1所示。 综合基态与激发态的组分及反应, 最后得到37个组分, 237个反应步的详细化学机理。

表1

Table 1

表1(Table 1)

表1 OH* /CH* 自由基化学反应动力学参数 Table 1 Kinetics parameters of OH* and CH* radicals reaction No. Reaction Ref. R1 CH+O2→ CO+OH* [11] R2 H+O+M→ OH* +M [12] R3 OH* → OH [13] R4 OH* +N2→ OH+N2 [13] R5 OH* +O2→ OH+O2 [13] R6 OH* +H2O→ OH+H2O [13] R7 OH* +H2→ OH+H2 [13] R8 OH* +CO2→ OH+CO2 [13] R9 OH* +CO→ OH+CO [13] R10 OH* +CH4→ OH+CH4 [13] R11 C2H+O2→ CO2+CH* [13] R12 C2H+O→ CO+CH* [13] R13 CH* → CH [13] R14 CH* +N2→ CH+N2 [13] R15 CH* +O2→ CH+O2 [13] R16 CH* +H2O→ CH+H2O [13] R17 CH* +H2→ CH+H2 [13] R18 CH* +CO2→ CH+CO2 [13] R19 CH* +CO→ CH+CO [13] R20 CH* +CH4→ CH+CH4 [13]

表1 OH* /CH* 自由基化学反应动力学参数 Table 1 Kinetics parameters of OH* and CH* radicals reaction

以实验室采用的同轴燃烧器作为数值计算的对象, 如图1所示, 火焰为对称性结构, 选取1/2截面进行数模计算, 计算区域沿轴向长度为100 mm, 沿径向总长为25 mm, 中心燃料出口为1 mm, 空气出口2 mm, 燃料和空气之间的壁厚2 mm, 且向上延伸5 mm, 具体二维计算区域如图2所示, 对表2和表3的工况进行数值模拟。

图1

Fig.1

	Figure Option ViewDownloadNew Window
	图1 化学发光测量系统Fig.1 Chemiluminesenes measurement system

图2

Fig.2

	Figure Option ViewDownloadNew Window
	图2 计算域示意图Fig.2 Calculate domain

表2

Table 2

表2(Table 2)

表2 扩散火焰工况 Table 2 Diffusion flame conditions 工况 vCH4/(m·s-1) vair/(m·s-1) 1 1 1 2 1 0.4 3 1 0.12

表2 扩散火焰工况 Table 2 Diffusion flame conditions

表3

Table 3

表3(Table 3)

表3 预混火焰工况 Table 3 Premix flame conditions 工况 vpre/ (m·s-1) vair/ (m·s-1) 局部当量 比φ 1 1 0.2 1.0 2 1 0.2 1.2 3 1 0.2 1.6

表3 预混火焰工况 Table 3 Premix flame conditions

燃料气体的进口边界条件为速度入口, 空气为压力入口, 底部右侧为壁面, 计算域的右边设置为壁面, 实验中采用玻璃罩, 左边界设置为对称轴, 上边界设置为压力出口边界条件, 操作压力为101 325 Pa。 为保证模拟的准确性, 所有网格均采用结构化网格, 网格质量均在0.95以上, 网格总数约为62万。

同轴射流燃烧器是用于产生射流扩散火焰的一种典型燃烧器, 结构简单, 主要分为中心射流喷嘴和同轴射流喷嘴两部分, 同样具有内外两个腔体, 但由于内外分别通入燃料和氧化剂, 不涉及混合, 此处两个腔体只作为集气腔, 保证喷嘴上游压力稳定。 贴壁面安装隔网, 在其中填充若干钢珠, 防止回火, 同时起到一定的整流作用, 保证出口处的气流速度均匀分布。 燃料由中心集气腔进入中心射流喷嘴, 再从射流喷嘴出口流出, 与同轴射流空气进行扩散燃烧。

由于化学发光十分微弱, 因此必须使用ICCD相机采集发光信号; 成像系统主要由ICCD相机、 镜头、 近摄镜、 滤波片构成, 相机可探测的波段为190~1 080 nm, 包含研究所需的300~500 nm波段, 并且曝光时间、 转换速率等参数均较好的满足试验需求。 为拍摄不同激发态粒子化学发光的空间分布, 在镜头前加装对应波段的滤波片, 过滤掉其他波段信号的干扰。 因主要研究OH^* 和CH^* 两种激发态粒子, 对应的中心波长分别为308和431 nm。 详细试验参考胡悦的甲烷-空气预混火焰中化学发光定量化测量放热率的试验研究^[14]。

图3为扩散火焰工况2下OH^* 和CH^* 的二维分布实验结果和模拟结果的对比, (a)和(b)分别是OH^* 摩尔质量分布的实验结果和模拟结果, (c)和(d)分别是CH^* 摩尔质量分布的实验结果和模拟结果。 从图中可以看出模拟结果与实验结果有相同的变化趋势以及很好的吻合性。 OH^* 分布主要分为三个燃烧区域, 离开喷嘴处摩尔质量最大的第一燃烧区, 然后向火焰下游延伸的第二燃烧区, 此区域范围较大, 第三燃烧区摩尔质量逐渐减小, 此区域范围较小。 CH^* 分布主要为两个燃烧区域, 离开喷嘴的第一燃烧区域摩尔质量最大, 然后向火焰下游扩散形成第二燃烧区。 模拟结果得到摩尔分数稍大, 这是由于燃烧机理相对于实际燃烧有简化, OH^* 和CH^* 的生成和淬灭略有偏差, 但该偏差对模拟结果影响较小, 可以认为本数值模拟能够准确地反映激发态粒子OH^* 和CH^* 的分布以及燃烧情况。

图3

Fig.3

	Figure Option ViewDownloadNew Window
	图3 OH* 和CH* 实验结果与模拟结果 (a): OH* 实验结果; (b): OH* 模拟结果; (c): CH* 实验结果; (d): CH* 模拟结果Fig.3 Distributions of OH* /CH* by measurement and computation (a): Distribution of OH* by measurment; (b): Distribution of OH* by computation; (c): Distribution of CH* by measurment; (d): Distribution of CH* by computation

由图4(a, b, c)中可以看出三种工况下中心CH₄速度不变, 逐渐增大空气流速时, 反应热和OH^* 和CH^* 的分布趋势一致, OH^* 和CH^* 质量分数增大时, 反应热也增大, 当激发态粒子OH^* 和CH^* 的质量分数达到峰值的时候, 反应热也达到峰值。 空气速度分别是0.12, 0.4和1.0 m·s^-1时, 反应热的第一峰值和激发态粒子的峰值分别位于0.003 1, 0.002 75和0.002 4 mm, 燃烧区域逐渐向中心轴移动, 由此可以看出OH^* 和CH^* 能标识反应热。

与激发态粒子分布特性不一样的是反应热呈现双峰现象。 反应热第一峰值的变化与OH^* 和CH^* 相似, 随着燃烧反应的进行, 反应热迅速增加达到第一波峰, 随后迅速减小, 但减小到离开激发态粒子的反应区域后与OH^* 和CH^* 的分布不相同, OH^* 和CH^* 是迅速减小到零, 但反应热产生拖尾现象, 释放热量逐渐增加到第二峰值后再缓慢减小。

反应热达到第一峰值后迅速减小, 然后出现第二个波峰后再逐渐减小为零, 为了解此现象将对化学反应速率进行分析。

图5(a, b, c)为不同空气速度下反应热、 OH^* 、 CH^* 生成速率的轴向分布图, 从图中可知CH₄速度不变, 逐渐增大空气速度, 不同工况下反应热和激发态粒子OH^* 和CH^* 的生成速率分布趋势一样, 反应热存在第一峰值和第二峰值呈现双峰现象。

图4

Fig.4

	Figure Option ViewDownloadNew Window
	图4 反应热与OH* /CH* 质量分数径向分布图 (a): 空气速度0.12 m·s-1; (b): 空气速度0.4 m·s-1; (c): 空气速度1.0 m·s-1Fig.4 Reaction heat and OH* /CH* mass fraction profiles along radius (a): Air velocity 0.12 m·s-1; (b): Air velocity 0.4 m·s-1; (c): Air velocity 1.0 m·s-1

图5

Fig.5

	Figure Option ViewDownloadNew Window
	图5 反应热与OH* /CH* 反应速率轴向分布图 (a): 空气速度0.12 m·s-1; (b): 空气速度0.4 m·s-1; (c): 空气速度1.0 m·s-1Fig.5 Reaction heat and OH* /CH* chemical reactionrateat profiles along axial (a): Air velocity 0.12 m·s-1; (b): Air velocity 0.4 m·s-1; (c): Air velocity 1.0 m·s-1

CH^* 的主要生成反应机理是C₂H+O₂=CH^* +CO₂(R11)和C₂H+O=CH^* +CO(R12)。 离开喷嘴出口CH₄迅速反应生成C₂H, 此时甲烷和空气接触区域内氧气含量高, 与C₂H反应生成CH^* 。 随着氧气参与反应后逐渐转变成游离氧, 在反应区域内游离氧的含量增加, 游离氧是R12的反应物, 因此反应速率增加, 此刻两者反应速率都达到最大, 释放的热最多, 反应热达到第一峰值。

OH^* 的主要生成反应机理是CH+O₂=OH^* +CO(R1)和H+O+M=OH^* +M(R2), 离开喷嘴后燃烧区域内氧气含量高, 但是CH含量很少, 因此R1反应速率很小。 随着火焰不断向下游方向传播, 游离氢和游离氧的质量分数不断增加, 两者迅速结合反应生成OH^* , 使R2的反应速率增加, 从而反应热增加, 达到第二波峰, 因此反应热呈现双峰现象。

从图6(a, b, c)中可以看出, 在三种工况下中心预混气体和空气速度不变, 逐渐增大中心预混气体的当量比, 反应热和激发态粒子OH^* 和CH^* 的分布趋势一样, 在激发态粒子OH^* 和CH^* 质量分数增大时, 反应热也增大, 当激发态粒子OH^* 和CH^* 的质量分数达到峰值的时候, 反应热也达到峰值, 由此可以表明激发态粒子OH^* 和CH^* 能标识反应热。

图6

Fig.6

	Figure Option ViewDownloadNew Window
	图6 反应热与OH* /CH* 质量分数径向分布图 (a): 局部当量比ε =1.0; (b): 局部当量比ε =1.2; (c): 局部当量比ε =1.4Fig.6 Reaction heat and OH* /CH* mass fraction profiles along radius (a): Part equivalenc ε =1.0; (b): Part equivalenc ε =1.2; (c): Part equivalenc ε =1.4

图7(a, b, c)为不同空气流速下反应热、 OH^* /CH^* 生成速率的轴向分布图, 从图中可见中心预混气体和空气速度不变, 逐渐增大中心预混气体的当量比, 反应热和激发态粒子OH^* 和CH^* 的生成速率分布一样。 当激发态粒子的反应速率达到最大值时, 反应热也达到最大值。

图7

Fig.7

	Figure Option ViewDownloadNew Window
	图7 反应热与OH* /CH* 反应速率径向分布图 (a): 空气速度0.12 m·s-1; (b): 空气速度0.4 m·s-1; (c): 空气速度1.0 m·s-1Fig.7 Reaction heat and OH* /CH* chemical reactionrateat profiles along radius (a): Air velocity 0.12 m·s-1; (b): Air velocity 0.4 m·s-1; (c): Air velocity 1.0 m·s-1

本研究是模拟预混火焰燃烧, 中心预混气体的当量比大于1, 因此CH^* 的主要生成反应机理是C₂H+O₂=CH^* +CO₂(R11), OH^* 的主要生成反应机理是 CH+O₂=OH^* +CO(R1)。 H+O+M=OH^* +M(R2)和C₂H+O=CH^* +CO(R12)的反应物为游离氧, 反应速率较R1和R11的反应速率小。 激发态生成反应速率最大的位置也是在反应热最大的位置, 由此可见在激发态粒子生成速率最大时, 释放的热量最多, OH^* 和CH^* 能标识反应热。