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.

First Principle Calculation on Multi-GPU Machines

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

First-principle nonlocal projector potential calculation on GPU cluster

The Analysis of a Plane Wave Pseudopotential Density Functional Theory Code on a GPU Machine

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

Large Scale Plane Wave Pseudopotential Density Functional Theory Calculations on GPU Clusters

Fast plane wave density functional theory molecular dynamics calculations on multi-GPU machines

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

GPU Tuning for First-Principle Electronic Structure Simulations

GPU_PBTE: an Efficient Solver for Three and Four Phonon Scattering Rates on Graphics Processing Units

Several Optimized Algorithms for Ultra-soft Psedupotential DFT Molecular Dynamics Simulation

Many-core Acceleration of the First-Principles All-Electron Quantum Perturbation Calculations.

Multi-GPU Hybrid Programming Accelerated Three-Dimensional Phase-Field Model in Binary Alloy

GPU Acceleration of Numerical Atomic Orbitals-Based Density Functional Theory Algorithms within the ABACUS package

Gpu-Based Parallel Calculation Method For (1)Molecular Weight Distribution Of Batch Free Radical Polymerization

Efficient molecular dynamics simulations with many-body potentials on graphics processing units

An alternative GPU acceleration for a pseudopotential plane-waves density functional theory code with applications to metallic systems

Accelerating Relativistic Exact-Two-Component Density Functional Theory Calculations with Graphical Processing Units

Parallel Computing of First-principles Based Quantum Transport Simulations

OpenCL-accelerated first-principles calculations of all-electron quantum perturbations on HPC resources

High Performance Computing for First-Principles Kohn-Sham Density Functional Theory Towards Exascale Supercomputers