基于FLUENT的库区涌浪数值模拟-水利水运工程学报.PDF

第3期 水 利 水 运 工 程 学 报 No.3 摇 2014年6月 HYDRO鄄SCIENCE AND ENGINEERING Jun. 2014 基于FLUENT 的库区涌浪数值模拟 邓成进, 袁秋霜, 侯延华, 贾摇 巍 (中国水电顾问集团西北勘测设计研究院有限公司, 陕西 西安摇 710065) 摘要:基于流体计算软件FLUENT,模拟某水电站库区近坝变形体可能失稳后下滑引起库区水面变化过程,分 析初始涌浪形成以及涌浪在对岸爬坡和涌浪沿岸传播的过程;研究挡水建筑物对库区涌浪沿岸传播的影响,得 出初始涌浪高度,以及对岸、沿岸、坝址处的最大浪高,并与潘家铮法估算结果进行对比分析。 分析结果表明, 数值模拟能较好反映波浪爬坡和沿岸传播过程,真实模拟库区水体相互作用;由于库区涌浪运动受大坝建筑物 阻挡作用,库区水面的反复震荡和涌浪叠加,会形成更高的涌浪。 计算的初始涌浪及库区各处的最大涌高更符 合实际情况,可为近坝库区的工程设计及涌浪灾害的预防提供参考。 关摇 键摇 词:库区;滑坡体;爬坡过程;最大浪高;涌浪叠加;挡水建筑物;数值模拟 中图分类号:TV697;O242摇 摇 摇 文献标志码:A摇 摇 摇 文章编号:1009-640X(2014)03-0084-08 库区存在的滑坡、崩滑体等不良地质体在失稳后高速滑入库中,产生的高速波浪足以对沿岸及下游建筑 [1] 物和居民生命财产构成巨大威胁。 1963年发生在意大利瓦依昂水库的滑坡 、1982 年发生的鸡扒子滑坡、 [2] 1985年发生的新滩滑坡 、2003年发生的千将坪滑坡等库区滑坡失稳后激起的巨大涌浪,均造成了巨大的 人员伤亡和财产损失。 因此,研究库区滑坡涌浪受到国内外学者的广泛关注。 但是,由于库区涌浪受多种因 素影响且十分复杂,涌浪形成的边界条件和初始条件难以明确定义,滑坡涌浪的计算还没有一种通用的方 [3-5] 法。 目前,研究库区滑坡涌浪的手段主要有物理试验、经验计算公式和数值模拟等 。 其中,经验公式法 主要有:Noda 法、潘家铮法、中国水利水电科学研究院经验公式法、美国土木工程学会法、Slinger鄄land [6] Voight 公式等。 物理模型试验也取得一定的经验,殷坤龙等 采用物理试验方法研究了三峡水库库岸滑坡 问题。 随着计算机及数值理论的迅速发展,采用FVM,FEM,VOF等方法能较好地模拟流体飞溅、融合等复 [7-9] 杂自由表面现象,因此,基于流体力学的涌浪数值模拟有了一定程度的研究 。 然而,现有数值模拟方面 的研究主要集中在模拟滑块垂直入小水槽内,而对在宽广水域中的涌浪问题研究较为欠缺,也很少有尝试通 过数值模拟研究挡水大坝对库区涌浪传播的影响。 本文采用FLUENT 软件建立某库区近坝库岸变形体失稳后产生涌浪的数值模型,模拟滑体滑动后库区 涌浪形成及传播的过程;计算分析了大坝对库区涌浪传播过程的影响,得出初始涌浪高度、涌浪沿岸的传播 高度以及涌浪到达大坝时的爬坡高度,并将其结果与潘家铮法估算的涌浪高度进行对比分析,为工程设计提 供参考。 1摇 基于FLUENT软件的数值模拟方法 FLUENT软件是一种基于有限体积法发展起来的通用商业流体计算动力学软件包,用于模拟从不可压 缩到高度可压缩范围内的复杂流体运动,在水工泄水建筑物水力学计算方面应用广泛。 1.1摇 控制方程 对于二维的瞬态流,其基本控制方程包括:连续性方程和以速度和压力为变量的动量守恒方程。 摇 摇 收稿日期:2013-10-29 摇 摇 作者简介:邓成进(1986-),男,湖北随州人,助理工程师,硕士研究生,主要从事水工设计及边坡数值模拟