Stokes Flow

Stokes flow (also termed creeping flow or low-Reynolds-number hydrodynamics) is the asymptotic regime of continuum fluid mechanics in which viscous shear forces dominate over fluid inertia. Named after the Anglo-Irish physicist and mathematician Sir George Gabriel Stokes, who formulated the governing theory in his 1851 monograph on the motion of pendulums in viscous media, Stokes flow emerges whenever the Reynolds number is much smaller than unity (Re = ρUL/μ ≪ 1). This condition is satisfied when characteristic length scales L are microscopic, flow velocities U are very slow, or the dynamic viscosity μ of the medium is exceptionally high.

Because the non-linear advective acceleration term (u · ∇u) of the Navier-Stokes equations becomes negligible in the Re → 0 limit, the governing equations reduce to a linear, instantaneous, time-reversible partial differential system. While planetary and synoptic flows in geophysical fluid dynamics operate at Re > 10⁹ (with the notable exception of mantle convection, where Re ≈ 10⁻²⁰), Stokes flow governs the physical world of microbiology, cell biology, aerosol sedimentation, colloidal suspension rheology, and lab-on-a-chip microfluidics.

Consider the incompressible Navier-Stokes equations scaled by a characteristic length L, velocity U, advective time scale L/U, and viscous pressure scale μU/L:

Taking the creeping-flow limit Re = ρUL/μ → 0 eliminates both the Eulerian local acceleration and the quadratic advective acceleration on the left-hand side. In dimensional variables, the velocity field u(x, t) and pressure field p(x, t) satisfy the linear Stokes equations:

where f represents external body forces per unit volume (such as gravity ρg or electrostatic forces).

Biharmonic Structure and Harmonic Pressure

Taking the divergence of the Stokes momentum equation and applying the incompressibility condition ∇ · u = 0 shows that in the absence of divergent body forces (∇ · f = 0), the pressure field is harmonic and satisfies Laplace's equation:

Taking the curl of the Stokes momentum equation eliminates pressure and shows that the vorticity vector ω = ∇ × u also satisfies Laplace's equation (∇²ω = 0)—the steady diffusive limit of the vorticity equation when advective transport and vortex stretching vanish. Applying the curl operator a second time (or introducing a two-dimensional or axisymmetric Stokes streamfunction ψ) demonstrates that the velocity field and streamfunction satisfy the fourth-order biharmonic equation:

Linearity, Superposition, and the Oseen Tensor (Stokeslet)

Because the Stokes equations are strictly linear, complex low-Reynolds-number flows can be constructed by linear superposition of fundamental singularity solutions. The Green's function of the Stokes equations—representing the velocity field u(r) generated at position r by a concentrated point force F acting at the origin in an unbounded fluid—is the Stokeslet, governed by the Oseen tensor J(r):

where r = |r|, r̂ = r/r is the unit radial vector, and I is the identity tensor. Crucially, velocity disturbances in Stokes flow decay algebraically as 1/r—an exceptionally long-range hydrodynamic interaction that couples the motion of suspended macromolecules, polymers, and swimming cells across distances many times their own radii.

Viscous Drag on a Rigid Sphere

In 1851, Stokes solved the axisymmetric biharmonic streamfunction equation ∇⁴ψ = 0 around a rigid sphere of radius R translating at constant speed U through an unbounded quiescent fluid, subject to the no-slip boundary condition at the sphere surface (r = R). In spherical coordinates (r, θ), the velocity components and pressure perturbation are:

Integrating the normal pressure stress (−p cos θ) and tangential viscous shear stress (τ_rθ sin θ) over the spherical surface S = 4πR² reveals that normal pressure stress contributes exactly one-third (2πμRU) of the total resistance—often called form drag—while tangential skin friction contributes the remaining two-thirds (4πμRU). Summing both components yields Stokes' law for viscous drag:

Expressed in terms of the dimensionless drag coefficient C_D = F_D / (½ρU² · πR²) and the diameter-based Reynolds number Re_D = 2ρRU/μ, Stokes' law takes the hyperbolic form C_D = 24 / Re_D. Notably, in two spatial dimensions, no steady creeping-flow solution exists for uniform flow past an infinite cylinder that satisfies both the no-slip condition at the cylinder wall and uniform velocity at infinity (Stokes' paradox, resolved by Carl Wilhelm Oseen in 1910 by retaining linearized far-field convective inertia).

Terminal Sedimentation Velocity and Brownian Diffusion

When a spherical particle of density ρ_p and radius R settles under gravity g (or centrifugal acceleration ω²r in an analytical ultracentrifuge) through a fluid of density ρ_f and viscosity μ, it accelerates until the downward net gravitational weight balances the upward Archimedes buoyancy force and Stokes drag:

This quadratic dependence on particle radius (U_sed ∝ R²) explains why cloud droplets (R ≈ 10 μm, U_sed ≈ 1.2 cm s⁻¹) remain suspended by weak atmospheric updrafts, why clinical hematology uses the erythrocyte sedimentation rate (ESR) to detect inflammation-induced rouleaux aggregation (which increases effective R), and how Svedberg sedimentation coefficients (1 S = 10⁻¹³ s) separate ribosomal subunits (e.g., 30S and 50S prokaryotic subunits) during ultracentrifugation.

For sub-micrometer particles and globular proteins, thermal agitation from surrounding solvent molecules drives Brownian motion. In 1905, Albert Einstein and William Sutherland equated the hydrodynamic friction coefficient ζ = F_D/U = 6πμR from Stokes' law to thermal kinetic energy k_BT, yielding the Stokes–Einstein relation for the translational diffusion coefficient D:

Purcell's Scallop Theorem

Because time t does not appear explicitly in the steady Stokes equations (−∇p + μ∇²u = 0), the velocity field responds instantaneously to boundary boundary motions and obeys kinematic time-reversibility: reversing the velocity of the boundaries (u_boundary → −u_boundary) reverses the velocity field everywhere in the fluid (u → −u) and reverses the sign of the pressure perturbation (p − p_∞ → −(p − p_∞)) along the exact same streamlines. Consequently, the trajectory of a swimmer in Stokes flow depends solely on the sequence of geometric configurations (shapes) it assumes, and is completely independent of the rate (fast or slow) at which it moves through that sequence.

In his classic 1977 lecture "Life at Low Reynolds Number", physicist Edward M. Purcell illustrated this constraint with the Scallop Theorem: a single-hinged reciprocal swimmer (a "scallop" with only one internal degree of freedom that opens its shell slowly and snaps it shut rapidly) cannot achieve any net translational displacement at Re ≪ 1. Whatever tiny distance the scallop advances while snapping shut is exactly erased when it reopens its shell along the time-reversed geometric path.

Breaking Time-Reversal Symmetry in Biological Swimmers

To achieve net propulsion in Stokes flow, biological microorganisms and synthetic microrobots must execute non-reciprocal deformations that enclose a finite area in configuration space (requiring at least two independent mechanical degrees of freedom):

1. Rotary Helical Flagella in Prokaryotes: Peritrichous bacteria such as Escherichia coli and Salmonella enterica propel themselves using rigid, left-handed helical flagellar filaments (~20 nm diameter, 5–10 μm long) driven at their base by a reversible rotary molecular motor powered by the transmembrane proton motive force (or Na⁺ gradient in marine Vibrio species) at 100–300 Hz. Because rotation of a chiral helix in a single sense is non-reciprocal under time reversal, resistive force theory (Gray & Hancock, 1955) shows that the anisotropy between the normal drag coefficient (ζ_⊥) and tangential drag coefficient (ζ_∥, where ζ_⊥ / ζ_∥ ≈ 1.8–2.0 for a slender filament) converts rotary torque directly into axial thrust along the helix axis.
2. Flexible Flagellar Bending Waves and Metachronal Cilia in Eukaryotes: Eukaryotic cells—including human spermatozoa, the biflagellate green alga Chlamydomonas reinhardtii, and the ciliated protozoan Paramecium—utilize the 9+2 microtubule axoneme powered by ATP-dependent dynein motor proteins. Rather than rigid rotation, dynein sliding generates Either a unidirectional propagating sinusoidal wave from base to tip (which breaks time-reversal symmetry by fixing a wave propagation vector k) or an asymmetric two-phase ciliary beat cycle consisting of an extended, high-drag effective (power) stroke followed by a flexed, low-drag recovery stroke sweeping close to the cell surface. In human respiratory epithelia, coordinated metachronal waves of cilia propel the periciliary mucus layer upward at ~4–10 mm min⁻¹ to clear inhaled pathogens and particulates.

Why can't a fish-like flapping motion work at the scale of a bacterium?

A fish swims at Re ≈ 10⁴–10⁶ by shedding momentum-carrying vortices into its wake and coasting on its own inertia, governed by the Euler fluid equations and Bernoulli's principle. If a bacterium of length L = 2 μm stops rotating its flagella at initial speed U₀ = 30 μm s⁻¹ in water, its coasting deceleration distance is d_coast ≈ ρ_pU₀R² / (9μ) ≈ 10⁻¹² m (0.01 Å—less than a tenth of the diameter of a hydrogen atom), stopping within τ ≈ 0.3 microseconds. Without inertia, a rigid reciprocal flap merely pushes fluid forward on the upstroke and pulls the exact same fluid backward on the downstroke.

What is Stokes' paradox in two dimensions?

In three dimensions, the velocity disturbance around a translating sphere decays as 1/r, allowing the Stokes equations to satisfy both the no-slip boundary condition (u = 0 at r = R) and a uniform freestream velocity (u → U_∞ as r → ∞). In two dimensions (flow past an infinite circular cylinder), the radial solutions of the biharmonic equation contain a logarithmic term (ln r) that diverges as r → ∞. Carl Wilhelm Oseen resolved this in 1910 by showing that at distances r ≳ L/Re (the Oseen length), convective inertia (ρU_∞ · ∇u) becomes comparable to viscous diffusion (μ∇²u) and must be retained in the far field, yielding a drag per unit cylinder length of F'_D ≈ 4πμU / [½ − γ_E − ln(Re / 8)], where γ_E ≈ 0.5772 is the Euler–Mascheroni constant.

How does Stokes flow apply to microfluidic lab-on-a-chip devices?

Microfluidic channels typically have hydraulic diameters of L ≈ 10–100 μm and flow speeds of U ≈ 0.1–10 mm s⁻¹, yielding Re ≈ 10⁻³–1. Because the flow is strictly laminar and governed by the linear Stokes equations (and Hagen–Poiseuille flow inside straight conduits), two parallel fluid streams brought together in a Y-junction do not mix by turbulent eddies; instead, they flow side by side and exchange solutes solely through slow molecular diffusion (D = k_BT/(6πμR)) across the interface. Rapid microfluidic mixing therefore requires chaotic advection via staggered herringbone grooves that fold streamlines in three dimensions.

1. Stokes, G. G. (1851). "On the effect of the internal friction of fluids on the motion of pendulums." Transactions of the Cambridge Philosophical Society, 9, 8–106. doi:10.1017/CBO9780511702266.002
2. Purcell, E. M. (1977). "Life at low Reynolds number." American Journal of Physics, 45(1), 3–11. doi:10.1119/1.10903
3. Hapel, J., & Brenner, H. (1983). Low Reynolds Number Hydrodynamics: with special applications to particulate media. Martinus Nijhoff Publishers, Springer. doi:10.1007/978-94-009-8352-6
4. Lauga, E., & Powers, T. R. (2009). "The hydrodynamics of swimming microorganisms." Reports on Progress in Physics, 72(9), 096601. doi:10.1088/0034-4885/72/9/096601