Exponential Integrators for Resistive Magnetohydrodynamics: Matrix-free Leja Interpolation and Efficient Adaptive Time Stepping

Abstract

We propose a novel algorithm for the temporal integration of the resistive magnetohydrodynamics (MHD) equations. The approach is based on exponential Rosenbrock schemes in combination with Leja interpolation. It naturally preserves Gauss’s law for magnetism and is unencumbered by the stability constraints observed for explicit methods. Remarkable progress has been achieved in designing exponential integrators and computing the required matrix functions efficiently. However, employing them in MHD simulations of realistic physical scenarios requires a matrix-free implementation. We show how an efficient algorithm based on Leja interpolation that only uses the right-hand side of the differential equation (i.e., matrix free) can be constructed. We further demonstrate that it outperforms Krylov-based exponential integrators as well as explicit and implicit methods using test models of magnetic reconnection and the Kelvin–Helmholtz instability. Furthermore, an adaptive step-size strategy that gives excellent and predictable performance, particularly in the lenient- to intermediate-tolerance regime that is often of importance in practical applications, is employed.

Combining DPG in space with DPG time-marching scheme for the transient advection–reaction equation

“…For a p-order method, the cost of sum-factorized matrix assembly is O(p 2d+1 ); so-called matrix-free approaches can compute the action of the matrix on a vector at a cost of O(p d+1 ) (see [3,Table 1]). Recently, Deka and Einkemmer have developed a matrix-free exponential integrator for resistive MHD [8]. Given the potential cost savings in high dimensions, we would like to develop a matrix-free exponential integrator analogous to the one in the present work.…”

A robust and conservative dynamical low-rank algorithm

“…This is, in particular, true for the examples we consider in this paper, where significantly more grid points are required in the v direction compared to the spatial directions. We also need to compute the right-hand side of equation (10). We have…”

A comparison of Leja- and Krylov-based iterative schemes for Exponential Integrators

“…This step size controller, or the cost controller from here on, has been optimised to incur minimum possible computational cost by choosing a smaller step size, than the one chosen by the traditional controller. The cost controller was first proposed for implicit integrators [14] and was later shown to be even more beneficial for exponential integrators [11,12]. This is because exponential integrators, in principle, are able to take much larger step sizes than implicit ones.…”