A simplified Cauchy-Kowalewskaya procedure for the implicit solution of generalized Riemann problems of hyperbolic balance laws

The Cauchy-Kowalewskaya (CK) procedure is a key building block in the design of solvers for the Generalised Rieman Problem (GRP) based on Taylor series expansions in time. The CK procedure allows us to express time derivatives in terms of purely space derivatives. This is a very cumbersome procedure, which often requires the use of software manipulators. In this paper, a simplification of the CK procedure is proposed in the context of implicit Taylor series expansion for GRP, for hyperbolic balance laws in the framework of [Journal of Computational Physics 303 (2015) 146-172]. A recursive formula for the CK procedure, which is straightforwardly implemented in computational codes, is obtained. The proposed GRP solver is used in the context of the ADER approach and several one-dimensional problems are solved to demonstrate the applicability and efficiency of the present scheme. An enhancement in terms of efficiency, is obtained. Furthermore, the expected theoretical orders of accuracy are achieved, conciliating accuracy and stability.

[1]  Michael Dumbser,et al.  Very high order PNPM schemes on unstructured meshes for the resistive relativistic MHD equations , 2009, J. Comput. Phys..

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

[3]  Rémi Abgrall,et al.  Multidimensional HLLC Riemann solver for unstructured meshes - With application to Euler and MHD flows , 2014, J. Comput. Phys..

[4]  Claus R. Goetz,et al.  Approximate solutions of generalized Riemann problems: The Toro-Titarev solver and the LeFloch-Raviart expansion , 2012 .

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

[6]  Michael Dumbser,et al.  Comparison of solvers for the generalized Riemann problem for hyperbolic systems with source terms , 2012, J. Comput. Phys..

[7]  Eleuterio F. Toro,et al.  Reformulations for general advection-diffusion-reaction equations and locally implicit ADER schemes , 2014, J. Comput. Phys..

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

[9]  Dinshaw S. Balsara.  Multidimensional HLLE Riemann solver: Application to Euler and magnetohydrodynamic flows , 2010, J. Comput. Phys..

[10]  Dinshaw S. Balsara,et al.  An efficient, second order accurate, universal generalized Riemann problem solver based on the HLLI Riemann solver , 2018, J. Comput. Phys..

[11]  G. Sod.  A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws , 1978 .

[12]  Eleuterio F. Toro,et al.  Pollutant transport by shallow water equations on unstructured meshes: Hyperbolization of the model and numerical solution via a novel flux splitting scheme , 2016, Journal of Computational Physics.

[13]  I. Men’shov,et al.  The generalized problem of breakup of an arbitrary discontinuity , 1991 .

[14]  P. Raviart,et al.  An asymptotic expansion for the solution of the generalized Riemann problem. Part 2 : application to the equations of gas dynamics , 1989 .

[15]  P. Raviart,et al.  An asymptotic expansion for the solution of the generalized Riemann problem Part I: General theory , 1988 .

[16]  Michael Dumbser,et al.  ADER discontinuous Galerkin schemes for aeroacoustics , 2005 .

[17]  Claus-Dieter Munz,et al.  ADER: A High-Order Approach for Linear Hyperbolic Systems in 2D , 2002, J. Sci. Comput..

[18]  Eleuterio F. Toro,et al.  Towards Very High Order Godunov Schemes , 2001 .

[19]  J. Falcovitz,et al.  A second-order Godunov-type scheme for compressible fluid dynamics , 1984 .

[20]  Michael Dumbser,et al.  A family of HLL-type solvers for the generalized Riemann problem , 2017, Computers & Fluids.

[21]  Eleuterio F. Toro,et al.  ADER: Arbitrary High Order Godunov Approach , 2002, J. Sci. Comput..

[22]  Randall J. LeVeque,et al.  A study of numerical methods for hyperbolic conservation laws with stiff source terms , 1990 .

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

[24]  Eleuterio F. Toro,et al.  Implicit, semi-analytical solution of the generalized Riemann problem for stiff hyperbolic balance laws , 2015, J. Comput. Phys..

[25]  S. Osher,et al.  Uniformly High-Order Accurate Nonoscillatory Schemes. I , 1987 .

[26]  Michael Dumbser,et al.  A new efficient formulation of the HLLEM Riemann solver for general conservative and non-conservative hyperbolic systems , 2016, J. Comput. Phys..

[27]  Michael Dumbser,et al.  Fast high order ADER schemes for linear hyperbolic equations , 2004 .

[28]  M. Dumbser,et al.  An arbitrary high-order discontinuous Galerkin method for elastic waves on unstructured meshes - I. The two-dimensional isotropic case with external source terms , 2006 .

[29]  Eleuterio F. Toro,et al.  ADER: Arbitrary-Order Non-Oscillatory Advection Schemes , 2001 .

[30]  Eleuterio F. Toro,et al.  Solvers for the high-order Riemann problem for hyperbolic balance laws , 2008, J. Comput. Phys..

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

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

[33]  Dinshaw S. Balsara.  A two-dimensional HLLC Riemann solver for conservation laws: Application to Euler and magnetohydrodynamic flows , 2012, J. Comput. Phys..

[34]  Gino I. Montecinos,et al.  A strategy to implement Dirichlet boundary conditions in the context of ADER finite volume schemes. One-dimensional conservation laws , 2016, ArXiv.

[35]  Jiequan Li,et al.  Hyperbolic balance laws: Riemann invariants and the generalized Riemann problem , 2007, Numerische Mathematik.

[36]  I. S. Men'shov.  Increasing the order of approximation of Godunov's scheme using solutions of the generalized riemann problem , 1990 .

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

[38]  E. F. Toro,et al.  Primitive, Conservative and Adaptive Schemes for Hyperbolic Conservation Laws , 1998 .

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

[40]  D. Balsara,et al.  A cell-centered polynomial basis for efficient Galerkin predictors in the context of ADER finite volume schemes. The one-dimensional case , 2017 .

[41]  Michael Dumbser,et al.  A Novel Solver for the Generalized Riemann Problem Based on a Simplified LeFloch–Raviart Expansion and a Local Space–Time Discontinuous Galerkin Formulation , 2016, J. Sci. Comput..

[42]  Eleuterio F. Toro,et al.  Arbitrarily Accurate Non-Oscillatory Schemes for Nonlinear Scalar Conservation Laws with Source Terms II , 2002 .