Comparison of time splitting and backsubstitution methods for integrating Vlasov's equation with magnetic fields

Abstract

The standard approach for integrating the multidimensional Vlasov equation using grid based, conservative schemes is based on a time splitting approach. Here, we show that although the truncation error is of second order, time splitting can introduce systematic heating of the plasma. We introduce a backsubstitution method, which not only avoids this deficiency but also is computationally less expensive. The general approach is demonstrated in conjunction with Boris' scheme for evaluating the characteristics.

to appear in Comp. Phys. Comm

Cited by in corpus (15)

• Kinetic Vlasov Simulations of collisionless magnetic Reconnection
• Stable interpolation with isotropic and anisotropic Gaussians using Hermite generating function
• A massively parallel semi-Lagrangian solver for the six-dimensional Vlasov-Poisson equation
• Numerical simulations of the Fourier transformed Vlasov-Maxwell system in higher dimensions --- Theory and applications
• Coupled Vlasov and two-fluid codes on GPUs
• Temperature gradient driven heat flux closure in fluid simulations of collisionless reconnection
• A Parallel Low-Rank Solver for the Six-Dimensional Vlasov-Maxwell Equations
• Multi-Moment Advection scheme for Vlasov simulations
• An Energy Conserving Vlasov Solver That Tolerates Coarse Velocity Space Resolutions: Simulation of MMS Reconnection Events
• The muphyII Code: Multiphysics Plasma Simulation on Large HPC Systems
• Fluid Simulations of Three-Dimensional Reconnection that Capture the Lower-Hybrid Drift Instability
• Enhanced conservation properties of Vlasov codes through coupling with conservative fluid models
• Convergence of splitting methods on rotating grids for the magnetized Vlasov equation
• Optimal Landau-type closure parameters for two-fluid simulations of plasma turbulence at kinetic scales
• Multiphysics simulations of collisionless plasmas