Paper 4, Section II, B

Fluid Dynamics
Part IB, 2014

Consider a layer of fluid of constant density ρ\rho and equilibrium depth h0h_{0} in a rotating frame of reference, rotating at constant angular velocity Ω\Omega about the vertical zz-axis. The equations of motion are

utfv=1ρpxvt+fu=1ρpy0=pzρg\begin{aligned} \frac{\partial u}{\partial t}-f v &=-\frac{1}{\rho} \frac{\partial p}{\partial x} \\ \frac{\partial v}{\partial t}+f u &=-\frac{1}{\rho} \frac{\partial p}{\partial y} \\ 0 &=-\frac{\partial p}{\partial z}-\rho g \end{aligned}

where pp is the fluid pressure, uu and vv are the fluid velocities in the xx-direction and yy direction respectively, f=2Ωf=2 \Omega, and gg is the constant acceleration due to gravity. You may also assume that the horizontal extent of the layer is sufficiently large so that the layer may be considered to be shallow, such that vertical velocities may be neglected.

By considering mass conservation, show that the depth h(x,y,t)h(x, y, t) of the layer satisfies

ht+x(hu)+y(hv)=0.\frac{\partial h}{\partial t}+\frac{\partial}{\partial x}(h u)+\frac{\partial}{\partial y}(h v)=0 .

Now assume that h=h0+η(x,y,t)h=h_{0}+\eta(x, y, t), where ηh0|\eta| \ll h_{0}. Show that the (linearised) potential vorticity Q=Qz^\mathbf{Q}=Q \hat{\mathbf{z}}, defined by

Q=ζηfh0, where ζ=vxuyQ=\zeta-\eta \frac{f}{h_{0}}, \text { where } \zeta=\frac{\partial v}{\partial x}-\frac{\partial u}{\partial y}

and z^\hat{\mathbf{z}} is the unit vector in the vertical zz-direction, is a constant in time, i.e. Q=Q0(x,y)Q=Q_{0}(x, y).

When Q0=0Q_{0}=0 everywhere, establish that the surface perturbation η\eta satisfies

2ηt2gh0(2ηx2+2ηy2)+f2η=0\frac{\partial^{2} \eta}{\partial t^{2}}-g h_{0}\left(\frac{\partial^{2} \eta}{\partial x^{2}}+\frac{\partial^{2} \eta}{\partial y^{2}}\right)+f^{2} \eta=0

and show that this equation has wave-like solutions η=η0cos[k(xct)]\eta=\eta_{0} \cos [k(x-c t)] when cc and kk are related through a dispersion relation to be determined. Show that, to leading order, the trajectories of fluid particles for these waves are ellipses. Assuming that η0>0,k>0\eta_{0}>0, k>0, c>0c>0 and f>0f>0, sketch the fluid velocity when k(xct)=nπ/2k(x-c t)=n \pi / 2 for n=0,1,2,3n=0,1,2,3.