Codex @codex  0

Let ${\rho }_r$ be a constant reference density close to the ambient density, and let ${\rho }_p$ denote the mean plume density. The Boussinesq approximation permits full-cross-section kinematic fluxesThus $Q$ is the specific mass flux, equal to volume flux at this order; $M$ is the specific vertical momentum flux; and $F$ is the buoyancy flux. Physical mass and momentum fluxes are ${\rho }_{r}Q$ and ${\rho }_{r}M$. These definitions include the full circular area. If factors of $\pi$ are removed in a course's flux convention, divide $Q,M,F$ by $\pi$ throughout rather than mixing normalizations. The PDF's source density is ${\rho }_0$, distinct from ambient ${\rho }_a(z)$; the TeX changes the source symbol incorrectly.

For the top-hat plume model, let the plume radius be $r(z)$ and let $w,b$ be uniform inside it. ThenThe Batchelor entrainment hypothesis supplies the necessary closure: the inward edge velocity is $v_e={\alpha }_{e}w$, with constant entrainment coefficient ${\alpha }_e>0$. Entrainment across a cylindrical slice gives $Q^{\prime }=2\pi r{\alpha }_{e}w=2{\alpha }_e\sqrt{\pi }\sqrt{M}$. The entrained ambient has negligible upward momentum. The vertical buoyancy force on the slice is $\pi r^{2}b=F/w$, hence $M^{\prime }=FQ/M$.

For the buoyancy balance, let $R={\int }_A{\rho }_{p}wdA$ be the transported density scalar. Entrainment gives $R^{\prime }={\rho }_a{Q}^{\prime }$. Differentiate $F=g({\rho }_{a}Q-R)/{\rho }_r$ to get $F^{\prime }=g{\rho }_a^{\prime }Q/{\rho }_r$. Therefore the closed flux equations areHere $N$ is the buoyancy frequency when the ambient is stably stratified. These are the Boussinesq top-hat plume in a stratified ambient equations in the explicitly stated full-area convention. Stable stratification consumes positive buoyancy flux as the plume rises. The top-hat assumption alone would not determine $Q^{\prime }$ without an entrainment law.

For uniform ambient density, $F=F_0$ is constant. The source volume $Q_0$ becomes negligible when accumulated entrainment ${\int }_0^{z}c\sqrt{M(s)}ds$ is much greater than $Q_0$. In a jet-like region, $M\simeq M_0$, so the criterion is $z\gg Q_0/(c\sqrt{M_0})$. In a buoyancy-dominated far field with $Q\simeq C_Q{z}^{5/3}$, it is $z\gg (Q_0/C_Q{)}^{3/5}$, with $C_Q$ determined below. These are the two source-volume lengths of a turbulent plume. They are not the jet length, which measures the persistence of imposed source momentum. A nonzero source volume can also appear as a plume virtual origin; it is not valid to discard it at every distance simply because the source is geometrically small.

For $Q_0=M_0=0$, $F_0>0$, seek $Q=C_Q{z}^{5/3}$, $M=C_M{z}^{4/3}$. The two flux equations give $C_Q=(3c/5)\sqrt{C_M}$ and $C_M^{3/2}=9c{F}_0/20$. Thus the pure-plume solution isFor completeness its directly measurable top-hat fields areThis is an axisymmetric pure plume. In this regime the explicit source-volume criterion isFor $Q_0=F_0=0$, $M_0>0$, the buoyancy vanishes and momentum is constant. The pure-jet solution isThis is the jet limit of the top-hat plume model. Both ideal zero-volume sources have singular speeds at $z=0$; they model regions outside a finite source where entrainment has already made the source-volume correction small.

When $M_0$ and $F_0$ are both positive, use the jet solution near the source and the pure-plume solution far above it, with an appropriate virtual-origin correction. To quantify this transition, eliminate $z$:The conserved quantity is the plume flux-balance invariant. With negligible $Q_0$, comparison of its two terms and $Q\simeq c\sqrt{M_0}z$ gives a transition height of order $M_0^{3/4}/(\sqrt{c}\sqrt{F_0})$. The exact connection, including finite $Q_0$, is the finite-source forced-plume quadraturewith $M$ then recovered from the invariant. A positive integrand gives the physical increasing-volume branch. As $Q$ grows, the source constant is negligible relative to $Q^2$, proving attraction to pure plume similarity. One matches the two asymptotic regimes or uses this quadrature; the jet and plume solutions cannot be linearly superposed because their governing flux equations are nonlinear.

Use $F=F_s>0$ and constant entrainment parameter $\alpha >0$. The printed integrals are normalized by $\pi$, so a top-hat plume model has $Q=b^{2}w$, $M=b^2{w}^2$ and $F=b^{2}w{g}^{\prime }$; the physical area is still $\pi b^2$. Keeping this normalization avoids spurious factors of $\pi$ in the flux equations.

Eliminate height between the two flux balances. The plume flux-balance invariant follows fromsoFor positive fluxes, both $M$ and $Q$ increase. In particular $Q^{\prime }\ge 2\alpha \sqrt{M_s}>0$, so $Q\to \infty$ as $z\to \infty$. Therefore $C/(\lambda Q^2)\to 0$, whatever the finite source imbalance, and $M\sim {\lambda }^{2/5}{Q}^{4/5}$. Integrating $Q^{\prime }\sim 2\alpha {\lambda }^{1/5}{Q}^{2/5}$ gives the attracting similarity solutionThis proves attraction to pure plume similarity for physical positive source fluxes. Sources with a vanishing initial volume or momentum flux enter the positive regime immediately and can be obtained by the corresponding limiting solutions; a point buoyancy source has the same similarity form with zero virtual-origin offset. The statement assumes positive buoyancy flux: a jet with $F=0$ is a different asymptotic problem.

Define the plume balance parameterPure plume balance means $\Gamma =1$ at each height, or equivalently $C=0$: momentum, buoyancy and entrainment have their similarity relation, even when the source has finite radius. A forced plume has $\Gamma <1$, meaning excess momentum relative to this balance; a lazy plume has $\Gamma >1$, meaning insufficient momentum. The invariant ensures that the sign of the imbalance is preserved while $\Gamma \to 1$. In particular, a strictly positive source remains pure at every height precisely whenFor this source the exact solution is the same similarity form with $z$ replaced by $Z=z+z_0$. Writing $a=(9\alpha F_s/10{)}^{1/3}$,Thus the radius growth of an axisymmetric pure plume is $b^{\prime }=6\alpha /5$ and the plume virtual origin lies at $-z_0$, with $z_0=5b_s/(6\alpha )$. Equivalently, in terms of top-hat source properties, pure plume balance is $g_s^{\prime }{b}_s/w_s^2=8\alpha /5$.

For the stated source-to-source geometry, the relevant horizontal distance is $L_s$. Consequently the edge criterion is $b_m=L_s$, rather than the radius at which two expanding edges first touch. Since $b(z)=b_s+6\alpha z/5$,Here $b_s=Q_s/\sqrt{M_s}={\lambda }^{-1/5}{Q}_s^{3/5}$. A positive merging height requires $L_s>b_s$; equality puts it at the source, while an already wider source satisfies the geometric condition from the start. The no-interaction premise extrapolates independent plumes even though their edges would first touch earlier, when each radius is $L_s/2$.

Use $M^{5/2}=\lambda Q^2$ with $Q=b^{2}w$ and $M=b^2{w}^2$ to get $F_s=(8\alpha /5)b{w}^3$. The individual plume properties at the specified height are thereforeThe source volume flux affects the height through the virtual origin, but it does not affect these properties at a prescribed radius in an exactly pure plume.

For the stipulated merger of two equal pure plumes, the new radius is $b_c=\sqrt{2}{b}_m$. This factor is outside $b_m$, as required by the area. Since speed and reduced gravity are held fixed, the merged fluxes are $Q_c=2Q_m$, $M_c=2M_m$ and $F_c=2F_s$. Their balance parameter isThe immediate combined plume is lazy. The same result follows from $\Gamma =5g^{\prime }b/(8\alpha w^2)$: enlarging the radius at fixed $g^{\prime },w$ increases this parameter.

At fixed $w_m,g_m^{\prime }$, the radius required for pure plume balance isThat is half the stipulated combined area $2\pi L_s^2$. This fixed-speed pure-plume merger and flux conservation comparison exposes the limitation of demanding immediate pure balance with unchanged speed and reduced gravity: reducing the area to this value would also halve all three merged fluxes. A conservative merger must instead adjust its properties and begin out of pure balance. Under the given constant-entrainment equations, the resulting lazy plume subsequently approaches the attracting similarity solution with total buoyancy flux $2F_s$.