1、SHIP ENGINEERING 船 舶 工 程 Vol.45 No.2 2023 总第 45 卷,2023 年第 2 期 99 文丘里管对洗涤塔流场影响的数值模拟文丘里管对洗涤塔流场影响的数值模拟 魏 凯1,2,龚泽儒1,李 冬1,邹 鹏1,朱元清2(1烟台龙源电力技术股份有限公司,山东烟台 264006;2哈尔滨工程大学 烟台研究院,山东烟台 264006)摘 要:摘 要:为了研究文丘里管对洗涤塔内部流场的影响,利用 Fluent 软件分别对文丘里 U 型塔与普通U 型塔的内流场特性进行数值模拟,为洗涤塔的结构优化提供依据。结果表明:文丘里管对洗涤塔速度场、温度场以及压力场影响较大;文丘里
2、U 型塔内反应温度较低,有利于气液之间高效反应;文丘里 U型塔流场在高负荷下分布更加均匀,“冲壁”现象减弱;文丘里 U 型塔压力损失始终高于普通 U 型塔,并且在文丘里管中,喉部与渐扩段压力损失较高,约占文丘里管压力总损失的 87%。关键词:关键词:洗涤塔;文丘里管;数值模拟;速度场;温度场;压力损失 中图分类号:中图分类号:U671.99 文献标志码:文献标志码:A 【DOI】10.13788/ki.cbgc.2023.02.13 Numerical Simulation on Influence of Venturi Tube on Flow Field of Scrubber Tower
3、 WEI Kai1,2,GONG Zeru1,LI Dong1,ZOU Peng1,ZHU Yuanqing2(1.Yantai Longyuan Power Technology Co.,Ltd.,Yantai 264006,Shandong,China;2.Yantai Research Institute,Harbin Engineering University,Yantai 264006,Shandong,China)Abstract:In order to study the influence of Venturi tube on the internal flow field
4、of the scrubber,the Fluent software is used to carry out numerical simulation research on the characteristics of the Venturi U-type column and the ordinary U-type column in order to provide basis for the structural optimization of scrubber.The results show that the Venturi tube has a great influence
5、 on the velocity field,temperature field and pressure field of the scrubbing tower,and the reaction temperature in the Venturi U-type tower is lower,which is conducive to the efficient reaction between gas and liquid.The Venturi U-type tower of the flow field is more uniformly distributed under high
6、 load,and the phenomenon of rushing the wall is weakened.The pressure loss of the Venturi U-type tower is always higher than that of the ordinary U-type tower,and the pressure loss of throat and diffuser is higher in Venturi tube,accounting for about 87%of the total pressure loss of Venturi tube.Key
7、 words:scrubber tower;Venturi tube;numerical simulation;velocity field;temperature field;pressure loss 0 引言引言 由于船舶柴油机排放的废气对全球环境造成了严重影响,国际海事组织(International Maritime Organization,IMO)制定了严格的 NOx和 SOx排放标准1-2。为了能够同时满足NOx和 SOx的排放限值,可采用低硫油联合选择性催化还原(Selective Catalytic Reduction,SCR)技术、SCR 联合废气洗涤技术、湿法氧化脱硫脱
8、硝技术等。考虑船舶内空间有限以及经济效益,利用洗涤塔进行湿法氧化脱硫脱硝技术成为目前的研究热点之一3。目前,洗涤塔已经广泛应用于电厂与船舶的尾气处理领域中,国外已经对其进行了大量的试验与模拟仿真研究。BROWN 等4对洗涤塔内喷嘴的结构与位置进行试验与数值模拟,研究结果表明,三喷嘴洗涤器中洗涤液消耗量少、脱除效率高,但会产生洗涤液聚集现象,可能会造成设备腐蚀;WANG 等5采用ANSYS CFX 软件建立了不同喷淋层组合方式下洗涤塔气液两相流场的三维数学模型,根据流场分析,喷雾高度越高的液滴在吸收塔中移动距离越长,气液接触时间越长,越有利于气液之间传热。收稿日期:2022-04-05;修回日期
9、:2022-08-02 基金项目:国家自然科学基金资助项目(U1906232)作者简介:魏 凯(1998),男,硕士研究生。研究方向:船舶废气污染物处理。船舶动力、推进装置和辅机设备 100 在国内,也有不少专家学者对洗涤塔内部流场特性进行研究分析。郭浩等6模拟了不同喷雾条件对洗涤塔内液滴的蒸发和流动特性的影响,研究结果表明:喷嘴锥角较大时,阻力较小,烟气易穿透喷雾区;喷淋温度高时,液滴在洗涤塔内分布会变差;贾振宇等7研究了喷淋液滴直径对塔内气液两相流的影响,结果表明:喷淋液滴直径较小难以贯穿整个洗涤塔,会导致塔内速度分布不均匀,影响脱除效率;CHEN 等8对导流洗涤塔的性能进行了试验和数值研
10、究,观察到导流板可以具有废气整流效果,促进了洗涤塔流场均匀性,解决了由于单一入口造成的偏流问题。洗涤塔可分为 I 型洗涤塔和 U 型洗涤塔,I 型塔结构相对简单,占地面积小,但难以满足多污染物处理的需求;U 型塔是在 I 型洗涤塔之前加装文丘里洗涤器,可有效去除颗粒物并降低废气温度9。洗涤塔内部流场的特性能够直接影响洗涤塔的脱除效率,目前对于洗涤塔的研究多集中于喷淋策略和内部结构等因素对塔内部流场的影响,而文丘里管对于洗涤塔内部流场影响的研究相对较少。本文分别建立带文丘里管U 型塔和普通U 型塔的三维模型,利用Fluent 仿真软件对其进行数值模拟分析,重点研究文丘里管对于洗涤塔内部流场流动特
11、性的影响,为洗涤塔的优化研究提供依据和思路。1 物理模型物理模型 为了提高运算速度和精度,对模型加以简化,去掉外支撑板等结构,分别建立如图 1 所示的文丘里 U型塔与普通U 型塔的三维模型。在2 种 U 型塔中,烟气依次经过第一喷淋层、第二喷淋层以及第三喷淋层。第一喷淋层采用单喷嘴布置,第二喷淋层和第三喷淋层均为 5 个喷嘴。为了显示文丘里管不同区域的温度及压力损失变化情况,在文丘里管U 型塔中沿文丘里管轴向选取10 个截面,即 y1y10,在普通U 型塔同一位置也选取 10 个截面,将简化后的两模型导入ICEM CFD 中采用非结构化网格划分。经过网格无关性验证,两模型最终网格数量约79 万
12、。图1 洗涤塔物理模型 2 数学模型及假设数学模型及假设 2.1 数学模型及边界条件数学模型及边界条件 洗涤塔中涉及柴油机高温废气与洗涤液间的两相流动,本文将气相视为连续相,液相视为离散相,采用欧拉-拉格朗日方法计算气液两相间的质量、能量和能动量传递过程,可以较好地实现两相之间的耦合。选取在模拟流体旋转流动性更加精确的RNG k-模型来获取洗涤塔内气液两相流动规律,具体方程如下。质量守恒方程:)_)_)_0uvwtxyz+=(1)式中:为密度;t 为时间;u、v 和 w 为 x、y 和 z 方向上的速度分量。能量守恒方程:Tppp()()()()TuTvTwTtxyzkTkTkTSx cxx
13、cyx cz+=+|(2)式中:T为系统温度;cp为比热容;k为传热系数;ST为黏性耗散。动量守恒方程:()div()yxxxzxxupuuFtxxyz+=+(3)()div()xyyyzyyupvuFtyxyz+=+(4)第三喷淋层 第二喷淋层 第一喷淋层第二喷淋层第三喷淋层第一喷淋层(a)文丘里 U 型塔(b)普通U 型塔魏凯等,文丘里管对洗涤塔流场影响的数值模拟 101 ()div()yzxzzzzupwuFtzxyz+=+(5)式(3)式(5)中:p 为流体微元上压强;xx、yx和 zx为分子黏性作用产生的黏性应力分量;Fx、Fy和 Fz为 x、y 和 z 3 个方向上的单位质量力。R
14、NG k-湍流方程:)_bM()iiijkikkkkktxxxGGYS+=+|+-+(6))_)_213b2()iiijjktxxCGC GCRkkx+=+|+|-(7)式(6)和式(7)中:i、j和k为常数;和i分别为层流与湍流的黏性系数;Gk为由层流速度梯度产生的湍流动能;Gb为由浮力产生的湍流动能;和k分别为k方程和方程的湍流Prandtl数;YM为湍流能;Sk为湍流耗散率;C1、C2和C3为经验系数;R为湍流能量耗散率方程的源项。在Fluent仿真计算中,系统入口采用速度入口边界条件,出口采用自由流出口边界条件,壁面采用无滑移、无热传导的边界条件,采用基于压力求解器的Simple算法,
15、对流差分格式采用2阶迎风格式。本文考虑了两型洗涤塔在WD10柴油机不同负荷下的工作情况,具体柴油机工作参数见表1。表1 WD10 柴油机工作参数 参数 100%负荷 75%负荷 50%负荷25%负荷功率/kW 140 105 70 35 转速/(r/min)1 500 1 365 1 200945排气体积/(m3/h)643 526 448 348排气温度/376 350 294 251排气质量流量/(kg/h)268 313 292 244 2.2 合理假设合理假设 本文着重关注文丘里管对塔内流场特性变化的影响,为了提高模拟计算效率,综合考虑多种因素,做如下合理假设:1)假设烟气在洗涤塔内处
16、于稳态流动;2)忽略喷嘴、喷淋管路等内部构件对流场的影响;3)忽略洗涤液之间发生的碰撞、破碎以及液滴蒸发现象;4)不考虑气、液之间的化学反应,洗涤液用液态水代替;5)将洗涤塔内浆液池液面视为壁面。3 结果与分析结果与分析 3.1 文丘里管对洗涤塔温度场的影响文丘里管对洗涤塔温度场的影响 图2为100%负荷下两型洗涤塔的温度云。由图2可看出,柴油机高温废气进入洗涤塔后,首先经过第一喷淋层洗涤,废气温度均有一定的下降。由图3可知:废气在文丘里管中温度下降约23,而在普通U型塔同一位置处,废气温度只下降约16,这是由于文丘里管喉部截面突然收缩,导致废气流速增快,湍动能增大,废气流与洗涤液充分混合,换热比较剧烈,温度下降较为明显;而在普通U型塔同一位置处,只有距离喷嘴较近位置,废气温度有所降低,靠近管壁处废气温度没有变化。随着废气流动,部分洗涤液会附着在壁面上,在两种洗涤塔中,废气与洗涤液换热效果都减弱,温度变化趋于平缓。图2 100%负荷下两型洗涤塔温度云图 图3 100%负荷下文丘里管废气温度随高度变化关系 当废气进入到洗涤塔主塔后,由于主塔内喷嘴较多、喷淋流量增加,气液之间换热更加剧烈,