煤粉燃烧过程的数值模拟
水泥回转窑内煤粉燃烧过程的数值模拟_王乃帅
水泥回转窑内煤粉燃烧过程的数值模拟王乃帅 温 治 楼国锋 刘训良 郑坤灿 张 欣(北京科技大学 北京 100083)摘 要:本文应用标准kε−湍流模型、随机颗粒轨道模型、即混即燃模型、P1辐射模型以及多块非均匀结构化网格划分技术,对配有四风道燃烧器的水泥窑内煤粉燃烧过程进行了数值仿真,研究了多种操作参数对窑内燃烧工况的影响规律,提出了相应的优化操作制度。
关键词:水泥回转窑 煤粉燃烧 数值模拟Numerical Simulation of Pulverized Coal Combustion Procedurein Cement Rotary KilnWang Nai-shuai、Wen Zhi、Lou Guo-feng、Liu Xun-liang、Zheng Kun-can、Zhang Xin (University of Science and Technology Beijing, Beijing 100083)Abstract: This paper presents a 3-D numerical study on the effects of various operational parameters on flame feature and temperature profile in cement rotary kiln with CFD software which combined the standard K-Epsilon turbulent model, random track mode, P1 radiation model and unstructured grid. A four-air channel coal burner and a cement rotary kiln are investigated. The optimal operational parameters are put forward.Key words: cement rotary kiln; pulverized coal combustion; numerical simulation1 引言水泥窑是一种以燃料燃烧、高温传热、生料反应以及生料输送为主要功能的水泥熟料生产设备,煅烧生料所需的热量来自于燃烧器向窑内直接喷煤燃烧[1]。
高炉内煤粉燃烧速度场的数值模拟
2 0 1 4年 2月
包
钢
科
技
Vo 1 . 4 0, No . 1 F e b r u a r y, 2 0 1 4
S c i e n c e a n d T e c h n o l o g y o f B a o t o u S t e e l
i t y i f e l d a n d t r a c k o f p a r t i c l e i n t h e p r o c e s s o f p u l v e i r z e d c o a l c o mb u s t i o n a r e a n a l y z e d a n d s t u d i e d . T h e s i mu l a t e d r e s u l t s i n - d i c a t e t h a t t h e r e s u l t s w i t h mo me n t u m e q u a t i o n . s t nd a a r d d o u b l e—e q u a t i o n t u r b u l e n c e mo d e l a n d d i s c r e t e p h a s e mo d e l a r e
Ab s t r a c t : I n t h e p a p e r , t h e 3 D ma t h e ma t i c a l mo d e l f o r p u l v e i r z e d c o a l c o mb u s t i o n i n t h e r a c e w a y o f b l a s t f u r n a c e i s e s — t a b l i s h e d w i t h f i n i t e v o l u me me t h o d a n d c o mme r c i l a s o f t wa r e F L U ENT . Wt h d i f f e r e n t c o a l r a t i o s a n d o x y g e n r a t i o s , t h e v e l o c -
水泥回转窑煤粉燃烧过程的数值模拟
第36卷第7期 娃酸盐通报Vol.36 No.7 2017 年7 月_________________BULLETIN OF THE CHINESE CERAMIC SOCIETY_________________July,2017水泥回转窑煤粉燃烧过程的数值模拟赵良侠,边永欢,刘仁平,罗晓,王山辉,巧呈然(河北科技大学环境科学与工程学院,石家庄050018)摘要:为了研究水泥回转窑N O%浓度与温度的变化规律,对其内煤粉燃烧过程进行数值模拟,模拟结果与监测数据 相对吻合。
研究结果表明:当煤粉喷入量为一次风量的0.5〜1.0倍时,水泥回转窑冷却段、煅烧段、过渡段、出口 N O,浓度分别为〇、1〇〇〜1300 m g/m3、200 〜1300 m g/m3、100 〜1000 m g/m3,温度分别为1100 〜1300 K、1400 〜1900 K、1600〜1750 K、1300〜1700 K;煤粉喷射量为一次风量的0.6〜0. 8倍时,水泥回转窑内N O,浓度较低、温度较高,较为适合实际水泥生产工艺。
关键词:水泥;回转窑;数值模拟;NO;温度中图分类号:X511 文献标识码:A 文章编号:1001-1625 (2017) 07-2220-06Numerical Simulation of Pulverized Coal CombustionProcess in the Cement Rotary KilnZHAO Liang-xia ,BIAN Yong-huan,LIU Ren-ping, LUO Xiao, WANG Shan-hui ,DIA0 Cheng-ran(School of Environmental Science and Engineering, Hebei University of Science and Technology, Shijiazhuang 050018 , China)Abstract:In order to study the variation laws of N0X concentrations and temperatures in cement rotary kilns, combustion processes of pulverized coals were simulated inside them. Meanwhile, the results were relatively consistent with monitoring data. The study results indicate that the N0X concentrations are respectively 0, 100-1300 m g/m3, 200-1300 mg/m3, 100-1000 m g/m3and the temperatures are respectively 1100-1300 K,1400-1900 K,1600-1750 K,1300 〜1700 K in cooling zones, sintering zones, transition zones and outlet when pulverized coals are 0. 5-1.0times of primary airs. Whats more, the condition which pulverized coals are injected as 0. 6-0. 8 times of primary airs is relatively suited for actual cement technologies, because the N0X concentration is lower and the temperature is higher than other conditions in the cement rotary kiln.Key words :cement ; rotary kiln ; numerical simulation ;N O; temperature1引言现阶段国内以颗粒物、S02、N0p汞化物和臭氧为特征的区域性复合型大气污染问题较为突出[1],严重制约着社会经济的可持续发展、威胁人类健康。
链条炉煤粉燃烧的数值模拟
链条炉煤粉燃烧的数值模拟摘要:利用Fluent软件对国内某热电厂75t/h链条锅炉燃烧时炉膛内部煤粉燃烧和温度分布情况进行数值模拟。
首先提出煤的燃烧模型,选择k-ε模型对湍流进行模拟,分析炉膛热烟气流动以及温度分布,为SNCR设计提供理论支持。
关键词:煤粉燃烧;数值模拟;SNCR;温度1 前言王春波通过对某钢厂300MW锅炉燃烧过程进行数值模拟,对比发现掺烧高炉煤气时炉内温度水平有明显下降,且随着掺烧量的增加而加剧,掺烧高炉煤气产生烟气量增多,炉膛出口烟速明显增加,煤粉颗粒停留时间减小,煤粉很难燃烧完全[1]。
同时NO的生成量随高炉煤气掺烧量的增加而明显减少。
朱梅通过数值模拟的方法,设计了炉膛高速喷嘴的结构,并且通过模拟结果的比较,寻找到一种更适合退火炉的结构,使得炉膛内温度分布更加均匀[2]。
通过模拟的方法可以合理的调整喷口结构参数,对于喷口的设计更加方便,减少了大量的投资。
刘亮通过实际现场的改进,采用一种链条炉粉煤复合燃烧增容节能改造方案,在链条炉前后墙分别布置一次风二次风喷口将煤粉喷入燃烧。
通过链条炉炉排的燃烧与炉膛内增加煤粉的复合燃烧方式,将炉排上煤产生的大量热量来加热由煤粉喷口喷出的煤粉。
煤粉受热被点燃,释放出大量的热又反过来加热炉排上的煤块,加快煤块的着火速度,提高炉排上煤块的燃尽率。
并且通过改进链条炉煤粉的制粉系统的同时,改变了省煤器结构,增加其受热面积,进一步降低了有炉膛排出的烟气的热损失[3]。
张海军通过对同一种钝体喷口在不同的湍流模型下的模拟结果,结合 PIV 测量方法下的实验进行验证,通过比较发现在不同湍流模型下的模拟结果存在比较大的差别。
通过模拟比较,发现k-ε模型和 RNG k-ε模型模拟的结果与实际过程更加符合[4]。
徐榕在数值模拟的基础上分别使用 k-ε,RNG k-ε和代数应力(ASM)三种湍流模型进行模拟环形燃烧室中的三维流场分布图,通过模拟结果的比较,得出了三种湍流模型模拟的结果大致趋势是一致的。
煤粉燃烧过程的数值模拟
煤粉燃烧过程的数值模拟Ryoichi Kurose京都大学高级研究院流体科学与工程学院机械工程与科学系Hiroaki Wata nabe and Hisao Makino中央研究所的电力行业能源工程研究实验室摘要煤炭是一种能够满足电力进一步需求的重要能源资源,而且煤炭比其他化石燃料的储量丰富得多。
在燃煤发电厂,改善对环境污染物如NOx,SOx及包括未燃尽的碳粒在内的灰粒的含量的控制技术十分重要。
随着计算机性能的显着提高,人们强烈希望计算流体动力学(CFD)成为一种工具,成为一种研发和设计这种合适的煤粉燃烧的燃烧炉膛和燃烧器的工具。
这次审查的重点是突出我们的CFD 研究的最新进展,即煤粉在燃烧中的平均雷诺数纳维斯托克斯(RANS)的模拟和大涡模拟(LES)的最新进展,及未来的一些前景。
关键词:煤粉燃烧,数值模拟,平均雷诺数纳维斯托克斯模拟,大涡模拟1.介绍煤炭是一种能够满足电力进一步需求的重要能源资源,而且煤炭比其他化石燃料的储量丰富得多。
在燃煤发电厂,改善对环境污染物如NOx,SOx及包括未燃尽的碳粒在内的灰粒的含量的控制技术十分重要。
为了实现这些目标和要求,了解煤粉燃烧机理和先进的燃烧技术的发展十分必要。
然而,由于煤粉燃烧是一个非常复杂的现象,其中最高的火焰温度超过1500C,以及某些物质难以进行测量,如一些原子团种类和一些高活性固体颗粒,因此在燃烧过程中的煤粉燃烧机理至今没有得到很好的解释。
而且由于研发过程包含许多步骤,因此,新的燃烧炉膛和燃烧器的发展需要较高的成本和较长的时间。
随着计算机性能的显著提升,煤粉燃烧领域的计算流体动力学正在被研发。
在这种方法中,电脑解决了燃烧领域的控制方程式,这使它能够提供温度和化学物质种类分布的详细信息和在整个燃烧空间中煤粉颗粒的行为,而上述那些通过实验是不能得到的。
此外,此种方法有助于在相对较低的成本条件下重复审查任意条件下的煤粉的流场和各种参数。
因此,强烈地希望计算流体动力学(CFD)能够成为燃烧炉炉膛和燃烧器研发和设计的一种工具。
高炉三维气固湍流和煤粉燃烧过程数值模拟
立并发展 了高炉风 1回旋区湍流气 固两相 流动 和煤粉燃烧 的三 维数 学模型 。用所建模 型分别对冷 : 7
态模型 内气 固两相流动和某企业 70m 5 3高炉风 1回旋 区 内的 气固两相 三维流动与 煤粉 燃烧进行 了 : 7 数值模拟。采 用三维激光相位 多普勒分析仪( D ) P A 对冷 态模 型 内气固两相流 场进行 了测量 , 实验 结 果与冷态两相 流动的模拟结果基本一致。热 态模拟结果给 出了气相温度和组 分浓度分布 , 模拟结果 与实验测量结果较吻合 . 示 了风 口回旋 区内气 固两相流动和煤粉燃烧的基本性质和特 点。 揭 关 键词 : 高炉风 口回旋 区 ; 流 ; 固两相 流 ; 粉 燃烧 ; 湍 气 煤 数值模 拟 中图分类 号 : F5 8 T 3 文章编 号 :0 5— 8 0 2 0 ) 2— 1 3— 5 10 9 3 ( 0 7 0 0 6 0
aye)m aue ns e ea are n a e e e ntecl l e sl n s d t.T e l r esrmet.G n rl gemet s d t e a ua dr ut adt t aa h z im bw h c t e s e
n me ia i lt n r s l r n a c r a c t x e me t a u e n s T e su y r v a s u rc smu ai e u t ae i c o d n e wi e p r n a me s r me t . h t d e e l l o s h i l
Nu e ia i u a i n fTh e - i e so a a - a tce m rc lS m l t so r e d m n in lG sp ril o Tu b ln o a d Pu v rz d Co lCo b si n i a tF n c r u e tFlw n l e ie a m u to n Bl s ur a e
煤粉锅炉富氧燃烧的数值模拟研究
煤粉锅炉富氧燃烧的数值模拟研究摘要:本文旨在通过数值模拟方法,来研究煤粉锅炉富氧燃烧技术。
本文采用CFD软件ANSYS FLUENT计算了不同条件下煤粉燃烧的温度、火花流体机理、热质传输特性等热力学参数。
通过分析模拟结果发现,富氧燃烧技术可以降低出口温度和增加热效率,有效改善煤粉锅炉的工况,增强其运行安全性。
关键词:煤粉锅炉;富氧燃烧;CFD;数值模拟;温度;热效率正文:煤粉锅炉是一种多用途的重要能源装置,具有经济性、安全性、灵活性和可操作性等优点,因此广泛应用于能源、化工及其他行业。
然而,煤粉锅炉的存在也伴随着一些问题,尤其是因过热、超温和污染的问题而带来的火灾和运行不安全的风险。
为了解决这些问题,研究者提出了富氧燃烧技术,以改善锅炉运行工况。
本文旨在使用CFD软件ANSYS FLUENT,通过数值模拟研究煤粉锅炉富氧燃烧技术。
我们首先采用RNG k-ε模型,建立定常空气-煤粉层燃烧模型,并采用双态壁界条件计算温度控制。
然后,通过不同的氧气比例高低,模拟火花流动模式和热质传输特性的变化。
最后,我们分析了结果数据,发现富氧燃烧技术可以降低出口温度和增加热效率,有效改善煤粉锅炉的工况,增强其运行安全性。
本研究的结果可以为以后的工程实践提供参考,帮助改善煤粉锅炉的工况,增强其运行安全性。
本文的研究也应当引起我们的重视,为了更好地实现富氧燃烧技术的运用,需要对对流传热、扩散热及燃烧效率等有影响因素进行详细分析,这[]将会为改善煤粉锅炉工况和节能提出新的思路。
此外,结构参数和操作参数等都会对煤粉锅炉富氧燃烧技术的效果产生不同程度的影响,今后可以以此为依据,进行更加详细的研究。
同时,借助于现代计算技术和模拟工具,如果能够进一步提高模拟精度,就可以及时根据实际情况,准确地指导实际工况的运行。
值得注意的是,本文的研究只采用了单一煤粉锅炉富氧燃烧技术,而在实际应用中,存在多种煤粉组合,必须根据具体组合来进行研究。
煤自燃过程的实验及数值模拟研究
煤自燃过程的实验及数值模拟研究
煤层自燃严重影响着煤炭工业发展,给矿井生产带来极大安全隐患。
由于实际条件下的煤自燃过程很难描述清楚,使得煤层自然发火预测预报技术的发展受到严重制约。
本文根据现场实际条件,抓住煤体氧化放热和环境散热这对主要矛盾,根据多孔介质渗流力学和传热传质学理论,建立了煤自燃过程数学模型。
依据模型确定出煤自燃过程数值模拟所需的关键参数为煤的耗氧速度和放热强度及松
散煤体内的氧气扩散系数和渗透系数,通过实验对其进行了研究和测试。
提出了实际条件下煤体放热强度和耗氧速度的计算方法,并通过实验得出了煤体粒度影响函数。
设计建造了国内最大的(装煤量15t)煤自然发火实验台,首次实现了煤自然发火全过程的实验模拟,掌握了煤自燃高温点的发生、发展、变化过程及停止供氧后高温点的降温规律。
根据实验台条件,建立了数学模型,通过数值模拟,解决了煤自然发火实验模拟条件单一的问题。
针对现场应用,提出了实际条件下漏风强度的测算方法及煤自燃危险区域的判定准则,通过煤自燃过程的数值模拟研究,发展了综放面采空区及巷道自然发火预测理论。
采用实验测定的关键参数,结合现场可测参数,并考虑这些参数的实际变化情况,运用理论与实践相结合的方法,通过数值模拟再现了综放面采空区和巷道松散煤体自燃的发生及发展过程,首次解决了实际条件下煤自然发火条件、地点和时间的判定及预测问题。
该理论成果经受了现场多次煤层自燃火灾预测预报和防灭火工作的检验,在该成果的指导下,成功地预防和扑灭了几十次矿井煤层火灾。
煤粉燃烧过程中NOx排放的数值模拟资料精
万千瓦 增加了 3000 万千瓦 按照国务院的要求 “十五”后三年 每年要新增
装机 2500 万千瓦 其中约 60 ~70 1500 — 1750 万千瓦 以上为火电机组
煤粉的燃烧将会产生大量的粉尘 SO2 和氮氧化物等有害污染物 其中燃煤 排放的 NOx 占总氮氧化物排放的 67 长期以来 我国对 NOx 污染问题的重视
旋风炉
1000~1500mg/N m3
燃油炉
600~1400mg/N m3
2002 年 1 月 30 日国务院第 54 次常务会议通过了新的收费制度 排污费 征收使用管理条例 2003 年 7 月 1 日起连同 排污费征收标准管理办法 简 称 收费标准 等配套规章一起施行 收费标准 的实行对火电行业形成了较 大的压力 排污收费标准规定 每排放 1kg 的 NOx 收费 0.63 元 上海市电站总 装机容量 9612MW 大部分为燃煤机组 燃煤锅炉造成了严重的环境污染 如 1998
II
上海交通大学硕士学位论文
NUMERICAL SIMULATION ON NOx EMISSIONS IN COAL COMBUSTION PROCESS
ABSTRACT
Pollutants from coal combustion are a major environmental problem since they produce large amount of dust, sulfur dioxide and nitrogen oxides. The formation of nitrogen oxide (NOx) in combustion system is a significant pollutant in the environment, and the control of NOx emissions is a world-wide concern as the utilization of fossil fuel continues to increase. A lot of researches on the control of nitrogen oxides has been done, however, so far many real flames still cannot be described and calculated precisely due to the complex interaction between reaction and diffusion occurring in coal combustion. So it is always important to develop the theoretical and experimental studies for the combustion.
Fluent软件模拟计算煤粉燃烧的机理及其模型实现的方式
2010 年第 11 期
丁历威, 等: Fluent 软件模拟计算煤粉燃烧的机理及其模型实现的方式
33
跟踪颗粒数=1×(Number of Tries) (1) 2.2 Group 类型
Group 类型是指煤粉颗粒按照直径分成多股 煤粉流, 每股煤粉流之间直径、质量、 流量、 位置 都各不同, 然后每股煤粉流按照 1 个 Single 类型 来处理, 所以该类型下 Fluent 跟踪的颗粒数为: 跟踪颗粒数=(Number of Group)×(Number of Tries)
(2) 2.3 Surface 类型
Surface 类型是指煤粉颗粒根据入口面网格分 成多股煤粉流。 每 1 个面网格相当于 1 个 Group, 所以该类型 Fluent 跟踪的颗粒数为:
跟踪颗粒数=(Number of Face)×(Number of Group)×(Number of Tries) (3)
率很慢, 扩散过来的氧量远远超过反应所需的氧 量, 这时燃烧速率与反应有关。
(2)当温度超过 1 400℃时, 由于化学反应速 率常数随着温度升高急剧增大, 致使反应所消耗 的氧量远远超过扩散过来的氧量, 这时燃烧速率 主要由氧气的扩散速度决定。
(3)当温度在 1 000~1 400℃之间 , 燃烧速率 受到反应速率和氧气扩散速率双重影响。 1.4 灰分的吸热或者冷却(Fourth Law)
烧模型实现的方式。在详细说明煤粉单个颗粒燃烧每一步过程的含义和切换的条件后, 通过自编用户自
定义函数(UDF)深入研究了煤粉燃烧模型中颗粒跟踪数计算的方法、 过程之间的切换、 过程与规 则之
间的关系、 颗粒生命周期内最多调用自定义规则的数目以及过程规则的调用机理, 并用示意图的方式
煤粉颗粒群着火和燃烧过程的数值模拟
煤粉颗粒群着火和燃烧过程的数值模拟赵云华, 何玉荣, 陆慧林, 刘文铁, 沈志恒, 叶校圳(哈尔滨工业大学能源科学与工程学院, 哈尔滨150001)摘要: 建立了一维非稳态球形煤粉颗粒团的群燃烧模型.数值模拟煤粉颗粒团的着火和燃烧过程,获得了颗粒团燃烧火焰随时间的变迁.分析了煤粉颗粒团内部参数和外部环境参数对颗粒团着火和燃烧的影响.随着颗粒团内煤粉浓度的增加,颗粒团的均相着火延迟先减小后增加.增加煤粉颗粒尺寸和降低外部温度都会明显延迟均相着火.环境氧气含量的增加会减小着火延迟,同时增加颗粒团的燃烧速率.模拟计算和文献试验结果的变化趋势相吻合.关键词: 颗粒群燃烧; 煤粉颗粒团; 着火中图分类号: TK16 文献标志码: A 文章编号: 100628740 (2007) 022*******N um erical Analysis of I gnition and Combustion ofPulverized Coal P article C loudZH AO Y un2hua , HE Y u2rong , LU Hui2lin , L I U Wen2tie , SHE N Zhi2heng , YE X iao2zhen ( S chool of Energ y S cience and Engineering , Har b in Institu te of Techn olog y , Harb in 150001 , C hina)Abstract : A transient g rou p combustion m od el w as presented for the p u lverized coal particles in sp herical cloud. H om og eneous and heterog eneous reactions w ere m o d eled to take account of the interactions o f particles. The variation of flame stru cture of cloud was o b tained. The effects o f param eters on cloud ig nition and combustion w ere analyzed. The hom og eneous ig nition d elay w as de2 creased , an d then increased w ith the increase of concentration of par ticles. The increase of p u lverized coal size and d ecrease of amb ient g as temperature d elay hom og eneous ig nition. The d elay time w as redu ced , and the combustion rate w as increased w ith the increase of oxyg en concentration. S imulation resu lts w ere in ag reement w ith p ub lished experimental find ings.K ey w ords : g rou p combustion ; p u lverized coal cloud ; ig nition煤粉颗粒的着火和燃烧过程对于煤粉锅炉燃烧效率、排放特性和运行控制等具有重要影响,成为煤燃烧领域的重要研究方向之一. 煤粉颗粒的加热着火和燃烧与煤质特性和结构以及加热条件等密切相关. 单煤粉颗粒的加热、挥发分析出和燃烧过程中煤粉颗粒间反应过程的相互作用对于着火和燃烧行为的影响可以忽略不计.然而,随着煤粉燃烧技术的发展( 如浓淡煤粉燃烧技术、超细煤粉再燃技术等) ,由于气体湍流和颗粒碰撞、颗粒表面力作用等内在因素和外部因素(如煤粉浓缩器) 等,使得在气流中出现局部高浓度煤粉颗粒,形成低颗粒浓度的稀相区和高颗粒浓度的密相区.在富集煤粉颗粒的密相区,煤粉颗粒之间物质扩散、化学反应等相互作用将直接影响煤粉颗粒的着火、燃烧、燃尽和排放特性 1 ,因此需要考虑高数密度煤粉颗粒群燃烧特性的影响. Annam alai 等 2 建立了稳态的焦炭颗粒团燃烧模型,分析了颗粒团浓度对着火过程的影响. Du 等 3 建立了一维圆柱形煤粉颗粒群燃烧模型, 考虑了煤颗粒相的热解挥发、焦炭异相燃烧和气相燃烧,模拟分析了颗粒群的非稳态均相着火和异相着火.姚强等 4 以单煤粉颗粒为基准,采用有限气体空间的假设建立颗粒群着火模型,研究辐射传热条件下煤粉颗粒群非稳态着火过程,预测了不同条件下的着火方式、着火时间和着火温度. 盛昌栋等 5 采用了假设颗粒团内部具有相同气体温度和颗粒温度的集总参数法,收稿日期: 20062 032 22 .作者简介: 赵云华(1981 —) ,男,博士研究生,zha oyunhua @gmail . c om.·108 · 燃 烧 科 学 与 技 术 第 13 卷 第 2 期研究了受辐射加热条件下煤粉浓度对颗粒群 着 火 方 式 、着火时间和着火温度的影响 .笔者以煤粉颗粒团作为研究对象 ,考虑颗粒团内部气体和固相温度的差异 ,以及煤粉颗粒挥发分析出和焦炭燃烧反应过程 ,建立一维非稳态球形颗粒团燃 烧模型 ,模拟煤粉颗粒团的着火和燃烧过程 ,分析了煤 粉浓度 、环境温度和煤粉颗粒直径对煤粉颗粒团着火 和燃烧的影响 .的挥发分析出率 m ´V 和焦炭的异相反应消耗率 m ´C 共同决定 ,即m ´P = m ´V + m ´C煤粉颗粒温度 T p 由热平衡计算得到 ,即(4)dT p = m ´V q V + m ´C q C + hA p ( T g - T p ) +(5)m p c p d t 4 4ζp σA p ( T w - T p ) 式中 :方程右边第 1 、第 2 项分别为热解挥发和焦炭反 应放热 , q V 为挥发吸热率 , q C 为焦炭燃烧放热率 ; 右 边第 3 项为对流换热 ,换热系数 h = N u λ/ d p ,λ为气相 导热系数 , A p 指煤粉颗粒表面积 ; 右边最后一项为辐 射换热 ,ζp 为煤粉黑度 ,σ为波尔兹曼常数 , T w 为辐射 边界温度 .假设煤粉颗粒燃烧过程中挥发分析出只影响煤粉 颗粒的密度 ,直径保持不变 ;煤粉的异相反应只影响煤煤粉颗粒团的群燃烧模型1 考虑一球形颗粒团的半径为 R C ,其中均匀分布着直径为 d p 、数密度为 n 的煤粉颗粒. 整个颗粒团处于 高温环境中 ,初始时刻颗粒团内的颗粒和气体温度均 为 T 0 ,环境温度为 T ∞ . 在对流传热和环境热辐射共同 作用下 ,分 ;挥发分与周围空气混合并发生氧化反应. 学反应速率达到一定程度 ,便产生着火燃烧 ,团燃烧火焰 .为简化分析 ,现假设如下 : ( 1) 颗粒团为球形对称结构 ; (2) 煤粉颗粒为点源 , 燃烧过 程 空 间 位 (3) 理想气体 ,全场等压 、无黏性 . 对于一维燃烧和传热 过程 ,煤粉颗粒群燃烧过程可用守恒方程表示.111 气相守恒方程质量守恒方程4 ,6 ] ,即(6) (7)[ 3 ]挥发分析出率采用两步竞争反应模型 ,即471336 ×105 m ´V , I = - αI 317 ×10 ex p - m d a f RT(8) 521511 ×1013m ´V , II = - αII 1146 ×10 ex p - m d a fRTm ´V , I m ´V , II m ´daf = + (9)5ερ 1 5 αIαII2ερ ) + 2 ( r u = S m r 5 r( ) 15 t 式中 :αI 和 αII 为热解常数 ; m da f 为尚未分解的原煤质量.物质组分守恒方程5 (ερY i ) 1 5 1 5 2ερ ) 2114 焦炭异相反应模型在模拟计算中 ,煤粉颗粒的异相反应过程为C + O 2 →C O 22C + O 2 →2C O C + C O 2 →2C O+ 2 ( r uY i r 5 r = 2 r 5 r r ερD · 5 t 5 Y i(2)(A ) (B ) ( C )+ S i5 r能量守恒方程 5 (ερh T ) 1 5 12ερ + ( r uh T= · r 2 5 r r 2) 化学反应速率按表 16 —8 计算. 各物质的总反应速率为u A u B u C 5 2ερD 5 h T+ + πd 2ρ m ´ = - (10) ( ) + S 3r C β β β p ph T 5 r 5 rCA CB CC2m ´O ,s = - ( u A + u B )πd p ρp(11)式中 :ρ为气相密度 ;ε为颗粒团空隙率 ; u 为混合气 体径向速度 ; r 为距颗粒团中心的距离 ; D 为气体扩散系数 ; Y i 为组份 i 质量分数 ; h T 为混合气体热焓 ; S m 、21 +βCA2m ´CO ,s = -u A - u C πd p ρp (12)βCA2其中 ,βCA 、βCB 和 βCC 为化学反应当量比 ,计算反应速率S i 和 S h 分别为煤粉颗粒相引起的气相质量 、物质组 T[ 6 ]所用的颗粒表面的气体浓度 Y O ,s 和 Y CO ,s 分别为 分和能量的增加率. 112 固相控制方程煤粉颗粒燃烧过程中质量变化率 m ´P 由热解产生22m ´O , s m ´O ,s m ´ c p p 2 2(13)Y O ,s = - + ex p Y O + m ´pπd pN u λpm ´p2 22007 年 4 月 赵云华等 :煤粉颗粒群着火和燃烧过程的数值模拟 ·109 ·3m ´CO , s响 ,本文采用 Du 等 的着火定义 : 局部气相温度超过 两边相邻层气体的温度 ,则认为发生均相着火 ;局部的 颗粒相温度超出它周围气相的温度 ,则认为发生异相 着火. 图 1 为煤粉团燃烧模型.m ´CO ,sm ´ c 2p pY CO ,s = - + 2ex p (14)Y + m ´p πd N uλ COm ´22p p p 其中 , c p 和λp 分别为煤粉颗粒的比热容和导热系数.115 气相反应假设煤粉颗粒热解挥发物质为 CH 4 和 C 3 H 8 ,气相 化学反应为2CH 4 + 3O 2 →2C O + 4 H 2 O 2C 3 H 8 + 7O 2 →6C O + 8 H 2 O 2C O + O 2 →2C O 2 2C O 2 →2C O + O 2其中 ,化学反应速率见表 1 所示.(D )( E ) ( F )( G )表 1 气相化学反应速率计算模型图 1 煤粉团燃烧模型2 图 2 表示颗粒团燃烧过程中气相温度的变化 . 由 图可见 ,颗粒团内部气相温度在前 16 m s 平缓升高 ,约在 1715 m s 时出现局部峰值 ,峰值点位置位于颗粒团 外 . 随后 ,气相温度逐渐升高并形成燃烧火焰 ,在 18 m s 之后 ,气相燃烧火焰分离成两个 ,一个穿透到颗粒团内 部 ,另一个则向外围传播 . 最终 ,内部火焰消失 ,仅有外 部火焰维持并控制着颗粒团的燃烧 .22 CH O g4 2 C H O g 3 8 2 H O O g2 2 2116 初始 、边界条件和计算方法在 t = 0 时刻气相 Y i = Y i ,0 , T g = T g ,0 , r ≤R CY i = Y i , ∞ , T g = T ∞ , r > R Cm p = m p ,0 ,ρp =ρp ,0 , d p = d p ,0 ,T p = T p ,0 , T w = T ∞r →∞: Y i = Y i , ∞ , T = T ∞ , T w = T ∞煤粉颗粒相 边界条件r = 0 : 5 T = 5 Y = u = 0 5 r 5 r g k图 2 气相温度的变化本文采用 Du 等 3 数据作为基本工况参数 : T =∞ 1 500 K , T g ,0 = T p ,0 = 300 K , Y O ,0 = 0123 , R C = 012252图 3 表示初始的单个火焰结构 . 由图可见 ,火焰的 位置在颗粒团表面附近 ,是由热解挥发物质和氧气的 剧烈反应诱发的 ,火焰的形成将快速加热火焰中的煤粉颗粒 ,从而产生高挥发分析出率 . 图 4 为双峰火焰结 构 . 由于初始火焰中氧气的耗尽和挥发分的大量析出 , 单个火焰最终被分离成两个火焰峰 . 内部火焰向煤粉 颗粒团内传播 ,一方面消耗颗粒团内部的氧气 ,另一方 面高温火焰快速加热内部煤粉颗粒 ; 而外部火焰则向 外界有氧的环境中扩散 . 图 5 为挥发分扩散控制的燃 烧火焰结构 . 此时 , 内部火焰因为 氧 气 的 耗 尽 已 经 消 失 ,整个颗粒团的燃烧主要表现为挥发分的扩散燃烧 .- 5cm , d p = 64μm ,ρD = 5 ×10 kg/ (m ·s) ;煤质成分 : w C =84163 % , w H = 5147 % , w N = 2113 % , w O = 7177 % ; 煤粉颗粒数密度 n = 61417 ×103cm - 3( 煤粉空气质量比为31) ,颗粒密度 ρp = 1 300 kg/ m . 方程 (1) ~(3) 采用隐式 上风格式离散求解 ,而源项的求解采用显式方法.2 计算结果和讨论211 颗粒团燃烧过程煤粉着火是一个复杂的过程 ,受到许多因素的影反应过程化学反应速率/ (m ol ·(cm 3 ·s ) - 1 )A B C D E F Gu A = 11225 ×103Y O ,s exp ( - 91977 ×104/ RT p ) u B = 11813 ×103Y O ,s exp ( - 11089 ×105/ RT p ) u C = 71351 ×103Y O ,s exp ( - 1138 ×105/ RT p ) u D = 219 ×1012 Y 017 Y018 exp ( - 24 358/ T ) u E = 816 ×1011 Y011 Y1165 exp ( - 15 098/ T )u F = 21239 ×1012 Y C OY 015 Y 0125 exp ( - 20 131/ T )u G = 510 ×108Y CO exp ( - 20 131/ T g )·110 ·燃烧科学与技术第13 卷第2 期图6 为剩余焦炭燃烧的火焰结构. 此时,热解挥发物质燃烧完毕,内部焦炭颗粒被氧化成C O ,颗粒团外形成C O 的扩散燃烧火焰.212 煤粉浓度的影响图7 表示不同煤粉浓度颗粒团的着火时间曲线.当煤粉浓度较低时,均相着火时间随着浓度的增加而减小;而当煤粉浓度较高时,均相着火时间随浓度的增加而增大. 因为在低浓度时,颗粒浓度的增加会产生更多的挥发分,从而更容易引发均相着火;但是当颗粒浓度增加到一定程度,固体颗粒团的温度升高将会吸收更多的热量, 从而延迟挥发分的析出和均相着火. Du等3进行了煤粉颗粒群的着火实验研究,得到了煤粉颗粒群的均相着火延迟时间.本文模拟获得的均相着火延迟时间与Du 等的实验结果具有相同的变化规律,两者定量上的差异主要来自于Du 等的圆柱模型与本文球型模型的不同以及实验中颗粒群与气流间存在相对速度等因素. 比较图7 中颗粒团均相和异相着火时间,可以看出,对于颗粒/ 空气质量比小于5 的稀疏煤粉颗粒团先发生异相着火,而密实的煤粉颗粒团则先发生均相着火. 盛昌栋等9 实验研究了煤粉气流着火方式与煤粉浓度的关系,结果也表明随着煤粉浓度的增加,着火方式由异相着火向均相着火过渡. 图7中,当过渡到颗粒团先发生均相着火后,异相着火时间有较大突变. 因为均相着火发生后,颗粒团周围会形成挥发分燃烧火焰,导致局部气体温升大于颗粒异相反应造成的煤粉颗粒温升,从而使得异相着火时间有较大延迟.图3初始火焰结构的物质组分和温度分布图4双峰火焰结构的物质组分和温度分布图5 挥发分火焰结构的物质组分和温度分布图7 煤粉浓度对着火延迟的影响213 颗粒尺寸的影响图8 表示相同煤粉颗粒浓度下不同煤粉粒径对煤粉团燃烧的影响. 由图可见,均相着火时间随着内部颗粒直径的减小而急剧减小,因为随着颗粒直径减小,颗粒团内的比表面积将增加,辐射传热和对流换热也随之加快,颗粒温度迅速升高并释放出挥发分,从而更快达到均相着火. 图中整个颗粒团的质量变化曲线基本相似,其变化可大致划分为以热解挥发为主的前一段和之后的焦炭燃烧段. Ayling 和Sm ith10 实验研究直径为45 μm 和25μm 的焦炭粒子流在1 415 K 的空气环图6 焦炭火焰结构的物质组分和温度分布2007 年 4 月 赵云华等 :煤粉颗粒群着火和燃烧过程的数值模拟 ·111 ·果 . 由图可见 ,两者的变化趋势是一致的 . 图 11 表示相同煤粉浓度条件下 ,氧气质量分数对颗粒团燃烧过程 的影响. 由图可知 ,随着氧气质量分数的增加 ,整个颗 粒团的热解挥发和焦炭燃烧都有明显加快.境中的燃烧速率 ,结果表明焦炭颗粒群的燃烧速率基 本相同. 图中计算结果也表明 ,整个颗粒团燃烧速率几乎不受内部煤粉颗粒尺寸的影响 .图 8 煤粉颗粒尺寸对着火和燃烧的影响图 10 氧气质量分数对着火的影响214 环境温度的影响图 9 表示相同煤粉浓度下环境温度对颗粒团燃烧的影响. 由图可见 ,随着环境温度的增加均相着火明显 提前. 但是整个颗粒团的燃烧速率并不随着外界温度 的变化而变化. Annam alai 2研究了外界温度对焦炭颗 粒团燃烧的影响 ,发现颗粒团的燃烧速率随着外界温度的增加而增加 ,同时作者指出了气流对颗粒团的吹 散作用 、温度对氧气扩散的作用等物理效应是促使燃 烧速率增加的主要原因 . 然而 ,在本文的模拟研究中 , 煤粉颗粒被固定 ,ρD 也被处理为常数 ,因此 ,模拟计算 结果 未 能 反 映 环 境 温 度 变 化 对 颗 粒 团 燃 烧 速 率 的 影响.图 11 氧气质量分数对群燃烧的影响3 结 语本文建立了一维 非 稳 态 煤 粉 颗 粒 团 的 群 燃 烧 模 型 ,数值模拟了煤粉颗粒团的燃烧过程 ,获得了着火的产生和整体火焰的变化 . 研究表明 ,当颗粒和空气质量 比小于 5 时 ,煤粉颗粒团的着火先是异相着火 ,而当质 量比大于 5 时 ,煤粉颗粒团的着火先是均相着火 . 随着 燃烧过程的进行 ,初始的燃烧火焰被分离 ,一个向颗粒 团内部传播 ,另一个向外传播 ,最终内部火焰消失 ,外 部火焰控制着整个颗粒团的燃烧 .随着煤粉团的颗粒浓度增加 ,煤粉颗粒团的均相 着火延迟先减小后增加 . 增加煤粉颗粒尺寸和降低外 部温度会明显延迟均相着火 ,但对颗粒团的燃烧速率 影响不大 . 环境氧气浓度的增加会减小着火延迟 ,同时 明显的加快颗粒团的燃烧速率. 图 9 外界温度对着火和燃烧的影响215 氧气质量分数的影响图 10 表示相同煤粉浓度下氧气质量分数对着火 时间的影响 . 由图可见 ,随着氧气质量分数的增加 ,均 相着火和异相着火都相应的会提前 . 但是 ,氧气质量分 数的增加对着火的影响幅度没有颗粒团内部颗粒直径 和外部温度的影响明显. 图 10 中同时给出 Saito 等11 对 45 μm 煤粉颗粒流在 1 200 ℃环境中的着火实验结参考文献 :1韩才元 ,徐明厚 ,周怀春 ,等 . 煤粉燃烧 M . 北京 : 科学出·112 · 燃 烧 科 学 与 技 术第 13 卷 第 2 期Fuel ,2003 ,82 :8932907.版社 ,2001.Han C aiyu an , Xu Ming hou , Zhou jet J . Hu aichun , et al . 7Westbrook C K , Dryer F L . C hemical kinetic m od eling of hydro 2 Pulver ized C oal C ombustion M . Beijing : S cience Press ,2001 (in Chinese ) .Annamalai K , R amaling am S C. G ro u p combustion of char/ car 2 Prog En ergy C o mbust Sci , 1984 , 10 : 12carb on combustion J . 57.2 b on particlesJ .C ombustion and Flame ,1987 ,70 :3072332.8Pre P , Hemati M , Marchand B. S tud y on natural g as combus 2 tion in flu id ized b eds :M o d eling and experimental valid ation J . Chemical E ngineering Science ,1998 ,53 :287122883.盛昌栋 ,齐 宏 ,徐明厚 ,等 . 煤粉气流着火方式与煤粉浓度的关系 J . 电站系统工程 ,1995 ,11 (3) :31237.Sheng C hang d ong , Qi H ong , Xu Ming hou , et al . Relationsh ip b etw een the ig nition m o d es and the coal concentrations in coal cloudsJ . P ower S ystem Engine ering , 1995 , 11 ( 3) : 31237 ( i n C hinese ) .3 Du X Y , G opalakrishnan C , Annam alai K. Ig nition and com 2 Fuel ,1995 ,74 (4) :4872494.bustion of coal particle streamJ .4姚 强 ,周俊虎 ,涂建华 ,等 . 煤粉群非稳态统一着火模型J . 浙江大学学报 ,1996 ,30 (4) :4462454.Y ao Qiang , Zhou J unhu , Tu Jianhu a , et al . An uns tead y un i 2 versal ig nition m od el of the p u lverized coal g ro u p J . J o urnal o f Zhejiang University ,1996 ,30 (4) :4462454 (in C hinese ) .盛昌栋 ,袁建伟 ,徐明厚 ,等 . 受辐射加热的煤粉颗粒群着火模型 J . 燃烧科学与技术 ,1996 ,2 (1) :38245.9 510 Ayling A B , S mith I W. Measured temperatures of burning p u l 2verized 2f u el particles and the nature of primary reaction produ ct Sheng C hangd ong , Yu an Jianw ei , Xu Ming hou , et al . m od el of p u lverized 2coal cloud heated b y rad iation J . Ig nition J o urnalC ombustion and Flame ,1972 ,18 :1732184.J .o f C o mbustion Science and Technology , 1996 , 2 ( 1 ) : 38245 ( in C hinese ) .G u o Y C , C han C K , Lau K S. Nu merical stud ies of p u lverized coal combustion in a tubu lar coal combustor w ith slanted oxyg en11 S aito M , S ad akata M , S ato M , et al . C ombustion rates of p u l 2verized coal particles in hig h 2temperatu re/ hig h 2oxyg en co ncen 2 C ombustion and Flame ,1991 ,87 :1212.6tration atm osp here J .。
煤粉颗粒团燃烧过程的数值模拟
第19卷 第3期2010年9月 矿 冶M IN ING &M ETALLURGYVol 119,No 13S epte m ber 2010文章编号:1005-7854(2010)03-0045-04煤粉颗粒团燃烧过程的数值模拟刘向军,林超(北京科技大学热能系,北京100083)摘 要:煤粉的燃烧在工程实际中很多情况下都是聚团燃烧,颗粒内部传热传质机理复杂。
本文根据多孔介质的传热传质学理论,引入对流项,建立了气流场中煤粉团燃烧的三维数学模型,并对不同孔隙率的煤粉团的燃烧过程进行了对比计算,计算结果合理,揭示了孔隙率对煤粉团着火燃烧的影响,为深入研究煤粉团燃烧特性打下基础。
关键词:煤粉颗粒团;燃烧;数值模拟中图分类号:TK16 文献标识码:ANUM ER I CAL SI M ULATION OF THE COAL CLUSTER CO M BUSTION PROCESSLIU X iang-j u n,LIN Chao(Ther m al Engineering D e part m ent ,University o f Science&Technology B eijing,Beijing 100083,Ch i n a)ABSTRACT :M ost o f the coa-l gas co m bustion process acti n g i n a dense partic l e -gas t w o -phase fl o w fiel d is coa l cluster co m bustion,w h ic h is d ifferent fr o m tha t o f a sing le coa-l partic l e co m bustion because of the co m plex inner -cluster heat and m ass transfer 1I n this paper ,a 3D m athe m atica lm ode lo f a coa l c l u ster co m bustion i n hot air fl o w ,w hich based on the heat and m ass conservation la w of porous m edia ,is established 1The governing equations i n -clude heat transfer equation and gas transfer equati o n 1The detailed heating -up and co m busti o n processes of a coa l cluster are obtained and t h e effects of the porosity are also ana l y sed 1A ll t h ese results are reasonab le and o ffer theo -retical bases for futher st u dy i n g t h e co mbusti o n pr ocess of coa-l c l u sters 1KEY W ORD S :coal cluster ;co m bustion ;nu m erical si m ulati o n 收稿日期:2010-04-28作者简介:刘向军,博士,教授,长期从事流动与燃烧过程的基础研究。
水平浓淡煤粉燃烧器内气固两相流的数值模拟
水平浓淡煤粉燃烧器内气固两相流的数值模拟水平浓淡煤粉燃烧器内气固两相流的数值模拟水平浓淡煤粉燃烧器内气固两相流的数值模拟是一个涉及流体力学、热力学和化学反应等多个学科的复杂课题。
它涉及到了流体在燃烧过程中的传热、传质以及相变等多个物理过程,而且在燃烧过程中,气固两相流的相互作用更是复杂多变。
其数值模拟需要考虑到多种因素,如湍流模型、燃烧模型以及颗粒运动模型等。
在这篇文章中,我将从基础概念开始,逐步深入,探讨水平浓淡煤粉燃烧器内气固两相流的数值模拟。
1. 什么是水平浓淡煤粉燃烧器内气固两相流?水平浓淡煤粉燃烧器内气固两相流是指在燃烧器内,煤粉和空气以及燃烧产物之间同时存在的流动状态。
其中,煤粉和空气混合形成气固两相流,而在燃烧过程中,燃烧产物也会与气体形成两相流动。
在水平浓淡煤粉燃烧器内,气固两相流的流动状态复杂多变,既有湍流现象,又有颗粒间的相互作用,因此需要进行数值模拟来更好地理解和控制这一过程。
2. 数值模拟的基本原理数值模拟是利用计算机对实际物理过程进行数值求解,以获得系统的流动信息、温度分布、物质转移等相关数据。
在水平浓淡煤粉燃烧器内气固两相流的数值模拟中,需要考虑到多个物理过程,例如流体的密度、粘度、热传导系数等,颗粒的运动状态、燃烧产物的生成和传输等。
数值模拟需要建立相应的数学模型,以描述和求解这些物理过程。
3. 湍流模型在数值模拟中的应用在水平浓淡煤粉燃烧器内,流体的运动状态往往处于湍流状态,因此需要采用湍流模型来描述流体的湍流运动。
常用的湍流模型包括k-ε模型和RANS模型等,在数值模拟中,选择合适的湍流模型对于准确描述气固两相流的运动状态至关重要。
4. 燃烧模型在数值模拟中的应用燃烧是水平浓淡煤粉燃烧器内最为重要的物理过程之一,燃烧模型的选取直接影响到数值模拟的准确性。
常用的燃烧模型包括简化化学反应模型、进一步细化的半简化模型以及详细化学动力学模型等。
在数值模拟中,需要选择适合燃烧过程特点的模型,以准确预测燃烧产物的生成和传输过程。
水平浓淡煤粉燃烧器内气固两相流的数值模拟
水平浓淡煤粉燃烧器内气固两相流的数值模拟在燃烧工程领域,水平浓淡煤粉燃烧器内气固两相流的数值模拟一直是一个备受关注的研究课题。
煤粉燃烧作为一种重要的能源利用方式,其优化设计和运行参数对于提高燃烧效率、减少污染排放具有重要意义。
水平浓淡煤粉燃烧器内气固两相流的数值模拟可以通过计算机模拟对流场、燃烧特性和热力学参数进行分析和预测,为煤粉燃烧工程的设计和优化提供理论依据和技术支持。
1. 水平浓淡煤粉燃烧器内气固两相流的数值模拟概述水平浓淡煤粉燃烧器是煤粉燃烧系统中的重要部件,其内部气固两相流动态特性对于燃烧效率和环境排放具有重要影响。
数值模拟是一种有效的研究方法,通过建立数学模型和求解数学方程来描述和预测水平浓淡煤粉燃烧器内的气固两相流动特性,从而为燃烧器的设计和优化提供科学依据。
2. 水平浓淡煤粉燃烧器内气相流动的数值模拟在水平浓淡煤粉燃烧器内,气体流动是影响燃烧效率和稳定性的重要因素。
数值模拟可以通过求解雷诺平均湍流模型(RANS)方程和离散相模型(DPM)方程来描述气相流动的特性,包括速度场、压力场和湍流特性等,从而揭示燃烧器内部气相流动的规律和规律。
3. 水平浓淡煤粉燃烧器内固相流动的数值模拟除了气相流动外,煤粉燃烧器内的固相流动也是十分复杂的,研究固相流动对优化燃烧过程至关重要。
数值模拟可以通过求解颗粒流体动力学(PFD)方程来描述固相颗粒的运动和燃烧过程,其中包括颗粒的输运、碰撞和燃烧过程,为燃烧器内固相流动的规律和规律提供重要信息。
4. 水平浓淡煤粉燃烧器内气固两相流的耦合数值模拟气固两相流动是水平浓淡煤粉燃烧器内最为复杂的部分,气相流动和固相流动之间存在多种相互作用和耦合关系。
数值模拟可以通过耦合求解气相流动和固相流动的方程来综合分析气固两相流动的特性,包括颗粒的输运、燃烧和热力学参数的耦合关系,为水平浓淡煤粉燃烧器内气固两相流的整体特性提供全面的理论支持。
5. 水平浓淡煤粉燃烧器内气固两相流的数值模拟在煤粉燃烧工程中的应用水平浓淡煤粉燃烧器内气固两相流的数值模拟在煤粉燃烧工程中具有重要的应用价值,可以为燃烧器的设计和运行参数优化提供重要的理论和技术支持。
煤粉颗粒团燃烧过程的数值模拟
煤粉颗粒团燃烧过程的数值模拟Chapter 1: Introduction- Background and significance of the research on coal powder particle group combustion process- Objectives and significance of the current study- Brief overview of the research methods and techniques utilized Chapter 2: Literature Review- Overview of the current state of research on coal powder particle group combustion process- Key studies and findings in the field- Gaps and limitations that exist in the current literature- Theoretical framework and models that have been developed Chapter 3: Methodology- Research design and approach- Description of the numerical simulation technique utilized and its key assumptions- Data collection and analysis methods- Case studies and scenarios investigatedChapter 4: Results and Discussion- Presentation of the numerical simulation results- Analysis and interpretation of the findings- Comparison with existing theoretical models and experimental results- Insights gained and implications of the findingsChapter 5: Conclusion and Future Research Directions- Summary of the research findings and their practical implications- Limitations of the study and areas for further research- Recommendations for improving the numerical simulation technique and expanding the scope of the study.第1章节:引言背景和研究的重要性煤粉颗粒群燃烧过程一直是煤燃烧领域的一个重要研究课题。
基于ANSYS CFX的旋转炉内煤粉燃烧数值模拟
第 2期
韩
颂
马 蓉
张浩浩 等
基于 A N S Y S C F X 的旋转炉 内煤粉燃烧数 值模 拟
明显 、 炉 内流 动在 两侧 呈现 较 强的 回流 、 中心形 成射 流 , 该技 术 为研 究 炉 内流动 与 反 应过 程提 供 了
新方法。
关键词 : A N S Y S C F X; 旋转 炉 ; 煤粉 燃烧 ; 数值 模 拟
中 图分类 号 : T D 7 1 2 . 6
i n c e n t e r .Nu me i r c l a s i mu l a t i o n t e c h n i q u e p r o v i d e s n e w me ho t d f o r he t s t u d y o f l f o w a n d r e a c t i o n p r o c e s s i n f u na r c e . Ke y wo r d s: ANS YS C F X;r o t a r y f u na r c e ;p u l v e iz r e d c o a l c o mb u s t i o n;n u me ic r l a s i mu l a t i o n
文献标 识 码 : A
文章 编号 : 1 6 7 1— 7 4 9 X( 2 0 1 5 ) 0 2— 0 0 5 3— 0 3
N唧 e r i c a l s i mu l a t i 0 n 0 f p ul v e r i z e d c o a l c o mb u s t i 0 n i n r 0 t a r v f u r na c e b a s e d 0 n ANS YS CFX
燃煤锅炉燃烧过程流场及温度场数值模拟
毕业设计(论文)题目燃煤锅炉燃烧过程流场及温度场数值模拟系别动力工程系专业班级环境工程班学生姓名指导教师二○一三年六月摘要数值模拟是以电子计算机为手段,通过数值计算和图像显示的方法,达到对工程问题和物理问题乃至自然界各类问题研究的目的。
课题涉及到三维燃烧过程,并带有两相流。
综合考虑,我选择了目前应用比较广泛的FLUENT软件作为数值模拟的工具。
本文对锅炉炉膛计算域通过GAMBIT软件构建三维框架结构,从而对其进行网格划分,确定合适的数学物理模型,设置边界条件,选用适当的变量和参数,对炉膛燃烧进行三维数值模拟,得出炉膛内流场与温度场分布。
最后经过简单的处理,将模拟结果以图片或图表的形式进行直观的展示。
通过对模拟结果的观察分析得出合理的结论,并分析不足之处。
改变燃尽风风速大小,选择30m/s、40m/s、49m/s及60m/s三种燃尽风速,研究燃尽风风速对炉内混合特性和炉内温度场的影响。
结果表明:燃尽风口风速增大时,炉内气流的旋转强度随之增强,燃尽风的穿透程度随之加强,相对容易穿透到炉膛中心,从而使得烟气与煤粉的混合加剧,有利于增加煤炭燃烧的效率;在一定条件下,随着燃尽风速的增加,炉膛中心的高温区域面积增加,而且相对集中;随着燃尽风速的增加,锅炉烟气出口的温度降低;燃尽风风速为49m/s时炉内燃烧状况最佳。
关键词:流场;温度场;数值模拟;燃尽风NUMERICAL SIMULATION ON FLOW FIELD AND TEMPERATURE FIELD OF THE COMBUSTION PROCESS IN THE FIRED PULVERIZED-COAL BOILERAbstractNumerical simulation uses electronic computers as the means. To achieve the purpose of engineering problems and physical problems as well as the nature of various problems, it uses the method of the numerical calculation and image shows. The topic relates to the three-dimensional combustion process and the two-phase flow. Considered, I chose the FLUENT as the tool for numerical simulation.In this paper, establishing the three dimensional frame construction with GAMBIT, carrying on the grid division, then selecting the appropriate model of mathematics and physics and the suitable parameter and the variable, setting up the boundary condition, making three-dimensional numerical simulation of furnace combustion, receiving the distributions of flow field and temperature field in the furnace. After simple processing, we can show the result by making the pictures or diagrams. Making a conclusion from the results and finding out the inadequacies of the results. Changing the size of velocity of over fired air,choose 30,40,49 and 60 meters per second,then discuss what will happen about mixing characteristics of the furnace and temperature field.Study results indicate that As burnout air speed increases, the rotation of the furnace air flow intensity increases, and the degree of penetration strength increases, it’s easy to penetrate into the center of the furnace relatively, so that the mixture of flue gas and coal increases, the increasing efficiency of coal combustion is also in favor.Under certain conditions, as the velocity of over fired air increases,the temperature of furnace center area increases, more concentrated;and the temperature of the boiler flue gas outlet slso reduced;the best overfire air velocity for optimum combustion furnace is 49 meters per second.Keywords:Flow Field;Temperature Field;Numerical Simulation;Over Fired Air目录摘要 ............................................................................................... 错误!未定义书签。
煤粉燃烧数值原理
T 表示的是颗粒沿着其运
′ u′ p u p (t + s)
0
(u ′ )
p
2
ds
(2 13) T 值越大就表明颗粒在流动过程中处于湍
积分时间与颗粒的湍流扩散率正比
流状态时间越长 由上式可知颗粒的湍流扩散率可写作 ui′u′j T 对于在流动区域中具有良好跟踪性 相间滑移速度接近于零 的细小颗粒 颗
ρp 为颗粒密度 Re =
d p 为颗粒直径
Re 为相对雷诺数 颗粒雷诺数
ρd p | up − u | µ
曳力系数 C D 可采用如下的表达式 CD = a1 + a2 a 2 + Re Re (2 12) 当流动
对于球形颗粒 状态为湍流时
在一定的雷诺数范围内 上式中的 a1 , a 2 , a 3 为常数 [13] 使用流体的时均速度 u 通过轨迹方程
3
上海交通大学硕士学位论文
参数的评价仍是可靠的 需要指出的是 明 这些基础前提有可能是不准的 而其正确性又得不到直接证
因此在整个模型的应用中 必须小心的结合实验数据 运用数值模拟方法来 从另一个角度来看 虽然燃烧过程中包含着复杂的
估价模型的精确性与可靠性
微观过程 但其宏观特性却呈现出明显的规律性 包括宏观的温度场 速度场 浓度场 行的 数值模拟的工作可以获得一般实验方法无法得到的信息 包括求出燃烧室的 温度分布和壁面热流分布,分析其热工况情况 括各股气流的混合 流速 湍流和回流情况 的反应过程经历 分析积灰 反应混合情况 求解气相流畅 分析流动工况 包 传热传质 流动等特性 这就表明用数学方法来描述这些过程是切实可
因此可
湍流模型 气体 其表达
的两个输运方程
(2 1) Sφ 是由气相引
燃煤锅炉催化燃烧的数值模拟
燃煤锅炉催化燃烧的数值模拟煤粉催化燃烧技术能够有效转变煤粉燃烧反应中的动力学特征,以此提高煤粉燃烧效率,本文就主要分析了燃煤锅炉催化燃烧的数值模拟,以供借鉴。
煤炭是我国的重要能源,我国能源生产消费中7成左右为煤炭资源。
现阶段,能源生产结构转型升级明显,受到气候变化等因素的影响,我国能源结构更为合理,煤炭在总能源使用中的占比有所减少。
但煤炭资源消耗速度明显快于其发展速度,相关部门需增大煤炭资源利用率。
1 建立双切圆锅炉燃烧模型本文以某热电厂的双切圆锅炉为例,建立双切圆锅炉的三维立体空间几何模型与二维平面数字模型,具体内容如下所述。
1.1建立模型1.1.1三维立体空间几何模型主燃烧器是双切圆锅炉的重要组成部分,是为锅炉提供燃烧反应场所和热力储备的。
与普通锅炉相比,双切圆锅炉的主燃烧器的尺寸规格更小,但是内部结构十分复杂。
若想建立单独且完整的主燃烧器模型,必定会增加模型网格数量,进一步加大模型网格划分难度。
由于双切圆锅炉的内外部结构具有一定的对称性,为简化锅炉结构,减少模型网格数量。
在绘制锅炉的主燃烧器时,结合主燃烧器的结构特点,使用水力直径相同的长方形替代锅炉的一次风进口和二次风进口,以便更容易的识别。
根据图1(锅炉的区域划分图)可知,将整个双切圆锅炉的炉膛区域划分为炉膛下部、炉膛中部、炉膛上部三个重要组成部分。
图2是锅炉内部立体空间方位图,其中X+、Y+与Z+分别代表水平方向、平面垂直方向和立体垂直方向三个维度。
1.1.2数学模型本文采用的数学模型如下:(1)基本控制方程组连续性方程:式中:代表气流密度;代表煤粉中水分蒸发量;S代表脱挥发引起的质量变化。
动量守恒方程:式中:p代表炉膛静压力;μ代表动力粘度,τij 代表应力张量;ρgi代表重力体积力;Fi代表两相之间体积力。
能量守恒方程式中:E代表流体微团总能;K1代表湍流导热系数;Keff代表有效导热系数;Jf 代表组分,J1的扩散通量;(τij)sff代表热剪切应力;Sh代表能量源项;h代表静态焓值;mf代表组分,J1的质量分数,Tref=298.15k。
煤粉工业锅炉燃烧数值模拟简介
通过对新型高效煤粉工业锅炉燃烧器和煤粉工业锅炉炉膛进行数值模拟,对炉膛内气固两相流动、温度场分布、各气相组份浓度场分布、发热率分布、煤粉颗粒的燃尽情况、污染物生成量等进行预测,并考察各输入量的变化对炉内燃烧的影响;通过将数值模拟的结果与实际运行数据进行验证,对存在的问题进行分析,并对优化运行提出建议。
CFD数值模拟概述
CFD是计算流体力学的简写(Computational Fluid Dynamics),其基本的定义是通过计算机进行数值计算和图像显示,分析包含流体流动、热传导、化学反应等现象的系统。近年来,CFD有了很大的发展,所有涉及到流体流动、热交换、分子运输等现象的问题,
几乎都可以通过计算流体力学的方法进行分析和模拟。CFD不仅作为一个研究工具,而且还作为设计工具在水利工程、土木工程、环境工程、食品工程、海洋结构工程、工业制造等领域发挥作用。
C. 验证及优化:根据数值模拟计算的结果,对燃烧器结构以及锅炉受热面布置提出优化建议。
数值模拟的主要研究内容
A.冷态单向流数值模拟:模仿冷态试验过程,对燃烧器进行相似模化,考察燃烧器内部的流场分布情况和燃烧器外部回流区的形成情况;
B.热以及调节各输入量(如风速、风温、送粉量、挡板开度等)对着火特性的影响;
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
煤粉燃烧过程的数值模拟
Ryoichi Kurose
京都大学
高级研究院流体科学与工程学院机械工程与科学系
Hiroaki Watanabe and Hisao Makino
中央研究所的电力行业能源工程研究实验室
摘要
煤炭是一种能够满足电力进一步需求的重要能源资源,而且煤炭比其他化石燃料的储量丰富得多。
在燃煤发电厂,改善对环境污染物如NOx,SOx及包括未燃尽的碳粒在内的灰粒的含量的控制技术十分重要。
随着计算机性能的显着提高,人们强烈希望计算流体动力学(CFD)成为一种工具,成为一种研发和设计这种合适的煤粉燃烧的燃烧炉膛和燃烧器的工具。
这次审查的重点是突出我们的CFD 研究的最新进展,即煤粉在燃烧中的平均雷诺数纳维斯托克斯(RANS)的模拟和大涡模拟(LES)的最新进展,及未来的一些前景。
关键词:煤粉燃烧,数值模拟,平均雷诺数纳维斯托克斯模拟,大涡模拟1.介绍
煤炭是一种能够满足电力进一步需求的重要能源资源,而且煤炭比其他化石燃料的储量丰富得多。
在燃煤发电厂,改善对环境污染物如NOx,SOx及包括未燃尽的碳粒在内的灰粒的含量的控制技术十分重要。
为了实现这些目标和要求,了解煤粉燃烧机理和先进的燃烧技术的发展十分必要。
然而,由于煤粉燃烧是一个非常复杂的现象,其中最高的火焰温度超过1500℃,以及某些物质难以进行测量,如一些原子团种类和一些高活性固体颗粒,因此在燃烧过程中的煤粉燃烧机理至今没有得到很好的解释。
而且由于研发过程包含许多步骤,因此,新的燃烧炉膛和燃烧器的发展需要较高的成本和较长的时间。
随着计算机性能的显著提升,煤粉燃烧领域的计算流体动力学正在被研发。
在这种方法中,电脑解决了燃烧领域的控制方程式,这使它能够提供温度和化学物质种类分布的详细信息和在整个燃烧空间中煤粉颗粒的行为,而上述那些通过实验是不能得到的。
此外,此种方法有助于在相对较低的成本条件下重复审查任意条件下的煤粉的流场和各种参数。
因此,强烈地希望计算流体动力学(CFD)能够成为燃烧炉炉膛和燃烧器研发和设计的一种工具。
直接数值模拟(DNS),大涡模拟(LES)和平均雷诺纳维斯托克斯模拟(RANS)是有关湍流模型中燃烧领域中的典型的计算流体动力学(CFD)方法。
直接数值模拟(DNS),它能通过在这些领域中设置低于最小漩涡的数值网格空间来直接解决流场、化学物质的浓度和温度的控制方程,在上述提及的方法中,它具有最高的数值精确度。
虽然它能有效地应用于基础研究,但是它在实际水平的燃烧领域应用却很难,因为它需要设置大量网格点和高负荷计算的计算机。
相反,平均雷诺纳维斯托克斯模拟(RANS)在实际应用中最常用。
因为这种方法通过在时间和置换雷诺切应力加以平均解决了控制方程和湍流标量通量的湍流模型。
它可以大大减少网格点的数量和计算机负载。
然而,平均雷诺纳维斯托克斯模拟有几个问题,如湍流模型的选择和其中包含的参数值的确定的困难,以及在预测不稳定的湍流运动中存在一定的劣势。
因此,最近大涡模拟收到了重视,这种方法直接解决了较大旋流的控制方程,剩余的小漩涡使用模型计算。
这是一种空间平均法,它具有一定的优势,即对不稳定的湍流运动进行了评估和模型中包含的参数数量减少了。
与平均雷诺纳维斯托克斯模拟相比,尽管大涡模拟(LES)对计算机有较高的负荷,但是如果考虑到计算机性能的改进取得进展,但是很可能在不久的将来LES会被应用到实际领域中。
这次审查的重点是突出我们的CFD研究的最新进展,即煤粉在燃烧中的平均雷诺数纳维斯托克斯(RANS)的模拟和大涡模拟(LES)的最新进展,及未来的一些前景。
在第2部分中,介绍了和现在数值模拟有关的煤粉燃烧试验炉膛和燃烧器以及煤粉燃烧锅炉的简明概念。
一些LES和RANS数值模拟结果分别显示在第3部分和第4部分。
最后,本次审查在第5部分下结论。
2.燃煤电站锅炉和基本测试实验设施及其数值模拟在这次审查中,主要关注的燃煤锅炉的燃烧场数值模拟。
锅炉是一个将热量转化为蒸汽热能的系统。
在通常的煤粉燃烧锅炉中,煤粉由安装的燃烧器提供,它安装在垂直方向的三个位置,如图1所示。
锅炉的燃烧器的数量为10-40个,每个燃烧器的给煤率最多为10吨/小时。
根据燃烧系统的不同,电站锅炉一般分为两类:即相对布置在前墙和后墙的燃烧器的锅炉(见图1)和布置在锅炉的四个角的四角切圆燃烧器的锅炉。
对于这些锅炉,相对于了解水平方向的火焰相互影响而言,了解垂直方向的火焰相互影响更为重要,因为火焰弯曲向上的主流和垂直方向的浮力使得火焰之间彼此影响。
为了从根本上研究火焰相互作用,中央电力行业研究所(CRIEPI)采用多燃烧器的煤粉燃烧试验炉膛,在Yokosuka实验室的试验炉中,垂直方向安装有
三个燃烧器,每个燃烧器具有约100千克/小时的燃烧能力,如图2所示。
多燃烧器炉膛紧挨着烟气处理装置,比如选择性催化脱硝装置,静电除尘装置,石膏石灰石湿式脱硫装置,这很类似电站锅炉。
此炉的高度,水平宽度,深度分别为11米,0.9 米和1.9米。
助燃空气通过燃烧器和分级布置的风口注入炉内,风口位于从燃烧平台至其以下2米的范围内。
此外,日本电力工业中央研究所的Yokosuka研究实验室也采用单燃烧器的燃煤试验炉(见图3),以进一步了解单一火焰的详细结构。
煤燃烧容量大约为100千克/小时。
单燃烧器锅炉是一个圆柱形锅炉,煤粉燃烧的火焰不弯曲且不相互扰动,从而可以检测其燃烧过程。
该炉的直径为0.85米,长度为8米。
助然空气通过燃烧器和分级布置在燃烧器以下3米范围内的风口注入炉内。
通过比较这两个试验炉之间的燃烧特性,希望获得在设计和运行实际电站锅炉方面有用的信息。
安装在上述试验炉中的燃烧器是日本电力工业中央研究所高级低污染高性能燃烧器(CI-α),它的燃煤能力大约为100千克/小时。
燃煤和一次风共同送入炉膛,二次风和三次风通过一次风的周围边界供应。
日本电力工业中央研究所高级低污染高性能燃烧器(CI-α)是这样设计的,它通过一次风的直接运动和二三次风的旋转运动(见图4)而产生的再循环流动来促进燃烧器区域的混合。
再循环流动延长了煤粉颗粒在燃烧器出口附近高温区域的停留时间,并且加速了挥发分物质的演变和焦炭反应的进程。
因此,有效地减少了飞灰中未燃尽碳的数量,但是在这个区域中氮氧化物的浓度也增加了。
在再循环区域之后存在着减小的火焰,然后氮氧化物立即减少变为氮气。
而且,分级燃烧的方法促进了氮氧化物的减少效果,在这种方法中,一些和燃烧器分开的燃烧空气通过安装在炉膛后部的注射风口供应。
有关日本电力工业中央研究所高级低污染高性能燃烧器(CI-α)和其性能的详细描述能够在我们以前的报纸中找到。
4.3燃煤领域的挑战
最近,在带有CI-α燃烧器的单一燃烧器炉膛中,大涡模拟(LES)已经应用到实际的燃煤领域中。
图21显示了在流动方向的平面上瞬时轴向速度,燃气温度和氧气浓度的预测分布。
能观察到在靠近燃烧器的中央区域形成了旋转的再循环流动及它的大小和动态特性随时间的变化。
同时,这种不稳定流动行为强烈影响了点火特性。
图22显示了平均时间烟气温度和氧气浓度的轴向分布在计算和实验中的比较。
定性预测分布和实验结果一致。
因此,不像平均雷诺数纳维斯托克斯(RANS)模拟,大涡模拟(LES)能够精确地捕捉这种不稳定燃烧流运动,它的这种特性具
有极大地提高数值精确度的可能。
此外,大涡模拟(LES)较之于平均雷诺数纳维斯托克斯(RANS)模拟有一个优点,就是它的模型参数数量较少。
然而,对于燃煤锅炉的大涡模拟(LES),仍然存在着困难的挑战,就是它不但涉及湍流燃烧的数学模拟和污染物的排放,而且涉及计算成本。
5.结论
在此综述中,有关平均雷诺数纳维斯托克斯(RANS)模拟和大涡模拟(LES)的煤粉燃烧计算流体动力学最新进展得到了讨论。
可以这样说,有关流场的预测,大涡模拟(LES)比平均雷诺数纳维斯托克斯(RANS)模拟更有优势,因此,大涡模拟(LES)希望在不久的未来成为新燃烧炉膛和燃烧器的发展和设计的有用工具。
然而,煤粉的燃烧时一个复杂的现象,在这种过程中包括各式各样的过程,比如煤粉的液化,挥发分物质的燃烧,在颗粒表面的燃烧,以及在气流中煤粉的分散行为及他们之间的相互作用。
因此,煤粉燃烧行为还没有被完全理解和模仿。
尤其是最近煤的液化模型能够强烈地影响数值的精确度。
湍流燃烧模型的提高也是必须的。
更为复杂的湍流燃烧模型,在这种模型中,小火焰单元模型等详细作用信息应当被考虑,它常用来气态的和喷雾液体的燃烧,这种模型是可取的。
此外,大涡模拟(LES)的更大更复杂领域比如在电站锅炉中的应用在数值稳定性、准确性和成本方面受到挑战。