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 Algorithms for Computing Matrix Function Vector Products with Applications to Exponential Integrators and Fractional Diffusion.

Parallelized Implementation of the Finite Particle Method for Explicit Dynamics in GPU

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

High Performance Matrix Multiplication on General Purpose Graphics Processing Units

GPU-Acceleration of Parallel Unconditionally Stable Group Explicit Finite Difference Method

An adaptive finite element multigrid solver using GPU acceleration

Fast hardware-aware matrix-free algorithms for higher-order finite-element discretized matrix multivector products on distributed systems

Acceleration of Approximate Matrix Multiplications on GPUs

An Efficient Acceleration of Solving Heat and Mass Transfer Equations with the Second Kind Boundary Conditions in Capillary Porous Radially Composite Cylinder Using Programmable Graphics Hardware

Solving the Caputo Fractional Reaction-Diffusion Equation on GPU

Fast hardware-aware matrix-free algorithm for higher-order finite-element discretized matrix multivector products on distributed systems

Finite Element Integration with Quadrature on the GPU

Accelerating Iterative Linear Solvers Using Multiple Graphical Processing Units

Hybrid GPU–CPU Efficient Implementation of a Parallel Numerical Algorithm for Solving the Cauchy Problem for a Nonlinear Differential Riccati Equation of Fractional Variable Order

Performant low-order matrix-free finite element kernels on GPU architectures

A New Sparse Matrix Vector Multiplication GPU Algorithm Designed for Finite Element Problems

GPU Implementation of Conjugate Gradient Method

Accelerating Dissipative Particle Dynamics Simulations on GPUs: Algorithms, Numerics and Applications

ADVANCES IN GRAPHICS PROCESSING UNITS' APPLICATIONS TO THE COMPUTATION OF LARGE-SCALE MECHANICAL PROBLEMS

Generalized Gpu Acceleration For Applications Employing Finite-Volume Methods

Employing Graphics Processing Unit Technology, Alternating Direction Implicit Method And Domain Decomposition To Speed Up The Numerical Diffusion Solver For The Biomedical Engineering Research