 Solution

Use the flux normalization of the supplied equations, with the common factor $\pi$ suppressed. Put $q=w{b}^2$, $m=w^2{b}^2$ and $B_0=g^{\prime }q$. The top-hat plume model reduces toEliminate height to obtain $dm/dq=B_{0}q/(2\alpha m^{3/2})$. Integrating with zero source volume and momentum flux gives the pure plume balanceEquivalently, seek $b=dz$ and $w=s{z}^{-1/3}$. The volume balance gives $d=6\alpha /5$, and the momentum and buoyancy balances give $g^{\prime }=(4s^2/3)z^{-5/3}$ and $B_0=4d^2{s}^3/3$. ThusA finite source replaces the singular point origin by a plume virtual origin and source conditions. The Boussinesq approximation requires the fractional density difference to be small, so density can be replaced by ambient density in inertia and entrainment while retaining the small difference in buoyancy. In this solution it requires $g^{\prime }/g\ll 1$. The approximation necessarily fails sufficiently near the ideal point source, where $g^{\prime }$ diverges, and can also fail for strongly heated or otherwise large-density-contrast sources. The similarity solution is appropriate farther away once entrainment has diluted the anomaly, assuming a uniform ambient, constant entrainment coefficient and negligible confinement, heat loss and stratification.