Unconfined Dense Plunging Jets Used for Brine Disposal from Desalination Plants

Laboratory experiments were conducted to measure entrained air bubble penetration depth and dilution of a dense vertical unconfined plunging jet to evaluate its performance as an outfall to dilute brine from desalination plants as well as a means to aerate water column. Experiments involved neutrally buoyant or dense plunging jets discharging in quiescent receiving water. The density difference between effluent and receiving water, the plunging jet length (height above water surface), and the receiving water salinity were varied in the experiments. Observed penetration depth for neutrally buoyant jets was somewhat greater than previously reported, and increased modestly with jet density. Increasing density also resulted in an increasing number of fine bubbles descending together with the dense plume. These observations can help guide the design of plunging jets to mitigate anoxic conditions in the water column when brine is introduced to a receiving water body, as with seawater desalination.


Introduction
This study is motivated by the potential use of plunging jets to dispose of reject desalination brine into Arabian Gulf (hereafter referred to as the Gulf) waters, in an effort to mitigate environmental impacts caused by brine discharges from Kuwait's desalination plants. Brine discharge from desalination plants using reverse osmosis (RO) and multistage flash (MSF) distillation technologies may produce depressions in the dissolved oxygen (DO) level, and increased salinity and sea surface temperature (SST) [1,2]. In turn, these could potentially cause degradation of Kuwait's marine diversity [3,4] both at the plant intakes and discharge outfalls [5] and further offshore [6,7].
With the Gulf Cooperation Council (GCC) countries' desalination capacity estimated at 20 million m 3 /day of freshwater in 2012 [8], and with an estimated 2 m 3 brine discharged per m 3 of produced freshwater [9], the total volume of brine discharged into the Gulf from the GCC alone could be as high as 40 million m 3 /day. In the Gulf, 88% of the desalination plants operated by GCC countries use MSF and <10% use RO. Most of these desalination plants are located on the coast, and discharge brine into the near shore after mixing it partially with ambient seawater [3,10]. In Kuwait, contaminants from RO plants are discharged to the storm drains that eventually go to the Arabian Gulf; while from MSF plants, brine is discharged directly to the Arabian Gulf [11][12][13][14][15]. Very shallow coastal bathymetry, increased brine discharge, large distance from the limited flushing from the Indian Ocean, and reduced fresh water discharge from the nearby Shatt Al-Arab have made Kuwait's marine waters the densest in the Gulf [16].
Reject brine discharges pose a risk of creating anoxic zones near discharge locations [7]. Therefore, it is worth exploring alternative engineering solutions that could combine the dilution of reject brine with reaeration from the surface. Figure 1 considers a pedagogic progression of four options for brine discharge: a) a bottom staged diffuser where dilution is achieved by jet momentum; b) a sinking dense plume caused by discharging brine from the water surface; c) unconfined plunging jets caused by discharging brine above the water surface; and d) confined plunging jets in which the release is shrouded by a downcomer. Figure 1a shows a staged diffuser scheme that might be used to dilute thermal or wastewater effluent discharge in the receiving water body. The dilution (S) resulting from a staged diffuser of length L is well studied and is dependent on the ambient current u a , as well as the nozzle velocity u 0 , as follows [17][18][19]: where H is the receiving water depth and Q 0 is the effluent discharge flowrate. Indian Ocean, and reduced fresh water discharge from the nearby Shatt Al-Arab have made Kuwait's marine waters the densest in the Gulf [16]. Reject brine discharges pose a risk of creating anoxic zones near discharge locations [7]. Therefore, it is worth exploring alternative engineering solutions that could combine the dilution of reject brine with reaeration from the surface. Figure 1 considers a pedagogic progression of four options for brine discharge: a) a bottom staged diffuser where dilution is achieved by jet momentum; b) a sinking dense plume caused by discharging brine from the water surface; c) unconfined plunging jets caused by discharging brine above the water surface; and d) confined plunging jets in which the release is shrouded by a downcomer. Figure 1a shows a staged diffuser scheme that might be used to dilute thermal or wastewater effluent discharge in the receiving water body. The dilution ( ) resulting from a staged diffuser of length is well studied and is dependent on the ambient current , as well as the nozzle velocity , as follows [17][18][19]: where is the receiving water depth and is the effluent discharge flowrate. Picking a worst case ambient velocity of zero and using typical values for a desalination plant discharge of = 1 m 3 /s (~20 million gallon per day), = 4 m/s, and = 3 m, the required diffuser length to achieve a target dilution of 20 (e.g., to reduce the excess salinity from ~40 psu to ~2 psu above ambient) is ≈ 175 m (Equation (1)). Figure 1b shows a schematic of a line source of buoyancy released on the surface of a receiving water body, e.g., via a floating pipe. The discharge can be thought of as an upside down staged diffuser where the discharge is located on the water surface, and there is no horizontal momentum. For a weak ambient current (normalized ambient velocity ⁄ ≤ 0.1, where = Δ / , is the buoyancy flux per unit length), and without ambient stratification, the dilution for a buoyant line source is [20,21]: Picking a worst case ambient velocity of zero and using typical values for a desalination plant discharge of Q 0 = 1 m 3 /s (~20 million gallon per day), u 0 = 4 m/s, and H = 3 m, the required diffuser length to achieve a target dilution of 20 (e.g., to reduce the excess salinity from~40 psu to~2 psu above ambient) is L ≈ 175 m (Equation (1)). Figure 1b shows a schematic of a line source of buoyancy released on the surface of a receiving water body, e.g., via a floating pipe. The discharge can be thought of as an upside down staged diffuser where the discharge is located on the water surface, and there is no horizontal momentum. For a weak ambient current (normalized ambient velocity u 3 a /b 0 ≤ 0.1, where b 0 = Q 0 ∆ρ 0 g/ρ w L, is the buoyancy flux per unit length), and without ambient stratification, the dilution for a buoyant line source is [20,21]: where q 0 = Q 0 /L is the initial brine flowrate per unit length, ∆ρ 0 is the effluent density difference, ρ w is the ambient water density, g is the gravitational acceleration, H is the ambient water depth, and c 1 is a proportionality constant reported as 0.27 by [20] for a finite length line source. To achieve a target near field dilution of 20, Equation (2) was used to predict the required diffuser length as L ≈ 225 m. Figure 1c depicts a line of dense plunging jets, similar to the outfall in Figure 1b, except that the ports are located at a fixed elevation above the surface, e.g., connected to a manifold supported by a pier. The plunging jet reactor concept has been in use for several decades, as a means of achieving high mass transfer rates by entraining gas bubbles into a liquid, at low capital and operating costs [22][23][24]. The typical objective is to contact the two phases to promote mass transfer. Plunging jet reactors are widely employed in a variety of processes, such as gas-liquid reactors in the chemical industry, aerobic wastewater treatment, air pollution abatement, froth flotation, and fermentation [25]. Figure 1d depicts an array of confined plunging liquid jet reactors (CPLJR), consisting of vertical confining tubes that are partially immersed in the liquid pool. The CPLJR is utilized to try to improve the mixing and mass transfer rates from the ambient air into the liquid. This can be achieved by: (1) increasing the jet penetration depth and (2) confining the bubbly jet region to increase the contact time between the gas and liquid [22,23]. Table 1 (modified from [26,27]) shows the oxygen transfer efficiencies (defined as the mass of oxygen introduced to a reactor system per kWh of energy expenditure) reported for a range of aeration technologies. As seen from Table 1, plunging jet reactors (both confined and unconfined) are capable of achieving oxygen efficiencies comparable to other conventional aeration technologies. Because plunging jets were historically focused on enhancing mass transfer from the gas phase to a single liquid phase, the vast majority of previous studies are performed with neutrally buoyant jets, i.e., the jet and the receiving water are of the same density [31][32][33]. Some statistical studies have been performed for plunging jets of neutral buoyancy [34], and numerical and experimental work has been performed for a plunging jet of a smaller density than the ambient [35], but there is a lack of literature available for plunging jets with higher density than the receiving water.
There is also a lack of available measurements of the dilution capacity created by a plunging jet when discharged either unconfined or confined into a receiving water body. The main objectives of this paper, therefore, are to assess the effect of density difference (created by brine jets of varying salinity) on the behavior of a plunging jet, as well as measure the dilution capability of a dense plunging jet to see if it is comparable to the dilution capability of submerged or surface discharges. This study offers a schematic model for the calculation of the increase in dissolved oxygen concentration in the water column as a result of the use of an unconfined plunging jet. These estimates may inform the design of plunging jet reactors for desalination brine discharge that are optimized to attain maximum brine dilution and DO dissolution from the atmosphere. This study focuses on unconfined plunging jet reactors (Figure 1c), but future effort will include confined plunging jet reactors.

Plunging Jet Behavior
The above discussion treated a row of buoyant jets as a line source. Of course, a diffuser is comprised of individual jets, so the analysis starts with a single (unconfined) liquid jet. As shown in Figure 2, the jet plunges into an open liquid pool, creating a conical downflow dispersion of fine bubbles and a surrounding upflow of larger coalesced bubbles. Depending on the density difference between the effluent and receiving water, the fine bubbles may or may not reach the bottom. For cases in which the fine bubbles reach the bottom, they spread radially outwards with the diluted effluent flow, and then rise, forming an annulus wider than the annulus of large bubbles. The observed behavior of neutrally buoyant and dense plunging jets is described below.
Processes 2020, 8, x FOR PEER REVIEW 4 of 18 this paper, therefore, are to assess the effect of density difference (created by brine jets of varying salinity) on the behavior of a plunging jet, as well as measure the dilution capability of a dense plunging jet to see if it is comparable to the dilution capability of submerged or surface discharges. This study offers a schematic model for the calculation of the increase in dissolved oxygen concentration in the water column as a result of the use of an unconfined plunging jet. These estimates may inform the design of plunging jet reactors for desalination brine discharge that are optimized to attain maximum brine dilution and DO dissolution from the atmosphere. This study focuses on unconfined plunging jet reactors (Figure 1c), but future effort will include confined plunging jet reactors.

Plunging Jet Behavior
The above discussion treated a row of buoyant jets as a line source. Of course, a diffuser is comprised of individual jets, so the analysis starts with a single (unconfined) liquid jet. As shown in Figure 2, the jet plunges into an open liquid pool, creating a conical downflow dispersion of fine bubbles and a surrounding upflow of larger coalesced bubbles. Depending on the density difference between the effluent and receiving water, the fine bubbles may or may not reach the bottom. For cases in which the fine bubbles reach the bottom, they spread radially outwards with the diluted effluent flow, and then rise, forming an annulus wider than the annulus of large bubbles. The observed behavior of neutrally buoyant and dense plunging jets is described below. When the vertical jet issues from the nozzle, the jet descends through the headspace over a vertical jet length . There is a minimum jet descent velocity, known as the onset velocity, above which the headspace air begins to be entrained by the water jet. The onset velocity is dependent on the water density, water dynamic viscosity, surface tension between the air and water phase, turbulent velocity fluctuation, turbulent length scale, and the angle of the jet [25,36]. So far, there are no consistent predictions for the onset velocity; however, for vertical water jets, the onset velocities reported are in the range of 1.5-2.5 m/s [33,37,38].
Al-Anzi et al. [22] describes that increasing the jet length increases the amplitude of the disturbances on the surface of a rough jet, because these disturbances have a longer time to grow since their inception at the nozzle exit. The increased roughness, in turn, increases the rate of air entrainment into the jet prior to the jet hitting the water surface.
After breaking the surface of water, the impact momentum of the jet is responsible for carrying the two-phase jet/air mixture downward. For jets that have entrained air bubbles upon plunging into the receiving water body, the air bubbles descend a certain maximum distance (the penetration depth, denoted ) before their buoyancy drives them back to the water surface (Figure 2a,b). Numerous When the vertical jet issues from the nozzle, the jet descends through the headspace over a vertical jet length L j . There is a minimum jet descent velocity, known as the onset velocity, above which the headspace air begins to be entrained by the water jet. The onset velocity is dependent on the water density, water dynamic viscosity, surface tension between the air and water phase, turbulent velocity fluctuation, turbulent length scale, and the angle of the jet [25,36]. So far, there are no consistent predictions for the onset velocity; however, for vertical water jets, the onset velocities reported are in the range of 1.5-2.5 m/s [33,37,38].
Al-Anzi et al. [22] describes that increasing the jet length increases the amplitude of the disturbances on the surface of a rough jet, because these disturbances have a longer time to grow since their inception at the nozzle exit. The increased roughness, in turn, increases the rate of air entrainment into the jet prior to the jet hitting the water surface.
After breaking the surface of water, the impact momentum of the jet is responsible for carrying the two-phase jet/air mixture downward. For jets that have entrained air bubbles upon plunging into the receiving water body, the air bubbles descend a certain maximum distance (the penetration depth, denoted H p ) before their buoyancy drives them back to the water surface (Figure 2a,b). Numerous authors have observed the bubble penetration depth of neutrally buoyant plunging jets in water [25,33,38,39]. Bin [25] fit the following empirical relations for H p (in m) as a function of the nozzle velocity V 0 (in m/s) and the nozzle diameter d 0 (in m): Processes 2020, 8, 696 Note that the quantity V 0 d 0 is related to the initial kinematic momentum of the jet as (3) and (4) are dimensional as they incorporate a constant value for the density difference between the air and water phases. In this section, the dependence of H p is reframed using the jet kinematic momentum upon impact,

The proportionality constants in Equations
; the density ratio of the jet and receiving water body ρ b ρ w ; and gravitational constant (g). Here, V j = V 2 0 + 2gL j (after [36]) and d j = V 0 d 0 2 /V j are the jet velocity and diameter at the water surface. For neutrally buoyant jets, ρ b ρ w = 1 and by dimensional analysis, the penetration depth of the bubbles in the receiving water should follow the relation: where α is a proportionality constant determined to be α = 4.0 in the section below. For dense jets, ρ b ρ w > 1, and the following is presumed: There exists a regime of dense brine plunging jets where the jet liquid buoyancy is sufficient to temporarily overcome the bubble's positive buoyancy, pushing the bubbles to the bottom of the receiving water body (Figure 2c). In this case, the penetration depth of the bubbles is greater than the receiving water depth (H). The bubbles, after reaching the bottom, travel radially outwards along with the dense jet liquid, and, once they are sufficiently far from the sinking plume core, they rise passively towards the water surface.
The plunging jet itself, whether or not it entrains air bubbles, is diluted as it mixes with the ambient receiving water body. Measurements of near field dilution of a plunging brine jet were compared with available literature (described below) for a dense jet released just below the water surface, or the inverted situation of a positively buoyant jet released at the bottom.
Chanson et al. [36] expressed the dilution at elevation z of a round momentum jet discharging vertically upward into an unstratified ambient using the expression: where l Q = Q 0 /M 0.5 0 is the volume length scale. For a buoyancy dominated discharge, Tian et al. [21] developed a prediction for the dilution of a buoyant freshwater or low salinity wastewater discharge released just above the bottom of the sea floor into a saline ambient. The dilution at depth H is given by: where B 0 = ∆ρ 0 ρ w gQ 0 is the initial buoyancy flux. Fischer et al. [40], Section 9.

Materials and Methods
The experimental setup for observing brine plunging jets is presented here. A water tank, at the Parsons Laboratory at MIT, of dimensions 4.9 m × 1.2 m × 0.6 m, filled with fresh or salt water to H = 0.5 m depth, acts as the receiving water body. A submersible pump located in a separate brine reservoir (~140 L capacity) is able to deliver a steady liquid flowrate of up to 0.17 L/s downwards perpendicular to the water surface to form the vertical plunging jet. The discharge is dyed with Rhodamine WT and its concentration is measured using a Turner Designs Cyclops 7F model fluorometer.
In order to determine the dilution and aeration capability of these dense plunging jets, the following parameters were varied: the plunging jet length L j (the distance of the nozzle above the water surface); the density of the brine made with sodium chloride (expressed as ∆s, the difference in salinity between the effluent and the receiving water); and the density of the receiving water. The main experimental parameters were as follows: ∆s = 0, 10, 40 psu and L j = 0, 0.15, 0.30, 0.60 m, resulting in a suite of 12 experimental runs with receiving water salinity (s w ) of zero. An additional set of 3 experiments with ∆s = 0, 10, 40 psu and L j = 0.30 m were conducted with s w = 40 psu. In all of the experiments, a straight cylindrical nozzle was used with an internal diameter of d 0 = 10 mm, and a liquid flowrate of Q 0 = 163 cm 3 /s (V 0 = 2.1 m/s). Table 2 shows the experimental parameters for all the experiments. The jet Froude number Fr j = V j 2 /gd j is also listed in Table 2.
The above experimental parameters can be scaled to prototype values using Froude scaling. Assuming a desalination plant with brine flowrate of 1 m 3 /s using 6 jets to discharge, a single plunging jet will have a flowrate of 0.167 m 3 /s. The volume ratio between the field prototype and lab is, thus, Q r = 0.167 m 3 /s 163 cm 3 /s = 1024. By Froude scaling, with similar density differences between model and prototype, this is equivalent to a length scale ratio of L r = Q 0.4 r = 16 between the field prototype and the lab. Thus, the receiving water depth of 50 cm and nozzle diameter of 10 mm in the lab scale to prototype water depth of 8 m and prototype nozzle diameter of 16 cm. The lab jet lengths of 0, 15, 30, and 60 cm scale to 0, 2.4, 4.8, and 9.6 m, respectively, in the field.
Measurements were made of: (1) the penetration depth of the plunging jet (H p , the maximum depth to which the entrained gas bubbles descend within the water column); and (2) dye concentration on the tank bottom at a radial distance from the source (R) of 0.5 m (i.e., equal to one water depth). For some experiments, additional dye concentration measurements were made at the impact point on the tank bottom (i.e., R = 0). The measured concentrations were used to calculate the dilution on the tank bottom at R = 50 cm (S R ) and dilution at the impact point on the tank bottom (S H ). S R was not measured for experiments with ∆s = 0 because the plume took a long time to reach the measurement location. H p was not measured for experiments with L j = 60 cm.

Measurement of Bubble Penetration Depth
For the measurement of penetration depth and the width of the bubble plume, a camera was mounted to the side of the tank to image the plunging jet behavior below the water surface at 60 frames per second. Since the maximum depth of the bubbles fluctuates with time (with a period of about 1 s and amplitude of about 2 cm), an average of all of the observed maximum depths and maximum bubble widths over the duration of each experiment was computed, using a total of~5000 frames per experimental run. Figure 3 shows snapshots of movies taken of laboratory runs with L j = 30 cm. It can be observed from Figure 3 that small and large bubbles alike are plunged down to the penetration depth. For neutrally buoyant and moderately dense jets (∆s = 0 and 10 psu), the bubbles were observed to plunge to about the same depths (range of 0.18-0.19 m; Table 2), similar to the predicted value of 0.19 m using Equation (3) [25]. However, for brine with ∆s = 40 psu, while large bubbles behaved similarly to the neutrally buoyant plunging jet, smaller bubbles were pushed downwards to the bottom of the tank, spreading radially along the tank bottom before rising again up to the surface, resulting in larger H p and plume widths. The penetration depth and width of the bubble plume were determined using image processing as follows. Each image of an experimental run was split into 3 separate images, each comprised of the red, green, and blue pixel intensities. MATLAB was used to apply an image filter that uses the Sobel operator to create a new image that exaggerates the edges present in each of the 3 images. The pixel intensities for each row of pixels were scanned from left to right, in the portion of the image that the plume is located (i.e., ignoring any bottom features as well as the bubbles on the water surface). Prior to each experimental run, an image of a ruler placed in the plane of the bubble plume centerline was taken in order to determine the conversion of image pixels to centimeters without the effect of parallax. The widest extent of the image edges (predominantly using the green image) was taken as , the width of the bubble plume. The deepest extent of the image edge was taken as the bubble penetration depth . See Figure 4 for an example of bubble plume width determination. Experimental results for the penetration depth are shown in Table 2.   The penetration depth and width of the bubble plume were determined using image processing as follows. Each image of an experimental run was split into 3 separate images, each comprised of the red, green, and blue pixel intensities. MATLAB was used to apply an image filter that uses the Sobel operator to create a new image that exaggerates the edges present in each of the 3 images. The pixel intensities for each row of pixels were scanned from left to right, in the portion of the image that the plume is located (i.e., ignoring any bottom features as well as the bubbles on the water surface). Prior to each experimental run, an image of a ruler placed in the plane of the bubble plume centerline was taken in order to determine the conversion of image pixels to centimeters without the effect of parallax. The widest extent of the image edges (predominantly using the green image) was taken as W, the width of the bubble plume. The deepest extent of the image edge was taken as the bubble penetration depth H p . See Figure 4 for an example of bubble plume width determination. Experimental results for the penetration depth are shown in Table 2.  The penetration depth and width of the bubble plume were determined using image processing as follows. Each image of an experimental run was split into 3 separate images, each comprised of the red, green, and blue pixel intensities. MATLAB was used to apply an image filter that uses the Sobel operator to create a new image that exaggerates the edges present in each of the 3 images. The pixel intensities for each row of pixels were scanned from left to right, in the portion of the image that the plume is located (i.e., ignoring any bottom features as well as the bubbles on the water surface). Prior to each experimental run, an image of a ruler placed in the plane of the bubble plume centerline was taken in order to determine the conversion of image pixels to centimeters without the effect of parallax. The widest extent of the image edges (predominantly using the green image) was taken as , the width of the bubble plume. The deepest extent of the image edge was taken as the bubble penetration depth . See Figure 4 for an example of bubble plume width determination. Experimental results for the penetration depth are shown in Table 2.

Measurement of Dilution
For the measurement of dilution, dye concentration was measured using a fluorometer. Each experiment was run for a duration of 5-6 min with the fluorometer recording concentrations every second. Each experiment was also replicated (number of replicates = 2-5) and the reported dilution is the average of all the measurements after concentrations appeared to be quasi-steady. Measured dilutions are reported in Table 2. The variation in dilution indicated in Table 2 shows the standard deviation of all the concentration measurements. It can be seen from Table 2 that the measured dilutions are comparable to the predicted dilutions for jets discharged just below the water surface.  (5) with α = 4.0 b calculated using Equation (7). c calculated using Figure 8.7, Section 9.2.3 from [40].

Results
The measurements of H p in Table 2 are plotted against M j /g 1/3 in Figure 5, and compared with experiments from earlier literature (compiled by [25,33]) for the four experiments for which fine bubbles were not pushed to the bottom of the tank. The experimental data from this study (with ρ b ρ w > 1) is observed to be somewhat higher than the fitted curve for ρ b ρ w = 1. However, there were insufficient experiments with significant variation in the parameter ρ b /ρ w to determine the form of function ϕ in Equation (6). Due to large variation, and the relatively small value of M j /g 1/3 for the experiments, our curve fitting was limited to M j /g 1/3 < 0.125 m. From Figure 5, the combined results yield a proportionality constant of α = 4.0, which was used for predictions of H p using Equation (5). The presence of the bubbles in the bulk fluid serves to decrease the buoyancy of the descending plume. The bulk buoyancy of the descending buoyant jet is given by: where Δ is the density difference between the air and the ambient water, and is the entrained gas flowrate.
An estimate of the maximum effect that the bubbles may have on the bulk buoyancy is to estimate the bulk jet buoyancy when all of the brine and the entrained bubbles are part of a single phase. The entrained gas volume rate for a plunging jet is dependent on jet momentum, jet length, and perimeter (which, in turn, is dependent on the jet diameter). The empirical relation for an unconfined jet plunging vertically downward into an unconfined water body is [25,41] Substituting Equation (10) into Equation (9) yields: Using values from the current experiments (for which was measured), the magnitude of the second term in the bracket of Equation (11) is in the range of 6 to 49 times the first term, i.e., the entrained air buoyancy is an order of magnitude greater than the buoyancy of the descending plume. This would indicate that if the entrained bubbles and the plunging brine jet were to act as a single phase, the jet would not sink. However, as the observations show the jets to sink, it is likely that only a small fraction of the entrained bubbles stays in a single phase with the plunging brine. The relative magnitudes of the positive buoyancy of the bubbles and negative buoyancy of the plume can also be used to explain the observed penetration depths. For experiments with ∆ = 10 psu, the positive buoyancy due to the bubbles is much larger than the negative buoyancy of the plume. Thus, the The presence of the bubbles in the bulk fluid serves to decrease the buoyancy of the descending plume. The bulk buoyancy of the descending buoyant jet is given by: where ∆ρ g is the density difference between the air and the ambient water, and Q g is the entrained gas flowrate. An estimate of the maximum effect that the bubbles may have on the bulk buoyancy is to estimate the bulk jet buoyancy when all of the brine and the entrained bubbles are part of a single phase. The entrained gas volume rate for a plunging jet Q g is dependent on jet momentum, jet length, and perimeter (which, in turn, is dependent on the jet diameter). The empirical relation for an unconfined jet plunging vertically downward into an unconfined water body is [25,41]: Substituting Equation (10) into Equation (9) yields: Using values from the current experiments (for which H p was measured), the magnitude of the second term in the bracket of Equation (11) is in the range of 6 to 49 times the first term, i.e., the entrained air buoyancy is an order of magnitude greater than the buoyancy of the descending plume. This would indicate that if the entrained bubbles and the plunging brine jet were to act as a single phase, the jet would not sink. However, as the observations show the jets to sink, it is likely that only a small fraction of the entrained bubbles stays in a single phase with the plunging brine. The relative magnitudes of the positive buoyancy of the bubbles and negative buoyancy of the plume can also be used to explain the observed penetration depths. For experiments with ∆s = 10 psu, the positive buoyancy due to the bubbles is much larger than the negative buoyancy of the plume.
Thus, the bubbles are not carried to the bottom of the tank for such cases. However, for cases with ∆s = 40 psu, the negative buoyancy of the plume is large enough to carry the bubbles to the bottom.
As discussed in the previous section, the current experimental parameters can be scaled to prototype parameters using Froude scaling. But this scaling does not apply to bubble size which depends on other parameters such as surface tension. The bubble size in the field can be estimated using a Weber number We = ρ w V 2 j d j /σ scaling. (Here, σ is the surface tension and ρ w is the density of water). If the ratio of bubble size to nozzle diameter scales as We −0.6 [42], then the ratio of bubble diameter in the field to that in the lab is L r −0.2 (equal to 0.57 with L r = 16). Thus, the bubble diameter in the field is expected to be somewhat smaller than the bubble diameter observed in the lab. Figure 6 shows the dilutions measured at R = 50 cm (S R ) plotted against ∆s. These measurements are expected to be greater than S H,pred because the plume is allowed to spread and mix with the ambient water as a radial gravity current after impacting the bottom of the tank. S R increases slightly with increase in ∆s, and shows a non-monotonic dependence on L j . As seen in Figure 6, dilution reduces as L j is increased from 0 to 15 cm, but increases as L j is increased further (for both values of ∆s). As shown in Table 2, S H,pred increases with increase in L j because the jet momentum used in the prediction increases. However, this trend is not replicated in S R because the jets lose some of the momentum upon impact at the water surface. The effect of receiving water salinity (s w ) can also be seen in Figure 6, with no significant difference seen in dilution for runs with s w = 0 and 40 psu (L j = 30 cm). bubbles are not carried to the bottom of the tank for such cases. However, for cases with ∆ = 40 psu, the negative buoyancy of the plume is large enough to carry the bubbles to the bottom. As discussed in the previous section, the current experimental parameters can be scaled to prototype parameters using Froude scaling. But this scaling does not apply to bubble size which depends on other parameters such as surface tension. The bubble size in the field can be estimated using a Weber number = ⁄ scaling. (Here, is the surface tension and is the density of water). If the ratio of bubble size to nozzle diameter scales as . [42], then the ratio of bubble diameter in the field to that in the lab is . (equal to 0.57 with = 16). Thus, the bubble diameter in the field is expected to be somewhat smaller than the bubble diameter observed in the lab. Figure 6 shows the dilutions measured at = 50 cm ( ) plotted against ∆ . These measurements are expected to be greater than , because the plume is allowed to spread and mix with the ambient water as a radial gravity current after impacting the bottom of the tank. increases slightly with increase in ∆ , and shows a non-monotonic dependence on . As seen in Figure 6, dilution reduces as is increased from 0 to 15 cm, but increases as is increased further (for both values of ∆ ). As shown in Table 2, , increases with increase in because the jet momentum used in the prediction increases. However, this trend is not replicated in because the jets lose some of the momentum upon impact at the water surface. The effect of receiving water salinity ( ) can also be seen in Figure 6, with no significant difference seen in dilution for runs with = 0 and 40 psu ( = 30 cm). Dilutions measured at = 0 are plotted in Figure 7. As expected, the dilutions measured at = 0 are smaller than measurements at = 50 cm because they do not account for mixing after the plume hits the bottom (when it spreads radially). For = 0, the measured dilutions are lower than the predicted values at = for a buoyant jet. This is likely due to the presence of a layer of diluted brine at the bottom which limits the entrainment of ambient water. For = 60 cm, the measured dilution is higher than that of a pure plume (discharge with zero momentum), but lower than the theoretical value for a jet with = 60 cm. This can be explained by: 1) the presence of layer of diluted brine; and 2) loss of momentum of the plunging jet after it hits the water surface. Note that the dilution of a pure plume is independent of (as it has no momentum). Dilutions measured at R = 0 are plotted in Figure 7. As expected, the dilutions measured at R = 0 are smaller than measurements at R = 50 cm because they do not account for mixing after the plume hits the bottom (when it spreads radially). For L j = 0, the measured dilutions are lower than the predicted values at z = H for a buoyant jet. This is likely due to the presence of a layer of diluted brine at the bottom which limits the entrainment of ambient water. For L j = 60 cm, the measured dilution is higher than that of a pure plume (discharge with zero momentum), but lower than the theoretical value for a jet with L j = 60 cm. This can be explained by: 1) the presence of layer of diluted brine; and 2) loss of momentum of the plunging jet after it hits the water surface. Note that the dilution of a pure plume is independent of L j (as it has no momentum).

Aeration of Water Column
Al Anzi et al. [22] proposed the use of plunging liquid jet reactors, arranged in a diffuser array, to promote the aeration of the underlying water via entrainment of air bubbles from above the water surface. In this section, results from the current study are used to estimate the amount of aeration achievable with the bubbles entrained by plunging brine jets that are in contact with a moving water column. Consider the scenario in Figure 8, where a jet array of length comprising jets discharges vertically downward into a water column of depth where the ambient current is . As the jets penetrate the water surface, they collectively entrain oxygen from the atmosphere into the water column at an oxygen flux rate of = , where is the gas flowrate per unit length of the jet array.

Aeration of Water Column
Al Anzi et al. [22] proposed the use of plunging liquid jet reactors, arranged in a diffuser array, to promote the aeration of the underlying water via entrainment of air bubbles from above the water surface. In this section, results from the current study are used to estimate the amount of aeration achievable with the bubbles entrained by plunging brine jets that are in contact with a moving water column. Consider the scenario in Figure 8, where a jet array of length L comprising n jets discharges vertically downward into a water column of depth H where the ambient current is u a . As the jets penetrate the water surface, they collectively entrain oxygen from the atmosphere into the water column at an oxygen flux rate of Q g = q g L, where q g is the gas flowrate per unit length of the jet array.

Aeration of Water Column
Al Anzi et al. [22] proposed the use of plunging liquid jet reactors, arranged in a diffuser array, to promote the aeration of the underlying water via entrainment of air bubbles from above the water surface. In this section, results from the current study are used to estimate the amount of aeration achievable with the bubbles entrained by plunging brine jets that are in contact with a moving water column. Consider the scenario in Figure 8, where a jet array of length comprising jets discharges vertically downward into a water column of depth where the ambient current is . As the jets penetrate the water surface, they collectively entrain oxygen from the atmosphere into the water column at an oxygen flux rate of = , where is the gas flowrate per unit length of the jet array. The increase in dissolved oxygen concentration in the water column δ CDO (mass/volume) can be calculated as the oxygen flux J divided by the water flux, or: where A = Nπd 2 is the total surface area of N air bubbles (diameter d) in contact with the water column, k the oxygen mass transfer coefficient, c s the saturation concentration of oxygen in water, and c the dissolved oxygen concentration in water.
In turn, the number of bubbles in contact with the water can be estimated as equal to the bubble input rate . N multiplied by the contact time (t). If the bubbles are all of the same size, the bubble number rate in is given by πd 3 6 . N = Q g = q g L. The contact time of the bubbles in the water column can be estimated by t = H p /u b , where u b is the bubble rise velocity. Substituting the above relations into Equation (12) yields: Based on Equation (13) above, the aeration of the water column depends on the following four factors, each embodied by a variable present in the expression. Each of these four factors are discussed quantitatively for the plunging jet below.
A. The gas entrainment due to the plunging jet impinging across the water surface (q g ). B. The bubble size distribution (and bubble velocity) due to the plunging jet (u b and d).
C. Depth of bubble penetration (H p ). D. Mass transfer coefficient (k).

A. Gas Entrainment
The empirical relation for gas entrainment rate for an unconfined plunging jet is given by Equation  (10) can be applied to predict their gas entrainment.
Miwa et al. [39] compiled many predictive formulas for the air entrainment rate for neutrally buoyant plunging jets, and also offered an empirical relation for the air entrainment rate as a function of the Laplace length scale La = σ/∆ρ gb g and the Weber number. However, the current experimental parameters are outside the range of applicability of their results.

B. Bubble Size Distribution and Rise Velocity
The size distribution of the bubbles created by the plunging jet determines both the total interfacial surface area across which the oxygen mass transfer can occur, and the contact time of the bubbles as they ultimately rise back to the surface from their maximum penetration depth.
The bubble shape depends on the Reynolds number, Eötvös number Eo = g∆ρ al d 2 /σ (ratio of buoyancy and surface tension), and the Morton number M = gµ 4 ∆ρ gb /ρ b 2 σ 3 (combination of buoyancy, surface tension, and viscosity, from [43]. For experiments presented in this work, the expected regime encountered for plunging jet entrained bubbles is the ellipsoidal bubble shape regime. Based on qualitative observations of size distributions of bubbles generated by freshwater, seawater, and saltwater plunging jets, seawater and saltwater plunging jets may generate finer bubbles than freshwater plunging jets because of (1) lesser entrainment of large-sized bubbles, (2) greater entrapment of fine bubbles, (3) less bubble coalescence in seawater, and (4) surface tension causing finer bubbles to be more stable in the water column. These effects may cause the fine bubbles to be more likely to be pushed to the bottom of the tank in the case of seawater and saltwater plunging jets [44]. However, the experiments described by Chanson et al. [44] only used liquids that plunged into the receiving water of the same density.
Although it is expected that the plunging jet will produce a size distribution of bubbles, a wide range of sizes of air bubbles rising in contaminated (i.e., not pure) water (expected to behave similarly to air bubbles in seawater) rise at a fairly constant velocity of roughly 20 cm/s, especially within the ellipsoidal shape range. For the current analysis, u b =20 cm/s is assumed, independent of the bubble size.

C. Depth of Bubble Penetration
Equation (4a) can be used to estimate the depth of the bubble penetration. H p affects aeration because there is an increased contact time between the two phases before the bubbles ultimately rise back to the surface, which increases the total mass transfer. A larger bubble penetration depth is, therefore, expected to increase the aeration of the water column.

D. Mass Transfer Coefficient
The mass transfer coefficient of oxygen from the gas phase to the water phase (k) can be found in literature (e.g., [43]). Empirical functions of k depend on the molecular diffusivity of oxygen in water D O 2,w , bubble diameter (d), the density difference of the bubble and continuous phases (∆ρ g /ρ w ), and the bubble shape regime. For ellipsoidal bubbles, the mass transfer coefficient (from [43]) is given by: where units of k and d are m/s and cm, respectively, D O 2 ,w = 2 × 10 −9 m 2 /s and ∆ρ gl /ρ w ≈ 1.

E. Example of Aeration Prediction
To estimate the increase of dissolved oxygen for a typical plunging jet, the following estimates of the parameters were used as typical for a brine discharge offshore of Kuwait [45]: , yielding q g = 2.2 × 10 −3 m 2 /s. (11) results in δc DO = 2.7 mg/L. The estimated dissolved oxygen concentration increase is therefore significant, more than half of the DO sag of 6.8 − 2 = 4.8 mg/L.

Substituting the above values into Equation
The estimation of δc DO from Equation (11) clearly depends on many uncertain parameters. Therefore, a preliminary sensitivity analysis has been presented in Table 3. The sensitivity analysis showed that δc DO is strongly dependent on the bubble size generated by the plunging jets. For example, halving the bubble size increases δc DO by a factor of 4 from the base case-resulting in a predicted δc DO greater than saturation concentration of oxygen in water, suggesting that super saturation of oxygen may be achieved if the plunging jet comprises of very fine bubbles. Table 3. Sensitivity of dissolved oxygen concentration change (δc DO ) to a variety of parameters. For all cases, the following quantities were kept constant: H = 3 m; c s = 6.8 mg/L (25 • C, s = 35 psu); L = 100 m; Q 0 = 1 m 3 /s; D O 2 ,w = 2 × 10 −9 m/s 2 . Bold values indicate a change from the base case. Equation (4a) was used to determine the penetration depth and Equation (9) was used to determine the air entrainment rate. Asterisks denote that the dissolved oxygen introduced is greater than saturation.  As shown in Table 3, the factors decreasing aeration the most are the factors which decrease the contact time between air and the water column (an increased ambient current, a large nozzle diameter, and a short jet length). Conversely, a decrease in the bubble size greatly increases the overall contact area between the phases, and therefore increases the expected aeration. Using a box model that estimates the aeration potential of an array of plunging jet in the water column, the jets are predicted to introduce 2.7 mg/L of DO.

Conclusions
This study presented experimental observations of vertical unconfined plunging brine jets discharging in quiescent freshwater or saltwater where the plunging jets are of the same or higher density than receiving water. The measurements from this study showed an increase in bubble penetration depth with an increase in jet density difference. Bubble plume widths were generally higher for bubbles plunging in a saltwater ambient. Plunging jets were shown to achieve dilutions that are comparable to dilutions for jets or plume discharging just below the water surface. Measured dilutions showed no effect of receiving water salinity, increased slightly with increase in brine salinity, and showed a non-monotonic trend with increasing jet length.
There are many factors that may affect the efficiency of plunging jets for aerating the water column, namely, gas entrainment across the water surface, bubble size distribution, bubble penetration, and mass transfer coefficient. Each of the effects described above were quantified and a box model case study was presented with predictions of the resulting aeration from different scenarios. The study concludes that generating finer bubbles within the water column will have the greatest effect on aeration of the water column.
The findings from the paper can be used to estimate design parameters for the scale up of plunging jets for discharge of desalination brine with enhanced dissolved oxygen, and to inform further work on the bubble dynamics of the plunging jet within the water column. Further experimental measurements of bubble size distribution and water column dissolved oxygen could be performed to confirm the mass transfer and aeration capability of these plunging jets (single and arranged in a multiple diffuser assemble) at the laboratory scale.