A High Order Time Splitting Method Based on Integral Deferred Correction for Semi-Lagrangian Vlasov Simulations

Abstract

Semi-Lagrangian schemes with various splitting methods, and with different reconstruction/interpolation strategies have been applied to kinetic simulations. For example, the order of spatial accuracy of the algorithms proposed in {[Qiu and Christlieb, J. Comp. Phys., 2010]} is very high (as high as ninth order). However, the temporal error is dominated by the operator splitting error, which is second order for Strang splitting. It is therefore important to overcome such low order splitting error, in order to have numerical algorithms that achieve higher orders of accuracy in both space and time. In this paper, we propose to use the integral deferred correction (IDC) method to reduce the splitting error. Specifically, the temporal order accuracy is increased by $r$ with each correction loop in the IDC framework, where $r=1,2$ for coupling the first order splitting and the Strang splitting, respectively. The proposed algorithm is applied to the Vlasov-Poisson system, the guiding center model, and two dimensional incompressible flow simulations in the vorticity stream-function formulation. We show numerically that the IDC procedure can automatically increase the order of accuracy in time. We also investigate numerical stability of the proposed algorithm via performing Fourier analysis to a linear model problem.

References in corpus (2)

• Energy-conserving discontinuous Galerkin methods for the Vlasov-Ampère system
• High order operator splitting methods based on an integral deferred correction framework

Cited by in corpus (12)

• A high order semi-Lagrangian discontinuous Galerkin method for Vlasov-Poisson simulations without operator splitting
• On the convergence of spectral deferred correction methods
• A WENO-based Method of Line Transpose Approach for Vlasov Simulations
• High-order partitioned spectral deferred correction solvers for multiphysics problems
• Quantification of mixing in vesicle suspensions using numerical simulations in two dimensions
• Numerical Study of the Two-Species Vlasov-Ampère System: Energy-Conserving Schemes and the Current-Driven Ion-Acoustic Instability
• Arbitrary-order time-accurate semi-Lagrangian spectral approximations of the Vlasov-Poisson system
• A Semi-Lagrangian Spectral Method for the Vlasov-Poisson System based on Fourier, Legendre and Hermite Polynomials
• Conservative Multi-Dimensional Semi-Lagrangian Finite Difference Scheme: Stability and Applications to the Kinetic and Fluid Simulations
• A high order conservative semi-Lagrangian discontinuous Galerkin method for two-dimensional transport simulations
• Implicit-Explicit Integral Deferred Correction Methods for Stiff Problems
• Comparison of semi-Lagrangian discontinuous Galerkin schemes for linear and nonlinear transport simulations