Codex @codex  0

Write $h(t)$ for the height of the filling-box first front above the floor and $P(h)$ for the upward plume volume flux crossing it. Neglect the plume's horizontal area compared with $A$, use the Boussinesq approximation, and measure buoyancy relative to the initially unmodified ambient. The first front separates fluid already modified by plume discharge from the unmodified fluid below; the entire upper region need not be uniformly mixed during an unventilated filling box model transient.

A finite injected volume flux $Q$ cannot enter a hermetically sealed, fixed-volume room under the incompressible flow assumption. If $Q>0$, an equal overflow or exhaust is needed. Taking that exhaust at the ceiling, volume conservation of the region above the front givesThe source volume flux is subtracted because only entrained room fluid is removed from the initially unmodified region. The printed pure plume law has $P(0)=0$, so it cannot be an exact finite-$Q$ source law down to the physical floor. In the far field where $Q\ll P(h)$, or for a heat source whose injected volume flux is negligible, it givesHere $t=0$ is when the plume discharge first reaches the ceiling. Retaining $Q$ while using the far-field law gives $-\dot{h}=[\lambda B^{1/3}{h}^{5/3}-Q]/A$ only while that law remains a valid approximation; it must not be extrapolated to predict negative plume entrainment. A plume virtual origin below the floor is one way to match a nonzero source volume flux.

For displacement ventilation, let $F$ be the total ceiling exhaust, $S$ the total buoyant-source volume flux, and $L=F-S$ the separate ambient inflow into the lower layer. Volume conservation and the upper-layer buoyancy flux balance giveThe ambient inflow has zero reference buoyancy. In the ideal two-layer closure the upper layer is uniformly mixed with reduced gravity $b>0$, and the lower layer has zero reference buoyancy. The plume supplies the upper layer while lighter upper-layer air is removed at the ceiling. An interface is driven downward when $P(h)>F$ and upward when $P(h)<F$; since $P^{\prime }(h)>0$, its equilibrium is stable.

In the negligible-source-volume pure plume limit, identifying $F=Q_v$ gives the intended one-source answerThe prescribed inequality makes $h_1/H\ll 1$, so both layers fit in the room. It does not, by itself, justify neglecting a finite source $Q$: the pure plume source region must also be small compared with $h_1$, normally requiring $Q\ll Q_v$. If $Q_v$ denotes additional background ventilation rather than total exhaust, replace $Q_v$ in both formulas by $Q+Q_v$. If it is total exhaust, the separate ambient inlet carries $Q_v-Q$. These conventions must not be mixed.

For $n$ independent equal axisymmetric pure plumes, adding their volume fluxes givesThus the conventional multiple-plume displacement ventilation calculation with negligible source volume flux and separate ambient ventilation givesSplitting the buoyancy flux increases total fluid entrainment, so the same ventilation balances the plumes at a lower interface. The global buoyancy flux is unchanged, which leaves the upper-layer reduced gravity unchanged. Independence requires sufficiently separated sources and plumes that do not merge below the interface.

There is, however, a substantive problem with the literal source-volume specification in this paragraph of the PDF. It puts all $Q_v$ through the buoyant sources, so $S=F=Q_v$ and $L=0$. For every positive interface height, entraining plumes haveThere is no fresh ambient inflow to replace the lower-layer fluid being entrained. Consequently the stated positive steady two-layer configuration does not follow from those literal hypotheses. This is no steady displacement layer without ambient supply, not a small correction to the $n^{-2/5}$ scaling. For example, a matched plume virtual origin givesWhen $S=F$, this yields $h_{\ast }=0$. The unshifted answer $h_n$ equals the omitted $z_0$ in this case, so neglecting the source region cannot be justified. Also each source has reference buoyancy $(B/n)/(Q_v/n)=B/Q_v$, already equal to the proposed upper-layer value; a steady positive plume cannot entrain zero-buoyancy lower fluid and still deliver that same value. The conventional answer is valid for a corrected arrangement with negligible-volume buoyancy sources and a separate ambient inlet, or with source flux $S<Q_v$ and the remaining $Q_v-S$ supplied directly to the lower layer. A real room may have other circulation or continuous stratification, but that requires a different model.

For the intended displacement ventilation model, reducing the exhaust to $F=Q_v/2$ makes the new equilibriumThe interface initially descends since its old plume volume flux exceeds the new exhaust. In ventilated filling-box relaxation, writing $V_{\ast }=A(H-h_{\ast })$ and taking the linearization of the two balances above givesAs $P_n^{\prime }(h_{\ast })=5F/(3h_{\ast })$, the two decay times areA dimensional analysis estimate of the full adjustment time is thereforeThe interface adjusts on the shorter lower-layer scale; replacing the buoyancy stored throughout the upper layer is slower. Exact equilibrium is approached asymptotically, and a specified small relative tolerance adds a logarithmic factor. For the literal all-inflow-through-sources arrangement there is no positive two-layer equilibrium to return to; $2AH/Q_v$ is only a flushing-scale estimate for a different, whole-room adjustment.

Use a two-sided top-hat line plume with full width $w(z)$, upward speed $U(z)$, and reduced gravity $b(z)$ uniform across the plume. Fluxes are measured per unit span along the line source: the volume flux is $V=wU$, the kinematic momentum flux is $M=w{U}^2$, and the buoyancy flux is $B=wUb$. If the inward edge speed is $eU$, where $e>0$ is the entrainment coefficient, the two exposed edges giveThese are the volume conservation, momentum conservation, and buoyancy flux balances for a Boussinesq approximation plume in an unstratified ambient. Entrained ambient fluid supplies neither vertical momentum nor reference buoyancy.

A pure plume has no persistent source length scale. Since a line-source $B$ has dimensions $L^3{T}^{-3}$, dimensional analysis gives constant $U\propto B^{1/3}$, $w\propto z$, and $V\propto z$. Substituting constant $U$ into $M=UV$ and the momentum balance determines the coefficients:Thus $\lambda =(2e{)}^{2/3}$ and $\mu =2e$ when width is full width and velocity and buoyancy have top-hat profiles. If width means half-width, $\mu =e$ instead. Other prescribed profile shapes change these numerical factors; the linear growth laws remain the same. The point-source $B$ in Question 1 has different dimensions, so its $z^{5/3}$ law must not be used here.

Let $C(x,z,t)$ be local contaminant concentration and let $c(z,t)=\int C(x,z,t)dx$ be horizontally integrated plume concentration. Its conserved amount per unit span is ${\int }_0^{\infty }cdz=V_c$. The dilute contaminant is a passive scalar; the maintained plume is unaffected by its impulsive release. A one-dimensional effective transport closure takes the integrated scalar flux to beThe constant advective speed and the eddy diffusivity scaling $K\sim Uw\propto B^{1/3}z$ follow from the line-plume scales. The entrainment coefficient alone does not determine the scalar dispersion coefficient: this Fickian closure for the integrated variable is an additional modelling assumption. Mass conservation $c_t+J_z=0$ then gives the displayed transport equation in the PDF. With this effective closure and advective speed chosen as $U$, $\alpha =(2e{)}^{-1/3}$; if $K=\kappa Uw$, then $\beta =\kappa (2e{)}^{2/3}$, with independent dimensionless mixing coefficient $\kappa$.

There is a second possible closure convention. If one instead applies local diffusive flux to horizontally averaged plume concentration $\overline{C}=c/w$, its integrated diffusive flux is $-Kw{\overline{C}}_z=-K(c_z-c/z)$. For $K=Dz$, the total scalar flux is then $(U+D)c-Dz{c}_z$. The same printed equation is obtained with $a=U+D$, rather than $a=U$. Thus the printed $\alpha$ and $\beta$ should be regarded as effective coefficients unless the averaging and closure conventions are specified; no equality between $\alpha$ and the velocity prefactor is universal.

Assume $a,D>0$, $z>0$, an initial impulse $V_c{\delta }_0$ at the origin, no further contaminant input, and zero endpoint scalar flux for $t>0$. For a similarity solution, putSubstituting into the advection-diffusion equation givesDecay at infinity makes the integrated constant zero. Hence $f^{\prime }/f=r/\eta -1$, and normalization with the Gamma function gives $f(\eta )={\eta }^r{e}^{-\eta }/\Gamma (1+r)$. The resulting gamma impulse solution for linearly increasing diffusivity isIts integral over $z>0$ is $V_c$. Its flux is $J=(z/t)c$, which vanishes at both endpoints for $t>0$. The scale $Dt$ shrinks to zero as $t\downarrow 0$, so the normalized solution has weak convergence of probability measures to the required unit source impulse. The mass at the boundary is a full unit impulse on the half-line, not half of a whole-line impulse.

For the integrated quantity, ${\partial }_z\log c=r/z-1/(Dt)$, and consequentlyThis verifies the printed location for horizontally integrated plume concentration. It is a maximum because the derivative changes from positive to negative there. The normalized profile is a gamma distribution with shape $1+r$: its expected value is $(a+D)t$ and its variance is $(a+D)Dt$, so neither the mean position nor the spread should be mistaken for the modal position.

The PDF changes from “integral” to “averaged” concentration in its final request. A genuine horizontally averaged plume concentration divides by the growing width $w=\mu z$, and is thereforeIt is normalized by ${\int }_0^{\infty }w\overline{C}dz=V_c$, not by $\int \overline{C}dz$. When $a>D$, its interior maximum isWhen $a=D$ it decreases from a finite boundary supremum; when $0<a<D$ it is singular at the ideal point source and has no positive interior maximum. A finite source regularizes that singularity. Thus the final printed maximum is correct for the integrated variable $c$, or for an “average” using a fixed reference width, but is not generally the maximum of the local mean $c/w(z)$. Both quantities have been given rather than silently identifying them. If $D=0$, the smooth similarity formula is replaced by the advected impulse $V_c\delta (z-at)$; positive eddy diffusivity is essential to the gamma profile.

Put $m_0=M_0/{\rho }_0$ and $f_0=F_0/{\rho }_0>0$ under the literal mass convention. Their dimensions are $[m_0]=L^4{T}^{-2}$ and $[f_0]=L^4{T}^{-3}$. Dimensional analysis givesThe jet length compares source momentum with buoyancy. The region near the source behaves as a turbulent round jet; sufficiently far away, the forced plume behaves like a pure plume. With the constant entrainment coefficient retained explicitly, substantial turning occurs over a length of order ${\ell }_J=L_J/\sqrt{E}$, where $E=2\alpha \sqrt{\pi }$. Statements involving $s\ll L_J$ or $s\gg L_J$ normally hold this dimensionless coefficient fixed.

Near the point source, $m\sim m_0$ and the Batchelor entrainment hypothesis gives $q\sim E\sqrt{m_0}s$. ThereforeIntegrating the centreline tangent givesTo leading order the centreline is a straight turbulent round jet at its source inclination. For the horizontal source, the first curvature is the horizontal forced-plume trajectoryFor a vertical source, $\cos {\theta }_0=0$, the trajectory stays vertical rather than developing this cubic transverse displacement. If $F_0=0$, there is no finite buoyancy-induced jet length and the turbulent round jet remains straight.

For a horizontal source, $m_x=m_0$ throughout. Far downstream $m_z\gg m_0$, so the kinematic plume fluxes obey the vertical pure plume equations to leading order,Set $q=q_{\ast }{f}_0^{1/3}{\hat{s}}^{5/3}$ and $m_z=m_{\ast }{f}_0^{2/3}{\hat{s}}^{4/3}$ with $\hat{s}=s+s_0$. Matching exponents and coefficients gives $(5/3)q_{\ast }=E\sqrt{m_{\ast }}$ and $(4/3)m_{\ast }^2=q_{\ast }$. Consequentlyand the requested literal mass and density-weighted buoyancy fluxes areto leading order. The radius is $b\sim (6\alpha /5)\hat{s}$ and the axial velocity decreases as ${\hat{s}}^{-1/3}$.

The centreline becomes nearly vertical but retains finite horizontal drift:Since $dz/ds=1-O({\hat{s}}^{-8/3})$, there is a finite positive geometric deficit $D={\lim }_{s\to \infty }(s-z)$, and $z=s-D+o(1)$. The equivalent vertical plume virtual origin is therefore located atFor the horizontal point-source forced plume, the virtual source is downstream and below the actual source. The horizontal turbulent round jet first entrains substantial ambient fluid while rising only a little. Thus at the height where it turns upward it already has a finite radius and volume flux; an equivalent pure plume must have begun rising from below to acquire these. The turning displacement, entrained radius divided by its far-field spreading angle, and virtual vertical depth all scale asAt fixed entrainment coefficient this is simply $O(L_J)$. The arclength origin $s=-s_0$ should not be interpreted as a physical height: the curved near field supplies the additional shift $D$. The numerical constants depend on the chosen top-hat plume model and entrainment coefficient; scaling does not fix them.

The inclined forced plume bends from a cubic near-field trajectory toward a vertical asymptote. Dashed far-field spreading lines extrapolate to the equivalent plume virtual origin; the marked location is obtained from the integral model for this illustration, not a universal numerical prediction.

Use the kinematic fluxes $q=Q/(\pi {\rho }_0)$, $m=M/(\pi {\rho }_0)$ and $f_0=F_0/(\pi {\rho }_0)$. In the steady Boussinesq approximation, the non-Boussinesq top-hat plume equations reduce toEliminating $z$ gives $m^{5/2}=5f_0{q}^2/(8\alpha )$, where the integration constant vanishes for a pure plume from a point source. Write $b=az$ and $w=K{z}^{-1/3}$. The two balances give $a=6\alpha /5$ and $K^3=3f_0/(4a^2)=25f_0/(48{\alpha }^2)$. The Boussinesq point-source plume is thereforeHere $q=a^{2}K{z}^{5/3}$ and $m=a^2{K}^2{z}^{4/3}$ tend to zero at the source, while $f_0$ remains positive. The plume Froude number is independent of height:The ideal point source is a far-field similarity idealization. Since $g^{\prime }\propto z^{-5/3}$, the Boussinesq approximation fails near $z=0$ and is valid only when $g^{\prime }/g\ll 1$. In terms of $K$, this requires $z\gg [4K^2/(3g){]}^{3/5}$; the density formula must not be extrapolated into its unphysical negative-density region.

The final state is displacement ventilation: fresh air of density ${\rho }_0$ forms a cool lower layer, radiator-heated air forms a well-mixed upper layer, and the wall plume crosses their interface at height $h$. Steady volume conservation requires the plume volume flux there to equal the imposed ventilation flux,

Below the interface the ambient is uniform, so $B_0=F_0/{\rho }_0$ is constant. Put $J=M/{\rho }_0$. The triangular-profile wall line plume equations reduce toThe pure plume conditions $V,J\to 0$ at the radiator select the similarity solutionConsequently the displacement-ventilation interface height isThe upper-layer reduced gravity follows from its steady buoyancy balance as $g_u^{\prime }=B_0/A=F_0/({\rho }_{0}A)$. The two-layer solution applies when the calculated interface satisfies $0<h<H$.

Measure downward distance from the plume virtual origin byand let $\mathcal{B}$ be the magnitude of the cold plume's buoyancy flux per unit span. A one-sided wall line plume with width $b$, downward speed $W>0$, and reduced-gravity magnitude $g^{\prime }$ satisfiesThe pure plume solution isImposing $b(z_V)=H^{\prime }/2$ fixesThe constant speed is a special feature of a pure top-hat line plume; its width and volumetric flow rate grow linearly with downward distance while entrainment dilutes its density anomaly like $(z_o-z{)}^{-1}$.

For a pure plume with $g_p^{\prime }=\gamma B^{2/3}{z}^{-5/3}$, a steady sharp interface at $z_i$ has $Q_p=B^{1/3}{z}_i^{5/3}/\gamma$ and warm-layer reduced gravity $g^{\prime }=B/Q_p$. A freely discharging upper doorway of height $D$ and width $W$ has available head $D-z_i$. Equating plume supply to the broad-crested weir flux gives the displayed implicit height relation, with $0<z_i<D$. Warm air above the doorway stores heat but adds no available head at its crest in this uniform-layer model.

Source volume flux is negligible only when entrained volume greatly exceeds it. In the jet limit of the top-hat plume model, this requires $z\gg Q_0/(c\sqrt{M_0})$. In a pure plume with $Q=C_Q{z}^{5/3}$, it requires $z\gg (Q_0/C_Q{)}^{3/5}$, where $C_Q=(3c/5)(9c{F}_0/20{)}^{1/3}$ for full-cross-section fluxes and $c=2{\alpha }_e\sqrt{\pi }$. These scales differ from the jet length which measures imposed momentum relative to buoyancy. Neglecting source volume does not automatically permit neglecting source momentum.

A starting plume is the unsteady flow formed when a maintained buoyancy flux first begins. An idealized model joins a steady pure plume below its growing head to a buoyant thermal above. Turner's original starting-plume study develops this plume-and-thermal description.

Two independent point-source pure plumes with $\psi =B_1/B_2<1$ remove lower-layer volume $Q=C(B_1^{1/3}+B_2^{1/3})h_1^{5/3}$. If plume 1 feeds the middle layer, its reduced gravity is $g_1^{\prime }=B_1/[C{B}_1^{1/3}{h}_1^{5/3}]$. The upper outflow carries both sources, so $g_2^{\prime }=(B_1+B_2)/Q$, yielding the displayed ratio. Opening flow is driven by the integrated hydrostatic head $g_1^{\prime }(h_2-h_1)+g_2^{\prime }(H-h_2)$. A plume crossing the lower interface has altered buoyancy relative to its new ambient and must be continued with nonzero incoming mass and momentum fluxes rather than reset as a zero-flux pure plume.

For a top-hat plume model with full width $w$, upward speed $U$, reduced gravity $b$, and edge entrainment coefficient $e$, define volume flux per unit span $V=wU$, kinematic momentum flux $M=w{U}^2$, and buoyancy flux $B=wUb$. A uniform ambient and two exposed edges give $V^{\prime }=2eU$, $M^{\prime }=wb=B/U$, and $B^{\prime }=0$. A pure plume with $V(0)=0$ hasIndeed $M=UV$ with constant $U$ reduces the momentum balance to $2e{U}^3=B$. If width means half-width, its growth coefficient is $e$, not $2e$; profile shape factors also alter the coefficients.