Parallel-in-Time High-Order Multiderivative IMEX Solvers

In this work, we present a novel class of high-order time integrators for the numerical solution of ordinary differential equations. These integrators are of the two-derivative type, i.e., they take into account not only the first, but also the second temporal derivative of the unknown solution. While being motivated by two-derivative Runge–Kutta schemes, the methods themselves are of the predictor-corrector type, constructed in such a way that time-parallelism through pipelining is possible. This means that predictor and corrector steps can be computed simultaneously on different processors. Two variants are shown, the second-one being of Gauss–Seidel-type and hence having lower storage requirements. It turns out that this second variant is not only low-storage, but also performs better in numerical simulations. The algorithms presented can be cast as implicit two-step-two-derivative-multi-stage methods for which we present a detailed mathematical convergence analysis. Subsequently, numerical results are shown, demonstrating first the algorithms’ ability to cope with very stiff equations and showing up to eighth order of accuracy. Following, the parallel-in-time capabilities are illustrated. We conclude that the class of methods is very well-suited if the time parallelization of very stiff equations is an option.