轴流式水泵非定常湍流数值模拟若干关键问题.pdfVIP

·泵及榘装置· 轴流式水泵非定常湍流数值模拟的若干关键问题 王福军，张玲，黎耀军，张志民 (中国农Nk大学水利与土木工程学院北京100083；中国水利水电科学研究院北京 100044) 摘要；为准确捕捉轴流式水泵不稳定流场特征，采用霄诺时均法和大涡模拟方法对轴流泵非定常湍流进行了数值 模拟，并就相关问题进行了深入研究。研究发现，与其他湍流分析方法相比较，大涡模拟方法在轴流泵非定常流场 分析巾具有更高的计算精度：同时将轴流泵进水流道与泵段一起进行流场计算后。发现叶轮进口的流动是极不均匀 的，靠近弯肘形流道内侧的速度明显大于其他部位，这一现象与通常假定的轴流泵进口流动是均匀的设计与分析理 论有很大区别；研究还发现，为了获得轴流泵流场内不同压力脉动频率成份，脉动计算的采用样时间应不少于8个 旋转周期；轴流泵内部的水压力脉动幅值在叶轮附近区域沿半径逐渐增大，而在导叶出口之后沿半径逐渐减小。 关键字：轴流泵；进水流道；非定常计算；压力脉动；采样时间 O前言 随着南水北调工程及国家中部地区大型排涝泵站更新改造工程的建设，大型轴流式水泵在我国 的应用越来越广，随之而来的运行稳定性问题日益突出，其中最严重的是机组在某些工况下出现较 大的水压力脉动和相应的水力激振Ⅲ。 日本茬原公司忙1和德国KSB公司p1分别通过实验研究和理论计算，发现轴流泵装置系统的水压力 脉动主要来自于进水流道不稳定流动、叶轮与导叶的动静耦合，以及叶尖二次流。而目前对这种复 杂的非稳定流场的研究还很不深入。随着计算流体动力学(CFD)的发展，三维非定常湍流分析方 法为揭示复杂流动机理提供了一种新的途径。基于雷诺时均法(RANS)及大涡模拟法(LES)的湍 流处理模式已开始应用于流体机械非稳定流场预测及优化设刊Ⅷ。但对轴流式水泵而言，现有的流 场分析多数基于定常理谢71。而定常计算很难揭示泵内压力脉动和二次流，以及由此引发的结构振动 和空化加剧等现象佴1，因此，急需开展轴流泵非定常湍流分析。 本文采用不同湍流模型对轴流泵三维不稳定流场进行了非定常计算，预测了泵内水流压力脉动，并 利用快速傅立叶变换对压力脉动频谱特征进行了分析，本文拟就其中的某些技术关键做一总结和探讨。 1湍流模型与计算方法的选择 水力机械中的湍流运动可用连续方程和Navier-Stokes方程来描述阳川： 望+v．(p跖)：o (1) 一 拼 旦生刍堕+v．(pⅣ即)：,oF—vp+必u (2) 0l 式中，p为水的密度，Ⅳ为速度矢量，，为质量力，P为压力，·∥为动力粘度。式(2)左边 第一项为由于运行非楚常性而引起的局部惯性力，第二项为由于运动的非均匀性而引起的变位惯性 力，右边三项分别为质量力、压力和粘性力。通常可将水流按不可压缩处理，故方程可以简化。 对于湍流的处理方式．也即式(2)中最后一项的计算。是决定流场数值模拟精度与效率的关键， 支持计划． 2009‘辛日^，《女≈＆mⅢt0饪新技术、新P∞交m^0le』《 H前辛婴打雷诺时均法和大涡模拟法州种…。各种^q仅方程授删以及雷辞；应力模+¨是衙讲时均法 的几种小同}B式，小文选抒标准t—s双^程模型”-为雷潇叫均法的代表，连同太描模拟方法一计对 轴漉泵0l的非稳定流功进行计算。 计算对象是一叶轮直径为300nm·的轴流泵，叶片数Z-6，导叫叶片数z=II，慢忱漉量口= 350协，扬程Ⅳ=1Iflm转速肟I450dmin。计算时考虑r进水流道的影响将进水流通勺叶轮 导af 弯臂等一起柞=为计算区域，如阿1所H；。 。圈 目l#沉“《H“#《≈ 采片j有限体积法对控制方程在字瑚域i一进行离散““，控制方靴中的糠项和扩敞项应用二阶中心 差分格式，对流项f“用二阶迎风格式，在时间域r．采用二阶争隐式{}}式进行离散。每个时问步内 采用PISO片法“o