Paper 4, Section II, 36B

Waves
Part II, 2015

The shallow-water equations

ht+uhx+hux=0,ut+uux+ghx=0\frac{\partial h}{\partial t}+u \frac{\partial h}{\partial x}+h \frac{\partial u}{\partial x}=0, \quad \frac{\partial u}{\partial t}+u \frac{\partial u}{\partial x}+g \frac{\partial h}{\partial x}=0

describe one-dimensional flow over a horizontal boundary with depth h(x,t)h(x, t) and velocity u(x,t)u(x, t), where gg is the acceleration due to gravity.

Show that the Riemann invariants u±2(cc0)u \pm 2\left(c-c_{0}\right) are constant along characteristics C±C_{\pm}satisfying dx/dt=u±cd x / d t=u \pm c, where c(h)c(h) is the linear wave speed and c0c_{0} denotes a reference state.

An initially stationary pool of fluid of depth h0h_{0} is held between a stationary wall at x=a>0x=a>0 and a removable barrier at x=0x=0. At t=0t=0 the barrier is instantaneously removed allowing the fluid to flow into the region x<0x<0.

For 0ta/c00 \leqslant t \leqslant a / c_{0}, find u(x,t)u(x, t) and c(x,t)c(x, t) in each of the regions

 (i) c0txa (ii) 2c0txc0t\begin{gathered} \text { (i) } \quad c_{0} t \leqslant x \leqslant a \\ \text { (ii) }-2 c_{0} t \leqslant x \leqslant c_{0} t \end{gathered}

explaining your argument carefully with a sketch of the characteristics in the (x,t)(x, t) plane.

For ta/c0t \geqslant a / c_{0}, show that the solution in region (ii) above continues to hold in the region 2c0tx3a(c0t/a)1/32c0t-2 c_{0} t \leqslant x \leqslant 3 a\left(c_{0} t / a\right)^{1 / 3}-2 c_{0} t. Explain why this solution does not hold in 3a(c0t/a)1/32c0t<x<a.3 a\left(c_{0} t / a\right)^{1 / 3}-2 c_{0} t<x<a .