 Solution

Interpret ${\rho }_2$ as the ambient density: the PDF's phrase assigning it a “depth” is a typographical error. Take $0<{\rho }_1-{\rho }_2\ll {\rho }_2$, and use a common reference density in inertia while retaining the density difference in buoyancy. This is the Boussinesq approximation, with reduced gravity $g^{\prime }=g({\rho }_1-{\rho }_2)/{\rho }_{\mathrm{ref}}>0$. For the shallow water approximation, the lower-layer depth must be small compared with longitudinal variation scales, vertical accelerations must be negligible and pressure hydrostatic. The deep ambient is approximated as a hydrostatic pressure reservoir away from the nose; assume a uniform velocity over each lower-layer cross-section and neglect drag and entrainment.

Write the prismatic triangular area as $A(h)=k{h}^2$ with fixed $k>0$, so the width at height $z$ is $2kz$. The excess hydrostatic pressure force divided by reference density isThe prismatic-channel shallow water equations are thereforeFor positive depth, divide the first equation by $2kh$ and use it in the second to obtainThe coefficient matrix in variables $(h,u)$ is $\left(\begin{array}{cc}u& h/2\\ g^{\prime }& u\end{array}\right)$. Its two distinct real characteristic speeds areThus the system is strictly hyperbolic for $h>0$ and $g^{\prime }>0$, though it degenerates at zero depth. Since $dc/dh=c/(2h)$, direct substitution in the two equations givesHence $R_{\pm }=u\pm 4c$ are the Riemann invariants conserved along their respective characteristic curves. The factor four, rather than the rectangular-channel factor two, comes from the triangular relation $A\propto h^2$.