新建75th高压高温水煤浆锅炉炉内冷态空气动力场 燃烧过程及NOx排放试验和数值模拟

《新建75th高压高温水煤浆锅炉炉内冷态空气动力场 燃烧过程及NOx排放试验和数值模拟》由会员分享，可在线阅读，更多相关《新建75th高压高温水煤浆锅炉炉内冷态空气动力场 燃烧过程及NOx排放试验和数值模拟（86页珍藏版）》请在装配图网上搜索。

1、聊旺大学硕士学位论文周平摘要摘要水煤浆是一种低污染、高效率、流动性能强的代油清洁燃料。为了降低化工产品的生产成本,提高企业的经济效益，中国石油吉林石化分公司委托浙江大学 设计并由南通锅炉厂制造了7t/l高压高温水煤浆锅炉,该锅炉在国内属于第一 台。锅炉冷态试验和燃烧调整试验表明,锅炉的设计是成功的，为进一步推广水煤浆锅炉在石化企业中的应用起到了示范作用。本文对四角切圆燃烧水煤浆锅炉炉内的流动、传熟、燃烧过程及氮氧化物 （X)的生成过程作了详细的数值研究，并尽可能结合实炉试验，对理论进行验 证。本文同时应用数值研究方法模拟了四角切圆锅炉一二次风大小变化、大小同 心切圆布置方式、一三角相切二四角对

2、冲、二次风配风方式改变,二次风反切等 对炉内流场特性的影响,为优化锅炉设计、改造、运行提供了及其有意义的理论 上的指导。关键词：高压高温，水煤浆,冷态试验，切向燃烧锅炉,数值模拟，NOx86 / 86浙江大学硕士学位论文周平ABSACToalwae surry i te of lowpolution,hfficiency，perfe liidty ul as sutitutio fooiIn orr iminish the cot animrve the benei fehemcal entrrise,one75h i【sh emprare and ih ressure coa.atr sur

3、ry iler as created b ntong iler Mauacory， which wa mmiioned y te gnicythtical Fcry o Jilin Petrchemial Co,addeiged Zhejn niverity。The oir oigialin Chinae col mel exrimnt and comuston adjust xperientidicate ht the desgnof t oeis sucesulndfntion s ademonsrtionf urtrin th aplaton of coalwate slurr iner

4、ochmical cororaions。In thispaer，etaledattention ha beenpaid on the modelingfor low，at trans，cobusi,rdon，NOx rucion nd oeetifare o on theo y xperiensManwhe smedtale umicleserces e on teiluene tthe aerdynic fied togho changng magnitude f riary airad the oayair，adoptig homocent two angentialcicles，tang

5、nially-fred in umbrone and numberthr comerfafired i nmetwo an numb fur oer，chang th ir distrion ofthscondry airandh evre wil of uppersecndry air,Which pvid hretical basis forbole design，otima eatin,reovatinecKeywods：HPandHT,co-water sluy,old moe exim,taentily fred ier，nmricalsmulaio,N浙江大学硕士学位论文周平第一章

6、绪论第一章绪论1。1引言水煤浆是由6075不同粒径分布的煤、3左右的水和约l的添加剂制成的煤 水混合物，是七十年代石油危机发展起来的一种新型低污染代油燃料。它既保持了煤炭原有 的物理特性，又具有石油一样的流动性和稳定性,被称为液态煤炭产品。水煤浆技术包括水煤浆制备、储运、燃烧等关键技术，是N涉及多门学科的系统技术。水煤浆具有燃烧效率高，污染物排放低等特点,可用于电站锅炉、工业锅炉和工业窑炉代油、代气和代煤燃烧， 亦可作为气化原料.用于生产合成氨、合成甲醇等。水煤浆技术是我国现行阶段适宜的代油、环保、节能技术.用煤制取清洁燃料,【;l煤代油.是我国能源长期稳定的发展战略和现实选择。吉林石化公司现

7、有燃油蒸汽有机锅炉1 l台，2001年消耗燃料油总量3。7万吨年。由于 目前燃料油价格高,导致燃油锅炉产生的蒸汽成本提高同时增加了化工产品的生产成本。 由于水煤浆具有以上等优点，并且燃烧效率、锅炉效率能达到燃油同等水平,且燃烧调节方 便，运行稳定可靠，贮存、运输比燃油简单、安全，便于管理，适宜于管道运输.因此，以水煤浆替代燃油，改变有机锅炉的燃料结构在技术上是可行、可靠的，经济上是合理的。中国石油吉林石化分公司为了降低化工产品的生产成本，提高企业的经济效益，新建了75t高压高温水煤浆锅炉工程，以满足蒸汽有机锅炉由燃油到燃水煤浆的燃料结构改变的 需要,以提升企业的综合竞争能力,适应未来市场发展的

8、要求.经初步预算:燃油锅炉的产汽成本在013元吨左右;水煤浆锅炉的产汽成本在0 90元吨左右.新建t/h高压高温水煤浆锅炉，年可节省燃料油、汽39万吨。该水煤浆锅炉由吉林石化有机合成厂和南通锅炉厂委托浙江大学设计,其容量和设计参数在国内尚属首台.该锅炉于03年10月安装完毕，之后进行了锅炉冷态空气动力场试验，并于003年2月和20年4月份2个阶段完成了燃烧调整试验。热态试验进行了全烧油、油气混烧、浆气混烧、全烧浆工况的性能测试,笔者参与了全部调试过程,试验结果达到了合同规定的技术要求和技术指标.锅炉经燃烧调试后，锅炉运行情况良好，并一至持续运行至今。总之。吉林石化公司有机合成厂国内首台新建75

9、讹高压高温水煤浆锅炉的设计和调试是成功的，锅炉负荷和蒸汽参数均符合设计技术指标的要求，并满足机组安全经济运行要求。 为进步优化锅炉运行,本文将从现场试验和计算机数值模拟相结合的方式研究四角切 圆燃烧锅炉炉内空气动力场特性三I及水煤浆燃烧过程特性、污染物排放特性等，用以指导水煤浆锅炉设计和实际运行。1四角切圆燃烧锅炉的优点、目前存在的主要问题四角燃烧锅炉广泛应用于燃煤和燃油以及水煤浆锅炉.通常直流式燃烧器布置在炉膛四角其出口气流几何轴线切于炉膛中心的假想切圆，造成气流在炉内强烈旋转。炉膛四周 是强烈的螺旋上升气流,在中心则是速度很低的微风区.四角切向燃烧锅炉由于其燃料的适应性及风粉混合均匀等特点

10、决定其是我国电站应用 最广、最成熟的燃烧方式，应用的锅炉容量从小到大不等，使用从优质煤到劣质煤几乎所有 的煤种 。对于水煤浆,更是普遍采用了四角切圆的燃烧方式.多年运行实践表明，切向燃烧锅炉具有以下优点阱:.。 （1)四角射流着火后相交，相互点燃,使煤粉着火稳定是煤粉着火稳定性较好的炉型。 由此可见，切向燃烧方式是以整个炉膛为单元来组织燃烧的，故燃烧器的燃烧工况和整个炉浙江大学硕士学位论文周平第一章绪论膛的空气动力特性关系十分密切。（2)由于四股射流在炉膛内相交后强烈旋转，湍流的热量、质量和动量交换十分强烈,故能加速着火后燃料的燃烬程度.（3)四角切向射流有强烈的湍流扩散和良好的炉内空气动力结

11、构,旋转火焰气流横扫过大部分炉膛容积，使炉膛火焰充满系数很高，整个炉膛形成一个特大型旋流燃烧室,燃烧充 分,炉内热负荷较均匀。（）旋转火焰气流使燃料、空气和烟气强烈混合,并随旋转气流呈螺旋状上升,直到炉膛出口,延长了煤粉粒子在炉内的行程,有利于煤粉的燃烬,减小了飞灰损失，经济性好。 （5)切圆燃烧时每角均由多个一、二、三次风喷嘴所组成,负荷的变化调节灵活，对煤种的适应性强，控制和调节手段也较多。 （6)炉内结构简单，便于大容量锅炉的布置.（7)采用摆动式直流燃烧器时运行中改变上、下摆动角度即可调整火焰中心高度,改 变炉膛出口烟温，达到调节过热蒸汽温度的目的。(）便于实现分段送风组织分段燃烧,从

12、而抑制q的排放.旋转火焰使炉内火焰温度 比较均匀,不易出现局部高温区，使得炉膛四面水冷壁的热负荷也较均匀,可避免出现膜态 沸腾。也有利于降低NO的生成.尽管四角切向燃烧锅炉在我国乃至世界范围内已得到了广泛的应用，在实际的运行中仍存在着许多问题。四角切向燃烧锅炉存在如下主要问题jJ：()污染物的生成锅炉燃烧所产生的污染物中，NOx对人类环境危害较大,会造成大面积酸雨,还会破坏臭氧层,如何有效地控制NO生成不仅已成为人们普遍关注的焦点，还 是能源发展面I临的重大障碍之一。并且随着四角切向燃烧锅炉容量的增加，断面热负荷越高。炉内温度也较高。煤粉的挥发份释放和焦炭燃烧速率也较大，使燃料型哝的比例较高，

13、同 时热力型NO的份额也大幅度提高。本文将对酞的生成进行数值模拟。(2）结渣锅炉结渣问题的存在非常普遍。在炉膛内煤粉燃烧过程中,烟气中夹带的熔 化或部分熔化的颗粒碰撞在炉壁上,水冷壁或管壁上被冷却凝固形成结渣。结渣的形态主要 是粘稠状的熔融沉淀物，结渣主要出现在辐射受热面上。结渣轻则影响传热，迫使锅炉降负 荷运行,降低锅炉的效率；重则导致非计划停炉或造成重大安全事故，是危及锅炉安全经济 运行的一大难题。结渣的内因是燃料的灰分特性;外因是炉内的燃烧工况。从内因上讲，我国的燃料政 策是电站锅炉尽量燃用劣质燃料。统计表明,我国大机组所用煤质约有半数属易结渣类型: 而且电厂燃用煤质多变,实际燃用煤质与

14、设计煤种相差较大；炉膛结构特性，燃烧器的结构 布置与燃料不匹配也容易出现结渣问题。燃用劣质煤的锅炉为了保证燃烧稳定，通常在设计上选择较高的热负荷和采用较保守的稳燃措施，这虽可取得较好的稳燃效果.但却导致结 渣。从外因上讲，锅炉运行方式不当、速度场和温度场分布不合理、火焰偏斜、冲墙、局部 出现还原性气氛也是导致结渣频繁发生的原因。锅炉结渣已经成为阻碍锅炉技术向前发展的巨大障碍.本文将结合炉内空气动力场特性及燃烧过程说明此问题。（3）烟温偏差近些年部分大容量火电机组非计划停运事故的统计分析表明:锅炉事故 最多的是四管爆破。即水冷壁、省煤器、过热器和再热器的管子爆破。过热器和再热器爆管所占的比例最大

15、，又因过热器工质的压力最高，温度最高，所以事故扩大性更为严重。超温 爆管的发生主要来自两方面的原因：锅内过程和炉内过程。由于锅内过程的设计是针对炉内 过程进行的,寻找并解决引起爆管的炉内过程方面的因素是解决超温爆管问题的基础。切圆燃烧锅炉水平烟道中存在着不同程度的烟速、烟温偏差,而且随着锅炉容量的增大偏差也增大，速度和温度偏差是造成高温对流过热器超温爆管的直接原因。国内部分大 机组锅炉运行情况表明，其水平烟道左右两侧的烟温偏差都在10”C以上，高者甚至可以达2浙江大学硕士学位论文周平 第一章绪论到270C左右。过热器和再热器超温爆管问题目趋严重,国产大容量电站锅炉及一些总体质 量上很好的进口大

16、容量电站锅炉上都发生过受热面超温问题.这类问题在技术上有一定的复杂性,其研究对于锅炉的安全性和经济性具有重要的意义.本文将对此作出定性分析。()低负荷稳燃效果差,飞灰可燃物也较高,特别对于劣质煤层、贫煤、无烟煤等，在 负荷低时必须投油助燃,而多数四角锅炉在燃烧稳定性、着火性能上都不理想，且飞灰中可 燃物含量过高。着火和燃烧的稳定性除了与煤质有关外,还与燃烧的组织有密切的关系。燃 烧器形式和布置、锅炉结构和运行方式（主要是风粉配合），要保证炉内合理的空气动力场和温度场。.3数值模拟的主要优缺点发现问题、分析问题并尽可能地去解决它,这无疑是人类认识自然和改造自然的一般 规律。本文正是基于这样的一个

17、思路，尝试用数值模拟的方法对锅炉冷态空气动力场、燃烧 特性及污染物排放进行预报，以达到为锅炉机组的优化设计和优化运行提供技术指导的目的。数值模拟具有如下优点14：（1）耗费少在许多应用中，计算机的计算费用要比相应实验费用低好几个数量级.当需 要研究的问题变得更复杂时，这个因素所起的作用更加重要.(）速度快计算研究能以非常快的速度进行。只要模型合理。计算所花费的时间将比试 验所用的时间少。(3)完整的资料一个问题的计算机解能给出详细和完整的资料。能提供整个计算区域内的所有相关变量的值(例如速度、压力、温度、浓度、湍流强度等）。试验研究不可能测出整个区域内所有变量的分布，即使完成了实验，在很大程度

18、上还需要计算机解去完善补充 实验资料。(4)模拟真实条件的能力在理论计算中，很容易模拟真实条件，没有必要求助于小尺寸和冷流动模型。虽然数值模拟具有以上优点，但计算的使用价值受数学模型和计算方法是否合理所限制。若使用了不适当的数学模型，即使是一个完全满意的数值方法也不能得到有价值的结果. 锅炉的试验研究，即使规模不太,从设计到实施都要消耗大量的人力和物力. 随着计 算机内存容量的不断增大，运算速度的加快，数值模拟技术与试验研究日益结合起来。数值 模拟是对试验研究成果的修正和补充,并且能够提供试验研究无法实施的结果。比如受到试 验条件和测试手段的限制,克服试验过程中难以避免的某些系统误差和偶然误差

19、：对设计方案或改造方案进行初步论证,减少试验工作的盲目性和重复性，降低工作量和资金费用。4电站锅炉冷热态数值模拟的研究现状对于电站锅炉炉内过程的数值预报,涉及计算流体力学、燃烧学、传热学、气固两相流、化学反应动力学等诸多领域.借助这些学科和计算机科学的发展，计算模拟炉内过程逐步成为一种强有力的研究手段。141研究气固两相流的两大经典方法炉内燃烧特性的数值模拟,首先是气固两相的计算流体力学问题。涉及气固两相流，则 不得不提到研究气固两相流的两大经典方法：欧拉法和拉格朗日法。前者从空间场的角度对流体现象进行描述，研究某一时刻质点群形成的空间场量;后者从时间场的角度对流体现象 进行描述,研究某一质点

20、在不同的运动时刻形成的时间场量。两种方法在本质上是相通的，浙江大学硕士学位论文周平第一章绪论可以相互转换。处理气固两相流时,欧拉法将固相当作连续相:拉格朗日法将固相当作离散 相，通过在拉格朗日坐标系下对大量颗粒的轨迹跟踪统计其物理量。两种方法在处理不同问 题时具有不同的优缺点。欧拉法由于可以方便地解决颗粒间的碰撞问题因而在处理稠密两相流动中得到充分的应用。它可以得到整个流场任何一个位置的运动情况，并且引入了相的容积分额概念，将固相和气相在数学上都当作相互渗透的连续相处理；拉格朗日法则是跟踪某 个湍流微团(或者是某个颗粒)，研究其在不同时刻内运动轨迹和状态的变化,它由于可以方 便地计入颗粒在运动

21、过程中沿程的物理化学变化，因而在处理有化学反应的气固两相流中得 到长足地发展，它需要假定固相占有很低的容积份额(固相的质量大于气相是容许的)。所以 在对整个流场作研究的时候，一般都采用欧拉方法,在求解运动轨迹的时候用拉格朗日方法. 现在有许多学者致力于这两种方法的结合。应用它们各自的优点更反映实际地去模拟各种气 固两相流动问题.14.2气相湍流流动模型综述处理湍流流动大体上有三种方法：直接模拟（DietNumeical Simulation，简称DNS）、雷 诺时均法和大涡模拟。直接模拟的方法可以获得复杂几何形状体内的高雷诺数湍流流动N。S 方程完整的时间相关精确解,但所需要的计算机的速度和内

22、存按现有的计算机水平是无法想 象的，在可预见的将来实现起来也不是容易的。在无法直接模拟的情况下，有两种替代性的 方法，一是雷诺时均法一是大涡模拟。这两种方法的特点是不直接模拟小尺度的湍流脉动。 无论是雷诺时均法还是大涡模拟。它们在控制方程中都引入了附加项。为保证控制方程的封 闭,需要对附加项进行模化:补充新的方程将附加项与控制方程中原有的项联系起来。雷诺时均法通过对NS方程时均化,求解平均流动量,通过采用的湍流模型描述所有尺度的湍流。这种处理方法大大减少了计算的工作量。在实际工程问题中被广泛采用。雷诺时 均法引入了雷诺应力项，根据处理附加雷诺应力项方法的不同，可以分为两大类:ousinq 模型

23、和雷诺应力模型.Boussiesq模型是将雷诺应力与平均速度梯度联系起来.像Salartllmras模型例、k一占模型1、k一占模型的改进型(NG k一占川、Realizble k一占嘲）、k一国模型驯和km模型的改进型（SSTk一)都属于此类。这种模型的优点在于计算湍 流粘性系数。时的代价很低。缺点是。各向同性的假定与实际不是很吻合。雷诺应力模型 (RM）”uJll【I“需要求解雷诺应力张量每一项所对应的输运方程。通常情况下基于Bouseq 假设”的湍流模型就可以取得很好的计算结果，如果采用雷诺应力模型的话，计算结果可 能会有改进，但相比与计算花费而言，将是得不偿失。当然，对于各向异性占主导

24、的湍流流动,雷诺应力模型的优越性表现无余。比如高强度的旋流和应力诱发的次级流动,采用雷诺应力模型是比较合适的,根据求解雷诺应力的方程又可分为雷诺应力微分方程模型 (Dffeal Strss de，简称DSM）、代数应力模型(Aebral Sts Mode,简称ASM)。 AM模型是DSM的简化模型.典型的代数应力模型有Craflaner模型4】、Shh模型115】【16、 Rodi模型【17。大涡模拟是近几十年来发展的一种具有一定应用前景的数值计算方法，其基本思想是: 紊流流动是由许多大小不同尺寸的涡旋组成,大尺度的涡旋对平均流动的影响比较大，各种 变量的紊流扩散、质量、动量和能量的交换以及雷

25、诺应力的产生都是通过大尺度涡旋来实现的而小尺度涡旋主要对耗散起作用。通过耗散脉动来影响各种变量。不同的流场形状和边 界条件对大涡有较大影响，使它具有明显的各相不均匀性。在高雷诺数下，小涡旋近似于均匀各相同性,受边界条件影响较小，有较大的共同性。因此可把紊流中大涡旋和小涡旋分开 处理，大涡旋通过修正的NS方程【lq直接求解，小涡旋可通过亚网格尺度模型,建立与大尺 度涡旋的关系对其进行模拟.亚网格尺度模型主要有:动力与代数亚网格尺度模型 、方浙江大学硕士学位论文周平第一章绪论程亚网格尺度模型“、动力双系数亚网格尺度模型”、可压缩流动力亚网格尺度模型o“。大涡模拟还处于工业应用的起步阶段.工程应用中

26、一般采用雷诺时均法，其中以k一占 双方程模型以及它的各种修正模型应用最为广泛。14.3气固两相流动模型综述流动时间（停留时间)f=叫扩散驰豫时间=；.瓦/18平均运动驰豫时间o。=TrllRe；3口；l脉动时间=上：兰碰撞时间r，=“j=(c厅#)。：）其中RP=lDupl以”根据上述特征时间的比值所组成的各个准则，可以定性地判断几种极端情况下的流动类型:T,f,1为无滑移流(平衡流）、TrI/，口为强滑移流（冻结流）、TrzT0 1为扩散冻结流（可忽略颗粒湍流扩散)、叶口为扩散平衡流(颗粒扩散等于流体组分扩 散)、0lfP 0 为稀疏悬浮流、rlf口为稠密悬浮流.可以看出r一/口,以）（以）

27、,其中以为颗粒相的体积分数。表1。4颗粒相模型;黧兰蒸麓篓竺兰翌竺兰竺 絮擞磬亲繁不考虑有嚣无c=q=五=o,扩散冻磊 无滑篓磬州鬈部分地考虑蒜嚣欧拉引h 2黼,2，j鬲历厚-一各种模型都是由不同角度对真实过程所作的近似和简化。不同模型使用范围大致取决于颗粒相相对浓度(体积分数或质量载荷比）和几种特征时间之比。这可以定性地用表1。4。 和图1。1来加以表示。浓度极低（丸极小）和相对滑移量很大（f，.f，很大)时为颗粒冻结流简化模型区。砟和。/r，稍大些的区域为单颗粒动力学简化模型区，这时颗粒浓度仍相当低,相对滑移量仍相当大（颗粒尺寸较大和颗粒Reyos数较大)，因而颗粒对流体作用可忽略，但已考

28、虑流体对颗粒的作用。处于中等的砟和f，f，值的大部分区域是完整两相流模拟区其中又可又可分成较大的,叶值与较小的。值（即扩散相对较小和颗粒浓浙江大学硕士学位论文周平第一章绪论度偏小）的确定轨道模型区，较小的r 区。随机轨道模型区中有一部分和拟流体模型区重叠，后种情况意味着，在足够大的妒。及足够小的.f，的范围内，两类模型均适用。当颗粒浓度屯相当大而肛r又较小时，即稠密两相流而又滑移量不大时，为拟流体模型区,这时轨道模型不适用.在极小的r肛，（滑移量趋于0）的范围内则为无滑移或单流体简化模型区。即两相平衡流区。衰1.42颗粒相模型适用范围圈I.颗粒相模型适用范围圈示墅垫塑苎翌生!塞立!二l(王量王

29、l!。冻结流 一 极大 一 一最早期（0年代至60年代初期)的单颗粒动力学模型是假设颗粒对流场无影响，考察已 知流场(通常是均匀的速度场和温度场）中颗粒的时平均或对流运动的轨道，以及颗粒速度 及温度沿轨道的变化。这是一种极端的简化模型。N60年代后期及7年代初期提出了两相流的单流体模型或无滑移模型241，即颗粒相与流体相达到动力平衡及热平衡，认为空间到处颗粒与流体时均速度与温度相等，颗粒扩散则 相当于流体组分的扩散（扩散平衡)，把颗粒与流体作为统一的流体来加以研究,这是另外 一种极端的简化模型.由60年代质期开始,.So（苏绍礼)开始提出用多相连续介质即颗粒群的小滑移模 型”来描述两相流，其中

30、对稀疏悬浮流虽忽略颗粒对流体的作用。然而对颗粒则由单颗粒的描述过渡到颗粒群拟流体概念，只考虑颗粒群集体速度与温度，承认颗粒与流体间的滑移， 但认为颗粒的滑移与其扩散溧移是同一现象的两个方面,可称为“小滑移模型,这是建立 较完整的两相流模型的开端。由7年代中期以后逐渐发展了更完整的两相流模型,即完整地考虑相间速度与温度滑 移(包括非起源于扩散的大滑移），颗粒扩散,完整地考虑相间耦合，包括质量，动量及能 量间耦合，特别是颗粒对流体的作用，以及颗粒反应的模型,其中又主要分成颗粒相轨道(包 括不考虑相间紊流扩散的颗粒群轨道模型和考虑相间紊流扩散的随机轨道模型)模型（前者般采用单元内颗粒源法(Parti

31、le urce in Cell thod”,简称PSI法）求解,后者一般按照Monte-Carlo法求解)及颗粒相拟流体模型两大类。颗粒群轨道模型的基本思想是：在Lagagan坐标系下考察各组颗粒群沿各自轨道运 动、质量损失及能量变化过程，在欧拉坐标系下处理气相场，由颗粒群沿轨道中阻力、蒸发、 挥发或燃烧产生颗粒群速度、温度、尺寸变化作为质量、动量及能量源加入气相场中成为气 相与颗粒相之间相互耦合.在用Lagrana法描述颗粒运动时,由于只考虑颗粒沿自身轨迹 互不干扰的运动，并认为沿轨迹颗粒数目不变,这就意味不考虑颗粒的湍流扩散效应。因此 必须进行修正。一类是认为的加入漂移速度鲫或漂移力的概念

32、，另一类是用随机的方法来处理颗粒的运动随机轨道法,是在计算颗粒的随机轨道时考虑气相脉动随机速度对颗粒运 动的影响。采用随机轨道模型。一般假定气相速度脉动各向同性并符合高斯分布，以气相湍 流脉动速度的平均平方根值为基础,给予不同的随机数以达到模拟颗粒的湍流扩散特性.浙 江大学岑可法,樊建人等认为一般的随即轨道模型未能考虑到各个湍流涡团的不同的脉动频6浙江大学硕士学位论文周平第一章绪论率、频谱、振幅和方向性而仅以气相湍流脉动速度的平均平方根为基础而提出了脉动频谱随机轨道模型【28l,用随机傅里叶级数来模拟气流的湍流脉动速度。 颗粒拟流体模型的基本思想就是把颗粒群看成是与流体互相渗透的拟流体它可以较

33、完整的考虑颗粒相的各种湍流输运，可以用统一的数值解法处理两相流问题，模拟结果能够 给出颗粒相空间分布的详尽信息,并且易于和试验数据对照。发展拟流体模型的关键是颗粒相湍流模型。到目前为止，建立双流体模型颗粒相基本方程及其紊流模型主要有两种途径:一是用雷诺平均建立颗粒连续、动量、能量颗粒素流动能和雷诺应力方程，并采用类似单项紊流模型的方法加以封闭；另一途径是，应用概率密度函 数建立概率密度（Pd)输运方程。进行统计平均或质量加权平均,建立紊流两相流动的 连续、动量、能量雷诺应力或紊流动能方程，并加以封闭，构成两相紊流模型。1970年，Hiz-he基于局部各相同性紊流脉动，提出以颗粒追随流体脉动概念

34、为基础的代数模型(HinzeThn公式），简称。模型。把气相k一占模型与爿.模型合称为一占一一。模型29J。按此模型颗粒的湍动能k。恒小于流体的湍动能，且颗粒尺寸越大， k.k越小。对此，198年，LZou（周力行）提出了建立颗粒湍动能七。的方程,构成 一占一k。模型pw.该模型认为颗粒湍流动能不仅取决于当地气体湍流的作用而且取决于其 本身的对流，剪切产生和扩散，因此一定情况下在流场的某些区域颗粒脉动比气体的强.由 于两相紊流是各向异性的，为此周力行等人从两相各自的瞬时及时均方程出发，导出了两相各自的雷诺应力输运方程和两相脉动速度关联的输运方程而建立了统一二阶矩（UM)模型。对于第二种途径，李

35、勇、周力行等通过研究紊流两相流动的瞬时量与脉动的Pd瀚运方 程,建立了两相紊流k-只矿模型p2J，该模型对流体相紊流采用一占模型而用df模 型来计算出颗粒相各方程的脉动关联相；而后又提出了气相雷诺应力方程模型和颗粒概率密度输运方程相结合而构成的一蹦两相湍流模型3:此外文献口4应用Pd蹦运方程，建立双流模型基本方程把联合概率密度函数(Pdf)辅运方程用于气相和颗粒相模拟两相紊流 流动，封闭两相紊流模型。44气相湍流燃烧（反应)模型综述燃烧现象是物理与化学过程相互作用的结果。由于燃烧反应的速率不仅与化学动力学因 素有关,而且还与反应体系内质量和能量的输运有密切关系.这种输运取决于系统的宏观运 动(

36、可燃混气的流动）以及系统内微观运动(分子热运动）所引起的热传导和扩散现象。在 实际燃烧过程中，燃烧速率往往受控于某些物理因素。因为在燃烧系统中燃料和氧化剂的浓 度足够大，气温温度又很高若混合均匀，则化学反应速率可能很大,但如果混合不均匀， 混合或扩散速率相当慢时,则燃烧速率实际上取决于扩散速率。这样.许多物理因素如混合 气体的速度分布、流动情况、燃烧设备形状和尺寸、紊流与局部扰动情况、分子扩散与涡旋扩散情况以及系统内热量分布等。都对燃烧速率有很大影响。对于混合气扩散速度比化学反 应速率低得多的燃烧过程。称为扩散燃烧。相反对于混合均匀的气体.如果化学反应速率 比扩散速率小得多的燃烧过程.称之为动

37、力燃烧.若化学反应速率与扩散速率相当，则可称 之为扩散动力燃烧。扩散燃烧主要受扩散、流动和其他物理混合过程的控制.而动力燃烧主 要受化学动力因素支配扩散动力燃烧则是化学动力学和物理混合因素同时起支配作用口”。一般来说，凡是燃料和氧化剂未预先混合的燃烧过程称为扩散燃烧,该火焰称为扩散火浙江大学硕士学位论文周平 第一章绪论焰。Saidig于1971年提出了预估紊流扩散火焰的k-g模型,后来演变为k一占一g模 型圳,其主要思想是,采用紊流k-s模型来描述紊流特性,燃烧按简单化学反应系统和快 速反应假设进行简化，建立厂和g,混合分数，是脉动的t是一个随机变量,可采用一个适当的概率密度函数(f)来描述厂

38、的脉动性质，如果P(f）已知,则各标量的时均值都可求出.一般P厂)是未知的,确定方法有两种:一是直接建立以P（m）输运方程，找出矽的分布形式,即概率密度分布函数输运方程模型田;另一种方法是根据对紊流脉动的设想 事先给出Pd分布形式,即简单概率函数模型或简化的pd模型。典型的pdf函数有：城 墙式(Batement）或方齿形分布、截尾正态分布概率密度分布函数（ClipedNormal Disiuion）、Bta函数形式、联合概率密度函数等。燃料与氧化剂在进入火焰区之前已经均匀混合的火焰。称之为预混火焰.关于紊流预混 燃烧反应速率计算模型，最早是Spalding提出的旋涡破碎模型唧（Edd Bra

39、k-Up Model,简 称EBU模型)，其基本思想是：把湍流燃烧区考虑成末燃气徽团和已燃气微团的混合物；化 学反应在这两种微团的交界面上发生；认为化学反应速率取决于未燃气微团在湍流作用下破 碎成更小微团的速率；认为破碎速率与湍流脉动动能衰变的速率成正比。由于该模型微考虑温度对化学反应速率的影响，Mason等人又提出TEBUrrhnius模型坩。17年,agnssen提出了涡团耗散模型I*0（EddyDisipainoel，简称M模型），其基本思想是：当气流 涡团因耗散而变小时，分子之间碰撞机会增多,反应才容易进行并迅速完成，故化学反应速 率很大程度上受紊流影响，而且反应速率还取决于涡团中包含

40、燃料、氧化剂和产物中浓度值 最小的一个。该模型反应速率取决于紊流脉动衰变速率“七，并能自动选择成分来控制速率. 因此该模型既能用于预混火焰，也能用于扩散火焰。但它只适用于快速反应的化学反应，为 此,9年anusen等人有提出了修正的ED模型,其主要思想是：当反应物气团能以 分子尺度进行混合，并且温度又较高时，化学反应才有可能发生。在紊流中，分子尺度混合 与紊流中最小尺度有关,这个混台过程与紊流能量的消耗过程相似。在紊流中涡团耗散不是 均匀分布而是集中在高度拉伸变形区域，在这些耗散区内充满了其特征尺寸与Kolmgor微尺度同数量级的小涡团，称这区域为细微结构(FieStructes）。1981年

41、，agnussn又提出了涡团耗散概念模型1421(EdyDssipaionConcepMol)。该模型的出发点是化学反 应由两个区组成：一是在紊流微细结构中分子间很容易碰撞,发生化学反应，此时化学反应受化学动力控制；二是在微细结构的周围较大涡团区域内，含有不同成分的涡团必须混合后 才能进行化学反应，混合所需时间比化学反应时间长，因此受混合速率控制.EDC模型较适 用于那些反应速率较低的化学反应。在紊流燃烧模型研究中，为了使基本方程得到封闭，需要模拟平均化学反应速率.还可 采用紊流燃烧关联矩模型m】(平均反应速率输运方程模型）、概率密度函数反应速率模型ml 或二阶矩-概率密度模型m1等有限反应速

42、率模型.45辐射换热模型综述辐射换热过程有其特殊性.燃烧装置中火焰的热辐射与火焰中介质的温度以及介质 的辐射吸收、散射能力有关，而介质的辐射、散射能力与辐射波长有关。燃烧空间中任意一 点对空间中其它任一点都有辐射作用，燃烧装置的壁面通常对辐射有反射作用。炉膛内的温度很高，炉内高温火焰和水冷壁之间的传热以辐射换热为主。占总换热量的0%左右。常用的计算辐射换热的方法有热流法（热通量法)、区域法、蒙特卡洛法(又称 概率模拟法)、离散传播方法；还有用数论方法模拟辐射换热。自从19年Shuer提,mx2通量模型【46】以来，目前已经发展到能应用于多维空间的四通 量、六通量模型【47l。浙江大学硕士学位论

43、文周平第一章绪论区域法通常又称otle法 J。其主要原理是将求解空间划分为若干体积区域，将四周壁 面划分为若干面积区，并假定任一区域(体积区域或面积区域）的温度和辐射性质都均匀一 致。区域法在原理上是好的。但由于计算工作量巨大,实际上只能把燃烧室分成较少的大区 域.目前区域法未能广泛应用于燃烧室传热过程的工程计算中。蒙特卡洛法(MonteCaro Mehd)也称概率模拟法，这是一种求定积分的数学方法,具 体应用到辐射换热的计算中I，便是对各微元体的吸收情况作概率模拟计算,而不需对各能束射线沿途计算辐射能量被逐步地吸收和透过的情况。其基本思想是：把燃烧室分成若干 个气体微元和面积微元，假定各区域

44、所具有的辐射能由很多条能束组成，每条能束具有一定的能量，依次跟踪每一条能束，直到被完全吸收为止，若跟踪足够多的能束历程,就可得到 具有统计意义的结果，从而求得各微元间的辐射换热量。蒙特卡洛法应用较多。我国也在这 方面作了大量的和卓有成效的研究工作，取得了可喜成绩。这种方法的计算结果有一定的统 计误差。离散传播法(scete Trasfe Metho ）的主要思想是考虑边界网格面作为辐射的吸 收和发射源，将边界网格面上向上半球空间发射的辐射能离散成有限束能束,这些能柬穿过 内部网格被介质吸收和散射后.到达另外的边界面上被吸收。在各边界网格面上进出的辐射 能达到平衡。数论法其主要思想为：求解高维多

45、重积分,首先用数论方法构造出高维立方体中的一般 分布点集贯,然后用事先选定的点集贯上的被积函数织构成的单和来逼近多重积分。在被积函数适合一定的条件时，逼近的误差主阶是与积分的维数无关的.所谓一致分布点集贯,又 称数论网格,是用数论方法构造的具有一定规律的点的集合。其实质就是把一个复杂的多重积分化为一个高维立方体上的多重积分.数论方法能够消除蒙特卡洛法难以克服的随机误 差，此外，由于数论法能用较少的能柬数给出结果，蒙特卡洛法需要更多的能束数,所以单从计算方法本身来讲,数论方法也有其自身的优越性.1.。污染物排放的数值模拟综述本文主要考虑NO的数值模拟.其主要模型可分为四类pq: (）简化动力学模

46、型。此类模型没有考虑紊流影响,适用于层流燃烧.(2）Arheni-E模型。此类模型简单，并考虑紊流影响，目前在工程问题上应用较多,但计算精度不够理想，需加以改进。（3)f模型。目前研究较多是二维联合概率密度函数，此类模型有些假设不太合理， 如用两个pdf的乘积（尸（r）(m）或P口）尸(厂）来代替联台概率密度这意味着温度与 浓度无关；在局部平衡pdf模型中假设反应度石与紊流脉动无关，这于实际情况不相符,计算工作量大。(4)关联矩模型。用二阶矩封闭法直接求时均反应速率的模型，比较简单，但引入温度指数项的级数展开近似,误差较大。为了避免这些不足之处。统一二阶矩模型直接导出反应 速率系数和浓度脉动关

47、联矩。并用统一二阶矩输运方程求解，避免了pal积分，这样较为合理,计算工作量也可能比df模型小些。1.5本文研究的主要内容及其意义吉林石化新建75讹高压高温水煤浆锅炉的设计在国内尚属首次,通过锅炉冷态试验和热态燃烧调整试验证明,锅炉各项性能基本达到了预期的目标。为了优化锅炉运行，针对四 角切圆燃烧锅炉存在的问题,本文通过数值模拟，对炉内冷态空气动力场、燃烧过程、NO排放进行了详细的模拟.主要内容包括:浙江大学硕士学位论文周平第一章绪论（1）锅炉冷态空气动力场试验： （2）试验工况配风下，锅炉冷态空气动力场的数值模拟； (3）通过数值模拟，研究改变锅炉配风大小，配风方式对炉内空气动力场特性的影响

48、； （）锅炉热态燃烧调整试验； (）锅炉8负荷下全烧水煤浆炉内燃烧、传热的数值模拟；（6)污染物o排放的数值模拟：(7)将部分数值计算结果与试验结果进行对比，以达到最终指导实践的目的,完成计算 机辅助优化数值试验在锅炉中的应用;(8）对模拟炉内过程的方法和模型进行概括。通过以上试验和数值研究，有助于了解新建7t四角切圆燃烧水煤浆锅炉空气动力场 特性、燃烧特性，用于指导四角切圆燃烧水煤浆锅炉的设计、运行方面存在的问题.本课题通过应用锅炉冷态现场试验、热态现场试验和数值模拟相结合的方法,目的明确，为锅炉设计和运行提供指导。16本章小结本章首先概述了吉林石化新建75h高温高压水煤浆锅炉工程项耳的意义

49、。指出了本课题的由来，总结了四角切圆燃烧锅炉的主要优缺点以及采用数值模拟方法优化锅炉设计及运 行的利与弊，然后对涉及计算流体力学、燃烧学、传热学、气固两相流，化学反应动力学等诸多领域的各种数学模型进行了综述，包括气相湍流流动模型、两相流动模型、辐射换热模 型、污染物排放模型等。最后，概括了本文的主要研究内容及其意义.参考文献：【】钱力庚,330W电站对冲锅炉炉内过程数值模拟和实验研究及四角锅炉变负荷研究及炉内过程通用程 序的设计与研究,浙江大学博士学位论文，2001年月【2l岑可法等.燃烧流体力学。水利电力出版社，1991年月【1李文艳，电站锅炉煤粉燃烧过程及结渣的数值模拟，华北电力大学博士学

50、位论文，02年12月f4】 V.帕坦卡，传热和流动的数值方法,安徽科学技术出版社,18年2月。【5 J PSpalrand SAlram，A ne。equtinrbulenmoel for rodynamc flw，ecnicalReprtIA2039。ncan Isitut fAronaut anronautics,19【6lB.Laudrand D.Blding，Lectures i MathematiaodeIs ofmbnlenc。adeicPr，LondonEnland,97f7J DChoudhury Inrutio to th Rermizutio Goup Method and

51、TubleceModeng.FluentIne Tehnil Meorndum T-107，199【8】-HShih,WW.i,AShai,ZYang,dJhu Nw k一占Eddy-Viscsity Model for Hgh Reynold Number Turuent FosMdel Delpmnnd alidatio.ots luis，240）:2278,15【9】D CiloxTudleneModelingorCFD.DCWIndutres，Inc，Laaad，llfoma，19【1】 M。Gibson andB ELundeGod Effecto Presur Fltutionsn

52、th Alanosheri Bonday LayerJ Flid Mch。86：491511978【l】B Eaude,SeconMoment Cosue:PretanFture？InmrJeatFlui Flow.10H）:282。300，t989，【12】B ELaundr，.JRcc，and WRodi。Progrein th Developmet o ynolds。Strs TurbulenceCsrc .luid Mc。68（3）:3-6,pril975【3】 O.Ha，rblen，Mcraw-Hll blishng Co,New Yrk。197【4 rafLT.，Launde,B

53、E。and Suga.，Development ad Application of a CubcEdViscosit MIf Ibulece，int JHea FluidFlow，01 7o.2,P108。1996【15】Shi，t H，Zhu，Jand Lo,W，Modeln ofTubulnSwirlngFlos,NASATM 13119971】ih，TH。,urbulecedlDeopenttICM，IAA 9833.198【7Rodi，。，A New gbric Relaion fr Calculatingthe Rnols tsas,ZAnge。Math Mech.V0.5浙江大学硕

54、士学位论文周平第一章绪论p210，975【8】Leonrd,A，Eergy aade in Larg Eddy Smlt fTubulence luid los，avGophys,l 68， 37，174【1】Smagoinsky，J，GenelCiculmon Eprimnts with te Pimitive Equmion I,he asc Experimn,eWether Rv，ol 91, 3,P 99，1963【201 im，。W，A Nw ynaic OneEation ubgridScale Moe forLage Edy Simulaions，AIAA 9035，1995【I

55、】arghini,F nd Pimelli,，Scaesimub鲥 Scss ModefoLare-Eddy iulat，IA9805998.【22】Ravikah，VIL ad Plecher,R H,rrge dd Smuifthe urulent FlowPastA Backard Facing Sep。IAA 200004,2000【3】周力行，湍流两相流动与燃烧的数值模拟，清华大学出版社，91年8月。【24】Sldg，DB,Basiuatio ofFlud Mcaic，H/66,art 10，1976【25soo,sL，FluidDynas ofMliphase Systems,Ba

56、ael，1967。【26】Cwe, t，The Partice-Souron-Celodcl岛 asrple Flow，SM75-WA7。2,1975.127】Smth，P。J,eal.PaperWSS-7,Wetentatscto,obstionIn。tare，1990【28】岑可法，樊建人，PaiuteSciece and Tecnoogy，V01.6,No，198【29】Zou。L ，They a NuroeriM Mdelin of TurblentGasPatile Flows andComustin，icePres（Beijng）an R Pres（Fovda，S），hap，pp5

57、6163。993。【30】hon， xan Huang，X Qienceinhin，。33,N,198。【31Zho,X。,iao，C。 ad hnt,PmenterSyru。ONum。Mth For Mlh.low,LakeTahoe,NV，PP。3073，1994. （】李勇、周力行,k一占一Pdf两相湍流模型和台阶喉方气粒两相流动的模拟,工程热物理学 报，Vo 17，No.2,P234,1【33】周力行、李勇，旋流两相流动的DSM一尸咿两相湍流模型。工程热物理学报。V0.17,No.2,P.22，1999.【4】Zhu，Mad Bray,K。C。，PDF Trt Equation fo

58、ToPhase eactve los ad Spry，Comb adFlame，V01。22，P 3138，20.【3】赵坚杼,燃烧的数值模拟，科学出版社,002年6月。【36】Saiding。DB，onentraio ltuations in a ound urlat FreeJt,Chcal Eineer Scince.V01.6，P9510，197137】r,KN。dMoss，I。，A UBiffdSatisica Modelof th rtixad Tfoln FlametAsrnautc，Vl ，No-4，P9。19938】Spding,D.，Mixin and Cemica Reacinn Sendy ofned Tubuent Flames，13h Sypoim(nt.）o Cmbstion，p649,71.【3】Mason，Had Sadin，DB。，Peiin of eaco Rates Tbulen Premixe BndyLaer Flo，Pmc Europen Symp on ombtion，01，1【4】Magu