煤燃烧中NOx生成的模拟课件

,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,6 煤燃烧中NOx 生成的模拟,氮的氧化物,NOx,主要是,NO和NO,2,，其次是,N,2,O，N,2,O,3,，,N,2,O,4,和N,2,O,5,1 燃烧过程中,NOx,的生成机理,NOx,生成的三种主要方式,热力,NOx,它是助燃空气中的,N,2,在高温下氧化而生成的氮氧化物，其生成机理最先是原苏联科学家,Zeldovich,等提出的,燃料,NOx,，系指燃料中的有机氮化合物在燃烧过程中氧化生成的氮氧化物,快速,NOx,，指碳化氢系燃料在燃烧时分解所产生的中间产物和,N,2,反应生成的氮氧化物,上述三种氮氧化物的组成随燃料含氮量不同而有差别。对于燃煤，通常燃料,NOx,占70%-85%，热力,NOx,占15%-25%，其余为少量的快速,NOx,热力NOx,高温下空气中的N,2,氧化的产物，其主要反应如下,ON,2,NON,（136）,NO,2,NOO,(137),NOHNOH,(138),以上反应的活化能很大，且生成速度与燃烧温度的关系密切,当燃烧温度低于1500时，热力NOx生成量极少,控制热力NOx 的关键在于降低燃烧温度水平，避免局部高温，同时降低氧气浓度,快速NOx,碳氢化合物燃烧过程中分解的CH，CH,2,和C,2,等基团破坏了空气中的N,2,分子键，并经反应生成HCN，NH和N等原子基团，它们再与O，OH等基团反应生成NO,快速NOx只有在富燃的情况下，即碳氢化合物较多，氧浓度相对较低时才发生,燃料NOx,影响因素多，反应复杂，因而至今对其反应机理的认识还不够完善,有文献认为，从燃料N向NO的转换是由两个互相竞争的过程所决定的。这两个过程是：由燃料N在高温下分解生成含有N原子的中间生成物I（主要是N、CN、HCN、和NHi等化合物），然后I和含有氧原子的反应物R（如O、OH、O,2,等）反应生成NO，或者和NO反应而使之还原为N,2,，即,IRNO,燃料N,INON,2,DeSoete,等人认为燃料N高温分解主要生成HCN，HCN氧化后生成NO，同时NO又被HCN还原生成N,2,，即,char,O,2,NO N,2,燃料N HCN,NON,2,煤粉燃烧之燃料NOx的生成机理更为复杂,既存在挥发分的单相燃烧，又有焦碳的多相燃烧，因而这时的燃料NOx应包括挥发分中的N生成的NO和残留在焦碳中的N生成的NO这两个部分,挥发分中的氮化合物主要有HCN和NH,3,，两者的比例不仅取决于煤的挥发分，而且与N和碳氢化合物的结合状态等化学性质有关,在典型的煤粉燃烧条件下，挥发分NOx约占燃料NOx的60-80%,相比挥发分NOx，焦碳NOx的生成机理更为复杂。焦碳NOx的生成速率与焦碳中的N含量、氧浓度和温度、煤颗粒孔隙结构及颗粒反应表面积等因素有关。同时，焦碳表面和CO等还原性气体对已生成的NO会产生还原分解作用，使燃料NOx减少,2.燃烧过程中NOx生成的模拟,1)基元反应模型,NOx在燃烧过程中生成复杂，有的研究者就忽略实际燃烧过程的特征，强调NOx生成的平衡反应，考虑体系内所有可能的基元反应和反应物质，按平衡方程求解,忽略了实际的湍流燃烧过程的特征，这样的处理尚未让众多的燃烧工作者接受,2)NOx生成的PDF输运方程模型,燃烧过程中NOx的生成不仅和本身的详细化学反应机理有关，而且和湍流流场及其之间的相互作用有关。诸如湍流场、温度场、燃料种类、氧量和辐射换热等，特别是化学反应机理的非线性和非平衡性，都会极大地影响污染物的最终排放量,通常的湍流模型在模拟化学反应时，一般是假设化学反应速率是线性的或极快的或极慢的，以获得平均反应速率的简化表达式,实际燃烧过程中，化学反应是非线性的、其反应速率是有限的，其中NOx的生成和详细的化学反应机理密切相关，难于用平均的总包反应速率来描述,PDF输运方程模型可以把与湍流输运及化学反应速率有关的项都用封闭的形式表达，无需模拟，任何复杂的化学反应机理都可精确计算，因而它在污染物生成过程的模拟中得到重视和应用,求解PDF输运方程所用的Monte Carlo统计方法在流体颗粒数选取较大时需要很大的计算工作量和存储空间,k-,湍流模型与化学热力学参数的标量联合PDF方程模型,考虑到计算工作量，其计算过程中选用的随机样本数并不大，这限制了PDF输运方程模型和Monte Carlo模拟对实际更为复杂的湍流反应流的计算,为了兼顾PDF输运方程模型的优点，克服其以上不足，目前实用的主要是给定PDF的简化PDF模型和简化PDF与局部瞬时平衡相结合的模型.前者仅适用于双组分反应的简单情况，后者则可适用于复杂的扩散燃烧,简化PDF与局部瞬时平衡相结合的模型(同轴射流射入突扩燃烧室),同轴射流射入突扩燃烧室时，其中心一次风射流为燃料与空气，环缝射流为二次风。混合物分数可定义为：,（139）,其中m,p,为来自一次风的流体质量，m,s,是来自二次风的流体质量。混合物分数f表示任一时刻在任一位置上的混合程度，即局部瞬时当量比。任何其它守恒标量可表示为：,f,p+(1-f),s,(140),其中,p,与,s,分别是一次风及二次风中的,值。由于在一般的反应中任何元素既不可能产生也不以消灭，所以元素k的局部瞬时质量分数“b,k,”也是一个恒标量。,当各元素的扩散系数都相等时，可有,b,k,=f b,kp,+(1-f)b,ks,(141),对于无辐射、无导热损失的绝热系统，气体焓也是守恒标量。在导热系数和元素扩散系数相等时，有,h=f h,p,+(1-f)h,s,(142),混合物分数的时均值 和脉动均方值g可由5节介绍的方法求得。从已求得的 ，g以及假设的PDF我们可以找到每个当地位置上的PDF。若对任何标量,，,(f)函数是已知的，则可得到,的时均值：,（143）,（144）,其中a,p,和a,s,是p(f)在f=1和f=0时的值；,p,和,s,分别是一次风和二次风中的,值。使用局部瞬时平衡的概念，能够找到瞬时温度、密度及组分浓度随焓和元素组分而变化的函数关系(当压力变化甚小时)：,T=T(b,k,h),=,(b,k,h),Ys=Ys(b,k,h),使用化学平衡计算可以找到上述函数关系。,对于绝热体系，焓与元素质量分数仅是混合物分数的函数，即,b,k,=b,k,(f),h=h(f),则有,T=T(f),=,(f),Ys=Ys(f),所以对这种情况，我们可以直接使用p(f)获得时均温度和时均气体密度和时均组分质量分数 和 ；,对非绝热系统，需求解能量方程找出焓h。这时有,T=T(f,h),=,(f,h),Ys=Ys(f,h),其时均值 、应决定于,（144a）,但是，联合PDFp(f,h)是难于预先给定的。有一种替代方法是假设,h=hf+hr,hf=hf(f),并令,p(f,h)=p(f),由此可得,(144),3),DeSoete,模型,在该模型中，所有中间组分都设为HCN，煤中N的释放速率RcoalHCN正比于煤粉热解及煤焦燃烧时的质量衰减速率Sp，即,（145）,式中f,N,为煤中氮的质量分量；M,HCN,、M,N2,分别为HCN及N,2,的分子量（下面各式类同）。,HCN氧化后生成NO，其反应速率为：,(146),式中,g,、T,g,分别为混合气的密度和温度；Yi表示组分i的质量分数；b为氧的反应级数。由DeSoete模型给定，可表达为氧浓度的分段线性函数。,同时，所生成的NO被HCN还原生成N,2,，其反应速率为,(147),对热力NOx可采用Zeldovich机理计算，其反应速率为：,（148）,根据上述NOx生成的机理模型，要考虑煤燃烧过程中NOx的生成，只需在数学模型中考虑两种组分的质量分数Y,HCN,和Y,NO,的输运方程，而其形式与气相标量（如混合分数f）输运方程一样，只是其源项分别由上述NOx 生成机理给定，即,(149),(150),其它化学组分（如O,2,、N,2,等）的质量分数可由混合分数f的值根据所采用的挥发分燃烧模型而确定。,4)扩展的DeSoete模型,DeSoete模型中仅考虑了燃料NOx里面的挥发分NOx中的一部分（即经由HCN产生的NOx）以及热力NOx,事实上，挥发分中还会释放NH,3,等组分，它们对NOx 的生成也有贡献,焦炭NOx 在很多情况下也占有总的NOx中的相当比例,再考虑近年来研究较多的燃料分级(Reburning)减少NOx排放的技术，则煤燃烧过程中NOx生成与排放的机理可用图5示意图表示（忽略快速NOx）,图5 煤燃烧过程中NOx 生成模型,在上述模型中，挥发分释放组分HCN及NH,3,的速率为：,（151）,（152）,以上两式中，,为挥发分中释放的组分HCN的质量分数；S,devo,l,为煤的热解速率（kg.m,-3,.s,-1,）。,在贫燃条件下，过程1和过程3生成与分解NO的速率为：,(153),(154),式中M,gas,为混合气的分子量；其余符号参见前述DeSoete模型中所述。,在富燃条件下，过程1和过程3生成与分解NO的速率为：,（155）,(156),过程2生成NO的速率为：,（157）,过程4消耗NO的速率为：,（158）,过程5系采用燃料分级技术降低NOx 生成的机制。此处假定再燃燃料为主要组分是CH,4,的天然气，该过程分解NO的速率为：,(159),式中,(160),其中 为CH4燃烧过程的化学当量比。,过程6是焦炭NOx生成的过程。在计算该过程中NO的生成速率时，假定焦炭NO产生于颗粒表面，其值正比于焦炭的燃烧速率（式161b）；同时，生成于颗粒表面的NO会发生还原反应（式161c）；颗粒表面NO的质量源就可以通过上述两式得到，进而求得焦炭表面的NO浓度（式161a）。总的NO生成率根据焦炭表面的NO浓度及其从颗粒表面向周围环境的扩散流由下式计算：,（161）,式中焦炭表面的NO浓度Y,NOsup,如下计算：,(161a),焦炭NO的产生率为,（161b）,其中S,char,为焦炭的燃烧速率。,焦炭表面NO的还原速率为,（161c）,其中A,E,为焦炭颗粒外表面积；m,char,为焦炭质量；VOL为其体积；P为总压力。,而焦炭表面NO向周围环境扩散系数为,（161d）,其中Sh为Sherwood数，取为2；D,NO,gas,为NO在燃烧产物中的双元扩散系数；d,p,为颗粒直径。,热力NO的生成（过程7）速率则由下式确定：,（162）,式中Y,0,为与Y,O2,及温度T,g,有关的量；k,i,为对应各反应的化学动力学常数。,有了上述各个子过程的反应速率，按照前述DeSoete模型中类似的方法，即可求解各组分（NO、HCN、NH,3,等）控制方程的源项，进而求得其空间分布。,由于NOx生成与分解反应的复杂性，目前无论是采用,DeSoete,模型，还是扩展的,DeSoete,模型，在计算燃烧过程NOx 生成特性时，基本上都是采用所谓的后置处理方法（,Post-Processor,）,即认为其化学反应耗氧、生成热及气相参数不影响主燃烧反应,事实上,它们之间是有强烈的相互作用的。但如果考虑这一相互作用，则整个炉内气固两相流动、传热和煤粉燃烧的计算都将作相应的改变,这一过程非常复杂，加之采用后置处理方法所求得的NOx生成特性基本满足工程需要，所以至今仍然以这种简单方法为主,