6.4 Stokes drag

When a small body moves slowly through a viscous fluid, the resistance it feels is set entirely by viscosity, and takes a simple closed form. This lesson derives the drag on a sphere at low Reynolds number — Stokes drag — and the friction coefficient that the Einstein relation of the next lesson depends on.

Drag in the inertia-free limit

A sphere of radius aa moving at velocity U\mathbf{U} through a fluid of viscosity μ\mu, slowly enough that the Reynolds number Re=ρUa/μ1\mathrm{Re} = \rho U a/\mu \ll 1, feels a drag force

Fdrag  =  6πμaU.\mathbf{F}_\text{drag} \;=\; -6\pi\mu a\,\mathbf{U}.

This is the central numerical result of low-Reynolds-number hydrodynamics. It follows from solving the linear Stokes equation around a no-slip sphere; the factor 6π6\pi is exact for a rigid spherical surface. Three features deserve emphasis. The drag is linear in velocity — double the speed, double the force — unlike the high-speed drag that grows as U2U^2. It is linear in size aa, not in cross-sectional area a2a^2, because at low Reynolds number the disturbance the sphere creates reaches out many radii into the fluid. And it depends on μ\mu but not on the fluid density, since inertia plays no role.

The form of the result can be read off from dimensions alone: the only quantities available are μ\mu, aa, and UU, and the single combination with units of force is μaU\mu a U. Solving the flow merely fixes the dimensionless prefactor at 6π6\pi.

Fdrag = 6πμaU = 18.85

The streamlines around a Stokes-flow sphere are *fore-aft symmetric* — the perturbation to the uniform flow is the same upstream and downstream. The total drag, integrated around the sphere, is exactly F = 6πμaU. Doubling the radius doubles the drag; doubling the velocity doubles the drag — a linear-response regime entirely controlled by viscosity.

The streamlines around the sphere are perfectly fore-aft symmetric — a hallmark of Stokes flow, where the time-reversible equation cannot tell upstream from downstream. There is no wake, no separation, and no pressure recovery asymmetry; the entire force comes from viscous shear and the gentle pressure field draped over the sphere.

The friction coefficient and terminal velocity

Writing the drag as Fdrag=γU\mathbf{F}_\text{drag} = -\gamma\,\mathbf{U} defines the friction coefficient

γ  =  6πμa,[γ]=Ns/m.\gamma \;=\; 6\pi\mu a, \qquad [\gamma] = \text{N}\cdot\text{s/m}.

A particle pulled by a steady external force FF — gravity, say — accelerates only until drag balances the force, then coasts at the terminal velocity Uterm=F/γU_\text{term} = F/\gamma. For a sphere settling under gravity in a fluid, balancing buoyancy-corrected weight against Stokes drag gives a settling speed proportional to a2a^2: larger particles fall far faster, which is why fine silt stays suspended for days while sand drops in seconds.

This same friction coefficient was the lever in Millikan’s oil-drop experiment, where the terminal velocity of a falling droplet fixed its radius and hence — through the balancing electric field — the elementary charge. And it is exactly the γ\gamma that ties drag to diffusion in the next lesson: the same molecular collisions that resist a particle’s steady motion also drive its random jitter, and the Einstein relation makes that connection quantitative.