A high order time splitting method based on integral deferred correction for semi-Lagrangian Vlasov simulations

JOURNAL OF COMPUTATIONAL PHYSICS

卷:267

Article

Christlieb, Andrew¹  Guo, Wei²  Morton, Maureen³  Qiu, Jing-Mei² 

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

[2] Univ Houston, Dept Math, Houston, TX 77204 USA

[3] Stark State Coll, Dept Math, North Canton, OH 44270 USA

关键词: Vlasov-Poisson system;    Guiding center model;    Operator/dimensional splitting;    Semi-Lagrangian method;    WENO reconstruction;    Integral deferred correction;   

DOI  :  10.1016/j.jcp.2014.02.012

来源: Elsevier

【 摘 要 】

Semi-Lagrangian schemes with various splitting methods, and with different reconstruction/interpolation strategies have been applied to kinetic simulations. For example, the order of spatial accuracy of the algorithms proposed in Qiu and Christlieb (2010) [29] is very high (as high as ninth order). However, the temporal error is dominated by the operator splitting error, which is second order for Strang splitting. It is therefore important to overcome such low order splitting error, in order to have numerical algorithms that achieve higher orders of accuracy in both space and time. In this paper, we propose to use the integral deferred correction (IDC) method to reduce the splitting error. Specifically, the temporal order accuracy is increased by r with each correction loop in the IDC framework, where r = 1, 2 for coupling the first order splitting and the Strang splitting, respectively. The proposed algorithm is applied to the Vlasov-Poisson system, the guiding center model, and two dimensional incompressible flow simulations in the vorticity stream-function formulation. We show numerically that the IDC procedure can automatically increase the order of accuracy in time. We also investigate numerical stability of the proposed algorithm via performing Fourier analysis to a linear model problem. (C) 2014 Elsevier Inc. All rights reserved.

【 授权许可】

Free