内部输水性壁面波形除雾器性能研究

彭玉凤, 陈 倬, 朱家骅

内部输水性壁面波形除雾器性能研究

彭玉凤, 陈 倬, 朱家骅

(四川大学 化学工程学院, 四川 成都 610065)

利用内部输水性材料作为波形除雾器壁面有利于液体传输，从而提升除雾性能。对该除雾器雾滴捕集率、液体传输速率及夹带抑制机理进行了研究，建立了预测雾滴捕集率的半经验模型和模拟计算方法，并得到实验验证，据此建立了内部输水性壁面降膜厚度与材料输水特性关联式，解释了材料内部输水性能对该除雾器的夹带抑制机理。模型计算表明，常规除雾器气速超过临界值4.82 m×s-1时即会导致液膜剪切产生夹带，而内部输水性壁面除雾器接近该气速时壁面尚处于无液膜状态，气速超过临界值19.52 m×s-1才开始产生液沫夹带，表明内部输水壁面能有效降低液膜厚度，且液滴在壁面完全铺展，抑制夹带发生，从而强化除雾器处理能力、提高分离效率。

输水；雾滴；除雾器；液膜

1 前言

波形除雾器能有效分离气流中的雾滴，被广泛应用于化工、能源、海水淡化等工业领域[1-2]。其优点之一是压降低，有利于采用高气速以提高除雾器的处理能力，但不利之处是波折板表面聚集的雾滴会被气流二次夹带使捕集率随之下降[3]，反而导致除雾器总体性能下降[4]。为了抑制二次夹带，一种方法是在波折板壁面安装排水钩[5]，但除雾器压降会急剧增加，且使设备结构复杂，使用效果不佳[6]。另一种方法是波折板表面复合多孔泡沫层[7]，减轻液体表面与气流的直接剪切作用，但液体在复杂多孔结构中传输阻力大，排液能力小，尤其在波折面水平放置、气液逆流工况下，雾滴二次夹带仍然不可避免。

由此可见，抑制二次夹带不仅要避免气液逆流，更重要的是及时移除捕集的雾沫、避免或降低其在壁面的聚集量。本研究采用内部输水性布料构建了波折面直立的气液错流波形除雾器(图1(a)、(b))，除雾器壁面捕集到的液体在毛细力和重力联合作用下通过布料纤维内部传输、从底部排出，力图使壁面保持无液膜。目前除雾器壁面输水过程研究[7-9]尚未见气液错流下壁面内部输水的报道，有必要对此类除雾器内部输水性壁面二次夹带抑制机理、输水能力和除雾器负荷的预测模型进行系统的研究。

图1 气液错流波形除雾器实验装置图

根据图1(c)所示实验装置获取的雾沫捕集数据分析了该过程机理，建立了雾滴捕集率预测模型；结合布料输水性能参数测试，推导了内部输水性壁面膜厚预测关联式，结合二次夹带临界气速判定模型，与常规壁面除雾器对比绘制了气速-膜厚曲线图，定量表达了内部输水性壁面夹带抑制机理，并可为合理选择除雾器操作气速提供理论参考。

2 实验与模拟

2.1 雾滴捕集

稳定工况下除雾器波形折板壁面捕集的雾滴量即为壁面内部输水量及外部降膜输水量的总和。为获得雾滴捕集率与壁面降膜膜厚之间的定量关系，首先需要根据雾滴捕集率实验数据建立雾滴捕集率预测模型。

2.1.1 雾滴捕集实验研究

本研究设计的气液错流波形除雾器实验装置及其几何结构如图1所示。表1中列出其详细几何参数。实验过程中，鼓风机提供0.5～1.0 kPa(表压)的循环空气经计量后进入喷雾段，携带喷雾器产生的细雾滴通过混合段均化，进入波形除雾器之前在入口测量段用马尔文帕纳科Spraytec喷雾粒度仪对气流中的雾滴进行表征，其粒径分布如图2所示，入口雾滴的粒径p在1～30 μm，索特平均直径3,2为7.04 μm。在除雾器迎流面用等动力学采样管[10]将含雾气流引入装有200 mL纯水的洗气瓶中，雾滴(浓度为0.01 mol×L-1的KCl稀溶液)带入的电解质增量Δe与洗气瓶内液体电导率变化量Δ满足标定的线性关系。实验在20 ℃下进行，采用同批配制的KCl稀溶液喷雾，并使用等动力学采样管、等容量采样方法。该除雾器共有120级波折板，每间隔10级波折板设置采样点，通过测量洗气瓶中电导率变化量即可从标定曲线得到该次采样气流中的雾滴质量分数。根据所测进口和出口处雾滴质量分数in和out，通过下式计算除雾器雾滴平均质量捕集率：

表1 波形除雾器结构尺寸

图2 除雾器入口雾滴的粒径分布

2.1.2 雾滴捕集模拟研究

电导率测量法获得的是除雾器总的质量捕集率，作为一种验证手段，借助图2所示粒径分布数据通过加权平均，可以检验数值模拟方法获得的不同粒径雾滴的粒级捕集率，这对除雾器性能表征具有重要意义。模拟基于以下合理化假设[11-13]：

(1) 除雾器通道由多个等间距的单通道并行排列而得，通道高度远大于通道间距与波折长度，因此将三维多通道多相流简化成二维单通道；

(2) 实验中气流Mach数远小于0.1，视为不可压缩流；

(3) 雾滴被视为球形液滴，只受流体曳力作用，且忽略液滴之间碰撞；

(4) 液滴接触壁面即认为捕集，不考虑液滴破碎及溅射造成的夹带影响。

使用ICEM CFD创建二维几何模型，对10级波折板(=10)的除雾通道进行结构化网格划分，并对近壁面附近区域利用边界层网格进行加密处理。网格划分完毕后导入Fluent 2019 R2进行数值计算。在模拟中，气相流场基于欧拉法求解，由于其流动为湍流，采用了Realizable-湍流模型[10,19]进行求解，用SIMPLE算法求解控制方程，并采用二阶迎风法将气体流动方程离散化。各变量的收敛标准设置为10-5。液滴轨迹则基于拉格朗日法进行计算，采用DPM(discrete phase model)模型模拟液滴运动。

根据被捕集的单一粒径雾滴质量流量q,d,t与其进口质量流量q,d,i之比，即可确定该粒径雾滴模拟捕集率：

2.2 内部输水性能实验

内部输水性能直接影响被捕集雾滴在布料壁面上的传输和成膜，由此影响捕集壁面的降膜及雾沫夹带。为研究不同布料的内部输水性能，设计了图3所示的布料内部输水速率测量实验，采用长宽为30 cm×10 cm的3种布料(苎麻布料、亚麻布料、棉布料)在室温下进行实验。

图3 布料内部输水实验装置图

实验中，水槽水位由外置恒液位槽控制，布料一端浸入水槽中，液体在毛细力与重力联合作用下于布料内部传输，在布料另一端汇聚成液滴滴落到锥形瓶中并由电子天平称量。在实验过程中(0 ～ 3.5 h)，布料表面未观察到液膜存在，表明液体转移全部通过布料内部传输完成。根据所测宽度为的布料在不同时间段Δ内传输液滴质量Δ，可计算出单位宽度布料内部的输水速率Δt：

3 实验结果与讨论

3.1 雾滴捕集实验

图4为波折板级数=10时粒级捕集率随气速的变化。可以看出，随气速增大，雾滴惯性越大，捕集率增大；粒径增大也使捕集率上升，当粒径达到20 μm时，雾滴接近全捕集，这是因为粒径越大的雾滴惯性越大，越不容易随气流方向改变运动方向，越容易被捕集。

图4 N =10时捕集率随气速的变化

波形除雾器为多级波折板的串联结构，Jackson等[14-16]认为，每级波折板对同一粒径的雾滴捕集率d,b均相等。因此除雾器单一粒径捕集率d,t随级板数变化的计算式为

模拟结果为单一粒径雾滴的捕集率，试验中利用马尔文激光粒度测量仪可得到某一粒径雾滴的体积分数(见图2)。由于雾滴密度为常数，其质量分数d等于其体积分数，采用权重加和方式，计算已知雾滴粒径分布的总质量捕集率：

该模拟计算加权统计值与实验结果具有同一基准下对比意义。

图5为不同气速下除雾器捕集率的实验值和模拟值的对比图，结果表明：在不同工况下，模拟值和实验值符合良好，相对误差均小于13.7%，说明本研究采用的计算模型可以较好地反映雾滴捕集情况。由图可见，增加气流速度和波折板级数均可提高除雾器捕集率。

图5 实验数据和模拟数据对比图

3.2 内部输水性能

如前所述以单位时间、单位布料宽度的输水速率Δt表示布料内部输水性能，3种布料在0 ～ 3.5 h内单位宽度输水质量累计值的实验检测结果如图6所示。由图6可见，干布料接触纯水初期，布料处于液体润湿扩散过程，前15 min输水量为0；随着传输时间延长，布料内部输水通道越来越通畅，输水量持续增加，输水量曲线斜率在45 min后基本保持恒定，输水质量与传输时间呈线性关系，说明形成内部输水通路后，输水速率Δt即保持稳定。由此可见，对于内部输水性布料，在除雾器长期稳态工况下，在不含堵塞性杂质的雾滴捕集过程中，布料内部液体传输速率Δt为恒定值，因此本研究设定传输稳定后的输水速率即为该布料的内部输水速率D。

研究选择了3种实验材料中内部输水速率最大的苎麻布料为除雾器装置内部的雾滴捕集壁面材料，其D=3.84×10-5kg×(m×s)-1。

图6 布料的内部输水质量随时间的变化

4 输水过程及夹带抑制机理分析

4.1 捕集模型的建立

除雾器壁面输水过程将影响雾沫夹带，而其壁面内部及外部传输的液体来自雾滴的捕集，二者互为因果，从机理分析角度更需要建立显式的雾滴捕集率计算模型。经典的Burkholz[17]雾滴捕集理论将粒径p与气流速度这2个关键参数归纳为液滴斯托克斯数，单级波折板雾滴捕集率d,b是和波折板偏转角的乘积：

式中：d为雾滴密度，kg×m-3；g为气体黏度，Pa×s。

但该模型未考虑湍流下液滴颗粒的混合，其对大的斯托克斯数范围预测偏差较大，本研究借鉴Wikinson模型[18]，引入湍流影响修正因子c，并利用前述模拟结果获得c的拟合表达式c=2.718×(4.44612+1)-0.6，由此得到新的单级波折雾滴捕集率预测模型：

该模型可计算不同气速下不同粒径的雾滴捕集率，为方便后续代入输水模型且不失代表性，借鉴文献[8]的做法，以雾滴粒径范围对应的索特平均粒径3,2代入式(8)计算得到的捕集率代表按粒径分布加权平均值，与本研究实验结果对比，其误差在±20% 以内，这说明在本实验条件下，采用该简化模型进行壁面输水过程机理分析是合理的。

4.2 输水及夹带模型

对于常规除雾器，随着雾滴不断被捕集，液滴或汇聚形成局部液膜沿壁面下降；或汇聚形成较大液滴，悬挂在壁面上，直到超过临界厚度后开始滑落[19]，故输水过程中壁面液膜厚度不均匀，如图7(a)所示。而内部输水性壁面由于其材料可浸润，被捕集的雾滴可立即被壁面吸收并铺展，若形成液膜也会是均匀的，如图7(b)所示。

图7 雾滴在不同壁面的输水过程示意图

基于雾滴捕集率计算模型式(8)，被捕集雾滴汇聚的体积流量(m3×s-1)：

式中：为板间距，m；为板高，m；in为入口气流中雾滴质量浓度，kg×m-3；雾滴捕集率为粒径3,2的代表值。

稳定工况下除雾器波形折板壁面捕集的雾滴量,t为壁面内部输水量,d,i(m3×s-1)及外部降膜输水量,d,s(m3×s-1)的总和：

基于内部输水性能实验测得的布料内部输水速率D，水平长度的壁面内部输水体积流量q,d,i为

根据Nusselt液膜理论，壁面液膜厚度与壁面降膜输水量q,d,s的关系为：

式中：为重力加速度，m×s-2。联立式(9)～(12)，得到内部输水性壁面波形除雾器液膜厚度关联式：

根据Wang等[20]夹带临界判定模型得到夹带临界气速gc与液膜厚度的关系：

式中：为波折板弯折处的曲率半径，m。

结合式(13)、(14)得到，壁面材料内部输水能力越强，壁面液膜厚度越薄，故发生夹带时的临界气速越高。

对于常规除雾器壁面，悬挂液滴处的局部膜厚大于降膜区域平均膜厚，更容易被气流剪切。液体与壁面之间的液固接触角越大，越易形成液滴悬挂，悬挂的液滴汇聚增大至脱落直径后沿壁面滑落。闵敬春等[21]研究了竖直平壁上液滴沿壁面下滑时的液滴临界厚度c(m)，可由式(15)得到

式中：l为液滴表面张力，N×m-1；A和R分别为前进和后退接触角，rad。

根据式(13)～(15)，绘制出不同壁面的液膜厚度与气速的关系图，如图8所示。图中a线为内部输水性壁面上的液膜厚度，其液膜为完全铺展状态；b线为假设液膜在常规非浸润壁面上为完全铺展状态时的膜厚。c线为悬挂液滴脱落临界厚度，d线为夹带临界线，e线为壁面悬挂液滴被夹带临界线。

在实际除雾过程中，雾滴在常规除雾器壁面上为液滴、液膜共存状态，即液膜厚度处于b、c线之间的A、B区域。在区域A内，其壁面上液膜和悬挂液滴的厚度均小于临界值，因此在区域A内不会发生夹带；当气速超过夹带临界气速(图中为4.82 m×s-1)时，气流的剪切力与壁面上悬挂液滴的黏性力无法平衡，此时在B区域内发生夹带，图中表现为d线(夹带临界线)穿越了B区域。

图8 除雾器膜厚与气速关系图

对于内部输水性壁面除雾器，当气速不超过降膜临界气速(图中为4.35 m×s-1)时，雾滴捕集量不会大于其壁面内部输水量，其壁面液膜厚度为0，完全不会发生夹带。气速超过降膜临界气速后，壁面开始出现降膜，随着气速增加，雾滴捕集量增加，壁面液膜厚度随之增加，气速超过夹带临界气速(理论计算值为19.52 m×s-1)后，才会产生夹带。相比于常规除雾器，苎麻布料的内部输水性能使液膜在捕集壁面上均匀铺展，故其壁面液膜厚度更薄，非夹带区域更大。

5 结论

(1) 雾滴捕集实验结果表明，采用提高气速、增加级数的措施可有效提高雾滴捕集率；基于惯性分离机理，建立了用于雾滴捕集率预测的半经验模型，模型计算值与实验值相比，误差在±20% 之内，能较好地预测不同工况下的雾滴捕集率。

(2) 通过对内部输水性波形除雾器壁面的输水过程进行分析，建立了膜厚预测关联式：不同内部输水性材料在稳定输水过程中具有不同的内部输水速率，壁面材料内部输水能力越强，壁面降膜厚度越薄。

(3) 内部输水壁面能有效降低液膜厚度，且液滴能在壁面完全铺展，从而抑制夹带发生。模型计算显示，常规除雾器当气速超过夹带临界气速4.82 m×s-1时即可能产生液沫夹带，而内部输水性壁面除雾器在气速4.35 m×s-1时尚处于无液膜状态，其临界夹带气速提高到19.52 m×s-1，可显著强化除雾器处理能力。

[1] Mao F, Tian R, Chen Y,Re-entrainment in and optimization of a vane mist eliminator [J]. Annals of Nuclear Energy, 2018, 120: 656-665.

[2] Venkatesan G, Kulasekharan N, Muthukumar V,Regression analysis of a curved vane demister with Taguchi based optimization [J]. Desalination, 2015, 370: 33-43.

[3] Song J, Hu X. A mathematical model to calculate the separation efficiency of streamlined plate gas-liquid separator [J]. Separation and Purification Technology, 2017, 178: 242-252.

[4] Azzopardi B J, Sanaullah K S. Re-entrainment in wave-plate mist eliminators [J]. Chemical Engineering Science, 2002, 57(17): 3557-3563.

[5] Galletti C, Brunazzi E, Tognotti L. A numerical model for gas flow and droplet motion in wave-plate mist eliminators with drainage channels [J]. Chemical Engineering Science, 2008, 63(23): 5639-5652.

[6] Kavousi F, Behjat Y, Shahhosseini S. Optimal design of drainage channel geometry parameters in vane demister liquid–gas separators [J]. Chemical Engineering Research and Design, 2013, 91(7): 1212-1222.

[7] Xu Y, Yang Z, Zhang J. Study on performance of wave-plate mist eliminator with porous foam layer as enhanced structure. Part II: Experiments [J]. Chemical Engineering Science, 2017, 171: 662-671.

[8] Xu J, Hrnjak P. Impinging oil separator for compressors [J]. International Journal of Refrigeration, 2020, 119: 110-118.

[9] Ruiz J, Cutillas C G, Kaiser A S,. Experimental study on pressure loss and collection efficiency of drift eliminators [J]. Applied Thermal Engineering, 2018, 149: 94-104.

[10] 张伟, 吴鑫宇, 李睿. 废气低浓度颗粒物采样常见问题的梳理 [J]. 环境与发展, 2020, 32(5): 126-127.

Zhang W, Wu X Y, Li R. Combating common problems in sampling low concentration particulates of exhaust gas [J]. Environment and Development, 2020, 32(5): 126-127.

[11] 郝雅洁, 刘嘉宇, 袁竹林, 等. 除雾器内雾滴运动特性与除雾效率 [J]. 化工学报, 2014, 65(12): 4669-4677.

Hao Y J, Liu J Y, Yuan Z L,Movement characteristics of droplets and demisting efficiency of mist eliminator [J]. CIESC Journal, 2014, 65(12): 4669-4677.

[12] 林其聪, 刘欣, 周翔, 等. 折线型与流线型除雾器性能的数值模拟与分析 [J]. 中氮肥, 2013(1): 34-35.

LIN Q C, LIU X, ZHOU X,Numerical simulation on performances of droplet eliminator [J]. M-Sized Nitrogenous Fertilizer Progress, 2013(1): 34-35.

[13] Liu Y, Qu Z. Numerical investigation of moisture separators with corrugated plates [J]. Energy Procedia, 2017, 105: 1501-1506.

[14] Jackson S, Calvert S. Entrained particle collection in packed beds [J]. AIChE Journal, 1966, 12(6): 1075-1078.

[15] Guan L, Yuan Z, Yang L,. Numerical study on the penetration of droplets in a zigzag demister [J]. Environmental Engineering Science, 2016, 33: 35-43.

[16] Wang Y I, James P W. The calculation of wave-plate demister efficiencies using numerical simulation of the flow field and droplet motion [J]. Chemical Engineering Research and Design, 1998, 76(8): 980-985.

[17] Bürkholz A. Droplet separation [M]. New York: VCH Publishers, 1989.

[18] Wilkinson D. Optimizing the design of waveplates for gas-liquid separation [J]. Proceedings of the Institution of Mechanical Engineers, Part E: Journal of Process Mechanical Engineering, 1999, 213(4): 265-274.

[19] Kim S W, Kim J, Park S S,Enhanced water collection of bio-inspired functional surfaces in high-speed flow for high performance demister [J]. Desalination, 2020, 479: 114314.

[20] Wang B, Tian R. Investigation on flow and breakdown characteristics of water film on vertical corrugated plate wall [J]. Annals of Nuclear Energy, 2019, 127: 120-129.

[21] 闵敬春, 彭晓峰, 王晓东. 竖壁上液滴的脱落直径 [J]. 应用基础与工程科学学报, 2002, 10(1): 57-62.

Min J C, Peng X F, Wang X D. Departure diameter of a drop on a vertical plate [J]. Journal of Basic Science and Engineering, 2002, 10(1): 57-62.

Study on wave-plate demister with internal transport wall

PENG Yu-feng, CHEN Zhuo, ZHU Jia-hua

(School of Chemical Engineering, Sichuan University, Chengdu 610065, China)

Wave-plate demisters made of internal transport materials can facilitate demist performance by water transportation. Droplet separation, water transportation and re-entrainment prevention mechanism were investigated by establishing corresponding models. A semi-empirical droplet separation efficiency prediction model was established and validated by experiment data, and the relationship between falling-film thickness and material properties was obtained. The re-entrainment prevention mechanism of the internal transport wall was analyzed. The results show that traditional demisters may form re-entrainment when gas velocity exceeds critical value of 4.82 m×s-1. However, the wave-plate demisters with internal transport wall have no liquid film under such conditions. Re-entrainment happened at gas velocity over 19.52 m×s-1. Droplets can be completely spread on internal transport wall, and the thickness of the liquid film can be effectively reduced to prevent re-entrainment. Therefore, demister performance can be improved.

liquid transport; mist; demister; liquid film

1003-9015(2021)06-0979-07

TQ 028.8

A

10.3969/j.issn.1003-9015.2021.06.004

2021-03-17；

2021-05-29。

国家科技支撑计划(2013BAC12B01)。

彭玉凤（1996-），女，四川内江人，四川大学硕士生。

朱家骅，E-mail：jhzhu@scu.edu.cn