宽高比对矩形微细通道中CH4Air火焰稳定性影响的数值分析.docVIP

基金项目：国家自然科学基金资助项目 中国工程热物理学会 燃烧学 学术会议论文 编号：184260 宽高比对矩形微细通道中CH4/Air火焰稳定性影响的数值分析 袁自力，范爱武* （华中科技大学能源与动力工程学院，武汉 430074） (Tel: Email: faw@hust.edu.cn) 摘 要 对CH4/air预混气在高度为4mm、宽高比为1~4的矩形通道中的燃烧过程进行了三维数值模拟。结果显示：对于方形截面通道，火焰在吹熄前不发生倾斜；宽高比为2和3时，随着进气速度的增大，火焰发生倾斜；而宽高比增大到4时，火焰在吹熄前出现了振荡。火焰吹熄极限随着宽高比的增大，先增大、后减小，在宽高比为3时获得最大值0.66 m/s。分析表明，稳燃极限在宽高比较小的通道中主要受散热损失的影响，而在宽高比较大的情况下则主要受水动力学不稳定性的影响。本工作可为矩形截面微燃烧器的优化设计提供理论指导。 关键词 微尺度燃烧；截面宽高比；火焰吹熄极限；散热损失；水动力学不稳定性 0 前言 微机电技术由微传感器、微电路、微动力系统等组成[1-2]。碳氢化合物燃料的能量密度远大于传统化学电池，且来源充足[3]，是微机电技术的理想能源。因此，微尺度燃烧受到了广泛关注。然而，微小尺度下的燃烧稳定性受到几个方面的影响。首先，微通道特征尺度缩小导致散热损失比例急剧增大[4]。其次，气体混合物在微通道中的驻留时间非常短暂[5]，燃料中的化学能来不及充分释放，从而降低火焰温度和燃烧效率。另外，燃烧室的比表面积大，也增加了壁面对自由基的捕获和销毁几率，产生壁面淬熄现象[6]，降低了火焰根部的稳定性。以上几个因素综合起来，使得微小尺度下的稳燃极限显著减小。目前常用的微尺度稳燃方法有：设计热循环结构，如瑞士卷燃烧器[7]；构建回流区，如凹腔和钝体稳燃器[8-10]；催化剂助燃[11]。 对于直通道内的燃烧来说，影响火焰稳定性的具体参数包括：壁面厚度及材料[12]，进气速度和当量比[13]，管径或微通道间距[14]以及通道截面形状[15, 16]等等。Norton和Vlachos [12, 13]通过数值模拟发现壁面厚度和热物性严重影响微通道的回热效果（热循环效应），并指出：进气速度过小会导致产热过少，散热损失比例增大；进气速度过大则减小了反应物驻留时间。Wan等[14]通过数值模拟，揭示了火焰分裂极限随通道间距增加呈现非单调变化规律。Pizza[15]使用三维模型研究了甲烷在不同直径的微圆管中的火焰动态。Di[16]使用三维模型计算了丙烷在方形和圆形截面微通道内的可燃速度范围。 为了减小计算工作量，许多学者[12-14]往往忽略微通道宽度方向（或周向）对火焰稳定性的影响，使用二维模型进行模拟研究。最近，Pashchenko[17]通过对比研究指出，只有当微圆管的长度相对于管径足够大时，二维模型才能得到与三维模型相近的结果。文献调研表明，目前尚无关于矩形微细通道的宽高比对火焰稳定性影响的研究报道。因此，本文使用三维模型对CH4/air 预混气体在高度为4mm的狭窄矩形通道中的燃烧过程进行数值模拟，考察截面宽高比（1、2、3、4）对火焰稳定性的影响，并从散热损失和流场特性两个方面进行了分析。 1 计算方法 矩形微细通道的横截面如图1所示。图中O点为通道进口的中心位置，取为坐标原点。通道宽度、高度、长度方向分别设为X、Y、Z方向。通道高度H、长度L和壁厚固定为4 mm、40 mm和1 mm。宽高比（W/H）分别取1、2、3、4四个值。 图1 矩形微细通道的横截面示意图 模拟使用的燃料为甲烷，氧化剂为空气，预混气当量比固定为1。通道壁面材料采用耐热石英玻璃，其密度、比热、热导率以及壁面发射率ε分别为2650 kg/m3、750 J/kg·K、1.05 W/m·K和0.92。进口最大雷诺数（Re）约为400，故采用层流模型。CH4/Air的燃烧反应机理采用Bilger等[18]提出的包含18种组分、58个基元反应的C1机理。基元反应速率采用层流有限速率模型进行计算。气体混合物假设符合理想气体状态方程，其热力学参数与输运参数来自CHEMKIN数据库。 边界条件设置如下：（1）采用速度入口条件，预混气初始温度为300 K，压力为1个大气压。（2）采用压力出口边界条件。（3）燃烧器的外壁面与环境之间的换热同时考虑辐射和自然对流两种方式，单位面积的散热速率为：