Stokesian dynamics^[1] 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 ${\displaystyle N}$ rigid particles of radius ${\displaystyle a}$ suspended in an incompressible Newtonian fluid of viscosity ${\displaystyle \eta }$ and density ${\displaystyle \rho }$, 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:

${\displaystyle \mathbf {m} {\frac {d\mathbf {U} }{dt}}=\mathbf {F} ^{\mathrm {H} }+\mathbf {F} ^{\mathrm {B} }+\mathbf {F} ^{\mathrm {P} }.}$

In the above equation ${\displaystyle \mathbf {U} }$ is the particle translational/rotational velocity vector of dimension 6N. ${\displaystyle \mathbf {F} ^{\mathrm {H} }}$ is the hydrodynamic force, i.e., force exerted by the fluid on the particle due to relative motion between them. ${\displaystyle \mathbf {F} ^{\mathrm {B} }}$ is the stochastic Brownian force due to thermal motion of fluid particles. ${\displaystyle \mathbf {F} ^{\mathrm {P} }}$ 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.^[2]

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:

${\displaystyle \mathbf {F} ^{\mathrm {H} }=-\mathbf {R} _{\mathrm {FU} }(\mathbf {U} -\mathbf {U} ^{\infty })+\mathbf {R} ^{\mathrm {FE} }:\mathbf {E} ^{\infty }.}$

Here, ${\displaystyle \mathbf {U} ^{\infty }}$ is the velocity of the bulk shear flow evaluated at the particle center, ${\displaystyle \mathbf {E} ^{\infty }}$ is the symmetric part of the velocity-gradient tensor; ${\displaystyle \mathbf {R} _{\mathrm {FU} }}$ and ${\displaystyle \mathbf {R} _{\mathrm {FE} }}$ are the configuration-dependent resistance matrices that give the hydrodynamic force/torque on the particles due to their motion relative to the fluid (${\displaystyle \mathbf {R} _{\mathrm {FU} }}$) and due to the imposed shear flow (${\displaystyle \mathbf {R} _{\mathrm {FE} }}$). Note that the subscripts on the matrices indicate the coupling between kinematic (${\displaystyle \mathbf {U} }$) and dynamic (${\displaystyle \mathbf {F} }$) quantities.

One of the key features of Stokesian dynamics is its handling 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 ${\displaystyle O(N^{3})}$ operations where N is the number of particles in the system (usually a periodic box). Recent advances have reduced the computational cost to about ${\displaystyle O(N^{1.25}\,\log N).}$^[3][4]

The stochastic or Brownian force ${\displaystyle \mathbf {F} ^{\mathrm {B} }}$ arises from the thermal fluctuations in the fluid and is characterized by:

${\displaystyle \left\langle \mathbf {F} ^{\mathrm {B} }\right\rangle =0}$

${\displaystyle \left\langle \mathbf {F} ^{\mathrm {B} }(0)\mathbf {F} ^{\mathrm {B} }(t)\right\rangle =2kT\mathbf {R} _{\mathrm {FU} }\delta (t)}$

The angle brackets denote an ensemble average, ${\displaystyle k}$ is the Boltzmann constant, ${\displaystyle T}$ is the absolute temperature and ${\displaystyle \delta (t)}$ is the delta function. The amplitude of the correlation between the Brownian forces at time ${\displaystyle 0}$ and at time ${\displaystyle t}$ results from the fluctuation-dissipation theorem for the N-body system.

• Immersed boundary methods
• Stochastic Eulerian Lagrangian methods