A study on conserving invariants of the Vlasov equation in semi-Lagrangian computer simulations

Abstract

The semi-Lagrangian discontinuous Galerkin method, coupled with a splitting approach in time, has recently been introduced for the Vlasov–Poisson equation. Since these methods are conservative, local in space and able to limit numerical diffusion, they are considered a promising alternative to more traditional semi-Lagrangian schemes. In this paper we study the conservation of important physical invariants and the long-time behaviour of the semi-Lagrangian discontinuous Galerkin method. To that end we conduct a theoretical analysis and perform a number of numerical simulations. In particular, we find that the entropy is non-decreasing for the discontinuous Galerkin scheme, while unphysical oscillations in the entropy are observed for the traditional method based on cubic spline interpolation.