Abstract:Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in physical space. The computational cost of evolving high-dimensional variables often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we developed techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We present theoretical hardware compute and communication bounds required for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU. We demonstrate that the multi-GPU Vlasov solver implementation VCK-GPU simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits, leaving minimal room for theoretical performance improvements. This resulted in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. We then demonstrate how the improved compute performance can be used to explore configurations which were previously computationally infeasible such as resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

The fast VFP code for solution of the Vlasov-Fokker-Planck equation

Higher-order moment convergent method in weakly anisotropic plasma and the NLVFP code for solution of the 0D-2V Vlasov-Fokker-Planck equation

A fast solver for the spatially homogeneous electron Boltzmann equation

A Vlasov-Fokker-Planck-Landau code for the simulation of colliding supersonic dense plasma flows

A Splitting Fourier Pseudospectral Method for Vlasov–Poisson–Fokker–Planck System

High Fidelity simulations of the multi-species Vlasov equation in the electro-static, collisional-less limit

A Parallel Low-Rank Solver for the Six-Dimensional Vlasov-Maxwell Equations

Vlasov-Fokker-Planck-Maxwell Simulations for Plasmas in Inertial Confinement Fusion

Vlasov methods in space physics and astrophysics

Relativistic Vlasov-Maxwell modelling using finite volumes and adaptive mesh refinement

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic Solvers with Adaptive Mesh in Phase Space for Low-Temperature Plasmas

A Fokker-Planck Code for Laser-Produced Plasmas

Relativistic Vlasov Code Development for High Energy Density Plasmas

SPACE: 3D Parallel Solvers for Vlasov-Maxwell and Vlasov-Poisson Equations for Relativistic Plasmas with Atomic Transformations

An efficient implicit boundary integral solver for the Vlasov-Maxwell system

Porting the grid-based 3D+3V hybrid-Vlasov kinetic plasma simulation Vlasiator to heterogeneous GPU architectures

A Massively Parallel Performance Portable Free-space Spectral Poisson Solver

A performance portable implementation of the semi-Lagrangian algorithm in six dimensions

Ultrahigh performance three-dimensional electromagnetic relativistic kinetic plasma simulation