Well-Balanced High-Order Centred Schemes for Non-Conservative Hyperbolic Systems. Applications to Shallow Water Equations with Fixed and Mobile Bed

This paper concerns the development of high-order accurate centred schemes for the numerical solution of one-dimensional hyperbolic systems containing non-conservative products and source terms. Combining the PRICE-T method developed in [Toro E, Siviglia A. PRICE: primitive centred schemes for hyperbolic system of equations. Int J Numer Methods Fluids 2003;42:1263–91] with the theoretical insights gained by the recently developed path-conservative schemes [Castro M, Gallardo J, Pares C. High-order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products applications to shallow-water systems. Math Comput 2006;75:1103–34; Pares C. Numerical methods for nonconservative hyperbolic systems: a theoretical framework. SIAM J Numer Anal 2006;44:300–21], we propose the new PRICE-C scheme that automatically reduces to a modified conservative FORCE scheme if the underlying PDE system is a conservation law. The resulting first-order accurate centred method is then extended to high order of accuracy in space and time via the ADER approach together with a WENO reconstruction technique. The well-balanced properties of the PRICE-C method are investigated for the shallow water equations. Finally, we apply the new scheme to the shallow water equations with fix bottom topography and with variable bottom solving an additional sediment transport equation.

[1]  Michael Dumbser,et al.  Finite volume schemes of very high order of accuracy for stiff hyperbolic balance laws , 2008, J. Comput. Phys..

[2]  Michael Dumbser,et al.  A unified framework for the construction of one-step finite volume and discontinuous Galerkin schemes on unstructured meshes , 2008, J. Comput. Phys..

[3]  Michael Dumbser,et al.  Building Blocks for Arbitrary High Order Discontinuous Galerkin Schemes , 2006, J. Sci. Comput..

[4]  Eleuterio F. Toro,et al.  PRICE: primitive centred schemes for hyperbolic systems , 2003 .

[5]  Eleuterio F. Toro,et al.  CENTERED DIFFERENCE SCHEMES FOR NONLINEAR HYPERBOLIC EQUATIONS , 2004 .

[6]  E. Toro,et al.  Solution of the generalized Riemann problem for advection–reaction equations , 2002, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[7]  I. Toumi.  A weak formulation of roe's approximate riemann solver , 1992 .

[8]  Valerio Caleffi,et al.  High-order balanced CWENO scheme for movable bed shallow water equations , 2007 .

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

[10]  Michael Dumbser,et al.  Arbitrary high order non-oscillatory finite volume schemes on unstructured meshes for linear hyperbolic systems , 2007, J. Comput. Phys..

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

[12]  Carlos Parés Madroñal,et al.  Numerical methods for nonconservative hyperbolic systems: a theoretical framework , 2006, SIAM J. Numer. Anal..

[13]  J. Glimm.  Solutions in the large for nonlinear hyperbolic systems of equations , 1965 .

[14]  Enrique D. Fernández-Nieto,et al.  Sediment transport models in Shallow Water equations and numerical approach by high order finite volume methods , 2008 .

[15]  A. Siviglia,et al.  Quasi-Conservative Formulation of the One-Dimensional Saint-Venant — Exner Model , 2008 .

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

[17]  G. Parker.  Surface-based bedload transport relation for gravel rivers , 1990 .

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

[19]  M. Altinakar,et al.  St. Venant-Exner Equations for Near-Critical and Transcritical Flows , 2002 .

[20]  R. Müller,et al.  Formulas for Bed-Load transport , 1948 .

[21]  Michael Dumbser,et al.  Quadrature-free non-oscillatory finite volume schemes on unstructured meshes for nonlinear hyperbolic systems , 2007, J. Comput. Phys..

[22]  Manuel Jesús Castro Díaz,et al.  High order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products. Applications to shallow-water systems , 2006, Math. Comput..

[23]  Eleuterio F. Toro,et al.  Centred TVD schemes for hyperbolic conservation laws , 2000 .

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

[25]  T. Hou,et al.  Why nonconservative schemes converge to wrong solutions: error analysis , 1994 .

[26]  Yulong Xing,et al.  High-Order Well-Balanced Finite Difference WENO Schemes for a Class of Hyperbolic Systems with Source Terms , 2006, J. Sci. Comput..

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

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

[29]  Manuel Jesús Castro Díaz,et al.  Why many theories of shock waves are necessary: Convergence error in formally path-consistent schemes , 2008, J. Comput. Phys..