Abstract:Harnessing the power of modern multi-GPU architectures, we present a massively parallel simulation system based on the Material Point Method (MPM) for simulating physical behaviors of materials undergoing complex topological changes, self-collision, and large deformations. Our system makes three critical contributions. First, we introduce a new particle data structure that promotes coalesced memory access patterns on the GPU and eliminates the need for complex atomic operations on the memory hierarchy when writing particle data to the grid. Second, we propose a kernel fusion approach using a new Grid-to-Particles-to-Grid ( G2P2G ) scheme, which efficiently reduces GPU kernel launches, improves latency, and significantly reduces the amount of global memory needed to store particle data. Finally, we introduce optimized algorithmic designs that allow for efficient sparse grids in a shared memory context, enabling us to best utilize modern multi-GPU computational platforms for hybrid Lagrangian-Eulerian computational patterns. We demonstrate the effectiveness of our method with extensive benchmarks, evaluations, and dynamic simulations with elastoplasticity, granular media, and fluid dynamics. In comparisons against an open-source and heavily optimized CPU-based MPM codebase [Fang et al. 2019] on an elastic sphere colliding scene with particle counts ranging from 5 to 40 million, our GPU MPM achieves over 100x per-time-step speedup on a workstation with an Intel 8086K CPU and a single Quadro P6000 GPU, exposing exciting possibilities for future MPM simulations in computer graphics and computational science. Moreover, compared to the state-of-the-art GPU MPM method [Hu et al. 2019a], we not only achieve 2x acceleration on a single GPU but our kernel fusion strategy and Array-of-Structs-of-Array ( AoSoA ) data structure design also generalizes to multi-GPU systems. Our multi-GPU MPM exhibits near-perfect weak and strong scaling with 4 GPUs, enabling performant and large-scale simulations on a 1024 3 grid with close to 100 million particles with less than 4 minutes per frame on a single 4-GPU workstation and 134 million particles with less than 1 minute per frame on an 8-GPU workstation.

Petascale turbulence simulation using a highly parallel fast multipole method on GPUs

Parallelized Implementation of the Finite Particle Method for Explicit Dynamics in GPU

Implementation of the moving particle semi-implicit method for free-surface flows on GPU clusters.

The Implementation of the Three-Dimensional Unified Gas-Kinetic Wave-Particle Method on Multiple Graphics Processing Units

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

A GPU-ready pseudo-spectral method for direct numerical simulations of multiphase turbulence

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Multiple-GPU accelerated high-order gas-kinetic scheme for direct numerical simulation of compressible turbulence

Toward Exascale: Overview of Large Eddy Simulations and Direct Numerical Simulations of Nuclear Reactor Flows with the Spectral Element Method in Nek5000

Method for scalable and performant GPU-accelerated simulation of multiphase compressible flow

Method for portable, scalable, and performant GPU-accelerated simulation of multiphase compressible flow

Efficient Magnetohydrodynamic Simulations on Graphics Processing Units with CUDA

Accelerating Dissipative Particle Dynamics Simulations on GPUs: Algorithms, Numerics and Applications

Strong Scaling for Numerical Weather Prediction at Petascale with the Atmospheric Model NUMA

A Peta-Scalable CPU-GPU Algorithm for Global Atmospheric Simulations

Quantum-Inspired Fluid Simulation of 2D Turbulence with GPU Acceleration

Gpu Speed-Up For The Implicit Navier-Stokes Solver

A massively parallel and scalable multi-GPU material point method

A FFT-based finite-difference solver for massively-parallel direct numerical simulations of turbulent flows

Efficient Computation of Large-Scale Statistical Solutions to Incompressible Fluid Flows

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