Structure preserving numerical methods for the Vlasov equation

TL;DR

Long-time Vlasov–Poisson integration must handle high-dimensional, nonlinear dynamics across widely separated timescales while preserving physically relevant invariants. The paper studies semi-Lagrangian discontinuous Galerkin methods, compares them with cubic spline interpolation, and reports two-stream-instability simulations indicating better entropy behavior for sLdG.

• Problem

  Vlasov–Poisson simulation is challenging because it may involve six-dimensional phase space, nonlinearity, and physical timescales much larger than the plasma period, while long-time integration requires preserving important invariants.

• Method

  The paper studies semi-Lagrangian discontinuous Galerkin methods for Vlasov–Poisson and compares them with cubic spline interpolation using characteristic-based advection.

• Results

  Numerical simulations indicate that cubic spline interpolation decreases entropy, whereas the semi-Lagrangian discontinuous Galerkin scheme does not show this deficiency; higher-order approximations also provide significant advantages.

• Takeaways & Limitations

  The study evaluates numerical methods by tracking mass, momentum, energy, positivity, entropy, and L2-norm behavior over long integrations.

• Takeaways & Limitations

  The report considers the electrostatic Vlasov–Poisson model, with the potential used to determine the electric field.