Efficient well-balanced hydrostatic upwind schemes for shallow-water equations

The proposed work concerns the numerical approximations of the shallow-water equations with varying topography. The main objective is to introduce an easy and systematic technique to enforce the well-balance property and to make the scheme able to deal with dry areas. To access such an issue, the derived numerical method is obtained by involving the free surface instead of the water height and this produces the scheme well-balanced independently from the numerical flux function associated with the homogeneous problem. As a consequence, we obtain an easy well-balanced scheme which preserves non-negative water height. When compared with the well-known hydrostatic reconstruction, the presented topography discretization does not involve any max function known to introduce some numerical errors as soon as the topography admits very strong variations or discontinuities. A second-order MUSCL accurate reconstruction is adopted. The proposed hydrostatic upwind scheme is next extended for considering 2D simulations performed over unstructured meshes. Several 1D and 2D numerical experiments are performed to exhibit the relevance of the scheme.

[1]  Yulong Xing,et al.  High order finite difference WENO schemes with the exact conservation property for the shallow water equations , 2005 .

[2]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[3]  Marica Pelanti,et al.  Approximation of the hydrostatic Navier-Stokes system for density stratified flows by a multilayer model: Kinetic interpretation and numerical solution , 2011, J. Comput. Phys..

[4]  M. Thanh,et al.  The Riemann problem for the shallow water equations with discontinuous topography , 2007, 0712.3778.

[5]  Jostein R. Natvig,et al.  Well-balanced finite volume schemes of arbitrary order of accuracy for shallow water flows , 2006, J. Comput. Phys..

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

[7]  W. Thacker.  Some exact solutions to the nonlinear shallow-water wave equations , 1981, Journal of Fluid Mechanics.

[8]  N. Gouta,et al.  A finite volume solver for 1D shallow‐water equations applied to an actual river , 2002 .

[9]  P. Raviart,et al.  Numerical Approximation of Hyperbolic Systems of Conservation Laws , 1996, Applied Mathematical Sciences.

[10]  Frédéric Lagoutière,et al.  Stability of reconstruction schemes for scalar hyperbolic conservations laws , 2008 .

[11]  Olivier Delestre,et al.  SWASHES: a compilation of shallow water analytic solutions for hydraulic and environmental studies , 2011, 1110.0288.

[12]  Olivier Delestre.  Simulation du ruissellement d'eau de pluie sur des surfaces agricoles. (Rain water overland flow on agricultural fields simulation) , 2010 .

[13]  Christophe Berthon,et al.  Robustness of MUSCL schemes for 2D unstructured meshes , 2006, J. Comput. Phys..

[14]  Shi Jin,et al.  Two Interface-Type Numerical Methods for Computing Hyperbolic Systems with Geometrical Source Terms Having Concentrations , 2005, SIAM J. Sci. Comput..

[15]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

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

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

[18]  Philippe G. LeFloch,et al.  A Godunov-type method for the shallow water equations with discontinuous topography in the resonant regime , 2011, J. Comput. Phys..

[19]  Shi Jin,et al.  An efficient method for computing hyperbolic systems with geometrical source terms having concentrations ∗ , 2004 .

[20]  C. Berthon,et al.  Stability of the MUSCL Schemes for the Euler Equations , 2005 .

[21]  Manuel J. Castro,et al.  WELL-BALANCED NUMERICAL SCHEMES BASED ON A GENERALIZED HYDROSTATIC RECONSTRUCTION TECHNIQUE , 2007 .

[22]  Fabien Marche,et al.  A Positive Preserving High Order VFRoe Scheme for Shallow Water Equations: A Class of Relaxation Schemes , 2008, SIAM J. Sci. Comput..

[23]  Christophe Chalons,et al.  Relaxation approximation of the Euler equations , 2008 .

[24]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

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

[26]  E. Toro,et al.  Restoration of the contact surface in the HLL-Riemann solver , 1994 .

[27]  Carlos Parés,et al.  A Q-SCHEME FOR A CLASS OF SYSTEMS OF COUPLED CONSERVATION LAWS WITH SOURCE TERM. APPLICATION TO A TWO-LAYER 1-D SHALLOW WATER SYSTEM , 2001 .

[28]  Jean-Marc Hérard,et al.  On the use of symmetrizing variables for vacuums , 2003 .

[29]  B. Larrouturou.  How to preserve the mass fractions positivity when computing compressible multi-component flows , 1991 .

[30]  Emmanuel Audusse,et al.  A Fast and Stable Well-Balanced Scheme with Hydrostatic Reconstruction for Shallow Water Flows , 2004, SIAM J. Sci. Comput..

[31]  L. Gosse.  A well-balanced flux-vector splitting scheme designed for hyperbolic systems of conservation laws with source terms☆ , 2000 .

[32]  Carlos Parés,et al.  On the well-balance property of Roe?s method for nonconservative hyperbolic systems , 2004 .

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

[34]  Steve Bryson,et al.  Well-balanced positivity preserving central-upwind scheme on triangular grids for the Saint-Venant system , 2011 .

[35]  Jean-Marc Hérard,et al.  Un schma simple pour les quations de Saint-Venant , 1998 .

[36]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .

[37]  T. Gallouët,et al.  Some approximate Godunov schemes to compute shallow-water equations with topography , 2003 .

[38]  T. Morales de Luna,et al.  On a shallow water model for the simulation of turbidity currents , 2009 .

[39]  F. Marche.  A Simple Well-Balanced Model for Two-Dimensional Coastal Engineering Applications , 2008 .

[40]  Thierry Gallouët,et al.  On a rough Godunov scheme , 1996 .

[41]  Christophe Berthon,et al.  Why the MUSCL–Hancock Scheme is L1-stable , 2006, Numerische Mathematik.

[42]  Shi Jin,et al.  AN EFFICIENT METHOD FOR COMPUTING HYPERBOLIC SYSTEMS WITH GEOMETRICAL SOURCE TERMS HAVING CONCENTRATIONS ∗1) , 2004 .

[43]  Thierry BuÄard.  A sequel to a rough Godunov scheme: application to real gases , 2000 .

[44]  Pilar García-Navarro,et al.  Efficient construction of high‐resolution TVD conservative schemes for equations with source terms: application to shallow water flows , 2001 .

[45]  F. Bouchut.  Nonlinear Stability of Finite Volume Methods for Hyperbolic Conservation Laws: and Well-Balanced Schemes for Sources , 2005 .

[46]  Jean-Marc Hérard,et al.  Some recent finite volume schemes to compute Euler equations using real gas EOS , 2002 .

[47]  Giovanni Russo,et al.  Central schemes for conservation laws with application to shallow water equations , 2005 .

[48]  J. Greenberg,et al.  Analysis and Approximation of Conservation Laws with Source Terms , 1997 .

[49]  Chi-Wang Shu,et al.  On positivity preserving finite volume schemes for Euler equations , 1996 .

[50]  B. Perthame,et al.  A kinetic scheme for the Saint-Venant system¶with a source term , 2001 .

[51]  B. Perthame,et al.  Relaxation of Energy and Approximate Riemann Solvers for General Pressure Laws in Fluid Dynamics , 1998 .

[52]  Randall J. LeVeque,et al.  A class of approximate Riemann solvers and their relation to relaxation schemes , 2001 .

[53]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[54]  Jean-Frédéric Gerbeau,et al.  Derivation of viscous Saint-Venant system for laminar shallow water , 2001 .