Abstract
In this paper, we present a semi-Lagrangian (SL) method based on a non-polynomial function space for solving the Vlasov equation. We find that a non-polynomial function based scheme is suitable to the specifics of the target problems. To address issues that arise in phase space models of plasma problems, we develop a weighted essentially non-oscillatory (WENO) scheme using trigonometric polynomials. In particular, the non-polynomial WENO method is able to achieve improved accuracy near sharp gradients or discontinuities. Moreover, to obtain a high-order of accuracy in not only space but also time, it is proposed to apply a high-order splitting scheme in time. We aim to introduce the entire SL algorithm with high-order splitting in time and high-order WENO reconstruction in space to solve the Vlasov-Poisson system. Some numerical experiments are presented to demonstrate robustness of the proposed method in having a high-order of convergence and in capturing non-smooth solutions. A key observation is that the method can capture phase structure that require twice the resolution with a polynomial based method. In 6D, this would represent a significant savings.
Original language | English |
---|---|
Pages (from-to) | 116-142 |
Number of pages | 27 |
Journal | Communications on Applied Mathematics and Computation |
Volume | 5 |
Issue number | 1 |
DOIs | |
Publication status | Published - Mar 2023 |
Bibliographical note
Publisher Copyright:© 2021, Shanghai University.
Keywords
- High-order splitting methods
- Non-polynomial basis
- Semi-Lagrangian methods
- Vlasov equation
- Vlasov-Poisson system
- WENO schemes