基于半开式叶轮离心泵内部气液两相不稳定流动特性研究

本文主要研究内容

由于半开式叶轮离心泵内部流动复杂，本文首先对纯水工况下半开式叶轮离心泵进行定常和非定常计算，并对外特性和内部流动状况进行分析，具体内容如下：
（1）计算1000 r/min、1200r/min、1450r/min时的流量-扬程，流量-效率曲线。
（2）进行定常计算，分析半开式离心泵中间截面的压力分布、湍动能分布，以及叶顶间隙的压力分布、流线分布等。
（3）进行非定常计算，不仅要对叶片出口处的监测点进行压力脉动的分析，对于蜗壳断面处也是不容忽视的，此外还要观测其中的流量变化。
气液两相流工况下，基于欧拉两相流模型和Musig模型进行稳态数值计算，并通过气液两相流实验对模型的准确性进行了相应的验证。其中，MUSIG模型加入了粒子形变过程，不仅考虑了入口气体粒子直径，还考虑了气体粒子的聚并和溃灭。
（1）在1000r/min，12.9m3/h（0.75Qd），以及不同含气率下计算两种模型的外特性，将两种模型计算的性能结果分别与实验结果进行比较，对比哪种模型可以更好的反映半开式叶轮离心泵内的流动状况。
（2）为了掌握离心泵内气泡变化引起的气液两相流内部流动规律，分析了不同含气量的半开式叶轮离心泵的气体分布、压力分布、流线分布、和湍动能分布等。
因此，本文还基于MUSIG模型进行了非定常数值计算。具体内容如下：
（1）分析不同时刻半开式叶轮离心泵内气体分布，并与可视化实验数据进行对比。
（2）分析不同含气率下半开式叶轮离心泵压力脉动，并将出口处的压力脉动与实验结果进行对比；

泵内气液两相流气液两相流

气泡在泵内轨迹
许多学者研究了泵中气液两相流的流型，并将其分为四种流动类型：聚合物气泡流、气囊流、均匀气泡流和分离流。

图2. 1离心泵气液两相流流型
Fig.2. 1 Flow pattern in centrifugal pump under gas-liquid two-phase
每个流型的特征和描述如下：
(1)均匀泡状流
当入流含气率比较小时，气体在离心泵内主要以小气泡的形式存在，沿着叶轮流道均匀分散，气泡大多呈球形。
(2)聚合泡状流
随着进气量的增加，气泡的数量更多也更大，使得气泡的间距就减小了，随之而来的是气泡间的相互作用更加明显，在叶轮流道内可观察到气泡的集结和湮灭。
(3)气囊流
随着入流含气率的进一步增加，叶轮内的气泡聚并加剧，形成一些大气泡。大气泡占据了叶轮流道的大面积区域且运动缓慢，限制了流动。不过，这种气泡结构的气液两相界面变形剧烈且不稳定，导致泵内部流动的不稳定，引发“喘振”。
(4)分离流
随着入流含气率的继续增加，之前保持不动的大气泡开始径向扩大尺寸，并逐渐占据整个叶轮流道，从而导致气液相分流，气相形成连续相，几乎导致泵发生空转，这种现象一般称为“气锁”，同时在叶轮通道的进口发现了小气泡聚集，出口区域由也于不稳定性导致大气泡尾部局部破裂，形成具有大量气泡群和强烈回流的区域。

纯水条件下半开式叶轮离心泵实验验证
流量-扬程无量纲化分析
图3.5（a）是叶轮离心泵在一定条件下的扬程曲线，具体条件为转速是、1200r/min和1000r/min时的扬程曲线，从图中可以看出：将转速设置为不同数值的情况下，流量范围与转速呈正比例趋势，在转速不断减小的过程当中，模型泵可运行的流量范围呈规律性减小[39]。
将其无量纲化得到图图3.5（b）所示的扬程系数曲线。从图中可以看出这三种曲线到最后基本上都是处于重合的状态，根据这个图像我们可以得知，相似定律在这种离心泵内，这种情况下依然是适用的。

（a）扬程曲线

（b）扬程系数曲线
图3.5不同转速下半开式叶轮离心泵扬程
Fig.3.5 Heads of semi-open impeller centrifugal pumps at different rotational speeds

基于MUSIG模型不同含气率下流场特性
基于本章4.2节分析半开式叶轮离心泵内流形变化的原因，本节就半开式叶轮离心泵中间截面压力梯度、液相速度、湍动能进行分析。

压力分布分析

半开式叶轮离心泵中间截面压力的分布如图4.10所示，叶轮流道内总压力显示出不对称性，观察图片，我们不难发现，在流体受到叶轮旋转的功的作用下，流道的半径在不断增大，叶轮流道的压力也随之增大，在隔舌位置存在压力突变。在叶片的进口处，吸力面有低压区存在。叶轮的进口这一位置的压力是最小的，在压差的作用下易于吸引并聚集气泡。由于叶轮区域流体的动能转化为蜗壳区域流体的压能，因此蜗壳压力最大可以达到123.5Mpa。在蜗壳出口处和隔舌处，可发现水压可以达到127Mpa。
同时，随着入流含气率的增大，压力分布均呈现逐步下降的趋势，靠近蜗壳的壁面高压区逐渐减小，蜗壳出口管处的高压区面积也逐渐减小。这是由于气相粒子之间聚合与破碎作用消耗能量较高，入口含气率越高压力脉动越剧烈造成的压降现象。

图4.10 1000r/min时不同含气率下半开式叶轮离心泵水体中间截面的压力分布
Fig.4.10 Pressure distribution of the middle section of the water body of the semi-open impeller centrifugal pump with different inlet gas volume fractions at 1000r/min

液相速度分析

图4.11呈现的是在不同的气体含量占比之下，叶轮水体中间部位截面的液相速度的分布情况。从这张图中我们可以得知，半开式叶轮离心泵的气体速度沿叶轮的半径方向呈现增加的趋势，位于叶轮出口的最大速度为7m/s，气体在流道中在压差的作用下，并且可以观察到叶轮流道中有大尺度漩涡和回流现象，这种现象随着含气率的增大越发明显。

图4.11 1000r/min时不同含气率下半开式叶轮叶顶间隙中间截面的液相速度
Fig.4.11 The liquid phase velocity of the middle section of the semi-open impeller under different inlet gas volume fractions at 1000r/mi

不同含气率下气泡分布随时间变化

数值计算结果分析
本节对转速为1000r/min，流量为12.9m3/h（0.75Qd）进行非定常计算，时间步长为，计算总时间为
图5.1为不同进口含气率下叶轮旋转1/6个周期内不同时刻叶轮中截面上的气体分布（T为叶片从某一位置旋转到相邻叶片位置的时间）。可以看出不同时刻下泵内气泡主要集中于叶轮流道出口附近。在含气率GVF=1%时，叶轮流道的气体主要集中于叶轮出口靠近压力面处，在圆周方向呈现非均匀分布，受蜗壳螺隔舌的影响，螺旋线起点下游附近叶轮内的气体不能及时排出，因此靠近隔舌下游叶轮流道内的含气率最高，随着叶轮的旋转，叶轮内流体逐渐从叶轮出口沿蜗壳螺旋线向出口运动，隔舌上游附近叶轮流道内气体体积分数相对较少。蜗壳内气体体积分数最大值出现在蜗壳第八断面附近，叶轮内的气体在隔舌附近排出叶轮[58]。
当含气率GVF=3%时，叶片中间部位靠近吸力面出开始附着有少量气泡，随着叶轮的旋转，叶轮流道内的部分气体在蜗壳第五断面附近开始往下游方向从蜗壳内流出，隔舌附近出现了气体的附着，影响蜗壳喉部过流性能。
当含气率GVF=4.5%时，叶片吸力面的附着气体增加且从出口延伸至叶片进口，占据叶轮流道主要面积的部分就是气体。在叶轮旋转过程中，隔舌下游附近叶轮出口的气体呈聚集状态，当叶轮向螺旋线下游运动时，叶轮出口的气体逐渐分散并往蜗壳内流动。在叶轮从当前位置旋转到相邻叶片位置的过程中，气体将叶轮内部的几乎每一个流道都阻断了，导致在当下情况下离心泵的性能急剧下降、大打折扣。蜗壳内部，气体在第七断面附近聚集严重，从第七断面至蜗壳出口，流道内始终被大量气体占据，存在断流的危险。

（a）IGVF=1%

（b）IGVF=3%

（c）IGVF=4.5%
图5.1不同含气率下不同时刻气体体积分数云图
Fig.5.1 Contour map of air volume fraction at different times under different inlet gas volume fractions

气液两相流下动静干涉作用分析

数值计算结果分析
图5.3为不同含气率下叶片出口处的压力脉动时频域图。在四个含气率下，压力脉动波形相似，压力波动呈非对称分布[62]。
在叶轮的叶片出口这一位置，压力脉动的主频全部是叶频和它的倍频，对这一部位的压力脉动起到重要影响作用的还是叶轮和蜗壳在动与静交替下的干涉。对比纯水工况下叶片出口的压力脉动频谱图可以发现，当泵内通入空气后，叶片出口压力脉动呈现出宽频特性，高频段的压力幅值显著增加[63]。叶频处压力脉动幅值随着入流含气率的增大出现先增大后减小的规律，主要原因是当含气率增加到4.5%时，叶片表面几乎被气体附着，且叶轮流道及出口出现大量气泡形成的气团，叶片对水的作用减弱，使得出口处水流与蜗壳隔舌的干涉作用减小[64]。从二倍叶频到九倍叶频，压力脉动幅值呈现增大的趋势。随着含气率的增加，压力脉动逐渐减小。由此可见，离心泵对入口含气率十分敏感，通入气体后极易改变叶轮内的流场分布，诱发高频脉动与结构振动。

图5.3不同含气率下叶片出口处时频域信息
Fig.5.3 Time-frequency domain information at the outlet of blade under different inlet gas volume fractions
图5.4为半开式叶轮离心泵蜗壳流道各监测点不同含气率时域图，整体来看，压力波动的周期性在入流气体含量占比为4.5% 时与入流气体含量占比为1%时差。当入流含气率比较小时（1%)，各监测点处压力值比较大，而且周期性比较好；当入口体积含气率入流含气率（IGVF）增大到2%以上，各监测点处压力值开始变小。对比不同含气率下的压力脉动可看出：随着含气率的增大，半开式叶轮离心泵蜗壳内相同位置处的压力脉动强度呈现减小的趋势[65]。

图5.4半开式叶轮离心泵蜗壳流道各监测点不同含气率时域图
Fig.5.4 Time domain diagram of different inlet gas volume fractions at each monitoring point of the volute channel of semi-open impeller centrifugal pump
图5.5为半开式叶轮离心泵蜗壳流道各监测点不同含气率频域图，当入流含气率比较小时（1%)，各监测点处压力值比较大;当入口体积含气率入流含气率增大到2%以上，各监测点处压力值开始变小。对比不同含气率下的压力脉动可看出:随着含气率的增大，半开式叶轮离心泵蜗壳内相同位置处的压力脉动强度呈现减小的趋势。产生低频脉动的原因可能是由两相混合流体流出叶轮后撞击蜗壳内壁产生了流动分离，蜗壳内侧气相的堆积使得液相流动受阻，气相分布广泛几乎占据了整个蜗壳流道，并在靠近叶轮出口的区域形成较多漩涡，这使得流动更加复杂并出现低频脉动能量分量[66]。

图5.5半开式叶轮离心泵蜗壳流道各监测点不同含气率频域图
Fig.5.5 The frequency domain diagram of the different inlet gas volume fractions at each monitoring point of the volute channel of the semi-open impeller centrifugal pump

研究总结

本文的设计采用的是数值模拟的手法，并进行对比实验。对于纯水入流和含有气体入流这两种情况下，半开式叶轮离心泵内部的流动情况进行深层次的探究。研究所得的结果比较客观公正且具体的呈现出了纯水入流和含有气体入流这两种情况下半开式叶轮离心泵内部流动情况的不稳定特点和引发种种现象的原因及性质。经过整理以下是本文全部的探究内容以及结论：
1.在控制变量，以纯水入流为基础情况，对半开式叶轮离心泵进行研究，分别实现非定常计算和定常计算。计算的结果呈现出：
（1）叶顶间隙中间截面的压力分布也是随着流量的增大而逐渐增大的，与离心泵中间截面的压力分布十分相似。湍动能分布较大的区域为叶轮流道和蜗壳隔舌附近以及蜗壳的扩散段，随着流量的增大，整体湍动能逐渐减小，局部湍动能较高的地方也逐渐消失。同时，随着流量的增大，叶顶间隙内流动状况得到改善，叶顶间隙涡逐渐减弱。
（2）叶轮与蜗壳的动态与静态的交替干涉导致了流道处于不同的位置下所具备的压力的分布情况是有差异的。因此，对该半开式叶轮离心泵叶轮和蜗壳流道的压力脉动进行了分析，发现：小流量工况下，由于叶轮内流动分离强烈，压力脉动复杂。蜗壳隔舌附近的压力脉动最为剧烈。且随着流量的增加，压力脉动幅值减弱。
2.本文接着运用两种模型分别进行了半开式叶轮离心泵在转速为1000r/min，12.9m3/h及不同入流含气率下的定常数值模拟，通过与模型泵外特性实验和可视化实验结果进行对比，以及分析半开式叶轮离心泵内部流场变化规律，得到了如下结论：
（1）对于此半开式叶轮离心泵，两种模型计算的外特性结果均与实验值相近，但Musig模型计算的外特性结果与实验测得的外特性结果匹配更为良好，且与可视化实验结果更加吻合。说明本章使用并分析的Musig模型、进口边界条件设置的方式是合理的。
（2）导致离心泵的性能大打折扣的最大的影响因素就是入流气体含量占比，入流含气率的变化直接影响叶轮内流型，流型的改变导致泵性能发生变化。入流的气体含量占比不断地上升，叶轮内部气泡的流动形式也会随之改变，均匀泡状流、聚合泡状流、气囊流和分离流动这几种情况都会出现。基于Musig模型分析了不同含气率下的流场特性，气相粒子之间聚合与破碎作用消耗能量较高，入口含气率越高压力脉动越剧烈会造成压降现象。叶轮流道中有大尺度漩涡和回流现象，这种现象随着含气率的增大越发明显。除此之外，还发现了湍动能如何分布最大的影响因素就是叶轮内部是如何分布的。哪个地方的气体含量多，哪个地方的湍动能就越大。这个现象说明叶轮流道内流动不稳定就是因为聚集了大量的气体，使得耗散了更多的流体粘性，进而使得叶轮内部流体在进行能量交换以及能量传递的时候出现了极大的问题，因此模型离心泵的运行性能就大大折扣了。
3.本文最后运用Musig模型进行了半开式叶轮离心泵在转速为1000r/min，12.9m3/h及不同入流含气率下的非定常数值模拟，通过对比半开式叶轮离心泵不同时刻的可视化照片，以及分析半开式叶轮离心泵内部流场变化规律，得到了如下结论：
（1）小含气率下，气泡才是气体在叶轮和蜗壳内部的主要存在形式。入流气体含量占比越来越大，气泡的数量越来越多，体积越来越大，就导致每个气泡与每个气泡之间的间距越来越小，因此气泡与气泡之间的相互作用就更为突出，在离心泵中就会观察到气泡有团聚现象出现。气泡的密度与液体的密度相比，气泡的密度要小于液体的密度，因此气泡的聚集地在叶片吸力面的侧面，由此向出口的部位流动。当气泡到达蜗壳隔舌地方的时候，流动到此的气泡会有一部分是停滞回旋的状态，而另一部分气泡只会随着液体的流动，回到叶轮流到的内部。距离蜗壳隔舌远的那一部分气泡则会随着流动的液体和浮力流出离心泵。
（2）在半开式叶轮离心泵中，叶频和它对应的倍频是组成压力脉动主频的主体，在叶片出口、监测蜗壳的各个监测位置以及离心泵出口这些位置周围，对于压力脉动起到重要影响作用的是叶轮与蜗壳的动态与静态的相互干涉，入流气体含量占比越大，压力脉动幅值就越小。主要原因是当含气率增加到4.5%时，叶片表面几乎被气体附着，且叶轮流道及出口，蜗壳处出现大量气泡形成的气团，叶片对水的作用减弱，使得水流与蜗壳隔舌的干涉作用减小。

研究展望

本文对半开式叶轮离心泵气液两相流下进行的数值计算与实际流动情况仍然存在定偏差，由于计算对实际流态进行了一定简化，没有考虑气体的可压缩性、气体与液体之间复杂的作用力以及气体的真实形态变化等。对于半开式叶轮离心泵内的气液两相流研究，今后还需要以下拓展研究：
（1）在进行半开式叶轮水体划分时，将叶顶间隙层水体与叶轮水体分开划分，在对其用网格的手法进行划分的操作时，将这两个部分按照结构网格划分出来，可能会提高计算的准确性。
（2）在进行可视化实验时，应使用微距镜头，应用的特殊方法进行图像处理[70]，以提供气泡直径分布和轨迹的定量测量。
（3）由于环形流动条件会严重影响泵的性能，因此，基于CFD的数值模拟必须同时考虑气泡和环形流动条件，才能获得更加准确的结果。