水轮机补气条件下压力脉动的数值分析论文.pdfVIP

水轮机补气条件下压力脉动的数值分析· 钱忠东杨建东 (武汉大学，武汉430072) 【摘要】本文对某电站原型混流式水轮机进行三维非定常多相湍流模拟，对大轴中心孔补气前后水轮机各过流 部件内的压力脉动进行计算。并与试验结果相比较，探讨了水轮机补气减振的机理。分析结果表明：大轴补气能 够减小尾水管断面的压差，从而降低尾水管低频压力脉动的振幅，但对振动频率影响不明显；补气孔入口与转轮 之间存在动静干扰。在转轮区域形成新的振源，增强了转轮前的中频压力脉动。 【关键词】混流式水轮机补气压力脉动数值模拟 近年来，混流式水轮发电机组运用越来越广泛，并向大容量、高比转速、高效率的方向发展，但同 时水轮机在稳定运行方面的问题也越来越突出，包括：运行时间不长。叶片就出现裂纹【·1；部分负荷工况 下，尾水管偏心涡带导致水流的低频压力脉动，从而造成水轮机出力摆动加大刚；导水机构形成的卡门涡、 叶道涡等造成转轮叶片的疲劳破坏[习。水轮机大轴中心孔补气是常用的减振措施，国内外学者对此进行了 大量的试验研究。取得了显著的减振效果啕。 本文以某电站原型混流式水轮机为对象，通过三维非定常多相湍流计算，对大轴补气前后水轮机蜗 竟、导水机构、转轮以及尾水管内的压力脉动进行预测，采用F兀’变换对压力脉动的幅值和频率进行分 忻。探讨了水轮机补气减振的机理。 I数值计算方法 为了对气体相的运动及相间作用进行较为精确的描述。计算中采用了不可压缩流的双流体模型[71。 液相控制方程为： 连续方程 砉(a帆)+V·(稍-)=o (1) 动量方程 砉(啪蚕t)+V‘(明谕)=一叩p+V。乃+明；+瓦 (2) 气相控制方程为： 连续方程 軎(嗽)+V。(讲。)=o (3) 动量方程 砉(口∥。)+v．(嘶驴一咖+v．％+啊+九 (4) 式中乃=一凡，为相间作用力。 由于水轮机内部的流动为高雷诺数流动，本文采用RNG(RenormalizationGroup)k-s模型封闭方程网， 时间项离散采用二阶全隐格式，扩散项采用中心差分格式，对流项采用二阶迎风差分格式，应用PISO算 法求解离散方程。 在混流式水轮机中，转轮与上游的活动导叶和下游的尾水管形成了两级动静干扰，在本计算中。采用 滑移网格技术解决固定部件与旋转部件之间的信息传递。从而模拟两级动静干扰流。 ·基金项目：水资源与水电工程科学国家重点实验室开放研究基金资助项目(2004c006)。 2＆6 第一届水力发电技术国际会议论文集(第一卷) 2计算对象 图1为某电站原型水轮机各部件示意图。水轮机参数如 下：水轮机叶片数gB=15，转轮直径D，=9．7m，转速／,／ 7．854#s，活动导叶数Z高24，上游水位85m，下游水位21m。 计算中采用适应性强的非结构化四面体网格划分蜗 壳、转轮、尾水管等空间形状复杂的部件，对导水机构 采用非结构化六面体网格划分。整个计算域网格节点数 538059。单元数1972557。 3计算结果与分析 图1水轮机结构示意 3．1压力计算面和压力计算点位置 为了便于分析引水部件的压力脉动情况，分别在蜗壳进口、活动导叶前、转轮前布置了三个压力计算 点，各点位置如图2所示。为了分析补气条件下尾水管内的涡带运动情况，在尾水管锥管段布置了压力计 算平面，并在该计算平面的上游侧和下游侧布置了两个压力计算点，计算平面及各点位置如图3所示。 点 图2引水部件压力计算点位置 图3尾水管压力计算平面及压力计算点位置 3．2压力脉动计算结果及分析 本文选择导叶开度铂=190工况进行分析，分别计算了大轴中心孔补气量为进El流量的0．25％、0．5％和 0．75％的压力脉动。计算时间为6个周期，时间步长为周期的1／100。在每一个时刻，分别计算各