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-Accelerated Flow Simulations on Unstructured Grids Using a Multi-colored Gauss-Seidel Method

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

GPU-Accelerated Computations for Supersonic Flow Modeling on Hybrid Grids

Optimization and acceleration of flow simulations for CFD on CPU/GPU architecture

Development of GPU-based Parallel CFD Solver for Unstructured Meshes

OpenACC acceleration of an unstructured CFD solver based on a reconstructed discontinuous Galerkin method for compressible flows

Acceleration of a Production-Level Unstructured Grid Finite Volume CFD Code on GPU

The research of discontinuous galerkin p-multigrid solver

CPU/GPU Heterogeneous Parallel CFD Solver and Optimizations

Analysis of Parallel Multigrid Methods in Real-Time Fluid Simulation

Accelerating CFD simulation with high order finite difference method on curvilinear coordinates for modern GPU clusters

A GPU cluster optimized multigrid scheme for computing unsteady incompressible fluid flow

A scalable multi-GPU method for semi-implicit fractional-step integration of incompressible Navier-Stokes equations

A GPU accelerated three-dimensional ghost cell method with an improved implicit surface representation for complex rigid or flexible boundary flows

3D Parallel Multigrid Methods for Real-Time Fluid Simulation

A Multi-Block Viscous Flow Solver Based on GPU Parallel Methodology

The Simulation of Three-Dimensional Flow by Using GPU-based MPS Method

Implementation of the Moving Particle Semi-Implicit Method on GPU

A Cfd Solver On Graphical Processing Unites For Turbulence Simulations

An adaptive finite element multigrid solver using GPU acceleration

Numerical Simulation of Three-Dimensional Violent Free Surface Flows by GPU-Based MPS Method