Hydrodynamic stability

Hydrodynamic stability studies whether small disturbances to a fluid base state decay, remain neutral, or grow. Linear normal modes reduce this question to eigenvalues or growth rates.

A global hydrodynamic mode is a separated solution $a(x)e^{-i\omega t}$ satisfying the complete spatial boundary conditions of a linearized flow problem. Its eigenvalue depends on the spatial domain and variable coefficients, rather than solely on a local dispersion relation.

A spatially nonuniform system is globally unstable when a normal mode satisfying its global spatial boundary conditions has positive temporal growth rate. This differs from purely local temporal amplification or convective wave-packet instability. A locally absolutely unstable region can be necessary without being sufficient for a growing global mode of finite spatial extent.

For the whole-line linear complex Ginzburg-Landau equation with real positive diffusion $\gamma$ and $\mu (x)={\mu }_0-\nu {ϵ}^2{x}^2$, removing drift and scaling $\xi =\sqrt{ϵ}(\nu /\gamma {)}^{1/4}x$ gives the Hermite differential equation. Decay at both infinities quantizes the spectral parameter to $2j+1$. The strongest mode is $j=0$, giving global instability iff ${\mu }_0>U^2/(4\gamma )+ϵ\sqrt{\gamma \nu }$. The confinement correction vanishes in the slow-profile limit, recovering the homogeneous absolute wave-packet instability threshold.

A critical layer is a region where the leading wave phase speed coincides with the local shear-flow velocity and an outer approximation degenerates. Small viscosity or other corrections must then be retained on a distinguished transverse scale. The width depends on the order of vanishing of $U-c_0$ and on the governing mode equation; simple-shear and stationary-profile layers have different balances.

For a monotone inviscid parallel shear flow, let a regular neutral Rayleigh equation for inviscid shear flow mode have $U(0)=c_0$, $U^{\prime \prime }(0)=0$, $d=U^{\prime }(0)>0$, and $p={\varphi }_0(0)\ne 0$. A perturbation $c=c_0+\epsilon c_1$ with $\mathrm{Im}{c}_1>0$ resolves the singularity on the inner variable $\eta =y/\epsilon$. The first outer correction has ${\varphi }_1=a+\beta y\log |y|+b_{\pm }y+\cdots$, with $\beta =c_1{U}^{\prime \prime \prime }(0)p/d^2$. The inner particular solution contains $(d\eta -c_1)\log (d\eta -c_1)$. Its complex logarithm has argument tending to zero at positive infinity and to $-\pi$ at negative infinity, giving the displayed matching jump.

Let $I=\int {\varphi }_0^{2}dy$, $J=\mathrm{PV}\int U^{\prime \prime }{\varphi }_0^2/(U-c_0{)}^{2}dy$, and $M=U^{\prime \prime \prime }(0){\varphi }_0(0{)}^2/U^{\prime }(0{)}^2$. The Wronskian identity for the first-order Rayleigh equation for inviscid shear flow, combined with neutral-mode critical-layer matching, gives $2k_0{k}_{1}I=-c_1(J+i\pi M)$. For real $k_1$, its imaginary part is $2k_0{k}_{1}I\pi M/(J^2+{\pi }^2{M}^2)$. The Cauchy principal value is necessary when the outer forcing has a simple real pole.

For the Prandtl normal-mode equation with viscosity scaled to one and a profile having a nondegenerate stationary point, $U-c_0\sim (y-y_c{)}^2/2$. A piecewise outer mode proportional to $U-c_0$ has a second-derivative jump. Balancing the advective and third-derivative terms smooths that jump over width ${\alpha }^{-1/4}$, with velocity amplitude ${\alpha }^{-1/2}$ and a phase-speed correction of the same order. Matching selects the eigenvalue through a third-order inner problem.

A stationary point in a non-monotone velocity profile can support a quarter-power critical layer whose phase-speed correction has positive imaginary part. The resulting normal mode grows at a rate proportional to the square root of the wavenumber. Such unbounded high-frequency amplification demonstrates linear instability of the frozen layer; nonlinear evolution and precise well-posedness consequences require separate analysis.

For a wave with zonal phase velocity $c$ in a basic shear $U(z)$, a critical level is a height $z_c$ at which $U(z_c)=c$. The Doppler-shifted material frequency then vanishes. This means equality of the laboratory velocities; the $U$ and $-c$ terms cancel in the wave frame. It does not mean $U=-c$. A semi-infinite Eady edge wave with zero interior potential vorticity anomaly can remain regular through its steering height.

Baroclinic instability converts potential energy of a stratified sheared basic flow into growing disturbances. In the Eady model, coupled Boundary Rossby waves provide a useful mechanism: their induced velocities displace boundary buoyancy anomalies, and suitable wavelengths permit phase locking. A single neutral boundary wave does not by itself imply exponential instability.

A high-Reynolds number long-wave mode near a wall has a viscous layer of thickness $\delta =(\alpha R{)}^{-1/3}$, an order-one region following the base profile, and a pressure-decay region of thickness $1/\alpha$. The lower Airy function solution fixes velocity displacement and fluid pressure. The middle Rayleigh equation for inviscid shear flow transports these to the outer constant-velocity region. Matching its decaying pressure produces an Airy wall-layer matching determinant. For order-one amplitudes, constant leading middle pressure requires ${\alpha }^2\ll \delta$; the distinguished balance $\alpha \sim \delta$ gives $\alpha \sim R^{-1/4}$.

With the Airy displacement integrals, the lower pressure-to-velocity ratio is $-iq\mathrm{Ai}\prime (-qC)/I_0(C)$. The decaying upper mode gives that ratio as $\alpha /\delta$. Hence for $\alpha =k{R}^{-1/4}$ the dispersion relation is $k^{4/3}{I}_0(C)+iq\mathrm{Ai}\prime (-qC)=0$. Writing it without division remains meaningful if an individual moment vanishes. The wavespeed is $c=\delta C$ and a positive wavenumber grows temporally when $\mathrm{Im}C>0$.

For an incompressible parallel flow with base velocity $U(y)$ and perturbation proportional to $e^{i\alpha (x-ct)}$, divergence of the linearized momentum equations gives $p\prime \prime -{\alpha }^{2}p=-2i\alpha U\prime v$. At a no-slip boundary condition, streamwise momentum gives $p(0)=u\prime \prime (0)/(i\alpha R)$. The equation encodes how base shear generates perturbation fluid pressure.

Relative to the kinetic-energy inner product, integration by parts reverses incompressible base transport, transposes the velocity-gradient coupling, and preserves viscous diffusion. With reverse time $\tau =-t$, the adjoint base state is evaluated at $-\tau$. The Leray-Helmholtz projection or adjoint pressure enforces divergence-free adjoint velocity. The curl of a cross product and vorticity cross-product identity give the equivalent form ${\partial }_{\tau }v=\Omega \times v-\nabla \times (U\times v)-\nabla \pi +{\mathrm{Re}}^{-1}\Delta v$.

For a constant reference density, fluid perturbation kinetic energy is proportional to $\Vert u{\Vert }^2/2$ under the displayed $L^2$ inner product. An adjoint operator for perturbation evolution must be defined relative to this inner product, including boundary conditions and the projection onto incompressible flow.

For the linear complex Ginzburg-Landau equation with real $U,c_d,\mu$ and positive real diffusion coefficient scaled to one, the maximum temporal growth rate is $\mu$ and the absolute growth rate is $\mu -U^2/[4(1+c_d^2)]$. For $U>0$, $\mu <0$ is stable, $0<\mu <{\mu }_a$ has convective hydrodynamic instability, and $\mu >{\mu }_a$ has absolute hydrodynamic instability. Both equality boundaries are marginal in exponential rate.

On a ray $x/t=V$, a localized wave packet is governed by a physically selected saddle point satisfying $d\omega /dk=V$ in its analytic dispersion relation. The exponent's real growth is $\mathrm{Im}({\omega }_{\star }-V{k}_{\star })$. The fixed-frame case $V=0$ gives the absolute growth rate. This ray formulation makes the frame dependence of convective hydrodynamic instability explicit.

A convective hydrodynamic instability amplifies an advected localized disturbance while its amplitude decays at any fixed position. The distinction from absolute hydrodynamic instability depends on the observation frame.

An absolute hydrodynamic instability makes the response to a localized impulse grow at a fixed spatial point in a specified frame. This differs from convective hydrodynamic instability, which amplifies a moving wave packet while its response decays at fixed positions. For analytic dispersion relations, the relevant absolute frequency is associated with an accessible zero-group velocity saddle, with the spatial-branch selection checked rather than assumed for an arbitrary dispersion relation.

For ${\eta }_t+\mathbf{c}\cdot \nabla \eta =\mu \eta +\nabla \cdot (\Gamma \nabla \eta )$ with symmetric positive definite $\Gamma$, the dispersion relation is $\omega =\mathbf{c}\cdot \mathbf{k}+i(\mu -{\mathbf{k}}^T\Gamma \mathbf{k})$. Its zero-group velocity saddle is ${\mathbf{k}}_0=-i{\Gamma }^{-1}\mathbf{c}/2$. The anisotropic heat kernel has fixed-position exponential rate $\mu -{\mathbf{c}}^T{\Gamma }^{-1}\mathbf{c}/4$, so this accessible saddle gives the actual threshold. At fixed speed $q$, the threshold ranges from $q^2/(4d_{\max })$ to $q^2/(4d_{\min })$, where $d_{\min },d_{\max }$ are the eigenvalues of $\Gamma$. Absence of absolute hydrodynamic instability in every direction requires the upper bound $\mu \le q^2/(4d_{\max })$.

The absolute frequency is the generally complex angular frequency at the physically selected absolute wavenumber. Its imaginary part is the fixed-position exponential growth rate under the convention $e^{-i\omega t}$. Marginal zero exponential growth may still carry an algebraic prefactor.

A time-independent ray Hamiltonian conserves the stationary-observer absolute frequency. A horizontally uniform background also conserves horizontal wavenumber. The intrinsic frequency $\hat{\omega }=\omega -kU(z)$ nevertheless changes when the ray moves through shear. These follow from the Hamiltonian ray-tracing equations and are essential to identifying a critical level of an internal gravity wave.

The absolute growth rate is the exponential growth rate of a localized impulse at a fixed position in the chosen observation frame. Under the convention $e^{i(kx-\omega t)}$, it is the imaginary part of the physically selected absolute frequency. Positive values give absolute hydrodynamic instability; negative values can still coexist with convective hydrodynamic instability.

An absolute wavenumber is a complex wavenumber at the selected zero-group velocity saddle of a dispersion relation. For the linear complex Ginzburg-Landau equation, $k_0=-U/[2(c_d-i)]$; the positive real part of the diffusion coefficient makes the Gaussian saddle accessible from the real Fourier contour.

The Orr-Sommerfeld equation governs the wall-normal velocity of normal mode disturbances to an incompressible viscous parallel shear flow. Here $D=d/dy$, ${\kappa }^2={\alpha }^2+{\beta }^2$, and $Re$ is the Reynolds number. It follows by eliminating pressure from the linearized Navier-Stokes equations. At a rigid no-slip wall, $\hat{v}=D\hat{v}=0$. In the inviscid limit it reduces to the Rayleigh equation for inviscid shear flow.

For temporal normal-mode stability of a viscous parallel shear flow, any unstable three-dimensional mode maps by the Squire transformation to an unstable two-dimensional mode at a no larger Reynolds number. Hence the first temporal modal instability occurs among two-dimensional disturbances. The theorem does not exclude three-dimensional transient growth or forced responses.

For positive streamwise wavenumber $k$, replace a three-dimensional Orr-Sommerfeld equation by a two-dimensional one with the displayed parameters. Both the transverse operator and the coefficient $1/(ikR)$ remain unchanged, preserving the eigenfunction and phase-speed eigenvalue. This proves the modal comparison in Squire's theorem.

For nondimensional constant shear $U=y$ and streamwise wavenumber one, put $W=(D^2-1)\hat{v}$. The Orr-Sommerfeld equation gives $W^{\prime \prime }=[1+iRe(y-c)]W$. The displayed substitution transforms it into the Airy ordinary differential equation, because $(e^{i\pi /6}{)}^3=i$. The scale $R{e}^{-1/3}$ balances shear advection against viscous diffusion.

The equation $(D^2-1)\hat{v}=W$ is inverted by the displayed formula: differentiating twice produces the endpoint contribution $W(y)$, since $\sinh 0=0$ and $(\sinh {)}^{\prime }(0)=1$. If $W=a_1\mathrm{Ai}(Z)+a_2\mathrm{Bi}(Z)$, there are four independent constants $C_{+},C_{-},a_1,a_2$. This reconstructs the velocity from the Airy reduction of the Orr-Sommerfeld equation in constant shear; four no-slip and no-flux conditions determine the boundary eigenproblem.

The Squire equation governs wall-normal vorticity $\eta ={\partial }_{z}u-{\partial }_{x}w$ in the same normal mode convention as the Orr-Sommerfeld equation. It follows by taking ${\partial }_z$ of the streamwise momentum equation minus ${\partial }_x$ of the spanwise momentum equation. At a rigid no-slip wall, $\hat{\eta }=0$.

A Squire mode has zero wall-normal velocity and nonzero wall-normal vorticity, satisfying the homogeneous Squire equation. For real base velocity, positive finite Reynolds number, and zero boundary terms,for a nontrivial mode in a finite channel with no-slip walls. This follows by multiplying the homogeneous equation by $\overline{\eta }$, integrating by parts and taking the real part. Such modes are damped even when the coupled velocity-vorticity system can exhibit transient growth from non-normal modes.

For real base-flow velocity and positive Reynolds number, multiply the homogeneous Squire equation by the conjugate mode and integrate with zero endpoint values. Taking real parts gives the displayed negative temporal growth rate. This proves decay of homogeneous Squire modes; the coupled shear-flow system can still be unstable or show transient growth because wall-normal velocity can force the Squire variable.

An Orr-Sommerfeld mode has nonzero wall-normal velocity satisfying the Orr-Sommerfeld equation. Its wall-normal vorticity satisfies the accompanying forced Squire equation. The triangular coupling separates the velocity eigenproblem from the vorticity forcing; the vorticity need not vanish in three dimensions.

Inertial instability occurs when displaced fluid parcels amplify their displacement by conserving absolute momentum in a rotating flow. For a uniformly stratified basic state on an f-plane, a necessary condition is that the Ertel potential vorticity have the opposite sign to the Coriolis parameter.

For a parallel inviscid base flow $U(z)$ and disturbance streamfunction $\varphi (z)e^{i\alpha (x-ct)}$, the Rayleigh equation is

For $U=-1$ below $z=-1$, $U=z$ between $z=-1$ and $z=1$, and $U=1$ above $z=1$, decaying normal modes solve ${\varphi }^{\prime \prime }-k^2\varphi =0$ in each region. Pressure matching at a piecewise-linear shear interface gives continuity of $\varphi$ and $(U-c){\varphi }^{\prime }-U^{\prime }\varphi$. The two corner conditions yield the displayed dispersion relation. At $k=1/2$ the wave speeds are $c=\pm i/e$, so one mode grows exponentially with rate $1/(2e)$.

An axisymmetric Stokes streamfunction disturbance $\varphi (r)e^{ik(z-ct)}$ of an axial basic flow obeys the displayed ordinary differential equation. Regularity requires $\varphi =O(r^2)$ at the axis, and an impermeable wall requires $\varphi (a)=0$ for nonzero $k$. Multiplication by ${\varphi }^{\ast }/r$ and integration by parts shows that a growing real-wavenumber mode requires $(U^{\prime }/r{)}^{\prime }$ to change sign; if it vanishes identically, the positive energy integral rules out such a mode. The condition concerns this cylindrical vorticity gradient rather than $U^{\prime \prime }$ alone.

When $U^{\prime }/r$ is constant and $U\ne c$, the axisymmetric inviscid pipe stability equation reduces to ${\varphi }^{\prime \prime }-{\varphi }^{\prime }/r-k^2\varphi =0$. Axis regularity selects $r{I}_1(kr)$ and impermeability gives $I_1(ka)=0$. These imaginary-wavenumber disturbances are spatially evanescent potential fields, not bounded whole-pipe temporal normal modes; their radial equation supplies no eigenvalue condition on $c$.

Decaying normal modes of this continuous piecewise-linear profile satisfy ${\varphi }^{\prime \prime }-k^2\varphi =0$ between the corners. Pressure matching at a piecewise-linear shear interface yields an odd varicose mode of a planar jet with real speed $c=(1-e^{-2k})/(2k)$ and an even sinuous mode of a planar jet whose speeds satisfyIts reduced discriminant factors as $(2k-3-e^{-2k}-4e^{-k})(2k-3-e^{-2k}+4e^{-k})$. The second factor is positive for $k>0$, while the first is strictly increasing and crosses zero once. Thus exactly the long-wave even branch is exponentially unstable, for $0<k<k_c$, where $2k_c-3-e^{-2k_c}-4e^{-k_c}=0$.

For a normal mode of an inviscid parallel flow, the perturbation pressure is $p=-(U-c){\varphi }^{\prime }+U^{\prime }\varphi$. At an interface where $U$ is continuous but $U^{\prime }$ jumps, continuity of normal velocity and pressure gives continuity of $\varphi$ and the displayed flux. Consequently $(U-c)[{\varphi }^{\prime }]=[U^{\prime }]\varphi$. A sinuous mode of a planar jet may therefore have a cusp at a symmetric velocity-profile corner: even parity alone does not force its one-sided derivative to vanish.

For a smooth inviscid parallel shear profile with an exponentially growing mode, let $U_s$ be its velocity at an inflection point. The real and imaginary parts of the integrated Rayleigh equation for inviscid shear flow imply $\int U^{\prime \prime }(U-U_s)|\psi {|}^2/|U-c{|}^2=-\int (|{\psi }^{\prime }{|}^2+k^2|\psi {|}^2)<0$. Consequently $U^{\prime \prime }(U-U_s)<0$ somewhere in the channel. This strengthens Rayleigh's inflection-point theorem but is still only a necessary condition, not a sufficient instability criterion.

For inviscid plane Couette flow $U=z$ between impermeable walls at $z=\pm 1$, the nonzero real streamwise wavenumber $\alpha$ gives no nontrivial smooth discrete vertical-velocity modes. Instead, the vorticity variable $q=(D^2-{\alpha }^2)w$ obeys $({\partial }_t+i\alpha z)q=0$. Localized vorticity sheets give generalized eigenfunctions with real frequencies ${\omega }_{\xi }=\alpha \xi$. Their closure forms a neutral continuous spectrum.

For constant shear $U=\lambda y$, the wall-normal Laplacian of velocity is advected with phase $e^{-i\alpha \lambda yt}$. Inverting $D^2-k^2$, $k^2={\alpha }^2+{\beta }^2$, with its Dirichlet Green function gives the displayed initial-value solution. Impermeable inviscid walls impose zero normal velocity, not zero tangential velocity or no-slip wall vorticity.

In unbounded Couette flow with velocity $(y,0)$, linear vorticity is transported exactly as $\omega (x,y,t)={\omega }_0(x-yt,y)$. For a nonzero streamwise Fourier mode $\alpha$ and sufficiently localized smooth initial data, inversion by the one-dimensional modified Helmholtz Green function gives the displayed pointwise decay of the streamfunction and velocity, although vorticity does not decay in magnitude. The mechanism is phase mixing: shearing produces increasingly rapid transverse oscillations, whose elliptic velocity reconstruction suppresses them. Smoothness and decay of the initial streamfunction without derivative control do not alone imply these algebraic rates.

Let $G_a$ be the Dirichlet Green function of $D^2-a^2$ on $[-1,1]$, with $a=|\alpha |$:It is continuous and has derivative jump one, so $(D^2-a^2)G_a=\delta (z-\xi )$. The Dirac delta multiplication identity gives $(z-\xi )\delta (z-\xi )=0$, proving it is a generalized eigenfunction of the inviscid Couette continuous spectrum. Integrating $G_a(z,\xi )q_0(\xi )e^{-i\alpha \xi t}$ reconstructs the evolving vertical velocity from its initial vorticity.

Across a jump of $U^{\prime }$ in a constant-density parallel flow with continuous $U$, the normal mode vertical velocity and pressure are continuous. The linearized momentum equation gives $\hat{p}=\rho [(U-c){\hat{w}}^{\prime }-U^{\prime }\hat{w}]/(ik)$, yielding the displayed conditions. Equivalently, away from $U=c$, $(U-c)[{\hat{w}}^{\prime }]=[U^{\prime }]\hat{w}$. Derivative continuity is generally incorrect when the background vorticity jumps. This is the constant-density, continuous-velocity specialization of jump conditions for stratified inviscid shear flow.

An inviscid parallel shear flow has base velocity $\mathbf{u}=U(y){\mathbf{e}}_x$. Two-dimensional normal-mode perturbations satisfy the Rayleigh equation for inviscid shear flow.

A uniform-density inviscid shear flow with $U=U_{0}z/d$ for $|z|<d$ and $U=\pm U_0$ outside has two vorticity jumps but no walls. Decaying normal modes satisfy $w^{\prime \prime }-k^{2}w=0$ in each region. Continuity of $w$ and pressure matching at a piecewise-linear shear interface gives a two-by-two determinant whose diagonal entries are $1-2kd\pm 2kdc/U_0$ and off-diagonal entries $e^{-2kd}$. Its vanishing gives the displayed dispersion relation. The Kelvin-Helmholtz instability band is $0<kd<{\alpha }_s$, where $2{\alpha }_s-1=e^{-2{\alpha }_s}$ has its unique positive root above $1/2$, approximately $0.63923227$. Short waves decouple the two edges and remain neutral; long waves approach the vortex sheet instability with $c=\pm i{U}_0$.

This three-layer shear flow has $U=\Delta Uz/h$ and density ${\rho }_0$ between $z=\pm h/2$. Above and below, the velocity is the corresponding constant $\pm \Delta U/2$, and the mass density is ${\rho }_0\mp \Delta \rho /2$. The density drops and vorticity jumps coincide at both interfaces. Its hydrodynamic stability involves coupled gravity-vorticity interface waves, unlike a profile with uniform shear extending through all three layers.

For the bounded-shear three-layer profile with density contrasts minus one, zero and plus one, the localized-mode matching determinant is $[a(1-c{)}^2-(1-c)-J][a(1+c{)}^2-(1+c)-J]-b^2(1-c^2{)}^2=0$, where $a=\alpha [1+\coth (2\alpha )]$ and $b=\alpha \mathrm{csch}(2\alpha )$. Dividing by $a^2-b^2=2\alpha a$ gives a quartic in $c$ and a quadratic in $c^2$. Its constant term is $\{[2\alpha -(1+J){]}^2-e^{-4\alpha }(1+J{)}^2\}/(4{\alpha }^2)$, whose negative sign proves the unstable band of a three-layer stratified shear flow.

For the three-layer stratified piecewise-linear shear flow, set $\alpha =kh/2$. The dispersion relation is quadratic in squared nondimensional phase velocity. Its constant term is $[2\alpha -(1+J)(1+e^{-2\alpha })][2\alpha -(1+J)(1-e^{-2\alpha })]/(4{\alpha }^2)$. Throughout the displayed band this is negative, so one squared speed is negative and gives an exponentially growing normal mode.

The Hazel model is a smooth stratified shear flow whose base velocity has a hyperbolic-tangent profile and whose squared buoyancy frequency has a hyperbolic-secant-squared profile. The equal-width nondimensional example is $U=\tanh z$ and $N^2=J{\mathrm{sech}}^{2}z$. For $J\ge 0$, its gradient Richardson number is $J{\cosh }^{2}z$, with minimum $J$.

The zero-phase velocity ansatz $w=(\mathrm{sech}z{)}^k(\tanh z{)}^{1-k}$ solves the Taylor–Goldstein equation away from its critical level of an internal gravity wave exactly when $J=k(1-k)$. Let $y=\tanh z$ and $a=1-k$. Then $w^{\prime }/w=a/y-y$ and $w^{\prime \prime }/w=a(a-1)/y^2-a-1+2y^2$, so the differential-equation residual divided by $w$ is $[J-k(1-k)](y^{-2}-1)$. The curve has maximum $1/4$ at $k=1/2$. The endpoint $k=0$ does not decay at infinity; fractional powers at the critical level require a branch and regularity qualification.

For $0<k<1$, a neutral mode of the equal-width Hazel model behaves as $w\sim z^{1-k}$ near $z=0$. Its derivative diverges as $z^{-k}$, so it is not a classical continuously differentiable mode across the critical level of an internal gravity wave. The horizontal velocity is proportional to $w^{\prime }/k$; its local kinetic energy is finite only for $k<1/2$, because ${\int }_0^{ϵ}{z}^{-2k}dz$ converges precisely then. One possible neutral-mode convention is the boundary value from $U-c$ with $c_i\downarrow 0$, which fixes the phase of the power on the negative-$z$ side. The Miles–Howard theorem concerns growing modes with nonreal $c$ and does not rule out such singular neutral limits.

Let $\zeta =z/L\in [-1,1]$, $h=L_s/L\in (0,1)$, and normalize velocity by its maximum magnitude. The profile is $U=\zeta /h$ for $|\zeta |<h$ and $U=\mathrm{sgn}\zeta$ outside. Each region has exponential vertical-velocity modes; the two walls and vorticity-jump matching for an inviscid shear flow determine their coupling. With $\alpha =kL$, $q=\alpha h$, $X=\tanh q$, $Y=\tanh [\alpha (1-h)]$, $H=q(1+XY)$ and $J=q(X+Y)$, elimination gives $c^2=(H-Y)(J-XY)/(HJ)$, with $c$ normalized by the velocity scale.

The bounded piecewise-linear shear layer has exponentially growing modes exactly when $h<1/2$. Its squared phase speed tends to $2h-1$ as $kL\to 0$ and to $1$ as $kL\to \infty$. For $h<1/2$, sufficiently long nonzero waves are therefore unstable. Conversely $J-XY>0$, and for $h\ge 1/2$, $q\ge \alpha (1-h)>Y$, so $H-Y>0$ for every positive wavenumber and $c^2>0$. This proves the threshold without assuming monotonicity of the dispersion curve.

An inviscid parallel shear flow can possess an exponentially growing normal mode only if $U^{\prime \prime }$ changes sign somewhere in the flow.

If an inviscid shear-flow mode has complex wave speed $c$ and $U_{-}\le U\le U_{+}$, then

For a nonreal regular Rayleigh equation for inviscid shear flow mode vanishing at the walls, put $F=\hat{w}/(U-c)$. The equation becomes $[(U-c{)}^2{F}^{\prime }{]}^{\prime }-k^2(U-c{)}^{2}F=0$. Multiplication by $\overline{F}$ and integration by parts gives the displayed identity. With $Q=|F^{\prime }{|}^2+k^2|F{|}^2$, its imaginary and real parts give $\int UQ=c_r\int Q$ and $\int U^{2}Q=|c{|}^2\int Q$. Since $(U-U_{\min })(U_{\max }-U)\ge 0$, these imply Howard's semicircle theorem.

Waves travelling in opposite directions relative to their local background currents can share a laboratory phase velocity. Their interaction can then lock their phases and extract energy from the mean shear, producing growing disturbances. In a layered flow, the coupling is often exponentially weak when the interface separation exceeds the waves' decay length. A dispersion relation for two density interfaces in uniform shear gives a simple explicit example.

For stable density jumps in the three-layer stratified piecewise-linear shear flow, the upper interface's intrinsically left-going gravity-vorticity interface wave and the lower interface's intrinsically right-going wave have opposite laboratory phase velocities. They coincide at zero when $1+J=2\alpha$. At large wavenumber, interaction is of order $e^{-2\alpha }$ and the unstable band of a three-layer stratified shear flow has width $4\alpha e^{-2\alpha }+O(\alpha e^{-6\alpha })$ around this resonance.

Individually decaying nonorthogonal modes can interfere constructively and produce temporary energy growth. This mechanism is absent for an orthogonal eigenbasis of a normal evolution operator.

The Orr mechanism amplifies a leading two-dimensional shear disturbance as the background flow swings its wavefronts toward alignment. Conserved perturbation vorticity implies velocity amplitude proportional to the inverse wavevector magnitude. This produces algebraic or transient energy amplification without an exponentially growing normal mode. Viscosity limits the benefit of choosing an initially tightly wound disturbance.

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

A streamwise-independent cross-stream disturbance in an inviscid parallel shear flow has time-independent wall-normal velocity. It advects the base streamwise momentum, producing the displayed linear growth of streamwise velocity and quadratic growth of perturbation kinetic energy. This is algebraic amplification, not an exponentially unstable normal mode; it exemplifies transient growth from non-normal modes.

The energy-stability threshold is the parameter value below which a positive perturbation energy decreases monotonically for every disturbance. It can be lower than the linear instability threshold when non-normal transient growth is possible.

Symmetric instability occurs in rotating stratified flow when perturbations aligned with the basic current can extract energy from vertical shear. It is excluded when the Ertel potential vorticity has the stable sign; for constant gradients this gives a condition relating stratification, rotation, and horizontal buoyancy gradient.

The Eady model describes quasi-geostrophic baroclinic instability in a uniformly stratified layer with constant vertical shear, zero interior potential-vorticity gradient, and two rigid horizontal boundaries.

When both rigid boundaries of a constant-shear Eady model have the same small slope $\alpha$, define $a=N^2\alpha /(f_0\Lambda )-1$, $M=N|k|D/|f_0|$ and $C=c/(\Lambda D)$. Zero-interior-PV modes have vertical structure $A\cosh (Mz/D)+B\sinh (Mz/D)$. The boundary conditions $(z/D-C){\hat{\psi }}_{z/D}+a\hat{\psi }=0$ at both endpoints give the displayed dispersion relation. For $a\ge 0$ its square root is real and there is no exponential baroclinic instability. The same-sign boundary-pseudomomentum criterion is necessary, not sufficient.

In a constant-shear Eady model with slopes $\alpha$ and $\gamma$, set $A_0=\Lambda -N^2\alpha /f_0$ and $A_D=\Lambda -N^2\gamma /f_0$. The interior potential vorticity obeys $({\partial }_t+\Lambda z{\partial }_x)q^{\prime }=0$. In the zero-interior-PV subspace, integrating ${\psi }_x^{\prime }{q}^{\prime }$ by parts and using the two buoyancy boundary equations proves conservation of the displayed boundary quadratic form. Opposite signs of $A_0,A_D$ make it definite and exclude growing normal modes. For arbitrary nonzero interior PV, its derivative is $(N^2/f_0^3)\int {\psi }_x^{\prime }{q}^{\prime }dV$ and need not vanish. The restriction to zero interior PV is automatic for a growing smooth mode with complex phase speed.

For constant shear $U=\Lambda z$ between impermeable horizontal walls at $0,D$, with constant buoyancy frequency and no interior potential-vorticity gradient, define $m=NK/|f_0|$. The boundary conditions yieldThe bracket is negative for $0<mD<2.399\dots$, giving baroclinic instability for nonzero zonal wavenumber. For large separation the boundary waves decouple; the lower-wave speed tends to $\Lambda /m$, matching the dispersion relation of a semi-infinite Eady edge wave.

A uniformly stratified quasi-geostrophic approximation fluid with constant vertical shear occupies $z>0$, has a flat rigid lower wall, and has decaying perturbations aloft. On an f-plane, the basic interior potential vorticity is constant. Its smooth edge-wave modes have zero interior anomaly and propagate along the lower-boundary buoyancy gradient. The absence of an upper wall removes the pair of interacting edge waves that makes the finite-depth Eady model unstable.

A growing complex phase speed cannot equal the real basic velocity. The linear interior potential-vorticity equation then forces zero anomaly everywhere, decay fixes one exponential, and the wall fixes the real speed $c=\Lambda /m$. Thus there is no growing smooth decaying edge-wave eigenmode. This excludes exponential normal-mode instability without claiming that every initial disturbance lacks transient amplification or a singular neutral component.

A lower-boundary Eady edge wave in the semi-infinite Eady model has a smooth exponentially decaying vertical structure and a real frequency. It is supported by the horizontal boundary buoyancy gradient while its interior potential vorticity anomaly vanishes. Its speed equals the basic velocity at one penetration depth.

For basic shear $U(z)=\Lambda z$ with a resting lower boundary, the rigid-boundary buoyancy condition for quasi-geostrophic waves and the decaying vertical exponential giveThis is real. The edge wave propagates in the same laboratory direction as the basic flow above its resting lower wall. If the basic buoyancy is $M^{2}y+N^{2}z$, then $\Lambda =-M^2/f_0$.

The critical level of the regular edge wave $c=\Lambda /m$, $U=\Lambda z$, is $z_c=1/m$ for $\Lambda \ne 0$. At that height $U=c$, so mean advection and the phase term oppose in the wave frame. A literal condition $U=-c$ instead gives the unphysical height $-1/m$ below the lower wall.

A smooth zero-interior-potential vorticity mode hasDecay for $z\to \infty$ gives $\hat{\psi }=C{e}^{-mz}$. Its penetration depth $1/m$ increases with horizontal wavelength and decreases with buoyancy frequency. Singular neutral interior sheets are a distinct continuous-spectrum class.