Axisymmetric pure plume

Consider a steady, axisymmetric pure plume in a uniform ambient fluid, using the Boussinesq approximation and uniform top-hat cross-sectional profiles. If $r(z)$ is its radius, $U(z)>0$ its vertical velocity and $g^{\prime }(z)>0$ its reduced gravity, define the physical fluxesThey are respectively volumetric flow rate, momentum flux per reference mass density, and buoyancy flux. Only $B$ is conserved in this unstratified, source-free model. Fluid entrainment adds ambient volume to $Q$, and buoyancy increases $M$.

For constant entrainment coefficient $\alpha >0$, the Batchelor entrainment hypothesis gives an inward boundary speed $\alpha U$. Volume, vertical momentum and buoyancy balances therefore giveTo find a pure similarity solution, put $Q\propto z^p$ and $M\propto z^q$ relative to the plume virtual origin. The first balance gives $p-1=q/2$; the second gives $q-1=p-q$. Thus $p=5/3$, $q=4/3$, andIn particular $r=Q/\sqrt{\pi M}$ grows linearly while $U=M/Q$ decreases as $z^{-1/3}$. A finite source shifts the virtual-origin coordinate. Dividing all three fluxes by $\pi$ changes their balance coefficients, so pi-normalized and physical-flux formulas must not be mixed.

For pi-normalized top-hat fluxes, $b=Q/\sqrt{M}$. The pure similarity solution therefore has $b=(6\alpha /5)(z+z_0)$, so the radius grows linearly with slope $6\alpha /5$. The corresponding plume virtual origin is $-z_0$ and a finite source has $z_0=5b_s/(6\alpha )$. This law translates geometric plume-contact criteria into heights.

A constant surviving buoyancy flux sets the amplitudes of a pure plume: $Q\propto F^{1/3}$ and $M\propto F^{2/3}$. Relative-flux perturbations have the neutral ratio $(\delta Q/Q,\delta M/M,\delta F/F)\propto (1,2,3)$. This is the zero-eigenvalue limit of growing mode of power-law plume similarity; nonzero endpoint stratification can still create cumulative logarithmic flux loss.

A steady top-hat plume model supplied with pi-suppressed kinematic buoyancy flux $f_0=b^{2}w{g}^{\prime }>0$ but zero source volume flux and momentum flux hasThe full-cross-section buoyancy flux is $\pi f_0$. The Batchelor entrainment hypothesis and momentum balance yield these constants. Its Froude number satisfies $w^2/(g^{\prime }b)=5/(8\alpha )$. The Boussinesq approximation requires $g^{\prime }/g\ll 1$ and therefore fails at the ideal point source itself.