A massively parallel semi-Lagrangian solver for the six-dimensional Vlasov-Poisson equation

03/01/2019
by   Katharina Kormann, et al.
0

This paper presents an optimized and scalable semi-Lagrangian solver for the Vlasov-Poisson system in six-dimensional phase space. Grid-based solvers of the Vlasov equation are known to give accurate results. At the same time, these solvers are challenged by the curse of dimensionality resulting in very high memory requirements, and moreover, requiring highly efficient parallelization schemes. In this paper, we consider the 6d Vlasov-Poisson problem discretized by a split-step semi-Lagrangian scheme, using successive 1d interpolations on 1d stripes of the 6d domain. Two parallelization paradigms are compared, a remapping scheme and a classical domain decomposition approach applied to the full 6d problem. From numerical experiments, the latter approach is found to be superior in the massively parallel case in various respects. We address the challenge of artificial time step restrictions due to the decomposition of the domain by introducing a blocked one-sided communication scheme for the purely electrostatic case and a rotating mesh for the case with a constant magnetic field. In addition, we propose a pipelining scheme that enables to hide the costs for the halo communication between neighbor processes efficiently behind useful computation. Parallel scalability on up to 65k processes is demonstrated for benchmark problems on a supercomputer.

READ FULL TEXT
research
03/03/2017

Parallel energy-stable phase field crystal simulations based on domain decomposition methods

In this paper, we present a parallel numerical algorithm for solving the...
research
02/29/2020

Convergence estimates of a semi-Lagrangian scheme for the ellipsoidal BGK model for polyatomic molecules

In this paper, we propose a new semi-Lagrangian scheme for the polyatomi...
research
10/15/2019

An explicit semi-Lagrangian, spectral method for solution of Lagrangian transport equations in Eulerian-Lagrangian formulations

An explicit high order semi-Lagrangian method is developed for solving L...
research
07/24/2023

Semi-Lagrangian Scheme with Arakawa Splitting for Gyro-kinetic Equations

The gyro-kinetic model is an approximation of the Vlasov-Maxwell system ...
research
09/23/2022

A mass-conserving sparse grid combination technique with biorthogonal hierarchical basis functions for kinetic simulations

The exact numerical simulation of plasma turbulence is one of the assets...
research
04/29/2021

Distributed Multigrid Neural Solvers on Megavoxel Domains

We consider the distributed training of large-scale neural networks that...
research
09/12/2019

PittPack: An Open-Source Poisson's Equation Solver for Extreme-Scale Computing with Accelerators

We present a parallel implementation of a direct solver for the Poisson'...

Please sign up or login with your details

Forgot password? Click here to reset