导叶翼型优化对水泵水轮机“S”特性的影响.pdf

第 44 卷第 8 期2023 年 8 月哈 尔 滨 工 程 大 学 学 报Journal of Harbin Engineering UniversityVol.44.8Aug.2023导叶翼型优化对水泵水轮机“S”特性的影响李琪飞1,2,黄腾1,郑非1,辛路1,谢耕达1(1.兰州理工大学 能源与动力工程学院,甘肃 兰州 730050;2.甘肃省流体机械及系统重点实验室,甘肃 兰州 730050)摘 要:为了改善水泵水轮机“S”特性,本文以某水泵水轮机模型机为研究对象,采用 SST k-湍流模型开展水泵水轮机全流道三维数值计算。通过优化活动导叶翼型,得出 n11-Q11模拟曲线和机组“S”特性改善情况,并对优化活动导叶翼型后的机组的无叶区进行了压力脉动分析。结果表明:优化活动导叶翼型对机组的“S”特性存在一定改善效果;翼型优化后机组的效率虽然有小幅度下降,但是整体效率依然保持在 90%以上;无叶区压力脉动主要是由活动导叶区域出流,转轮叶片工作面发生冲击导致的分流以及转轮叶片与活动导叶组成的动静叶栅干涉效应引起的。改型后的活动导叶降低了无叶区的压力脉动幅值,在一定程度上提升了水泵水轮机的安全系数。关键词:水泵水轮机;“S”特性曲线;活动导叶;翼型设计;数值计算;动静叶栅干涉;安全系数DOI:10.11990/jheu.202106034网络出版地址:https:/ 文献标志码:A 文章编号:1006-7043(2023)08-1354-07Effect of guide wane airfoil optimization on the S characteristics of pump turbines LI Qifei1,2,HUANG Teng1,ZHENG Fei1,XIN Lu1,XIE Gengda1(1.School of Energy and Power Engineering,Lanzhou University of Technology,Lanzhou 730050,China;2.State Key Laboratory of Fluid Machinery and Systems of Gansu Province,Lanzhou 730050,China)Abstract:In this study,a pump turbine model was taken as the research object,and the SST k-turbulence model was used to conduct the three-dimensional numerical calculation of the entire channel of a pump turbine to improve the S characteristics of pump turbines.The n11-Q11 simulation curve was obtained,and the S characteristics of the unit were improved by optimizing the airfoil of the movable guide vane,and pressure pulsation analysis was per-formed on the vaneless zone of the unit after the airfoil of the movable guide vane was optimized.Results showed that the optimized movable guide vane airfoil improved the S characteristics of the unit.Although the efficiency of the unit was slightly reduced after airfoil optimization,the overall efficiency of the unit remained above 90%.Pres-sure pulsation in the vaneless zone was mainly caused by the divergence of the flow from the movable guide vane ar-ea and the working surface of the runner blade due to impact as well as interference between the runner vane and the active guide vane cascades composed of dynamic and static vane cascades.The modified movable guide vane reduced the amplitude of pressure pulsation in the vaneless zone.This effect improved the safety factor of the pump turbine to a certain extent.Keywords:water pump turbine;S characteristic curve;movable guide vane;airfoil design;numerical calculation;interference between the runner vane and the active guide vane cascades;safety factor收稿日期:2021-06-10.网络出版日期:2023-05-05.基金项目:国家自然科学基金项目(52066011);西华大学教育部重点实验室开放课题(szjj2019-025).作者简介:李琪飞,男,副教授.通信作者:李琪飞,E-mail:lqfy .水泵水轮机在“S”特性区会出现转速对应多个流量的非设计工况,因而机组在“S”特性区内的运行非常不稳定1-3。机组启动阶段,多种工况之间的瞬态变化可能导致机组并网时间变长甚至并网失败;机组甩负荷后,由于引水管道内依然存在大量的水,会引起转轮高速旋转,并使机组进入反水泵工况4-6;当机组在“S”特性显著的工况下运行时会导致机组无叶区强烈的压力脉动,引起机组发生强烈震动7-9。水泵水轮机在“S”特性区运转极其不稳定,转速值对应多个流量值有可能使得机组在多个工况间频繁转变,这种情况会导致机组运转严重失稳10-13;预开导叶可以改善水泵水轮机进入不稳定的“S”区,但是有可能增大机组无叶区压力脉动值,从而影响水泵水轮机机组安全稳定运行14-16。第 8 期李琪飞,等:导叶翼型优化对水泵水轮机“S”特性的影响本文采用 SST k-湍流模型,以国内某型抽水蓄能电站水泵水轮机为研究对象,通过优化活动导叶翼型改善水泵水轮机“S”特性以及无叶区压力脉动。1 设计流程及模型建立1.1 活动导叶设计流程 由于水泵水轮机的导水机构为圆柱形导水机构,因此可以通过改变活动导叶的翼型型线对现有的活动导叶改型。本文参考水轮机优秀水力设计模型、水轮机设计手册导叶叶型断面尺寸表以及原有翼型设计数据对活动导叶翼型进行改型,在保持原有导叶弦长不变的情况下设计出新的活动导叶。通过改变活动导叶翼型改变入水方向与速度从而改善水泵水轮机机组“S”特性。改型前后导叶翼型如图1 所示。注:D0为 638.619 mm,a 为 11.036 mm,a1为 2.85 mm,b 为 11.247 mm,b1为 5.011 mm,c 为 11.145 mm,c1为6.903 mm,d 为 10.273 mm,d1为 7.714 mm,L1为 47.385 mm,e 为 9.272 mm,e1为 8.238 mm,d0为 16.099 mm,m 为7.468 mm,m1为 5.891 mm,k 为 0.912 mm,f 为 4.722 mm,r 为 6.016 mm,L 为 89.897 mm,L2为 42.611 mm图 1 改型前后导叶翼型示意Fig.1 Schematic diagram of the modified anterior and posterior guide vane airfoils从理论分析,通过优化前后速度三角形,U=nD/60,V1 m=Q/A,A 指叶片进口边过流面积,同时通过 CFDpost测出叶片进口安放角的大小,通过叶片进口安放角,绝对流速 V 以及轴面流速 V1 m得出速度三角形。如图 2 所示在假定圆周速度方向和大小不变化条件下,如果转轮的进口水流角(水流相对速度与圆周速度之间的夹角)和叶片进口安放角 相等时,无脱流现象,其水力损失最小。因此,水流冲角 越大,水力损失越大。活动导叶转动角度不变,活动导叶分布圆直径增大,绝对速度受到牵连速度的影响发生变化,绝对速度减小,来流液体与转轮叶片冲角 减小,有效过流流量增大,转轮水力损失相应减小。图 2 改型前后导叶翼型速度三角形对比及叶片进口速度三角形Fig.2 Comparison of the front and rear guide vane airfoil velocity triangles of the modification1.2 网格划分及模型建立 本文研究对象为水泵水轮机模型,过流部件由蜗壳、固定导叶、活动导叶、转轮和尾水管组成,示意如图 3 所示,具体参数如表 1 所示。为了保证数值计算结果可行可靠,本文网格划分采用商业软件 ANSYS 的子功能 ICEM 进行全流道六面体网格划分。经过网格无关性验证,最终采用网格质量为 0.3 左右,网格总数为 610 万左右的网格进行计算。网格划分结果如图 4 所示。图 3 模型水泵水轮机计算区域Fig.3 Model pump turbine calculation area表 1 模型水泵水轮机几何参数Table 1 Geometric parameters of model pump turbine参数名称数值叶片数/个9活动导叶/个20转轮高压侧直径/mm473.6蜗壳进口直径/mm315导叶高度 b0/mm66.72固定导叶数/个20导叶高度/mm66.72转轮低压侧直径/mm300尾水管出口直径/mm6602 试验台简介及数值计算方法 为了验证湍流模型的可靠性以及建模与数值计算的可行性,本文拟将数值计算结果与实验结果进5531哈 尔 滨 工 程 大 学 学 报第 44 卷行对比。试验台示意如图 5 所示。本文模型试验采用 H=30 m 恒定水头对模型机进行试验,利用转矩仪以及流量计测出转轮转速以及流量并与计算结果进行对比。图 4 局部网格示意Fig.4 Local grid diagram注:1.偏流器,2.喷嘴,3.低压箱,4.测功电机,5.扭矩测量系统,6.水泵水轮机 7.高压箱,8.支架,9.流量计,10.封闭系统回路管,11.水泵 12.敞开系统回水管,13.称重传感器,14.称筒 15.水冷系统,16.换向管路图 5 水泵水轮机试验平台Fig.5 Water pump turbine test platform由于 SST k-湍流模型可以有效捕捉近壁面的流动,尤其是对水泵水轮机这种多导叶、多叶片的复杂几何模型有较强的准确性以及适应性,因此本文选取 SST k-湍流模型开展数值模拟。流体介质设为常温水,壁面采用无滑移壁面边界条件;进口和出口分别设置为质量流量进口和自由出流;静止域与旋转域之间的数据传递依靠 INTERFACE 边界条件。采用 SIMPLE C 速度压力耦合算法,残差值设置为 10-6,时间步长设置为 0.000 467 3 s。3 计算结果及分析3.1 可靠性验证 本文研究选取模型水泵水轮机活动导叶 a0=33 mm 进行数值计算可靠性验证。选取7 个工况点对其进行定常数值计算。将数值计算的结果进行单位转速 n11、单位流量 Q11的转换分别为:n11=nD2/H(1)Q11=Q/(D22H)(2)通过换算所得 Q11和 n11,绘制出 n11-Q11特性曲线。利用换算所得 n11-Q11特性曲线与试验曲线进行对比,其结果如图 6 所示。通过对比,两者具有很高的吻合度,误差值保持在 6%以内,满足工程研究要求。因此,本文数值计算所选用的模型具有较高的可靠性。图 6 实验与模拟结果对比Fig.6 Comparison of experimental and simulation results3.2 特性曲线分析 根据前文的换算结果,对原有导水机构导叶替换为优化后的翼型导叶模型水泵水轮机进行数值模拟,分别得到数据。并绘制出 n11-Q11曲线,如图 7所示。图 7 翼型优化对“S”特性影响Fig.7 Influence of airfoil optimization on“S”character-istics从图 7 的试验曲线可以看出,在导水机构导叶为原有翼型的情况下出现了很明显的“S”特性区。由于“S”特性的存在导致机组在运行的过程中及其的不稳定。将原有导水机构中的导叶替换为优化翼型的导叶后,机组“S”特性得到了一定的改善。因此,通过优化导叶翼型来改善水泵水轮机“S”特性可行。3.3 S1流面流态分析 为了更好展现出水泵水轮机流场特性,本文选取转轮部位 S1流面作为进一步分析的研究对象。将转轮上冠下环的等距面选取作为 S1流面。利用后处理软件 CFD-Post 做出 S1流面的速度流线分布图,如图 8 所示。6531第 8 期李琪飞,等:导叶翼型优化对水泵水轮机“S”特性的影响图 8 不同工况下 S1流面流线图Fig.8 Streamline diagram of S1 stream surface under different working conditions 由图 8 可以看出,147.9 kg/s 流量工况下,不管活动导叶翼型优化与否,在无叶区均出现了较明显的高速水环,但是将活动导叶翼型优化后,无叶区的高速水环相较优化前有了一定的削弱,从转轮区域的叶道流场分析,147.9 kg/s 流量工况下,装配原有活动导水机构机组的叶道内分布了大量的涡,且基本集中在距转轮入口处,然而优化活动导叶翼型后,转轮区域流态较优化前有了较好的改善,部分叶道内未出现明显的旋涡流动。从 325.3 kg/s 以及414.078 kg/s 的 2 个工况点的流面看,优化活动导叶翼型前后流态相似,并且在活动导叶及转轮流域内均未出现明显的旋涡。根据以上分析,优化活动导叶翼型后,由于翼型骨线相较优化前发生了变化,活动导叶区域出流角也随之改变,从而使得无叶区高速水环在一定程度上被破坏,进而改善了转轮入流量,这也是导致机组“S”特性得到一定改善的其中一个原因。3.4 机组效率计算 为了研究采用优化翼型活动导叶的机组的效率,以导叶开度 a0=33 mm 工况下模型水泵水轮机为研究对象,以选用优化活动导叶翼型后的水泵水轮机为模型,在流量特性曲线上选取 Q1=750 kg/s,Q2=650 kg/s,Q3=530 kg/s,Q4=300 kg/s 共 4 个工况点进行数值计算,计算出机组效率,如图 9 所示。由图 9 可以看出,将活动导叶翼型优化后,机组效率出现了微小幅度的下降,但针对可逆式机组小流量工况的过渡过程,其效率的微小幅度下降是在可允许范围内的。而且,在机组效率为小幅度下降的前提下,对机组“S”特性的改善应该占据主导地位,以保证机组在小流量工况下的稳定运行。显然优化活动导叶翼型对改善机组“S”特性取得了较好的效果。因而本文活动导叶翼型的优化是有效的。图 9 机组各工况点效率Fig.9 The efficiency of the unit operating point3.5 S1流面湍动能分析 由于将原有活动导叶全部替换为优化翼型导叶之后,机组整体效率出现了微小幅度下降。而转轮又是机组主要的能量转化过流部件,因此分析转轮部位 S1流面湍动能分布是有必要的。由图10 可以看出,湍动能区域出现在无叶区附近影响了叶片进口来流,最大湍动能值为32.44 kg/J。改型后最大湍动能值增大为 44.24 kg/J,均集中在无叶区,但无叶区附近的湍动能区域面积有所减小但7531哈 尔 滨 工 程 大 学 学 报第 44 卷其强度变高,整体出现了一些能量耗散较大的区域,比较 Q4工况,改型前后湍动能分布区域较为相似,均出现无叶区附近,只是改型后环形无叶区的湍动能面积相应减小,活动导叶出口区域能量损失降低,从整体来看,将活动导叶翼型优化后,无叶区以及转轮区域均出现了一定的能量耗散,导致机组整体效率出现了微小幅度的下降。而造成能量损失主要的原因有:1)活动导叶和转轮叶片之间产生了动静干涉效应,导致无叶区出现了能量耗散;2)优化翼型后的活动导叶附着流动不佳,在靠近导叶进口端就出现了流动剥离,从而使得在活动导叶区域就出现了大量的涡,且由于迟滞效应这些涡会一直存在到转轮区域才会逐渐消散;3)进口流量与转轮转速的匹配对湍动能会产生较大影响,从而造成一定的能量损失。这些原因共同导致了机组效率出现微小幅度下降。图 10 不同工况点下 S1流面的湍动能分布Fig.10 Turbulent kinetic energy distribution of the flow surface S1 at different operating points3.6 无叶区压力脉动频域分析 水泵水轮机的运行过程中不稳定性主要来源于机组内部水力振动,而机组的水力振动又主要来自于无叶区的压力脉动17。为了合理的解释无叶区压力脉动,引入无量纲参数 H/H 表示压力脉动的程度为:H/H=Pi-PP(3)式中:H/H 为振幅相对值,%;Pi为压力检测值,Pa;P 为监测点平均压力值。为了细致地了解到导叶翼型的优化对附近流域的影响,在无叶区选取了 4 个监测点 P1P4,进行无叶区压力检测,如图 11 所示。图 11 压力监测点示意Fig.11 Schematic diagram of pressure monitoring points图 12 展示了转轮区叶片区域内各监测点在水轮机的压力脉动时域图与频域瀑布图。此工况虽然8531第 8 期李琪飞,等:导叶翼型优化对水泵水轮机“S”特性的影响压力脉动复杂但依旧存在周期性的波动变化,改型后的压力波动 260 kPa 附近波动且相比原型波动幅度降低。频域图可以看出在监测得时间点内可以看出压力脉动多以低频率为主。整体频率成分相对简单,第 1 主频受到动静干涉影响降低从而更加接近其叶片转动频率。图 12 P1与 P3时域图与频域图Fig.12 P1 and P3 time domain plots and frequency domain plots 因为来流在进入转轮流道时与转轮叶片工作面发生冲击并产生了冲击分流,这些冲击分流使得无叶区出现了不稳定的压力脉动;另一方面,由转轮和活动导叶组成的动静叶栅干扰也是无叶区出现压力脉动的一个主要原因,转轮不断扫过由流体流经活动导叶组成的静叶栅所形成的尾迹,并且不断击打着这些尾迹,从而使得无叶区出现了脉动幅值周期性变化的压力脉动。改型后,无叶区内单位时间内涌入了更多的无撞击的有效水流,不仅最大程度上破坏了无叶区的高速水环,也减少了对转轮区域的直接冲击作用,压力幅值均有所下降,对于水泵水轮机这种经常变工况运行的水力机械来说存在一定是益处的。4 结论 1)优化翼型后的活动导叶对无叶区高速水环造成一定的破坏,使得转轮叶道部位入流状况得到一定程度改善,从而对机组“S”特性起到一定的改善作用。2)在活动导叶转动角度不变的条件下,优化导叶翼型后,入流角增大,来流液体与转轮叶片冲角减小,使得有效过流面积增大。这是改善内流特性的根本原因。3)水泵水轮机的振动主要来自于无叶区压力脉动,受动静干涉的影响,无叶区振动会引起转轮叶片疲劳损坏。优化后,压力脉动振幅出现不同程度的下降,水流自导叶出口到叶片尾端,压力脉动强度分布稳定,脉动频域成分简化,提升了水泵水轮机的平稳运行能力。参考文献:1 杨建东,曾威,杨威嘉,等.水泵水轮机飞逸稳定性及其与反 S 特性曲线的关联J.农业机械学报,2015,46(4):59-64.YANG Jiandong,ZENG Wei,YANG Weijia,et al.Runa-way stabilities of pump-turbines and its correlations with S characteristic curvesJ.Transactions of the Chinese soci-ety for agricultural machinery,2015,46(4):59-64.2 刘君,张继成,刘立峰.抽水蓄能电站调节保证设计的控制参数取值分析J.西北水电,2015(6):64-68.LIU Jun,ZHANG Jicheng,LIU Lifeng.Analysis on selec-tion of control parameters for design of regulating guarantee of pumped storage power plantJ.Northwest hydropower,2015(6):64-68.9531哈 尔 滨 工 程 大 学 学 报第 44 卷3 李琪飞,蒋雷,李仁年,等.水泵水轮机反水泵工况区压力脉动特性分析J.水利学报,2015,46(3):350-356.LI Qifei,JIANG Lei,LI Rennian,et al.Study of pump-turbines pressure fluctuations at reverse pump modeJ.Journal of hydraulic engineering,2015,46(3):350-356.4 李仁年,刘殿兴,董志强,等.水泵水轮机“S”形区全流道数值模拟J.排灌机械工程学报,2013,31(5)401-405.LI Rennian,LIU Dianxing,DONG Zhiqiang,et al.Nu-merical analysis of flow in pump-turbine operating in S-shape regionJ.Journal of drainage and irrigation ma-chinery engineering,2013,31(5)401-405.5 ZENG Wei.Extreme water-hammer pressure during one-af-ter-another load shedding in pumped-storage stationsJ.Renewable energy,2016,99:35-44.6 CHEN Zhang.Influence mechanism of geometric character-istics of water conveyance system on extreme water hammer during load rejection in pumped storage plantsJ.Ener-gies,2019,12(15):2854.7 HASMATUCHI V,FARHAT M,ROTH S,et al.Experi-mental evidence of rotating stall in a pump-turbine at off-design conditions in generating modeJ.Journal of fluids engineering,2011,133(5):1.8 XIA Linsheng,CHENG Yongguang,YOU Jianfeng,et al.Mechanism of the S-shaped characteristics and the runaway instability of pump-turbinesJ.Journal of fluids engineer-ing,2017,139(3):031101.9 ZUO Z,LIU S,SUN Y,et al.Pressure fluctuations in the vaneless space of High-head pump-turbinesa reviewJ.Renewable and sustainable energy reviews,2015,41:965-974.10 SUH J W,KIM S J,YANG H M,et al.A comparative study of the scale effect on the S-shaped characteristics of a pump-turbine unitJ.Energies,2021,14(3):525.11 ZHANG Wenwu.Pressure fluctuation and flow instability in S-shaped region of a reversible pump-turbineJ.Re-newable energy,2020,154:826-840.12 CHITRAKAR Sailesh.Review on numerical techniques applied in impulse hydro turbinesJ.Renewable energy,2020,58(6):843-859.13 张金凤,赖良庆,陈圣波,等.基于改进粒子群算法的水泵水轮机多目标优化J.华中科技大学学报(自然科学版),2021,49(3):86-92.ZHANG Jinfeng,LAI Liangqing,CHEN Shengbo,et al.Multi-objective optimization of pump turbine based on im-proved partical swarm optimization algorithmJ.Journal of Huazhong University of Science and Technology(nature science edition),2021,49(3):86-92.14 李琪飞,龙世灿,周峰,等.导叶不同开度下水泵水轮机内流特性分析J.排灌机械工程学报,2020,38(2)133-138.LI Qifei,LONG Shican,ZHOU Feng,et al.Analysis of internal flow characteristics of pump-turbine under differ-ent openings of misaligned guide vanesJ.Journal of drainage and irrigation machinery engineering,2020,38(2)133-138.15 XIAO Y X,ZHANG J,LUO Y Y,et al.Numerical analysis of misaligned guide vanes effect pressure oscillations in a prototype pump turbineJ.IOP conference series materials science and engineering,2013,52(5):234-246.16 李琪飞,张正杰,权辉,等.水泵水轮机“S”特性区的流动机理及研究进展J.热能动力工程,2017,32(9):1-6.LI Qifei,ZHANG Zhengjie,QUAN Hui,et al.Flow mechanism and research progress in“S”characteristic re-gion of pump-turbineJ.Journal of engineering for ther-mal energy and power,2017,32(9):1-6.17 王小龙,刘德民,刘小兵,等.水泵水轮机无叶区内部流场及压力脉动分析J.水力发电学报,2021,40(4):59-72.WANG Xiaolong,LIU Demin,LIU Xiaobing,et al.Analysis on flow structures and pressure pulsation in vane-less space of reversible pump turbineJ.Journal of hy-droelectric engineering,2021,40(4):59-72.本文引用格式:李琪飞,黄腾,郑非,等.导叶翼型优化对水泵水轮机“S”特性的影响J.哈尔滨工程大学学报,2023,44(8):1354-1360.LI Qifei,HUANG Teng,ZHENG Fei,et al.Effect of guide wane airfoil optimization on the S characteristics of pump turbinesJ.Journal of Harbin En-gineering University,2023,44(8):1354-1360.0631