Non-Boussinesq top-hat plume equations

For a turbulent plume in a homogeneous ambient of density ${\rho }_0$, define actual sectional mass flux, momentum flux and density-weighted buoyancy flux by $Q=\pi b^2\rho w$, $M=\pi b^2\rho w^2$ and $F=\pi b^{2}wg({\rho }_0-\rho )$. With non-Boussinesq Batchelor entrainment $v_e=\alpha w\sqrt{\rho /{\rho }_0}$, their balances areHere $Q^2/M$ is sectional mass, $QF/M$ is sectional density-deficit buoyancy, and the entrained ambient initially has zero vertical momentum. The last equation follows by subtracting mass conservation from ${\rho }_0$ times volume conservation. Scase, Caulfield, Dalziel and Hunt's time-dependent plume model derives the equivalent system with the common factor $\pi$ suppressed.