Flow Topology in the Confluence of an Open Channel with Lateral Drainage Pipe

The purpose of this paper is to develop design guidelines for flood control channel height in the vicinity of the confluence of a submerged drainage pipe and a flood control channel. The water exchange in the confluence of an open channel with a lateral drainage pipe produces unique hydraulic characteristics, ultimately affecting the water surface elevation in the channel. An accurate prediction of the water surface elevation is essential in the successful design of a high-velocity channel. By performing several experiments, and utilizing a numerical model (FLOW-3D), this study investigated the impact of submerged lateral drainage pipe discharges into rectangular open channels on flow topology in the confluence hydrodynamics zone (CHZ). The experiments were conducted in different flume and junction configurations and flow conditions. Moreover, the simulations were performed on actual size channels with different channel, pipe, and junction configurations and flow conditions. The flow topology in the CHZ was found to be highly influenced by the junction angle, as well as the momentum ratios of the channel flow and the pipe flow. The findings of this study were used to develop conservative design curves for channel confluences with lateral drainage pipe inlets. The curves can be used to estimate water surface elevation rise in different channel and pipe configurations with different flow conditions to determine the channel wall heights required to contain flows in the vicinity of laterals.


Introduction
Confluences of open channels and lateral drainage pipes are important elements in the hydraulic networks of man-made canals [1]. They are necessarily encountered at the end of urban sewage or agricultural subsurface drainage networks, where they release water [2]. Flow through lateral drainage pipes to open channels impacts channel capacity and disrupts the natural flow of water by creating turbulence in the vicinity of the inlet pipe, as well as influencing flow characteristics in the confluence hydrodynamics zone (CHZ) (Figure 1). Even though these confluences represent a critical component of drainage system geometry, they have received surprisingly little attention from researchers and engineers. This is mainly due to the relatively large number of parameters involved, along with the complex flow features occurring in the CHZ [3]. As a result, no comprehensive data set has been compiled that describes the 3D flow field within the CHZ [4]. The complexity of an open channel confluence stems from flow mixing, secondary circulation, post-confluence flow separation, contraction, and backwater effects. These effects, in turn, result in a large number of parameters required to adequately describe flow patterns and turbulent flow structures, due to flow merging [5]. In a CHZ, the flow topology is governed by parameters such as confluence geometry (angle, channel widths, bed discordance presence, and cross-sectional shape) and momentum ratio, as well as the ratio between the incoming flows, etc. [6].
Studying the combination of flows in open channel confluences with lateral drainage pipes has a direct application in drainage system open channel network design. One important parameter in such design is estimating the rise in flow depth, due to the lateral flow effect, and in a CHZ this effect can be significant [7]. An understanding of these flow conditions is necessary for the hydraulic design of the channel walls (height and length) that are required to contain these flows in the vicinity of laterals [8]. In this regard, the present study seeks to answer the following questions:


What would be the impact of submerged lateral drainage pipe discharges into rectangular open channels on the flow topology (including maximum water height) in a CHZ for different channel widths, flow rates, Froude numbers, confluence angles, inlet pipe diameters, and lateral flows to the main channel flow ratios?  Under what circumstances would the maximum increase in water surface elevation be confined to the desired values (e.g., 10%, or 20%)?
The first and only study to date on the hydraulic behavior of a pipe outlet into an open channel, allowing the combination of free-surface and pressurized flows, was conducted by Nedelec and Gay in 2008 [9]. This limited experimental study explored the hydraulic properties of subsurface drainage outlets into an open-air stream. This experimental study had several limitations, one of which was that the flow in channels during floods is super-critical, while the study was confined to sub-critical flow. Another limitation was that the angle of the inlet pipe was 90° with respect to the flume, while some flood control agencies recommend smaller pipe entry angles. Additionally, the results of the small-sized flume experiments could not be extrapolated to actual flood control channel dimensions.
Recognizing the knowledge gap, this study addresses the limitations identified in the literature review, utilizing experimental and numerical studies, and contributes to a better knowledge of flow topology in the CHZ of an open channel and a pipe. Additionally, this study can provide a benchmark experimental dataset for the validation of future numerical models. The experimental studies were performed for sub-critical and super-critical flow conditions with different junction angles. Results from these experiments were used to validate the numerical model (FLOW-3D) simulations. Finally, the validated numerical model was implemented to simulate actual flood control channel dimensions, in order to propose recommendations for the design of channels with lateral pipe inlets. These recommendations are based on conservative design curves, developed on the basis of the momentum principle, according to the relationship between the relative rise in The complexity of an open channel confluence stems from flow mixing, secondary circulation, post-confluence flow separation, contraction, and backwater effects. These effects, in turn, result in a large number of parameters required to adequately describe flow patterns and turbulent flow structures, due to flow merging [5]. In a CHZ, the flow topology is governed by parameters such as confluence geometry (angle, channel widths, bed discordance presence, and cross-sectional shape) and momentum ratio, as well as the ratio between the incoming flows, etc. [6].
Studying the combination of flows in open channel confluences with lateral drainage pipes has a direct application in drainage system open channel network design. One important parameter in such design is estimating the rise in flow depth, due to the lateral flow effect, and in a CHZ this effect can be significant [7]. An understanding of these flow conditions is necessary for the hydraulic design of the channel walls (height and length) that are required to contain these flows in the vicinity of laterals [8].
In this regard, the present study seeks to answer the following questions: • What would be the impact of submerged lateral drainage pipe discharges into rectangular open channels on the flow topology (including maximum water height) in a CHZ for different channel widths, flow rates, Froude numbers, confluence angles, inlet pipe diameters, and lateral flows to the main channel flow ratios? • Under what circumstances would the maximum increase in water surface elevation be confined to the desired values (e.g., 10%, or 20%)?
The first and only study to date on the hydraulic behavior of a pipe outlet into an open channel, allowing the combination of free-surface and pressurized flows, was conducted by Nedelec and Gay in 2008 [9]. This limited experimental study explored the hydraulic properties of subsurface drainage outlets into an open-air stream. This experimental study had several limitations, one of which was that the flow in channels during floods is super-critical, while the study was confined to sub-critical flow. Another limitation was that the angle of the inlet pipe was 90 • with respect to the flume, while some flood control agencies recommend smaller pipe entry angles. Additionally, the results of the small-sized flume experiments could not be extrapolated to actual flood control channel dimensions.
Recognizing the knowledge gap, this study addresses the limitations identified in the literature review, utilizing experimental and numerical studies, and contributes to a better knowledge of flow topology in the CHZ of an open channel and a pipe. Additionally, this study can provide a benchmark experimental dataset for the validation of future numerical models. The experimental studies were performed for sub-critical and super-critical flow conditions with different junction angles. Results from these experiments were used to validate the numerical model (FLOW-3D) simulations. Finally, the validated numerical model was implemented to simulate actual flood control channel dimensions, in order to propose recommendations for the design of channels with lateral pipe inlets.
These recommendations are based on conservative design curves, developed on the basis of the momentum principle, according to the relationship between the relative rise in channel water height, lateral discharge to channel flow ratio, and confluence angle. Figure 2 shows the flow chart of this study.
Hydrology 2020, 7, x FOR PEER REVIEW 3 of 51 channel water height, lateral discharge to channel flow ratio, and confluence angle. Figure 2 shows the flow chart of this study.

Experiments
The flume in the Fluid Dynamics Laboratory of the University of Nevada, Las Vegas (UNLV) is 14.63 m (48') long and 0.81 m (32") wide, with a height of 0.61 m (24"). In order to obtain super-critical flows in the flume, considering the availability of pumps at UNLV, the flume had to be narrowed to get reasonable flow velocities, water depths, and desired Froude numbers. Therefore, to provide the desired flow conditions, as well as to evaluate the effect of channel width on the main channel flow disturbances, at first, the flume was narrowed to 0.61 m (24") (wide channel) and then to 0.41 m (16") (narrow channel). A 5.1-cm (2") diameter inlet pipe was installed through the flume wall 7.62 m (25') downstream of the flume entrance to provide lateral flow to the channel. Flow through the inlet was provided through a separate variable-speed pump fed by water from the flume sump ( Figure 3).

Experiments
The flume in the Fluid Dynamics Laboratory of the University of Nevada, Las Vegas (UNLV) is 14.63 m (48 ) long and 0.81 m (32") wide, with a height of 0.61 m (24"). In order to obtain super-critical flows in the flume, considering the availability of pumps at UNLV, the flume had to be narrowed to get reasonable flow velocities, water depths, and desired Froude numbers. Therefore, to provide the desired flow conditions, as well as to evaluate the effect of channel width on the main channel flow disturbances, at first, the flume was narrowed to 0.61 m (24") (wide channel) and then to 0.41 m (16") (narrow channel). A 5.1-cm (2") diameter inlet pipe was installed through the flume wall 7.62 m (25 ) downstream of the flume entrance to provide lateral flow to the channel. Flow through the inlet was provided through a separate variable-speed pump fed by water from the flume sump ( Figure 3).   After setting up the flume, the experiments were conducted in both wide and narrow channel configurations; 30, 45, and 90° junction angles; with lateral flow to main channel flow ratios of 2.5% and 5% for the wide channel, and 2.5%, 5%, and 7.5% for the narrow channel; and different Froude numbers, from 0.72 (sub-critical) to 3.27 (super-critical). Figure 6 shows how lateral flows were fed to the main channel through the pipe. Table 1 summarizes the experimental scenarios.   After setting up the flume, the experiments were conducted in both wide and narrow channel configurations; 30, 45, and 90° junction angles; with lateral flow to main channel flow ratios of 2.5% and 5% for the wide channel, and 2.5%, 5%, and 7.5% for the narrow channel; and different Froude numbers, from 0.72 (sub-critical) to 3.27 (super-critical). Figure 6 shows how lateral flows were fed to the main channel through the pipe. Table 1 summarizes the experimental scenarios.   After setting up the flume, the experiments were conducted in both wide and narrow channel configurations; 30, 45, and 90° junction angles; with lateral flow to main channel flow ratios of 2.5% and 5% for the wide channel, and 2.5%, 5%, and 7.5% for the narrow channel; and different Froude numbers, from 0.72 (sub-critical) to 3.27 (super-critical). Figure 6 shows how lateral flows were fed to the main channel through the pipe. Table 1 summarizes the experimental scenarios. After setting up the flume, the experiments were conducted in both wide and narrow channel configurations; 30 • , 45 • , and 90 • junction angles; with lateral flow to main channel flow ratios of 2.5% and 5% for the wide channel, and 2.5%, 5%, and 7.5% for the narrow channel; and different Froude Hydrology 2020, 7, 57 5 of 53 numbers, from 0.72 (sub-critical) to 3.27 (super-critical). Figure 6 shows how lateral flows were fed to the main channel through the pipe. Table 1 summarizes the experimental scenarios.  In this study, the ultrasonic GREYLINE AVFM 5.0 sensor was mounted on the invert of the flume to measure the flow rate, flow velocity, and water surface elevation in the flume. The sensor mounts on the bottom of a channel with a stainless-steel mounting bracket and a single screw into the bottom of the channel. This device uses a submerged ultrasonic sensor to measure the velocity, flow rate, and water surface elevation. The sensor is completely sealed with no orifices or ports, and resists fouling, corrosion, and abrasion. The AVFM 5.0 can measure forward flow velocity up to 6 m/s (20 ft/s) and reverse flow up to 1.5 m/s (5 ft/s). The submerged velocity/level sensor measures flow in partially full and surcharged pipes and channels with pressure up to 68.95 kPa (10 psi) [10].

Measurement Equipment
In this study, the ultrasonic GREYLINE AVFM 5.0 sensor was mounted on the invert of the flume to measure the flow rate, flow velocity, and water surface elevation in the flume. The sensor mounts on the bottom of a channel with a stainless-steel mounting bracket and a single screw into the bottom of the channel. This device uses a submerged ultrasonic sensor to measure the velocity, flow rate, and water surface elevation. The sensor is completely sealed with no orifices or ports, and resists fouling, corrosion, and abrasion. The AVFM 5.0 can measure forward flow velocity up to 6 m/s (20 ft/s) and reverse flow up to 1.5 m/s (5 ft/s). The submerged velocity/level sensor measures flow in partially full and surcharged pipes and channels with pressure up to 68.95 kPa (10 psi) [10].

Numerical Modeling
Conventional 1D numerical modeling treats the junctions by means of internal boundary conditions for the shallow water equations by applying mass continuity combined with a simplified momentum or energy conservation at the junctions [1,4,8,[11][12][13][14]. It is readily apparent that such boundary conditions and simplified mathematical models are incapable of fully describing the highly turbulent 3D flow features present at a CHZ. The difficulty of adequately describing this flow with simplified mathematical models leads to the necessity of using a 3D computational fluid dynamics (CFD) code to visualize and understand the flow conditions in a CHZ. Numerical modeling is a less time-consuming and expensive way to obtain the various flow parameters needed for engineering design. Once a 3D model is validated, numerical experiments can be designed and conducted for the systematic examination of the standalone effects of different primary controls upon the flow structure development. A 3D CFD model solves the full 3D form of the Navier-Stokes equations, with (Reynolds-averaging (RANS)) or without temporal averaging (direct numerical simulation (DNS)); hence, it is capable of adequately reproducing the evolving primary and secondary flow patterns and turbulent flow structures.
Computer-based simulations are now the dominant tool for understanding and visualizing turbulent flow structures [15]. However, controlled flow visualization experiments are still necessary to direct, develop, and validate the numerical simulations now dominant in the field [16]. Considering the large number of scenarios and limitations in setting up different laboratory experiments, a numerical model (FLOW-3D) was used in this study to simulate the flow topology for different channel configurations, lateral inlet pipe sizes and angles, main channel flow, and lateral flow conditions. FLOW-3D is a general-purpose computational fluid dynamics (CFD) software. It employs specially developed numerical techniques to solve motion equations for fluids to obtain transient, three-dimensional solutions to multi-scale, multi-physics flow problems. An array of physical and numerical options allows users to apply FLOW-3D to a wide variety of fluid flow phenomena [17].
In this study, the simulations were performed on: Three different channel widths of 3.05 (10 ), 7.62 (25 ), 15.24 m (50 ); with three initial water heights in the channel of 1.22 (4 ), 2.13 (7 ), 3.05 m (10 ); three junction angles of 30, 45, 90 • ; and three inlet pipe diameters of 0.46 (1.5 ), 0.91 (3 ), 1.52 m (5 ); with lateral flow to main channel flow ratios of 2.5%, 5.0%, 10% (based on the flow rates in the 3.05 m (10 ) channel); and, three Froude numbers of 0.8, 1.2, 2.0. In all scenarios, the channel bed slope was 0.05%. Table 2 summarizes the simulation scenarios. It should be noted that the lateral flow proportions in the 7.62 (25 ) and 15.24 m (50 ) channels were based on the 3.05 m (10 ) channel, since 2.5%, 5.0%, and 10% of the flow rate in the 7.62 (25 ) and 15.24 m (50 ) channels would result in very high and uncommon velocities in the lateral inlet pipe. According to the Hydrologic Criteria and Drainage Design Manual (HCDDM) of the Clark County Regional Flood Control District (CCRFCD), the maximum flow velocity in a lateral drainage pipe into a flood control channel should not exceed 10.67 m/s (35 ft/s) [18]. For this reason, the lateral pipe flow proportions in all scenarios in this study were based on the 3.05 m (10 ) channel. However, pipe flow velocities up to 15.24 m/s (50 ft/s) were considered, regarding the fact that higher flow velocities may occur under critical conditions. The governing equations on fluid flow are mass continuity and momentum equations. The general mass continuity equation is: where V F is the fractional volume open to flow, ρ is the fluid density, R DIF is a turbulent diffusion term, and R SOR is a mass source. The velocity components (u, v, w) are in the coordinate directions (x, y, z) or (r, R SOR , z). A x is the fractional area open to flow in the x-direction, while A y and A z are similar area fractions for flow in the y and z directions, respectively. The coefficient R depends on the choice of the coordinate system. When Cartesian coordinates are used, R is set to unity and ξ is set to zero. The first term on the right side of Equation (1), is a turbulent diffusion term: where the coefficient v ρ is equal to S c µ ρ , in which µ is the coefficient of momentum diffusion (i.e., the viscosity); and S c is a constant, whose reciprocal is usually referred to as the turbulent Schmidt number. This type of mass diffusion only makes sense for turbulent mixing processes in fluids having a non-uniform density. The last term, R SOR , on the right side of Equation (1) is a density source term that can be used, for example, to model mass injection through porous obstacle surfaces. Compressible flow problems require the solution of the full density transport equation, as stated in Equation (1). For incompressible fluids, ρ is a constant, and Equation (1) reduces to the incompressibility condition: 8 of 53 For problems in which the propagation of acoustic pressure waves is important, but the fluid may otherwise be treated as incompressible, the density time derivative is approximated by: where c 2 is the square of the sound speed, and p is the pressure. This approximation is valid in the range δρ ρ < 0.1. With this approximation, the modified continuity equation then becomes:

Momentum Equations
The equations of motion for the fluid velocity components (u, v, w) in the three coordinate directions are Navier-Stokes equations with some additional terms: In these equations, (G x , G y , G z ) are body accelerations, ( f x , f y , f z ) are viscous accelerations, (b x , b y , b z ) are flow losses in porous media or across porous baffle plates, and the final terms account for the injection of mass at a source represented by a geometry component. The term U w = (u w , v w , w w ) in Equation (5) is the velocity of the source component, which will generally be non-zero for a mass source in a general moving objects model. The term U s = (u s , v s , w s ) is the velocity of the fluid at the surface of the source, relative to the source itself. It is computed in each control volume as: where dQ is the mass flow rate, ρQ is the fluid source density, dA is the area of the source surface in the cell, and n is the outward normal to the surface. When δ = 0.0, in Equation (5), the source is of the stagnation pressure type, whereas if δ = 1.0, the source is of the static pressure type. For a variable dynamic viscosity µ, the viscous accelerations are: where In the above expressions, the terms wsx, wsy, and wsz are wall shear stresses. If these terms are omitted, there is no wall shear stress because the remaining terms contain the fractional flow areas (A x , A y , A z ), which vanish at the walls. The wall stresses are modeled by assuming a zero tangential velocity on the portion of any area closed to flow. Mesh and moving obstacle boundaries are an exception because they can be assigned non-zero tangential velocities. In this case, the allowed boundary motion corresponds to a rigid body translation of the boundary parallel to its surface. For turbulent flows, a law-of-the-wall velocity profile is assumed near the wall, which modifies the wall shear stress magnitude.

Fluid Interfaces and Free-Surfaces
Fluid configurations are defined in terms of a volume of fluid (VOF) function, F (x, y, z, t). This function represents the volume of fluid #1 per unit volume and satisfies the equation: where The diffusion coefficient is defined as v F = c F µ ρ , where c F is a constant whose reciprocal is sometimes referred to as a turbulent Schmidt number. This diffusion term only makes sense for the turbulent mixing of two fluids whose distributions are defined by the F function. The term F SOR corresponds to the density source R SOR in Equation (1); F SOR is the time change rate of the fluid #1 volume fraction, associated with the mass source for fluid #1. The interpretation of F depends on the type of problem being solved. Incompressible problems must involve either a single fluid with a free-surface or two fluids and no free-surfaces. For a single fluid, F represents the volume fraction occupied by the fluid. Thus, fluid exists where F = 1, and void regions correspond to locations where F = 0. "Voids" are regions without fluid mass that have a uniform pressure assigned to them. Physically, they represent regions filled with a vapor or gas, whose density is insignificant with respect to the fluid density. Two-fluid problems may be composed of either two incompressible fluids, or one incompressible and one compressible fluid. F represents the volume fraction of the incompressible fluid component in either case, which is referred to as fluid #1. The complementary region, with the volume fraction 1 − F, represents fluid #2, and may have either a constant density, or its density may be computed from the compressible fluid equation-of-state. More details regarding FLOW-3D governing equations can be found in the FLOW-3D user's manual [17].

Water Height Measurements
Changes in flow proportions and junction angles resulted in measurable changes in the flow pattern. Figures 7-9 provide examples of water height measurements in the wide channel configuration, with and without lateral flow in a 90 • junction, which were used to calculate the maximum water height ratios in the flume. Other flume test results are available in Appendix A.
According to Figure 7, since the flow rates used for all scenarios were almost similar, in cases where the flow velocity was higher, the water height was lower. For instance, the water height for the Fr = 0.8 scenario was the highest, while for the Fr = 3.1, it was the lowest. Higher fluctuations were observed at the upstream of the flume in higher Froude numbers. This was due to the angled plate that was installed at the entrance of the flume (see Figure 4) to narrow the channel, which caused the formation of waves. However, the impacts diminished further downstream of the flume.  Figure 8 shows that as the lateral flow entered the channel, water heights in all scenarios were impacted either upstream or downstream of the pipe. The results showed that for Fr = 0.8, 1.0, and 1.2, the maximum water height occurred upstream of the pipe, while for Fr = 1.6, 2.0, 2.4, and 3.1, the maximum water height occurred downstream of the pipe. This was because the flow momentum was higher in the higher Froude numbers, which pushed the jet of water through the pipe further downstream. High peaks in the flume's entrance should be neglected as the highest water heights, since their formation was due to the angled wall installed upstream of the flume, which created local waves.  Figure 8 shows that as the lateral flow entered the channel, water heights in all scenarios were impacted either upstream or downstream of the pipe. The results showed that for Fr = 0.8, 1.0, and 1.2, the maximum water height occurred upstream of the pipe, while for Fr = 1.6, 2.0, 2.4, and 3.1, the maximum water height occurred downstream of the pipe. This was because the flow momentum was higher in the higher Froude numbers, which pushed the jet of water through the pipe further downstream. High peaks in the flume's entrance should be neglected as the highest water heights, since their formation was due to the angled wall installed upstream of the flume, which created local waves. Furthermore, comparing Figures 8 and 9, as the lateral flow proportion increased, water height peaks increased as well. According to Figure 9, the water height peak for Fr = 0.8, 1.0, 1.2, and 1.6 occurred upstream of the pipe, while for Fr = 2.0, 2.4, and 3.1, it occurred downstream of the pipe. Compared with the 2.5% lateral flow, more water height peaks were observed upstream of the pipe.

Water Height Ratios
Experimental test results were used to calculate the maximum water height ratios in the channel, by dividing the maximum water height in the channel with lateral flow, by the water height in the channel without lateral flow. In Figures 10-15, the maximum water height ratios for both wide and narrow channels are presented. Moreover, Tables 3-8 present maximum water height ratios, along with the flow rate and velocity in the flume and the pipe for each scenario. In the wide channel, the flow rate in the flume was maintained at the maximum value that could be achieved in the lab, but different flow rates were used in the narrow channel. For both channel configurations, to achieve certain flow velocities, and consequently, certain Froude numbers in the flume, the flume angle was adjusted.   Figure 8 shows that as the lateral flow entered the channel, water heights in all scenarios were impacted either upstream or downstream of the pipe. The results showed that for Fr = 0.8, 1.0, and 1.2, the maximum water height occurred upstream of the pipe, while for Fr = 1.6, 2.0, 2.4, and 3.1, the maximum water height occurred downstream of the pipe. This was because the flow momentum was higher in the higher Froude numbers, which pushed the jet of water through the pipe further downstream. High peaks in the flume's entrance should be neglected as the highest water heights, since their formation was due to the angled wall installed upstream of the flume, which created local waves.

90° Junction
Furthermore, comparing Figures 8 and 9, as the lateral flow proportion increased, water height peaks increased as well. According to Figure 9, the water height peak for Fr = 0.8, 1.0, 1.2, and 1.6 occurred upstream of the pipe, while for Fr = 2.0, 2.4, and 3.1, it occurred downstream of the pipe. Compared with the 2.5% lateral flow, more water height peaks were observed upstream of the pipe.

Water Height Ratios
Experimental test results were used to calculate the maximum water height ratios in the channel, by dividing the maximum water height in the channel with lateral flow, by the water height in the channel without lateral flow. In Figures 10-15, the maximum water height ratios for both wide and narrow channels are presented. Moreover, Tables 3-8 present maximum water height ratios, along with the flow rate and velocity in the flume and the pipe for each scenario. In the wide channel, the flow rate in the flume was maintained at the maximum value that could be achieved in the lab, but different flow rates were used in the narrow channel. For both channel configurations, to achieve certain flow velocities, and consequently, certain Froude numbers in the flume, the flume angle was adjusted.  According to Figure 10, the maximum water height ratios in the 5% lateral flow were higher than in the 2.5%. Moreover, in the 2.5% case, the maximum ratio occurred for Fr = 1.2, and for higher Froude numbers, the water height peak decreased. This was because the water jet through the pipe could reach the outer wall of the flume, and the water height peak was created by the combined effects of splashing, along with the incoming flow from the flume upstream. As the Froude number increased, the jet was deflected and could not reach the outer wall of the flume to create flow splashing. For Froude numbers smaller than 1.2, the water height in the channel was higher, compared with other Froude numbers, and the water height around the pipe was less impacted by the lateral flow through the submerged pipe. This resulted in lower maximum water height ratios in the channel. According to Figure 10, the maximum water height ratios in the 5% lateral flow were higher than in the 2.5%. Moreover, in the 2.5% case, the maximum ratio occurred for Fr = 1.2, and for higher Froude numbers, the water height peak decreased. This was because the water jet through the pipe could reach the outer wall of the flume, and the water height peak was created by the combined effects of splashing, along with the incoming flow from the flume upstream. As the Froude number increased, the jet was deflected and could not reach the outer wall of the flume to create flow splashing. For Froude numbers smaller than 1.2, the water height in the channel was higher, compared with other Froude numbers, and the water height around the pipe was less impacted by the lateral flow through the submerged pipe. This resulted in lower maximum water height ratios in the channel. Similarly, for the 5% lateral flow, in Fr = 1.2, the maximum water height ratio increased compared with the smaller Froude numbers. However, the maximum ratio occurred for Fr = 2. This can be justified by the fact that as the water height for similar flow rates is lower for higher Froude numbers, and as the 5% lateral flow jet is strong enough to create a blockage in the water path, this blockage could cause a hydraulic jump that increases the water height in the confluence.

• Junction
In the narrow channel ( Figure 11), similar to the wide channel ( Figure 10), as the lateral flow proportion increased, maximum water height ratios increased in the flume as well. In the 2.5% lateral flow case, the maximum water height ratio occurred for Fr = 0.95, but for 5% and 7.5% lateral flow proportions, the maximum occurred for Fr = 1.45. Similarly, for the 5% lateral flow, in Fr = 1.2, the maximum water height ratio increased compared with the smaller Froude numbers. However, the maximum ratio occurred for Fr = 2. This can be justified by the fact that as the water height for similar flow rates is lower for higher Froude numbers, and as the 5% lateral flow jet is strong enough to create a blockage in the water path, this blockage could cause a hydraulic jump that increases the water height in the confluence.
In the narrow channel ( Figure 11), similar to the wide channel ( Figure 10), as the lateral flow proportion increased, maximum water height ratios increased in the flume as well. In the 2.5% lateral flow case, the maximum water height ratio occurred for Fr = 0.95, but for 5% and 7.5% lateral flow proportions, the maximum occurred for Fr = 1.45.
numbers, and as the 5% lateral flow jet is strong enough to create a blockage in the water path, this blockage could cause a hydraulic jump that increases the water height in the confluence.
In the narrow channel ( Figure 11), similar to the wide channel ( Figure 10), as the lateral flow proportion increased, maximum water height ratios increased in the flume as well. In the 2.5% lateral flow case, the maximum water height ratio occurred for Fr = 0.95, but for 5% and 7.5% lateral flow proportions, the maximum occurred for Fr = 1.45.   In the 45 • junction, the results show that similar to the 90 • junction, higher lateral flows correspond to higher maximum water height ratios in the flume. According to Figure 12, in the wide channel, for the 2.5% lateral flow, the maximum water height ratio occurred for Fr = 2, while for the 5% lateral flow, the maximum occurred for Fr = 1.6. In the 45° junction, the results show that similar to the 90° junction, higher lateral flows correspond to higher maximum water height ratios in the flume. According to Figure 12, in the wide channel, for the 2.5% lateral flow, the maximum water height ratio occurred for Fr = 2, while for the 5% lateral flow, the maximum occurred for Fr = 1.6. In the narrow channel, and for the 45° junction ( Figure 13), when the water was sub-critical in the channel, the maximum water height ratios were almost the same for all scenarios. For supercritical flow conditions, the maximum values did not show a distinct pattern. However, the maximum water height ratios for 2.5%, 5%, and 7.5% lateral flow proportions occurred at Fr = 1.74, 2.01, and 2.01, respectively. In the narrow channel, and for the 45 • junction ( Figure 13), when the water was sub-critical in the channel, the maximum water height ratios were almost the same for all scenarios. For super-critical flow conditions, the maximum values did not show a distinct pattern. However, the maximum water height ratios for 2.5%, 5%, and 7.5% lateral flow proportions occurred at Fr = 1.74, 2.01, and 2.01, respectively.     In the 30 • junction, and the wide channel ( Figure 14), for both 2.5% and 5% lateral flow proportions, the maximum water height ratio occurred for Fr = 1.2. In the narrow channel ( Figure 15), the maximum water height ratios for 2.5%, 5%, and 7.5% lateral flow proportions occurred at Fr = 1.74, 2.01, and 1.74, respectively.  Based on the results for the 90°, 45°, and 30° junctions, both in the wide and narrow channel configurations, it can be concluded that, as the lateral pipe angle and lateral flow proportion increase, the maximum water height ratios in the channel also increase.  Based on the results for the 90°, 45°, and 30° junctions, both in the wide and narrow channel configurations, it can be concluded that, as the lateral pipe angle and lateral flow proportion increase, the maximum water height ratios in the channel also increase.  Based on the results for the 90 • , 45 • , and 30 • junctions, both in the wide and narrow channel configurations, it can be concluded that, as the lateral pipe angle and lateral flow proportion increase, the maximum water height ratios in the channel also increase.

Geometry
Geometry is constructed in FLOW-3D by assembling solid geometric objects to define the flow region for a simulation. The flow geometry is then embedded in the computational grid by the preprocessor using a technique called FAVOR TM , an acronym for fractional area/volume obstacle representation. This method computes open area fractions on the cell faces, along with the open volume fraction, and reconstructs the geometry based on these parameters. This approach offers a simple and accurate way to represent complex surfaces in the domain without requiring a body-fitted grid. Figure 16 shows a 3D view of the channel and its components, modeled in FLOW-3D. A roughness height (absolute roughness) of 0.001 m was assigned to all components, as the absolute roughness of ordinary concrete.

Meshing
FLOW-3D uses a structured mesh that may be defined in either a cartesian or cylindrical coordinate system. It is possible to define multiple mesh blocks (multi-block gridding) to create more efficient meshes when modeling complex flow domains. In this study, finer grids were used near the confluence, and coarser grids were used further away from the junction, where the influence of the confluence diminishes ( Figure 17). This kind of grid design has the potential to capture the near confluence flow behavior efficiently, while reducing computation cost. Table 9 presents the grid size for each mesh block. Figure 17. Meshing the channel, using four different grid sizes. Table 9. Grid sizes in each mesh block in the channel.

Meshing
FLOW-3D uses a structured mesh that may be defined in either a cartesian or cylindrical coordinate system. It is possible to define multiple mesh blocks (multi-block gridding) to create more efficient meshes when modeling complex flow domains. In this study, finer grids were used near the confluence, and coarser grids were used further away from the junction, where the influence of the confluence diminishes ( Figure 17). This kind of grid design has the potential to capture the near confluence flow behavior efficiently, while reducing computation cost. Table 9 presents the grid size for each mesh block. Table 9. Grid sizes in each mesh block in the channel.

Mesh Block
Grid Size (cm) 0.7 coordinate system. It is possible to define multiple mesh blocks (multi-block gridding) to create more efficient meshes when modeling complex flow domains. In this study, finer grids were used near the confluence, and coarser grids were used further away from the junction, where the influence of the confluence diminishes ( Figure 17). This kind of grid design has the potential to capture the near confluence flow behavior efficiently, while reducing computation cost. Table 9 presents the grid size for each mesh block. Figure 17. Meshing the channel, using four different grid sizes. Table 9. Grid sizes in each mesh block in the channel.

Favor
FAVOR TM is a very powerful method for incorporating geometry effects into the governing equations, but like all discrete methods, it is affected by the resolution of the computational grid. This is because the preprocessor generates area fractions for each cell face in the grid by determining which corners of the face are inside of a defined geometry. If all four corners of a cell face are inside the geometry, then the entire face is defined to be within the geometry. Similarly, if all corners lie outside, then the entire face is assumed to be outside the geometry. When some face corners are inside a geometry and some are outside, the intersection of the geometry with face edges is computed. Area fractions are then computed from these intersection points, assuming straight-line connections between the intersection points within the face. The straight-line assumption introduces a small error

Favor
FAVOR TM is a very powerful method for incorporating geometry effects into the governing equations, but like all discrete methods, it is affected by the resolution of the computational grid. This is because the preprocessor generates area fractions for each cell face in the grid by determining which corners of the face are inside of a defined geometry. If all four corners of a cell face are inside the geometry, then the entire face is defined to be within the geometry. Similarly, if all corners lie outside, then the entire face is assumed to be outside the geometry. When some face corners are inside a geometry and some are outside, the intersection of the geometry with face edges is computed. Area fractions are then computed from these intersection points, assuming straight-line connections between the intersection points within the face. The straight-line assumption introduces a small error in the fractional area when the geometry boundary is curved inside the cell. The approximation is consistent with the other assumptions in the equation development, and improves as the grid resolution is refined [17]. Figure 18 shows a 3D view of the favored channel. According to the figure, all components are properly captured in the mesh space. Moreover, Figure 19 shows how the pipe was modeled as a hole in a pipe complement attached to the inner wall of the flume, along with its favored view for 45 • and 30 • junctions.
Hydrology 2020, 7, x FOR PEER REVIEW 17 of 51 in the fractional area when the geometry boundary is curved inside the cell. The approximation is consistent with the other assumptions in the equation development, and improves as the grid resolution is refined [17]. Figure 18 shows a 3D view of the favored channel. According to the figure, all components are properly captured in the mesh space. Moreover, Figure 19 shows how the pipe was modeled as a hole in a pipe complement attached to the inner wall of the flume, along with its favored view for 45° and 30° junctions.

Boundary Conditions
The equations governing the motion of fluid flow can be categorized as initial-boundary value problems, meaning that the solution must be known at the initial time and boundaries to be found. This paradox is managed by assuming the solution at the boundaries (setting the boundary conditions) and initial time (initial conditions), and using these assumed values to solve the governing PDE inside of the domain. This implies that the solution to the problem is fundamentally defined by what is assumed at the boundaries. Therefore, it is quite important that the solution assumed at the boundary is a good approximation at this location. Figure 20 shows the boundary conditions on each mesh block face for the entire channel volume. The boundary conditions on the inner and outer walls, as well as and the channel bottom, were "wall." The boundary condition at the top of each mesh block was "symmetry." At the channel entrance, between the mesh blocks, and at the end of the channel, the boundary conditions were "volume flow rate," "grid overlay," and "outflow," respectively.

Boundary Conditions
The equations governing the motion of fluid flow can be categorized as initial-boundary value problems, meaning that the solution must be known at the initial time and boundaries to be found. This paradox is managed by assuming the solution at the boundaries (setting the boundary conditions) and initial time (initial conditions), and using these assumed values to solve the governing PDE inside of the domain. This implies that the solution to the problem is fundamentally defined by what is assumed at the boundaries. Therefore, it is quite important that the solution assumed at the boundary is a good approximation at this location. Figure 20 shows the boundary conditions on each mesh block face for the entire channel volume. The boundary conditions on the inner and outer walls, as well as and the channel bottom, were "wall". The boundary condition at the top of each mesh block was "symmetry". At the channel entrance, between the mesh blocks, and at the end of the channel, the boundary conditions were "volume flow rate", "grid overlay", and "outflow", respectively.
Hydrology 2020, 7, x FOR PEER REVIEW 18 of 51 The wall boundary condition applies the no-slip condition at the boundary, as well as a zerovelocity condition normal to the boundary. The symmetry boundary condition applies a zerogradient condition at the boundary, as well as a zero-velocity condition normal to the boundary. The volume flow rate boundary condition applies the specified flow rate at the boundary. The grid overlay boundary condition applies the solution from the restart source simulation, as a boundary condition in a restart simulation. Finally, the outflow boundary condition uses a Sommerfeld radiation condition to dynamically estimate the conditions at the boundary.

Initial Conditions
Initial conditions define the starting conditions in a simulation. As mentioned in the mesh boundary conditions, the initial state of the solution for transient fluid flow problems must be known in order to find a solution. Further, in a manner similar to what is done with boundary conditions, the initial conditions are assumed, approximating the true state at time t = 0. The accuracy of the initial conditions is somewhat less important than that of the boundary conditions because their effect is reduced as time progresses. Because of this effect, it is common to define the initial fluid geometry with reasonable accuracy, but to assume that the pressure and velocity fields are uniform. This The wall boundary condition applies the no-slip condition at the boundary, as well as a zero-velocity condition normal to the boundary. The symmetry boundary condition applies a zero-gradient condition at the boundary, as well as a zero-velocity condition normal to the boundary. The volume flow rate boundary condition applies the specified flow rate at the boundary. The grid overlay boundary condition applies the solution from the restart source simulation, as a boundary condition in a restart simulation. Finally, the outflow boundary condition uses a Sommerfeld radiation condition to dynamically estimate the conditions at the boundary.

Initial Conditions
Initial conditions define the starting conditions in a simulation. As mentioned in the mesh boundary conditions, the initial state of the solution for transient fluid flow problems must be known in order to find a solution. Further, in a manner similar to what is done with boundary conditions, the initial conditions are assumed, approximating the true state at time t = 0. The accuracy of the initial conditions is somewhat less important than that of the boundary conditions because their effect is reduced as time progresses. Because of this effect, it is common to define the initial fluid geometry with reasonable accuracy, but to assume that the pressure and velocity fields are uniform. This typically produces good results within a relatively short time, where the solution is affected by physical initial conditions.
In this study, the initial conditions were: An empty channel at time zero; a certain water height and velocity upstream of the channel, depending on the desired Froude number; and a certain flow rate through the pipe, depending on the desired flow proportion relative to the flow rate in the channel.

Other Considerations in Setting up the Model
In setting up the model in FLOW-3D, the reference pressure was set to 1 atm; the reference temperature was 20 • C; the fluid used was water at 20 • C; and the finish time for most of the simulations was 20 s unless longer simulation times were required for the flow in the channel to reach steady-state conditions. Moreover, the renormalized group (RNG) model was used as the turbulence model.

Convergence and Validation
Fluid flow problems are highly nonlinear. Only through imposing restrictive conditions can the governing Navier-Stokes equations be solved analytically. As a result, CFD solutions must be calculated iteratively. By their nature, iterative solution methods require a convergence criterion that is used to decide when the iterations can be terminated. In FLOW-3D, all relaxation and convergence criteria are selected by the program itself. Further, all selections are adjusted dynamically by the program to follow the development of the solution. One such case would be the use of a very large convergence criteria and no over-relaxation, as a way to reach steady-state conditions with less CPU time. Figures 21 and 22 show that by the end of the simulation time (for one of the simulated scenarios, as an example), the solution converged and the time-step size was within the stability limits, respectively.
Hydrology 2020, 7, x FOR PEER REVIEW 19 of 51 scenarios, as an example), the solution converged and the time-step size was within the stability limits, respectively.    Using CFD allows researchers to understand flow behavior and quantify important flow parameters, given that the CFD solver has been properly validated using observations. Once validated, the 3D numerical model can be used to investigate the flow characteristics at the junction.
The results of three experiments with different Froude numbers (0.8, 2.0, 3.1) were used to validate the FLOW-3D model. These three cases were selected to cover both the sub-critical and super-critical conditions in the flume. The flume wall was made of plexiglass sheets attached at certain lengths (1.03 m), making a total of 13 columns. The simulation and experimental results of water heights at each column are presented in Figures 23-25.
Using the Nash-Sutcliffe efficiency function (NS) (Equation (16)), the goodness of fit between the observed and simulated values were evaluated and are presented in Table 10. In Equation (16), "H" is a variable, such as the flow surface elevation; "m" and "s" stand for measured and simulated; the bar indicates the average; and "i" is the ith measured or simulated value. The NS function has a range of −∞ to 1. NS = 1 corresponds to a perfect match of simulated values to the observed data. The values between 0 and 1 indicate that the simulated and observed values are in close agreement, whereas values less than 0 show that the model has no predictive power [19,20].
Good agreement was found between predicted water surface levels and the experimental measurements across the stream, throughout the flume, as shown in Using CFD allows researchers to understand flow behavior and quantify important flow parameters, given that the CFD solver has been properly validated using observations. Once validated, the 3D numerical model can be used to investigate the flow characteristics at the junction.
The results of three experiments with different Froude numbers (0.8, 2.0, 3.1) were used to validate the FLOW-3D model. These three cases were selected to cover both the sub-critical and super-critical conditions in the flume. The flume wall was made of plexiglass sheets attached at certain lengths (1.03 m), making a total of 13 columns. The simulation and experimental results of water heights at each column are presented in Figures 23-25.
Using the Nash-Sutcliffe efficiency function (NS) (Equation (16)), the goodness of fit between the observed and simulated values were evaluated and are presented in Table 10. In Equation (16), "H" is a variable, such as the flow surface elevation; "m" and "s" stand for measured and simulated; the bar indicates the average; and "i" is the ith measured or simulated value. The NS function has a range of −∞ to 1. NS = 1 corresponds to a perfect match of simulated values to the observed data. The values between 0 and 1 indicate that the simulated and observed values are in close agreement, whereas values less than 0 show that the model has no predictive power [19,20].
Good agreement was found between predicted water surface levels and the experimental measurements across the stream, throughout the flume, as shown in Figures 23-25. The high NS values obtained indicated that the model was able to simulate flow profiles in an open channel with lateral flow, with high accuracy. The small differences between simulation and experimental results can be attributed to errors in reading water heights in the lab, due to flow turbulence and flume waves. lateral flow, with high accuracy. The small differences between simulation and experimental results can be attributed to errors in reading water heights in the lab, due to flow turbulence and flume waves.  can be attributed to errors in reading water heights in the lab, due to flow turbulence and flume waves.

Surface Water Elevations and Conservative Design Curves
The incoming tributary flow caused a rise of water surface elevation, and for sub-critical flow conditions, resulted in a back-water effect towards the upstream ( Figure 26). Additionally, larger junction angles allowed the tributary flow to penetrate further into the confluence, which in some cases resulted in flow splashing on the outer wall of the channel (Figure 27). After the mixing between the tributary and mainstream, the confluence cross-sectional flow field gradually recovered to that of a normal open channel, at the downstream of the confluence.

Surface Water Elevations and Conservative Design Curves
The incoming tributary flow caused a rise of water surface elevation, and for sub-critical flow conditions, resulted in a back-water effect towards the upstream ( Figure 26). Additionally, larger junction angles allowed the tributary flow to penetrate further into the confluence, which in some cases resulted in flow splashing on the outer wall of the channel (Figure 27). After the mixing between the tributary and mainstream, the confluence cross-sectional flow field gradually recovered to that of a normal open channel, at the downstream of the confluence.

Surface Water Elevations and Conservative Design Curves
The incoming tributary flow caused a rise of water surface elevation, and for sub-critical flow conditions, resulted in a back-water effect towards the upstream (Figure 26). Additionally, larger junction angles allowed the tributary flow to penetrate further into the confluence, which in some cases resulted in flow splashing on the outer wall of the channel (Figure 27). After the mixing between the tributary and mainstream, the confluence cross-sectional flow field gradually recovered to that of a normal open channel, at the downstream of the confluence.  Among the 648 scenarios simulated, in 42 cases, as a result of the interaction between channel flow and lateral flow through the pipe, waves were formed in the channel, and in other cases the impact was local. Waves were not observed in 15.24 m (50 )-wide channels, since the channel flow momentum could overcome the pipe flow momentum, and the jet could not reach the outer wall. Among 42 cases that showed the formation of waves in 3.05 (10 ) and 7.62 m (25 ) channels, 30 cases occurred in 3.05 m (10 )-channels and 12 cases in 7.62 m (25 )-channels. Moreover, in 27 of the 42 cases, the flow was sub-critical, and in the other 15, it was super-critical. These findings indicate that narrower channel and pipe confluences, with lower Froude numbers, are more prone to the formation of waves in the channels.
The simulation results and conservative design curves for the three channel widths, and based on the Froude numbers of the flow in the channels, are presented in Figures 28-33. Detailed FLOW-3D results are available in Appendix B. In the following figures, W is the channel width; H 0 is the initial water height in the channel without lateral flow; H is the maximum water height in the channel with lateral flow; Q is the flow rate in the channel; V is the flow velocity in the channel; q is the flow rate in the pipe; v is the flow velocity in the pipe; and θ is the junction angle. As the two flows merged in the confluence, the flow rate and velocity in the channel and pipe were the two factors impacting the maximum water height in the channel. Accordingly, the momentum ratios of the channel and the pipe ( QV + qv(cos θ) qv(sin θ) ), could be related to the H/H 0 , as shown in Figure 28, Figure 30, and Figure 32. Among the 648 scenarios simulated, in 42 cases, as a result of the interaction between channel flow and lateral flow through the pipe, waves were formed in the channel, and in other cases the impact was local. Waves were not observed in 15.24 m (50')-wide channels, since the channel flow momentum could overcome the pipe flow momentum, and the jet could not reach the outer wall. Among 42 cases that showed the formation of waves in 3.05 (10') and 7.62 m (25') channels, 30 cases occurred in 3.05 m (10')-channels and 12 cases in 7.62 m (25')-channels. Moreover, in 27 of the 42 cases, the flow was sub-critical, and in the other 15, it was super-critical. These findings indicate that narrower channel and pipe confluences, with lower Froude numbers, are more prone to the formation of waves in the channels.
The simulation results and conservative design curves for the three channel widths, and based on the Froude numbers of the flow in the channels, are presented in Figures 28-33. Detailed FLOW-3D results are available in Appendix B. In the following figures, W is the channel width; H0 is the initial water height in the channel without lateral flow; H is the maximum water height in the channel with lateral flow; Q is the flow rate in the channel; V is the flow velocity in the channel; q is the flow rate in the pipe; v is the flow velocity in the pipe; and θ is the junction angle. As the two flows merged in the confluence, the flow rate and velocity in the channel and pipe were the two factors impacting the maximum water height in the channel. Accordingly, the momentum ratios of the channel and the pipe (                The following design guidelines are conservative. The guidelines provide conservative designs for extreme flow conditions in the channel and the pipe. Accordingly, the designer needs to first establish the expected extreme flow conditions in the channel and the pipe. Then, using the guidelines with the expected extreme flow conditions in the channel and the pipe as functional design criteria, the design will accommodate these extreme flow conditions. The uncertainty resides in the estimation of the expected extreme flow conditions in the channel and the pipe and not in the design guidelines. Moreover, it should be mentioned that the presence of impurities such as sediments in water will impact the unit weight (density) of the water and the roughness of the materials impacts the momentum of the water. The design guidelines are based on the relative momentums of the flows in the channel and the pipe. Since the momentum is a function of the unit weight (density), the guidelines include the effect of the presence of impurities in water and the roughness of the materials.

Conclusions and Recommendations
This study utilized a hybrid experimental-numerical approach to investigate the flow topology in the confluence of rectangular open channels with submerged lateral drainage pipe discharges to develop design guidelines for flood control channel height. The scenarios covered a range of channel widths, flow rates, Froude numbers, junction angles, inlet pipe diameters, and lateral flows to main channel flow ratios. Clear water at 20 °C with no impurities was used both in the experiments and numerical modeling.
The experimental results showed that as the junction angle and the lateral flow proportions increased, the maximum water height ratios in the channel also increased. Moreover, the numerical ranges. The following design guidelines are conservative. The guidelines provide conservative designs for extreme flow conditions in the channel and the pipe. Accordingly, the designer needs to first establish the expected extreme flow conditions in the channel and the pipe. Then, using the guidelines with the expected extreme flow conditions in the channel and the pipe as functional design criteria, the design will accommodate these extreme flow conditions. The uncertainty resides in the estimation of the expected extreme flow conditions in the channel and the pipe and not in the design guidelines.
Moreover, it should be mentioned that the presence of impurities such as sediments in water will impact the unit weight (density) of the water and the roughness of the materials impacts the momentum of the water. The design guidelines are based on the relative momentums of the flows in the channel and the pipe. Since the momentum is a function of the unit weight (density), the guidelines include the effect of the presence of impurities in water and the roughness of the materials.

Conclusions and Recommendations
This study utilized a hybrid experimental-numerical approach to investigate the flow topology in the confluence of rectangular open channels with submerged lateral drainage pipe discharges to develop design guidelines for flood control channel height. The scenarios covered a range of channel widths, flow rates, Froude numbers, junction angles, inlet pipe diameters, and lateral flows to main channel flow ratios. Clear water at 20 • C with no impurities was used both in the experiments and numerical modeling.
The experimental results showed that as the junction angle and the lateral flow proportions increased, the maximum water height ratios in the channel also increased. Moreover, the numerical modeling results showed that as the ratio of channel flow momentum to pipe flow momentum ( QV + qv(cos θ) qv(sin θ) ) increased, H/H 0 in the channel decreased. Additionally, narrower channel and pipe confluences, with lower Froude numbers, were found to be more prone to the formation of waves in the channel.
The results indicated that the momentum ratios of 500 and 1000 could be considered as critical numbers. Among almost all scenarios, in the cases where the value of QV + qv(cos θ) qv(sin θ) was greater than 500, the maximum increase in channel water height was 20%. Similarly, in cases where QV + qv(cos θ) qv(sin θ) was greater than 1000, the maximum water height increase in the channel was 10%. Additionally, in cases where the QV + qv(cos θ) qv(sin θ) values were greater than 500, the impacts on water height changes in the channel were local (at the pipe), while channel wave impacts were observed in low momentum ratios. Furthermore, for the channel's sub-critical and super-critical conditions, the waves propagated upstream and downstream of the channel, respectively.
Based on the findings of this study, if junctions are to be constructed at angles larger than 30 • , it is recommended that the momentum ratio between the channel and the pipe flow ( QV + qv(cos θ) qv(sin θ) ) be greater than 500, to limit the impacts of the lateral drainage pipe flow to the main channel flow around the pipe (local); it is also recommended to limit the maximum increase in the water height to 20%.
The design guidelines presented in this study are conservative, and they provide conservative designs for extreme flow conditions in the channel and the pipe. The designer has to first establish the expected extreme flow conditions in the channel and the pipe. Then, using the guidelines with the expected extreme flow conditions in the channel and the pipe, the design will accommodate these extreme flow conditions. The uncertainty resides in the estimation of the expected extreme flow conditions in the channel and the pipe and not in the design guidelines.
For future extension of this research, it is recommended to study: (1) The impact of flow through submerged lateral pipes into other channel and pipe confluence configurations (e.g., trapezoidal and triangular channels), different pipe sizes, and junction angles; (2) Froude numbers and lateral flow proportions outside the range covered in this study; (3) water with debris and impurities, under different environmental conditions (e.g., in the presence of wind); (4) the penetration length of the jet through the pipe into the channel, and its impact on the flow characteristics and flow topology in the channel.

Appendix B
In Tables A1-A9, Figure A45. Water height ratios in a 90 • junction for Fr = 3.27 in the narrow channel.