A WENO-based Method of Lines Transpose approach for Vlasov simulations

JOURNAL OF COMPUTATIONAL PHYSICS

卷:327

Article

Christlieb, Andrew¹  Guo, Wei²  Jiang, Yan² 

[1] Michigan State Univ, Dept Computat Math Sci & Engn, E Lansing, MI 48824 USA

[2] Michigan State Univ, Dept Math, E Lansing, MI 48824 USA

关键词: Method of Lines Transpose;    Implicit Runge-Kutta methods;    Weighted essentially non-oscillatory methodology;    High order accuracy;    Vlasov-Poisson;    Demensional splitting;    Positivity-preserving;   

DOI  :  10.1016/j.jcp.2016.09.048

来源: Elsevier

【 摘 要 】

In this paper, a high order implicit Method of Lines Transpose (MOLT) method based on a weighted essentially non-oscillatory (WENO) methodology is developed for onedimensional linear transport equations and further applied to the Vlasov-Poisson (VP) simulations via dimensional splitting. In the MOLT framework, the time variable is first discretized by a diagonally implicit strong-stability-preserving Runge-Kutta method, resulting in a boundary value problem (BVP) at the discrete time levels. Then an integral formulation coupled with a high order WENO methodology is employed to solve the BVP. As a result, the proposed scheme is high order accurate in both space and time and free of oscillations even though the solution is discontinuous or has sharp gradients. Moreover, the scheme is able to take larger time step evolution compared with an explicit MOL WENO scheme with the same order of accuracy. The desired positivity-preserving (PP) property of the scheme is further attained by incorporating a newly proposed high order PP limiter. We perform numerical experiments on several benchmarks including linear advection, solid body rotation problem; and on the Landau damping, two-stream instabilities, bump-on-tail, and plasma sheath by solving the VP system. The efficacy and efficiency of the proposed scheme is numerically verified. (C) 2016 Elsevier Inc. All rights reserved.

【 授权许可】

Free