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.

Accelerating the simulation of kinetic shear Alfvén waves with a dynamical low-rank approximation

Accelerating the simulation of kinetic shear Alfvén waves with a dynamical low-rank approximation

A Gyrokinetic Simulation Model for Low Frequency Electromagnetic Fluctuations in Magnetized Plasmas

A stable multiplicative dynamical low-rank discretization for the linear Boltzmann-BGK equation

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

Interpolatory dynamical low-rank approximation for the 3+3d Boltzmann-BGK equation

Splitting scheme for gyro-kinetic equations with Semi-Lagrangian and Arakawa substeps

A Gyrokinetic Electron and Fully Kinetic Ion Plasma Simulation Model

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

Accelerating Kinetic Simulations of Electrostatic Plasmas with Reduced-Order Modeling

Robust and conservative dynamical low-rank methods for the Vlasov equation via a novel macro-micro decomposition

Low-rank approximation in the numerical modeling of the Farley-Buneman instability in ionospheric plasma

High-order Discretization of a Gyrokinetic Vlasov Model in Edge Plasma Geometry

A Deterministic Dynamical Low-rank Approach for Charged Particle Transport

Low rank compression in the numerical solution of the nonequilibrium Dyson equation

Finite-mass Fluid Electron Simulation of Kinetic and Inertial Alfvén Waves in a Sheared Magnetic Field

Nonlinear gyrokinetic equations for low-frequency electromagnetic waves in general plasma equilibria

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

An Extended Hybrid Magnetohydrodynamics Gyrokinetic Model For Numerical Simulation Of Shear Alfven Waves In Burning Plasmas

Asymptotic-Preserving and Energy Stable Dynamical Low-Rank Approximation

Low rank approximation for the numerical simulation of high dimensional Lindblad and Riccati equations