Velocity Verlet algorithm for dissipative-particle-dynamics-based models of suspensions

A velocity Verlet algorithm for velocity dependent forces is described for modeling a suspension of rigid body inclusions. The rigid body motion is de termined from the quaternion-based scheme of Omelyan [Comput. Phys. 12, 97 (1998)]. An iterative method to determine angular velocity in a self-consis tent fashion for this quaternion-based algorithm is presented. This method is tested for the case of liquid water. We also describe a method for evalu ating the stress tensor for a system of rigid bodies that is consistent wit h the velocity Verlet alogorithm. Results are compared to the constraint-ba sed rattle algorithm of Anderson [J. Comput. Phys. 52, 24 (1993)]. [S1063-6 51X(99) 13203-3].