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.

Implementation of a Lattice Boltzmann Method for Large Eddy Simulation on Multiple GPUs

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

The Implementation of the Three-Dimensional Unified Gas-Kinetic Wave-Particle Method on Multiple Graphics Processing Units

Direct Numerical Simulation and Large Eddy Simulation on a Turbulent Wall-Bounded Flow Using Lattice Boltzmann Method and Multiple GPUs

Multi-GPU Performance of Incompressible Flow Computation by Lattice Boltzmann Method on GPU Cluster.

Direct Numerical Simulation of the Wall-Bounded Turbulent Flow by Lattice Boltzmann Method Based on Multi-GPU

GPU implementation of lattice Boltzmann method for flows with curved boundaries

Accelerated Lattice Boltzmann Simulation Using GPU and OpenACC with Data Management

The Immersed Boundary-Lattice Boltzmann Method Parallel Model for Fluid-Structure Interaction on Heterogeneous Platforms

An efficient flux-reconstructed lattice boltzmann flux solver for flow interaction of multi-structure with curved boundary

Accelerating Large-Scale CFD Simulations with Lattice Boltzmann Method on a 40-Million-core Sunway Supercomputer.

Numerical Simulation of a 2D Electrothermal Pump by Lattice Boltzmann Method on GPU

Application of Lattice Boltzmann Methods for the Multiphase Fluid Pipe Flow on Graphical Processing Unit

Accelerating Unstructured Large Eddy Simulation Solver with GPU

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

Performance engineering for the Lattice Boltzmann method on GPGPUs: Architectural requirements and performance results

High-Performance Numerical Simulation of Jet in Cross-Flow Based on Lattice Boltzmann Method

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

Multi-GPU thermal lattice Boltzmann simulations using OpenACC and MPI

Three-Dimensional Numerical Simulation of Droplet Evaporation Using the Lattice Boltzmann Method Based on GPU-CUDA Accelerated Algorithm

Implementation of a direct-addressing based lattice Boltzmann GPU solver for multiphase flow in porous media