An exponential integrator for the drift-kinetic model
arXiv:1705.09923 · doi:10.1016/j.cpc.2017.11.003
Abstract
We propose an exponential integrator for the drift-kinetic equation in cylindrical geometry. This approach removes the CFL condition from the linear part of the system (which is often the most stringent requirement in practice) and treats the remainder explicitly using Arakawa's finite difference scheme. The present approach is mass conservative, up to machine precision, and significantly reduces the computational effort per time step. In addition, we demonstrate the efficiency of our method by performing numerical simulations in the context of the ion temperature gradient instability. In particular, we find that our numerical method can take time steps comparable to what has been reported in the literature for the (predominantly used) splitting approach. In addition, the proposed numerical method has significant advantages with respect to conservation of energy and efficient higher order methods can be obtained easily. We demonstrate this by investigating the performance of a fourth order implementation.
References in corpus (8)
- A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov-Poisson equations
- Convergence analysis of Strang splitting for Vlasov-type equations
- Convergence analysis of a discontinuous Galerkin/Strang splitting approximation for the Vlasov--Poisson equations
- A splitting approach for the Kadomtsev--Petviashvili equation
- A conservative discontinuous Galerkin scheme for the 2D incompressible Navier--Stokes equations
- A study on conserving invariants of the Vlasov equation in semi-Lagrangian computer simulations
- Towards optimal explicit time-stepping schemes for the gyrokinetic equations
- Mixed semi-Lagrangian/finite difference methods for plasma simulations
Cited by in corpus (7)
- Exponential methods for solving hyperbolic problems with application to kinetic equations
- An accurate and time-parallel rational exponential integrator for hyperbolic and oscillatory PDEs
- A -mode integrator for solving evolution equations in Kronecker form
- Efficient adaptive step size control for exponential integrators
- LeXInt: Package for Exponential Integrators employing Leja interpolation
- Exponential Integrators for Resistive Magnetohydrodynamics: Matrix-free Leja Interpolation and Efficient Adaptive Time Stepping
- Post-Moore Technologies for Plasma Simulation: A Community Roadmap