Multidimensional Vlasov–Poisson Simulations with High-order Monotonicity- and Positivity-preserving Schemes
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577 (Japan)
- Department of Mathematical Science and Advanced Technology, Japan Agency for Marine-Earth Science and Technology, 3173-25, Syowa-machi, Kanazawa-ku, Yokohama, Kanagawa 236-0001 (Japan)
- Department of Physics, The University of Tokyo, Bunkyo, Tokyo 113-0033 (Japan)
We develop new numerical schemes for Vlasov–Poisson equations with high-order accuracy. Our methods are based on a spatially monotonicity-preserving (MP) scheme and are modified suitably so that the positivity of the distribution function is also preserved. We adopt an efficient semi-Lagrangian time integration scheme that is more accurate and computationally less expensive than the three-stage TVD Runge–Kutta integration. We apply our spatially fifth- and seventh-order schemes to a suite of simulations of collisionless self-gravitating systems and electrostatic plasma simulations, including linear and nonlinear Landau damping in one dimension and Vlasov–Poisson simulations in a six-dimensional phase space. The high-order schemes achieve a significantly improved accuracy in comparison with the third-order positive-flux-conserved scheme adopted in our previous study. With the semi-Lagrangian time integration, the computational cost of our high-order schemes does not significantly increase, but remains roughly the same as that of the third-order scheme. Vlasov–Poisson simulations on 128{sup 3}×128{sup 3} mesh grids have been successfully performed on a massively parallel computer.
- OSTI ID:
- 22875625
- Journal Information:
- Astrophysical Journal, Vol. 849, Issue 2; Other Information: Country of input: International Atomic Energy Agency (IAEA); ISSN 0004-637X
- Country of Publication:
- United States
- Language:
- English
Similar Records
Optimized 1d-1v Vlasov-Poisson simulations using a filtered Fourier-Hermite spectral discretization
A Semi-Lagrangian Spectral Method for the Vlasov–Poisson System Based on Fourier, Legendre and Hermite Polynomials