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.

Venice : A multi-scale operator-splitting algorithm for multi-physics simulations

Venice: a multi-scale operator-splitting algorithm for multi-physics simulations

Multi-physics simulations using a hierarchical interchangeable software interface

Space-Time-Parallel Poroelasticity Simulation

Enhancing the description of multi-time-scale geophysical phenomena: Incorporating finite time Scale separation and feedback, illustrated with the case of a 1D variable of interest

Operator splitting for coupled linear port-Hamiltonian systems

Intelligent Time-Scale Operator-Splitting Integration for Chemical Reaction Systems

An adaptive multiphysics model coupling vertical equilibrium and full multidimensions for multiphase flow in porous media

Scaling up Probabilistic PDE Simulators with Structured Volumetric Information

Scalable and fast heterogeneous molecular simulation with predictive parallelization schemes

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

Commutator-based operator splitting for linear port-Hamiltonian systems

Nonlocal multicontinua with Representative Volume Elements. Bridging separable and non-separable scales

A framework for multi-scale modelling

Effectiveness of the operator splitting for solving the atmospherical shallow water equations

On Splitting-Based Numerical Methods for Convection-Diffusion Equations

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

High order finite volume schemes with IMEX time stepping for the Boltzmann model on unstructured meshes

Multiscale Simulation Based on the Meshless Local Petrov-Galerkin (mlpg) Method

Multi-scale S-fraction reduced-order models for massive wavefield simulations

A unified operator splitting approach for multi-scale fluid-particle coupling in the lattice Boltzmann method