ENO and WENO Schemes with the Exact Conservation Property for One-Dimensional Shallow Water Equations

In this paper we propose new finite difference numerical schemes for hyperbolic conservation law systems with geometrical source terms. In the development of the new schemes we use the essentially nonoscillatory (ENO) and weighted ENO (WENO) reconstruction, developed by Harten, Osher, Engquist, Chakravarthy, Shu, and Jiang, and the idea of the balancing between the flux gradient and the source term, introduced by Bermudez and Vazquez. Actually, the new schemes are ENO and WENO schemes with the source term decomposed, i.e., the ENO and WENO reconstruction is applied not only to the flux but to a combination of the flux and the source term. In particular, when new schemes are applied to the shallow water equations the new schemes verify the exact conservation property (C-property). We present the algorithm, the proof of the exact C-property, and numerical results for several test problems.

[1]  Randall J. LeVeque,et al.  Wave Propagation Methods for Conservation Laws with Source Terms , 1999 .

[2]  Guang-Shan Jiang,et al.  A High-Order WENO Finite Difference Scheme for the Equations of Ideal Magnetohydrodynamics , 1999 .

[3]  M. Vázquez-Cendón Improved Treatment of Source Terms in Upwind Schemes for the Shallow Water Equations in Channels with Irregular Geometry , 1999 .

[4]  D. Causon,et al.  The surface gradient method for the treatment of source terms in the shallow-water equations , 2001 .

[5]  Chi-Wang Shu,et al.  Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy , 2000 .

[6]  Z. Xin,et al.  The relaxation schemes for systems of conservation laws in arbitrary space dimensions , 1995 .

[7]  Philippe Montarnal,et al.  Real Gas Computation Using an Energy Relaxation Method and High-Order WENO Schemes , 1999 .

[8]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[9]  S. Osher,et al.  Uniformly High-Order Accurate Nonoscillatory Schemes. I , 1987 .

[10]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[11]  L. Margolin,et al.  MPDATA: A Finite-Difference Solver for Geophysical Flows , 1998 .

[12]  Patrick Jenny,et al.  Rankine-Hugoniot-Riemann Solver Considering Source Terms and Multidimensional Effects , 1998 .

[13]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[14]  Shi Jin,et al.  A steady-state capturing method for hyperbolic systems with geometrical source terms , 2001 .

[15]  Jean-Antoine Désidéri,et al.  Upwind schemes for the two-dimensional shallow water equations with variable depth using unstructured meshes , 1998 .

[16]  J. Greenberg,et al.  A well-balanced scheme for the numerical processing of source terms in hyperbolic equations , 1996 .

[17]  Alfredo Bermúdez,et al.  Upwind methods for hyperbolic conservation laws with source terms , 1994 .

[18]  Pilar García-Navarro,et al.  Flux difference splitting and the balancing of source terms and flux gradients , 2000 .