A well-balanced positivity preserving second-order scheme for shallow water flows on unstructured meshes

We consider the solution of the Saint-Venant equations with topographic source terms on 2D unstructured meshes by a finite volume approach. We first present a stable and positivity preserving homogeneous solver issued from a kinetic representation of the hyperbolic conservation laws system. This water depth positivity property is important when dealing with wet-dry interfaces. Then, we introduce a local hydrostatic reconstruction that preserves the positivity properties of the homogeneous solver and leads to a well-balanced scheme satisfying the steady-state condition of still water. Finally, a formally second-order extension based on limited reconstructed values on both sides of each interface and on an enriched interpretation of the source terms satisfies the same properties and gives a noticeable accuracy improvement. Numerical examples on academic and real problems are presented.

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

[2]  Benoît Perthame,et al.  Maximum principle on the entropy and second-order kinetic schemes , 1994 .

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

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

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

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

[7]  Rémi Abgrall,et al.  On essentially non-oscillatory schemes on unstructured meshes: analysis and implementation , 1994 .

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

[9]  Laurent Gosse,et al.  A Well-Balanced Scheme Using Non-Conservative Products Designed for Hyperbolic Systems of Conservati , 2001 .

[10]  Ramaz Botchorishvili,et al.  Equilibrium schemes for scalar conservation laws with stiff sources , 2003, Math. Comput..

[11]  P. Souganidis,et al.  Existence and stability of entropy solutions for the hyperbolic systems of isentropic gas dynamics in Eulerian and Lagrangian coordinates , 1998 .

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

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

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

[15]  Alexander Kurganov,et al.  On a hybrid finite-volume-particle method , 2004 .

[16]  B. Perthame Second-order Boltzmann schemes for compressible Euler equations in one and two space dimensions , 1992 .

[17]  Marie-Odile Bristeau,et al.  Transport of pollutant in shallow water using kinetic schemes , 2001 .

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

[19]  E. Audusse,et al.  A multilayer Saint-Venant model: Derivation and numerical validation , 2005 .

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

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

[22]  F. Saleri,et al.  A new two-dimensional shallow water model including pressure effects and slow varying bottom topography , 2004 .

[23]  Chaowei Hu,et al.  No . 98-32 Weighted Essentially Non-Oscillatory Schemes on Triangular Meshes , 1998 .

[24]  B. Perthame,et al.  A new model of Saint Venant and Savage–Hutter type for gravity driven shallow water flows , 2003 .

[25]  J. Hervouet A high resolution 2-D dam-break model using parallelization , 2000 .

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

[27]  Enrique D. Fernández Nieto,et al.  An Entropy-correction free solver for non-homogeneous Shallow Water equations , 2003 .

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

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

[30]  Emmanuel Audusse,et al.  Transport of Pollutant in Shallow Water A Two Time Steps Kinetic Method , 2003 .

[31]  Benoît Perthame,et al.  A variant of Van Leer's method for multidimensional systems of conservation laws , 1994 .

[32]  Emmanuel Audusse,et al.  A Multilayer Saint-Venant Model , 2004 .

[33]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[34]  Marie-Odile Bristeau,et al.  Boundary Conditions for the Shallow Water Equations solved by Kinetic Schemes , 2001 .

[35]  J. Craggs Applied Mathematical Sciences , 1973 .

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

[37]  B. Perthame,et al.  A new kinetic scheme for Saint-Venant equations applied to debris avalanches , 2002 .

[38]  Kyōto Daigaku. Sūgakuka Lectures in mathematics , 1968 .

[39]  Sebastian Noelle,et al.  Shock waves, dead zones and particle-free regions in rapid granular free-surface flows , 2003, Journal of Fluid Mechanics.