Past exam of the mathematics course of the University of Cambridge / 2017 / ib / Paper 2 / 7D (Fluid Dynamics)

For constant fluid density $\rho$, steady Euler equations for an inviscid fluid under a conservative force per unit mass $-\nabla \Phi$ read $(\mathbf{u}\cdot \nabla )\mathbf{u}=-\nabla p/\rho -\nabla \Phi$. Here $\mathbf{u}$ is velocity, $p$ pressure and $\Phi$ potential per unit mass. Dotting with $\mathbf{u}$ givesThis is the Bernoulli equation. The constants can differ between streamlines. If density varies, the same derivation for a barotropic fluid replaces $p/\rho$ by ${\int }^{p}d{p}^{\prime }/\rho (p^{\prime })$; constant density is the appropriate water approximation here.

For the draining cone, the radius at height $h$ is $h$ and the horizontal area is $\pi h^2$. Both the free surface and the jet are at atmospheric pressure; their gravitational potential difference is approximately $gh$. The speed of the descending surface is smaller than the jet speed by the area ratio ${ϵ}^2/h^2$, so it can be neglected when $h\gg ϵ$. The Bernoulli equation then gives Torricelli's law, $U\simeq \sqrt{2gh}$. Conservation of mass givesHence the leading drainage time of a conical tank isRemoving the tip places the aperture at height of order $ϵ$, and the quasi-steady calculation fails in the last layer of that depth. Integrating down only to $h=O(ϵ)$ changes the displayed leading time by relative order $(ϵ/h_0{)}^{5/2}$; replacing the actual head $h-ϵ$ by $h$ introduces a larger, but still vanishing, relative correction of order $ϵ/h_0$. The formula is the ideal inviscid estimate under the question's stated approximation, not an exact law for the final transient or a real contracted jet.