DOI: 10.11817/j.issn.1672-7207.2018.10.015

混流泵启动过程瞬态流场的涡动力学分析

李伟¹，季磊磊¹，张扬¹，施卫东²，杨勇飞¹

(1. 江苏大学 国家水泵及系统工程技术研究中心，江苏 镇江，212013；

2. 南通大学 机械工程学院，江苏 南通，226019)

摘要：为了深入分析混流泵启动过程的瞬态流动结构，研究启动过程叶轮内部能量分布特性及其对瞬态性能的影响，基于正则化螺旋度法提取瞬态流场涡核，对启动过程进口段、叶轮和导叶段内部流动进行涡结构分析，并运用过流断面诊断法对混流泵启动过程内部流动进行诊断。研究结果表明：进口观测截面流场的涡核结构受叶轮叶片数的影响较大，涡核总体呈现从分散到集中的演化过程；随着转速增加，叶轮内涡结构正向和反向涡交替变化，并在转速稳定后流动逐渐趋于稳定；导叶内的涡结构在启动初期呈非对称性分布，当转速稳定后，涡结构区域逐渐减少并呈现规律性分布，流体流动趋于稳定。在混流泵启动过程中，随着叶轮旋转加速，总压流随之迅速增大，叶轮对流体做功，流体获得的能量迅速增加；由于流体惯性，加速末期流体获得了大于稳态转速下的能量，这种瞬态效应的外部体现就是在加速末期泵装置获得了瞬时冲击扬程。

关键词：混流泵；启动过程；瞬态流场；涡结构；涡动力学

中图分类号：TH313             文献标志码：A         文章编号：1672-7207(2018)10-2480-10

Vortex dynamics analysis of transient flow field at starting process of mixed-flow pump

LI Wei¹, JI Leilei¹, ZHANG Yang¹, SHI Weidong², YANG Yongfei¹

(1. National Research Center of Pumps and Pumping System, Jiangsu University, Zhenjiang 212013, China;

2. College of Mechanical Engineering, Nantong University, Nantong 226019, China)

Abstract: In order to further analyze the transient flow structure of the mixing pump at the start-up process and research the energy distribution in the impeller during the start up period as well as its effect on the transient characteristic of mixed-flow pumps, vortex core of transient flow field was extracted based on the regularized helical method. The vortex structure in the pump inlet, impeller and the guide vane was then analyzed and diagnosed through the information on the flow section method. The results show that the vortex structure on the section of the inlet is greatly affected by the blade number of impeller and the vortex concentrates from different directions. With the increase of the rotating speed of impeller, the positive and negative vortexes appear alternately and tend to be stable when the speed reaches the rated speed. The vortex structure in the guide vane is distributed asymmetrically; the vortex region shrinks and is regularly distributed when the speed becomes stable. During the start up period of mixed-flow pump, the total pressure increases quickly with the increase of the rotating speed. Due to the work of impeller, the energy of fluid increases rapidly. As a result of the inertial force of fluid, the energy of fluid is higher than that in stable working condition. The external appearance of such a transient effect is that there is a shock head at the end of the accelerating process.

Key words: mixed-flow pump; starting process; transient flow field; vortex structure; the vortex dynamics

混流泵在南水北调工程、国防军工等领域都有广泛应用^[1-3]。混流泵启动过程是一种特殊的瞬态加速过程，可为特殊的应用场合提供瞬时流体动力^[4]，但混流泵在启动过程中的瞬态水力效应容易引起冲击负载、水力激振和空化破坏等负面影响^[5-6]：因此，研究混流泵启动过程的瞬态水力特性，探索混流泵启动过程内部非定常流动结构，对深入地了解混流泵启动过程的瞬态水力效应，改进混流泵瞬态工况设计具有重要意义。传统的流场分析通常是通过流场的压力分布状况来判断流场的质量^[7-9]，但是利用这些常规流动参数只能从宏观上对流场的优劣进行判断，无法充分揭示水力设计方法对流场潜在的流体动力学影响。涡动力学研究方法是在20世纪80年代中期被提出并开始被大量采用^[10-12]，在研究涡量产生、发展及其与物体壁面相互作用以及在湍流发生、演化和流体控制中的作用等方面起到了主导性和基础性作用。涡动力学主要是通过分析任意垂直于流动方向的截面上总压流分部积分项的分布状况，找到导致总压增益不良的动力学根源及其发生位置，为设计提供直接的理论和诊断依据。WU等^[13-16]在该领域进行了大量研究，基于涡动力学对机翼表面的不佳流动结构进行了分析和预测。ZHU等^[17]基于涡动力学研究了机器鸟翅膀的旋转幅度和旋转时间对飞行的状态影响，发现适当地选取这2个参数可以控制机器鸟在水平方向和竖直方向上的速度。徐朝晖^[18]利用BVF方法诊断了高速离心泵叶片表面流场，在出现裂纹的相应位置上捕捉到BVF的剧烈变化。李凤超等^[19]引入边界涡量动力学理论对导叶和桨叶叶型进行全流道三维湍流模拟，分析了导叶和桨叶表面上边界涡量流分布状态，找到了局部流动不良的部位。樊红刚等^[20-21]采用涡动力学方法对可逆式转轮表面流场进行了BVF分析，为转轮的设计提供了诊断依据。张梁等^[22]将涡动力学理论用于混流式水轮机转轮内部流动研究，根据其内部涡核分布状况，预测了小流量工况下叶片吸力面上的流动分离。但是在混流泵研究方面，尤其是混流泵启动过程瞬态特性研究方面尚未见到相关基于涡动力学的研究报道。本文作者采用涡动力学方法对混流泵启动过程内部瞬态流场演化过程进行了非定常分析和诊断，对比分析试验和数值模拟下混流泵启动过程的外特性，基于正则化螺旋度法对进口流场、叶轮内部流场、导叶流场的涡结构演化过程进行诊断，分析沿流线截面上总压流的分布，以便为改进瞬态工况混流泵设计提供了直接的诊断依据。

1  研究模型

1.1  计算模型

本文研究的混流泵模型主要参数为：流量Q=380 m³/h，扬程H=6 m，转速n=1 450 r/min，比转速n_s=480。叶片数Z=4，导叶叶片数Z_d=7。通过Pro/E软件对混流泵模型的进口段、叶轮、导叶、蜗室和出口段分别进行建模，装配后得到混流泵全流道的三维实体造型，计算区域为从泵的进口段到环形蜗室出口段的整个装置段。混流泵模型如图1所示。

图1  混流泵模型

Fig. 1  Mixed-flow pump model

1.2  网格划分

在网格划分时，首先对混流泵各部分进行离散化，然后，分别对各部分进行网格划分，最后在ANSYS CFX中进行装配。考虑到叶顶间隙相对叶轮的尺寸很小，为保证后期计算的精准性，间隙内部必须提供足够的网格单元数和有效节点数以及轮缘间隙向叶轮内部网格均匀过渡，采用六面体网格进行划分。在叶轮和导叶处分别采用J/O型拓扑结构和H/O型拓扑结构，进口段采用Y-block拓扑形式，并进行周期阵列形成完整的过流通道。通过对每根拓扑线上节点的控制，使得网格按照各个壁面的曲率均匀变化，并对叶轮进行加密处理。

图2  全流道计算区域网格

Fig. 2  Computational area mesh of whole flow passage

1.3  控制方程和边界条件

在众多关于泵启动数值计算的研究中，已经证实standard k-ε模型能够较好地模拟其启动过程的内部瞬态流动特性^[23-24]，因此，本文也选用该模型进行数值研究，假设混流泵从启动开始，其内部流动即为湍流，以时均N-S方程作为基本控制方程，调用standard k-ε双方程湍流模型，采用二阶精度迎风格式，以基于微元中心有限体积法空间离散方式，实现压力速度的耦合求解。

数值计算中将泵划分为静止区域与旋转区域，对旋转区域做整体加速运动，旋转区域与静止区域之间通过滑移交界面进行连接，采用多重坐标系算法，在旋转区域加速旋转过程中，将旋转域与静止域之间的交界面选择Transient Frozen Rotor模式，交界面之间的数据通过差值方式进行传递。在数值模拟时，为保证数值计算的准确性，将加速时间设置与实验中的实际情况一致(混流泵启动加速过程实际时间为1.35 s)，采用匀加速启动。将试验流量随时间的变化曲线用三角函数进行拟合，并将此函数用CEL表达式写入CFX中作为出口流量的变化条件。将进口设置为opening，参考压力为1.01×10⁵ Pa。由于转轮室壁面附属于旋转区域，而转轮室壁面为绝对静止状态，因此，将其设置为 The Counter Wall。静止区域壁面设置为No Slip Wall。数值计算时主要通过调用CEL表达式来控制泵的加速运动过程及对变量进行监测，假设泵的加速过程为线性加速。计算介质为常温清水，密度1 t/m³，动力黏度1.0 mPa·s，并考虑重力影响。在进行非定常计算时，取总时间步长(Total time)为5 s，每个时间步长(Time step)取0.000 413 8 s，即每转内经历100 步。在每个时间步长内取最大迭代次数为2 000 次，以保证在每个时间步长内都绝对收敛。残差收敛精度设置为10^-4。

1.4  网格无关性验证

本文对混流泵在设计工况下稳定运行时的网格数进行了无关性检验。采用同样的网格拓扑结构，通过改变拓扑线条上的网格节点数目，并调整相应节点，使得网格质量保持一致。分别以多种不同数目的结构网格为载体进行计算，采用相同控制方程和边界条件以控制网格数量的单一变量。当网格数量达到110.01万左右时，再增加网格数量所得的计算扬程变化很小，误差在±5%以内，符合网格无关性检验要求。获得不同数目网格单元下的计算扬程如图3所示。

图3  计算网格数下的扬程对比

Fig. 3  Comparison of head with different mesh elements

2  数值方法的试验验证

2.1  试验装置和试验方案

为了验证混流泵启动过程瞬态数值模拟的准确性，搭建了如图4所示的实验装置系统来测量混流泵启动过程的瞬态外特性。

图4  试验装置图

Fig. 4  Experimental setup

由于电机启动过程中瞬时电流冲击过大，考虑到试验装置的安全性，采用15 kW的变频电机启动，输入电压为380 V，转速为1 500 r/min，频率为50 Hz。试验中采用自主开发的瞬态同步触发器测量瞬态转速，采用精度为0.2级的ZJ型转矩转速测量仪测量扭矩，采用精度为0.5级的LWGY型祸轮流量计测量流量，采用采样频率为17.4 kHz、测量精度为0.5%FS的MPM480型高频压力传感器测量动态扬程，由HSJ2010水力机械综合测试仪采集相关数据。

在启动实验开始前，先进行混流泵稳态工况实验，调节出口阀门，待流量稳定在设计工况点Q=380 m³/h之后，保持阀门开度不变，关闭混流泵机组。当实验管路内流体重新处于静止状态时，先启动性能参数测量仪，并设置变频器的启动时间为1 s，再启动混流泵机组。后续处理时，通过判断高频压力脉动信号中的突然波动为零时刻基准。由于压力传感器采用频率为17.4 kHz，由此导致的时间误差可以忽略不计。当转速稳定在1 450 r/min、流量稳定在380 m³/h时停机，待管路内流体重新稳定，进行上述3次重复性试验。

2.2  瞬态外特性对比

数值计算的瞬态外特性预测结果与实验瞬态外特性结果对比曲线如图5所示。

图5  瞬态预测结果与实验结果对比

Fig. 5  Comparison of experiment and transient calculation

由图5可知：虽然变频器设置的启动时间1 s，但在实际测量时，混流泵从转速为0 r/min加速到1 450 r/min的时间为1.35 s。在混流泵启动过程中，试验转速近似呈匀加速上升趋势，与数值模拟结果相一致。数值计算扬程与实验测量扬程保持了较好的一致性，随着转速到达最大值，扬程也立即到达最大值，并均出现1个瞬时冲击扬程，但流量随时间的变化滞后于扬程的变化。当启动过程结束时，转速逐渐趋于稳定，试验扬程和瞬态计算扬程均随着时间的增加呈下降趋势，在3 s后扬程逐渐趋于稳定值，下降幅度为1 m左右。在转速稳定阶段，数值计算结果高于实验测试扬程0.3 m左右。从预测扬程与实验扬程的相似程度来看，两者误差较小、变化趋势基本一致，说明混流泵启动过程的数值计算方法具有较高准确性。

3  混流泵启动过程涡运动结构

3.1  基于正则化螺旋度法涡核提取方法

考虑计算资源的限制，本文采用正则化螺旋度法分析混流泵内部流场的涡结构。采用螺旋度法可以准确地反映主涡与二次涡流动，描述旋涡形态以及变化规律。

正则化螺旋度法使用正则化螺旋度H_n提取涡核。正则化螺旋度H_n定义为速度与涡量的点积除以速度的模与涡量的模的乘积，如下式所示。

              (1)

式中：H_n为正则化螺旋度，[-1，1]；v为速度矢量；ω为涡量矢量。在涡核区域，速度矢量方向与涡量矢量方向近于平行，正则化螺旋度H_n趋于±1。正则化螺旋度H_n的符号表明涡旋转的方向，以流动方向为正方向，若H_n为正，则涡旋转方向为逆时针方向；若H_n为负，则涡旋转方向为顺时针方向。

3.2  涡结构辨识结果

在混流泵进口，由于叶轮加速旋转的影响，压力分布起伏较大，容易形成旋涡。利用正则化螺旋度法研究进口段的诱导旋涡形态以及变化规律。在距离叶轮进口10 mm设立了观测面，均垂直于旋转轴。图6所示为混流泵启动过程(从叶轮出口向叶轮进口方向观察，下同)10 mm观测平面的正则化螺旋度H_n云图。

从图6可以看出：受到叶轮叶片数的影响，混流泵启动过程观测平面内在不同时刻涡核结构始终有4个相似区域沿周向均匀分布。在启动初期，观测平面上涡核分布较为分散，旋涡旋转方向与叶轮旋转方向一致。在0.48 s时，涡核结构迅速增加并向圆心方向集中。随着转速继续增大，涡核结构由圆心向外扩张，分布于截面圆周上。当叶轮转速逐步稳定时，涡核区域又向圆心方向集中，截面圆心并没有出现旋涡结构。

获得叶轮YZ截面(如图7所示)在启动过程不同时刻的正则化螺旋度H_n云图，如图8所示。

从图8可以看出：在混流泵启动过程初始阶段的0.18 s时，叶轮进口截面流道内出现正向涡结构，同时，在靠近叶轮叶片区域出现反向涡结构；随着转速的增加，在0.48 s时，叶轮流道内的正向涡结构强度减弱并且区域减小，但在靠近叶片吸力面附近，涡结构强度增加并且区域增大。在0.78 s时，YZ截面已经位于叶轮叶片尾部，流道内靠近轮毂区域出现大块反向涡集中区，同时，在靠近叶轮端壁的流道内，正向涡结构强度增加。当混流泵转速继续增大并达到最高时，叶轮流道内基本被正向涡所占据，但强度不高。

图6  启动过程进口段内距离叶轮10 mm处H_n云图

Fig. 6  H_n diagrams of inlet section in start-up process at 10 mm

图7  叶轮YZ截面

Fig. 7  YZ section of impeller

为了分析导叶内部涡核结构形态以及变化规律，分别在导叶进口、中部和出口分别截取3个截面，截面之间的距离为60 mm。图9所示为进口截面上的正则化螺旋度H_n云图。从图9可以看出：导叶进口截面流场涡结构呈现明显的非周期性，完全区别于进口段和叶轮内的涡结构，这可能是叶轮和导叶动静干涉或者启动初始阶段流体突然获得能量所导致。混流泵启动过程初始阶段，在0.18 s时，正向涡结构区域在导叶流道内间隔出现，并占据大部分流道。随着转速的增加，在0.48 s时，流道内涡结构减少并伴随反向涡的产生。随着转速进一步增加，2种涡结构区域均减少，正向涡区域收缩至导叶工作面和轮毂附近。当转速继续增加并趋于稳定时，截面内涡强度缓慢减弱，沿圆周方向在靠近端壁附近产生了非对称性的、不连续的正向涡和反向涡区域。

图10所示为启动过程导叶中间截面的正则化螺旋度H_n云图。从图10可以看出：在混流泵启动过程初始阶段，导叶中间截面内涡结构的非周期性依然明显。在0.48 s时，导叶流道内被大量正向涡占据并且呈无规则分布；随着转速的增加，这些正向涡结构集中于导叶轮毂附近，并且在导叶流道内呈现一定的周期性。随着转速继续增加并趋于稳定，正向涡结构强度逐渐减弱并且区域缓慢缩小，在靠近端壁区，出现了周期性的反向涡结构。

图11所示为启动过程导叶出口截面的正则化螺旋度H_n云图。从图11可以看出：在混流泵启动初始阶段，由于湍流尚未充分发展，导叶出口截面的涡结构也是非周期性的，在导叶3个流道内正向涡和反向涡交替出现，呈现无规则性。在0.48 s时，各流道中部均有正向涡存在，靠近轮毂处均出现反向涡结构。随着转速进一步增加，在到0.78 s时，流道内正向涡区域增加，反向涡结构减少，正向涡区域逐步呈现一定的周期性分布；随着转速趋于稳定，导叶内正向涡结构区域形态趋于一致，并在流道靠近端壁处产生了小范围的反向涡结构。

图8  启动过程叶轮YZ截面H_n云图

Fig. 8  H_n diagrams of YZ section in start-up process

图9  启动过程导叶进口截面H_n云图

Fig. 9  H_n diagrams of guide vane inlet section in start-up process

图10  启动过程导叶中部截面H_n云图

Fig. 10  H_n diagrams of guide vane middle section in start-up process

图11  启动过程导叶出口截面H_n云图

Fig. 11  H_n diagrams of guide vane outlet section in start-up process

4  混流泵启动过程叶轮内部涡结构诊断

4.1  过流断面诊断方法

流体动量方程的积分形式为

         (2)

式中：V为控制容积；D/Dt为随体导数；u为速度矢量；f为体积力；S为控制体的面积；τ为空间变量、时间变量t和面元方向n的函数，τ=τ(x，y，n)。这说明存在二秩张量T(x,t)，使得τ=τ(x，y，n)=n·T(x，y，n)，则式(2)可化为

             (3)

对式(3)点乘速度矢量u，考虑T的对称性，根据雷诺输运定理，在大雷诺数Re下惯性力远大于黏性力，有

          (4)

式中：ΩM_z为叶轮施加给流体的轴功率；K为总动能之和，P和D分别为整个控制容积所作的压缩功和耗散功，分别表示如下：

            (5)

              (6)

式中：Φ为熵增引起的耗散率。式(4)中的G表示流体经过流道后能量的增加过程，可以表示为

          (7)

其中：W为流道的过流断面；u_l为沿流线方向的速度；U为轴向速度；S_in为流道的进口断面；p^*和p^*_∞分别表示为：

            (8)

           (9)

式中：为静压力，并令。

利用式(4)计算叶轮施加给流体的轴功率时，对于启动过程的非定常流动，K/t由启动过程的加速度、叶轮和流道几何形状共同决定。流体在流道进口的涡量和张量很小，在流道内计算得到的P和D较P_u小很多数量级，因此，利用式(4)计算叶轮叶片施加给流体的轴功率时，P_u是主要参数，称P_u为总压流，从客观上反映了流道中流体能量变化过程。

4.2  叶轮内过流断面涡诊断

在靠近叶轮水体段进口、中部和出口，沿着流体流动方向取3个与叶轮流道近乎垂直的过流断面(如图12所示)，并分别对各个断面做总压流积分，获得混流泵启动过程不同时刻总压流积分的变化曲线，如图13所示。

从图13可以看出：在混流泵启动过程中，随着叶轮旋转加速，P_u随之迅速增大，由于叶轮对流体不断做功，流体获得的能量迅速增加；在t=1.38 s时，叶轮加速基本完成，P_u出现了最大值随后其值略有减小，并逐渐趋于一个稳定值。上述现象可以理解为：在叶轮加速到额定转速的过程中，由于流体惯性，其在加速末期获得的流体能量大于稳态转速下流体获得能量。这也是混流泵启动过程外特性研究中加速末期泵装置获得瞬时冲击扬程的本质原因。从图13还可以看出：从叶轮进口到叶轮出口，在加速初期各个截面P_u变化不大，当t=1.08 s之后，P_u曲线呈现1个明显的增大过程，瞬态效应凸显。

图12  叶轮3个截面

Fig. 12  Three sections of impeller

图13  启动过程叶轮各个截面总压流积分P_u变化曲线

Fig. 13  P_u curves of each section in start-up process

图14  启动过程叶轮内3个截面H_n云图

Fig. 14  H_n diagrams of impeller three sections in start-up process

图14所示为3个截面的正则化螺旋度H_n云图。从图14可以看出：在混流泵启动初期，叶轮流道内反向涡结构随着时间的增加明显减少，随着加速结束，流场趋于稳定，由于流体惯性，在1.68 s时略有回升。这和上述总压流积分曲线相对应，说明由于流体惯性，在叶轮转速上升过程中，叶轮流道内正向涡结构(反向涡结构)并不是一直增加(减少)的，而是与时间尺度相关，在某个时刻存在1个最大值(最小值)。随着叶轮内流场结构趋于稳定，叶轮流道内正向涡结构占据整个叶轮流道，只在叶轮出口附近有较小区域的反向涡结构存在。

5  结论

1) 混流泵启动过程中，叶轮流道内正向涡和反向涡交替出现，当转速稳定后，正向涡占据主导作用，叶轮做功和扰动效应明显。但由于流体惯性，随着叶轮旋转加速，叶轮流道内正向涡结构(反向涡结构)并不是一直增加(减少)的，而是和时间尺度相关，在某个时刻存在1个最大值(最小值)。

2) 在加速初期，叶轮各截面总压流变化不大，但随着转速不断升高，在1.08 s以后，总压流沿叶轮流道出现急速增加后逐渐降低并趋于稳定的过程，在外部能量特性表现为加速末期动态扬程出现峰值，随后又回落到稳态工况的水平，瞬态效应凸显。

3) 本文利用涡动力学理论，诊断了混流泵启动过程三维瞬态流场的涡核结构，分析了混流泵叶轮内总压流变化的原因及其对启动性能的影响，为研究瞬态工作水泵的流体能量分布特性和叶轮做功情况提供了一种直观有效的分析方法。

参考文献：

[1] BING H, CAO S, TAN L, et al. Effects of meridional flow passage shape on hydraulic performance of mixed-flow pump impellers[J]. Chinese Journal of Mechanical Engineering, 2013, 26(3): 469-475.

[2] KIM S, LEE K Y, KIM J H, et al. High performance hydraulic design techniques of mixed-flow pump impeller and diffuser[J]. Journal of Mechanical Science and Technology, 2015, 29(1): 227-240.

[3] LI W, ZHOU L, SHI W D, et al. PIV experiment of the unsteady flow field in mixed-flow pump under part loading condition[J]. Experimental Thermal and Fluid Science, 2017, 83(4): 191-199.

[4] DAZIN A, CAIGNAERT G, BOIS G. Transient behavior of turbomachineries: applications to radial flow pump startups[J]. Journal of Fluids Engineering, 2007, 129(11): 1436-1444.

[5] 吴大转, 焦磊, 王乐勤. 离心泵启动过程瞬态空化特性的试验研究[J]. 工程热物理学报, 2008, 29(10): 1682-1684.

WU Dazhuan, JIAO Lei, WANG Leqin. Experimental study on cavitation performance of a centrifugal pump during starting period[J]. Journal of Engineering Thermophysics, 2008, 29(10): 1682-1684.

[6] 季磊磊, 李伟, 施卫东, 等. 导叶式混流泵内部非定常流动特性数值模拟[J]. 农业机械学报, 2016(S1): 155-162.

JI Leilei, LI Wei, SHI Weidong, et al. Numerical simulation of unsteady flow characteristics in mixed-flow pump with guide vanes[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016(S1): 155-162.

[7] 季磊磊, 李伟, 施卫东, 等. 叶片数对混流泵内部非定常压力脉动特性的影响[J]. 排灌机械工程学报, 2017, 35(8): 666-673.

JI Leilei, LI Wei, SHI Weidong, et al. Influence of different blade numbers on unsteady pressure pulsations of internal flow field in mixed-flow pump[J]. Journal of Drainage and Irrigation Machinery Engineering, 2017, 35(8): 666-673.

[8] ZHOU L, SHI W D, CAO W D, et al. CFD investigation and PIV validation of flow field in a compact return diffuser under strong part-load conditions[J]. Science China Technological Sciences, 2015, 58(3): 405-414.

[9] LI W, JIANG X, PANG Q, et al. Numerical simulation and performance analysis of a four-stage centrifugal pump[J]. Advances in Mechanical Engineering, 2016, 8(10): 1-8.

[10] SHAO C, ZHOU J, GU B, et al. Experimental investigation of the full flow field in a molten salt pump by particle image velocimetry[J]. Journal of Fluids Engineering, 2015, 137(10): 104501.

[11] WU Jiezhi, WU Jianming. Interactions between a solid surface and a viscous compressible flow field[J]. Journal of Fluid Mechanics, 1993, 254: 183-211.

[12] WU X H, WU J Z, WU J M. Effective vorticity-velocity formulations for three-dimensional incompressible viscous flows[J]. Journal of Computation Physics, 1995, 112: 68-82.

[13] WU J Z, WU X H, WU J M. Streaming vorticity flux from oscillating walls with finite amplitude[J]. Physics Fluids A, 1993, 5(8): 1933-1938.

[14] WU J Z, MA H Y, ZHOU M D. Vorticity and vortex dynamics[M]. Berlin: Springer-Verlag, 2006: 1-782.

[15] WU J Z, TRAMEL R W. A vorticity dynamics theory of three-dimensional flow separation[J]. Physics of Fluids, 2002, 12(8): 1932-1954.

[16] WU J C. Elements of vorticity aerodynamics[M]. Shanghai: Shanghai Jiaotong University Press, 2018: 1-147.

[17] ZHU L, GUAN H, WU C. Vorticity dynamics and control of self-propelled flying of a three-dimensional bird[J]. Acta Aerodynamica Sinica, 2016, 34(2): 204-213.

[18] 徐朝辉. 高速离心泵内全流道三维流动及其流体诱发压力脉动研究[D]. 北京: 清华大学能源与动力工程系, 2004: 1-163.

XU zhaohui. The analysis of three-dimensional flow in high-speed pump and research of its fluid-induced pressure fluctuation[D]. Beijing: Tsinghua University. Department of Energy and Power Engineering, 2004: 1-163.

[19] 李凤超, 樊红刚, 王正伟, 等. 贯流式水轮机桨叶涡动力学优化设计[J]. 清华大学学报(自然科学版), 2011, 51(6): 836-839.

LI Fengchao, FAN Honggang, WANG Zhengwei, et al. Optimum design of runner blades of a tubular turbine based on vorticity dynamics[J]. Journal of Tsinghua University(Science and Technology), 2011, 51(6): 836-839.

[20] 杨琳, 樊红刚, 陈乃祥. 基于涡动力学的可逆转轮双向流动诊断及优化设计[J]. 清华大学学报(自然科学版), 2007, 47(5): 686-690.

YANG Lin, FAN Honggan, CHEN Naixiang. Bidirectional flow diagnosis to optimize the design of a pump-turbine runner using vorticity dynamics theory[J]. Journal of Tsinghua University(Science and Technology), 2007, 47(5): 686-690.

[21] 樊红刚, 陈乃祥, 杨琳. 可逆转轮三维流动的涡动力学诊断研究[J]. 水力发电学报, 2007, 26(3): 124-128.

FAN Honggang, CHEN Naixiang, YANG Lin. Three dimensional flow diagnosis of the pump turbine runner based on the dynamic vorticity[J]. Journal of Hydroelectric Engineering, 2007, 26(3): 124-128.

[22] 张梁, 刘树红, 张乐福, 等. 混流式转轮内部流场的涡动力学分析[J]. 水力发电学报, 2007, 26(6): 106-110.

ZHANG Liang, LIU Shuhong, ZHANG Lefu. Vorticity dynamics analysis of flow field in Francis runner[J]. Journal of Hydroelectric Engineering, 2007, 26(6): 106-110.

[23] 李伟. 斜流泵启动过程瞬态非定常内流特性及实验研究[D]. 镇江: 江苏大学流体机械工程技术研究中心, 2012: 1-159.

LI Wei. Experimental study and numerical simulation on transient characteristics of mixed-flow pump during starting period[D]. Zhenjiang: Jiangsu University. Research Center of Fluid Machinery Engineering and Technology, 2012: 1-159.

[24] 李伟, 季磊磊, 施卫东, 等. 混流泵非均匀轮缘间隙流场数值计算[J]. 农业机械学报, 2016, 47(10): 66-72.

LI Wei, JI Leilei, SHI Weidong, et al. Numerical calculation of internal flow field in mixed-flow pump with non-uniform tip clearance[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(10): 66-72.

(编辑  杨幼平)

收稿日期：2017-10-10；修回日期：2017-12-04

基金项目(Foundation item)：国家自然科学基金资助项目(51679111，51579118，51409127)；国家重点研发计划项目(2017YFC0403703)；江苏省重点研发计划项目(BE2015119，BE2015001-4，BE2016319，BE2017126)；江苏省自然科学基金资助项目(BK20161472)；江苏高校优势学科建设工程资助项目(PAPD)(Projects(51679111, 51579118, 51409127) supported by the National Natural Science Foundation of China; Project (2017YFC0403703) supported by the National Key Research and Development Program of China; Projects(BE2015119, BE2015001-4, BE2016319, BE2017126) supported by the Key Research and Development Program of Jiangsu Province; Project(BK20161472) supported by the Natural Science Foundation of Jiangsu Province; Project(PAPD) supported by the Priority Academic Program Development of Jiangsu Higher Education Institutions)

通信作者：李伟，研究员，博士生导师，从事流体机械及工程研究；E-mail：lwjiangda@ujs.edu.cn