300MW燃煤锅炉SNCR过程的数值模拟.pdf
第 4 2卷第 4期 2 0 1 2年 7月 东 南 大 学 学 报(自 然 科 学 版 ) J O U R N A LO FS O U T H E A S TU N I V E R S I T Y( N a t u r a l S c i e n c eE d i t i o n ) V o l . 4 2 N o . 4 J u l y2 0 1 2 d o i 1 0 . 3 9 6 9 / j . i s s n . 1 0 0 1- 0 5 0 5 . 2 0 1 2 . 0 4 . 0 2 2 3 0 0MW 燃煤锅炉 S N C R过程的数值模拟 姜 敏 金保癉 周英贵 王晓佳 ( 东南大学能源热转换及其过程测控教育部重点实验室,南京 2 1 0 0 9 6 ) 摘要利用 F L U E N T软件对 3 0 0MW 燃煤锅炉的选择性非催化还原( S N C R ) 过程进行了模拟. 重点分析尿素溶液的喷射位置、 液滴喷射速度、 雾化平均粒径对 S N C R过程脱硝效率以及漏氨 量的影响, 从而得到实际运行过程中最佳运行条件. 结果显示, 增加喷射速度能加快蒸发过程; 雾 化粒径越大越有利于还原剂的平均分布, 但也会延长蒸发时间; 在第 1 , 2层同时喷射尿素溶液, 并控制雾化平均粒径为 0 3m m和液滴喷射速度为 3 0m/ s 时, 能够获得较好的脱硝效果. 将不 同氨氮摩尔比模拟结果与实际过程进行对比, 结果表明两者能够较好地吻合. 关键词数值模拟; 选择性非催化还原; 雾化; 氮氧化物 中图分类号X 5 1 1 文献标志码A 文章编号 1 0 0 1- 0 5 0 5 ( 2 0 1 2 ) 0 4 0 6 9 1 0 5 N u me r i c a l s i mu l a t i o no f S N C Ri na3 0 0MW p u l v e r i z e d c o a l b o i l e r J i a n gMi n J i nB a o s h e n g Z h o uY i n g g u i Wa n gX i a o j i a ( K e yL a b o r a t o r yo f E n e r g yT h e r m a l C o n v e r s i o na n dC o n t r o l o f Mi n i s t r yo f E d u c a t i o n ,S o u t h e a s t U n i v e r s i t y ,N a n j i n g2 1 0 0 9 6 ,C h i n a ) A b s t r a c t Mo d e l i n ga n dn u m e r i c a l s i m u l a t i o nf o rS N C R( s e l e c t i v en o n c a t a l y t i cr e d u c t i o n )i na 3 0 0MW p u l v e r i z e d c o a l b o i l e r w e r e c o n d u c t e db a s e do nt h e F L U E N Ts o f t w a r e .Ma i nf a c t o r s i n c l u d i n gt h e i n j e c t i o np o s i t i o no f u r e a s o l u t i o n ,d r o p l e t v e l o c i t y ,d r o p l e t s i z e a n dt h e i r i n f l u e n c e o nN O x r e m o v a l e f f i c i e n c ya n da m m o n i as l i pw e r ea n a l y z e d ,t h u st h eb e s t o p e r a t i o n a l c o n d i t i o n sw e r eo b t a i n e d .T h e r e s u l t s i n d i c a t e t h a t i n c r e a s i n gt h e i n j e c t i o ns p e e dc a na c c e l e r a t e t h e e v a p o r a t i o np r o c e s s . T h eg r e a t e r t h e d r o p l e t s i z e i s , t h e b e t t e r t h e r e d u c t a n t d i s t r i b u t e s , a n dt h e l o n g e r t h e e v a p o r a t i o nt i m e i s .I n j e c t i n gu r e a s o l u t i o no nb o t ht h ef i r s t a n ds e c o n df l o o r w i t ham e a na t o m i z i n gp a r t i c l es i z eo f 0 3m ma n dav e l o c i t yo f 3 0m/ s s h o w s a ni d e a l r e s u l t .S i m u l a t i o na n dm e a s u r e dv a l u e s u n d e r d i f f e r e n t m o l a r r a t i o s o f N H 3/ N Oa r ec o m p a r e d ,w h i c hs h o w s ag o o da g r e e m e n t . K e yw o r d s n u m e r i c a l s i m u l a t i o n ;s e l e c t i v en o n c a t a l y t i cr e d u c t i o n ;a t o m i z i n g ;n i t r i co x i d e 收稿日期 2 0 1 2 0 3 1 2 . 作者简介姜敏( 1 9 8 8 ) , 女, 硕士生; 金保癉( 联系人) , 男, 教授, 博士生导师, b s j i n @s e u . e d u . c n . 基金项目科技部中欧国际合作资助项目( 2 0 1 0 D F A 6 1 9 6 0 ) . 引文格式姜敏, 金保癉, 周英贵, 等. 3 0 0MW 燃煤锅炉 S N C R过程的数值模拟[ J ] . 东南大学学报 自然科学版, 2 0 1 2 , 4 2 ( 4 ) 6 9 1 6 9 5 . [ d o i 1 0 . 3 9 6 9 / j . i s s n . 1 0 0 1- 0 5 0 5 . 2 0 1 2 . 0 4 . 0 2 2 ] 燃煤电厂的 N O x排放是重要的污染源. 针对 N O x的排放, 目前电站广泛采用的控制方法主要 有选择性非催化还原技术( S N C R ) 和选择性催化 还原技术( S C R ) 等. S N C R技术相对于 S C R技术, 投资和运行成本低. 通过在炉膛中喷入还原剂的方 法, 有利于旧电厂的改造, 比较适合中国国情. 但 S N C R技术的还原效果受到多方面的制约, 如反应温度、 混合程度、 还原剂的停留时间、 氨氮摩 尔比、 烟气组分、 N O x的初始浓度等 [ 1 ]. 在实验室条 件下, 能达到8 0 % ~ 9 0 %的还原率, 但是由于还原剂 与 N O x混合不理想等原因, 在实际的运行过程中, 脱硝率一般不会超过5 0 %. S N C R过程非常复杂, 因 实验过程中无法测出炉内的温度场及浓度场, 给研 究工作带来了很多不便, 而数值模拟可以详细地反 映出这一过程. N g u y e n等[ 23 ]以 F L U E N T软件为平 台, 利用简化的还原反应对管式炉中的 S N C R过程 进行了模拟, 表明数值模拟的方法具有一定的指导 性. 但由于规模较小, 其结果对锅炉实际运行帮助有 限. 本文在前人研究的基础上, 使用 F L U E N T软件 对大型3 0 0 MW 燃煤锅炉进行三维全尺寸模拟, 针 h t t p / / j o u r n a l . s e u . e d u . c n 对燃烧温度、 N O x生成及 S N C R过程进行预测, 研究 尿素溶液喷射位置、 喷射速度、 雾化粒径以及摩尔比 对脱硝效果的影响, 确定其最佳运行参数. 1 模型及计算方法 1 1 数学模型 炉膛中的燃烧模拟已较为成熟, 采用 R N Gκ ε 模型预测炉膛内烟气的流动状态, 该模型相对于标 准 κ ε模型考虑了旋流的影响. 采用离散相模型处 理煤粉, 模拟煤粉运动采用拉格朗日随机轨道模 型. 由于炉膛内温度较高, 辐射传热占主要地位, 采 用 P 1模型模拟辐射传热. 焦炭热解采用双竞争反 应热解模型, 挥发分燃烧采用混合分数概率密度函 数 ( p r o b a b i l i t yd e n s i t yf u n c t i o n , P D F ) 法, 焦炭燃烧 采用动力 扩散模型. 利用 S I MP L E算法, 采用一阶 迎风格式对控制方程进行求解. 由于燃烧过程中产 生的 N O x浓度非常小, 对煤粉的燃烧过程几乎没 有影响, 所以预测 N O x生成采用燃烧后处理方法. 由于快速型 N O x的生成量较少, 所以本文只考虑 燃料型及热力型 N O x的生成. 其中, 假设 9 0 %的挥 发分氮转变为 H C N , 其余则转变为 N H 3, 而焦炭氮 直接氧化为 N O [ 4 ]. 选用尿素作为 S N C R过程的还原剂, 质量分 数为 5 0 %的尿素溶液与稀释水混合后经雾化喷入 炉膛中, 经过水的蒸发、 尿素分解及 H N C O的水解 3个步骤后与 N O x反应 [ 56 ] 1 )水的蒸发 C O ( N H 2)2( a q →)C O ( N H2) 2( s 或 l )+ H2O ( g ) ( 1 ) 2 )尿素分解 C O ( N H 2)2( s 或 l →)N H3( g )+ H N C O ( g ) ( 2 ) 3 )H N C O的水解 H N C O ( g )+ H 2O ( g →)N H3( g )+ C O2( g ) ( 3 ) 湍流与化学反应的耦合选用涡耗散概念 ( E D C ) 模型. 本文中尿素的分解动力学参数采用 由 A o k i 等[ 7 ]提出的分解反应模型, 该化学反应动 力学模型与尿素瞬间分解的假设相比, 更接近于实 验值[ 8 ]. 对于 N O x与还原剂反应的过程, B r o u w e r 等[ 9 ]提出了简化的 7步反应模型, N g u y e n等[ 3 ]又 对该模型进行了改进. 本文中采用的尿素分解及简 化还原反应的动力学参数如表 1所示. 表 1 尿素分解及简化还原反应的动力学参数 序号反应方程指前因子 A / ( m3s - 1K- bk m o l- 1) 温度指数 b活化能 Ea/ ( J k m o l - 1) 1C O ( N H2)→ 2 N H3+ H N C O 1 2 7 1 0 4 06 5 0 1 0 7 2C O ( N H2)2+ H2→O2 N H3+ C O26 1 3 1 0 1 0 08 7 8 1 0 7 3N H3+→N ON2+ H2O+ H4 2 4 1 0 5 5 3 03 5 0 1 0 8 4N H3+ O→ 2 N O+ H2O+ H 3 5 0 1 0 2 7 6 55 2 4 1 0 8 5H N C O+→MH+ N C O+ M2 4 0 1 0 1 4 0 8 52 8 5 1 0 8 6N C O+→N ON2O+ C O1 0 0 1 0 1 0 0- 1 6 3 1 0 6 7N C O+→O HN O+ C O+ H1 0 0 1 0 1 3 00 8N2O+→O HN2+ O2+ H2 0 0 1 0 1 2 04 1 9 1 0 7 9N2O+→MN2+ O+ M1 6 0 1 0 1 4 02 1 6 1 0 8 注 反应 1和反应 2为尿素分解反应[ 7 ]; 反应 3~ 反应 9为简化还原反应[ 3 ]; M为反应的中间体. 1 2 边界条件及初始条件 计算原型为3 0 0MW 四角切圆煤粉炉, 炉膛宽 1 2 4 6 8m, 深 1 4 0 4 8m, 总高 5 2 1m, 具体参数如 图 1所示. 根据锅炉在 B MC R ( 最大连续蒸发量) 下的实际运行情况, 确定模型进口一次风、 二次风 速度分别为 2 5 , 4 4m/ s , 温度分别为 3 4 8 , 6 0 6K . 煤 粉的 速 度 及 温 度 与 一 次 风 相 同, 给 煤 总 量 为 3 6 6 7 k g / s , 煤质分析如表2 所示. 锅炉水平烟道中布置有大量的过热器、 再热器 管道, 若按照实际状况来模拟, 划分的网格数会较 多, 对计算机的要求较高, 且耗费大量时间. 所以本 文将过热器、 再热器管群简化为一定厚度的屏, 并 且适当减少屏的数量, 这样网格数会相应减少. 同 图 1 四角切圆煤粉炉各示意图( 单位 m) 296东南大学学报( 自然科学版) 第 4 2卷 h t t p / / j o u r n a l . s e u . e d u . c n 时考虑到计算机的硬件限制, 将燃烧及 S N C R部 分分开计算, 并将炉膛出口的速度、 温度、 浓度场导 入 S N C R计算域的进口. 简化后燃烧部分的网格 数为 1 4 3 1 0 6, S N C R部分的网格数为 7 8 1 05. 尿素溶液喷射方式采用 s o l i d c o n e 类型, 锥形 角为1 4 , 粒径分布符合 r o s i n r a m m l e r , 最小粒径为 0 0 3m m, 最大粒径为 0 8m m, 稀释后的尿素溶液 质量分数为 5 %. 在实际运行时, 炉膛的温度及烟 气流动基本稳定. 确定喷射位置后, S N C R过程的 脱硝效率与尿素溶液的喷射速度 u 、 颗粒粒径 d及 氨氮摩尔比 r 有关. 计算时, 改变其中一个变量, 其 他的按照基准参数设定( 见表 3 ) . 表 2 燃料的元素分析及工业分析 元素分析/ %工业分析/ % w ( Ca r)w ( Ha r)w ( Oa r)w ( Na r)w ( S a r) A a r V a r Ma r 低位发热量/ ( k J k g - 1) 5 6 0 63 38 0 10 6 10 5 62 1 8 62 3 4 09 6 02 1 2 0 0 表 3 计算参数 参数u / ( ms - 1) d / m mr 基准参数3 00 31 2 变化范围1 5~ 4 00 1~ 0 70 6~ 2 . 0 2 结果与讨论 燃烧的模拟首先进行连续相的流动计算, 然后 与离散相进行耦合, 最后加入传热及燃烧反应. 图 2为炉内的温度场纵向剖面. 可以看出, S N C R合 适的温度窗口大约在标高 4 0~ 4 5m之间. 图 2 锅炉纵向剖面温度分布( 单位 K) 2 1 各因素对 S N C R过程的影响 2 1 1 还原剂喷入点位置的影响 还原剂的喷入位置很大程度上影响到还原效 果, 温度过高, 还原剂易被氧化; 温度过低, 造成还 原剂的逃逸, 对下游设备造成腐蚀. 吕洪坤等[ 1 0 ]对 电站锅炉 S N C R过程进行了试验研究, 得出最佳 脱硝炉膛温度约为 12 2 3K , 选择在高温侧喷射能 够减少 N H 3的逃逸, 但同时也会有较多 N2O生 成. 本文选择在平均温度 1 4 4 2 , 1 2 0 0 , 1 1 0 0K 的平 面上布置喷嘴, 标高分别为 3 5 , 4 0和 4 2m, 由下往 上称为第 1层, 第 2层和第 3层. 其中第 1层在前 墙、 左墙和右墙上共布置 2 1个喷嘴, 第 2 , 3层只在 前墙布置 9个喷嘴. 为了研究喷入点位置对脱硝的 影响, 本文分别对 4种工况进行了研究比较, 分别 为第 1层、 第 2层、 1 , 2层组合及 1 , 2 , 3层组合, 设 定每个喷嘴喷射量相同, 模拟结果如图 3 ( a ) 所示. 图 3 各因素对脱硝效率及漏氨的影响 396第 4期姜敏, 等 3 0 0MW 燃煤锅炉 S N C R过程的数值模拟 h t t p / / j o u r n a l . s e u . e d u . c n 由图 3 ( a ) 可得, 单独投运第 1层获得了最大 的脱硝效率及最小的漏氨量. 炉膛三面布置喷枪的 方式使得溶液更易到达炉膛中心, 有利于烟气与尿 素溶液的均匀混合. 尿素溶液进入炉膛经水分蒸发 分解后再与 N O x反应, 粒径0 3m m的液滴在炉膛 中的停留时间约为 0 5~ 0 6s , 而烟气在折焰角处 的烟气速度为 7~ 1 0m/ s , 假设液滴的速度与烟气 流速相同, 则 0 3m m的尿素溶液液滴在炉膛中运 动 3 5~ 6m后开始分解成氨气和 H N C O , 后者再 与 N O x混合发生还原反应. 所以第 1层虽然喷射 点温度较高, 但大部分液滴在反应区间内蒸发完并 与 N O x反应, 获得了较好的脱硝效果. 第 2层与第 3层喷射的液滴大部分在温度区间外完成蒸发过 程, 因温度过低反应不完全而发生氨逃逸. 单独投运第 1层, 虽能够获得较好的脱硝效 果, 但由于温度高, 会生成大量的 N 2O ; 单独投运 第 2层, 效果较差; 1 , 2层组合以及 1 , 2 , 3层组合 脱硝效果较好, 但是 3层全部投运引起氨的逃逸较 多, 对下游设备有影响. 综上所述, 1 , 2层组合效果 最佳. 2 1 2 液滴喷射速度的影响 液滴喷射速度从 2 0m/ s 增加到 4 0m/ s , 对还 原率及漏氨的影响如图 3 ( b ) 所示. 由图可见, 随着 速度的增加, 还原率有所增加, 漏氨量减少. 喷射速 度对脱硝效果的影响主要表现在对混合和蒸发 2 个过程. 一方面增大液滴喷射初速度有利于增加液 滴的轨迹长度, 使得液滴更易穿透到主流区深处, 使得烟气与还原剂的混合更均匀; 另一方面, 速度 的增加有利于气流与液滴之间的传热传质, 蒸发速 度加快, 蒸发时间变短, 颗粒在反应区间外完成蒸 发的概率减少[ 1 1 ]. 当液滴喷射速度为 2 0m/ s 时, 液滴的平均停留时间为 0 9 2s ; 而速度增加到 4 0 m/ s 时, 缩短到 0 8 4s , 停留时间大大减少. 以上 2 个因素导致还原率随液滴喷射速度的增大而增大. 但是液滴喷射速度大于 3 0m/ s 后的还原率变化并 不明显, 原因可能是喷射速度对蒸发的影响变小. 而喷射速度越大, 耗能则越高. 所以液滴喷射初速 度选择 3 0m/ s 较合适. 2 1 3 液滴雾化粒径的影响 粒径小的液滴, 由于刚性差, 穿透能力小, 不能 到达炉膛中心, 只能还原近壁面部分的氮氧化物, 导致还原率较小; 相反粒径大的液滴, 刚性强, 能够 到达炉膛中心, 但是蒸发过程需要的时间较长, 可 能在反应区间内未蒸发完, 导致还原剂在合适的温 度区间内停留时间过短, 反应不完全, 出口的氨含 量增加, 还原率降低. 图 3 ( c ) 的模拟结果也证明了 这一点. 随着液滴粒径的增大, 漏氨量呈现上升的 趋势, 而平均粒径为 0 3m m的雾化条件, 穿透能 力好且在反应窗口的停留时间长, 获得了最佳的还 原率和较低的漏氨量. 2 2 模拟结果比较 在不同氨氮摩尔比下模拟脱硝效率, 由图 4可 见, 摩尔比越大, 喷入的还原剂越多, 还原率增大, 在氨氮摩尔比大于 1 8后还原率趋于稳定. 由模拟 值与电厂实测值拟合出的二次曲线可见, 当氨氮摩 尔比为 0 6时, 还原率相差最大, 相对误差小于 1 0 %. 随着氨氮摩尔比的增加, 2条曲线逐渐接近, 模拟结果与实测结果较为吻合[ 1 2 ]. 氨氮摩尔比对 还原率的影响, 一方面是由还原剂量的变化引起 的; 另一方面, 液流量的变化也有一定作用. 氨氮摩 尔比越大, 尿素溶液流量也越大. 增加液流量, 延长 了蒸发时间. 在速度不变的情况下, 液滴更容易进 入炉膛中心, 增加混合, 有利于还原率的提高. 但 是, 液流量的增加会使温度降低, 影响了再热器、 过 热器的传热效率. 由此可见, 利用数值计算的方法 模拟 S N C R过程, 其结果对实际运行一定有理论 指导作用. 图 4 模拟值与实测值的比较 3 结论 1 )通过对 4种组合方式脱硝效果对比, 研究 了尿素溶液的喷射位置对还原效果的影响. 尿素溶 液的整个蒸发过程必须在温度窗口内完成, 才能获 得更好的还原效果. 2 )尿素溶液的喷射速度对颗粒的蒸发及混合 过程有较大影响, 速度越大, 蒸发所需时间越短, 还 原剂与 N O x的混合越充分, 脱硝效果更好. 但是当 速度大于 3 0m/ s 后, 效果变得不明显, 实际运行过 程中可以选用 3 0m/ s 的喷射速度. 3 )雾化粒径越大, 颗粒刚性越强, 还原剂能到 496东南大学学报( 自然科学版) 第 4 2卷 h t t p / / j o u r n a l . s e u . e d u . c n 达烟气主流区, 增强混合过程, 但是蒸发时间也越 长. 颗粒可能在反应区间外蒸发完, 导致氨和 H N C O在温度窗口的停留时间短. 颗粒粒径平均为 0 3m m, 刚性强, 在温度窗口有足够的停留时间, 能够获得良好的脱硝效果. 4 )通过数值模拟的方法来研究 S N C R过程的 影响因素, 确定最佳运行条件, 对实际运行过程有 一定指导作用. 参考文献 ( R e f e r e n c e s ) [ 1 ]J a v e dM T , I r f a nN , G i b b s BM.C o n t r o l o f c o m b u s t i o n g e n e r a t e dn i t r o g e no x i d e sb ys e l e c t i v en o n c a t a l y t i cr e d u c t i o n[ J ] .J o u r n a lo fE n v i r o n m e n t a lMa n a g e m e n t , 2 0 0 7 , 8 3 ( 3 ) 2 5 1 2 8 9 . [ 2 ]N g u y e nTDB , K a n gT , L i m Y , e t a l .A p p l i c a t i o no f u r e a b a s e dS N C Rt oa m u n i c i p a l i n c i n e r a t o r o n s i t e t e s t a n dC F Ds i m u l a t i o n[ J ] .C h e m i c a l E n g i n e e r i n gJ o u r n a l , 2 0 0 9 , 1 5 2 ( 1 ) 3 6 4 3 . [ 3 ]N g u y e nTDB , L i m Y , K i m S , e t a l .E x p e r i m e n t a n d c o m p u t a t i o n a lf l u i d d y n a m i c s( C F D) s i m u l a t i o n o f u r e a b a s e ds e l e c t i v e n o n c a t a l y t i c r e d u c t i o n( S N C R )i na p i l o t s c a l ef l o w r e a c t o r[ J ] .E n e r y&F u e l , 2 0 0 8 , 2 2 ( 6 ) 3 8 6 4 3 8 7 6 . [ 4 ] C h o i CR , K i m CN .N u m e r i c a l i n v e s t i g a t i o no nt h e f l o w , c o m b u s t i o na n dN Oxe m i s s i o nc h a r a c t e r i s t i c si na 5 0 0MWet a n g e n t i a l l yf i r e dp u l v e r i z e d c o a l b o i l e r [ J ] . F u e l , 2 0 0 9 , 8 8 ( 9 ) 1 7 2 0 1 7 3 1 . [ 5 ]B i r k h o l dF , Me i n g a s t U , Wa s s e r m a n nP , e t a l .Mo d e l i n g a n ds i m u l a t i o no f t h e i n j e c t i o no f u r e a w a t e r s o l u t i o nf o r a u t o m o t i v e S C RD e N O x s y s t e m s [ J ] . A p p l i e dC a t a l y s i s B E n v i r o n m e n t a l , 2 0 0 7 , 7 0 ( 1 / 2 / 3 / 4 ) 1 1 9 1 2 7 . [ 6 ]吕洪坤, 杨卫娟, 周俊虎, 等.尿素溶液高温热分解特 性的实验研究 [ J ] .中国电机工程学报, 2 0 1 0 , 3 0 ( 1 7 ) 3 5 4 0 . L H o n g k u n , Y a n gWe i j u a n , Z h o uJ u n h u , e t a l .I n v e s t i g a t i o no nt h e r m a ld e c o m p o s i t i o no nc h a r a c t e r i s t i c so f u r e as o l u t i o nu n d e r h i g ht e m p e r a t u r e[ J ] .P r o c e e d i n g s o f t h eC S E E , 2 0 1 0 , 3 0 ( 1 7 ) 3 5 4 0 .( i nC h i n e s e ) [ 7 ]A o k i H ,F u j i w a r a T ,Mo r o z u m i Y ,e t a l .Me a s u r e m e n t o f u r e at h e r m a l d e c o m p o s i t i o nr e a c t i o nr a t ef o r N Os e l e c t i v en o n c a t a l y t i cr e d u c t i o n[ C ] / / F i f t hI n t e r n a t i o n a l C o n f e r e n c e o nT e c h n o l o g i e s a n dC o m b u s t i o nf o r aC l e a n E n v i r o n m e n t .L i s b o n ,P o r t u g a l , 1 9 9 9 1 1 5 1 1 8 . [ 8 ]R o t aR , A n t o s D , Z a n o e l o , e t a l .E x p e r i m e n t a l a n d m o d e l i n ga n a l y s i s o f t h e N OxO U Tp r o c e s s [ J ] .C h e m i c a l E n g i n e e r i n gS c i e n c e , 2 0 0 2 , 5 7 ( 1 ) 2 7 3 8 . [ 9 ]B r o u w e r J , H e a pM P ,P e r s h i n gDW, e t a l .Am o d e l f o r p r e d i c t i o no f s e l e c t i v en o n c a t a l y t i cr e d u c t i o no f n i t r o g e no x i d e s b ya m m o n i a , u r e a , a n dc y a n u r i ca c i dw i t h m i x i n gl i m i t a t i o n s i nt h ep r e s e n s eo f C O[ C ] / / T w e n t y s i x t hS y m p o s i u m ( I n t e r n a t i o n a l ) o nC o m b u s t i o n / T h e C o m b u s t i o nI n s t i t u t e .N a p l e s , I t a l y , 1 9 9 6 2 1 1 7 2 1 2 4 . [ 1 0 ]吕洪坤, 杨卫娟, 周俊虎, 等.电站锅炉选择性非催化 还原脱硝实验研究 还原区温度、尿素溶液喷射 体积流量的影响 [ J ] .浙江大学学报 工学版, 2 0 0 9 , 4 3 ( 9 ) 1 6 5 5 1 6 6 0 . L H o n g k u n , Y a n gWe i j u a n , Z h o uJ u n h u , e t a l . E x p e r i m e n t a l r e s e a r c ho nS N C Ra p p l i e di np o w e r p l a n t b o i l e r i m p a c t o ft e m p e r a t u r ea n du r e as o l u t i o ni n j e c t i o n v o l u m ef l o wr a t e [ J ] .J o u r n a l o f Z h e j i a n gU n i v e r s i t y E n g i n e e r i n gS c i e n c e , 2 0 0 9 , 4 3 ( 9 ) 1 6 5 5 1 6 6 0 .( i n C h i n e s e ) [ 1 1 ]王海涛, 杨卫娟, 周俊虎, 等.液滴在高温气流中蒸发 混合特性计算分析 [ J ] .浙江大学学报 工学版, 2 0 1 1 , 4 5 ( 5 ) 8 7 8 8 8 4 . Wa n gH a i t a o , Y a n gWe i j u a n , Z h o uJ u n h u , e t a l .C a l c u l a t i o na n da n a l y s i so ne v a p o r a t i o na n dm i x i n gc h a r a c t e r i s t i c s o f d r o p l e t s i nh i g ht e m p e r a t u r e f l u e [ J ] .J o u r n a l o f Z h e j i a n gU n i v e r s i t y E n g i n e e r i n gS c i e n c e , 2 0 1 1 , 4 5 ( 5 ) 8 7 8 8 8 4 .( i nC h i n e s e ) [ 1 2 ]李芳芳, 吴少华, 曹庆喜, 等.湍流混合限制下选择性 非催化还原过程的数值模拟 [ J ] .中国电机工程学 报, 2 0 0 8 , 2 8 ( 2 3 ) 4 6 5 2 . L i F a n g f a n g ,WuS h a o h u a , C a oQ i n g x i , e t a l .N u m e r i c a ls i m u l a t i o n o fs e l e c t i v e n o n c a t a l y t i c r e d u c t i o n p r o c e s s w i t ht u r b u l e n c e m i x i n gi n f l u e n c e [ J ] . P r o c e e d i n g s o f t h e C S E E , 2 0 0 8 , 2 8 ( 2 3 ) 4 6 5 2 . ( i nC h i n e s e ) 596第 4期姜敏, 等 3 0 0MW 燃煤锅炉 S N C R过程的数值模拟