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

Abstract

We present a scalable dissipative particle dynamics simulation code, fully implemented on the Graphics Processing Units (GPUs) using a hybrid CUDA/MPI programming model, which achieves 10-30 times speedup on a single GPU over 16 CPU cores and almost linear weak scaling across a thousand nodes. A unified framework is developed within which the efficient generation of the neighbor list and maintaining particle data locality are addressed. Our algorithm generates strictly ordered neighbor lists in parallel, while the construction is deterministic and makes no use of atomic operations or sorting. Such neighbor list leads to optimal data loading efficiency when combined with a two-level particle reordering scheme. A faster in situ generation scheme for Gaussian random numbers is proposed using precomputed binary signatures. We designed custom transcendental functions that are fast and accurate for evaluating the pairwise interaction. The correctness and accuracy of the code is verified through a set of test cases simulating Poiseuille flow and spontaneous vesicle formation. Computer benchmarks demonstrate the speedup of our implementation over the CPU implementation as well as strong and weak scalability. A large-scale simulation of spontaneous vesicle formation consisting of 128 million particles was conducted to further illustrate the practicality of our code in real-world applications.

Cited by in corpus (12)

• Strong scaling of general-purpose molecular dynamics simulations on GPUs
• Multiscale Universal Interface: A Concurrent Framework for Coupling Heterogeneous Solvers
• A dissipative particle dynamics method for arbitrarily complex geometries
• OpenRBC: A Fast Simulator of Red Blood Cells at Protein Resolution
• TorchMD-Net 2.0: Fast Neural Network Potentials for Molecular Simulations
• GPU-accelerated Red Blood Cells Simulations with Transport Dissipative Particle Dynamics
• A GPU-accelerated package for simulation of flow in nanoporous source rocks with many-body dissipative particle dynamics
• A High-Throughput Solver for Marginalized Graph Kernels on GPU
• Mirheo: High-Performance Mesoscale Simulations for Microfluidics
• Accelerating Particle-in-Cell Monte Carlo Simulations with MPI, OpenMP/OpenACC and Asynchronous Multi-GPU Programming
• A portable code for dissipative particle dynamics (DPD) simulations with additional specific interactions
• Mesoscopic modeling of heptane: A surface tension calculation