1、文章编号：1009 444X（2022）04 0347 06气体燃烧室中热流场与生成 NOx浓度场的数值模拟研究张诺楠，郭韵（上海工程技术大学 机械与汽车工程学院,上海 201620）摘要：对于一种燃气燃烧器，选取合适的燃烧室尺寸进行匹配，进而构成一个完整的计算区域.天然气燃料和助燃空气通过该燃气燃烧器喷入燃烧室实现完整燃烧过程.通过简化燃气燃烧器完成几何建模，使用热流与燃烧分析软件对计算区域进行数值模拟，选用 Realizable k-湍流模型、ED 燃烧模型、DO 辐射模型及热力型 NOx生成模型，探究过量空气系数 对该燃气燃烧器燃烧过程的影响，得出其温度场及生成 NOx浓度场的相关云图及
2、分布情况并进行分析.当过量空气系数 =1.1 时，NOx浓度最大.研究结果可为相关研究实验和设计低 NOx燃烧器提供参考依据.关键词：燃气燃烧器；氮氧化物；数值模拟；过量空气系数中图分类号：TK17 文献标志码：ANumerical simulation of heat flow field and NOx concentrationfield in gas combustorZHANGNuonan，GUOYun（School of Mechanical and Automotive Engineering,Shanghai University of Engineering Science,
3、Shanghai 201620,China）Abstract：For a kind of gas burner,the appropriate size of the combustion chamber is selected to match,andthen a complete calculation area was formed.Natural gas fuel and combustion air were injected into thecombustion chamber through the gas burner to realize the complete combu
4、stion process.By simplifying the gasburner to complete geometric modeling,heat flow and combustion analysis software were used for numericalsimulation of the calculation area,and the Realizable k-turbulence model,ED combustion model,DOradiation model and thermal NOx generation model were selected to
5、 explore the influence of excess aircoefficient on the combustion process of the gas burner.The cloud image and distribution of temperaturefield and NOx concentration field were obtained and analyzed.When the excess air coefficient =1.1,the NOxconcentration reached the maximum.The research results c
6、an provide reference for relevant researchexperiments and design of low NOx burners.Key words：gas burner；nitrogen oxide；numerical simulation；excess air coefficient 收稿日期：2022 05 02基金项目：国家自然科学基金项目资助（51606116）；上海市科委长三角科技联合攻关领域项目资助（19195810800）作者简介：张诺楠（1996 ），女，在读硕士，研究方向为能源装备与过程控制.E-mail:通信作者：郭韵（1976 ），
7、女，副教授，博士，研究方向为清洁能源.E-mail: 第 36 卷 第 4 期上 海 工 程 技 术 大 学 学 报Vol.36 No.42022 年 12 月JOURNAL OF SHANGHAI UNIVERSITY OF ENGINEERING SCIENCEDec.2022 工业锅炉是氮氧化物工业源的重要组成之一，其在化石燃料利用中应用较多，产生的氮氧化物排放量较多.随着天然气锅炉数量越来越多1，其产生的氮氧化物量不容忽视，因此有效降低天然气锅炉产生的污染物 NOx排放量2刻不容缓.燃烧器是燃气锅炉的重要燃烧设备之一，通过将燃气燃料和助燃空气一起送入到燃烧室内部，进行后续燃烧过程，其对
8、燃烧过程有很大影响.对燃烧器燃烧性能的改善非常重要，它会影响产量、绿色环保、能源消耗等指标.目前燃烧技术研究方向主要聚集在低污染燃烧技术、低氮燃烧技术和燃烧器低氮燃烧技术.在现有研究基础上，对于污染物的减少重点还是要合理选择燃烧方式3.现今采用的大多数方法是设计或者改进优化燃气锅炉燃烧器的相关结构，与低氮燃烧技术相结合，对燃烧性能进行改良，使燃烧热效率得到提高，最终降低大气环境污染4.Su 等5对旋流燃烧器工业炉内过量空气比对燃烧性能的影响进行数值模拟.选用相关模型对进气温度为 1 273 K 的天然气进行模拟.研究过量空气比对炉内 NO 排放、温度和 CO 分布的影响，发现当过量空气比从 1
9、.05 增加到 1.25 时，热态 NO排放从 5106增加到 70106.Ariwibowo 等6采用计算流体力学（Computational Fluid Dynamics，CFD）方法研究非预混天然气空气在工业炉内的燃烧特性.何建军等7将一种新型低氮燃烧器用于燃气锅炉，结合多种低氮技术，烟气循环率达到20%40%，有效降低氮氧化物生成量，炉膛出口NOx浓度小于 30 mg/m3.史丹君等8采用计算颗粒流体力学（Computational Particle Fluid Dynamics，CPFD）进行低氮燃烧数值模拟，结合空气分级和低过量空气燃烧方法，使锅炉炉膛温度场分布均衡，O2浓度降低，
10、最终取得很好的 NOx降低效果.本研究对燃气燃烧器进行几何建模，匹配合适 的 燃 烧 室 模 型，用 ICEM 进 行 网 格 划 分，用ANSYS Fluent 2021 软件对计算区域模型进行数值模拟，选择过量空气系数作为可变条件，研究不同过量空气系数下温度场及氮氧化物 NOx浓度场的分布，探究过量空气系数对燃烧过程的影响.1 燃烧器结构设计燃气燃烧器使用三维软件 UG 按照等比例进行三维几何结构建模，整体长度为 820 mm，直径为 270 mm，空气与燃料的进口直径分别为 270 mm与 40 mm，燃 烧 器 头 部 内 径 为 140 mm，外 径 为210 mm，宽度为 50 m
11、m，分流管道内径为 10 mm.燃气燃料通过燃料入口进入燃料管道进入燃烧器内部，燃烧器头部有 4 个分流管，燃气燃料经过第一次的分级，分流流向 4 个管道，燃烧器头部的形状为圆形，形成一个环形腔体；当燃气燃料在整个环形腔体均匀充满后再向外发出喷射.燃烧器头部的环形腔体分别有 24 个一级燃料喷嘴与 3 个二级燃料喷嘴，喷嘴直径均为 8 mm，其三维结构模型如图 1 所示.燃料管道空气入口燃料入口燃料器头部旋流器分流管图 1 燃气燃烧器三维模型简图Fig.1 Three-dimensional model diagram of gas burner 2 模型建立 2.1 计算区域确定为使燃气燃料
12、的燃烧过程更加清晰明了且利于分析，需要匹配恰当的燃烧室模型，故建立一个圆柱体燃烧区域进行模拟燃气燃烧室，与原型燃气燃烧器组合在一起构建完整的计算区域，如图 2所示.燃烧室尺寸见表 1.为减少整体结构数值模拟的复杂性，燃气燃烧器减少了旋流器结构，并对主管道与燃料管道进行相应简化，其头部保留，包括燃料入口、空气入口及一二级燃料喷嘴，保存全部边界条件，故不会影响模拟结果的精准性.炉膛燃烧室烟气出口燃气燃烧器图 2 完整的流体计算区域示意图Fig.2 Complete schematic diagram of fluid computing area 348 上 海 工 程 技 术 大 学 学 报第
13、36 卷 2.2 控制方程自然界中的任何流体流动和化学反应过程都应遵循相对应的定律.在本研究的数值模拟中，均需要用到以下三大守恒方程9.质量守恒方程表达式为t+xi(ui)=Sm（1）其中，Sm的表达式为tvoldxdydz+AdA=0（2）动量守恒方程表达式为(ui)t+(uiuj)xi=pxi+ijxi+gi+Fi（3）能量守恒方程表达式为t(cpT)+xi(cpuiT)=xiprcpTxi+wsQs（4）式中：为微元流体的密度，kg/m3；t 为时间变量，s；xi为 i 方向的空间坐标分量，m；ui为 i 方向上的流体速度分量，m/s；uj为 j 方向上的流体速度分量，m/s；p 为微元
14、流体所受到的静压力，Pa；ij为应力张量，m/s；gi为 i 方向上的重力体积力；Fi为 i 方向上的外部体积力；cp为微元流体的定压比热容，J/(kgK)；T 为开尔文温度，K；为黏度，Pas；为导热系数，W/(mK)；wsQs为其他形式导入的能量.2.3 求解模型本研究利用 ANSYS 进行数值模拟，采用的求解模型见表 2.2.4 其他设置初始设置完成后，对计算迭代次数进行设置.选择 SIMPLE 算法.在数值模拟10 11过程中，对计算结果残差等进行监控，残差值设定见表 3.表 3 残差参考值设定Table 3 Residual reference value setting 残差监控D
15、OEnergy其他项设定值106106103 2.5 网格划分及无关性验证 2.5.1 网格划分在保证计算区域参数及工况设定一致情况下，对整体几何结构模型进行网格划分以及部分区域网格加密，选取 4 种不同的网格数量，见表 4.表 4 4 种网格数量与网格节点数量Table 4 Four kinds of grids numbers andgrid nodes numbers 网格编号网格数量网格节点数量1516 34590 88121 192 578203 88131 532 122261 44442 010 254342 451 2.5.2 网格无关性验证燃烧室沿 Y 轴方向火焰中部的温度与
16、速度分布如图 3 所示.Z 轴方向为燃烧室中心轴方向.由图可知，4 种不同网格数量的温度与速度分布曲线的变化趋势近似，说明本研究采用的数值模拟计算方法具有稳定性与可重复性.4 种网格数量的模拟结果有较小的差距，推测当网格数量再增加时，可以忽视对最终计算结果的影响.考虑计算机硬件性能、数值模拟计算结果的精度以及时间成本，最终选取网格数量 1 192 578 的模型进行模拟，最终网格划分如图 4 所示.3 过量空气系数对燃烧过程的模拟 3.1 过量空气系数燃烧反应中，过量空气系数 是反映燃料与空气比例的重要参数，过大还是过小都不利于燃烧过程的进行.在保持燃料燃烧充分的基础上，降低 表 1 燃烧室模型相关尺寸Table 1 Combustion chamber model related dimensions 半径/m长度/m半径长度比出口半径/m出口角度/（）0.41.80.2220.17550 表 2 采用的求解模型Table 2 Adoption of solving model 求解模型采用形式湍流模型Realizable k-模型燃烧模型涡耗散ED模型辐射模型DO辐射模型NOx形成