Next Article in Journal
Variations in Canopy Cover and Its Relationship with Canopy Water and Temperature in the Miombo Woodland Based on Satellite Data
Previous Article in Journal
Advances in Flood Early Warning: Ensemble Forecast, Information Dissemination and Decision-Support Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Mohammad Nazari-Sharabian
1,
Moses Karakouzian
1,* and
Donald Hayes
2
1
Department of Civil and Environmental Engineering and Construction, University of Nevada, Las Vegas, NV 89154, USA
2
U.S. Army Engineer Research and Development Center, Vicksburg, MS 39180, USA
*
Author to whom correspondence should be addressed.
Hydrology 2020, 7(3), 57; https://doi.org/10.3390/hydrology7030057
Submission received: 5 July 2020 / Revised: 11 August 2020 / Accepted: 12 August 2020 / Published: 15 August 2020

Abstract

:
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.

1. 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 channel water height, lateral discharge to channel flow ratio, and confluence angle. Figure 2 shows the flow chart of this study.

2. Materials and Methods

2.1. 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).
Moreover, three inlet connections at 30, 45, and 90° were constructed to feed lateral flow to the channel at desired angles. Figure 4 and Figure 5 show the wide channel and narrow channel configurations, respectively.
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.

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].

2.2. 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.

2.2.1. Governing Equations

The governing equations on fluid flow are mass continuity and momentum equations. The general mass continuity equation is:
V F ρ t + x ( ρ u A x ) + R y ( ρ v A y ) + z ( ρ w A z ) + ξ ρ u A x x = R D I F + R S O R
where V F is the fractional volume open to flow, ρ is the fluid density, R D I F is a turbulent diffusion term, and R S O R is a mass source. The velocity components ( u , v , w ) are in the coordinate directions ( x , y , z ) or ( r ,   R S O R , 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:
R D I F = x ( v p A x ρ x ) + R y ( v ρ A y R ρ y ) + z ( v ρ A z ρ z ) + ξ v ρ A x x ρ x
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 S O R , 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:
x ( u A x ) + R y ( v A y ) + z ( w A z ) + ξ u A x x = R S O R ρ
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:
ρ t 1 c 2 p t
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:
V F ρ c 2 p t + u A x x + R v A y y + w A z z + ξ u A x x = R S O R ρ

2.2.2. 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:
u t + 1 V F { u A x u x + v A y R u y + w A z u z } ξ A y v 2 x V F          = 1 ρ p x + G x + f x b x R S O R ρ V F ( u u w δ u s ) v t + 1 V F { u A x v x + v A y R v y + w A z v z } + ξ A y u v x V F          = 1 ρ ( R p y ) + G y + f y b y R S O R ρ V F ( v v w δ v s ) w t + 1 V F { u A x w x + v A y R w y + w A z w z }          = 1 ρ p z + G z + f z b z R S O R ρ V F ( w w w δ w s )
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:
U s = d Q ρ Q d A n
where d Q is the mass flow rate, ρ Q is the fluid source density, d A 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:
ρ V F f x = w s x { x ( A x τ x x ) + R y ( A y τ x y ) + z ( A z τ x z ) + ξ x ( A x τ x x A y τ y y ) } ρ V F f y = w s y { x ( A x τ x y ) + R y ( A y τ y y ) + z ( A z τ y z ) + ξ x ( A x + A y ) τ x y } ρ V F f z = w s z { x ( A x τ x z ) + R y ( A y τ y z ) + z ( A z τ z z ) + ξ x ( A x τ x z ) }
where
τ x x = 2 μ { u x 1 3 ( u x + R v y + w z + ξ u x ) }
τ y y = 2 μ { R v y + ξ u x 1 3 ( u x + R v y + w z + ξ u x ) }
τ z z = 2 μ { w z 1 3 ( u x + R v y + w z + ξ u x ) }
τ x y = μ { v x + R u y ξ v x }
τ x z = μ { u z + w x }
τ y z = μ { v z + R w y }
In the above expressions, the terms w s x , w s y , and w s z 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.

2.2.3. 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:
F t + 1 V F [ x ( F A x u ) + R y ( F A y v ) + z ( F A z w ) + ξ F A x u x ] = F D I F + F S O R
where
F D I F = 1 V F { x ( v F A x F x ) + R x ( v F A y R F y ) + z ( v F A z F z ) + ξ v F A x F x }
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 S O R corresponds to the density source R S O R in Equation (1); F S O R 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].

3. Results and Discussion

3.1. Flume Test Results

3.1.1. Water Height Measurements

Changes in flow proportions and junction angles resulted in measurable changes in the flow pattern. Figure 7, Figure 8 and Figure 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.
Furthermore, comparing Figure 8 and Figure 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.

3.1.2. 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 Figure 10, Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15, the maximum water height ratios for both wide and narrow channels are presented. Moreover, Table 3, Table 4, Table 5, Table 6, Table 7 and Table 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.

90° Junction

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.
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.

45° Junction

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 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.

30° Junction

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.

3.2. Numerical Model Setup

3.2.1. 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 FAVORTM, 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.

3.2.2. 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.

3.2.3. Favor

FAVORTM 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.

3.2.4. 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.
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.

3.2.5. 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.

3.2.6. 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.

3.2.7. 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. Figure 21 and Figure 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.
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 Figure 23, Figure 24 and Figure 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 Figure 23, Figure 24 and Figure 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.
NS = 1   i | H m   H s | i 2 i | H m , i H ¯ m | i 2

3.3. Simulation Results

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 Figure 28, Figure 29, Figure 30, Figure 31, Figure 32 and Figure 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 ( QV   +   qv ( cos θ ) qv ( sin θ ) ), could be related to the H/H0, as shown in Figure 28, Figure 30, and Figure 32.
In all three charts, in Figure 28, Figure 30, and Figure 32, the results show that the highest H/H0 values were observed for W = 3.05 m (10′), and as the channel width increased to 7.62 (25′), and then to 15.24 m (50′), the points, representing the maximum water height ratios in the channel (H/H0), accumulated at the lower parts of the charts. The results indicate that QV   +   qv ( cos θ ) qv ( sin θ ) values of 500 and 1000 could be considered as critical numbers. Among almost all scenarios, in cases where the value of QV   +   qv ( cos θ ) qv ( sin θ ) was greater than 500, the maximum increase in water height was 20%. Similarly, in cases where QV   +   qv ( cos θ ) qv ( sin θ ) was greater than 1000, the maximum increase in water height was 10%. Additionally, in cases where the QV   +   qv ( cos θ ) qv ( sin θ ) values were greater than 500, impacts on the channel flow were local, while wave impacts were observed in channels in smaller QV   +   qv ( cos θ ) qv ( sin θ ) 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.

4. 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/H0 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.

Author Contributions

Conceptualization, M.N.-S., M.K., and D.H.; methodology, M.N.-S., M.K., and D.H.; software, M.N.-S.; validation, M.N.-S.; resources, M.N.-S. and M.K.; writing—original draft preparation, M.N.-S.; writing—review and editing, M.K. and D.H.; supervision, M.K. and D.H.; funding acquisition, M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Clark County Regional Flood Control District (CCRFCD), and the publication fees for this article were supported by the UNLV University Libraries Open Article Fund.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Water height ratios in a 30° junction for Fr = 0.8 in the wide channel.
Figure A1. Water height ratios in a 30° junction for Fr = 0.8 in the wide channel.
Hydrology 07 00057 g0a1
Figure A2. Water height ratios in a 45° junction for Fr = 0.8 in the wide channel.
Figure A2. Water height ratios in a 45° junction for Fr = 0.8 in the wide channel.
Hydrology 07 00057 g0a2
Figure A3. Water height ratios in a 90° junction for Fr = 0.8 in the wide channel.
Figure A3. Water height ratios in a 90° junction for Fr = 0.8 in the wide channel.
Hydrology 07 00057 g0a3
Figure A4. Water height ratios in a 30° junction for Fr = 1.0 in the wide channel.
Figure A4. Water height ratios in a 30° junction for Fr = 1.0 in the wide channel.
Hydrology 07 00057 g0a4
Figure A5. Water height ratios in a 45° junction for Fr = 1.0 in the wide channel.
Figure A5. Water height ratios in a 45° junction for Fr = 1.0 in the wide channel.
Hydrology 07 00057 g0a5
Figure A6. Water height ratios in a 90° junction for Fr = 1.0 in the wide channel.
Figure A6. Water height ratios in a 90° junction for Fr = 1.0 in the wide channel.
Hydrology 07 00057 g0a6
Figure A7. Water height ratios in a 30° junction for Fr = 1.2 in the wide channel.
Figure A7. Water height ratios in a 30° junction for Fr = 1.2 in the wide channel.
Hydrology 07 00057 g0a7
Figure A8. Water height ratios in a 45° junction for Fr = 1.2 in the wide channel.
Figure A8. Water height ratios in a 45° junction for Fr = 1.2 in the wide channel.
Hydrology 07 00057 g0a8
Figure A9. Water height ratios in a 90° junction for Fr = 1.2 in the wide channel.
Figure A9. Water height ratios in a 90° junction for Fr = 1.2 in the wide channel.
Hydrology 07 00057 g0a9
Figure A10. Water height ratios in a 30° junction for Fr = 1.6 in the wide channel.
Figure A10. Water height ratios in a 30° junction for Fr = 1.6 in the wide channel.
Hydrology 07 00057 g0a10
Figure A11. Water height ratios in a 45° junction for Fr = 1.6 in the wide channel.
Figure A11. Water height ratios in a 45° junction for Fr = 1.6 in the wide channel.
Hydrology 07 00057 g0a11
Figure A12. Water height ratios in a 90° junction for Fr = 1.6 in the wide channel.
Figure A12. Water height ratios in a 90° junction for Fr = 1.6 in the wide channel.
Hydrology 07 00057 g0a12
Figure A13. Water height ratios in a 30° junction for Fr = 2.0 in the wide channel.
Figure A13. Water height ratios in a 30° junction for Fr = 2.0 in the wide channel.
Hydrology 07 00057 g0a13
Figure A14. Water height ratios in a 45° junction for Fr = 2.0 in the wide channel.
Figure A14. Water height ratios in a 45° junction for Fr = 2.0 in the wide channel.
Hydrology 07 00057 g0a14
Figure A15. Water height ratios in a 90° junction for Fr = 2.0 in the wide channel.
Figure A15. Water height ratios in a 90° junction for Fr = 2.0 in the wide channel.
Hydrology 07 00057 g0a15
Figure A16. Water height ratios in a 30° junction for Fr = 2.4 in the wide channel.
Figure A16. Water height ratios in a 30° junction for Fr = 2.4 in the wide channel.
Hydrology 07 00057 g0a16
Figure A17. Water height ratios in a 45° junction for Fr = 2.4 in the wide channel.
Figure A17. Water height ratios in a 45° junction for Fr = 2.4 in the wide channel.
Hydrology 07 00057 g0a17
Figure A18. Water height ratios in a 90° junction for Fr = 2.4 in the wide channel.
Figure A18. Water height ratios in a 90° junction for Fr = 2.4 in the wide channel.
Hydrology 07 00057 g0a18
Figure A19. Water height ratios in a 30° junction for Fr = 3.1 in the wide channel.
Figure A19. Water height ratios in a 30° junction for Fr = 3.1 in the wide channel.
Hydrology 07 00057 g0a19
Figure A20. Water height ratios in a 45° junction for Fr = 3.1 in the wide channel.
Figure A20. Water height ratios in a 45° junction for Fr = 3.1 in the wide channel.
Hydrology 07 00057 g0a20
Figure A21. Water height ratios in a 90° junction for Fr = 3.1 in the wide channel.
Figure A21. Water height ratios in a 90° junction for Fr = 3.1 in the wide channel.
Hydrology 07 00057 g0a21
Figure A22. Water height ratios in a 30° junction for Fr = 0.72 in the narrow channel.
Figure A22. Water height ratios in a 30° junction for Fr = 0.72 in the narrow channel.
Hydrology 07 00057 g0a22
Figure A23. Water height ratios in a 45° junction for Fr = 0.72 in the narrow channel.
Figure A23. Water height ratios in a 45° junction for Fr = 0.72 in the narrow channel.
Hydrology 07 00057 g0a23
Figure A24. Water height ratios in a 90° junction for Fr = 0.72 in the narrow channel.
Figure A24. Water height ratios in a 90° junction for Fr = 0.72 in the narrow channel.
Hydrology 07 00057 g0a24
Figure A25. Water height ratios in a 30° junction for Fr = 0.95 in the narrow channel.
Figure A25. Water height ratios in a 30° junction for Fr = 0.95 in the narrow channel.
Hydrology 07 00057 g0a25
Figure A26. Water height ratios in a 45° junction for Fr = 0.95 in the narrow channel.
Figure A26. Water height ratios in a 45° junction for Fr = 0.95 in the narrow channel.
Hydrology 07 00057 g0a26
Figure A27. Water height ratios in a 90° junction for Fr = 0.95 in the narrow channel.
Figure A27. Water height ratios in a 90° junction for Fr = 0.95 in the narrow channel.
Hydrology 07 00057 g0a27
Figure A28. Water height ratios in a 30° junction for Fr = 1.45 in the narrow channel.
Figure A28. Water height ratios in a 30° junction for Fr = 1.45 in the narrow channel.
Hydrology 07 00057 g0a28
Figure A29. Water height ratios in a 45° junction for Fr = 1.45 in the narrow channel.
Figure A29. Water height ratios in a 45° junction for Fr = 1.45 in the narrow channel.
Hydrology 07 00057 g0a29
Figure A30. Water height ratios in a 90° junction for Fr = 1.45 in the narrow channel.
Figure A30. Water height ratios in a 90° junction for Fr = 1.45 in the narrow channel.
Hydrology 07 00057 g0a30
Figure A31. Water height ratios in a 30° junction for Fr = 1.74 in the narrow channel.
Figure A31. Water height ratios in a 30° junction for Fr = 1.74 in the narrow channel.
Hydrology 07 00057 g0a31
Figure A32. Water height ratios in a 45° junction for Fr = 1.74 in the narrow channel.
Figure A32. Water height ratios in a 45° junction for Fr = 1.74 in the narrow channel.
Hydrology 07 00057 g0a32
Figure A33. Water height ratios in a 90° junction for Fr = 1.74 in the narrow channel.
Figure A33. Water height ratios in a 90° junction for Fr = 1.74 in the narrow channel.
Hydrology 07 00057 g0a33
Figure A34. Water height ratios in a 30° junction for Fr = 2.01 in the narrow channel.
Figure A34. Water height ratios in a 30° junction for Fr = 2.01 in the narrow channel.
Hydrology 07 00057 g0a34
Figure A35. Water height ratios in a 45° junction for Fr = 2.01 in the narrow channel.
Figure A35. Water height ratios in a 45° junction for Fr = 2.01 in the narrow channel.
Hydrology 07 00057 g0a35
Figure A36. Water height ratios in a 90° junction for Fr = 2.01 in the narrow channel.
Figure A36. Water height ratios in a 90° junction for Fr = 2.01 in the narrow channel.
Hydrology 07 00057 g0a36
Figure A37. Water height ratios in a 30° junction for Fr = 2.32 in the narrow channel.
Figure A37. Water height ratios in a 30° junction for Fr = 2.32 in the narrow channel.
Hydrology 07 00057 g0a37
Figure A38. Water height ratios in a 45° junction for Fr = 2.32 in the narrow channel.
Figure A38. Water height ratios in a 45° junction for Fr = 2.32 in the narrow channel.
Hydrology 07 00057 g0a38
Figure A39. Water height ratios in a 90° junction for Fr = 2.32 in the narrow channel.
Figure A39. Water height ratios in a 90° junction for Fr = 2.32 in the narrow channel.
Hydrology 07 00057 g0a39
Figure A40. Water height ratios in a 30° junction for Fr = 2.80 in the narrow channel.
Figure A40. Water height ratios in a 30° junction for Fr = 2.80 in the narrow channel.
Hydrology 07 00057 g0a40
Figure A41. Water height ratios in a 45° junction for Fr = 2.80 in the narrow channel.
Figure A41. Water height ratios in a 45° junction for Fr = 2.80 in the narrow channel.
Hydrology 07 00057 g0a41
Figure A42. Water height ratios in a 90° junction for Fr = 2.80 in the narrow channel.
Figure A42. Water height ratios in a 90° junction for Fr = 2.80 in the narrow channel.
Hydrology 07 00057 g0a42
Figure A43. Water height ratios in a 30° junction for Fr = 3.27 in the narrow channel.
Figure A43. Water height ratios in a 30° junction for Fr = 3.27 in the narrow channel.
Hydrology 07 00057 g0a43
Figure A44. Water height ratios in a 45° junction for Fr = 3.27 in the narrow channel.
Figure A44. Water height ratios in a 45° junction for Fr = 3.27 in the narrow channel.
Hydrology 07 00057 g0a44
Figure A45. Water height ratios in a 90° junction for Fr = 3.27 in the narrow channel.
Figure A45. Water height ratios in a 90° junction for Fr = 3.27 in the narrow channel.
Hydrology 07 00057 g0a45

Appendix B

In Table A1, Table A2, Table A3, Table A4, Table A5, Table A6, Table A7, Table A8 and Table A9, detailed FLOW-3D simulation results are presented. The last column in each table shows whether the impact of the lateral flow on the main channel flow was local (L), or it created waves (W) that propagated upstream (for sub-critical flow condition in the channel (Fr = 0.8)) or downstream (for super-critical flow conditions in the channel (Fr = 1.2, 2)).
Table A1. FLOW-3D simulation results for W = 3.05 m (10′) and Fr = 2.
Table A1. FLOW-3D simulation results for W = 3.05 m (10′) and Fr = 2.
W (m)H0 (m)Fr (-)Q Chan. (m3/s)V Chan. (m/s)D Pipe (m)q Pipe (m3/s)v Pipe (m/s)θ Junc. (deg.) Q V + q v ( cos θ ) q v ( sin θ ) H H 0 Impact
3.053.052101.6310.940.462.54015.489028.281.16L
3.053.052101.6310.940.462.54015.484540.991.12L
3.053.052101.6310.940.462.54015.483058.281.1L
3.053.052101.6310.940.9110.16315.48907.071.76L
3.053.052101.6310.940.915.0837.749028.241.27L
3.053.052101.6310.940.912.5403.8790113.101.17L
3.053.052101.6310.940.9110.16315.484510.991.25L
3.053.052101.6310.940.915.0837.744540.941.2L
3.053.052101.6310.940.912.5403.8745160.951.16L
3.053.052101.6310.940.9110.16315.483015.861.17L
3.053.052101.6310.940.915.0837.743058.221.16L
3.053.052101.6310.940.912.5403.8730227.941.15L
3.053.052101.6310.941.5210.1635.579019.631.68L
3.053.052101.6310.941.525.0832.799078.461.29L
3.053.052101.6310.941.522.5401.3990314.181.15L
3.053.052101.6310.941.5210.1635.574528.751.35L
3.053.052101.6310.941.525.0832.7945111.951.24L
3.053.052101.6310.941.522.5401.3945445.311.15L
3.053.052101.6310.941.5210.1635.573040.981.32L
3.053.052101.6310.941.525.0832.7930158.641.23L
3.053.052101.6310.941.522.5401.3930630.081.15L
3.052.13259.529.150.461.4889.079040.371.17L
3.052.13259.529.150.461.4889.074558.101.17L
3.052.13259.529.150.461.4889.073082.481.17L
3.052.13259.529.150.915.9529.079010.091.89W
3.052.13259.529.150.912.9764.549040.371.32L
3.052.13259.529.150.911.4882.2790161.491.19L
3.052.13259.529.150.915.9529.074515.271.26L
3.052.13259.529.150.912.9764.544558.101.18L
3.052.13259.529.150.911.4882.2745229.391.18L
3.052.13259.529.150.915.9529.073021.921.18L
3.052.13259.529.150.912.9764.543082.481.18L
3.052.13259.529.150.911.4882.2730324.721.18L
3.052.13259.529.151.525.9523.269028.041.84W
3.052.13259.529.151.522.9761.6390112.151.33L
3.052.13259.529.151.521.4880.8290448.591.2L
3.052.13259.529.151.525.9523.264540.651.49L
3.052.13259.529.151.522.9761.6345159.601.28L
3.052.13259.529.151.521.4880.8245635.411.19L
3.052.13259.529.151.525.9523.263057.811.43L
3.052.13259.529.151.522.9761.6330226.031.3L
3.052.13259.529.151.521.4880.8230898.921.17L
3.051.22225.716.920.461.2867.839017.661.38W
3.051.22225.716.920.460.6433.929070.651.23L
3.051.22225.716.920.461.2867.834525.981.26L
3.051.22225.716.920.460.6433.9245100.911.2L
3.051.22225.716.920.461.2867.833037.061.2L
3.051.22225.716.920.460.6433.9230143.031.2L
3.051.22225.716.920.912.5713.929017.661.93W
3.051.22225.716.920.911.2861.969070.651.39W
3.051.22225.716.920.910.6430.9890282.601.17L
3.051.22225.716.920.912.5713.924525.981.53W
3.051.22225.716.920.911.2861.9645100.911.36L
3.051.22225.716.920.910.6430.9845400.661.25L
3.051.22225.716.920.912.5713.923037.061.5L
3.051.22225.716.920.911.2861.9630143.031.36L
3.051.22225.716.920.910.6430.9830566.931.25L
Table A2. FLOW-3D simulation results for W = 7.62 m (25′) and Fr = 2.
Table A2. FLOW-3D simulation results for W = 7.62 m (25′) and Fr = 2.
W
(m)
H0
(m)
Fr (-)Q Chan. (m3/s)V Chan. (m/s)D Pipe (m)q Pipe (m3/s)v Pipe (m/s)θ Junc. (deg.) Q V + q v ( cos θ ) q v ( sin θ ) H H 0 Impact
7.623.052254.0910.940.462.54015.489070.691.1L
7.623.052254.0910.940.462.54015.4845100.981.06L
7.623.052254.0910.940.462.54015.4830143.121.03L
7.623.052254.0910.940.9110.16315.489017.661.38W
7.623.052254.0910.940.915.0837.749070.611.19L
7.623.052254.0910.940.912.5403.8790282.771.1L
7.623.052254.0910.940.9110.16315.484525.981.15L
7.623.052254.0910.940.915.0837.7445100.861.11L
7.623.052254.0910.940.912.5403.8745400.901.1L
7.623.052254.0910.940.9110.16315.483037.061.08L
7.623.052254.0910.940.915.0837.7430142.961.07L
7.623.052254.0910.940.912.5403.8730567.281.06L
7.623.052254.0910.941.5210.1635.579049.071.32W
7.623.052254.0910.941.525.0832.7990196.151.18L
7.623.052254.0910.941.522.5401.3990785.481.1L
7.623.052254.0910.941.5210.1635.574570.391.2L
7.623.052254.0910.941.525.0832.7945278.401.15L
7.623.052254.0910.941.522.5401.39451111.841.1L
7.623.052254.0910.941.5210.1635.573099.861.17L
7.623.052254.0910.941.525.0832.7930394.041.15L
7.623.052254.0910.941.522.5401.39301572.691.1L
7.622.132148.819.150.461.4889.0790100.931.07L
7.622.132148.819.150.461.4889.0745143.741.01L
7.622.132148.819.150.461.4889.0730203.601.01L
7.622.132148.819.150.915.9529.079025.231.36W
7.622.132148.819.150.912.9764.5490100.931.14L
7.622.132148.819.150.911.4882.2790403.731.07L
7.622.132148.819.150.915.9529.074536.691.11L
7.622.132148.819.150.912.9764.5445143.741.09L
7.622.132148.819.150.911.4882.2745571.971.07L
7.622.132148.819.150.915.9529.073052.201.11L
7.622.132148.819.150.912.9764.5430203.601.08L
7.622.132148.819.150.911.4882.2730809.201.06L
7.622.132148.819.151.525.9523.269070.091.43W
7.622.132148.819.151.522.9761.6390280.371.14L
7.622.132148.819.151.521.4880.82901121.481.07L
7.622.132148.819.151.525.9523.2645100.131.22L
7.622.132148.819.151.522.9761.6345397.501.12L
7.622.132148.819.151.521.4880.82451587.011.06L
7.622.132148.819.151.525.9523.2630141.921.2L
7.622.132148.819.151.522.9761.6330562.471.12L
7.622.132148.819.151.521.4880.82302244.701.06L
7.621.22264.286.920.461.2867.839044.161.17L
7.621.22264.286.920.460.6433.9290176.631.07L
7.621.22264.286.920.461.2867.834563.451.08L
7.621.22264.286.920.460.6433.9245250.791.02L
7.621.22264.286.920.461.2867.833090.041.01L
7.621.22264.286.920.460.6433.9230354.981.01L
7.621.22264.286.920.912.5713.929044.161.39W
7.621.22264.286.920.911.2861.9690176.631.17L
7.621.22264.286.920.910.6430.9890706.501.08L
7.621.22264.286.920.912.5713.924563.451.27L
7.621.22264.286.920.911.2861.9645250.791.15L
7.621.22264.286.920.910.6430.98451000.141.07L
7.621.22264.286.920.912.5713.923090.041.25L
7.621.22264.286.920.911.2861.9630354.981.14L
7.621.22264.286.920.910.6430.98301414.731.06L
Table A3. FLOW-3D simulation results for W = 15.24 m (50′) and Fr = 2.
Table A3. FLOW-3D simulation results for W = 15.24 m (50′) and Fr = 2.
W
(m)
H0
(m)
Fr (-)Q Chan. (m3/s)V Chan. (m/s)D Pipe (m)q Pipe (m3/s)v Pipe (m/s)θ Junc. (deg.) Q V + q v ( cos θ ) q v ( sin θ ) H H 0 Impact
15.243.052508.1510.940.462.54015.4890141.381.01L
15.243.052508.1510.940.462.54015.4845200.941.01L
15.243.052508.1510.940.462.54015.4830284.491.01L
15.243.052508.1510.940.9110.16315.489035.331.23L
15.243.052508.1510.940.915.0837.7490141.221.07L
15.243.052508.1510.940.912.5403.8790565.521.01L
15.243.052508.1510.940.9110.16315.484550.961.07L
15.243.052508.1510.940.915.0837.7445200.721.02L
15.243.052508.1510.940.912.5403.8745800.761.01L
15.243.052508.1510.940.9110.16315.483072.381.02L
15.243.052508.1510.940.915.0837.7430284.171.01L
15.243.052508.1510.940.912.5403.87301132.761.01L
15.243.052508.1510.941.5210.1635.579098.131.24L
15.243.052508.1510.941.525.0832.7990392.281.08L
15.243.052508.1510.941.522.5401.39901570.881.03L
15.243.052508.1510.941.5210.1635.5745139.771.14L
15.243.052508.1510.941.525.0832.7945555.771.08L
15.243.052508.1510.941.522.5401.39452222.551.05L
15.243.052508.1510.941.5210.1635.5730197.981.12L
15.243.052508.1510.941.525.0832.7930786.291.06L
15.243.052508.1510.941.522.5401.39303143.481.01L
15.242.132297.619.150.461.4889.0790201.871.05L
15.242.132297.619.150.461.4889.0745286.481.01L
15.242.132297.619.150.461.4889.0730405.471.01L
15.242.132297.619.150.915.9529.079050.471.25L
15.242.132297.619.150.912.9764.5490201.871.11L
15.242.132297.619.150.911.4882.2790807.471.09L
15.242.132297.619.150.915.9529.074572.371.09L
15.242.132297.619.150.912.9764.5445286.481.06L
15.242.132297.619.150.911.4882.27451142.931.02L
15.242.132297.619.150.915.9529.0730102.671.07L
15.242.132297.619.150.912.9764.5430405.471.04L
15.242.132297.619.150.911.4882.27301616.671.01L
15.242.132297.619.151.525.9523.2690140.191.3L
15.242.132297.619.151.522.9761.6390560.741.11L
15.242.132297.619.151.521.4880.82902242.961.07L
15.242.132297.619.151.525.9523.2645199.251.17L
15.242.132297.619.151.522.9761.6345794.011.07L
15.242.132297.619.151.521.4880.82453173.031.02L
15.242.132297.619.151.525.9523.2630282.101.17L
15.242.132297.619.151.522.9761.63301123.211.08L
15.242.132297.619.151.521.4880.82304487.661.02L
15.241.222128.566.920.461.2867.839088.311.08L
15.241.222128.566.920.460.6433.9290353.251.06L
15.241.222128.566.920.461.2867.8345125.891.03L
15.241.222128.566.920.460.6433.9245500.571.01L
15.241.222128.566.920.461.2867.8330178.361.01L
15.241.222128.566.920.460.6433.9230708.231.01L
15.241.222128.566.920.912.5713.929088.311.23L
15.241.222128.566.920.911.2861.9690353.251.06L
15.241.222128.566.920.910.6430.98901413.001.03L
15.241.222128.566.920.912.5713.9245125.891.15L
15.241.222128.566.920.911.2861.9645500.571.04L
15.241.222128.566.920.910.6430.98451999.281.02L
15.241.222128.566.920.912.5713.9230178.361.14L
15.241.222128.566.920.911.2861.9630708.231.03L
15.241.222128.566.920.910.6430.98302827.731.01L
Table A4. FLOW-3D simulation results for W = 3.05 m (10′) and Fr = 1.2.
Table A4. FLOW-3D simulation results for W = 3.05 m (10′) and Fr = 1.2.
W
(m)
H0
(m)
Fr (-)Q Chan. (m3/s)V Chan. (m/s)D Pipe (m)q Pipe (m3/s)v Pipe (m/s)θ Junc. (deg.) Q V + q v ( cos θ ) q v ( sin θ ) H H 0 Impact
3.053.051.260.976.560.461.5249.299028.261.2L
3.053.051.260.976.560.461.5249.294540.971.2L
3.053.051.260.976.560.461.5249.293058.251.2L
3.053.051.260.976.560.916.0979.29907.071.78L
3.053.051.260.976.560.913.0484.649028.261.43L
3.053.051.260.976.560.911.5242.3290113.041.18L
3.053.051.260.976.560.916.0979.294510.991.26L
3.053.051.260.976.560.913.0484.644540.971.23L
3.053.051.260.976.560.911.5242.3245160.861.15L
3.053.051.260.976.560.916.0979.293015.861.22L
3.053.051.260.976.560.913.0484.643058.251.2L
3.053.051.260.976.560.911.5242.3230227.811.15L
3.053.051.260.976.561.526.0973.349019.631.84L
3.053.051.260.976.561.523.0481.679078.501.43L
3.053.051.260.976.561.521.5240.8490314.001.19L
3.053.051.260.976.561.526.0973.344528.751.57L
3.053.051.260.976.561.523.0481.6745112.021.33L
3.053.051.260.976.561.521.5240.8445445.061.17L
3.053.051.260.976.561.526.0973.343040.981.54L
3.053.051.260.976.561.523.0481.6730158.731.32L
3.053.051.260.976.561.521.5240.8430629.731.17L
3.052.131.235.715.490.461.78510.889010.101.28L
3.052.131.235.715.490.460.8935.449040.381.2L
3.052.131.235.715.490.461.78510.884515.281.2L
3.052.131.235.715.490.460.8935.444558.111.2L
3.052.131.235.715.490.461.78510.883021.921.2L
3.052.131.235.715.490.460.8935.443082.501.2L
3.052.131.235.715.490.913.5715.449010.101.85L
3.052.131.235.715.490.911.7852.729040.381.36L
3.052.131.235.715.490.910.8931.3690161.541.15L
3.052.131.235.715.490.913.5715.444515.281.33L
3.052.131.235.715.490.911.7852.724558.111.23L
3.052.131.235.715.490.910.8931.3645229.451.12L
3.052.131.235.715.490.913.5715.443021.921.27L
3.052.131.235.715.490.911.7852.723082.501.22L
3.052.131.235.715.490.910.8931.3630324.811.11L
3.052.131.235.715.491.523.5711.969028.041.87L
3.052.131.235.715.491.521.7850.9890112.181.35L
3.052.131.235.715.491.520.8930.4990448.711.15L
3.052.131.235.715.491.523.5711.964540.661.66L
3.052.131.235.715.491.521.7850.9845159.641.31L
3.052.131.235.715.491.520.8930.4945635.581.13L
3.052.131.235.715.491.523.5711.963057.821.63L
3.052.131.235.715.491.521.7850.9830226.091.31L
3.052.131.235.715.491.520.8930.4930899.161.13L
3.051.221.215.434.150.461.5439.41904.411.65L
3.051.221.215.434.150.460.7724.709017.661.2L
3.051.221.215.434.150.460.3862.359070.621.08L
3.051.221.215.434.150.461.5439.41457.241.1L
3.051.221.215.434.150.460.7724.704525.971.07L
3.051.221.215.434.150.460.3862.3545100.881.03L
3.051.221.215.434.150.461.5439.413010.561.03L
3.051.221.215.434.150.460.7724.703037.041.03L
3.051.221.215.434.150.460.3862.3530142.981.03L
3.051.221.215.434.150.911.5432.359017.661.78L
3.051.221.215.434.150.910.7721.189070.621.28L
3.051.221.215.434.150.910.3860.5990282.501.09L
3.051.221.215.434.150.911.5432.354525.971.52L
3.051.221.215.434.150.910.7721.1845100.881.25L
3.051.221.215.434.150.910.3860.5945400.511.08L
3.051.221.215.434.150.911.5432.353037.041.48L
3.051.221.215.434.150.910.7721.1830142.981.25L
3.051.221.215.434.150.910.3860.5930566.721.07L
Table A5. FLOW-3D simulation results for W = 7.62 m (25′) and Fr = 1.2.
Table A5. FLOW-3D simulation results for W = 7.62 m (25′) and Fr = 1.2.
W
(m)
H0
(m)
Fr (-)Q Chan. (m3/s)V Chan. (m/s)D Pipe (m)q Pipe (m3/s)v Pipe (m/s)θ Junc. (deg.) Q V + q v ( cos θ ) q v ( sin θ ) H H 0 Impact
7.623.051.2152.436.560.461.5249.299070.661.01L
7.623.051.2152.436.560.461.5249.2945100.921.01L
7.623.051.2152.436.560.461.5249.2930143.051.01L
7.623.051.2152.436.560.916.0979.299017.661.27W
7.623.051.2152.436.560.913.0484.649070.661.07L
7.623.051.2152.436.560.911.5242.3290282.631.01L
7.623.051.2152.436.560.916.0979.294525.981.06L
7.623.051.2152.436.560.913.0484.6445100.921.01L
7.623.051.2152.436.560.911.5242.3245400.691.01L
7.623.051.2152.436.560.916.0979.293037.061.01L
7.623.051.2152.436.560.913.0484.6430143.051.01L
7.623.051.2152.436.560.911.5242.3230566.981.01L
7.623.051.2152.436.561.526.0973.349049.071.26L
7.623.051.2152.436.561.523.0481.6790196.271.17L
7.623.051.2152.436.561.521.5240.8490785.071.01L
7.623.051.2152.436.561.526.0973.344570.391.14L
7.623.051.2152.436.561.523.0481.6745278.571.06L
7.623.051.2152.436.561.521.5240.84451111.261.01L
7.623.051.2152.436.561.526.0973.343099.871.13L
7.623.051.2152.436.561.523.0481.6730394.271.05L
7.623.051.2152.436.561.521.5240.84301571.881.01L
7.622.131.289.285.490.461.78510.889025.241.04L
7.622.131.289.285.490.460.8935.4490100.981.01L
7.622.131.289.285.490.461.78510.884536.701.01L
7.622.131.289.285.490.460.8935.4445143.801.01L
7.622.131.289.285.490.461.78510.883052.221.01L
7.622.131.289.285.490.460.8935.4430203.691.01L
7.622.131.289.285.490.913.5715.449025.241.21W
7.622.131.289.285.490.911.7852.7290100.981.06L
7.622.131.289.285.490.910.8931.3690403.911.01L
7.622.131.289.285.490.913.5715.444536.701.04L
7.622.131.289.285.490.911.7852.7245143.801.02L
7.622.131.289.285.490.910.8931.3645572.211.01L
7.622.131.289.285.490.913.5715.443052.221.01L
7.622.131.289.285.490.911.7852.7230203.691.01L
7.622.131.289.285.490.910.8931.3630809.541.01L
7.622.131.289.285.491.523.5711.969070.121.21W
7.622.131.289.285.491.521.7850.9890280.491.05L
7.622.131.289.285.491.520.8930.49901121.961.01L
7.622.131.289.285.491.523.5711.9645100.171.15L
7.622.131.289.285.491.521.7850.9845397.671.04L
7.622.131.289.285.491.520.8930.49451587.691.01L
7.622.131.289.285.491.523.5711.9630141.981.13L
7.622.131.289.285.491.521.7850.9830562.711.04L
7.622.131.289.285.491.520.8930.49302245.661.01L
7.621.221.238.604.150.461.5439.419011.041.25W
7.621.221.238.604.150.460.7724.709044.161.09L
7.621.221.238.604.150.460.3862.3590176.621.01L
7.621.221.238.604.150.461.5439.414516.611.18L
7.621.221.238.604.150.460.7724.704563.451.05L
7.621.221.238.604.150.460.3862.3545250.791.01L
7.621.221.238.604.150.461.5439.413023.811.01L
7.621.221.238.604.150.460.7724.703090.041.01L
7.621.221.238.604.150.460.3862.3530354.981.01L
7.621.221.238.604.150.911.5432.359044.161.27L
7.621.221.238.604.150.910.7721.1890176.621.12L
7.621.221.238.604.150.910.3860.5990706.501.04L
7.621.221.238.604.150.911.5432.354563.451.2L
7.621.221.238.604.150.910.7721.1845250.791.09L
7.621.221.238.604.150.910.3860.59451000.141.02L
7.621.221.238.604.150.911.5432.353090.041.17L
7.621.221.238.604.150.910.7721.1830354.981.08L
7.621.221.238.604.150.910.3860.59301414.731.01L
Table A6. FLOW-3D simulation results for W = 15.24 m (50′) and Fr = 1.2.
Table A6. FLOW-3D simulation results for W = 15.24 m (50′) and Fr = 1.2.
W
(m)
H0
(m)
Fr (-)Q Chan. (m3/s)V Chan. (m/s)D Pipe (m)q Pipe (m3/s)v Pipe (m/s)θ Junc. (deg.) Q V + q v ( cos θ ) q v ( sin θ ) H H 0 Impact
15.243.051.2304.836.560.461.5249.2990141.301.01L
15.243.051.2304.836.560.461.5249.2945200.831.01L
15.243.051.2304.836.560.461.5249.2930284.331.01L
15.243.051.2304.836.560.916.0979.299035.331.1L
15.243.051.2304.836.560.913.0484.6490141.301.02L
15.243.051.2304.836.560.911.5242.3290565.201.01L
15.243.051.2304.836.560.916.0979.294550.961.03L
15.243.051.2304.836.560.913.0484.6445200.831.01L
15.243.051.2304.836.560.911.5242.3245800.311.01L
15.243.051.2304.836.560.916.0979.293072.381.01L
15.243.051.2304.836.560.913.0484.6430284.331.01L
15.243.051.2304.836.560.911.5242.32301132.131.01L
15.243.051.2304.836.561.526.0973.349098.131.1L
15.243.051.2304.836.561.523.0481.6790392.501.02L
15.243.051.2304.836.561.521.5240.84901570.001.01L
15.243.051.2304.836.561.526.0973.3445139.771.03L
15.243.051.2304.836.561.523.0481.6745556.081.02L
15.243.051.2304.836.561.521.5240.84452221.321.01L
15.243.051.2304.836.561.526.0973.3430197.981.05L
15.243.051.2304.836.561.523.0481.6730786.731.02L
15.243.051.2304.836.561.521.5240.84303141.731.01L
15.242.131.2178.545.490.461.78510.889050.481.02L
15.242.131.2178.545.490.460.8935.4490201.921.01L
15.242.131.2178.545.490.461.78510.884572.391.01L
15.242.131.2178.545.490.460.8935.4445286.561.01L
15.242.131.2178.545.490.461.78510.8830102.691.01L
15.242.131.2178.545.490.460.8935.4430405.571.01L
15.242.131.2178.545.490.913.5715.449050.481.11L
15.242.131.2178.545.490.911.7852.7290201.921.03L
15.242.131.2178.545.490.910.8931.3690807.681.01L
15.242.131.2178.545.490.913.5715.444572.391.02L
15.242.131.2178.545.490.911.7852.7245286.561.01L
15.242.131.2178.545.490.910.8931.36451143.241.01L
15.242.131.2178.545.490.913.5715.4430102.691.01L
15.242.131.2178.545.490.911.7852.7230405.571.01L
15.242.131.2178.545.490.910.8931.36301617.101.01L
15.242.131.2178.545.491.523.5711.9690140.221.12L
15.242.131.2178.545.491.521.7850.9890560.891.04L
15.242.131.2178.545.491.520.8930.49902243.571.01L
15.242.131.2178.545.491.523.5711.9645199.311.06L
15.242.131.2178.545.491.521.7850.9845794.221.02L
15.242.131.2178.545.491.520.8930.49453173.891.01L
15.242.131.2178.545.491.523.5711.9630282.181.06L
15.242.131.2178.545.491.521.7850.98301123.521.02L
15.242.131.2178.545.491.520.8930.49304488.871.01L
15.241.221.277.164.150.461.5439.419022.071.19L
15.241.221.277.164.150.460.7724.709088.281.04L
15.241.221.277.164.150.460.3862.3590353.121.01L
15.241.221.277.164.150.461.5439.414532.211.01L
15.241.221.277.164.150.460.7724.7045125.851.01L
15.241.221.277.164.150.460.3862.3545500.391.01L
15.241.221.277.164.150.461.5439.413045.871.01L
15.241.221.277.164.150.460.7724.7030178.291.01L
15.241.221.277.164.150.460.3862.3530707.971.01L
15.241.221.277.164.150.911.5432.359088.281.14L
15.241.221.277.164.150.910.7721.1890353.121.05L
15.241.221.277.164.150.910.3860.59901412.481.02L
15.241.221.277.164.150.911.5432.3545125.851.05L
15.241.221.277.164.150.910.7721.1845500.391.02L
15.241.221.277.164.150.910.3860.59451998.551.01L
15.241.221.277.164.150.911.5432.3530178.291.02L
15.241.221.277.164.150.910.7721.1830707.971.01L
15.241.221.277.164.150.910.3860.59302826.691.01L
Table A7. FLOW-3D simulation results for W = 3.05 m (10′) and Fr = 0.8.
Table A7. FLOW-3D simulation results for W = 3.05 m (10′) and Fr = 0.8.
W
(m)
H0
(m)
Fr (-)Q Chan. (m3/s)V Chan. (m/s)D Pipe (m)q Pipe (m3/s)v Pipe (m/s)θ Junc. (deg.) Q V + q v ( cos θ ) q v ( sin θ ) H H 0 Impact
3.053.050.840.664.380.462.03312.39907.071.33W
3.053.050.840.664.380.461.0176.209028.261.17L
3.053.050.840.664.380.462.03312.394510.991.1L
3.053.050.840.664.380.461.0176.204540.971.1L
3.053.050.840.664.380.462.03312.393015.861.1L
3.053.050.840.664.380.461.0176.203058.251.1L
3.053.050.840.664.380.914.0666.20907.071.71W
3.053.050.840.664.380.912.0333.109028.261.39W
3.053.050.840.664.380.911.0171.5590113.041.19L
3.053.050.840.664.380.914.0666.204510.991.36L
3.053.050.840.664.380.912.0333.104540.971.31L
3.053.050.840.664.380.911.0171.5545160.861.18L
3.053.050.840.664.380.914.0666.203015.861.32L
3.053.050.840.664.380.912.0333.103058.251.29L
3.053.050.840.664.380.911.0171.5530227.811.18L
3.053.050.840.664.381.524.0662.239019.631.7W
3.053.050.840.664.381.522.0331.129078.501.4W
3.053.050.840.664.381.521.0170.5690314.001.2L
3.053.050.840.664.381.524.0662.234528.751.58W
3.053.050.840.664.381.522.0331.1245112.021.35W
3.053.050.840.664.381.521.0170.5645445.061.2L
3.053.050.840.664.381.524.0662.233040.981.54W
3.053.050.840.664.381.522.0331.1230158.731.35W
3.053.050.840.664.381.521.0170.5630629.731.2L
3.052.130.823.813.660.462.38114.51902.521.72W
3.052.130.823.813.660.461.1917.269010.091.36W
3.052.130.823.813.660.460.5953.639040.361.14L
3.052.130.823.813.660.462.38114.51454.571.08L
3.052.130.823.813.660.461.1917.264515.271.08L
3.052.130.823.813.660.460.5953.634558.071.08L
3.052.130.823.813.660.462.38114.51306.781.05L
3.052.130.823.813.660.461.1917.263021.911.05L
3.052.130.823.813.660.460.5953.633082.451.05L
3.052.130.823.813.660.912.3813.639010.091.76W
3.052.130.823.813.660.911.1911.819040.361.41W
3.052.130.823.813.660.910.5950.9190161.431.16L
3.052.130.823.813.660.912.3813.634515.271.47W
3.052.130.823.813.660.911.1911.814558.071.32L
3.052.130.823.813.660.910.5950.9145229.291.16L
3.052.130.823.813.660.912.3813.633021.911.42W
3.052.130.823.813.660.911.1911.813082.451.31L
3.052.130.823.813.660.910.5950.9130324.591.15L
3.052.130.823.813.661.522.3811.319028.031.73W
3.052.130.823.813.661.521.1910.6590112.101.4W
3.052.130.823.813.661.520.5950.3390448.411.17L
3.052.130.823.813.661.522.3811.314540.631.7W
3.052.130.823.813.661.521.1910.6545159.541.35W
3.052.130.823.813.661.520.5950.3345635.151.17L
3.052.130.823.813.661.522.3811.313057.781.64W
3.052.130.823.813.661.521.1910.6530225.941.35W
3.052.130.823.813.661.520.5950.3330898.551.17L
3.051.220.810.282.770.461.0286.26904.421.72W
3.051.220.810.282.770.460.5143.139017.671.34L
3.051.220.810.282.770.460.2571.579070.691.15L
3.051.220.810.282.770.461.0286.26457.251.18L
3.051.220.810.282.770.460.5143.134525.991.18L
3.051.220.810.282.770.460.2571.5745100.971.09L
3.051.220.810.282.770.461.0286.263010.571.18L
3.051.220.810.282.770.460.5143.133037.081.18L
3.051.220.810.282.770.460.2571.5730143.111.09L
3.051.220.810.282.770.911.0281.579017.671.65W
3.051.220.810.282.770.910.5140.789070.691.38W
3.051.220.810.282.770.910.2570.3990282.761.17L
3.051.220.810.282.770.911.0281.574525.991.63W
3.051.220.810.282.770.910.5140.7845100.971.34W
3.051.220.810.282.770.910.2570.3945400.881.17L
3.051.220.810.282.770.911.0281.573037.081.55W
3.051.220.810.282.770.910.5140.7830143.111.31W
3.051.220.810.282.770.910.2570.3930567.241.17L
Table A8. FLOW-3D simulation results for W = 7.62 m (25′) and Fr = 0.8.
Table A8. FLOW-3D simulation results for W = 7.62 m (25′) and Fr = 0.8.
W
(m)
H0
(m)
Fr (-)Q Chan. (m3/s)V Chan. (m/s)D Pipe (m)q Pipe (m3/s)v Pipe (m/s)θ Junc. (deg.) Q V + q v ( cos θ ) q v ( sin θ ) H H 0 Impact
7.623.050.8101.664.380.462.03312.399017.661.14L
7.623.050.8101.664.380.461.0176.209070.651.01L
7.623.050.8101.664.380.462.03312.394525.981.01L
7.623.050.8101.664.380.461.0176.2045100.911.01L
7.623.050.8101.664.380.462.03312.393037.061.01L
7.623.050.8101.664.380.461.0176.2030143.031.01L
7.623.050.8101.664.380.914.0666.209017.661.4L
7.623.050.8101.664.380.912.0333.109070.651.16L
7.623.050.8101.664.380.911.0171.5590282.601.04L
7.623.050.8101.664.380.914.0666.204525.981.14L
7.623.050.8101.664.380.912.0333.1045100.911.08L
7.623.050.8101.664.380.911.0171.5545400.661.01L
7.623.050.8101.664.380.914.0666.203037.061.09L
7.623.050.8101.664.380.912.0333.1030143.031.07L
7.623.050.8101.664.380.911.0171.5530566.931.01L
7.623.050.8101.664.381.524.0662.239049.061.41L
7.623.050.8101.664.381.522.0331.1290196.251.17L
7.623.050.8101.664.381.521.0170.5690785.001.04L
7.623.050.8101.664.381.524.0662.234570.381.28L
7.623.050.8101.664.381.522.0331.1245278.541.15L
7.623.050.8101.664.381.521.0170.56451111.161.04L
7.623.050.8101.664.381.524.0662.233099.861.26L
7.623.050.8101.664.381.522.0331.1230394.231.13L
7.623.050.8101.664.381.521.0170.56301571.731.04L
7.622.130.859.553.660.462.38114.51906.311.25L
7.622.130.859.553.660.461.1917.269025.231.11L
7.622.130.859.553.660.460.5953.6390100.921.01L
7.622.130.859.553.660.462.38114.51459.921.01L
7.622.130.859.553.660.461.1917.264536.681.01L
7.622.130.859.553.660.460.5953.6345143.721.01L
7.622.130.859.553.660.462.38114.513014.351.01L
7.622.130.859.553.660.461.1917.263052.191.01L
7.622.130.859.553.660.460.5953.6330203.571.01L
7.622.130.859.553.660.912.3813.639025.231.38L
7.622.130.859.553.660.911.1911.8190100.921.13L
7.622.130.859.553.660.910.5950.9190403.671.03L
7.622.130.859.553.660.912.3813.634536.681.18L
7.622.130.859.553.660.911.1911.8145143.721.09L
7.622.130.859.553.660.910.5950.9145571.871.01L
7.622.130.859.553.660.912.3813.633052.191.17L
7.622.130.859.553.660.911.1911.8130203.571.08L
7.622.130.859.553.660.910.5950.9130809.061.01L
7.622.130.859.553.661.522.3811.319070.081.38L
7.622.130.859.553.661.521.1910.6590280.321.14L
7.622.130.859.553.661.520.5950.33901121.301.04L
7.622.130.859.553.661.522.3811.3145100.111.3L
7.622.130.859.553.661.521.1910.6545397.441.13L
7.622.130.859.553.661.520.5950.33451586.751.03L
7.622.130.859.553.661.522.3811.3130141.891.29L
7.622.130.859.553.661.521.1910.6530562.381.13L
7.622.130.859.553.661.520.5950.33302244.321.02L
7.621.220.825.712.770.461.0286.269011.051.3L
7.621.220.825.712.770.460.5143.139044.201.16L
7.621.220.825.712.770.460.2571.5790176.821.04L
7.621.220.825.712.770.461.0286.264516.631.11L
7.621.220.825.712.770.460.5143.134563.521.06L
7.621.220.825.712.770.460.2571.5745251.061.02L
7.621.220.825.712.770.461.0286.263023.831.04L
7.621.220.825.712.770.460.5143.133090.141.03L
7.621.220.825.712.770.460.2571.5730355.371.01L
7.621.220.825.712.770.911.0281.579044.201.35L
7.621.220.825.712.770.910.5140.7890176.821.15L
7.621.220.825.712.770.910.2570.3990707.281.05L
7.621.220.825.712.770.911.0281.574563.521.28L
7.621.220.825.712.770.910.5140.7845251.061.12L
7.621.220.825.712.770.910.2570.39451001.241.03L
7.621.220.825.712.770.911.0281.573090.141.27L
7.621.220.825.712.770.910.5140.7830355.371.13L
7.621.220.825.712.770.910.2570.39301416.291.01L
Table A9. FLOW-3D simulation results for W = 15.24 m (50′) and Fr = 0.8.
Table A9. FLOW-3D simulation results for W = 15.24 m (50′) and Fr = 0.8.
W
(m)
H0
(m)
Fr (-)Q Chan. (m3/s)V Chan. (m/s)D Pipe (m)q Pipe (m3/s)v Pipe (m/s)θ Junc. (deg.) Q V + q v ( cos θ ) q v ( sin θ ) H H 0 Impact
15.243.050.8203.324.380.462.03312.399035.331.03L
15.243.050.8203.324.380.461.0176.2090141.301.01L
15.243.050.8203.324.380.462.03312.394550.961.01L
15.243.050.8203.324.380.461.0176.2045200.831.01L
15.243.050.8203.324.380.462.03312.393072.381.01L
15.243.050.8203.324.380.461.0176.2030284.331.01L
15.243.050.8203.324.380.914.0666.209035.331.18L
15.243.050.8203.324.380.912.0333.1090141.301.05L
15.243.050.8203.324.380.911.0171.5590565.201.01L
15.243.050.8203.324.380.914.0666.204550.961.04L
15.243.050.8203.324.380.912.0333.1045200.831.01L
15.243.050.8203.324.380.911.0171.5545800.311.01L
15.243.050.8203.324.380.914.0666.203072.381.01L
15.243.050.8203.324.380.912.0333.1030284.331.01L
15.243.050.8203.324.380.911.0171.55301132.131.01L
15.243.050.8203.324.381.524.0662.239098.131.17L
15.243.050.8203.324.381.522.0331.1290392.501.04L
15.243.050.8203.324.381.521.0170.56901570.001.01L
15.243.050.8203.324.381.524.0662.2345139.771.1L
15.243.050.8203.324.381.522.0331.1245556.081.03L
15.243.050.8203.324.381.521.0170.56452221.321.01L
15.243.050.8203.324.381.524.0662.2330197.981.09L
15.243.050.8203.324.381.522.0331.1230786.731.01L
15.243.050.8203.324.381.521.0170.56303141.731.01L
15.242.130.8119.073.660.462.38114.519012.611.16L
15.242.130.8119.073.660.461.1917.269050.451.03L
15.242.130.8119.073.660.460.5953.6390201.791.01L
15.242.130.8119.073.660.462.38114.514518.841.01L
15.242.130.8119.073.660.461.1917.264572.341.01L
15.242.130.8119.073.660.460.5953.6345286.371.01L
15.242.130.8119.073.660.462.38114.513026.961.01L
15.242.130.8119.073.660.461.1917.2630102.621.01L
15.242.130.8119.073.660.460.5953.6330405.301.01L
15.242.130.8119.073.660.912.3813.639050.451.18L
15.242.130.8119.073.660.911.1911.8190201.791.07L
15.242.130.8119.073.660.910.5950.9190807.141.01L
15.242.130.8119.073.660.912.3813.634572.341.08L
15.242.130.8119.073.660.911.1911.8145286.371.03L
15.242.130.8119.073.660.910.5950.91451142.471.01L
15.242.130.8119.073.660.912.3813.6330102.621.08L
15.242.130.8119.073.660.911.1911.8130405.301.02L
15.242.130.8119.073.660.910.5950.91301616.011.01L
15.242.130.8119.073.661.522.3811.3190140.131.17L
15.242.130.8119.073.661.521.1910.6590560.511.06L
15.242.130.8119.073.661.520.5950.33902242.061.02L
15.242.130.8119.073.661.522.3811.3145199.171.12L
15.242.130.8119.073.661.521.1910.6545793.691.04L
15.242.130.8119.073.661.520.5950.33453171.751.01L
15.242.130.8119.073.661.522.3811.3130281.991.12L
15.242.130.8119.073.661.521.1910.65301122.761.05L
15.242.130.8119.073.661.520.5950.33304485.851.01L
15.241.220.851.402.770.461.0286.269022.091.2L
15.241.220.851.402.770.460.5143.139088.361.09L
15.241.220.851.402.770.460.2571.5790353.441.02L
15.241.220.851.402.770.461.0286.264532.241.07L
15.241.220.851.402.770.460.5143.1345125.961.04L
15.241.220.851.402.770.460.2571.5745500.851.01L
15.241.220.851.402.770.461.0286.263045.911.02L
15.241.220.851.402.770.460.5143.1330178.451.01L
15.241.220.851.402.770.460.2571.5730708.621.01L
15.241.220.851.402.770.911.0281.579088.361.24L
15.241.220.851.402.770.910.5140.7890353.441.08L
15.241.220.851.402.770.910.2570.39901413.781.02L
15.241.220.851.402.770.911.0281.5745125.961.19L
15.241.220.851.402.770.910.5140.7845500.851.07L
15.241.220.851.402.770.910.2570.39452000.381.01L
15.241.220.851.402.770.911.0281.5730178.451.16L
15.241.220.851.402.770.910.5140.7830708.621.01L
15.241.220.851.402.770.910.2570.39302829.291.01L

References

  1. Gurram, S.K.; Karki, K.S.; Hager, W.H. Subcritical junction flow. J. Hydraul. Eng. 1997, 123, 447–455. [Google Scholar] [CrossRef]
  2. Best, J.L. Flow Dynamics and Sediment Transport at River Channel Confluences. Ph.D. Thesis, University of London, London, UK, 1985. [Google Scholar]
  3. U.S. Army Engineer District. Supercritical Flow at Open Channel Junction; Report No. 2-100; U.S. Army Engineer District: Los Angeles, CA, USA, 1975. [Google Scholar]
  4. Taylor, E.H. Flow characteristics at rectangular open-channel junctions. Trans. ASCE 1944, 109, 893–902. [Google Scholar]
  5. Mignot, E.; Bonakdari, H.; Knothe, P.; Liperne Kouyi, G.; Bessette, A.; Riviere, N.; Bertrand-Krajewski, J.L. Experiments and 3D simulations of flow structures in junctions and their influence on location of flowmeters. Water Sci. Technol. 2012, 66, 1325–1332. [Google Scholar] [CrossRef] [PubMed]
  6. Constantinescu, G.; Miyawaki, S.; Rhoads, B.; Sukhodolov, A. Numerical evaluation of the effects of planform geometry and inflow conditions on flow, turbulence structure, and bed shear velocity at a stream confluence with a concordant bed. J. Geophys. Res. 2014, 119, 2079–2097. [Google Scholar] [CrossRef]
  7. Webber, N.B.; Greated, C.A. An investigation of flow behaviour at the junction of rectangular channels. ICE Proc. 1966, 34, 321–334. [Google Scholar] [CrossRef]
  8. Hsu, C.; Wu, F.; Lee, W. Flow at 90° equal-width open-channel junction. J. Hydraul. Eng. 1998, 124, 186–191. [Google Scholar] [CrossRef]
  9. Nedelec, Y.; Gay, B. Experimental study of a right-angled end junction between a pipe and an open channel. J. Hydraul. Eng. 2008, 134. [Google Scholar] [CrossRef]
  10. GREYLINE AVFM 5.0 User’s Manual (2017) [Online]. Available online: http://www.greyline.com/pdf/AVFM%205.0%20Manual.pdf (accessed on 15 February 2018).
  11. Lin, J.; Soong, H. Junction losses in open channel flows. Water Resour. Res. 1979, 15, 414–418. [Google Scholar] [CrossRef]
  12. Garcia-Navarro, M.; Saviron, J. Numerical simulation of unsteady flow at open channel junctions. J. Hydraul. Res. 1992, 30, 595–609. [Google Scholar] [CrossRef]
  13. Mays, L.W. Stormwater Collection Systems Design Book; McGraw-Hill: New York, NY, USA, 2011. [Google Scholar]
  14. Yen, B.C. Hydraulics of sewers. Adv. Hydrosci. 1986, 14, 1122. [Google Scholar]
  15. Nazari-Sharabian, M.; Nazari-Sharabian, A.; Karakouzian, M.; Karami, M. Sacrificial Piles as Scour Countermeasures in River Bridges—A Numerical Study using FLOW-3D. Civ. Eng. J. 2020, 6. [Google Scholar] [CrossRef]
  16. Karami, M.; Kabiri-Samani, A.; Nazari-Sharabian, M.; Karakouzian, M. Investigating the Effects of Transient Flow in Concrete-Lined Pressure Tunnels and Developing a New Analytical Formula for Pressure Wave Velocity. Tunn. Undergr. Space Technol. 2019, 91, 102992. [Google Scholar] [CrossRef]
  17. Flow Science, Inc. Santa Fe, NM, USA. FLOW-3D® Version 12.0 User’s Manual (2018) [Online]. Available online: https://www.flow3d.com (accessed on 3 February 2019).
  18. Hydrologic Criteria and Drainage Design Manual (HCDDM) of the Clark County Regional Flood Control District (CCRFCD). Available online: http://www.ccrfcd.org/hcddm.htm (accessed on 3 February 2019).
  19. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models. Part 1: A discussion of principles. J. Hydrol. 1970, 103, 282–290. [Google Scholar] [CrossRef]
  20. Kenworthy, S.; Rhoads, B. Hydrologic control of spatial patterns of suspended sediment concentration at a stream confluence. J. Hydrol. 1995, 168, 251–263. [Google Scholar] [CrossRef]
Figure 1. The confluence of an open channel with a lateral drainage pipe.
Figure 1. The confluence of an open channel with a lateral drainage pipe.
Hydrology 07 00057 g001
Figure 2. The flow chart of this study.
Figure 2. The flow chart of this study.
Hydrology 07 00057 g002
Figure 3. The 2” diameter inlet pipe installed on the flume wall.
Figure 3. The 2” diameter inlet pipe installed on the flume wall.
Hydrology 07 00057 g003
Figure 4. Wide channel configuration (top view).
Figure 4. Wide channel configuration (top view).
Hydrology 07 00057 g004
Figure 5. Narrow channel configuration (top view).
Figure 5. Narrow channel configuration (top view).
Hydrology 07 00057 g005
Figure 6. Lateral flow to the main channel (back view).
Figure 6. Lateral flow to the main channel (back view).
Hydrology 07 00057 g006
Figure 7. Water heights in the wide channel, without lateral flow.
Figure 7. Water heights in the wide channel, without lateral flow.
Hydrology 07 00057 g007
Figure 8. Water heights in the wide channel, with 2.5% lateral flow.
Figure 8. Water heights in the wide channel, with 2.5% lateral flow.
Hydrology 07 00057 g008
Figure 9. Water heights in the wide channel, with 5.0% lateral flow.
Figure 9. Water heights in the wide channel, with 5.0% lateral flow.
Hydrology 07 00057 g009
Figure 10. Maximum water height ratios in a 90° junction (wide channel).
Figure 10. Maximum water height ratios in a 90° junction (wide channel).
Hydrology 07 00057 g010
Figure 11. Maximum water height ratios in a 90° junction (narrow channel).
Figure 11. Maximum water height ratios in a 90° junction (narrow channel).
Hydrology 07 00057 g011
Figure 12. Maximum water height ratios in a 45° junction (wide channel).
Figure 12. Maximum water height ratios in a 45° junction (wide channel).
Hydrology 07 00057 g012
Figure 13. Maximum water height ratios in a 45° junction (narrow channel).
Figure 13. Maximum water height ratios in a 45° junction (narrow channel).
Hydrology 07 00057 g013
Figure 14. Maximum water height ratios in a 30° junction (wide channel).
Figure 14. Maximum water height ratios in a 30° junction (wide channel).
Hydrology 07 00057 g014
Figure 15. Maximum water height ratios in a 30° junction (narrow channel).
Figure 15. Maximum water height ratios in a 30° junction (narrow channel).
Hydrology 07 00057 g015
Figure 16. 3D view of the channel components and mesh planes.
Figure 16. 3D view of the channel components and mesh planes.
Hydrology 07 00057 g016
Figure 17. Meshing the channel, using four different grid sizes.
Figure 17. Meshing the channel, using four different grid sizes.
Hydrology 07 00057 g017
Figure 18. 3D favored view of the channel.
Figure 18. 3D favored view of the channel.
Hydrology 07 00057 g018
Figure 19. (a,c): Pipe defined as a hole that passes through the channel wall and the pipe complement; (b,d): Favored view of the lateral inlet pipe.
Figure 19. (a,c): Pipe defined as a hole that passes through the channel wall and the pipe complement; (b,d): Favored view of the lateral inlet pipe.
Hydrology 07 00057 g019
Figure 20. Boundary conditions in the simulated channel.
Figure 20. Boundary conditions in the simulated channel.
Hydrology 07 00057 g020
Figure 21. The convergence of the numerical modeling solution.
Figure 21. The convergence of the numerical modeling solution.
Hydrology 07 00057 g021
Figure 22. Time-step information of the simulation.
Figure 22. Time-step information of the simulation.
Hydrology 07 00057 g022
Figure 23. Water heights in the flume and based on FLOW-3D results, for Fr = 0.8.
Figure 23. Water heights in the flume and based on FLOW-3D results, for Fr = 0.8.
Hydrology 07 00057 g023
Figure 24. Water heights in the flume and based on FLOW-3D results, for Fr = 2.0.
Figure 24. Water heights in the flume and based on FLOW-3D results, for Fr = 2.0.
Hydrology 07 00057 g024
Figure 25. Water heights in the flume and based on FLOW-3D results, for Fr = 3.1.
Figure 25. Water heights in the flume and based on FLOW-3D results, for Fr = 3.1.
Hydrology 07 00057 g025
Figure 26. Back-water effect in a channel with sub-critical flow.
Figure 26. Back-water effect in a channel with sub-critical flow.
Hydrology 07 00057 g026
Figure 27. Splashing on the outer wall of a channel with super-critical flow and 90° junction angle.
Figure 27. Splashing on the outer wall of a channel with super-critical flow and 90° junction angle.
Hydrology 07 00057 g027
Figure 28. Simulation results for Fr = 2 in the channel.
Figure 28. Simulation results for Fr = 2 in the channel.
Hydrology 07 00057 g028
Figure 29. Conservative design curves for Fr = 2 in the channel.
Figure 29. Conservative design curves for Fr = 2 in the channel.
Hydrology 07 00057 g029
Figure 30. Simulation results for Fr = 1.2 in the channel.
Figure 30. Simulation results for Fr = 1.2 in the channel.
Hydrology 07 00057 g030
Figure 31. Conservative design curves for Fr = 1.2 in the channel.
Figure 31. Conservative design curves for Fr = 1.2 in the channel.
Hydrology 07 00057 g031
Figure 32. Simulation results for Fr = 0.8 in the channel.
Figure 32. Simulation results for Fr = 0.8 in the channel.
Hydrology 07 00057 g032
Figure 33. Conservative design curves for Fr = 0.8 in the channel.
Figure 33. Conservative design curves for Fr = 0.8 in the channel.
Hydrology 07 00057 g033
Table 1. Experimental scenarios.
Table 1. Experimental scenarios.
Flume Width (m)Flume Slope (%)Flow Rate (m3/s)Fr (-)Junction Angle (degree)Lateral Flow Proportion (%)
0.41 (16″)0.000.130.7230, 45, 902.5, 5.0, 7.5
0.420.120.9530, 45, 902.5, 5.0, 7.5
0.680.141.4530, 45, 902.5, 5.0, 7.5
0.970.141.7430, 45, 902.5, 5.0, 7.5
1.680.142.0130, 45, 902.5, 5.0, 7.5
2.150.132.3230, 45, 902.5, 5.0, 7.5
2.600.132.8030, 45, 902.5, 5.0, 7.5
3.510.133.2730, 45, 902.5, 5.0, 7.5
0.61 (24″)0.000.180.8030, 45, 902.5, 5.0
0.440.181.030, 45, 902.5, 5.0
0.560.181.230, 45, 902.5, 5.0
0.890.181.630, 45, 902.5, 5.0
1.670.182.030, 45, 902.5, 5.0
2.220.182.430, 45, 902.5, 5.0
3.330.183.130, 45, 902.5, 5.0
Total number of tests: 120
Table 2. Simulation scenarios.
Table 2. Simulation scenarios.
Width (m)Initial Water Height (m)Fr (-) Inlet Pipe Diameter (m)Junction Angle (Degree)Lateral Flow Proportion (%)
3.05 (10′)1.22 (4′)0.8, 1.2, 2.00.46 (1.5′), 0.91 (3′)30, 45, 902.5, 5, 10
2.13 (7′)0.8, 1.2, 2.00.46 (1.5′), 0.91 (3′), 1.52 (5′)30, 45, 902.5, 5, 10
3.05 (10′)0.8, 1.2, 2.00.46 (1.5′), 0.91 (3′), 1.52 (5′)30, 45, 902.5, 5, 10
7.62 (25′)1.22 (4′)0.8, 1.2, 2.00.46 (1.5′), 0.91 (3′)30, 45, 902.5, 5, 10
2.13 (7′)0.8, 1.2, 2.00.46 (1.5′), 0.91 (3′), 1.52 (5′)30, 45, 902.5, 5, 10
3.05 (10′)0.8, 1.2, 2.00.46 (1.5′), 0.91 (3′), 1.52 (5′)30, 45, 902.5, 5, 10
15.24 (50′)1.22 (4′)0.8, 1.2, 2.00.46 (1.5′), 0.91 (3′)30, 45, 902.5, 5, 10
2.13 (7′)0.8, 1.2, 2.00.46 (1.5′), 0.91 (3′), 1.52 (5′)30, 45, 902.5, 5, 10
3.05 (10′)0.8, 1.2, 2.00.46 (1.5′), 0.91 (3′), 1.52 (5′)30, 45, 902.5, 5, 10
Total number of simulations: 648
Table 3. Experimental results for the 90° junction in the wide channel.
Table 3. Experimental results for the 90° junction in the wide channel.
Fr
(-)
Q Flume (m3/s)V Flume (m/s)θ Junc.
(Deg.)
2.5% Lateral Flow5% Lateral Flow
q Pipe (m3/s)v Pipe (m/s)H/H0q Pipe (m3/s)v Pipe (m/s)H/H0
0.80.181.37900.00452.231.050.00914.451.1
10.181.55900.00452.231.160.00914.451.26
1.20.181.65900.00452.231.460.00914.451.63
1.60.181.86900.00452.231.290.00914.451.57
20.182.10900.00452.231.220.00914.451.74
Table 4. Experimental results for the 90° junction in the narrow channel.
Table 4. Experimental results for the 90° junction in the narrow channel.
Fr
(-)
Q Flume(m3/s)V Flume(m/s)θ Junc. (deg.)2.5% Lateral Flow5% Lateral Flow7.5% Lateral Flow
q Pipe (m3/s)v Pipe (m/s)H/H0q Pipe (m3/s)v Pipe (m/s)H/H0q Pipe (m3/s)v Pipe (m/s)H/H0
0.720.13221.16900.00341.651.140.00653.261.170.00994.911.20
0.950.11581.34900.00281.431.300.00572.871.400.00884.301.48
1.450.13681.89900.00341.681.120.00683.381.580.01025.061.73
1.740.13852.13900.00341.711.130.00683.411.290.01055.121.71
2.010.13992.35900.00341.741.160.00713.441.320.01055.181.26
2.320.12942.53900.00311.581.240.00653.201.290.00964.791.53
2.80.13222.87900.00341.651.200.00653.261.600.00994.911.87
3.270.12803.17900.00311.581.290.00653.171.860.00964.752.50
Table 5. Experimental results for the 45° junction in the wide channel.
Table 5. Experimental results for the 45° junction in the wide channel.
Fr
(-)
Q Flume (m3/s)V Flume (m/s)θ Junc.
(Deg.)
2.5% Lateral Flow5% Lateral Flow
q Pipe (m3/s)v Pipe (m/s)H/H0q Pipe (m3/s)v Pipe (m/s)H/H0
0.80.181.37450.00452.231.030.09754.451.06
10.181.55450.00452.231.140.09754.451.14
1.20.181.65450.00452.231.160.09754.451.2
1.60.181.86450.00452.231.190.09754.451.3
20.182.10450.00452.231.220.09754.451.28
Table 6. Experimental results for the 45° junction in the narrow channel.
Table 6. Experimental results for the 45° junction in the narrow channel.
Fr
(-)
Q Flume (m3/s)V Flume (m/s)θ Junc.
(Deg.)
2.5% Lateral Flow5% Lateral Flow7.5% Lateral Flow
q Pipe (m3/s)v Pipe (m/s)H/H0q Pipe (m3/s)v Pipe (m/s)H/H0q Pipe (m3/s)v Pipe (m/s)H/H0
0.720.13221.16450.00341.651.120.00653.261.130.00994.911.13
0.950.11581.34450.00281.431.160.00572.871.160.00884.301.15
1.450.13681.89450.00341.681.040.00683.381.160.01025.061.08
1.740.13852.13450.00341.711.190.00683.411.130.01055.121.17
2.010.13992.35450.00341.741.110.00713.441.160.01055.181.21
2.320.12942.53450.00311.581.120.00653.201.120.00964.791.24
2.80.13222.87450.00341.651.400.00653.261.330.00994.911.53
3.270.12803.17450.00311.581.360.00653.171.290.00964.751.36
Table 7. Experimental results for the 30° junction in the wide channel.
Table 7. Experimental results for the 30° junction in the wide channel.
Fr
(-)
Q Flume (m3/s)V Flume (m/s)θ Junc.
(Deg.)
2.5% Lateral Flow5% Lateral Flow
q Pipe (m3/s)v Pipe (m/s)H/H0q Pipe (m3/s)v Pipe (m/s)H/H0
0.80.181.37300.00452.231.030.0094.451.03
10.181.55300.00452.231.100.0094.451.10
1.20.181.65300.00452.231.120.0094.451.12
1.60.181.86300.00452.231.100.0094.451.10
20.182.10300.00452.231.110.0094.451.11
Table 8. Experimental results for the 30° junction in the narrow channel.
Table 8. Experimental results for the 30° junction in the narrow channel.
Fr
(-)
Q Flume (m3/s)V Flume (m/s)θ Junc.
(Deg.)
2.5% Lateral Flow5% Lateral Flow7.5% Lateral Flow
q Pipe (m3/s)v Pipe (m/s)H/H0q Pipe (m3/s)v Pipe (m/s)H/H0q Pipe (m3/s)v Pipe (m/s)H/H0
0.720.13221.16300.00341.651.130.00653.261.130.00994.911.13
0.950.11581.34300.00281.431.130.00572.871.180.00884.301.18
1.450.13681.89300.00341.681.120.00683.381.170.01025.061.08
1.740.13852.13300.00341.711.190.00683.411.190.01055.121.24
2.010.13992.35300.00341.741.160.00713.441.210.01055.181.21
2.320.12942.53300.00311.581.290.00653.201.290.00964.791.29
2.80.13222.87300.00341.651.400.00653.261.330.00994.911.40
3.270.12803.17300.00311.581.200.00653.171.290.00964.751.27
Table 9. Grid sizes in each mesh block in the channel.
Table 9. Grid sizes in each mesh block in the channel.
Mesh BlockGrid Size (cm)
12
21.5
31
4 (Pipe)0.7
Table 10. FLOW-3D validation results.
Table 10. FLOW-3D validation results.
Froude NumberNash–Sutcliffe Efficiency
0.80.87
2.00.83
3.10.85

Share and Cite

MDPI and ACS Style

Nazari-Sharabian, M.; Karakouzian, M.; Hayes, D. Flow Topology in the Confluence of an Open Channel with Lateral Drainage Pipe. Hydrology 2020, 7, 57. https://doi.org/10.3390/hydrology7030057

AMA Style

Nazari-Sharabian M, Karakouzian M, Hayes D. Flow Topology in the Confluence of an Open Channel with Lateral Drainage Pipe. Hydrology. 2020; 7(3):57. https://doi.org/10.3390/hydrology7030057

Chicago/Turabian Style

Nazari-Sharabian, Mohammad, Moses Karakouzian, and Donald Hayes. 2020. "Flow Topology in the Confluence of an Open Channel with Lateral Drainage Pipe" Hydrology 7, no. 3: 57. https://doi.org/10.3390/hydrology7030057

APA Style

Nazari-Sharabian, M., Karakouzian, M., & Hayes, D. (2020). Flow Topology in the Confluence of an Open Channel with Lateral Drainage Pipe. Hydrology, 7(3), 57. https://doi.org/10.3390/hydrology7030057

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop