Exponential Integrators for Resistive Magnetohydrodynamics: Matrix-free Leja Interpolation and Efficient Adaptive Time Stepping
arXiv:2108.13622 · doi:10.3847/1538-4365/ac5177
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.
arXiv commentAccepted for publication in ApJS; 20 pages, 13 figures
References in corpus (15)
- Cosmological Hydrodynamics with Adaptive Mesh Refinement: a new high resolution code called RAMSES
- PLUTO: a Numerical Code for Computational Astrophysics
- Athena: A New Code for Astrophysical MHD
- A High Order Godunov Scheme with Constrained Transport and Adaptive Mesh Refinement for Astrophysical MHD
- CASTRO: A New Compressible Astrophysical Solver. I. Hydrodynamics and Self-Gravity
- A Parallel Mesh-Adaptive Framework for Hyperbolic Conservation Laws
- The CRONOS Code for Astrophysical Magnetohydrodynamics
- Numerical study of Kelvin-Helmholtz instability and its impact on synthetic emission from magnetized jets
- On the performance of exponential integrators for problems in magnetohydrodynamics
- Magnus integrators on multicore CPUs and GPUs
- A splitting approach for the Kadomtsev--Petviashvili equation
- An accurate and time-parallel rational exponential integrator for hyperbolic and oscillatory PDEs
- An exponential integrator for the drift-kinetic model
- An adaptive step size controller for iterative implicit methods
- A modern resistive magnetohydrodynamics solver using C++ and the Boost library