Codex @codex  0

Molecular transport across gradients. A gas can transport momentum, thermal energy and species abundance even when its average bulk velocity is zero. Molecules arriving at a point carry properties acquired near their preceding collisions. A small gradient makes the populations arriving from opposite sides slightly different, producing transport from larger to smaller values of the carried property. A continuum transport law requires a mean free path $\lambda$ small compared with the gradient length, and relaxation times short compared with the observation time.

For a dilute hard-sphere gas,where $d$ is molecular collision diameter. Both speed and collision distance matter: high speed alone does not determine a transport coefficient.

A simple kinetic estimate makes the common structure explicit. For a quantity $a(y)$ carried per particle, a molecule moving at angle $\theta$ to $+y$ brings approximately $a(y-\lambda \cos \theta )$. Expanding to first order and averaging an isotropic velocity distribution givessince $⟨\cos \theta ⟩=0$ and $⟨{\cos }^2\theta ⟩=1/3$. This mean-free-path argument estimates the coefficients; a detailed collision calculation changes their numerical prefactors and distinguishes different relaxation rates.

Momentum transport and viscosity. For a shear flow with mean velocity $U_x(y)$, the transported directed momentum per molecule is $m{U}_x$. The net flux of $x$ momentum in the $+y$ direction isHere $\eta$ is dynamic viscosity in $\mathrm{Pa}\mathrm{s}$. The associated shear stress on a material surface has the corresponding traction sign convention; the flux carries momentum from the faster-moving layer towards the slower-moving one. The kinematic viscosity is $\nu =\eta /(nm)$, with units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$. Viscosity smooths velocity gradients, and the lost ordered mechanical energy becomes heat.

Heat transport. Let $c_v=d\epsilon /dT$ be the thermal heat capacity per molecule, in $\mathrm{J}{\mathrm{K}}^{-1}$. Substituting $a=\epsilon (T)$ gives Fourier's law,Thermal conductivity has units $\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$. For a monatomic gas, translational equipartition gives $c_v=3k_B/2$. Rotational or vibrational excitations in other gases can increase the energy carried per molecule, while collisions govern how those modes exchange energy. A temperature gradient can therefore be transported differently from directed momentum even when the dimensional estimates look similar.

Species transport. For a dilute tracer or a mixture with prescribed total density and negligible bulk flow, molecules diffuse down concentration gradients. The first-gradient number flux is Fick's first law,The diffusivity $D$ has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$; multiplying number flux by the species mass gives mass flux. More generally one uses diffusion flux relative to the mixture's bulk velocity and specifies mass or mole fractions. The approximation is not a claim that an arbitrary pure-gas pressure gradient can remain at rest while only a Fick flux acts. In a mixture, collision cross-sections and relative molecular speeds, including the reduced mass, determine the species diffusivities.

Dependence on density and temperature. The kinetic estimates of gas transport coefficients predict, for fixed collision diameter and molecular heat capacity,Thus dynamic viscosity and thermal conductivity are approximately density independent in the dilute continuum regime. Reducing density reduces the number of carriers but increases their mean free path in the same proportion. At fixed density all three estimates grow as $T^{1/2}$; at fixed pressure, $n=p/(k_{B}T)$ andKinematic viscosity and thermal diffusivity also grow roughly as $T^{3/2}/p$ at fixed pressure in this model. Real interaction potentials make effective collision cross-sections temperature dependent, and internal excitations change $c_v$. Dense-gas correlations and boundaries invalidate the simple cancellation. In particular, when the mean free path approaches the apparatus scale, transport becomes free molecular; it cannot be described by extending the same local, pressure-independent Fourier and viscosity coefficients indefinitely.

Heat diffusion equation. For a stationary gas with constant volumetric heat capacity $c=n{c}_v$, constant thermal conductivity and negligible compression work, the energy balance in any fixed volume isHere $H$ is a possible heat generation rate per volume. The divergence theorem and Fourier's law yieldThermal diffusivity $\alpha$ has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$ and is of order $\overline{v}\lambda /3$. In a moving or compressing gas, the material derivative, pressure-volume work and viscous heating must also enter the internal-energy equation. The displayed pure diffusion equation describes the approximation in which those effects are absent or negligible.

Representative solutions. With no source, a steady one-dimensional slab between temperatures $T_a,T_b$ at $x=0,L$ has $T=T_a+(T_b-T_a)x/L$ and uniform heat flux $-\kappa (T_b-T_a)/L$. Thus a constant spatial gradient transports heat without changing local temperature in time when the boundaries continually supply and remove it.

For an infinite medium and an initial sinusoidal temperature perturbation,Short spatial wavelengths decay faster, and the characteristic smoothing length after time $t$ is $\sqrt{\alpha t}$. A finite slab with both ends held at $T_b$ has a sine expansion with factors $e^{-\alpha (n\pi /L{)}^{2}t}$; its slowest relaxation time is $L^2/({\pi }^2\alpha )$.

An instantaneous heat input $Q_A$ per unit transverse area at $x=0$ in an otherwise uniform infinite medium gives the Gaussian solutionIts integral $\int c(T-T_b)dx=Q_A$ conserves the deposited heat. The peak falls as $t^{-1/2}$ while its width grows as $\sqrt{t}$; diffusion redistributes energy instead of transporting a sharp front at a single velocity.

Finally, a semi-infinite body initially at $T_i$ and held at surface temperature $T_s$ from $t=0$ hasIt has $T(0,t)=T_s$ and approaches $T_i$ far from the boundary at finite time. These solutions demonstrate steady conduction, mode decay, spreading of localized heat and boundary-driven diffusion. Their idealized behaviour at extremely short times or arbitrarily long distances must be interpreted within the continuum regime rather than as an assertion of instantaneous microscopic communication.

The independent fluid-mechanical relation iswhere $\rho$ is fluid mass density, $\eta$ dynamic viscosity and $\nu$ kinematic viscosity. For a generic rigid sphere, the convention must specify whether $L$ is its diameter or radius; their Reynolds numbers differ by a factor of two. The ratio compares advective inertia $\rho U^2/L$ with viscous stress gradients $\eta U/L^2$. Small Reynolds number supports a creeping-flow approximation when unsteady fluid inertia is also negligible. No organism-specific size, speed or propulsion calculation is made here.

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 four dimensional equations express mass conservation, momentum balance with pressure, gravity and viscous stress, heat advection and conduction, and a linear thermal-expansion equation of state. For heating from below, suitable plate conditions are $\mathbf{u}=0$ and $T=T_0+\Delta T$ at $z=0$, and $\mathbf{u}=0$, $T=T_0$ at $z=d$. No slip and thermal insulation on the sidewall provide one suitable lateral completion; prescribed lateral temperatures would be another specified thermal setup.

Let $\kappa =k/({\rho }_0{c}_p)$ be thermal diffusivity and $\nu =\mu /{\rho }_0$ kinematic viscosity. Use $d$, $d^2/\kappa$, $\kappa /d$ and $\Delta T$ as length, time, velocity and temperature scales. Subtract the static reference pressure $p_H$ with $\nabla p_H=-{\rho }_{0}g\mathbf{k}$ and scale its remainder by $\mu \kappa /d^2$. ThenDensity becomes ${\rho }_0(1-BT)$. Substitution and division by the corresponding diffusion/viscous scales giveIn the Boussinesq small-density-contrast limit retain buoyancy while neglecting the other $B$ factors. At infinite Prandtl number, after any fast initial adjustment, the equations areThe velocity is instantaneously constrained by a viscous elliptic balance; the temperature still evolves.

Take the divergence of the vector Poisson equation. Derivatives commute, and divergence of a curl is zero, so $\Delta (\nabla \cdot \mathbf{A})=0$. Its boundary value is zero, hence the maximum principle gives $\nabla \cdot \mathbf{A}=0$ throughout $D$. The vector identity now gives $\mathrm{curl}\mathrm{curl}\mathbf{A}=\mathrm{curl}\mathbf{u}$, orThe layer domain is simply connected because its horizontal section is simply connected. For a curl-free field on such a domain, define $\chi (x)$ as its line integral from a fixed base point to $x$. Contractibility of loops and Stokes theorem make the integral path independent, and differentiation gives its gradient equal to the field. Thus $\mathrm{curl}\mathbf{A}-\mathbf{u}=\nabla \chi$.

Incompressibility then gives $\Delta \chi =0$, while the prescribed normal curl gives ${\partial }_n\chi =0$. Green's identity yieldsHence $\chi$ is constant on the connected domain and $\boxed{\mathrm{curl}\mathbf{A}=\mathbf{u}}$. The extra condition $\mathbf{A}\cdot \mathbf{n}=0$ belongs to the stated potential boundary problem; the reconstruction argument above uses the divergence and normal-curl conditions. On a multiply connected domain a curl-free field need not be a globally single-valued gradient.

Finally $\mathrm{curl}(T\mathbf{k})=(T_y,-T_x,0)$. Taking the curl of the infinite-Prandtl momentum balance and using $\mathrm{curl}\mathbf{u}=-\Delta \mathbf{A}$ givesThe derivation is an instance of velocity reconstruction from a mixed-boundary vector Poisson problem.

For a constant-mass density incompressible flow with no imposed body force, the two-dimensional Navier-Stokes equation isA conservative body force, if retained, may be incorporated into a modified pressure. Using scales $L,U_0,L/U_0$ and $\rho U_0^2$ givesHere $\nu =\mu /\rho$ is the molecular kinematic viscosity. At high Reynolds number, inertia and the pressure gradient dominate on the outer length scale, while molecular viscous diffusion is small there. Thin boundary layers, narrow shear flows and fine turbulent scales can still have sufficiently large velocity gradients for viscous effects to matter. High Reynolds number permits instability and turbulence, but does not alone prove that every flow is turbulent. Molecular viscous dissipation remains essential at the small end of the turbulent cascade.

Let $\delta \ll L$ be a representative mixing-layer thickness, ${ϵ}_m=\delta /L$, and let $U_s$ be a representative streamwise velocity. Use the scalesThe transverse velocity scale follows from the mean continuity equation. For turbulent transverse transport to balance streamwise inertia,With this ordering, the streamwise eddy-stress divergence isThus transverse diffusion of streamwise linear momentum is retained, while streamwise diffusion and the $V_x$ contribution are smaller by ${ϵ}_m^2$. Molecular kinematic viscosity is neglected where ${\nu }_T\gg \nu$.

In the transverse momentum equation, inertia and the retained eddy stresses have scale ${ϵ}_m{U}_s^2/L$, whereas an order-one transverse pressure variation would give $P_y/\rho \sim U_s^2/({ϵ}_{m}L)$. Hence the leading dimensionless transverse balance is $P_y=0$. The resulting boundary layer equations for the planar turbulent mixing layer areHere ${\nu }_T={\mu }_T/\rho$. The effective pressure convention explained in part (a)(iii) is understood. Since the outer streams have uniform horizontal velocities, their inviscid momentum equations give $P_x=0$; choosing the far-field gauge pressure as zero gives $P=0$ throughout this approximation.

For a steady two-dimensional incompressible boundary layer with tangential velocity $u(x,y)$, normal velocity $v(x,y)$, outer tangential velocity $U(x)$, and kinematic viscosity $\nu$,The pressure is uniform across the layer, and the outer Euler equation supplies $-p_x/\rho =UdU/dx$.

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.

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]

In his 1883 paper, Osborne Reynolds investigated how the character of pipe flow depends jointly on velocity, pipe diameter and viscosity. His final dye apparatus drew water from a large glass tank through a horizontal glass tube with a smooth trumpet-shaped entrance. A separate reservoir supplied strongly colored water just upstream of the entrance; a valve below the tank controlled the outflow. The glass tube was immersed in the tank, so its internal dye band could be watched. Reynolds repeated the trials with different tube diameters and water temperatures. ^[ref]

The volumetric flow rate was measured independently of the appearance of the dye. Reynolds timed a one-inch fall of the tank's water level, corresponding to a calibrated discharged volume. In modern notation, if this volume is $V$ and the elapsed time is $\Delta t$, then $Q=V/\Delta t$ and $U=4V/(\pi D^2\Delta t)$. He determined the glass tubes' internal diameters by filling them with water. For a uniform bore, a known filled length $\ell$ and measured water volume $V_{\mathrm{tube}}$ determine $D=\sqrt{4V_{\mathrm{tube}}/(\pi \ell )}$; checking bore uniformity matters because the same volume would only give an average cross-section in a tapered tube. Thermometers at different tank levels checked both the temperature and its uniformity. Reynolds compared the viscosity at different temperatures using Poiseuille's results and checked his own water by separate experiments of that type. Thus the Reynolds number estimate did not come from judging the dye band or from assuming that transition occurred at a preset value. ^[ref]

A practical repetition can instead collect a known outflow volume in a measuring vessel and time its collection, or weigh the collected water and divide its mass by mass density. Measure the internal diameter, rather than the external diameter, and obtain kinematic viscosity at the measured water temperature from independently established properties or a separate viscosity measurement. Record the dye appearance and observation position alongside $Q$, $D$, temperature and the resulting Reynolds number. These are alternative measurement procedures following $Q=V/\Delta t$, not additional historical claims about Reynolds' dye apparatus.

Begin with settled water and a small dye supply, then open the outlet valve in small steps. Reynolds saw a thin straight band at sufficiently low speeds. With a tank that had not fully settled, the band could wander while retaining its coherence. At higher discharge it could abruptly spread into the surrounding water downstream of the entrance. Short illumination by an electric spark resolved curls within what otherwise looked like a colored cloud. This connected the visible breakup with internal irregular motion rather than merely the appearance of an opaque pipe. ^[ref]

The controls are part of the result. Reynolds used a smooth entrance and let the tank rest to reduce inlet disturbances. He found that a heated tank could retain slow circulation as it cooled, so measuring one temperature was insufficient: thermal uniformity and residual motion also mattered. His wire-coil trial broke up the band below the previous onset speed, but he explicitly cautioned that once the dye was mixed he could not tell from that dye alone whether the motion became orderly again farther along the tube. In a repetition, keep the dye supply small enough to act approximately as a passive scalar, avoid a strong injection jet, and record the inlet condition and downstream observation distance. If also testing the Hagen-Poiseuille equation, allow the profile to develop and measure pressure drop between positions within the developed region; do not treat the entrance and outlet losses as part of its uniform pressure gradient. The dye experiment by itself does not establish fully developed flow. ^[ref]

The experiment's usefulness is that it compares a visible change in the motion with independently measured scales. Repeating it with several diameters and temperatures, while controlling disturbances, tests the organizing role of the Reynolds number. Its limitations are equally instructive: dye is a tracer of transport, diffusion still acts in laminar flow, and a local dye observation cannot by itself prove that turbulence persists throughout an arbitrarily long pipe.

The molecular Schmidt number is the ratio of kinematic viscosity to scalar diffusivity. A large value allows a passive scalar to develop structure below the viscous smoothing scale of the velocity field. It is distinct from the turbulent Schmidt number, which compares effective turbulent transport coefficients, and from the Schmidt number used for quantum entanglement.

For a steady zero-mass-flux background and constant viscous transport response exponent, this change removes the first derivative in the radial diffusion operator. If ${\Sigma }_0\propto r^{-p}$, the kinematic viscosity scales as $r^{p-1/2}$, so the transformed diffusivity scales as $x^{2p-3}$. It is constant at $p=3/2$.

For a homogeneous Newtonian fluid with constant kinematic viscosity, a steady axisymmetric purely azimuthal flow obeys $v_{\varphi }^{\prime \prime }+v_{\varphi }^{\prime }/r-v_{\varphi }/r^2=0$. Writing $v_{\varphi }=r\Omega$ gives $(r^3{\Omega }^{\prime }{)}^{\prime }=0$, so $\Omega =A+B/r^2$. A domain containing the rotation axis requires $B=0$ for regularity; rotating walls in an annulus can sustain $B\ne 0$. An arbitrary time-dependent profile instead evolves by ${\partial }_t\Omega =\nu ({\Omega }^{\prime \prime }+3{\Omega }^{\prime }/r)$ within this ansatz.

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.

With constant kinematic viscosity and zero inner viscous torque in an accretion disk, $\nu \Sigma \sqrt{r}=(6\pi {)}^{-1}{\int }_0^r{F}_{\mathrm{acc}}(r^{\prime })/\sqrt{r^{\prime }}d{r}^{\prime }$. For the steady photoevaporating-disk mass flux, the density is $({\dot{M}}_d-{\dot{M}}_w)/(3\pi \nu )$ inside $r_g$ and the displayed logarithmic profile outside. At equal feeding and wind rates the inner density vanishes, and the outer normalized profile rises from zero with zero slope and tends to one. Its near-edge expansion is $(r/r_g-1{)}^2/8$.

A surface Ekman layer is the stress-driven boundary correction below a rotating fluid’s upper boundary. For upward $z\le 0$, $f>0$, constant kinematic viscosity and imposed complex wind stress $T$, its steady velocity is $W_E=T{e}^{\lambda z}/(\rho \nu \lambda )$, $\lambda =(1+i)/{\delta }_E$. The underlying geostrophic flow is independent of $T$ unless an additional boundary velocity condition couples them.

For a layer rotating at angular speed $\Omega$, of depth $d$ and kinematic viscosity $\nu$, the Taylor number compares the Coriolis force to viscous damping through the squared ratio $(2\Omega d^2/\nu {)}^2$. Definitions in other geometries may include different geometric factors; this is the plane-layer convention.

When the applied inner torque dominates the advective angular-momentum flux, the constant-kinematic viscosity inner profile resembles a decretion disk. This resemblance does not alone imply outward mass transport.

The Townsend-Betchov cascade cartoon organizes inertial range vorticity into sheets by bi-axial strain, then uses sheet instability and roll-up to generate finer vortex tubes. Dissipative tube cores balance stretching with kinematic viscosity. Such differing mechanisms across scales challenge literal identical-step, memoryless cascade models without contradicting a mean energy flux.

Although the local viscous dissipation is $2\nu S_{ij}{S}_{ij}$, increasingly intense velocity gradients can compensate for vanishing kinematic viscosity. A Burgers vortex illustrates this mechanism: its radius shrinks while its excess integrated viscous dissipation remains finite. Such a local model does not by itself prove the zeroth law of turbulence for a whole flow.