一维浅水流动方程的godunov 格式求解 - 水科学进展.pdfVIP

第 14 卷 第 4 期 水 科 学 进 展 Vol14 ,No4 2003 年 7 月 　　 ADVANCES IN WATER SCIENCE 　 Jul , 2003 　 一维浅水流动方程的 Godunov 格式求解 1 ,2 2 2 潘存鸿 , 林炳尧 , 毛献忠 ( 1 上海大学理学院 , 上海 　200436 ; 2 浙江省水利河口研究院 , 浙江 杭州 　310020) ( 摘要 : 以准确 Riemann 解为基础 , 建立了求解一维非平底浅水流动方程的 Godunov 格式 , 用“水位方程法 Water Level Formulation , WLF) ”求解 Riemann 解 ,结合中心差分和 Riemann 解离散底坡项 , 保证了计算格式的和谐性 。经算例验 证 , 方法健全 、通用 , 且分辨率高 。 关 　键 　词 : 一维浅水流动方程 ; Godunov 格式 ; Riemann 解 ; 源项 ( ) 中图分类号 : TV 1312 　　　文献标识码 : A 　　　文章编号 : 2003 0443007 自从 1959 年 Godunov 在博士论文中提出了利用 Riemann 解求解双曲型方程的格式以来 , 40 余年中 , 已经不 断改进 。Van Leer 等发展了二阶精度 Godunov 格式 , 使得该格式更有吸引力 。因为 Godunov 格式具备模拟大梯 度流动和自动捕捉激波的能力 , 在计算流体力学中得到了广泛的应用[ 1] 。 1981 年 , Marshall 和 Mendez[2 ] 应用 Riemann 解求解浅水流动方程 , 开辟了此方程新的求解途径 。最近 10 年 来 , 以Riemann 解为基础的 Godunov 型格式求解浅水流动方程的研究成果逐渐增多 , 从一维发展到二维 , 成为 求解浅水大梯度流动的主要方法[3 ] 。 如果不考虑底坡和摩阻作用 , 则无源项的浅水流动方程类似于空气动力学中的欧拉方程 , 早期工作大多针 [3～12 ,20 ] [9 ] 对齐次浅水流动方程 , 直接移植空气动力学中成熟的 Godunov 格式 。在此基础上 , 许为厚等 应用欧拉 拉格朗日统一坐标系准确模拟了带滑移线的浅水间断流动 。 事实上 , 浅水流动模拟的关键是河床起伏的模拟 。应用常规的分步法求解非平底浅水方程时 , 计算结果不 具和谐性 , 即在静水条件下 , 流速不为零 , 且水位不再为常数 。近几年已有学者致力于非平底浅水流动方程和