高炉鼓风脱湿系统除雾器挡板结构两相流场数值模拟

罗雨慧　谢建中　左可

摘要：

对高炉鼓风脱湿系统除雾器的除雾机理进行了详细分析，在此基础上参考实际运行参数，采用计算流体力学（CFD）方法对除雾器挡板结构内的气液两相流动进行了数值模拟，得到了液滴的运动轨迹以及液滴质量浓度、压力、速度和旋涡分布情况，并对除雾器挡板结构内部两相流场进行了深入分析.结果表明：鼓风脱湿系统除雾器能够对携雾气流中的液滴进行有效的分离，从而保证鼓风脱湿系统连续可靠地运行.研究对除雾器的优化设计和运行具有指导意义.

关键词：

除雾器； 鼓风脱湿； 两相流； 压降； 流场

中图分类号： TQ 051.8文献标志码： A

Abstract：

Based on the analysis of the demisting mechanism of the demister in the blower dehumidification system of a blast furnace，the computational fluid dynamics（CFD） method was used to simulate numerically the twophase flow of the gas and liquid in the baffle structure of a demister with the actual operation conditions as a reference.The detailed distributions of the pressure，velocity，vortex，and droplet concentration were provided.The simulation also showed the tracks of the droplets.In addition，an indepth analysis and research on the twophase flow were conducted.The results showed that the demister could separate the droplets from the airflow effectively，so as to ensure that the blower dehumidification system could run continuously and reliably.This study is helpful to the optimization and operation of the demister.

Keywords：

demister； blast dehumidification； two phase flow； pressure drop； flow field

钢铁企业高炉鼓风站的电能消耗十分巨大，其中，鼓风脱湿系统中除雾器阻损对鼓风机的电能消耗影响较大[1].高炉鼓风系统设置的除雾器主要用于清除由鼓风冷凝时产生的游离雾滴，提高鼓风质量.

本文拟对鼓风脱湿系统中丝网与挡板组合式除雾器的分离机理进行详细分析，在此基础上对除雾器挡板结构内部气液两相流动进行数值模拟，得到流道内液滴的运动轨迹及压力、速度、旋涡分布，为除雾器的进一步优化设计提供指导.

1除雾器的结构及工作原理

在脱湿除雾器中，挡板主要起导流作用，究其原因在于多层丝网结构的除雾效率很高（根据现场测试数据，对于粒径大于1 μm的粒子，分离效率高达99%；对于粒径大于2 μm的粒子，分离效率高达99.6%）.携雾气流流过丝网结构后，其中所夹带的液滴几乎全部被捕集除去.液滴在重力作用下沿丝网的丝径向向下运动，同时继续吸附气体中夹带的雾滴，凝聚变大的雾滴滴落在挡板上，沿挡板流入水槽中，进而排出除雾器外，实现除雾.

本文所研究的鼓风脱湿系统中除雾器布置在脱湿器内二级冷却器之后，为压板固定的框架式8层×8块钢丝网结构（7块斜挡板）.除雾器除雾构件如图1所示.在除雾器中，携雾气流流过钢丝网层实现气液分离，而被丝网捕集分离和与挡板撞击分离的液滴则经挡板导流引入水槽中，进而被排出除雾器外.

携雾气流在除雾构件内流动，由于流线偏折，在惯性作用下，液滴不能随气流偏转而撞击到挡板上，其中：动量较小的液滴粘附在挡板表面被捕获；动量较大的大液滴撞击挡板表面发生溅射，产生多个小液滴.聚集在挡板表面的水滴受三种力的作用，即拉力、重力和表面张力.当重力占主导地位时，液滴在重力作用下沿挡板表面流入集液槽排出.气流冲刷挡板表面的液膜，将其卷起、带走.撞击在挡板表面的液滴由于自身的动量过大而破裂、飞溅，均可能导致雾沫的二次夹带.

2模型建立及网格划分

携雾气流在除雾器内的流动实际上是一种可压缩的黏性流体的三维、非定常的复杂流动.对这一实际流动情况无法采用数学形式精确描述.本文从既能较好反映实际情况又力求模型构造简单角度出发，在合理误差范围内对流场作简化：

（1） 由于气流速度较小，马赫数远小于0.1，故可把气体视为不可压缩气体处理[2]；

（2） 考虑在实际的稳定工作条件下，流动参数与时间的关系及气流的振荡对流场的影响可忽略，故将流动视为定常流动；

（3） 采用冷态条件，过程中不考虑温度的影响，不考虑气液传热和液体蒸发等现象；

（4） 由于液滴粒径很小，故可作球形处理，考虑流动中其重力和气相对液滴的曳力；同时假定其在运动中直径不变，不考虑液滴之间的碰撞、聚合等现象，且忽略蒸发、摩擦、撕裂及热效应的影响，不考虑气液两相之间的任何能量交换；

（5） 液滴无溅射，不考虑壁面反弹、液膜形成与撕裂，忽略二次带水的影响；

（6） 液滴接触挡板，即认为被捕集；液滴到达除雾器出口时，即认为液滴逃逸.

本文采用计算流体力学（CFD）应用软件Fluent进行数值计算，应用其前处理软件Gambit生成网格，并参考以往的数值计算结果和经验[3]，整个计算区域采用非均匀的网格布置方式.为了提高整体网格质量，网格划分时先对除雾器部分进行网格加密，采用内部面将除雾器所在体与其他体分离.

将模型分割为四个计算区域：方管段、除雾器、渐缩段和圆管段.先划分除雾器的相关网格，在挡板结构上生成线网格，采用相似边界软连接，再采用混合结构网格从已有边界网格生成挡板结构体网格，最后生成除雾器的体网格.除雾器简化三维模型如图2所示.其四个计算区域的网格划分均采用混合结构网格，共生成75 965个节点，353 251个混合形式网格.除雾器网格划分如图3所示.

3数值计算方法及边界条件

在除雾器内流场的数值模拟中，通道内包括互相之间交换质与能的三相——气相、气流夹带的液滴相和挡板表面的液膜相.本文采用两相流模型，只考虑气相和液滴相.

对于连续相（气相），由雷诺数可知流动为湍流.本文采用SIMPLE算法进行压力-速度耦合，采用有限体积法对算例进行离散处理.压力采用Standard离散格式，动量、湍流动能、湍流耗散率均采用二阶迎风离散格式，以获得较准确的解.

对于离散相（液滴相），由于除雾器中气流内液滴的体积分数小于10%，故可忽略颗粒之间的相互作用、颗粒对气相的影响、颗粒的运动轨迹[3-5].本文选用基于LagraianEulerian法的DPM（discrete phase model）模型，按拉格朗日方法对各个颗粒方程进行积分求解.

交替求解离散相与连续相的控制方程，直到两者均收敛从而实现离散相与连续相双向耦合.

携雾气流中液滴流动轨迹计算时只需考虑阻力和重力的影响[6].

气相作为连续相在欧拉坐标系中描述；液滴相作为离散相在拉格朗日坐标系中描述[7].初始条件及边界条件分别为：

（1） 连续相（气相）

介质：空气，密度为1.205 kg·m-3，动力黏度（20℃）为18.1 μPa·s.

进口条件：给定气流速度uy=1.2 m·s-1，ux=0，假设其在进口截面为均匀分布，湍流度为0.05[8].

出口条件：自由出流.

壁面条件：无滑移，绝热.

（2） 离散相（液滴相）：

介质：水，呈细小液滴状，密度为1.0×103 kg·m-3，平均液滴直径为5 μm，鼓风含湿量为10 g·Nm-3.

进口条件：设定液滴速度与气流速度相同[9]，喷射类型选为表面，使液滴在进口截面均匀分布.

壁面条件：选择捕集类型，即不考虑反弹，液滴触及壁面即认为被捕集，不考虑二次夹带效应.

当计算连续性残差、速度残差、湍流动能残差、湍流耗散率残差均降至10-3，且入口与出口气相流量相差小于3%时认为计算收敛[10].

4流场分布及分析

4.1液滴质量浓度分布

图4为计算得到的液滴质量浓度c的分布.由图可知：在捕集液滴的除雾构件区域及起导流作用的挡板附近液滴质量浓度明显较高，而在整个除雾构件中第四至第六块斜挡板之间区域的液滴质量离效果已足以保证系统连续、可靠地运行.

4.2液滴颗粒的运动轨迹

采用Fluent软件对除雾器进行模拟，得到的连续相和离散相的运动轨迹如图5所示，其中图5（b）为计算得到的除雾器内液滴的随机运动轨迹.由图5（b）可知，液滴运动轨道主要集中在除雾器除雾构件中的中间五块挡板部分区域，这对于除雾器的工程优化设计有着重要的指导意义.此外，液滴主要在除雾器中被捕集，小部分进入渐缩段与壁面发生碰撞而被分离.被捕集的液滴主要源自携雾气流的主流两侧，由于受到气流速度梯度影响，其在流线偏折时所受离心惯性作用较大，容易发生偏转，撞击在壁面上而被捕集.故除雾器入口处液滴参数的设计对除雾效率的提高具有较大影响.

由图5（b）亦可知，鼓风脱湿系统中的丝网与挡板组合式除雾器可有效地实现气液分离.

4.3压力及速度场分析

除雾器中压力和速度分布云图如图6所示.由图可知，方管段和圆管段压力、速度分布较为均匀.携雾气流流过除雾器挡板结构的过程中未出现显著的压降，而流过除雾器出口处的渐缩段时发生了明显的压降，降幅达142 Pa左右.其产生的主要原因是流体的流通面积减小，并且截面形状突然改变，流体流线被迫发生改变[11].由于旋流作用较强，在渐缩段方圆突变截面后出现了明显的低压区，在渐缩段后方圆突变截面处及其后一段距离内亦出现了显著的压降.

4.4旋涡分布

脱湿除雾器内的湍流动能和湍流耗散率分布如图7所示.在整个脱湿除雾器中，湍流耗散最强烈的区域为渐缩段后一段距离.在除雾器除雾构件中，中间五块挡板部分区域的湍流动能最大，湍流耗散强烈，是除雾器内实现气液分离的关键区域，携雾气流流过除雾器通道时，气流主要流道在该区域.由于惯性力的作用，气流中液滴的跟随性变差，速度迟豫时间延长，所以易于碰到壁面而被捕集[12].在第三、四块挡板下部区域的高压低速区是除雾器挡板结构中湍流动能最大、湍流耗散最强烈的区域.该区域除雾效率较高.气体流过除雾器后在最上及最下处竖直挡板后易形成回流区，此处亦可能产生角涡，加剧了流场的扰动，增大了系统压力损失.上部竖直挡板附近流域局部放大图如图8所示.

5结论

本文运用两相流模型对鼓风脱湿系统除雾器挡板结构的主要性能进行了数值模拟，得到了液滴质量浓度、压力、速度及旋涡分布情况，展示了其三维内部流场，揭示了实验手段难以获得的数据和现象.计算结果表明：鼓风脱湿系统除雾器能够有效实现携雾气流中液滴的分离.模拟结果对除雾器结构设计、降低除雾器内流动阻力具有一定的参考性意义.

参考文献：

[1]张永红，袁熙志，罗冬梅，等.我国钢铁行业节能降耗现状与发展[J].工业炉，2013，35（3）：12-17.

[2]周光炯，严宗毅，许世雄，等.流体力学[M].2版.北京：高等教育出版社，2006.

[3]李鹏飞，徐敏义，王飞飞.精通CFD工程仿真与案例实战[M].北京：人民邮电出版社，2011.

[4]孙秀君，孙海鸥，姜任秋.油气分离器内油滴轨迹的数值模拟[J].应用科技，2006，33（10）：69-72.

[5]李文艳，徐妍，申林艳.数值模拟技术在湿法脱硫除雾器优化设计上的应用[J].热力发电，2007，29（5）：10-15.

[6]韩占忠，王敬，兰小平.FLUENT——流体工程仿真计算实例与应用[M].2版.北京：北京理工大学出版社，2010.

[7]郭鹏宇，杨震.除雾器通道内二维两相流场的数值模拟[J].电力科学与工程，2005（2）：31-33.

[8]王福军.计算流体动力学分析[M].北京：清华大学出版社，2004.

[9]赵毅，华伟，王亚君，等.湿式烟气脱硫塔中折线型挡板除雾器分离效率的数值模拟[J].动力工程学报，2005，25（2）：293-297.

[10]陈凯华，宋存义，李强，等.湿法烟气脱硫系统中折板式除雾器性能的数值模拟[J].环境工程学报，2007，1（7）：91-96.

[11]CRANE工程部编.流体流经阀门、管件和管道的流体计算TP410[M].北京：化学工业出版社，2013.

[12]RUDINGER G.气体-颗粒流基础[M].张远君译.北京：国防工业出版社，1986.