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.

Assume $U>0$, $k=\alpha >0$, constant $N$, negligible rotation, dissipation and reflection from infinity, small $k{\eta }_0$ and $m{\eta }_0$, and no wave overturning. The terrain is stationary, so laboratory frequency is zero and intrinsic frequency is $\hat{\omega }=-kU$. Propagation requires $kU<N$. PutThe kinematic boundary condition at the linearized bed is $w(x,0)=U{\eta }_x=-kU{\eta }_0\sin (kx)$. The upward radiation condition selects the positive $m$ branch because $c_{gz}=-\hat{\omega }m/(k^2+m^2)>0$. ThusThe internal-wave polarization and every other perturbation follow from this vertical velocity and the linear equations.

In the frame moving with the fluid, $x^{\prime }=x-Ut$, the phase is $k{x}^{\prime }+mz+kUt$, andThe intrinsic energy ray slopes upward and upstream, $dz/d{x}^{\prime }=-k/m$, along a constant-phase line of an internal gravity wave. At fixed intrinsic frequency the spatial wave equation is hyperbolic, with the two characteristic curve line slopes $dz/d{x}^{\prime }=\pm k/m$; the radiation condition selects the upward ray of the appropriate horizontal propagation branch.

In the terrain frame, phase lines remain fixed and the pattern phase velocity is zero. A transported packet hasIts energy path slopes upward and downstream, $dz/dx=m/k$, and need not coincide with a stationary phase line. One should not add $(U,0)$ to the chosen normal intrinsic phase velocity and call that the normal ground-frame phase velocity: tangential motion along a phase line is arbitrary; zero laboratory frequency makes the entire crest pattern stationary. The sketch distinguishes stationary crests, intrinsic characteristics and transported energy.

If $kU>N$, $m$ is imaginary and the bounded field is an evanescent wave; at $kU=N$ the vertical group velocity vanishes and an ordinary upward radiation interpretation degenerates.

For $\omega >N$, the dispersion relation gives $k^2=-m^2{\omega }^2/({\omega }^2-N^2)$. Put $\kappa =m\omega /\sqrt{{\omega }^2-N^2}>0$. Boundedness as $x\to +\infty$ selects $k=i\kappa$, rather than the exponentially growing root. The internal-wave polarization relations now yieldAgain take real parts and set $v=0$. This is an evanescent wave, decaying over distance ${\kappa }^{-1}$; there is no propagating far-field internal gravity wave at this frequency. The pressure amplitude $P=-i\omega u/\kappa$ is in quadrature with the horizontal velocity, so the time-averaged horizontal pressure-work flux vanishes. The equality $\omega =N$ is a singular limiting case and is outside the two frequency ranges considered here.

For $0<\omega <N$, put $\alpha =\omega /\sqrt{N^2-{\omega }^2}$. The dispersion relation permits $k=\pm \alpha m$. On the positive-frequency branch, for $k>0$,The radiation condition therefore selects $k=\alpha m$, since energy must travel into the fluid. Substitution into the internal-wave polarization relations gives the complex fieldsThe physical fields are their real parts; $v=0$. The vertical component of group velocity is negative, so this positive-$m$ plane internal gravity wave carries energy down and right, although its phase propagates in a different direction.

In the Boussinesq approximation, density departures from a constant reference density ${\rho }_0$ are small enough to neglect in inertia and mass continuity, but their gravitational effect is retained. Write the total buoyancy acceleration as $b=-g(\rho -{\rho }_0)/{\rho }_0$, choose a hydrostatic background $\bar{b}(z)$, and set $\sigma =b-\bar{b}$. Thus $\sigma$ is the buoyancy perturbation, positive for a parcel lighter than the background at its present height. The squared buoyancy frequency isStable stratification has $N^2>0$. A small vertical displacement $\xi$ of a parcel conserving its buoyancy gives $\sigma =-N^2\xi$, explaining the restoring acceleration and the name of $N$.

Let $P=p^{\prime }/{\rho }_0$ be pressure with the background hydrostatic pressure removed, $\mathbf{u}=(u,v,w)$, and $D/Dt={\partial }_t+\mathbf{u}\cdot \nabla$ the material derivative. The ideal, nonrotating Boussinesq equations areThe last equation is equivalent to $Db/Dt=0$. There is neither viscosity nor buoyancy diffusion. The reference density occurs in the pressure normalization, while the small density anomaly occurs only through buoyancy.

For the two-dimensional disturbance, linearization about rest givesA plane internal gravity wave proportional to $e^{i(kx+mz-\omega t)}$ therefore satisfiesThese internal-wave polarization relations and the dispersion relation determine the individual responses in parts i and ii. A radiating response has no incoming internal gravity wave from $x=+\infty$; an evanescent response is bounded there. The harmonic solutions describe the response after excluding arbitrary initial transients.

For the localized oscillatory boundary, define $\alpha =\omega /\sqrt{N^2-{\omega }^2}$ and $I(z)=ϵ{\int }_0^{\infty }{\hat{u}}_0(m)e^{imz}dm=f(z)+ig(z)$. Both signs of vertical wavenumber must radiate into $x>0$. Their horizontal wavenumber is $k=\alpha |m|$, so the positive-$m$ contribution involves $Z_1=z+\alpha x$ and the negative-$m$ contribution involves $Z_2=z-\alpha x$. Consequently the complex horizontal velocity isTaking its real part provesThe positive-$m$ internal-wave polarization has $w=-\alpha u$ and the negative-$m$ one has $w=+\alpha u$. HenceIn particular, $2f(z)=ϵu_0(z)$, so the boundary condition is recovered exactly. This is two-beam radiation from a localized oscillating boundary.

Take a smooth localized profile, for example $u_0(z)=e^{-z^2/H^2}$. Its effective support is a neighbourhood of width a few $H$; increasing that neighbourhood makes the tails arbitrarily small. The boundary displacement is $ϵu_0(z)\sin \omega t/\omega$ to first order, so at $t=0$ it is flat and moving at maximum speed. At that instant,There are two narrow oblique beams centered on $z=-\alpha x$ and $z=+\alpha x$. Away from their overlap, the velocity is parallel to the relevant beam, in the directions $(1,-\alpha )$ and $(1,+\alpha )$. The figure uses $\alpha =3/4$ and measures velocities in units of $ϵ$.

At other phases the $g$ terms contribute. For the positive-frequency convention above, $g$ is the Hilbert transform of $f$; it generally has longer tails than $f$. Thus even when the boundary is nearly stationary, the interior need not be stationary or confined to the bright cores shown here. Half a period later all velocities reverse. The beam directions remain the same, but the instantaneous velocity profile and its quadrature component change.

A short pulse has a broad frequency spectrum. Frequencies below $N$ produce outgoing internal gravity wave packets with frequency-dependent ray slopes $\pm \omega /\sqrt{N^2-{\omega }^2}$; frequencies above $N$ contribute an evanescent near-boundary response. Moreover the group velocity depends on wavelength as well as direction. The radiated packet therefore spreads into a dispersive fan or wake, with different frequencies traveling at different angles and different scales separating along those directions. After the forcing ends the propagating packet travels away, rather than maintaining the two permanently driven monochromatic beams.

Choose the positive-frequency branch with $k>0$, let $K=|\mathbf{k}|$, and write $\mathbf{k}=K(\cos \theta ,\sin \theta )$. The internal gravity wave dispersion relation is $\omega =Nk/K=N\cos \theta$. The phase velocity and group velocity areThey are perpendicular: the phase advances normal to the crests while energy travels along them. For a general signed $k$, the positive frequency is $N|k|/K$ and the horizontal component of group velocity changes sign with $k$.

If $\omega >N$, there is no real wave vector satisfying this relation. For prescribed real horizontal wavenumber, $m^2=k^2(N^2/{\omega }^2-1)<0$, so a bounded boundary response is an evanescent wave, rather than a freely propagating plane wave.

For a single incompressible plane wave, internal-wave polarization gives ${\mathbf{u}}_1\cdot \mathbf{k}=0$. Its oscillatory displacement ${\mathbf{\xi }}_1$ is parallel to ${\mathbf{u}}_1$ and therefore also perpendicular to $\mathbf{k}$. Every spatial derivative of the plane-wave field is proportional to a component of $\mathbf{k}$, soThere is no leading Stokes drift for that single wave. The fluid parcels oscillate along straight lines within phase planes and never sample a spatial change of the wave along their displacement. This conclusion does not prohibit cross terms between different wave vectors.

Assume $\kappa ,\omega >0$ and impose upward group velocity in $z>0$. The two positive-frequency components then have wave vectors ${\mathbf{k}}_1=(\kappa ,m)$ and ${\mathbf{k}}_2=(-\kappa ,m)$ with $m<0$. Equal lengths are automatic, and orthogonality gives $m^2={\kappa }^2$. HenceThe given orthogonal-wave condition is therefore possible only at this forcing frequency. The upward radiation condition excludes the $m=+\kappa$ branch.

Use superposition with vertical amplitude $W$ in each component and internal-wave polarization:The two sets of crests have slopes $+1$ and $-1$; their energy propagates upward and outward at angles of $45^{\circ }$, while their phases travel downward. Their interference makes fixed horizontal bands of vertical-velocity nodes, at $\cos \kappa x=0$, and horizontal-velocity nodes, at $\sin \kappa x=0$.

Put $A=\kappa x$, $B=\kappa z+\omega t$. The first-order displacement relative to mean parcel labels isThe horizontal displacement correction vanishes identically, whileConsequently the formal cross-wave Stokes drift isIt has alternating downward and upward columns, no depth decay in this ideal infinite plane-wave field, and zero horizontal-period-averaged vertical transport. Each wave alone has zero drift; the interference terms create this result.

This calculation is a kinematic second-order consequence of the linear field, not a complete prediction of mean transport. In particular, an impermeable oscillating lower boundary cannot have that nonzero mean normal parcel velocity. Its first-order displacement is ${\zeta }_b=(2W/\omega )\cos (\kappa x)\sin \omega t$. Expanding the exact kinematic boundary condition on $z={\zeta }_b$ givesThis Eulerian mean velocity cancels the computed normal Stokes drift at the material boundary. Away from it, the second-order mean-flow and buoyancy equations must also be satisfied. Thus the alternating drift is a valid formal cross-wave effect, but treating it alone as persistent transport with zero Eulerian mean would be physically incomplete. A prescribed permeable boundary would require its own second-order flux specification instead.

Use the nonrotating linear Boussinesq approximation with reference density ${\rho }_0$ and $N^2=-(g/{\rho }_0)d\bar{\rho }/dz>0$:Let $K_h^2=k^2+l^2>0$, $K^2=K_h^2+m^2$ and $E=e^{i(kx+ly+mz-\omega t)}$. The horizontal momentum equations give $u^{\prime }=k{p}^{\prime }/({\rho }_0\omega )$ and $v^{\prime }=l{p}^{\prime }/({\rho }_0\omega )$. Continuity requires $k{u}^{\prime }+l{v}^{\prime }+m{w}^{\prime }=0$. Thus the internal-wave polarization isPhysical perturbations are the real parts. The density is in temporal quadrature with the vertical velocity, because buoyancy responds to the vertical fluid displacement. Substitution into the vertical momentum equation determines the internal gravity wave dispersion relation,A nontrivial propagating wave requires a nonzero horizontal wave vector and a nonzero frequency. The formulas are not to be divided by $K_h^2=0$: incompressibility excludes a nonzero oscillatory vertical velocity for a purely vertical wave vector.

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.