320期应用力学学报Vol.32No.02015年00月.PDFVIP

第 卷 第 期 应 用 力 学 学 报 32 0 Vol.32 No.0 年 月 2015 00 CHINESEJOURNAL OF APPLIED MECHANICS Xxx. 2015 基金项目：三峡大学人才科研启动基金 （KJ2014B031）；国家自然科学基金 收稿日期：2016-11-22 修回日期：2017-01-23 第一作者简介：吴泽艳，男，1976 年生，博士，三峡大学，讲师；研究方向——计算流体力学。 通讯作者：杨忠勇，E-mail: ayong0710@163.com 第0 期 吴泽艳，等：基于简单WENO-RKDG 的溃坝问题数值模拟 2 第0 期 吴泽艳，等：基于简单WENO-RKDG 的溃坝问题数值模拟 3 文章编号：1000- 4939(2015) 00-0000-00 基于简单WENO-RKDG 的溃坝问题数值模拟 吴泽艳 杨忠勇 (三峡大学水利与环境学院 443002 宜昌) 摘要：对溃坝问题水流间断面的高精度、高分辨率数值模拟是水动力学的重要内容。简单加权本 质无振荡 （WENO ）限制器由 “问题单元”及其相邻单元的解重构 “问题单元”的解，从而抑制 数值解的非物理振荡，能够很好地模拟间断问题。本文详细介绍了简单WENO 限制器的基本原理 和过程。将简单WENO 限制器-Runge-Kutta 间断Galerkin 方法应用于二维浅水控制方程的求解中， 对二维矩形明渠中大坝瞬间全溃、局部溃塌所致的水流运动进行了数值模拟，并将数值计算结果 与理论分析进行了比较，证实了方法的可行性。计算结果表明简单WENO 限制器-RKDG 方法能够 很好地模拟溃坝波的演进过程。 关键词：溃坝; 浅水波方程；Runge-Kutta 间断Galerkin；简单WENO 限制器；数值模拟 中图分类号：TV131.4；O241.82 文献标识码：A 面运动的VOF方程进行离散建立了溃坝水流运动的 1 引 言 三维数学模型并对急弯河道进行了数值模拟计算； [7] [8] [9] 袁美等 、夏军强等 、宋利祥等 对浅水波方程, 大坝失事或毁坏会对下游地区人民生命财产安 在二维非结构网格上给出了一种非振荡有限体积方 全造成灾难性的破坏，并对大坝下游的生态环境造 法，对二维溃坝问题进行了数值模拟；张春泽等[10] 成不良影响，因此对溃坝洪水进行数值模拟与分析 运用计算二维明渠非恒定流的格子Boltzmann模 具有重要的理论与现实意义。多年来，国内外水利 型 ，利用格子 Boltzmann 算法耦合度低和 工程师和研究工作者做了很多工作。 CUDA(ComputeUnified DeviceArchitecture)线程 早期对溃坝问题的高精度、高分辨率数值模拟 级并行的优势，在GPU(Graphic Processing Unit) [2,3] [11] 方法综述可见文献[1]。近期，魏文礼等 先后采 上实现了该模型