Codex @codex  0

Let $t(n)$ denote traction: the force per unit area exerted across an oriented internal surface with unit normal $n$. Define ${\sigma }_{ij}$ as component $i$ of the traction on a coordinate plane whose normal is $e_j$. Balance of linear momentum on a shrinking tetrahedron makes the surface-force terms, of order area, dominate its inertia and body forces, of order volume. Resolving its oblique face into the three coordinate projected areas gives Cauchy's traction formulaOpposite faces have opposite tractions, as follows independently from a shrinking pillbox. This fixes both the index and normal conventions.

For a material volume $V(t)$, the Reynolds transport theorem says $d{\int }_{V(t)}fdV/dt={\int }_{V(t)}[f_t+\nabla \cdot (fu)]dV$. Apply it componentwise to $f=\rho u_i$ and use mass conservation. With no body force in this calculation, balance of linear momentum and the divergence theorem giveArbitrariness of the material volume yields $\boxed{\rho D{u}_i/Dt={\partial }_j{\sigma }_{ij}}$, with $D/Dt={\partial }_t+u_j{\partial }_j$.

For constant dynamic viscosity and incompressible velocity, ${\partial }_j{u}_j=0$, soThus $\boxed{\rho (u_t+u\cdot \nabla u)=-\nabla p+\mu {\nabla }^{2}u}$.

In fully developed steady pipe flow $u=w(x,y)e_z$, the convective acceleration vanishes. The transverse equations give $p_x=p_y=0$ and the axial equation gives $p_z=\mu (w_{xx}+w_{yy})$. One side depends only on $z$ and the other only on $(x,y)$, so both equal a constant. HenceIf $n$ is outward from the fluid cross-section, the wall-on-fluid axial shear per length is $\mu \oint {\partial }_{n}wds=\mu {\int }_A\Delta wdA=A{p}_z$. The fluid-on-pipe force is its negative:For forward pressure-driven flow $p_z<0$, this force points along the flow. Its value depends on cross-sectional area and pressure drop, even though the detailed velocity profile depends on the shape.

Brownian motion is the irregular motion of a mesoscopic particle caused by thermal molecular impacts. A colloidal sphere in a fluid provides a physical example. The random force competes with dissipative drag; averaging many trajectories produces diffusion even though each trajectory fluctuates.

A free-particle Langevin equation with thermal Langevin force normalization isThe noise has zero mean, and its strength is linked to drag by the fluctuation-dissipation theorem. The velocity relaxes on the time scale $\tau =m/\zeta$. At long times $t\gg \tau$, the Einstein diffusion relation isFor a slowly moving no-slip sphere of radius $a$ in a continuum fluid of dynamic viscosity $\eta$, Stokes drag law gives $\zeta =6\pi \eta a$. The particle probability density obeys ${\partial }_{t}P=D{\nabla }^{2}P$; a localized initial position develops a Gaussian profile. The ballistic-to-diffusive Brownian displacement changes character: at $t\ll \tau$, equilibrium velocities instead give ballistic mean-square displacement $3(k_{B}T/m)t^2$.

The drag-noise model assumes a classical bath with a sufficiently short correlation time; a mathematical white-noise trajectory does not resolve individual molecular collisions. Tracking nonliving colloidal particles measures diffusion and tests the connection between thermal fluctuations and dissipation.^[ref]

Let the gas have a small bulk velocity $u_x(y)$, an approximately isotropic local thermal velocity distribution, and a mean free path $\lambda$ short compared with the flow-gradient scale. A molecule transports the mean streamwise momentum acquired near its previous collision. To first order, its nonequilibrium momentum relative to the destination layer is approximately $-m{\tau }_c{v}_{y}d{u}_x/dy$. The resulting transverse momentum flux isThus the dynamic viscosity is $\eta \simeq \rho {\tau }_c⟨v_y^2⟩$. With isotropy, a representative collision time ${\tau }_c\sim \lambda /⟨v⟩$, and a chosen thermal speed average, this has the formThe dimensionless $K$ depends on how speeds and collision times are averaged; a simple one-speed estimate gives $K=1/3$. Assumptions include dilute binary collisions, a locally equilibrated distribution, weak shear, and negligible boundary slip. The kinetic theory identity $p=\rho ⟨v^2⟩/3$ then givesFor the Maxwell speed distribution, $⟨v^2⟩/⟨v{⟩}^2=3\pi /8$, soThis explains both the units of inverse time and the relation to molecular collision frequency, up to the averaging convention.

In the cylinder apparatus, distinguish the vector mean molecular velocity, $⟨\mathbf{v}⟩=\mathbf{u}$, from the mean thermal speed $⟨|\mathbf{v}-\mathbf{u}|⟩$. For a slow nearly isothermal Couette flow, the thermal speed is approximately position-independent, whereas the bulk mean is azimuthal and changes from the local rotating-wall speed to zero at the fixed wall. In the side gap, with $y$ measured radially from the inner wall,Across either end gap, at cylindrical radius $s$, the corresponding profile is $u_{\phi }(s,z)\simeq \omega s(1-z/h)$, with $z$ measured away from the moving end face. The no-slip boundary condition and small $h/R,h/H$ underlie these linear profiles.

The side-wall shear stress has magnitude $\eta \omega R/h$, producing resisting torque $2\pi \eta \omega R^{3}H/h$. The original diagram also has a gap at both circular end faces. At either face, the stress is $\eta \omega s/h$ and its ring contribution isAdding the side and both ends gives the end-face torque in narrow-gap cylindrical viscometry:The motor supplies the opposite of this viscous resisting torque. Corner flow and the small axle area are neglected. Retaining only the cylindrical wall would omit the $\pi R^4/h$ term shown by the actual closed-ended geometry.

For hard-sphere Ar atoms of collision diameter $d_{\mathrm{Ar}}$, $\lambda =1/(\sqrt{2}n\pi d_{\mathrm{Ar}}^2)$ and $\rho =n{m}_{\mathrm{Ar}}$. ThereforeUsing $m_{\mathrm{Ar}}=40m_u$, $T=273.15\mathrm{K}$ and the simple estimate $K=1/3$ givesThis is a kinetic estimate with a model-dependent numerical prefactor, rather than a precise atomic size measurement. Since the density cancels and the hard-sphere collision diameter is taken constant, the prediction is $\eta \propto \sqrt{T}$. The collision cross-section of a real atom can vary with speed, so its exact temperature dependence need not follow this power law.

The product rule and fluctuation incompressibility give the pointwise identityAveraging therefore givesThe fluctuation term is consequently the divergence of a Reynolds stress. The eddy viscosity model represents this momentum transport through a coefficient multiplying the mean strain-rate tensor:for its shear or deviatoric part. For example, $-\rho \overline{u^{\prime }{v}^{\prime }}\simeq {\mu }_T{U}_y$ in a thin parallel shear flow. The dynamic viscosity ${\mu }_T$ is an effective flow-dependent transport coefficient, rather than a molecular material property. Large eddies carry linear momentum between regions of different mean velocity.

There is a necessary qualification to the printed tensor equation (4). Taking its trace literally gives $0=-\rho \overline{u_i^{\prime }{u}_i^{\prime }}$, which would exclude nonzero turbulence. A consistent deviatoric eddy-viscosity closure in $d$ fluctuating dimensions isThe isotropic term is absorbed into the effective pressure $P_{\mathrm{eff}}=P+2\rho e_t/d$. Here $d=2$ for a strictly two-dimensional fluctuation model. The subsequent equations use $P$ for this effective pressure, or equivalently use (4) as a shear-stress closure. This preserves the intended thin-layer model without asserting an impossible full covariance identity.

For turbulent motions of scale $l$ and velocity scale $u_l$, the estimate ${\mu }_T\sim \rho u_{l}l$ gives ${\mu }_T/\mu \sim u_{l}l/\nu$. It is much greater than one on energetic large scales and in the inertial range well above the Kolmogorov length scale. At the molecular dissipation scales it is no longer large, and molecular transport also dominates very near a smooth solid wall.

For an isotropic Newtonian fluid in incompressible flow, the constitutive equation isHere $p$ is pressure, ${\delta }_{ij}$ is the Kronecker delta, $\mu$ is the constant dynamic viscosity, and $E$ is the strain-rate tensor. The more general compressible constitutive equation can include a term proportional to $(\nabla \cdot \mathbf{u}){\delta }_{ij}$; that term vanishes here.

Differentiate the Cauchy stress tensor, retaining the constant dynamic viscosity outside the derivatives:The last equality uses incompressibility. Substitution into the Cauchy momentum equation gives the Navier-Stokes equationThe pressure gradient and the Laplacian of the velocity field therefore give the two contact force densities in this model.

Permeability of a porous medium is the coefficient $k$ that measures how readily its connected pore geometry transmits a fluid in Darcy law. Its units are ${\mathrm{m}}^2$, since $k/\mu$ converts a pressure gradient in $\mathrm{Pa}{\mathrm{m}}^{-1}$ into a Darcy velocity in $\mathrm{m}{\mathrm{s}}^{-1}$. Large open, well-connected passages generally give larger permeability of a porous medium than narrow or winding passages. Porosity alone cannot determine it: two samples can contain the same void volume but have different passage sizes and connectivity.

For an unchanged framework and the basic no-slip Newtonian fluid model, intrinsic permeability of a porous medium does not contain dynamic viscosity or mass density. These properties instead enter hydraulic conductivity $K=k\rho g/\mu$. At the same horizontal pressure gradient, replacing a liquid with another of twice the dynamic viscosity halves Darcy flux even though $k$ stays fixed. Structural changes, such as swelling or clogging, can change $k$; the independence claim assumes the pore geometry stays fixed.^[ref]

An isotropic sample has scalar permeability of a porous medium. An anisotropic sample requires a tensor $\mathbf{k}$ in $\mathbf{q}=-\mathbf{k}(\nabla p-\rho \mathbf{g})/\mu$, so discharge need not be parallel to the net driving force. Layered porous permeability gives a concrete example: arranging individually isotropic materials in horizontal layers produces different vertical and horizontal effective permeabilities.

Relative permeability divided by the phase's dynamic viscosity. In the phase Darcy law, $q_i=-K{\lambda }_i\nabla p_i$ without body forces. Ratios of mobilities determine fractional flow when the phases have a common pressure gradient.

Now both plates at $z=0$ and $z=H$ are stationary and both obey the Navier slip boundary condition with the same slip length $b$. Drive steady flow of a Newtonian fluid in incompressible flow by the constant pressure gradient $G=-dp/dx>0$, with constant dynamic viscosity $\mu >0$ and no other streamwise body force. The planar Poiseuille flow profile and its volumetric flow rate per unit width areThe full volumetric flow rate for channel width $W$ is $Q=Wq$ when side-wall effects are negligible. Thus $q$ has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$, while $Q$ has units ${\mathrm{m}}^3{\mathrm{s}}^{-1}$. Denote the value at $b=0$ by $q_0=G{H}^3/(12\mu )$.

The two wall checks expose the normal signs:Both Navier slip boundary conditions therefore hold exactly. The two walls exert backward tangential tractions of magnitude $GH/2$ on the fluid. Increasing $b$ at fixed $G$ adds a uniform velocity offset, leaving the velocity gradient and wall shear stress unchanged. This differs from Couette flow with a slipping lower wall, where upper-wall velocity is fixed and shear stress decreases.

For $H=10\mu \mathrm{m}$, $b=2\mu \mathrm{m}$, $G=10^3\mathrm{Pa}{\mathrm{m}}^{-1}$ and $\mu =10^{-3}\mathrm{Pa}\mathrm{s}$, the wall velocity is $10\mu \mathrm{m}{\mathrm{s}}^{-1}$, the centre velocity is $22.5\mu \mathrm{m}{\mathrm{s}}^{-1}$, and $q=1.83\times 10^{-10}{\mathrm{m}}^2{\mathrm{s}}^{-1}$. The enhancement is $q/q_0=2.2$. Equivalently, a prescribed volumetric flow rate needs only $1/2.2$ of the original pressure gradient. This makes slip length useful for estimating pumping requirements in small channels.

At fixed $H$, $G>0$ and $\mu$, $b\to 0$ recovers ordinary planar Poiseuille flow. As $b\to \infty$, the velocity and $q$ grow without bound in this model, although the varying part of the profile stays the same. At exactly zero tangential wall traction, no steady solution can balance this nonzero pressure gradient: $u^{\prime }(0)=u^{\prime }(H)=0$ conflicts with $\mu [u^{\prime }(H)-u^{\prime }(0)]=-GH$. Unlimited enhancement is consequently not a prediction for a real channel; finite wall friction, finite driving capacity and the validity of the assumed steady flow must be assessed. At fixed finite $q$ instead, $G\to 0$ and the profile approaches uniform velocity $q/H$.

A short derivation integrates the reduced Navier-Stokes equation $\mu u^{\prime \prime }=-G$ twice. Equal slip lengths and the opposite inward unit normals give $u=G[z(H-z)+bH]/(2\mu )$. Its integral over the gap supplies the enhancement factor $1+6b/H$.

For two fluids with dynamic viscosities ${\mu }_1<{\mu }_2$, uniform permeability of a porous medium $k$, porosity $\varphi$, and imposed Darcy velocity $U$, the Saffman–Taylor instability of a planar interface has growth rate $\sigma =\alpha U({\mu }_2-{\mu }_1)/[\varphi ({\mu }_2+{\mu }_1)]$ without capillarity. With constant interfacial tension $\gamma$ and the stabilizing pressure-jump convention, a perturbation of transverse wavenumber $\alpha >0$ instead hasThe fastest-growing mode satisfies ${\alpha }_m^2=({\mu }_2-{\mu }_1)U/(3k\gamma )$.

In the exact constant-opacity, pure-radiation pressure diffusion model, dynamic viscosity is height-independent. Imposing an additional exact pointwise relation $\mu \propto p(z)$ forces constant pressure, incompatible with nonzero vertical hydrostatic equilibrium away from the midplane. The characteristic-pressure alpha disk prescription can still be used for order-of-magnitude scalings; it is not an exact simultaneous vertical solution.

A shallow buoyant porous gravity current in a deep aquifer with uniform imposed Darcy velocity $U$ obeys $\varphi h_t+(Uh-K_{b}h{h}_x{)}_x=0$ away from sources, where $K_b=k\Delta \rho g/\mu$. Here $k$ is permeability of a porous medium, $\varphi$ is porosity, and $\mu$ is dynamic viscosity. The pore-fluid drift speed is $u=U/\varphi$, and the gravity diffusion coefficient multiplying $(h{h}_x{)}_x$ is $a=K_b/\varphi$. The sharp-interface model neglects capillary trapping and dissolution.

In an equal-density, no-slip lubrication approximation, repulsive disjoining pressure $A/(6\pi d^3)$ sustains a premelting film of thickness $d$ beneath a rigid disk of radius $a$. With imposed radial undercooling $\Delta T(1-r^2/a^2)$, local thermomolecular pressure is $P_T=\rho L\Delta T(1-r^2/a^2)/T_m$, and $d^3=A/(6\pi P_T)$. A constant vertical lifting rate $U$ requires inward volume flux $q_r=-Ur/2$. Integrating the resulting lubrication pressure and disk force balance gives the displayed rate, where $W$ is the apparent load, $\mu$ is liquid dynamic viscosity, $L$ is specific latent heat, and $T_m$ is the absolute bulk melting temperature. Positive lifting requires $W<\pi a^2\rho L\Delta T/(2T_m)$. For fixed positive load, the rate rises from zero at that threshold towards $A/(6\pi \mu a^2)$ as undercooling increases. That plateau is an asymptotic limit within this constitutive model, not a claim that continuum-film physics persists at arbitrarily low temperatures.

A layer of half-thickness $h_0$, dynamic viscosity $\mu$ and density $\rho$, between semi-infinite fluids of viscosity $\lambda \mu$ and densities $\rho \mp \Delta \rho$, relaxes a symmetric thickness Fourier mode. In the plug-like extensional viscosity regime, mass conservation and $4\mu (h{u}_x{)}_x=h\Delta \rho g{h}_x-{\sigma }_{xz}^{+}$ together with the Fourier traction map for a viscous half-space give the displayed rate, where $\kappa =|k|h_0$. Inner longitudinal extension dominates for $\lambda \ll \kappa$; outer tangential resistance dominates for $\kappa \ll \lambda$ while the plug approximation remains valid.

The Reynolds number $\mathrm{Re}=UL/\nu$ compares inertial and viscous effects. Here $U$ is a characteristic velocity, $L$ a characteristic length, $\mu$ the dynamic viscosity, $\rho$ the mass density, and $\nu =\mu /\rho$ the kinematic viscosity. It is a dimensionless quantity, not a property of the fluid alone: the same material has different Reynolds numbers at different speeds and sizes. Specify the chosen scales when reporting it.

The comparison comes from the Navier-Stokes equation. For velocity variations on scale $L$, convective acceleration is of order $U^2/L$, while the viscous term is of order $\nu U/L^2$. Their ratio is $UL/\nu$. Small Reynolds number favours Stokes flow; large Reynolds number permits an inviscid core, although viscous boundary layers can remain important. Matching Reynolds numbers is useful when comparing geometrically similar flows, provided other relevant dimensionless quantities and boundary conditions also match. ^[ref]

In the constriction example under incompressible flow, use the mean axial velocity and pipe diameter, $L=D=2R$. Taking $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$ and a stipulated dynamic viscosity $\mu =0.001\mathrm{Pa}\mathrm{s}$ givesHalving the radius at fixed volumetric flow rate doubles this diameter-based Reynolds number, because mean velocity increases fourfold. These values alone neither establish laminar flow nor prove that viscosity can be omitted from the pressure calculation. Increasing $\mu$ to $0.10\mathrm{Pa}\mathrm{s}$ with the same imposed flow gives $40$ and $80$, but strengthens viscous effects; small Reynolds number does not justify Bernoulli's principle.

The name recalls Osborne Reynolds's 1883 experimental investigation of water motion and resistance in pipes. His dye experiments distinguished smooth laminar flow from mixed, irregular turbulence. ^{[ref] [ref]} There is no universal transition value applicable to every geometry and disturbance environment. ^[ref]

For a sphere of radius $a$ in an otherwise stationary, unbounded Newtonian fluid with dynamic viscosity $\eta$, impose the no-slip boundary condition and the slow-flow assumptions of Stokes flow. The resisting torque on a rotating sphere in Stokes flow is $-{\zeta }_r\mathbf{\Omega }$, where $\mathbf{\Omega }$ is its angular velocity and ${\zeta }_r=8\pi \eta a^3$. Thermal rotational diffusion at absolute temperature $T$ then givesHere $k_B$ is the Boltzmann constant; the translational expression is the Stokes–Einstein relation. The ratio has units of inverse length squared, so comparing the numerical values of $D_r$ and $D_t$ without their units is meaningless. At fixed temperature and dynamic viscosity, decreasing radius tenfold increases $D_r$ a thousandfold but $D_t$ only tenfold. This strong size dependence explains why small particles can lose directional memory extremely quickly.

Choose $T=298\mathrm{K}$ and $\eta =0.001\mathrm{Pa}\mathrm{s}$ as illustrative model inputs. Using $k_B=1.380649\times 10^{-23}\mathrm{J}{\mathrm{K}}^{-1}$ gives $k_{B}T=4.11433402\times 10^{-21}\mathrm{J}$ and the following predictions:These are bulk three-dimensional sphere calculations, explicitly not a fit to the confined ellipsoid experiment below. Walls, nonspherical shape and departure from the no-slip boundary condition can change rotational resistance. A visible surface mark is still necessary to read the sphere's turning in an image, even when its rotational diffusion coefficient is well defined physically.

Take identical rigid spheres of radius $a$, particle density ${\rho }_p$ and displaced liquid density ${\rho }_f$, at uniform absolute temperature $T$. Let $z$ increase upwards from a reference layer and write $\Delta \rho ={\rho }_p-{\rho }_f>0$. The buoyancy correction gives the effective gravitational massThe actual grain mass is ${\rho }_{p}V$; the displaced liquid subtracts ${\rho }_{f}V$ from the effective gravitational mass. In a dilute, noninteracting particle suspension, the Boltzmann distribution therefore givesHere $n$ is number density and $k_B$ is the Boltzmann constant. Equivalently, the dilute osmotic pressure $\Pi =n{k}_{B}T$ supports the grains above it: $d\Pi /dz=-n{m}_{b}g$. Dividing by $n$ and integrating yields the same exponential function. No dynamic viscosity appears in this equilibrium argument. If ${\rho }_p<{\rho }_f$, the population instead increases upwards; at equal densities this gravitational model predicts uniform $n$.

A microscope counts a finite optical slab, not an infinitesimal plane. For equal viewed area $A$ and slab thickness $\delta$, an ideal uniform detection window centred at $z$ givesThus equal optical volumes preserve the exponential function count ratio exactly, even for finite $\delta$, provided the windows lie wholly inside the accessible chamber and have the same detection response. Unequal depth response or overlapping a boundary can corrupt this inference. Repeated observations estimate the mean count; closely spaced visits need not be independent because grains remain in view.

For an original chosen model, set $a=0.250\mu \mathrm{m}$, ${\rho }_p=1200\mathrm{kg}{\mathrm{m}}^{-3}$, ${\rho }_f=1000\mathrm{kg}{\mathrm{m}}^{-3}$, $T=298\mathrm{K}$ and $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$. These are illustrative inputs, not Jean Baptiste Perrin's measurements. Using $k_B=1.380649\times 10^{-23}\mathrm{J}{\mathrm{K}}^{-1}$ givesIf equal observation effort gives expected count $400$ in the reference layer, the expected counts $20$, $40$ and $60\mu \mathrm{m}$ higher are approximately $214.3$, $114.8$ and $61.49$. Fractions describe expected values over observations, not fractional grains. Using the actual grain mass instead of its buoyancy-corrected value would make the inferred Avogadro constant six times too small for these chosen densities.

For a steady incompressible flow of a Newtonian fluid with constant density $\rho$ and constant dynamic viscosity $\mu >0$, neglecting inertia in the Navier-Stokes equation givesThe first equation is local force balance; the second is mass conservation. The term $\mu {\nabla }^2\mathbf{u}$ is the divergence of the viscous part of the Newtonian fluid stress tensor. The body force $\mathbf{f}$ has units of force per unit volume, $\mathrm{N}{\mathrm{m}}^{-3}$, and equals $\rho \mathbf{g}$ for uniform gravity. Appropriate boundary data for velocity or traction, together with suitable distant conditions for an exterior flow, complete the problem; a rigid surface normally imposes the no-slip boundary condition.

With uniform gravity, write $p=p_h+p^{\prime }$ where the hydrostatic pressure satisfies $\nabla p_h=\rho \mathbf{g}$. Then the disturbance equation is $-\nabla p^{\prime }+\mu {\nabla }^2\mathbf{u}=0$. This removes gravity from the fluid disturbance equation, while the hydrostatic pressure still supplies buoyancy on an immersed body. Omitting it from the total force balance would give the wrong Stokes settling velocity.

Stokes flow is the limit of fluid mechanics in which inertia is negligible: pressure gradients, viscosity and any body forces balance at each instant. A small Reynolds number can arise from a small object, a slow velocity, or a large dynamic viscosity. A microscopic sphere moving through water is a useful example: the surrounding fluid continually dissipates the supplied power through viscous dissipation.

The picture is a meridional slice of the three-dimensional axisymmetric flow around a sphere, with the sphere held fixed and the distant fluid moving rightward at speed $U$. The no-slip boundary condition makes both components of velocity vanish on the surface. The streamlines have fore-aft reflection symmetry; reversing the incident velocity reverses their arrows. The underlying field is the translating sphere in Stokes flow solution viewed in the stationary-sphere frame.

For an incompressible flow of a Newtonian fluid with constant dynamic viscosity $\mu$, the steady Stokes equation isHere $p$ is pressure, $\mathbf{u}$ is velocity, and $\mathbf{f}$ is body force per unit volume. The pressure gradient enforces incompressible flow; it need not be small just because inertia is small. At speed $U$ and length $L$, the motion-induced pressure gradient has typical scale $\mu U/L^2$, corresponding to a disturbance pressure difference of order $\mu U/L$. See Shreyas Mandre's fluid dynamics notes, sections 9.1.1–9.1.2.

The controlling Reynolds number iswhere $\rho$ is density and $\nu$ is kinematic viscosity. The condition ${\mathrm{Re}}_L\ll 1$ suppresses advection of momentum on the length $L$. Rapid changes require an additional check of the viscous diffusion time: small Reynolds number alone does not eliminate the time derivative in the Unsteady Stokes equation.

Two rewards of this limit are the exact Stokes's law for the drag on a sphere and kinematic reversibility of Stokes flow. The first predicts settling velocities and connects to the Stokes–Einstein relation; the second explains why a microswimmer needs a nonreciprocal stroke, as expressed by the scallop theorem. Despite viscous dissipation, reversing the driving can retrace the motion: an especially striking consequence of Linearity of Stokes flow.

George Gabriel Stokes obtained the sphere drag and its terminal velocity in his 1851 investigation of how viscosity changes physical pendulum motion. He also used the slow fall of small water droplets to discuss clouds. The paper was read on 9 December 1850; 1851 is its publication date. See Stokes, On the effect of the internal friction of fluids on the motion of pendulums, original 1851 printing, especially articles 42–43; the formulas are also accessible in the collected-paper reprint, articles 42–43.

Take downward as positive and let gravity have magnitude $g$. For a rigid sphere of density ${\rho }_p$ in otherwise stationary fluid of density ${\rho }_f$, its weight minus buoyancy is $\frac{4}{3}\pi a^3({\rho }_p-{\rho }_f)g$. 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 $a^3$ while linear drag at fixed speed grows as $a$; their balance gives the settling-velocity radius scaling $U_s\propto a^2$. 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 ${\rho }_f=1000\mathrm{kg}{\mathrm{m}}^{-3}$, $\mu =1.00\times 10^{-3}\mathrm{Pa}\mathrm{s}$, ${\rho }_p=2500\mathrm{kg}{\mathrm{m}}^{-3}$, $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$, and $a=10.0\mu \mathrm{m}$. These are specified rounded inputs rather than a claim of exact water properties at every temperature. They giveThe diameter is $20.0\mu \mathrm{m}$, and both the net downward weight and the upward drag force are $6.16\times 10^{-11}\mathrm{N}$. At this terminal velocity, falling $1.00\mathrm{mm}$ takes $3.06\mathrm{s}$. At $a=20.0\mu \mathrm{m}$ the prediction is $1.308\mathrm{mm}{\mathrm{s}}^{-1}$ with ${\mathrm{Re}}_d=0.05232$; at $a=50.0\mu \mathrm{m}$ it is $8.175\mathrm{mm}{\mathrm{s}}^{-1}$ with ${\mathrm{Re}}_d=0.8175$, 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].

For an isolated rigid sphere of radius $a$ moving at velocity $\mathbf{U}$ relative to a distant uniform Newtonian fluid, the drag force exerted by the fluid isThe assumptions are steady incompressible flow, constant dynamic viscosity $\mu$, an unbounded surrounding fluid, the no-slip boundary condition, and a small diameter-based particle Reynolds number, ${\mathrm{Re}}_d=2a{\rho }_{f}U/\mu \ll 1$. Integrating the surface traction in translating-sphere Stokes flow gives the factor $6\pi$: the disturbance pressure contributes $2\pi \mu aU$ to the drag magnitude and the viscous shear stress contributes $4\pi \mu aU$. See Richard Fitzpatrick, Axisymmetric Stokes Flow Around a Solid Sphere, equations 10.120–10.129.

The result is linear drag: doubling $U$ doubles the drag force, and doubling $a$ at fixed $U$ also doubles it. The required power is $6\pi \mu a{U}^2$, 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 $4\pi \mu a$ 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.