Past exam of the mathematics course of the University of Cambridge / 2015 / iii / Paper 74 / 1

Use the inviscid, nonrotating Boussinesq approximation, with a stable background mass density $\bar{\rho }(z)$ and reference density ${\rho }_{\ast }$. Write $B=-g{\rho }^{\prime }/{\rho }_{\ast }$ for the perturbation buoyancy and $\pi =p^{\prime }/{\rho }_{\ast }$ for the kinematic pressure. The buoyancy frequency is $N_0^2=-(g/{\rho }_{\ast }){\bar{\rho }}_z>0$. Dropping products of perturbations in the Boussinesq equations gives the Linearized Boussinesq equationsThe last equation expresses incompressible flow; the second follows by advecting the background mass density gradient. Neglecting rotation and viscosity is part of this internal gravity wave model.

For a plane internal gravity wave proportional to $\exp [i({\mathbf{k}}_h\cdot {\mathbf{x}}_h+mz-\omega t)]$, put $\kappa =|{\mathbf{k}}_h|>0$ and $K^2={\kappa }^2+m^2$. Eliminating the horizontal velocity, pressure and buoyancy from the linear equations givesFor example, the horizontal momentum and continuity equations give $\pi =-\omega mw/{\kappa }^2$; substituting $B=-i{N}_0^{2}w/\omega$ into vertical momentum yields the displayed dispersion relation. Thus an IGW has $0<|\omega |\le N_0$ in this model.

On the positive-frequency branch, the phase velocity normal to a constant-phase plane and the group velocity areConsequently $({\mathbf{k}}_h,m)\cdot {\mathbf{c}}_g=0$, and phase and group velocity are perpendicular. Equivalently, the dispersion relation is homogeneous of degree zero in the wave vector, so differentiating with respect to its scale proves the same orthogonality. Energy travels with the group velocity, along the phase planes. At the degenerate limit $m=0$, $\omega =N_0$ and the group velocity vanishes; the orthogonality statement then has this limiting interpretation.

For a two-dimensional internal gravity wave with fixed horizontal wave number $k\ne 0$, the vertical velocity amplitude satisfiesThis follows from the Linearized Boussinesq equations just as in part (a), allowing the background buoyancy frequency to vary. The WKB approximation gives locallyIt requires slowly varying stratification compared with the local vertical wavelength, in particular $|m_z|\ll m^2$, and a small-amplitude wave in the inviscid, nonrotating regime. Propagation requires $N>|\omega |$; where $N<|\omega |$, the local vertical wave number is imaginary and the solution is evanescent. Turning levels and abrupt changes need a separate connection calculation. Finite viscosity, background flow, rotation or nonlinear effects require modifications of this model.

The ray tracing equations use the local dispersion relation as a Hamiltonian:The stationary, horizontally uniform background conserves the frequency and horizontal wave number. In two dimensions the group velocity gives $dx/dz=-m/k$. For the specified buoyancy frequency and $\omega =1$, $m^2=256k^2{z}^2$. Hence for $z>0$ the two geometric internal-wave rays through the specified point areEach curve has both directed propagation senses, giving four local directed rays. Their formal intersections with $z=0$ are $x=-8$ and $x=8$.

At $z=0$, $N=\omega$ and $m=0$: this is a double-zero internal-wave turning level, since $m^2$ touches zero without changing sign. The ray becomes vertical geometrically, but its group velocity tends to zero. Taking $k>0$ without loss of generality, the downward branch hasThe formal travel time diverges logarithmically as $z\to 0$. Meanwhile $|m_z|/m^2=1/(16|k|z^2)$ grows, so the WKB approximation fails before that limit; its predicted $|m{|}^{-1/2}$ amplitude cannot be extrapolated to infinity. The spatial ray reaches a finite limiting position, but finite-wavelength behaviour near the zero must be found beyond WKB. There is no evanescent half-space supplied by this particular $N(z)$, so ordinary simple-turning-point reflection cannot be assumed solely from the local ray construction.

Use the spatial internal-wave ray tracing model, with $dx/dz=\pm 16z$, and reflect the rays so that the frequency and tangential wave number are preserved. At a vertical wall the horizontal energy direction reverses; at a horizontal wall the vertical energy direction reverses. On the sloping wall the reflection depends on internal-wave slope criticality. The three-reflection orbit below encounters that wall at a height above $3/16$, so the wall reverses the horizontal energy direction there.

First take the packet to leave the bottom upward and to the left, with $0<x_0<6$. The four segments of the loop areThe left-wall reflection is at $z=\sqrt{x_0/8}$ and the top reflection at $x=8-x_0$. The slope intersection solves $8z_s^2+3z_s+x_0-10=0$, giving the internal-wave ray return mapThus the next formal bottom reflection is at $(R(x_0),0)$. On $0<x_0<6$, the slope heights are between approximately $0.544$ and $0.946$, consistent with the assumed reflection type.

An internal-wave attractor is a stable periodic ray, so solve $R(x_{\ast })=x_{\ast }$. This gives $z_s=2/3$, and thereforeThe successive vertices areThis is clockwise circulation in the $(x,z)$ plane. Differentiating the internal-wave ray return map showsNeighbouring bottom intersections converge to the periodic orbit, geometrically focusing the energy beam.

For completeness, the opposite upward launch direction has the reversed three-reflection itinerary: slope, top, left. In the interval where this itinerary is valid, putThe three-reflection reverse map is $R^{-1}$, with domain $R(0)<x_0<R(6)$. Its fixed point is the same, but $R_{\mathrm{rev}}^{\prime }(x_{\ast })=41/23>1$. Counterclockwise packets defocus from that orbit. Except for the exactly periodic reverse ray, they eventually leave this itinerary. In the full geometric reflection construction, omitted-side or additional-slope reflections redirect generic rays into the clockwise attracting branch. Tangency and exactly critical reflections are singular exceptions requiring finite-wavelength or dissipative treatment.

The bottom reflection is a formal spatial-ray idealization: part (b) shows that the WKB approximation itself fails near $z=0$, so the diagram does not imply arrival and reflection there in finite WKB travel time. In a real fluid, viscous attenuation of an internal-wave beam, finite wavelength, internal-wave breaking and ensuing turbulence and mixing limit the beam width and energy density. Unlimited focusing is a property of the ideal geometric model, not a prediction for the resolved physical wave.