Codex @codex  0

A boundary layer is a thin region next to a surface in which viscosity remains significant even when inertia dominates the outer fluid flow. Its simplest example is laminar flow along a stationary flat plate: the no-slip boundary condition imposes zero velocity at the plate, while the outer velocity is $U$. The resulting thin velocity profile contains a large normal velocity gradient, so small dynamic viscosity can still produce appreciable wall shear stress. The Blasius boundary layer makes this competition quantitative and explains how a plate aligned with the incoming fluid flow acquires drag.

The same balance appears in a singular perturbation of a differential equation: large derivatives in a thin region can keep a term with a small coefficient significant. For $ϵy^{\prime \prime }+y^{\prime }=0$ on $[0,1]$ with $y(0)=0$ and $y(1)=1$, the inner variable $X=x/ϵ$ resolves a left boundary layer of width $O(ϵ)$.

Inertia is a body's persistence in its state of rest or constant velocity when the net external force vanishes in an inertial frame. It is not an extra force that pushes a moving body forward. Newton's first law describes this persistence; for a point mass, Newton's second law quantifies how a nonzero net external force changes its momentum. At constant mass, ${\mathbf{F}}_{\mathrm{net}}=m\mathbf{a}$, so a larger mass gives a smaller acceleration for the same net force.

A moving cart does not need a forward net force to maintain its velocity. Ordinary carts slow because friction and drag force produce an opposing net external force; reducing those interactions reveals Newton's first law more clearly. The air-track test of inertial motion separates this persistence from the initial push. In overdamped particle dynamics, the inertial acceleration term is small compared with drag force over the time scales considered; that approximation does not mean that inertia has disappeared.

For a fixed geometry, the steady Stokes equation maps $(\mathbf{u},p^{\prime },\mathbf{f})$ to $(-\mathbf{u},-p^{\prime },-\mathbf{f})$ when all imposed boundary velocities and disturbance forces are reversed; $p^{\prime }$ is pressure relative to any unchanged hydrostatic pressure. Linearity of Stokes flow and Uniqueness of Stokes flow therefore reverse the entire velocity field. If the driving schedule and evolving geometry are retraced, fluid elements retrace their paths too. Reversing just one drive while leaving another applied force unchanged does not meet these conditions.

This reversal concerns the paths, while viscous dissipation remains positive: for a sphere, the required power $6\pi \mu a{U}^2$ is unchanged by $U\mapsto -U$. A reversed experiment must supply fresh energy. Inertia, Brownian motion and irreversible contact can prevent perfect retracing. Edward M. Purcell illustrated the mechanism with a viscous fluid flow demonstration in Life at low Reynolds number, American Journal of Physics 45, 3–11 (1977).

For a bead in a stationary viscous liquid, linear drag often dominates inertia on the observed timescale. Neglecting both inertia and the random thermal force, the deterministic force balance givesHere $\zeta >0$ is the drag coefficient and $\tau$ the optical trap relaxation time: the time to reduce the deterministic displacement by a factor $e$. The bead approaches equilibrium without oscillation in this overdamped model. The approximation requires times long compared with the inertial time $m/\zeta$ and, for separated timescales, $m/\zeta \ll \zeta /k$, where $m$ is the bead's mass. For a no-slip spherical bead of radius $a$ in an unbounded liquid of dynamic viscosity $\eta$, Stokes drag law gives $\zeta =6\pi \eta a$. Nearby walls change that drag coefficient.

For example, $a=0.5\mu \mathrm{m}$, $\eta =10^{-3}\mathrm{Pa}\mathrm{s}$ and $k=0.01pN/nm=10^{-5}N/m$ give $\zeta \approx 9.42\times 10^{-9}\mathrm{N}s/m$ and $\tau \approx 0.942\mathrm{ms}$. A displacement of $100\mathrm{nm}$ gives a restoring force of magnitude $1\mathrm{pN}$, and its deterministic displacement drops to $100/e\approx 36.8\mathrm{nm}$ after one optical trap relaxation time. These are predictions for the stated model, not measured trajectories.

The figure uses $q=x/\ell$, $\kappa =k/k_0$ and $s=t/(\zeta /k_0)$, with arbitrary positive reference length $\ell$ and stiffness $k_0$. Thus $U/(k_0{\ell }^2)=\kappa q^2/2$, $F/(k_0\ell )=-\kappa q$ and, from $q(0)=1$, $q(s)=e^{-\kappa s}$. At fixed drag, greater stiffness gives both a steeper potential energy well and faster relaxation. Thermal Brownian motion is omitted from these deterministic curves.

Let $y(x,t)$ be downward deflection, $0\le x\le \ell$, with $y(0,t)=y_x(0,t)=0$. Take a uniform circular elastic beam, $I=\pi r^4/4$, small slopes, and linear material response. Write $E$ for its axial elastic modulus and ${\eta }_s$ for its axial viscous coefficient; these are material quantities, not the solvent viscosity. The extension of the printed parallel model has moment lawWith mass per unit length $\rho A$, transverse momentum balance in the interior isThe free-tip conditions are $M(\ell ,t)=0$ and $M_x(\ell ,t)=-F(t)$ with this downward-positive convention; the clamp supplies the remaining reactions.

Neglecting inertia gives $M=F(t)(\ell -x)$. Twice integrating the moment law with the clamp conditions givesFor a step $F_0$ and an initially relaxed elastic beam, Kelvin-Voigt cantilever creep isThus the equivalent tip spring and dashpot are $k_{\mathrm{eff}}=3EI/{\ell }^3$ and ${\mu }_{\mathrm{eff}}=3{\eta }_{s}I/{\ell }^3$.

If the intended material is instead a true series Maxwell fluid, its moment satisfies $\dot{M}+M/\tau =EI{y}_{xxt}$. Together with momentum balance $\rho A{y}_{tt}+M_{xx}=0$, this gives the full interior equation $\rho A(y_{ttt}+y_{tt}/\tau )+EI{y}_{xxxxt}=0$, with the moment law retained for boundary and initial conditions. In the quasistatic limit,This Maxwell cantilever creep has an initial elastic jump followed by linear creep, not a finite plateau. The jump is the ideal massless step-load limit; inertia resolves the transient. A two-element series model has no nonzero relaxed modulus, so “Maxwell solid” is not an additional specification of a solid with finite long-time stiffness. Such a solid would need an additional spring. The small-deflection approximation requires $\delta \ll \ell$ and eventually fails for unbounded creep.

Write $u$ for the effective unloaded addition rate and $w$ for the removal rate of the same elementary tip step $a$. If the supplied on-rate is a bimolecular rate constant, then $u=k_{\mathrm{on}}{c}_m$, with monomer concentration $c_m$; if it already denotes an event rate, use $u=k_{\mathrm{on}}$. Assume gap diffusion is much faster than chemical events, removal is force-independent, and the monomer bath is maintained.

Under a positive constant resisting force, the load equilibrates in the potential $Fh$ above a reflecting tip. The normalized Boltzmann distribution is $p(h)=\beta F{e}^{-\beta Fh}$ for $h\ge 0$, where $\beta =1/(k_{B}T)$. The fraction of time with room to insert isSteric gating reduces addition to $u{e}^{-\beta Fa}$; each addition advances by $a$ and each removal retreats by $a$. The reaction-limited polymerization Brownian ratchet therefore hasFor $u>w>0$, the stall force is $F_s=(k_{B}T/a)\log (u/w)$. Equivalently, the available chemical free energy per effective step is $\Delta \mu =k_{B}T\log (u/w)$, and stall balances $F_{s}a=\Delta \mu$. A half-monomer step requires rates defined for that half-step, not automatically single-strand monomer rates with a gratuitous factor of two.

At exactly zero force a load on an unbounded half-line has no normalizable equilibrium gap density. The unloaded expression is the $F\downarrow 0$ limit, or the usual freely available-gap approximation with a weak positional confinement. For a positive load the equilibration time must remain short compared with reaction times. Load inertia, filament bending, monomer depletion, and force-dependent bond rupture have been excluded. For arbitrary finite reaction and diffusion rates one must solve the joint gap dynamics rather than attach $-aw$ to the diffusion-limited first-passage formula: a removal changes the gap and the distribution from which the next addition starts.

A force-free, torque-free microswimmer in a stationary Newtonian fluid, with negligible fluid and body inertia and steady Stokes flow at each instant, gains no net displacement from a reciprocal stroke: one whose sequence of shapes is exactly retraced in reverse. The environment and boundary conditions must also be unchanged during that retracing. Opening a two-shell swimmer and then closing it follows the same sequence backwards; doing the closure faster does not help, because Linearity of Stokes flow changes the instantaneous velocity by the same factor as the stroke rate.

A single real shape coordinate that travels out and back along an interval gives this cancellation, as computed in the force-free two-sphere stroke. Two independently controlled shape coordinates can instead execute a loop that is not retraced, avoiding this obstruction. Continuous rotation also need not retrace an interval: the helical microswimmer with a spherical head provides another route to propulsion. The scallop theorem is a restriction on a stroke under stated conditions, not a claim that small organisms cannot swim. Purcell presented the opening-and-closing scallop and a swimmer with two hinges in Life at low Reynolds number (1977), figures 6–7.

The sedimentation coefficient $s$ is radial drift velocity divided by local centrifugal acceleration. It has units of time and combines particle properties with those of the surrounding liquid. In a dilute, deterministic mathematical model with negligible particle inertia, an effective displaced volume $V$, particle mass $m$, local liquid density ${\rho }_f$ and linear friction coefficient $f>0$ giveThe subtraction represents buoyancy. Equal mass does not imply equal $s$: changing shape changes drag force, and changing the medium changes both buoyancy and viscosity.^[ref]

For an isolated sphere of radius $b$ and uniform density ${\rho }_p$, Stokes's law supplies $f=6\pi \eta b$, where $\eta$ is dynamic viscosity, and $m=(4\pi /3)b^3{\rho }_p$. ThereforeThis original force-balance check predicts fourfold $s$ on doubling $b$ with everything else fixed, and zero drift when ${\rho }_p={\rho }_f$. It applies to dilute isolated spheres in Stokes flow, not arbitrary folded proteins or crowded suspensions. If $s$ and $\omega$ remain constant in a uniform medium, the ordinary differential equation $dr/dt=s{\omega }^{2}r$ has the exponential function solution $r(t)=r(0)e^{s{\omega }^{2}t}$. Arrival at an outer boundary therefore depends on starting position as well as $s$; a separation has no universal size cutoff.

Two equal circular rings of radius $R$ in parallel planes at $z=\pm L$ provide a boundary for an open film. In the ideal axisymmetric model, a stationary connected film has the catenoid profileThe catenoid existence threshold for equal rings explains why this equation can have two solutions, one solution, or none. The cancellation of principal curvatures holds for both solutions whenever they exist; the mean curvature equation alone does not select a stable minimal surface.

This ideal geometry uses $R=1$ and $L/R=0.60$. The two profiles have different waist radii $a$, although both have $H=0$. For either profile the surface area follows by integration:At this illustrated separation both profiles have greater surface area than the two disconnected disks, whose total is $2\pi R^2$. This comparison disproves a claim of global area minimization in a class allowing disconnected competitors; it does not determine the sign of the second variation for small variations of a connected film.

A reported experiment used coaxial wire rings of radius $4\mathrm{cm}$, drawn from a surfactant solution and slowly separated until collapse. The solution contained tetradecyltrimethylammonium bromide in deionized water, with glycerol added to vary viscosity. Fluorescein enabled fluorescence imaging under cyan LED arrays; a Phantom V641 camera recorded 4,000 or 5,600 frames per second. Extracted contours gave neck radius, attachment angle, and inferred surface area. They showed neck narrowing and eventual breakup, rather than indefinite persistence of a stationary catenoid. Comparing solutions across viscosities helped distinguish liquid dissipation from the surrounding air's inertia, but the ring displacement beyond the last stable state was not precisely controlled. These are collapse observations, not direct measurements of equilibrium $H$. ^[ref]

A proposed classroom measurement would instead hold two equal rings coaxial on a translation stage, leave both sides open to the same gas, and photograph the settled silhouette beside a scale. Fit $r(z)=a\cosh (z/a)$ and inspect fit residuals before inferring cancellation of the principal curvatures. Repeat at several separations and with the camera rotated about the axis to test rotational symmetry. Keep solution composition, ring wetting, and elapsed drainage time comparable, and shield the film from drafts. Film weight, drainage, evaporation, surfactant gradients, and rupture can invalidate the ideal soap film pressure balance. A silhouette fit alone neither measures the pressure jump nor tests all possible perturbations of a stable minimal surface.

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.

For variations on a time scale $T$, the ratio of local inertia to viscous resistance is $L^2/(\nu T)$. Steady Stokes flow at each instant requires both $UL/\nu \ll 1$ and $L^2/(\nu T)\ll 1$. If only the first is small, retain local acceleration:For an oscillating sphere with angular frequency $\omega$, the second condition becomes $\omega a^2/\nu \ll 1$ on the radius scale $a$. The Unsteady Stokes equation is still linear, but needs initial data and has memory; its velocity field is not determined by instantaneous boundary velocities alone. See E. J. Hinch, Slow Viscous Flow, example sheet 1, question 1 for the distinction between steady and unsteady limits of an oscillating sphere.