High-order multiderivative IMEX schemes

From MaRDI portal
Publication:2227687



Abstract: Recently, a 4th-order asymptotic preserving multiderivative implicit-explicit (IMEX) scheme was developed (Sch"utz and Seal 2020, arXiv:2001.08268). This scheme is based on a 4th-order Hermite interpolation in time, and uses an approach based on operator splitting that converges to the underlying quadrature if iterated sufficiently. Hermite schemes have been used in astrophysics for decades, particularly for N-body calculations, but not in a form suitable for solving stiff equations. In this work, we extend the scheme presented in Sch"utz and Seal 2020 to higher orders. Such high-order schemes offer advantages when one aims to find high-precision solutions to systems of differential equations containing stiff terms, which occur throughout the physical sciences. We begin by deriving Hermite schemes of arbitrary order and discussing the stability of these formulas. Afterwards, we demonstrate how the method of Sch"utz and Seal 2020 generalises in a straightforward manner to any of these schemes, and prove convergence properties of the resulting IMEX schemes. We then present results for methods ranging from 6th to 12th order and explore a selection of test problems, including both linear and nonlinear ordinary differential equations and Burgers' equation. To our knowledge this is also the first time that Hermite time-stepping methods have been applied to partial differential equations. We then discuss some benefits of these schemes, such as their potential for parallelism and low memory usage, as well as limitations and potential drawbacks.


The author presents a class of multiderivative IMEX (Implicit-Explicit) schemes of arbitrary order, focusing on orders 6--12 following [\textit{J. Schütz} and \textit{D. C. Seal}, Appl. Numer. Math. 160, 84--101 (2021; Zbl 1459.65159)] on the 4th-order method of this family. The computations in these schemes are highly parallelizable, and may be advantageous on current and future parallel architectures. In addition to this, many high-order multistep and Runge-Kutta methods require storing a number of intermediate stages of the solution, while the methods presented in this work only require storing values at the beginning and end of a time step. Another feature, such high-order schemes offer advantages when one aims to find high-precision solutions to systems of differential equations containing stiff terms, which occur in physics. Some numerical tests are presented to show the efficiency of the schemes proposed in this work.



Cites work



Describes a project that uses

Uses Software






This page was built for publication: High-order multiderivative IMEX schemes

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2227687)