Exponential methods for solving hyperbolic problems with application to kinetic equations

Abstract

The efficient numerical solution of many kinetic models in plasma physics is impeded by the stiffness of these systems. Exponential integrators are attractive in this context as they remove the CFL condition induced by the linear part of the system, which in practice is often the most stringent stability constraint. In the literature, these schemes have been found to perform well, e.g., for drift-kinetic problems. Despite their overall efficiency and their many favorable properties, most of the commonly used exponential integrators behave rather erratically in terms of the allowed time step size in some situations. This severely limits their utility and robustness.Our goal in this paper is to explain the observed behavior and suggest exponential methods that do not suffer from the stated deficiencies. To accomplish this we study the stability of exponential integrators for a linearized problem. This analysis shows that classic exponential integrators exhibit severe deficiencies in that regard. Based on the analysis conducted we propose to use Lawson methods, which can be shown not to suffer from the same stability issues. We confirm these results and demonstrate the efficiency of Lawson methods by performing numerical simulations for both the Vlasov-Poisson system and a drift-kinetic model of a ion temperature gradient instability.

References in corpus (17)

• A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov-Poisson equations
• Structure and structure-preserving algorithms for plasma physics
• A Hamiltonian splitting for the Vlasov-Maxwell system
• The Leja method revisited: backward error analysis for the matrix exponential
• Convergence analysis of Strang splitting for Vlasov-type equations
• A comparison of semi-Lagrangian discontinuous Galerkin and spline based Vlasov solvers in four dimensions
• Convergence analysis of a discontinuous Galerkin/Strang splitting approximation for the Vlasov--Poisson equations
• High performance computing aspects of a dimension independent semi-Lagrangian discontinuous Galerkin code
• On the performance of exponential integrators for problems in magnetohydrodynamics
• A strategy to suppress recurrence in grid-based Vlasov solvers
• A splitting approach for the Kadomtsev--Petviashvili equation
• Semi-Lagrangian Vlasov simulation on GPUs
• An exponential integrator for the drift-kinetic model
• A study on conserving invariants of the Vlasov equation in semi-Lagrangian computer simulations
• An adaptive step size controller for iterative implicit methods
• An asymptotic preserving scheme for the relativistic Vlasov--Maxwell equations in the classical limit
• Exponential Integrators on Graphic Processing Units

Cited by in corpus (5)

• An accurate and time-parallel rational exponential integrator for hyperbolic and oscillatory PDEs
• A $\mu$-mode integrator for solving evolution equations in Kronecker form
• LeXInt: Package for Exponential Integrators employing Leja interpolation
• Efficient adaptive step size control for exponential integrators
• An efficient stochastic particle method for high-dimensional nonlinear PDEs