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.

Hybrid Decomposition Method in Parallel Molecular Dynamics Simulation Based on SMP Cluster Architecture

Cluster-Based parallel simulation for large scale molecular dynamics in microscale thermophysics

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

Parallel Computation Of Meshfree Methods For Extremely Large Deformation Analysis

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

An Efficient Domain Decomposition Method in Multi-physics Simulation of Package Systems

A New Dynamical Domain Decomposition Method for Parallel Molecular Dynamics Simulation on Grid

A new hierarchical parallel algorithm on molecular dynamics simulation

MPI+X:Massive Parallelization and Dynamic Load Balance of a Production-level Unstructured DSMC Solver

Efficient Parallel Implementation of Molecular Dynamics with Embedded Atom Method on Multi-core Platforms

Scalable and fast heterogeneous molecular simulation with predictive parallelization schemes

Parallel approach of three-dimensional phase-field model of binary alloy in MPI+OpenMP environment

Optimization techniques for parallel force-decomposition algorithm in molecular dynamic simulations

Multi-GPU Hybrid Programming Accelerated Three-Dimensional Phase-Field Model in Binary Alloy

Simulation computation of MPI-based parallel PF-LBM 3-D dendritic growth model

Parallel SPH modeling using dynamic domain decomposition and load balancing displacement of Voronoi subdomains

Hybrid PSM/FDM method for seismic wavefield simulation

Employing multi-GPU power for molecular dynamics simulation: an extension of GALAMOST

Parallelization of 3D Explicit Material Point Method Using OpenMP

Parallel Algorithm for Particle-Grid Dual Discretization

Development of a Hybrid Parallel MCV-based High-Order Global Shallow-Water Model