富氧燃烧中气体辐射模型对燃烧与换热数值模拟的影响

吴晓峰,范卫东

(上海交通大学 机械与动力工程学院,上海 200240)

摘 要:富氧燃烧过程中,由于使用再循环烟气代替空气中N2作为稀释剂,烟气中存在大量CO2和H2O。CO2和H2O作为非极性三原子分子,具有N2没有的辐射能力,导致富氧燃烧中气体辐射特性发生变化。在数值模拟过程中,气体辐射模型是一个重要的子模型。前人提出多种修改后适用于富氧燃烧的气体辐射模型,但不同气体辐射模型在不同富氧燃烧工况数值模拟中的影响尚未有统一研究。为了研究不同炉型下,气体燃烧和煤粉燃烧中气体辐射模型对燃烧换热模拟结果的影响,通过编程,将一种考虑CO影响的气体辐射模型以及文献中的6种典型气体辐射模型耦合入数值模拟计算。结果表明,在气体富氧燃烧中,气体辐射模型影响了火焰结构。同时,燃烧温度分布有所变化,不同模型结果之间差别最高可到500 K。气体与壁面之间的辐射换热受到影响。气体辐射模型对炉膛中心火焰区域影响较大,而对非火焰区域影响较小。在煤粉富氧燃烧过程中,当有效辐射层厚度在0.3 m左右时,如在100 kW下行炉中,气体辐射模型对煤粉燃烧数值模拟结果几乎没有影响。这可能是由于颗粒辐射在辐射换热计算中占主导地位。而当有效辐射层厚度在16 m左右时,如1 000 MW塔式炉中,气体辐射模型对炉内切圆燃烧火焰温度以及组分浓度影响较大,温度差别可到100 K左右。而气体辐射模型对炉膛中心模拟结果没有影响。

关键词:富氧燃烧;煤燃烧;气体辐射;灰气体加权模型

中图分类号:TK11TK16

文献标志码:A

文章编号:1006-6772(2021)02-0150-11

收稿日期:2020-09-06;责任编辑:白娅娜

DOI:10.13226/j.issn.1006-6772.CCUS20090601

移动阅读

基金项目:国家重点研发计划资助项目(2018YFB0605301)

作者简介:吴晓峰(1994—),男,江苏苏州人,博士研究生,从事富氧燃烧方面研究。E-mail:xfengwu@qq.com。

通讯作者:范卫东,教授,从事高效清洁燃烧方面研究。E-mail:wdfan@sjtu.edu.cn

引用格式:吴晓峰,范卫东.富氧燃烧中气体辐射模型对燃烧与换热数值模拟的影响[J].洁净煤技术,2021,27(2):150-160.WU Xiaofeng,FAN Weidong.Effect of gas radiation model on the numerical simulation of combustion and heat transfer in oxy-fuel combustion[J].Clean Coal Technology,2021,27(2):150-160.

Effect of gas radiation model on the numerical simulation of combustion and heat transfer in oxy-fuel combustion

WU Xiaofeng,FAN Weidong

(School of Mechanical and Power Engineering,Shanghai Jiao Tong University,Shanghai 200240,China)

Abstract:In oxy-fuel combustion,flue gas is mainly composed of CO2 and H2O because the recycled flue gas is used as diluent instead of N2 in the air. CO2 and H2O have an intensive radiation ability as non-polar triatomic molecule while N2 has no radiation ability,which leads to the change of gas radiation characteristics in oxy-fuel combustion. In the numerical simulation,gas radiation model is an important sub model. Several modified gas radiation models have been proposed for applying to oxy-fuel combustion,but the effect of different gas radiation model on the numerical simulation of oxy-fuel combustion in different conditions has not been unified. In order to research the effect of gas radiation model on the simulation results of combustion heat transfer in different sizes of furnace with gas or coal combustion,one new gas radiation model considering the effect of CO and other six typical models in the literatures were coupled into the numerical simulation by user defined function programming. The results show that:in the oxy-gas combustion,the flame structure is affected by gas radiation models. Meanwhile,the combustion temperature distribution changes and the difference between the results of different models can reach 500 K. Radiation heat transfer between gas and wall is also controlled by gas radiation models. The effect of gas radiation models on flame zone is rather large while it is ignored in the non-flame zone. In the oxy-coal combustion,when the effective radiation layer thickness is around 0.3 m such as in a 100 kWth down-fired furnace,gas radiation models almost have no effect on the coal combustion numerical simulation results. This may be because particle radiation plays a dominant role in the radiation heat transfer calculation. When the path length is around 16 m such as in a 1 000 MW tower type boiler,gas radiation models have a large effect on the combustion temperature and species concentration. The temperature difference can reach 100 K. However,the gas radiation model has no effect on the simulation results in the non-flame zone in middle of furnace.

Key words:oxy-fuel combustion;coal combustion;gas radiation;weighted sum of gray gas model

0 引 言

全球变暖是世界上主要环境问题之一[1]。在导致全球变暖的各种温室气体之中,CO2对全球温升影响最大,约70%的CO2排放来自化石燃料的燃烧[2]。未来20年,化石能源依然是我国能源供应的主要来源[3]。因此,为了降低碳排放,碳捕捉、利用与封存(Carbon capture,utilization and storage,CCUS)技术应运而生。CCUS不仅能将CO2捕捉封存,更可以将CO2“变废为宝”,被认为是应对全球气候变化的关键技术之一[4]

富氧燃烧作为一种燃烧中捕集CO2的技术,被认为是可应用于未来实际电站锅炉的CCUS技术。富氧燃烧中,空气分离器生产的纯氧与再循环烟气混合代替空气作为助燃剂送入炉膛。由于烟气中大量CO2代替空气中氮气作为稀释剂,当采用干式循环富氧燃烧时,炉膛出口排烟中CO2浓度理论上可以高达95%[5-7]。但当前阶段富氧燃烧的电厂仍处于设计研究阶段。对富氧燃烧的研究主要集中于实验室研究和CFD研究[8]。利用CFD工具对富氧燃烧锅炉炉膛内的温度、流场、传热情况进行预测是有效手段[9-15]。由于富氧燃烧中使用CO2代替空气中N2作为稀释剂,与常规空气燃烧相比,富氧工况下的燃烧特性和传热特性有很大不同。适用于空气燃烧的反应动力学、辐射模型等子模型由于其局限性还不能直接应用于富氧燃烧模拟中[16-17],需要针对富氧燃烧进行相应改进。辐射模型是燃烧模拟中的重要子模型之一,决定了颗粒以及混合气体在炉膛内的辐射换热量,进而影响了整个炉内的燃烧温度以及壁面热流。常规空气燃烧中,烟气成分主要以N2为主,N2作为非极性双原子分子不具有辐射能力。但富氧燃烧中的烟气以H2O和CO2为主,大量CO2作为稀释剂代替了常规空气燃烧中的N2。CO2作为极性的三原子分子,具有较强的辐射能力。在湿式循环富氧燃烧中,H2O含量可以提高到50%。H2O作为极性三原子分子,其辐射能力也很强。相比于空气燃烧,富氧燃烧中炉内气体辐射传热特性有所改变,火焰总发射率大大提高。在数值模拟中,考虑到计算开销,需使用简化的辐射模型来计算气体的辐射特性参数。灰气体加权模型(Weighted sum of gray gases model,WSGGM)是一种简单的气体辐射模型。在合理误差范围内,该模型计算速度非常快,易与计算流体力学(Computational fluid dynamics,CFD)软件耦合[18],被认为是最适用于三维空间数值模拟的气体辐射模型[19]

Kangwanpongpan等[20]在富氧工况下优化了WSGGM的参数,并在100 kW富氧煤燃烧炉膛上,比较了改进后的模型与Fluent软件的默认模型(Smith模型[21])在流速、氧浓度和温度预测精度上的差别。结果表明,改进后的辐射模型在流速和氧浓度的预测没有变化,但降低了烟气温度,更符合试验结果。Guo等[13]研究了改进后气体辐射模型对温度与组分浓度预测的影响,发现气体辐射对火焰结构、温度以及组分浓度影响较小,但颗粒辐射影响较大。冯艳[22]研究了Yin等[23]和Johansson等[24]提出的改进气体辐射模型对煤粉燃烧的影响,发现Yin的模型和Smith模型类似,而与Johansson模型有所差异。在理论计算上,Johansson等[24]、Guo等[18]、Kwangwanpongpan等[25]和Bordbar等[26]已证明了各自改进模型的准确性。但不同学者提出的富氧燃烧下气体辐射改进模型对三维燃烧数值模拟的影响以及不同气体辐射模型在各种工况下对燃烧模拟的影响尚未有统一结论。现有气体辐射模型均以CO2和H2O构成的富氧燃烧烟气为出发点进行计算。由于CO2存在,富氧燃烧中有强烈的焦炭气化反应,造成大量CO生成。若采取氧分级燃烧技术,主燃区内呈现富燃料燃烧,导致主燃区和还原区内CO出现。因此,在富氧燃烧炉膛内,CO是一种不可忽视的气体组分。一方面,CO作为非极性的双原子分子,本身存在一定的辐射能力。另一方面,作为CO2的前驱物,CO会占据一部分辐射气体体积,也会降低CO2的分压比,导致烟气辐射特性变化。但目前鲜见考虑CO气体辐射的模型。本文整理了近年来多种基于富氧燃烧改进的气体辐射模型,同时提出了一种考虑富氧燃烧中CO影响的改进模型。通过将这些模型耦合入CFD进行数值模拟,研究了在下行炉和实炉的气体和煤粉富氧燃烧过程中,不同气体辐射模型对燃烧和传热模拟结果的影响。

1 气体辐射特性

1.1 气体辐射特性计算模型

对于气体辐射特性的计算模型,目前主要应用有线性模型、谱带模型和总体模型[27-28]。不同气体辐射模型及分类如图1所示。

图1 气体辐射计算模型的分类

Fig.1 Categories of calculation models of gas radiation

线性模型一般指逐线性积分法(Line by line,LBL)。该方法是目前最准确的气体辐射特性计算方法,可作为其他方法的基准[29]。但采用逐线性积分法时,需知道气体分子每条谱线的详细光谱特性参数,包括谱线位置、谱线强度、谱线半宽、谱线跃迁能级能量等一系列参数。该方法计算量巨大,不适合工程问题求解。谱带模型建立谱线辐射强度与谱线参数的关系式,可从理论上推出气体在一定波数范围内的总辐射强度或发射率与谱线参数、谱线分布参数的关联式[29]。但谱带模型不能计算气体的光谱辐射特性,只能计算气体某一波长范围或全波长的辐射特性。谱带模型方法分为窄谱带模型和宽谱带模型两类。总体模型方法中主要包括k-分布方法和灰气体加权模型。k-分布方法在不同波段内计算不同波数的吸收系数,进行高斯积分,使计算效率大大提高。该方法可以和任意求解辐射传递方程的方法结合[27-28]

1.2 WSGGM及其改进模型

WSGGM模型最早由Hottle和Sarofim提出[30],后来Smith等[21]发展并拟合了当前使用最广泛的WSGGM模型参数。WSGGM模型是利用几种假想的灰气体以及一种透明气体去模拟实际气体的非灰辐射特性,并认为每种假想灰气体的吸收系数与波数无关。利用WSGGM模型,可通过以下公式计算混合气体的发射率ε

(1)

其中,Ng为假想的灰气体数量;αi为第i个假想灰气体的发射率权重因子;Ki为第i个假想灰气体的吸收系数;pa为吸收气体的总压;L为路径长度。α0为透明气体的权重因子,其吸收系数K0=0。

WSGGM模型避免了谱带模型的缺陷,可以与任意辐射传递方程求解方法结合,在工程上应用较广泛。Smith等[21]给出了一系列WSGGM多项式系数,利用Fletcher-Powell算法来求解多变量约束的非线性问题。基于指数宽谱带模型(Exponential wide-band mode,EWBM),Smith使用3种灰气体模拟混合气体的吸收系数。根据CO2或H2O分压不同,选择不同的关联系数计算混合气体的总发射率。该模型广泛使用于空气燃烧,成为了Fluent中的默认模型,用于计算混合气体的吸收系数。

在富氧燃烧中,许多学者提出了改进后的WSGGM。Yin[23,31]基于EWBM模型,利用4种灰气体给出了一种适用于空气和富氧多工况的WSGGM关联系数,基于CO2和H2O分压和总压力,给出了11种工况下WSGGM模型关联系数。Johansson等[24]使用统计窄谱带模型(Statistical narrow band model,SNBM)计算CO2和H2O混合气体发射率并提出了基于4种灰气体的WSGGM关联系数。但基于富氧燃烧中干、湿循环2种工况,只给出了0.125和1两种不同H2O/CO2摩尔比工况下的拟合系数。Guo等[18]结合了全光谱k-分布模型和WSGGM的特点,从k-分布中直接拟合了改进后WSGGM模型的权重因子和吸收系数。一些学者使用HITEMP2010数据库计算出混合气体总体发射率,用于WSGGM模型的拟合。Kangwanpongpan等[25]基于HITEMP2010光谱数据库,拟合了新的WSGGM关联系数,并给出了7组适用于不同H2O/CO2摩尔比工况下的关联系数。Bordbar等[26]使用5种灰气体拟合了一组新的WSGGM关联系数,该组关联系数适用于整个H2O/CO2摩尔比范围。

本文基于HITEMP2010光谱数据库,考虑了富氧燃烧中CO对气体辐射特性的影响,重新拟合了适用于浓度0~45% CO工况下的WSGGM关联系数。在富氧燃烧中,分级燃烧技术会引起主燃区的不完全燃烧,同时烟气中大量CO2引起的气化反应也会生成大量CO。因此,考虑CO影响,改进的WSGGM格式为

(2)

(3)

(4)

(5)

(6)

其中,Tr为相对温度;TTref分别为气体温度和参考气体温度;Mr为H2O和CO2分压比;P(H2O)、P(CO2)分别为H2O和CO2分压;bi,jci,j,kdi,k为多项式系数;jk分别为多项式项数。拟合后的关联系数见表1。

本文研究的7种WSGGM模型适用范围见表2,提出的模型命名为New-WSGGM。这些气体辐射模型均被重新编写成用户自定义函数(User defined function,UDF),耦合入Fluent软件进行数值模拟计算。

表1 考虑CO影响后的WSGGM关联系数

Table 1 Correlation coefficients of WSGGM considering effect of CO

关联系数ijk=0k=1k=2k=3k=4ci,j,k100.127 90.363 4-0.386 30.135 1-0.015 4110.242 8-0.306 40.601 7-0.242 10.029 412-0.433 80.196 9-0.456 60.184 7-0.022 3130.212 6-0.095 60.164 5-0.062 00.007 114-0.034 10.018 0-0.022 20.007 5-0.000 8200.164 21.029 8-0.811 30.256 4-0.027 7210.120 9-3.257 62.486 7-0.762 80.080 2220.386 12.684 4-2.063 40.633 5-0.066 623-0.347 8-0.845 60.666 2-0.206 60.021 8240.070 10.091 8-0.074 70.023 4-0.002 5300.171 20.042 40.015 7-0.013 30.002 0310.112 9-0.520 30.501 5-0.182 60.021 832-0.228 50.377 7-0.486 50.195 6-0.024 5330.091 6-0.080 30.166 0-0.073 80.009 734-0.011 20.002 6-0.019 90.009 8-0.001 3400.426 8-1.471 61.346 5-0.454 00.050 9410.183 33.512 5-3.563 91.248 7-0.143 042-0.693 6-2.644 33.065 4-1.119 50.130 8430.387 20.824 3-1.089 80.410 9-0.048 644-0.065 8-0.093 30.139 9-0.054 00.006 5di,k1—-0.075 710.413 8-9.338 23.330 4-0.395 02—0.003 20.116 2-0.107 40.038 0-0.004 43—8.971 1137.797 9-124.378 942.468 8-4.846 34—0.127 40.647 1-0.572 20.211 2-0.025 9

表2 本文研究的7种WSGGM模型

Table 2 Seven WSGGMs studied in this paper

文献温度/K参考温度/K分压比有效辐射层厚度/m灰气体数量基准模型Smith[21]600~2 400—1,20.001~103EWBMYin[23]500~3 000—0.125~40.001~604EWBMJohansson[24]500~2 5001 2000.125,10.01~603或4SNBMGuo[18]600~2 5001 5000.05~20.05~504FSKMKangwanpongpan[25]400~2 5002 0000.125~40.001~604LBLBordbar[26]300~2 4001 2000.01~40.01~604LBLNew-WSGGM400~3 0001 2000.01~40.01~604LBL

2 炉膛结构与计算模型

2.1 IFRF 0.8 MW富氧天然气燃烧炉膛

IFRF 0.8 MW富氧天然气燃烧炉膛如图2所示。IFRF炉膛截面尺寸为1.05 m×1.05 m,长度为3.44 m。天然气从燃烧器中心管喷入炉膛,中心管直径为16 mm。中心管周围布置一圈二次风,内径28 mm,外径36 mm。纯氧通过二次风管送入炉膛。尾部连接一个直径为500 mm的烟囱用于排烟。炉内有效辐射层厚度为0.8 m。具体炉膛结构可以参考文献[31-32]。采用结构化的六面体网格进行网格划分,炉膛中心采用网格加密以适应火焰处剧烈的温度梯度和浓度梯度。整体网格数量约83万。

图2 IFRF 0.8 MW富氧天然气燃烧炉膛

Fig.2 IFRF 0.8 MW oxy-natural gas combustion furnace

天然气由甲烷、乙烷等多种碳氢燃料组成,计算时简化为一种C1.122H4.244碳氢燃料。不同轴向距离的炉膛壁温不同,壁面温度分布Tw以多项式形式给出,即

Tw(z)=1 700.6+212.59z-46.669z2

(7)

采用离散坐标法(DO)求解辐射传输方程。湍流气相燃烧模型采用涡耗散概念模型(Eddy dissipation concept,EDC),采用细化后的Jones & Lindstedt 4步反应机理[33]见式(8)~(11)。相应的动力学参数见表3,其中A为指前因子,(kmol/m3)1-n/s;Ea为反应活化能,J/kmol;n为反应级数。

CH4+0.5O2CO+2H2

(8)

CH4+H2OCO+3H2

(9)

H2+0.5O2H2O,

(10)

CO+H2OCO2+H2

(11)

表3 天然气燃烧使用的四步化学反应动力学参数

Table 3 four-step reactions kinetic parameters in oxy-natural gas combustion

反应A/((kmol·m-3)1-n·s-1)Ea/(J·kmol-1)n(8)4.4×10111.26×108[CH4]0.5[O2]1.25(9)3.0×1081.26×108[CH4] [H2O](10)5.69×10111.465×108[H2] [O2]0.5(forward)(11)2.75×1098.36×107[CO] [H2O](forward)

2.2 100 kWth半工业规模炉膛

图3 100 kWth富氧燃烧炉膛结构

Fig.3 100 kWth oxy-fuel combustion furnace

该100 kWth半工业规模炉膛是一个垂直炉膛,来自德国亚琛工业大学传热传质研究所,主要用于富氧煤燃烧的相关研究[20,34]。100 kWth富氧燃烧炉膛结构如图3所示。炉膛内径0.4 m,长度2.1 m。炉内有效辐射层厚度为0.3 m。燃烧器由3个圆环喷口组成。其中一次风携带煤粉颗粒送入炉膛,风量17.6 kg/h,风温40 ℃。二次风为旋流风,风量为26.6 kg/h,旋流数设置为1.2,风温60 ℃。三次风用于扫气,风量为1.5 kg/h,风温60 ℃。在炉膛入口壁面设置了分级风,分级风风量为54.9 kg/h,风温900 ℃。

煤样采用莱茵河褐煤,煤质分析见表4。给煤量为6.5 kg/h,此次试验燃烧功率为40 kW。详细的炉膛结构和运行参数参见文献[33]。网格划分采用结构化的六面体网格,炉膛中心处加密,整体网格数量约98万。

表4 莱茵河褐煤煤质分析

Table 4 Properties of the Rhine lignite

工业分析/%MarAarVarFCar元素分析/%CarHarOarNarSar8.404.1046.6040.9067.404.2414.700.860.30

计算采用有限速率涡耗散模型(Finite rate eddy dissipation model,FREDM)进行湍流气相燃烧的模拟。采用多表面燃烧反应进行颗粒相的燃烧与气化模拟。详细反应机理[34]见式(13)~(18),相应的化学反应动力学参数见表5。湍流模型采用Realizable k-ε模型。采用DO模型进行辐射传输方程的求解。

(12)

CO+0.5O2CO2

(13)

CO2CO+0.5O2

(14)

H2+0.5O2H2O,

(15)

Cchar+0.5O2CO,

(16)

Cchar+CO22CO,

(17)

Cchar+H2OCO+H2

(18)

表5 煤粉燃烧化学反应机理动力学参数

Table 5 Reaction kinetic parameters of oxy-coal combustion

序号A/((kmol·m-3)1-n·s-1)Ea/(J·kmol-1)n(12)3.05×1075.065×107[CxHyOlNmSn] [O2](13)2.238×1064.187×107[CO] [O2]0.25 [H2O]0.5(14)1.095×10133.282×108[CO2] [O2]-0.25 [H2O]0.5(15)1×1088.36×106[H2] [O2](16)0.0057.4×107[O2](17)0.006 351.62×108[CO2](18)0.001 921.47×108[H2O]

2.3 1 000 MW塔式炉富氧燃烧

图4 某1 000 MW塔式锅炉结构

Fig.4 1 000 MW tower type boiler

基于某1 000 MW塔式锅炉,假想将其改造成富氧燃烧锅炉,如图4所示。锅炉宽度和深度均为23.16 m,高度为113.4 m。炉内有效辐射层厚度为16.21 m。燃烧器分为4组,3组主燃烧器和1组分离式燃尽风喷口。主燃烧器采用一、二次风相间分布。详细的锅炉结构和运行参数见文献[35]。富氧燃烧改造时,氧浓度保持21%不变,将空气中N2替换为CO2。数值模拟时,湍流模型采用标准的k-ε模型。采用DO模型求解辐射传输方程。挥发分热解采用双方程竞争模型。挥发分的气相湍流燃烧以及焦碳颗粒的氧化、气化反应与2.2节一致。采用网格分区划分,将整个炉膛分为灰斗区域、主燃区和受热面区域3部分。网格采用六面体结构化网格,主燃区网格进行加密,整体网格数量约190万。

3 结果与讨论

3.1 气体辐射模型在气体燃烧中的影响

不同气体辐射模型下炉膛截面的温度场分布如图5所示。Smith模型为Fluent软件自带的默认模型,修改后的模型采用本文提出的New-WSGGM模型。由图5可知,不同的气体辐射模型对天然气着火燃烧产生影响。Smith模型下,燃料流外侧着火较快。而在Modified气体辐射模型下,天然气气流外侧着火较慢一点。同时火焰高温区更加集中。

图5 气体辐射模型修改前后炉内温度场比较

Fig.5 Temperature distribution before and after the modification of gas radiation model

炉膛中心轴线处,气体温度和气体吸收系数随中心轴线距离的变化如图6所示。可知不同气体辐射模型对气体温度和气体吸收系数影响较大。由图6(a)可知,在高温火焰处,不同模型计算的火焰温度差别最高可达到500 K。由图6(b)可知,这些区域的吸收系数差别较大,吸收系数较高的模型导致气体温度降低。气体燃烧过程中,不同气体辐射模型计算导致的不同气体吸收系数对火焰的燃烧温度和火焰结构产生影响。当轴线距离大于1 m后,不同气体辐射模型计算的气体吸收系数有所差别,但气体温度差距不明显。从图7壁面辐射热流可以看出,辐射热流受气体吸收系数影响很大。气体燃烧过程中,与壁面之间的辐射换热完全由气体辐射特性决定。尽管气体辐射模型对气体温度的影响不大,但细小的差别对辐射换热产生了巨大影响。

图6 炉膛中心轴线处的气体温度和吸收系数

Fig.6 Gas temperature and absorption coefficient in the center line of the furnace

图7 炉膛壁面轴线处辐射热流

Fig.7 Radiation heat flux in the axial line of the furnace wall

z=0.22 m处径向距离上试验结果与模拟结果对比如图8所示。可知不同气体辐射模型在径向距离上的计算结果都与试验结果规律吻合较好。具体数值上的差别可能是入口气流速度分布、湍流模型以及化学反应机理上模拟不精确性导致。由于火焰集中在炉膛中心,模拟结果在炉膛中心处受气体辐射模型影响较大。在远离火焰至炉膛壁面处,辐射模型对温度或组分浓度的影响十分微小。

图8 z=0.22 m处径向距离模拟结果与试验结果对比

Fig.8 Comparison of simulation results and experimental results in the radial direction of z=0.22 m

3.2 气体辐射模型在煤燃烧半工业炉膛中的影响

选取另一个与IFRF燃气炉膛有效辐射层厚度类似的富氧燃煤锅炉,详细研究气体辐射模型在煤粉燃烧中的影响。气体辐射模型修改前后炉膛内的火焰情况如图9所示。可知煤燃烧过程中,气体辐射模型对火焰形状影响较小,只在高温火焰的尾迹处略有不同。在径向距离上,模拟结果与试验结果的比较如图10所示。温度与氧浓度均能较好符合在径向的变化规律。模拟精度受边界条件、湍流模型、反应机理等影响,导致结果与试验数据不够吻合。

图9 气体辐射模型修改前后炉内温度场

Fig.9 Temperature distribution before and after the modification of gas radiation model

图10 z=0.05 m处径向距离上模拟结果与试验结果比较

Fig.10 Comparison of simulation results and experimental results in the radial direction of z=0.05 m

图11 炉膛壁面轴向壁面辐射热流

Fig.11 Radiation heat flux in the axial line of the furnace wall

不同气体辐射模型下,壁面轴线处的辐射热流密度如图11所示。对比图10、11可知,不同气体辐射模型在富氧煤燃烧中对气体温度、组分以及辐射换热影响较小。不同模型的计算结果差别可以忽略不计。炉膛中心轴线处,气体吸收系数和颗粒吸收系数的变化情况如图12所示。可知气体吸收系数先降后升,在火焰高温区吸收系数较低。整体上看,气体吸收系数在0.49 m-1左右。而颗粒吸收系数则完全不同,在炉膛中心轴线处的颗粒吸收系数在0~20 m-1。不同气体辐射模型的数值模拟过程中,颗粒的发射率均设置为常数(0.9)。可以看出,富氧煤燃烧过程中,颗粒吸收系数远大于气体的吸收系数。在辐射换热中,颗粒的辐射特性起决定性作用。不同气体辐射模型计算的气体吸收系数会改变气体辐射换热特性。但由于颗粒存在,气体辐射特性参数的改变很难影响气体和颗粒与壁面之间的换热量。同时,在与颗粒耦合计算中,也很难改变颗粒的燃烧特性和辐射特性参数。因此,当有效辐射层厚度较小时,由于颗粒的高辐射特性参数,导致气体辐射模型几乎不会影响富氧煤燃烧过程中的燃烧和换热特性。

图12 炉膛中心轴线处的气体吸收系数和颗粒吸收系数

Fig.12 Gas and particle absorption coefficients in the center line of furnace

3.3 气体辐射模型在1 000 MW实炉中的影响

对一台1 000 MW塔式锅炉进行富氧燃烧改造,氧浓度保持21%不变。相比100 kW富氧煤燃烧炉膛,1 000 MW的塔式锅炉中有效辐射层厚度增大到16.21 m。不同气体辐射模型下,炉内中心截面处的温度分布以及炉膛壁面处的辐射热流分布如图13所示。可知气体辐射模型对炉内整体温度分布影响很小,对炉膛壁面辐射热流影响也较小,但燃尽风后炉膛壁面处的吸热情况略有不同。

图13 气体辐射模型修改前后的炉内温度场和壁面辐射热流分布

Fig.13 Temperature and wall radiation heat flux distribution before and after the modification of gas radiation model

典型一、二次风平面中心线处的温度和CO2浓度分布分别如图14、15所示。可知在不同的气体辐射模型下,炉膛中心处的温度分布和CO2浓度分布均未发生变化。而在靠近壁面的切圆部分,即炉内高温火焰处,不同模型对温度和组分影响产生较大变化。在燃烧最剧烈的火焰部分,温度和组分浓度变化最剧烈,气体辐射模型对温度和浓度的计算产生影响。不同模型之间,最高温度偏差在二次风界面中心线可以达到100 K,最大浓度偏差在二次风界面中心线也可以达到5%。

图14 典型温度分布

Fig.14 Typical temperature distribution

图15 典型CO2浓度分布

Fig.15 Typical CO2 concentration distribution

典型一次风平面中心线处的气体吸收系数如图16(a)所示。在大炉膛空间内,气体吸收系数很小,约0.042 m-1,而变化区间在0.004 m-1左右。不同气体辐射模型最终计算的气体吸收系数区别很小。相同位置的颗粒吸收系数如图16(b)所示,可知颗粒吸收系数约在6 m-1变化,该数值仍远高于气体吸收系数。

与3.2节半工业规模小炉膛不同的是,在大炉膛内,气体辐射模型对炉内气体辐射特性的影响很小,但会导致最终计算的高温火焰温度和组分出现较大差别。这说明,当有效辐射层厚度较大时,气体辐射模型对高温火焰区的温度、组分浓度等燃烧特性产生较大影响,这在有效辐射层厚度较小时无法体现。

图16 典型一次风截面中心线处的气体吸收系数和颗粒吸收系数

Fig.16 Gas and particle absorption coefficients in the center line of cross section of primary flow

4 结 论

本文提出了一种考虑CO影响的气体辐射模型,并将文献中其他6种改进后的气体辐射模型通过UDF方式耦合入数值模拟计算。在下行炉和实炉中,研究了气体和煤粉富氧燃烧过程中,气体辐射模型对数值模拟结果的影响。主要结论如下:

1)富氧气体燃烧中,气体辐射模型会影响气体燃料流的着火情况,改变火焰温度分布,温度差别最高达到500 K。同时影响了气体与壁面之间的辐射换热。气体辐射模型对炉膛轴线处的高温火焰影响较大,而对靠近壁面的非火焰区域几乎没有影响。

2)在下行炉的煤粉富氧燃烧过程中,颗粒吸收系数远大于气体吸收系数,在辐射传热中占主导地位。当有效辐射层厚度较小时(0.3 m),不同气体辐射模型在燃烧数值模拟过程中,对温度、组分浓度以及换热方面的预测影响可以忽略不计。

3)在实炉中的煤粉富氧燃烧中,有效辐射层厚度较大(16.21 m)时,气体辐射模型对炉内切圆燃烧高温火焰区域影响很大,而对炉膛中心处的计算结果没有影响。应用不同气体辐射模型,火焰温度差别可到100 K,组分浓度差别可达5%。

4)在富氧气体燃烧数值模拟以及大炉膛内富氧煤粉燃烧数值模拟中,需要考虑气体辐射模型对数值模拟结果的影响。

参考文献(References):

[1] 刘宇,曹江,朱声宝.挑战全球气候变化——二氧化碳捕集与封存[J].前沿科学,2010,4(1):40-51.

LIU Yu,CAO Jiang,ZHU Shengbao. Challenging climate change:Carbon dioxide capture and storage[J]. Frontier Science,2010,4(1):40-51.

[2] 潘一,梁景玉,吴芳芳,等.二氧化碳捕捉与封存技术的研究与展望[J].当代化工,2012,41(10):1072-1075,1078.

PAN Yi,LIANG Jingyu,WU Fangfang,et al. Research and prospect of the carbon dioxide capture and storage technology[J]. Contemporary Chemical Industry,2012,41(10):1072-1075,1078.

[3] 何铮,李瑞忠.未来20年中国能源需求预测[J]. 当代石油石化,2016,24(9):1-8.

HE Zheng,LI Ruizhong. Prediction of China's energy consumption in the future 20 year[J]. Petroleum & Petrochemical Today,2016,24(9):1-8.

[4] 张卫东,张栋,田克忠.碳捕集与封存技术的现状与未来[J].中外能源,2009,14(11):7-14.

ZHANG Weidong,ZHANG Dong,TIAN Kezhong. Carbon capture and sequestration technology[J]. Sino-global Energy,2009,14(11):7-14.

[5] HEES J,ZABRODIEC D,MASSMEYER A,et al. Experimental investigation into the influence of the oxygen concentration on a pulverized coal swirl flame in oxy-fuel atmosphere[J]. Fuel,2019,240:64-74.

[6] JIANG X,HUANG X,LIU J,et al. NOx Emission of fine- and superfine- pulverized coal combustion in O2/CO2 atmosphere[J]. Energy & Fuels,2010,24:6307-6313.

[7] ZHONG S,ZHANG F,PENG Z,et al. Roles of CO2 and H2O in premixed turbulent oxy-fuel combustion[J]. Fuel,2018,234:1044-1054.

[8] WALL T,STANGER R,LIU Y. Gas cleaning challenges for coal-fired oxy-fuel technology with carbon capture and storage[J]. Fuel,2013,108:85-90.

[9] EDGE P,GHAREBAGHI M,IRONS R,et al. Combustion modelling opportunities and challenges for oxy-coal carbon capture technology[J]. Chemical Engineering Research and Design,2011,89:1470-1493.

[10] EDGE P,GUBBA S R,MA L,et al. LES modelling of air and oxy-fuel pulverised coal combustion—impact on flame properties[J]. Proceedings of the Combustion Institute,2011,33:2709-2716.

[11] EDGE P J,HEGGS P J,POURKASHANIAN M,et al. An integrated computational fluid dynamics-process model of natural circulation steam generation in a coal-fired power plant[J]. Computers & Chemical Engineering,2011,35:2618-2631.

[12] GUO J,HU F,JIANG X,et al. Experimental and numerical investigations on heat transfer characteristics of a 35 MW oxy-fuel combustion boiler[J]. Energy Procedia,2017,114,481-489.

[13] GUO J,HU F,JIANG X,et al. Effects of gas and particle radiation on IFRF 2.5 MW swirling flame under oxy-fuel combustion[J]. Fuel,2020,263:116634.

[14] GUO J,LIU Z,HUANG X,et al. Experimental and numerical investigations on oxy-coal combustion in a 35 MW large pilot boiler[J]. Fuel,2017,187:315-327.

[15] GUO J,LIU Z,WANG P,et al. Numerical investigation on oxy-combustion characteristics of a 200 MWe tangentially fired boiler[J]. Fuel,2015,140:660-668.

[16] WALL T,LIU Y,SPERO C,et al. An overview on oxyfuel coal combustion:State of the art research and technology development[J]. Chemical Engineering Research and Design,2009,87:1003-1016.

[17] WALL T F. Combustion processes for carbon capture[J]. Proceedings of the Combustion Institute,2007,31:31-47.

[18] GUO J,LI X,HUANG X,et al. A full spectrum k-distribution based weighted-sum-of-gray-gases model for oxy-fuel combustion[J]. International Journal of Heat and Mass Transfer,2015,90:218-226.

[19] BECHER V,GOANTA A,SPLIETHOFF H. Validation of spectral gas radiation models under oxyfuel conditions - Part C:Validation of simplified models[J]. International Journal of Greenhouse Gas Control,2012,11:34-51.

[20] KANGWANPONGPAN T,CORRA da Silva R,KRAUTZ H J. Prediction of oxy-coal combustion through an optimized weighted sum of gray gases model[J]. Energy,2012,41:244-251.

[21] SMITH T F,SHEN Z F . Evaluation of coefficients for the weighted sum of gray gases model,ASME and American Institute of Chemical Engineers[C]. 20th National Heat Transfer Conference,1982,104:602-608.

[22] 冯艳.富氧气氛下辐射模型对煤粉燃烧数值模拟影响特性的研究[D]. 马鞍山:安徽工业大学,2017:73-76.

FENG Yan. Study on the influence of radiation model on the numerical simulation of pulverized coal combustion in oxygen-enriched atmosphere[D].Maanshan:Anhui University of Technology,2017:73-76.

[23] YIN C,JOHANSEN L C R,ROSENDAHL L A,et al. New weighted sum of gray gases model applicable to computational fluid dnamics(CFD) modeling of oxy-fuel combustion:Derivation,validation,and implementation[J]. Energy & Fuels,2010,24:6275-6282.

[24] JOHANSSON R,ANDERSSON K,LECKNER B,et al. Models for gaseous radiative heat transfer applied to oxy-fuel conditions in boilers[J]. International Journal of Heat and Mass Transfer,2010,53:220-230.

[25] KANGWANPONGPAN T,FRANÇA F H R,CORRA da Silva R,et al. New correlations for the weighted-sum-of-gray-gases model in oxy-fuel conditions based on HITEMP 2010 database[J]. International Journal of Heat and Mass Transfer,2012,55:7419-7433.

[26] BORDBAR M H,WCEL G,HYPPNEN T.. A line by line based weighted sum of gray gases model for inhomogeneous CO2-H2O mixture in oxy-fired combustion[J]. Combustion and Flame,2014,161:2435-2445.

[27] 齐永刚. 富氧燃烧过程炉膛烟气辐射特性与辐射传热的数值模拟[D].武汉:华中科技大学,2011.

QI Yonggang,Numerical simulation of radiative heat transfer and radiative properties of gases in oxy-fuel conditions [D]. Wuhan:Huazhong University of Science and Technology,2011.

[28] 姚放,加压富氧下煤粉燃烧过程的数值模拟[D]. 北京:北京交通大学,2014.

YAO Fang,Numerical simulation of the pressurized oxy-coal combustion process[D]. Beijing:Beijing Jiao Tong University,2014.

[29] 谈和平,夏新林,刘林华,等. 红外辐射特性与传输的数值计算:计算热辐射学[M]. 哈尔滨:哈尔滨工业大学出版社,2006:69-86.

TAN Heping,XIA Xinlin,LIU Linhua,et al,Numerical calculation of infrared radiation characteristics and transmission:Calculation of thermal radiation science[M]. Harbin:Harbin Institute of Technology Press,2006:69-86.

[30] HOTTEL H C,SAROFIM A F. Radiative transfer[M]. New York:McGraw-Hill Book Company,1967.

[31] YIN C,ROSENDAHL L A,K☞R S K . Chemistry and radiation in oxy-fuel combustion:A computational fluid dynamics modeling study[J]. Fuel,2011,90:2519-2529.

[32] YIN C . Prediction of air-fuel and oxy-fuel combustion through a generic gas radiation property model[J]. Applied Energy,2017,189:449-459.

[33] JONES W P,LINDSTEDT R P. Global reaction schemes for hydrocarbon combustion[J].Combustion and Flame,1988,73:233-249.

[34] TOPOROV D,BOCIAN P,HEIL P,et al. Detailed investigation of a pulverized fuel swirl flame in CO2/O2 atmosphere[J]. Combustion and Flame,2008,155:605-618.

[35] WU X,FAN W,LIU Y,et al. Numerical simulation research on the unique thermal deviation in a 1 000 MW tower type boiler[J]. Energy,2019,173:1006-1020.