1. 引言
水泵水轮机在运行过程中,安全稳定性至关重要。活动导叶起到控制流量的作用,同时,使水流在流经转轮前形成一定的环量。水泵水轮机运行过程中,活动导叶受到水力矩的作用。若水力矩使导叶有自开趋势,则不利于机组安全运行。另外,若某一导叶水力矩过大,则导叶保护装置生效,剪断销及导叶摩擦装置将起到保护作用,使得其它导叶仍可以正常动作,防止进一步破坏。因此,在设计阶段,了解导叶的水力矩特性对导叶保护装置等部件的设计十分重要。
目前,可以通过计算流体力学(CFD)的方法来预测导叶的水力矩特性。对非同步导叶情形下机组内部流场的研究 [1] [2] [3] [4] 表明,非同步导叶可降低处转轮外其他流域内压力脉动的峰峰值。在启动过程空载工况下会导致低频压力脉动的频率发生改变。飞逸工况下非同步导叶的应用破坏了水流的轴对称性,使得转轮受到非均匀的水力激励,从而引发振动。对导叶水力矩的实验研究 [5] 表明,导叶开度逐渐增大时,不同导叶的水力矩因数变化一致。启用非同步导叶后,非同步导叶与同步导叶开度相差较大时,不同导叶水力矩因数相差较大。Christophe等人 [6] 采用实验及数值模拟的方式研究了同步导叶及非同步导叶的水力矩,对于同步导叶,导叶开度较小时,导叶水力矩使得导叶有自关趋势。当导叶开度较大时,导叶水力矩使得导叶有自开趋势。而非同步导叶将造成其自身及相邻导叶的水力矩突增。Fan等人 [7] 对水轮机工况开停机过程中导叶水力矩进行了预测,发现停机过程中小开度区导叶水力矩变化较快,是导致振动和噪声的主要原因。但上述对于导叶水力矩的研究,均假定活动导叶出口为一固定静压值,未考虑到转轮对导叶的动静干涉作用。李琪飞等人 [8] 通过实验与数值模拟的方式,研究了水泵水轮机在水轮机工况下的导叶水力矩特性。研究针对模型水泵水轮机,并建立了水泵水轮机的全流道模型。研究发现导叶水力矩会被导叶附近的静压影响,导叶附近水压力的周期性变化导致了导叶水力矩随时间的变化。
目前,针对导叶水力矩的非定常特性研究较少,同一工况下不同导叶水力矩的变化规律尚未得到深入研究。本文采用数值模拟的方法,针对原型水泵水轮机稳定发电工况下的导叶水力矩特性进行了研究,分析了导叶水力矩的非定常特性。
2. 物理模型及数值模拟方法
本文计算采用的原型水泵水轮机如图1所示。转轮进水边直径D1 = 2.72 m,出水边直径D2 = 2.066 m,转轮叶片数Z = 7,固定导叶数Zsv = 20,活动导叶数Zgv = 20,导叶高度b0 = 1.112 m,活动导叶分布圆直径Dgv = 3.0875 m,转轮额定转速为200 rpm。将机组的导叶进行编号,如图2所示。
对上述模型进行离散,得到机组网格如图3所示。为确保计算结果的准确性,基于相似的网格拓扑结构,改变单元尺度,生成了3套网格,节点数分别为150万、350万和670万。选取额定工况进行网格无关性验证。不同节点数下机组外特性如图4所示。当网格数超过350万之后,机组的出力及效率趋于稳定。引入水力矩因数
其中,
为导叶水力矩因数;
为导叶水力矩,N∙m;
为水的密度,kg/m3;g为重力加速度,m/s2;H为水头,m;
为转轮进水边直径,m。
图5对比了不同节点数下导叶水力矩因数。可见节点数超过350万之后,导叶水力矩因数相差不大。
(a)
(b)
Figure 3. The grid of pump turbine. (a) The mesh of flow path, (b) The mesh of guide vane
图3. 导叶编号示意图。(a) 机组全流道网格。(b) 活动导叶域网格
![](//html.hanspub.org/file/15-2340543x19_hanspub.png)
Figure 4. External characteristics of the unit
图4. 不同节点数下机组外特性
![](//html.hanspub.org/file/15-2340543x20_hanspub.png)
Figure 5. The factor of hydraulic torque of guide vane
图5. 不同节点数下导叶水力矩因数
综合机组外特性及导叶水力矩计算结果,选取350万节点数的网格进行计算。
本文采用SST k-ω湍流模型进行计算。计算过程中在蜗壳进口给定总压边界条件,在尾水管出口给定静压边界条件。
3. 结果分析
图6显示了不同导叶开度下导叶水力矩因数的时均值。即使在同一导叶开度的情况下,每个导叶的水力矩大小并不相同。这主要是上游来流的不均匀性引起的。随着导叶开度的增加,导叶水力矩因数降低,且额定开度下,部分导叶水力矩的值由正变为负,形成自开趋势。
图7显示了额定开度(25˚)下机组相邻三个导叶的水力矩因数时域图及频域图。图中T为转轮旋转一周所需时间,fn为转频。从时域图中可见,每个导叶水力矩均呈周期性变化,且周期性明显。对导叶水力矩进行傅里叶变换,得到水力矩因数的频域图。可见水力矩因数的频率成分为转轮对导叶的动静干涉频率及其整数倍。这反映了动静干涉对导叶水力矩的影响。导叶水力矩的周期性变化主要是由于动静干涉引起的。另外,各导叶水力矩的主频及次频均为14倍转频及21倍转频,这表明影响导叶水力矩的主要压力脉动频率成分为14倍及21倍转频。
从时域图中可以发现,两个相邻导叶的水力矩之间存在时间差。假设#i导叶水力矩因数随时间变化的规律为
则#(i + 1)导叶水力矩因数随时间变化的规律可表示为
因此,对各导叶水力矩因数进行傅里叶变换,提取各导叶主频对应的相位
,求得相邻导叶之间的
相位差
。并将其换算成时间差
,其中f为主频。对于次频成分也做类似处理。结
果如图8所示。相邻导叶间水力矩时间差为转轮旋转周期的1/20,是转轮叶片扫过一个导叶所需时长。转轮顺时针旋转(从上往下看),因此转轮叶片按导叶编号顺序依次扫过各导叶,先扫过的导叶水力矩先发生变化,后扫过的导叶水力矩后发生变化。这进一步反映了转轮旋转对导叶水力矩的影响。
![](//html.hanspub.org/file/15-2340543x26_hanspub.png)
Figure 6. Mean value of the factor of hydraulic torque of guide vane
图6. 导叶水力矩因数时均值
![](//html.hanspub.org/file/15-2340543x27_hanspub.png)
Figure 7. The hydraulic torque factor of different guide vanes (left: time domain diagram; right: frequency domain diagram)
图7. 不同导叶水力矩因数(左:时域图,右:频域图)
图9为一个转轮旋转周期内#1导叶水力矩因数随时间的变化。从0.875T到0.935T这段时间内,导叶水力矩先上升后下降,直至变为负值。在这一时间段内,选取了8个时刻,如图9中8个点所示。定义相对压力
为
![](//html.hanspub.org/file/15-2340543x29_hanspub.png)
Figure 8. Time difference of hydraulic torque between adjacent guide vane. (a) Main frequency; (b) Sub frequency
图8. 相邻导叶间水力矩时间差。(a) 主频成分;(b)次频成分
![](//html.hanspub.org/file/15-2340543x30_hanspub.png)
Figure 9. Hydraulic torque factor of #1 guide vane
图9. #1导叶水力矩因数
其中
为尾水管出口压力,H为水头,
为水的密度,g为重力加速度。
图10为8个时刻下的导叶中截面相对压力分布。t = 0.875T时,#1导叶处于转轮两个叶片之间。随
![](//html.scirp.org/file/15-2340543x34_hanspub.png)
Figure 10. Relative pressure distribution in the cross section of the guide vane area
图10. 活动导叶区域中截面相对压力分布
着转轮叶片的旋转,#1导叶吸力面靠近头部的位置压力逐渐降低。受到转轮叶片头部压力面高压区的影响,#1导叶吸力面靠近尾部的位置压力逐渐升高。因此#1导叶水力矩上升。t = 0.91T时,导叶水力矩出现极大值。其后叶片继续扫过#1导叶,此时#1导叶主要受到转轮叶片背面低压区的影响,导叶吸力面靠近尾部的位置压力逐渐降低,导叶水力矩减小。t = 0.935T时,低压区压力极低,且接近#1导叶尾部,离导叶轴较远,贡献了较大的反向水力矩,因此在该时刻下水力矩为负值。
4. 结论
本文采用数值模拟方法研究了水泵水轮机在正常发电情况下的导叶水力矩特性。得到的主要结论如下:
随着导叶开度的增大,导叶水力矩的数值逐渐减小,并有变为负值的趋势。导叶开度较大时,机组运行过程中导叶水力矩在正值与负值之间变化,导叶轴上所受力矩不断反转,容易引起材料疲劳;
各导叶水力矩呈周期性变化,变化频率为转轮对导叶的动静干涉频率及其倍频。动静干涉现象是引起导叶水力矩变化的主要原因;
各个导叶水力矩变化趋势相同,但相邻导叶间水力矩变化存在时间差,为转轮旋转周期的1/20。这表明转轮叶片在旋转过程中依次扫过各导叶,使得导叶水力矩先后产生相似变化。
NOTES
*通讯作者。