Past exam of the mathematics course of the University of Cambridge / 2024 / iii / Paper 345 / 2

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.

Let $\theta$ be the angle made by the group velocity ray with the horizontal. The internal gravity wave dispersion relation givesEach sawtooth face changes height by $2h_0$ over horizontal distance ${\lambda }_T/2$, so its slope magnitude isThe internal-wave slope criticality criterion is $s_T<\tan \theta$. Hence reflection is subcritical forIn the corresponding internal-wave ray tracing sketch, every incident ray meets one planar face and leaves it into the fluid at the same angle $\theta$ to the horizontal. The reflected ray is steeper than either face, so it clears the sawtooth rather than running into an adjacent corner.

For a face of signed slope $s$, conservation of frequency and tangential wavenumber gives, on the branch relevant to an incident downward-right ray,At $s\mu =1$, the reflected wavenumber diverges and the reflected group velocity becomes tangent to the face. On a supercritical face the denominator changes sign: horizontal propagation reverses and rays from the two faces are directed towards a sawtooth corner. Successive reflections therefore focus energy and shorten the wavelength.

The inviscid ray pattern cannot persist indefinitely. Near-critical focusing amplifies gradients until kinematic viscosity, mass diffusivity, nonlinear wave steepening, and wave breaking matter; a real corner is also rounded on some finite scale. These effects replace the singular ray construction by dissipative boundary layers, mixing, and a finite-width reflected beam.

The maximum boundary slope is $s_{\max }=h_0{k}_T$, so subcritical internal-wave reflection requiresor equivalently

Take $k,m>0$ without loss of generality, putand write the incident field as the imaginary part of ${\tilde{w}}_i{e}^{i(kx+mz-\omega t)}$. The kinematic boundary condition on $z=h(x)$ isFor a flat boundary the reflected wave is $-{\tilde{w}}_i\sin (kx-mz-\omega t)$. Expanding the boundary condition in a Taylor expansion about $z=0$ creates the topographic sidebands of an internal gravity wave. With upward-radiating vertical wavenumber $-\mu |p_j|$, their complex amplitudes through second order areThus the general compact result isFor the convenient nondegenerate case $k>2k_T$, all displayed sideband wavenumbers are positive. Defining ${\Theta }_j=p_{j}x-\mu p_{j}z-\omega t$, the same answer is the explicitly real formulaValidity requires a linear incident wave, an inviscid uniformly stratified bulk, an outgoing-radiation condition, strict separation from critical slopes, and small boundary excursions for every retained mode, in particular $h_0{k}_T\ll 1$ and $\mu h_0{\max }_{|j|\le 2}|p_j|\ll 1$. A vanishing $p_j$ is a degenerate zero-horizontal-wavenumber case and must be treated by taking the corresponding zero-amplitude limit rather than dividing by $p_j$.

If $h_0{k}_T>\tan \theta$, each period contains two critical points satisfying $|h_x|=\tan \theta$, with supercritical intervals around the steepest parts. Internal-wave ray tracing sends neighbouring reflected rays towards caustics attached to those critical points; rays can reverse horizontal direction in the supercritical intervals and intersect rays reflected elsewhere on the sinusoid. At equality the inviscid reflected wavelength collapses locally. The physical pattern is therefore a set of intense finite-width beams and mixing regions once viscosity, scalar diffusion, and wave breaking regularize the ray caustics.