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.

Investigation of Single Bubble Behavior under Rolling Motions Using Multiphase MPS Method on GPU

Numerical Study of Bubble Rising in Quiescent Liquid under Rolling Motion Conditions Using MPS – MAFL Method

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

Numerical Study on the Single Bubble Rising Behaviors under Rolling Conditions

An Improved Computational Fluid Dynamics-Discrete Bubble Model for Multiphase Bubble Flow Simulation with a Wide Grid/bubble Size Ratio

Numerical Study of Laminar Flow and Friction Characteristics in Narrow Channels under Rolling Conditions Using MPS Method

NUMERICAL SIMULATION OF BUBBLE DYNAMIC UNDER OCEAN CONDITIONS

Investigation of the bubble deflection using graphics processing unit accelerated multiple relaxation time lattice Boltzmann method: The effect of bubble size and liquid viscosity

Numerical Simulation on Bubble Growth and Departure by MPS

Numerical Simulation of Taylor Bubble Formation in Micro-channel by MPS Method

A Multiphase MPS Method for Bubbly Flows with Complex Interfaces

Numerical Analysis of Bubble Rising Behavior in a Liquid Metal Using MPS

Simulation of Single Bubble Dynamic Process in Pool Boiling Process under Microgravity Based on Phase Field Method

Numerical Simulation on Single Taylor Bubble Rising in LBE Using Moving Particle Method

Numerical Analysis of Single Pulsating Bubble Interaction with a Free Surface in a Viscous Compressible Liquid

An Improved Multiphase Moving Particle Semi-implicit Method in Bubble Rising Simulations with Large Density Ratios

Current Achievements on Bubble Dynamics Analysis Using MPS Method

Numerical Study of Bubble Rising and Coalescence Characteristics under Flow Pulsation Based on Particle Method

Numerical Investigation on Bubble Dynamics During Flow Boiling Using Moving Particle Semi-Implicit Method

NUMERICAL INVESTIGATION ON TURBULENCE MODELS IN BUBBLY FLOW USING EULER-LAGRANGE APPROACH

GPU accelerated fast multipole boundary element method for simulation of 3D bubble dynamics in potential flow