Codex @codex  0

Use the steady axisymmetric top-hat plume model, with plume radius $b(z)$, uniform mean upward speed $w(z)>0$ and reduced gravity $g^{\prime }(z)=g[{\rho }_a(z)-{\rho }_p(z)]/{\rho }_0$. The Boussinesq approximation uses constant reference density ${\rho }_0$ in inertia and retains the small density difference only in buoyancy. Assume a quiescent hydrostatic ambient, a slender fully turbulent plume, negligible molecular transport and an entrainment velocity $\alpha w$, with constant positive entrainment coefficient $\alpha$. Define the ambient buoyancy frequency by $N^2=-(g/{\rho }_0)d{\rho }_a/dz$; positive $N^2$ is stable stratification.

The integrated volume, momentum and buoyancy balances of a Boussinesq top-hat plume in a stratified ambient areThe first follows from side entrainment, the second from the integrated buoyancy force, and the third accounts for the changing ambient density. Common factors of $\pi$ are cancelled. Equivalently, the requested radius, speed and reduced-gravity equations areThese equations stop describing an upward light plume once $w$ vanishes, and their light-plume interpretation requires $g^{\prime }>0$.

For nonzero $C$, seek plume similarity in power-law stratification in the form $b=\beta z$, $w=A{z}^n$ and $g^{\prime }=G{z}^m$. The first balance gives $\beta =2\alpha /(n+2)$. Momentum gives $m=2n-1$ and $G=2(n+1)A^2$. The buoyancy balance then requires $3n=2+n+p$, soChoosing $A>0$ gives the exact power-law family, including its normalization:These are the amplitudes and physical branches of a power-law stratified plume. No source flux was specified, so an arbitrary finite-source solution is not determined by $C,p$ alone; the displayed family is an exact similarity solution, rather than the general solution for arbitrary initial fluxes.

With the specific flux convention $Q=b^{2}w$, $M=b^2{w}^2$ and $B=b^{2}w{g}^{\prime }$, their values areThe physical volume flux is $\pi Q$, mass flux is ${\rho }_0\pi Q$, momentum flux is ${\rho }_0\pi M$, and kinematic buoyancy flux is $\pi B$. Thus the “specific mass flux” is the mass flux divided by reference density, with or without the explicitly stated factor $\pi$.

The parameter restrictions follow directly from positivity. Positive radius needs $p>-6$. Positive $g^{\prime }$ and real positive $w$ further require $p>-4$ and $-C/(3p+8)>0$, givingFor stable $C>0$, buoyancy flux decreases with height, while volume and momentum fluxes increase. However, $B\propto z^{4+3p/2}$ diverges at $z=0$ throughout this interval. Also the ambient density gradient and plume speed are singular there. Thus this stable similarity branch can describe a region away from a virtual source, but cannot represent a finite-buoyancy-flux point source extending all the way to $z=0$. For ordinary regular stable stratifications, including positive $p$, a finite-source rising plume instead loses buoyancy and eventually reaches neutral buoyancy and a finite rise limit; the pure-power ansatz does not supply such a solution.

For unstable $C<0$ and $p>-8/3$, the buoyancy flux grows with height and tends to zero at the origin. The plume draws its increasing buoyancy from an unstable ambient rather than a finite positive source buoyancy flux. Such an ambient must be maintained or treated locally, since it is itself convectively unstable. If $p\le -1$, its density profile is singular at the source, while for $p>-1$ its increasing density contrast eventually violates the Boussinesq approximation at sufficiently large height. These are formal local similarity solutions with corresponding physical limits. For $C<0$ and $-6<p<-4$, the algebra also gives a real positive speed but negative $g^{\prime }$: this is a heavy upward flow, not the specified light plume.

At nonzero $C$, the values $p=-6,-4,-8/3$ do not give finite coefficients. Respectively they would require zero entrained volume-flux growth, zero momentum-flux growth despite nonzero buoyancy, or constant buoyancy flux despite nonzero ambient stratification. They are not omitted regular endpoints.

When $C=0$, the ambient is unstratified regardless of the nominal value of $p$. The nontrivial positive-buoyancy point-source solution is the Boussinesq point-source plume:Here $B_0>0$ is the constant specific buoyancy flux with $\pi$ removed. Its $Q\propto z^{5/3}$ and $M\propto z^{4/3}$ vanish at the point source, while $B=B_0$. Its near-source singularity is the usual ideal point-source limit; a finite vent replaces that limit by finite source data or a virtual origin.

For the steady positive-velocity branch, conservation of buoyancy flux gives $bhu=F_0$. Substituting into the other two balances givesUse a power-law ansatz $h=C{x}^p$, $u=D{x}^q$. The entrainment balance forces $p=1$ and $C(1+q)=\alpha$. The momentum powers are $x^{2q}$, $x^{-q}$ and $x^{-q}$, so the nonzero inertial branch has $q=0$. Therefore the constant-speed entraining gravity current on a slope isIts positive speed requires $\tan \theta >\alpha /2$ for $\theta <\pi /2$. In the shallow, strictly hyperbolic range this meansThe same steady integral balances have a formal vertical-plume limit at $\theta =\pi /2$, although their time-dependent system then loses strict hyperbolicity. At the lower threshold the assumed positive-flux, finite-speed solution fails.

The line source is singular: $h,uh$ and the momentum flux tend to zero, but $b$ diverges as $x\downarrow 0$. Thus this is a downstream ideal-source solution, valid only far enough away that $b/g\ll 1$ and the Boussinesq approximation and top-hat plume model hold. It does not describe the immediate source region.

Write $G=g^{\prime }$ for the positive reduced gravity of dense fluid. With $z$ upwards, $G_z>0$ places denser fluid above lighter fluid and permits convective turbulence. A displacement over the radius scale $d$ samples a buoyancy difference of order $d{G}_z$; balancing kinetic energy per mass with buoyant work over $d$ givesThe buoyancy-gradient mixing-length closure then gives an eddy diffusivity $K\sim ud=d^2\sqrt{G_z}$. We use the paper's unit-prefactor convention for this dimensional closure and its effective cross-sectional area $d^2$; order-one mixing constants and the cylinder's geometric factor $\pi$ are suppressed in its displayed evolution equation.

The downward reduced-gravity transport magnitude and its convergence are different quantities:The printed expression labelled a flux is actually the second quantity, a flux convergence per unit height. The dimensions distinguish them: $\mathcal{F}$ has units $L^4{T}^{-3}$, whereas $F_{\mathrm{conv}}$ has units $L^3{T}^{-3}$. A constant positive gradient has nonzero transport but zero convergence, providing a direct counterexample to identifying the derivative with the transport itself.

The upward total buoyancy flux is $QG-\mathcal{F}$. Conservation therefore yieldsThis flux convergence in turbulent buoyancy mixing obtains the displayed evolution equation with the correct transport interpretation.

For the arrested cubic buoyancy profile, an undisturbed region below $-H$ has $G=0$ and no turbulent buoyancy flux. The joining condition is $G(-H)=0$ andIn a steady state, integrating the conservation equation from that front gives $QG=d^4(G_z{)}^{3/2}$. For a nonzero mixed region,and henceThe associated gradient is $Q^2(z+H{)}^2/(9d^8)$, so it matches the required zero-gradient state continuously. We take the mixed interval to start immediately above $-H$; adding a further zero interval would merely relocate the arrest front.

The source buoyancy flux is $B_s=Q_s{g}_s^{\prime }$, with $g_s^{\prime }=g\Delta \rho /\rho$ in the Boussinesq approximation. The stated top-source convention means that this entire flux enters the downward turbulent transport there: $\mathcal{F}(0)=B_s$. Since steady transport also satisfies $\mathcal{F}(0)=QG(0)$,This is the extent in the paper's dimensional-closure normalization. It presupposes positive $Q$; with no upflow there is no finite steady arrest depth. For a mixture of the two original fluids, $G(0)\le g_s^{\prime }$ is necessary, so an extrapolation to $Q_s/Q>1$ cannot represent a physical concentration profile under this top-flux boundary model. A dilute-source regime avoids appreciable source-volume corrections.

Assuming the same scalar-mixing coefficient for a passive tracer, its eddy diffusivity in the mixed region isIt vanishes at the arrest front and increases linearly towards the source. A different turbulent scalar diffusivity ratio would multiply this result by its corresponding constant.

The geometric normalization can be restored without changing the reasoning. With actual area $S=\pi d^2$ and $K=C_m{d}^2\sqrt{G_z}$, put $R=S{C}_m{d}^2$. Then $G=Q^2(z+H{)}^3/(27R^2)$, $H=3R^{2/3}{B}_s^{1/3}/Q$, and $K=Q(z+H)/(3S)$. These formulas show which numerical prefactors are convention-dependent.

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$.

At a fixed height, first determine the stored quantities per unit vertical height and per unit length along the line source. With $\Delta \rho =\rho -{\rho }_0$, the triangular-profile line plume has the useful profile integralThe geometric volume is $\mathcal{A}=2b$; the mass is $\mathcal{D}={\rho }_0{I}_0+\Delta \rho I_1=b(\rho +{\rho }_0)$. The upward momentum is $\mathcal{P}=W({\rho }_0{I}_1+\Delta \rho I_2)$, while the buoyancy force is $\mathcal{B}=-g\Delta \rho I_1$. ThusHere positive buoyancy corresponds to a rising light plume. If instead buoyancy is defined kinematically, divide $\mathcal{B}$ by ${\rho }_0$; the density-weighted convention here is the one matching the printed equations.

Multiply each local density by the vertical velocity to obtain the fluxes through a horizontal section, still per unit source length:In particular, the stored momentum equals the mass flux $Q$, and the buoyancy flux is $F=g({\rho }_{0}V-Q)$. The different triangular-profile factors in storage and flux must not be replaced by top-hat coefficients.

Use the steady Boussinesq top-hat plume in a stratified ambient, with negligible ambient vertical momentum and no externally added heat or particles. The Batchelor entrainment hypothesis supplies inward velocity $\alpha w$ across circumference $2\pi b$. Volume conservation gives $\pi Q_z=2\pi b\alpha w$, and $bw=\sqrt{M}$, henceThe momentum flux increases by the integrated buoyancy force. Writing $g^{\prime }=g\Delta \rho /{\rho }_o$, the force per reference density and per $\pi$ is $g^{\prime }{b}^2$. Since $F=g^{\prime }Q$ and $b^2=Q^2/M$, this givesFor the buoyancy flux, entrained material has the local environmental mass density. Let ${\rho }_p$ be the plume's averaged density. Its scalar mass conservation is $({\rho }_{p}Q{)}_z={\rho }_a{Q}_z$. Since $F=(g/{\rho }_o)({\rho }_a-{\rho }_p)Q$, differentiating cancels both entrainment terms and givesIn stable density stratification, ${\rho }_a^{\prime }<0$ and a rising plume loses buoyancy flux. These are kinematic fluxes with the common factor $\pi {\rho }_o$ removed, so no additional factor of $\pi$ belongs in the three differential equations.

Define $K=N_s^2{z}_s^{-\beta }>0$. The buoyancy frequency definition gives $F_z=-KQ{z}^{\beta }$. The powers in a plume similarity in power-law stratification satisfySolving givesA rising positive-buoyancy turbulent plume entrains, so $Q_z>0$; its momentum flux increases, so $M_z>0$; its buoyancy flux decreases in the stable ambient, so $F_z<0$. With positive coefficients, this requires $q>0$, $m>0$ and $f<0$. The strongest restrictions areIn this interval $1<q<5/3$, $0<m<4/3$ and $-2<f<0$.

Coefficient matching makes the positive similarity solution explicit:soIt is an ideal scaling for $z>0$, often a far-field description. As $z\to 0$, the prescribed ambient stratification and this model's buoyancy flux are singular; a finite physical source would require a regularized near-source region.

Expand the square root and quotient in the logarithmic-height plume equations:The linear stability analysis givesIn particular, the second row contains $-2m{M}^{\prime }$, since the reciprocal of $\hat{M}$ contributes another $-M^{\prime }$.

For the growing mode of power-law plume similarity, the characteristic polynomial isWithin the physical interval $q,m>0$ and $f<0$, $P(0)<0$ while $P(\lambda )\to +\infty$ as $\lambda \to +\infty$. There is a positive real eigenvalue, so the similarity equilibrium is unstable with increasing height. Perturbations behave as $e^{\lambda \xi }\propto z^{\lambda }$; the result is not stability in physical time. Thus this exact positive-buoyancy similarity trajectory is not a generic attracting far field.

At the upper endpoint,The limiting polynomial is $\lambda (\lambda +1)(\lambda +10/3)$, reproducing the three stated eigenvalues. The positive eigenvalue tends to zero, with $\lambda \sim -2f/3$. The two negative modes decay as $z^{-1}$ and $z^{-10/3}$; the remaining mode is neutral at the limiting linear system. Its eigenvector is $(1,2,3)$ in relative-flux coordinates.

The neutral buoyancy-flux mode of a pure plume has a simple interpretation: a constant surviving buoyancy flux fixes the amplitudes of a pure plume, with $Q\propto F^{1/3}{z}^{5/3}$ and $M\propto F^{2/3}{z}^{4/3}$. A small fractional change in $F$ produces changes in $Q,M$ in the ratio $1:2:3$. It is retained rather than damped out.

There is a qualification at the exact endpoint. For fixed nonzero $K$ and $\beta =-8/3$, taking $Q\propto z^{5/3}$ gives $F_z\propto -z^{-1}$, not zero. The flux coefficients above also diverge as $f\to 0$. Thus the constant-flux powers are the formal limiting scaling, not an exact constant-buoyancy solution in a nonzero endpoint stratification; cumulative ambient effects can require logarithmic corrections.

Seek a steady downslope solution with $u>0$ constant. The volume conservation equation reduces to $h_x=E$, so the constant-speed entraining gravity current on a slope hasThe conserved buoyancy flux per unit width is $B=u{g}^{\prime }h$, givingSubstitute these expressions into the momentum equation. Its advective term is $E{u}^2$ and its hydrostatic term is $\frac{1}{2}E(B/u)\cos \theta$. The balance is thereforewhich givesThe assumed slope inequality makes ${\beta }_{\ast }>0$. The increasing thickness dilutes reduced gravity while leaving $g^{\prime }h$ and the buoyancy flux constant. The denominator includes both bottom drag and the momentum needed to accelerate newly entrained ambient fluid.

For a steady no-source particle-laden gravity current, $B_x=-v_s{g}^{\prime }<0$. Simply entraining particle-free ambient fluid increases volume and dilutes particle volume fraction; it cannot create particle buoyancy flux.

A steady increase can occur if sediment resuspension or bed erosion adds more suspended particle buoyancy than settling removes, if a tributary or lateral particle source supplies a positive flux, or if the buoyancy per unit particle amount increases through a change of the surrounding fluid's density. Changes of particle properties or another buoyancy source can also alter the flux, but they go beyond the fixed-particle, fixed-reference model. In a model with a source $S_B$,Thus the buoyancy production in a particle-laden current requires $S_B>v_s{g}^{\prime }$ in a steady current. In an unsteady current, $B_x$ can be positive even without external addition: local depletion has ${\partial }_t(g^{\prime }h)<-v_s{g}^{\prime }$, so more buoyancy leaves a section than enters it while stored particles are being released. An instantaneous spatial flux increase is not by itself proof of particle production.