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.

Parallelization with load balancing of the weather scheme WSM7 for heterogeneous CPU-GPU platforms

A Parallel Simulator for Massive Reservoir Models Utilizing Distributed-Memory Parallel Systems

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

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

Optimizing the Weather Research and Forecasting Model with OpenMP Offload and Codee

Accelerating the 3D Euler Atmospheric Solver Through Heterogeneous CPU-GPU Platforms

AutoWM: a Novel Domain-Specific Tool for Universal Multi-/Many-core Accelerations of the WRF Cloud Microphysics

Heterogeneous acceleration algorithms for shallow cumulus convection scheme over GPU clusters

Parallel SnowModel (v1.0): a parallel implementation of a distributed snow-evolution modeling system (SnowModel)

Accelerating Lagrangian transport simulations on graphics processing units: performance optimizations of Massive-Parallel Trajectory Calculations (MPTRAC) v2.6

Solving the Global Atmospheric Equations Through Heterogeneous Reconfigurable Platforms.

A Peta-Scalable CPU-GPU Algorithm for Global Atmospheric Simulations

Efficiency and scalability of fully-resolved fluid-particle simulations on heterogeneous CPU-GPU architectures

A Flexible Patch-Based Lattice Boltzmann Parallelization Approach for Heterogeneous GPU-CPU Clusters

Refactoring and Optimizing WRF Model on Sunway TaihuLight

Unleashing the Performance Potential of CPU-GPU Platforms for the 3D Atmospheric Euler Solver.

Accelerating atmospheric physics parameterizations using graphics processing units

GPU acceleration of MPAS microphysics WSM6 using OpenACC directives: Performance and verification

Model Coupling between the Weather Research and Forecasting Model and the DPRI Large Eddy Simulator for Urban Flows on GPU-accelerated Multicore Systems

Towards a GPU-Parallelization of the neXtSIM-DG Dynamical Core

Solving global shallow water equations on heterogeneous supercomputers