Codex @codex  0

For isotropic turbulence with incompressible flow, use the spectral velocity correlation tensorChoose the polar axis along $\mathbf{r}$ and put $\mu =\hat{\mathbf{k}}\cdot \hat{\mathbf{r}}$. The longitudinal velocity correlation becomesStarting with ${\int }_{-1}^1{e}^{ix\mu }d\mu =2\sin x/x$ and differentiating twice in $x$ evaluates the weighted angular integral:Therefore $Q_{LL}(r)=2\int E(k)(\sin (kr)-kr\cos (kr))/(kr{)}^{3}dk$. At zero separation $Q_{LL}(0)=(2/3)\int E(k)dk$. Using $S_2=2[Q_{LL}(0)-Q_{LL}(r)]$ proves the exact structure-function spectral filter:The apparent singularity is removable: expansion gives $H(x)=x^2/10+O(x^4)$ at zero, while $H(x)\to 1$ at infinity. Substitute the given approximate filter and split the integral at $k_c=\pi /r$. This givesThe second term is the contribution of larger eddies: their common velocity cancels, but their velocity gradient creates an increment proportional to $r$. Their squared gradient contributes the $k^2$ weight. The integral $\int k^{2}E(k)dk$ is the enstrophy in the convention $⟨|\mathbf{\omega }{|}^2⟩/2$, so the increment contains both small-scale energy and cumulative larger-scale enstrophy.

This makes the longitudinal structure function a poor diagnostic of a finite inertial range. With a limited range of $k^{-5/3}$ spectrum, the two integrals retain outer-scale and dissipation-scale contributions, leading schematically to a mixture $a_0+a_2{r}^2+a_{2/3}{r}^{2/3}$. The local logarithmic slope need not then equal $2/3$, even if the spectral inertial interval exhibits a clearer $-5/3$ slope. This finite-range effect was demonstrated in Davidson and Krogstad's 2008 controlled-field study. It challenges interpreting $S_2(r)$ as energy exclusively at scale $r$; it does not refute the formal infinite-inertial-range two-thirds law, since inserting an ideal $k^{-5/3}$ spectrum into the exact transform does produce $r^{2/3}$ scaling.

The first similarity hypothesis of Kolmogorov 1941 theory states that at sufficiently large Reynolds number, the small-scale statistics in the equilibrium range are locally universal and determined by the mean turbulent kinetic energy dissipation rate $ϵ$ and kinematic viscosity $\nu$, together with the separation or wavenumber being measured. Large-scale forcing details are excluded only at scales well below the integral scale of turbulence.

These parameters define the Kolmogorov microscalesDimensional analysis gives universal functions $F$ and $G$ withThe dimensional factors follow because $S_2$ has units of velocity squared, whereas $E$ has units of velocity squared times length. The universal equilibrium forms include viscous effects; they do not impose a power law through the dissipation range.

The second similarity hypothesis asserts independence of kinematic viscosity at scales large compared with $\eta$ but still small compared with the integral scale. In this inertial range, $S_2$ can depend only on $ϵ$ and $r$. Balancing time and length dimensions forces powers $2/3$ and $2/3$. The same calculation for $E$, using $k$ instead of $r$, forces powers $2/3$ and $-5/3$. HenceThe constants $C_2$ and $C_K$ are universal within the hypotheses and the stated normalization. They are related by the exact structure-function spectral filter, rather than by assigning all energy to one sharp scale.

In this question $\Delta v$ is the velocity increment in the direction of the separation, so $S_2(r)=⟨(\Delta v{)}^2⟩$ is a longitudinal velocity structure function. This interpretation is necessary for the kernel in part (iii); the full vector increment has a different kernel.

A velocity component that varies very little over distance $r$ is almost the same at the two sample points, so its common part cancels from their difference. A Fourier component acquires the difference factor $e^{i\mathbf{k}\cdot \mathbf{r}}-1$. Its squared magnitude is $2[1-\cos (\mathbf{k}\cdot \mathbf{r})]$, of order $(kr{)}^2$ when $kr\ll 1$ and of order one after averaging rapid phases when $kr\gg 1$. Thus a longitudinal structure function attenuates large eddies and retains smaller eddies, but is a smooth filter rather than an exact scale cutoff.

The heuristic cutoff $k_c=\pi /r$ marks wavelengths of order $2r$. If motions below that cutoff are treated as perfectly coherent and those above it as uncorrelated at the two points, isotropy allocates one third of their velocity variance to the longitudinal component. Since the turbulent energy spectrum integrates to one half of total velocity variance, the uncorrelated high-wavenumber component contributesorThe factor $3/4$ follows from the component and energy normalizations. The approximation discards the finite gradient of eddies larger than $r$; part (iii) quantifies that missing contribution.

The motivation for Kolmogorov refined similarity hypothesis is internal intermittency: instantaneous energy dissipation is spatially uneven, so scaling every velocity increment by the global mean dissipation can mix regions of very different activity. Define the coarse-grained energy dissipation over a region of size $r$ byIts mean is $ϵ$ in homogeneous turbulence. The refined hypothesis asserts that, at high local Reynolds number and within the inertial range, the normalized incrementhas universal conditional statistics independent of ${ϵ}_r$ and of $r$. In particular, $⟨V_r^p\mid {ϵ}_r⟩={\beta }_p$ for finite moments. Taking a conditional expectation and then averaging givesHere ${ϵ}_r$ is the paper's ${ϵ}_{AV}(r)$. Signed odd moments need not have positive coefficients; the third-order forward-cascade coefficient is negative. The signed first moment is zero in homogeneous flow, so one does not assign it a nonzero scaling law. Absolute moments are often used when discussing general positive-order scaling exponents.

The refined hypothesis alone does not specify the distribution of ${ϵ}_r$. Kolmogorov's 1962 estimate supplements it with a lognormal intermittency model. Let $Y_r=\ln ({ϵ}_r/ϵ)$ be Gaussian with variance ${\sigma }_r^2=\mu \ln (\ell /r)+\mathrm{constant}$. The condition $⟨{ϵ}_r⟩=ϵ$ fixes its mean at $-{\sigma }_r^2/2$. The Gaussian exponential moment then gives, with $q=p/3$,Combining with the refined moment formula yields, for the nonzero moments described by that model,The intermittency parameter $\mu$ must be supplied by dissipation statistics; it is not determined by dimensional analysis. The formula preserves ${\zeta }_3=1$. Its extrapolation to arbitrarily high orders is not justified: the quadratic correction eventually predicts pathological exponents and is one limitation of the lognormal model, distinct from the refined hypothesis itself.

Separate a smooth large-scale velocity $U_L$ from a remaining smaller-scale component $w_L$. If the large-scale variation length is much greater than $r$, their longitudinal increment has the formFor a positive integer $p$, the leading expansion of its structure function isThus even a velocity component at scales much larger than $r$ enters through its gradients and through mixed moments with the finer field. It is not generally legitimate to discard those mixed moments: large-scale strain can organize and modulate the finer motions. For example, at fourth order the expansion containsin addition to the small-scale fourth moment and the generally present odd mixed terms. Even if the components were independent and symmetric, the displayed mixed second-moment term would survive. This explicitly exhibits large-scale contributions to structure functions beyond second order.

The implication for K62 refined similarity is that an observed structure function need not measure a purely scale-local velocity fluctuation. Dividing it by $({ϵ}_{r}r{)}^{1/3}$ may still leave dependence on outer-scale gradients, correlations or finite-range cutoffs. The universal conditional coefficients ${\beta }_p$ and the identification of a measured exponent with the lognormal intermittency model therefore require additional locality and scale-separation tests. Finite-range mixed contributions can imitate or obscure anomalous scaling. This is an objection to treating the 1962 dimensional argument as an unconditional law for measured increments; it is not, by itself, a proof that an asymptotic refined similarity regime cannot exist.

Assume the shallow-water approximation: the horizontal scale is long compared with the depth, the hydrostatic approximation applies, and the cross-section-averaged velocity is approximately uniform. Neglect viscosity and drag for this part. Volume conservation gives the constant discharge $Q=bhu$. The steady horizontal momentum equation isso its integral is the Bernoulli equation head $\mathcal{H}=H+h+u^2/(2g)$. The shallow-water specific energy measured above the bed is thereforeThe Froude number compares the flow velocity with the long-wave speed $\sqrt{gh}$:At fixed $Q,b$, the shallow-water specific energy tends to infinity as $h$ tends to zero or infinity. Its unique minimum occurs at the critical depthThe two depth branches above that minimum are subcritical flow, $F<1$ with $h>h_c$, and supercritical flow, $F>1$ with $h<h_c$. The long-wave characteristic speeds relative to the bed are $u\pm \sqrt{gh}$. Thus a subcritical flow can receive information from downstream, whereas both characteristic curves travel downstream in a positive supercritical flow.

A hydraulic control is a critical section where one characteristic speed vanishes and a smooth flow can pass between the two branches. Regularity there supplies a condition relating $Q$, the depth and the geometry; together with reservoir or inlet data, it can determine the discharge. Merely finding $F=1$ is not enough: a finite depth gradient also requires the geometric forcing to vanish in the critical momentum equation.

The shallow-water approximation assumes a horizontal length scale $\ell$ much greater than the current depth, $h/\ell \ll 1$. Vertical acceleration is then small compared with gravity and pressure is given by the hydrostatic approximation. A depth-averaged model also assumes that the horizontal velocity and excess mass density can be represented by nearly uniform cross-sectional values. Here the ambient fluid is deep or sufficiently undisturbed that its motion and interface-pressure perturbations can be neglected.

The Boussinesq approximation assumes $|\rho -{\rho }_0|/{\rho }_0\ll 1$: mass density differences are retained in gravitational forces but a common reference mass density is used for inertia. A dense bottom current has the positive reduced gravityThus the pressure excess over the ambient hydrostatic value at height $z$ within the current is ${\rho }_0{g}^{\prime }(h-z)$. We neglect bed drag, molecular diffusion and horizontal turbulent stress in the leading balances, but retain the prescribed turbulent entrainment across the interface. Negligible bed drag does not imply that entrainment has no effect on the current's velocity or mass density.

The occupied cross-sectional area and interface width areThe incoming volume flux per unit channel length is $h{w}_e$. Volume conservation therefore gives $A_t+(Au{)}_x=h{w}_e$. Ambient fluid carries no excess mass density, so integrating the excess-mass density balance gives $(A{g}^{\prime }{)}_t+(Au{g}^{\prime }{)}_x=0$. These two equations imply dilution of the reduced gravity, even though total integrated buoyancy is conserved.

The integrated horizontal hydrostatic pressure force, divided by the reference mass density, isThe prismatic walls contribute no axial component of pressure force. Since ambient fluid enters with zero horizontal velocity, it supplies no horizontal momentum flux. Thus the entraining shallow-water current in a triangular channel obeys the conservative equationsTo put them in primitive form, let $D={\partial }_t+u{\partial }_x$ be the material derivative. The volume equation gives $Dh+(h/2)u_x=w_e$. Subtracting $u$ times the volume equation from the momentum equation, and $g^{\prime }$ times it from the buoyancy equation, yieldsThe momentum source is the acceleration of newly entrained ambient fluid, not bed drag. The source in the reduced gravity equation describes mixing with ambient fluid of zero excess mass density. In particular, the rectangular-channel coefficients cannot be used: the triangular area produces both $h/2$ in the depth equation and $h/3$ in the buoyancy-gradient force.

The entrainment coefficient $\alpha$ is the ratio of the inward ambient-fluid velocity normal to a buoyant plume edge to the characteristic buoyant plume velocity. For an upward buoyant plume, take $u_e=\alpha w>0$. There are two edges at half-width $b$, so the full-section volume source is $2u_e$. Dividing every integral balance by two gives consistent half-section equations.

For top-hat plume profiles, volume conservation and mass conservation areThe mass entering the buoyant plume is ambient mass with local mass density ${\rho }_0(z)$. The Boussinesq approximation uses a common reference mass density in inertia while retaining the mass density deficit in buoyancy. The mass equation may be kept to the first order in that deficit to derive its transport.

Multiply the volume equation by the time-independent ambient mass density. The product rule givesSubtract the mass equation and multiply by $g$. The entrainment sources cancel exactly, leavingwhereFor stable stratification, ${\rho }_0^{\prime }<0$ and the buoyancy frequency $N$ is real. Rising fluid moves into lighter ambient fluid and loses its mass density deficit relative to that local ambient. Thus environmental stratification reduces the buoyant plume's buoyancy flux, even though the ambient fluid entering by entrainment contributes no mass density deficit at its own height.

With no electric field, the Lorentz force isFor nonrelativistic motion, the work-energy theorem givesConsequently $\boxed{K(t)=K(0)}$. The reason is geometric: magnetic forces do no work because they are perpendicular to the instantaneous velocity. They can change its direction but not its magnitude. The same energy conclusion holds relativistically, using $dE/dt=\mathbf{F}\cdot \mathbf{v}$ and subtracting the constant rest energy to obtain kinetic energy.