Past exam of the mathematics course of the University of Cambridge / 2002 / iii / Paper 52 / Section C

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.

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.

Let ${\mathcal{B}}_0>0$ be the source buoyancy flux per unit length. The conversion from the heat-source strength depends on its units: if $F$ is a kinematic temperature-volume flux, ${\mathcal{B}}_0=g\beta F$; if it is a physical heat power per unit length, ${\mathcal{B}}_0=g\beta F/({\rho }_0{c}_p)$, where $c_p$ is specific heat capacity at constant pressure. This explicit convention avoids treating heat power as a buoyancy flux without conversion.

For uniform ambient $N^2=0$, so $\mathcal{B}={\mathcal{B}}_0$. The ideal pure line-source limit has negligible initial volume and momentum transport. The flux equations admit $w=w_{\ast }$ constant, $q=2\alpha w_{\ast }z$ and $b=\alpha z$. Substituting into the momentum equation gives $2\alpha w_{\ast }^2=2bG$, and the conserved buoyancy flux then fixes $2\alpha w_{\ast }^3={\mathcal{B}}_0$. ThusThis verifies all three differential equations, not just their dimensional exponents. The plume grows in width while maintaining constant upward velocity; its temperature anomaly decreases like $z^{-1}$. The zero-height singularity is the ideal point-source limit. A finite source or nonzero momentum input introduces a plume virtual origin or a forced near-source region; the heat strength alone cannot uniquely determine that finite-source region.

Let the closed box have ceiling height $H$ above the source. The printed width $2R$ alone does not specify this height, so it must appear as an additional geometrical parameter. Take $t=0$ here when the plume has first reached the ceiling and its thin outflow has spread across the box. The leading filling box model assumes a negligible source volume, insulated boundaries, fast horizontal spreading, slow evolution compared with the plume transit time, negligible diffusion, and $b\ll R$. Below the interface the ambient still has its original temperature, so the local plume flux remains $q(z)=2\alpha w_{\ast }z$ from part (c).

The ceiling outflow creates a warm upper region. Fluid entrainment draws the original lower fluid into the plume, making that upper region deepen. The filling-box first front is the lower material boundary of fluid already modified by plume discharge, separating it from the initially unmodified ambient; it need not bound a uniformly mixed upper layer. With negligible occupied plume area, volume conservation at its height $h(t)$ givesThereforeIf source switch-on is the time origin instead, replace $t$ by elapsed time after the initial transit and spreading stage; the quasi-steady model does not resolve that stage.

The buoyancy jump requires distinguishing ambient fluid at the front from the plume crossing it. Just after the first ceiling outflow, the ambient buoyancy immediately above the nascent front isOutside the narrow plume, ambient temperature is transported without diffusion. If $G_a$ is its buoyancy relative to the original fluid, its equation isThe first front follows precisely these material trajectories, so the temperature of the fluid immediately above it is conserved. The initially unmodified fluid below has $G_a=0$. Thus the ambient jump isBy contrast, the contemporaneous plume immediately below the front has $G_p(h,t)=w_{\ast }^2/h(t)$, which grows as the front descends. That plume value is not the ambient jump across the material first front. Equating the two would inadvertently impose a uniformly mixed upper layer and discard the stable density stratification of the standard filling box model. Diffusion or extensive interfacial mixing can smooth and alter the ideal jump.

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.

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.

Normalize all fluxes by $\pi$, rather than including that factor in only some of them. The specific plume volume flux, momentum flux and buoyancy flux areThe ambient transport is $-Q$ and its specific momentum flux isPut $s=b^2/R^2=Q^2/(R^{2}M)$, with $0<s<1$. The total specific momentum flux isThe first differential equation follows immediately from relative-speed fluid entrainment:For the second, differentiate $J=M/(1-s)$ and use $s^{\prime }=2s{Q}^{\prime }/Q-s{M}^{\prime }/M$. Since $J^{\prime }=b^2{g}^{\prime }=B_{s}Q/M$, this yieldsThis is the requested two-equation system, with $B=B_s$ and zero total volume flux as the additional conservation conditions. Near the pure source, $Q,M\to 0$ and $s\to 0$. Eliminating height in that limit givesThe zero integration constant selects a pure plume; a nonzero imposed source momentum would select a different branch and generally a different turning height. The leading radius is $b\sim 6\alpha z/5$, which recovers the unconfined pure-plume limit near the source.

It is helpful to solve using the nonsingular total flux $J$ and volume flux $Q$ before recovering $M$. The relation $J=M/(1-Q^2/(R^{2}M))$ is quadratic:The plus sign selects the small-area pure-source branch. Its discriminant vanishes when $J=4Q^2/R^2$, givingThereforeThe specific upward plume and ambient momentum-flux contributions are equal, not opposite. The ambient's transported momentum has the opposite sign from its transported volume, so their product remains positive.

To show that the pure branch actually reaches this fold in finite height and that it is a genuine pole, introduceThen $s=q^2/m$ andThe pure-source condition is $j\sim (5/8{)}^{2/5}{q}^{4/5}$. Before the fold, $s<1/2$ and $m>2q^2$. Thus for every fixed positive $q_1$, the derivative $dj/dq$ is bounded above by a constant times $q^{-2}$ on $q\ge q_1$. If the branch continued for arbitrarily large $q$, $j$ would remain bounded, contradicting $j\ge 4q^2$. Hence it reaches $j=4q^2$ at a finite $q_o>0$. Its height is finite because $d\zeta /dq=(1-s)/(2\sqrt{m})$ is finite away from the source and behaves like a constant times $q^{-2/5}$ at the source, an integrable singularity.

For the numerator, write $m=j(1-s)$. ThenAt the fold $j=4q_o^2$, so $q_o^3\ge 5/256$. On the other hand the dimensionless numerator in the $m$ equation at $s=1/2$ isIt is strictly negative because $5/256>1/(64\sqrt{2})$. The denominator $1-2s$ approaches zero from above. Consequently $dM/dz\to -\infty$ at finite $z_o$, although the fluxes and velocities themselves remain finite. This proves the confined-plume momentum fold rather than merely noticing a possible zero denominator.

The dimensionless pure-source branch uniquely determines the remaining amplitudes. They can be specified without an elementary antiderivative by integrating the displayed scalar equation for $j(q)$ until $j=4q^2$, and settingA direct numerical integration gives $q_o\simeq 0.281610$, ${\zeta }_o\simeq 0.400049$, so the dimensional results for this specified entrainment convention areA convention using $\alpha w$ instead of $\alpha (w-v)$ would produce different numerical amplitudes; the relative-speed convention was stated in part (a).

Physically $z_o$ is a turnover or breakdown height for a slender steady plume surrounded by uniform counterflow. The opposing streams then occupy equal areas and have equal speeds; the single-valued small-area branch cannot continue smoothly. Above it one expects lateral spreading, stronger recirculation, mixing and unsteady buoyant structures rather than a continuation of the same slender top-hat column. A detailed flow or upper-layer evolution requires restoring pressure and stress dynamics and, eventually, the modified ambient stratification. The model does not predict an infinite physical velocity, nor does it require that all buoyant fluid forever remain below $z_o$.