管状火焰燃烧器中甲烷富氧燃烧分析.docVIP

基金项目：国家自然科学基金资助项目（No 中国工程热物理学会 燃烧学 学术会议论文 编号：184002 管状火焰燃烧器中甲烷富氧燃烧分析 马康，初庆钊，石保禄*，唐国荣，王宁飞 （北京理工大学 宇航学院，北京 100081） (Tel: Email: shibaolu@bit.edu.cn) 摘 要: 采用OpenFOAM对急速混合管状火焰燃烧器中以CO2为稀释剂的甲烷富氧燃烧进行大涡模拟，分析了流速、CH、OH、温度及CO分布特征，并与冷流流场对比，探讨了急速混合管状火焰燃烧器中层流火焰的流场特征，构建了火焰结构。结果表明，冷流条件下的湍流流动在强旋流燃烧器中转化为稳定的层流火焰，温度沿径向均匀分布，且切向入口附近壁面温度较低，有助于保护燃烧器；通过调节流量增加入口切向速度，燃料与氧化剂掺混加强，可大幅降低CO排放量。 关键词: OpenFOAM；管状火焰；富氧燃烧；大涡模拟 0前言 随着经济社会的发展，人类对能源的需求日益增长[1]。传统化石燃料燃烧方式带来全球变暖和空气质量污染等问题。富氧火焰可以提高燃烧效率并减少NOx排放，被广泛应用于工业生产和生活中[2]。近年来，Satoru Ishizuka[3]等人的将燃料和氧化剂分别切向通入燃烧室，燃料与氧化剂在燃烧室内快速混合和燃烧，形成稳定的管状层流火焰。这种急速混合管状旋流燃烧器具有火焰稳定性强，无回火，有利于燃料和氧化剂的掺混等特点，被用于开展甲烷富氧燃烧。 现今，常见的甲烷富氧燃烧的氧化剂是空气和氧气，由于空气中含有约71%的氮气，大量的氮气并不参加反应，带走了大量的热，同时可能会在反应过程中生成有害的氮氧化合物，燃烧效率不高且对环境有害[4]。近年来，国内外研究者们提出了以CO2为稀释剂的富氧燃烧（Oxy-fuel Combustion）技术，一方面有利于减少热力型NOx的排放，另一方面，燃烧产物主要为CO2和H2O，有利于提高CO2的回收和循环利用效率，减少污染，节约成本[5]。由于实验仪器的限制，对于急速混合管状火焰中组分、污染物、速度场的研究分析还很欠缺。为了更好地利用CCS[6]（CO2 Capture and Sequestration）技术，进一步地研究富氧燃烧的燃烧特性，有必须要开展管状旋流燃烧器中甲烷富氧燃烧数值计算。目前，湍流流动的数值模拟方法有三类，即直接模拟方法，大涡模拟方法和Reynolds时均方程的模拟方法。直接模拟方法对每个网格求解瞬态Navier-Stokes方程，其计算量与雷诺数的三次方成正比，对计算资源需求巨大，目前难以实现对高雷诺数复杂湍流的数值模拟。其中Reynolds时均方程模拟方法在工程流动计算中应用范围较广，虽然能在较低的计算资源要求下实现高雷诺数复杂流动的数值模拟，但是精度较差。上世纪70年代提出的大涡模拟[7]通过滤波方程过滤出较大空间尺度和时间尺度的涡，对大尺度的涡直接求解，小尺度的涡采用模拟的方法。因此大涡模拟的精度和计算量介于直接模拟和雷诺平均模拟之间，流场解析度较高，能更真实地模拟流场特征，在当前的计算资源下有更大的优势。LES方法在十多年前就已经是CFD领域研究和应用的热点[8-9]。LES是通过建立数学滤波函数，从湍流瞬时运动方程中将尺度比滤波函数尺度小的涡过滤掉，从而分解出描述大涡流场的运动方程，而被过滤掉的小涡对大涡运动的影响，则是通过建立亚格子尺度模型（SGS）[10-14]来实现。 本文采用了入口宽度为1mm的管状旋流燃烧器，结合开源软件OpenFOAM，借助大涡模拟，对甲烷富氧燃烧进行了数值计算。在相同氧气浓度为0.5，当量比为1条件下，首先对Q=3m3/h情况下冷态流场和燃烧流场，组分和温度变化情况等进行了具体分析。随后对不同入口总流量(1，3，6m3/h)的工况分析，研究CO2稀释燃烧的火焰的燃料掺混，出口附近污染物排放。 1计算方法 1.1 计算仿真模型 实验所用的管状火焰旋流燃烧器的结构如图1所示，右侧为燃烧器整体结构，左侧为入口喉部截面。燃烧器主要由灰色部分的不锈钢结构，四个均匀分布的切向入口，石英观察窗和内径16mm，长L=75mm的石英玻璃管构成，石英观察窗和玻璃管用于观察和记录火焰结构，为了与实验参数形成对比验证算法，三维计算模型几何参数与该实验燃烧器完全一致，三维网格划分如图2所示，网格总数约250万。 图 SEQ 图_ \* ARABIC 1 管状火焰旋流燃烧器 Figure