For the pressure-neglecting confined plume with compensating return flow, define specific plume momentum flux and area fraction . Then and
The branch issuing from a pure plume reaches a fold at , where plume and ambient have equal areas, opposite velocities and equal specific momentum fluxes. The numerator is nonzero on this branch, so diverges although remain finite. This signals failure of the steady uniform-profile description, not infinite physical velocity.
For an axisymmetric top-hat plume model in uniform ambient density, use full kinematic fluxes and . The equations , conserve . Solve this invariant for , then integrate to obtain the displayed exact quadrature. It determines the transition from the near-source jet limit of the top-hat plume model to the distant pure plume. The two regimes cannot be added by linear superposition because the flux equations are nonlinear.
Forced plume Created 2026-10-05 Updated 2026-10-06
A forced plume has both imposed momentum flux and positive buoyancy flux. Its near-source region behaves like a jet, while its far field approaches a pure plume. The crossover is measured by a jet length.
Merger of two equal pure plumes Created 2026-10-06 Updated 2026-10-14
If two equal pure plumes merge while keeping their common velocity and reduced gravity, conservation of cross-sectional area gives radius times the individual radius. All three integral fluxes double. The plume balance parameter increases by , so the immediate combined plume is a lazy plume. Similarity attraction may restore pure balance farther downstream, but instantaneous pure balance does not follow from adding two pure plumes.
A constant surviving buoyancy flux sets the amplitudes of a pure plume: and . Relative-flux perturbations have the neutral ratio . This is the zero-eigenvalue limit of growing mode of power-law plume similarity; nonzero endpoint stratification can still create cumulative logarithmic flux loss.
For an ideal-gas heated pure plume with entrainment proportional to , the excess heat flux implies conserved kinematic buoyancy flux . Put , . Their balances are and , identical to the pi-suppressed Boussinesq pure-plume flux equations. With and , , and . The radius is and reduced gravity is . At density is half ambient; far above this height the Boussinesq approximation becomes accurate.
Consider entraining plumes supplied with total source volume flux , a separate lower-layer ambient inflow , and ceiling exhaust . Write , where is total fluid entrainment from the lower layer. Volume conservation gives . A steady positive lower layer therefore requires . If and for every , no positive steady displacement-ventilation interface height exists. In a shifted pure plume approximation matched to source flux at , this follows immediately from for . This is a statement about the ideal two-layer model, not a prohibition of all possible stratification in a real ventilated room.
For a rising light plume use 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 are
Here 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 , where is the dimensionless entrainment coefficient. Volume added per unit height is . With , the integrated plume balances in uniform ambient density are
Seek a pure plume with no persistent source volume or momentum scale. Equating powers of gives and ; constant buoyancy gives . Equating coefficients yields, relative to the ideal source or a fitted plume virtual origin,
Thus , and . A finite nozzle or nonzero source momentum requires a near-source transition rather than this singular point-source solution at arbitrarily small .
In a stably stratified ambient, entrained fluid is progressively lighter with increasing height. Relative to that ambient, the plume's buoyancy flux satisfies . 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 , the balances give a useful exact relation:
for an ideal zero-momentum source. The rising branch begins at , attains its largest momentum at , and reaches its terminal zero-momentum limit at . Entrainment and decreasing buoyancy make the associated rise interval finite.
An estimate follows already from the unstratified pure plume flux: , so . Hence
The 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 and . The sketch distinguishes the neutral level from the maximum rise.
Figure 1.
Entrainment, loss of buoyancy, momentum overshoot and lateral intrusion of a stratified plume
. Details.
The ratio compares the reference constant-stratification rise length with the height at which the actual buoyancy frequency equals . Equivalently, the actual gradient sampled around that reference height is . Positive means stratification strengthens upwards; negative means it weakens upwards.
For a finite-rise branch at fixed , increasing makes the gradient stronger at each fixed normalized height , so decreases. A small means a larger region of weaker-than-reference stratification and permits relatively greater rise. For , the same comparison reverses: increasing weakens the sampled stratification, increases , and can permit escape when the upper stratification decays sufficiently fast. At the ratio is one. Comparison at is not an exact equality for nonzero , because the plume samples a whole range of heights, not only .
Using the pure plume as a scaling estimate before appreciable buoyancy loss,
when and the lower integral is regular. This gives the preceding monotonic trends, but not a universal numerical height. In particular, sufficiently negative powers are singular at a literal point source at . A finite source elevation or a regularized near-source stratification must then be specified, as in part (d); no finite positive source buoyancy can absorb an infinite near-source buoyancy loss. Once such a finite source is used, has an integrable far-field tail for . A sufficiently buoyant plume can then retain positive buoyancy and rise without bound. The exact decaying-buoyancy infinite-rise branch in part (d) is another, specially matched possibility.
Write for the height of the warm-layer interface above the heat source. A negligible-volume point heat source supplies a pure plume with conserved buoyancy flux . The given plume buoyancy implies its volume flux at the interface is
At steady displacement ventilation, the plume feeds the warm layer at the same rate it empties through the upper vent: . The upper-layer buoyancy balance gives , so
Eliminate and to get
The height equation has one solution with . Its height is independent of source strength in this ideal pure-plume model, while its warm-layer density contrast increases as . The sketch shows cool replacement air below, an entraining plume through that region, and the upper mixed layer feeding the ceiling vent.
Figure 1.
Plume-fed displacement ventilation through two vents and through a hydraulically controlled upper doorway
. Details.
Normalize all fluxes by , rather than including that factor in only some of them. The specific plume volume flux, momentum flux and buoyancy flux are
The ambient transport is and its specific momentum flux is
Put , with . The total specific momentum flux is
The first differential equation follows immediately from relative-speed fluid entrainment:
For the second, differentiate and use . Since , this yields
This is the requested two-equation system, with and zero total volume flux as the additional conservation conditions. Near the pure source, and . Eliminating height in that limit gives
The 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 , which recovers the unconfined pure-plume limit near the source.
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 . Let be the upper-layer reduced gravity. Steady volume conservation makes the magnitudes of inlet and outlet flow equal to .
If the opening pressure drops are and , the orifice laws and hydrostatic pressure difference give
Define the reduced area by
Then the throughflow in this displacement ventilation model is
Lower-layer conservation requires . Global buoyancy conservation gives . Eliminating yields
This uniquely determines an interface height between zero and 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.
Above , plume 2 entrains middle-layer fluid rather than the original ambient. That fluid is already warm and light, so the plume's buoyancy relative to its surroundings is reduced. Its pre-interface density contrast and zero-source pure plume formula cannot simply be continued unchanged through the interface.
A consistent model treats the middle layer as a uniform new ambient, assumes a sharp interface and weak overall density differences, and carries the incident volume and momentum fluxes continuously across it. The relative buoyancy flux immediately above the interface is
Continue an entraining plume model from these nonzero incident fluxes, using buoyancy measured relative to the middle layer. In that uniform layer the relative buoyancy flux is conserved, while entrainment raises the absolute buoyancy carried relative to the outside environment. This is generally a finite-flux or forced plume continuation, not a newly reset point-source similarity solution. A virtual-origin approximation can be used only after matching its fluxes appropriately.
At its matching level , the plume density reaches the upper-layer density; it spreads and detrains, with possible finite-momentum overshoot. The model neglects mutual plume entrainment, wall heat transfer, interface thickness and detailed intrusion dynamics, and assumes each ambient layer is well mixed. In the equal-source limit , the continuation is a neutrally buoyant momentum-driven flow, emphasizing why the upper interface loses its independent density meaning.
The Batchelor entrainment hypothesis assumes that the mean ambient inflow normal to the edge of a turbulent turbulent plume is proportional to its representative axial velocity, . With top-hat plume model fields and circumference , the added volume per height is . It is an integral turbulent mixing closure; it does not assert that every instantaneous eddy has the same inward speed. Use constant positive , uniform resting ambient, Boussinesq approximation, negligible drag and no buoyancy loss.
Since , . Eliminate between the volume and momentum balances:
The zero-source-flux condition gives , the pure plume balance. Substitution in and integration gives
Recovering the top-hat fields from , and yields
The paper's ratio is therefore
This is the constant inverse-normalized ratio used here, not the catalog's normalized plume balance parameter, which is one for a pure plume. The zero-flux point source is a limiting similarity origin: the speed and density anomaly diverge there, so a physically finite source must regularize the near field. The expressions are valid in the region where the Boussinesq approximation and turbulent closure apply.
Differentiate and use both balances:
Substituting gives
For , the radius grows faster than the pure-plume slope and its growth rate decreases toward that slope as increases. For , it grows at that slope exactly. For , it grows more slowly, possibly initially contracting, and its slope increases toward . Since and , all branches satisfy as .
Initial necking requires
Thus a sufficiently lazy plume initially narrows as buoyancy accelerates its fluid faster than fluid entrainment supplies additional volume. It still has increasing volume and momentum fluxes; shrinking radius does not imply volume loss. At equality the initial slope is zero, and for the slope later becomes positive once . This gives a finite minimum radius before the eventual pure plume regime.
Use and constant entrainment parameter . The printed integrals are normalized by , so a top-hat plume model has , and ; the physical area is still . Keeping this normalization avoids spurious factors of in the flux equations.
Eliminate height between the two flux balances. The plume flux-balance invariant follows from
so
For positive fluxes, both and increase. In particular , so as . Therefore , whatever the finite source imbalance, and . Integrating gives the attracting similarity solution
This proves attraction to pure plume similarity for physical positive source fluxes. Sources with a vanishing initial volume or momentum flux enter the positive regime immediately and can be obtained by the corresponding limiting solutions; a point buoyancy source has the same similarity form with zero virtual-origin offset. The statement assumes positive buoyancy flux: a jet with is a different asymptotic problem.
Define the plume balance parameter
Pure plume balance means at each height, or equivalently : momentum, buoyancy and entrainment have their similarity relation, even when the source has finite radius. A forced plume has , meaning excess momentum relative to this balance; a lazy plume has , meaning insufficient momentum. The invariant ensures that the sign of the imbalance is preserved while . In particular, a strictly positive source remains pure at every height precisely when
For this source the exact solution is the same similarity form with replaced by . Writing ,
Thus the radius growth of an axisymmetric pure plume is and the plume virtual origin lies at , with . Equivalently, in terms of top-hat source properties, pure plume balance is .
For the stated source-to-source geometry, the relevant horizontal distance is . Consequently the edge criterion is , rather than the radius at which two expanding edges first touch. Since ,
Here . A positive merging height requires ; equality puts it at the source, while an already wider source satisfies the geometric condition from the start. The no-interaction premise extrapolates independent plumes even though their edges would first touch earlier, when each radius is .
Use with and to get . The individual plume properties at the specified height are therefore
The source volume flux affects the height through the virtual origin, but it does not affect these properties at a prescribed radius in an exactly pure plume.
For the stipulated merger of two equal pure plumes, the new radius is . This factor is outside , as required by the area. Since speed and reduced gravity are held fixed, the merged fluxes are , and . Their balance parameter is
The immediate combined plume is lazy. The same result follows from : enlarging the radius at fixed increases this parameter.
At fixed , the radius required for pure plume balance is
That is half the stipulated combined area . This fixed-speed pure-plume merger and flux conservation comparison exposes the limitation of demanding immediate pure balance with unchanged speed and reduced gravity: reducing the area to this value would also halve all three merged fluxes. A conservative merger must instead adjust its properties and begin out of pure balance. Under the given constant-entrainment equations, the resulting lazy plume subsequently approaches the attracting similarity solution with total buoyancy flux .
Expand the square root and quotient in the logarithmic-height plume equations:
The linear stability analysis gives
In particular, the second row contains , since the reciprocal of contributes another .
For the growing mode of power-law plume similarity, the characteristic polynomial is
Within the physical interval and , while as . There is a positive real eigenvalue, so the similarity equilibrium is unstable with increasing height. Perturbations behave as ; the result is not stability in physical time. Thus this exact positive-buoyancy similarity trajectory is not a generic attracting far field.
At the upper endpoint,
The limiting polynomial is , reproducing the three stated eigenvalues. The positive eigenvalue tends to zero, with . The two negative modes decay as and ; the remaining mode is neutral at the limiting linear system. Its eigenvector is in relative-flux coordinates.
The neutral buoyancy-flux mode of a pure plume has a simple interpretation: a constant surviving buoyancy flux fixes the amplitudes of a pure plume, with and . A small fractional change in produces changes in in the ratio . It is retained rather than damped out.
There is a qualification at the exact endpoint. For fixed nonzero and , taking gives , not zero. The flux coefficients above also diverge as . Thus the constant-flux powers are the formal limiting scaling, not an exact constant-buoyancy solution in a nonzero endpoint stratification; cumulative ambient effects can require logarithmic corrections.
Write for the height of the filling-box first front above the floor and for the upward plume volume flux crossing it. Neglect the plume's horizontal area compared with , 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 cannot enter a hermetically sealed, fixed-volume room under the incompressible flow assumption. If , an equal overflow or exhaust is needed. Taking that exhaust at the ceiling, volume conservation of the region above the front gives
The source volume flux is subtracted because only entrained room fluid is removed from the initially unmodified region. The printed pure plume law has , so it cannot be an exact finite- source law down to the physical floor. In the far field where , or for a heat source whose injected volume flux is negligible, it gives
Here is when the plume discharge first reaches the ceiling. Retaining while using the far-field law gives 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 be the total ceiling exhaust, the total buoyant-source volume flux, and the separate ambient inflow into the lower layer. Volume conservation and the upper-layer buoyancy flux balance give
The ambient inflow has zero reference buoyancy. In the ideal two-layer closure the upper layer is uniformly mixed with reduced gravity , 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 and upward when ; since , its equilibrium is stable.
In the negligible-source-volume pure plume limit, identifying gives the intended one-source answer
The prescribed inequality makes , so both layers fit in the room. It does not, by itself, justify neglecting a finite source : the pure plume source region must also be small compared with , normally requiring . If denotes additional background ventilation rather than total exhaust, replace in both formulas by . If it is total exhaust, the separate ambient inlet carries . These conventions must not be mixed.
For independent equal axisymmetric pure plumes, adding their volume fluxes gives
Thus the conventional multiple-plume displacement ventilation calculation with negligible source volume flux and separate ambient ventilation gives
Splitting 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 through the buoyant sources, so and . For every positive interface height, entraining plumes have
There 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 scaling. For example, a matched plume virtual origin gives
When , this yields . The unshifted answer equals the omitted in this case, so neglecting the source region cannot be justified. Also each source has reference buoyancy , 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 and the remaining 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 makes the new equilibrium
The interface initially descends since its old plume volume flux exceeds the new exhaust. In ventilated filling-box relaxation, writing and taking the linearization of the two balances above gives
As , the two decay times are
A dimensional analysis estimate of the full adjustment time is therefore
The 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; is only a flushing-scale estimate for a different, whole-room adjustment.