基于CFD的船舶脱硫系统废水导流装置建模与仿真研究

北极星水处理 讯:摘 要：船用海水脱硫技术是目前世界上技术最成熟、使用最为广泛的 SOx排放控制技术，该技术采用导流稀释的形式处理工艺过程中产生大量酸性较强的脱硫废水。文章采用常用的商用计算流体力学软件，通过对经过导流装置后排出废水的浓度、速度及 pH值进行数值计算，以验证该装置高效的稀释作用。结果表明，该导流装置具有高效的稀释作用，排出废水满足法规排放要求。

关键词：船用脱硫；废水排放；法规；计算流体力学

引言

挖泥船 EGC 系统洗涤废水通过导流装置后与舷外海水进行的混合稀释过程，验证了在船舶静止状态下，距离舷外洗涤废水排口 4m 处海水 pH 值的达标情况。

1 湍流流动模型

1.1 脱硫废水导流装置

1.2 基本控制方程

流体的传热和流动过程都遵守质量守恒、动量守恒和能量守恒这三个基本物理规律，在数值求解流动过程中，基于以上三个基本控制方程由统一的控制方程来计算：

通常采用有限体积法（FVM）对上述基本控制方程进行离散化处理。CFD 计算软件最通用的离散方法为有限体积法（FVM），先将计算域划分为不重叠但有限的控制体积，再对每个控制体的控制方程进行积分计算，得出一组离散的非线性代数方程，最后联立求解得出最终结果。

1.3 湍流流动模型

目前的湍流数值模拟方法可以分为直接数值模拟方法和非直接数值模拟方法，其中主流的计算方法包含：直接数值模拟（DNS）、大涡模拟 （LES） 和雷诺时平均法（RANS）。其中，DNS 方法是对纳维-斯托克斯（N-S）方程进行直接求解，计算结果相对准确，计算量非常巨大，对内存空间和计算速度要求非常高，目前常用于低 Re 数的简单流动计算；LES 方法计算精度稍低，采用特定滤波函数将涡按尺度大小进行分离，并分别按照不同的模型进行求解，这种方法对内存和 CPU 速度要求仍比较高，无法直接应用工程；RANS 方法的核心是不直接求解瞬时的 Navier-Stokes 方程，而是求解时均化的 RANS 方程，极大降低计算工作量，是目前使用最为广泛的湍流数值模拟方法。

基于船舶 EGC 系统废水通过导流装置后具备强烈旋流的流场特征，综合考虑计算量、计算时间以及计算精度，相比于标准 k-ε 双方程模型，本论文采用更适合模拟弯曲流线流动的可实现 k-ε 双方程湍流模型对废水通过导流装置后的湍流流动进行模拟。可实现 k-ε 双方程湍流模型的输运方程如下：

2 络模型与边界条件

2.1 数值计算 格划分

2.2 数值计算边界条件

3 模拟结果及其分析

3.1 速度分布

图 4（a）给出了舷外大空间中心截面的速度分布图。可以看出，排水管道内废水进入舷外后，受到舷外海水的缓冲作用，速度迅速降低至 0.1m/s 左右。同时，如图 4（b）所示，由于导流装置的存在，废水经过导流叶片后产生强烈的旋流，在舷外近排水管道区域形成海水回流，后方回流海水在排水管道出口附近对废水进行快速地混合与稀释，从而提高废水 pH 值。

3.2 温度分布

图 5 所示为舷外大空间中心截面的温度分布图。从温度分布可以看出，高温废水进入舷外后与舷外低温海水进行充分的混合，废水温度在很短的时间和距离内便与周围海水达到了同一水平，从另一个方面也体现了导流装置对增强废水与海水之间相互混合地积极作用。

3.3 pH 值度分布

图 6（a）和（b）为舷外大空间中心截面的 pH 值分布图以及中心线 pH 值变化曲线。从 pH 值分布云图可以看出，废水通过导流装置后产生强烈的旋流，在舷外附近形成了一个低压区，导致后方高 pH 值海水持续回流，对排水管道出口处的低 pH 值废水进行充分地混合稀释，使得当地废水 p H 值快速上升。从中心线 p H 值变化曲线可以看出，舷外 4m 处的废水 pH 值已经达到 7 左右，大于 6.5，符合IMO 对船舶 EGC 脱硫废水的排放要求。

4 结论

本论文采用CFD数值仿真方法对某耙吸式挖泥船EGC系统废水通过导流装置后与舷外海水进行的混合稀释过程进行了模拟，通过分析舷外废水的流场、温度场及p H值分布，论证了在船舶静止状态下，舷外4m处海水pH值的达标情况。

仿真证明，排水管道内装配的导流装置能够有效地增强废水与舷外海水的混合稀释，某耙吸式挖泥船 EGC 系统所产生的废水（流量为 320m3/h，p H 值为 5.5）的 pH 值能够在 4m 内迅速提升到 7 左右，完全满足 IMO 对船舶 EGC系统废水排放指标的要求。

原标题:基于CFD的船舶脱硫系统废水导流装置建模与仿真研究

全国能源信息平台联系电话：010-65367702，邮箱：hz@people-energy.com.cn，地址：北京市朝阳区金台西路2号人民日 社