利用燃煤电厂尾部烟道的烟气余热来实现脱硫废水的喷雾蒸发是实现其零排放的有效途径，以国内某燃煤电厂330MW火力机组的烟道为研究对象，利用DPM模型对雾化液滴群在高温烟道内的蒸发及流动特性进行了研究，考察了不同雾化嘴角情况下液滴碰壁情况、不同负荷下液滴的蒸发情况，研究结果表明：在50%、75%、100%烟气负荷工况下，烟气温度越高、烟气速度越快，雾化液滴群完全蒸发所需时间越少，液滴最大蒸发时间在2.85～3.36s之间。在单烟道结构的最佳喷嘴雾化锥角为65°情况下，越靠近烟道内侧，涡的尺寸越大，越有利于促进喷嘴区的局部液滴群不断向其他区域扩散。
中国是以煤炭为主要能源的国家，2017年燃煤火力发电量占全年总发电量的67%。发电过程中煤炭燃烧产生的二氧化硫排放问题尤为引人关注，在一定的气象条件下产生复杂的化学反应，是形成雾霾和酸雨的重要前驱体。石灰石-石膏湿法烟气脱硫工艺应用最广，然而，循环浆液将持续富集来自烟气及脱硫剂中的重金属元素和氯离子，从而产生高浓度的脱硫废水，废水直接排放对环境产生负面影响。
如采用常规工艺进行废水零排放处理，则高浓度氯离子的腐蚀性对设备材质要求很高，造价昂贵。使用喷嘴将脱硫废水雾化为液滴群并喷入空气预热器至电除尘器间的烟道内，利用高温烟气与常温废水的传热作用实现脱硫废水的零排放，有投资少、工艺流程短、去除重金属离子、建设工期短、维护成本低等特点，被推荐为实现脱硫废水零排放的可行性技术。针对脱硫废水液滴群在烟气中蒸发与流动特性的优化是实现脱硫废水烟道蒸发零排放的关键。
目前，国内外对于脱硫废水烟道蒸发工艺的研究主要集中在以数值模拟的方式研究脱硫废水蒸发特性、流动特性两方面，同时，伴以一定的工程或实验数据作为参照。张子敬等研究认为喷雾液滴群蒸发特性受到液滴加热升温(传热过程)和喷雾液滴群在烟气中的扩散(传质过程)两方面的共同作用。Strotos G等建立了单个液滴在高温燃气中蒸发、运动过程的数学模型，获得了不同燃气温度和速度下液滴的蒸发规律。
冉景煜等对不同物性液滴在低温烟气环境中的运动，以及受热和蒸发过程中的传热传质特性进行了理论分析。李明波等通过计算流体动力学软件Fluent，对空气预热器出口至电除尘器入口段烟道内的烟气流动情况进行了模拟。
Laín等的以拉格朗日湍流颗粒分散体模型的建立为基础，提出携带稀薄粒子的气流在一定条件下，假设粒子为球体，只考虑曳力和重力作用。Young等应用离散多组分(DMC)燃料液滴模型对多组分燃料喷雾的蒸发进行了数值模拟。Pinto等研究了双流体喷嘴的喷雾干燥，成功地预测了干燥时间和最终含水量随着初始液滴直径变化的趋势。
晋银佳等等提出深度过滤脱硫废水预处理工艺，将脱硫废水在雾化蒸发前进行固液深度分离预处理以解决硫废水中悬浮颗粒物堵塞问题。
国内外学者已对液滴蒸发的机理进行了深入研究，重点考察烟气温度、速度、液滴直径、液滴速度对蒸发的影响，但是，不同雾化嘴角对脱硫废水蒸发的影响尚未有明确的解释，文中结合国内某燃煤电厂330MW机组空气预热器至电除尘器间的烟道中喷雾蒸发实现脱硫废水零排放工程实践，数值模拟不同烟气负荷和不同喷雾锥角对脱硫废水喷雾蒸发流动特性的影响。
1 方法与模型
脱硫废水在烟道中喷雾蒸发属于典型的气液两相流流动，在数值模拟中以空气为连续相，以喷雾液滴为离散相，主要考虑连续相和离散相之间的相间运动和相互作用。首先，建立烟道的物理模型，根据连续相和离散相方程，以确定的边界条件进行相应数值模拟计算。
1. 1 物理模型
图1所示为空气预热器与电除尘器之间烟道和尺寸的物理模型。烟道分为入口段、下弯头、竖直段烟道、上弯头、异型弯头和水平烟道6个部分。采用ANSA软件对烟道几何模型进行网格划分，该烟道模型结构简单，流场结构均匀，在计算速度上采用有明显优势的全六面体网格，生成总网格数为200万。
经检验，该模型网格EquiSize Skew值在0～0.4之间的网格数占98.09%，网格划分质量较高。采用网格数分别为200万和300万和400万的网格进行无关化验证，对竖直段烟道内的6个点进行速度监测，3种网格计算结果相差不大，为了节省计算资源，选择网格数量约为200万的网格进行模拟，如图2所示。
1. 2 数学模型
1. 2. 1 连续相方程
在气液两相流动中，尽管控制方程独立，两相却是相互耦合的。液滴作为质量源、动量源和能量源被引入到气相方程中，并通过这些源项影响气相流场，气相流场又反过来通过其速度场、温度场、压力场等来影响液滴的本身状态。下列方程为气相控制方程，其表达式分别如下。
连续性方程：
2 结果与分析
2. 1 烟气负荷对液滴群蒸发及运动过程的影响
脱硫废水在锅炉尾部烟道中的雾化及流动、蒸发过程可分为初始阶段和稳态阶段。初始阶段，常温液滴群作为吸热蒸发的分布热汇，充分吸收烟气流的余热，所吸收的热量大部分用于液滴群温度的升高，同时，在烟气速度的影响下，该阶段的液滴群速度不断增大;在很短时间内，雾化液滴群即达到稳态阶段，此时，液滴群被烟气加热到稳定值，吸收的所有热量都用于液滴群蒸发，液滴群速度与来流烟气速度一致。
液滴群的蒸发效果主要由以下参数共同决定：气相温度、传输特性、液相温度、运动速度以及气液两相的传热、传递效率，分别选取330MW机组50%、75%、100%烟气负荷工况下，3种不同烟温( 120.3、125.1、128.9℃)及烟速( 9.19、11.56、14.64m/s)的气相条件对脱硫废水蒸发及流动特性的影响作定量分析，并结合传热传质理论加以解释。
图3显示了50%、75%、100%3种不同烟气负荷下，以不同的雾化锥角进行喷雾，运动液滴最大蒸发时间T的模拟结果， T值随烟气负荷的增加呈现近乎相同的线性下降趋势。随着负荷的增加烟气量增加，烟道烟气温度降低减少，蒸发时间减少，其中，50%、75%及100%烟气负荷工况运动液滴最大蒸发时间T分别在3.07～3.36s、2.85～3.04s和2.57～2.80s范围内。
选取喷嘴雾化锥角65°配置下的各烟气负荷颗粒运动轨迹，如图4所示。
液滴群颗粒皆能蒸发完全，100%烟气负荷对应的最大蒸发时间最短，50% 烟气负荷对应的最大蒸发时间最长，由此可见，对于相同粒径的液滴，气体环境温度越高、烟气速度越快，液滴群的汽化速率越高、蒸发效果越好。
其中，由于100%负荷下烟气速度相较于75%和50% 负荷时更快，则脱硫废水颗粒衰减后的速度仍然较快，若烟道长度不足，仍有蒸发不完全的可能性，从图中可看出，烟气速度的变化对液滴最大完全蒸发时间的影响较小，故在单烟道结构中，烟气温度对蒸发效果起主导作用。
若烟气温度升高，则气液两相的温差增大，气体环境向液滴群的传热增强，从而使液滴表面蒸发及传质扩散速率不断增大，因此，液滴温度持续升高，其到达临界蒸发温度所需时间变短，液滴自喷入烟道至完全蒸发的停留时间随烟气温度升高而逐渐减少。
2. 2 雾化锥角对液滴群蒸发及运动过程的影响
为定量分析雾化锥角对雾化液滴群流动特性的影响，定义被烟道壁面捕捉的液滴数量占液滴颗粒总数比为A0。A0值可反映出脱硫废水喷雾蒸发结晶后，在烟道内壁积灰的可能性大小。
图5显示了在20°、35°、50°、65°、80°、95°6种不同雾化锥角下在单烟道壁面被捕捉的液滴数量分数的模拟结果，A0值随雾化锥角的变化呈现近乎相同的先平稳下降、后明显上升趋势。
图5表明：在雾化锥角由20°至50°增加的过程中，A0值变化相对平稳，由于雾化角过小，液滴蒸发速度较慢，易撞击顶部水平烟道;当雾化锥角增加至65°，烟道捕捉的液滴数达到最小值，说明65°雾化锥角在烟道内壁积灰可能性最小;雾化锥角由65°至95°继续增大的过程中，A0值呈明显增加趋势，此时，由于雾化角过大液滴易撞击竖直烟道，但雾化锥角大于90°后，增加速率有所放缓，且有下降趋势，随着雾化角的增大，液滴蒸发速度变快，液滴碰壁的可能性变小。
当喷嘴雾化锥角过小时，相同工况下液滴蒸发较慢。当液滴进入水平烟道时，由于液滴的直径相对较大，随流能力也就越弱，液滴越撞击水平烟道形成积灰。当喷嘴雾化锥角过大时，液滴容易直接撞击竖直烟道形成积灰。因此，存在1个最佳的雾化锥角使液滴的碰壁数量最小，经过验证当雾化锥角为65°时撞击烟道的液滴数量最小。
单烟道结构75%烟气负荷工况下，最佳雾化锥角65°时，对于脱硫废水蒸发及流动特性的定量及烟道截面速度矢量图，如图6所示。
由图6可知，在喷雾蒸发的初始阶段，传质扩散及蒸发速率较快，喷雾对烟气的剪切卷吸形成了一个较大的不规则的涡。
由于烟道内侧的烟气体积流量较大，喷嘴截面沿烟气流动方向1m处烟气以较快的速度冲入对墙，造成其上部有较大压强差而形成回流，故越靠近烟道内侧，涡的形态越大，有利于促进喷嘴区的局部液滴群不断向其他区域扩散。随着蒸发及传质扩散的进一步均匀化，喷雾蒸发进入稳态阶段，烟道通流截面涡增大，截面涡的形态逐渐规则化，速度矢量场趋于稳定。
3 结 论
1 ) 50%、75%、100%3种烟气负荷工况下，在单烟道壁面被捕捉的液滴数量分数随雾化锥角的增加皆呈现先平稳下降、后明显上升趋势。
2)在20°、35°、50°、65°、80°、95°6种不同雾化锥角下运动液滴最大蒸发时间值随烟气负荷的增加呈现近乎相同的线性下降趋势。在最佳喷嘴雾化锥角65°配置下，对于相同粒径的液滴，气体环境温度越高、烟气速度越快，液滴群的汽化速率越高、蒸发效果越好。其中，烟气速度的变化对液滴最大完全蒸发时间的影响较小，烟气温度对蒸发效果起主导作用。脱硫废水喷雾后形成的液滴群可在烟道中完全蒸发。
3 )最佳雾化锥角配置下的速度矢量图显示，越靠近烟道内侧，涡的尺寸越大，有利于促进喷嘴区的局部液滴群不断向其他区域扩散;喷雾蒸发初始阶段的传质扩散及蒸发速率较快，速度矢量图呈现出一个较大的不规则的涡形态;喷雾蒸发稳态阶段烟道通流截面涡增大、形态逐渐规则化，速度矢量场趋于稳定。
