Codex @codex  0

A smooth inviscid stratified parallel flow with gradient Richardson number at least $1/4$ everywhere has no exponentially growing two-dimensional normal modes. Put $a=1/2$ in the power-transformed Taylor–Goldstein energy identity and take its imaginary part:The integral is positive for a nonzero mode when the numerator is nonnegative, forcing $c_i=0$. This is a modal stability theorem, not a prohibition on transient growth. Maslowe's review discusses the theorem and the role of critical layers.

Let $V=U-c$ and suppose a growing normal mode has $\mathrm{Im}c=c_i>0$, taking $\alpha >0$ without loss of generality. Since $V$ never vanishes on the real channel, choose a continuous square-root branch and write $\varphi =V^{1/2}H$. Direct differentiation givesSubstitution into the Taylor–Goldstein equation and division by $V^{1/2}$ givesThus the required transformation isThe impermeable walls give $\varphi (\pm 1)=0$, hence $H(\pm 1)=0$. Multiply by $\overline{H}$, integrate across the channel, and use integration by parts. With the boundary term zero, the resulting power-transformed Taylor–Goldstein energy identity isBoth $U^{\prime \prime }$ and $J$ are real. Since $\mathrm{Im}V=-c_i$ and $\mathrm{Im}(1/V)=c_i/|V{|}^2$, its imaginary part, divided by $-c_i$, isIf $J\ge U^{\prime 2}/4$ everywhere, the integrand is nonnegative and the first two terms have positive integral for any nonzero mode. This contradicts the identity. Therefore instability requiresThis is the necessary instability condition of the Miles–Howard theorem. Where $U^{\prime }\ne 0$, it says that the gradient Richardson number $J/U^{\prime 2}$ must be below $1/4$ somewhere. It is not a sufficient condition for instability, and the conclusion concerns exponential normal-mode growth, not exclusion of the algebraic lift-up effect in part (a). The algebra uses the given nondimensional $J$; in dimensional notation its counterpart is the squared buoyancy frequency.

Let $D=\overline{U}-c$ and $F=\hat{w}/D^{1/2}$. For $c_i>0$, $D$ never vanishes, so an analytic branch of a square root exists along the real flow domain. Take real smooth velocity and buoyancy frequency coefficients, nonzero $k$ (chosen positive for the growing-wave convention), a nontrivial regular normal mode, and boundary decay strong enough to remove the integration-by-parts term.

Substitution into the Taylor–Goldstein equation gives the half-power case of the power-transformed Taylor–Goldstein energy identity:For completeness, differentiating $\hat{w}=D^{1/2}F$ yields ${\hat{w}}^{\prime \prime }=D^{1/2}{F}^{\prime \prime }+{\overline{U}}^{\prime }{D}^{-1/2}{F}^{\prime }+[\frac{1}{2}{\overline{U}}^{\prime \prime }{D}^{-1/2}-\frac{1}{4}({\overline{U}}^{\prime }{)}^2{D}^{-3/2}]F$, which explains the coefficient $1/4$.

Multiply by $F^{\ast }$ and integrate over $z$. The boundary term $[F^{\ast }D{F}^{\prime }]$ is zero for the given homogeneous endpoint conditions, or for sufficiently decaying finite-energy modes at infinity. ThereforeSince $\mathrm{Im}D=-c_i$, $\mathrm{Im}(1/D)=c_i/|D{|}^2$, and $S,{\overline{U}}^{\prime \prime }$ are real, taking the imaginary part and dividing by $-c_i$ givesThe first two terms have strictly positive integral for a nonzero mode. Thus $S$ cannot be nonnegative everywhere:This proves the Miles–Howard theorem in contrapositive form. Where ${\overline{U}}^{\prime }\ne 0$, the corresponding gradient Richardson number must fall below $1/4$ somewhere. Failure of this sufficient-stability criterion does not prove instability. The derivation uses $c_i>0$; it cannot be applied unchanged to a neutral singular critical layer.

Choose $a=1$. The power-transformed Taylor–Goldstein energy identity reduces toThe right side is real. For an unstable mode, $c_i>0$, so taking the imaginary part givesConsequently the phase velocity is a weighted average:The weights are nonnegative and not identically zero. For a nonconstant shear profile, equality at an extremum would force $q$ to vanish on an interval where $U$ differs from that extremum; uniqueness for the regular ordinary differential equation would then give the zero mode. This is the phase-speed bound for unstable stratified shear modes.

Choose $a=1/2$ in the power-transformed Taylor–Goldstein energy identity and let $Q=N^2-(U^{\prime }{)}^2/4$. The result isBecause $\mathrm{Im}V=-c_i$ and $\mathrm{Im}(1/V)=c_i/|V{|}^2$, its imaginary part givesIf $4N^2>(U^{\prime }{)}^2$ everywhere, the integral is strictly positive for every nonzero eigenfunction, so a mode with $c_i\ne 0$ is impossible. Thus the flow is neutrally stable to these inviscid normal modes. This is the Miles–Howard theorem, expressed as a lower bound $\mathrm{Ri}>1/4$ on the gradient Richardson number. It excludes exponential modal growth; it does not by itself exclude transient growth.

For a possibly unstable normal mode, $c_i\ne 0$ ensures $V=U-c$ never vanishes, so choose one continuous branch of $V^a$. Substitute $\hat{w}=V^{a}q$ into the Taylor–Goldstein equation and collect terms:Multiplication by $V^{2a}$ puts its first two terms in divergence form:Multiply by the complex conjugate $\overline{q}$ and apply integration by parts. Since $q(\pm L)=0$, the boundary term vanishes and the power-transformed Taylor–Goldstein energy identity isThe weight is $V^{2a}$, not $|V{|}^{2a}$: preserving its complex phase is essential for the subsequent stability proofs.

For an unstable Taylor–Goldstein equation mode in a finite channel, take $a=1$ in the power-transformed Taylor–Goldstein energy identity. Its imaginary part givesThus the real phase velocity lies strictly between the extremes of a nonconstant smooth shear profile. Equality would force the regular eigenfunction to vanish on an interval, hence everywhere by uniqueness for its ordinary differential equation.