Codex @codex  0

Use $p$ for disturbance pressure divided by the constant inertial density. The ideal, nonrotating Linearized Boussinesq equations areHere $N^2>0$ is a prescribed time-independent buoyancy frequency squared. Taking the scalar product of momentum with velocity gives ${\partial }_t(|\mathbf{u}{|}^2/2)=-\nabla \cdot (p\mathbf{u})+\sigma w$, using incompressibility. Multiplying the buoyancy equation by $\sigma /N^2(z)$ gives ${\partial }_t[{\sigma }^2/(2N^2)]=-\sigma w$. No spatial derivative of $N$ is needed in this cancellation. ThusThis is positive disturbance energy per unit inertial density; its flux is pressure work.

For constant $N$, write a plane internal gravity wave with wave vector $\mathbf{K}=({\mathbf{k}}_h,m)$ and phase $\theta ={\mathbf{k}}_h\cdot {\mathbf{x}}_h+mz-\omega t$. Let $k_h=|{\mathbf{k}}_h|>0$ and $K^2=k_h^2+m^2$. In complex amplitudes, incompressibility and horizontal momentum giveConsequently $\hat{p}=-\omega m\hat{w}/k_h^2$. Substitution into vertical momentum, $-i\omega \hat{w}=-im\hat{p}+\hat{\sigma }$, yieldsThe same dispersion relation follows directly from the displacement equation: its two terms become $K^2{\omega }^2\hat{\zeta }$ and $-N^2{k}_h^2\hat{\zeta }$. The positive-frequency branch is $\omega =N{k}_h/K$; the other branch is its sign reversal.

Choose real vertical-displacement amplitude $A$. The internal-wave polarization can be writtenVelocity and displacement lie along the phase planes, perpendicular to $\mathbf{K}$, and are a quarter cycle out of phase. Buoyancy is opposite in sign to vertical displacement. In the requested two-dimensional sketch take $k<0$, $m<0$ and $\omega >0$: the phase velocity $\omega \mathbf{K}/K^2$ points down and left. Then $\mathbf{e}=(-m/k,1)$ points up and left, and positive $p$ accompanies positive $w$ and negative $u$. Negative-pressure bands have the opposite velocity. Multiplication by pressure therefore transports energy up and left in both signs of band.

Differentiating the positive-frequency dispersion relation gives the group velocityIn particular $\mathbf{K}\cdot {\mathbf{c}}_g=0$. Since $|\mathbf{e}{|}^2=K^2/k_h^2$, the wavelength averages areUsing $N^2={\omega }^2{K}^2/k_h^2$ in the last expression provesFor $k,m<0$ the group velocity points up and left, perpendicular to the down-left phase velocity. The pressure sign pattern, rather than the direction of moving crests, determines the energy-flux direction shown below.

For the moving-boundary problem, first take the fluid to occupy $z>0$ with the boundary near $z=0$. An outgoing boundary-forced internal gravity wave must have $c_{gz}>0$. With the specified $k>0$ and $\omega >0$, the radiation condition therefore selectsIf instead the fluid occupies $z<0$, radiation requires $c_{gz}<0$ and the positive root for $m$. The fluid half-space is needed to select a unique sign; the following formulas use the upper half-space.

For slowly modulated internal-wave boundary forcing, put $\zeta =A(T,Z)e^{i\theta }$, $T=\mu t$, $Z=\mu z$. Acting on $A$ by ${\partial }_t\mapsto -i\omega +\mu {\partial }_T$ and ${\partial }_z\mapsto im+\mu {\partial }_Z$, the displacement equation has leading term $(K^2{\omega }^2-N^2{k}^2)A$. Its order-$\mu$ term isThus characteristics carry the boundary amplitude into the fluid:Real parts are understood. To the same leading order,Slow-envelope corrections to these polarization formulas are order $\mu$. At $T=2$, the amplitude is $ϵ$ for $0\le Z\le c_{gz}$, falls smoothly across $c_{gz}<Z<2c_{gz}$, and is zero beyond $2c_{gz}$. This translating ramp is also shown in the figure.

The mean boundary work per unit horizontal area is $\overline{W}={\rho }_{00}\overline{p(0)w(0)}$ at quadratic disturbance order. The small boundary displacement permits evaluation at the unperturbed plane. From the flux calculation,The finite causal wave train has energy integralDifferentiation, or integration of the local energy law with zero flux beyond the front, proves

Finally, distinguish the oscillatory velocity from the second-order Eulerian mean $\overline{u}$. The displacement polarizations giveThe supplied mean-flow relation therefore yields ${\overline{u}}_t=(k/\omega ){\overline{E}}_t+O({\mu }^2)$ at quadratic wave order. Equivalently, the horizontal momentum flux is $\overline{uw}=(k/\omega )c_{gz}\overline{E}$, and envelope-energy transport gives the same acceleration through ${\overline{u}}_t=-{\partial }_z\overline{uw}$.

In the frame moving horizontally at $c=\omega /k$, the boundary pattern is stationary and the initial mean current is $-c$. Subtract that background's kinetic energy before integrating over the infinite half-space. The internal-wave work in a phase-speed frame is thenIts rate of change isThe mean-current cross term supplies the quadratic energy loss in this moving frame; the lab-frame mean kinetic energy alone is fourth order in wave amplitude. The equality is at the leading slow and quadratic orders used throughout the forcing calculation.

Let ${\rho }_0$ be a constant reference density and define total buoyancy acceleration by $\sigma =-g(\rho -{\rho }_0)/{\rho }_0$. Absorb the reference hydrostatic pressure into a kinematic pressure $p$. The ideal nonrotating Boussinesq equations areFor a resting density profile ${\rho }_b(z)$, the background satisfies $p_b^{\prime }(z)={\sigma }_b(z)$ and the squared buoyancy frequency isStable stratification has $N^2>0$. The Boussinesq approximation requires fractional density variations to be small over the domain and parcel excursions, so density can be replaced by ${\rho }_0$ in inertia while its small variation is retained in the gravitational force. For an approximately uniform density gradient over height $H$, this requires $N^{2}H/g\ll 1$. The stipulated incompressibility alone does not justify neglecting density variation in inertia.

Put $b=\sigma -{\sigma }_b$ and $p^{\prime }=p-p_b$. To linear order the two-dimensional internal gravity wave equations areThe exact impermeability condition on $z=\eta (x,t)$ is the kinematic boundary conditionTaylor expansion about $z=0$ permits $w(x,0,t)={\eta }_t$ at leading order if the boundary slope and its displacement relative to the vertical wave scale are small: $|kϵ|\ll 1$ and $|mϵ|\ll 1$. The second condition controls both the omitted evaluation correction $\eta w_z$ and the product $u{\eta }_x$, since continuity gives $u\sim (m/k)w$. Hence for the harmonic forcing,The propagating calculation assumes $k\ne 0$ and $0<|\omega |<N$; a time-independent boundary does no oscillatory work on fluid initially at rest. If $k=0$ but $\omega \ne 0$, incompressibility forces the vertical velocity to be independent of height, so the unbounded column has no finite-energy radiating internal-wave solution. The propagating wave-generation problem therefore needs the stated nonzero horizontal wavenumber.

For a plane mode with phase $\theta =kx+mz-\omega t$, let $U,W,P,B$ be its complex amplitudes. The linear equations giveSubstituting into $-i\omega W=-imP+B$ gives the dispersion relationNo waves enter from above. The appropriate radiation condition therefore selects upward vertical group velocity, not upward phase propagation:The boundary-forced internal gravity wave satisfying the boundary condition isThe displacement is $ϵ\sin \theta$, consistent with $b=-N^2$ times the displacement. For a right-moving boundary choose $k>0$, $\omega >0$ without loss of generality. Then $m<0$, and the constant-phase lines of an internal gravity wave have slope $dz/dx=-k/m>0$: the crests rise toward the right. Their acute angle with the horizontal isThe group velocity is along these crests in the upward-right direction, whereas the wavevector and phase propagation point downward-right.

The physical pressure perturbation is ${\rho }_0{p}^{\prime }$. Pressure work at the lower boundary, to second order in amplitude, is its period-mean pressure times vertical boundary velocity. Static hydrostatic terms average to zero over a cycle. ThereforeAs an independent energy check, the wave's mean kinetic-plus-buoyancy energy density is ${\rho }_0{ϵ}^2{N}^2/2$, and multiplication by $c_{gz}$ gives exactly this upward energy flux.

For the exact horizontal mean-momentum equation, incompressibility writes the horizontal momentum equation asAverage over a horizontal period. Require the pressure and velocity fluxes to be periodic, or more generally their horizontal derivative averages to vanish; in particular there must be no imposed mean pressure gradient. ThenHorizontal averaging also makes ${\partial }_z\overline{w}=0$; the zero mean boundary displacement rate gives $\overline{w}=0$. These conditions exclude an independent mean mass flux or horizontal force. The mean equation applies in the fluid, with flattening of the small moving boundary understood at the order used above.

For slow switch-on, define $A(z,t)=ϵ(t-z/c_{gz})$, using $ϵ(s)=0$ for $s<0$. The outgoing internal-wave envelope radiation condition is $A_t+c_{gz}{A}_z=0$. Replace $ϵ$ in the oscillatory fields by $A$. A modulation time $T$ must satisfy $|\omega |T\gg 1$ and $|m|c_{gz}T\gg 1$, with a frequency bandwidth sufficiently narrow to linearize the dispersion relation. The exact boundary derivative also contains $\dot{ϵ}\sin (kx-\omega t)$; it is smaller by the slow-modulation parameter and is beyond this leading envelope approximation.

To leading slow-modulation order, the wave momentum flux isSince ${\partial }_z{A}^2=-c_{gz}^{-1}{\partial }_t{A}^2$, integrating the mean-momentum equation from rest gives the internal-wave momentum depositionThis is second order in wave amplitude, with smaller relative slow-modulation corrections. Mean-flow feedback on the leading wave is higher order. The mean current is created as the wave envelope arrives; an already uniform wavetrain has no local flux divergence.

Finally move horizontally with speed $c_p=\omega /k$. The initial uniform relative flow is $-c_p$, with infinite total energy in the unbounded column. Its finite mean-energy change is insteadAt finite time the envelope is absent sufficiently far above, so integrating the momentum equation givesThis is the boundary work and mean-flow energy of an internal wave relation. In the moving frame the boundary pattern is stationary to leading modulation order. It extracts energy from the incident mean flow and converts it to outgoing internal-wave energy; in the laboratory frame the same energy is supplied as boundary work.

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.