Abstract:The Vlasov-Poisson equations describe the evolution of a collisionless plasma, represented through a probability density function (PDF) that self-interacts via an electrostatic force. One of the main difficulties in numerically solving this system is the severe time-step restriction that arises from parts of the PDF associated with moderate-to-large velocities. The dominant approach in the plasma physics community for removing these time-step restrictions is the so-called particle-in-cell (PIC) method, which discretizes the distribution function into a set of macro-particles, while the electric field is represented on a mesh. Several alternatives to this approach exist, including fully Lagrangian, fully Eulerian, and so-called semi-Lagrangian methods. The focus of this work is the semi-Lagrangian approach, which begins with a grid-based Eulerian representation of both the PDF and the electric field, then evolves the PDF via Lagrangian dynamics, and finally projects this evolved field back onto the original Eulerian mesh. In particular, we develop in this work a method that discretizes the 1+1 Vlasov-Poisson system via a high-order discontinuous Galerkin (DG) method in phase space, and an operator split, semi-Lagrangian method in time. Second-order accuracy in time is relatively easy to achieve via Strang operator splitting. With additional work, using higher-order splitting and a higher-order method of characteristics, we also demonstrate how to push this scheme to fourth-order accuracy in time. We show how to resolve all of the Lagrangian dynamics in such a way that mass is exactly conserved, positivity is maintained, and high-order accuracy is achieved. The Poisson equation is solved to high-order via the smallest stencil local discontinuous Galerkin (LDG) approach. We test the proposed scheme on several standard test cases.

A Numerical Method For Solving The One-Dimensional Vlasov-Poisson Equation In Phase Space

A Novel Method for Solving the Linearized 1D Vlasov--Poisson Equation

A Splitting Fourier Pseudospectral Method for Vlasov–Poisson–Fokker–Planck System

A discontinuous Galerkin method for the Vlasov-Poisson system

A numerical integration scheme for vectorised phase-space of one-dimensional collision-free, electrostatic systems

Conservative stabilized Runge-Kutta methods for the Vlasov-Fokker-Planck equation

A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov-Poisson equations

A Vlasov description of the Euler equation

A numerical method for computing time-periodic solutions in dissipative wave systems

A Characteristic Mapping Method for Vlasov-Poisson with Extreme Resolution Properties

1D Vlasov-Poisson Equations with Electron Sheet Initial Data

Anti-symmetric and Positivity Preserving Formulation of a Spectral Method for Vlasov-Poisson Equations

Smooth Solution To The One-Dimensional Inhomogeneous Non-Automorphic Landau-Lifshitz Equation

6D phase space collective modes in Vlasov-Maxwell system

Asymptotically stable Particle-in-Cell methods for the magnetized Vlasov--Poisson equations in orthogonal curvilinear coordinates

A massively parallel semi-Lagrangian solver for the six-dimensional Vlasov-Poisson equation

Asymptotic-preserving Schemes for Kinetic-Fluid Modeling of Disperse Two-Phase Flows

Numerical study of an anisotropic Vlasov equation arising in plasma physics

Solving Vlasov-Maxwell Equations By Using Hamiltonian Splitting

Comparison of two different integration methods for the (1+1)-Dimensional Schrödinger-Poisson Equation