竖井旋流式泄洪洞数值模拟.pdfVIP

水 利 学 报 2003 年08 月 SHUILI XUEBAO 第8 期 文章编号：0559-9350 （2003）08-0058-06 竖井旋流式泄洪洞数值模拟 1 1 1 1 张晓东 ，刘之平 ,高季章 ，王晓松 (1.中国水利水电科学研究院 水力学研究所，北京 100038) 摘要：本文采用k- ε双方程紊流模型及基于水气两相流的VOF 方法，数值模拟竖井旋流泄洪洞复杂的水流运动，得 到自由水面，流速及压力等水力要素的分布规律。数值模拟结果与已有模型试验数据的对比分析表明，采用的数值 方法能够较好地模拟竖井旋流泄洪洞泄流过程。此外，利用数值模拟的结果对竖井的泄流能力和消能率进行了验证。 关键词：竖井旋流；消能；数值模拟；紊流模型；VOF 方法 中图分类号：TV135.2 文献标识码：A 与传统的竖井溢洪道相比，竖井旋流式泄洪洞是一种新型消能工，由于水流绕竖井轴心旋转呈贴壁流 动，在中央形成稳定空腔，既保持水流流态的稳定又形成良好的通气条件，避免了气囊失稳现象的出现， 在适当的设计体型下，还能有效地消除部分水流能量。对于将导流洞改建为泄洪洞的工程，竖井旋流式消 能工具有其独特的优势，经济效益显著。 目前对旋流式消能工性能的研究方法，主要采用物理模型试验和数值模拟。模型试验多为常压下进行， 通过对旋流式内消能工的流态、压力、流速、泄流能力、消能率等方面的分析，探讨其消能机理。由于竖 井旋流的复杂性，用数值模拟法研究旋流式内消能工的文献较少。本文结合沙牌水电站竖井式泄洪洞的试 验结果[1]，对所研究体型进行数值模拟，其结果与模型试验结果基本相符。 1 旋流式泄洪洞的数值计算 1.1 工程简介 沙牌旋流竖井泄洪洞由压力短进水口、引水道、涡室、竖井和导流洞组成。1995 年中 国水利水电科学研究院水力学所曾就其体型优化进行了系列模型试验，所推荐的体型如图1 所示。 本文根据推荐的体型确定计算域为引水道上游断面至桩号 0+235m 处导流洞断面，不包括压力短进水 口。就上游水位1866m，下游导流洞明流无压出流工况进行数值模拟。计算域及计算域离散网格见图2。 1.2 数值计算 1.2.1 紊流模型 采用 k- ε双方程紊流模型，针对复杂曲线自由水面和边界条件，利用非结构网格来剖 分三维计算区域；引入适用于分层两相流的VOF 方法求解曲线自由水面。k- ε双方程模型对雷诺应力各分 量采用各向同性假设，对于各向异性较强的螺旋水流的模拟不如雷诺应力模型或代数应力模型精确，但根 据文献［2］对阶梯溢流坝的计算经验，采用VOF 方法由k- ε模型计算比较符合实际情况。标准k- ε紊流模 型的连续方程、动量方程和k、ε方程分别为：连续方程： ∂ρ ∂ρui + 0 （1） t x ∂ ∂ i 收稿日期：2003-03-28 作者简介：张晓东(1976-)，男，山东寿光人，在读博士生，主要研究方向为高速水力学与泄洪消能。 58 水 利 学 报 2003 年08 月 SHUILI XUEBAO 第8 期 图1 竖井结构布置(单位:m) 图2 竖井计算网格示意 动量方程：