自激振荡燃烧的模拟和数值计算.pdf

中国工程热物理学会 燃烧学 学术会议论文 编号：084265 自激振荡燃烧的模拟和数值计算 姚兆普，朱 民 清华大学热能工程系，热科学与动力工程教育部重点实验室，100084 (Tel:010 Email: zhumin@tsinghua.edu.cn) 摘要：燃烧室中的燃烧稳定性直接关系到燃气轮机的噪声，效率及污染物的产生。振荡燃烧是由压力 波动和非稳定燃烧之间的相互耦合导致系统共振引起的。本文利用热声学的相关理论，构建用于自激 振荡燃烧数值计算的边界模型，并将计算结果与实验结果相比较，两者吻合良好。通过对燃烧区典型 位置的分析，获得自激振荡燃烧周期内流场中相关参数的变化规律和耦合关系，这些结果对于研究振 荡燃烧的机理，发展抑制振荡燃烧的方法具有重要意义。 关键词：自激振荡燃烧，热声耦合，数值计算。 0 前言 燃烧室作为燃气轮机的核心部件之一，其发展直接影响到燃气轮机甚至整个动力系 统的性能。液体或气体燃料在燃烧室中的燃烧状况直接关系到燃气轮机的噪声，效率及 污染物的产生。振荡燃烧是燃烧室的工作过程中经常出现的一种非稳定燃烧状况。造成 振荡燃烧的原因是多种多样的，例如燃烧室内压力的脉动，燃烧区放热的脉动，燃料流 量随工作状态的波动以及进气流场的不稳定性。从理论上讲，振荡燃烧是由压力波动和 非稳定燃烧之间的相互耦合导致系统共振引起的。振荡燃烧不但会造成燃烧效率下降， 污染物生成率的提高，严重地时候甚至会对燃烧室造成结构损坏。因此，研究燃烧室振 荡燃烧的有关现象和本质，找到振荡燃烧的有关规律，进而发展抑制振荡燃烧的方法， 在工业的实际应用中具有重要意义。 对于燃烧振荡的最初认识来源于 Rijke 在 1859 年完成的一个实验。在实验中，Rijke 观察到了一个两端开口的圆管中，被火焰激发的空气振动现象。瑞利从此类现象出发， 发展了压力脉动和火焰放热率脉动之间的耦合准则。可以写作如下形式： 1 G p =Ω t q Ω t dtdΩ （1） ∫ ∫ ( , ) ( , ) 0 Ω T T 其中p 和q 分别为压力脉动和放热率脉动。T 为描述的时间范围，Ω为描述的流 场范围。方程(1)表示当压力脉动和放热率脉动为同相时，燃烧系统将可能出现振荡；反 之，燃烧系统将是稳定的。 国内外的有关学者已经作了很多振荡燃烧的实验研究和分析计算相关工作。例如 Bloxsidge[1]构建了一个关于不稳定燃烧的火焰模型并将其和在一个气体燃烧室中得到 的实验数据相比较。Schuller[2]和 Noiray[3]对于环形燃烧室预混火焰的稳定性进行了分 资助项目：国家重点基础研究发展计划（973 计划）(No. 2007CB210106) 析。他们提出了一个用于解释此种条件下不稳定燃烧的机理，认为不稳定燃烧主要是由 于火焰对于流场振荡的响应与燃烧室的声学特性相互耦合造成的。Uhm 和 Acharya[4]对 于高速射流中的液滴的不稳定燃烧问题进行了实验研究。他们在高速射流的基础上添加 人为扰动，使整个流场产生周期性振荡，并对此进行控制，抑制不稳定燃烧现象。Zhu[5] 研究了具有实际工业背景的燃气轮机燃烧室的自激振荡燃烧情况。 本文将对特定实验下振荡燃烧的形成机理和相关特性进行研究[6]，主要采用理论分 析和数值计算的研究方法，并和实验的结果进行比较。通过热声学有关原理推导出适合 于数值计算的自激振荡边界条件的相关模型，并将数值计算结果和实验结果进行对比。 1 边界模型 实验装置如图 1 所示，其中位置 1 为堵塞板，位置 1 到位置 2 为进气段，位置 3 为 燃烧段入口，其出口为位置 5，通向大气环境。通过三维非定常 CFD 模拟该系统的振荡 现象