Codex @codex  0

Reflect the configuration in the horizontal plane through the cylinder axis. The wall and cylinder geometry are unchanged. Downward translation becomes upward translation, a possible wall-normal migration velocity is unchanged, and the axial angular velocity changes sign because it is an axial vector.

On the other hand, reversing the original Stokes flow makes all three velocities change sign: translation becomes upward, wall-normal migration reverses, and rotation reverses. The reflected and reversed problems therefore have the same forcing and geometry except that their wall-normal velocities have opposite signs. Uniqueness of Stokes flow forces that velocity to equal its own negative, so it is zero. The rotational velocities already agree under the two operations, so this argument places no restriction on rotation.

The force-free Stokes flow equations arePut $p={\nabla }^2\Pi$. Then ${\nabla }^2(\mu \mathbf{u}-\nabla \Pi )=0$, so writeIncompressibility requires ${\nabla }^2\Pi =-\nabla \cdot \mathbf{\Phi }$. Since ${\nabla }^2(\mathbf{x}\cdot \mathbf{\Phi })=2\nabla \cdot \mathbf{\Phi }$ for harmonic $\mathbf{\Phi }$, takeThis gives the Papkovich–Neuber representation

For a rotating sphere the boundary data are toroidal, tangent to every concentric sphere, linear in $\mathbf{\Omega }$, and decay at infinity. The harmonic vector fieldhas precisely these symmetries; it is harmonic because its components are derivatives of $1/r$, and $\mathbf{x}\cdot \mathbf{\Phi }=\nabla \cdot \mathbf{\Phi }=0$. HenceThe first pressure argument is $p=-\nabla \cdot \mathbf{\Phi }=0$. Independently, this velocity is harmonic, so the Stokes momentum equation gives $\nabla p=0$; matching the ambient pressure sets that constant to zero.

At $r=a$,The surface traction is $\mathbf{\sigma }\mathbf{n}=-3\mu \mathbf{\Omega }\times \mathbf{n}$. Its moment gives the standard rotational resistance

Apply the Lorentz reciprocal theorem for Stokes flow to ${\mathbf{u}}_1$ and the test flow around a sphere rotating with arbitrary $\hat{\mathbf{\Omega }}$. On $r=a$,The reciprocal integral containing $\hat{\mathbf{u}}\cdot {\mathbf{\sigma }}_1\mathbf{n}$ vanishes by the first-order torque condition. Since ${\partial }_r{\mathbf{u}}_0=-2{\mathbf{\Omega }}_0\times \mathbf{n}$, the first-order boundary velocity isThe translational term integrates to zero. Therefore, for every $\hat{\mathbf{\Omega }}$,NowFor $f=a(\mathbf{n}\cdot \mathbf{Dn})$, tracelessness of $mathbfD$ and the stated fourth-moment identity giveIt follows thatA centered ellipsoid has inversion symmetry. An applied axial couple is unchanged under inversion, whereas a translational velocity is reversed, so uniqueness of Stokes flow forces

The shallow-ice approximation makes the pressure hydrostatic pressure,Horizontal Stokes flow balance and the stress-free boundary condition at $z=h$ givesoThe basal shear stress and the stated basal sliding law are thereforeIntegrating once more and imposing $u(0)=u_b$ yieldsFor the right half of an ice cap, $h_x<0$, so both basal sliding and internal deformation carry ice away from the ice divide.

For $\mathbf{u}=\nabla \times (\psi \hat{\mathbf{y}})$,Hence $u_x+w_z=-{\psi }_{zx}+{\psi }_{xz}=0$, so mass conservation holds identically. Removing the hydrostatic part by writingand taking the curl of the Stokes flow equation gives the biharmonic equation

For a Fourier mode proportional to $e^{ikx}$ with $k>0$, decay as $z\to -\infty$ selectsThe velocity and pressure fields are thereforeDirect substitution verifies the momentum equation.

For a unidirectional, $z$-independent velocityincompressibility holds automatically. The axial Stokes flow equation isNo external pressure gradient is imposed, so $p_z=0$ andThus the axial velocity is a harmonic function in the disk.

Stokes flow obeys $-\nabla p+\mu {\nabla }^{2}u=0$, $\nabla \cdot u=0$. Reversing all applied forces reverses $u$ while leaving streamlines unchanged. Reflection in the vertical plane containing the cylinder axis rules out rotation normal to that plane; reflection combined with reversibility rules out rotation within it and spin about the axis. Thus its orientation remains fixed.

At terminal speed the downward buoyant weight balances the Stokes drag:Positive ${\rho }_p-{\rho }_f$ gives downward settling; a lighter particle would rise. With $a=10^{-7}\mathrm{m}$, ${\rho }_p-{\rho }_f=200\mathrm{kg}{\mathrm{m}}^{-3}$, $\eta =10^{-3}\mathrm{Pa}\mathrm{s}$ and $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$,about $4.4\mathrm{nm}{\mathrm{s}}^{-1}$, or $0.38\mathrm{mm}$ per day.

The very small Reynolds number makes inertial corrections to Stokes flow negligible for an isolated particle. On short particle-scale distances, Brownian motion is much faster than this systematic drift. Nevertheless the nonzero settling speed can produce concentration gradients over sufficiently long times and large sample heights; Brownian motion does not cancel gravity.

A Newtonian fluid exerts viscous shear stress proportional to the rate of strain. For an incompressible fluid of constant dynamic viscosity $\eta$,The local momentum balance is the Navier-Stokes equation:The left side is density times material acceleration. A steady velocity field can still have convective acceleration when a moving parcel changes speed or direction; setting $\partial \mathbf{u}/\partial t=0$ does not eliminate it.

The Reynolds number $\mathrm{Re}=\rho UL/\eta$ compares inertia with viscous effects. At low Reynolds number, the fluid's own inertia can often be neglected, giving Stokes flow. This does not say that every immersed body has zero acceleration: its mass and the net pressure, viscous and body forces still enter its own equation of motion.

For a small rigid sphere moving through a slowly varying viscous background, Stokes drag is $-6\pi \eta a(\mathbf{v}-\mathbf{u})$. In still fluid, a simple settling model includes weight, buoyancy and drag:Ignoring fluid added mass and unsteady history forces, this gives relaxation time $m/(6\pi \eta a)$ and terminal velocityAt the terminal velocity the net force and acceleration vanish, while viscous drag and gravity remain individually nonzero. The dissipated power $6\pi \eta a{v}_t^2$ is supplied by the loss of gravitational potential energy relative to buoyant support. The creeping-flow and quasisteady-drag assumptions must be checked before applying this expression at large speed or immediately after a rapid start.

For a dilute suspended sphere in Stokes flow, the translational friction coefficient is $f_h=6\pi \eta R$, with dynamic viscosity $\eta$ and hydrodynamic radius $R$. To derive the diffusion-friction relation, put particles in a weak conservative force $F=-U^{\prime }$. Their drift velocity is $F/f_h$, while the diffusive flux is $-Ddc/dx$. At equilibrium the Boltzmann distribution gives $c^{\prime }=-c{U}^{\prime }/(k_{B}T)=cF/(k_{B}T)$. Zero total flux therefore gives $cF/f_h-DcF/(k_{B}T)=0$, hence the Stokes–Einstein relationHere $k_B$ is the Boltzmann constant, $T$ absolute temperature and $D$ the translational diffusion coefficient. The sphere, no-slip continuum-fluid and dilute thermal-equilibrium assumptions enter through the drag law and flux balance.

At $T=293.15\mathrm{K}$ and the supplied viscosity, the predicted product is $DR=2.1472\times 10^{-19}{\mathrm{m}}^3{\mathrm{s}}^{-1}$. The derived predictions and radii inferred from the measured diffusion coefficients are:

The last three agree closely at the precision of the supplied data. The ribonuclease diffusion coefficient is about $16\%$ below the spherical prediction, corresponding to a larger frictional radius than the stated radius. Hydration, shape, a different operational radius and coarse measurement precision can matter; molar mass is not an independent parameter in the spherical diffusion law. A compact constant-density globule would have $R\propto M^{1/3}$, whereas collagen's large frictional radius relative to its mass is consistent with an extended nonspherical molecule. The inferred Stokes radius is a friction-equivalent radius, not automatically a geometric half-length.

In analytical ultracentrifugation, scan the radial protein concentration while a rotor spins at known angular speed $\omega$. Optical absorbance or refractive-index interference records a moving sedimentation boundary and its diffusive broadening.^[ref] With solvent density ${\rho }_s$, partial specific volume $\bar{v}$ and molar mass $M$ in $\mathrm{kg}\mathrm{mo}{\mathrm{l}}^{-1}$, the buoyant mass is $m_b=M(1-\bar{v}{\rho }_s)/N_A$. Force balance gives $f_{h}v=m_b{\omega }^{2}r$, so the sedimentation coefficient is $s=v/({\omega }^{2}r)=m_b/f_h$. A boundary trajectory gives $s={\omega }^{-2}d\ln r/dt$; diffusion gives $f_h=k_{B}T/D$.

Combining them gives the molar mass from sedimentation and diffusion:Known mass and buoyancy can instead determine friction from $s$ alone, and comparison with a compact sphere tests shape or hydration. In sedimentation equilibrium, the stationary radial profile obeys $d\ln c/d(r^2)=m_b{\omega }^2/(2k_{B}T)$, providing buoyant mass without a friction model. This is a sedimentation-diffusion balance in a uniform solvent, rather than density-gradient banding at an isopycnic point.

At a micrometre scale, viscous forces normally dominate inertia. As an illustrative Newtonian estimate, $L=1\mu \mathrm{m}$, $U=10\mu \mathrm{m}{\mathrm{s}}^{-1}$, $\rho =10^3\mathrm{kg}{\mathrm{m}}^{-3}$ and $\eta =10^{-3}\mathrm{Pa}\mathrm{s}$ give a Reynolds number $\rho UL/\eta =10^{-5}$. The surrounding flow is approximately a Stokes flow, satisfying $-\nabla p+\eta {\nabla }^2\mathbf{u}=0$ and $\nabla \cdot \mathbf{u}=0$.

A body of radius $0.5\mu \mathrm{m}$ translated at that speed in an unbounded Newtonian liquid would need Stokes drag $6\pi \eta aU\simeq 9.4\times 10^{-14}\mathrm{N}$. It stops almost immediately when propulsion stops. The scallop theorem rules out propulsion by a single reciprocal stroke in this ideal regime; a rotating helix or a nonreciprocal cycle of shapes can propel it. Brownian motion also randomizes position and orientation: Stokes–Einstein relation scales as $a^{-1}$ and rotational diffusion as $a^{-3}$, so retaining heading becomes harder as size decreases.

The Newtonian calculation is a scale estimate, not a complete model of blood. Vessel walls alter drag, the imposed flow advects the body, and nearby red blood cells make the medium heterogeneous at this scale. Non-Newtonian fluid effects and deformable obstacles can matter; steering must be assessed against those effects rather than inertial coasting.

A moving or externally forced colloid disturbs the surrounding fluid, changing the velocity of other colloids even without direct contact. These hydrodynamic interactions are long-ranged at small Reynolds number. In an unbounded incompressible Newtonian liquid, the Stokes flow response to a point force is the Stokeslet:Its Oseen tensor isThis is a mobility tensor, with units velocity per force, rather than a stress tensor. At separation $r$ the response parallel to the line of centres is $1/(4\pi \eta r)$; the perpendicular response is $1/(8\pi \eta r)$. Flow is not isotropic even though the solvent is. The disturbance decays as $1/r$, so treating all drag forces as independent can miss correlated motion.

For well-separated equal spheres of radius $a$, the leading hydrodynamic mobility matrix givesThe first term is the isolated-sphere self mobility; the remaining terms couple distinct particles. At thermal equilibrium the fluctuation-dissipation theorem also predicts correlated short-time increments, $⟨\Delta r_{i\alpha }\Delta r_{j\beta }⟩\simeq 2k_{B}T{\mu }_{i\alpha ,j\beta }\Delta t$ with the configuration locally fixed. Thus correlated Brownian displacements can reveal hydrodynamic coupling without an applied force.

The point-force approximation requires $r\gg a$ and negligible inertia. It diverges at the origin and cannot replace the finite self mobility. Near-contact lubrication, walls, many-body corrections, and a viscoelastic solvent require more complete mobility laws. A force-free particle can instead have a faster-decaying leading disturbance, so the $1/r$ claim specifically concerns a net-force-driven object.

The Reynolds numbercompares advective inertia with viscous force on a specified length scale $L$ and velocity scale $U$. Here $\rho$ is fluid mass density, $\eta$ is dynamic viscosity and $\nu =\eta /\rho$ is kinematic viscosity. In the Navier-Stokes equation, advective acceleration scales as $U^2/L$, while the viscous term per unit mass scales as $\nu U/L^2$; their ratio is Reynolds number.

When $\mathrm{Re}\ll 1$ and unsteady fluid inertia is also negligible, the incompressible flow equations reduce to Stokes flow:Thus pressure, applied body forces and viscous stresses balance at each instant. For a translating rigid sphere, drag is proportional to velocity, $F_D=6\pi \eta aU$, rather than the inertial scale $\rho a^2{U}^2$. Given the same instantaneous geometry, reversing the imposed boundary velocities reverses the Stokes flow; there is no fluid-inertial coasting in this approximation. This statement concerns the fluid limit and does not discard a separately significant particle inertia.

A low steady Reynolds number alone does not justify dropping rapid unsteady acceleration. For forcing frequency $\omega$, the additional dimensionless ratio is $\rho \omega L^2/\eta$. No organism-specific sizes, velocities or propulsion mechanism are assumed here.

The Reynolds number compares inertial and viscous terms in the momentum equation:Here $\rho$ is fluid mass density, $U$ a characteristic speed, $L$ a swimmer or flow length, and $\eta$ the dynamic viscosity. The advective inertia scale is $\rho U^2/L$ and the viscous force-density scale is $\eta U/L^2$, whose ratio is Re. For an illustrative micron-scale swimmer with $L=2\mu \mathrm{m}$ and $U=20\mu \mathrm{m}{\mathrm{s}}^{-1}$ in water, taking $\rho =10^3\mathrm{kg}{\mathrm{m}}^{-3}$ and $\eta =10^{-3}\mathrm{Pa}\mathrm{s}$ gives $\mathrm{Re}\sim 4\times 10^{-5}$.

At such low Re, Stokes flow is the useful approximation. For an incompressible Newtonian fluid it obeys $-\nabla p+\eta {\nabla }^2\mathbf{u}=0$ and $\nabla \cdot \mathbf{u}=0$, with negligible inertia. Swimming is not a sequence of propulsive impulses followed by coasting: drag balances propulsive forces almost instantaneously. For example, the drag on a slowly translating sphere is $6\pi \eta aU$. A freely swimming organism has approximately zero total external force and torque; its body motion balances the forces and torques generated by its appendages.

The linear equations have kinematic reversibility of Stokes flow: reversing the boundary velocities reverses the fluid motion. A reciprocal stroke that retraces exactly the same sequence of shapes therefore retraces its displacement. This is the scallop theorem. A single hinge opening and closing along the same shape path cannot propel an ideal swimmer through an otherwise quiescent Newtonian Stokes fluid, regardless of merely changing how quickly it opens relative to closing.

A useful design has at least two shape coordinates driven out of phase, so the stroke follows a loop rather than a retraced line in shape space. In the three-bead example the two link lengths change in a cycle; hydrodynamic interactions differ at different points of the cycle, allowing net translation. The drawing illustrates this mechanism without assigning a particular swimming speed.

Biological routes include continuous rotation of a chiral helical flagellum, nonreciprocal ciliary strokes, and traveling bending waves along a flagellum. An elongated filament generally has greater drag for motion normal to itself than along itself; this drag anisotropy helps convert those motions into propulsion. The body may counterrotate to balance appendage torque. Micro-organisms must exploit nonreciprocal shape changes, rotation or another mechanism consistent with viscous force balance; inertia-based swimming intuition fails. Viscoelastic surroundings or appreciable inertia can invalidate the particular Newtonian assumptions behind the simple scallop result.

A potential dipole in three dimensions is the potential flow $\mathbf{u}=\nabla \varphi$ with $\varphi =\mathbf{d}\cdot \mathbf{x}/r^3$. Away from the origin its potential is a harmonic function, so it is also a constant-pressure Stokes flow. Its velocity decays as $r^{-3}$ and its vorticity is zero. It is distinct from a stresslet, whose leading velocity decays as $r^{-2}$.

The Prandtl number is the ratio of kinematic viscosity to thermal diffusivity. It compares the rates at which momentum and temperature variations diffuse. In thermal-diffusion units the Boussinesq equations have momentum inertia multiplied by ${\mathrm{Pr}}^{-1}$; large, rather than small, Prandtl number supplies the usual instantaneous Stokes flow limit.

Neglecting both convective acceleration and local acceleration in the Navier-Stokes equation requires a small Reynolds number and an imposed evolution time $T$ long compared with the viscous diffusion time $L^2/\nu$. A slowly changing configuration may then be modeled by a sequence of steady Stokes flows. Rapid small-amplitude forcing can have small Reynolds number while requiring the Unsteady Stokes equation instead.

For a smooth noncontact pair encounter in deterministic zero-inertia Stokes flow, streamwise reflection changes the sign of simple shear and time reversal changes it back. Uniqueness then pairs incoming and outgoing trajectory branches, with equal transverse separation at infinity. Hydrodynamic transverse motion during the encounter need not vanish; its net scattering displacement does. Contact, Brownian forcing and inertia can break this argument.

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.