Urban development near rivers has led to the increased release of municipal and industrial wastewater into water bodies. Although several studies have explored wastewater movement in rivers, further investigation is needed to understand the impact of effluent discharge on river health [1], [2]. The discharge of effluent into flowing water, termed jet in crossflow, has applications in industry, hydraulics, and the environment [3], [4], [5], [6].
The hydrodynamics of a jet in crossflow is inherently complex, especially when the jet is released in areas of the river with complex structures, such as a confluence channel junction. To the best of the authors’ knowledge, the mixing of a jet in a confluence zone has never been investigated. This study aims to investigate effluent discharges released in confluence areas to examine the effect of the junction on the outfall. First, the hydrodynamics of both jets in crossflow and confluences are reviewed briefly in the following sections.
Jet in crossflow
When a less-dense discharge is introduced into a denser fluid, it gives rise to a positive buoyant jet. Specifically, a positive buoyant surface jet is formed when discharged near the free surface. This phenomenon is commonly observed in both natural and industrial settings, such as the release of heated cooling water [7] and discharge of treated wastewater into water bodies [8]. This method is the most traditional and economical method for releasing pollutants into rivers or coastal areas. Due to their limited dilution factor, usually ten or below, surface discharges tend to pose greater environmental risks compared to submerged alternatives [9].
The dilution factor is a measure of the amount of pollution that has been diluted or reduced in concentration. The dilution factor is calculated using the following equation.
$$\:Dilution\:Factor=\frac{Volume\:of\:initial\:Concentration}{Volume\:of\:final\:Concentration}$$
1
In Eq. (1), "Volume of Initial Solution" refers to the volume of the concentrated pollution before dilution, and "Volume of Final Solution" refers to the total volume of the diluted pollution after dilution.
The trajectory of the jet is established by connecting the points with the maximum time-averaged concentration in each cross section (Fig. 1). The angle 𝜃 represents the inclination angle of the jet, signifying the direction angle of the instantaneous velocity along the central axis in relation to the crossflow direction (𝜃𝑜 is the initial angle, while 𝜃 denotes angles along the trajectory). Theoretically, 𝜃 reaches zero degrees at the end of the trajectory.
The length-scale parameters for buoyant surface jet classifications were introduced by Jirka et al. [10], and Chu and Jirka [11] provided an initial classification for buoyant surface jets in crossflows. A comprehensive classification system by Jones et al. [9] categorized surface jets into four main classes as either free jets, shoreline-attached jets, wall jets, or upstream intruding plumes, using the parameters defined in Tables 1 and 2 where a0 is the cross-sectional area of the jet, g is the gravitational acceleration, ρ is the water density, and Δρ0 is the difference in discharge density (Δρ0 = ρa − ρ0).
Table 1
Dynamic length scales of jets
Discharge length scale | \(\:{L}_{Q}={Q}_{0}/{M}_{0}^{1/2}\) |
Jet-to-plume length scale | \(\:{L}_{M}={M}_{0}^{3/4}/{J}_{0}^{1/2}\) |
Jet-to-crossflow length scale | \(\:{l}_{m}={M}_{0}^{1/2}/\stackrel{-}{{u}_{a}}\) |
Plume-to-crossflow length scale | \(\:{l}_{b}={J}_{0}/{\stackrel{-}{{u}_{a}}}^{3}\) |
Table 2
Flow characteristics of jets
Jet flow rate | \(\:{Q}_{0}=\stackrel{-}{{u}_{0}}{a}_{0}\) |
Jet momentum flux | \(\:{M}_{0}=\stackrel{-}{{u}_{0}}{Q}_{0}\) |
Jet buoyancy flux | \(\:{J}_{0}={Q}_{0}\:{g}_{0}^{{\prime\:}}\) |
Discharge buoyant acceleration | \(\:{g}_{0}^{{\prime\:}}=g(\varDelta\:{\rho\:}_{0}/{\rho\:}_{0})\) |
A crossflow buoyant jet can be categorized into four distinct asymptotic regimes based on whether buoyancy, momentum, or crossflow primarily influences the jet's mixing characteristics (Fig. 1) [12]:
-
Momentum Dominated Near-Field (y/lm < < 1): This regime occurs when the initial jet momentum significantly influences the jet discharge.
-
Momentum Dominated Far-Field (y/lm > > 1): In this regime, the crossflow effects become significant in the jet discharge.
-
Buoyancy Dominated Near-Field (y/lb < < 1): This regime is characterized by the initial plume buoyancy playing a crucial role in the plume discharge.
-
Buoyancy Dominated Far-Field (y/lb > > 1): The crossflow effects become prominent in the plume discharge.
where y represents the cross-channel distance, and lm and lb are the lengths of the near-field for a momentum-dominated jet or buoyancy-dominated plume, respectively (Fig. 1).
Early investigations into surface thermal discharge patterns were initiated by Dunn et al. [13], followed by computational modeling using mathematical models and turbulence closure techniques [14]. In the 1960s, research focused on surface buoyant jets, prompted by concerns over the environmental impact of heated cooling water discharges from power plants. Notable contributions during this period were laboratory experiments and integral jet models [10].
A surface jet discharging into a crossflow transition along its trajectory. In the "zone of flow establishment" immediately following the nozzle, the flow undergoes a transformation into a sheared jet-like velocity distribution. In this region, buoyancy has a minimal impact, and the flow is primarily governed by the momentum of the jet, conforming to a Gaussian distribution. As the momentum diminishes, the jet shifts towards the crossflow, entering a plume-like condition. Mixing in this zone is influenced by buoyancy effects and secondary flows, which collectively form a near-field zone. In the far field, mixing is predominantly governed by buoyant spreading or passive diffusion (Fig. 1) [8], [9].
Researchers have explored surface jets to enhance understanding of these phenomena. Anwar [15] conducted experiments to assess the vertical and horizontal isotherms of warm surface jets (73 ◦C) discharged from a rectangular channel into a cold crossflow (15 ◦C). The study revealed a strong depth dependence of the jet axis, and the mutual interaction between the buoyant surface jet and crossflow generated a zone with overshoot velocity at the outer edge of the jet. Additionally, a slow-moving zone was observed on the lee side of the buoyant jet, along with a return flow closer to the inner wall (near the discharge point) of the flume. Gawad et al. [16] investigated the effects of horizontal surface and submerged jets in crossflows, highlighting similar Gaussian distributions in the excess jet velocity and concentration in both the lateral and vertical directions for cases without strong bank reattachment and recirculation.
Gharavi et al. (2020) [8] conducted experimental measurements of the 3-D velocity field to enhance the understanding of the flow structure of surface jets under crossflow conditions. Stereoscopic Particle Image Velocimetry (PIV) was employed to measure the instantaneous spatial and temporal velocity distributions in a measurement zone downstream of the jet discharge point. This study identified the formation and evolution of vortices in the flow structure. This vortex in the surface jet in crossflow exhibited similarities to one half of the counter-rotating vortex pair (CVP) typical of submerged jets in crossflows. It was inferred that the water surface behaved as a plane of symmetry (Fig. 2).
Confluence Hydrodynamics
Channel confluences are hydrological systems in which two or more channels converge, resulting in complex flow structures. These structures can be challenging to understand and model due to the complexity of the flow physics involved. Over the years, researchers have attempted to gain a better understanding of the behavior of open-channel confluences using experimental methods [17], [18], [19], [20], [21], numerical methods [31], [32], [33], [34], and field experiments [26], [27], [28], [29]. Best (1987) [30] proposed a generalized model for river confluences that included six different zones (Fig. 3). These zones are characterized by different flow patterns and velocities. The first zone is the stagnation zone, where water from the tributary branch river is impeded by the main river, causing a build-up of water. The second zone is the deflection zone, where the flow from the branch river changes direction to align with that of the main river channel. The third zone is the separation zone, which is the area immediately after the junction close to the inner bank of the main river. The velocity in this zone is usually low, yet highly turbulent with flow recirculation. The fourth zone is the maximum velocity zone, which occurs in the contraction or acceleration zone between the shear layer and outer and inner banks of the main river. The fifth zone is the flow recovery zone, in which the flow begins to recover from the turbulence generated in the shear plane. The sixth and final zone is the shear plane, which occurs at the branch river entrance to the main river and is the primary source of turbulence generation due to differences in their velocity fields.
After the confluence junction, various river characteristics change significantly, including alterations in water surface and bed level. Additionally, a range of phenomena can be observed, such as the transport of sediment, dispersion of pollutants, erosion and deposition processes, and river migration [31].
Several key elements, such as junction angle, ratio of discharge or momentum between the main channel and tributary flows, bed level (accordant or discordant), channel geometry, and any differences in density between the incoming flows, can affect the morphodynamics of a confluence [32], [33], [34], [35]. While channel confluences have been studied for several decades, it is evident that examination of these details continues to accelerate [36], [37], [38], [39]. Significant research efforts have focused on investigating the phenomenon of flow mixing at confluences [27], [40], [41]. Furthermore, there has been considerable interest in studying the relationship between sediment transport and the morphology of river systems [29], [42], [43].
According to previous studies, confluences can be categorized based on their width-to-depth ratio as either small (W/H < 10), medium (10 < W / H < 50), or large scale (W / H > 50) [29], [33]. Recently, researchers have explored the possibility of applying findings from small laboratory and field experiments to large rivers, thus this scaling categorization should be considered [44]. It has been suggested that for large width-to-depth ratios, the influence of secondary currents is diminished due to increased roughness caused by the shape of the riverbed. Consequently, it may be challenging to extend observations from small to large-scale confluences [33].
Objective
Despite numerous investigations into confluence channels or jets in crossflow, there is still a gap in existing studies. To the best of the authors’ knowledge, no prior study has explored the dynamics of positively buoyant surface jets in confluence regions. Recognizing this research gap, our study aims to provide a comprehensive investigation of the behavior and characteristics of positively buoyant surface jets in the complex context of confluence regions. We inspect not only the fundamental features of these jets but also their concentration patterns and trajectories within the intricate dynamics of river confluences. Different jet in cross flow categories that can occur when a buoyant surface jet is discharged into a confluence are examined, by considering variable jet location, momentum, and buoyancy. Specifically, laboratory flume experiments were conducted in an open-channel confluence with injection of an effluent jet into either the tributary or main channel upstream of the confluence zone. For various jet and confluence flow scenarios, laser-induced fluorescence was utilized to evaluate the mixing of the effluent within the confluence zone. Through a systematic exploration of this uncharted territory, we seek to contribute valuable insights into the understanding of the flow patterns, turbulence, concentration dynamics, and trajectory behaviors associated with positively buoyant surface jets in confluence areas.