求解二维浅水方程的一种高分辨率有限体积法.pdfVIP

维普资讯 2006年 12月 应用数学与计算数学学报 第 20卷 第2期 Dec．，2006 COMM ．ONAPPL．MATH．AND COMPUT V01．20 No．2 求解二维浅水方程的一种高分辨率有限体积法 窦红 汪继文 1．安徽大学数学与计算科学学院，合肥，230039 2．安徽大学计算智能与信号处理教育部重点实验室，合肥，230039 摘要。本文将一种vanAlbada型可微的限制器函数引入到二维浅水方程的求解 中， 发展了一种求解二维浅水方程的有 限体积法．数值实验结果表明，该方法不仅计算精度 高，而且较其它求解二维浅水方程的高精度有限体积法，在数值解的收敛性能方面大有改 善． 关键词t浅水方程，有限体积法，限制器函数 1．引 言 浅水方程的求解是计算数学、计算流体力学的一个重要研究内容．近年来，求解浅 水方程的有限体积法 (FVM)越来越受到人们的重视．FVM具有如下一些优点 【】： FVM从控制方程的守恒形式出发，利用 Riemann解算子，可以准确地捕捉到激波； FVM可以采用无结构阿格，比有限差分法更能适应复杂的求解区域；FVM的计算比 有限元方法的计算代价低，更容易实现． 在二阶和二阶以上精度的FVM中，为了消除数值振荡，都要在变量的重构过程中 引入限制器 (Limiter)函数，近十多年来，不少学者设计了各种各样的限制器，在一维 问题中，几个著名的限制器有 van—Leer，Superbee，Minmod，vanAlbada等限制器 [1】_在 二维问题中，Barth和Jesperson[】给出的限制器及其不同的变形被广泛使用．限制器 的使用能有效地消除虚假振荡，但也同时产生了两个负面问题，一是会阻止数值解收 敛到稳定状态 I引．对此，我们在求解浅水方程的数值计算中也深有体会；二是会降低 光滑极值点附近解的精度．这些现象在多维问题中更加严重．为消除限制器的这些缺 点，人们进行了不少努力，但Jawahar和Kamath4[】认为结果都不尽如人意，为此，[41 中给出了一种 VailAlbada型可微的限制器函数，发展了一种求解 Euler方程的高分辨 率FVM．4【】中的实验结果表明，该限制器的使用取得了理想的效果． 和 Aliabadi[51 将这种限制器与间断Galerkin方法相结合，发展了一种求解空气动力学问题的高分辨 率方法．但使用这种限制器的求解浅水方程的高分辨率方法尚未见报道 本文中，我 们也采用 4【]中的限制器，发展了一种求解二维浅水方程的FVM．基lF求解二维浅水 方程的斜水跃问题，我们对这种新发展的FVM与已有的有代表性的高精度 FVM在 精度与收敛性等方面进行了比较．数值实验结果表明，本文方法不仅消除了振荡，计 本文 2006年 3月6日收到． 安徽省高校 自然科学研究项 目(2006KJ028B) 维普资讯 应用数学与计算数学学报 20卷 算精度高，而且有效地解决了目前求解浅水方程常用的高精度 FVM数值解不收敛或 难以收敛到稳定状态的问题． 2．控制方程与数值方法 水波运动是一种在重力作用下有 自由面的水运动．水一般视为不可压流体， 浅水 波理论又假定水的速度沿深度方向不变．在假定 自由面压力分布为常数， 并忽略表面 风切应力，Coriolis力和内部水平切应力时，描写二维瞬变 (transient)自由面流体运 +V ·F=S， (1) 这里 U是守恒变量向量，S是源项，包括