1. Introduction
A round buoyant jet is the flow out of a pipe or circular nozzle that discharges into a quiescent or moving volume of fluid of different density. A vertical jet is positively buoyant if the direction of flow initially is that of the effective gravity, due to the density difference between jet and ambient fluid at the elevation of discharge; otherwise, it is called a negatively buoyant jet or fountain. A positively buoyant jet rises and mixes up to the surface of a surrounding motionless fluid of higher uniform density. A (momentum-dominated) fountain rises and mixes in an ambient of uniform, lower density to an elevation where the momentum flux of the mixed fluid vanishes. Then, the flow reverses and returns to the elevation of injection.
In a density-stratified, stable, calm ambient, a vertical jet that is positively buoyant at the elevation of injection entrains and mixes with the surrounding fluid up to an elevation where the average jet density is equal to that of the ambient fluid (neutral buoyancy). From that point on, the vertical kinematic momentum flux (inertia) of the jet will force it to a higher elevation, while ambient fluid still entrains and mixes with jet fluid, until the vertical momentum flux vanishes. Then, the jet flow reverses and returns to the elevation of neutral buoyancy, where it spreads horizontally. A vertical fountain (jet negatively buoyant at injection elevation) in a stable, density-stratified ambient will rise and mix with ambient up to the elevation where its kinematic momentum flux vanishes, under the influence of reverse buoyancy. Then, the direction of motion reverses and the fluid returns to an elevation of neutral density where it spreads horizontally. The maximum steady-state elevation attained by a vertical jet that is initially positively buoyant or a fountain in a density-stratified ambient is defined as “the terminal height of rise” (THR), and the mean elevation of neutral density where it spreads horizontally is defined as the “spreading height” (SH).
Atmosphere and oceans are both density-stratified. When studying the initial dilution of turbulent buoyant jets which discharge in the atmosphere or in a body of water, the ambient is assumed to be motionless, since the timescale of the initial evolution of the flow is rather short. Subsequently, it may be considered as a moving body of fluid for the study of the dispersion in the horizontal spreading field which is formed.
Following the pioneering work of Morton, Taylor, and Turner [
1], the mechanics of positively buoyant jets and fountains in a calm, linearly density-stratified ambient has been the subject of theoretical and experimental investigations for several decades. Vertical fountains in a uniform ambient fluid have also been studied in experiments [
2,
3]. Interest has been focused on the terminal rise height reached by the flow, as well as the mean spreading elevation where the mixed fluid subsequently expands horizontally. Direct applications of the results of such studies are related to the design of outfalls for the disposal of treated municipal wastewater and brine from desalination plants into the sea, as well as to the dispersion of smoke from factory stacks in the atmosphere. Further applications are the study of geophysical flows such as the motion of smoke plumes above fires, volcano eruptions, black smokers in the deep ocean, and magma chambers. Note that, after the eruptions of the Eyjafjallajökull volcano in Iceland (April–May 2010), ash and gases spread at an elevation around 7000 m in a very large area, causing great problems in the aviation industry for a long time.
Earlier studies of round buoyant jets and fountains in a linearly density-stratified ambient that focused on the terminal height of rise (THR) are summarized in [
4]. In cases where the initial momentum of the flow dominates the buoyancy force at the injection point, the flow is characterized as momentum-driven or as jet-like. If the initial buoyancy force dominates the momentum at the injection point, the flow is characterized as buoyancy-driven or as a plume. The fountain flow is jet-like since it is a consequence of the initial jet momentum. Fischer et al. [
4] using the experiments in [
5,
6] proposed that the normalized THR c
j (to be defined in
Section 2) for positively buoyant, momentum-driven flows (jets) is c
j = 3.80, as also presented in [
7]. In (positively buoyant) plumes, from the experiments in [
1,
8,
9], [
4] proposed the same constant c
p ≈ 3.80 (c
p is defined in
Section 2) for the normalized THR, while [
7] recommended a constant equal to 5.
Experiments on positively buoyant, vertical, round jets that discharge into a linearly density-stratified fluid were implemented by [
10]. The authors classified the flow into two asymptotic cases, jet-like flows and plumes, on the basis of the dimensionless parameter
MN/
B, where
M and
B are the initial specific (per unit mass) momentum and buoyancy fluxes and
N is the buoyancy frequency of the density stratification to be defined in
Section 2. Experiments in the regime 0.80 <
MN/
B < 2.40 regarding the THR were also reported [
11] from measurements. Similar experiments were reported for a wide range of
MN/
B by [
12,
13], the latter [
13] including data for the mean spreading elevation of the buoyant jet. Experiments were reported by [
14] for pure momentum-driven and negatively buoyant jets (fountains) in a linearly, density-stratified ambient. The authors [
14] also solved the steady-state equations of motion numerically and found the predictions to be congruent with their experiments.
Recently, the spreading height of buoyant jets in a linearly density-stratified ambient was measured [
15], concluding that the normalized spreading height is 1.53 when
MN/
B > 7 and 2.775 when
MN/
B < 7. Measurements of the velocity field of buoyant jets in linearly stratified ambient were performed [
16] using particle image velocimetry (PIV). The terminal rise height in plume-like flows was measured, along with the turbulence parameters that were evaluated for modeling.
The model proposed in [
1] was extended [
17] to compute the THR and SH in geophysical flows on Earth, as well as on other planets. An attempt was made to connect the rise height to a plume length scale (
Lp defined in next chapter), where both are normalized by the plume diameter, via the coefficient of entrainment that is variable in the range 0.05 to 0.16. This normalization does not comply with dimensional arguments presented in the next section.
The scope of the present investigation is to provide a systematic new set of data regarding the THR and SH of buoyant jets and fountains in a linearly density-stratified ambient and compare it to earlier experiments. Following the dimensional arguments (see next section), the experiments extend to a range of the dimensionless parameter
MN/
B that includes the two asymptotic regimes of initially jet-like flows and plumes, as well as the transition in between. Emphasis is given to the spreading height (SH) that is a very important parameter in the dispersion of buoyant jets in a stratified ambient, where available data are limited. An integral model using an improved entrainment coefficient function derived in
Appendix A from the entrainment equations is applied, in an effort to match the results obtained from the experiments. The entrainment function obtains a reduced coefficient if compared to that of a simple jet, when the buoyancy acts against the movement, in accordance with the results in [
18,
19].
2. Dimensional Analysis—Review of Earlier Experiments
A round vertical (positively) buoyant jet of density
ρo and uniform velocity
W discharges vertically into an infinite volume of calm, linearly density-stratified fluid. The initial jet volume flux
Q, specific (per unit mass) momentum
M, and buoyancy
B fluxes with dimensions L
3/T, L
4/T
2, and L
4/T
3, respectively, are calculated as follows:
where
ρa is the density of the ambient fluid at the jet nozzle elevation. The buoyancy frequency
N of the linearly density-stratified ambient with dimensions T
−1 is computed from
Let us assume fully developed turbulent flow and adopt the Boussinesq approximation (density differences do not exceed 5%). The terminal (steady-state) height of rise
Z is related to the initial flow parameters with a functional relationship
Z =
Z(
Q,
M,
B,
N). Far from the ejection point, the jet volumetric flow rate is much greater than the initial one, such that
Q may be neglected [
4]; thus,
Z =
Z(
M,
B,
N). Two length scales
Lj and
Lp may be defined, one using the initial jet momentum flux and
N and the other using the initial buoyancy flux and
N [
1,
4,
9,
20], as follows:
The jet is characterized initially as momentum (inertia)-driven or as buoyancy-driven, from the ratio of the inertia to buoyancy forces that is the initial densimetric Froude number [
4]
Fro =
W/(
g’
oD)
1/2,
g’
o = (
ρa −
ρo)/
ρo. Thus, the flow can be classified into four categories on the basis of the initial momentum and buoyancy fluxes
B and
M, respectively (or on
Fro), as well as the (initial) density difference between jet
ρο and ambient fluid
ρa at nozzle elevation: (i) jet-like flow when the momentum flux is dominant and
ρο ≤
ρa or
Fro is very large, (ii) plume-like flow, when the buoyancy flux is dominant and
ρο ≤
ρa or
Fro→1, (iii) buoyant jet (forced plume) flow when momentum and buoyancy fluxes affect the flow equally and
ρο ≤
ρa, and (iv) fountain flow when the momentum flux is dominant and
ρο ≥
ρa (
B < 0).
- (i)
A jet-like discharge (
B ≈ 0,
ρο ≤
ρa) will rise up to an elevation
Z =
Z(
M,
N) that is only a function of the initial specific momentum flux
M and the buoyancy frequency
N, thus leading to a dimensionless elevation,
which must be constant.
- (ii)
A plume (
M ≈ 0,
ρο ≤
ρa), will rise up to a dimensionless elevation,
which must also be constant.
- (iii)
A buoyant jet with initial kinematic momentum and buoyancy fluxes
M and
B, respectively, of similar strength that discharges into an ambient fluid with linear density stratification of buoyancy frequency
N will rise to an elevation
Z that is a function of
M,
B, and
N. Then, the normalized rise height will attain a functional form,
where
MN/
B is a characteristic dimensionless parameter. In asymptotically jet-like flows (
M >>
B),
MN/
B >> 1 and Equation (4) is valid, while, in plumes,
MN/B < 1 and Equation (5) is valid. If
MN/
B >> 1, dividing both sides of Equation (4) by
Lp obtains
meaning that the dimensionless jet rise
Z/
Lp is proportional to (
MN/
B)
1/4. Similarly, in plumes (
MN/
B < 1),
The analysis presented above is also valid if we substitute the terminal rise height Z with spreading elevation ZS, giving ZS/Lj = constant for jet-like flows and ZS/Lp = constant for plumes.
- (iv)
A fountain is the flow when the initial momentum flux is large but the buoyancy flux
B is negative (opposite to the direction of motion). In a uniform, calm ambient, the normalized terminal height of rise of a turbulent fountain is constant [
21],
if the length scale
lM =
M3/4/|
B|
1/2 [
2,
4] is large, or alternatively when the initial jet Richardson number
Rio =
Q|
B1/2|/
M5/4 <<1 [
21]. For fountains (large
lM) in a linearly density-stratified ambient, the asymptotic normalized terminal height of rise
Z is a function of
M and
N only and, therefore,
Z/
Lj takes a limiting value as shown from Equation (4). Thus, for large values of the parameter
MN/|
B|,
Z/
Lj = constant, regardless of the sign of
B. Rearranging the terms in Equation (4) gives
When the stratification parameter
N is small (
N→0,
MN/
B << 1) the ambient fluid asymptotically can be considered as homogeneous, and the turbulent fountain will reach a terminal rise height so that Equation (9) is fulfilled, i.e.,
Note that, in the case of a fountain, the spreading height (SH) elevation obeys the dimensional arguments made for the terminal rise height (THR).
An average dilution at elevation
z is the ratio
S =
μ(
z)/
Q of the local jet volume flux
μ(
z) to the initial one
Q. Following the same procedure as in the previous paragraphs, if the momentum is initially dominant (
MN/
B >> 1, jet-like flow), the normalized volume flux may be obtained as
If the buoyancy is dominant (
MN/
B < 1, plume-like flow),
where
dj and
dp are constants. The results from the asymptotic analysis are summarized in
Table 1.
Data from earlier experiments regarding the THR and SH were normalized and plotted together to justify the findings from dimensional analysis and evaluate the constants. The constants of proportionality from different experiments are shown in
Table 2. The parameters presented in this section are used for comparison with present data. The descriptions “ortho” and “round” for the data by Wong and Wright [
10] refer to the shape of the horizontal cross-section of the dispersion tank they employed, which was orthogonal and circular, respectively, for the two different sets of experiments they implemented.
The normalized THR
Z/
Lj and
Z/
Lp are plotted versus the dimensionless parameter
MN/
B in
Figure 1. From
Figure 1a, it is evident that, when
MN/
B > 10, the normalized THR
Z/
Lj is constant, in accordance with dimensional analysis, with an average value around
cj ≈ 3.50 according to all authors but [
14]. From
Figure 1b, it is clear that, when
MN/
B < 1, the normalized THR
Z/
Lp is a constant, congruent with dimensional analysis regarding initially plume-like flows, with an average constant from the data presented of
cp ≈ 4.35, greater than those measured in [
1,
11,
16]. Note that [
1,
14] added the distance
zo of virtual origin to the THR, estimated to be located 5 cm and 1 cm upstream from the pipe outlet, respectively. Data from other authors did not involve
zo, including Fan’s [
5] data, from which we subtracted 6.2 jet diameters.
3. Experiments—The Integral Model
Experiments were performed in a tank made of 1.25 cm thick Lucite with horizontal dimensions 0.90 m × 0.60 m and 0.80 m deep. An overflow around the top was used for excess fluid removal. The jet fluid was dyed with food dye, and the flow was illuminated with two 1000 W flood lights through a 5 mm thick panel of milky-white plexiglass, placed behind the rear window to diffuse the light uniformly. In the front and back (opposite) Lucite windows, square grids of 5 cm × 5 cm were drawn. For a vertical jet placed at the mid-plane of the tank, the terminal height of rise (THR) and mean spreading height (SH) were the average values read in the front and rear grid [
22]. The linear density gradient was obtained using the “two-tank technique” [
23]. Two similar tanks connected with a pipe were used, the first filled with fresh water and the other with an aqueous salt (NaCl) solution. Heavy fluid withdrawn from the continuously stirred second tank at small flow rates was pumped toward two floating surface spreaders positioned inside the dispersion tank. In this manner, lighter fluid was continuously added on the surface of the underlying heavier one, until the dispersion tank was filled up to the elevation of the overflow.
Four nozzles with diameters 0.5, 0.75, 1, and 1.5 cm were used to create the buoyant jets. They were mounted at the lower end of a jet plenum made of a 200 mm long piece of polyvinyl chloride (PVC) pipe with an inner diameter of 40 mm. A 2 cm long cylinder of sponge and a subsequent 6 cm long honeycomb section were set inside the jet plenum, to destroy any large eddies launched at the upstream end of it from the supply pipe and reduce the jet turbulence intensity. The inner side of all nozzles (inside the jet plenum) was rounded in order to provide smooth transition from the 40 mm pipe diameter to the size of the nozzle diameter, thus producing a uniform velocity profile at the exit [
24,
25].
When the stratification tank was filled up to the overflow, the surface spreaders were removed. The density profile was measured 20 min later waiting for the tank fluid to settle, using a calibrated microscale conductivity probe [
26] that was mounted on a vertical, 60 cm long, linear positioning device, the accuracy of which was 0.1 mm. Calibration of the conductivity probe was performed before each experiment using samples of salt water, the density of which was measured with a hydrometer of 2% accuracy. The conductivity probe was lowered abruptly, while the density and position were recorded simultaneously at a rate of 200 Hz with a 12 bit data acquisition card and Labview
® software. From the density profile (a typical one is shown in
Figure 2), the density gradient d
ρ/d
z was evaluated and the buoyancy frequency
N was computed. The jet plenum was placed at the top of the stratification tank pointing downward. The nozzle elevation was located at least 15 cm below the free surface, to ensure that the density gradient was linear throughout the whole depth where the jet flow developed. For positively buoyant jets at nozzle elevation, the jet fluid was heavier than the ambient, while, for negatively buoyant ones, the jet fluid was lighter than the ambient.
The jet water supply was through a constant head tank with an overflow, positioned 2 m above the free surface of the dispersion tank. It was filled with water of known density from a supply tank on the laboratory floor, via a submersible pump at flow rates higher than that of the jet, in order to ensure overflow. The jet flow rate was measured via a calibrated rotameter in the range 7–55 cm
3/s with 2.5% accuracy. All experiments were planned according to the dimensional analysis discussed earlier. Experiments were performed in the regime 0.10 <
MN/
B < 1000. The buoyancy frequency parameter
N2 varied from 0.058 to 0.16 s
−2, the Reynolds number varied between 1870 and 8710 (the transition to turbulence occurred within five jet diameters from the nozzle, [
24,
27]), and the initial jet Richardson number varied between 0.003 and 0.406.
The experiments were video-taped at 25 fps with a three-CCD (charge-coupled device) camera, positioned 4 m far from the dispersion tank. The terminal rise height and spreading elevation were evaluated from the time series of the movie frames obtained, using Autocad
® software and the two grids in the front and back transparent Lucite panels of the dispersion tank. The THR fluctuated around an average value, while the SH did not vary essentially throughout the experiment. Approximately 60 frames (an example is shown in
Figure 3) sampled every 1 s were used to evaluate an average of THR and SH, before the laterally expanding “wastefield” reached the tank walls. Initial and measured flow parameters in buoyant jets and fountains are shown in
Table 3.
Note that
N was computed as an “average over depth” least-square line from density profiles such as in
Figure 2. This means that deviations from the “linear density” due to possible flow rate changes of the pump supplying the surface spreaders would produce local density gradients that affect both THR and SH; hence, repeatability of the experiment was impossible. Therefore, our experiments were planned and executed by varying only the dimensionless parameter
MN/
B.
The equations of motion of the mean flow in round vertical turbulent buoyant jets that discharge in a density-stratified fluid with buoyancy frequency
N, assuming top-hat distributions for the time-averaged velocity and density excess or deficiency, are written [
4] as follows:
where
μ,
m, and
β are the local volume, specific momentum, and buoyancy fluxes, and
α is the local entrainment coefficient. The initial conditions at the jet origin
zo are
An estimate of the virtual origin distance
zo from the nozzle can be obtained from the theory by List and Imberger [
29]. In a buoyant jet, the ratio
Cp =
μ/
zm1/2 is defined to be the constant jet width parameter [
30] found to be equal to 0.27. Substituting the local volume and momentum fluxes in the equation above with the initial values, the jet virtual origin can be estimated as
zo =
Q/(
CpM1/2) = 3.28
D. The system of Equation (14) can be integrated starting at
zo with the initial conditions in Equation (15).
5. Discussion
From
Figure 4, it is evident that positively buoyant jets and fountains converge asymptotically to the same THR
Z/
Lj ≈ 3.50 and SH
ZS/
Lj ≈ 2 when
MN/
B > 10. The results regarding THR are congruent with earlier measurements [
5,
10,
12,
13,
16] but higher than the 3.0 reported in [
14], as shown in
Table 2. The results regarding SH are also higher than the 1.50 reported in [
14,
16] and lower than the 2.38 in [
13], as shown in
Table 2 as well. Analyzing the data in [
16] using dimensional arguments from
Section 2, the normalized SH elevation must be corrected to 1.89, a value that is near the measurements of the present study.
In the same figure, one may note that the data by [
14] regarding fountains in a linearly, density-stratified fluid for
MN/
B >> 1 give asymptotic dimensionless constants of 2.88 for THR and 1.52 for SH, which are also lower from those in the present investigation. There are also differences in the dimensionless THR and SH between the present measurements and those in [
14] when the parameter
MN/
B is small, the values of the latter being systematically lower than those measured in the present experiment. According to dimensional analysis presented in
Section 2, as
MN/
B→0 when
N→0 (uniform ambient density) and the jet momentum is driving the flow, the normalized THR
Z/
lM→
C [
21] or
Z/
Lj =
C(
MN/
B)
1/2. Recent experiments [
21,
31] indicated that
C ≈ 2, while [
32] reported a constant
C = 1.70. In
Figure 4a, the line
Z/
Lj = 2(
MN/
B)
1/2 is plotted for comparison, which is closer to the present data than those reported in [
14]. In the same graph, we show four points from the unpublished data of [
12] that are near the present findings.
When the flow is buoyancy-driven (
B > 0) using dimensional arguments in
Section 2, analysis of the present data in
Table 3 shows that the normalized THR
Z/
Lp ≈ 4.50 and SH
ZS/
Lp ≈ 3.0 when
MN/
B < 1. The results regarding the normalized THR are congruent to those in [
10,
12] but higher than those in [
1,
13,
16]. The results regarding the normalized SH are very close to 2.80 [
15] and 3.0 [
13], thus confirming their validity.
Numerical predictions using the entrainment coefficient derived from mass and momentum conservation equations led to accurate predictions regarding the spreading elevation and average dilution if compared to our measurements, as seen in
Figure 4 and
Figure 5, when the flow is buoyancy-driven (
MN/
B < 1). In this regime the predicted THR is a little lower than our measurements. In jet-like flow and fountain regimes (
MN/
B >> 1), however, the predictions are lower than our measurements (3.06 vs. 3.50) for THR and much lower (1.06 vs. 2) for SH. The normalized average dilution in this regime turned out to be lower than the measured minimum dilution [
10], a result that is against the physics of flow.
Nevertheless, in this regime, the results from integral models employed by others provided even lower values of normalized parameters. For example, the computed normalized THR in plumes [
10] was 3.55 vs. measured 4.45 and that in jets was 3.06 vs. 3.60. The computed normalized dilutions in [
10] were 0.42 vs. measured 0.68 in jets and 0.70 vs. measured 0.80 in plumes, results that are lower than minimum dilutions measured and contradictory to physics. Jirka [
33] computed the normalized THR and SH using CorJet software and found values of 3.30 and 1.50, respectively, values that are lower than those measured in jets. To estimate THR, in addition to the elevation
Z of the THR where momentum flux vanishes, √2
b (where
b is the local jet width) was added in [
34]; thus, (
Z + √2
b)/
Lj = 3.30.
6. Conclusions
Turbulent vertical buoyant jets and fountains in a linearly density-stratified fluid were investigated experimentally and numerically. Data from earlier investigations were collected, analyzed, and put together according to dimensional analysis, leading to asymptotic constants. A positively buoyant jet (B ≥ 0) is momentum-driven if the dimensionless parameter MN/B > 10. In this regime, the normalized terminal height of rise (THR) Z/Lj and mean spreading elevation (SΗ) ZS/Lj take constant values around 3.5 and 2, respectively. The flow is buoyancy-driven if MN/B < 1, and, in this regime, the normalized THR Z/Lp and SΗ ZS/Lp take values around 4.50 and 3, respectively. The regime 1 < MN/B < 10 is considered to be the transition from plume-like to jet-like flows, where the initial specific momentum and buoyancy fluxes M and B are of equal importance for the development of the flow.
A negatively buoyant jet (
B < 0) or fountain is momentum-driven if
MN/
B > 10. The normalized THR
Z/
Lj and SΗ
ZS/
Lj take asymptotically the same values as momentum-driven positively buoyant flows, which are around 3.50 and 2, respectively. For
MN/
B < 10, the flow is buoyancy-driven, and the normalized THR
Z/
lM takes an asymptotic value around 2 as
MN/
B→1, which is compatible with dimensional analysis and experiments of vertical fountains in a uniform ambient [
21,
31], while the SH is that of the nozzle elevation (
ZS = 0) as
MN/
B→1.5.
Both THR and SΗ were obtained numerically from the equations of motion, using variable entrainment coefficient α that is a function of the local Richardson number of the flow. In buoyancy = driven flows (MN/B < 1), the computed THR and SH were found to be congruent to the measured ones. In momentum-driven flows (MN/B > 10) and fountains (B < 0), both THR and SΗ obtained numerically were not compatible with measurements.
The average dilution obtained numerically in buoyancy-driven jets (
MN/
B < 1) was in accordance with earlier dilution measurements [
10]. However, the computed average dilution in the jet-like flow regime (
MN/
B > 10) was far from the measured minimum dilution [
10], although the experimental data showed substantial scatter. We believe that velocity and concentration measurements using laser light-based tomography of the flow have to be made in order to understand the physics of dispersion of buoyant jets in a density-stratified ambient.
The entrainment coefficient function that was derived from the equations of motion of vertical buoyant jets in a uniform ambient seems to work quite well, especially in regimes where the discharges are plume-like initially or if
MN/
B < 1. Since it is related to the squared local Richardson number of the flow, when the buoyancy acts against the movement, the computed entrainment coefficient is lower than that measured in jets, a result that were previously thoroughly discussed [
18,
19,
34].