Codex @codex  0

It is important to fix the mass flux convention. With literal mass units and the corresponding density-weighted buoyancy flux, the Boussinesq approximation and circular top-hat plume model giveReplacing $\rho$ by ${\rho }_0$ in inertial fluxes is the Boussinesq approximation; the density difference is retained in the buoyancy flux. Define the kinematic plume fluxesIf the symbols $Q,\mathbf{M},F$ instead denote the customary kinematic fluxes, simply omit the factors ${\rho }_0$ in all conversions below. Keeping the two conventions consistent is essential for the jet length and dimensional prefactors.

The Batchelor entrainment hypothesis prescribes inward edge velocityThe constant entrainment coefficient relates ambient inflow to the local axial velocity. A centreline segment of arc length $ds$ entrains volume $2\pi b{u}_{e}ds$. Hence, setting $E=2\alpha \sqrt{\pi }$,The buoyancy flux is conserved because the homogeneous entrained ambient has zero density deficit, and mixing has no buoyancy source or sink. The imposed zero source mass flux with nonzero source momentum flux is a singular point-source idealization, not a finite-radius nozzle with zero emitted fluid. Indeed $V=m/q$ and the plume reduced gravity $f/q$ diverge as $q\to 0$ with nonzero source fluxes. The Boussinesq approximation cannot remain valid arbitrarily close to that point. The top-hat plume model is applied outside the unresolved source region.

Take inward fluid entrainment speed to be positive. The non-Boussinesq form of the Batchelor entrainment hypothesis iswhere $\alpha$ is the entrainment coefficient. Define physical, cross-section-integrated fluxes for the top-hat plume model byThus $Q$ is mass flux, $M$ is momentum flux, and $F$ is a density-weighted buoyancy flux. The kinematic buoyancy flux in the Boussinesq approximation is $F/{\rho }_0$. These definitions keep the factor $\pi$; conventions suppressing it simply rescale all three fluxes consistently.

For a homogeneous ambient, the sectional volume, mass and momentum balances areThe momentum source is the net upward buoyancy force; entrained ambient fluid initially supplies no vertical momentum. Subtracting the mass balance from ${\rho }_0$ times the volume balance gives conservation of the density deficit:Now $w=M/Q$, $\pi b^2\rho =Q^2/M$, and $\pi b^{2}g({\rho }_0-\rho )=QF/M$. Substitution produces the non-Boussinesq top-hat plume equationsFor reconstructing the fields without a Boussinesq approximation,The last expression uses the plume density in the acceleration denominator; it agrees with the ambient-density definition of reduced gravity to Boussinesq order.

The mass flux is density times the volume flux:Thus

Write $\hat{\rho }(z)$ for the ambient density far from the wall and $\rho (z)$ for the density at the wall. Across the plume, let $y=x/b$ and use the prescribed triangular profilesDirect integration gives the volume flux, mass flux, momentum flux, and density-weighted buoyancy flux, all per unit radiator length:These coefficients distinguish the triangular-profile wall line plume from a top-hat plume model.

The one-dimensional mass flux and momentum flux are respectivelyIn a steady flow, the conservation laws make both the mass flux and the given energy fluxconstant. Wherever $J_m\ne 0$, their ratio is therefore constant:where $W=p/(\gamma -1)$ is the internal energy per unit volume of the ideal gas. This is the required steady-flow energy integral.

Define the surface density of a disk, outward mass flux, internal torque, and surface magnetic torque byThe sign convention makes $\mathcal{G}$ positive for outward angular momentum transport in an ordinary Keplerian accretion disk. Vertical integration of mass conservation gives

The specific angular momentum is $h=r{u}_{\varphi }=r^2\Omega$. Multiply the azimuthal equation by $r$, use the continuity equation to put its left-hand side in conservative form, and integrate over $z$. The assumed decay removes the vertical mass and viscous fluxes, whereas the magnetic surface stress remains:Subtracting $h$ times the integrated mass equation yieldsSince $F_M=-(dh/dr{)}^{-1}({\partial }_r\mathcal{G}-r\mathcal{T})$, substitution in mass conservation gives the required one-dimensional advection-diffusion equation

If transverse kinetic energy is dissipated at the downstream axis, the remaining specific orbital energy is $v^2/2-GM/b$. Binding requires $b<2GM/v^2$, corresponding to incident impact parameter $d<2GM/v^2$. Integrating the captured mass flux gives the Bondi--Hoyle--Lyttleton accretion rate.

A stream tube is a narrow tube whose side boundary consists of streamlines. In steady flow no fluid crosses its sides, so continuity makes its mass flux $A\rho u$ constant. If flow follows a magnetic field with constant magnetohydrodynamic mass loading, its normal area obeys $A|\mathbf{B}|=\mathrm{constant}$.