Rotational and Translational Diffusion of Confined Rigid Bodies by Steven Delong, Florencio Balboa, Blaise Delmotte, Brennan Sprinkle and Aleksandar Donev ([email protected]) Courant Institute of Mathematical Sciences.
This package contains several python codes to run simulations of rigid bodies made out of rigidly connected blobs, and confined above a a single wall (floor). These codes can compute the mobility of complex shape objects, solve mobility or resistance problems for suspensions of many bodies or run deterministic or stochastic dynamic simulations.
For the theory behind the numerical methods consult the references:
-
Brownian Dynamics of Confined Rigid Bodies, S. Delong, F. Balboa Usabiaga, and A. Donev, The Journal of Chemical Physics, 143, 144107 (2015). DOI arXiv
-
Hydrodynamics of suspensions of passive and active rigid particles: a rigid multiblob approach F. Balboa Usabiaga, B. Kallemov, B. Delmotte, A. Pal Singh Bhalla, B. E. Griffith, and A. Donev, Communications in Applied Mathematics and Computational Science, 11, 217 (2016). DOI arXiv
-
Brownian dynamics of condined suspensions of active microrollers, F. Balboa Usabiaga, B. Delmotte and A. Donev, The Journal of Chemical Physics, 146, 134104 (2017). DOI arXiv
-
Large Scale Brownian Dynamics of Confined Suspensions of Rigid Particles, B. Sprinkle, F. Balboa Usabiaga, N. Patankar and A. Donev, The Journal of Chemical Physics, 147, 244103 (2017) DOI arXiv.
Several example scripts for simulating immersed rigid bodies near a single wall are present in subfolders.
For usage see doc/README.md.
- doc/: documentation.
- body/: it contains a class to handle a single rigid body.
- boomerang/: older stochastic example from [1], see documentation
boomerang/README.md. - sphere/: the folder contains an example to simulate a sphere
whose center of mass is displaced from the geometric center
(i.e., gravity generates a torque), sedimented near a no-slip wall
in the presence of gravity, as described in Section IV.C in [1].
Unlike the boomerang example this code does not use a rigid
multiblob model of the sphere but rather uses the best known
(semi)analytical approximations to the sphere mobility.
See documentation
doc/boomerang.txt. - many_bodyMCMC/: Markov Chain Monte Carlo code for rigid bodies.
- mobility/: it has functions to compute the blob mobility matrix M and the product Mf using CPUs or GPUs, see [2].
- multi_bodies/: codes to run many-body simulations, based on [3] (minimally-resolved active rollers) and primarily on [4] (general many-particle case).
- quaternion_integrator/: it has a small class to handle quaternions and the schemes to integrate the equations of motion, see [1] and [4].
- stochastic_forcing/: it contains functions to compute the product M^{1/2}z necessary to perform Brownian simulations, see [3] and [4].
- utils.py: this file has some general functions that would be useful for general rigid bodies (mostly for analyzing and reading trajectory data and for logging).