A high-order relaxation method with projective integration for solving nonlinear systems of hyperbolic conservation laws

Abstract : We present a general, high-order, fully explicit relaxation scheme which can be applied to any system of nonlinear hyperbolic conservation laws in multiple dimensions. The scheme consists of two steps. In a first (relaxation) step, the nonlinear hyperbolic conservation law is approximated by a kinetic equation with stiff BGK source term. Then, this kinetic equation is integrated in time using a projective integration method. After taking a few small (inner) steps with a simple, explicit method (such as direct forward Euler) to damp out the stiff components of the solution, the time derivative is estimated and used in an (outer) Runge–Kutta method of arbitrary order. We show that, with an appropriate choice of inner step size, the time step restriction on the outer time step is similar to the CFL condition for the hyperbolic conservation law. Moreover, the number of inner time steps is also independent of the stiffness of the BGK source term. We discuss stability and consistency, and illustrate with numerical results (linear advection, Burgers' equation and the shallow water and Euler equations) in one and two spatial dimensions.
Document type :
Journal articles
Complete list of metadatas

https://hal.archives-ouvertes.fr/hal-01706186
Contributor : Pauline Lafitte-Godillon <>
Submitted on : Saturday, February 10, 2018 - 7:22:19 PM
Last modification on : Thursday, February 21, 2019 - 10:31:47 AM

Links full text

Identifiers

Citation

Pauline Godillon-Lafitte, Ward Melis, Giovanni Samaey. A high-order relaxation method with projective integration for solving nonlinear systems of hyperbolic conservation laws. Journal of Computational Physics, Elsevier, 2017, 340, pp.1 - 25. ⟨10.1016/j.jcp.2017.03.027⟩. ⟨hal-01706186⟩

Share

Metrics

Record views

84