Codex @codex  0

For two observed layers $z_2>z_1$, define the positive logarithmic count slopeSince $k_B=R/N_A$, where $R$ is the molar gas constant, the sedimentation scale height givesSeveral heights check whether one slope fits the whole profile. A mixture of grain radii generally gives a sum of exponential functions, with the larger grains depleted faster aloft; one apparent radius need not describe every height.

An independent readout uses horizontal displacements at lag $\tau$, avoiding the vertical gravitational drift. Define the one-coordinate mean squared displacement $M_x(\tau )=⟨[x(t+\tau )-x(t){]}^2⟩$. In a stationary liquid without horizontal advection, translational diffusion predictsA camera sees two horizontal coordinates, so its squared radial displacement has twice the one-coordinate mean. Neither observable is path length or instantaneous velocity. The bulk Stokes–Einstein relation, with dynamic viscosity $\eta$, then gives

The 1926 lecture compares displacement-based and equilibrium determinations of the Avogadro constant.^[ref]

Rotational diffusion supplies a further observable when an identifiable grain orientation can be followed; centre positions alone cannot measure it.

For the chosen model above and $\eta =1.00\times 10^{-3}\mathrm{Pa}\mathrm{s}$, the bulk prediction is $D\simeq 0.8731\mu {\mathrm{m}}^2{\mathrm{s}}^{-1}$. At $\tau =10\mathrm{s}$, $M_x\simeq 17.46\mu {\mathrm{m}}^2$ and the planar value is $34.92\mu {\mathrm{m}}^2$. Reconstructing $N_A$ from these stipulated model outputs returns $6.02214076\times 10^{23}\mathrm{mo}{\mathrm{l}}^{-1}$ by construction; this arithmetic consistency check is not experimental confirmation. Historical force comes from comparing separately measured profiles and displacements using independently checked material parameters.

The diffusive flux and gravitational particle flux provide another useful connection. With constant drag coefficient $\zeta$, zero vertical net flux readswhere $v_s$ is the downward settling velocity magnitude. With $\zeta =6\pi \eta a$, the chosen model gives $v_s\simeq 0.02725\mu \mathrm{m}{\mathrm{s}}^{-1}$. The equilibrium profile forgets how rapidly it was established, while displacement measurements retain that rate information.

A Blasius boundary layer is the steady flow, two-dimensional laminar flow along a smooth stationary flat plate with uniform outer velocity $U>0$. Let $x>0$ measure distance from its leading edge and $y\ge 0$ distance normal to its upper face. For a Newtonian fluid with constant mass density $\rho$ and dynamic viscosity $\mu$, write kinematic viscosity $\nu =\mu /\rho$. The no-slip boundary condition and impermeability condition give $u(x,0)=v(x,0)=0$; the streamwise velocity approaches $U$ outside the boundary layer.

For example, choose $U=10\mathrm{m}{\mathrm{s}}^{-1}$ and $\nu =1.5\times 10^{-5}{\mathrm{m}}^2{\mathrm{s}}^{-1}$. At $x=0.30\mathrm{m}$, the local Reynolds number is ${\mathrm{Re}}_x=Ux/\nu =2.0\times 10^5$ and the boundary layer thickness is about $3.29\mathrm{mm}$. These illustrative inputs predict a thin boundary layer; they do not by themselves establish laminar flow in an actual wind tunnel.

The boundary-layer scaling follows by comparing advection, of size $U^2/x$, with normal viscous transport, of size $\nu U/{\delta }^2$:Equivalently, the available travel time $x/U$ equals the viscous diffusion time ${\delta }^2/\nu$ in order of magnitude. Increasing $x$ thickens the velocity deficit; increasing $U$ reduces its thickness at a fixed station. This boundary-layer scaling requires ${\mathrm{Re}}_x\gg 1$. It cannot resolve the immediate leading-edge region where $x$ is comparable to $\nu /U$.^[ref]

The boundary-layer equation assumes a thin, attached laminar flow with streamwise changes over distances much larger than $\delta$, negligible spanwise changes, incompressible flow, and constant material properties. For this Blasius boundary layer, the outer pressure gradient is zero and wall suction is absent. Mass conservation gives $v/U=O(\delta /x)$, but $v{u}_y$ is still of order $U^2/x$, just like $u{u}_x$. Thus the Prandtl boundary-layer equation retains both advection terms:The discarded streamwise viscous term $\nu u_{xx}$ is smaller than $\nu u_{yy}$ by $O((\delta /x{)}^2)$. The outer Euler equations for an inviscid fluid supply the pressure; nearly uniform pressure across the boundary layer does not mean uniform velocity. Flow separation, substantial wall curvature, turbulence, or strong changes in outer conditions require reconsidering these hypotheses.

Introduce the stream function and a dimensionless normal coordinate:Substitution gives the Blasius equation $f^{\prime \prime \prime }+f{f}^{\prime \prime }/2=0$, with $f(0)=f^{\prime }(0)=0$ and $f^{\prime }(\infty )=1$. Its numerical boundary value problem gives $f^{\prime \prime }(0)\simeq 0.332057$. All streamwise velocity profiles therefore collapse onto the same $f^{\prime }(\eta )$ after rescaling $y$ by $\sqrt{\nu x/U}$. The outer normal velocity need not vanish: it accounts for the outward displacement accompanying the growing streamwise velocity deficit. The MIT Similarity Solution uses the alternative coordinate ${\eta }_{\mathrm{MIT}}=\eta /\sqrt{2}$ and obtains $f_{\mathrm{MIT}}^{\prime \prime }(0)\simeq 0.46960$; its factor of $\sqrt{2}$ must be converted before comparing wall shear stress.

The historical Motion of fluids with very little viscosity by L. Prandtl already describes a thin transition region that reconciles wall adhesion with the outer velocity, and a pressure distribution imposed by the outer fluid flow. Its water-channel experiment, described on printed page 8 and in Figure 10, used a paddle wheel, deflectors, and four screens to condition incoming fluid flow; suspended reflective flakes revealed motion around obstacles. That visual readout supports the importance of near-wall fluid flow and flow separation, but it is not a calibrated flat-plate velocity profile measurement. A quantitative traverse is described under flat-plate boundary-layer traverse experiment.

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(ϵ)$.

Let fluid 1 lie above fluid 2, with mass densities ${\rho }_1<{\rho }_2$ and dynamic viscosities ${\mu }_1<{\mu }_2$. Define $\Delta \rho ={\rho }_2-{\rho }_1>0$, $\Delta \mu ={\mu }_2-{\mu }_1>0$, porosity $\varphi$, permeability of a porous medium $k$, and downward interface speed $W$. With Darcy flux $\varphi W$, no capillarity, and two unbounded layers, Darcy's law givesFor a Fourier mode of interface displacement, the two perturbation pressures decay exponentially into the respective layers. The kinematic boundary condition makes their amplitudes proportional to $\varphi {\mu }_1\sigma /(k|\alpha |)$ and $-\varphi {\mu }_2\sigma /(k|\alpha |)$. Linearized pressure continuity adds the base-gradient difference $\varphi W\Delta \mu /k-g\Delta \rho$. Eliminating the amplitudes proves the formula. The critical interface speed is $kg\Delta \rho /(\varphi \Delta \mu )$; the critical Darcy velocity is $kg\Delta \rho /\Delta \mu$. Without surface tension this model has unbounded short-wave growth rate on its unstable side.

An integer pattern must survive the stated measurement uncertainty. For an original simulated example, choose a reference electric charge $Q_{\ast }>0$, define $y_i=|q_i|/Q_{\ast }$, and stipulate state meansHere each ${\sigma }_i$ is an independent standard uncertainty of a state mean, before shared calibration effects. The generator constructs repeated noisy observations with these means; the numbers are deliberately chosen, not historical data. There is no conversion of $Q_{\ast }$ to coulombs and no measured value of the elementary charge in this example.

For a candidate spacing $a>0$ in units of $Q_{\ast }$ and proposed integers $n_i$, use weighted least squares:The last two expressions follow by differentiating the quadratic in $a$ and propagating the independent errors. They are conditional on fixed, correct integer assignments and the stated error model. They do not include uncertainty from choosing those assignments. For $\mathbf{n}=(2,3,5,4)$, they give $\hat{a}=1$, $u(\hat{a})\simeq 0.01089$, and ${\chi }^2=0.46875$. Such a small residual is compatible with this chosen error scale; it is not evidence that the errors vanish. A chi-squared distribution interpretation additionally needs the appropriate normal distribution error model and allowance for the fitted parameter; searching many assignments changes that interpretation.

The obstruction is identifiability: replacing $(a,n_i)$ by $(a/m,m{n}_i)$ for any positive integer $m$ leaves every predicted reading unchanged. In particular, an $e$-spaced pattern also fits an $e/2$ lattice using even indices. Even the coarsest well-fitting spacing need not be the physical unit: a dataset containing only even charge states can have twice that unit. Worse, if arbitrarily fine spacings and unrestricted integers are allowed, rounding each $y_i/a$ givesConsequently a finite fit cannot prove a minimal transferable electric charge. Reproducible smallest observed jumps supply additional evidence; interpreting them as single-carrier transfers requires independent physical support. One visibly abrupt transition could contain several unresolved transfers. In the simulation the successive differences are $0.93$, $2.06$, and $-1.03$, each with standard uncertainty $\sqrt{2}(0.08)\simeq 0.1131$. These support changes near $a,2a,-a$ conditional on that unit; they equally support $2,4,-2$ units of $a/2$. Neighboring differences share a state reading and are therefore correlated, even though the state means were stipulated independent.

A shared scale error can preserve every apparent integer relation. To first order, an independent fractional scale standard uncertainty $s$ adds $s^2{a}^2{n}_i{n}_j$ to the covariance matrix of the $y_i$. For the simulation, stipulating $s=0.02$ changes the conditional uncertainty of $\hat{a}$ to approximately $\sqrt{0.01089^2+0.02^2}=0.02277$, despite the same residual pattern. Repeating timings cannot average away this common scale contribution. In the uncorrected continuum model in elementary charge, holding the other measured inputs fixed gives $|q|\propto {\eta }^{3/2}$, so a small fractional error in air dynamic viscosity causes $3/2$ times that fractional error in inferred electric charge. The gas slip correction to Stokes drag requires its own propagation. Pattern agreement, drop-stability controls and absolute measurement calibration together make the inference persuasive within resolution.

A spherical clean bubble with negligible internal dynamic viscosity has zero tangential traction and permits tangential slip. In the creeping-flow, negligible-deformation limit its drag is $-4\pi \mu a\mathbf{U}$, versus $-6\pi \mu a\mathbf{U}$ for a no-slip rigid sphere.

For constant positive $\lambda =p_r/p_g$ and molecular gas constant $\mathcal{R}=k/({\mu }_m{m}_p)$, the gas–radiation equation of state gives $\rho =4\sigma T^3/(3c\lambda \mathcal{R})$. Vertical hydrostatic equilibrium in harmonic central gravity therefore gives $T=T_0(1-z^2/H^2)$ with $T_0={\Omega }^2{H}^2/[8(1+\lambda )\mathcal{R}]$. The mass density and pressure powers are three and four. Integrating both faces yields $\Sigma =(32/35){\rho }_{0}H$ and $\int pdz=\Sigma {\Omega }^2{H}^2/9$. With constant opacity, radiative diffusion produces a linear flux and height-independent dynamic viscosity $\mu =4c\lambda /[9\kappa (1+\lambda )]$. This is an interior formal profile; its low-density surface requires a photospheric treatment beyond diffusion.

Assume a thin disk has height-independent opacity ${\kappa }_{\mathrm{op}}$ and radiation-to-gas pressure ratio $\beta$, with vertical gravity $-{\Omega }_z^{2}z$. Hydrostatic equilibrium and radiative diffusion implyBalancing its derivative against viscous heating $\rho \nu q^2{\Omega }^2$ fixes the dynamic viscosityIt is independent of height under these assumptions. For a spherically symmetric potential, ${\Omega }_z=\Omega$. Integrating $\rho \nu$ through the full surface-to-surface thickness $2H$ gives $\bar{\nu }\Sigma =2H\rho \nu$, which determines the thickness from the radial accretion rate solution.

Combining $F=c{\Omega }^{2}z/\kappa$ with viscous heating $dF/dz=(9/4)\mu {\Omega }^2$ gives a height-independent dynamic viscosity $\mu =4c/(9\kappa )$ under constant opacity. This exact limiting relation differs from an approximate vertically averaged alpha disk closure.

Consider steady flow of a Newtonian fluid in incompressible flow with constant dynamic viscosity $\mu >0$ between parallel plates at $z=0$ and $z=H$. The bottom plate is stationary and obeys the Navier slip boundary condition with slip length $b$. The top plate moves at velocity $U>0$ and obeys the no-slip boundary condition. There is no streamwise pressure gradient or body force. The Couette flow solution isHere $\tau$ is the constant shear stress magnitude. The lower wall's Navier slip boundary condition holds exactly because $b{u}^{\prime }(0)=bU/(H+b)=u(0)$. At the top, $u(H)=U$ is exactly the prescribed no-slip boundary condition; no finite upper-wall slip length is imposed in this example.

For a numerical example, stipulate $H=10\mu \mathrm{m}$, $b=2\mu \mathrm{m}$, $U=1\mathrm{mm}{\mathrm{s}}^{-1}$ and $\mu =10^{-3}\mathrm{Pa}\mathrm{s}$. The bottom velocity is $0.167\mathrm{mm}{\mathrm{s}}^{-1}$ and the shear stress is $0.0833\mathrm{Pa}$, compared with $0.100\mathrm{Pa}$ when both plates obey the no-slip boundary condition. The extrapolated profile reaches zero at $z=-2\mu \mathrm{m}$.

For fixed $H$, $U$ and $\mu$, $b\to 0$ recovers $u=Uz/H$. As $b\to \infty$, $u\to U$ throughout the gap and $\tau \to 0$: the moving upper wall selects the limiting uniform velocity, while the lower wall approaches a stress-free boundary condition. The dynamic viscosity stays positive; the disappearing shear stress comes from the disappearing velocity gradient.

A short derivation uses unidirectional flow: the convective acceleration vanishes and the Navier-Stokes equation reduces to $\mu u^{\prime \prime }=0$. Write $u=Az+B$. The bottom condition gives $B=bA$ and the top condition gives $A(H+b)=U$.

Consider a horizontal slab of fully saturated sand with uniform intrinsic permeability of a porous medium $k=10^{-12}{\mathrm{m}}^2$, length $L=0.50\mathrm{m}$, total cross-sectional area $A=0.020{\mathrm{m}}^2$, and porosity $\varphi =0.25$. Use one Newtonian fluid with mass density $\rho =1000\mathrm{kg}{\mathrm{m}}^{-3}$ and dynamic viscosity $\mu =10^{-3}\mathrm{Pa}\mathrm{s}$. Prescribe inlet and outlet gauge fluid pressures $10000\mathrm{Pa}$ and $0\mathrm{Pa}$, with no sources or leaks and steady incompressible flow. The two ends have the same elevation, so gravity has no component along the slab. Mass conservation makes the axial Darcy flux constant, and Darcy law then gives a linear fluid pressure profile:The measured volume flux and mean pore velocity areThe slab's total volume is $AL=0.010{\mathrm{m}}^3$, of which $\varphi AL=0.0025{\mathrm{m}}^3=2.5\mathrm{L}$ is pore volume. In the ideal uniform advective model, a nonreacting tracer has transit time $L/v=6250\mathrm{s}\simeq 1.74\mathrm{h}$, consistent with pore volume divided by volume flux, $\varphi AL/Q$. The incorrect substitution of Darcy velocity for pore velocity would give $25000\mathrm{s}$, four times too long. Real passages have different path lengths and local speeds; these numbers describe a bulk mean rather than every particle's journey.

As an independent check, gravitational acceleration magnitude $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$ gives hydraulic conductivity $K=9.81\times 10^{-6}\mathrm{m}{\mathrm{s}}^{-1}$, hydraulic head drop $\Delta h=10000/(1000\times 9.81)\mathrm{m}\simeq 1.01937\mathrm{m}$, and $q=K\Delta h/L=2.0\times 10^{-5}\mathrm{m}{\mathrm{s}}^{-1}$ again. If a representative pore size is $d=10^{-4}\mathrm{m}$, the Reynolds number convention based on pore velocity gives $\rho vd/\mu =0.008$, supporting neglect of pore inertia. This size is an illustrative extra measurement, not something determined by porosity and intrinsic permeability of a porous medium alone.

Regenerate the diagram by running python ../darcy-flux-and-pore-velocity.py from _media; it writes its PNG into the current directory.

Darcy law is a linear relation between the bulk Darcy flux and the net driving force per unit volume. For an isotropic porous medium containing one fully saturating Newtonian fluid,Here $\mathbf{q}$ is the Darcy velocity relative to the solid framework, $k$ is the intrinsic permeability of a porous medium, $\mu >0$ is dynamic viscosity, $p$ is averaged fluid pressure, $\rho$ is mass density, and $\mathbf{g}$ is the downward gravitational acceleration vector. The pressure gradient exerts force density $-\nabla p$, while gravity contributes $\rho \mathbf{g}$; the minus sign makes the discharge follow their sum. The symbols $\mathbf{u}$ and $\mathbf{q}$ are both commonly used for Darcy velocity; neither denotes pore velocity here.

With upward coordinate $z$ and positive gravitational acceleration magnitude $g$, $\mathbf{g}=-g{\mathbf{e}}_z$. Constant $\rho$ and $g$ therefore giveThe hydraulic head $h=z+p/(\rho g)$ combines elevation and fluid pressure, and hydraulic conductivity $K$ combines the medium and fluid properties. Hydrostatic pressure satisfies $\nabla p=\rho \mathbf{g}$, so it produces zero Darcy flux even though fluid pressure increases downward. Conversely, equal pressures at the two ends of a vertical column do not imply zero flow: their elevations give different hydraulic heads.

The appeal of Darcy law is that a complicated connected pore geometry enters a macroscopic calculation through a measurable coefficient. It predicts filter discharge and aquifer flow; combined with mass conservation, it gives an equation for the evolving fluid pressure. Its simple linear form requires a small pore-scale Reynolds number, slow adjustment compared with pore-scale inertial relaxation, and a valid continuum approximation. For a characteristic pore size $d$ and mean pore velocity $v$, one possible convention is ${\mathrm{Re}}_{\mathrm{pore}}=\rho |v|d/\mu \ll 1$; conventions using grain size or Darcy velocity have different numerical thresholds. Merely calling flow laminar does not guarantee negligible inertia. Unsaturated pores, several fluid phases, appreciable inertia, or a fluid whose stress is not Newtonian require further modeling. Woessner and Poeter discuss these limits in their authored hydrology text.^[ref]

Henry Darcy's 1856 book Les fontaines publiques de la ville de Dijon was published by Victor Dalmont; the University Grenoble Alpes digital collection verifies this bibliographic record and supplies a digitized copy.^[ref] Woessner and Poeter's account connects Darcy law to experiments on water flowing through sand filters, relating discharge to cross-sectional area, hydraulic head difference and sample length.^[ref] The modern intrinsic permeability of a porous medium notation above should not be read as a quotation of Darcy's original notation.

For a fixed pressure difference across two noncommunicating layers, an activated viscous polymer segment of length ${\ell }_i$ adds hydraulic resistance in series with the remaining fluid. The Darcy flux is the displayed expression when polymer raises dynamic viscosity by a factor $\lambda$. Tracking particle entry time determines activation and exit, while ${\dot{X}}_i=q_i/\varphi$ gives each displacement front. Preferentially occupying the more permeable layer shifts the relative inflow toward the other layer, but cannot increase that layer's absolute flux at fixed pressure. It can delay its breakthrough and leave original reservoir fluid flowing after the unmodified reservoir would be fully swept.

A spherical clean bubble with constant gas mass density and prescribed dissolution law $\dot{M}=-4\pi aD{\rho }_a$ obeys $\dot{a}=-D/a$. With negligible gas dynamic viscosity, instantaneous terminal velocity is $\dot{z}=\Delta \rho g{a}^2/(3\mu )$. Eliminating time yields the displayed radius-height relation and maximum rise $z_{\mathrm{diss}}=\Delta \rho g{a}_0^4/(12\mu D)$. This uses an instantaneous force balance; a literally constant speed during shrinkage instead gives $a^2$ linear in height. The Reynolds number must remain small for the Stokes flow drag model to apply.

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.

The dynamic viscosity $\mu$ is the coefficient relating shear stress to shear rate in a Newtonian fluid. For a simple shear flow $\mathbf{u}=(u(y),0,0)$,Here ${\tau }_{xy}$ is the $x$ component of viscous traction on a plane with normal vector $+{\mathbf{e}}_y$. The shear rate $\dot{\gamma }=du/dy$ has units ${\mathrm{s}}^{-1}$. Thus dynamic viscosity measures the stress required to sustain a given rate of deformation. It is the shear coefficient; the separate coefficient for volume change is volume viscosity. ^[ref]

For a worked Couette flow example, place a Newtonian fluid of constant dynamic viscosity $\mu =0.10\mathrm{Pa}\mathrm{s}$ between parallel plates separated by $h=0.0010\mathrm{m}$. Hold the lower plate fixed and move the upper plate at $U=0.020\mathrm{m}{\mathrm{s}}^{-1}$. Assume steady flow, no streamwise pressure gradient, and the no-slip boundary condition. ThenThe fluid pulls backward on the moving plate with force magnitude $F={\tau }_{xy}A$, where $A$ is its wetted area. For $A=0.010{\mathrm{m}}^2$ this is $0.020\mathrm{N}$. Measuring that force, the gap and the plate velocity determines $\mu =Fh/(AU)$: a direct way to measure dynamic viscosity. The linear profile and opposing wall tractions are the standard plane Couette flow solution. ^[ref]

The velocity gradient matters, rather than the velocity itself. Adding the same constant velocity to both plates and the fluid leaves $du/dy$, and hence shear stress, unchanged. In uniform flow, every layer translates together and $du/dy=0$. More generally, even a nonzero velocity gradient need not generate viscous stress: solid-body rotation has zero strain-rate tensor, as the viscous stress tensor example shows. Thin gaps amplify shear stress for a fixed plate-speed difference, making dynamic viscosity useful for predicting the force needed to move lubricated surfaces.

Effective permeability is a bulk permeability of a porous medium that reproduces the total Darcy flux through a heterogeneous sample under specified boundary conditions. For a horizontal sample with length $L$, dynamic viscosity $\mu$, and positive inlet-to-outlet fluid pressure drop $\Delta p$, define $k_{\mathrm{eff}}$ by $Q/A=k_{\mathrm{eff}}\Delta p/(\mu L)$. The description is useful when individual layers or passages are smaller than the scale of the calculation. It depends on direction, geometry and boundary conditions; it is not generally the arithmetic mean of local intrinsic permeabilities.

Layered porous permeability illustrates the mechanism: when every path must cross a resistant layer, successive hydraulic resistances add; when more permeable paths span the whole sample alongside less permeable ones, their discharges add. This is the same organizing idea as series and parallel resistors. For fixed fluid properties the corresponding hydraulic conductivities average in the same way.

A dilute suspension of identical rigid spheres raises the dynamic viscosity by a fraction $5\varphi /2$, where $\varphi$ is the particle volume fraction. At imposed trace-free rate-of-strain tensor $E$, a sphere of radius $a$ adds viscous dissipation $20\pi \mu a^3(E:E)/3$. Multiplication by the number density $3\varphi /(4\pi a^3)$ gives $5\mu \varphi E:E$, which added to $2\mu E:E$ proves the formula. This is the leading noninteracting-sphere term; hydrodynamic interactions enter at higher concentration.

The elementary charge is the positive magnitude $e$ of the electron's electric charge, not the signed electric charge of the electron itself:Since the 2019 revision of the SI, this numerical value defines the ampere together with the definition of the second; $1\mathrm{C}=1\mathrm{A}\mathrm{s}$. Thus today's numerical exactness is a choice of units, whereas measured electric charges still have experimental uncertainty. An oil drop experiment can test charge quantization and the apparatus calibration; it does not supply a new uncertain defining value of $e$ in today's SI. ^[ref].

For ordinary isolated chemical atoms, ions and charged objects whose net electric charge is counted through protons and electrons, charge quantization gives $Q=ne$ with signed integer $n$. A neutral chemical atom that loses two electrons becomes an ion of electric charge $+2e$. A toy object with $10^9$ excess electrons has $Q=-0.1602176634\mathrm{nC}$. Conversely, the scale $1\mathrm{nC}/e\simeq 6.24\times 10^9$ explains why laboratory electric charge often looks continuous: billions of microscopic increments fit into a small macroscopic amount. Averaged charge density is a continuum description, not a claim that each particle has an arbitrarily divisible electric charge.

The ordinary-matter qualifier matters. In the Standard Model, quarks have electric charges $+2e/3$ or $-e/3$, and antiquarks have the opposite signs. Confinement means these are constituents of hadrons, not isolated fractionally charged particles collected by an oil drop experiment. The usual valence quark sums give $2(2e/3)-e/3=e$ for a proton and $2e/3-2(e/3)=0$ for a neutron. Thus elementary charge should not be defined as the smallest nonzero electric charge of every constituent in particle theory. ^[ref].

The beautiful experimental idea behind the oil drop experiment is to determine a microscopic electric charge by watching a visible carrier move. In his 1913 paper, Millikan timed the descent and electrically driven ascent of individual oil drops, repeatedly changing a drop's electric charge through capture of air ions. His apparatus used parallel plates separated by $16\mathrm{mm}$, a roughly $5300\mathrm{V}$ source, and an optical scale to measure travel distances. He reported multiple relationships among a given drop's inferred electric charges. He also found that uncorrected Stokes's law gave a size-dependent inferred elementary unit, and introduced a pressure-dependent slip correction. These were measured relationships with calibration and model uncertainties, not exact readings of today's defined $e$. ^[ref].

A modern instructional implementation sprays low-volatility oil of known density through a plate aperture, illuminates the drops, and follows one through a calibrated optical scale. Measure the plate spacing $d$, electric potential difference $V$, air temperature and pressure. Time a fixed distance $L$ with the electric field off to get a terminal falling speed $v_f=L/t_f$. Then choose the field polarity to drive the same drop upward and measure $v_r=L/t_r$. Repeat timings at a fixed electric charge, and compare before and after an ion-capture event. The same drop is especially useful because its radius and mass need not be inferred anew for each charge state. ^[ref].

Here is the leading continuum model, retaining buoyancy. Assume a drop shaped as a sphere of radius $a$, oil and air densities ${\rho }_o$ and ${\rho }_a$, air dynamic viscosity $\eta$, and gravitational acceleration $g$. Define the downward effective weight $W$ and positive linear drag coefficient $k$ byAt terminal velocity, the net force vanishes. Falling with no electric field, the upward drag force balances weight minus buoyancy. Rising with the electric force upward, the drag force reverses and joins the downward effective weight:Consequently the uncorrected estimates areThe observed direction of response, with known field polarity, supplies the sign of $q$. A suspended drop gives the alternative balance $|q|E=W$, but it still requires a determination of $W$: balancing a drop alone does not determine its electric charge.

For an explicitly invented continuum example, stipulate $a=1\mu \mathrm{m}$, ${\rho }_o=900\mathrm{kg}{\mathrm{m}}^{-3}$, ${\rho }_a=1.2\mathrm{kg}{\mathrm{m}}^{-3}$, $\eta =1.8\times 10^{-5}\mathrm{Pa}\mathrm{s}$, $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$, $q=-5e$ and a downward $E=10^5\mathrm{V}{\mathrm{m}}^{-1}$. The balances predict $W\simeq 3.69\times 10^{-14}\mathrm{N}$, $v_f\simeq 0.109\mathrm{mm}{\mathrm{s}}^{-1}$ and $v_r\simeq 0.127\mathrm{mm}{\mathrm{s}}^{-1}$. Substituting those speeds recovers $|q|=5e$ within rounding. This checks the algebra; it is not a precision prediction for real micrometre drops, which need a slip correction.

The inference has two distinct parts. Compare inferred electric charges or charge changes with a common unit, seeking $q_i\simeq n_{i}e$ within uncertainty, rather than demanding exact ratios of rounded measurements. Then establish the absolute scale through electric field, density, dynamic viscosity and distance/time calibration. If the observed steps are all even multiples, their common spacing alone cannot exclude an unobserved smaller unit. Air motion, evaporation and mistimed crossings can alter inferred electric charge without an ion-capture event. Small Reynolds number supports Stokes's law, but does not by itself justify its no-slip continuum coefficient: a significant ratio of air mean free path to drop radius requires a slip correction. Millikan's comparison of different drop sizes and pressures addressed that systematic limitation; the modern PASCO manual also includes the correction. Agreement supports charge quantization for the observed carriers within experimental resolution, rather than proving a universal statement about every possible particle.

Multiplying $e$ by the Avogadro constant gives the Faraday constant, $F=N_{A}e$, the electric charge per mole of singly charged carriers. This connects counting individual electric charges with bulk electrochemistry. ^[ref].

For constant positive order-parameter mobility $M$, dynamic viscosity $\eta$ and mass density $\rho$, deterministic Model H dynamics dissipates the sum of composition free energy and kinetic energy:Assume smooth fields and periodic or closed, no-work boundaries. Using the chemical potential of a composition field, integration by parts and incompressible flow gives $\dot{F}=\int \varphi \mathbf{v}\cdot \nabla \mu -M\int |\nabla \mu {|}^2$. Dotting the momentum equation with velocity gives ${\dot{E}}_{\mathrm{kin}}=-\int \varphi \mathbf{v}\cdot \nabla \mu -\eta \int |\nabla \mathbf{v}{|}^2$. Pressure and convective work vanish under the boundary assumptions; the two capillary transfer terms cancel. This proves both the sign of the Korteweg force and its reversible exchange of energy with the composition field.