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.

Massively Parallel Computation of Accurate Densities for N-body Dark Matter Simulations using the Phase-Space-Element Method

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

Massively Parallel Multiphysics Solver for Simulating Large-Scale Antenna Array and Other 3-D Structures

High performance massively parallel direct N-body simulations on large GPU clusters

Supermassive Black Hole Binaries in High Performance Massively Parallel Direct N-body Simulations on Large GPU Clusters

Fast emulation of cosmological density fields based on dimensionality reduction and supervised machine-learning

Astrophysical particle simulations with large custom GPU clusters on three continents

How to Build and Use Special Purpose PC Clusters in Stellar Dynamics

Astrophysical Particle Simulations with Custom GPU Clusters

Parallel local time stepping for rigid bodies represented by triangulated meshes

Acceleration of Hybrid MPI Parallel NBODY6++ for Largen-Body Globular Cluster Simulations

Phase-Space Delaunay Tessellation Field Estimator

Up to 700k GPU Cores, Kepler, and the Exascale Future for Simulations of Star Clusters Around Black Holes.

A Stable Finite-Volume Method for Scalar-Field Dark Matter

Large Scale Plane Wave Pseudopotential Density Functional Theory Calculations on GPU Clusters

GPGPU Acceleration of All-Electron Electronic Structure Theory Using Localized Numeric Atom-Centered Basis Functions

From Newton to Einstein –N-Body Dynamics in Galactic Nuclei and SPH Using New Special Hardware and Astrogrid-D

Astrophysical SPH Simulations with raceSPH Library

An Efficient Multi-GPU Implementation for Linear-Response Time-Dependent Density Functional Theory

Up to 700k GPU cores, Kepler, and the Exascale future for simulations of star clusters around black holes

Cornerstone: Octree Construction Algorithms for Scalable Particle Simulations