Codex @codex  0

For a flow satisfying incompressibility, this antisymmetric tensor writes advection and vortex stretching in conservative form. With its first index as transport direction, $(\nabla \cdot T{)}_j={\partial }_i{T}_{ij}$.

For differentiable homogeneous incompressible flow, $⟨\mathrm{tr}{A}^3⟩=0$. One divergence representation is $\mathrm{tr}{A}^3={\partial }_i[u_j({\partial }_j{u}_k)({\partial }_k{u}_i)-u_i\mathrm{tr}{A}^2/2]$. Splitting the velocity gradient tensor into strain and rotation gives $⟨{\omega }_i{\omega }_j{S}_{ij}⟩=-4⟨abc⟩$. It connects mean vortex stretching to the cubic strain statistic, without requiring isotropy.

Enstrophy measures squared vorticity: a local density is ${\omega }^2/2$, and its spatial integral or ensemble mean gives a global measure. In homogeneous incompressible flow, viscous dissipation satisfies $ϵ=\nu ⟨{\omega }^2⟩$. Vortex stretching can produce enstrophy in three dimensions; kinematic viscosity destroys it.

For isotropic incompressible flow, directional averaging gives $⟨({\partial }_x{u}_x{)}^3⟩=(8/105)⟨\mathrm{tr}{S}^3⟩=-(2/35)⟨{\omega }_i{\omega }_j{S}_{ij}⟩$. Positive mean vortex stretching therefore implies negative longitudinal derivative skewness. Statistical homogeneity alone does not determine the sign in a fixed direction.

Write $\mathbf{\omega }=\nabla \times \mathbf{u}$ and use the vector identity $(\mathbf{u}\cdot \nabla )\mathbf{u}=\nabla (|\mathbf{u}{|}^2/2)-\mathbf{u}\times \mathbf{\omega }$. Taking the curl of the momentum equation eliminates all pure gradients, including gravity and the $\nabla (\nabla \cdot \mathbf{u})$ viscous contribution. Before invoking barotropy, the pressure term givesThis is baroclinic vorticity generation. For a barotropic fluid, $p=p(\rho )$, so $\nabla p$ is parallel to $\nabla \rho$ and this source vanishes. Constant $\nu$ permits the curl and Laplacian to commute. The result isNo incompressibility assumption was used. Since $\nabla \cdot \mathbf{\omega }=0$, an equivalent form isThe terms represent vortex stretching/tilting, compression/dilatation and viscous diffusion. In particular, removing the baroclinic source does not make the vorticity of a three-dimensional parcel constant. A spatially varying entropy-dependent $K$ in $p=K{\rho }^{\gamma }$ need not define a single barotropic law; the common-pressure-law assumption matters.

The quasi-geostrophic approximation isolates slowly evolving, nearly balanced motions in a rotating, stably stratified fluid. Ocean eddies and mid-latitude atmospheric disturbances often have this character: a leading Coriolis force–pressure balance suppresses rapid horizontal acceleration, while hydrostatic balance suppresses rapid vertical acceleration. Small ageostrophic flow nevertheless supplies the vertical motion that changes the leading vorticity. The approximation filters fast inertia-gravity waves while retaining nonlinear transport, vertical shear and Rossby waves.

Use a reference density ${\rho }_0$ and remove the background hydrostatic pressure. Let $\pi =p^{\prime }/{\rho }_0$ and let the upward buoyancy perturbation be $b=-g{\rho }^{\prime }/{\rho }_0$. A vertically varying reference density stratification has $N^2(z)>0$. Under the Boussinesq approximation, the inviscid and adiabatic equations areHere $D/Dt={\partial }_t+{\mathbf{u}}_H\cdot {\nabla }_H+w{\partial }_z$. Density variations are small compared with ${\rho }_0$ but retained in the vertical buoyancy force. This is a useful local, incompressible description; a compressible atmosphere needs the corresponding weighted or pressure-coordinate version.

Let horizontal and vertical lengths be $L,H$, speed be $U$, and take the slow time $L/U$. Take $f_0>0$ for definiteness; use $|f_0|$ in the scale estimates in the opposite hemisphere. The principal small parameter is the Rossby number $\mathrm{Ro}=U/(f_{0}L)\ll 1$ at a latitude where $f_0\ne 0$. For the hydrostatic, strongly stratified regime, $H/L\ll 1$ and $f_0/N\ll 1$. Retaining both horizontal vorticity and vortex stretching gives the Burger numberThe geostrophic balance pressure scale is $f_{0}UL$ and the buoyancy perturbation scale is $B=f_{0}UL/H$. Hence $B/(N^{2}H)=\mathrm{Ro}/\mathrm{Bu}\ll 1$: buoyancy perturbations produce small displacements relative to the background stratification scale. The ageostrophic flow is $O(\mathrm{Ro}U)$ and incompressibility then gives $w=O(\mathrm{Ro}UH/L)$. Vertical acceleration is smaller than hydrostatic terms by order ${\mathrm{Ro}}^2(f_0/N{)}^2$ in this scaling. Weak forcing or diffusion can be added at the slow order, but neither is present here.

Define the quasi-geostrophic streamfunction by $\psi =\pi /f_0$. Leading geostrophic balance and hydrostatic balance implyDifferentiating horizontally and vertically gives the thermal-wind balanceThus vertical shear is determined by horizontal buoyancy gradients. The homogeneous Taylor–Proudman theorem is recovered only when those gradients vanish. Stable stratification allows sloping density surfaces and different geostrophic flow at different levels, without abandoning leading balance.

The Prandtl ratio of scales makes that modification quantitative. Horizontal and vertical terms in balanced potential-vorticity inversion have sizes $\psi /L^2$ and $(f_0^2/N^2)\psi /H^2$. Equality givesFor constant $N$, $Z=(N/|f_0|)z$ makes the inversion operator an ordinary three-dimensional Laplacian. Equivalently, an interior balanced disturbance of horizontal size $L$ naturally penetrates a depth of order $|f_0|L/N$. With $N\gg |f_0|$ this is a thin, vertically sheared structure, rather than the arbitrarily tall homogeneous columns suggested by strict Taylor–Proudman theorem behaviour. This is an aspect-ratio statement, not the diffusivity Prandtl number.

To allow a weak meridional variation of the Coriolis parameter, use the beta plane $f=f_0+\beta y$, with $\beta L/f_0=O(\mathrm{Ro})$, or equivalently $\beta L^2/U=O(1)$. The small correction to the leading balance is dynamically important on the slow time. Let $\xi ={\nabla }_H^2\psi$ be vertical relative vorticity and let $D_g={\partial }_t+{\mathbf{u}}_g\cdot {\nabla }_H$. At the first nontrivial order, the vertical vorticity equation and buoyancy equation reduce toRelative-vorticity stretching and tilting are higher order here; planetary-vorticity stretching $f_0{w}_z$ must remain. Eliminate $w$ to expose a materially conserved balanced scalar. Because ${\mathbf{u}}_{g,z}=(-b_y/f_0,b_x/f_0)$, the commutator term ${\mathbf{u}}_{g,z}\cdot {\nabla }_{H}b$ is exactly zero. For $N=N(z)$ this givesAdding the two equations proves the quasi-geostrophic potential-vorticity equationAn additive constant such as $f_0$ changes no dynamics. The horizontal streamfunction advection bracket can express this as $Q_t+{\psi }_x{Q}_y-{\psi }_y{Q}_x=0$. Although small departures from balance generated this evolution equation, the leading geostrophic flow itself advects $Q$ nonlinearly. The stretching term couples horizontal motion to vertical displacement of density surfaces.

This gives the practical advection-and-inversion description. Start with interior three-dimensional quasi-geostrophic potential vorticity, advect it with ${\mathbf{u}}_g$, and recover $\psi$ fromThe stratified quasi-geostrophic inversion is elliptic for positive $N^2$. It must be supplied with lateral and vertical boundary conditions and any circulation or mean-flow constraints; interior potential vorticity alone is not sufficient. At rigid horizontal lids, $w=0$ makes boundary buoyancy satisfy $D_{g}b=0$, and $b=f_0{\psi }_z$ supplies the vertical derivative data for inversion. Boundary buoyancy perturbations can therefore support flow even when the interior anomaly is zero. Once $\psi$ is found, its horizontal derivatives give velocity, its vertical derivative gives buoyancy, and the thermodynamic equation diagnoses $w$. This balanced evolution excludes the independent fast-wave initial data present in the full equations.

Finally, a gradient in $f$ supplies a restoring mechanism even in an otherwise uniform fluid. Linearizing about rest with constant $N$, the quasi-geostrophic potential-vorticity equation becomesSubstitution of $\psi \propto e^{i(kx+ly+mz-\omega t)}$ gives the Rossby wave dispersionFor $\beta >0$, its zonal phase speed is westward relative to a resting background. A meridionally displaced parcel retains potential vorticity, so it acquires a relative-vorticity anomaly that induces the velocity tending to return the disturbance. Background shear and stratification modify the full potential-vorticity gradient; they can alter propagation or permit instability. On a constant f-plane the planetary contribution vanishes, but gradients of relative vorticity, boundary buoyancy or topography can still support balanced wave motion. Close enough to the equator that $f_0$ vanishes, this mid-latitude quasi-geostrophic approximation and its scaling must be replaced.

The basic quasi-geostrophic streamfunction $\bar{\psi }=-\Lambda yz$ gives $\bar{u}=\Lambda z$, $\bar{v}=0$ and constant quasi-geostrophic potential vorticity $\bar{Q}=f$. DefineBecause the background has no potential-vorticity gradient, the linear potential-vorticity equation is simplyIts characteristics have constant $y,z,x-\Lambda zt$, so the general solution is $q^{\prime }=\mathcal{F}(y,z,x-\Lambda zt)$. In particular, $q^{\prime }=F(y,z)G(x-\Lambda zt)$ is the requested arbitrary separated family; one product is not the most general initial field.

Linearizing the vertical velocity requires retaining advection of ${\bar{\psi }}_z=-\Lambda y$ by $v^{\prime }={\psi }_x^{\prime }$. ThereforeThe impermeable boundary at $z=0$ imposesThis is conservation of boundary buoyancy, not merely a condition that ${\psi }_z^{\prime }$ vanish.

For zero interior potential vorticity, write ${\psi }^{\prime }=\mathrm{Re}\{\varphi (z,t)e^{i(kx+ly)}\}$. With $K_h=\sqrt{k^2+l^2}>0$, stratified quasi-geostrophic inversion givesThe decaying solution isNo time dependence has been assumed to obtain this vertical structure. Applying the boundary condition gives $-\mu \dot{A}-i\Lambda kA=0$, henceThus the Eady edge wave travels eastward with horizontal phase velocity $\omega /k=\Lambda /\mu$ for $k\ne 0$.

The propagation mechanism can be made explicit using boundary parcels. Thermal wind gives the background buoyancy gradient ${\bar{b}}_y=f{\bar{\psi }}_{zy}=-f\Lambda$. A small meridional parcel displacement $d$ produces boundary buoyancy perturbation $b^{\prime }=-d{\bar{b}}_y=f\Lambda d$. Meanwhile the decaying potential-vorticity inversion gives $b^{\prime }=f{\psi }_z^{\prime }=-f\mu {\psi }^{\prime }$ at the boundary. Consequently ${\psi }^{\prime }=-\Lambda d/\mu$ there, and the parcel equation $d_t=v^{\prime }={\psi }_x^{\prime }$ becomesA displaced buoyancy contour therefore induces a velocity field that advances that displacement eastward. No interior potential-vorticity gradient is needed because the boundary buoyancy supplies the invertibility data.

Vortex stretching is what sets the induced flow's vertical reach. In the interior, the zero-potential vorticity condition requiresBoundary isopycnal displacement decreases with height as $e^{-\mu z}$, so fluid columns are stretched or compressed over depth ${\mu }^{-1}=f/(N{K}_h)$. Conservation of potential vorticity makes this stretching balance the relative vorticity. The resulting inversion factor $b^{\prime }=-f\mu {\psi }^{\prime }$ gives the propagation speed: stronger density stratification reduces penetration and speed, while larger $f$ increases them. Zero interior potential vorticity therefore does not imply zero interior flow.

For imposed boundary vertical velocity, retain the same decaying structure and writeThe boundary relation becomes the forced Eady edge wave equationThe integrating factor $e^{i\omega t}$ gives, for ${\omega }_0\ne \omega$,The arbitrary homogeneous coefficient $B$ sets the initial disturbance. With $A(0)=0$, this is $A=C(e^{-i{\omega }_{0}t}-e^{-i\omega t})/[i(\omega -{\omega }_0)]$.

At resonance, ${\omega }_0=\omega$, integration instead yieldsAn initially undisturbed fluid has ${\psi }^{\prime }=Ct{e}^{-\mu z}\cos (kx+ly-\omega t)$. The imposed pumping reinforces the same boundary normal mode every cycle, with no dissipative loss, so the linear amplitude grows secularly. The divergence describes resonance of the undamped linear model; sufficiently large wave amplitude invalidates the assumed small-disturbance approximation.

Potential-vorticity inversion is the diagnostic recovery of the flow from its advected vortical scalar, using the balance relation and specified boundary conditions. Its usefulness is a separation of tasks: transport the scalar by the current velocity, invert to find the new velocity, and repeat. The inverse is spatially nonlocal even though the transport law is local. The four models differ chiefly in the operator being inverted and the information needed at boundaries.

In two-dimensional vortex dynamics, planar incompressibility permits $\mathbf{u}=(-{\psi }_y,{\psi }_x)$, and vertical vorticity is $\zeta ={\nabla }^2\psi$. The curl of inviscid momentum givesThere is no vortex stretching in strictly planar motion, so vorticity is materially conserved. The two-dimensional potential-vorticity inversion is the Poisson equation ${\nabla }^2\psi =\zeta$. On the plane, with suitable integrability and far-field behavior, its Green's function is $(2\pi {)}^{-1}\log |\mathbf{x}-{\mathbf{x}}^{\prime }|$, givingDifferentiation gives the planar Biot-Savart kernel. Here ${\psi }_h$ is a harmonic function fixed by boundary, circulation and far-field data. A rigid impermeable wall is a streamline, so $\psi$ is constant on it; constants or circulations around distinct holes require specification. The vorticity alone cannot distinguish two flows differing by a permitted harmonic contribution, such as a uniform current on an unbounded plane. A constant added to $\psi$ has no physical effect.

In nonrotating layerwise-two-dimensional vortex dynamics, strong stable density stratification constrains motion to almost horizontal layers. In the ideal independent-layer limit, vertical velocity vanishes at leading order, $(u,v)=(-{\psi }_y,{\psi }_x)$ andhold separately at each fixed $z$. Height is a parameter: the horizontal inverse Laplacian contains no vertical derivative. With background buoyancy $N^{2}z$, the leading Ertel potential vorticity is $Q=\mathbf{\omega }\cdot \nabla b\simeq N^2\zeta$, so the materially advected normalized scalar is $Q/N^2\simeq \zeta$. The same horizontal Poisson equation is inverted layer by layer, using each layer's boundary data. This is an asymptotic constrained model, not arbitrary three-dimensional Euler motion with vertical velocity set to zero: its validity requires small Froude number and controlled vertical scales. At finite buoyancy frequency, the weak vertical motion and tilted isopycnals restore coupling. The inversion may be independent by layer while the reconstructed fields vary with height. The horizontal Euler equations for an inviscid fluid determine the pressure gradient, while hydrostatic balance relates its vertical derivative to buoyancy.

In shallow-water quasi-geostrophic inversion, begin with the exact shallow-water potential vorticity $(f+\zeta )/h$. Linearizing its anomaly about depth $H$ gives, after multiplying by $H$, $\zeta -f_0\eta /H$. For small Rossby number, geostrophic balance gives $\psi =g\eta /f_0$, while $\zeta ={\nabla }_h^2\psi$. Keeping the beta plane background gradient givesThus $q-\beta y$ is inverted with the modified Helmholtz equation $({\nabla }_h^2-L_D^{-2})\psi =q-\beta y$, under the specified boundary and far-field conditions. The velocity follows from $\psi$ and surface elevation from $\eta =f_0\psi /g$. A horizontal Fourier mode of the anomaly satisfiesThe Rossby deformation radius therefore limits the reach of the balanced response. Layer stretching, represented by $-\psi /L_D^2$, is absent from ordinary planar vorticity inversion. Unlike the pure Poisson equation, the finite-deformation operator also responds to the constant part of $\psi$, which specifies mean elevation rather than an arbitrary streamfunction offset. Specified bottom topography adds a known source to $q$; it is subtracted along with the background before inversion.

In stratified quasi-geostrophic inversion, continuously varying vertical displacement replaces the shallow layer's single stretching term. For positive $N^2(z)$ and constant nonzero $f_0$, the three-dimensional quasi-geostrophic potential vorticity isGeostrophic balance and hydrostatic balance give velocity $(-{\psi }_y,{\psi }_x)$ and buoyancy perturbation $b=f_0{\psi }_z$. Invert the anisotropic elliptic boundary value problem using interior $q$, prescribed top and bottom buoyancy as Neumann boundary conditions, and the lateral, far-field and global-flow data. Boundary buoyancy evolves by its own material conservation when the boundary is impermeable. It cannot be discarded: the Eady edge wave in the preceding solution has zero interior potential vorticity anomaly but nonzero boundary buoyancy and a nonzero induced flow.

For constant $N,f_0$, the stretched coordinates for quasi-geostrophic inversion $Z=Nz/|f_0|$ turn the inversion operator into the three-dimensional Laplacian. Its infinite-space Green's function is $-1/(4\pi |\mathbf{X}-{\mathbf{X}}^{\prime }|)$ in $(x,y,Z)$ coordinates. Hence a horizontal disturbance of scale $L$ communicates vertically over scale $|f_0|L/N$, explicitly coupling layers. In a bounded domain with homogeneous difference data, multiplying the difference of two inversions by its streamfunction and applying integration by parts givesThis proves uniqueness up to a constant when the relevant boundary terms vanish; a mean or pressure reference fixes the remaining constant. Compatible boundary fluxes are needed for a purely Neumann boundary-value problem.

Across all four models, inversion provides the induced flow, while scalar transport supplies time evolution. The independent horizontal inverse in the first two models becomes a finite-deformation inverse in shallow water and a vertically coupled elliptic inverse in continuous stratification. The principle applies within the stated balance or constrained model. In the unrestricted Boussinesq equations or shallow water equations, gravity waves and divergent motion can share the same potential vorticity, so the vortical scalar alone does not reconstruct every possible flow.

For constant density, divide the incompressible Euler equations by $\rho$. The vector identity $(\mathbf{u}\cdot \nabla )\mathbf{u}=\nabla (|\mathbf{u}{|}^2/2)-\mathbf{u}\times \mathbf{\omega }$ converts them into ${\partial }_t\mathbf{u}=\mathbf{u}\times \mathbf{\omega }-\nabla (p/\rho +|\mathbf{u}{|}^2/2)$. Taking the curl eliminates the gradient. Expand the curl of the cross product, use $\nabla \cdot \mathbf{u}=0$ and $\nabla \cdot \mathbf{\omega }=0$, and obtainThis is the vorticity equation, including vortex stretching. Without constant density there would in general be an additional baroclinic term; the displayed source equations use the constant-density interpretation.

For the specified velocity, direct differentiation gives $\mathbf{\omega }=2\Omega (t)\mathbf{k}$. It is spatially uniform, while $(\mathbf{\omega }\cdot \nabla )\mathbf{u}=2\beta \mathbf{\omega }$. Thus ${\Omega }^{\prime }=2\beta \Omega$ andIf the positive $z$ direction agrees with the initial vorticity, $2\Omega (0)={\omega }_0$ as written in the question. For opposite rotation the initial component is signed and the vector has the corresponding minus sign.

For a material particle, $d(x^2+y^2)/dt=2x(-\beta x-\Omega y)+2y(-\beta y+\Omega x)=-2\beta (x^2+y^2)$, while $dz/dt=2\beta z$. The initial material circle therefore evolves intoThe azimuthal rotation changes the labels of points on the circle but not the circle as a set. Its counterclockwise circulation, viewed from positive $z$, isfor the positive-orientation convention. The strain contributes no tangential velocity. Thus circulation stays constant although vorticity grows and the enclosed material area shrinks, illustrating Kelvin's circulation theorem for a material loop in constant-density inviscid flow with no nonconservative body force.

The displayed evolution law uses the usual assumptions of incompressible flow, constant density and conservative body force. With $\mathbf{\omega }=\nabla \times \mathbf{u}$, the Euler equation isUse $(\mathbf{u}\cdot \nabla )\mathbf{u}=\nabla (|\mathbf{u}{|}^2/2)-\mathbf{u}\times \mathbf{\omega }$ and take a curl. Gradients have zero curl, andHere $\nabla \cdot \mathbf{\omega }=0$ because it is a curl, and $\nabla \cdot \mathbf{u}=0$ by incompressibility. The vorticity equation is thereforeThe material derivative combines local change with advection by fluid particles. The right-hand side is vortex stretching and tilting: spatial variation of velocity along a vortex line changes its strength and direction. Along each particle path this is a homogeneous linear equation for vorticity. For a smooth velocity field its solution with zero initial vorticity remains zero by uniqueness. Thus an initially irrotational flow remains irrotational while these assumptions and smoothness hold.

For a plane flow $\mathbf{u}=(u(x,y,t),v(x,y,t),0)$, the vorticity is $\mathbf{\omega }=(0,0,\zeta )$ and $(\mathbf{\omega }\cdot \nabla )\mathbf{u}=\zeta {\partial }_z\mathbf{u}=0$. Hence $D\zeta /Dt=0$. Every particle retains its initial value, so an initially uniform value givesInviscid motion alone is insufficient for the printed simplified law: compressible variable-density motion generally adds $-\mathbf{\omega }\nabla \cdot \mathbf{u}+\nabla \rho \times \nabla p/{\rho }^2$. Those terms vanish under the assumptions used above.

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.

Put $W=⟨{\omega }^2⟩$, $C=⟨|\nabla \times \mathbf{\omega }{|}^2⟩$ and $M=⟨G^3⟩$, where $G={\partial }_x{u}_x$. The small-$r$ Taylor series giveIn the Kármán-Howarth equation, the coefficient of $r^2$ on the left is $-\dot{W}/30$. On the right, the triple-correlation term contributes $7M/6$ and the viscous term contributes $\nu C/15$. HenceComparing with the mean enstrophy balance identifies the mean vortex stretching:Also $S_2=2(u^2-F)=W{r}^2/15+\cdots$, whereas differentiability gives $S_2=⟨G^2⟩r^2+\cdots$. Thus $⟨G^2⟩=W/15$. Homogeneity makes $⟨G⟩=0$, so the longitudinal velocity-gradient skewness is $S_0=M/(W/15{)}^{3/2}$. Substitution givesPositive mean vortex stretching therefore requires negative skewness. A centered Gaussian distribution has zero third moment, so the derivative statistics cannot be Gaussian. This conclusion uses the physically positive mean stretching, not homogeneity alone.

The probability density function has mean zero but a longer or stronger negative tail, representing intermittent compression. The illustration uses a standardized Gaussian mixture distribution solely to sketch this sign of skewness; it is not turbulence simulation data or a fitted universal distribution.

For unforced incompressible flow, multiply the vorticity equation by ${\omega }_i$. The antisymmetric part of the velocity gradient tensor drops out of ${\omega }_i{\omega }_j{\partial }_j{u}_i$, leavingStatistical homogeneity removes mean divergences. Integration by parts, together with $\nabla \cdot \mathbf{\omega }=0$, givesThe last equality holds for the averages, not pointwise. It follows from the usual curl/divergence identity after the divergence terms have averaged to zero. Vortex stretching produces enstrophy, while viscous diffusion destroys it.

Let $T=\ell /u$, $ϵ\sim u^3/\ell$, and use the Kolmogorov microscales $\eta =({\nu }^3/ϵ{)}^{1/4}$ and ${\tau }_{\eta }=(\nu /ϵ{)}^{1/2}$. Since $⟨{\omega }^2⟩=ϵ/\nu$, a slowly changing mean enstrophy has a time derivative of order $ϵ/(\nu T)$. The usual small-scale estimate givesThe relative imbalance is consequently ${\tau }_{\eta }/T\sim (\nu /u\ell {)}^{1/2}={\mathrm{Re}}^{-1/2}$, yieldingThe estimate assumes an established high-Reynolds-number cascade and integral-time evolution of the mean, not a rapid initial transient. In a forced flow the exact balance includes $⟨\mathbf{\omega }\cdot \nabla \times \mathbf{f}⟩$; forcing smooth on the integral scale of turbulence is smaller than the leading small-scale terms.

The high-Reynolds-number balance gives positive mean vortex stretching. The Betchov relation therefore gives $⟨abc⟩<0$. To connect this to a fixed longitudinal velocity derivative, additionally assume isotropic turbulence, or average over the direction of that derivative. For a unit vector $\mathbf{n}$, the derivative is $d={\mathbf{n}}^{T}S\mathbf{n}$. Uniform directional averaging, with $\mathrm{tr}S=0$, givesFor example these identities follow by averaging products of direction cosines on a sphere. Hence in isotropic turbulenceand division by the positive derivative variance to the power $3/2$ gives negative longitudinal velocity-gradient skewness. Homogeneity alone does not fix the orientation statistics and does not establish this sign for every chosen coordinate direction.

Order the principal strain rates as $a\ge b\ge c$. Except at degeneracies, $a>0$ and $c<0$, so $abc<0$ corresponds to $b>0$: two extensional directions and one compressive direction, or bi-axial strain. The negative average of $abc$ says that such configurations dominate this cubic amplitude-weighted statistic over axial-extension configurations with $b<0$. In the usual physical interpretation, comparably intense strain events are therefore biased toward bi-axial extension.

An unweighted frequency claim needs more information. As a counterexample to deducing it solely from the cubic average, assign strain eigenvalues $(1,1,-2)$ to 10% of samples and $(-0.1,-0.1,0.2)$ to 90%, with randomly oriented eigenvectors. Then $⟨abc⟩=-0.1982<0$, despite bi-axial extension occurring in only 10% of the samples. This example addresses what the moment inequality can infer, not a construction of a complete turbulent velocity field. A measured frequency preference is a further statistical observation.

The geometric significance is that a fluid element stretches in two directions and contracts in the third, becoming sheet-like while preserving volume. Distributed vorticity is swept into thin layers and its component along an extensional direction is amplified. The associated shear layers can subsequently roll up and produce vortex tubes. Thus the weighted strain bias helps explain sheet formation in the energy cascade even though the most intense dissipative structures can be tubes.

Assume inviscid, unforced evolution, so shallow-water potential vorticity is conserved. If a column has initial depth $H_1$ and relative vorticity ${\zeta }_1$, moving to $H_2>H_1$ at fixed Coriolis parameter givesIts absolute vorticity is amplified by vortex stretching. A column initially without relative vorticity develops cyclonic spin: positive in the Northern Hemisphere and negative in the Southern Hemisphere. More generally the increase is in the direction of its original absolute vorticity, not necessarily the direction of $f$ if the initial spin is unusually strong.

Conservation of shallow-water potential vorticity with fixed relative vorticity givesFor ordinary planetary-vorticity-dominated columns, moving towards the equator decreases the magnitude of absolute vorticity, so the column must become shallower and spread horizontally. This is the compression counterpart of vortex stretching. For $\zeta =0$, the depth ratio is $f_2/f_1$, in either hemisphere before reaching the equator. The displayed general formula is essential if the fixed relative vorticity is large: the sign of the depth change then depends on the original absolute vorticity, and conservation cannot continue through zero absolute vorticity with positive finite depth and nonzero conserved potential vorticity.

Let $\mathbf{\zeta }=\nabla \times \mathbf{u}$ and $\mathbf{\omega }=\mathbf{\zeta }+2\Omega {\mathbf{e}}_z$ be the relative vorticity and absolute vorticity. The nonlinear acceleration identity rewrites the momentum equation asTaking the curl removes the pressure, kinetic energy gradient and shearing-sheet tidal potential. ThusBoth $\nabla \cdot \mathbf{u}$ and $\nabla \cdot \mathbf{\omega }$ vanish. The vorticity equation is consequentlyFor the absolute-vorticity flux tensor $T_{ij}=u_i{\omega }_j-{\omega }_i{u}_j$, use the flux convention $(\nabla \cdot \mathbf{T}{)}_j={\partial }_i{T}_{ij}$, with the first index denoting transport direction. This givesThe scalar $\theta$ is a buoyancy displacement variable, so dimensional consistency gives it units of length; the printed potential-temperature terminology presupposes this normalization. Its curl source directly forces only the $y,z$ components of absolute vorticity. Hence radial vorticity conservation in a shearing sheet takes the local flux formIts volume integral is constant when the net boundary flux vanishes. However $D_t{\omega }_x=\mathbf{\omega }\cdot \nabla u_x$ still permits vortex stretching; the radial component is not generally a materially conserved scalar in three dimensions. Using divergence on the second index of the stated tensor would reverse the transport sign, so the tensor convention matters.

For the curl of a cross product, the Cartesian contraction isThe product rule then givesFor the fluid equation, replace the advective term by $\mathbf{w}\times \mathbf{u}+\nabla (|\mathbf{u}{|}^2/2)$ using the preceding identity. Take its curl. Both gradient terms have zero curl, while $\nabla \cdot \mathbf{w}=\nabla \cdot (\nabla \times \mathbf{u})=0$ and incompressibility gives $\nabla \cdot \mathbf{u}=0$. The vorticity equation is consequentlyThe left side is the material derivative of vorticity; the right is vortex stretching. It follows here because the given force per mass is a gradient, in addition to the flow being inviscid and incompressible.

For a material fluid volume, mass conservation yields ${\rho }_t+\nabla \cdot (\rho \mathbf{u})=0$. Constant mass density reduces this to $\nabla \cdot \mathbf{u}=0$. For inviscid flow, the Cauchy stress tensor is $-pI$; its net surface force is the volume integral of $-\nabla p$ by the divergence theorem. Newton's second law for each material volume therefore gives $\rho D\mathbf{u}/Dt=-\nabla p$. Since $D/Dt={\partial }_t+\mathbf{u}\cdot \nabla$, these are exactly the two stated equations.

In two dimensions they readDifferentiate the second momentum equation in $x$ and subtract the $y$ derivative of the first. The pressure terms cancel. For $\omega =v_x-u_y$, collecting terms givesIncompressibility removes the last term, givingThus each fluid parcel retains its scalar vorticity as it moves. In particular an initially irrotational flow remains irrotational in this smooth two-dimensional inviscid evolution. There is no vortex stretching in this two-dimensional setting.

Flow visualization and the idealization. A streamline is everywhere tangent to the instantaneous velocity field at one fixed time. A streakline is the locus, at one observation time, of particles that previously passed through a specified injection point. A marked particle's own trajectory is a pathline. These coincide geometrically in steady flow, but need not do so in an unsteady draining flow.

The intended ideal liquid is approximately incompressible and inviscid: viscosity and associated shear stresses are negligible on the scales and times being considered, the Mach number is small, and the density is effectively constant. This does not require zero vorticity; the initially rotating state is already a counterexample. It also does not justify ignoring viscosity in an arbitrarily thin boundary layer or at a real drain.

Spiralling flow and angular momentum. In an axisymmetric inviscid flow, pressure and gravity exert no torque about the vertical axis. A parcel therefore conserves axial specific angular momentum,The late, shallow remnant in the intended draining picture is supplied by parcels initially near the outer rim, so $r_0\simeq R$. The slow-drain angular-momentum amplification then givesThe assumption about the parcels' origin is necessary: torque conservation by itself does not assign $r_0=R$ to every parcel at every time.

A continuous dye feed near the outer free surface produces inward spirals that wind increasingly rapidly near the drain. The figure uses a late quasisteady model, $u_r=-D/r$ and $u_{\varphi }=C/r$ with positive constants $D,C$, giving $\varphi ={\varphi }_0-(C/D)\ln (r/r_0)$. Its streaklines and streamlines approximately coincide over an observation interval short compared with the overall draining time. Strictly, individually marked particles at the initial surface trace pathlines, not a continuously injected streakline.

The free-surface depression. Let $z=H(r)$ be the surface. Slow radial draining makes radial inertial corrections small compared with the azimuthal centripetal acceleration. The radial Euler equation and vertical hydrostatic equilibrium giveThus $g{H}^{\prime }(r)=u_{\varphi }^2/r=C^2/r^3$, with $C={\omega }_0{R}^2$. Integrating from $r$ to the rim givesFor $r\ll R$ this becomesThe surface rises outward because the inward pressure gradient force must sustain the circular motion. The dimensions are metres, and the $1/r^2$ profile is not valid inside the finite drain or vortex core.

The terrestrial estimate. Use the vertical component of Earth's rotation as the ideal initial angular speed, ${\omega }_0\simeq {\Omega }_E\sin \lambda$. At Cambridge, taking $\lambda \simeq 52.2^{\circ }$ and ${\Omega }_E=7.292\times 10^{-5}{\mathrm{s}}^{-1}$, this gives ${\omega }_0\simeq 5.76\times 10^{-5}{\mathrm{s}}^{-1}$. Hence, for $R=1\mathrm{m}$, $r=0.01\mathrm{m}$ and $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$,Taking the whole terrestrial angular speed rather than its local vertical projection would give about $2.7\mu \mathrm{m}$. The Cambridge estimate uses the projection appropriate to the vertical vortex. This very small head assumes terrestrial rotation supplies the only initial swirl; imposed residual rotation would change it quadratically.

Curl of the azimuthal velocity. For axisymmetric $u_{\varphi }=r\omega (r)$ with no vertical dependence, the rotational part of the vorticity is vertical:The ideal quasisteady radial sink has no curl outside its core. The source's symbol $\Omega$ denotes this vertical vorticity component, not the fluid's angular velocity. In general its magnitude is the absolute value of the bracket.

The two limiting states. For initial uniform rotation, $\omega ={\omega }_0$, whereas the late approximation has $\omega =C/r^2$. ThereforeThe late flow can have large angular velocity and nonzero circulation while being locally irrotational. The drain core is excluded from the domain of the $1/r$ velocity field.

What is actually conserved. Inviscid does not mean that each parcel's vorticity magnitude is fixed in three dimensions. For an incompressible barotropic liquid with conservative body forces, the vorticity equation isVertical compression of a material column changes its vorticity through vortex stretching, even while its volume is conserved by horizontal expansion. In the columnar shallow-layer limit, vorticity dilution in a shallow draining layer givesThus a remnant much shallower than its original column can have much smaller vertical vorticity. At finite depth there is generally a residual value; the exactly irrotational late profile is a leading limiting approximation, not an exact parcelwise conversion of nonzero vorticity to zero in finite time.