轴流泵不稳定流场压力脉动特性研究.pdf

水 利 学 报 年 月 第 卷 第 期 文章编号 轴流泵不稳定流场的压力脉动特性研究 王福军 张 玲 张志民 中国农业大学水利与土木工程学院 北京 中国水利水电科学研究院 北京 摘要 流场压力脉动是影响大型轴流泵运行稳定性的关键因素 本文采用时间相关的瞬态流分析理论及大涡模拟 方法研究轴流泵内部非定常流动 得到了不同工况下泵内水压力脉动结果 通过与实测扬程和功率对比 证明本 文所提出的方法可较准确地反映泵的流动特征 研究表明 轴流泵内最大压力脉动发生在叶轮进口前 压力脉动 频率主要受叶轮转频控制 在叶轮进口与出口处 从轮毂到轮缘压力脉动逐渐增大 而在导叶中间及导叶出口处 结果正好相反 偏离最优工况越远 脉动的相对振幅越大 在 流量工况下泵内压力脉动约为最优工况的 倍 关键词 轴流泵 不稳定流场 压力脉动 大涡模拟 中图分类号 文献标识码 转动的叶片与静止的导叶间的相对运动 偏离最优工况时吸水室水流圆周运动 局部空化及二次流 等因素 都可能导致泵内水压力随时间不断快速变化 即出现压力脉动 压力脉动规律可能是周期性 的 也可能是随机的 压力脉动的最大危害是使机组的结构振动加剧 即出现水力激振 同时还可引发 进一步的局部空化 甚至在某些情况下引起共振 水力激振问题是大型泵站水泵机组所面临的突出问 题之一 大型轴流泵在我国长江中下游地区的排涝泵站以及正在建设的南水北调东线工程中广为使 用 而到目前为止 国内外对轴流泵压力脉动特性的研究几乎是空白 因此 出于提高轴流泵运行稳 定性 保证机组在正常安装高程下不发生空化的需要 对轴流泵压力脉动特性的研究显得尤为重要 对于泵内非稳定流场的压力脉动特性研究 主要有实验和数值计算两种手段 由于成本过高 采用 实验手段监测泵内 特别是转动的叶片上 压力脉动通常比较困难 此外 压力脉动及水力激振特性经 常不具相似性 使得在实验室内采用模型泵进行压力脉动实验研究的意义大打折扣 因此 多数情况 下只能在水泵外特性实验结果的辅助下进行内部流场的三维数值模拟 随着计算流体动力学 技术的发展 水力机械流场计算目前已发展到全三维 黏性 非定常阶段 国内外已有相当深入的研究成 果 但到目前为止 对轴流泵的流场模拟还基本停留在定常计算阶段 为此 本文拟采用大涡模拟 方法对轴流泵三维非定常湍流进行数值模拟 以期揭示泵内压力脉动特征 轴流泵大涡模拟计算模型 除近壁面区域外 轴流泵内的水流一般处于湍流状态 湍流包含一系列大小不一的涡团 大涡模 拟的思想是将比网格尺度大的涡团通过瞬时 方程直接计算出来 而小尺度涡对大尺度涡 运动的影响则通过一定的模型在瞬时 方程中体现出来 为此 使用大涡模拟前 需要先 建立数学滤波函数 以便从湍流瞬时运动方程中将尺度小的涡滤掉 形成大涡流场运动方程 这里 采 收稿日期 基金项目 北京市自然科学基金项目 国家自然科学基金项目 教育部新世纪优秀人才支持计划项目 作者简介 王福军 男 河北唐山人 博士 教授 主要从事水力机械研究 用如下方式来构造滤波函数 ? 式中 为实际流动区域中的空间坐标 为滤波后的大尺度空间上的空间坐标 为滤波函数 为控制体积所占几何空间的大小 为控制体所在的计算域 用式 表示的滤波函数处理瞬时状态下的 方程及连续方程 可得到 式 和式 构成了大涡模拟方法所使用的控制方程组 式中带有上划线的量为滤波后的场变量 为亚格子尺度应力 它体现了小尺度涡的运动对所求解的运动方程的影响 大涡模拟计算的关键是如何表示上述应力 即如何选取或构造亚格子尺度应力模型 本文选取文 献 中的方法来确定亚格子尺度应力 式中 是亚格子尺度的湍动黏性系数 在