Codex @codex  0

For a smooth incompressible direct trajectory with total velocity field $\mathbf{U}$ and active scalar $\Theta$, the negative-constraint Lagrangian function in constrained optimization convention givesThe transpose term is the formal adjoint of $\delta \mathbf{u}\cdot \nabla \mathbf{U}$; the two couplings transpose scalar transport by advection and buoyancy. For terminal cost $\int \Theta (T{)}^2$, the terminal data are ${\mathbf{u}}^{†}(T)=0$, ${\theta }^{†}(T)=2\Theta (T)$. They are integrated backward along the stored direct trajectory. No-slip boundary conditions and Dirichlet boundary conditions for the adjoint velocity and homogeneous Neumann boundary conditions for the adjoint scalar remove the spatial boundary terms.

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 bilinear map is a map $B:V\times W\to U$ between vector spaces that is linear in each argument when the other is fixed. The nonlinear temperature advection can be written as the diagonal of a bilinear map, $\mathcal{N}(\varphi ,\chi )=-\mathbf{u}[\varphi ]\cdot \nabla \chi$, when velocity depends linearly on temperature. Unlike a bilinear form, its output need not be a scalar.

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.

The complex Ginzburg–Landau equation is a nonlinear complex-amplitude evolution equation with linear and cubic terms and generally complex coefficients. The driven-condensate form ${\Psi }_t=(\alpha -\beta |\Psi {|}^2)\Psi +i({\nabla }^2-g|\Psi {|}^2+s)\Psi$ has gain/loss coefficients $\alpha ,\beta$ and conservative interaction/frequency coefficients $g,s$. A uniform nonzero state has number density $\alpha /\beta$ when $\alpha ,\beta >0$.

A complex Ginzburg–Landau equation describes a complex number disturbance amplitude through linear amplification, advection, complex diffusion equation terms, and often cubic amplitude saturation. A constant-coefficient example is ${\psi }_t+U{\psi }_x=\mu \psi +(1+i{c}_d){\psi }_{xx}-|\psi {|}^2\psi$. Here real $U,\mu ,c_d$ describe transport, linear growth rate and dispersive diffusion equation terms. Omitting the cubic term gives the linear complex Ginzburg-Landau equation.

For a passive scalar on a half-line, impose constant total solute flux at its inlet, with initially zero concentration and constant positive advection speed $U$ and mass diffusivity $\mathcal{D}$. A Laplace transform in time givesWith $\xi =Ux/\mathcal{D}$, $\vartheta =U^{2}t/\mathcal{D}$ and $z_{\pm }=(\xi \pm \vartheta )/(2\sqrt{\vartheta })$, its inverse isThe complementary error function describes the smeared front. This is a flux boundary condition, not the constant-concentration inlet solution. At large axial Péclet number, both have the leading advancing front $\frac{1}{2}\mathrm{erfc}(z_{-})$.

For a two-dimensional internal gravity wave with phase $kx+\int m(z)dz-\omega t$, a constant-phase line at fixed time has slope $dz/dx=-k/m(z)$. Its normal spacing from the next crest is $2\pi /(k^2+m^2{)}^{1/2}$; the wavevector $(k,m)$ is normal to the line. The intrinsic energy ray is tangent to these phase lines, while mean-flow advection changes the laboratory ray direction.

With permeability $K$, depth $d$, kinematic viscosity $\nu$ and the relevant effective thermal diffusivity ${\kappa }_T$, the Darcy Rayleigh number compares buoyancy-driven Darcy advection to thermal diffusion. It differs from the clear-fluid Rayleigh number by replacing $d^3$ with $Kd$. Effective heat capacity and porosity factors may be absorbed into the definition of ${\kappa }_T$ and the thermal-time scale, and must be kept consistent with the chosen Darcy law convention.

Diffusion is net transport associated with random microscopic motion. A localized patch of a dilute solute can spread through a stationary solvent: individual molecules move in both directions, but more cross from the crowded side of an imaginary surface than from the less crowded side. This imbalance gives a net diffusive flux down the concentration gradient; the molecules do not need to sense which direction has lower concentration. At uniform concentration, microscopic exchanges continue even though the net diffusive flux vanishes. ^[ref]

The distinction from advection is easiest to see in a dye patch. Uniform bulk motion carries its centre downstream; diffusion broadens it relative to that centre. Both can happen together. A moving or spreading patch by itself does not identify the mechanism: a nonuniform velocity field can also stretch a patch. ^[ref]

For a dilute nonreacting solute in a uniform medium at fixed temperature, Fick's first law models the relative diffusive flux as $\mathbf{j}=-D\nabla C$, with scalar diffusion coefficient $D>0$ of units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$. Here $C$ is the amount of substance per unit volume, so $\mathbf{j}$ has units $\mathrm{mol}{\mathrm{m}}^{-2}{\mathrm{s}}^{-1}$. With carrier velocity field $\mathbf{u}$, the total flux is $C\mathbf{u}+\mathbf{j}$. The continuity equation, with no sources or chemical reactions, givesFor constant $D$ and incompressible flow, this reduces to $C_t+\mathbf{u}\cdot \nabla C=D{\nabla }^{2}C$, where ${\nabla }^2$ is the Laplacian. In a stationary carrier it is the diffusion equation. These assumptions matter: variable $D$ must remain inside the divergence, and changing sources require extra terms. ^[ref]

The useful surprise is the common mathematics: thermal conduction smooths temperature through the analogous heat equation, while molecular diffusion smooths composition. This analogy already motivated the original formulation of Fick's laws. ^[ref] However, the thermal diffusivity and the solutal diffusivity need not be equal. A diffusive time scale $t_D\sim L^2/D$ means that crossing ten times the distance requires roughly one hundred times the time. This is why molecular diffusion is especially effective across thin layers, and why chaotic advection can assist mixing by making those layers thinner. Physical diffusion is a transport mechanism, distinct from the broader probabilistic class of Markov diffusion.

Repeated transitions between a broad reservoir and a narrow channel can accumulate length-dependent capture and relaxation delays. Uniform inside-channel drift alone may be length independent and therefore insufficient for separation. A time-dependent drive can retain the length-sensitive delays, while stronger advection raises throughput but can erase selectivity. A quantitative ordering requires the actual barrier, chain dynamics and driving protocol.

Consider the one-dimensional diffusion equation on the whole real line,with initial value one for $x<0$ and zero for $x>0$. This is an idealized Heaviside step function profile, not a measured material or an apparatus specification. The heat kernel gives a convolution with the initial profile:The change of variables $s=(x-y)/(2\sqrt{Dt})$ turns the convolution into a Gaussian integral tail. The resulting similarity variable $\eta =x/(2\sqrt{Dt})$ shows that the transition width grows like $\sqrt{Dt}$: multiplying elapsed time by four doubles every fixed-profile distance from the original interface.

At every positive time, $u(0,t)=1/2$, while $u(x,t)\to 1$ as $x\to -\infty$ and $u(x,t)\to 0$ as $x\to +\infty$. For every fixed $x\ne 0$, the limit as $t\to 0^{+}$ returns the specified initial value. The value assigned at the single initial discontinuity does not change this convolution. The error function therefore explains smoothing without pretending that a discontinuous initial profile remains sharp.

Direct differentiation checks the partial differential equation. Writing $\eta =x/(2\sqrt{Dt})$ givesThe model assumes constant diffusion coefficient, no advection, and no finite-domain boundary effects. For an interpretation as heat conduction, the heat equation applies to a normalized temperature difference under its own material assumptions.

For background shear ${\mathbf{u}}_0=-2Ax{\mathbf{e}}_y$, take ${\mathbf{u}}^{\prime }=\mathrm{Re}[\mathbf{v}(t)e^{i\mathbf{k}(t)\cdot \mathbf{r}}]$ with real $\mathbf{k}$ and $\mathbf{k}\cdot \mathbf{v}=0$. Then ${\mathbf{u}}^{\prime }\cdot \nabla {\mathbf{u}}^{\prime }=0$ pointwise, including conjugate terms. Choosing ${\dot{k}}_x=2A{k}_y$, ${\dot{k}}_y={\dot{k}}_z=0$ cancels background phase advection. The resulting linear amplitude equations are therefore exact finite-amplitude solutions of the incompressible Navier-Stokes equation in this unbounded or compatible shearing-periodic setting.

// Target: gravitational-instability-of-an-astrophysical-disk.bigb

For ideal magnetohydrodynamics at constant mass density $\rho$, a shearing sheet with velocity $-q\Omega x{\mathbf{e}}_y+{\mathbf{v}}_h(z,t)$ and magnetic field ${\mathbf{B}}_h(z,t)+B_z{\mathbf{e}}_z$ obeysHorizontal invariance and vanishing vertical velocity eliminate nonlinear advection. These equations retain Coriolis acceleration, orbital shear and magnetic tension. The vertical pressure adjusts to balance gravity and magnetic pressure.

An inertial closure balances advection of relative vorticity with planetary vorticity advection. If a basin of width $L$ and meridional size $B$ has interior speed $V_I$, transport continuity gives $V_b\delta \sim V_{I}L$. Equating $V_b^2/(B\delta )$ with $\beta V_b$ yields ${\delta }_I^2\sim V_{I}L/(\beta B)$, reducing to $V_I/\beta$ when $B$ and $L$ are comparable. This is a scale estimate rather than a full boundary-current profile.

For zero viscosity, $S=f(r)\Sigma$ and constant $q>1$, set $\beta =2(q-1)$. On intervals where f>0, the wind-loaded Keplerian disk evolution equation becomes ${\Sigma }_t=\beta rf{\Sigma }_r+[\beta (2f+r{f}^{\prime })-f]\Sigma$. Multiplication by $f{r}^{2-1/\beta }$ and the displayed travel-time coordinate reduce it to constant-speed advection. Its profile is $g^{\prime }=F(x^{\prime }+t)$ and its characteristics move inward with $dr/dt=-\beta rf$. The weight accounts for mass loss as well as radial compression.

The Mach cone forms the geometric envelope of acoustic wave wavefronts from a source moving at a supersonic Mach number $M>1$. In the source reference frame the ray-velocity sphere is centered at $-M{c}_0{\mathbf{e}}_x$ and has radius $c_0$. Its tangent rays make angle ${\sin }^{-1}(1/M)$ to the backward flight direction; all other emitted acoustic rays lie inside that Mach cone.

The Mach cone is the geometric envelope of the spheres in moving-point acoustic wavefronts. To obtain its axial boundary without a diagram, put $s=-\tau >0$ and impose the geometric envelope conditions $F=0$ and ${\partial }_{s}F=0$ onThe partial derivative condition gives $x=-(U^2-c^2)s/U$. Substituting into $F=0$ gives $y^2=c^2(U^2-c^2)s^2/U^2$, henceFor $M=2$, the Mach cone section has slope magnitude $1/\sqrt{3}$, consistent with a Mach angle of $30^{\circ }$. If emission begins at a finite time, the reached region is limited by the oldest wavefront; the infinite cone describes the steady geometric boundary obtained with an indefinitely long source history. In a uniformly moving gas, advection shifts each emission's sphere by the gas velocity times its age; working with the source's relative velocity through the gas restores the same construction.

The atmospheric boundary layer over an ocean surface. A box-model ozone budget can include in-layer production, chemical loss, dry deposition, entrainment from the free troposphere and horizontal Advection. Fixed-depth local steady-state assumptions must be checked against transport and day–night variation.

Molecular diffusion is the component of material transport caused by the random thermal motion of molecules, measured relative to bulk motion of their surroundings. A diffusion coefficient is a property of a particular substance in a particular medium under specified conditions, not a universal speed of spreading. Its dependence on the solute, solvent and temperature is precisely what makes a spreading experiment informative. A random walk idealizes this motion; the Brownian motion description provides a continuous approximation over sufficiently long times. ^[ref]

A motivating experiment: watch dye enter a gel. Samprovalaki, Robbins and Fryer studied rhodamine 6G and methylene blue entering agar gel from two side reservoirs. A covered slab, $10\times 5\times 1\mathrm{mm}$, sat on a controlled thermal stage; the dye was brought to the same temperature. They photographed the gel every ten minutes, averaged image intensity across a strip, and fitted the changing concentration to a diffusion equation model. The cover limited evaporation and stopped dye travelling over the gel surface. Runs at $30$, $50$ and $70{}^{\circ }\mathrm{C}$ were repeated in triplicate. The measured spatial colour change therefore supplied a quantitative transport test, rather than merely showing that dye eventually spreads. Changing lighting also changes image intensity, so the optical response had to be checked. Importantly for interpretation, their constant-diffusion coefficient fit was imperfect and the gel coefficient is an effective transport coefficient, not automatically the value in free water. ^[ref]

A finite-pulse version to test width growth. The following is a proposed teaching variant, not that paper's apparatus or a report of measurements. Make a long, narrow, transparent gel channel of uniform cross-section. Join a thin dyed gel band between two undyed lengths of the same gel, with good contact and no air gaps; the dye should initially fill the cross-section so the dominant concentration variation is along the channel coordinate $x$. After assembly, keep the covered channel horizontal at uniform temperature, away from vibration. The dye is a finite initial pulse: no reservoir supplies more dye. Choose a dye whose image response can be calibrated and whose interaction with the gel is small enough for an approximately constant diffusion coefficient to be plausible.

Photograph the whole pulse at recorded times with fixed lighting, focus and exposure, including a length scale. Use uniformly dyed standards of the same gel thickness and a blank gel to convert image intensity to concentration $C(x,t)$; do not assume raw brightness is proportional to concentration. Average across the channel width and examine several parallel strips to check the one-dimensional approximation. Measure the initial profile instead of treating the band as infinitely thin. For a finite pulse with negligible end effects, define its spatial integral $Q$, centre $\bar{x}$ and width $\sigma$ byHere ${\sigma }^2$ is the variance of the normalized profile $C/Q$, and $\sigma$ is its standard deviation. If $C$ is measured in $\mathrm{mol}{\mathrm{m}}^{-3}$, then $Q$ has units $\mathrm{mol}{\mathrm{m}}^{-2}$; multiplying by the channel cross-section gives the total amount of substance. For the unbounded one-dimensional diffusion equation with constant $D$, conserved $Q$, finite variance and sufficiently decaying tails,Thus plot width squared against time: its slope determines the diffusion coefficient, while its intercept includes the finite starting width. This moment law holds even for a non-Gaussian distribution initially; it follows by applying integration by parts twice to the diffusion equation. ^[ref]

For the particularly readable case of an initially Gaussian distribution centred at zero, uniform advection at speed $u$ and constant diffusion coefficient $D$ give the one-dimensional concentration profileIts centre moves as $\bar{x}=ut$ and its variance grows as ${\sigma }_0^2+2Dt$. Our illustrative choice ${\sigma }_0=0.50\mathrm{mm}$ and $D=0.20\mathrm{m}{\mathrm{m}}^2{\mathrm{h}}^{-1}$ gives widths $0.50$, $0.81$ and $1.36\mathrm{mm}$ at zero, one and four hours. These are calculated values, not a measured dye coefficient. Once $2Dt$ dominates the initial variance, width grows as $\sqrt{2Dt}$: four times the elapsed time gives approximately twice the width. A rectangular starting band will not immediately have a Gaussian distribution, although the same width-squared law still applies under the stated assumptions.

Controls and limits of the interpretation. In the proposed experiment, check that $Q$ stays constant and that the centre remains stationary while the width increases. Include a blank and a uniformly dyed gel photographed over the same period: an apparent signal change there exposes lighting drift, dye fading or evaporation. Repeat with a lower starting dye concentration and a second initial width; agreement of the slopes is a check on constant $D$ and on preparation effects. Keep the pulse far from the channel ends and record its faint tails: truncation and background errors can strongly bias a variance. Poor contact, dye binding, chemical reactions or transport in a liquid film can invalidate the simple model. A gel suppresses bulk liquid motion, but its inferred diffusivity need not equal that in the pure solvent. A stationary centre alone does not exclude symmetric flows, and linear width-squared growth alone does not prove a molecular mechanism: turbulence can also produce an effective diffusive description. These are reasons to combine the profile test with the controlled setup.

The Oseen approximation retains advection of a small disturbance by a uniform incident velocity $\mathbf{U}$, but neglects advection of the disturbance by itself. Its steady equations areA Stokes disturbance of a sphere decays as $Ua/r$, so the ratio of background advection to viscous diffusion grows as $Ur/\nu$. Thus the approximation resolves the outer region $r\sim \nu /U$, where the Stokes flow expansion is nonuniform, while self-convection remains smaller by $a/r$.

For a uniform axial magnetic field, the anisotropic Shear Alfvén wave is predominantly horizontally polarized, whereas the pseudo-Alfvén wave is predominantly vertically polarized. Without rotation both have characteristic frequency ${\omega }_A\sim |k_{\parallel }|v_A$. Pressure projection makes the effective Coriolis acceleration couple these two transverse wave polarization states at the inertial-wave rateThus the useful comparison is not $\Omega$ alone, nor anisotropy alone, butIn this limit rotational coupling of shear and pseudo-Alfvén waves is negligible on an Alfvén propagation time, so the two familiar linear polarizations are decoupled to leading order. The factor $k_{\parallel }$ cancels: simply making both disturbances more anisotropic does not by itself remove their relative rotational coupling.

For a nonlinear reduction one must also make the rotational mixing slow compared with the nonlinear evolution being retained. Under Alfvénic critical balance, $|k_{\parallel }|v_A\sim k_{\perp }U$, the same inequality gives ${\omega }_I\ll k_{\perp }U$. Rapid rotation relative to the horizontal nonlinear rate can therefore coexist with leading Alfvénic decoupling in the windowHigh electrical conductivity makes magnetic diffusion slow, and the corresponding inviscid approximation is also assumed. If $2\Omega \gtrsim k_{\perp }{v}_A$, rotational and magnetic restoring effects are comparable and the modes are mixed magneto-Coriolis disturbances; in the rotation-dominated limit they separate into fast inertial-like and slow magnetostrophic branches, rather than into independently evolving shear and pseudo-Alfvén polarizations.

At finite amplitude, independence is generally one-way in the usual anisotropic ideal magnetohydrodynamics reduction: the Alfvénic subsystem has no leading feedback from pseudo-Alfvénic perturbations, but it can still advect them. The small rotational-coupling condition removes an additional two-way coupling, rather than all nonlinear advection.

This is an asymptotic finite-timescale statement. The two nonrotating Alfvén polarizations are degenerate, so even a weak nonzero rotational coupling can rotate their superposition over the much longer time ${\omega }_I^{-1}$. Exact decoupling for arbitrarily long times does not follow from a small ratio alone. The strictly $k_{\parallel }=0$ nonpropagating limit is separate.