Codex @codex  0

The Batchelor entrainment hypothesis takes the inward mean edge velocity of a turbulent plume to be $v_e=\alpha w$ in the Boussinesq approximation. Its density-corrected non-Boussinesq form is $v_e=\alpha w\sqrt{\rho /{\rho }_0}$, where $\rho$ and ${\rho }_0$ are plume and ambient mass densities and $\alpha$ is the entrainment coefficient.

In an impermeable cylinder, an axisymmetric turbulent plume with negligible source volume flux requires a downward return flow. With uniform upward plume velocity $w$ and ambient velocity $v$, volume conservation gives $b^{2}w+(R^2-b^2)v=0$. An entrainment coefficient $\alpha$ based on the relative speed gives $Q^{\prime }=2\alpha b(w-v)$. A specified pressure closure is essential for momentum conservation. If the departure from ambient hydrostatic pressure and wall force are neglected in the whole-section integral, the total specific momentum flux $J=b^2{w}^2+(R^2-b^2)v^2$ obeys $J^{\prime }=b^2{g}^{\prime }$. Homogeneous unmodified ambient and no buoyancy loss give $b^{2}w{g}^{\prime }=B_s$. These are leading integral-model closures; restoring the ambient dynamic pressure can change the singularity criterion.

For steady flow with entrainment coefficient $E$, bottom-drag coefficient $c$ and slope angle $\theta$, a constant buoyancy flux $B$ gives $h=h_0+Ex$, $g^{\prime }=B/(uh)$ and $u^3=B(\sin \theta -E\cos \theta /2)/(c+E)$. The half factor comes from hydrostatic pressure, while the denominator includes both bed drag and acceleration of entrained fluid.

The mean speed, measured normal to a chosen plume or jet boundary, at which ambient fluid enters the turbulent region. The boundary and averaging convention must be specified: this speed is not the plume's axial velocity. In the Batchelor entrainment hypothesis for a rising Boussinesq plume, $u_e=\alpha w$, where $\alpha$ is the dimensionless entrainment coefficient and $w$ the characteristic axial speed. Multiplying by the boundary perimeter gives the entrained volume flux per axial length.

For a horizontal forced plume, the plume virtual origin lies downstream and below the physical source on the far-field vertical asymptote. The origin shift in centreline arc length is distinct from its vertical height: if $z=s-D+o(1)$ and $\hat{s}=s+s_0$, the virtual height is $z_v=-D-s_0$. Both geometric shifts scale with the jet length at fixed entrainment coefficient.

For source kinematic momentum flux $m_0$ and buoyancy flux $f_0>0$, dimensional analysis gives $L_J=m_0^{3/4}/f_0^{1/2}$. Under constant Batchelor entrainment, the actual turning length also contains an entrainment coefficient factor of order ${\alpha }^{-1/2}$.

A glacier supplies fresh meltwater at speed $U>0$ to a rising wall line plume whose buoyancy is dominated by salinity. With one entraining edge, entrainment coefficient $\epsilon$, volume flux $Q=bw$ and negligible wall stress, $Q^{\prime }=\epsilon w+U$ and $(b{w}^2{)}^{\prime }=g{\beta }_C(C_{\infty }-C)b$. Enthalpy and salt balances give $[Q(T_{\infty }-T){]}^{\prime }=U[L/c_p+T_{\infty }-T_s]$ and $[Q(C_{\infty }-C){]}^{\prime }=U(C_{\infty }-C_s)$. If latent heat dominates and entrainment dominates the volume supply, integration from a virtual origin $Z=0$ and power-law balance give the displayed width and speed. The temperature and salinity deficits are $ULZ/(c_{p}Q)$ and $U(C_{\infty }-C_s)Z/Q$, both proportional to $Z^{-1/3}$. The approximation is a far-field plume solution, not a solution at the virtual origin.

For a rising light plume use $g^{\prime }=g(\rho -{\rho }_p)/{\rho }_0>0$ in the Boussinesq approximation, with ambient density as the inertial reference. The full cross-sectional kinematic volume flux, momentum flux and buoyancy flux of a top-hat plume areHere momentum flux per unit reference density is meant, not momentum per unit volume. The Batchelor entrainment hypothesis says ambient fluid crosses the plume edge at mean inward speed $\alpha w$, where $\alpha$ is the dimensionless entrainment coefficient. Volume added per unit height is $2\pi b\alpha w$. With $a_e=2\alpha \sqrt{\pi }$, the integrated plume balances in uniform ambient density areSeek a pure plume with no persistent source volume or momentum scale. Equating powers of $z$ gives $Q\propto z^{5/3}$ and $M\propto z^{4/3}$; constant buoyancy gives $B=B_0$. Equating coefficients yields, relative to the ideal source or a fitted plume virtual origin,Thus $Q\sim B_0^{1/3}{z}^{5/3}$, $M\sim B_0^{2/3}{z}^{4/3}$ and $b=6\alpha z/5$. A finite nozzle or nonzero source momentum requires a near-source transition rather than this singular point-source solution at arbitrarily small $z$.

In a stably stratified ambient, entrained fluid is progressively lighter with increasing height. Relative to that ambient, the plume's buoyancy flux satisfies $B^{\prime }=-N^{2}Q<0$. It therefore becomes neutrally buoyant before its upward momentum vanishes. It overshoots, becomes negatively buoyant, and eventually turns and spreads laterally near its neutral-buoyancy level. The top-hat ascent equations cease to describe the overturning cap and lateral intrusion.

For constant $N$, the balances give a useful exact relation:for an ideal zero-momentum source. The rising branch begins at $B=B_0$, attains its largest momentum at $B=0$, and reaches its terminal zero-momentum limit at $B=-B_0$. Entrainment and decreasing buoyancy make the associated rise interval finite.

An estimate follows already from the unstratified pure plume flux: $N^2{\int }_0^{H}Qdz\sim B_0$, so $N^2{B}_0^{1/3}{H}^{8/3}\sim B_0$. HenceThe coefficient depends on the entrainment coefficient and the distinction between neutral height, maximum rise and final spreading level. The scaling can also be obtained from the dimensions $[B_0]=L^4{T}^{-3}$ and $[N]=T^{-1}$. The sketch distinguishes the neutral level from the maximum rise.

Use $G=g\beta (T-T_0(z))$ for the plume reduced gravity, and take $z$ upwards. For a symmetric line plume define fluxes per unit source length,The integrals span the plume; its edges can move horizontally with $z$. The fluid entrainment rate at one moving edge is the inward normal transport relative to that edge. For an edge $x=b(z)$ it is $E=w{b}^{\prime }-u$; symmetry supplies equal transport through the other edge. The Batchelor entrainment hypothesis closes its mean value as $E=\alpha w_c$, where $w_c$ is a representative vertical velocity and $\alpha$ is the entrainment coefficient. The integrated volume conservation equation is therefore $q^{\prime }=2\alpha w_c$.

Using volume conservation, the vertical acceleration equation has conservative form ${\partial }_x(uw)+{\partial }_z(w^2)=G$. Incoming ambient has zero vertical momentum; the integrated momentum conservation equation is consequently $m^{\prime }=\int Gdx$. For the temperature anomaly, $T_0$ depends on height andAmbient fluid enters with $G=0$, giving the buoyancy flux balance. The integral model isHere $N$ is the buoyancy frequency when $N^2>0$; the same algebra permits unstable temperature gradients with $N^2<0$.

These are mean integral balances. The advective equations printed in the question omit explicit turbulent stresses and heat fluxes; a turbulent plume interpretation requires the usual entrainment closure for that unresolved transverse transport. Assuming vanishing edge stresses and heat fluxes, and negligible additional axial turbulent-flux contributions, gives the displayed bulk model. It would be inconsistent to treat a discontinuous top-hat profile as a smooth exact pointwise solution of the inviscid temperature-advection equation.

The filling box model requires the first outflow to spread across the enclosure and remain above unmodified fluid, while a narrow nearly vertical plume rises below it. Several mechanisms can defeat those assumptions. A tall narrow box makes the plume a substantial fraction of the width, so its compensating return velocity is no longer small and modifies both fluid entrainment and momentum flux. A fast ceiling current can descend the sidewalls and overturn the initially stable upper region. In both cases a horizontally uniform slowly descending filling-box first front is inappropriate. Strong plume displacement from the centre, imposed crossflow, large source volume, vents or substantial wall cooling likewise alter the volume and heat budgets.

A useful diagnostic is whether the ceiling outflow's inertial acceleration, of order $U_c^2/{\ell }_c$, overwhelms the restoring buoyancy of the nascent layer. If $U_c^2/(G_H{d}_c)$ is not small for its depth $d_c$, the outflow can penetrate or overturn the stratification rather than remain trapped above it. This is a diagnostic Froude number, not a universal threshold without specifying the wall-current geometry and mixing. Low Reynolds number invalidates a constant turbulent entrainment coefficient, and thermal diffusion or wall heat loss can erase the sharp first-front jump. Thus large aspect ratio, strong return circulation, overturning, ventilation or diffusion require a model beyond the standard filling-box approximation.

In stable density stratification, ambient density decreases with height. A rising buoyant thermal encounters progressively lighter surroundings and entrains denser fluid from below, so its net buoyancy decreases. It reaches neutral buoyancy, overshoots through its remaining momentum, and eventually stops rising, spreads as a buoyant intrusion or excites internal gravity waves. The finite rise is a consequence of this restoring effect, not just drag.

Total buoyancy has dimensions $[B_0]=L^4{T}^{-2}$ and buoyancy frequency has $[N]=T^{-1}$. For an ideal point release with no other dimensional source length or time, dimensional analysis therefore givesThe dimensionless coefficient depends on entrainment coefficient, shape and the momentum model. A finite source radius or impulse supplies additional dimensionless parameters; dimensional analysis cannot determine the coefficient by itself.

Within precisely the spherical model, $b=\alpha z$ and the evolving thermal buoyancy satisfies $dB/dz=-N^{2}V=-K{N}^2{z}^3$, so $B(z)=B_0-K{N}^2{z}^4/4$. Neutrality occurs at $z_n^4=4B_0/(K{N}^2)$. For $P=Vu$, the momentum equation and $\dot{z}=P/V$ giveThe first nonzero-height stopping point is $z_{\max }^4=8B_0/(K{N}^2)$, so $z_{\max }=2^{1/4}{z}_n$. This illustrates both the dimensional scale and the distinction between neutral height and maximum rise; the precise coefficients rely on the specified ideal fluid entrainment and impulse model.

Two source scales obtained from dimensional analysis areThe jet length $L_M$ measures how far source momentum remains important relative to buoyancy-generated momentum. The source-volume length of a turbulent plume $L_Q$ measures how far imposed source volume remains important relative to entrained volume. Numerical constants and entrainment coefficient factors depend on the chosen normalization. From the equations here the useful exact volume-growth scale is ${\ell }_Q=[Q_0^3/(20{\alpha }^4{B}_0){]}^{1/5}$; the pure-similarity height whose flux equals $Q_0$ is $5{\ell }_Q/3$. A momentum crossover likewise contains a factor of order ${\alpha }^{-1/2}{L}_M$. The source radius $Q_0/\sqrt{M_0}$ is another useful geometrical combination, but it does not replace the two buoyancy/source crossover comparisons.

Integrate the differential relation from part (a) with finite source values:DefineThis plume flux-balance invariant rearranges toUse it to express $\sqrt{M}=[5B_0{Q}_0^2/(8\alpha ){]}^{1/5}[(Q/Q_0{)}^2-c{]}^{1/5}$, and substitute into $Q^{\prime }=2\alpha \sqrt{M}$. The coefficient satisfies $(2\alpha /Q_0{)}^5[5B_0{Q}_0^2/(8\alpha )]=20{\alpha }^4{B}_0/Q_0^3$, givingWith positive source fluxes, $c<1$. For $c<0$ the source has excess momentum relative to pure plume balance, so it is a momentum-dominated forced plume. For $c=0$ it is already in pure plume balance, and the exact solution is the pure similarity translated to a finite plume virtual origin. For $0<c<1$ the source has deficient momentum relative to its buoyancy and volume fluxes, a lazy plume. Here “forced” in the broad sense merely means finite source fluxes; this sign classification distinguishes its jet-like and lazy branches. At large height, finite $c{Q}_0^2/Q^2$ vanishes and every positive branch approaches pure plume balance.

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.

The Batchelor entrainment hypothesis makes the inward entrainment velocity proportional to the local plume speed:where $\alpha >0$ is the dimensionless entrainment coefficient. For an upward plume $w>0$, this is $u_e=\alpha w$. The hypothesis represents turbulent ingestion of ambient fluid; it is a closure rather than a consequence of mass conservation.

Let $A=b^2$, suppressing the common factor $\pi$ from cross-sectional integrals. In the Boussinesq approximation, use a constant reference density ${\rho }_r$ in the leading-order mass conservation equation. Dividing that equation by ${\rho }_r$ gives volume conservation:The factor two is the ratio of the perimeter $2\pi b$ to the suppressed area factor $\pi$. Density differences are retained in buoyancy, even though they are neglected in the leading inertial and volume terms.

To derive the density-deficit balance, first leave the small density differences explicit in the supplied mass equation. Because the ambient density ${\rho }_0(z)$ is time independent, the product rule givesThus the ambient mass entrained through the perimeter cancels exactly against its contribution to the volume budget, yielding the required buoyancy balance. It is understood to the accuracy of the underlying Boussinesq approximation; the volume equation is not an additional exact constraint on a fully variable-density fluid.

For later use define reduced gravity $g^{\prime }=g({\rho }_0-\rho )/{\rho }_r$ and the squared buoyancy frequency $N^2=-g{\rho }_{0z}/{\rho }_r$. With $B=A{g}^{\prime }$ and $Q=Aw$, the balance isA stably stratified ambient has $N^2>0$, so upward motion consumes the plume's positive density deficit; an unstable ambient reverses this sign.

Use pi-normalized kinematic plume fluxes: physical volume, momentum and buoyancy fluxes are $\pi Q$, $\pi M$ and $\pi F$. For a top-hat plume model in a homogeneous ambient, constant entrainment coefficient $\alpha$ givesEliminating height between the first two equations yields $d(M^{5/2})/dQ=5FQ/(4\alpha )$. Hence the plume flux-balance invariant is $M^{5/2}-5F{Q}^2/(8\alpha )$. Pure plume balance means that this invariant vanishes, including at a finite-volume source:It does not mean that $Q_s$ and $M_s$ must individually be zero. The balanced solution has a plume virtual origin below the physical source:where $z_v=[Q_s/(C{F}_s^{1/3}){]}^{3/5}>0$. For given nonzero source fluxes, the displayed balance is a compatibility condition; they cannot be chosen independently and still describe an exactly balanced source.

The entrainment coefficient $\alpha$ is the ratio of the inward ambient-fluid velocity normal to a buoyant plume edge to the characteristic buoyant plume velocity. For an upward buoyant plume, take $u_e=\alpha w>0$. There are two edges at half-width $b$, so the full-section volume source is $2u_e$. Dividing every integral balance by two gives consistent half-section equations.

For top-hat plume profiles, volume conservation and mass conservation areThe mass entering the buoyant plume is ambient mass with local mass density ${\rho }_0(z)$. The Boussinesq approximation uses a common reference mass density in inertia while retaining the mass density deficit in buoyancy. The mass equation may be kept to the first order in that deficit to derive its transport.

Multiply the volume equation by the time-independent ambient mass density. The product rule givesSubtract the mass equation and multiply by $g$. The entrainment sources cancel exactly, leavingwhereFor stable stratification, ${\rho }_0^{\prime }<0$ and the buoyancy frequency $N$ is real. Rising fluid moves into lighter ambient fluid and loses its mass density deficit relative to that local ambient. Thus environmental stratification reduces the buoyant plume's buoyancy flux, even though the ambient fluid entering by entrainment contributes no mass density deficit at its own height.

The Boussinesq approximation retains small density differences only in the gravitational buoyancy force and uses a constant reference density in inertia and volume conservation. The entrainment coefficient $\alpha$ specifies the inward velocity of ambient fluid relative to the thermal boundary, $u_e=\alpha W$. Turbulent interface motions engulf ambient water and mix it into the coherent rising volume.

For a spherical buoyant thermal, $V=4\pi a^3/3$. To leading order in the dilute gas fraction, expansion of the gas does not dominate the total-volume balance, so fluid entrainment givesBut $dV/dt=4\pi a^{2}da/dt$ and $dz/dt=W$. Therefore, wherever the thermal is rising,This growth law for a dilute bubbly thermal is independent of stratification under the constant-$\alpha$ model. Stratification changes the speed and eventual stopping height; it does not enter this volume-per-distance relation. In real flows the entrainment coefficient itself can vary with stratification and flow structure.

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.

Assume the ambient fluid has zero downslope velocity, no bed drag, and an entrainment coefficient $\alpha >0$. The entraining shallow-water layer receives volume at rate $\alpha |u|$, but neither excess buoyancy nor downslope momentum from that ambient fluid. Consequently volume conservation, conservation of excess buoyancy, and momentum conservation giveThe pressure contribution is the integral of the hydrostatic pressure from part (a). The last source is the downslope component of the gravitational force.

With $D={\partial }_t+u{\partial }_x$, the equivalent equations areThe last term is the momentum required to accelerate newly entrained fluid. For variables $(h,u,b)$, the matrix of the spatial derivatives isThus it is a strictly hyperbolic system for $h,b>0$ and $0\le \theta <\pi /2$. At $\theta =\pi /2$ all three characteristic speeds coincide and the matrix is not diagonalizable. The exclusion of that angle presumes the physical range with $\cos \theta \ge 0$: allowing an overhanging channel with $\cos \theta <0$ would also destroy hyperbolicity.