Parallel implementation of isothermal and isoenergetic Dissipative Particle Dynamics using Shardlow-like splitting algorithms

JP Larentzos and JK Brennan and JD Moore and M Lisal and WD Mattson, COMPUTER PHYSICS COMMUNICATIONS, 185, 1987-1998 (2014).

DOI: 10.1016/j.cpc.2014.03.029

A parallel implementation of the Shardlow splitting algorithm (SSA) for Dissipative Particle Dynamics (DPD) simulations is presented. The isothermal and isoenergetic SSA implementations are compared to the DPD version of the velocity-Verlet integrator in terms of numerical stability and performance. The integrator stability is assessed by monitoring temperature, pressure and total energy for both the standard and ideal DPD fluid models. The SSA requires special consideration due to its recursive nature resulting in more inter-processor communication as compared to traditional DPD integrators. Nevertheless, this work demonstrates that the SSA exhibits stability over longer time steps that justify its regular use in parallel, multi-core applications. For the computer architecture used in this study, a factor of 10-100 speedup is achieved in the overall time-to-solution for isoenergetic DPD simulations and a 15-34 speedup is achieved for the isothermal DPD simulations. 2014 Elsevier B.V. All rights reserved.

Return to Publications page