浅析轴流式水轮机模型飞逸过程三维湍流数值模拟.pdf

水 利 学 报 2010年 2月 SHUILI XUEBAO 第4l卷 第2期 文章编号：0559—9350(2010)02—0233—06 轴流式水轮机模型飞逸过程三维湍流数值模拟 周大庆 ，吴玉林 ，刘树红 (1．河海大学 能源与电气学院，江苏 南京 210098；2．清华大学 水沙科学与水利水电512程国家重点实验室 ，北京 100084) 摘要：通过作者建立的非稳态过程三维湍流数值模拟方法，模拟了轴流式水轮机模型飞逸过程，并进行了其三维 数值湍流计算 ，获得 了机组到达最大转速时所需时间、最大飞逸转速值、飞逸过程 内特性流场演变规律 以及转 速 、流量、力矩与测点静压随时间变化的曲线 。计算结果表 明，所得最大飞逸转速与试验值非常接近 ，误差在 1．5％以内；所得力矩及流量等外特性参数变化过程均符合高 比转速水轮机理论特性 ；当转速增大到某一数值之 后 ，尾水管水力稳定性急剧恶化是引起机组振动的主因。 关键词 ：轴流式水轮机模型；飞逸过程 ；三维非稳态湍流数值模拟 中图分类号：TK733-3 文献标识码：A 1 研究背景 当水轮发 电机组处于飞逸状态时，由旋转体 中的质量不平衡以及机组过流部件 中出现的水力不 稳定现象所产生的动载将大大增加 ，其所引发 的振动对水轮发 电机组会造成极大的威胁 。根据对水 轮机边界的处理方式不同，当前常用 的数值解法有 以下两大类 ：(1)外特性数值解法 ，以完整的水轮 机综合特性 曲线或全特性曲线作为求解过渡过程的边界条件 ，联立机组转动部分运动方程式和一 元非恒定管流基本方程组 ，然后通过特征线法求解 ；(2)内特性数值解法 ，只需 已知水轮机及其装 置的几何参数和基本结构参数以及过渡过程计算 的初始条件 ，根据 叶片式水力机械的广义基本方程 式建立不稳定工况下水轮机动态力矩和动态有效水头表达式 ，在调节元件运动规律 已知的情况下 ， 与有压非恒定管流基本方程组联立求解 。上述两种方法在解决包括 引水系统在 内的整个水力系统过 渡过程问题 中发挥着重要作用 一，但它们在精细反映水轮机 内流特性规律方面存在一定不足。因此 本文通过建立一种非定常三维湍流数值模拟方法 ，以研究水轮机的飞逸过程。该模型可详细揭示机 组转速 、流量 以及相关测点压力随时间变化规律 ，获得到达最大转速时所需时间，及最大转速数值 大小 ，从而为机组结构设计和运行管理提供参考 ，并为寻求有效的防飞逸措施及降低飞逸最大转速 提供科学依据。 2 飞逸过程三维湍流数值计算方法 2．1 控制方程 针对飞逸过程中水轮机运行特点，转轮室区域采用 以下流动控制方程 连续性方程 ： (一 )=0 收稿 日期 ：2008—11-26 基金项 目：国家 自然科学基金资助项 目90410019)；河海大学 自然科学基金青年项 目 作者简介 ：周大庆 (1976一)，男，浙江泰顺人，讲师，主要从事流体机械及工程研究。E-mail：zhoudaqing@hhu．edu．cn - - - — — 233 -·—·— 式中： 为水流绝对运动速度 ；五为转轮区网格运动速度。 动量方程为 +(一)鲁=一吉罢+91ZZli (2) 式 中： 为单位体积力 ；P为密度 ；P为压强 ；l，为水体运动黏性系数。 水轮机飞逸过程 ，是一个不断加速的过程 ，转轮区网格运动速度 由转轮旋转平衡方程确定， ．， (3) 式 中： 为合力矩 ，本计算忽略机械摩擦力矩与风阻力矩 ，单位为N·nl；．，为机组转动惯量 ，单位为 ·m2；cc，为转轮旋转角速度 ，单位为rad／s。 则在笛卡尔坐标系下，转轮区网格运动速度分量为：