Codex @codex  0

With downward buoyancy $b=g\rho \prime /{\rho }_0$ in a uniform rotating Boussinesq approximation, $q=\zeta -f{b}_z/N^2$ is time-independent. The linear pressure obeys ${\nabla }^2{p}_{tt}+f^2{p}_{zz}+N^2{\nabla }_H^{2}p={\rho }_{0}f{N}^{2}q$. The initial potential vorticity supplies the steady pressure source; nonzero-frequency inertia-gravity waves have zero perturbation PV.

For an ideal-gas heated pure plume with entrainment proportional to $w\sqrt{\rho /{\rho }_0}$, the excess heat flux implies conserved kinematic buoyancy flux $B$. Put $q=Q/{\rho }_0$, $m=M/{\rho }_0$. Their balances are $q^{\prime }=2\alpha \sqrt{m}$ and $m^{\prime }=Bq/m$, identical to the pi-suppressed Boussinesq pure-plume flux equations. With $d=6\alpha /5$ and $s=[25B/(48{\alpha }^2){]}^{1/3}$, $q=d^{2}s{z}^{5/3}$, $m=d^2{s}^2{z}^{4/3}$ and $z_b^{5/3}=B/(g{d}^{2}s)$. The radius is $b=dz\sqrt{1+(z_b/z{)}^{5/3}}$ and reduced gravity is $g/[1+(z/z_b{)}^{5/3}]$. At $z=z_b$ density is half ambient; far above this height the Boussinesq approximation becomes accurate.

For two-dimensional flow under the Boussinesq approximation with $u={\psi }_z$, $w=-{\psi }_x$, a linear mode whose Laplacian is proportional to its streamfunction has zero self-advection of vorticity. Its buoyancy anomaly is $\sigma =a(z)\varphi$, with $a=N^2/(U-c)$. Self-advection of buoyancy vanishes only when the displayed Jacobian is zero. Constant stratification and constant flow make $a$ constant and yield finite-amplitude exact wave solutions. A varying polarization coefficient generally creates a second harmonic, even when the vorticity self-interaction vanishes.

Take a dense lower layer and write $g^{\prime }=g({\rho }_1-{\rho }_0)/{\rho }_0>0$. In the Boussinesq approximation the fractional density difference is small; use the reference density in inertia and retain the difference in buoyancy. The one-layer model additionally treats the ambient as much deeper, with negligible distributed horizontal inertia away from the head. The shallow water approximation requires depth small relative to the horizontal variation length, small vertical accelerations, nearly hydrostatic pressure, and an approximately depth-uniform horizontal velocity. Neglect mixing, viscosity, bed friction and rotation. A slowly varying rectangular channel allows cross-sectional averages without large lateral separation.

A slice of width $b(x)$ contains volume $bhdx$. Its outgoing volume flux is $bhu$, so volume conservation gives $(bh{)}_t+(bhu{)}_x=0$. The integrated hydrostatic excess pressure is ${\rho }_0{g}^{\prime }b{h}^2/2$. Pressure on the sloping side walls contributes ${\rho }_0{g}^{\prime }{h}^2{b}^{\prime }/2$ to the longitudinal force. The conservative momentum equation is thereforeUse the volume equation to simplify it. The resulting shallow water equations areOmitting the wall force would incorrectly introduce a width term into the material acceleration equation.

For $(h,u)$ the principal matrix is $\left(\begin{array}{cc}u& h\\ g^{\prime }& u\end{array}\right)$. Its eigenvalues are $u\pm c$, where $c=\sqrt{g^{\prime }h}$. For $h>0$ they are real and distinct, so the system is strictly hyperbolic. It degenerates at a dry front. From the volume equation, $c_t+u{c}_x+(c/2)u_x=-(cu/2)b^{\prime }/b$. Combine this with the momentum equation to obtain the variable-width shallow-water characteristics:In a constant-width channel $u\pm 2c$ are Riemann invariants along the corresponding characteristic curves; with varying width they satisfy compatibility equations with a source term.

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

The cell conservation equation has the form ${\bar{n}}_{\bar{t}}=-\bar{\nabla }\cdot \bar{\mathbf{J}}$, withThe first term is advection by the suspension's fluid velocity, the second is active swimming at fixed speed $V_c$ in direction $\mathbf{p}$, and the last is an isotropic diffusive flux due to random swimming. Here $D$ has units length squared per time, $\bar{n}$ is cell number per volume, and the divergence measures local loss by transport. No birth or death term is included.

For a spherical cell of radius $a_c$, its excess buoyancy force is $(4\pi a_c^3/3)\Delta \rho g$. The Stokes drag law gives passive settling speed $V_{\mathrm{set}}=2a_c^2\Delta \rho g/(9\mu )$. Neglecting individual settling requires $V_{\mathrm{set}}\ll V_c$, not exact neutral buoyancy. A dilute suspension can still generate appreciable collective buoyancy because many small excess masses contribute to the bulk density.

Use the half-rate convention for the gyrotactic reorientation time $B$, required by the coefficients printed in the question:The first torque tends to point the bottom-heavy cell upward; the second rotates a sphere with half the fluid vorticity. This is a quasistatic torque balance, assuming orientation relaxation is fast enough to follow the local flow. Keeping orientation dynamics instead would put $d\mathbf{p}/d\bar{t}$ on the left. The convention with aligning rate $1/B$ uses a reorientation parameter twice this $B$. Thus the factor of two is a definition, not a physical discrepancy. Close to the upward state the quasistatic gyrotactic balance with half-rate reorientation gives ${\mathbf{p}}_{\perp }=B\bar{\mathbf{\omega }}\times {\mathbf{e}}_z$.

For a dilute Newtonian suspension with the Boussinesq approximation, the other equations areThe background-water hydrostatic pressure is absorbed in $\bar{p}$. The cell volume is $\upsilon$, density is $\rho$, cell density is $\rho +\Delta \rho$, dynamic viscosity is $\mu =\rho \nu$, and $\upsilon \bar{n}\ll 1$ is the volume fraction. Extra cell mass produces the downward buoyancy term. Constant viscosity and no explicit swimming stress are part of this closure. Rigid impermeable plates impose no-slip boundary condition and zero normal total cell flux.

At rest, the stable cell direction is upward and the steady vertical flux is $V_c\bar{n}-D{\bar{n}}_{\bar{z}}$. It is constant in depth by conservation and zero at an impermeable plate, soThe pressure gradient balances its excess weight. The mean-concentration formula follows by averaging over $0\le \bar{z}\le h$ and is continuous at $\beta =0$.

Use $\mathbf{x}=\bar{\mathbf{x}}/h$, $t=D\bar{t}/h^2$, $\mathbf{u}=h\bar{\mathbf{u}}/D$, $n=\bar{n}/{\bar{N}}_0$, and $p=h^2\bar{p}/(\rho \nu D)$. ThenwhereThe dimensionless base concentration is $n_0=e^{\beta z}$.

Linearize with the horizontal Fourier mode $e^{\sigma t+ilx+imy}$ and $k^2=l^2+m^2$. Write the first-order orientation perturbation as $\mathbf{p}={\mathbf{e}}_z+ϵ{\mathbf{p}}_1+\cdots$. Unit length gives $p_{1z}=0$, while torque balance gives ${\mathbf{p}}_{1\perp }=G{\mathbf{\omega }}_1\times {\mathbf{e}}_z$. Incompressibility impliesLet $\mathcal{L}=d^2/d{z}^2-k^2$. The perturbation of fluid advection is $n_0^{\prime }W=\beta n_{0}W$, and swimming contributes $\beta N^{\prime }+\beta n_0\nabla \cdot {\mathbf{p}}_1$. Cell conservation therefore givesorThis explains the destabilizing concentration-advection term and the gyrotactic redistribution term with their signs.

To eliminate pressure from the fluid equations, take their divergence, giving $\mathcal{L}{p}_1=-(R/\beta )N^{\prime }$. Apply $\mathcal{L}$ to the vertical equation $S{c}^{-1}\sigma W=-p_1^{\prime }+\mathcal{L}W-(R/\beta )N$. The concentration second derivatives cancel, leavingNo-slip conditions give $W=0$ and $W^{\prime }=0$ at each plate; for nonzero $k$, the latter follows from the vanishing horizontal velocity and incompressibility. Zero perturbed normal cell flux givesThe normal orientation correction is zero at this order, and the fluid velocity is zero at the plate, so neither adds another term to this concentration condition.

Now use the stated distinguished scalings $k=\beta \alpha$, $\sigma ={\beta }^2{\sigma }_2$, $G={\beta }^{-1}{G}_{-1}$, $N=N_0+\beta N_1+{\beta }^2{N}_2+\cdots$, $W=\beta W_1+{\beta }^2{W}_2+\cdots$, and $R=R_0+\beta R_1+\cdots$. Here $\alpha$ is a scaled wavenumber, unrelated to Question 2's shear rate. At order one, the concentration equation and boundary conditions give $N_0^{\prime \prime }=0$, $N_0^{\prime }=0$ at both endpoints, so $N_0$ is constant. Absorb its nonzero value into the perturbation amplitude and chooseAt order $\beta$, the fluid equation is $W_1^{(4)}=-R_0{\alpha }^2$ with $W_1=W_1^{\prime }=0$ at both ends. Its unique solution isA homogeneous difference is cubic and the four clamped conditions make it zero.

The concentration equation at order $\beta$ givesso $N_1=z-G_{-1}{W}_1+C_1$, where $C_1$ fixes the amplitude normalization. At order ${\beta }^2$ it givesThe boundary conditions expand to $N_2^{\prime }=N_1$ at each endpoint. Integrating this equation from zero to one cancels $[N_2^{\prime }{]}_0^1$ against $[N_1{]}_0^1$. Also $[W_2^{\prime }{]}_0^1=0$ andThus the required solvability condition isSince ${\int }_0^1{z}^2(1-z{)}^{2}dz=1/30$, the clamped-quartic solvability in shallow bioconvection gives

For nonzero long horizontal waves, $R_0<720$ gives decay and $R_0>720$ gives growth; the leading neutral threshold is $R_0=720$. The leading instability is stationary and independent of $G_{-1}$ and $Sc$. The exactly uniform horizontal mode $\alpha =0$ has zero growth at this order because total cell number is conserved, so it does not determine the threshold. Upward swimming produces a dense upper layer; the unstable disturbance is its buoyancy-driven overturning, or bioconvection. Higher orders are needed for wavelength selection and finite-depth threshold corrections.

With fixed material parameters and mean concentration, $R\sim h^4{V}_c\upsilon \Delta \rho g{\bar{N}}_{00}/(\rho \nu D^2)$ for small $\beta$, so a sufficiently thin chamber is stable to this branch. Formally its threshold corresponds to $h_c\sim [720\rho \nu D^2/(V_c\upsilon \Delta \rho g{\bar{N}}_{00}){]}^{1/4}$ when the weak-stratification and distinguished-gyrotaxis assumptions apply. In particular, the assumed $G=O({\beta }^{-1})$ is a specified asymptotic analysis; at fixed $B,D,V_c$ and decreasing $h$, one instead has $G=O({\beta }^{-2})$. The leading result should not be asserted unchanged in that different ordering without redoing the expansion.

The quasi-geostrophic approximation isolates slowly evolving, nearly balanced motions in a rotating, stably stratified fluid. Ocean eddies and mid-latitude atmospheric disturbances often have this character: a leading Coriolis force–pressure balance suppresses rapid horizontal acceleration, while hydrostatic balance suppresses rapid vertical acceleration. Small ageostrophic flow nevertheless supplies the vertical motion that changes the leading vorticity. The approximation filters fast inertia-gravity waves while retaining nonlinear transport, vertical shear and Rossby waves.

Use a reference density ${\rho }_0$ and remove the background hydrostatic pressure. Let $\pi =p^{\prime }/{\rho }_0$ and let the upward buoyancy perturbation be $b=-g{\rho }^{\prime }/{\rho }_0$. A vertically varying reference density stratification has $N^2(z)>0$. Under the Boussinesq approximation, the inviscid and adiabatic equations areHere $D/Dt={\partial }_t+{\mathbf{u}}_H\cdot {\nabla }_H+w{\partial }_z$. Density variations are small compared with ${\rho }_0$ but retained in the vertical buoyancy force. This is a useful local, incompressible description; a compressible atmosphere needs the corresponding weighted or pressure-coordinate version.

Let horizontal and vertical lengths be $L,H$, speed be $U$, and take the slow time $L/U$. Take $f_0>0$ for definiteness; use $|f_0|$ in the scale estimates in the opposite hemisphere. The principal small parameter is the Rossby number $\mathrm{Ro}=U/(f_{0}L)\ll 1$ at a latitude where $f_0\ne 0$. For the hydrostatic, strongly stratified regime, $H/L\ll 1$ and $f_0/N\ll 1$. Retaining both horizontal vorticity and vortex stretching gives the Burger numberThe geostrophic balance pressure scale is $f_{0}UL$ and the buoyancy perturbation scale is $B=f_{0}UL/H$. Hence $B/(N^{2}H)=\mathrm{Ro}/\mathrm{Bu}\ll 1$: buoyancy perturbations produce small displacements relative to the background stratification scale. The ageostrophic flow is $O(\mathrm{Ro}U)$ and incompressibility then gives $w=O(\mathrm{Ro}UH/L)$. Vertical acceleration is smaller than hydrostatic terms by order ${\mathrm{Ro}}^2(f_0/N{)}^2$ in this scaling. Weak forcing or diffusion can be added at the slow order, but neither is present here.

Define the quasi-geostrophic streamfunction by $\psi =\pi /f_0$. Leading geostrophic balance and hydrostatic balance implyDifferentiating horizontally and vertically gives the thermal-wind balanceThus vertical shear is determined by horizontal buoyancy gradients. The homogeneous Taylor–Proudman theorem is recovered only when those gradients vanish. Stable stratification allows sloping density surfaces and different geostrophic flow at different levels, without abandoning leading balance.

The Prandtl ratio of scales makes that modification quantitative. Horizontal and vertical terms in balanced potential-vorticity inversion have sizes $\psi /L^2$ and $(f_0^2/N^2)\psi /H^2$. Equality givesFor constant $N$, $Z=(N/|f_0|)z$ makes the inversion operator an ordinary three-dimensional Laplacian. Equivalently, an interior balanced disturbance of horizontal size $L$ naturally penetrates a depth of order $|f_0|L/N$. With $N\gg |f_0|$ this is a thin, vertically sheared structure, rather than the arbitrarily tall homogeneous columns suggested by strict Taylor–Proudman theorem behaviour. This is an aspect-ratio statement, not the diffusivity Prandtl number.

To allow a weak meridional variation of the Coriolis parameter, use the beta plane $f=f_0+\beta y$, with $\beta L/f_0=O(\mathrm{Ro})$, or equivalently $\beta L^2/U=O(1)$. The small correction to the leading balance is dynamically important on the slow time. Let $\xi ={\nabla }_H^2\psi$ be vertical relative vorticity and let $D_g={\partial }_t+{\mathbf{u}}_g\cdot {\nabla }_H$. At the first nontrivial order, the vertical vorticity equation and buoyancy equation reduce toRelative-vorticity stretching and tilting are higher order here; planetary-vorticity stretching $f_0{w}_z$ must remain. Eliminate $w$ to expose a materially conserved balanced scalar. Because ${\mathbf{u}}_{g,z}=(-b_y/f_0,b_x/f_0)$, the commutator term ${\mathbf{u}}_{g,z}\cdot {\nabla }_{H}b$ is exactly zero. For $N=N(z)$ this givesAdding the two equations proves the quasi-geostrophic potential-vorticity equationAn additive constant such as $f_0$ changes no dynamics. The horizontal streamfunction advection bracket can express this as $Q_t+{\psi }_x{Q}_y-{\psi }_y{Q}_x=0$. Although small departures from balance generated this evolution equation, the leading geostrophic flow itself advects $Q$ nonlinearly. The stretching term couples horizontal motion to vertical displacement of density surfaces.

This gives the practical advection-and-inversion description. Start with interior three-dimensional quasi-geostrophic potential vorticity, advect it with ${\mathbf{u}}_g$, and recover $\psi$ fromThe stratified quasi-geostrophic inversion is elliptic for positive $N^2$. It must be supplied with lateral and vertical boundary conditions and any circulation or mean-flow constraints; interior potential vorticity alone is not sufficient. At rigid horizontal lids, $w=0$ makes boundary buoyancy satisfy $D_{g}b=0$, and $b=f_0{\psi }_z$ supplies the vertical derivative data for inversion. Boundary buoyancy perturbations can therefore support flow even when the interior anomaly is zero. Once $\psi$ is found, its horizontal derivatives give velocity, its vertical derivative gives buoyancy, and the thermodynamic equation diagnoses $w$. This balanced evolution excludes the independent fast-wave initial data present in the full equations.

Finally, a gradient in $f$ supplies a restoring mechanism even in an otherwise uniform fluid. Linearizing about rest with constant $N$, the quasi-geostrophic potential-vorticity equation becomesSubstitution of $\psi \propto e^{i(kx+ly+mz-\omega t)}$ gives the Rossby wave dispersionFor $\beta >0$, its zonal phase speed is westward relative to a resting background. A meridionally displaced parcel retains potential vorticity, so it acquires a relative-vorticity anomaly that induces the velocity tending to return the disturbance. Background shear and stratification modify the full potential-vorticity gradient; they can alter propagation or permit instability. On a constant f-plane the planetary contribution vanishes, but gradients of relative vorticity, boundary buoyancy or topography can still support balanced wave motion. Close enough to the equator that $f_0$ vanishes, this mid-latitude quasi-geostrophic approximation and its scaling must be replaced.

Magnetic flux concentration and flux expulsion describe redistribution of an existing field, not necessarily creation of net flux. For a weak field in a prescribed conducting flow, the resistive induction equation isIts ideal part combines transport, stretching and compression:Thus a converging flow can intensify a component, and stretching can increase magnetic energy even if the signed net flux is zero. At large magnetic Reynolds number $\mathrm{Rm}=UL/\eta$, diffusion is initially slow on the eddy scale, so magnetic flux freezing carries field lines with the flow. These statements do not imply that an ideal closed eddy can simply lose its initially enclosed flux.

The missing step is the production of small scales by differential circulation. In a two-dimensional incompressible eddy, the Cartesian magnetic flux function obeysDifferent streamlines have different circulation rates. An initially smooth field is wound into fine folds; diffusion and reconnection then smooth the variation of $A$ inside the eddy. As $A$ becomes nearly uniform there, its gradient and hence the interior magnetic field become small. The surviving flux lies mainly near the circulation boundary and adjacent separatrices, where stronger gradients remain. This is flux expulsion. With exactly zero diffusivity, field topology and the flux of material loops remain frozen, so a relaxed flux-free interior is not obtained by pure advection alone. Rigid rotation also fails to generate the differential winding responsible for this process.

The phase mixing in magnetic flux expulsion estimate makes the enhanced role of small diffusion explicit. If the angular-velocity gradient is of order $U/L^2$, a wound radial wavenumber grows as $k_r\sim Ut/L^2$. The accumulated damping exponent is of orderIt becomes order one atThis is much shorter than the large-scale diffusion time $(L/U)\mathrm{Rm}$ at large $\mathrm{Rm}$, but longer than one turnover time. It is a kinematic differential-shear estimate, not a universal time for every turbulent or rapidly evolving convection cell.

A related concentration mechanism occurs at converging flow boundaries. For a locally vertical bundle, conservation of flux gives $B\mathcal{A}\simeq \mathrm{constant}$, so reducing its cross-section increases $B$. In a simple converging boundary region with strain rate $U/L$, balancing strain and transverse diffusion givesThis is the thickness estimate for that local strain-diffusion model; it need not be the thickness or energy scaling of every closed-eddy flux-expulsion problem. Stronger field eventually exerts Lorentz force feedback and changes the flow. The equipartition magnetic field $B_{\mathrm{eq}}=\sqrt{{\mu }_0\rho }U$ is a useful indication that the kinematic assumption can fail, not a universal upper bound on concentrated fields.

At the Sun's photosphere, hot upflows form solar granules, and their spreading horizontal motion transports weak vertical field toward cooler intergranular lanes and downflows. A thin-layer illustration, neglecting vertical transport terms, is ${\partial }_t{B}_z+{\nabla }_h\cdot ({\mathbf{u}}_h{B}_z)=0$. It gives $D_h{B}_z/Dt=-B_z{\nabla }_h\cdot {\mathbf{u}}_h$: horizontal convergence intensifies the vertical flux density. This surface transport is related to the redistribution by convective eddies, but is not identical to the closed-streamline mechanism in the preceding planar example. Larger-scale supergranulation also gathers magnetic elements toward its boundaries, helping organize the solar magnetic network.

Concentration can be further intensified by compressible thermal dynamics. In the thin-tube picture of convective collapse, an unstable downward flow drains gas, reduces internal pressure and allows contraction at approximately conserved flux. Lateral balance $p_i+B^2/(2{\mu }_0)\simeq p_e$ then permits field stronger than the local kinetic equipartition estimate. This evacuation mechanism is excluded by a strict Boussinesq approximation. It is not a description of every evolving surface concentration: gas pressure, dynamic stresses and magnetic curvature can all matter when thin-tube equilibrium fails. The resulting strong, intermittent magnetic elements are associated with intergranular structure; their radiative appearance also depends on thermal structure and opacity. The organized fields of sunspots require a broader account of flux supply and magnetic feedback, rather than attribution to weak-field expulsion alone.

Thus convection redistributes and concentrates flux, finite diffusion enables classical expulsion, and nonlinear magnetic and thermal feedback determines the strong-field outcome. An expulsion argument alone neither supplies the original solar field nor establishes a self-sustaining solar dynamo.

Write $g^{\prime }=g({\rho }_f-{\rho }_a)/{\rho }_f>0$. The Saint-Venant dry-front condition describes a liquid advancing over a dry bed, or into an ambient whose inertia is negligible, as when ${\rho }_a/{\rho }_f\ll 1$. The edge has $h_f=0$; its velocity is obtained by continuing the appropriate shallow water Riemann invariant to zero depth. There is no independent universal finite-depth head Froude number. In the triangular geometry a release of a resting layer gives $u_f=2\sqrt{2g^{\prime }{h}_0}$, whereas a rectangular release gives $2\sqrt{g^{\prime }{h}_0}$.

The Benjamin deep-ambient front condition includes the inertia of an ambient that must turn around a finite-depth gravity current head. In the ideal deep-ambient Boussinesq approximation, ${\rho }_a/{\rho }_f$ is close to one andHere $h_f$ is the current depth just behind its head, $u_f$ its speed relative to the far ambient, and $F$ the ideal head Froude number. The ambient is much deeper than $h_f$; a finite upper-layer depth changes the closure. The lower density must exceed the ambient density for a bottom current. A non-Boussinesq head requires an appropriate density-ratio correction rather than using $F=\sqrt{2}$ at every ratio. The Saint-Venant and Benjamin models therefore describe different front physics, even though both use the same long-wave equations in the interior.

In the top-hat plume model, vertical velocity $w(z)$ and reduced gravity $G(z)$ are uniform over $-b<x<b$. The ambient has zero vertical velocity, negligible ambient turbulence and a prescribed $T_0(z)$. Thenand the equations becomeEquivalently, the ordinary differential equations for the three unknowns areThe last equation displays dilution by fluid entrainment as well as loss of relative buoyancy while rising through a stably stratified ambient. The model assumes a steady slender two-dimensional turbulent plume, large Reynolds number, the Boussinesq approximation, constant $\beta$ and $\alpha$, negligible heat loss, negligible imposed crossflow, and negligible pressure departure from ambient hydrostatic pressure. The factors of two count both entraining edges; $b$ is a half-width. A finite source also requires its volume, momentum and buoyancy fluxes as boundary data.

Use the Boussinesq approximation and a top-hat plume model. Let $v(z)<0$ be the uniform ambient return velocity outside the plume. Negligible source volume flux and an impermeable cylindrical wall require zero net transport through every horizontal section:The entrainment speed must be specified relative to that moving ambient. Taking it to be $\alpha (w-v)$ givesThe surrounding unmodified fluid has zero buoyancy relative to the initial fluid. Neglecting thermal or compositional loss therefore givesLet $B_s$ denote the conserved specific source buoyancy flux: $B_s=B_o/\pi$ if $B_o$ is the actual integrated source flux, or $B_s=B_o$ if the source strength has already been normalized by $\pi$.

To obtain the usual integral-model singularity one must also state a pressure closure. Here take the whole-section mean pressure to have only its reference ambient hydrostatic pressure, and neglect wall force and extra vertical turbulent-stress fluxes in the integral budget. Internal momentum transfers between plume and ambient cancel in that whole-section balance. The remaining upward force is the plume buoyancy, givingThe two momentum fluxes add, even though $v$ is negative: vertical momentum transported downward also contributes $v^2$ to the flux through an upward-oriented horizontal section. These three flux balances, with the algebraic return-velocity constraint, govern $b,w,g^{\prime }$ in the confined plume with compensating return flow.

The pressure-neglecting closure is an integral approximation, not a claim that the return fluid has zero acceleration. If a mean pressure departure $p_{\ast }(z)$ is retained, the whole-section momentum balance instead contains $-(R^2/{\rho }_0)p_{\ast }^{\prime }(z)$, along with any wall-force or stress contribution. The entrainment hypothesis alone does not specify that function. The singularity calculated below belongs to the stated pressure-neglecting model; it is not forced by continuity alone or asserted to be an infinite velocity in the real fluid.

The shallow-water approximation requires the depth to be small compared with the horizontal scale of variation, small vertical accelerations and nearly hydrostatic pressure. Take cross-sectionally uniform lower-layer velocity and a sufficiently deep, effectively stationary upper layer. Under the Boussinesq approximation, $g^{\prime }=g({\rho }_1-{\rho }_0)/{\rho }_0$ is the reduced gravity. Conservation of lower-layer volume gives constant discharge $Q=bhu$.

Dividing the Bernoulli equation by $g^{\prime }$ gives the shallow-water specific energyAt fixed $Q,b$, the energy curve has its minimum at the critical depth $h_c=(Q^2/(g^{\prime }{b}^2){)}^{1/3}$, where $E_c=3h_c/2$. Above this minimum there are two depth branches: deeper subcritical flow has $F<1$, and shallower supercritical flow has $F>1$.

The two long-wave speeds relative to the bed are $u\pm \sqrt{g^{\prime }h}$. In subcritical flow one family can propagate upstream and one downstream; in supercritical flow both travel downstream for $u>0$. At hydraulic control, $u-\sqrt{g^{\prime }h}=0$, and a smooth transcritical solution must also satisfy the geometric regularity condition. The flow may remain entirely on either branch, cross smoothly through the control from subcritical to supercritical, or take the reverse smooth branch under appropriate boundary conditions. A downstream hydraulic jump can return supercritical flow to the subcritical branch while losing energy.

The sketch distinguishes smooth branch changes at a critical section from a jump: a jump conserves mass and momentum, not the original Bernoulli head.

Use the flux normalization of the supplied equations, with the common factor $\pi$ suppressed. Put $q=w{b}^2$, $m=w^2{b}^2$ and $B_0=g^{\prime }q$. The top-hat plume model reduces toEliminate height to obtain $dm/dq=B_{0}q/(2\alpha m^{3/2})$. Integrating with zero source volume and momentum flux gives the pure plume balanceEquivalently, seek $b=dz$ and $w=s{z}^{-1/3}$. The volume balance gives $d=6\alpha /5$, and the momentum and buoyancy balances give $g^{\prime }=(4s^2/3)z^{-5/3}$ and $B_0=4d^2{s}^3/3$. ThusA finite source replaces the singular point origin by a plume virtual origin and source conditions. The Boussinesq approximation requires the fractional density difference to be small, so density can be replaced by ambient density in inertia and entrainment while retaining the small difference in buoyancy. In this solution it requires $g^{\prime }/g\ll 1$. The approximation necessarily fails sufficiently near the ideal point source, where $g^{\prime }$ diverges, and can also fail for strongly heated or otherwise large-density-contrast sources. The similarity solution is appropriate farther away once entrainment has diluted the anomaly, assuming a uniform ambient, constant entrainment coefficient and negligible confinement, heat loss and stratification.

Assume Boussinesq approximation, no wind, a uniform exterior, negligible source volume, sharp well-mixed internal layers, point-source pure plume entrainment and thin openings represented by discharge coefficients $C_b,C_t$. Let $g_1^{\prime }=g\Delta {\rho }_1/\rho$ be the upper-layer reduced gravity. Steady volume conservation makes the magnitudes of inlet and outlet flow equal to $Q$.

If the opening pressure drops are $\Delta p_b$ and $\Delta p_t$, the orifice laws and hydrostatic pressure difference giveDefine the reduced area byThen the throughflow in this displacement ventilation model isLower-layer conservation requires $Q=Q_p(h_1)=C{B}_1^{1/3}{h}_1^{5/3}$. Global buoyancy conservation gives $g_1^{\prime }Q=B_1$. Eliminating $Q,g_1^{\prime }$ yieldsThis uniquely determines an interface height between zero and $H$ for a given positive area parameter. Opening area controls the nondimensional height; source strength controls the density contrast and flow magnitude. Factors of two depend on the stated effective-area convention, which has been fixed explicitly here.

Let $S(z)=d{\rho }_s/dz$, $N^2=-gS/{\rho }_0>0$, and define the quasi-geostrophic streamfunction by $\psi =\tilde{p}/({\rho }_0{f}_0)$. Use horizontal scale $L$, vertical scale $D$, velocity scale $U$ and advective time $L/U$. The required quasi-geostrophic approximation has small Rossby number $\epsilon =U/(|f_0|L)$, $\beta L/|f_0|=O(\epsilon )$, and Burger number $N^2{D}^2/(f_0^2{L}^2)=O(1)$. Hydrostatic, thin-layer scaling has $D/L\ll 1$; density variations are small compared with ${\rho }_0$, as required by the Boussinesq approximation. The leading horizontal flow is nondivergent, with vertical velocity $w=O(\epsilon UD/L)$. Retain forcing of the same slow order as density advection, $r=O(U\tilde{\rho }/L)$; stronger forcing need not preserve this balanced scaling.

Leading geostrophic balance and the hydrostatic approximation giveWrite $D_g/Dt={\partial }_t-{\psi }_y{\partial }_x+{\psi }_x{\partial }_y$. The density equation becomesTaking the curl of the horizontal momentum equations at the next order and using continuity gives the stretching equationSubstitute the expression for $w$. Vertical differentiation commutes with the geostrophic material derivative acting on ${\psi }_z$ here: the extra Jacobian is $J({\psi }_z,{\psi }_z)=0$. Since $N$ depends only on height, this yields the diabatically forced quasi-geostrophic potential vorticity equationThe minus sign is appropriate to the specified density forcing; a positive upward buoyancy forcing is $R=-gr/{\rho }_0$ and gives $D_{g}q/Dt=f_0{\partial }_z(R/N^2)$.

Now take constant $N$ and linearize about rest. Put $s=f_0^2/N^2$ and $S_0=g{f}_0\mu r_0/({\rho }_0{N}^2)$. The forced equation isFor the harmonic heating response of quasi-geostrophic flow, represent the two traveling Fourier components with frequencies ${\omega }_{+}={\omega }_0$ and ${\omega }_{-}=-{\omega }_0$:Each component of the forcing has amplitude $S_0/2$, soEquivalently its vertical structure obeysThese are the time-periodic forced responses; unspecified initial data may add free Rossby waves. Selecting no incoming waves describes forcing switched on causally and then its persistent harmonic part.