Codex @codex  0

The drag coefficient measures a drag force against a specified scale of force. Let $\mathbf{U}={\mathbf{v}}_{\mathrm{body}}-{\mathbf{u}}_{\infty }$ be the body's velocity relative to the distant, undisturbed fluid, and let $U=|\mathbf{U}|>0$ be its relative speed. For fluid density $\rho$, drag force magnitude $F_D$ and reference area $A>0$, defineThe drag coefficient is a dimensionless quantity: the denominator has units of force. The reference area must be reported with it. For a sphere of radius $a$, the usual choice is the frontal area $A=\pi a^2$. Choosing its entire surface area $4\pi a^2$ instead divides the reported drag coefficient by four for the same physical drag force. Comparing coefficients without comparing their reference areas can therefore be misleading. ^[ref]

For a rigid sphere in steady flow through an unbounded Newtonian fluid with constant dynamic viscosity $\mu$, incompressible flow and the no-slip boundary condition, the creeping-flow limit requires the diameter-based particle Reynolds number ${\mathrm{Re}}_d=2a\rho U/\mu \ll 1$. The Stokes drag law gives $F_D=6\pi \mu aU$. Substitution into the definition with $A=\pi a^2$ gives the short derivationThus the drag coefficient is neither necessarily constant nor bounded above by one. The formula containing $U^2$ still describes linear drag here, because the drag coefficient varies as $1/U$. The relation $C_D=24/{\mathrm{Re}}_d$ is the leading small-Reynolds number approximation, not a measured curve valid at every Reynolds number. ^[ref]

For a worked example, stipulate dynamic viscosity $\mu =0.10\mathrm{Pa}\mathrm{s}$, density $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$, radius $a=10^{-4}\mathrm{m}$ and relative speed $U=10^{-3}\mathrm{m}{\mathrm{s}}^{-1}$. Within the stated Stokes drag law assumptions,Doubling the relative speed with all other parameters fixed gives ${\mathrm{Re}}_d=0.004$ and $C_D=6000$: the predicted drag force doubles to $3.76991118\times 10^{-7}\mathrm{N}$. Falsely holding the drag coefficient at $12000$ would instead quadruple the drag force. Both particle Reynolds numbers remain small, so the comparison is within the creeping-flow approximation.

The useful idea is to separate the force scale from the dependence on geometry and flow conditions. A drag coefficient can help compare experiments, but changing Reynolds number can change the coefficient itself; NASA Glenn explicitly stresses matching the relevant flow conditions. The elegant cancellation above also exposes a common mistake: a drag coefficient definition does not by itself impose a quadratic drag law.

Dynamic pressure is the bulk kinetic energy per unit volume associated with fluid motion at speed $U$ in a specified reference frame:Here $\rho =m/V$ is mass density. The word "pressure" identifies its units and its use as a pressure scale; $q$ is not an extra isotropic stress to add to the local static pressure in every surface traction. A stationary fluid can have large static pressure with $q=0$. Conversely, knowing $U$ and $\rho$ fixes $q$ without fixing static pressure. The definition remains useful in viscous fluid flow and compressible flow, even when the simplest Bernoulli equation is inapplicable. ^[ref]

For stipulated mass density ${\rho }_{\infty }=1.20\mathrm{kg}{\mathrm{m}}^{-3}$ and reference speed $U_{\infty }=20\mathrm{m}{\mathrm{s}}^{-1}$, $q_{\infty }=240\mathrm{Pa}$. At $40\mathrm{m}{\mathrm{s}}^{-1}$ with the same mass density, $q_{\infty }=960\mathrm{Pa}$. This quadratic scaling makes dynamic pressure a useful normalization for the lift coefficient, drag coefficient and pressure coefficient; it predicts four times a particular force only if its corresponding coefficient remains unchanged.

The connection with stagnation pressure requires hypotheses. For steady flow, constant mass density, inviscid flow, and gravity as the only body force, the Bernoulli equation gives, along one streamline,Here $p$ is static pressure, $g$ is gravitational acceleration, and $z$ is height. Bring that flow to rest without losses at the same height and the limiting static pressure is the stagnation pressure $p_0$:There must be no intervening pump supplying work or other energy input. Constancy along one streamline does not assert equal $B$ on all streamlines; irrotational flow in a connected region supplies that additional conclusion. Changing height requires retaining $\rho gz$. In compressible flow the stopping pressure difference generally differs from $\rho U^2/2$, and viscous dissipation can reduce recovered stagnation pressure. ^[ref]

A Pitot-static tube turns this relation into a measurement: its upstream-facing opening senses stagnation pressure, while separate side openings sample static pressure. Connect these to the positive and negative ports of a differential pressure transducer or manometer. Under the hypotheses above, the measured difference $\Delta p$ givesThus the stipulated $240\mathrm{Pa}$ and $1.20\mathrm{kg}{\mathrm{m}}^{-3}$ imply $20\mathrm{m}{\mathrm{s}}^{-1}$, not an independently measured speed. Align the probe with the incoming flow and keep its passages clear: blockage and pinched tubing invalidate the comparison. At very low speed, the difference can be smaller than instrument error; at high Mach number, compressibility and possibly a shock wave require a different reduction. ^[ref]

A useful control is differential-pressure calibration. For a digital manometer, zeroing establishes the reading at equal port pressures; it does not establish the gain away from zero. The manufacturer's Series 475 Handheld Digital Manometer: Specifications - Installation and Operating Instructions by Dwyer Instruments, Inc. directs the user to vent both ports to the atmosphere before zeroing. Its stated accuracy is a fraction of full scale, not of the reading, and depends on instrument temperature. Select a suitable range and record the instrument variant and units; the manual covers several ranges.^[ref]

An independent physical check can use a liquid-column manometer: apply a stable pressure difference to both instruments and compare the digital reading with the hydrostatic pressure inferred from the vertical level difference. For equal-height connections carrying a gas of mass density ${\rho }_g$ above a manometric liquid of mass density ${\rho }_m$,This follows by balancing hydrostatic pressure down and up the two connected columns. Use vertical height, even if a tube is inclined. Known mass densities, measured height and gravitational acceleration set the comparison scale. Check several differences across the working range and reverse the port connections to check the sign. This is a proposed independent check, not a claim that the cited teaching experiment performed a full instrument calibration. Recheck zero after reconnecting the sealed tubing; a zero adjustment alone cannot detect a blocked passage or an incorrect scale.

The dimensional linear friction coefficient relates drag force and velocity, $F_d=-\zeta v$, and has units of force divided by speed. It differs from the dimensionless aerodynamic drag coefficient. For an isolated sphere in Stokes flow, the Stokes drag law gives $\zeta =6\pi \eta a$.

Primary paper, in the English translation collected in Investigations on the Theory of the Brownian Movement: sites.pitt.edu/~jdnorton/lectures/Rotman_Summer_School_2013/Einstein_1905_docs/Einstein_Brownian_English.pdf

Sections 3–5 connect equilibrium drift and diffusion, the spherical-particle drag coefficient, and mean-squared displacement. The linked scan includes the 1926 translation reprinted in the 1956 Dover edition; the research paper itself is dated 1905.

A Reynolds number based on particle diameter and relative speed in the carrier fluid. Small values justify Stokes drag law; a high-Reynolds-number regime with nearly constant drag coefficient instead gives quadratic drag. The drag law must be evaluated in the appropriate regime for every species.

The filament bending modulus has units force times length squared. The perpendicular drag coefficient in this equation is drag per unit length divided by velocity, so its units are force times time per length squared. Hence $A/(\zeta \omega )$ has units length to the fourth power. Balancing oscillatory viscous drag $\zeta \omega h$ against elastic force density $Ah/{\ell }^4$ gives the elastohydrodynamic penetration lengthIncreasing frequency confines the response nearer the actuator. The chosen coefficient one in this characteristic length is a convention; the exact decay lengths acquire the root-dependent factors found below.

Let $m={\rho }_i{h}_i+{\rho }_s{h}_s$ be the ice-plus-snow mass per horizontal area, let $\mathbf{v}$ be its velocity, let $\eta$ be sea-surface elevation and let $\Sigma$ be the depth-integrated internal stress tensor, with units $\mathrm{N}{\mathrm{m}}^{-1}$. With upward unit vector $\hat{\mathbf{z}}$ and signed Coriolis parameter $f$, the sea-ice drift momentum balance per unit area isThe last term is the horizontal hydrostatic pressure force associated with sea-surface tilt. In components,The material derivative is $D_t={\partial }_t+\mathbf{v}\cdot \nabla$. Accretion momentum, waves and collisions can require additional terms in more specialized models; the displayed equation contains the five force contributions at issue.

Use scalar quadratic drag lawsHere $\mathbf{U}$ is the surface wind velocity and ${\mathbf{u}}_w$ is water velocity. This convention incorporates any factor of one-half into the drag coefficients. Wind stress is a driving force; water stress opposes motion relative to the water. The Coriolis force is perpendicular to ice velocity and does no work. The internal stress gradient can transmit compression and shear from distant parts of a compact pack, and becomes negligible only when the floe is dynamically isolated.

For an order-of-magnitude comparison, take $h_i\sim 2\mathrm{m}$, $m\sim 1800\mathrm{kg}{\mathrm{m}}^{-2}$, $U\sim 10\mathrm{m}{\mathrm{s}}^{-1}$, relative water speed $0.15$--$0.20\mathrm{m}{\mathrm{s}}^{-1}$, $|f|\sim 1.4\times 10^{-4}{\mathrm{s}}^{-1}$, ${\rho }_a\sim 1.3$, ${\rho }_w\sim 1025\mathrm{kg}{\mathrm{m}}^{-3}$, $C_a\sim 1.3\times 10^{-3}$ and $C_w\sim 5\times 10^{-3}$. Then wind stress is about $0.17\mathrm{N}{\mathrm{m}}^{-2}$, water stress about $0.12$--$0.21\mathrm{N}{\mathrm{m}}^{-2}$, and the Coriolis force about $0.04$--$0.05\mathrm{N}{\mathrm{m}}^{-2}$. Sea-surface slope $|\nabla \eta |\sim 10^{-6}$--$10^{-5}$ gives $mg|\nabla \eta |\sim 0.02$--$0.2\mathrm{N}{\mathrm{m}}^{-2}$. Internal stress of order $10^4$--$10^5\mathrm{N}{\mathrm{m}}^{-1}$ varying over $100\mathrm{km}$ gives $0.1$--$1\mathrm{N}{\mathrm{m}}^{-2}$. Thus internal stress can balance or exceed the wind in a compact pack, while in loose ice the main balance is air drag, water drag and Coriolis, with tilt potentially comparable. These are regime estimates, not a universal ranking.

An acceleration of $0.2\mathrm{m}{\mathrm{s}}^{-1}$ over a day gives only about $0.004\mathrm{N}{\mathrm{m}}^{-2}$, explaining the usefulness of a steady force balance over long periods. Over an hour the same change gives $0.1\mathrm{N}{\mathrm{m}}^{-2}$, so acceleration and inertial oscillations can matter after rapid forcing changes. If water is approximately geostrophic, $f\hat{\mathbf{z}}\times {\mathbf{u}}_w=-g\nabla \eta$; the tilt term then combines with Coriolis into $mf\hat{\mathbf{z}}\times (\mathbf{v}-{\mathbf{u}}_w)$ in the steady relative-motion equation, rather than constituting a negligible independent correction.

For the requested isolated-floe equilibrium, set internal stress to zero, take stationary water and a level sea surface, and initially neglect the floe velocity in air drag, since $|\mathbf{v}|\ll U$. Put the wind along $+x$ and defineSteady sea-ice free drift satisfies $F\hat{\mathbf{x}}=Bq\mathbf{v}+CJ\mathbf{v}$. Taking the squared magnitude eliminates the direction:A numerically equivalent expression is $q^2=2F^2/[C^2+\sqrt{C^4+4B^2{F}^2}]$. Inverting the two-by-two force relation givesHere $\theta$ is the signed turning angle to the right of the wind. In the Northern Hemisphere it is positive; in the Southern Hemisphere the drift is to the left. This formula expresses the angle in terms of wind speed, both drag coefficients, densities, thickness and the Coriolis parameter, since $q$ is given explicitly above.

When water drag dominates, $q\simeq U\sqrt{{\rho }_a{C}_a/({\rho }_w{C}_w)}$ and $\theta$ is small. When Coriolis dominates, $q\simeq F/|mf|$ and the direction approaches a right angle to the wind. Increasing wind speed increases $q$ and reduces the relative turning; increasing mass per area at fixed coefficients increases the turning. For the numerical example above, the formula gives $q\simeq 0.178\mathrm{m}{\mathrm{s}}^{-1}$ and $\theta \simeq 15.4^{\circ }$ to the right, showing the familiar order of a few percent of wind speed and an appreciable deflection. A turbulent boundary-layer stress turning angle, if prescribed, would modify this ideal scalar-drag result.

For completeness, retain the ice velocity in the air stress. Set $A={\rho }_a{C}_a$ and $w=|\mathbf{U}-\mathbf{v}|$. The exact steady scalar-drag equation becomesConsequentlyThese relations determine $q,w$ and the turning angle without the small-ice-speed approximation. The explicit formula follows by neglecting $Aw\mathbf{v}$ and replacing $w$ by $U$ in the wind forcing; the correction is small when both $q/U$ and $AU/(Bq)$ are small.

An iceberg has much larger mass per horizontal area than a thin ice floe. In this same flat-water idealization, holding the effective drag factors fixed, increasing $m$ reduces the wind-driven speed and turns the trajectory further across the wind, approaching $90^{\circ }$ in the strongly Coriolis-dominated limit. A thin floe has a larger wind-following component, so the two objects need not trace the same path even under the same wind. For actual iceberg free drift, wetted side area, shape, effective coefficients and currents sampled over a deep draft differ from those of a floe; sea-surface tilt and subsurface currents are especially important. Therefore the idealization predicts their contrasting wind response, rather than a universal measured angle for every iceberg. Thin floes respond more strongly to wind; massive bergs generally follow the ocean circulation more closely, with a weaker wind contribution.

To connect the prescribed interfacial drag coefficient to the energy supply, state the lumped interfacial power model for entrainment explicitly:It estimates the power transmitted by lid-driven flow to the mixing interface. Constant geometric and transfer factors are included in $c_D$ or $C_W$. It is the local stress-power interpretation required for the depth laws in this question, not an identity equating arbitrary lid shaft power to interfacial power. With the potential energy of two-layer entrainment from (a),In model S, $u_I=\Omega R$ and $B$ are constant, soThus the depth increases at a constant rate proportional to the inverse bulk Richardson number, with $h=h_0+2C_W{c}_D\Omega Rt/{\mathrm{Ri}}_B$, until the layer reaches $H$. Here the prescribed $u_I$ must be read as a characteristic speed. Literal pure solid-body rotation has zero radial velocity, so the PDF's parenthetical description of $u_I$ as radial cannot simultaneously be a strict component statement and the assumption $u_I=\Omega R$.

If the phrase about lid work is taken literally as shaft power $P_{\mathrm{lid}}=\Omega {\mathcal{T}}_{\mathrm{lid}}$, the interfacial stress alone does not specify the lid torque or its energy transfer to the interface. The preceding closure is therefore a visible model assumption, not a hidden consequence of solid-body rotation.

Take every finite positive drag coefficient $\gamma >0$ and defineThe coefficients of the characteristic polynomial are $c_3=2\gamma$, $c_2=K+{\gamma }^2$, $c_1=2L\gamma$ and $c_0=P{\gamma }^2$. The nontrivial Routh-Hurwitz criterion reduces exactly toThe printed strict chain $K>L>P>0$ makes every coefficient positive and both bracket terms positive, so it proves attraction for every positive drag coefficient.

There is, however, a boundary error in the asserted necessity. To have the bracket positive for every $\gamma >0$, its constant and slope as a function of ${\gamma }^2$ must be nonnegative, and cannot both vanish. Together with coefficient positivity this gives the exact drag-polynomial stability for all positive stopping rates:Necessity follows by taking arbitrarily small and arbitrarily large positive $\gamma$; if one bracket coefficient is negative, the inequality fails in the corresponding limit. Sufficiency follows because at least one coefficient is strictly positive. Thus $K=L>P$ and $K>L=P$ still give strict decay at each finite positive $\gamma$. They need not give a uniform stability margin in the zero- or infinite-drag limit.

In a Keplerian shearing sheet, $S=3\Omega /2$ andThe differences areHence $r>3$ establishes the requested strict chain and trapping conclusion. At $r=3$, however, $K=L={\Omega }^2$ and $P=9{\Omega }^2/16$, soAll coefficients remain positive. This is an explicit counterexample within the stated vortex model: $r=3$ also attracts particles for every finite positive drag coefficient, so the literal strict “if and only if” is false. The exact Keplerian all-positive-drag condition is $r\ge 3$. If $\gamma =0$ is included, the polynomial has zero and purely imaginary roots and there is no asymptotic attraction; “all drag coefficients” must mean positive drag. The local attraction claim presumes the particle remains in the patch, as is assured for sufficiently small perturbations of its centre.

Let the particle and carrier-fluid mass densities be ${\rho }_p$ and ${\rho }_f$, with identical density contrast for both species. At small particle Reynolds number, balance Stokes drag law against submerged weight:ConsequentlyAt large particle Reynolds number, use approximately constant drag coefficient $C_D$ and quadratic drag. Balancing $\frac{1}{2}{C}_D{\rho }_f\pi r^2{W}_s^2$ against the same submerged weight gives $W_s^2=8({\rho }_p-{\rho }_f)gr/(3C_D{\rho }_f)$, soBoth particles must lie in the constant-$C_D$ regime for this simple square-root ratio. More generally it contains an additional factor $\sqrt{C_D(r_2)/C_D(r_1)}$. The speeds are positive downward magnitudes; the signs of the shock velocities below depend on the chosen spatial coordinate.

Let $\Delta \rho ={\rho }_p-{\rho }_f>0$. A spherical grain has submerged weightNeglect lift and contact torque, and use a sliding static friction model with normal reaction $R=G$. The inertial quadratic drag iswhere the last equality fixes the convention $\tau ={\rho }_f{U}^2$. Equivalently $U$ is the bed shear velocity and $C_D$ is an effective drag coefficient referred to it. A literal grain-level flow speed can differ from the shear velocity; that conversion must then be absorbed into $C_D$.

Downstream sliding begins when $F_D>{\mu }_{s}R$. Defining the Shields parameter by $\Theta =\tau /(\Delta \rho gd)$, the threshold force balance isThe grain moves downstream above this threshold in the stated sliding model. Real grain motion can instead involve lift, rolling, irregular contacts, or viscous drag; the printed constant belongs to the particular inertial-drag convention and force balance above.

The normal contact force and resisting static friction balance the drag force and submerged weight at impending motion. The right panel uses locally bed-tangent drag force.

Use locally bed-tangent drag force, consistent with interpreting $\tau$ as bed shear stress. Let the bed rise downstream at angle $\alpha$ and retain the same effective drag coefficient. The normal reaction is $R=G\cos \alpha$, while the opposing downslope component of submerged weight is $G\sin \alpha$. At impending upslope motion,Thus the inclined-bed sediment threshold isFor $|\alpha |\ll 1$,An uphill slope increases the sediment entrainment threshold; a downhill slope lowers it, until spontaneous gravity-driven motion invalidates the resting-bed model.

There is a geometric convention to specify: if the drag force remains horizontal while the bed tilts, then $R=G\cos \alpha +F_D\sin \alpha$. The corresponding force balance gives insteadThe subsequent bed shear stress analysis uses the first, locally tangent convention. The two interpretations should not be mixed.

Descent equation. Take downward velocity $v$ and underwater depth $z$ as positive, with $z=0$ at the lake surface. For a sphere of radius $r$, the pebble mass is $m=4\pi r^3{\rho }_p/3$ and the frontal area is $A=\pi r^2$. Newton's second law includes downward gravitational weight, upward buoyancy and upward quadratic drag:ThusHere $b$ has units of inverse length. For arbitrary signed velocity, the last term is $-bv|v|$, so drag always opposes motion. The displayed downward form suffices for this trajectory.

Limiting speed. At terminal velocity, the net force and $\dot{v}$ vanish. For ${\rho }_p>{\rho }_w$,Use ${\rho }_p=2500\mathrm{kg}{\mathrm{m}}^{-3}$, ${\rho }_w=1000\mathrm{kg}{\mathrm{m}}^{-3}$, $r=0.020\mathrm{m}$, $C=0.7$ and $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$. ThenThe ${\rho }_p-{\rho }_w$ factor is essential: omitting buoyancy would overestimate the limiting speed.

For the numerical descent, assume the pebble is released from rest, neglect air drag, take $g$ and the water mass density as constant, and apply the given constant drag coefficient once it is submerged. Neglect the finite entry interval, entry energy losses and added mass, so its initial underwater speed isIt consequently slows down towards terminal velocity; a rest-starting underwater trajectory would describe a different initial condition.

An upward drag force $F_D\ge 0$ changes the variable-mass system equation toFor an ordinary drag force increasing with positive speed, the new curve lies below the no-drag curve. It still starts with slope $g$ if $F_D=0$ at rest. The dashed curve illustrates quadratic drag on a sphere with constant radial growth, using $F_D=\frac{1}{2}{C}_D{\rho }_{\mathrm{air}}\pi r^2{v}^2$ with constant drag coefficient $C_D$ and air mass density ${\rho }_{\mathrm{air}}$. The equation for these dimensionless quantities isThe figure uses $\eta =2$ solely to show the shape; it is not an additional datum of the paper.

There is a limitation to the usual terminal velocity sketch here. The examiner's commentary describes a curve approaching terminal velocity, but the paper maintains unbounded radius growth and specifies no drag law. A horizontal limit is justified if one additionally treats the radius as approximately fixed, or otherwise assumes a drag-to-mass ratio with a fixed limiting balance. With the growing radius and constant-$C_D$ quadratic drag just written, the late balance instead givesThis curve flattens in the sense that its slope tends to zero, but its speed remains unbounded. Thus a lower, more strongly flattening curve is warranted; a finite terminal speed requires an additional assumption. The no-growth approximation gives the familiar finite value $v_T=\sqrt{8\rho gr/(3C_D{\rho }_{\mathrm{air}})}$ for a specified fixed radius.^[ref]

Let $C(x)$ be the particle number density, $U(x)$ the conservative potential energy and $F=-dU/dx$ the applied force. In the overdamped regime a viscous drag coefficient $\zeta$ gives drift velocity $v_d=F/\zeta$. Combining this drift with Fick's first law gives the number fluxAt thermal equilibrium, the Boltzmann distribution is $C_{\mathrm{eq}}\propto e^{-U/(k_{\mathrm{B}}T)}$, so $d{C}_{\mathrm{eq}}/dx=F{C}_{\mathrm{eq}}/(k_{\mathrm{B}}T)$. Setting its net flux to zero yieldsThis is the Einstein diffusion relation. The derivation assumes a dilute equilibrium suspension, uniform temperature, linear drift response and a well-defined long-time diffusion coefficient. With a nonzero conservative force, the full transport equation includes drift; the force-free equation in the next section describes spreading without that drift.

Let $U(x)$ be the static potential and $\gamma$ the drag coefficient, so terminal Genetic drift under force $-U^{\prime }$ is $v=-U^{\prime }/\gamma$. The probability current is $J=-D{U}^{\prime }p/(k_{B}T)$ only after imposing the desired relation; instead start independently from $J=-D{p}^{\prime }-U^{\prime }p/\gamma$. At thermal equilibrium $p\propto e^{-U/(k_{B}T)}$ and detailed balance requires $J=0$. Substitution yieldsAt 298 K with $a=0.5\mu$m, this gives approximately $4.37\times 10^{-13}$ m${}^2$ s${}^{-1}$. The argument assumes an overdamped particle in an equilibrium thermal bath, rather than arbitrary driven diffusion.

The Stokes–Einstein relation connects thermal diffusion with dissipative mobility. For a dilute suspension of isolated rigid spheres of hydrodynamic radius $a$ in an unbounded Newtonian fluid of viscosity $\eta$,where $\zeta$ is translational drag coefficient, $D$ diffusion coefficient and $T$ absolute temperature. Stokes drag supplies $\zeta =6\pi \eta a$ under no-slip, continuum, low-Reynolds-number conditions. The Einstein diffusion relation $D=k_{B}T/\zeta$ is more general than that particular spherical drag law.

A short equilibrium argument explains the connection. Let a weak external potential $U$ exert force $\mathbf{F}=-\nabla U$, so drift velocity is $\mathbf{F}/\zeta$. The number flux is $\mathbf{j}=n\mathbf{F}/\zeta -D\nabla n$. A Boltzmann distribution $n\propto e^{-U/(k_{B}T)}$ has $\nabla n=n\mathbf{F}/(k_{B}T)$; requiring zero equilibrium flux fixes $D=k_{B}T/\zeta$. thermal fluctuations and viscous dissipation are therefore linked, rather than independent adjustable effects. The original theory also relates diffusion to measured random displacements.^[ref]

On times long compared with momentum relaxation, an isotropic diffusing sphere hassingle-particle tracking can measure the slope. Alternatively, dilute single-scattering dynamic light scattering measures a field correlation $g_1(q,t)=e^{-D{q}^{2}t}$; the corresponding intensity-correlation excess is proportional to $e^{-2D{q}^{2}t}$. The scattering wavenumber and the distinction between field and intensity decay must be retained.

A Stokes-Einstein experimental consistency test independently measures temperature and solvent viscosity, uses a narrow distribution of inert spheres, checks decay rate proportional to $q^2$ or displacement variance proportional to time, and extrapolates toward infinite dilution. Subtract drift, account for wall-enhanced drag, and verify that aggregation is absent. Tracking additionally needs spatial calibration and correction for finite exposure and localization precision. At finite concentration, collective diffusion is not automatically the single-particle self-diffusion in the equation. A hydrodynamic radius is not generally identical to a polymer radius of gyration; nonspherical particles, slip or viscoelastic solvents require a different mobility model.

There is a source error to resolve first. For a drag coefficient defined by $F_{\mathrm{drag}}=-\xi v$, the velocity relaxation time is $\tau =m/\xi$. The mass and drag coefficient in the PDF's exponential are interchanged: $\xi t/m$ is dimensionless, whereas $mt/\xi$ is not. The consistent mean relaxation law is therefore $⟨v(t)⟩=v_0{e}^{-t/\tau }$.

A mean relaxation law alone does not determine a variance. Use the usual thermal Langevin equation, with a bath at temperature $T$ and the fluctuation-dissipation relation for a Langevin particle:Here the calculation is for one component in Cartesian coordinates. With fixed initial position and velocity, integration gives the mean displacementIntegrating the stochastic part of the velocity once more gives the centered displacementThe delta correlation reduces its variance to a single integral:This is the integrated Ornstein-Uhlenbeck displacement for a prescribed initial velocity. Its short-time variance is $2D{t}^3/(3{\tau }^2)+O(t^4)$; at long times it is $2Dt+O(1)$. Thus the diffusion coefficient isIf the requested mean-squared fluctuation instead refers to the conventional equilibrium mean squared displacement, the initial velocity is also averaged over its thermal distribution, with $⟨v_0⟩=0$ and $⟨v_0^2⟩=k_{B}T/m$. Its contribution adds $D\tau (1-e^{-t/\tau }{)}^2$, givingThis has the ballistic limit $(k_{B}T/m)t^2$ and the same diffusive limit $2Dt$. For fixed nonzero $v_0$, the uncentered mean square instead equals the first boxed variance plus $v_0^2{\tau }^2(1-e^{-t/\tau }{)}^2$. In three independent spatial dimensions, the centered variances add, and the long-time result is $⟨|\mathbf{r}-⟨\mathbf{r}⟩{|}^2⟩\sim 6Dt$. Stating the initial-velocity convention is necessary to distinguish these finite-time formulas.

Quadratic drag has magnitude proportional to the square of the relative speed and points opposite to the relative velocity.

For an isotropic model of a body of fixed shape, let $\mathbf{U}={\mathbf{v}}_{\mathrm{body}}-{\mathbf{u}}_{\infty }$ denote its relative velocity through the distant fluid, with $U=|\mathbf{U}|$. A constant positive resistance $k$ givesHere density $\rho$, reference area $A$ and drag coefficient $C_D$ are held fixed over the modeled range of relative speeds. The vector formula makes the direction explicit: the drag force on the body opposes its motion relative to the fluid, and vanishes at $U=0$. For a general body, orientation and other flow conditions can affect the drag coefficient; treating it as a constant requires a regime in which that is an adequate approximation. ^[ref]

Compare two distinct models at a chosen relative speed $U_0>0$. A linear drag model has $F_D(U)/F_D(U_0)=U/U_0$, while a constant-coefficient quadratic drag model has $F_D(U)/F_D(U_0)=(U/U_0{)}^2$. At $U=2U_0$ their drag force ratios are two and four respectively.

This distinction matters when estimating terminal velocity or trajectories such as descent from rest under quadratic drag: doubling relative speed increases the resistance much more strongly in the quadratic drag model. An approximately constant drag coefficient can occur in another Reynolds number regime; its use must be justified there. Conversely, substituting the creeping-flow sphere coefficient $C_D=24/{\mathrm{Re}}_d$ into $F_D=\frac{1}{2}\rho U^{2}A{C}_D$ recovers linear drag, so writing this general formula does not make the resistance quadratic. Neither approximation should be extended across all Reynolds numbers.

The distinction already appears in the original investigation underlying the Stokes drag law, read on 9 December 1850: article 42 identifies resistance proportional to relative speed and distinguishes it from a contribution proportional to the square of relative speed. ^[ref]

The sediment entrainment threshold is the bed shear stress at which resting grains begin moving. An elementary force balance using spherical grains, tangential quadratic drag, and static friction gives ${\Theta }_{\mathrm{th},0}=4{\mu }_s/(3C_D)$ when the drag coefficient is defined relative to the bed shear velocity.