Next Article in Journal
Effects of Inflow Condition on RANS and LES Predictions of the Flow around a High-Rise Building
Previous Article in Journal
Experimental and Mathematical Tools to Predict Droplet Size and Velocity Distribution for a Two-Fluid Nozzle
Previous Article in Special Issue
Simplified Modelling of Inclined Turbulent Dense Jets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Vertical Round Buoyant Jets and Fountains in a Linearly, Density-Stratified Fluid

by
Panos N. Papanicolaou
1,* and
George C. Stamoulis
2
1
Applied Hydraulics Laboratory, School of Civil Engineering, National Technical University of Athens, 5 Heroon Polytechniou Street, 15780 Zografou, Greece
2
Hydromechanics and Environmental Engineering Laboratory, Department of Civil Engineering, University of Thessaly, Pedion Areos, 38334 Volos, Greece
*
Author to whom correspondence should be addressed.
Fluids 2020, 5(4), 232; https://doi.org/10.3390/fluids5040232
Submission received: 10 November 2020 / Revised: 30 November 2020 / Accepted: 2 December 2020 / Published: 4 December 2020
(This article belongs to the Special Issue Advances in Turbulent Buoyant Jets)

Abstract

:
Turbulent round buoyant jets and fountains issuing vertically into a linearly density-stratified calm ambient have been investigated in a series of laboratory experiments. The terminal (steady-state) height of rise and the mean elevation of subsequent horizontal spreading have been measured in positively buoyant jets (at source level), including pure momentum jets and plumes, as well in momentum-driven negatively buoyant jets (fountains). The results from experiments confirmed the asymptotic analysis that was based on dimensional arguments. The normalized terminal height and spreading elevation with respect to the elevation of injection of momentum-driven (positively) buoyant jets and fountains attained the same asymptotic values. The numerical results from the solution of entrainment equations, using an improved entrainment coefficient function, confirmed the results related to buoyancy dominant flows (plumes), while their predictions in momentum-driven flows were quite low if compared to measurements.

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 cj (to be defined in Section 2) for positively buoyant, momentum-driven flows (jets) is cj = 3.80, as also presented in [7]. In (positively buoyant) plumes, from the experiments in [1,8,9], [4] proposed the same constant cp ≈ 3.80 (cp 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 L3/T, L4/T2, and L4/T3, respectively, are calculated as follows:
Q = W = π D 2 4 W ,           M = Q W ,           B = g o Q ; g o = ρ a ρ o ρ o g ,
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
N 2 = g ρ o d ρ d z .
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:
L j = M 1 / 4 N 1 / 2   and   L p = B 1 / 4 N 3 / 4 .
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/(goD)1/2, go = (ρ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,
Z N 2 M 1 / 4 = Z L j = c j ,
which must be constant.
(ii)
A plume (M ≈ 0, ρορa), will rise up to a dimensionless elevation,
Z N 3 / 4 B 1 / 4 = Z L p = c p ,
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,
Z L p o r Z L j = f M B N ,
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
Z L p = c j L j L p = c j M B N 1 / 4 ,
meaning that the dimensionless jet rise Z/Lp is proportional to (MN/B) 1/4. Similarly, in plumes (MN/B < 1),
Z L j = c p L p L j = c p M B N 1 / 4 .
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],
Z l M = C ,
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
Z l M = c j L j l M = c j M N B 1 / 2 .
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.,
Z l M = C Z L j = C l M L j = C M N B 1 / 2   or   alternatively   Z L p = C M N B 3 / 4 .
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
μ ( Z ) N 1 / 2 M 3 / 4 = d j   or   μ ( Z ) N 5 / 4 B 3 / 4 = d j M N B 3 / 4 .
If the buoyancy is dominant (MN/B < 1, plume-like flow),
μ ( Z ) N 5 / 4 B 3 / 4 = d p   or   μ ( Z ) N 1 / 2 M 3 / 4 = d p M N B 3 / 4 ,
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ρ/dz 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 cm3/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:
d μ d z = 2 π α m 1 / 2 ,   d m d z = μ β m ,   d β d z = N 2 μ ,
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
μ ( z o ) = Q = π D 2 4 W ,   m ( z o ) = M = π D 2 4 W 2 ,   β ( z o ) = B = ρ a ρ o ρ o g Q .
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.28D. The system of Equation (14) can be integrated starting at zo with the initial conditions in Equation (15).

4. Results

4.1. Positively Buoyant Jet Experiments

The terminal height of rise and spreading height normalized with Lj are plotted in Figure 4a,b, respectively, versus the dimensionless parameter MN/B. For large values of MN/B, i.e., for momentum-driven flows initially, both normalized THR Z/Lj and SH ZS/Lj take constant values around 3.58 ± 0.12 and 1.94 ± 0.11, respectively. These results are congruent with earlier experiments [12,13], while the data spread of ZS/Lj is lower if compared to that of the data by [13]. When the dimensionless parameter MN/B is small, i.e., the flow is initially buoyancy-driven, the normalized THR Z/Lp and SH ZS/Lp are proportional to (MN/B)−1/4 as predicted from dimensional analysis. The normalized THR and SH take the values Z/Lp = 4.52 ± 0.34 and SH ZS/Lp = 2.95 ± 0.34, respectively, in agreement with results reported by the latter authors.

4.2. Negatively Buoyant Jet (Fountain) Experiments

The terminal height of rise and spreading height normalized with Lj are also plotted in Figure 4a,b, respectively, versus MN/B. For large values of MN/B, i.e., for initially momentum-driven flows, both THR and SH take asymptotically constant values Z/Lj ≈ 3.60 and ZS/Lj ≈ 2.0 as in buoyant jets. When MN/B is small, i.e., the flow is initially buoyancy-driven, Z/Lj ≈ 2(MN/B)1/2, while ZS/Lj→−∞. In fact, from Figure 4b, the present data indicate that ZS/Lj→0 when MN/B→1. Our negatively buoyant jet data are limited to MN/B > 1 because, at lower values, the spreading elevation coincides with the receptor boundary that is located upstream from the nozzle at the top layers of the dispersion tank.

4.3. Numerical Results: THR and SH, Average Dilution

The system of Equation (14) with initial conditions in Equation (15) is an initial value problem. The solution proposed is similar to those in [10,14]. Nevertheless, there are differences regarding the variable entrainment coefficient used in the present investigation and the virtual origin. In our computations, we made the assumption that the density gradient in the spreading (waste) field does not vary significantly if compared to the original one. The system of entrainment equations is an initial value problem and was solved using a fourth-order Runge–Kutta routine, starting from the jet virtual origin at elevation zo. The entrainment coefficient at distance z from the source was computed as a function of the local Richardson number Ri(z) = μ(z)β(z)1/2/m(z)5/4 (see Appendix A) from the equation below.
α = C p 2 π 1 + 1 2 R i 2 ( z ) C p .
In positively buoyant jets (B ≥ 0), the jet entrainment coefficient varies from αj = 0.072 to ap = 0.119 (considering a corrected Cp = 0.254 and zo = 3.49D), while, when the buoyancy reverses, Ri2(z) is negative; therefore, a < 0.072. The spreading height can be roughly computed as
N 2 = g ρ o Δ ρ ¯ Δ z Δ z = Z s Z = g ρ o Δ ρ ¯ N 2 = β ( Z ) N 2 μ ( Z ) ,
assuming no further mixing during the descending flow from THR to SH [19], where Δρ, μ(Z) and β(Z) are the average density difference, the specific mass, and buoyancy fluxes, respectively, at the THR.
Note that the findings from computations using top-hat, time-averaged velocity and density difference profiles are not any different from those using Gaussian distributions.
The results from computations regarding the normalized THR and SH are plotted in Figure 4 for both buoyant jets and fountains. For jet-like flows and fountains (MN/B >> 1), the numerical predictions are quite lower than our measurements, since the computed Z/Lj = 3.0 and ZS/Lj = 1.07 are lower than measured values of 3.50 and 2.0, respectively. For initially buoyancy-driven jets (MN/B < 1), the numerical predictions give Z/Lp = 4.25 and ZS/Lp = 3.1. The normalized THR is near the measured one (4.50), while the normalized SH is nearly the same as the predicted SH.
Measurements regarding the dilution of buoyant jets in a linearly stratified uniform ambient were reported by Wong and Wright [10], who measured the dilution of fluorescein dye in the spreading layer. They reported the normalized minimum dilution (that corresponds to maximum time-averaged concentration cm) SmQ/(B3/4/N5/4) versus MN/B. If the initial dye concentration of a jet is Co, then Sm = Co/cm. From tracer conservation equation QCo = μ(z)c (where c is the average dilution), the dimensionless minimum dilution can be expressed as
Q S m N 1 / 2 M 3 / 4 = Q ( C o / c m ) N 1 / 2 M 3 / 4 = ( c / c m ) μ ( z ) N 5 / 4 B 3 / 4 ,
where the ratio c/cm < 1, as the average over a cross-section dye concentration, is only a portion of the maximum time-averaged concentration.
In Figure 5, the normalized minimum dilution from experiments in [10] is plotted along with the computed normalized average dilution μN1/2/M3/4. From this figure, one may note several findings. In plumes (MN/B < 1), the normalized average dilution takes the value dp = 0.87, which is near the 0.84 suggested in [4] and higher than 0.80, an average value of the minimum dilution proposed by [10], as expected. From Equation (16), cm/c = 1.09, a value that is different from 1.40 proposed by [4] for round plumes in uniform ambient. From experiments in [10], the vertical distribution of concentration in the spreading field seems to be uniform. Therefore, an average concentration should be near the maximum, resulting in similar values of minimum and average dilutions. In jets (MN/B > 10), the computed dimensionless average dilution is dj = 0.52, a value that is lower than the minimum dilution 0.67 measured in [4], following the trend of the normalized THR and SH in this regime.

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/lMC [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, √2b (where b is the local jet width) was added in [34]; thus, (Z + √2b)/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].

Author Contributions

Conceptualization, P.N.P.; methodology, P.N.P. and G.C.S.; software, P.N.P.; validation, G.C.S.; investigation, P.N.P. and G.C.S.; data curation, G.C.S.; writing—review and editing, P.N.P.; supervision, P.N.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

Part of this work was presented at the Sixth International Symposium of Environmental Hydraulics, 23–25 June 2010, Athens, Greece. The assistance of Elias Pappas and Demetris Karaberopoulos is gratefully appreciated.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. The Entrainment of Round Vertical Turbulent Buoyant Jets Re-Revisited

Appendix A.1. Introduction

The use of the entrainment coefficient is necessary for computations of buoyant jets using one-dimensional modeling and employing the set of entrainment equations. List, in [4], proposed two equations, a linear one and an empirical exponential, where the asymptotic values of jet and plume entrainment coefficients aj and ap, respectively, are used along with the constant plume Richardson number Rip.
α = α j α j α p R i ( z ) R i p 2 and α = α j exp R i ( z ) R i p 2 ln α p α j .
The equation to the left is attributed to a derivation where the equation of energy was employed [35], while the equation to the right is empirical. Both produce the same asymptotic values for pure jets and plumes. Apparently, to employ these relationships, one must know a priori aj, ap and Rip. The entrainment coefficients for Gaussian distributions of velocity and density difference as suggested in [4] are aj = 0.0535 ± 0.0025 (0.075), aj = 0.0833 ± 0.0042 (0.118), and Rip = 0.557, the values in parentheses corresponding to top-hat distributions of velocity and density difference. In the next section, we derive an equation for the local entrainment coefficient, using the entrainment equations of a vertical, buoyant, axisymmetric jet.
When computing fountains or jets of reverse (with respect to the motion) buoyancy, it was shown [18,34] that the entrainment coefficient is lower than aj. This assumption was also made by [19], where the jet entrainment coefficient was reduced in order to attain numerically the values of measurements of parameters such as the terminal rise height of buoyant jets and fountains in linearly, density-stratified and in uniform ambient, respectively.

Appendix A.2. Derivation of Entrainment Coefficient in Round Buoyant Jets

The entrainment equations of motion of a fully developed, axisymmetric, Boussinesq, turbulent buoyant jet in a uniform and motionless fluid of greater density, using top-hat velocity and density deficit profiles, are written as follows [1,4]:
d μ d z = 2 π α m 1 / 2 ,   d m d z = μ β m ,   d β d z = 0 ,
where α is the local entrainment coefficient, μ, m, and β are the local specific (per unit mass) mass, momentum, and buoyancy fluxes, respectively, written in terms of the time-averaged streamwise top-hat velocity w, the jet width b, and the time-averaged density difference between jet and ambient fluid Δρ as
μ ( z ) = π b 2 w ,   m ( z ) = π b 2 w 2 ,   β ( z ) = π Δ ρ ρ o g b 2 w ,
and ρο is a reference density, such as the jet density at the source. In a motionless uniform ambient, the buoyancy flux is usually conserved (β = B = constant); thus, the momentum equation can be written as
d m d z = μ β m = β m 1 / 2 μ m 1 / 2 m 1 / 2 d m d z = β C p z ,
where Cp is the constant buoyant jet width parameter [29],
C p = μ z m 1 / 2 ,
evaluated from the measurements [30] shown in Figure A1.
Figure A1. Buoyant jet width parameter plotted against the dimensionless elevation (data points from Papanicolaou and List [30]).
Figure A1. Buoyant jet width parameter plotted against the dimensionless elevation (data points from Papanicolaou and List [30]).
Fluids 05 00232 g0a1
Hence, the solution of the ordinary differential Equation (A4) is
d m 3 / 2 d z = 3 2 β C p z m 3 / 2 = 3 4 β C p z 2 + C ,
where C is the integration constant. At z = 0, m = M; therefore, C = M3/2, and the normalized local specific momentum flux m(z)/M as a function of the normalized distance z/lM from the source is
m M = 1 + 3 4 C p z l M 2 2 / 3 ,   l M = M 3 / 4 B 1 / 2 .
From the momentum equation, one may obtain the local volume flux (β = constant) as
d m d z = μ β m μ = 1 β m d m d z = 1 2 β d m 2 d z .
After some algebra, Equations (A6) and (A7) give
μ Q R i o = μ β 1 / 2 M 5 / 4 = C p z l M 1 + 3 4 C p z l M 2 1 / 3 ,   R i o = Q B 1 / 2 M 5 / 4 .
Equation (A8) is a measure of a jet bulk (average) dilution S(z) = μ(z)/Q at distance z from the source. Equations (A6) and (A8) for the computation of the momentum and buoyancy fluxes, respectively, are congruent with the experiments [30], as shown in [36]. Substituting the momentum and volume fluxes from Equations (A6) and (A8) into the continuity equation, one can obtain the buoyant jet entrainment coefficient as follows:
α = 1 2 π m d μ d z = C p 2 π 1 + 1 2 C p z l M 2 1 + 3 4 C p z l M 2 1 .
The local Richardson number is obtained from R i ( z ) = μ β 1 / 2 / m 5 / 4 upon substitution from Equations (A6) and (A8); thus,
R i ( z ) = C p z l M 1 + 3 4 C p z l M 2 1 / 2 .
Therefore, from Equations (A9) and (A10), the entrainment coefficient in terms of the local Richardson number is written as
α = C p 2 π 1 + 1 2 R i 2 ( z ) C p ; R i 2 ( z ) = μ 2 β m 5 / 2 .

Appendix A.3. Verification of Entrainment Coefficient in Round Buoyant Jets

Equations (A6) and (A8) were plotted [36] for comparison with data [30] obtained in dimensionless form. The normalized computed volume flux was congruent with measurements, while the normalized computed momentum flux was somewhat higher from measurements. Using the time-averaged streamwise velocity profile [30] in computing of the local momentum flux, [36,37] indicated that turbulent intensity may contribute as much as 14% and 11%, respectively, of the momentum from time-averaged velocity to the flow. Increasing the momentum flux by an average of 12.5 %, from Figure A1, the adjusted Cp would read Cp = 0.27/1.131/2 = 0.254. Hence, the asymptotic values of flow parameters for vertical axisymmetric pure jets and plumes derived from the analytical solution in the previous section are as follows:
(i)
Momentum flux: When z / l M < < 1 (jet-like flow), Equation (A6) leads to m = M, while, when z / l M > 10 (plume-like flow), it leads to
m M 3 4 C p z l M 2 2 / 3 m 3 4 C p 2 / 3 β 2 / 3 z 4 / 3 = 0.331 β 2 / 3 z 4 / 3 .
(ii)
Volume flux: Similarly, when z / l M < < 1 , from Equation (A8), one may obtain
μ Q = C p R i o z l M μ Q = C p M 5 / 4 Q β 1 / 2 z β 1 / 2 M 3 / 4 μ Q = C p z l Q = C p 4 π z D = C p z l Q ; l Q = π D 4 ,
while, when z / l M > 10 , from Equation (A8) one may obtain
μ Q C p 4 / 3 R i o 3 4 1 / 3 z l M 5 / 3 μ 3 4 1 / 3 C p 4 / 3 β 1 / 3 z 5 / 3 = 0.146 β 1 / 3 z 5 / 3 .
(iii)
Richardson number: In jets (z/lM→0), from Equation (A10), R i ( z ) C p z / l M , while, in plumes, (z/lM→∞) R i ( z ) 4 C p / 3 = 0.58
(iv)
Entrainment coefficient: In momentum jets (lM→∞), α j = C p / 2 π = 0.072 , while, in plumes (z/lM→∞),
α = C p 2 π 1 + 1 2 C p z l M 2 1 + 3 4 C p z l M 2 1 C p 2 π 1 + 2 3 = 5 6 C p π = 0.119 .
The values compare well to those proposed by [4] and to those measured by [30,38], as shown in Table A1. Note that the above values are for top-hat velocity and concentration distributions and must be divided by √2 in order to obtain those for Gaussian distributions.
Table A1. Comparison between measured and computed parameters in vertical jets and plumes using our analytical solution (top-hat formulation).
Table A1. Comparison between measured and computed parameters in vertical jets and plumes using our analytical solution (top-hat formulation).
JetPlumeJetPlumeJetPlumePlume
C p m M m β 2 / 3 z 4 / 3 μ / Q z / l Q μ β 1 / 3 z 5 / 3 α j α p R i p
0.27010.3450.2700.1590.0760.1270.600
0.25410.3310.2540.1460.0720.1190.582
Fischer et al. [4] 0.2500.1500.0760.1180.557
Papanicolaou and List [30] 0.2900.2500.1400.0770.1240.630
Wang and Law [38] 0.0740.1240.584

Appendix A.4. Conclusions

The analysis presented in this Appendix regards a direct derivation of the entrainment coefficient for buoyant jets using the system of entrainment equations, i.e., the mass and momentum conservation. The entrainment coefficient is a function of the jet growth parameter (jet width) only, which is near the value 0.25z, where z is the elevation from the source. It is evaluated using the local Richardson number squared, such that its sign may be positive or negative if the flow is either positively buoyant or a fountain, respectively. This is congruent with derivations [18,34], as well as the assumption [19] that it should be reduced below the jet entrainment coefficient in flows where buoyancy has the opposite sign with respect to the direction of motion.
The local buoyant jet entrainment coefficient computed directly from the local Richardson number using Equation (A11) is analogous to the linear Equation (A1), since the entrainment coefficient is a linear function of Ri2(z), and quite different from the exponential one. The entrainment coefficient computed from Equation (A11) is independent of the sign of buoyancy with respect to the direction of motion. If the buoyancy has opposite sign (negative) from that of the direction of motion, then a is reduced and it takes positive values up to Richardson numbers greater than or equal to −2Cp. For local Richardson numbers lower than −2Cp, the entrainment coefficient becomes negative; therefore, it is meaningless and, in the computation of flows with negative or reversing buoyancy, should be considered zero. Equation (A11) derived here is better to be applied to flows with negative buoyancy instead of Equation (A1) on the left, since, in the latter, the asymptotic plume Richardson number of positively buoyant jets has to be employed, which does not appear in negative or reversing buoyancy. Note that, in negatively buoyant jets, the local Richardson number Ri2(z) has a limit at −∞.

References

  1. Morton, B.R.; Taylor, G.I.; Turner, J.S. Turbulent gravitational convection from maintained instantaneous sources. Proc. R. Soc. Lond. 1956, A234, 1–23. [Google Scholar]
  2. Turner, J.S. Jets and plumes with negative or reversing buoyancy. J. Fluid Mech. 1966, 26, 779–792. [Google Scholar] [CrossRef]
  3. Abraham, G. Jets with negative buoyancy in homogeneous fluid. J. Hydraul.Res. 1967, 5, 235–248. [Google Scholar] [CrossRef]
  4. Fischer, B.H.; List, J.E.; Koh, R.C.Y.; Imberger, J.; Brooks, N.H. Mixing in Inland and Coastal Waters; Academic Press, Inc.: New York, NY, USA, 1979. [Google Scholar]
  5. Fan, L.-N. Turbulent Buoyant Jets into Stratified or Flowing Ambient Fluids; Technical Report, KH-R-15; WM Keck Laboratory of WM Keck Laboratory of Hydraulics & Water Resources, California Institute of Technology: Pasadena, CA, USA, 1967. [Google Scholar]
  6. Fox, D.G. Forced plumes in a stratified fluid. J. Geophys. Res. 1970, 75, 6818–6835. [Google Scholar] [CrossRef]
  7. Chen, C.J.; Rodi, W. Vertical Turbulent Buoyant Jets—A Review of Experimental Data; Pergamon Press: Oxford, UK, 1980. [Google Scholar]
  8. Crawford, T.V.; Leonard, A.S. Observations of buoyant plumes in calm stably stratified air. J. Appl. Meteorol. 1962, 1, 251–256. [Google Scholar] [CrossRef]
  9. Briggs, G.A. Plume Rise; TID 25075; Clearinghouse of Federal Scientific and Technical Information: Springfield, VA, USA, 1969. [Google Scholar]
  10. Wong, D.R.; Wright, S.J. Submerged turbulent jets in stagnant linearly stratified fluids. J. Hydraul. Res. 1988, 26, 199–223. [Google Scholar] [CrossRef]
  11. Abraham, G.; Eysink, W.D. Jets issuing into fluid with density gradient. J. Hydraul. Res. 1969, 7, 145–175. [Google Scholar] [CrossRef]
  12. Papanicolaou, P.N.; Synolakis, C.E.; Hodge, D. The Maximum Height of Rise of Elliptical and Round Buoyant Jets in A Linear Density—Stratified Fluid; School of Engineering, University of Southern California: Los Angeles, CA, USA, 1990. [Google Scholar]
  13. Konstantinidou, K.; Papanicolaou, P.N. Vertical round and orthogonal buoyant jets in a linear density-stratified fluid. In Proceedings of the 30th IAHR Congress, Thessaloniki, Greece, 24–29 August 2003; Theme C. Volume 1, pp. 293–300. [Google Scholar]
  14. Bloomfield, L.J.; Kerr, R.C. Turbulent fountains in a stratified fluid. J. Fluid Mech. 1998, 358, 335–356. [Google Scholar] [CrossRef]
  15. Richards, T.S.; Aubourg, Q.; Sutherland, B.R. Radial intrusions from turbulent plumes in uniform stratification. Phys. Fluids 2014, 26, 036602. [Google Scholar] [CrossRef] [Green Version]
  16. Zhang, W.; He, Z.; Jiang, H. Scaling for turbulent viscosity of buoyant plumes in stratified fluids: PIV measurement with implications for submarine hydrothermal plume. Deep Sea Res. Part I 2017, 129, 89–98. [Google Scholar] [CrossRef] [Green Version]
  17. Carazzo, G.; Kaminski, E.; Tait, S. On the rise of turbulent plumes: Quantitative effects of variable entrainment for submarine hydrothermal vents, terrestrial and extra terrestrial explosive volcanism. J. Geophys. Res. 2008, 113. [Google Scholar] [CrossRef] [Green Version]
  18. Kaminski, E.; Tait, S.; Carazzo, G. Turbulent entrainment in jets with arbitrary buoyancy. J. Fluid Mech. 2005, 526, 361–376. [Google Scholar] [CrossRef] [Green Version]
  19. Papanicolaou, P.N.; Papakonstantis, I.G.; Christodoulou, G.C. On the entrainment coefficient in negatively buoyant jets. J. Fluid Mech. 2008, 614, 447–470. [Google Scholar] [CrossRef]
  20. Wood, I.R.; Bell, R.G.; Wilkinson, D.R. Ocean Disposal of Wastewater; World Scientific: Singapore, 1993. [Google Scholar]
  21. Papanicolaou, P.N.; Kokkalis, T.J. Vertical buoyancy preserving and non-preserving fountains, in a homogeneous calm ambient. Int. J. Heat Mass Transfer 2008, 51, 4109–4120. [Google Scholar] [CrossRef]
  22. Michas, S.N.; Papanicolaou, P.N. Horizontal round heated jets into calm uniform ambient. Desalination 2009, 248, 803–815. [Google Scholar] [CrossRef]
  23. Oster, G.; Yamamoto, M. Density gradient techniques. Chem. Rev. 1963, 63, 257–268. [Google Scholar] [CrossRef]
  24. Papanicolaou, P.N. A piston driven jet for the study of the zone of flow establishment. Exp. Fluids 1994, 17, 287–289. [Google Scholar] [CrossRef]
  25. Mi, J.; Nobes, D.S.; Nathan, G.J. Influence of jet exit conditions on the passive scalar field of an axisymmetric free jet. J. Fluid Mech. 2001, 432, 91–125. [Google Scholar] [CrossRef] [Green Version]
  26. Head, M.J. The Use of Miniature 4-Electrode Conductivity Probes for High Resolution Measurements of Turbulent Density or Temperature Variations in Salt Stratified Water Flows. Ph.D. Thesis, University of California San Diego, San Diego, CA, USA, 1983. [Google Scholar]
  27. Kwon, S.J.; Seo, I.W. Reynolds number effects on the behavior of a non-buoyant round jet. Exp. Fluids 2005, 38, 801–812. [Google Scholar] [CrossRef]
  28. Stamoulis, G. Vertical Round Positively and Negatively Buoyant Jets in a Linearly Density-Stratified Ambient. Diploma Thesis, University of Thessaly, Department of Civil Engineering, Volos, Greece, 2010. (In Greek). [Google Scholar]
  29. List, E.J.; Imberger, J. Turbulent entrainment in buoyant jets and plumes. ASCE J. Hyraul. Div. 1973, 99, 1461–1474. [Google Scholar]
  30. Papanicolaou, P.N.; List, E.J. Investigations of round vertical turbulent buoyant jets. J. Fluid Mech. 1988, 195, 341–391. [Google Scholar] [CrossRef]
  31. Papakonstantis, E.G.; Christodoulou, G.C.; Papanicolaou, P.N. Inclined negatively buoyant jets 1: Geometrical characteristics. J. Hydraul. Res. 2011, 49, 3–12. [Google Scholar] [CrossRef]
  32. Bloomfield, L.J.; Kerr, R.C. A theoretical model of a turbulent fountain. J. Fluid Mech. 2000, 424, 197–216. [Google Scholar] [CrossRef]
  33. Jirka, G.H. Integral model for turbulent buoyant jets in unbounded stratified flows. Part I: Single round jet. Env. Fluid Mech. 2004, 4, 1–56. [Google Scholar] [CrossRef]
  34. Carazzo, G.; Kaminski, E.; Tate, S. The route to self similarity in turbulent jets and plumes. J. Fluid Mech. 2006, 547, 137–148. [Google Scholar] [CrossRef] [Green Version]
  35. Priestley, C.H.B.; Ball, F.K. Continuous convection from an isolated sourse of heat. Q.J.R. Met. Soc. 1955, 81, 144–157. [Google Scholar] [CrossRef]
  36. Papanicolaou, P.N. Analytical solution in vertical buoyant jets. In Proceedings of the 11th HSTAM International Congress on Mechanics, Athens, Greece, 27–30 May 2016. [Google Scholar]
  37. Papanicolaou, P.N. Mass and Momentum Transport in a Turbulent Buoyant Vertical Axisymmetric Jet. Ph.D. Thesis, W.M. Keck Laboratory of Hydraulics and Water Resources, California Institute of Technology, Pasadena, CA, USA, 1984. Report No. KH-R-46. [Google Scholar]
  38. Wang, H.; Law, A.W.-K. Second-order integral model for a round turbulent buoyant jet. J. Fluid Mech. 2002, 459, 397–428. [Google Scholar] [CrossRef]
Figure 1. Normalized terminal rise height (a) Z/Lj and (b) Z/Lp versus the dimensionless parameter MN/B from earlier experiments.
Figure 1. Normalized terminal rise height (a) Z/Lj and (b) Z/Lp versus the dimensionless parameter MN/B from earlier experiments.
Fluids 05 00232 g001
Figure 2. Density profile in the stratified tank up to 37.5 cm deeper than nozzle elevation (Exp-14). The thick line is the measured density, while the darker line is the least square linear approximation ρ(z) = 1006.8 + 0.163z (z in cm).
Figure 2. Density profile in the stratified tank up to 37.5 cm deeper than nozzle elevation (Exp-14). The thick line is the measured density, while the darker line is the least square linear approximation ρ(z) = 1006.8 + 0.163z (z in cm).
Fluids 05 00232 g002
Figure 3. Typical image extracted from the videotape of an experiment of a positively buoyant jet in a linearly density-stratified ambient. Note the double grid in the front and rear Lucite panels of the dispersion tank. The jet plenum edge with nozzle is shown at the top of the photo.
Figure 3. Typical image extracted from the videotape of an experiment of a positively buoyant jet in a linearly density-stratified ambient. Note the double grid in the front and rear Lucite panels of the dispersion tank. The jet plenum edge with nozzle is shown at the top of the photo.
Fluids 05 00232 g003
Figure 4. (a) Normalized terminal rise height Z/Lj and (b) spreading elevation ZS/Lj (bottom) versus MN/B.
Figure 4. (a) Normalized terminal rise height Z/Lj and (b) spreading elevation ZS/Lj (bottom) versus MN/B.
Fluids 05 00232 g004
Figure 5. Normalized dilution plotted versus MN/B for positively buoyant jets from measurements of Wong and Wright [10] (squares and circles) and the computation using top-hat modeling.
Figure 5. Normalized dilution plotted versus MN/B for positively buoyant jets from measurements of Wong and Wright [10] (squares and circles) and the computation using top-hat modeling.
Fluids 05 00232 g005
Table 1. Summary of the asymptotic relationships resulting from dimensional analysis.
Table 1. Summary of the asymptotic relationships resulting from dimensional analysis.
MN/B >> 1
Jets and fountains Z L j = c j Z L p = c j M N B 1 / 4 Z l M = c j M N B 1 / 2 μ ( Z ) N 1 / 2 M 3 / 4 = d j
MN/B < 1
Plumes Z L j = c p M N B 1 / 4 Z L p = c p Z l M = c p M N B 3 / 4 μ ( Z ) N 5 / 4 B 3 / 4 = d p
Table 2. Evaluation of the dimensionless parameters from earlier experiments.
Table 2. Evaluation of the dimensionless parameters from earlier experiments.
AuthorZ/LjZS/LjZ/LpZS/LpdjdpComments
Morton, Taylor, and Turner [1] 1 3.81 ± 0.25 Plume
Fan [5]3.41 ± 0.13 MN/B > 10
Abraham and Eysink [11] 3.63 ± 0.20 MN/B = 0.80–2.40
Wong and Wright [10] for rectangular tank3.72 ± 0.17 4.50 ± 0.46 0.77 ± 0.140.78 ± 0.11MN/B > 10
MN/B < 1
Wong and Wright [10] for circular tank3.23 ± 0.09 4.28 ± 0.37 0.68 ± 0.030.84 ± 0.12MN/B > 10
MN/B < 1
Wong and Wright [10] for all data3.60 ± 0.27 4.46 ± 0.44 0.74 ± 0.120.81 ± 0.11MN/B > 10
MN/B < 1
Papanicolaou et al. [12]3.46 ± 0.16 4.60 ± 0.39 MN/B > 10
MN/B < 1
Bloomfield and Kerr [14] 1,23 ± 0.23 (3.15) *1.53 ± 0.10 (1.52) * B = 0
Bloomfield and Kerr [14] 1,32.88 ± 0.101.35 ± 0.11 B < 0
Konstantinidou and Papanicolaou [13]3.55 ± 0.162.38 ± 0.393.97 ± 0.292.99 ± 0.29 MN/B > 10
MN/B < 1
Richards et al. [15] 1.50 (1.89) * 2.7 (2.78) * MN/B = 0.4 ÷ 46
Zhang et al. [16] 3.47 ± 0.26 MN/B < 2.80
Fischer et al. [4]3.80 3.80 Proposed values
Chen and Rodi [7]3.80 5.00 Proposed values
1 Including virtual origin distance; 2 zero buoyancy flux at source; 3 fountain; * linear fit.
Table 3. Initial flow parameters and measured terminal height of the rise (THR) and spreading height (SH) from positively buoyant jet and fountain experiments [28].
Table 3. Initial flow parameters and measured terminal height of the rise (THR) and spreading height (SH) from positively buoyant jet and fountain experiments [28].
RunDQρaροN2ReRoZZS
(cm)(cm3/s)(gr/cm3)(gr/cm3)(s−2) (cm)(cm)
Positively buoyant jets
Exp-11.0051.261.00731.00920.130371200.02047.324.2
Exp-31.0030.251.00661.00660.135043040.00333.917.7
Exp-41.0021.251.00691.01000.130230230.06133.318.8
Exp-50.5015.241.00571.00620.154143380.00533.317.4
Exp-101.5025.751.00831.01400.107420550.18533.720.7
Exp-140.7527.251.00701.00910.158442400.01937.820.8
Exp-150.5028.751.00751.01000.127672380.00747.825.1
Exp-160.5028.751.00731.00870.126967110.00548.227.7
Exp-201.5021.251.00451.01420.145420630.29436.023.9
Exp-211.5019.751.00311.01500.080919180.35048.631.4
Exp-221.5019.751.00281.01880.070518730.40555.136.6
Exp-231.0024.251.00321.01010.074534500.07951.533.0
Exp-241.0024.251.00701.01000.150534500.05235.120.7
Exp-251.5028.751.00761.01190.152027270.14535.221.5
Exp-261.0033.251.00591.00950.117647310.04238.222.2
Exp-271.0025.751.00651.00840.147334060.03932.818.7
Exp-281.5027.251.00631.01510.119624030.21834.422.0
Exp-291.0036.251.00491.00600.154546790.02137.319.4
Fountains
Exp-21.0036.251.00661.00600.149250360.01538.920.7
Exp-60.5021.251.00711.00320.097860460.01243.324.8
Exp-70.5030.251.00650.99690.122687100.01347.624.8
Exp-80.7530.251.00810.99690.058458070.04034.010.0
Exp-91.0036.251.00890.99860.160443400.06528.49.8
Exp-111.0036.251.00741.00540.137743400.02834.316.7
Exp-120.7528.151.00950.99860.163444930.04228.512.5
Exp-130.5027.251.00980.99880.135263610.01641.420.0
Exp-171.0027.251.01050.99880.129731800.09219.40.3
Exp-181.0026.651.00981.00150.156535260.07922.34.0
Exp-191.0021.251.00580.99750.166228810.09917.91.0
Exp-301.0030.251.00490.99780.104140020.06525.37.6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Papanicolaou, P.N.; Stamoulis, G.C. Vertical Round Buoyant Jets and Fountains in a Linearly, Density-Stratified Fluid. Fluids 2020, 5, 232. https://doi.org/10.3390/fluids5040232

AMA Style

Papanicolaou PN, Stamoulis GC. Vertical Round Buoyant Jets and Fountains in a Linearly, Density-Stratified Fluid. Fluids. 2020; 5(4):232. https://doi.org/10.3390/fluids5040232

Chicago/Turabian Style

Papanicolaou, Panos N., and George C. Stamoulis. 2020. "Vertical Round Buoyant Jets and Fountains in a Linearly, Density-Stratified Fluid" Fluids 5, no. 4: 232. https://doi.org/10.3390/fluids5040232

APA Style

Papanicolaou, P. N., & Stamoulis, G. C. (2020). Vertical Round Buoyant Jets and Fountains in a Linearly, Density-Stratified Fluid. Fluids, 5(4), 232. https://doi.org/10.3390/fluids5040232

Article Metrics

Back to TopTop