Abstract:The moving particle semi-implicit (MPS) method performs well in simulating incompressible flows with free surfaces. Despite its applicability, the MPS method suffers from the fundamental instability problem and high computational cost in its practical application. Substantial research has been conducted on improving the stability and accuracy of the MPS method. Moreover, graphics processing units (GPUs), which are multi-processors that execute many three-dimensional geometric processes at high speed, provide unprecedented capability for scientific computations. However, the usage of a single GPU card is not sufficient for engineering applications that require several million particles that predict the desired physical processes, because the available memory space is insufficient. In this work, the dynamic stability (DS) algorithm and particle shifting (PS) algorithm have been used to overcome the instability and inaccuracies caused by tensile instability and non-uniform particle distribution, respectively. Based on the stable MPS method, a GPU-based MPS code that uses the compute unified device architecture (CUDA) language has been developed. An efficient neighborhood particle search is performed using an indirect method, and the matrix for the pressure Poisson equation (PPE) is assembled in parallel. Based on the single-GPU version, a multi-GPU MPS code has been developed. The approach uses a non-geometric dynamic domain decomposition method that provides homogeneous load balancing whereby different portions (subdomains) of the physical system under study are assigned to different GPUs. Communication between devices is achieved with the use of a message passing interface (MPI). Based on the neighborhood particle search, the techniques for building and updating the “halo” are described in detail. The speed-up of the single-GPU version is analyzed for different numbers of particles, and the scalability of the multi-GPU version is analyzed for different numbers of particles and different numbers of GPUs. Last, an application with more than 107 particles is presented to show the capability of the code in handling large-scale simulations.

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

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

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

Toward a First-Principles Integrated Simulation of Tokamak Edge Plasmas

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

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

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

Particle-in-cell plasma simulation on heterogeneous cluster systems

Particle-in-Cell Code for GPU Systems

A Study of Performance Portability in Plasma Physics Simulations

Optimizing BIT1, a Particle-in-Cell Monte Carlo Code, with OpenMP/OpenACC and GPU Acceleration

Ultrahigh performance three-dimensional electromagnetic relativistic kinetic plasma simulation

Pegasus: A New Hybrid-Kinetic Particle-in-cell Code for Astrophysical Plasma Dynamics.

Smilei : A collaborative, open-source, multi-purpose particle-in-cell code for plasma simulation

The muphyII Code: Multiphysics Plasma Simulation on Large HPC Systems

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

High-order Implicit Particle-in-cell Method for Plasma Simulations at Solid Densities.

Maximus: a Hybrid Particle-in-Cell Code for Microscopic Modeling of Collisionless Plasmas

Understanding Large-Scale Plasma Simulation Challenges for Fusion Energy on Supercomputers

VPIC 2.0: Next Generation Particle-in-Cell Simulations

SPACE: 3D Parallel Solvers for Vlasov-Maxwell and Vlasov-Poisson Equations for Relativistic Plasmas with Atomic Transformations