壁面参数对发动机微燃烧室氢气空气预混燃烧的影响.docVIP

壁面参数对发动机微燃烧室氢气空气预混燃烧的影响 　　摘要：根据微尺度燃烧室内氢气和空气的预混燃烧过程，利用STAR-CD计算软件建立微燃烧的物理和数学模型，在试验验证的基础上模拟壁面导热系数、壁厚和外壁面传热系数等壁面参数对氢气和空气预混合燃烧的影响。结果表明：靠近燃烧室入口处的火焰中心温度随导热系数减小而增加，出口处混合气体的温度随导热系数减小而降低；随着外壁面厚度增加，外壁面温度降低，壁面温度分布更均匀；随着外壁面传热系数增加，燃烧室内火焰温度降低，出口处的气体温度和外壁面温度相应降低。 　　关键词：微燃烧室；壁面参数；数值模拟；燃烧特性；氢气 　　中图分类号：TK42 文献标识码：A 文章编号：1674-1161（2017）09-0044-04 　　微燃烧器是微热光电系统的核心部件，其结构和性能对外壁面温度分布的均匀性有重要影响，进而影响光电能量的转换效率。近年来，国内外学者对微燃烧器的结构和性能进行了大量研究。受加工尺寸的限制，研究壁面参数对微燃烧的影响比较困难。为此，利用STAR-CD计算软件建立微燃烧的物理和数学模型，在试验验证的基础上模拟研究燃烧室壁面导热系数、壁厚和外壁面传热系数等壁面参数对微尺度燃烧室内氢气和空气预混合燃烧的影响，旨在为微燃烧室结构设计和性能改善提供一定的参考。 　　1 数值模型和方法 　　1.1 物理数学模型 　　研究的模拟对象（如图1所示），为内径0.2 mm、壁厚0.05 mm、长1.6 mm的微型圆柱管。采用STAR-CD前处理软件pro-star对微燃烧室进行网格划分。燃烧室结构较为规则，因此采用结构化网格形式，在简单计算区域内可以直接生成。利用网格特点可以构造各种算法。 　　1.2 数学模型 　　数学模型包括以下基本方程：质量守恒方程、动量守恒方程、化学组分守恒方程、流体能量守恒方程、壁面能量守恒方程。 　　1.3 流体流动模型 　　该燃烧室尺度较小，但研究的传热及流动问题仍属于连续介质区，适用Navier-Stokes方程及Fourier热传导定律，可运用商业通用软件STAR-CD进行数值模拟研究。在微圆管内通入氢气/空气预混气体燃烧，以空气物性参数作为整个混合气流的物性参数，认为是层流粘性流，可用理想气体定律计算流体密度。用物质性质的质量分数的加权平均数来计算流体混合物的比热、粘度、热传导率。 　　1.4 燃烧模型 　　燃烧反应模型很多。考虑到氢气和氧气的化学反应，选择适于层流预混合燃烧的Chemical Kinetics模型，其化学反应率遵循Arrhenius定律。 　　1.5 传热模型 　　在微尺度燃烧室内燃烧过程中，由于表面积体积比增大，因此壁面热损失不能忽略，必须综合考虑导热、对流和辐射散热产生的热损失。 　　1.5.1 导热 壁面导热损失采用Conduct and conjugate Heat Transfer模型计算。模拟计算时采用双精度，这样利于解决潜在的收敛问题。 　　1.5.2 对流换热 在计算中需要确定对流换热产生热损失的传热系数，计算公式为： 　　1.5.3 辐射换热 主要考虑燃烧室外壁面的辐射散热对微尺度燃烧室内燃烧的影响。在STAR-CD模拟中，如果计算区域中没有建立空气网格，燃烧室外壁面辐射散热可不设置辐射模型，但需要在固体外壁面边界条件中给出辐射散热的热阻。辐射传热系数hra的计算公式为： 　　1.6 边界条件和算法 　　在边界条件中，入口条件选取速度边界条件，需要给出预混合气在入口的流速、温度和气体混合物质量百分比。燃烧室的结构属于轴对称，为全面分析微燃烧室的火焰传播特性，采用三维数值模拟。数值模拟采用1/2圆柱作为计算区域，计算区域的底部边界为轴对称边界条件。出口压力边界条件固定为 　　1.013 25×105 Pa。燃烧室壁面采用无滑移、无物质流出表面壁面边界条件。 　　燃烧室壁面的热边界条件与2个主要传热因素有关：壁面散失到外界环境的热损失和墙壁内的热传导。壁面包括内壁面边界和外壁面边界。内壁面边界热条件选择为传导，不需要给出具体的导热热阻，STAR-CD会根据材料特性和壁面厚度自行计算。外壁面边界条件只需给出综合的传热热阻就可以确定散热量。 　　采用只用于稳态计算的SIMPLE算法求解上述控制方程。在计算过程中，收敛条件规定各个控制方程迭代计算的余量至少小于10-3。速度的差分格式采用重构格式（MARS，Monotone Advection and Reconstruction Scheme），?@样不仅可提高梯度剧烈变化时的计算精度，还能减少对网格的依赖性。 　　2 模型验证 　　所模拟的研究对象尺寸较小，已经超出了目前国内能加工的最小尺寸。为验证计算模型的准确性，与黄俊等开展的试验进