基于复合结构的两级空化发生器设计优化及数值模拟

2020-04-17杨兴林石园园闫利坤李自强

舰船科学技术 2020年2期

杨兴林，石园园，闫利坤，程林，李自强

(江苏科技大学能源与动力学院，江苏镇江 212000)

0 引言

水力空化现象是指液体内局部压力降低时，液体内部或液固交界面上蒸汽或气体。空泡的形成、发展和溃灭的过程[1–2]。利用空泡溃灭时释放的大量能量，水力空化可以应用在废水处理、强化工业过程及空化水射流等领域。同时在螺旋桨空化噪声，超空泡减阻技术的研究目前也很多[3–4]。空化发生器原理是当水流流过收缩装置时，过流断面急剧变化，水流被增速降压，当压强降到水的气化压强下时发生空化[5]。空化效果的强弱影响空化发生器的效率，确定出最优化的复合空化器的结构参数来强化空化的效果，设计出高效率的空化发生器装置对生产应用具有重要的意义。

目前，国内外很多学者都对水力空化装置进行了探讨研究。例如：王智勇等[6]针对不同压力、不同喉部直径和扩散段长度对文丘里管内部空化现象的影响进行数值模拟；章昱等采用标准k－ω 模型对孔板的水力空化装置进行数值模拟，探讨空化效应强度的影响因素问题；曾章美等[9]对孔板位置一定时，不同的文丘里管喉部的组合情况进行数值模拟，得出急型的复合结构的空化器空化效果更强。但这些研究一般集中于组合空化发生器中的文丘里管结构的研究，其空化效果还有待于提高。

本文结合文丘里管和孔板各自的优点，设计了一种把孔板嵌套在文丘里管喉部位置组成复合结构的两级空化发生器，在管内流量和过流截面积一定时，通过用多个小管径的孔道代替大直径的单文丘里管管道，分别在孔板孔道和文丘里管渐扩段发生2 次空化过程，两级空化发生器中液流经过了2 次加速，从而能够产生更多的空化泡，空化泡不断地产生和溃灭，大大提高了空化效应，以此来提高空化发生器整体的空化效果。两级空化发生器的空化面积和空化维持时间都要更长，从整体上来看，更有利于长时间的作业和强化。且此设计是一体式复合结构空化发生器，制造成本低，结构紧凑，更易运输和节省材料。且本文在所设计的复合空化发生器形式的基础上，分别针对3 组孔不同排布方式孔板与文丘里管喉部进行复合的装置，利用Fluent 软件分别对3 种空化发生器的气液两相流场进行数值计算，对比分析不同形式流场的空化效果，以获得相对较高效率的两级空化发生器结构的优化设计，此数值模拟的结果可为研究多级空化发生器结构设计提供参考信息。

1 数值理论计算

本文装置中的液流视为不可压缩的稳定流动，利用Fluent 软件采用标准K-ε 湍流模型和混合模型（misture mordel），压力梯度采用Standard 格式离散，动量方程和湍流输运方程采用绝对稳定的一阶迎风格式，压力速度耦合采用Simple 方式，在低压区引入空化模型的混合多项流模式对流场进行计算。

1.1 湍流模型

本文计算中，空化发生器中的水流为不可压缩的稳定流动，湍流模型选择应用较广泛的k-ε 两方程模型。标准k-ε 模型基于Wilcoxk-ε，在加入Reynolds 数以及各种流体特性的条件下改善得到的，所以本文求解采用标准k-ε 模型。标准k-ε 模型的连续性方程为：

动量方程为：

式中：ρ 为流体密度；μ 为时均速度；p 为流体压力；η 为流体动力粘度，为紊流应力。

湍流动能方程为：

式中：Gb为由浮力产生的湍动能；YM为可压缩湍流中阶梯的散布所生成的波动；Gk为由层流的速率差生成的湍动能。

湍流耗散方程为：

式中：C1ε，C2ε，C3ε为常数；σk，σε为两方程的Prandtll 数。

1.2 空化模型

考虑到本文研究的对象和体系，综合分析选择Schnerr and Sauer 空化模型。其空化模型表达式为：

2 计算模型及参数

2.1 数学模型

图1 为设计的复合两级空化发生器结构示意图。图中，入口和出口直径D1=20 mm，孔板直径D2=5 mm，孔板小孔直径d=1 mm，渐缩段长度L1=24 mm，孔板长度L2=10 mm，渐扩段长度L3=63 mm，入口锥角α1为40°，出口锥角α2为18°。图2为设计的3 种不同的孔板孔的排布方式，分别为a 环状排布、b 直线排布、c 中心孔辐射状排布。

2.2 网格划分及边界条件

图 1 复合结构空化发生器结构示意图Fig. 1 Schematic diagram of composite structure cavitation generator

图 2 孔板孔的排布方式Fig. 2 Arrangement of orifice holes

图 3 空化发生器的网格示意图Fig. 3 Schematic diagram of the cavitation generators

本文利用CFD 对三维模型进行网格划分，如图3所示。采用非结构性网格进行计算，对孔板孔道部分进行加密处理，通过网格检查，确保网格质量满足要求。本文选取25 ℃水和25 ℃的水蒸气为对象进行数值模拟，设置进口条件为速度等于0.5 m/s，出口条件为压力0 Pa，无滑移标准壁面函数边界条件。水的密度为998 kg/m3，表面张力为0.717 N/m，设定不可压缩相为水。

3 计算结果及分析

通过使用Fluent 软件对流场进行数值计算，使用Tecplot 软件进行后处理。

3.1 速度分布

图4 为3 种不同复合结构空化发生器的速度分布云图。由图可知，3 种空化发生器的速率都是先增大，在孔板孔口处速率达到最大。由于流体惯性，通过孔板孔道后流体依然具有较高的流速，从而形成射流，整个过程表现为速度分布由中心轴线向管壁转移的过程。而在文丘里管渐扩段的锥角影响下，水流的流动区域变大，使水的流动速率逐渐变小，在管道表面周围水流受到摩擦力作用流体流速开始逐渐降低。不同的孔板排布对空化发生器的速度场分布有显著的影响。3 种装置的核心区的平均速度是c＞a＞b，这是由于中心孔辐射状排布的空化发生器c 有中心孔，总体的速度更容易叠加，使流体受到的阻力较小，从而可以将更多地能量转化为空化效应的势能，所以核心区的速度相对最高，径向范围的速度场也最大；而复合结构空化发生器b 和发生器c 相比没有中心孔，所以整体的速度场分布范围和整体核心区的平均速度相对于c 范围略微减小；直线分布的空化发生器b 由于沿直线均匀分布，孔间隔较大，水流在流出孔道后的射流不集中，相对于其他两种装置核心区的水流速度相对较小。轴向方向速度矢量存在负值，说明有回流区，这是因为高速水流在文丘里管壁面阻挡所致。

图 4 空化发生器的速度云图Fig. 4 Cavitation generator speed cloud

3.2 湍动能分布

图5 为3 种复合结构空化发生器的湍动能云图分布。从图中可以观察到，水流在进入孔板孔道前的流场湍动能很小，在流进孔道后迅速增大后又减小，在文丘里管渐扩段部分，湍动能再次增大。空化产生期间空泡的生成和破灭导致了流场的振荡，湍动能也同时进一步改变，在此过程中湍动能有2 次幅值的变化，即发生了2 次空化过程。由图可知，不同的孔板排布方式对空化发生器的整体湍动能散布有显著的影响，其中空化发生器c 的湍动能相对于其他2 种空化发生的散布区域更广和湍动能更强。这和前面的速度场分布有类似的规律。

3.3 汽含率分布云图和散点图

图 5 空化发生器的湍动能云图Fig. 5 The kinetic energy cloud diagram of the cavitation generator

图 6 空化发生器的汽含率云图Fig. 6 Cavitation rate of the cavitation generator

图6 和图7 分别为3 种空化发生器的汽含率云图和分布散点图。可以清晰地看出汽含率的分布情况，其中图7 的横坐标是沿管轴的距离，孔板位于28～38 mm 范围内。可知，3 种两级复合空化发生器的汽含率分布形状基本相似，3 种复合结构空化器在Y 轴坐标上汽含率都呈现出2 次最高值，说明产生了2 次空化。汽含率较高点主要产生在孔道区域（28～35 mm）和文丘里管锥角扩散段范围内（38～64 mm），表明空化大部分是产生在这些区域中，这与前面分析过的速度场和湍动能分布规律相似，汽含率的最高点都产生在文丘里管锥形段的流场中，主要集中在50%～97%左右。从图6 和图7 可以看出，空化发生器c 的最高汽含率大概是99%，说明核心区的部分区域几乎完全空化，核心区的平均汽含率明显高于其他2 种空化发生器。而且对比3 种空化发生器的空化区域，很明显的空化发生器c 的空化区域范围最长，在28～35 mm 和38～64 mm 之间。说明不同的孔板分布结构对于复合式空化发生器的空化场有较显著的影响，3 种两级空化的复合空化发生器中带中心孔的辐射状排布的孔板复合结构空化发发生器的空化效果最好，这也和前面的速度场和湍动能的分析结果一致。