预混火焰不稳定性的大涡模拟.pdfVIP

中国工程热物理学 燃烧学 学术会议论文 编号：094019 预混火焰不稳定性的大涡模拟 刘联胜 刘晶 郭平平 谷岩 段润泽 河北工业大学能源与环境工程学院，天津，300132 联系电话Email: liuliansheng@hebut.edu.cn 摘要：本文在湍流燃烧大涡模拟的基础上，利用气动声学FW −H 方程对甲烷/空气预混火焰的不稳定 性进行了数值分析。局部化学当量比的周期性波动，造成燃烧室内出现周期性的压力振荡，随着化学 当量比的增大或减小，压力振荡幅值均有所提高；贫氧工况下，压力振荡的频率较低，贫燃工况下， 压力振荡频率较高。燃烧室内存在周期性的温度振荡和速度脉动；压力振荡与速度脉动趋于同频、同 相，而温度振荡在相位上稍有提前。流向涡涡量分布呈现周期性变化趋势，相对而言，贫氧工况下， 流向涡涡量较小、涡团尺度较大；贫燃工况下，流向涡涡量较大、涡团尺度较小。 关键词：预混火焰，热声不稳定性，化学当量比，压力振荡，涡量 0 前言 贫燃预混、部分预混和空气分级燃烧技术是目前燃气轮机等动力机械普遍采用的低 NO x 燃烧技术。由于上述燃烧过程在偏离化学恰当比情况下进行，而过低或过高的化学 当量比都将恶化燃料与空气的混合效果，加剧燃烧器下游速度分布和浓度分布的不均匀 性，从而使火焰结构发生显著变化，火焰出现非线性振荡现象。在某些条件下，火焰的 非线性振荡将诱发火焰热释放与燃烧室声振动之间发生激振，形成热声耦合振荡[1] 。文 献[2-7][15] 的实验研究表明，火焰热声耦合振荡是声压振动、热释放波动以及流场脉动相 互作用的结果，其诱发机理十分复杂，局部化学当量比、预混流场结构、火焰表面涡团 尺度及涡团脱落频率、燃烧室结构等参数的变化，均可能诱发低频火焰振荡（振荡频率 f ≤500Hz ），造成火焰的动力学失稳，使燃烧状况恶化。 随着燃烧理论和计算机技术的发展，数值模拟已经成为研究燃烧过程的一个重要手 段。文献[8]-[13]利用不同的湍流燃烧模型模拟了预混燃烧过程，发现模拟所得到流场结构、 污染物生成以及化学动力学特性等与实验诊断结果在趋势上基本吻合；通过将湍流燃烧 模型与声学模块的耦合求解，可以得到燃烧室内的声压振动情况。但是，由于计算模型 对壁面条件的简化，其结果尚无法准确反映壁面反射、折射作用对压力振荡的影响。 本文利用FLUENT6.2 对不同化学当量比φ 的甲烷/空气非稳态预混燃烧过程进行了 大涡模拟，并通过耦合求解气动声学FW −H 方程，分析了化学当量比φ 对预混火焰不稳 定性的影响规律。本文结果为预混燃烧热声耦合振荡控制策略的制定提供了基础数据和 理论依据。 1 物理模型 [14] 本文燃烧室的物理模型根据文献 的燃烧装置简化得到，实验装置及物理模型如图 1 所示。燃烧室内径为90mm，长度为200mm ，甲烷/空气预混气入口直径10mm， 国家自然科学基金资助项目）；河北省自然科学基金资助项目（E2009000079 ） 环形直流空气入口内径30mm、外径50mm，燃烧室出口直径为40mm 。本文采用 Gambit2.1 构建三维轴对称计算区域，如图2 所示，网格采用非结构化Pave 网格，网格 数10046 个，最小控制体积为4.67 ×10−10 m3 ，最大控制体积为6.4 ×10−8 m3 。 图1 实验装置及物理模型 图2 网格划分 边界条件如下：（1）入口条件：假设预混气体为理想气体，入口边界条件设为速度 入口。预混可燃气由中心入口以旋流形式进入燃烧室，其轴向速度和切向速度均为 5 m s ，径向速度为零；环形入口速度根据空气体积流量确定，直流送入。（2）壁面边 界条件：FLUENT6.2 默认的壁面材料为铝，并假设壁面处无速度滑移，且湍流脉动量 为零；壁面温度恒定为300K 。 2 数学模型 因此本文利用大涡模拟方法（LES