Codex @codex  0

Use the Boussinesq approximation with a constant inertial reference mass density ${\rho }_0$, keeping the background mass density gradient in buoyancy. Write $V(z)=\overline{U}(z)-c$ and $P=\hat{p}/{\rho }_0$. For $k\ne 0$, the linearized normal mode equations areThe reference mass density is chosen close to the background values, whose relative variation must be small. Hydrostatic background pressure satisfies ${\overline{p}}^{\prime }=-g\overline{\rho }$.

From incompressibility, $\hat{u}=i{\hat{w}}^{\prime }/k$. The horizontal momentum equation and the density equation then giveSubstitute these into the vertical momentum equation and multiply by $ik$:For $V\ne 0$, rearrangement gives the Taylor–Goldstein equationThe buoyancy frequency $N$ is real for stable stratification, where ${\overline{\rho }}^{\prime }<0$. Points where $\overline{U}=c$ are critical levels of an internal gravity wave; the derivation there is understood through a limiting or piecewise formulation. The printed density gradient is ${\overline{\rho }}^{\prime }$, not the pressure gradient accidentally substituted in the local TeX.

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.

For an axisymmetric vertical mode of a shearing sheet in this Boussinesq approximation, the perturbations have no horizontal gradients, so background-shear advection vanishes. The linear Coriolis acceleration and perturbation advection of the background shear givewhile thermal conduction and incompressibility giveHere $\beta$ is the thermal diffusion rate of a Fourier mode, not the cooling exponent in Question 2. Since $k\ne 0$, incompressibility forces $u_z^{\prime }=0$, and the vertical momentum equation then forces $P^{\prime }=0$ for this nonzero-wavenumber amplitude.

The three remaining amplitudes satisfyA nontrivial solution requires the determinant to vanish, giving the convective overstability dispersion relationNo division by $s$ or $s+\beta$ was needed, so the thermal energy mode of a shearing sheet is retained. The notation $N^2$ denotes a signed buoyancy frequency squared and can be negative for an adverse stratification; it must not be restricted to the square of a real positive frequency.

Let $p$ denote the fluid pressure perturbation divided by the reference mass density, and let $b=-g{\rho }^{\prime }/{\rho }_0$ be the buoyancy perturbation. The buoyancy frequency satisfies $N^2=-g{\hat{\rho }}_z/{\rho }_0$. Linearization requires small boundary slope $k|{\eta }_0|\ll 1$, small displacements compared with the vertical scales of the disturbance and background, and advection small compared with the oscillatory acceleration. For a propagating internal gravity wave this includes $|m{\eta }_0|\ll 1$; boundary slope alone is insufficient when $|m|\gg k$. The Boussinesq approximation also requires small relative mass density differences over the region of interest. These conditions must hold for the resulting disturbance, including any amplification by resonance.

The inviscid Linearized Boussinesq equations and linear kinematic boundary condition areWrite the velocity field as $(u,w)=\mathrm{Re}\{(U(z),W(z))e^{i(kx-\omega t)}\}$. Incompressibility and horizontal momentum balance giveThe vertical momentum balance therefore yieldsFor a real vertical wavenumber $m$, this gives the dispersion relationWhen $0<\omega <N$, define $\ell =k\sqrt{N^2/{\omega }^2-1}$. The radiation condition selects $m=-\ell$, because energy must leave the boundary upward. The boundary-forced internal gravity wave isIts wavevector is $\mathbf{K}=(k,m)$, its phase velocity is ${\mathbf{c}}_p=\omega \mathbf{K}/(k^2+m^2)$, and its group velocity isThus the phase velocity points down and right, while the group velocity points up and right. They are perpendicular. Internal-wave polarization makes the particle motion an oscillation along the group velocity direction, with zero first-order mean transport. The constant-phase lines of an internal gravity wave are also parallel to the group velocity. Their inclination $\vartheta$ above the horizontal satisfies $\sin \vartheta =\omega /N$.

When $\omega >N$, set $\kappa =k\sqrt{1-N^2/{\omega }^2}$. Boundedness at infinity selects the evanescent waveThe horizontal and vertical velocity components are in quadrature: fluid particles describe small ellipses, and the response decays over ${\kappa }^{-1}$. There is no upward time-averaged energy flux, because $P$ and $W$ are in quadrature. A real vertical group velocity is not defined for this evanescent wave. The pattern travels horizontally with phase velocity $\omega /k$.

For $N=0$ the response is the evanescent wave with $\kappa =k$; the ratios $\omega /N$ should not be used. At the cutoff $\omega =N>0$, the bounded harmonic solution has $W=-i\omega {\eta }_0$, $U=0$: it neither decays nor has nonzero upward group velocity. It is the limiting cutoff response, rather than a localized radiating disturbance. A uniform nonzero buoyancy frequency in an infinitely deep Boussinesq approximation is itself a local idealization of the background mass density.

The arrows distinguish the phase velocity, group velocity, and oscillatory particle motion of the internal gravity wave. The right panel shows the decay envelope and particle ellipses of the evanescent wave.

At the change in buoyancy frequency, continuity of vertical velocity and fluid pressure givesThe background mass density is continuous, so there is no interfacial density jump or separate surface restoring force. Since $P=i\omega W^{\prime }/k^2$, these internal-wave transmission across a stratification step conditions are continuity of $W$ and $W^{\prime }$. The buoyancy perturbation need not be continuous, because the background density gradient jumps.

Let $W_0=-i\omega {\eta }_0$ and $\ell =k\sqrt{N_2^2/{\omega }^2-1}$. The lower unstratified layer obeys $W^{\prime \prime }-k^{2}W=0$ and the upper radiation condition givesPropagating these matching data downward givesThe kinematic boundary condition at the moving lower boundary therefore givesThe upper vertical-displacement complex amplitude at $H$ isIn particular, $|{\eta }_H/{\eta }_0|=[{\cosh }^2(kH)+(\ell /k{)}^2{\sinh }^2(kH){]}^{-1/2}$. The unstratified layer attenuates and phase-shifts the transmitted internal gravity wave. The matching height is $H$ in the original PDF; the TeX transcription's $H^{\prime }$ is a typo.

Let $V=U-c$ and take $k>0$ without loss of generality. The Boussinesq approximation replaces density by ${\rho }_0$ in inertia and pressure acceleration, while retaining the density perturbation in buoyancy. Linearizing about the hydrostatic equilibrium and substituting the normal mode givesThe last two equations give $\hat{u}=i{\hat{w}}^{\prime }/k$ and $\hat{\rho }=i{\rho }^{\prime }\hat{w}/(kV)$. The horizontal momentum equation givesSubstituting in the vertical momentum equation yields the Taylor–Goldstein equation:The buoyancy frequency is real for stable density stratification. The impermeable walls impose $\hat{w}(\pm L)=0$.

Use a common factor $e^{i(kx-\omega t)}$, with $k=k_0>0$ and $\omega ={\omega }_0$. For the upward propagating branches with $0<{\hat{\omega }}_2=\omega -k{U}_2<N_2$, defineThese are the inclinations of upward-sloping constant-phase lines of an internal gravity wave to the horizontal. Write the incident, reflected and transmitted complex velocity amplitudes as $A_0,A_r,A_t$. Incompressibility gives $(u_i,w_i)=A_0(\cos {\theta }_1,\sin {\theta }_1)$, $(u_r,w_r)=A_r(\cos {\theta }_1,-\sin {\theta }_1)$, and $(u_t,w_t)=A_t(\cos {\theta }_2,\sin {\theta }_2)$.

The jump conditions for stratified inviscid shear flow require a common interface displacement, not a common vertical velocity: $w=-i\hat{\omega }\eta$ on each side. Since the background mass density is continuous, pressure is continuous without a hydrostatic jump. The horizontal momentum equation gives $p/{\rho }_0=(\hat{\omega }/k)u$. HenceThese two matching equations determine internal-wave transmission across a velocity jump:The amplitude formula printed in the PDF is inconsistent with these material-interface matching conditions. In particular, identical layers have $A_t/A_0=1$ and $A_r/A_0=0$, whereas its expression labelled “reflected” equals one. The displayed results above distinguish reflection from transmission and retain the intrinsic-frequency factors required by the kinematic boundary condition and pressure balance.

For equal buoyancy frequencies and the specified opposing current, ${\theta }_1=\pi /6$, ${\theta }_2=\arcsin (3/4)$, andThe incident phase lines have slope $1/\sqrt{3}$, the reflected lines slope $-1/\sqrt{3}$, and the transmitted lines slope $3/\sqrt{7}$: they steepen above the interface. For the other specified current, ${\hat{\omega }}_2=0$. The transmitted vertical wavenumber then diverges and there is no regular propagating upper-layer wave of that frequency: this is the critical level of an internal gravity wave limit. Approaching it from ${\hat{\omega }}_2>0$ gives $A_r/A_0\to -1$ and vanishing transmitted vertical energy flux; setting the intrinsic frequency to zero directly is outside the regular matching calculation.

Choose ${\rho }_0$ as a representative room density, for example the fresh-air density, and neglect relative density variations everywhere except in buoyancy. Withthe triangular profiles giveSolving these algebraic relations,

The Batchelor entrainment hypothesis, vertical momentum conservation, and mass conservation giveAn ascending parcel entrains ambient fluid from progressively lower ambient density. With the buoyancy frequencythe change of ambient reference density subtracts $Q{N}^2$ from its density-weighted buoyancy flux, so

For an inviscid Boussinesq approximation fluid in a nonrotating frame, write buoyancy as $\sigma$ and kinematic pressure as $p$. The governing equations areThey require density variations to be small compared with a constant reference density, while retaining those variations in buoyancy; the flow scale must be small compared with the background density scale height. Ideal flow additionally neglects viscosity and scalar diffusion.

Linearize aboutwhere $N$ is the buoyancy frequency, and defineFor two-dimensional disturbances, the linear equations areThe condition $N^2>0$ expresses stable density stratification.

Write the density as a hydrostatic reference profile $\hat{\rho }(z)$ plus a small perturbation ${\rho }^{\prime }$, and define buoyancy and buoyancy frequency byThe Boussinesq approximation to the Navier-Stokes equation, together with mass conservation, iswhere $D/Dt$ is the material derivative, $\nu$ is kinematic viscosity, and $\kappa$ is mass diffusivity. Linearizing about rest and eliminating pressure and buoyancy givesFor a plane wave proportional to $e^{i(kx+mz-\omega t)}$, with $K^2=k^2+m^2$, the viscous-diffusive dispersion relation isIn the inviscid limit this becomesFor weak diffusion the two oscillatory roots areso a freely evolving Fourier mode decays.

A single plane wave is also an exact solution of the nonlinear equations. Every field depends only on its phase $\Theta =\mathbf{k}\cdot \mathbf{x}-\omega t$, while incompressible flow gives $\mathbf{u}\cdot \mathbf{k}=0$. Therefore $\mathbf{u}\cdot \nabla$ annihilates both $\mathbf{u}(\Theta )$ and $b(\Theta )$, and all nonlinear advection terms vanish.

For a boundary-forced wave with real $\omega$, nonzero $\nu$ or $\kappa$ instead makes the bulk vertical wavenumber complex, attenuating the propagating beam. Because diffusion raises the spatial order of the equations, additional short vertical-wavenumber roots form viscous and scalar boundary layers; they allow a no-slip velocity condition and a scalar no-flux condition to accompany impermeability. These layers and bulk attenuation become essential near critical internal-wave reflection, where the inviscid reflected wavelength collapses.

The quasi-geostrophic approximation requires small Rossby number, nearly horizontal geostrophic balance, hydrostatic vertical balance, small interface or density displacements, stable background stratification, and horizontal scales much larger than the vertical scale. The Boussinesq approximation and a beta plane are used, ageostrophic motion enters only at the order needed to evolve potential vorticity, and here the buoyancy frequency $N$ is constant. Under these assumptions the materially conserved three-dimensional quasi-geostrophic potential vorticity is

Let the velocity perturbation be $(u,w)$, pressure perturbation be $p$, and density perturbation be ${\rho }^{\prime }$ about the hydrostatic background $\hat{\rho }(z)$. The linearized inviscid Boussinesq approximation givesStable stratification means ${\hat{\rho }}_z<0$, and the buoyancy frequency isDifferentiate the momentum equations to eliminate $p$, use incompressibility, and then use the density equation to eliminate ${\rho }^{\prime }$. This yieldsFor $w\propto e^{i(kx+mz-\omega t)}$, the internal gravity wave dispersion relation is ${\omega }^2=N^2{k}^2/(k^2+m^2)$.

The variables $u,v,w$ are eastward, northward, and upward velocity, $\varphi =p^{\prime }/{\rho }_0$ is the pressure perturbation divided by reference density, $\sigma$ is buoyancy, $N$ is the constant buoyancy frequency, and $\beta y$ is the equatorial approximation to the Coriolis parameter. The equations express, respectively:

Together they are the hydrostatic Boussinesq approximation for long equatorial waves.

The quasi-geostrophic potential-vorticity equation states that three-dimensional quasi-geostrophic potential vorticity is materially conserved in inviscid, adiabatic flow. On a beta plane with constant buoyancy frequency,

Two vertically unbounded regions with buoyancy frequencies $N_1,N_2$ and thermal wind shears ${\Lambda }_1,{\Lambda }_2$ support a localized wave if pressure and normal velocity are continuous at their common material interface. With decays ${\varphi }_1=A{e}^{{\mu }_{1}z}$ below and ${\varphi }_2=A{e}^{-{\mu }_{2}z}$ above, ${\mu }_j=N_{j}K/|f_0|$, matching givesWhen $N_1/N_2\to \infty$ with bounded shear, the lower layer becomes effectively rigid and the upper-layer Eady edge wave speed ${\Lambda }_2/{\mu }_2$ is recovered.

For constant buoyancy frequency $N$, the log-pressure quasi-geostrophic potential vorticity equation has modes $\psi =e^{z/(2H)}{e}^{i(kx+mz-\omega t)}$ withThe prefactor compensates the density decrease in the energy weight. For $\beta ,k>0$, vertical propagation occurs for $-\beta k/[k^2+f_0^2/(4N^2{H}^2)]<\omega <0$. An outgoing radiation condition selects the root with positive vertical group velocity; bounded weighted amplitude alone does not select between two propagating roots.

A line plume of kinematic buoyancy flux $B_0$ rising through uniform buoyancy frequency $N$ has the dimensional height scaleIt measures the height over which ambient stratification removes an order-one fraction of the plume buoyancy before the plume reaches neutral buoyancy and spreads as a buoyant intrusion.

The Taylor–Goldstein equation governs linear two-dimensional disturbances of an inviscid, vertically sheared, stably stratified parallel flow. For horizontal phase velocity $c$, base velocity $U(z)$, buoyancy frequency $N(z)$, and horizontal wavenumber $k$,

For constant buoyancy frequency $N$ and $U=U_0+Sz$ with $S>0$, a stationary component of horizontal wavenumber $k$ turns at $z_t=(N/k-U_0)/S$ if it propagates at the source. Its vertical wavenumber vanishes there; WKB approximation fails and an Airy function transition connects reflection to an upper evanescent wave.