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.

A streaming multi-GPU implementation of image simulation algorithms for scanning transmission electron microscopy

A Fast Algorithm for Scanning Transmission Electron Microscopy (STEM) Imaging and 4D-STEM Diffraction Simulations

A fast image simulation algorithm for scanning transmission electron microscopy

ToTEM: A Software for Fast TEM Image Simulation

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

ATOM 1.0:A GPU Powered Package for Electron Tomography Reconstruction

SuperDragon: A Heterogeneous Parallel System for Accelerating 3D Reconstruction of Cryo-Electron Microscopy Images

Real-time simulations of ADF STEM probe position-integrated scattering cross-sections for single element fcc crystals in zone axis orientation using a densely connected neural network

Single-Particle 3d Reconstruction From Cryo-Electron Microscopy Images On Gpu

GPU Acceleration of Monte Carlo Simulations for Polarized Photon Scattering in Anisotropic Turbid Media.

Graphics processing unit (GPU)-accelerated particle filter framework for positron emission tomography image reconstruction.

GPU-enabled FREALIGN: accelerating single particle 3D reconstruction and refinement in Fourier space on graphics processors.

Parallelizing cryo-EM 3D reconstruction on GPU cluster with a partitioned and streamed model

PCIe-based FPGA-GPU heterogeneous computation for real-time multi-emitter fitting in super-resolution localization microscopy

Realizing smart scanning transmission electron microscopy using high performance computing

GPU Accelerated 3D Tomographic Reconstruction and Visualization From Noisy Electron Microscopy Tilt-Series

CuQ-RTM: A CUDA-based Code Package for Stable and Efficient Q-compensated Reverse Time Migration

Iterative Reconstruction of Micro Computed Tomography Scans Using Multiple Heterogeneous GPUs

PhoenixMR: A GPU-based MRI simulation framework with runtime-dynamic code execution

Low cost, high performance GPU computing solution for atomic resolution cryoEM single-particle reconstruction.

Petascale XCT: 3D Image Reconstruction with Hierarchical Communications on Multi-GPU Nodes