溃坝水流数值模拟的研究进展.pdfVIP

第 17 卷 第 1 期 水 科 学 进 展 Vol17 ,No1 　2006 年 1 月 　 ADVANCES IN WATER SCIENCE 　 J an . , 2006 　 溃坝水流数值模拟研究进展 史宏达 , 刘 　臻 ( 中国海洋大学工程学院 , 山东 青岛　266071) 摘要 : 溃坝问题在水利工程的设计管理中具有重要地位 , 也是广大学者长期以来一直关注和研究的课题 。回顾和总 结了国内外对溃坝水流演进问题的研究进展 : 介绍了溃坝水流的数学模型及解析解法存在的困难 , 进而讨论了数 值解法的最新进展 ; 论述了求解溃坝水流一维问题的有限差分法 、近似黎曼解的 Godunov 格式法 、Boltzmann 法 、 KFVS 法和二维问题的 TVD 格式法 、间断有限元法 、有限体积法 、特征线法 , 并分析了各种方法的适用范围和优缺 点 , 及讨论了限制函数的使用 ; 介绍了利用自由水面追踪方法计算溃坝水流的研究进展 , 并根据 目前存在的不足和 实际工程的需要 , 提出了进一步研究的方向和发展趋势 。 关 　键 　词 : 溃坝水流 ; 数值模拟 ; 研究进展 ( ) 中图分类号 : TV139231 ; G35311 　　　文献标识码 : A 　　　文章编号 : 2006 在水利工程的设计和管理中预报坝体溃决这一灾害性水流现象十分重要 。溃坝计算是对水库和堤防的失事 影响做出定量估算 , 并合理确定水库或堤坝防洪设计标准以及避险措施的有效手段 。溃坝计算的主要内容是算 出溃坝坝址的流量和水位过程线 , 以及下游洪水演进过程中沿程各处的流量 、水位 、流速 、波前和洪峰到达的 时间等 。溃坝水流的构成复杂 , 通常包含激波 , 亚临界流 , 超临界流等区域 。通过数值解与试验数据比较 , 认 为浅水方程能够较好的描述溃坝水流 。问题最终归结为求解控制水流运动的非恒定流拟线性双曲型偏微分方程 组的有间断问题 。 1 　溃坝水流的控制方程 对于一维溃坝问题 , 瞬间全溃引起的不稳定流动可视为一维流动 , 如果假定为静水压力分布和小底坡 , 则 可用圣维南方程[ 1～4 ] 描述 。对于二维溃坝问题 , 在静压假定和忽略风应力和柯氏力的条件下 , 描述溃坝洪水演 进的二维控制方程为浅水方程[ 1～4 ] 。 2 　溃坝问题的理论解法 尽管数值计算方法在模拟溃坝水流运动方面取得了一定成功 , 但其效果及优劣通常需要将其与解析法所得 结果作比较来判断 。因此 , 解析解的重要意义是不容置疑的。1957 年 , Stoker 将坝址流态分为连续波流 、临界 [5 ] 流和不连续流三个流态 , 推导了矩形河谷和下游有水但起始流速为零情况的瞬间全溃坝址峰顶流量公式 。 1982 年 , 谢任之在吸收前人经验的基础上 , 综合连续波与间断波的解法 , 用抛物线概化河谷断面 , 去掉对下 游水深和流速的控制条件 , 推出了 “统一公式”, 可用于各种情况的瞬间全溃的坝址峰顶流量计算 , 并给出了 便于查用的表格[6 ]