High-order multiderivative time integrators for hyperbolic conservation laws

Abstract

Multiderivative time integrators have a long history of development for ordinary differential equations, and yet to date, only a small subset of these methods have been explored as a tool for solving partial differential equations (PDEs). This large class of time integrators include all popular (multistage) Runge-Kutta as well as single-step (multiderivative) Taylor methods. (The latter are commonly referred to as Lax-Wendroff methods when applied to PDEs.) In this work, we offer explicit multistage multiderivative time integrators for hyperbolic conservation laws. Like Lax-Wendroff methods, multiderivative integrators permit the evaluation of higher derivatives of the unknown in order to decrease the memory footprint and communication overhead. Like traditional Runge-Kutta methods, multiderivative integrators admit the addition of extra stages, which introduce extra degrees of freedom that can be used to increase the order of accuracy or modify the region of absolute stability. We describe a general framework for how these methods can be applied to two separate spatial discretizations: the discontinuous Galerkin (DG) method and the finite difference essentially non-oscillatory (FD-WENO) method. The two proposed implementations are substantially different: for DG we leverage techniques that are closely related to generalized Riemann solvers; for FD-WENO we construct higher spatial derivatives with central differences. Among multiderivative time integrators, we argue that multistage two-derivative methods have the greatest potential for multidimensional applications, because they only require the flux function and its Jacobian, which is readily available. Numerical results indicate that multiderivative methods are indeed competitive with popular strong stability preserving time integrators.

Journal of Scientific Computing

References in corpus (1)

• A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov-Poisson equations

Cited by in corpus (12)

• A Compact Fourth-order Gas-kinetic Scheme for the Euler and Navier-Stokes Solutions
• Finite Difference Weighted Essentially Non-Oscillatory Schemes with Constrained Transport for Ideal Magnetohydrodynamics
• A Family of High-order Gas-kinetic Schemes and Its Comparison with Riemann Solver Based High-order Methods
• A gradient-compression-based compact high-order gas-kinetic scheme on three-dimensional hybrid unstructured mesh
• Two-derivative deferred correction time discretization for the discontinuous Galerkin method
• An acoustic and shock wave capturing compact high-order gas-kinetic scheme with spectral-like resolution
• Jacobian-free explicit multiderivative Runge-Kutta methods for hyperbolic conservation laws
• High resolution compact implicit numerical scheme for conservation laws
• A single-step third-order temporal discretization with Jacobian-free and Hessian-free formulations for finite difference methods
• Functional-preserving predictor-corrector multiderivative schemes
• Multiderivative time integration methods preserving nonlinear functionals via relaxation
• Unconditionally stable higher order semi-implicit level set method for advection equations