二维溃坝波遇障碍物水流泥沙数值模拟99808.pdfVIP

水 　　利 　　学 　　报 　2005 年 5 月 SHUILI 　　XUEBAO 第 36 卷 　第 5 期 ( ) 文章编号 2005 二维溃坝波遇障碍物的水流泥沙数值模拟 白玉川,许栋 ,王玉琦 ,张梅亭 (天津大学 建筑工程学院 ,天津　300072) 摘要 :本文通过对浅水方程进行 Godunov 差分离散和对泥沙运动方程进行一阶迎风差分离散 ,建立了平面二维溃坝 水流泥沙数值模型。在算例中计算了溃坝波遇障碍物后的水流反射、绕射及泥沙冲淤 ,计算结果表明 ,溃坝波遇障 碍物后会产生明显的反射和绕射 ,并在障碍物前形成壅水、障碍物的侧面会发生剧烈的冲刷 ,障碍物前后则会发生 局部淤积。本模型能够较好地模拟溃坝间断水流的运动 ,结果符合水流泥沙运动的基本规律。 关键词 :溃坝 ;水流 ;泥沙 ; Godunov 格式 中图分类号:TV1392 文献标识码 :A 水库堤防的安全是水工建筑物设计和管理的核心问题。长久以来 ,水利科学工作者建立了许多实 用的溃坝水流数值计算模型 ,如 Godunov 差分格式、TVD 差分格式、MacCormack 差分格式、有限元模型 等 ,能够对一维、二维甚至三维的溃坝水流进行较好地模拟[1～6 ] 。然而 ,实际的溃坝问题是一个水流泥 沙相互作用的复杂过程 ,一方面 ,溃坝水流由于流速较大 ,会强烈地冲刷河床、河岸及口门 ,并在障碍物 附近及下游等处形成泥沙淤积;另一方面 ,河床、河岸及口门的剧烈变形也在不断地改变水流的特性 ,使 其向着输沙平衡的方向发展。因此 ,完全忽略河床变形的溃坝水流计算可能带来较大的计算误差[7 ,8 ] ; 此外 ,研究溃坝水流作用下的泥沙运动有利于更全面地考虑溃坝对下游河道及建筑物的影响 ,也可以促 ( ) 进对溃坝机理和溃坝过程的进一步研究 尤其是土坝 ,具有一定的实际意义。 本文在利用 Godunov 差分格式进行平面二维溃坝水流计算的基础上 ,加入了非恒定流输沙的计算 , 综合考虑了推移质和悬移质作用下的河床变形 ,并在计算过程中 ,不断根据上一时段泥沙计算的结果 , 调整水流计算中的水深 , 以反映河床变形对溃坝水流的影响。 1 　溃坝水流计算 11 　水流控制方程 　溃坝水流的控制方程为浅水方程 ,在平面二维计算中 ,其守恒型微分形式为 h hu hv + + = 0 t x y hu 2 1 2 huv + hu + g h + = F ( 1) t x 2 y hv huv 2 1 2 + + hv + g h = G t x y 2 h 式中 : h 为水深 ; u 、v 为水流速度 ; g 为重力加速度 ; F 、G