混流式水轮机变工况点数值研讨.pdfVIP

混流式水轮机变工况点数值研究* 逯鹏郭鹏程彭文季 罗兴铸 (西安理工大学水利水电学院，西安710048) 【摘要】基于N—s方程，采用有限体积法和标准k—e湍流模型，对某模型转轮进行CFD分析。计算工况点涵 盖了模型特性曲线上从小流量到大流量，不同单位转速共30个工况。通过计算，获取了不同工况下转轮内部流 态，详细分析了转轮叶片上的压力和速度矢量随单位转速和开度变化的分布规律。小开度下转轮的汽蚀性能较 优，但转轮叶片背面头部靠上冠处的脱流较严重，随开度增加，该区域流态有所改善。叶片正面进水边高压区 随着单位转速的变化在上冠和下环之间转移分离，低压区主要集中于出水边；而背面低压区出现在距出水边约 1／3弦长靠近下环的区域。 【关键词】混流式水轮机流态数值计算 1前言 转轮是水轮机过流部件的核心，水轮机性能的好坏直接影响到电站的安全稳定运行。由于转轮自身的 旋转特性，及复杂外形结构，很难用实验的方法测得转轮内部流态，分析它在不同工况下的流动规律难度 更大。随着计算流体力学方法和计算机硬件性能的提高，水力透平机械三维紊流的数值模拟变得愈来愈精 确。曹树良等基于N—S方程和标准的k一￡紊流模型[1]，采用贴体坐标和交错网格系统，通过自编三维 紊流计算程序，给出了模型转轮在三个典型工况下的计算结果，初步分析了转轮内部的汽蚀和能量特性。 高忠信等对混流式轮转内部四个典型工况的三维紊流流动进行了数值模拟【2J，获得了合理的流速和压力分 布，在此基础上计算了各个工况下的转轮的水力效率，并与试验结果进行了比较，取得了一致的结果。所 以在合适的离散方法、湍流模型、网格质量以及计算域和边界条件的选取的基础上，数值计算能够较为准 确的反映转轮内部真实流态。 本文基于N—s方程，采用有限体积法和标准k—e湍流模型，运用非结构化网格和SIMPLEC算法， 对模型转轮特性曲线上30个工况点进行CFD分析，得到在各工况下，转轮内部流场结构，总结了随着单 位转速和开度变化时，叶片上压力和速度分布规律。 2数值计算方法和计算模型 2．1控制方程 在笛卡尔坐标系中，水流运动规律可以用下述瞬态的连续性方程和动量方程描述，假设在流动过程中 没有热传导。 连续性方程 生Dt+V·(』D方)=o N—S方程 —+ 式中 u——速度5 p——密度； 卢——动力黏性系数5 ·国家自然科学基金资助项目。 第一届水力发电技术国际会议论文集(第二卷) 占——克罗内克系数； SM——动量源项； p——压力。 ^ ^．=} 水力机械内部水流运动可假设为不可压定常流动，则上式中P密度为常数，且宴：掣：0。将 at Ot 上述方程做时均化处理后，引人了雷诺应力项，使的方程组不封闭，故引入湍流模型。由于标准k—e湍 流模型比零方程和一方程模型有较大的改进，在科学及工程实际中得到非常广泛的应用，并对于水轮机内 部流态的模拟可以获得相对准确的数值解，它对于计算机硬件资源和计算时间的要求也较高级的湍流模型 低，故本文选用标准k一￡模型。控制方程的离散采用有限体积法。 2．2计算流域 考虑到计算机硬件资源和计算时间的限制，并为了获得相对准确的结果，本文的计算分两步进行。首 先计算由固定导叶、活动导叶和转轮及其出口延伸组成的单周期通道，得到转轮进口速度分布，将其作为 转轮和尾水管联合计算的进口条件。计算模型如图1所示。各部分网格节点数如表1所示，对于活动导叶 不同开度下网格质量一致。模型转轮直径0．34m，叶片数为13个，固定导叶和活动导叶数均为24。 图1计算域模型网格 (a)固定导叶、活动导叶和转轮的周期通道网格；(b)转轮尾水管通道网格 表1 各部分通道网格节点数 计算步 第一步周期通道模型