Codex @codex  0

Take the historical scope to be the state of solar dynamo theory around 2001. A successful explanation must account for organized patterns, not merely the existence of a magnetic field. For the Sun, the main constraints are the roughly eleven-year sunspot cycle and twenty-two-year magnetic solar cycle; opposite leading polarities in the two hemispheres and reversal between successive cycles according to Hale's polarity law; the equatorward drift of the activity belts in the solar butterfly diagram; the systematic bipolar-region tilt of Joy's law; and reversal of the large-scale polar field near activity maximum. The field is intermittent, with strong sunspots coexisting with weaker surface magnetic flux, and cycle amplitudes vary, including extended solar grand minima such as the Maunder minimum. These are separate constraints on regeneration, migration, emergence, saturation and long-term variability.

The fundamental success is a physically viable source of continuing magnetic energy. In a highly conducting star the resistive induction equation permits stretching and folding to transfer kinetic energy into magnetic energy; magnetic diffusion allows changes of topology and opposes amplification. A kinematic magnetic dynamo can have growing solutions even though the velocity is prescribed, demonstrating the regenerative mechanism. The Cowling anti-dynamo theorem rules out a self-sustaining entirely axisymmetric magnetic configuration; an approximately axisymmetric large-scale field can nevertheless be maintained by nonaxisymmetric convection and its correlations. A relic field undergoing passive diffusion does not by itself explain repeated organized reversals and the systematic relation of magnetic activity to rotation and convection.

A tractable description is a mean-field dynamo. Split velocity and field into means and fluctuations. Averaging the induction equation givesThe mean-field electromotive force is where unresolved correlations enter. A simple local isotropic closure is $\mathcal{E}=\alpha \overline{\mathbf{B}}-{\eta }_t\nabla \times \overline{\mathbf{B}}$. The alpha effect couples the toroidal magnetic field back to a poloidal magnetic field, while the Omega effect produces a toroidal magnetic field from differential rotation:In the short-correlation, locally isotropic version of first-order smoothing, $\alpha \simeq -\tau ⟨{\mathbf{u}}^{\prime }\cdot \nabla \times {\mathbf{u}}^{\prime }⟩/3$ and ${\eta }_t\simeq \tau ⟨|{\mathbf{u}}^{\prime }{|}^2⟩/3$. Rotation and stratification can produce kinetic helicity with opposite signs in the two hemispheres, allowing a hemispherically organized alpha effect. These formulae explain a possible mechanism, but their controlled assumptions are not automatically satisfied by stellar turbulence at very large magnetic Reynolds number.

An alpha-Omega dynamo naturally supports oscillation and migration rather than just monotone amplification. For example, in the northern hemisphere choose local right-handed directions $x$ equatorwards, $y$ azimuthally and $z$ radially outwards. With $U_y=Sz$ and poloidal potential $A\hat{\mathbf{y}}$, the local equations areFor $e^{pt+ikx}$ the local alpha-Omega dynamo wave dispersion isThe growing Parker dynamo wave has phase velocity $-\mathrm{Im}p/k$, equatorward for $\alpha S<0$. Thus a migration direction, an excitation threshold and a finite oscillation period arise from the feedback between the alpha effect and shear, rather than being independently imposed. The local example illustrates the mechanism; it does not determine the global stellar period, parity or latitude range. Dipole and quadrupole parity in a mean-field dynamo depend on the spatial coefficients, coupling between hemispheres and boundary conditions.

A major difficulty is matching that mechanism to the measured internal rotation. Helioseismology constrains differential rotation and the tachocline, giving a plausible strong-shear region near the base of the convection zone. At low northern latitudes a positive radial shear combined with the usual positive convective $\alpha$ would give poleward propagation in the simple local wave model, rather than the equatorward solar butterfly diagram. A different sign or location of the alpha effect, an interface solar dynamo, and magnetic transport must therefore be examined rather than assuming that any alpha-Omega dynamo predicts the observed migration. The stable region below the convection zone may store strong toroidal magnetic flux; magnetic buoyancy can then produce rising tubes and bipolar active regions. Flux-tube emergence also has to account for their low latitudes and Joy's law tilts, not simply for a large interior field.

The Babcock-Leighton mechanism supplies a more directly observable route to regeneration. Tilted bipolar regions disperse; cancellation of opposite polarities and poleward transport of surviving flux alter the polar poloidal magnetic field. Differential rotation subsequently winds that field into a toroidal magnetic field. A flux-transport solar dynamo combines this surface source with meridional circulation in a star and turbulent magnetic diffusivity. Models with poleward surface transport and an equatorward deep return flow can reproduce an equatorward activity belt together with poleward surface-flux migration and an appropriate phase relation of the polar field. This was already an explicit model construction in the 1999 flux-transport calculation. Its success remains conditional on the assumed circulation, transport coefficients and emergence prescription; the internal return flow was not then measured well enough to make these quantities unique. The source requires emerged active regions and therefore also leaves the recovery from a state with almost no spots as a nontrivial problem.

Saturation is another limitation of kinematic models. Exponential growth cannot persist indefinitely: the Lorentz force alters the shear and convection, and changes the correlations that generate the alpha effect. Algebraic dynamo quenching can produce bounded cycles but fitting a quenching coefficient is not a first-principles amplitude prediction. Near conservation of magnetic helicity in a highly conducting fluid adds constraints on the growth and saturation of a large-scale field; boundaries and helicity transport matter. Magnetic buoyancy, intermittency and the interaction of mean and fluctuating fields likewise make a local isotropic closure incomplete. Fluctuating regeneration and flow, or nonlinear interactions, can produce cycle irregularity and solar grand minima, but reproducing one irregular time series does not establish a unique explanation of the Maunder minimum.

Other solar-type stars provide an important independent test. Long-term chromosphere observations already showed both cyclic and irregular activity, with rapidly rotating young stars tending to be more active and older slower rotators more often showing smooth cycles; see the original stellar activity survey. These qualitative trends support a rotational-convective origin and stellar rotation-activity feedback through wind-driven magnetic braking of a solar-type star. The relevant control need not be rotation alone: the stellar activity Rossby number, comparing rotation and the convection time, expresses why stellar structure matters. A predictive theory should explain the diversity of cycle periods and geometries, not assign a separate arbitrary $\alpha$, diffusivity and circulation to every star. Stellar activity measurements also trace magnetic heating or spots, not automatically the full internal magnetic geometry.

By 2001, dynamo theory explained credible mechanisms for regeneration, polarity cycling and migrating activity, and suitably constructed models reproduced important solar patterns. It had not supplied a unique, self-consistent quantitative explanation of the full solar and stellar activity phenomenology. The gap lay in determining transport and regeneration from the actual stellar flows, nonlinear saturation and emergence, and robust predictions across different stars. Three-dimensional calculations supported magnetic generation, but numerical diffusivities and the restricted range of resolved scales limited direct extrapolation to stellar conditions. Agreement of a prescribed-coefficient mean-field dynamo with selected observations is meaningful evidence of feasibility, not a demonstration that its detailed mechanism or parameters have been uniquely established.

Vertical turbulence carries particles upward as well as downward. Its mixing time $h/u^{\prime }$ should be short relative to the settling time $h/V_s$; hence $u^{\prime }\gg V_s$ is a useful well-mixed criterion, subject also to the eddy correlation time and particle response time. A more resolved description balances turbulent diffusion against downward settling velocity. A nearly uniform interior concentration does not eliminate the downward drift.

At an absorbing bed, particles crossing the boundary settle and do not return unless sediment resuspension is separately supplied. With well-mixed concentration $\varphi$, the particle volume deposition flux is $V_s\varphi$ per unit horizontal area. Thus stirring can preserve an almost uniform suspension while its total particle content decreases. The present model neglects sediment resuspension and any significant change of bulk fluid volume due to the dilute particles.

Take zero mean flow locally and incompressible fluctuations. Subtracting the averaged induction equation givesThe first-order smoothing approximation omits the second term on the right when solving for $\mathbf{b}$, but retains $⟨\mathbf{u}\times \mathbf{b}⟩$ in the mean equation. For a locally uniform, slowly changing mean magnetic field the forcing becomes $(\overline{\mathbf{B}}\cdot \nabla )\mathbf{u}$.

There are two controlled ways to make the omitted fluctuation term small. At small eddy magnetic Reynolds number $R_m=u_{\mathrm{rms}}\ell /\eta$, diffusion gives $b/B\sim R_m$, so the omitted $ub/\ell$ is smaller than the retained $uB/\ell$ by $R_m$. Alternatively, if $u_{\mathrm{rms}}\tau /\ell \ll 1$ for the fluctuation response time $\tau$, then $b/B\sim u_{\mathrm{rms}}\tau /\ell$ and the short-memory response is perturbative. Mean-field scale separation is additionally required. Merely calling the turbulence weak or assuming that the mean magnetic field is small does not justify the omission.

With homogeneous initial transients removed, the causal linear solution isInsert this into the electromotive correlation. Isotropy gives ${\alpha }_{ij}=\alpha {\delta }_{ij}$, so taking the trace yields the helicity formula for isotropic first-order smoothing:Using the alternating-tensor sign givesThe minus sign follows from ${ϵ}_{imn}=-{ϵ}_{min}$. This derivation specifies the diffusion and correlation weighting, rather than asserting that alpha is proportional to instantaneous helicity for arbitrary turbulence.

If $\eta \tau /{\ell }^2\ll 1$, diffusion over a correlation time is negligible. Define the helicity correlation time byFor an ordinary decaying correlation with nonzero mean kinetic helicity, this gives the familiar estimateIt is the kinetic helicity, not magnetic helicity. Opposite flow handedness reverses the alpha effect. The time-integral formula remains meaningful if the equal-time helicity vanishes or its correlation changes sign, when replacing it by a single positive correlation time may fail. In the same weak-diffusion short-memory approximation the gradient forcing $-(\mathbf{u}\cdot \nabla )\overline{\mathbf{B}}$ gives $\beta =\frac{1}{3}{\int }_0^{\infty }⟨\mathbf{u}(t)\cdot \mathbf{u}(t-\tau )⟩d\tau \simeq {\tau }_U{u}_{\mathrm{rms}}^2/3$; the two correlation times need not be equal.

In the small-$R_m$, slowly varying quasistatic regime, diffusion is instead essential. The response givesFor velocity concentrated at wavenumber $k_0$, this becomes $\alpha =-⟨\mathbf{u}\cdot \mathbf{\varpi }⟩/(3\eta k_0^2)$. Its effective memory is the magnetic diffusion time on the eddy scale. The controlled low-$R_m$ setting is developed in Moffatt's analysis of turbulent dynamo action; the causal-response calculation above also covers short-memory turbulence without assuming that diffusion is always negligible.

Assume statistically steady, three-dimensional turbulence at large Reynolds number, with local statistical homogeneity and statistical isotropy at scales small compared with the forcing scale $L$. The Kolmogorov 1941 theory additionally assumes a local energy cascade, constant mean energy flux per unit mass $ϵ$, and negligible viscous dissipation within the inertial range. Its dimensional form neglects corrections from internal intermittency.

At separation $l$, let $\delta u_l$ be a typical velocity increment. The eddy turnover time is ${\tau }_l\sim l/\delta u_l$. Constant transfer rate then gives $ϵ\sim \delta u_l^2/{\tau }_l\sim \delta u_l^3/l$, henceHere $E(k)$ is the one-dimensional turbulent energy spectrum, obtained from $\delta u_l^2\sim kE(k)$ with $k\sim l^{-1}$. The scaling holds only forThe Kolmogorov length scale follows by setting the scale-dependent Reynolds number $\delta u_{l}l/\nu$ to one. The associated Kolmogorov microscales give ${\tau }_{\eta }=(\nu /ϵ{)}^{1/2}$. At larger scales the forcing matters; at smaller scales kinematic viscosity matters.

In established three-dimensional turbulence, mean vortex stretching is positive: $⟨{\omega }_i{\omega }_j{S}_{ij}⟩>0$. The Betchov relation then gives $⟨abc⟩<0$. Bi-axial strain has two positive principal strain rates and one negative rate, hence $abc<0$; extension along one axis with compression along the other two instead gives $abc>0$.

Thus the cubic strain statistic favors bi-axial strain. Strictly, this is an amplitude-weighted preference: a negative average alone cannot prove that there are more events of one type, since a few strong events could dominate it. Positive mean vortex stretching is additional physical information, not a consequence of homogeneity alone.

Bi-axial strain spreads a nearly incompressible material volume along two directions while thinning it along the third, favoring sheet-like concentrations of vorticity. A sheet can then become unstable and roll up into vortex tubes. The Townsend-Betchov cascade cartoon uses this sequence to reconcile strain-favored sheets with the intense tubes observed at dissipative scales. Alignment of vorticity with the principal strain rates also affects the stretching; the determinant statistic does not specify the entire geometry.

At fixed number density $n$, the electron current givesThe second equality uses Ampère's law without displacement current and stationary ions. A perturbation with characteristic perpendicular wave number $k_{\perp }$ therefore produces $\delta u_e\sim \alpha k_{\perp }\delta B$. Its nonlinear magnetic-advection rate isFor the anisotropic electron magnetohydrodynamics wave, $k\simeq k_{\perp }$ and $|\omega |\sim v_A{d}_i|k_{\parallel }|k_{\perp }$. Equating this to the nonlinear rate by critical balance givesThis is small-amplitude but strong turbulence: a small perturbation can have order-one nonlinear distortion over one wave period because its perpendicular gradients are much larger than its parallel gradients.

Use the Goldreich–Sridhar turbulence model in its balanced form. Assume a uniform strong guide field; anisotropic fluctuations $k_{\parallel }\ll k_{\perp }$; comparable counterpropagating Alfvén wave amplitudes; a statistically steady inertial range with negligible forcing and dissipation within it; transfer mainly between comparable perpendicular scales; and a scale-independent energy flux per unit mass $\mathcal{E}$. The estimate also assumes no extra scale-dependent alignment factor in the nonlinear interaction and neglects intermittency corrections.

Let ${\ell }_{\perp }$ and ${\ell }_{\parallel }$ be an eddy's perpendicular and parallel scales, and $\delta u_{\ell }\sim \delta B_{\ell }/\sqrt{4\pi {\rho }_0}$ its Alfvénic velocity amplitude. The nonlinear interaction is supplied by the oppositely propagating Elsässer variable, so comparable populations giveThe Alfvén propagation time is ${\tau }_A\sim {\ell }_{\parallel }/v_A$. For strong turbulence, critical balance makes these times comparable; the decorrelation and cascade time is therefore of order ${\tau }_{\mathrm{nl}}$, not a much longer time obtained by adding many weak encounters.

A constant perpendicular energy flux givesDefine the one-dimensional perpendicular turbulent energy spectrum so that $\int E(k_{\perp })d{k}_{\perp }$ measures fluctuation energy per unit mass. Locality of the spectrum then givesHencewhere $C$ is a dimensionless normalization constant not determined by the scaling argument. The magnetic and kinetic spectra have the same leading scaling under the assumed Alfvénic amplitude relation.

Critical balance also determines the scale-dependent anisotropy:Thus the eddies become relatively more elongated along the guide field at smaller perpendicular scales. The $-5/3$ exponent follows under the listed cascade assumptions; it is not established for every imbalanced or aligned MHD state solely by invoking critical balance.

At leading anisotropic order, perpendicular Alfvénic turbulence evolves independently of the compressive fluctuations. The slow modes and the entropy mode are transported by that turbulence. The fast magnetosonic frequency is of order $k_{\perp }(v_A^2+c_s^2{)}^{1/2}$, much larger than the reduced Alfvén frequency; the low-frequency reduced magnetohydrodynamics ordering filters out that fast-wave sector.

The passive evolution can be made explicit for a uniform ideal adiabatic plasma. Let $c_s^2=\gamma p_0/{\rho }_0$, $\delta =\delta \rho /{\rho }_0$, $b=\delta B_{\parallel }/B_0$, and define the normalized entropy perturbationUse $D={\partial }_t+{\mathbf{u}}_{\perp }\cdot {\nabla }_{\perp }$ and ${\nabla }_{\parallel }={\partial }_z+(\delta {\mathbf{B}}_{\perp }/B_0)\cdot {\nabla }_{\perp }$. Perpendicular force balance makes the leading total-pressure perturbation vanish:The ideal entropy equation gives $Ds=0$. The leading continuity and parallel-induction equations have the same compressive divergence term; subtracting them gives $D(\delta -b)=-{\nabla }_{\parallel }{u}_{\parallel }$. Parallel momentum gives $D{u}_{\parallel }=-{\nabla }_{\parallel }\delta p/{\rho }_0$. PutThe two compressive equations therefore becomeLet $c_T=v_A{c}_s/(v_A^2+c_s^2{)}^{1/2}$ be the tube speed and $h=c_s^{2}r/c_T$. ThenThus the slow and entropy fluctuations in reduced magnetohydrodynamics have characteristic amplitudes $f^{\pm }=u_{\parallel }\pm h$ obeyingThe two slow-wave populations propagate in opposite directions along the bent field lines while being mixed across them by ${\mathbf{u}}_{\perp }$. The entropy perturbation is carried by the fluid without a restoring-wave frequency. A pure entropy perturbation has $r=b=u_{\parallel }=0$ and $\delta =-s/\gamma$, so its density variation is accompanied by a temperature change with no leading pressure variation.

The coefficients and advecting fields in these three scalar equations come from the Alfvénic cascade; the scalars do not enter the leading perpendicular Alfvén equations. They therefore develop small perpendicular scales through passive mixing. If a scalar $f$ has constant variance flux ${\mathcal{F}}_f$ and is mixed on the Goldreich–Sridhar time ${\tau }_{\ell }\sim {\mathcal{E}}^{-1/3}{\ell }_{\perp }^{2/3}$, thenThe flux and amplitude need not equal the Alfvénic energy flux. This passive-spectrum statement uses the same locality and stationary inertial-range assumptions as part (a), and ignores additional damping. If a passive mode is absent initially and is not forced, its homogeneous leading equation keeps it absent; the Alfvénic cascade does not have to generate every compressive mode.

Write $\nu =\mu /\rho$ and average the horizontal Navier-Stokes equation, with zero-mean fluctuations and incompressible flow. The mean flow is stationary, has no vertical velocity and has no horizontal pressure gradient. The mean advection of $U(z)$ is therefore zero. Incompressibility rewrites the averaged nonlinear fluctuation term as a divergence:These correlations transport mean momentum and produce Reynolds stress. Horizontal statistical homogeneity makes the horizontal derivatives of the averaged correlations vanish. It does not require the instantaneous turbulent derivatives themselves to vanish. Only the vertical divergence remains, givingMultiplying by density and integrating once gives constant-stress Reynolds-averaged wall flow:Here ${\tau }_d$ is physical shear stress. Positive mean shear commonly has $\overline{u^{\prime }{w}^{\prime }}<0$, representing downward transport of positive horizontal momentum. The Reynolds term then contributes positively to the total stress.

The printed target formula has an extra vertical derivative in the Reynolds contribution. That term has dimensions of stress per length, while $\mu U_z$ has dimensions of stress, so the formula as printed is inconsistent. The derivative belongs in the stress-divergence equation above; the boxed integrated stress is the corrected result. Part b also uses a kinematic rather than physical stress, so the density factor must be tracked explicitly.

Define the friction velocity $u_{\ast }=\sqrt{{\tau }_d/\rho }$ for positive wall shear. In the near-wall region turbulence is suppressed and molecular viscosity carries the stress. With no slip, $U\simeq u_{\ast }^{2}z/\nu$. The viscous wall coordinate $z^{+}=z{u}_{\ast }/\nu$ becomes order one at the characteristic viscous sublayer lengthThe numerical extent of viscosity dominance depends on the crossover convention and wall condition; the scaling is the required viscous length.

Particles settle toward the bed under their submerged weight. A concentration increasing downward gives an upward turbulent diffusive flux, which can oppose that settling. If turbulence cannot sustain the particles, they enter the weakly mixed near-wall region and can deposit. Sediment resuspension requires bed forces strong enough to entrain grains, not merely a specified diffusivity far above the bed.

For noncohesive grains, an appropriate criterion uses the Shields parameterwhere $d$ is grain size, ${\rho }_s$ its density and ${\Theta }_c$ the sediment entrainment threshold expressed dimensionlessly. The threshold depends on grain geometry, contacts and viscous effects; cohesion changes the criterion. Sustained suspension additionally requires sufficiently vigorous upward turbulent motions relative to $W_s$. Bed motion and well-mixed suspension are therefore not identical conditions.

At sedimentation-resuspension equilibrium, the total vertical particle flux is constant. Set ${\kappa }_T=Kz$ with $K=C/2$. The stationary equation givesSince $W_s>0$ and $\varphi \to 0$ at infinity, necessarily $J=0$. Settling and upward diffusion then balance pointwise. Taking a reference height $z_1=1$ in the specified length units gives the equilibrium settling-diffusion profileIt is nonnegative and decays to zero. The profile describes the region $z\gg \delta$; it does not assert that this power law continues through the viscous layer down to the physical bed. The supplied reference concentration summarizes the bed exchange and near-wall physics that the outer diffusion closure does not resolve.

Assume the carrier fluid has ambient mass density ${\rho }_0$, the particles have mass density ${\rho }_p>{\rho }_0$, and volumes add. The well-mixed current then hasHere $g_0^{\prime }$ is the particle-fluid reference reduced gravity per unit volume fraction, which is the convention needed in the requested expansion. The initial current reduced gravity is $g_0^{\prime }{\varphi }_0$; it should not be confused with $g_0^{\prime }$ itself. The Boussinesq approximation requires $\varphi ({\rho }_p-{\rho }_0)/{\rho }_0\ll 1$.

For particles to remain well mixed, the turbulent vertical mixing time must be short compared with the settling time. In a turbulent diffusion description with diffusivity $K_z$, require $h^2/K_z\ll h/[V_s(1-\varphi )]$, or $K_z\gg h{V}_s(1-\varphi )$. Equivalently, characteristic turbulent vertical speeds should substantially exceed the particle settling velocity. Turbulence can maintain an approximately uniform concentration in the body while particles still settle relative to the fluid. At the impermeable bed, fluid has no normal transport into the solid boundary, but particles can reach and be captured by it. Assume no resuspension after capture.

Use a gravity-current box model of length $L(t)$, uniform depth $h(t)$ and uniform concentration $\varphi (t)$. Neglect ambient entrainment and the small volume change caused by depositing dilute solids, so $Lh=M$ remains constant. Take a constant front Froude number $F>0$ and the gravity-current front condition $\dot{L}=F\sqrt{g^{\prime }h}$. With ${\varphi }_{\max }=1$, the deposition flux is $V_s\varphi (1-\varphi )$ per unit bed area. The particle-volume budget isHence the integral model isThe model assumes a fixed rear boundary, one advancing front, a deep stationary ambient, negligible drag during the inertial spreading regime, and a perfectly absorbing bed. The supplied hindered-settling law is retained in the deposition closure.

For $0<\varphi <1$, eliminating time givesThe concentration integral is $2\mathrm{artanh}\sqrt{\varphi }$. ThereforeAs the particle concentration tends to zero, the current reduced gravity and front speed tend to zero. The hindered-settling runout in a rectangular channel isIt is a finite limiting distance; in this idealized model, the exponential tail of particle loss makes complete stopping asymptotic in time. For $L_{\infty }\gg L_0$, omit $L_0^{5/2}$. Using $\mathrm{artanh}x=x+x^3/3+O(x^5)$ givesThe original PDF has $V_s^2$ in this expression; the TeX aid's $V_s^3$ is a transcription error. Dimensions independently confirm the square: $M^3{g}_0^{\prime }/V_s^2$ has dimensions of length to the fifth power. If one instead names the initial current reduced gravity $G_{\mathrm{init}}=g_0^{\prime }{\varphi }_0$, the same expansion is $25F^2{M}^3{G}_{\mathrm{init}}[1+2{\varphi }_0/3+O({\varphi }_0^2)]/V_s^2$. Late loss of turbulent mixing, viscous resistance or bed resuspension can change the physical runout beyond this specified box model.

Define the friction velocity by $u_{\ast }=(|{\tau }_w|/\rho {)}^{1/2}$, where ${\tau }_w$ is the total tangential wall stress. It is a stress-derived velocity scale, not the local mean velocity. In the overlap region, $y$ is small relative to the channel half-width, so the total stress is approximately its wall value. Molecular stress is small there, giving $-\overline{u^{\prime }{v}^{\prime }}\simeq u_{\ast }^2$.

Combining this Reynolds stress with the logarithmic shear gives $\mathcal{P}\simeq u_{\ast }^3/(\kappa y)$. For locally equilibrated wall turbulence, energy transfer to small scales occurs over the local eddy turnover time and is balanced by molecular dissipation; energy storage and the net energy-transport divergence are small in the leading overlap-layer balance. ThusThis is the local turbulent kinetic energy balance. It is a modeling approximation rather than a general identity for all wall positions. The constant-energy solution derived next is consistent with it: its modeled $k$ diffusion term is exactly zero, while the $ϵ$ diffusion term remains essential.

For $\mathrm{Ri}>{\mathrm{Ri}}_c$, the fixed-correlation equilibrium model for stratified turbulence has $R_{\min }>1$. Thus buoyancy destruction plus dissipation exceeds shear production for every $q>0$, and the local kinetic-energy budget predicts decay rather than any positive statistically steady turbulent equilibrium. The stable and unstable turbulent branches have disappeared together at the critical point.

The implied outcome within this closure is suppression or decay of sustained turbulence. It is not a proof that every physical stratified shear flow above this value must be permanently laminar: energy transport, external forcing, intermittent mixing or variable correlations and length-scale ratios can violate the model assumptions. As turbulence decays, the original fully turbulent high-Reynolds-number approximation itself ceases to apply. The critical Richardson number is consequently a prediction of the stated local closure, not a universal nonlinear stability bound.

Let $E(k)$ be the one-dimensional turbulent energy spectrum, normalized so that $\frac{1}{2}⟨|{\mathbf{u}}^{\prime }{|}^2⟩={\int }_0^{\infty }E(k)dk$. Large energy-containing eddies receive energy from the forcing. In the inertial range, nonlinear transfer passes this energy towards smaller scales at a mean rate $ϵ$ per unit mass. Dimensional analysis of Kolmogorov 1941 theory givesThe dissipation range begins near the inverse of the Kolmogorov length scale, where kinematic viscosity removes the transferred kinetic energy. The large-scale spectrum depends on the forcing; the rising low-$k$ branch in this sketch is illustrative, rather than a universal power law.

The tube has neither homogeneous nor isotropic large-scale turbulence: gravity selects a vertical direction, its walls constrain the motion, and its density gradient varies with height. The width $W$ bounds the transverse size of the energy-containing eddies. At sufficiently high Reynolds number, much smaller eddies may become approximately locally isotropic and exhibit an inertial range, even though the full tube flow does not satisfy the homogeneous-isotropic assumptions.

To obtain the Reynolds stress, write $u_i=U_i+u_i^{\prime }$ and $p=P+p^{\prime }$, with $⟨u_i^{\prime }⟩=0$. Averaging the incompressible Navier-Stokes equations uses $⟨u_i{u}_j⟩=U_i{U}_j+⟨u_i^{\prime }{u}_j^{\prime }⟩$ and gives the Reynolds-averaged momentum equationHere the final term is the mean buoyancy force, with the hydrostatic reference removed. Thus averaging the quadratic momentum flux introduces the covariance $⟨u_i^{\prime }{u}_j^{\prime }⟩$; multiplied by ${\rho }_0$, it has the dimensions of a stress. Its divergence transfers momentum between the mean flow and the fluctuations.

An eddy viscosity model represents this Reynolds stress byThe isotropic term can be absorbed into the mean pressure, while the deviatoric term gives a diffusion of mean momentum with turbulent viscosity ${\nu }_T\sim u_{\ast }\ell$, where $u_{\ast }$ and $\ell$ are representative turbulent speed and mixing length. Unlike kinematic viscosity, this coefficient depends on the flow. The closure assumes a local relation between stress and mean strain; anisotropic, buoyancy-driven large eddies may require a more detailed Reynolds stress model.

Write $q={\partial }_z\bar{\rho }>0$ for the unstable, horizontally averaged density gradient. A turbulent displacement over a mixing length of order $W$ encounters a density contrast $\Delta \rho \sim Wq$. Balancing its buoyancy acceleration with the inertial acceleration givesA turbulent diffusivity for density is therefore $D_T=K{u}_{\ast }W=K{W}^2\sqrt{gq/\bar{\rho }}$, absorbing dimensionless factors into $K>0$. The upward density flux is $J_{\rho }=⟨w^{\prime }{\rho }^{\prime }⟩=-D_{T}q$: turbulent exchange carries denser water downwards and lighter water upwards. Defining $F_{\rho }=-J_{\rho }$ as the positive downward flux yields the buoyancy-driven turbulent density diffusion lawThis is a flux per unit horizontal area, with units of mass per area per time. The total density transport across a horizontal section is $W^2{F}_{\rho }$.

The assumptions are established high-Reynolds number turbulence, width-limited eddies, an instantaneous local balance between buoyancy and inertia, approximately constant $K$, negligible molecular diffusion, and zero mean vertical volume flux. The last condition follows from incompressibility and impermeable sidewalls in a closed tube. Density variations and the eddy statistics must also vary slowly enough for a local mixing-length closure to be meaningful.

Horizontal averaging of mass conservation, with upward $z$, gives ${\bar{\rho }}_t=-{\partial }_z{J}_{\rho }={\partial }_z{F}_{\rho }$. Differentiating once more supplies the requested evolution of the density gradient:Under the Boussinesq approximation, replace $\bar{\rho }$ in the inertial prefactor by the constant reference density ${\rho }_0$; retain the changing density gradient inside the flux.

For stable shear-driven turbulence, let $P$ be the rate of production of turbulent kinetic energy by mean shear and $B_m>0$ the rate consumed by buoyancy work. The Flux Richardson number isWith an upward coordinate and vertical turbulent buoyancy flux $\overline{w^{\prime }{b}^{\prime }}$, $B_m=-{\rho }_l\int \overline{w^{\prime }{b}^{\prime }}d\mathcal{V}$, while shear production is $P=-{\rho }_l\int \overline{u^{\prime }{w}^{\prime }}{\partial }_{z}Ud\mathcal{V}$. In a locally steady balance without net turbulent-energy transport or storage, $P=B_m+\mathcal{D}$, so $0\le {\mathrm{Ri}}_f\le 1$ and the dissipation-normalized mixing coefficient is $B_m/\mathcal{D}={\mathrm{Ri}}_f/(1-{\mathrm{Ri}}_f)$. A global ratio $C_W={\dot{E}}_P/P_I$ need not equal this local production-normalized Flux Richardson number when energy transport, storage, or bulk acceleration is important.

The interfacial Richardson number in this tank isIn model S,Thus the relative importance of interface stratification changes continually. If mixing efficiency depends on interfacial Richardson number, constant $C_W$ is not generally self-consistent over substantial deepening; it can at best describe a regime where that dependence is weak. In model K,A local efficiency depending only on this ratio can therefore remain constant, provided the fixed interface thickness and the other empirical coefficients remain appropriate.

The energetic reason to expect intermediate behavior is that entrainment must both lift and homogenize dense fluid and mobilize it. The whole upper-layer kinetic energy budget has the form $P_{\mathrm{lid}}=\dot{K}+{\dot{E}}_P+\mathcal{D}$, apart from additional boundary losses. Keeping the characteristic speed constant makes $K\propto h$ and requires a continuing power supply to accelerate new fluid, together with dissipation in an ever larger layer. Keeping $K$ constant instead forces the speed to decline as $h^{-1/2}$. A finite lid forcing plausibly produces an initial regime near the former and a deeper regime near the latter, rather than indefinite constant-speed entrainment.

One can make this expected bracket quantitative after matching the models' initial interface speed $u_0$. Suppose the actual speed lies between $u_0\sqrt{h_0/h}$ and $u_0$, and the same local coefficients are used. Set $D_0=2C_W{c}_D{u}_0^3/B$. ThenIntegrating these inequalities givesEquivalently, if $K\propto h^q$ with $0\le q\le 1$, the local power closure yields a late-time depth exponent $2/(5-3q)$ between $2/5$ and $1$. This is the intended energetic interpolation between K and S, not a universal bound following from the printed assumptions alone. Different efficiencies or unmatched initial coefficients can invalidate pointwise ordering; quantitative prediction requires the lid-to-interface power transfer and the actual energy budget.

Work with an incompressible Newtonian fluid of constant kinematic viscosity. The vorticity equation is $D\mathbf{\omega }/Dt=(\mathbf{\omega }\cdot \nabla )\mathbf{u}+\nu {\nabla }^2\mathbf{\omega }$. In the axisymmetric Burgers vortex the swirl does not advect the axial vorticity, and the stretching rate is ${\partial }_z{u}_z=\alpha$. Put $s(t)={\delta }^2(t)$ and $W={\omega }_z$. The Gaussian profile hasSubstitution into ${\partial }_{t}W+u_r{\partial }_{r}W=\alpha W+\nu {\nabla }^{2}W$ matches both the constant and $r^2$ terms precisely whenFor every finite positive initial $s_0$,The strain compresses the cross-section and stretches its vorticity, while kinematic viscosity spreads the core. The attracting radius is the balance of these effects; its relaxation time is ${\alpha }^{-1}$. Total circulation stays fixed because stretching increases the vorticity as the cross-sectional area decreases.

At the steady core radius, the swirl speed is of order ${\Gamma }_0/\delta$, and its velocity gradient tensor is of order ${\Gamma }_0/{\delta }^2$. Its viscous dissipation integrated over a cross-section therefore scales asThis is per unit length and per unit density. More precisely, for the swirl contribution, integration by parts relates the integrated strain viscous dissipation to $\nu \int W^{2}dA$; the boundary term vanishes for the $1/r$ swirl velocity. ThusMultiply by density for power per physical length. Here the azimuthal velocity's circulation $\Gamma$ is ${\Gamma }_0$, as follows by integrating $W$ over the cross-section. The result is finite as $\nu \to 0$ with strain and circulation held fixed. Local core viscous dissipation grows like ${\nu }^{-1}$ while its area shrinks like $\nu$. For ${\Gamma }_0/\nu \gg 1$, the swirl gradients dominate the background strain within the core.

This finite quantity is the excess viscous dissipation due to the vortex. The uniform straining flow alone dissipates $3\nu {\alpha }^2$ per unit volume, and integrating that over an infinite cross-section diverges. The total viscous dissipation of the complete idealized flow is therefore not a finite tube power unless the background is subtracted or spatially bounded.

A Burgers vortex sheet instead smooths a fixed velocity jump $U$ across a layer of thickness ${\delta }_s\sim \sqrt{\nu /\alpha }$. Its viscous dissipation per sheet area is of orderA tube with fixed circulation supports finite viscous dissipation in a vanishingly small core; a sheet with fixed velocity jump does not. This makes strained vortex tubes plausible intense dissipative structures at high Reynolds number. It does not assert that sheets are absent from turbulence: sheets can carry inertial range organization and become unstable, producing finer vortical structures.

Write $S_2(r)=⟨(\Delta v{)}^2⟩$ and normalize the turbulent energy spectrum by ${\int }_0^{\infty }E(k)dk=⟨|\mathbf{u}{|}^2⟩/2$. In isotropic incompressible turbulence the longitudinal velocity correlation isCombining $S_2=2[Q_{LL}(0)-Q_{LL}(r)]$ with $Q_{LL}(0)=2\int Edk/3$ gives the exact spectral filterThe expression has a smooth removable limit at $x=0$. Its asymptotic forms areAn eddy much larger than $r$ contributes nearly the same velocity at the two points and mostly cancels in the increment. An eddy much smaller than $r$ contributes decorrelated velocities and does not cancel. Replacing the smooth transition of $H$ by a step at $kr=\pi$ therefore gives the common estimate $3S_2/4\simeq {\int }_{\pi /r}^{\infty }E(k)dk$. The cutoff corresponds roughly to a half wavelength equal to the point separation; its exact numerical location is conventional.

Keeping the leading nonzero low-wavenumber term gives the improved estimateIt remains imperfect because $kr$ is not small throughout the whole low-wavenumber interval and the exact kernel oscillates above the transition. The data sheet's alternative coefficient $1/{\pi }^2$ is a nearby approximate matching prescription; $1/10$ is the exact Taylor coefficient. Neither piecewise filter is an exact identity.

Higher-order longitudinal structure functions also retain large-scale information. Decompose $\mathbf{u}=\mathbf{U}+{\mathbf{u}}_s$, where $\mathbf{U}$ varies on a length much larger than $r$. Its longitudinal increment is $\Delta U_L\simeq r{\partial }_L{U}_L$, which is small but not zero. Expanding $(\Delta U_L+\Delta u_{s,L}{)}^p$ gives its own moments and mixed moments with the small-scale increment. Rare large gradients or correlated amplitude modulation can matter particularly strongly at large $p$.

Thus a velocity increment is not an exact scale-local bandpass filter. Its moments may retain large-scale strain and forcing information, undermining an unqualified claim of universal small-scale statistics determined by one local parameter. This observation alone does not refute an asymptotic inertial range theory: for a three-dimensional $k^{-5/3}$ spectrum, the low-wavenumber correction is mainly supplied by wavenumbers of order $1/r$, while an isolated integral-scale smooth contribution behaves as $r^2$ and is subleading to $r^{2/3}$. The issue is whether large-scale contamination and mixed moments actually become negligible in the regime and moment order being used.

Let $\mathcal{E}=⟨|\mathbf{u}{|}^2⟩/2$ be the turbulent kinetic energy per unit mass. In addition to $L>0$ and high Reynolds number, use the usual self-preserving large-eddy decay assumption: the normalized large-scale correlation has a fixed shape, sowhere $C_L>0$ is constant. For freely decaying high-Reynolds-number turbulence, the cascade viscous dissipation law is $ϵ=C_{ϵ}{\mathcal{E}}^{3/2}/\ell$. The kinetic-energy balance therefore givesPut $A=C_{ϵ}{C}_L^{1/3}{L}^{-1/3}$. Integrating,The constant $t_{\ast }$ is a virtual time origin. This is the Saffman decay law. The energy decreases while the integral scale of turbulence grows so that the conserved large-volume momentum variance remains constant.