水平旋流式内消能泄洪洞数值模拟分析.pdfVIP

第 42卷 第 5期 2011年 3月 人 民 长 江 Vo1．42．No．5 Yangtze River Mar．， 2011 文章编号 ：1001—4179(2011)05—0078—04 水平旋流式内消能泄洪洞数值模拟分析 邬 海 波 ，赵 振 兴 (1．河海大学 水利水 电工程学院，江苏 南京 210098； 2．河海大学 力学与材料学院，江苏 南京 210098) 摘要 ：为 了探 究水平旋流式 内消能泄洪洞的消能机理及相关水力要素的变化情况，采用 RNGk—s双方程紊流 模型及基于水气两相流的VOF方法，并结合公伯峡水 电站的原型观测资料 ，得到 了洞 内流速和压强等水力要 素的分布规律。数值模拟结果与公伯峡水电站原型观测数据 的对 比分析表明，数值模拟方法能够较好地模拟 水平旋流泄洪洞的泄流过程 ；通过改变旋 流室入流 角度 ，得 出了压强 、通 气量和 消能率的变化规律 。 关 键 词 ：紊流模型 ；数值模拟 ；水平旋流 ；入 流角度 ；消能率 中图法分类号：TV651．3 文献标志码 ：A 水平旋流式 内消能泄洪洞作为一种新型的消能 工 ，有其 自身独特的优势 ，通过水流贴洞室壁高速旋转 流动并和空气混掺，可有效减免空化空蚀 ，同时利用导 流洞泄洪可提高经济效益 ，洞室中形成的空腔环流增 加了水流的紊动和水流之间的剪切作用 ，增大了消能 率。 目前旋流式消能工在模 型实验方面有很多研究 ， 图 1 水平旋流泄洪洞体型布置 (单位 ：m) 包括流速 、压强 、旋流室空腔直径 、旋流夹角、泄流量和 消能率等。 本文用数值模拟 的方法对其流场进行计算分析 ， 并通过与黄河公伯峡水 电站原型观测资料进行对 比验 证 ，结果与原型资料基本符合 。 1 工程概况 公伯峡水电站位于黄河干流上，装机容量 1500 MW，最大坝高 139m，总库容6．2亿m ，是国内第 1个 泄量超过 1000m ／s，水头超过 100m 由右岸导流洞 改建成水平旋流泄洪洞 的大型水利工程 。2006年的 图2 水平旋流洞计算网格 (单位 ：m) 原型观测表 明，工程获得圆满成功。 本文首先选取竖井段至旋流洞段进行计算 ，上游 2 数值模拟 水位选取2004m，下游为无压出流。接着选取竖井段 2．1 数值计算与原型 比较 至导流洞0+420m处作为另一计算体型，工程体型依 据公伯峡电站设计资料 ，见图 1，计算离散区域见图 2．1．1 紊流模型 2 采用RNGk—s双方程紊流模型，针对复杂曲线 自 收稿 日期 ：2011—01—24 作者简介：邬海波，男，硕士研究生，主要从事流体力学的研究工作。E—mail：wuhaibo8@163．corn 第 5期 邬海波 ，等：水平旋流武 内消能泄洪洞数值模 拟分析 79 由水面和边界条件 ，利用结构网格来剖分主要三维计 占计算区域体积的相对 比例。对于某个计算单元，若 F 算区域 ；引入适用于分层两相流的VOF方法求解 曲线 = 1，则表示该单元被流体全部充满；若 F=0，则表示 自由水面。k—s双方程模型对雷诺应力各分量采用各 它是一个空单元 ；若0F 1，则表示该单元部分充 向同性假设 ，对于各 向异性较强 的螺旋水流 的模拟不 满流体。显然 ，对于 自由表面问题 ，自由边界存在于第 如雷诺应力模型或代数应力模型精确 ，但通过对湍动 3种单元中。F的梯度可以用来确定 自由边界的法