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.

Graphics processing unit (GPU)-enhanced nonhydrostatic model with grid nesting for global tsunami propagation and coastal inundation

A Numerical Study of Tsunami Wave Impact and Run-Up on Coastal Cliffs Using a CIP-based Model

Computationally Efficient Modeling of Hydro-Sediment-morphodynamic Processes Using a Hybrid Local Time Step/global Maximum Time Step

Numerical Investigation on the Hydrodynamic Characteristics of Coastal Freak Wave Using a CIP-based Model

Development of Two-Dimensional Non-Hydrostatic Wave Model Based on Central-Upwind Scheme

Research on the single-host parallel computing with the local time step scheme for modeling of hydro-sediment-morphodynamic processes

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

A typhoon-induced storm surge numerical model with GPU acceleration based on an unstructured spherical centroidal Voronoi tessellation grid

Study on Multi-Time Scale Hydrodynamic Model Based on Local Time Stepping Scheme and GPUs and its Application in Urban Inundation

An Engineering-Oriented Shallow-water Hydro-Sediment-Morphodynamic Model Using the GPU-acceleration and the Hybrid LTS/GMaTS Method

A GPU-accelerated and LTS-based 2D hydrodynamic model for the simulation of rainfall-runoff processes

A well-balanced positivity preserving two-dimensional shallow flow model with wetting and drying fronts over irregular topography

A GPU-Accelerated and LTS-Based Finite Volume Shallow Water Model

2d Well-Balanced Positivity Preserving Shallow Flow Model With Wetting And Drying On Triangular Grids

A Well-Balanced Numerical Model for the Simulation of Long Waves over Complex Domains

Three-Dimensional Numerical Method for Simulating Large-Scale Free Water Surface by Massive Parallel Computing on a GPU

Modeling and Visualization of Tsunamis

Design and evaluation of an efficient high-precision ocean surface wave model with a multiscale grid system (MSG_Wav1.0)

A GPU Accelerated Finite Volume Coastal Ocean Model

GPU‐Accelerated Urban Flood Modeling Using a Nonuniform Structured Grid and a Super Grid Scale River Channel

Simulating Dam-Break Flooding with Floating Objects Through Intricate City Layouts Using Gpu-Based Sph Method