A "metric" semi-Lagrangian Vlasov-Poisson solver

Abstract

We propose a new semi-Lagrangian Vlasov-Poisson solver. It employs elements of metric to follow locally the flow and its deformation, allowing one to find quickly and accurately the initial phase-space position $Q(P)$ of any test particle $P$, by expanding at second order the geometry of the motion in the vicinity of the closest element. It is thus possible to reconstruct accurately the phase-space distribution function at any time $t$ and position $P$ by proper interpolation of initial conditions, following Liouville theorem. When distorsion of the elements of metric becomes too large, it is necessary to create new initial conditions along with isotropic elements and repeat the procedure again until next resampling. To speed up the process, interpolation of the phase-space distribution is performed at second order during the transport phase, while third order splines are used at the moments of remapping. We also show how to compute accurately the region of influence of each element of metric with the proper percolation scheme. The algorithm is tested here in the framework of one-dimensional gravitational dynamics but is implemented in such a way that it can be extended easily to four or six-dimensional phase-space. It can also be trivially generalised to plasmas.

arXiv comment32 pages, 14 figures, accepted for publication in Journal of Plasma Physics, Special issue: The Vlasov equation, from space to laboratory plasmas

References in corpus (13)

• Discreteness effects in simulations of Hot/Warm dark matter
• A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov-Poisson equations
• N-body simulations of gravitational dynamics
• Simulation techniques for cosmological simulations
• An adaptively refined phase-space element method for cosmological simulations and collisionless dynamics
• Discreteness effects in cosmological N-body simulations
• ColDICE: a parallel Vlasov-Poisson solver using moving adaptive simplicial tessellation
• Direct Integration of the Collisionless Boltzmann Equation in Six-dimensional Phase Space: Self-gravitating Systems
• Energy ejection in the collapse of a cold spherical self-gravitating cloud
• A cloudy Vlasov solution
• Scaling up Tides in Numerical Models of Galaxy- and Halo-Formation
• Arbitrarily high order Convected Scheme solution of the Vlasov-Poisson system
• Vlasov versus N-body: the Hénon sphere

Cited by in corpus (5)

• Large-scale dark matter simulations
• Cosmological Vlasov-Poisson Simulations of Structure Formation with Relic Neutrinos: Nonlinear Clustering and the Neutrino Mass
• Cosmological Vlasov-Poisson equations for dark matter: Recent developments and connections to selected plasma problems
• A $\delta f$ PIC method with Forward-Backward Lagrangian reconstructions
• Axion Structure Formation II: The Wrath of Collapse