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.

Co-design of a Particle-in-Cell Plasma Simulation Code for Intel Xeon Phi: A First Look at Knights Landing

Optimizing A Particle-In-Cell Code On Intel Knights Landing

The Plasma Simulation Code: A modern particle-in-cell code with load-balancing and GPU support

Toward a First-Principles Integrated Simulation of Tokamak Edge Plasmas

Optimization and Evaluation of VLPL-S on Knights Landing

The Plasma Simulation Code: A Modern Particle-in-cell Code with Patch-Based Load-Balancing.

Parallelization and I/O Optimization for a 3-D Plasma Simulation Code

Particle-in-cell plasma simulation on heterogeneous cluster systems

QuickPIC: a highly efficient fully parallelized PIC code for plasma-based acceleration

Particle-in-Cell Code for GPU Systems

Parallelization and Optimization of Laser-Plasma-Interaction Simulation

Modern Gyrokinetic Particle-In-Cell Simulation of Fusion Plasmas on Top Supercomputers

Exploiting Parallelism on Shared Memory in the QED Particle-in-Cell Code PICADOR with Greedy Load Balancing

A Particle-In-Cell Code Comparison for Ion Acceleration: EPOCH, LSP, and WarpX

PIFE-PIC: Parallel Immersed-Finite-Element Particle-In-Cell For 3-D Kinetic Simulations of Plasma-Material Interactions

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

A NVIDA kepler based acceleration of PIC method

Leveraging HPC Profiling & Tracing Tools to Understand the Performance of Particle-in-Cell Monte Carlo Simulations

Performance Analysis and Application of 2-D Plasma Parallel Particle-in-cell Simulations

MILC Code Performance on High End CPU and GPU Supercomputer Clusters

Photon-Plasma: a modern high-order particle-in-cell code