We have developed and tested a velocity-Verlet algorithm for a dissipative particle dynamics based model describing the motion of rigid body inclusions. A simple procedure for calculating the stress tensor contribution from the rigid body which is consistent with velocity-Verlet algorithm was given. The velocity-Verlet algorithm for DPD is less sensitive to variation in time step size than the Euler algorithm presented in the original DPD papers, thus significantly improving numerical accuracy at little computational cost. Although the original motivation of the paper was to improve upon the orginal DPD algorithm such that the DPD particles which represent a solvent and the rigid body motion are treated in a self consistent fashion, the numerical techniques presented in this paper should not be limited to DPD.
Nicos S. Martys would like to acknowledge support for this research from the National Institute of Standards and Technolgy Program on High-Performance Construction Materials and Systems.