Velocity and density measurements in forced fountains with negative buoyancy

In ﬂuid mechanics, fountains take place when a source ﬂuid is driven by its own momentum into a surrounding ambient ﬂuid, 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 ﬂow: we report the analysis of the turbulent kinetic energy, the velocity skewness, and the Reynolds stresses, including a quadrant analysis of the ﬂuctuating velocities. For some tests, the correlation between density and vertical velocity is investigated for both mean and ﬂuctuating values. We have quantiﬁed 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 ﬂuxes are present. The ability of acoustic Doppler current proﬁlers to identify the rise height of the fountain and to measure the velocity ﬁeld 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,3Fountains 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. 5Fountains 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. 6ther 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,9Hunt 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, Fr 0 ¼ w 0 = ffiffiffiffiffiffiffiffi g 0 0 r 0 p , where w 0 is the velocity at which fluid is ejected from the source, r 0 is the radial scale for the source, and g 0 0 is the buoyancy of the source fluid defined as g 0 0 ¼ gðq 0 À q a Þ=q a , where q 0 and q 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.
Furthermore, fountains may be regarded as laminar for source Reynolds numbers Re 0 Շ 120, turbulent for Re 0 տ 2000, and transitional for 120 Շ Re 0 Շ 2000, where Re 0 ¼ w 0 r 0 = 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 T at lower Fr 0 ), and they proposed Re T ¼ 75 Fr 0 þ 350 for Fr 0 > 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. 6fterward, 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,6Several 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. 169][20] Self-similarity is widely encountered in the study of buoyancy-and gravity-driven phenomena, including gravity currents and non-Newtonian flows. 21st 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. 13More 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. 23However, 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 Â 400 Â 800 mm 3 , as shown in Fig. 2. A vertical rigid tube with internal diameter D int ¼ 7:8 mm 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 6320 mm s À1 , 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 < N t , where the threshold value was taken as N t ¼ N e =3 and N e 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%. 24For 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 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 %20 Hz.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 %6 mm s À1 , 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 r 0 ¼ D int =2 ¼ 3:9 mm, the source fluid density was in the range 1020 q 0 1092 kg m À3 , and a value of source discharge Q ¼ 15 ml s À1 was used, which yields the source velocity w 0 ¼ Q=ðpr  TABLE II.Parameters of the experiments.Video indicates whether the video camera was used (active) or not (no).ADCP indicates (i) the number of acoustic Doppler current profilers deployed (ADCP1, ADCP2, or ADCP3), (ii) if the probes were in a fixed position (fixed) or they moved up and down (moving), and (iii) if the conductivity probe (Cond) was present.The variables r 0 , Q, w 0 , Re 0 , Fr 0 , and q 0 are the internal radius, the source fluid discharge, velocity, Reynolds number, Froude number, and density at the source section, respectively.The parameter Re T represents the Reynolds number threshold at which the fountain rise height is independent of Re (see Burridge, Mistry, and Hunt 12 ).

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 2 bss , the mean rise height peak z pe , and the mean rise height trough z tr .The analytical values of the defined quantities are as follows: where T is the acquisition time, z f ðtÞ the instantaneous value of the fountain height, bðz ss ; tÞ ¼ bðz ss ; tÞ when z f ðtÞ !z ss and bðz ss ; tÞ is the fountain half-width at z ¼ z ss .Furthermore: (i) z pe ðtÞ ¼ z f ðtÞ when z f ðtÞ !z ss þ r ss , (ii) z tr ðtÞ ¼ z f ðtÞ when z f ðtÞ z ss À r ss , and (iii) dt h ; dt pe , and dt tr are the total periods for which z f ðtÞ !z ss ; z f ðtÞ !z ss þ r ss , and z f ðtÞ z ss À r ss , respectively.The term r ss represents the standard deviation of the vertical fluctuation z ss over the acquisition time T. In addition, we define the magnitude of the vertical fluctuations as dz ss ¼ z pe À z tr .
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, z ss =r 0 , and the magnitude of the vertical fluctuations, dz ss (scaled both with the width of the forced fountains, dz ss =2 bss , and with the quasi-steady rise height, dz ss =z ss ), are well aligned.As found in BH2012, it suggests that (i) the forced fountains scale as z ss / r 0 Fr 0 (z ss ¼ 2:22 r 0 Fr 0 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 L , as St H ¼ f H T 0 and St L ¼ f L T 0 , respectively, where T 0 denotes an adequate timescale.A conventional timescale is represented by r 0 =w 0 , which results in a Strouhal number St / Fr À2 0 .A timescale proposed for forced fountains was w 0 =g 0 , for which the corresponding forced Strouhal number, St for , appears independent of Fr 0 .Another relevant timescale comes from the largeeddies length and velocity scales, i.e., / 2 bss =w ss , where w ss 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, St top .
Figure 5 shows the comparison between the values of the nondimensional found in BH2012 and in the present work.The results show a good overlap of the conventional Strouhal numbers St H and St L , as well as for the forced and top Strouhal numbers St for and St top (related to f H ).
In our experiments, we consider the quasi-steady rise height z ss as the vertical length scale and the fountain width bss for the horizontal length scale, while the vertical and horizontal velocities are nondimensional 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,28Here, 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  and extracted by the ADCP echo signal (red lines); panels (c) and (d) report the vertical velocity signal w(t); panels (e) and (f) report velocity power density spectra, S (vertical dashed lines represent the lower f L and higher f H peak frequencies).Upper panels refer to experiment 1 (q 0 ¼ 1030 kg m À3 ), lower panels refer to experiment 6 (q 0 ¼ 1089 kg m À3 ).

Physics of Fluids
interface (video data) and the velocity measured at the mean rise height (ADCP).We also see that fluctuations of w ADCP decrease with decreasing Fr 0 .Finally, the spectral analysis of signals is reported in panels (e) and (f) of Fig. 6.The dominant low frequency, f L , and the dominant high frequency, f H , of the fluctuations signal are calculated according to BH2013, and their value is compared with the power spectral density (PSD) of w video and w ADCP .The results qualitatively show that the dominant frequencies well represent the peak frequencies for w v as expected and also capture the main peaks in the spectra of the velocity w ADCP .
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 mm 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ðaÞ þv 0 ða; tÞ, where v(t, a) is the instantaneous velocity along the measuring direction a, while vðaÞ and v 0 ð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 %2 , which means that at the farther limit of the range (10 cm) the footprint radius is %3:5 mm 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, 2b u ðzÞ, equals the footprint diameter, d ADCP ðzÞ.Where d ADCP ðzÞ > 2b u ð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=z ss % 0: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=z ss % 1:2, and with a linear trend down to z=z ss % 0: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=z ss % 1: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=b ss % 65 (herein, the symbol ~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 0 ¼ v À v.The quantity v 0 rms 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 Շ z=z ss Շ 1; w 0 rms 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=z ss % 1: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 տ 1:1, with more intense fluctuations for higher densities.Figure 8(b) shows the map of u 0 rms , with a magnitude of the order of w 0 rms 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=z ss % 1.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=z ss % 1: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 03 represents the transport of w 0 w 0 by the turbulence itself.Moreover, skewness plays the same role in the equation for the evolution of turbulent kinetic energy (TKE).Hence when w 0 w 0 (and therefore skewness) is positive, both w 0 w 0 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=z ss % 1:1.Then it decreases upward (with a minimum at z % 1:4 À 1:5), and finally it increases up to null or also positive values (at z=z ss % 1:8).We infer that in the upper part of the fountain (z=z ss > 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=z ss < 1:1, which means that both w 0 w 0 and TKE are being transported upwards.For z=z ss > 1:1, s v is negative again (indicating transport in the downward direction).
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 Àu 0 w 0 , 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=z ss % 0: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. 31Then, conditionally sampling according to the quadrant gives the statistics of the events, as shown in Fig. 13.
The average shear stress for the ith 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 0 > 0 and w 0 > 0, quadrant 2 (Q2) for u 0 < 0 and w 0 > 0, quadrant 3 (Q3) for u 0 < 0 and w 0 < 0, and quadrant 4 (Q4) for u 0 > 0 and w 0 < 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 outdownward 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=z ss increases; we infer that in the region above the nozzle the transport tends to be outupward.
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 0 w 0 1 slightly larger than u 0 w 0 3 and maximum values at z=z ss ¼ 0:7 À 0:9.This is the area where the fountain widens and the momentum transport is mainly out-upward.The vertical profiles of u 0 w 0 2 and u 0 w 0 4 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 q ¼ q þ q 0 , where q is the mean (time-averaged) signal and q 0 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, q and q 0 rms , 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=z ss % 1), 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 and w 0 RMS , 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 q w and q 0 rms w 0 rms [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. 32The 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=z ss > 0:4.The comparison with traditional techniques is good, both as regards fluctuations of the interface and spectral analysis.

Physics of Fluids
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=z ss % 1: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=z ss % 1.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=z ss % 1:4.The transport of w 0 w 0 and TKE by the turbulence is mainly downward for z=z ss > 1:1, except for tests with the higher densities.Moreover, on the sides of the jet (for 0:2 < z=z ss < 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=z ss % 0: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 outdownward 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=z ss % 0:7 À 0: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 q 0 rms v 0 rms , 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.
u 0 w 0 Reynolds stress w vertical velocity w 0 velocity scale (source velocity) x horizontal axis z vertical axis z pe mean rise height peak z ss quasi-steady rise height z tr mean rise height trough dz ss magnitude of the vertical fluctuations q 0 source density

0: 3 Շ
FIG. 1. Scheme of a forced turbulent fountain: [(a) and (b)] at the initial stage and [(c) and (d)] during the subsequent steady state.
2 0 Þ ¼ 314 mm s À1 .The ranges of the nondimensional groups are 5:4 6 0:4 Fr 0 16:4 6 1:4 and 1030 6 70 Re 0 1180 6 80 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.

FIG. 4 .
FIG. 4. Rise height statistics in the present experiments (red filled diamonds) and in BH2012 (gray filled crosses): (a) non-dimensional quasi-steady rise height z ss =r 0 , symbols are experimental points and solid line represents the best-fitting relation by BH2012 z ss ¼ 2:46 r 0 Fr 0 (z ss ¼ 2:22 r 0 Fr 0 considering only the present experiments); (b) vertical fluctuations scaled with the fountain width dz ss =2 bss ; and (c) vertical fluctuations scaled with the quasi-steady rise height dz ss =z ss ).The vertical lines separate the fountains in VWWI (very weak-weak-intermediate), forced and highly forced regimes, following the classification by BH2012.

5 VFIG. 6 .
FIG. 5. Strouhal number (representing non-dimensional fluctuations peak frequencies) as a function of Fr 0 .Diamonds are experiments of the present work, crosses are experiments from BH2013.(a) Conventional Strouhal number St ¼ fr 0 =w 0 ; red diamonds refer to higher peak frequency f H , blue diamonds refer to lower peak frequency f L ; (b) forced Strouhal number St for ¼ f H w 0 =g 0 ; (c) top Strouhal number St top ¼ f H 4 bss =w ss ; the horizontal line indicate St top ¼ 1.The vertical line separates the fountains in intermediate and forced regimes, following the classification by BH2013.

FIG. 10 .
FIG. 10.(a) Two-component turbulent kinetic energy (TKE) obtained using measurements from ADCP1 and ADCP3.(b) Comparison of RMS vertical fluctuating velocities w 0 rms =w 0 at the nozzle axis between present (colored symbols) and literature experiments (black filled circles).

FIG. 14 .
FIG. 14.Average shear stress quadrant decomposed for ADCP1-ADCP3 measurements (axial position); values are non-dimensional with respect to w 2 0 .See the caption of Fig. 7 for details about symbols.

FIG. 15 .
FIG. 15.Average shear stress quadrant decomposed for ADCP2-ADCP3 measurements (non-axial position); values are non-dimensional with respect to w 2 0 .See the caption of Fig. 7 for details about symbols.

TABLE I .
Classification of fountains according to the source Froude number.

TABLE III .
Detection of the rise height of the turbulent fountain: comparison between video (subscript v) and acoustic Doppler current profilers (subscript D).Expt.indicates the number of the experiment from TableII, q 0 is the density at the source section, z ss and r zss are the mean and the standard deviation of the rise height, respectively.