High order finite difference WENO schemes with the exact conservation property for the shallow water equations

Shallow water equations with nonflat bottom have steady state solutions in which the flux gradients are nonzero but exactly balanced by the source term. It is a challenge to design genuinely high order accurate numerical schemes which preserve exactly these steady state solutions. In this paper we design high order finite difference WENO schemes to this system with such exact conservation property (C-property) and at the same time maintaining genuine high order accuracy. Extensive one and two dimensional simulations are performed to verify high order accuracy, the exact C-property, and good resolution for smooth and discontinuous solutions.

[1]  J LeVequeRandall.  Balancing Source Terms and Flux Gradients in High-Resolution Godunov Methods , 1998 .

[2]  Enrique D. Fernández Nieto,et al.  A family of stable numerical solvers for the shallow water equations with source terms , 2003 .

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

[4]  Randall J. LeVeque,et al.  Balancing Source Terms and Flux Gradients in High-Resolution Godunov Methods , 1998 .

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

[6]  Chi-Wang Shu.  Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[7]  Luka Sopta,et al.  WENO schemes for balance laws with spatially varying flux , 2004 .

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

[9]  Doron Levy,et al.  CENTRAL-UPWIND SCHEMES FOR THE SAINT-VENANT SYSTEM , 2002 .

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

[11]  P. Lax,et al.  On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws , 1983 .

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

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

[14]  Francis Filbet,et al.  Approximation of Hyperbolic Models for Chemosensitive Movement , 2005, SIAM J. Sci. Comput..

[15]  ShuChi-Wang,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes, II , 1989 .

[16]  Luka Sopta,et al.  ENO and WENO Schemes with the Exact Conservation Property for One-Dimensional Shallow Water Equations , 2002 .

[17]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[18]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

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

[20]  Giovanni Russo,et al.  Central Schemes for Balance Laws , 2002 .

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

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

[23]  Kun Xu,et al.  A well-balanced gas-kinetic scheme for the shallow-water equations with source terms , 2002 .

[24]  Luka Sopta,et al.  Balanced finite volume WENO and central WENO schemes for the shallow water and the open-channel flow equations , 2004 .