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.

Acceleration Of Three-Dimensional Tokamak Magnetohydrodynamical Code With Graphics Processing Unit And Openacc Heterogeneous Parallel Programming

Heterogeneous Programming and Optimization of Gyrokinetic Toroidal Code and Large-Scale Performance Test on TH-1A.

Multi-GPU thermal lattice Boltzmann simulations using OpenACC and MPI

Multi-GPU Hybrid Programming Accelerated Three-Dimensional Phase-Field Model in Binary Alloy

OpenACC Vs the Native Programming on Sunway TaihuLight: A Case Study with GTC-P

A GPU-Accelerated Modern Fortran Version of the ECHO Code for Relativistic Magnetohydrodynamics

OpenACC acceleration of an unstructured CFD solver based on a reconstructed discontinuous Galerkin method for compressible flows

Study on Acceleration of Three-Dimensional Method of Characteristics by GPU

An Efficient Heterogenous Parallel Algorithm of the 3D MOC for Multizone Heterogeneous Systems

Efficient Magnetohydrodynamic Simulations on Graphics Processing Units with CUDA

Particle-resolved thermal lattice Boltzmann simulation using OpenACC on multi-GPUs

OpenACC offloading of the MFC compressible multiphase flow solver on AMD and NVIDIA GPUs

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Multiple-GPU accelerated high-order gas-kinetic scheme on three-dimensional unstructured meshes

Performance Portable Graphics Processing Unit Acceleration of a High-Order Finite Element Multiphysics Application

Hybrid simulation of energetic particles interacting with magnetohydrodynamics using a slow manifold algorithm and GPU acceleration

Accelerating a three-dimensional MOC calculation using GPU with CUDA and two-level GCMFD method

Implicit Low-Order Unstructured Finite-Element Multiple Simulation Enhanced by Dense Computation using OpenACC

High Performance GPU Accelerated MuST Software

Accelerating an Iterative Eigensolver for Nuclear Structure Configuration Interaction Calculations on GPUs Using OpenACC