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 and optimization of <scp>GPU‐based</scp> parallel <scp>one‐step</scp> leapfrog <scp>ADI‐FDTD</scp> for far‐field scattering problems

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

A Subgridding Scheme Using Hybrid One-Step Leapfrog ADI-FDTD and FDTD Methods

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

Novel Parallel Meshing Technique Implemented into Parallel FDTD Simulation Using Function Language for Solving Electrically Large and Complex EMC Problems

A hybrid one‐step leapfrog ADI‐FDTD and subgridding approach based on a heterogeneous computing platform

Parallel Shooting and Bouncing Ray Method on Gpu Clusters for Analysis of Electromagnetic Scattering

FROM CPU TO GPU: GPU-BASED ELECTROMAGNETIC COMPUTING (GPUECO)

Mapping the Sbr and Tw-Ildcs to Heterogeneous Cpu-Gpu Architecture for Fast Computation of Electromagnetic Scattering

Parallel Transient Stability-Constrained Optimal Power Flow Using GPU as Coprocessor.

Parallel Fast Algorithm for Large-Scale Electromagnetic Scattering

EM Scattering at Rough Surface by Using Parallel FDTD Algorithm

Implementation of the Fdtd Method Based on Lorentz-Drude Dispersive Model on Gpu for Plasmonics Applications

Solving Electromagnetic Scattering Problems with Tens of Billions of Unknowns Using GPU Accelerated Massively Parallel MLFMA

A fine-grained parallel EMTP algorithm compatible to graphic processing units

Rapid simulation of elastic problems based on GPU

Generalized GPU Acceleration for Applications Employing Finite-Volume Methods.

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

An adaptive finite element multigrid solver using GPU acceleration

Massive parallelization and performance enhancement of an immersed boundary method based unsteady flow solver

Accelerating wave-propagation algorithms with adaptive mesh refinement using the Graphics Processing Unit (GPU)