混流泵兼具离心泵和轴流泵2种泵型的优点，可在较宽的流量区间内高效运行，广泛应用于水利、电力、水力推进等行业领域。混流泵在启动、停机、调转速、调阀等变工况条件下，内部会出现流动冲击、进出口回流、二次流、分离流、旋转失速等现象，将引起性能的剧烈变化，瞬态特性十分复杂。

近年来，许多学者相继开展了叶片泵启动过程的瞬态特性及其与内流场关联规律的相关研究。Tsukamoto等[1]对低比转速离心泵的启动过程进行研究，结果发现瞬态特性与准稳态假设的结果存在差异，泵无量纲扬程系数在启动早期远大于稳态假设值，之后迅速下降并低于稳态值。文[2-4]分别在不同启动加速度和不同管路系统中对离心泵进行启动试验，进一步验证了瞬态性能与稳态性能的差别。王乐勤等[5-6]发现启动过程中流量和扬程的增大滞后于转速的增加。在此基础上，胡征宇等[7]研究了不同阀门开度对离心泵启动特性的影响，结果发现阀门全关时实测扬程与理论扬程吻合较好，泵启动过程中一部分扬程被用作管路中流体加速。吴大转等[8-9]研究了不同启动加速度下离心泵的启动过程，结果发现较大的加速度能更好地抑制空化。

试验研究主要侧重泵启动过程的外特性，对泵启动过程的内流场研究主要是数值模拟。Li等[10-11]采用动态滑移区域法研究了离心泵启动中叶轮旋转加速引起的流场变形问题，结果表明叶片之间的主涡会在启动过程中占据流道造成扬程低于准稳态假设值。杨从新等[12]分析了离心泵启动过程叶片及蜗壳内的压力分布规律，发现启动初期受脉动作用的影响，泵内流场压力呈现进出口高、流道中间低等特点。陈红勋等[13]对ANSYS CFX软件进行了二次开发，根据输入驱动力矩自动更新转速曲线，所得计算结果与试验结果存在偏差，但变化趋势定性一致，验证了此种计算方法的有效性。文[14-15]通过数值计算分析了混流泵启动过程中压力场和速度场的变化过程，揭示了启动过程中大尺度涡堵塞流道是造成流量上升滞后的原因。文[16-17]对混流泵在不同进口压力和不同流量下的启动过程开展研究，采用高速摄影捕获了泵内部空化的组成及随时间变化的过程，结果发现转速加速到额定转速后叶顶泄漏涡空化首先出现，随后向流道内发展，进而形成叶道涡空化。

综上所述，关于泵启动过程的研究主要集中在外特性的试验测量和内流场的数值模拟，而对启动过程瞬态性能的理论建模和性能预测尚不深入。本文以混流泵为研究对象，建立了瞬态扬程与管路阻力耦合求解的方程，实现了混流泵瞬态特性预测，采用数值模拟验证理论预测结果的可靠性，为混流泵瞬态特性的优化控制奠定了基础。

以混流泵叶轮流道内控制体为研究对象，如图 1所示。图 1中，O为转轴上的一点，S为叶轮流道内控制体的封闭曲面，P为控制体内任意一点处的流体微团，r为点P到转轴的径向距离。dr为径向长度微元，dm=dr/sinγ为中间轴面流线的长度微元，γ为中间轴面流线与转轴的夹角。

基于不可压和无黏假设，对流道内控制体建立动量矩方程：

其中：T为动量矩，ρ为流体密度，t为时间。C为流体微团的绝对速度，V为体积，OP为流体微团的矢径，F、τ分别为作用在微团上的质量力和控制体封闭曲面S法线方向的主应力。

由Gauss定理，式(1)可写为

其中：Cn为闭合曲面S外法线方向上的速度，除叶轮进出口面外，其他面上的Cn均为0。故有

其中：T为T在转轴上的投影，Cm为轴面流速，Cu为周向流速，S1和S2分别代表叶轮的进出口面。

由式(3)可知，流道内流体的动量矩包含稳态项和瞬态项，其中稳态项可进一步表达为

瞬态项可进一步表达为

其中：ω为叶轮转速，β为叶片安放角，Iz为转动惯量。

因此，转轴上的动量矩大小可表示为

根据能量守恒定律，动量矩做功将转化为流体动能的增长量与静压能所做的功，即

其中：P为流体静压，C为绝对速度C的模。牵连速度U=ωr，相对速度W=Cm/sinβ。将式(6)代入式(7)，消去相同的项，可得到泵瞬态扬程Htransient，表示如下：

其中g为重力加速度。流体微元体积dV=2πrbψdm，其中: ψ为叶片排挤系数，b为过流断面线宽度。

将流体微元体积关系式代入式(8)，得到Htransient表示如下：

其中：Hsteady为稳态运行时的扬程，Q为流量。

式(9)表明，泵的瞬态扬程与稳态项扬程、由转速随时间变化引起的加速项扬程和由流量随时间变化引起的惯性项扬程有关。

启动加速过程中，管路系统的阻力由管路摩擦损失及管路流体加速的压差组成：

其中：φ为管路的损失系数，L、A和D分别为管路等效水力长度、等效水力面积和等效水力直径，λ为管路的阻力系数。当在稳态工况下运行时，管路系统阻力与泵扬程满足如下关系：

为得到流量和扬程随时间变化的规律，需对设计转速下的稳态扬程Hsteady, d进行多次函数拟合。本文选用6次多项式函数进行拟合，表示如下：

其中a0、a1、……、a6为多项式函数的系数。

启动过程中，转速在给定时间内加速至设计转速。为建立预测模型，假设泵在启动过程中满足相似定律，则有如下相似换算公式：

其中：n、nd分别为启动过程中的实际转速和设计转速。将式(13)代入式(12)，可得到基于相似换算的Hsteady，表示如下：

本文提出的启动过程扬程理论预测方法的计算流程如下：输入给定的转速加速变化曲线；由式(12)拟合得到流量-稳态扬程曲线；由式(9)计算瞬态扬程；由式(10)计算管路阻力；联立求解式(9)和(10)，得到dQ/dt；由dQ/dt得到Htransient与时间t的关系。

混流泵的设计参数为：流量Qd=0.03 m3/s，转速nd=1 000 r/min，扬程Hd=1.8 m。混流泵的主要几何参数为：叶轮叶片数为5，导叶叶片数为6，泵进口直径为150 mm，泵出口直径为180 mm。混流泵及其循环管路如图 2所示。

对混流泵的叶轮和导叶采用Turbogrid软件生成结构化网格，对管路采用ICEM软件生成六面体结构化网格，对弯管及阀门生成非结构化网格。网格无关性验证结果如表 1所示。

当网格总数达到3 216 432后，泵的扬程和效率变化不超过0.2%，满足网格无关性要求。因此，本文采用网格方案3开展研究。

采用k-ω SST湍流模型进行稳态和瞬态启动过程数值模拟。壁面设置为无滑移条件，对流项的空间离散采用高分辨率格式，瞬态项的时间离散采用二阶隐式。

泵稳态特性计算中，计算域包括进口延长段、混流泵、出口延长段。泵瞬态特性计算中，计算域包括混流泵、管路、阀门和水箱。进口边界给定为总压条件，出口边界给定为质量流量。叶轮进口与进口延长段，叶轮出口与导叶之间的动静交界面为Frozen-Rotor。启动工况下，叶轮转速由0上升至设计转速的加速时间为0.2 s，计算总时长为1 s，时间步长为0.000 2 s。动静交界面为Transient Rotor-Stator。

对混流泵不同流量下的稳态工况进行数值模拟，得到稳态扬程特性曲线，如图 3所示。混流泵在0.015~0.025 m3/s的小流量区，扬程存在马鞍区，这与比转速有关。同时采用6次函数对流量-扬程曲线进行拟合，如图中虚线所示，拟合后的残差平方和为0.034 4，说明拟合曲线与数值模拟结果吻合较好。

给定3组阀门开度A1、A2、A3，使泵在稳定运行工况下流量分别为0.8Qd、1.0Qd、1.5Qd。针对3组阀门开度，对启动过程混流泵的瞬态特性变化开展理论模型预测和数值模拟。

图 4为3种阀门开度下理论模型预测和数值模拟的流量变化曲线。结果表明：不同阀门开度下流量在加速初期快速增长，加速结束后缓慢增长并趋于稳定。小阀门开度A1下，数值模拟与理论模型预测的流量变化曲线几乎重合。随着阀门开度增大，两者在0.3 s后偏差增大，但在流量接近稳定阶段偏差减小。原因是理论预测模型中阻力系数为定值，而在实际流动中，阻力系数与流量和管道内的流动状态相关。

图 5为3种阀门开度下瞬态扬程式(9)中各项扬程的变化规律。结果表明，在转速加速阶段稳态项扬程的变化规律相同，在加速结束时达到最大值，随后迅速下降。在0.3 s后，随着阀门开度的增大，稳态项扬程开始逐渐下降，原因是阀门开度大，对应的流量大、扬程小。

3种开度下，加速项扬程的变化规律一致，在转速加速阶段存在恒定的加速项扬程，原因是转速加速为线性变化规律，加速度为恒定值。在0.2 s加速结束后加速项扬程变为0，说明该项仅仅与启动加速度和叶片参数有关。

3种开度下，惯性项扬程在转速加速阶段的变化规律相同，在加速结束时达到最大值。阀门开度大，管路系统阻力小，流量上升快，因此惯性项扬程更大。

将稳态项扬程、加速项扬程和惯性项扬程3者相加得到理论模型预测的瞬态扬程，并与数值模拟的结果进行对比，3种阀门开度下的结果如图 6所示。结果表明，理论模型预测结果与数值模拟结果基本一致，尤其在转速加速阶段吻合，验证了本文理论预测模型的准确性。t=0.2 s后，转速保持稳定，理论模型预测扬程与数值模拟扬程都逐渐下降；A1开度下理论预测扬程下降并保持稳定，与数值模拟扬程的偏差先逐渐增大直至基本不变；A2开度下理论预测值扬程先下降后缓慢增长，与数值模拟扬程逐渐接近；A3开度下理论预测扬程先下降后又上升，再急剧下降，在t=0.62 s附近与数值模拟扬程产生交点。

理论模型预测扬程与数值模拟扬程之间存在偏差的主要原因是：1) 各项扬程的结果表明，加速项扬程和惯性项扬程相比稳态项扬程小，稳态项扬程是影响预测结果准确性的重要因素。本文中稳态项扬程是由不同工况下扬程进行拟合得到的，与实际启动过程中扬程存在偏差。2) 启动过程主要对应小流量区间，在此区间泵内流态不稳定，使得数值模拟结果与实际启动过程中扬程存在偏差。3) 由于流量与扬程的相关性，理论模型预测流量与数值模拟流量的偏差会传递给扬程。

本文通过理论推导建立了混流泵启动过程瞬态性能理论预测模型，并通过数值模拟方法计算了混流泵和管路系统的三维瞬态流动，主要结论如下：

1) 基于混流泵瞬态扬程方程和管路系统阻力方程，建立了混流泵启动过程瞬态性能理论预测模型，通过给定转速加速曲线，得到启动过程中流量和扬程随时间的变化规律。理论模型预测结果与数值模拟结果吻合较好，验证了本文理论预测模型的准确性。

2) 转速加速过程中，混流泵流量随转速增大而快速增长。转速加速结束后，流量缓慢增长并最终趋于稳定。

3) 混流泵启动过程中，瞬态扬程可分为稳态项扬程、加速项扬程和惯性项扬程。转速加速过程中，加速项扬程影响较大。转速加速结束后，惯性项扬程影响较大。