Rigid body dynamics approach to Stokesian dynamics simulations of nonspherical particles
No Thumbnail Available
Date
2010-05-07
Authors
Journal Title
Journal ISSN
Volume Title
Publisher
American Institute of Physics
Abstract
We describe an algorithm for performing Stokesian dynamics (SD) simulations of suspensions of arbitrary shape rigid particles with hydrodynamic interactions, modeled as rigid groups of spheres, the hydrodynamic mobility matrix of which is accurately computable by several established schemes for spheres. The algorithm is based on Stokesian rigid body equations of translational and rotational motion, which we have derived by an approach formally analogous to that of Newtonian rigid body dynamics. Particle orientation is represented in terms of Euler parameters (quaternion of rotation). This rigid body SD algorithm (RBSDA) complements recently described constraint SD algorithms [ R. Kutteh, J. Chem. Phys. 119, 9280 (2003) ; R. Kutteh, Phys. Rev. E 69, 011406 (2004) ], over which it offers the same computational advantages in imposing total rigidity that the basic rigid body molecular dynamics (MD) algorithm offers over constraint MD algorithms. We show that SD simulation results generated with the RBSDA, in bounded and unbounded geometries, agree very well with those from experiment and other SD and non-SD methods, and are numerically identical to those from a constraint SD algorithm, HSHAKE. Finally, for completeness we also describe a third (additional to the constraint SD and rigid body SD approaches) more traditional approach for SD simulations of arbitrary shape rigid particles modeled as rigid groups of spheres. © 2010, American Institute of Physics
Description
Keywords
Computational fluid dynamics, Hydrodynamics, Molecular dynamics method, Suspensions, Algorithms, Spheres
Citation
Kutteh, R. (2010). Rigid body dynamics approach to Stokesian dynamics simulations of nonspherical particles. Journal of Chemical Physics, 132(17), 174107. doi:10.1063/1.3358330