溃坝水流模拟方法的对比分析-水利水运工程学报.pdfVIP

第3期 水 利 水 运 工 程 学 报 No.3 摇 2012年6月 HYDRO鄄SCIENCE AND ENGINEERING Jun. 2012 溃坝水流模拟方法的对比分析 1 1 2 肖摇 潇 , 吴时强 , 程永舟 (1. 南京水利科学研究院 水文水资源与水利工程科学国家重点实验室, 江苏 南京摇 210029;2. 长沙理工大 学, 湖南 长沙摇 410114) 摘要:采用FLUENT软件与光滑粒子流体动力学法(SPH)分别对溃坝水流进行数值模拟,对比分析溃坝水流 的流动形态、流速及水位变化.结果表明,两者较为一致地反映了溃坝水流的流动过程及流体的速度场,其中 SPH方法能精确地捕捉自由液面,所得自由水面附近结果更为准确;FLUENT所计算的流速、水位等数据更加精 确,能很好地体现水流的表面张力.分析两种方法各自的优缺点和适用性,使其在面对复杂的工程实际问题时 发挥更大的作用,促进溃坝洪水数学模型的应用与发展. 关摇 键摇 词:溃坝;SPH; 自由面;VOF模型 中图分类号:TV131摇 摇 摇 摇 文献标志码:A摇 摇 摇 摇 文章编号:1009-640X(2012)03-0037-07 溃坝洪水是指由于大坝等挡水物体突然溃决,发生水体突泄而形成的洪水.溃坝洪水波洪峰流量大、传 播速度快,破坏力极强,一旦发生将对下游人民生命财产安全带来巨大威胁.基于溃坝后果的严重性,溃坝洪 水波演进模拟一直是众多学者所关注的研究重点. [1] [2] J郾C郾Martin等 对在重力作用下直立方形水柱塌落水波进行了试验研究;P郾K郾Stansby等 试验研究了 [3] 前后具有不同高程比的两块水体溃决过程;M郾A郾Cruchaga等 对具有不同高程及不同密度的水体塌落过程 进行了试验分析研究.这些试验数据为数学模型的验证提供了依据.随着数值计算方法不断进步,数值模拟 手段越来越多地被应用于溃坝洪水波传播问题的研究中.常见的数值方法有有限元法、有限差分法、有限体 积法及近几十年发展起来的无网格方法,如粒子法、SPH法等.但这些数值方法由于离散方式、计算方法不 同,从而对模拟结果的影响也不同,因此,研究不同数值方法对溃坝洪水波演进模拟适用条件显得尤为重要. 计算流体力学软件FLUENT采用基于有限体积法对Navier鄄Stokes控制方程进行离散求解,并采用标准 [4] 紊流模型对方程进行封闭,自由面采用VOF方法处理.宋新远等 利用FLUENT模拟了二维滑坡涌浪;刑爱 [5] 国等 利用FLUENT分析了堰塞湖溃坝洪水演进.光滑粒子流体动力学法(SPH)是一种基于拉格朗日形式 的无网格法,通过粒子逼近得出SPH形式的Navier鄄Stokes方程,从而建立相应的数学模型.J郾J郾Monaghan [6] [7] 等 通过与SPH相类似的方法生成孤立波,且通过相关试验来研究缓坡上波的破碎现象.E郾Y郾M郾Lo等 将SPH法与大涡模相结合拟研究溃坝洪水波传播问题,并深入研究了孤立波的爬坡与破碎现象. [8] A郾Colagrossit等 详细研究溃坝水流与直立墙作用形成的复杂界面问题. 考虑到FLUENT软件是应用有限体积法的基于网格的模拟软件,而SPH方法是提出最早且发展比较完 善的无网格法,本文采用这两种方法分别模拟溃坝水流问题,比较两种方法模拟的溃坝水流流态及水力特 性,分析两种方法的优缺点,探讨两种方法模拟溃坝水流问题