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.

GPU-based Parallel Algorithm for Particle Contact Detection and Its Application in Self-Compacting Concrete Flow Simulations

A GPU-Based Parallel Algorithm for 2D Large Deformation Contact Problems Using the Finite Particle Method

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

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

GPU-accelerated Vector-Form Particle-Element Method for 3D Elastoplastic Contact of Structures

A Novel GPGPU-parallelized Contact Detection Algorithm for Combined Finite-Discrete Element Method

An efficient and robust GPGPU-parallelized contact algorithm for the combined finite-discrete element method

GPU-accelerated discrete element simulation of granular and gas-solid flows

A GPU-based DEM Approach for Modelling of Particulate Systems

Analysis of parallel spatial partitioning algorithms for GPU based DEM

A GPU-based DEM framework for simulation of polyhedral particulate system

GPU-accelerated DEM implementation with CUDA

A GPU Based Hybrid Material Point and Discrete Element Method (MPDEM) Algorithm and Validation

A GPU-based DEM Model for the Pebble Flow Study in Packed Bed: Simulation Scheme and Validation

Parallelized Combined Finite‐discrete Element (FDEM) Procedure Using Multi‐gpu with CUDA

Multi-GPUs DEM algorithm and its application in the simulation of granular materials

GPU-accelerated DEM implementation with CUDA.

A GPU-based DEM for Modelling Large Scale Powder Compaction with Wide Size Distributions

Computing Framework of Large-Scale Dilated Polyhedral DEM Particles Simulation Based on GPU

A GPU-based DEM-FEM Computational Framework for Tire-Sand Interaction Simulations

A GPU-based SPH Algorithm