气体燃烧室中热流场与生成NO_x浓度场的数值模拟研究_张诺楠.pdf

文章编号：1009 444X（2022）04 0347 06气体燃烧室中热流场与生成 NOx浓度场的数值模拟研究张诺楠，郭韵（上海工程技术大学 机械与汽车工程学院,上海 201620）摘要：对于一种燃气燃烧器，选取合适的燃烧室尺寸进行匹配，进而构成一个完整的计算区域.天然气燃料和助燃空气通过该燃气燃烧器喷入燃烧室实现完整燃烧过程.通过简化燃气燃烧器完成几何建模，使用热流与燃烧分析软件对计算区域进行数值模拟，选用 Realizable k-湍流模型、ED 燃烧模型、DO 辐射模型及热力型 NOx生成模型，探究过量空气系数 对该燃气燃烧器燃烧过程的影响，得出其温度场及生成 NOx浓度场的相关云图及分布情况并进行分析.当过量空气系数 =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,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 combustion 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 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 can 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 ），女，副教授，博士，研究方向为清洁能源.E-mail: 第 36 卷 第 4 期上 海 工 程 技 术 大 学 学 报Vol.36 No.42022 年 12 月JOURNAL OF SHANGHAI UNIVERSITY OF ENGINEERING SCIENCEDec.2022 工业锅炉是氮氧化物工业源的重要组成之一，其在化石燃料利用中应用较多，产生的氮氧化物排放量较多.随着天然气锅炉数量越来越多1，其产生的氮氧化物量不容忽视，因此有效降低天然气锅炉产生的污染物 NOx排放量2刻不容缓.燃烧器是燃气锅炉的重要燃烧设备之一，通过将燃气燃料和助燃空气一起送入到燃烧室内部，进行后续燃烧过程，其对燃烧过程有很大影响.对燃烧器燃烧性能的改善非常重要，它会影响产量、绿色环保、能源消耗等指标.目前燃烧技术研究方向主要聚集在低污染燃烧技术、低氮燃烧技术和燃烧器低氮燃烧技术.在现有研究基础上，对于污染物的减少重点还是要合理选择燃烧方式3.现今采用的大多数方法是设计或者改进优化燃气锅炉燃烧器的相关结构，与低氮燃烧技术相结合，对燃烧性能进行改良，使燃烧热效率得到提高，最终降低大气环境污染4.Su 等5对旋流燃烧器工业炉内过量空气比对燃烧性能的影响进行数值模拟.选用相关模型对进气温度为 1 273 K 的天然气进行模拟.研究过量空气比对炉内 NO 排放、温度和 CO 分布的影响，发现当过量空气比从 1.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浓度降低，最终取得很好的 NOx降低效果.本研究对燃气燃烧器进行几何建模，匹配合适 的 燃 烧 室 模 型，用 ICEM 进 行 网 格 划 分，用ANSYS Fluent 2021 软件对计算区域模型进行数值模拟，选择过量空气系数作为可变条件，研究不同过量空气系数下温度场及氮氧化物 NOx浓度场的分布，探究过量空气系数对燃烧过程的影响.1 燃烧器结构设计燃气燃烧器使用三维软件 UG 按照等比例进行三维几何结构建模，整体长度为 820 mm，直径为 270 mm，空气与燃料的进口直径分别为 270 mm与 40 mm，燃 烧 器 头 部 内 径 为 140 mm，外 径 为210 mm，宽度为 50 mm，分流管道内径为 10 mm.燃气燃料通过燃料入口进入燃料管道进入燃烧器内部，燃烧器头部有 4 个分流管，燃气燃料经过第一次的分级，分流流向 4 个管道，燃烧器头部的形状为圆形，形成一个环形腔体；当燃气燃料在整个环形腔体均匀充满后再向外发出喷射.燃烧器头部的环形腔体分别有 24 个一级燃料喷嘴与 3 个二级燃料喷嘴，喷嘴直径均为 8 mm，其三维结构模型如图 1 所示.燃料管道空气入口燃料入口燃料器头部旋流器分流管图 1 燃气燃烧器三维模型简图Fig.1 Three-dimensional model diagram of gas burner 2 模型建立 2.1 计算区域确定为使燃气燃料的燃烧过程更加清晰明了且利于分析，需要匹配恰当的燃烧室模型，故建立一个圆柱体燃烧区域进行模拟燃气燃烧室，与原型燃气燃烧器组合在一起构建完整的计算区域，如图 2所示.燃烧室尺寸见表 1.为减少整体结构数值模拟的复杂性，燃气燃烧器减少了旋流器结构，并对主管道与燃料管道进行相应简化，其头部保留，包括燃料入口、空气入口及一二级燃料喷嘴，保存全部边界条件，故不会影响模拟结果的精准性.炉膛燃烧室烟气出口燃气燃烧器图 2 完整的流体计算区域示意图Fig.2 Complete schematic diagram of fluid computing area 348 上 海 工 程 技 术 大 学 学 报第 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 为微元流体所受到的静压力，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 残差监控DOEnergy其他项设定值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 轴方向火焰中部的温度与速度分布如图 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形成模型热力型NOx模型第 4 期张诺楠 等：气体燃烧室中热流场与生成 NOx浓度场的数值模拟研究 349 NOx生成量，选择合适的过量空气系数范围.工况为满负荷，燃料及空气入口类型均为速度入口，空气温度设为 298 K，烟气出口边界采用pressure-outlet，燃烧室壁面条件为绝热.燃气体积流量为 100 m3/h，一二级燃料占比为 81，燃料管道进口速度为 50 m/s.使燃气燃烧器能在一种比较好的运行工况下运行，设置 5 种过量空气系数 进行模拟计算以及云图分析.3.2 模拟结果分析 3.2.1 温度场分析分别取 =1.0、1.1、1.2、1.3 和 1.4 时，轴向与径向温度云图如图 5 和图 6 所示.当 =1.0 时，火焰形状比较分散且火焰长度比较长，燃烧室内部的高温分布区域范围较大，燃烧反应过程较剧烈，且左下角角落有部分较高温区域，这是因为燃烧器的水平布置导致高温烟气在燃烧室的下边角落形成聚集，散热效果不够好.从轴向和径向截图的温度云图可以看出，随着 的增加，混合气体的体积以及流速增加，燃烧室中部上边的较高温区域被冲散，高温区域变得更加集中，火焰形状不再较为分散，火焰长度逐渐缩 0.400.20.4Y/m(a)温度分布曲线0.2Temperature/K8004001 2002 0001 600516 3451 192 5781 532 1222 010 2540.400.20.4Y/m(b)速度分布曲线0.2Velocity/(ms1)20100305040516 3451 192 5781 532 1222 010 254图 3 4 种不同网格数量的温度分布及速度分布曲线Fig.3 Four kinds of temperature distribution and velocitydistribution curves with different grid numbers 图 4 最终确定的网格划分示意图Fig.4 Final grid division diagram =1.2=1.1=1.02.98e+024.40e+025.82e+027.25e+028.67e+021.01e+031.15e+031.29e+031.44e+031.58e+031.72e+031.86e+032.00e+032.10e+03Contour1Static temperature/K2.98e+024.41e+025.84e+027.28e+028.71e+021.01e+031.16e+031.30e+031.44e+031.59e+031.73e+031.87e+032.02e+032.11e+03Contour1Static temperature/K2.98e+024.42e+025.86e+027.30e+028.75e+021.02e+031.16e+031.31e+031.45e+031.60e+031.74e+031.88e+032.03e+032.12e+03Contour1Static temperature/K=1.4=1.32.98e+024.41e+025.83e+027.26e+028.68e+021.01e+031.15e+031.30e+031.44e+031.58e+031.72e+031.87e+032.01e+032.10e+03Contour1Static temperature/K2.98e+024.42e+025.86e+027.30e+028.74e+021.02e+031.16e+031.31e+031.45e+031.59e+031.74e+031.88e+032.02e+032.12e+03Contour1Static temperature/K图 5 不同过量空气系数下计算域中心截面的轴向温度云图（Y=1 m）Fig.5 Axial temperature nephogram of central section of calculation domain under different excess air coefficients(Y=1 m)350 上 海 工 程 技 术 大 学 学 报第 36 卷短，左下角的较高温区域逐渐减小.燃烧室内部的平均温度偏高，并且随着 的逐渐增加先升高后降低；随着 逐渐增加，温度峰值的变化趋势先增加后减小，但增加和减小的幅度不是太大，都稳定在 2 100 K 左右，同时还能看到燃烧室径向的高温区域范围逐渐减小，且高温峰值也在逐渐减小.3.2.2 NOx浓度场分析大气污染物 NOx的主要来源之一就是工业锅炉燃烧.工业锅炉在进行燃烧之后，其排放的氮氧化物 NOx给人类带来极大的影响.随着人们对绿色环保理念的认可及国家对污染物排放标准的严格把控12，目前已形成一个相关的重要研究方向，即如何进行 NOx的控制排放及低 NOx燃烧理论与技术研究等.在探究不同过量空气系数对燃烧过程的影响时，除温度场分析，还需要对氮氧化物 NOx浓度场分析.本研究氮氧化物 NOx主要为热力型 NOx.图 7 为 取值不同时，燃烧室中心截面的 NOx浓度场分布云图.=1.2=1.1=1.00.00e+004.94e 069.88e 061.48e 051.98e 052.47e 052.97e 053.46e 053.95e 054.45e 054.94e 055.44e 055.93e 056.26e 05Contour1Mole fraction of pollutant no0.00e+005.04e 061.01e 051.51e 052.02e 052.52e 053.02e 053.53e 054.03e 054.54e 055.04e 055.54e 056.05e 056.38e 05Contour1Mole fraction of pollutant no0.00e+005.15e 061.03e 051.55e 052.06e 052.58e 053.09e 053.61e 054.12e 054.64e 055.15e 055.67e 056.19e 056.53e 05Contour1Mole fraction of pollutant no=1.4=1.30.00e+002.30e 064.61e 066.91e 069.22e 061.15e 051.38e 051.61e 051.84e 052.07e 052.30e 052.53e 052.77e 052.92e 05Contour1Mole fraction of pollutant no0.00e+002.31e 064.63e 066.94e 069.25e 061.16e 051.39e 051.62e 051.85e 052.08e 052.31e 052.54e 052.78e 052.93e 05Contour1Mole fraction of pollutant no图 7 不同过量空气系数下计算域中心截面的氮氧化物分布云图Fig.7 Clouds of nitrogen oxide distribution in the central section of calculation domain under different excess air coefficients 由图中可以看出，NOx高浓度的区域主要出现在高温区域，高温区域在上文中的温度分布云图也可以进行验证，其中温度云图中的高温区域的平均温度约为 1 400 K.因为天然气燃料和助燃空气进行激烈的燃烧反应，大量的燃烧热被释放出来，温度迅速升高.由于燃烧温度越高，空气中O2与 N2在高温下反应剧烈，热力型 NOx形成，所以氮氧化物主要存在于燃烧室的高温分布区域.为便于直观地进行数据分析，计算得出在不同工况下燃烧室出口的氮氧化物浓度，分布曲线如图 8 所示.从图中可以发现，当 为 1.01.1 时，NOx浓度呈上升趋势，当 为 1.11.4 时，NOx浓度开始逐渐降低.这是因为随着 从 1.0 增大到1.1 时，更多的空气进入燃烧室内，燃烧室的氧气含量迅速增多，更多的氧气参加生成氮氧化物的化学反应，造成燃烧室出口处的 NOx浓度提高.但 =1.2=1.1=1.0=1.4=1.3Contour1Static temperature/K2.09e+031.99e+031.85e+03 1.71e+031 57e+031 43e+031.29e+031.15e+031.01e+038.68e+027.27e+025 87e+024.46e+023.06e+02/KContour1Static temperature/K2.09e+031.99e+031.85e+03 1.71e+031 57e+031 43e+031.29e+031.15e+031.01e+038.68e+027.27e+025 87e+024.46e+023.06e+02/KContour1Static temperature/K2.00e+031.91e+031.78e+03 1.64e+031 51e+031 38e+031.24e+031.11e+039.74e+028.40e+027.06e+025.72e+024.38e+023.05e+02/KContour1Static temperature/K1.99e+031.91e+03 1.77e+031.64e+031.50e+031.37e+031.24e+031.10e+039.71e+028.38e+027.05e+025.71e+024.38e+023.04e+02/KContour1Static temperature/K2.00e+031.91e+031.78e+03 1.65e+031.51e+031.38e+031.24e+031.11e+039.75e+028.41e+027.07e+025.72e+024.38e+023.04e+02/K图 6 不同过量空气系数下计算域中心截面的径向温度云图（Z=1 m）Fig.6 Radial temperature nephogram of central section of calculation domain under different excess air coefficients(Z=1 m)第 4 期张诺楠 等：气体燃烧室中热流场与生成 NOx浓度场的数值模拟研究 351 随着 进一步增大，燃烧室内的燃烧过程虽然处于富氧状态，但是大量的冷空气喷入到燃烧室，其温度相对于燃烧室内部的高温来说很低，造成燃烧室的高温区域减小且温度降低，从而抑制氮氧化物生成，NOx浓度降低.即随着 的增加，燃烧室出口 NOx浓度的变化趋势是先增加后减小.当=1.1 时，NOx浓度最大，约为 37 106.的发展规律与文献 13 中的相关研究结论接近.1.01.21.11.4过量空气系数 1.3标况下 NOx 浓度 1062420263836322822183034图 8 不同过量空气系数下燃烧室出口的 NOx浓度分布曲线Fig.8 Distribution curve of NOx concentration at outlet ofcombustion chamber under different excess air coefficients 4 结语当过量空气系数 为 1.01.1 时，氮氧化物浓度逐渐增加，=1.0 时的 NOx浓度低于 =1.1时的 NOx浓度，但此时导致天然气燃烧不完全甚至燃烧不稳定的危险，还可能会引起燃烧效率降低；当 为 1.11.4 时，氮氧化物浓度逐渐减少，冷空气会随着 的增大进一步增多，进入过多的冷空气会吸收更多的燃烧热，导致燃烧温度降低，燃烧室出口的 O2浓度增大，过大的 也会导致燃烧效率降低.太大或太小都存在相应问题，所以在实际工程应用中，要选择合适的.参考文献：车光兰.天然气锅炉低氮氧化物燃烧技术研究J.工业加热，2020，49（3）：13 16.1 姜涌,夏明明,覃绍亮,等.热力型NOx的抑制J.电站系统工程，2005，21（2）：2.2 唐黔,卢先德,王树海,等.工业炉窑烟气污染及其治理(上)J.工业炉，1995（4）：42 47,58.3 CHEN C,CHENG S,GUO X.Hazard control of NOx inhot stoveJ.Journal of Iron and Steel Research(International)，2014，21（3）：306 311.4 SU Y X,WANG W H.Combustion performance and NOemission in industrial furnace under preheated air conditionwith different excess air ratioJ.Advanced MaterialsResearch，2011（402）：463.5 ARIWIBOWO T H,SAFITRA A G,SANTOSO H.Investigation on combustion characteristics ofnonpremixed combustion in an 85 MWe wall firing boilerunder various air-fuel ratio using CFDJ.2017International Electronics Symposium on EngineeringTechnology and Applications(IES-ETA)，2017：20 26.6 何建军,陆全华,王春华.用于燃气锅炉烟气回流内外双循环低氮燃烧器与燃气锅炉：中国，CN208859604UP.2019-05-14.7 史丹君,张建.循环流化床锅炉低氮燃烧的CPFD数值模拟J.锅炉技术，2020，51（1）：30 36.8 岑可法，姚强，骆仲泱，等.燃烧理论与污染控制M.北京:机械工业出版社,2004.9 刘效洲,朱光羽.小型天然气锅炉低氮燃烧技术的数值模拟和实验研究J.热能动力工程，2022，37（1）：132142.10 潘书婷.低NOx天然气燃烧器的数值模拟研究D.武汉:武汉科技大学,2021.11 苏涛.大气中氮氧化物的形成及防治J.科学咨询(决策管理)，2009（6）：43 44.12 王明雨.超低氮侧壁燃气燃烧器的数值模拟研究D.南京:南京师范大学,2020.13（编辑：韩琳）352 上 海 工 程 技 术 大 学 学 报第 36 卷