Towards a compact high-order method for non-linear hyperbolic systems, II. The Hermite-HLLC scheme

In a finite-volume framework, we develop an approximate HLL Riemann solver specific to weakly hyperbolic systems. Those systems are obtained by considering not only the variable but also its first spatial derivative, as unknowns. To this aim, we rely upon the theory of ''@d-shock waves'', newly developed in the scalar case. First, we demonstrate that the extended version of the HLLE scheme to weakly hyperbolic systems is compatible with the existence of Dirac measures in the solution. Then, we develop a specific Hermite Least-Square (HLSM) interpolation that enables to generate a high-order and compact scheme, without creating spurious oscillations in the reconstruction of the variable or its first derivative. Extensive numerical experiments make it possible to validate the method and to check convergence to entropy solutions. Relying upon those results, we construct a new HLL Riemann solver, suited for the extended one-dimensional Euler equations. For this purpose, we introduce the contribution of a contact discontinuity inside the definition of the solver. By using a formal analogy with the scalar study, we demonstrate that this solver tolerates the existence of ''@d-shock waves'' in the solution. Numerical experiments that follow help to validate some of the assumptions made to generate this scheme.

[1]  Hailiang Liu,et al.  Formation of δ-Shocks and Vacuum States in the Vanishing Pressure Limit of Solutions to the Euler Equations for Isentropic Fluids , 2003, SIAM J. Math. Anal..

[2]  Chi-Wang Shu.  Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[3]  J. F. Mckenzie,et al.  Interaction of Linear Waves with Oblique Shock Waves , 1968 .

[4]  D. Schaeffer,et al.  Riemann problems for nonstrictly hyperbolic 2×2 systems of conservation laws , 1987 .

[5]  Guy Capdeville,et al.  Towards a compact high-order method for non-linear hyperbolic systems. I: The Hermite Least-Square Monotone (HLSM) reconstruction , 2009, J. Comput. Phys..

[6]  Jianxian Qiu,et al.  Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous Galerkin method: one-dimensional case , 2004 .

[7]  V. M. Shelkovich,et al.  The Riemann problem admitting δ-, δ-shocks, and vacuum states (the vanishing viscosity approach) , 2006 .

[8]  V. M. Shelkovich,et al.  Dynamics of propagation and interaction of δ-shock waves in conservation law systems , 2005 .

[9]  E Weinan,et al.  Generalized variational principles, global weak solutions and behavior with random initial data for systems of conservation laws arising in adhesion particle dynamics , 1996 .

[10]  Hanchun Yang.  Riemann Problems for a Class of Coupled Hyperbolic Systems of Conservation Laws , 1999 .

[11]  Tong Zhang,et al.  Delta-Shock Waves as Limits of Vanishing Viscosity for Hyperbolic Systems of Conservation Laws , 1994 .

[12]  P. Roe,et al.  On Godunov-type methods near low densities , 1991 .

[13]  V. Shelkovich,et al.  δ-Shock waves as a new type of solutions to systems of conservation laws , 2006 .

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

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

[16]  Derek M. Causon,et al.  On the Choice of Wavespeeds for the HLLC Riemann Solver , 1997, SIAM J. Sci. Comput..

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