A low-rank projector-splitting integrator for the Vlasov--Poisson equation

TL;DR

This paper introduces a low-rank splitting integrator for the Vlasov--Poisson equation that significantly reduces computational costs by approximating the solution within a low-rank manifold and splitting the dynamics into lower-dimensional advection problems.

Contribution

The paper develops a dynamical low-rank approximation method with a splitting integrator for the Vlasov--Poisson equation, enabling efficient simulation in high-dimensional phase space.

Findings

Reduces computational effort in plasma simulations.

Successfully applied to 2D and 4D problems like Landau damping.

Maintains accuracy with lower-dimensional advection equations.

Abstract

Many problems encountered in plasma physics require a description by kinetic equations, which are posed in an up to six-dimensional phase space. A direct discretization of this phase space, often called the Eulerian approach, has many advantages but is extremely expensive from a computational point of view. In the present paper we propose a dynamical low-rank approximation to the Vlasov--Poisson equation, with time integration by a particular splitting method. This approximation is derived by constraining the dynamics to a manifold of low-rank functions via a tangent space projection and by splitting this projection into the subprojections from which it is built. This reduces a time step for the six- (or four-) dimensional Vlasov--Poisson equation to solving two systems of three- (or two-) dimensional advection equations over the time step, once in the position variables and once in the velocity variables, where the size of each system of advection equations is equal to the chosen rank. By a hierarchical dynamical low-rank approximation, a time step for the Vlasov--Poisson equation can be further reduced to a set of six (or four) systems of one-dimensional advection equations, where the size of each system of advection equations is still equal to the rank. The resulting systems of advection equations can then be solved by standard techniques such as semi-Lagrangian or spectral methods. Numerical simulations in two and four dimensions for linear Landau damping, for a two-stream instability and for a plasma echo problem highlight the favorable behavior of this numerical method and show that the proposed algorithm is able to drastically reduce the required computational effort.