Vlasov-Poisson system tackled by particle simulation utilising boundary element methods
Abstract
This paper presents a grid-free simulation algorithm for the fully three-dimensional Vlasov--Poisson system for collisionless electron plasmas. We employ a standard particle method for the numerical approximation of the distribution function. Whereas the advection of the particles is grid-free by its very nature, the computation of the acceleration involves the solution of the non-local Poisson equation. To circumvent a volume mesh, we utilise the Fast Boundary Element Method, which reduces the three-dimensional Poisson equation to a system of linear equations on its two-dimensional boundary. This gives rise to fully populated matrices which are approximated by the $\mathcal H^2$-technique, reducing the computational time from quadratic to linear complexity. The approximation scheme based on interpolation has shown to be robust and flexible, allowing a straightforward generalisation to vector-valued functions. In particular, the Coulomb forces acting on the particles are computed in linear complexity. In first numerical tests, we validate our approach with the help of classical non-linear plasma phenomena. Furthermore, we show that our method is able to simulate electron plasmas in complex three-dimensional domains with mixed boundary conditions in linear complexity.
- Publication:
-
arXiv e-prints
- Pub Date:
- November 2018
- DOI:
- 10.48550/arXiv.1811.03404
- arXiv:
- arXiv:1811.03404
- Bibcode:
- 2018arXiv181103404K
- Keywords:
-
- Mathematics - Numerical Analysis;
- 35Q83;
- 65Z05;
- 65N75;
- 65N38;
- 68W25
- E-Print:
- final draft post-refereeing