Lateral Shear Layer and Its Velocity Distribution of Flow in Rectangular Open Channels

The lateral velocity distribution of flow in the shear layer of open channel is required to many problems in river and eco-environment engineering, e.g. distribution of pollutant dispersion, sediment transport and bank erosion, and aquatic habitat. It is not well understood about how the velocity varies laterally in the wall boundary layer. This paper gives an analytical solution of lateral velocity distribution in a rectangular open channel based on the depth-averaged momentum equation proposed by Shiono & Knight. The obtained lateral velocity distributions in the wall shear layer are related to the two hydraulic parameters of lateral eddy viscosity (λ) and depth-averaged secondary flow (Γ) for given roughened channels. Preliminary relationships between the above two parameters and the aspect ratio of channel, B/H, are obtained from two sets of experimental data. The lateral width (δ) of the shear layer was investigated and found to relate to the λ and the bed friction factor (f), as described by Equation (26). This study indicates that the lateral shear layer near the wall can be very wide (δ/H = 14.6) for the extreme case (λ = 0.6 and f = 0.01).


Introduction
In open channel flow, like natural rivers and canals, flows are of highly three-dimensional characteristics owing to the influence of boundaries, varying roughness and non-uniform shapes.Even modeling flow in straight prismatic see www.flowdata.bham.ac.uk.The results of the lateral width of the shear layer indicate that the shear layer near the wall could be very wide for natural rivers when the λ is very large.
In the following sections, Section 2 provides theoretical background with the assumption of model and its analytical solution, followed by experimental data used for model parameter evaluation in Section 3. Section 4 shows the detailed comparison between the analytical model and experiment, along with the discussion on the lateral extent of shear layer.Some conclusions are given in Section 5.

Introduction
For steady flow in a prismatic open channel, the governing momentum equation in the stream wise direction may be combined with the continuity equation to give [8] [16]: where U, V, W are the time-averaged velocity components in the x (streamwise), y (lateral) and z (normal to bed) directions respectively (see Figure 1), S o is the X.Tang For steady flow in a rectangular channel, consider that and by integrating Equation (1) over the water depth (H), the depth-averaged momentum equation can be obtained as (Shiono & Knight [8]): where τ b is the bed shear stress, the overbar or subscript d denotes the depth-averaged value, and in which u and v are the fluctuating velocity component in x and y direction respectively.Through an eddy viscosity concept, the Reynolds stress, τ yx , can be described as: where λ is the dimensionless eddy viscosity coefficient and U * ( τ ρ = b ) is the local shear velocity.The bed shear stress, τ b , may be related to the depth-averaged velocity, U d , using the Darcy-Weisbach friction coefficient (f), which is Then, inserting Equations ( 5) and ( 7) into (2) gives: ( ) Experimental results by Shiono & Knight [8] demonstrate that in overbank flow the term of secondary flow, H(ρUV) d , varies almost linearly with y in the main and floodplain regions of a channel.Therefore, the lateral gradient of the secondary flow force may be approximated as where Γ is a dimensionless constant of secondary flow, which varies in different flow regions (Knight et al. [17]).Thus Equation ( 9) can be expressed by: which is the governing equation of depth-averaged flow in the streamwise direction.

Analytical Solution for d U 2
For given channels, when f, λ and Γ are fixed values, Equation ( 11) becomes a standard ordinary differential equation in terms of variable 2 d U .An analytical solution to Equation ( 11) can be obtained under appropriate boundary conditions.The detailed solution process is given by Knight et al. [17].The solutions for the channel with linearly varying bed can also refer to Knight & Shiono [18] and Tang & Knight [13] [19] for further details.
In a rectangular channel, which has the constant depth of flow (i.e.H = constant), the analytical solution of U d to Equation (11) has the following form: ( ) where ( ) where A 1 and A 2 are two integral constants, which can be obtained by applying the boundary conditions of flow continuity and the no-slip condition at the wall, as shown by Knight et al. [20] [21].For a simple rectangular channel of semi-width, b, see Figure 2, a single panel is used to represent half of the cross-section, and then the two A coefficients are obtained as: ( ) Thus for a rectangular channel, Equation ( 12) becomes ( ) where ( )

Calibration of Coefficients and Modelling Philosophy
Strictly speaking, Equation ( 15) suits for the flow region of constant depth where a constant Γ value exits.If Γ varies laterally, the cross section of a channel can be split into a discrete number of panels if necessary.Generally, panel junctions coincide with the place where the roughness of channel bed or flow depth is discontinuous.Thus, for a symmetric uniform rectangular channel, only half the channel needs to be modeled because of symmetry of the channel.If complex secondary flow cells present, then further division of the channel into more panels may be required [21].The friction coefficient is generally assumed to be constant in each panel, and it can be back calculated from ( ) Note that the variation of λ has great effect for open-channel flows due to lateral shear near wall [22].Hence, both the secondary flow parameter (Γ) and lateral eddy viscosity coefficient (λ) require calibration in modeling the lateral distribution of depth-averaged streamwise velocity.

Experimental Data in the Study
The data of lateral flow velocity used herein were obtained from Knight et al.
[14] and Atabay [15].These experiments have covered a wide range of aspect ratios (B/H), which vary from 1.0 to 16.0.The earlier experiments by Knight et al. [14] were conducted in a non-tilting 15 m long flume with a rectangular cross-section of width B =152 mm and a bed slope S o of 9.66 × 10 −4 .The more recent experiments by Atabay [15] were undertaken in a 22 m long flume with a channel width of 50 mm and a bed slope of 2.024 × 10 −3 .These data are also accessible from www.flowdata.bham.ac.uk.
In the present study, only half the channel was modeled due to the symmetry of flow, as shown in Figure 3, where only one panel is used for Equation (17).
The parameters adopted in the modeling are given in Table 1, where DWK and  AS denote the experimental runs by Knight et al. [14] and Atabay [15], respectively.In the modelling, ρ = 1000 kg/m 3 and g = 9.807 m/s 2 .

Prediction of Lateral Distributions of Ud
For given parameters of f, λ and Γ, Equation ( 17) can predict the lateral variation of depth-averaged velocity.Figure 4 and Figure 5 show the lateral distributions of U d between the prediction and experimental data.The results show that the prediction by Equation ( 17) agrees well with the experimental data for flows in rectangular open channels, which have a wide range of aspect ratios (B/H).
To evaluate the robustness of Equation ( 17 where E u,i is the error percentage of predicted velocity; u a,i and u e,i are the predicted and observed velocity at ith measured point (y), respectively.Therefore, the average error percentage of model for an experiment can be computed by ( ) where N is the total number of obervation in an experiment.
The values (E u ) of MAPE for all experiments are shown in Table 1, which shows that the prediction by Equation ( 17) has the averaged percentage error of 1% to 6%.
In this study, the optimum values of λ and Γ in are found to be well related to the aspect ratios of channel, as expressed by: ( ) where B is the width of channel, H is the flow depth, and R is the correlation coefficient.

Lateral Extent of Shear Layer
From the analytical solution of U d , given by Equation ( 17), it reveals that the depth-averaged velocity has little change in the major central part of channel, but it decreases rapidly toward zero at the wall, as shown in Figure 4 & Figure 5.
With increasing the aspect ratio (B/H) of channel, the extent of almost constant flow velocity is getting wider in the central channel, which can be treated as one-dimensional uniform flow.
When a channel is shallow and very wide, i.e. the ratio of B/H is very large, the value of cosh(γb) in Equation ( 18) can be very large (∞), so that the C valve in Equation ( 17) becomes a very small value (close to zero).Note that the secondary flow is weak in such a shallow channel, thus Γ is negligible (Knight et al. [21]).Therefore Equation ( 17) becomes: The obtained Equation ( 23) is the well-known velocity formula for one-dimensional ( ) Inserting Equation ( 18) into (24b) yields Since the parameter (γ), given by Equation ( 14), is a function of depth, Equation ( 25) is implicit for δ.However, δ can be obtained numerically by solving Equation (25), as graphically shown in Figure 6.This implies that when γb is about 5 the value of U d at y = b − δ is approximate to 0.99U ∞ .In other words, in the region away from the wall, where γδ = 5, the flow can be treated as uniform flow.
In order to estimate the shear layer width (δ) for a natural river, it needs to estimate the eddy viscosity coefficient (λ).When bed generated turbulence dominates the lateral momentum transfer, Cunge et al. [27] suggested that the typical value of λ is 0.5, which represents the value for the cross-sectional exchange coefficient in the range of 0.23 to 0.6.By taking the extreme case where λ = 0.6 and f = 0.01, Equation ( 26) gives δ/H ≈ 14.6.
In practical applications for open channel flows, within the boundary width (γδ = 5) Equation ( 17) can be used to predict the U d distribution in the lateral shear layer near the wall, whereas the velocity is calculated by Equation ( 23) in the central region of channel beyond the shear layer.In fact, the flow in the central part is uniform flow if the channel is very shallow (i.e.large aspect ratio B/H).
By taking the values of f and λ from Table 1, Equation (26) will establish a relationship between δ/H and B/H, as shown in Figure 7.The relationship can be described as:

Conclusions
The depth-averaged velocity distribution in lateral shear layer in a rectangular open channel can be described by Equation (17), which is analytically obtained based on the depth-averaged momentum equation proposed by Shiono & Knight.The analytical solution of Equation ( 17) agrees well with the experimental data with the aspect ratios ranging from 1 to 16.The percentage error in MAPE is 1% to 6%.
Based on Equation (17), the lateral width of the shear layer is found to be related to the dimensionless eddy viscosity coefficient (λ) and the bed friction (f), as given by Equation (26).The shear layer width (δ) near wall varies with the

Figure 1 .
Figure 1.Streamwise Reynolds stress in an open channel.
) against the experimental data, error analyses were carried out.MAPE (Mean Absolute Percentage Error) is used as a measure of the error percentage for predicted values of model against the measured values.The individual differences are called as residuals for the data sample that is used for estimation, and the residuals are known as estimation errors for the sample [23] [24] [25].The percentage of error in predicted discharge of each flow depth is calculated by, X. Tang DOI: 10.4236/jamp.2019.74056835 Journal of Applied Mathematics and Physics

Figure 4 .
Figure 4. U d distributions for the experiments by Knight et al. [14].

( 1 -
D) uniform flow, i.e.Darcy-Weisbach equation.In open channel flow, a lateral shear layer exists near the wall.The extent of lateral shear layer can be defined by a lateral distance, δ, away from the wall, for which ∞ is the velocity of 1-D flow, given by Equation (23), as illustrated in Figure3.The boundary distance δ in the shear layer can be oby = b − δ gives:

Table 1 .
Parameters adopted in the application of Equation (17).