Stokesian dynamics

Stokesian dynamics
is a solution technique for the Langevin equation, which is the relevant form of Newton's 2nd law for a Brownian particle. The method treats the suspended particles in a discrete sense while the continuum approximation remains valid for the surrounding fluid, i.e., the suspended particles are generally assumed to be significantly larger than the molecules of the solvent. The particles then interact through hydrodynamic forces transmitted via the continuum fluid, and when the particle Reynolds number is small, these forces are determined through the linear Stokes equations (hence the name of the method). In addition, the method can also resolve non-hydrodynamic forces, such as Brownian forces, arising from the fluctuating motion of the fluid, and interparticle or external forces. Stokesian Dynamics can thus be applied to a variety of problems, including sedimentation, diffusion and rheology, and it aims to provide the same level of understanding for multiphase particulate systems as molecular dynamics does for statistical properties of matter. For formula_1 rigid particles of radius formula_2 suspended in an incompressible Newtonian fluid of viscosity formula_3 and density formula_4, the motion of the fluid is governed by the Navier–Stokes equations, while the motion of the particles is described by the coupled equation of motion:
In the above equation formula_6 is the particle translational/rotational velocity 
vector of dimension 6N. formula_7 is the hydrodynamic force, i.e., force exerted by the fluid on the particle due to relative motion between them. formula_8 is the stochastic Brownian force due to thermal motion of fluid particles. formula_9 is the deterministic nonhydrodynamic force, which may be almost any form of interparticle or external force, e.g. electrostatic repulsion between like charged particles. Brownian dynamics is one of the popular techniques of solving the Langevin equation, but the hydrodynamic interaction in Brownian dynamics is highly simplified and normally includes only the isolated body resistance. On the other hand, Stokesian dynamics includes the many body hydrodynamic interactions. Hydrodynamic interaction is very important for non-equilibrium suspensions, like a sheared suspension, where it plays a vital role in its microstructure and hence its properties. Stokesian dynamics is used primarily for non-equilibrium suspensions where it has been shown to provide results which agree with experiments.
Hydrodynamic interaction.
When the motion on the particle scale is such that the particle Reynolds number is small, the hydrodynamic force exerted on the particles in a suspension undergoing a bulk linear shear flow is:
Here, formula_11 is the velocity of the bulk shear flow evaluated at the particle 
center, formula_12 is the symmetric part of the velocity-gradient tensor; formula_13 and formula_14 are the configuration-dependent resistance matrices that give the hydrodynamic force/torque on the particles due to their motion relative to the fluid (formula_13) and due to the imposed shear flow (formula_14). Note that the subscripts on the matrices indicate the coupling between kinematic (formula_6) and dynamic (formula_18) quantities.
One of the key features of Stokesian dynamics is its handing of the hydrodynamic interactions, which is fairly accurate without being computationally inhibitive (like boundary integral methods) for a large number of particles. Classical Stokesian dynamics requires formula_19 operations where "N" is the number of particles in the system (usually a periodic box). Recent advances have reduced the computational cost to about formula_20
Brownian force.
The stochastic or Brownian force formula_8 arises from the thermal fluctuations in the fluid and is characterized by:
The angle brackets denote an ensemble average, formula_24 is the Boltzmann constant, formula_25 is the absolute temperature and formula_26 is the delta function. The amplitude of the correlation between the Brownian forces at time formula_27 and at time formula_28 results from the fluctuation-dissipation theorem for the N-body system.