柴油机舷侧排气喷淋流动及传热模拟

2019-10-24王振吴炜杨先勇安一峰刘元春

中国舰船研究 2019年5期

王振，吴炜，杨先勇，安一峰，刘元春

中国舰船研究设计中心，湖北武汉430064

0 引言

舷侧排气是一种常见的排气方式，许多中、小型船舶的推进柴油机及大型船舶的发电柴油机采用了舷侧排气。针对舰船，可以将舷侧排气系统的排气口设置在水面附近，船体和波浪能够在一定程度上对排气起到遮蔽作用，有利于红外隐身，但由于排气温度较高（可达到500℃）［1-3］，若不采取降温措施，正对着排气口方位的舷侧排气口仍然是非常强烈的红外辐射源，甚至在某些特定的方位比烟囱内向上直排的排气红外辐射更加强烈，因此，有必要针对其排气降温技术进行研究。

海水喷淋冷却技术是一种降低船舶排气温度的高效且经济的方式，通过向排气管内喷射具有一定压力和流量的海水，利用海水气化潜热和海水与排气的掺混换热，可以显著降低排气温度。喷淋降温技术在压气机进气道降温［4］、火箭发动机尾焰流场降温［5］、高压喷水冷却［6］、排气冷却消声［7］等多个领域都有应用。但针对船舶柴油机舷侧排气喷淋效果及对排气管路阻力影响的研究尚不多见。

喷淋降温的效果评估是喷淋冷却技术应用的基础，随着计算流体动力学（CFD）技术的发展，数值计算方法被应用于对该类问题进行数值模拟和评估［8］。赖孝君等［9］对导弹舱室固体火箭发动机意外点火的喷雾降温开展了数值模拟。闫丹丽等［10］利用Mixture 多相流模型和沸腾传质模型对废热锅炉两相流区的沸腾传热进行了数值模拟。黄亮［11］利用FLUENT 软件对潜艇通气管状态下柴油机排气喷淋降温效果开展了仿真研究。

本文拟针对一种结构紧凑、适用于布置空间小的舷侧排气系统的多孔射流式排气喷淋装置，应用Mixture 多相流模型和蒸发冷凝模型对其排气管内的喷淋降温过程进行数值模拟，分析不同喷射流量对流动和传热特性的影响，为排气喷淋装置的设计和传热特性评估问题提供依据。

1 数值计算方法

1.1 柴油机排气喷水冷却过程物理模型

根据船舶柴油机舷侧排气的常见布置型式，在排气管末段，较长的直管段部分可能为水平、倾斜、竖直3 种状态。以某型船柴油发电机排气管的实际布置为例，选取竖直排气管段设计排气喷淋装置。为简化计算，不考虑出舷侧的弯管段，简化后的物理模型如图1 所示，图中H 为装置的高度，D 为排气管直径。排气管内为柴油机排出的具有一定流量的高温废气，排气管上开设多个喷淋孔，冷却水在压力作用下通过喷淋孔射入排气管内形成膜状喷雾流；与高温烟气混合后，冷却水发生蒸发，通过气化潜热吸收大量的热量，从而降低排气管内烟气的温度。由于大量冷却水气化，烟气中的水蒸气含量大幅度增加，最终将含有大量水蒸气的低温气体排出。

图1 多孔射流式排气喷淋装置原理图Fig.1 Schematic diagram of multiple jet spray device for exhaust

1.2 排气喷淋冷却气液两相流数学模型

排气喷淋过程涉及了两相流和相变传热，由于Mixture 模型可较好地处理高温、高压、可压缩性气、汽、液多相流动复杂问题，因此本文选择Mixture多相流模型及蒸发—冷凝模型来进行计算。

本问题中，气、液两相的流动速度都不高，故可假定气、液两相为不可压缩流体，其有限体积法的控制方程如下所示。

连续性方程：

动量方程：

能量方程：

蒸发—冷凝模型：

蒸发—冷凝模型是机理模型［12］，适用于Mixture 模型。其液—气质量输运可以通过蒸汽输运方程描述为

式中：下标“v”为水蒸气相；下标“l”为液相；αv为水蒸气相体积分数；ρv为水蒸气相密度，kg/m3；Vv为水蒸气相速度，m/s；m˙l→v和m˙v→l分别为蒸发和冷凝引起的质量传递速率，kg·s-1·m-3。

本文基于FLUENT 平台进行数值求解，其中湍流模型选用RNG k-ε模型，壁面处理采用增强壁面函数。

1.3 计算模型

1.3.1 几何模型

计算模型及计算域如图2 所示。计算域直径为排气管直径的10 倍，计算域长度为排气管长度的6.7 倍。坐标原点设在烟气入口的中心，其中坐标系z轴与排气管轴线重合，方向与重力加速度相反。鉴于实际海水中含有大量杂质，喷孔过小容易造成堵塞，综合考虑了防堵和冷却效果，将喷淋孔直径设计为2 mm。

1.3.2 网格划分

采用非结构化网格，分内、外2 个流域，网格总数约1.7×106，对排气管壁面及喷淋孔附近的网格进行了局部加密（图3）。

图2 计算模型及计算域（y=0 剖面）示意图Fig.2 Schematic diagram of computation model and computation domain（y=0 cross section）

图3 网格划分Fig 3 Mesh generation

1.3.3 边界条件

边界条件设置如表1 所示，其具体位置如图2所示。

表1 边界条件Table 1 Boundary conditions

2 结果与讨论

2.1 网格无关性分析及模型验证

为验证本文计算结果与网格数量的无关性，分别选取了网格数量为8.1×105，1.7×106，3.5×106的3 套网格进行计算。在除网格外的其他条件保持一致的情况下，比较了3 套网格条件下排气管出口的平均温度，结果如图4 所示。由图可见，计算结果基本保持不变。因此，本文计算选取的网格数量1.7×106能够满足网格无关性要求。

图4 不同网格数下排气管截面平均温度沿程分布曲线Fig.4 Distribution curves of average temperature along the cross section of exhaust pipe with different meshes

为验证本文数值模拟方法的准确性，采用该方法对文献［13］中的气液两相管喷淋冷却实验进行数值模拟，模拟得到的冷却后气体平均温度为312.65 K，与实验值较为接近，从而证明了本文方法的准确性（表2）。

表2 数值模拟与实验结果对比Table 2 Comparison between numerical simulation and experimental results

2.2 不同流量比下的温度场分布

图5 不同冷却流量比下y=0 剖面的温度分布Fig.5 Temperature contours of y=0 cross section with different cooling flux ratios

定义冷却水质量流量ql与排气质量流量qg之比为冷却流量比ξ，ξ=ql/qg。ξ值的大小表示冷却水量的相对量。本文分别取ξ=2.9，1.4，1.0，0.7 进行数值模拟，得到了流场和温度场，取y=0剖面的温度场进行分析，结果如图5 所示。由图可见，该装置存在锥形排气高温核心区，且随着ξ的增大，排气高温核心区的范围逐渐缩小。

2.3 沿程温度分布

为便于描述沿排气管轴向温度变化的情况，定义χ=l/l0，其中l为排气管某截面距排气进口的轴向距离，l0为排气管总长，χ为无量纲沿程距离。图6 为排气管截面平均温度随无量纲沿程距离的变化曲线，由图可看出，在χ=0.1～0.4，排气管截面平均温度急剧下降，之后排气管截面平均温度变化趋缓。χ=0.1～0.4 区域大致相当于喷孔区，这证明了喷淋降温的效果显著。此外，降温效果随着冷却流量比ξ的增大而变好。

图6 不同冷却流量比下排气管截面平均温度沿程分布曲线Fig.6 Distribution curves of average temperature along the cross section of exhaust pipe with different cooling flux ratios

2.4 阻力特性

为评估排气喷淋阻力特性对柴油机性能的影响，将不同冷却流量比ξ下排气管路总压损失与无喷淋、光滑排气管的理论沿程阻力损失估算值［14］进行了比较，结果如表3 所示。由表可见，喷淋引起的排气阻力并不会随着冷却流量比的增大而单调增大；喷淋引起的阻力损失与沿程阻力损失大致相当，一般柴油机排气管路总排气阻力为千帕量级，喷淋引起的总压损失几乎可以忽略不计，因而对柴油机性能影响很小。从表中还可以看出，当冷却流量比较大，喷淋所增加的动能大于气液碰撞等造成的压差阻力和摩擦阻力时，喷淋造成总压损失还可能为负值，也即喷淋反而降低了排气阻力。