High Order Semi-Lagrangian Discontinuous Galerkin Method Coupled with Runge-Kutta Exponential Integrators for Nonlinear Vlasov Dynamics

In this paper, we propose a semi-Lagrangian discontinuous Galerkin method coupled with Runge-Kutta exponential integrators (SLDG-RKEI) for nonlinear Vlasov dynamics. The commutator-free Runge-Kutta (RK) exponential integrators (EI) were proposed by Celledoni, et al. (FGCS, 2003). In the nonlinear transport setting, the RKEI can be used to decompose the evolution of the nonlinear transport into a composition of a sequence of linearized dynamics. The resulting linearized transport equations can be solved by the semi-Lagrangian (SL) discontinuous Galerkin (DG) method proposed in Cai, et al. (JSC, 2017). The proposed method can achieve high order spatial accuracy via the SLDG framework, and high order temporal accuracy via the RK EI. Due to the SL nature, the proposed SLDG-RKEI method is not subject to the CFL condition, thus they have the potential in using larger time-stepping sizes than those in the Eulerian approach. Inheriting advantages from the SLDG method, the proposed SLDG-RKEI schemes are mass conservative, positivity-preserving, have no dimensional splitting error, perform well in resolving complex solution structures, and can be evolved with adaptive time-stepping sizes. We show the performance of the SLDG-RKEI algorithm by classical test problems for the nonlinear Vlasov-Poisson system, as well as the Guiding center Vlasov model. Though that it is not our focus of this paper to explore the SLDG-RKEI scheme for nonlinear hyperbolic conservation laws that develop shocks, we show some preliminary results on schemes' performance on the Burgers' equation.

[1]  Elena Celledoni,et al.  Commutator-free Lie group methods , 2003, Future Gener. Comput. Syst..

[2]  Jennifer K. Ryan,et al.  Extension of a Post Processing Technique for the Discontinuous Galerkin Method for Hyperbolic Equations with Application to an Aeroacoustic Problem , 2005, SIAM J. Sci. Comput..

[3]  Zhengfu Xu,et al.  High order maximum principle preserving semi-Lagrangian finite difference WENO schemes for the Vlasov equation , 2013, J. Comput. Phys..

[4]  Jianxian Qiu,et al.  An h-Adaptive RKDG Method for the Two-Dimensional Incompressible Euler Equations and the Guiding Center Vlasov Model , 2017, J. Sci. Comput..

[5]  Wei Guo,et al.  A high order semi-Lagrangian discontinuous Galerkin method for Vlasov-Poisson simulations without operator splitting , 2017, J. Comput. Phys..

[6]  O. Pironneau On the transport-diffusion algorithm and its applications to the Navier-Stokes equations , 1982 .

[7]  Ilaria Perugia,et al.  An A Priori Error Analysis of the Local Discontinuous Galerkin Method for Elliptic Problems , 2000, SIAM J. Numer. Anal..

[8]  Elena Celledoni,et al.  Semi-Lagrangian Runge-Kutta Exponential Integrators for Convection Dominated Problems , 2009, J. Sci. Comput..

[9]  Wei Guo,et al.  A High Order Semi-Lagrangian Discontinuous Galerkin Method for the Two-Dimensional Incompressible Euler Equations and the Guiding Center Vlasov Model Without Operator Splitting , 2018, Journal of Scientific Computing.

[10]  Giovanni Russo,et al.  A High Order Multi-Dimensional Characteristic Tracing Strategy for the Vlasov–Poisson System , 2017, J. Sci. Comput..

[11]  Wei Guo,et al.  A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations , 2016, Journal of Scientific Computing.

[12]  Elena Celledoni,et al.  High Order Semi-Lagrangian Methods for the Incompressible Navier–Stokes Equations , 2012, J. Sci. Comput..

[13]  Magdi Shoucri A two-level implicit scheme for the numerical solution of the linearized vorticity equation , 1981 .

[14]  Chi-Wang Shu,et al.  The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems , 1998 .

[15]  David C. Seal,et al.  A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov-Poisson equations , 2010, J. Comput. Phys..

[16]  Jianxian Qiu,et al.  An h-Adaptive RKDG Method for the Vlasov–Poisson System , 2016, J. Sci. Comput..

[17]  Eric Sonnendrücker,et al.  Conservative semi-Lagrangian schemes for Vlasov equations , 2010, J. Comput. Phys..

[18]  Xiangxiong Zhang,et al.  On maximum-principle-satisfying high order schemes for scalar conservation laws , 2010, J. Comput. Phys..

[19]  Lorenzo Pareschi,et al.  Modeling and Computational Methods for Kinetic Equations , 2012 .

[20]  Todd Arbogast,et al.  A semi-Lagrangian finite difference WENO scheme for scalar nonlinear conservation laws , 2016, J. Comput. Phys..

[21]  B. K. Kometa On Semi-Lagrangian Exponential Integrators and Discontinuous Galerkin Methods , 2011 .

[22]  Douglas N. Arnold,et al.  Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems , 2001, SIAM J. Numer. Anal..