Non-Boussinesq pure-plume density transition ($\rho ={\rho }_0/[1+(z_b/z{)}^{5/3}]$)

For an ideal-gas heated pure plume with entrainment proportional to $w\sqrt{\rho /{\rho }_0}$, the excess heat flux implies conserved kinematic buoyancy flux $B$. Put $q=Q/{\rho }_0$, $m=M/{\rho }_0$. Their balances are $q^{\prime }=2\alpha \sqrt{m}$ and $m^{\prime }=Bq/m$, identical to the pi-suppressed Boussinesq pure-plume flux equations. With $d=6\alpha /5$ and $s=[25B/(48{\alpha }^2){]}^{1/3}$, $q=d^{2}s{z}^{5/3}$, $m=d^2{s}^2{z}^{4/3}$ and $z_b^{5/3}=B/(g{d}^{2}s)$. The radius is $b=dz\sqrt{1+(z_b/z{)}^{5/3}}$ and reduced gravity is $g/[1+(z/z_b{)}^{5/3}]$. At $z=z_b$ mass density is half ambient; far above this height the Boussinesq approximation becomes accurate.