A Simple Extension of the Osher Riemann Solver to Non-conservative Hyperbolic Systems

We propose a simple extension of the well-known Riemann solver of Osher and Solomon (Math. Comput. 38:339–374, 1982) to a certain class of hyperbolic systems in non-conservative form, in particular to shallow-water-type and multi-phase flow models. To this end we apply the formalism of path-conservative schemes introduced by Parés (SIAM J. Numer. Anal. 44:300–321, 2006) and Castro et al. (Math. Comput. 75:1103–1134, 2006). For the sake of generality and simplicity, we suggest to compute the inherent path integral numerically using a Gaussian quadrature rule of sufficient accuracy. Published path-conservative schemes to date are based on either the Roe upwind method or on centered approaches. In comparison to these, the proposed new path-conservative Osher-type scheme has several advantages. First, it does not need an entropy fix, in contrast to Roe-type path-conservative schemes. Second, our proposed non-conservative Osher scheme is very simple to implement and nonetheless constitutes a complete Riemann solver in the sense that it attributes a different numerical viscosity to each characteristic field present in the relevant Riemann problem; this is in contrast to centered methods or incomplete Riemann solvers that usually neglect intermediate characteristic fields, hence leading to excessive numerical diffusion. Finally, the interface jump term is differentiable with respect to its arguments, which is useful for steady-state computations in implicit schemes. We also indicate how to extend the method to general unstructured meshes in multiple space dimensions. We show applications of the first order version of the proposed path-conservative Osher-type scheme to the shallow water equations with variable bottom topography and to the two-fluid debris flow model of Pitman & Le. Then, we apply the higher-order multi-dimensional version of the method to the Baer–Nunziato model of compressible multi-phase flow. We also clearly emphasize the limitations of our approach in a special chapter at the end of this article.

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

[2]  Eleuterio F. Toro,et al.  On some fast well-balanced first order solvers for nonconservative systems , 2009, Math. Comput..

[3]  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..

[4]  M. J. Castro,et al.  ADER schemes on unstructured meshes for nonconservative hyperbolic systems: Applications to geophysical flows , 2009 .

[5]  M. Baer,et al.  A two-phase mixture theory for the deflagration-to-detonation transition (ddt) in reactive granular materials , 1986 .

[6]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

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

[8]  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..

[9]  G. Warnecke,et al.  The Riemann problem for the Baer-Nunziato two-phase flow model , 2004 .

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

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

[12]  G. D. Maso,et al.  Definition and weak stability of nonconservative products , 1995 .

[13]  Donald W. Schwendeman,et al.  The Riemann problem and a high-resolution Godunov method for a model of compressible two-phase flow , 2006, J. Comput. Phys..

[14]  M. J. Castro,et al.  FORCE schemes on unstructured meshes II: Non-conservative hyperbolic systems , 2010 .

[15]  Miltiadis V. Papalexandris,et al.  An exact Riemann solver for compressible two-phase flow models containing non-conservative products , 2007, J. Comput. Phys..

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

[17]  R. Abgrall,et al.  A Multiphase Godunov Method for Compressible Multifluid and Multiphase Flows , 1999 .

[18]  Rémi Abgrall,et al.  A comment on the computation of non-conservative products , 2010, J. Comput. Phys..

[19]  Long Le,et al.  A two-fluid model for avalanche and debris flows , 2005, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[20]  S. Osher,et al.  Upwind difference schemes for hyperbolic systems of conservation laws , 1982 .

[21]  R. Abgrall How to Prevent Pressure Oscillations in Multicomponent Flow Calculations , 1996 .

[22]  Eleuterio F. Toro,et al.  Exact solution of the Riemann problem for the shallow water equations with discontinuous bottom geometry , 2008, J. Comput. Phys..

[23]  Michael Dumbser,et al.  A new 3D parallel SPH scheme for free surface flows , 2009 .

[24]  Angela Ferrari,et al.  SPH simulation of free surface flow over a sharp-crested weir , 2010 .

[25]  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..

[26]  A. Stroud Approximate calculation of multiple integrals , 1973 .

[27]  Michael Dumbser,et al.  Well-Balanced High-Order Centred Schemes for Non-Conservative Hyperbolic Systems. Applications to Shallow Water Equations with Fixed and Mobile Bed , 2009 .

[28]  P. García-Navarro,et al.  On numerical treatment of the source terms in the shallow water equations , 2000 .

[29]  Sander Rhebergen,et al.  Discontinuous Galerkin finite element methods for hyperbolic nonconservative partial differential equations , 2008, J. Comput. Phys..

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

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

[32]  Eleuterio F. Toro,et al.  ADER schemes for three-dimensional non-linear hyperbolic systems , 2005 .

[33]  Marica Pelanti,et al.  A Roe-type scheme for two-phase shallow granular flows over variable topography , 2008 .

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

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

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

[37]  C. Parés Numerical methods for nonconservative hyperbolic systems: a theoretical framework. , 2006 .

[38]  E. Toro,et al.  Well-balanced high-order centered schemes on unstructured meshes for shallow water equations with fixed and mobile bed , 2010 .

[39]  S. Osher Riemann Solvers, the Entropy Condition, and Difference , 1984 .

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