交错网格下的浅水方程高分辨Gauss型格式.pdfVIP

第 13 卷 第 4 期 水 科 学 进 展 Vol13 ,No4 2002 年 7 月 　　 ADVANCES IN WATER SCIENCE 　　Jul . ,2002 　 交错网格下的浅水方程高分辨 Gauss 型格式 1 2 ,3 2 2 王如云 , 邱建贤 , 戴嘉尊 , 赵 　宁 ( 1 河海大学交通与海洋工程学院 , 江苏 南京 　210098 ; 2 南京航空航天大学 , 江苏 南京 　210016 ; 3 中国科学技术大学数学系 , 安徽 合肥 　230026) 摘要 : 在交错网格上 , 基于高精度的 Gauss 积分公式 , 针对浅水波方程设计了对模拟涌波具有高 分辨率的完全二阶精度的数值计算格式 。由于采用了交错网格 , 差分格式不需要解 Riemann 问 题 , 因此本文格式具有计算简单 、工作量少 、编程简便等特点 。另外 , 在一维单个方程时 , 本 文格式在 CFL (Courant Friedrich Lewy) 条件限制下为 TVD ( Total Variation Diminishing) 格式 , 在二维和 ( ) 三维情况下格式具有 MmB Maximum and Minimum Bounds Preserving 性质 。利用国家高性能计算中 ( ) 心 合肥 的曙光1000 型分布存储大规模并行机 ,对在交错网格下所构造的求解浅水方程的高分 辨差分格式进行了并行实现 , 几个算例的计算结果令人满意 。 关 　键 　词 : 交错网格 ; 浅水方程 ; 差分格式 ; Gauss 积分式 ( ) 中图分类号 : TV 1314 　　　文献标识码 : A 　　　文章编号 2002 0440306 溃坝波 、涌潮及滑坡涌浪等现象 , 虽然产生的物理背景不同 , 但它们的波形都有一个明显 ( ) [ 1] 的特点 , 即波前是一间段面 断波 ,此类波一般统称为涌波 。由于涌波波前的间断性 , 使得 其更具破坏性 。因此 , 对涌波的形成和传播过程的研究 , 具有重大的现实意义 。对于涌波的研 究 , 一般采用 Airy 在 1845 年假定水波的波长无限长 , 同时波高是有限值导出的浅水波方程 。 这不仅是因为涌波具有波长无限、波高有限的特点 , 而且是因为浅水波方程属双曲型偏微分方 程 , 本身具有间断解的特点[2 ] 。 关于双曲守恒型偏微分方程的数值求解 , 基于设计满足守恒性 、熵条件等要求的格式 , 使 得激波捕捉方法中的误差可局限在准确物理间断面的附近 ,