Multi-Moment Advection scheme for Vlasov simulations
arXiv:1105.3061 · doi:10.1016/j.jcp.2011.05.010
Abstract
We present a new numerical scheme for solving the advection equation and its application to Vlasov simulations. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, for better conservation of the information entropy. We have developed one- and two-dimensional schemes and show that they provide quite accurate solutions within reasonable usage of computational resources compared to other existing schemes. The two-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile for more than hundred rotation periods with little numerical diffusion. This is crucially important for Vlasov simulations of magnetized plasmas. Applications of the one- and two-dimensional schemes to electrostatic and electromagnetic Vlasov simulations are presented with some benchmark tests.
52 pages, 18 figures, accepted for the publication in Journal of Computational Physics
References in corpus (3)
- Electron Shock Surfing Acceleration in Multidimensions: Two-dimensional Particle-In-Cell Simulation of Collisionless Perpendicular Shock
- Kinetic Vlasov Simulations of collisionless magnetic Reconnection
- Comparison of time splitting and backsubstitution methods for integrating Vlasov's equation with magnetic fields