尾水管非定常流动模拟及不规则压力脉动预测.docVIP

尾水管非定常流动模拟及不规则压力脉动预测 桂中华1 唐 澍1 潘罗平1 韩凤琴2 1中国水利水电科学研究院，北京，100038 2华南理工大学电力学院，广东，广州 510640 摘要：尾水管内螺旋形涡带引起的压力脉动是造成混流式水轮机组振动的主要根源之一，直接威胁着机组的安全运行。本文提出了一种基于CFD技术的水轮机尾水管压力脉动数字化预测法，并利用此法对一大型混流式水轮机典型偏工况下尾水管内流动进行了长时间非定常计算，详细讨论该工况下尾水管内死水域与涡带的运动规律，预测了尾水管不规则压力脉动，压力脉动分析结果表明其波形、频率、相位与实际基本一致，证明预测方法是可行的。 关键词：尾水管 非定常流动 压力脉动预测 前言 水轮机尾水管是能量回收的重要部件，对机组的整体能量特性和稳定运行具有很大的影响，其内部水流从垂直方向转向水平方向，流动受离心力的作用而存在二次流，且过流断面沿流向存在扩散、收缩、再扩散的过程，其流动复杂，常常产生局部脱流和回流等现象。尤其在偏离最优工况运行时，进入尾水管的流动更加复杂，水流夹带着空化气泡在离心力的作用下形成同水流共同旋进的尾水管涡带，涡带在周期性非平衡因素的影响下产生偏心，这种偏心涡带大大降低水轮机效率,其诱发的压力脉动频率接近机组的某一个固有频率时,将会引起强烈共振，威胁机组运行的安全性。尾水管的压力脉动特性是水轮机振动与稳定性的重要评价指标，研究其特性对于解决水轮机稳定性具有直接意义。在过去的几年中，与混流式水轮机稳定性运行有关的压力脉动的研究取得了一定进展.文献[1]对模型水轮机进行了试验研究，分析了尾水管内部不同位置不同工况下压力脉动变化规律，为改善尾水管设计提供了参考。文献[2][3]对单个尾水管进行了内部非定常流动解析，仿真出了部分工况下尾水管内的涡带,但其计算时间较短。本文旨在通过对联合转轮的混流式水轮机尾水管进行长时间的非定常流解析，进而预测尾水管内死水域与涡带及其不规则的压力脉动，为提高机组运行稳定性提供可靠的技术保障。 压力脉动预测法 压力脉动预测法以CFD技术为基础,首先对水轮机全流道进行三维造型,并选择恰当网格尺度对水轮机各个部件进行网格划分，接着选用合适的湍流模型进行定常、非定常数值解析，最后提取测点的压力脉动值预测脉动变化趋势。具体方案如图1所示。 数值解析是该预测法的基石，在对尾水管流动进行解析时湍流模型的选取尤为重要。研究表明粘性涡模型（k -ε两方程模型）计算得不出由雷诺应力引起的二次回流，因为该模型采用了各向同性假设，在动能方程和耗散方程的模化过程中，将湍流粘性系数处理为各向同性的标量。这使得标准ε模型无法计入离心力、哥氏力对湍流场的影响,Ω模型本质上与k -ε模型 数值解析设置 3.1 计算域与网格 过去在计算资源不能满足的 这种解析策略虽然也能仿真出尾水管涡带运动，预测出压力脉动，但是还存在下列不足：1)虽然考虑了转轮出口流动对尾水管的影响，但这种方法只考虑了某一个时刻的影响，而实际上是转轮出口流动每时每刻都在影响着尾水管内的流动情况；2)转轮与尾水管间存在的干涉是相互的，上面的处理方法仅仅考虑转轮对尾水管的部分影响而根本没有涉及尾水管对转轮的影响。因此上面的解析方案不是最佳方案。但是一台混流式水轮机真机，几何尺寸非常大，把所有的过流通道，包括、固定导叶、活动导叶、转轮及尾水管一起进行整体解析，要想获得比较理想的结果，则要求网格的数量非常巨大，并采用周期性边界条件；(解析域)见图2。解析域各组件采用了六面体结构化网格（图3），解析总网格数近50万，其中固导活导10万、转轮叶片10万、尾水管30万。 3.2 边界条件及工况参数 本研究对象为一混流式水轮机，转轮直径为5.52m，固定导叶数12，活动导叶数24，叶片数13，机组水头57.2m。计算工况参数见表1。 n11(rpm) Q11(l/s) Q11 /Q11opt H(m) Q(m3/s) ag(m) 73.00 900 88.7 57.2 207.56 0.2913 进口条件：计算的入口设在蜗壳与固定导叶交界处，采用流量进口条件，流量可以从模型特性曲线相应工况点参数进行计算得到，方向结果计算得出。出口条件的压力值按平均静压为零给定。 100rev/min速度绕中心轴旋转。 交接面：活动导叶叶栅与转轮的交界面、转轮和锥管的交界面为滑移交界面，采用stage类型。 周期面：在导叶与转轮叶片两侧设置周期性边界。 计算步长：非定常解析的时间步长为?T=0.005s，为转轮转动周期的1/120，每计算4步保存一个数据，即数据采样周期为?Ts=0.02s。为了正确预测尾水管内流动的不规则行为，本文进行了物理时间80s的非定常计算。 解析结果分析 4.1 尾水管内死水域与涡带的形成 从计算结