The Effect of Variations of Flow from Tributary Channel on the Flow Behavior in a T-Shape Confluence

Channel confluences are of the common structures in fluid transport channels. In this study, a series of numerical simulations were performed, utilizing a 3D code to investigate the reaction of the flow parameters and vortical structure to the variations in flow discharge and its Froude number from both main channel and tributary branch in a T-shape junction. The code was calibrated with the experimental data. Parameters, including the velocity, the turbulence energy, stream surface profile, head losses, and the transverse flow motions, were considered in different situations. It was concluded that increasing the ratio of discharge of flow from side-channel to the main channel (Q*) increased the area and power of the recirculation zone, as well as the width of separation plate downstream of the confluence, while it reduced the area of the stagnation zone (or the wake vortex) within the side-channel. It was also indicated that increasing the discharge ratio from side-channel resulted in an increase in the upstream water level in the main channels, which was dependent on the upstream discharge. Record Type: Published Article Submitted To: LAPSE (Living Archive for Process Systems Engineering) Citation (overall record, always the latest version): LAPSE:2020.0883 Citation (this specific file, latest version): LAPSE:2020.0883-1 Citation (this specific file, this version): LAPSE:2020.0883-1v1 DOI of Published Version: https://doi.org/10.3390/pr8050614 License: Creative Commons Attribution 4.0 International (CC BY 4.0) Powered by TCPDF (www.tcpdf.org)


Introduction
A confluence is considered to be a point where two various flow characteristics join together and make a single flow, which causes significant turbulence characteristics. Confluences and diversions are the most common occurrences along both natural rivers and artificial channels. Different parameters, such as geometrical properties (i.e., channel cross-section, the angle between the branch and the main channel), and hydraulic parameters (i.e., Froude number of the flows, the discharge ratio of the branch to the main channel), have been found to be effective on the flow behavior and its effect on the loose bed of earth channels. The interaction between the main channel and the side inflow creates a 3D flow pattern with separation zones at both side and main channels. It also causes a significant increase in flow speed along with a drop in flow surface, which leads to a large amount of scouring at the opposite side of junction point and sedimentation adjacent to the inner wall of the main channel [1]. Generally, a comprehensive understanding of the sediment and junction or the interaction and behavior of diversion structures requires a thorough knowledge of the three-dimensional flow behavior around these structures [2].
Numerous investigations have been carried out to study the flow pattern through the junction of the channel, and there is a plethora of experimental, analytical, and numerical research on different aspects of flow behavior in such structures. From the very beginning research performed on confluences, various researchers, such as Taylor [3], Webber and Greated [4], Shumate [5], Mignot et al. [6], have considered the ratio of side-channel to main channel discharge as a key parameter of the flow behavior in channel confluences and diversions. Webber et al. [7] studied the pattern of the flow in a 90 • open channel via conducting a series of laboratory experiments. They considered the flow behavior in different ratios of side-channel discharge to the total channel varying from 0.08 to 0.9 and discussed the effect of discharge ratios on flow patterns in the channel. They provided significant data on the velocity domain, turbulence stresses, and water surface elevations. Liu et al. [8] experimentally studied the flow patterns and sediment transport in a confluence. They stated that increasing the ratio of tributary flow to total flow discharge restricted the development of the separation zone. They also suggested that the ratio of bedload transport and the scouring intensity relied mainly on the discharge ratio from the tributary channel. Zhang et al. [9] used 3D k-ω turbulence model and the Reynolds averaged Navier-Stokes equations to investigate the effect of the ratio of the flow in upstream of the main channel to the overall flow downstream on the separation zone's shape and size, the contraction coefficient, and cross-sectional mean flow angle for a 90 • open-channel with the same width. They concluded that larger discharge ratios of side-channel to the main channel might result in less uniformly distribution of flow and velocity profiles. They also stated that despite the variations in length and width of the recirculation zone in the main channel, its shape index remained almost constant. By utilizing a hybrid Reynolds averaged Navier Stokes (RANS)-based large eddy simulation (LES) model, Zeng and Li [10] simulated flow in an open-channel T-junction. Their results were far accurate and precise compared to the RANS results, which clearly indicated the preference of this model. Ting et al. [11] used a 3D K-ω model to simulate the flow for the 90 • open-channel junction with the same width and investigated the influence of the ratio of the discharge on the separation zone, the cross-sectional mean velocity of the flow, and the coefficient of the contraction. Riley and Rhoads [12] investigated both flow pattern and morphology of the bed; furthermore, they studied the impacts of the curvature of the channel and the angle of the confluence on the confluence morpho-dynamics. Mohammadiun et al. [1] performed a series of numerical simulations and studied the effect of different geometrical modifications on the flow behavior in a 90-degree junction. Modifications included adding converging and diverging arcs at corners of the junction, as well as the elimination of recirculation zones by making arc shape modifications in the main channel based on the idea of streamlinization. Their results indicated that the modifications successfully improved the flow pattern and eliminated the formation of recirculation areas. Liu et al. [13] investigated the flow-sediment interaction in confluences with different junction angles. Their results showed that the sediment transport rate increased with increasing the inflow ratio and junction angle. Biswal et al. [14] studied the hydraulics of flow in a right-angled confluence in the case of compound channels. They stated that the flow in such channels was highly different from that in simple channels. Sharifipour et al. [15] used a k-ω turbulence model to simulate the effect of tributary channel width on the flow characteristics in the main channel. It was indicated that increasing the ratio between tributary channel widths to main channel width from 0.5 to 2 might lead to a decrease in the effect of tributary flow on the characteristics of the main channel flow. Schindfessel et al. [16] numerically investigated the confluence flow patterns with a predominant fraction (90%<) of the total flow consisted of the tributary inflow. It was found that the tributary flow collided the opposing bank in the predominant tributary flow, triggering a recirculating eddy in the upstream channel of the confluence. Wang et al. [17] experimentally investigated the flow confluences on the natural rivers. They studied the flow junction impact on the stage-discharge relationship at the confluence. Sui and Huang [18] numerically simulated the separation zone. This quasi-2-D simulation was done in the high Reynolds number and using the LES model. The comparison between numerical results and the experimental measurements showed perfect acceptance. Ramos et al. [19] compared the accuracy of flat and rigid-lid techniques using the LES turbulence model for the simulation of free-surface flow in channel junctions. In another research, Ramos et al. [20] studied the effect of bed elevation discordance of the flow patterns and flow turbulent structure in a right-angled confluence.
Reviewing literature reveals that there is a significant amount of data on the different aspects of flow behavior in different geometrical and hydraulic conditions of junction's performance; yet, various researchers are working on different aspects of the flow and sediment transport within confluences. Several researchers have considered the effect of the ratio of tributary channel discharge to total channel discharge. They have stated this parameter as the governing parameter in sediment scour and flow structure at confluences. Nevertheless, different combinations of discharges from the tributary channel and main channel with a constant total Froude number-at downstream section after confluence-could be formed. The present study made efforts to explain the differences in flow patterns and their power induced by variations in discharges coming from the tributary channel and main channel in conditions with constant or varying total discharge. Sediment transport and sedimentation is a direct consequence of flow patterns and vortex structure above the bed and close to sensitive points. Hence, this paper made efforts to investigate the 3D behavior of flow at a T-shape junction in different ratios of side channels discharge to the main channel.

Laboratory Model
The present study utilized the experimental data of Webber et al.

Numerical Model
This study utilized the FLOW-3D numerical model as a solver of the Navier-Stokes equation to calculate fully 3D flow patterns at the channel junction.

Governing Equations
The governing equations comprised the momentum and continuity equations. The continuity equation, regardless of the density of the fluid in the form of Cartesian coordinates x, y, and z, was as follows: where u, v, and w are the components of velocity in the x, y, and z directions, respectively; A x , A y , and A z are the surface flow fractions in the x, y, and z, respectively; V F is the flow volume fraction; ρ is the fluid density; t is time, and R sor is the mass source. Equation (2) is the momentum equation in three dimensions: where G x , G y , and G z are the accelerations due to the gravity in the x, y, and z directions, respectively; f x , f y , and f z are the accelerations due to viscosity in the x, y, and z directions, respectively.

Turbulence Modeling
The utilized turbulence models in this study were the k-e, k-w, and k-ε renormalized group (RNG) models. The k-ε model is highly sophisticated, widely-used and consisted of two transport equations: one for turbulent kinetic energy k, and the other for its dissipation ε [20]. The 3D governing equations could be expressed as follows [21,22]: where µ t represents eddy viscosity, G b and G k are generations of the turbulent kinetic energy due to buoyancy and mean velocity gradients, respectively; Y M is the contribution of the fluctuating dilatation incompressible turbulence to the overall dissipation rate; C 1ε , C 2ε , C 3ε , and C µ are constants (1.44, 1.92, 0.09, and 0.09 respectively); σ k and σ ε are turbulent Prandtl numbers for k and ε (1.0 and 1.3, respectively). The RNG model used equations analogous to the equations for the k-ε model. However, constants that were found empirically in the standard k-ε model were derived explicitly in the RNG model. Generally, the RNG model has wider applicability than the standard k-ε model. In particular, the RNG model is known to present more accurately for low-intensity turbulence flows and flows with strong shear regions [2].

Boundary Conditions and Gridding
In this study, a structured non-uniform block of mesh was used to simulate the main and side channels. The meshes were more intense in the vicinity of the confluence to increase the accuracy of computations. Boundary conditions included inflow discharge for the main channel and the lateral channel entrance, constant level outflow for the main channel outlet, wall for bed and the right boundaries of the main channel, and symmetry for the top border. Other required wall conditions were produced using solid components. Given that using a single mesh domain caused two large empty domains at the sides of the lateral channel, these parts were eliminated from calculations using void or space components. Considering restrictions in the available processing power, the main mesh block was non-uniformly identified in order to increase the accuracy at the vicinity of the junction and reducing the number of mesh cells in less important areas. Considering the limitations in processing power (Core i7 CPU-turbo speed of up to 4.0 GHz), a trial and error process was performed to find a proper balance between the consumed simulation time and the resulted accuracy of the simulations. Different mesh sizes were tested; however, a severe increase in computational time-due to decreased mesh sizes-as well as the increased number of cells, restricted the use of meshes with side dimensions less than 4 mm. Thus, a mesh domain, with 17 m length, 0.6 m width, and 0.4 m height and consisting of 2,524,480 cells, was opted to simulate the channels. The mesh domain was composed of 320 cells in the x-direction, 161 in the y-direction, and 49 cells in the z-direction. Neglecting the void spaces, all active cells were 1,839,166. As mentioned above, the cell sizes were reduced near the walls and bed to 4 mm. Besides, the same cell size was identified at the corners of the junction and the space of flow mixing in order to increase the accuracy of the calculations. The plan and side view of the mesh domain, as well as the boundary conditions, are plotted in Figure 1.

Model Verification
A first-order expansion of the momentum advection, the implicit generalized minimal residual method (GMRES) with an allowed number of 120 internal iterations in each step, an initial time step of 1E-6, and the fractional area-volume obstacle representation (FAVOR) method for estimation of the free surface were used to estimate the momentum and pressure field. The verification of the model was performed using experimental data from Webber et al.

Model Verification
A first-order expansion of the momentum advection, the implicit generalized minimal residual method (GMRES) with an allowed number of 120 internal iterations in each step, an initial time step of 1E-6, and the fractional area-volume obstacle representation (FAVOR) method for estimation of the free surface were used to estimate the momentum and pressure field. The verification of the model was performed using experimental data from Webber et al.

Model Verification
A first-order expansion of the momentum advection, the implicit generalized minimal residual method (GMRES) with an allowed number of 120 internal iterations in each step, an initial time step of 1E-6, and the fractional area-volume obstacle representation (FAVOR) method for estimation of the free surface were used to estimate the momentum and pressure field. The verification of the model was performed using experimental data from Webber et al.

Simulation Scenarios
A total number of seven different conditions were simulated with different ratios of side-channel to main channel discharge. Table 1 presents the simulation conditions and discharge ratios. Various researchers have concluded that the ratio between discharges is a key parameter. Still, some studies have also stated that the cross-sectional area of channels is a dominant parameter, affecting the flow behavior. On the other hand, changing the width of either of channels with the same geometry will have a significant effect on the flow structure. Hence, considering the ratio between Froude number of flow in side-channel (Frsc) to Froude number of flow in the main channel (Frmc), as well as the Froude number of the overall flow after the confluence (FR*), as varying parameters in simulation scenarios would provide better understanding and could result in better and more generalized conclusions. Qsc is the discharge from the side-channel. Qmc is the discharge from the upstream main channel. Qt is the total discharge after confluence. Q* is ratio between discharge of tributary flow and the main channel flow. FR* is the ratio between Froude number of tributary flow and the main channel flow.

Simulation Scenarios
A total number of seven different conditions were simulated with different ratios of side-channel to main channel discharge. Table 1 presents the simulation conditions and discharge ratios. Various researchers have concluded that the ratio between discharges is a key parameter. Still, some studies have also stated that the cross-sectional area of channels is a dominant parameter, affecting the flow behavior. On the other hand, changing the width of either of channels with the same geometry will have a significant effect on the flow structure. Hence, considering the ratio between Froude number of flow in side-channel (Fr sc ) to Froude number of flow in the main channel (Fr mc ), as well as the Froude number of the overall flow after the confluence (FR*), as varying parameters in simulation scenarios would provide better understanding and could result in better and more generalized conclusions.

Velocity Distribution
The distribution of the stream-wise component of velocity (U) and the normal component of velocity (V) at x-y plane around the confluence structure and along the main channel are presented in Figures 4 and 5, respectively, for three different cases, including case (1) with Fr* = 0.28, case (4) with Fr* = 44, and case (7) with Fr* = 0.59. The results were presented in four different depths, including z = 0.01 m (near the bed), z = 0.08 m, z = 0.17 m, and z = 0.27, close to the water surface. The other cases were not shown here due to the similarity between results. As expected, there was a significant difference in maximum velocity between case 1 (Figure 4a-d) and case 2 (Figure 4e-h). This difference was a direct consequence of increased discharge from side-channel in case 4 compared to case 1, as well as the reduced width of the main channel after confluence due to the increase in the size of the recirculation zone, and agreed with the increased Froude number of total discharge. Unexpectedly, comparing case 1 to case 7 indicated that the velocity of flow after the confluence in case 7 was higher than that in case 1, especially near the bed and close to the surface, although the overall discharge in case (7) was less than that in case 1 (up to 7%). Besides, the high-velocity zone moved downstream. Yet, it could be seen that the velocity gradient toward the recirculation zone at the downstream corner of confluence was more intense in case (1) compared to the case (7), which was caused by increased discharge and total Froude number of the tributary inflow.      The same conclusion could be inferred from the distribution of the normal component of velocity in Figure 5. The V velocity gradient in case 7 was more intense than that in case (1) despite the reduced discharge and total Froude number of the main channel after confluence (Fr t = 0.31 for case 7) in comparison with other cases. It could be interpreted through the increased width of the recirculation zone in the channel in case (7) compared to other cases. Generally, it could be inferred that increasing the Froude number of the side-channel or the ratio between the Froude numbers of side-channel to the Froude number of the main channel might result in higher velocity gradients near the bed, which might result in increased sediment scour, as well as increased sedimentation due to the increased area of recirculation zone. Figure 6 presents the effect of variations in discharges from the side and main channel on the stream-wise and normal component of velocity for all test cases in the main channel after confluence at x = 7.91 m at the downstream edge of the junction and at x = 8.35 m (0.5* Wc downstream of the junction) close to the flow surface. Generally, increasing the total discharge in cases 1 to 4 resulted in an increased stream-wise velocity (U) throughout the channel cross-section at the junction point. But, decreased discharge of the main channel in cases 5 to 7 resulted in a reduced velocity profile. Yet, similarities between cases 1 to 6 and cases 2 to 5 could be seen due to the same total discharge after confluence (Figure 6a). It could be inferred that various combinations of side channel discharge (Q sc ) and that of main channel (Q mc ) could result in relatively similar stream-wise velocity distribution (Figure 6a,c), yet, the main difference, which caused variations in flow patterns, was caused by significant fluctuations in V. Considering the normal component of velocity (V) at the junction section (Figure 6b), it was seen that for cases 1 to 4, increasing the discharge from side-channel resulted in an increased V at the confluence. Besides, it was evident that the V velocity gradient decreased as the flow moved toward the opposite wall. Figure 6d indicates the cross-sectional variations of the normal component of velocity (V) at 0.5 Wc downstream of the confluence. Velocity gradients were significantly reduced compared to the junction section, yet there were distinct variations in different cases. Even though that the differences had no clear trend, yet, it could be inferred that by changing conditions from case 1 to case 7 (increasing the Fr*), the absolute value of V velocity close to the inner wall increased, which was equal to the increased width and length of recirculation zone.

Flow Surface Profile
Comparisons of the flow surface profiles for different cases at two transverse stations are presented in Figure 7. Generally, it was indicated that increasing Fr* might increase the upstream water level before the confluence section. It should be noted that the upstream water level for cases 6 and 7 was higher than that for case 1 (Figure 7a). A surface drop at the tail of the confluence section was observed, which indicated sharper declines for cases with higher total discharges. Comparing similar cases, 1 with 6 and 2 with 5 (cases with same total discharge), indicated that the upstream level for cases with higher Fr* (cases 5 and 6) was higher than that for cases with less Fr* (1 and 2) despite the similar discharges. Hence, it could be concluded that the upstream water level could be higher in cases with greater Fr* but same or less total discharge, which indicated the governing role of discharge ratio in flow behavior, which surpassed the discharge values.
The difference in the downstream water level after the confluence was more significant. Surface drop in cases with less Fr* was smaller and milder in profile, which suggested higher losses in cases with increased Fr*. Considering the normal component of velocity (V) at the junction section (Figure 6b), it was seen that for cases 1 to 4, increasing the discharge from side-channel resulted in an increased V at the confluence. Besides, it was evident that the V velocity gradient decreased as the flow moved toward the opposite wall. Figure 6d indicates the cross-sectional variations of the normal component of velocity (V) at 0.5 Wc downstream of the confluence. Velocity gradients were significantly reduced compared to the junction section, yet there were distinct variations in different cases. Even though that the differences had no clear trend, yet, it could be inferred that by changing conditions from case 1 to case 7 (increasing the Fr*), the absolute value of V velocity close to the inner wall increased, which was equal to the increased width and length of recirculation zone.

Flow Surface Profile
Comparisons of the flow surface profiles for different cases at two transverse stations are presented in Figure 7. Generally, it was indicated that increasing Fr* might increase the upstream water level before the confluence section. It should be noted that the upstream water level for cases 6 and 7 was higher than that for case 1 (Figure 7a). A surface drop at the tail of the confluence section was observed, which indicated sharper declines for cases with higher total discharges. Comparing similar cases, 1 with 6 and 2 with 5 (cases with same total discharge), indicated that the upstream level for cases with higher Fr* (cases 5 and 6) was higher than that for cases with less Fr* (1 and 2) despite the similar discharges. Hence, it could be concluded that the upstream water level could be higher in cases with greater Fr* but same or less total discharge, which indicated the governing role of discharge ratio in flow behavior, which surpassed the discharge values.

Streamlines and Flow Patterns
Generally, three significant vertical structures formed at a confluence, including a recirculation zone at the downstream corner of the junction adjacent to the inner wall of the main channel, a stagnation zone (or a wake vortex) at the inner corner of the tributary channel (Figure 8a,b), and the spiral movement and mixing of the tributary flow with the main channel flow (Figure 8c).  Figure 9 presents the streamlines for cases 1, 4, and 7 in the mentioned depths in the previous section. It was clear that the size of the recirculation zone at the main channel increased as the Froude ratio increased from case 1 to case 7, which was responsible for increased maximum velocity in case 7 compared to case 1 despite less discharge. The results were consistent with the findings of other The difference in the downstream water level after the confluence was more significant. Surface drop in cases with less Fr* was smaller and milder in profile, which suggested higher losses in cases with increased Fr*.

Streamlines and Flow Patterns
Generally, three significant vertical structures formed at a confluence, including a recirculation zone at the downstream corner of the junction adjacent to the inner wall of the main channel, a stagnation zone (or a wake vortex) at the inner corner of the tributary channel (Figure 8a,b), and the spiral movement and mixing of the tributary flow with the main channel flow (Figure 8c).

Streamlines and Flow Patterns
Generally, three significant vertical structures formed at a confluence, including a recirculation zone at the downstream corner of the junction adjacent to the inner wall of the main channel, a stagnation zone (or a wake vortex) at the inner corner of the tributary channel (Figure 8a,b), and the spiral movement and mixing of the tributary flow with the main channel flow (Figure 8c).  Figure 9 presents the streamlines for cases 1, 4, and 7 in the mentioned depths in the previous section. It was clear that the size of the recirculation zone at the main channel increased as the Froude ratio increased from case 1 to case 7, which was responsible for increased maximum velocity in case 7 compared to case 1 despite less discharge. The results were consistent with the findings of other  Figure 9 presents the streamlines for cases 1, 4, and 7 in the mentioned depths in the previous section. It was clear that the size of the recirculation zone at the main channel increased as the Froude ratio increased from case 1 to case 7, which was responsible for increased maximum velocity in case 7 compared to case 1 despite less discharge. The results were consistent with the findings of other researchers, such as [23] and [7]. With increasing the Froude ratio, loss of flow energy at the junction increased, as the momentum of tributary flow, which needed to be turned into the main channel flow, was increased. Consequently, higher energy drop might lead to a more significant separation zone [9]. Besides, an important dividing line starting from the upstream edge of the tributary channel was visible (where streamlines gather), which was located at different stations in different elevations of flow, indicating the section of shear flow at the mixing interface of tributary flow with the main stream. researchers, such as [23] and [7]. With increasing the Froude ratio, loss of flow energy at the junction increased, as the momentum of tributary flow, which needed to be turned into the main channel flow, was increased. Consequently, higher energy drop might lead to a more significant separation zone [9]. Besides, an important dividing line starting from the upstream edge of the tributary channel was visible (where streamlines gather), which was located at different stations in different elevations of flow, indicating the section of shear flow at the mixing interface of tributary flow with the main stream. A significant stagnation zone at the upstream corner of the confluence could also be observed in lower discharge ratios, which was vanished as the discharge from side-channel increased. However, it could be noted that the stagnation zone in lower discharge ratios was, in fact, a wake vortex and is plotted in Figures 8b and 9c,d. As stated, the width and length of the recirculation zone increased with increasing the Froude ratio. Figure 10 represents the cross-sectional variations of the recirculation zone in each case from the bed to the flow surface. It was indicated that the cross-sectional area of the recirculation zone increased from bed to water surface in all cases. Besides, it was shown that the overall and the width of the recirculation zone at all elevations increased with increasing the Froude ratio from case 1 to case 7. Figure 11 presents the length and width of the recirculation zone, as well as the shape index for each case and different Fr* values. It was indicated that the length and width of the recirculation area increased with different slopes, resulting in a slightly increased shape index with increasing the Fr*, which could be assumed horizontal. A significant stagnation zone at the upstream corner of the confluence could also be observed in lower discharge ratios, which was vanished as the discharge from side-channel increased. However, it could be noted that the stagnation zone in lower discharge ratios was, in fact, a wake vortex and is plotted in Figures 8b and 9c,d. As stated, the width and length of the recirculation zone increased with increasing the Froude ratio. Figure 10 represents the cross-sectional variations of the recirculation zone in each case from the bed to the flow surface. It was indicated that the cross-sectional area of the recirculation zone increased from bed to water surface in all cases. Besides, it was shown that the overall and the width of the recirculation zone at all elevations increased with increasing the Froude ratio from case 1 to case 7.

Turbulent Kinetic Energy
The distribution of the turbulent kinetic energy (TKE) at the transverse sections of x = 7.91 and x = 8.35 for all simulation cases are presented in Figure 12. It was indicated that a high positive gradient of TKE was formed at the downstream corner of the tributary channel. It could be seen that the maximum value of TKE was related to the case 4 with the maximum amount of total discharge and Froude number; for other cases, with decreasing the total flow discharge, the maximum values of TKE were reduced. It could be noted that the minimum value of TKE among all cases was related to case 1 and case 2. A comparison between case 1 and 2 to case 6 and 7 showed that the maximum value of TKE in cases 6 and 7 with less total discharge and higher Fr* was higher than that for cases 1 and 2, suggesting larger recirculation zones and higher ratios of turbulence despite less discharge (see Figure 12a). Considering the turbulent energy at x = 0.5 Wc downstream of the junction, a relatively similar trend with that of the previous transverse section could be observed.  Figure 11 presents the length and width of the recirculation zone, as well as the shape index for each case and different Fr* values. It was indicated that the length and width of the recirculation area increased with different slopes, resulting in a slightly increased shape index with increasing the Fr*, which could be assumed horizontal.

Turbulent Kinetic Energy
The distribution of the turbulent kinetic energy (TKE) at the transverse sections of x = 7.91 and x = 8.35 for all simulation cases are presented in Figure 12. It was indicated that a high positive gradient of TKE was formed at the downstream corner of the tributary channel. It could be seen that the maximum value of TKE was related to the case 4 with the maximum amount of total discharge and Froude number; for other cases, with decreasing the total flow discharge, the maximum values of TKE were reduced. It could be noted that the minimum value of TKE among all cases was related to case 1 and case 2. A comparison between case 1 and 2 to case 6 and 7 showed that the maximum value of TKE in cases 6 and 7 with less total discharge and higher Fr* was higher than that for cases 1 and 2, suggesting larger recirculation zones and higher ratios of turbulence despite less discharge (see Figure 12a). Considering the turbulent energy at x = 0.5 Wc downstream of the junction, a relatively similar trend with that of the previous transverse section could be observed.

Turbulent Kinetic Energy
The distribution of the turbulent kinetic energy (TKE) at the transverse sections of x = 7.91 and x = 8.35 for all simulation cases are presented in Figure 12. It was indicated that a high positive gradient of TKE was formed at the downstream corner of the tributary channel. It could be seen that the maximum value of TKE was related to the case 4 with the maximum amount of total discharge and Froude number; for other cases, with decreasing the total flow discharge, the maximum values of TKE were reduced. It could be noted that the minimum value of TKE among all cases was related to case 1 and case 2. A comparison between case 1 and 2 to case 6 and 7 showed that the maximum value of TKE in cases 6 and 7 with less total discharge and higher Fr* was higher than that for cases 1 and 2, suggesting larger recirculation zones and higher ratios of turbulence despite less discharge (see Figure 12a). Considering the turbulent energy at x = 0.5 Wc downstream of the junction, a relatively similar trend with that of the previous transverse section could be observed.

Conclusion
The present study investigated the effect of various combinations of flow from the tributary channel to the discharge from the main channel at a confluence. Seven different cases were simulated using a 3D CFD code, and the results were validated using the experimental data from the literature. Results indicated that the ratio between the Froude numbers of flows from the tributary channel to the flow of the main channel had a significant effect on the flow patterns and momentum exchange at the confluence.
Results indicated that the velocity gradient toward the recirculation was more intense in cases with a higher ratio of flow from the side-channel to the flow from the main channel, even in cases with less total discharge. It was due to the increased discharge and total Froude number of the tributary inflow. Besides, increasing the discharge ratio from side-channel resulted in an increase in the upstream water level in the main channel so that the upstream level for a flow with less total discharge but higher discharge ratio from tributary channel could be higher than that in cases with higher discharges but less portion coming from the tributary.
It could also be concluded that increasing the discharge from the tributary channel might result in higher velocity gradients near the bed, which, in turn, increased scour as well as sedimentation due to the expanded area of the recirculation zone. It also increased the dimensions of the recirculation zone as well as its power.
Finally, increasing the total discharge increased the turbulent energy within the domain, while for cases with constant total discharge, the values of TKE were higher for cases with less discharge from the main channel and more discharge coming from the tributary.

Conclusions
The present study investigated the effect of various combinations of flow from the tributary channel to the discharge from the main channel at a confluence. Seven different cases were simulated using a 3D CFD code, and the results were validated using the experimental data from the literature. Results indicated that the ratio between the Froude numbers of flows from the tributary channel to the flow of the main channel had a significant effect on the flow patterns and momentum exchange at the confluence.
Results indicated that the velocity gradient toward the recirculation was more intense in cases with a higher ratio of flow from the side-channel to the flow from the main channel, even in cases with less total discharge. It was due to the increased discharge and total Froude number of the tributary inflow. Besides, increasing the discharge ratio from side-channel resulted in an increase in the upstream water level in the main channel so that the upstream level for a flow with less total discharge but higher discharge ratio from tributary channel could be higher than that in cases with higher discharges but less portion coming from the tributary.
It could also be concluded that increasing the discharge from the tributary channel might result in higher velocity gradients near the bed, which, in turn, increased scour as well as sedimentation due to the expanded area of the recirculation zone. It also increased the dimensions of the recirculation zone as well as its power.
Finally, increasing the total discharge increased the turbulent energy within the domain, while for cases with constant total discharge, the values of TKE were higher for cases with less discharge from the main channel and more discharge coming from the tributary.