分享自:

水泵水轮机在飞逸工况下不稳定特性的数值研究

期刊:Renewable EnergyDOI:10.1016/j.renene.2021.01.063

该文档属于类型a(单篇原创研究报告)。以下是针对Jie Lu等人发表在《Renewable Energy》期刊上的研究的学术报告:


一、作者及研究背景

本研究由武汉大学水资源与水电工程科学国家重点实验室的Jie LuZhongdong Qian(通讯作者)联合韩国海洋大学(Korea Maritime and Ocean University)的Young-Ho Lee共同完成,发表于《Renewable Energy》2021年第169卷。研究聚焦水泵水轮机(pump turbine)在飞逸工况(runaway condition)下的非定常流动特性,旨在揭示其内部涡旋运动、压力脉动与熵生成(entropy generation)的关联机制。

学术背景

随可再生能源(风能、太阳能等)并网需求增长,抽水蓄能电站因其快速响应能力成为电网调频的关键设备。然而,水泵水轮机在涡轮模式运行时易出现S形特性曲线(S-shaped characteristic),导致飞逸工况下剧烈的压力脉动和机械振动,威胁机组安全。前人研究虽观察到该现象,但对非定常涡旋运动的能量耗散机制及其与压力波动的定量关系尚未阐明。本研究通过数值模拟与熵生成理论相结合,填补了这一空白。


二、研究流程与方法

1. 数值模型构建

  • 研究对象:中国湖北某抽水蓄能电站的原型水泵水轮机(额定转速250 rpm,飞逸转速302.3 rpm,叶片数9,导叶数20)。
  • 网格划分:采用ANSYS ICEM CFD生成结构化六面体网格(总数1145万),确保近壁面(导叶、转轮区域)细化,验证网格独立性(图6)。
  • 边界条件:螺旋蜗壳入口设为质量流量入口(mass flow inlet),尾水管出口为压力出口(pressure outlet),壁面无滑移(no-slip wall)。

2. 数值模拟方法

  • 湍流模型:采用SST k-ω模型(shear-stress transport k-omega)模拟非定常流动,时间步长为0.00055126 s(相当于转轮周期的1/360)。
  • 涡识别方法:基于Q准则(Q criterion)量化涡旋强度,定义Q>0区域为涡旋主导区(公式11)。
  • 熵生成分析:通过局部熵生成率(LEGR)解析能量耗散,区分直接耗散(粘性耗散)与间接耗散(湍流脉动耗散,公式5)。

3. 数据监测与验证

  • 监测点设置
    • 无叶区(vaneless region):导叶吸力面附近24个点,记录压力、Q值及LEGR信号(图15)。
    • 转轮区:叶片压力面与吸力面各21个点(图28)。
  • 实验验证:对比4个工况点的数值与实验数据(单位转速n11、单位流量q11),误差在可接受范围内(图7);压力脉动频谱分析验证了数值可靠性(表4)。

三、主要研究结果

1. 无叶区涡旋运动与能量耗散

  • 涡旋演化:无叶区涡旋以转轮同向(顺时针)运动(图12),形成纺锤形涡管(spindle-shaped vortex),其靠近导叶时诱发流动分离(图13c)和LEGR骤增(图14b-c)。
  • 特征频率:压力脉动、Q值与LEGR信号均出现1.4fn(fn为转频)的次主导频率(图21),与涡旋运动周期一致,表明涡旋是压力波动与能量耗散的主因。

2. 转轮区压力波动机制

  • 特征频率:转轮前缘监测点(如RP21)出现3.6fn的高幅压力脉动(图30),其与无叶区1.4fn的关系符合多普勒效应(Doppler effect):
    • 推导公式:f_runner = (v_wave - v_blade)/v_wave × f_vaneless × N_blade(N_blade=9),理论计算值3.6fn与实测值吻合(图31)。
  • 物理机制:无叶区涡旋产生的压力波(1.4fn)与转轮叶片相对运动,导致频率偏移(图33-34)。

3. 熵生成分布

  • 空间分布:转轮区熵生成占比最高(图9),主要由湍流脉动耗散主导(占比99.94%,表5),验证了飞逸工况下湍流耗散的关键作用。
  • 涡旋-LEGR关联:涡旋扩展至底环时,LEGR达峰值(图23-Ⅲ),表明涡旋尺度与能量耗散呈正相关。

四、结论与价值

科学价值

  1. 机制揭示:首次通过熵生成理论和涡识别方法,阐明了飞逸工况下涡旋运动→流动分离→能量耗散→压力波动的完整链条。
  2. 频率关系建模:发现无叶区1.4fn与转轮区3.6fn的线性关系,为后续振动抑制提供理论依据。

工程应用价值

  • 为水泵水轮机优化设计(如导叶布局、叶片倾角)提供数据支持,可减少S形区的运行风险。

五、研究亮点

  1. 方法创新:将LEGR分析与Q准则结合,定量解析能量耗散与涡旋的时空耦合关系。
  2. 发现特殊性:捕捉到1.4fn的次主导频率,修正了前人仅关注0.7fn(Hasmatuchi et al.)的局限,提出其谐波成分的物理意义(图24对比)。

六、其他亮点

  • 补充数据(附录)提供了详细的网格敏感性分析与实验对比,增强结果可信度。
  • 研究得到中国国家自然科学基金(51779186)和湖北省自然科学基金(2018CFA010)支持,体现其学术重要性。
上述解读依据用户上传的学术文献,如有不准确或可能侵权之处请联系本站站长:admin@fmread.com