用守恒元和解元法数值模拟二维溃坝洪水波.pdfVIP

水 利 学 报 2005 年 10 月 SHUILI XUEBAO 第 36 卷 第 10 期 文章编号: 0559-9350(2005)10-1224-06 用守恒元和解元法数值模拟二维溃坝洪水波 张永祥，陈景秋 (重庆大学 工程力学系，重庆 400030) 摘要: 本文采用守恒元和解元方法(简称CE/SE法)来求解二维浅水方程，该方法从守恒积分型方程出发，严格遵 守物理意义上的守恒律，通过设立守恒元和解元，从而提出一种新的数值模拟方法。用该方法计算得到的全溃的 溃坝洪水波的数值解和理论解与文献中的实验值符合的较好。在此基础上，进一步对二维局部溃坝洪水波的演进 过程进行模拟，结果表明，CE/SE法精度较高，稳定性较好，其溃坝洪水波的间断具有很高的分辨率，能很好的模 拟溃坝洪水波的运动过程。 关键词: 反拱水垫塘；拱座；整体稳定性 中图分类号: TV131.4 文献标识码:A 1 研究背景 大坝溃决后的洪水波会对滞洪区和下游造成很大危害。溃坝所造成的洪水波是前进长波，非线性特 征较强，对于其演进过程的研究，可为指导防汛抗洪提供依据[1]。在大坝全溃的情况下，溃坝洪水波的 流动特性和状态通常可以按一维流动来模拟，但是大坝溃决通常以局部溃坝为主，这时需按二维流动来 模拟。在数学表达上，这样的流动可以用二维圣维南浅水方程来描述，该方程是非线性的双曲型偏微分 方程[2] [3] ，只有特定的条件下，其一维情况才有理论解，而一般情况下理论解很难得到 ，因此，需要采 用数值方法来求解浅水方程。通常有3种基本方法来计算这种带有强烈的间断解的方程，即特征线法、 激波拟合法和激波捕捉法。在浅水方程的数值计算中，激波捕捉法是目前主要使用的方法。 [4] Chang S C 提出一种新的数值方法——时空守恒元和解元方法(简称CE/SE法)，该方法无论从概念 上还是从构造方法上都与传统的数值方法(有限元法、有限体积法、有限差分法等)不同。首先，其构造 思想简单、物理意义清晰、格式通用性好，它不用求解黎曼问题，也不用求解特征问题；其次，它把时 间与空间完全统一同等对待，并从守恒积分型方程出发，通过设立守恒元和解元，使局部和全局都严格 保证物理意义上的守恒律；第三，该方法把流场变量及其空间导数均作为变量，同时求解，这样与传统 的差分格式相比，在相同的网格点数的情况下，其格式精度可以达到更高，同时更便于精确地满足边界 条件。文献[5～8]对此方法进行了改进和优化，使此方法更利于离散与求解，同时易于编程。本文应用改 进的CE/SE法，通过自行开发程序对二维浅水方程进行离散求解，成功展示了溃坝洪水波的运动特性。 2 浅水方程的守恒元和解元的格式 二维圣维南浅水方程的守恒形式为 收稿日期: 2005-07-26 基金项目：国家自然科学基金资助项目；教育部重点科技项目(02058) 作者简介：张永祥(1978-)，浙江诸暨人，博士生，主要从事计算流体力学研究。E-mail:flasit@126.com 1 水 利 学 报 2005 年 10 月 SHUILI XUEBAO 第 36 卷 第 10 期 ∂qm ∂f m ∂g m + + Sm , m 1,2,3