电站进水口流场数值模拟.pdfVIP

1 电站进水口流场的数值模拟 陈云良，伍 超，叶 茂，李 静 四川大学高速水力学国家重点实验室（610065） email: liangyunchen@163.com 摘 要：本文采用三维k −ε 紊流模型模拟水电站进水口流场。对进水口前缘水面和水下流 态，以及喇叭形进口段流场进行了研究。计算采用原型比尺，并与 1：50 的模型试验作了对 比分析，结果表明：进水口前缘的回流只存在于水面，并未延伸至水下，进口的佛汝德数较 小，速度环量在回流区从外到内衰减较快，不具备形成立轴旋涡的条件；喇叭形进口段的水 流顺畅、进水匀称。 关键词：进水口，速度环量，立轴旋涡 1. 引言 保证进水口水流顺畅、进流匀称及流态平稳，并尽量减小水头损失，是水电站进水口设 计的主要目标。立轴旋涡是水电站进水口前常见的水力现象，它的存在将恶化进口流态，降 低过流能力，减少发电量，甚至造成水工建筑的振动和空化破坏。 立轴旋涡的机理复杂， 自由面与水工建筑物等边界条件多变，理论研究进展甚微[1]-[4], 现阶段主要还是采用实验方法进行研究。进水口立轴旋涡受边界条件的影响较大，这已是研 究者们的共识[5]- [7]，而在模型试验中，边界条件的模拟不可避免地有一定的近似性。另外 由于立轴旋涡流场的复杂性，相似律目前尚无统一的观点，多数情况下利用佛汝德数相似设 计，认为雷诺数超过一定极限值后，粘滞力的影响就可以忽略不计[8]，但吸气旋涡并非总是 与雷诺数无关。王水田[9]在某水电站进水口以佛汝德数相似设计的模型试验中只发生了表面 凹陷旋涡现象，而原型观测中却出现了明显的漏斗旋涡，说明实验结果受比尺效应的影响。 因此模型实验中未出现立轴旋涡并不等价于原型中不发生立轴旋涡现象。 与模型试验相比，数值计算有花费少、变方案快、信息完整、模拟能力强的优势，而且 不存在比尺效应的问题。但由于立轴旋涡现象的运动比较复杂，旋涡中心位置不确定，其尺 度一般较小，而速度梯度的变化大,对数值计算方法、局部网格的高度细化以及计算机的处 理能力都是一种巨大的挑战，尤其是自由水面的模拟。Lai[10]采用刚盖法近似处理自由水面， 对水泵进水口的模型进行了三维数值模拟。模型试验和数值计算相结合是研究立轴旋涡的有 效方式，在Wawushang水电站可行性研究的进水口水工模型试验中，库水位较低时有间歇性 的旋涡出现，初计阶段通过进水口位置的重选和进口体型的修正，作者对库区进水口流场进 行了三维数值模拟和 1：50 的模型试验。三维数模以分层两相流的k −ε紊流模型为理论基 础；运用 VOF 模型处理复杂的自由水面；网格剖分采用非结构化混合网格，以灵活地适应极 1本课题得到高等学校博士学科点专项科研基金（项目编号：20020610016）资助。 -1- 不规则的地形和复杂的进水口体型，较精确地模拟边界。通过对进水口原型的模拟，提供了 完善的流场数据。 2. CFD 模型 分层两相流的k −ε紊流模型的连续方程、动量方程和k 、ε方程可分别表示如下： ∂ρ ∂ρui 连续方程： + 0 （1） ∂t ∂x i ∂ρu p ⎡ ∂u ∂u ⎤ i ∂ ∂ ∂ i j 动量方程： ρu u µ µ