In fluid mechanics, fountains take place when a source fluid is driven by its own momentum into a surrounding ambient fluid, and it is counterbalanced by buoyancy. These phenomena are largely encountered in nature and human activities. Despite the numerous studies on the subject, few experimental data are available about the internal structure of turbulent fountains. Here, we present a set of laboratory experiments with the aim to (i) get direct velocity and density measurements of fountains in a controlled environment and (ii) obtain insights about the basic physics of the phenomenon. The results concern the characteristics of the mean and turbulent flow: we report the analysis of the turbulent kinetic energy, the velocity skewness, and the Reynolds stresses, including a quadrant analysis of the fluctuating velocities. For some tests, the correlation between density and vertical velocity is investigated for both mean and fluctuating values. We have quantified the momentum transport, which is mainly out-downward at the nozzle axis with peaks at the mean rise height, where also maximum levels of the buoyancy and mass fluxes are present. The ability of acoustic Doppler current profilers to identify the rise height of the fountain and to measure the velocity field is also discussed.

## I. INTRODUCTION

Turbulent fountains are generally defined as localized vertical flows of a source fluid into an ambient fluid with different density. The result is a jet with opposing buoyancy invested by a return flow.^{1}

The research activity on fountains is justified by their widespread occurrence and range of practical applications. One of the main examples is the role of fountains in heating and cooling within buildings (e.g., in air conditioning).^{2,3} Fountains form as cool air is mechanically injected upward through a floor-level cooling diffuser^{4} and as warm air is injected downward to form curtains, which are commonly used in tunnels and shop entrances as a means of segregating regions of fluid.^{5} Fountains in the built environment also include those that may form more naturally (e.g., during an enclosure fire).

The interest in fountains is also strong in geophysical sciences and environmental engineering. When a cloud tower is growing upwards into a dry environment, the evaporation of liquid water near the edge causes cooling and hence a buoyancy inversion; the result is the formation of heavier fluid, which drives the flow down again.^{6} Other examples of natural fountains are the evolution of volcanic eruption columns,^{7} and the replenishment of magma chambers in the earth's crust (through the cyclic intrusion of pulses of dense magma that give rise to fountain-like flows).^{8,9} Hunt and Burridge^{1} present a detailed review of many other applications that have been studied in the literature.

Different classes of fountains exist, and they can be defined depending on the source Froude number, $Fr0=w0/g0\u2032r0$, where *w*_{0} is the velocity at which fluid is ejected from the source, *r*_{0} is the radial scale for the source, and $g0\u2032$ is the buoyancy of the source fluid defined as $g0\u2032=g(\rho 0\u2212\rho a)/\rho a$, where *ρ*_{0} and *ρ _{a}* are the densities of the source and ambient fluid, respectively. A typical classification is the one proposed by Kaye and Hunt,

^{10}extended by Burridge and Hunt

^{11}and reported in Table I.

$0.3\u2272Fr0\u22721.0$ | Very weak fountains |

$1.0\u2272Fr0\u22722.0$ | Weak fountains |

$2.0\u2272Fr0\u22724.0$ | Intermediate fountains |

$Fr0\u22734.0$ | Forced and highly forced fountains |

$0.3\u2272Fr0\u22721.0$ | Very weak fountains |

$1.0\u2272Fr0\u22722.0$ | Weak fountains |

$2.0\u2272Fr0\u22724.0$ | Intermediate fountains |

$Fr0\u22734.0$ | Forced and highly forced fountains |

Furthermore, fountains may be regarded as laminar for source Reynolds numbers $Re0\u2272120$, turbulent for $Re0\u22732000$, and transitional for $120\u2272Re0\u22722000$, where $Re0=w0r0/\nu $ with the representative kinematic viscosity, *ν*, typically taken as that of the source fluid. In addition, Burridge, Mistry, and Hunt^{12} found that the threshold Reynolds number, *Re _{T}*, separating transitional to turbulent regime, is not constant and depends on the Froude number (with lower

*Re*at lower

_{T}*Fr*

_{0}), and they proposed $ReT=75\u2009Fr0+350$ for $Fr0>2$.

The present work focuses on forced fountains. The dynamics of such fountains is characterized by a first pulse of fluid, that is, a starting plume with a vortex-like front and nearly steady plume behind.^{6} Afterward, the plume broadens, comes to rest, and fell back. In fact, the fluid initially rises before the opposing buoyancy force arrests the flow and subsequently induces a returning counterflow (rise and fall behavior). Finally, the fountain settles down to a nearly steady state, with an up-flow in the center and a down-flow surrounding this. It is worth mentioning that the maximum distance from the source is reached by the first pulse. Figure 1 shows a schematic illustration of a forced fountain at the initial stage and during the subsequent steady state.

The approach to the subject is mainly experimental, with most of the literature studies regarding the steady, vertical, upwards injection of a heavy salt solution into a freshwater tank.^{2,6} Several variants have been proposed, including the use of aqueous potassium chloride (KCl) solutions^{13} and glycerol–water mixtures;^{14} also downward ejections of positively buoyant source fluid have been generated, with jets of heated water into cooler water^{15} or warm air into cool air.^{16} A first theoretical approach was proposed by the pioneering work of Morton, Taylor, and Turner,^{17} and self-similar solutions have been further developed since then.^{18–20} Self-similarity is widely encountered in the study of buoyancy- and gravity-driven phenomena, including gravity currents and non-Newtonian flows.^{21}

Most of the experimental data available in the literature are represented by measurements of the bulk flow, mainly obtained with image processing techniques.^{1} In the case of vertical turbulent jet with negative buoyancy, velocity profile measurements were obtained using a hot film anemometer,^{15} while the flow structures were visualized by means of a Particle Image Velocimetry (PIV) system for transient positive and negative buoyant fountains.^{13} More recently, PIV and planar laser induced fluorescence (PLIF) have been used to simultaneously measure the velocity and scalar concentration fields^{22} and to give a first description of the turbulent structure.^{23} However, additional experimental datasets and interpretations are needed to understand (i) the complex interaction between the upflow and the counterflow and (ii) the phenomena related to the momentum and mass exchanges.

In this work, we present experiments of forced fountains generated by the injection of dyed salt-water in homogeneous fresh water. Vertical and horizontal velocity profiles are acquired using an Acoustic Doppler Current Profiler (ADCP), and for some tests, also the vertical density profile is measured by means of a conductivity sensor. The aim of the paper is (i) to give further details about the mean flow and turbulence in fountains with negative buoyancy and (ii) to discuss the performance of the adopted instruments and techniques.

The paper is structured as follows. In Sec. II, the experiments are described: facility and instrumentation are illustrated along with the data processing and methodology. Experimental results and their discussion are reported in Sec. II B, including the analysis of the turbulent kinetic energy (*TKE*), the velocity skewness, the Reynolds stresses, and a quadrant analysis of the fluctuating velocities. Main conclusions are summarized in Sec. III.

## II. STUDY CASE AND METHODOLOGY

The present section describes the experimental setup and the experimental program. In addition, processing techniques and physical quantities are introduced, and some representative parameters of the tests are compared with the literature data for a better overview of the study case.

### A. Experimental facility and program

A series of experiments have been performed at the Hydraulics Laboratory of the University of Parma (Italy). As stated above, the experimental activity aimed to reproduce a vertical plume of a denser fluid which propagates into lighter fluid, and thus subjected to a negative buoyancy. In these tests, we used dyed brine for the denser fluid forming the plume and homogeneous fresh water for the ambient fluid.

The experimental apparatus consists of a square-section tank with dimensions $400\xd7400\xd7800\u2009mm3$, as shown in Fig. 2. A vertical rigid tube with internal diameter $Dint=7.8\u2009mm$ is fixed at the bottom of the tank, and it protrudes upward for a length of 300 mm. The tube is connected to an external pump, which allows the generation of fountains by injecting the salt water into the freshwater tank. A proportional-integral-derivative controller (PID) was used to control the flow rate, which was measured with a turbine flow-meter. The tank was filled with fresh water up to 600–650 mm before starting the experiments.

In a first set of experiments, a video camera with a resolution of 2 MP (1920 × 1080 pixels) was used to detect the interface between the ambient fluid (fresh-water) and the vertical plume (dyed salt-water), and hence to determine the main statistics of the fountain rise height. Before testing, a grid with known coordinates was inserted inside the tank and recorded for the extrinsic calibration of the camera, in order to transform the coordinate system from pixels to meters. During the test, the grid was removed to avoid disturbances to the flow.

A variable number of acoustic Doppler current profilers (DOP2000 by Signal Processing S.A.), hereinafter referred to as ADCPs, have been placed inside the tank to measure velocity in the vertical (*z*) and horizontal (*x*) directions. The ADCP averages data within control volumes (gates) at incremental distances from the probe, providing instantaneous velocity profiles with a rate ≈20 Hz. Notice that such a sampling frequency makes it possible to observe turbulent structures at a scale that is not affected by the viscosity and therefore that is substantially independent of the Reynolds number. The instrument we used is monostatic (i.e., it acts like a transceiver), and it is controlled by a computer, which allows the user to define a range of settings. We set an acoustic wave carrier with frequency of 8 MHz, a velocity measurements range of $\xb1320\u2009mm\u2009\u2009s\u22121$, and a spatial resolution of 1.5 mm, which determined the spatial range 0–100 mm starting from the probe. The estimated beam divergence angle is ≈2°, and the probe diameter is equal to 8 mm.

In order to filter the measurements, we disregarded velocity values with a number of echoes $N<Nt$, where the threshold value was taken as $Nt=N\xafe/3$ and $N\xafe$ is the time-averaged number of echoes (such a filtering is carried out independently for each value of *z*). In this way, we removed data with a poor backscatter, which could increase the experimental uncertainty. Notice that (i) the injected fluid was seeded with TiO_{2} particles, characterized by high sonic impedance, and (ii) a poor backscatter is associated with low (or null) tracer concentration in the ambient fluid. We calculated the Stokes number of the particles, *Stk*, which is an indicator of the fidelity of the flow tracers in turbulent flows, lower than 0.1, so the expected error due to the tracer is less than 1%.^{24} For this reason, no specific correction was applied to the velocity data. Vice versa, because the speed of sound depends on the density and temperature of the fluid, and to avoid errors of the order of 5%, we have corrected the position of the gate and the particle velocity using a model for density-bulk modulus-salinity suggested by Mackenzie.^{25}

A conductivity probe (Conduino) was installed together with the ADCP during some experiments. The primary sensor is represented by two pins (micro USB type B connectors) that work as electrodes spaced ≈0.2 mm. The volume of measurement is a cylinder of approximate height 4 mm and radius 2 mm, and the data rate is $\u224820\u2009Hz$. The voltage output is proportional to the fluid salinity which, in turn, gives the instant value of the density in a point. Further details on this type of instrumentation and its applications can be found in Petrolo and Longo.^{26}

The configurations of the different experimental conditions are shown in Fig. 3. For a first set of experiments [expts. 1–5, Fig. 3(a)], the video camera was active, and a single ADCP was installed at a fixed position measuring the vertical velocity profile above the inflow section. Then, for two experiments [expts. 6–7, Fig. 3(b)] an ADCP was mounted together with the Conduino on a traverse system, which continuously moved up and down in the tank during the experiments. This moving support covered the entire extension of the fountain with a velocity $\u22486\u2009mm\u2009s\u22121$, allowing to obtain the profile of both vertical velocity and density. The Conduino was vertical, aligned with the source of the fountain, while the ADCP was mounted by the side with an inclination of 20°. Finally, three ADCP were installed on the moving support [Fig. 3(c)], two of which measuring vertical velocities (ADCP1 above the inflow section and ADCP2 with a horizontal offset of 10 mm) and one measuring horizontal velocity (ADCP3). In particular, expts. 8–10 were realized using ADCP1 and ADCP3 in movement; expts. 11–14 still involved ADCP1 and ADCP3 but in a fixed position; and expts. 15–20 were performed using all three probes in movement.

A linear potentiometer was connected to the traverse system and used to measure the position in time of the probes. An external trigger was used to start the experiments and the data acquisition (video, ADCP and Conduino).

The main parameters of the experiments are listed in Table II. In our experiments, the internal radius was $r0=Dint/2=3.9\u2009mm$, the source fluid density was in the range $1020\u2264\rho 0\u22641092\u2009kg\u2009m\u22123$, and a value of source discharge $Q=15\u2009ml\u2009s\u22121$ was used, which yields the source velocity $w0=Q/(\pi r02)=314\u2009mm\u2009s\u22121$. The ranges of the non-dimensional groups are $5.4\u2009\xb1\u20090.4\u2264Fr0\u226416.4\u2009\xb1\u20091.4$ and $1030\u2009\xb1\u200970\u2264Re0\u22641180\xb180$ at the inflow section, which indicate that we are dealing with (highly) forced fountains, following the classification by Burridge and Hunt^{11} and Burridge and Hunt,^{27} herein BH2012 and BH2013, respectively. By adopting the relation proposed by Burridge, Mistry, and Hunt^{12} (reported in Sec. I), we find that all the fountains generated by the present activity can be considered turbulent, except for two tests (expts. 7 and 20), which are very close to the threshold and, in any case, far from the laminar conditions.

Expts. # . | Video . | ADCP . | Q $(ml\u2009s\u22121)$
. | r_{0} (mm)
. | w_{0} $(mm\u2009s\u22121)$
. | ρ_{0} $(kg\u2009m\u22123)$
. | Re_{0}
. | Re
. _{T} | Fr_{0}
. |
---|---|---|---|---|---|---|---|---|---|

1 | Active | 1 fixed | 15.3 | 3.9 | 320 | 1030 | 1180 | 1059 | 9.5 |

2 | Active | 1 fixed | 15.3 | 3.9 | 320 | 1051 | 1130 | 894 | 7.3 |

3 | Active | 1 fixed | 15.3 | 3.9 | 320 | 1069 | 1080 | 817 | 6.2 |

4 | Active | 1 fixed | 15.3 | 3.9 | 320 | 1092 | 1030 | 755 | 5.4 |

5 | Active | 1 fixed | 15.3 | 3.9 | 320 | 1089 | 1040 | 762 | 5.5 |

6 | No | 1 moving + Cond. | 15.0 | 3.9 | 314 | 1028 | 1160 | 1069 | 9.6 |

7 | No | 1 moving + Cond. | 15.0 | 3.9 | 314 | 1021 | 1180 | 1181 | 11.1 |

8 | No | 2 moving | 15.0 | 3.9 | 314 | 1060 | 1080 | 841 | 6.6 |

9 | No | 2 moving | 15.0 | 3.9 | 314 | 1050 | 1110 | 888 | 7.2 |

10 | No | 2 moving | 15.0 | 3.9 | 314 | 1040 | 1130 | 952 | 8.0 |

11 | No | 2 fixed | 15.0 | 3.9 | 314 | 1060 | 1080 | 841 | 6.6 |

12 | No | 2 fixed | 15.0 | 3.9 | 314 | 1050 | 1110 | 888 | 7.2 |

13 | No | 2 fixed | 15.0 | 3.9 | 314 | 1040 | 1130 | 952 | 8.0 |

14 | No | 2 fixed | 15.0 | 3.9 | 314 | 1030 | 1160 | 1045 | 9.3 |

15 | No | 3 moving | 15.0 | 3.9 | 314 | 1070 | 1060 | 805 | 6.1 |

16 | No | 3 moving | 15.0 | 3.9 | 314 | 1060 | 1080 | 841 | 6.6 |

17 | No | 3 moving | 15.0 | 3.9 | 314 | 1050 | 1110 | 888 | 7.2 |

18 | No | 3 moving | 15.0 | 3.9 | 314 | 1040 | 1130 | 952 | 8.0 |

19 | No | 3 moving | 15.0 | 3.9 | 314 | 1030 | 1160 | 1045 | 9.3 |

20 | No | 3 moving | 15.0 | 3.9 | 314 | 1020 | 1180 | 1201 | 11.4 |

Expts. # . | Video . | ADCP . | Q $(ml\u2009s\u22121)$
. | r_{0} (mm)
. | w_{0} $(mm\u2009s\u22121)$
. | ρ_{0} $(kg\u2009m\u22123)$
. | Re_{0}
. | Re
. _{T} | Fr_{0}
. |
---|---|---|---|---|---|---|---|---|---|

1 | Active | 1 fixed | 15.3 | 3.9 | 320 | 1030 | 1180 | 1059 | 9.5 |

2 | Active | 1 fixed | 15.3 | 3.9 | 320 | 1051 | 1130 | 894 | 7.3 |

3 | Active | 1 fixed | 15.3 | 3.9 | 320 | 1069 | 1080 | 817 | 6.2 |

4 | Active | 1 fixed | 15.3 | 3.9 | 320 | 1092 | 1030 | 755 | 5.4 |

5 | Active | 1 fixed | 15.3 | 3.9 | 320 | 1089 | 1040 | 762 | 5.5 |

6 | No | 1 moving + Cond. | 15.0 | 3.9 | 314 | 1028 | 1160 | 1069 | 9.6 |

7 | No | 1 moving + Cond. | 15.0 | 3.9 | 314 | 1021 | 1180 | 1181 | 11.1 |

8 | No | 2 moving | 15.0 | 3.9 | 314 | 1060 | 1080 | 841 | 6.6 |

9 | No | 2 moving | 15.0 | 3.9 | 314 | 1050 | 1110 | 888 | 7.2 |

10 | No | 2 moving | 15.0 | 3.9 | 314 | 1040 | 1130 | 952 | 8.0 |

11 | No | 2 fixed | 15.0 | 3.9 | 314 | 1060 | 1080 | 841 | 6.6 |

12 | No | 2 fixed | 15.0 | 3.9 | 314 | 1050 | 1110 | 888 | 7.2 |

13 | No | 2 fixed | 15.0 | 3.9 | 314 | 1040 | 1130 | 952 | 8.0 |

14 | No | 2 fixed | 15.0 | 3.9 | 314 | 1030 | 1160 | 1045 | 9.3 |

15 | No | 3 moving | 15.0 | 3.9 | 314 | 1070 | 1060 | 805 | 6.1 |

16 | No | 3 moving | 15.0 | 3.9 | 314 | 1060 | 1080 | 841 | 6.6 |

17 | No | 3 moving | 15.0 | 3.9 | 314 | 1050 | 1110 | 888 | 7.2 |

18 | No | 3 moving | 15.0 | 3.9 | 314 | 1040 | 1130 | 952 | 8.0 |

19 | No | 3 moving | 15.0 | 3.9 | 314 | 1030 | 1160 | 1045 | 9.3 |

20 | No | 3 moving | 15.0 | 3.9 | 314 | 1020 | 1180 | 1201 | 11.4 |

### B. Physical quantities and scales

The main quantities that characterize the rise height and the rhythm of the fountains are the quasi-steady rise height *z _{ss}*, the fountain width $2b\u0303ss$, the mean rise height peak $zpe\xaf$, and the mean rise height trough $ztr\xaf$. The analytical values of the defined quantities are as follows:

where *T* is the acquisition time, $zf(t)$ the instantaneous value of the fountain height, $b\u0303(zss,t)=b(zss,t)$ when $zf(t)\u2265zss$ and $b(zss,t)$ is the fountain half-width at $z=zss$. Furthermore: (i) $zpe(t)=zf(t)$ when $zf(t)\u2265zss+\sigma ss$, (ii) $ztr(t)=zf(t)$ when $zf(t)\u2264zss\u2212\sigma ss$, and (iii) $\delta th,\u2009\delta tpe$, and $\delta ttr$ are the total periods for which $zf(t)\u2265zss,\u2009zf(t)\u2265zss+\sigma ss$, and $zf(t)\u2264zss\u2212\sigma ss$, respectively. The term *σ _{ss}* represents the standard deviation of the vertical fluctuation

*z*over the acquisition time

_{ss}*T*. In addition, we define the magnitude of the vertical fluctuations as $\delta zss=zpe\xaf\u2212ztr\xaf$.

Figure 4 shows the comparison between the rise height statistics found in BH2012 and in the present work. The non-dimensional quasi-steady rise height, $zss/r0$, and the magnitude of the vertical fluctuations, $\delta zss$ (scaled both with the width of the forced fountains, $\delta zss/2\u2009b\u0303ss$, and with the quasi-steady rise height, $\delta zss/zss$), are well aligned. As found in BH2012, it suggests that (i) the forced fountains scale as $zss\u221dr0\u2009Fr0$ ($zss=2.22\u2009r0\u2009Fr0$ by fitting the present experiments), and (ii) the height of the vertical fluctuations is of the same order of the large-scale eddies at the fountains top and it is independent of *Fr*_{0}.

The fluctuations of the fountain top were found in BH2013 to be prominently bi-chromatic in the frequency domain. Thus, it is possible to define a Strouhal number of the higher peak frequency, *f _{H}*, and of the lower peak frequency,

*f*, as $StH=fH\u2009T0$ and $StL=fL\u2009T0$, respectively, where

_{L}*T*

_{0}denotes an adequate timescale. A conventional timescale is represented by $r0/w0$, which results in a Strouhal number $St\u221dFr0\u22122$. A timescale proposed for forced fountains was $w0/g\u2032$, for which the corresponding forced Strouhal number,

*St*, appears independent of

_{for}*Fr*

_{0}. Another relevant timescale comes from the large-eddies length and velocity scales, i.e., $\u221d2b\u0303ss/wss$, where

*w*is the root mean square (RMS) vertical velocity of the fountain top (calculated as the time derivative of the interface signal). The latter timescale is associated with the top Strouhal number,

_{ss}*St*.

_{top}Figure 5 shows the comparison between the values of the non-dimensional found in BH2012 and in the present work. The results show a good overlap of the conventional Strouhal numbers *St _{H}* and

*St*, as well as for the forced and top Strouhal numbers

_{L}*St*and

_{for}*St*(related to

_{top}*f*).

_{H}In our experiments, we consider the quasi-steady rise height *z _{ss}* as the vertical length scale and the fountain width $b\u0303ss$ for the horizontal length scale, while the vertical and horizontal velocities are non-dimensional with the source flow velocity

*w*

_{0}.

In this study, we are not interested in the early stage of the fountain (negatively buoyant jet), which has been widely studied in recent works.^{23,28} Here, we want to focus on the steady state of the forced fountain, when the up- and counter-flow are both present and interact with each other.

### C. Detection of the fountain interface

A well known method to extract the interface position in experimental fountains is through the use of a video camera,^{29} as we also did in our experiments. The image analysis we used mainly follows the same procedure reported by BH2012 and BH2013. A MATLAB script (i) extracts the pixel array above the inflow midsection at each instant, and (ii) concatenates successive arrays to build a resulting image, which represents the temporal evolution of the fountain interface along *z*. Moreover, the sharp density interface between the salt water and the fresh water is responsible for a net discontinuity (a “jump”) in the echoes number of the ADCP1 signal, which indicates the instantaneous position of the interface. Thus, we retrieved the time series of the interface fluctuations also by following the signal of the ADCP1 echoes in time (expts. 1–5). The video frames are also used to get the instantaneous vertical velocity of the interface, *w _{v}*. On the other hand, we used the ADCP1 measurements to extract the vertical velocity

*w*at

_{ADCP}*z*, that is the mean rise height.

_{ss}Since we have both video and ADCP measurements, we can compare the results of the two techniques in terms of temporal and spectral signals. As an example, we show the results for expts. 1 and 5, corresponding to a source flow density $\rho 0=1030\u2009kg\u2009m\u22123$ and $\rho 0=1089\u2009kg\u2009m\u22123$. Panels (a) and (b) in Fig. 6 show the temporal evolution of the interface. Both the positive and the negative peaks of the signals are well-individuated by the ADCP-extracted interface, especially in the case of the lower density, which is characterized by a greater amplitude (and period) of the fluctuations. The results of the comparison of the two techniques are also reported in Table III, for both the mean value and standard deviation of the rise height. The difference between the average values is in the range $1.3\u22122.9\u2009mm$, while for the standard deviation (which represents fluctuations) it is $0.2\u22120.6\u2009mm$. The discrepancies are within the experimental uncertainty. Panels (c) and (d) show the time series of the velocities, *w _{v}* and

*w*. The agreement between video and ADCP data is acceptable. Some discrepancies are present, and they can be explained by considering that the comparison is made between the velocity of the moving interface (video data) and the velocity measured at the mean rise height (ADCP). We also see that fluctuations of

_{ADCP}*w*decrease with decreasing

_{ADCP}*Fr*

_{0}. Finally, the spectral analysis of signals is reported in panels (e) and (f) of Fig. 6. The dominant low frequency,

*f*, and the dominant high frequency,

_{L}*f*, of the fluctuations signal are calculated according to BH2013, and their value is compared with the power spectral density (PSD) of

_{H}*w*and

_{video}*w*. The results qualitatively show that the dominant frequencies well represent the peak frequencies for

_{ADCP}*w*as expected and also capture the main peaks in the spectra of the velocity

_{v}*w*.

_{ADCP}Expt. no . | ρ_{0} (kg m^{−3})
. | $zss,v\xaf$ (mm) . | $zss,D\xaf$ (mm) . | $\sigma zss,v$ (mm) . | $\sigma zss,D$ (mm) . |
---|---|---|---|---|---|

1 | 1030 | 81.2 | 82.5 | 4.6 | 4.2 |

2 | 1051 | 60.1 | 63 | 3.3 | 2.7 |

3 | 1069 | 52.4 | 54.0 | 3.0 | 3.2 |

5 | 1089 | 48.4 | 50.0 | 3.0 | 3.4 |

Expt. no . | ρ_{0} (kg m^{−3})
. | $zss,v\xaf$ (mm) . | $zss,D\xaf$ (mm) . | $\sigma zss,v$ (mm) . | $\sigma zss,D$ (mm) . |
---|---|---|---|---|---|

1 | 1030 | 81.2 | 82.5 | 4.6 | 4.2 |

2 | 1051 | 60.1 | 63 | 3.3 | 2.7 |

3 | 1069 | 52.4 | 54.0 | 3.0 | 3.2 |

5 | 1089 | 48.4 | 50.0 | 3.0 | 3.4 |

The overall results suggest that, in order to trace the interface of two fluids with slightly different densities, the use of ADCP is comparable to the current detection methods, and it could be a good alternative to the use of a video camera. In particular, this avoids the storage of a large amount of data and the subsequent image processing.

## III. RESULTS AND DISCUSSION

The present section describes the results of measurements and data processing, with the aim of characterizing the flow field in the fountain (both in axis with the jet emission and along a vertical that is $10\u2009mm$ away from the same axis). The mean flow and turbulence are analyzed, and for tests 6–7, also the density profile and fluctuations are taken into account.

### A. The mean flow

The measured velocity is decomposed as $v(a,t)=v\xaf(a)+v\u2032(a,t)$, where *v*(*t*, *a*) is the instantaneous velocity along the measuring direction *a*, while $v\xaf(a)$ and $v\u2032(t,a)$ are the mean and fluctuating components, respectively. For horizontal velocity *a* =* x* and *v* =* u*, while for vertical velocity *a* =* z* and *v* =* w*. The mean velocity profile is obtained by time-averaging the ADCP signal

At each point along the vertical, the ADCP provides the measurement as an average on a circular footprint (disk), which slightly enlarges as the distance from the transducer increases. The divergence angle is $\u22482\xb0$, which means that at the farther limit of the range ($10\u2009cm$) the footprint radius is $\u22483.5\u2009mm$ bigger than the ultrasound source. On the opposite, the fountain widens away from the nozzle outlet (upwards). Notice that the nozzle diameter is approximately equal to the ADCP probe transceiver, thus there must be a point in the vertical where the fountain width, $2bu(z)$, equals the footprint diameter, $dADCP(z)$. Where $dADCP(z)>2bu(z)$, the measure is not reliable because it is the result of the interaction between upflow and counterflow in the footprint. This explains why we have limited the presentation of almost all the results to the lower limit of $z/zss\u22480.4$. This only applies for ADCP1.

Figure 7(a) shows the mean vertical velocity profile of ADCP1 that is aligned to the fountain source (i.e., axial measurement). Data collapse fairly well on a single curve, with null velocity at $z/zss\u22481.2$, and with a linear trend down to $z/zss\u22480.4$. Beyond this limit, the vertical velocities show a non monotonic profile when approaching the source inflow section; this is a non physical behavior that can be explained if we keep in mind the operating principle of the probe (described above). The comparison between present experiments and the results by Mizushina *et al.*^{15} is reported in Fig. 7(b): away from the source the agreement with literature data is within the experimental uncertainty.

Figure 7(c) shows the mean vertical velocity profile of ADCP2 (non-axial measurements). Velocities are slightly negative above $z/zss\u22481.1$, indicating the presence of the counterflow and/or of a current induced by the counterflow itself. Lower down, the behavior is strictly related to the density of the injected fluid. For higher densities, the measurement volumes are entirely within the counterflow, with negative velocity values decreasing downward. Vice versa, in the case of lower densities, the upflow widens more, and the probe registers positive velocities. Then, when the counterflow expands and invades the region next to the inlet pipe, the measures are negative again. This offers an indirect measurement of the shape of the counterflow.

Figure 7(d) shows the mean horizontal velocity registered by ADCP3, at fixed positions. These results show how far the fountain effects are felt in terms of induced currents and recirculation. Regardless of the density of the jet, the flow field extends at least up to $x/bss\u2248\xb15$ (herein, the symbol $\u0303$ over *b _{ss}* is omitted for simplicity). Moreover, results from expts. 8–10 and 15–20, with the probes moving up and down, allow reconstructing the horizontal velocity map. Figure 8(a) shows the results for test 16 (3 moving ADCP), which is representative for all ADCP3 data. Just above the nozzle the velocities are inward, due to the drag effect that draws ambient fluid from the surrounding areas. This favors mixing, even if the ambient fluid is only involved in the early stage of the fountain, while later it is the turn of counterflow fluid. In the upper part, there is a substantial symmetry, with the flow directed toward the outside of the fountain itself. Figures 8(b) and 8(c) show, respectively, the standard deviation and the skewness of the horizontal velocity, and they will be discussed below. Notice that Fig. 8 is a merge of two tests performed in the same experimental conditions, in order to cover the whole height of the fountain (Expts. 8 and 16).

### B. Turbulent kinetic energy

The root mean square (RMS) of the fluctuating velocity is defined as follows:

where the overline indicates herein the time average of the argument, and the fluctuating velocity is obtained by subtracting the mean value to the whole signal, $v\u2032=v\u2212v\xaf$. The quantity $vrms\u2032$ also represents the velocity standard deviation, and it is strictly related to the two-component turbulent kinetic energy (*TKE*), which is calculated as

Figures 9(a) and 9(b) report the standard deviation of the vertical velocity for ADCP1 and ADCP3. In the range $0.4\u2272z/zss\u22721,\u2009wrms\u2032$ is almost constant, with values between 0.1 and 0.16 for probe ADCP1, and between 0.8 and 0.12 for probe ADCP2. In the case of axial measurements, some tests (especially those characterized by the higher densities) present a peak at $z/zss\u22481.1$, where the amplitude of the fluctuations is maximum. As expected, the trend of the series is then slightly decreasing upwards. In the case of ADCP3 (non-axial measurements), a density-dependent trend can be observed for $z\u22731.1$, with more intense fluctuations for higher densities. Figure 8(b) shows the map of $urms\u2032$, with a magnitude of the order of $wrms\u2032$ and higher values in the inner part of the upflow. Figure 10 shows the profile of the turbulent kinetic energy calculated with measurements from ADCP1 and ADCP3. Data are slightly dispersed, but it is possible to observe a common trend with a maximum of TKE at $z/zss\u22481$. This suggests that the turbulence is mainly developed in the upper part of the fountain (at the mean rise height), where the flow from the nozzle collides with the flow generated by the periodic collapses of the plume (rise and fall behavior). Then, TKE decreases upwards and becomes almost null at $z/zss\u22481.4$.

### C. Skewness

The statistics of turbulence can also be characterized by the velocity skewness, which is an indicator of the probability density function (PDF) symmetry with respect to a Gaussian distribution. For a normal distribution, the skewness is zero. Negative values indicate that the signal distribution peak is shifted toward the right tail of the PDF, while positive values indicate that the signal distribution peak is shifted toward the left tail. The velocity skewness represents the third central moment of the velocity signal and it is calculated as

for vertical velocities. The indicator *s _{w}* gives also information on the structure of the flow field,

^{30}since the triple correlation $w\u20323\xaf$ represents the transport of $w\u2032w\u2032\xaf$ by the turbulence itself. Moreover, skewness plays the same role in the equation for the evolution of turbulent kinetic energy (TKE). Hence when $w\u2032w\u2032\xaf$ (and therefore skewness) is positive, both $w\u2032w\u2032\xaf$ and TKE are being transported upwards. Similar considerations can be made in the case of horizontal skewness,

*s*.

_{u}Figures 11(a) and 11(b) show the vertical skewness profiles for probe ADCP1 and ADCP2, respectively. In the case of axial measurements (ADCP1), *s _{w}* is almost null up to $z/zss\u22481.1$. Then it decreases upward (with a minimum at $z\u22481.4\u22121.5$), and finally it increases up to null or also positive values (at $z/zss\u22481.8$). We infer that in the upper part of the fountain ($z/zss>1.1$) the transport phenomena are mainly downwards due to the formation and action of the counterflow. The experiments that exhibit a positive value are those with the higher densities, for which the fluctuations seem to present an upward transport ability. In the case of non-axial measurements (ADCP2), the skewness is slightly positive in the range $0.4<z/zss<1.1$, which means that both $w\u2032w\u2032\xaf$ and TKE are being transported upwards. For $z/zss>1.1$,

*s*is negative again (indicating transport in the downward direction).

_{v}Figure 8(c) shows the map of the horizontal velocity skewness for test 16. In the inner area of the fountain, the scenario is quite varied, and no particular conclusion can be drawn. On the contrary, on the sides of the jet (in the areas enclosed in the dotted rectangles) a clear tendency to the outward transport can be observed, as a consequence of the progressive widening of the jet. Notice that the same information can be extracted from the maps referring to other tests.

### D. Reynolds stresses and quadrant analysis

Substituting the mean and fluctuating components of the velocity in the momentum equation yields the turbulent stress components, which arise from the fluctuations. For expts. 15–20, the relative position of ADCP1, ADCP2, and ADCP3 allowed us to find the overlapping measurement volume between ADCP1-ADCP3 and between ADCP2-ADCP3, and to calculate the fluctuating velocity correlations $\u2212u\u2032w\u2032\xaf$, which represent the Reynolds shear stress at the net of the fluid density (note that correlation and stress have opposite sign).

Figure 12 shows the calculated non-dimensional Reynolds stresses. Even with the significant data dispersion, the different development of the axial and non-axial terms is clear. At the nozzle axis (ADCP1-ADCP3, panel a), the shear stress profile seems to be negative on average; we also notice that the peaks are observed for higher densities and around the mean rise height, where mixing conditions are enhanced. On the contrary, the non-axial profile (ADCP1-ADCP3, panel b) is always positive and presents the largest values at $z/zss\u22480.7$, well below the mean rise height; this is particularly true for test with lower densities, for which the plume widens not far from the nozzle.

To give a more detailed description of the turbulence structure, Reynolds shear stresses contributions are categorized according to their origin and divided into four quadrants.^{31} Then, conditionally sampling according to the quadrant gives the statistics of the events, as shown in Fig. 13.

The average shear stress for the *i*th quadrant is

where *N* is the total number of events and *j* is the current sample number. The total shear stress is

In our experiments, we define quadrant 1 (Q1) for $u\u2032>0$ and $w\u2032>0$, quadrant 2 (Q2) for $u\u2032<0$ and $w\u2032>0$, quadrant 3 (Q3) for $u\u2032<0$ and $w\u2032<0$, and quadrant 4 (Q4) for $u\u2032>0$ and $w\u2032<0$.

Figure 14 shows the Reynolds shear stress from each quadrant using ADCP1-ADCP3 velocity correlations (in axis measurements). Q3 and Q4 show that the highest relative contributions with maximum values are reached around the rise height, indicating an out-downward transport of momentum (both to the right and to the left). On the other hand, Q1 and Q2 have similar profiles with values that are maximum near the nozzle and decrease as $z/zss$ increases; we infer that in the region above the nozzle the transport tends to be out-upward.

The average shear stresses for ADCP1-ADCP3 velocity correlations are reported in Fig. 15. In this case (and especially for low density tests), the larger shear stresses are observed in Q1 and Q3, with $u\u2032w\u2032\xaf1$ slightly larger than $u\u2032w\u2032\xaf3$ and maximum values at $z/zss=0.7\u22120.9$. This is the area where the fountain widens and the momentum transport is mainly out-upward. The vertical profiles of $u\u2032w\u2032\xaf2$ and $u\u2032w\u2032\xaf4$ are nearly constant, but with larger values in Q4.

### E. Density measurements

For expts. 6 and 7, a conductivity probe was used to measure the temporal evolution of the density vertical profile. Similarly to the velocity components, we split the density as $\rho =\rho \xaf+\rho \u2032$, where $\rho \xaf$ is the mean (time-averaged) signal and $\rho \u2032$ is the fluctuating component obtained by subtracting the mean part to the whole signal (see Sec. II C).

A map showing the vertical density evolution is reported in Fig. 16. Results suggest that the vertical profile is subjected to fluctuations, but it is relatively stable over time (stationary) thanks to the continuous flux of source fluid, which is injected and mixes with the surrounding ambient fluid.

Figure 17(a) reports the mean and fluctuating densities, $\rho \xaf$ and $\rho rms\u2032$, respectively. The density profile shows a maximum near the inflow and linearly decreases with the distance from the nozzle. The RMS value is nearly constant in the vertical and shows a maximum at the rise height elevation ($z/zss\u22481$), where flow starts to reverse its direction and most of the mixing takes place. For a better comprehension of the results, we report the vertical velocity profile for both the mean and fluctuating components $w\xaf$ and $wRMS\u2032$, respectively [Fig. 17(b)]. Comparing the density and velocity profile, we see that the mean components have a similar trend (nearly linear decrease with distance from the source inflow). On the other hand, the density RMS shows a maximum at the rise height, while the velocity RMS is slightly decreasing upwards (this is especially true for low density test). We also calculate the non-dimensional correlation quantities $\rho \xafw\xaf$ and $\rho rms\u2032wrms\u2032$ [Fig. 17(c)], which are related to both the buoyancy and momentum fluxes, i.e., to the stabilizing and acting forces of the turbulent fountains, respectively. Combining velocity, density, and salinity, we can retrieve crucial parameters in order to determine mixing condition, e.g., the total buoyancy flux at the midsection of the source inflow.^{32} The overall results suggest that the fluctuating correlations have higher values at the mean rise height, enhancing density fluxes and mixing.

## IV. CONCLUSIONS

The widespread occurrence of fountains, along with the scarcity of data concerning the turbulent structure of the flow field, makes novel laboratory investigations a key element for further advances in the subject.

In this framework, experiments on forced fountains have been carried out in a controlled environment at the University of Parma (Italy). Present activity includes (i) measurements of the vertical and horizontal velocities in different positions and (ii) density profiling for some of the tests. The analysis regards the mean and turbulent characteristics of the flow, and it includes details about the turbulent kinetic energy, the velocity skewness, the Reynolds stress, and the correlation between density and velocity. In order to have an idea of the variability in the radial direction, data have been collected both along axial and non-axial vertical profiles.

The acoustic Doppler current profiler (ADCP) has been demonstrated to provide reliable measurements of the rise height, with the advantage of reducing the amount of data to be processed (especially when compared to image analysis techniques). On the other hand, the ADCP returns an average on a footprint whose diameter depends on the size of the probe itself, making the measurement non-punctual. This also affects the axial measurements, which we have only studied for $z/zss>0.4$. The comparison with traditional techniques is good, both as regards fluctuations of the interface and spectral analysis.

The vertical profiles of the mean velocity collapse fairly well on a single curve (that is a straight line) for ADCP1, with null velocity at $z/zss\u22481.2$. In the case of ADCP2, it is possible to observe the effects of the counterflow on the average velocity profiles, with a behavior that is strictly related to the density of the injected fluid. For lower densities, the fountain has a larger horizontal spreading that is detected by the probe, thus offering an indirect measurement of the boundary between the upflow and the counterflow.

The turbulent kinetic energy shows a maximum at $z/zss\u22481$. This suggests that the turbulence is mainly developed at the mean rise height, where the mixing between the flow from the nozzle and the flow generated by the periodic collapses of the plume takes place. Then, the turbulent kinetic energy decreases upward and becomes almost null at $z/zss\u22481.4$. The transport of $w\u2032w\u2032$ and TKE by the turbulence is mainly downward for $z/zss>1.1$, except for tests with the higher densities. Moreover, on the sides of the jet (for $0.2<z/zss<0.75$) a clear tendency to the horizontal outward transport can be observed.

The Reynolds shear stress profiles are quite disperse, but a clear spatial variability can be observed: (i) at the nozzle axis, the maximum values tend to gather at the mean rise height especially for the higher densities; (ii) at non-axial position, the vertical profile is negative with peaks at $z/zss\u22480.7$. A quadrant analysis was performed to highlight the main contributors to the stresses and their transport directions. The most relevant results are that (i) at the nozzle axis, higher shear stresses are observed in Q3 and Q4, corresponding to an out-downward transport of momentum (both to the right and to the left), with a peak at the mean rise height for higher densities; (ii) at a non-axial position, the higher shear stresses are observed in Q1 and Q3, with peaks at $z/zss\u22480.7\u22120.9$ for lower densities, in the area where the fountain widens and the momentum transport is mainly out-upward.

For expts. 6 and 7, the density profile presents a nearly linear trend, decreasing with distance from the nozzle. The correlation $\rho \u2032rmsv\u2032rms$, as well as the density RMS, shows a maximum at the rise height, indicating high levels of the buoyancy and momentum fluxes.

In summary, the present work aims to give a contribution for a better understanding of forced turbulent fountains, providing novel laboratory data and measurement techniques for both velocity and density.

## ACKNOWLEDGMENTS

The authors thank Professor Sandro Longo for the stimulating discussions and the valuable suggestions and Diana Petrolo for providing useful scripts for data analysis.

## DATA AVAILABILITY

The data that support the findings of this study are available from the corresponding author upon reasonable request.

## Nomenclature

- ADCP
acoustic Doppler current profiler

- $b\u0303ss$
fountain half-width

- expts
experiments

*Fr*_{0}source Reynolds number

- $g0\u2032$
buoyancy

*Q*discharge

*r*_{0}radial scale (internal nozzle radius)

*Re*_{T}threshold Reynolds number

*Re*_{0}source Reynolds number

- RMS
root mean square

*s*_{u}skewness of horizontal velocity

*s*_{w}skewness of vertical velocity

*St*Strouhal number

*t*time

*TKE*Turbulent kinetic energy (two components)

*u*horizontal velocity

- $u\u2032w\u2032\xaf$
Reynolds stress

*w*vertical velocity

*w*_{0}velocity scale (source velocity)

*x*horizontal axis

*z*vertical axis

- $zpe\xaf$
mean rise height peak

*z*_{ss}quasi-steady rise height

- $ztr\xaf$
mean rise height trough

- $\delta zss$
magnitude of the vertical fluctuations

*ρ*_{0}source density