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.

GPN-based simulation of seismic wave propagation with hybrid PSM/FDM method

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.

Three-dimensional parallel hybrid PSM/FDM simulation of seismic wave propagation

Hybrid PSM/FDM method for seismic wavefield simulation

Simulation of the Seismic Response of 2D Sedimentary Basin with Hybrid PSM/FDM Method

3D wave equation prestack depth migration on Graphic Processing Unit

Implementation of the Moving Particle Semi-Implicit Method on GPU

Modeling Seismic Wave Propagation in Heterogeneous Medium Using Overlap Domain Pseudospectral Method

GPU-accelerated smoothed particle finite element method for large deformation analysis in geomechanics

Hybrid MPI and CUDA Parallelization for CFD Applications on Multi-GPU HPC Clusters

A multi-GPU and CUDA-aware MPI-based spectral element formulation for ultrasonic wave propagation in solid media

Analysis method for seismic response of high-rise structure based on CPU-GPU heterogeneous platform

Acceleration of a 2D/3D Finite-Discrete Element Code for Geomechanical Simulations Using General Purpose GPU Computing

Parallel seismic propagation simulation in anisotropic media by irregular grids finite difference method on PC cluster

GPU Based Two-Level CMFD Accelerating Two-Dimensional MOC Neutron Transport Calculation

Multidomain Finite-Difference and Chebyshev Pseudospectral Hybrid Method for Elastic Wave Simulation: Two-Dimensional Case

DEM–SPH Coupling Method for Landslide Surge Based on a GPU Parallel Acceleration Technique

Development of a GPGPU‐parallelized hybrid finite‐discrete element method for modeling rock fracture

Matrix‐based implementation and GPU acceleration of hybrid FEM and peridynamic model for hydro‐mechanical coupled problems

High Efficient Sedimentary Basin Simulations on Hybrid CPU-GPU Clusters