Past exam of the mathematics course of the University of Cambridge / 2008 / iii / Paper 81 / 1 / Solution

In the Boussinesq approximation, replace density by a constant reference density ${\rho }_0$ in inertia and incompressible flow, retaining the small density departure only in the gravitational force. Write $\rho =\bar{\rho }(z)+{\rho }^{\prime }$ and subtract the resting hydrostatic pressure. The upward buoyancy perturbation acceleration and squared buoyancy frequency areThus stable density stratification has $N^2>0$. Let $p$ henceforth denote the pressure anomaly divided by ${\rho }_0$. The ideal, nonrotating Boussinesq equations areThe buoyancy equation follows by materially conserving the total buoyancy $b_0(z)+\sigma$, where $b_0^{\prime }=N^2$; it does not assume constant $N$.

For an $x$-periodic flow of period $P$, define $\bar{a}=P^{-1}{\int }_0^{P}adx$. An infinite-domain average with vanishing end contributions works as well. The necessary pressure condition is no mean zonal pressure gradient, $\overline{p_x}=0$: periodic pressure satisfies it, whereas a pressure containing a term proportional to $x$ does not. Use incompressible flow to write the zonal momentum equation in conservative form,Averaging its zonal derivative terms yields the exact Reynolds stress relation

For two-dimensional small disturbances about rest, the Linearized Boussinesq equations are $u_t=-p_x$, $w_t=-p_z+\sigma$, ${\sigma }_t=-N^{2}w$, and $u_x+w_z=0$. Differentiate the two momentum equations to eliminate pressure: $(u_z-w_x{)}_t=-{\sigma }_x$. Differentiating in $x$ and using $u_{xz}=-w_{zz}$ gives ${\nabla }^2{w}_t={\sigma }_{xx}$. One further time derivative givesThe equation remains valid for $N=N(z)$. No vertical derivative of $N$ was taken: ${\partial }_x^2(N^{2}w)=N^2{w}_{xx}$.

For constant $N$ and $k\ne 0$, a plane internal gravity wave hasBoth signs give the prescribed lower-boundary velocity. The radiation condition excludes energy arriving from infinity. At positive frequency, the vertical group velocity isFor a source below the fluid, energy must propagate upward, so the outgoing solution has $m<0$, even though its vertical phase velocity $\omega /m$ is downward.

To obtain this selection causally, let $T=\mu t$ and $Z=\mu z$, and replace the constant amplitude by $a(Z,T)$. The method of multiple scales treats $(z,t,Z,T)$ as independent during differentiation: the carrier derivatives hold $Z,T$ fixed, while $a_T$ holds $Z$ fixed and $a_Z$ holds $T$ fixed. Only after differentiating do we restrict to $T=\mu t$, $Z=\mu z$. Acting on the complex wave, the operators becomeWith $K^2=k^2+m^2$, the wave equation has residualThe zeroth-order term is the dispersion relation; the first-order term vanishes precisely whenIf $a_0(T)$ is the switched-on boundary amplitude, with $a_0(T)=0$ for $T<0$, the outgoing envelope is $a(Z,T)=a_0(T-Z/b)$ for $b>0$. Its characteristics carry the boundary signal into initially undisturbed fluid. If $b<0$, characteristics instead arrive from the upper half-space; zero initial disturbances and no incoming signal cannot establish the proposed boundary-forced wavetrain. This is the internal-wave envelope radiation condition.

The leading incompressible flow polarization gives $u^{\prime }=-(m/k)w^{\prime }$. Hence, to second order in amplitude and leading order in the slow modulation,There is no meridional transport here. Envelope transport gives ${\partial }_t|a{|}^2=-b{\partial }_z|a{|}^2$, so an initially vanishing mean develops asThese are the leading expressions for internal-wave momentum deposition; higher slow-modulation corrections to the stress are beyond the displayed approximation. In particular, a spatially uniform fully established wavetrain has zero stress divergence, but the front has already produced its mean flow.

The mean wave energy density per unit mass follows from ${\sigma }^{\prime }=-i{N}^2{w}^{\prime }/\omega$:Thus $\bar{u}=(k/\omega )\mathcal{E}$. In the original frame the mean kinetic energy is ${\bar{u}}^2/2=O(a^4)$. In a frame translating at $c=\omega /k$, the undisturbed velocity is $-c$, and the change in mean kinetic energy is $-c\bar{u}+O(a^4)=O(a^2)$. Indeed $\mathcal{E}-c\bar{u}=0$ at this order behind the front. Mean-flow energy is therefore essential at the same order as wave energy in the translating frame. This follows also from the frame transformation of energy, $E\mapsto E-cP$ up to the fixed background constant, and is why conservation of momentum cannot be omitted from the moving-frame conservation of energy budget.