A performance comparison of semi-Lagrangian discontinuous Galerkin and spline based Vlasov solvers in four dimensions

JOURNAL OF COMPUTATIONAL PHYSICS

卷:376

Article

Einkemmer, Lukas^1,2 

[1] Univ Tubingen, Dept Math, Tubingen, Germany

[2] Univ Innsbruck, Dept Math, Tech Str 13, A-6020 Innsbruck, Austria

关键词: Semi-Lagrangian discontinuous Galerkin method;    Vlasov-Poisson equation;    Semi-Lagrangian methods;    Performance comparison;    Parallelization;   

DOI  :  10.1016/j.jcp.2018.10.012

来源: Elsevier

【 摘 要 】

The purpose of the present paper is to compare two semi-Lagrangian methods in the context of the four-dimensional Vlasov-Poisson equation. More specifically, our goal is to compare the performance of the more recently developed semi-Lagrangian discontinuous Galerkin scheme to cubic spline interpolation, which is widely used in Eulerian Vlasov simulation. To that end, we perform simulations for nonlinear Landau damping and a two-stream instability and provide benchmarks for the SeLaLib and SLDG codes, both on a workstation and using MPI on a cluster. In addition, we will present results for the graphic processing unit (GPU) implementation contained in SLDG. We find that the semi-Lagrangian discontinuous Galerkin scheme shows a moderate improvement in run time for nonlinear Landau damping and a substantial improvement for the two-stream instability. It should be emphasized that these results are markedly different from results obtained in the asymptotic regime, which favor spline interpolation. Thus, we conclude that the traditional approach of evaluating numerical methods is misleading, even for short time simulations. In addition, the absence of any global communication in the semi-Lagrangian discontinuous Galerkin method gives it a decisive advantage for scaling to more than 256 cores. (C) 2018 Elsevier Inc. All rights reserved.

【 授权许可】

Free