Stokes settling velocity Created 2026-09-28 Updated 2026-10-14
Figure 1.
Sphere settling in water: terminal speed grows as radius squared while Reynolds number grows as radius cubed
. Details.
Take downward as positive and let gravity have magnitude . For a rigid sphere of density in otherwise stationary fluid of density , its weight minus buoyancy is . At terminal velocity, this balances the Stokes drag law:
This signed downward settling velocity is negative for a buoyantly rising sphere and zero for equal densities. For fixed density contrast and dynamic viscosity, the submerged weight grows as while linear drag at fixed speed grows as ; their balance gives the settling-velocity radius scaling . Meanwhile,
Thus doubling the radius multiplies the predicted Stokes settling speed by four and the particle Reynolds number by eight. This makes the formula useful for separating particles by size, while also providing its own check against extrapolation to large particles.
For a computed water example, take , , , , and . These are specified rounded inputs rather than a claim of exact water properties at every temperature. They give
The diameter is , and both the net downward weight and the upward drag force are . At this terminal velocity, falling takes . At the prediction is with ; at it is with , which is no longer a convincing small-Reynolds number case. The figure labels these as predictions of the Stokes settling velocity formula and shows where its inertia check fails.
The calculation describes the eventual terminal velocity, not the release transient governed by particle acceleration and possibly the Unsteady Stokes equation. Walls, hydrodynamic interactions, departures from a rigid sphere, or a different surface condition require a different resistance. For sufficiently small particles, Brownian motion can obscure individual settling trajectories; the formula still specifies the gravitational drift within its assumptions. The balancing argument already appears in article 43 of [ref].
Stokes's law 2026-10-14
For an isolated rigid sphere of radius moving at velocity relative to a distant uniform Newtonian fluid, the drag force exerted by the fluid is
The assumptions are steady incompressible flow, constant dynamic viscosity , an unbounded surrounding fluid, the no-slip boundary condition, and a small diameter-based particle Reynolds number, . Integrating the surface traction in translating-sphere Stokes flow gives the factor : the disturbance pressure contributes to the drag magnitude and the viscous shear stress contributes . See Richard Fitzpatrick, Axisymmetric Stokes Flow Around a Solid Sphere, equations 10.120–10.129.
The result is linear drag: doubling doubles the drag force, and doubling at fixed also doubles it. The required power is , all supplied to viscous dissipation. Density enters the check of the particle Reynolds number, although it does not appear in the leading drag formula.
Nearby walls or other particles change the hydrodynamic resistance matrix and invalidate the isolated-sphere coefficient, even at very small Reynolds number. Finite Reynolds number introduces inertia corrections; the Oseen approximation resolves the distant region where the leading Stokes flow approximation ceases to be uniform. A clean gas bubble has tangential slip and instead obeys clean-bubble Stokes drag, with coefficient in its spherical, negligible-internal-viscosity limit. See H. K. Moffatt, Six Lectures on General Fluid Dynamics, section 3.4 for the sphere solution and its small-Reynolds number limitation. The continuum mechanics and boundary assumptions must also remain valid as the particle becomes very small.
Balancing this Stokes drag law against gravity and buoyancy gives the Stokes settling velocity of an isolated rigid sphere.
Consider two widely separated rigid spheres of equal radius , with separation , in Stokes flow. Apply a force to the first sphere and let the second be force-free and torque-free. The second sphere samples the rate-of-strain tensor of the first sphere's Stokeslet, . Its leading disturbance is the stresslet of a sphere in a uniform straining Stokes flow. Evaluating that disturbance back at the first sphere gives
The hydrodynamic mobility matrix decreases at fixed force, consistently with the increase of resistance caused by an additional rigid inclusion.
Surface energy 2026-10-14
For a liquid interface with constant surface tension , the mechanical surface energy is the energy associated with its surface area :
The additive constant does not affect the work needed to change surface area. For several interfaces with the same surface tension, means their total surface area; the two faces of a thin soap film must both be counted. [ref]
With assumed , creating of one interface requires of reversible work. In the sliding-rim example illustrated under surface tension, moving a rim outwards by creates across the two faces and requires . Independently, multiplying its holding force by the displacement gives the same work.
This makes the preference for a round isolated droplet precise. At fixed volume , the three-dimensional isoperimetric inequality gives
with equality for a sphere. Thus a sphere minimizes surface energy among sufficiently regular closed shapes enclosing the same volume, when surface tension is uniform and other energy contributions are negligible. This is an application of the isoperimetric inequality, rather than a proof of it. Walls, wetting, and gravity can change the shape by adding constraints or other energy contributions. [ref]
The differential rule integrates directly to because surface tension is constant. At fixed liquid volume, decreasing surface area can drive viscous fluid flow; the released surface energy can supply viscous dissipation. This mechanical bookkeeping is useful for predicting which deformations can release energy. If surface tension depends on temperature or composition, distinguish this mechanical area-work relation from the complete interfacial thermodynamic free energy and internal energy.
Surface area is the two-dimensional measure induced on the unit sphere in three dimensions. It is invariant under orthogonal rotations. It is not the restriction of ambient three-dimensional Lebesgue measure, under which the sphere has measure zero. Normalizing by the sphere area gives a probability measure.
Cover the sphere by complements of its two poles. Each open set is contractible, and their intersection deformation retracts to the equatorial sphere. The Poincaré lemma and homotopy invariance of de Rham cohomology make the two neighboring groups in the Mayer--Vietoris sequence for de Rham cohomology zero, so its connecting map is an isomorphism. Its inverse takes a top-degree form to the difference of its local primitives on the equator, up to the connecting-map sign convention.
Tangent plane 2026-10-14
For an embedded regular surface and , the affine tangent plane is
It is a set of points through , whereas the tangent space is a vector space of displacements with zero vector. Some texts also call that two-dimensional tangent space the tangent plane; here the translation by distinguishes the two uses. For a regular parametric surface with ,
where the partial derivatives are evaluated at . Regularity means that these tangent vectors are linearly independent, equivalently . Their cross product is a normal vector. [ref]
For an implicit regular surface , require a smooth function with . Then
Differentiating proves that every surface velocity lies in the displayed kernel of a linear map; the implicit function theorem gives a local graph with two independent coordinates, proving equality. The gradient is the normal vector with respect to the ambient dot product. [ref]
For the sphere , , and the affine tangent plane has equation . At the unit sphere's north pole , compare the plane with the displacement tangent space . The regular parametric surface , for , gives and , confirming the same answer.
Figure 1.
The sphere's tangent plane passes through the base point; its tangent space contains displacements
. Details.
This analytical illustration shows a unit sphere and a vertical cross-section. The smooth curves and give the two displayed tangent vectors at . The straight point for any nonzero is outside the sphere when , since . A tangent plane therefore captures first-order motion, not an exact patch of the regular surface.
The same idea makes local approximation useful. For the graph of a function with a smooth function, the affine tangent plane is . [ref] With continuous second partial derivatives, the vertical error is in Big O notation. On the unit sphere near its north pole, : the tangent plane loses quadratic bending while preserving all first-order velocities.
A failed regularity test needs interpretation. Replacing the sphere equation by makes everywhere on the same perfectly regular sphere, so this defining smooth function cannot recover its tangent plane. Similarly, the usual spherical parametric surface loses matrix rank at the poles, while the graph parametrization above remains regular there. A vanishing gradient or degenerate parametrization invalidates that calculation; it alone does not prove that the underlying geometry has no tangent plane.