复杂地形条件下weno-roe方法在浅水模拟中的改进-应用力学学报.pdfVIP

第28 卷 第3 期 应 用 力 学 学 报 Vol.28 No.3 2011 年6 月 CHINESE JOURNAL OF APPLIED MECHANICS Jun. 2011 文章编号：1000-4939(2011) 03-0249-05 复杂地形条件下WENO-Roe 方法 在浅水模拟中的改进 1 2 1 1 王党伟 余明辉 陈建国 吉祖稳 ( 中国水利水电科学研究院 100048 北京) 1 (武汉大学水资源与水电工程科学国家重点实验室 430072 武汉) 2 摘要：在分析浅水方程与二维空气运动方程差异的基础上，在不改变原有浅水方程形式的前提下， 提出了“局部水位法” 离散连续方程，并针对动量方程中底坡项提出了更具普遍意义的离散方法， 保证了浅水方程离散后的平衡性。通过不规则地形下潮流以及混合流流动的模拟，得到的计算结 果符合物理实际，与精确解符合良好，最大相对误差不超过 4% ，验证了此方法在复杂地形上的 平衡性，同时本方法又具有良好的间断捕捉能力与稳定性。 关键词：WENO-Roe 方法；浅水方程；复杂地形；局部水位法 中图分类号：TV131.2 文献标识码：A 出现的动量方程中底坡项离散的不平衡性进行了细 1 引 言 致的研究。目前尚缺乏对于直接采用Roe 方法离散 连续方程可能出现不合理结果的研究。本文从水流 传统的数值方法模拟包含激波的水流运动如 和空气流动的物理机制出发，通过两者的对比，分 溃坝、溃堤水流时，结果往往会出现非物理震荡或 析了Roe 方法直接用于模拟浅水运动而导致不合理 波形严重坦化失真等现象。鉴于平面二维水流与气 结果的根源，提出了有效的解决方案。 体运动方程方程的相似性，可以将空气动力学中用 于捕捉激波的数值方法移植到浅水模拟中，文献[1] 2 控制方程与总体离散方法 较早进行了这方面的尝试，采用 TVD 格式模拟溃 坝波并取得了成功。由于 TVD 格式的局限性，其 基本方程采用守恒形式的控制方程。 U F G 精度较低，WENO 插值方法提供了一条较为满意的 ∂ ∂ ∂ (1) + + S (U) t x y 解决途径，通过其特定的插值方式可以将原有的数 ∂ ∂ ∂ 值方法精度大大提高[2] 。文献[3]采用WENO 方法插 其中 T S [0 gh(s =− f ) gh(s − f )] ， 值界面两侧变量值，配合Roe 的FDS 格式计算界面 x x y y