Next Article in Journal
Use of Factor Analysis (FA), Artificial Neural Networks (ANNs), and Multiple Linear Regression (MLR) for Electrical Conductivity Prediction in Aquifers in the Gallikos River Basin, Northern Greece
Next Article in Special Issue
Hydrological Modeling in Agricultural Intensive Watershed: The Case of Upper East Fork White River, USA
Previous Article in Journal
Estimating the Precipitation Amount at Regional Scale Using a New Tool, Climate Analyzer
Previous Article in Special Issue
Simulation of Dam Breaks on Dry Bed Using Finite Volume Roe-TVD Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reinvestigating the Parabolic-Shaped Eddy Viscosity Profile for Free Surface Flows

ECAM-EPMI, LR2E-Lab, Quartz-Lab (EA 7393), 13 Boulevard de l’Hautil, 95092 Cergy-Pontoise, France
Hydrology 2021, 8(3), 126; https://doi.org/10.3390/hydrology8030126
Submission received: 19 July 2021 / Revised: 14 August 2021 / Accepted: 17 August 2021 / Published: 20 August 2021

Abstract

:
The flow in rivers is turbulent. The main parameter related to turbulence in rivers is the eddy viscosity, which is used to model a turbulent flow and is involved in the determination of both velocities and sediment concentrations. A well-known and largely used vertical distribution of eddy viscosity in free surface flows (open channels and rivers) is given by the parabolic profile that is based on the logarithmic velocity profile assumption and is valid therefore only in the log-law layer. It was improved thanks to the log-wake law velocity profile. These two eddy viscosities are obtained from velocity profiles, and the main shortcoming of the log-wake profile is the empirical Coles’ parameter. A more rigorous and reliable analytical eddy viscosity model is needed. In this study, we present two analytical eddy viscosity models based on the concepts of velocity and length scales, which are related to the exponentially decreasing turbulent kinetic energy (TKE) function and mixing length, namely, (1) the exponential-type profile of eddy viscosity and (2) an eddy viscosity based on an extension of von Karman’s similarity hypothesis. The eddy viscosity from the second model is Re * -independent, while the eddy viscosity from the first model is Re * -dependent (where Re * is the friction Reynolds number). The proposed analytical models were validated through computation of velocity profiles, obtained from the resolution of the momentum equation and comparisons to experimental data. With an additional correction function related to the damping effect of turbulence near the free surface, both models are similar to the log-wake-modified eddy viscosity profile but with different values of the Coles’ parameter, i.e., Π = 0.2 for the first model and Π = 0.15 for the second model. These values are similar to those found in open-channel flow experiments. This provides an explanation about the accuracy of these two analytical models in the outer part of free surface flows. For large values of Re * ( Re * > 2000), the first model becomes Re * independent, and the two coefficients reach asymptotic values. Finally, the two proposed eddy viscosity models are validated by experimental data of eddy viscosity.

1. Introduction

Determination of velocity distribution in open-channel flows and rivers is a topic of high interest and is involved in different practical applications [1,2,3,4,5,6,7]. However, the hydrodynamic in rivers and open-channel flows is strongly influenced by turbulence [8,9,10,11]. Different experimental studies were conducted to better understand the effect of turbulence on the streamwise velocity distribution [12,13,14,15,16].
Open-channel flow studies are considered an important preliminary step to investigate more complex river flows. The main interest in laboratory investigations of free surface turbulent flows is related to the experimental conditions that are chosen to be in agreement with the assumptions related to the models. Therefore, experiments in laboratory flumes allowed analytical models of turbulence to be developed that are in the form of analytical solutions, semi-theoretical or empirical relationships [17,18,19,20]. These models were proposed for mean velocities, turbulent kinetic energy (TKE), mixing length, eddy viscosity with their link to flow parameters (flow depth, friction velocity) [11,18]. The developed analytical models are mostly two-dimensional and for uniform flows, while free surface turbulent flows in rivers are generally strongly three-dimensional and non-uniform. Even if the assumptions are far from the real-life conditions in rivers, and the developed analytical models cannot account for the full complexity of turbulent river flows, these models present high practical interest. For example, measured data in the central part of rivers are well described by the analytical expressions obtained from laboratory open-channel flows [21].
In open-channel flows and rivers, in the classical two-layer approach, the flow is divided into two regions: an inner ( ξ < 0.2 ) and an outer ( ξ > 0.2 ) region (where ξ = y / h is the ratio of the distance from the bed y to flow depth h ). For smooth open-channel flows, log law allows accurate description of mean streamwise velocities U + ( y + ) in the logarithmic layer ( 30 < y + < 0.2 R e * or 30 / R e * < ξ < 0.2 , where in wall unit y + = y   u * / ν ,   u * is the friction velocity, ν the kinetic viscosity, and R e * = h   u * / ν is the friction Reynolds number) [11]. For river flows, Franca and Lemmin [22] found from analysis of experimental data from a field study of extremely rough, three-dimensional river flows, that in more than 65% of the profiles, the log law can be applied up to ξ = 0.4 , while above this value, mean velocities show deviations from the logarithmic profile.
In the outer region, in addition to the simple power law [23,24], the log-wake law is largely used. It is an extension of the log law by adding the Coles’ wake function, which contains the Coles’ parameter Π [25,26]. However, this method is empirical and no physical-based approach is available to determine the Coles’ parameter Π [27]. The value of Π is therefore not universal. From experiments in zero-pressure-gradient boundary layers, Cebeci and Smith [28] found that Π increases with the Reynolds number and becomes Π = 0.55 at high Reynolds numbers. For open-channel flows over smooth beds, the following values were found: Π increases from zero with Re * and becomes Π ≈ 0.2 for Re * = 2000 [18], Π ≈ 0.08 [29], Π = 0.1 [30], Π = 0.3 [31], Π = 0.45 for velocity distribution with dip-phenomenon in narrow open channels [32].
However, the more rigorous method for the prediction of velocity profiles is based on the resolution of the momentum equation. This method needs a reliable model for the eddy viscosity, which is the main parameter related to turbulence used in free surface flows since it is involved in the determination of both velocities and sediment concentrations (through the sediment diffusivity, i.e., the product of the eddy viscosity by the inverse of the turbulent Schmidt number) [33,34]. A well-known and largely used vertical distribution of eddy viscosity in free surface flows (open channels and rivers) is given by the parabolic profile [11]. This profile is based on the logarithmic velocity profile assumption and is valid therefore only in the log-law layer. It was improved thanks to the log-wake law velocity profile. These two eddy viscosities are obtained from velocity profiles. As for velocity profiles in the outer region, the main shortcoming of the log-wake eddy viscosity profile is the non-universal Coles’ parameter Π.
An analytical eddy viscosity model is therefore needed to predict velocity profiles. In this study, we present two eddy viscosity models based on the concepts of velocity and length scales, which are related, respectively, to the exponentially decreasing turbulent kinetic energy (TKE) function [11] and mixing length, namely, (1) the exponential-type profile of eddy viscosity [35,36] and (2) an eddy viscosity based on an extension of von Karman’s similarity hypothesis [37,38,39]. An additional correction is used in order to account for the damping effect of turbulence near the free surface. The proposed analytical models are validated through computation of velocity profiles, obtained from the resolution of the momentum equation, and comparisons to experimental data. This study aims to provide an explanation and a theoretical foundation to the empirical well-known eddy viscosity profiles.

2. Literature for Eddy Viscosity Models for Open-Channel Flows

2.1. Parabolic Eddy Viscosity

The widely used eddy viscosity ( ν t ) formulation is the parabolic profile given by [40,41]
ν t ( y ) = κ u * y ( 1 y h ) ,
where y is the vertical distance from the bed, κ is the von Karman constant, h the flow depth, and u * the friction or shear velocity. Equation (1) is based on shear stress, which decreases linearly with distance from channel bed y and a logarithmic velocity profile [33].

2.2. Log-Wake-Modified Eddy Viscosity Profile

When used for the resolution of the momentum equation, the parabolic eddy viscosity profile (1) is unable to predict accurately velocities outside the log layer [36]. In order to improve the parabolic eddy viscosity (Equation (1)), it is corrected in accordance with Coles’ log-wake law for velocities as [18]
ν t ( y ) = κ u * y ( 1 y h ) 1 + π   Π   y h sin ( π y h )
where Π is the Coles’ parameter. In Equation (2), the eddy viscosity (1) is corrected by dividing the parabolic eddy viscosity profile by the term 1 + π Π ( y / h ) sin ( π y / h ) of the log-wake velocity profile.

2.3. Mixing Length and Mixing Velocity

In order to predict the velocity profile over the entire flow depth, it is more suitable to define the eddy viscosity from the concepts of velocity and length scales, which are here given by mixing velocity and mixing length as
ν t = w m l m
From the parabolic profile given by Equation (1), mixing length and mixing velocity are given, respectively, by l m = κ y 1 y / h and w m = u * 1 y / h . The eddy viscosity of Equation (2) allows for the following expression:
l m ( y ) = κ y 1 y / h 1 + π Π ( y / h ) sin ( π y / h )
Equation (4) consists of a correction of the “parabolic” mixing length κ y 1 y / h by dividing it by the term 1 + π Π ( y / h ) sin ( π y / h ) as in Equation (2). The related mixing velocity is given by w m = u * 1 y / h [42].

3. Proposed Eddy Viscosity Models for Free Surface Flows

The eddy viscosity is related to turbulent kinetic energy (TKE) as
ν t = C µ 1 4 k l m
where k is the turbulent kinetic energy (TKE), and C µ = 0.09 .

3.1. Mixing Velocity from TKE Profile

A semi-theoretical function for TKE is given by [11]:
  k ( ξ ) = D k u * 2 e 2 C k ξ
where ξ = y / h , C k and D k are empirical constants, D k = 4.78 , and C k = 1 [11]. Equation (6) was validated by direct numerical simulation DNS data [43]. Instead of the mixing velocity given in the above section by w m = u * 1 y / h , the shape of mixing velocity should be supported by turbulence intensity measurements [38] and in agreement with the TKE formulation of Equation (6).
From Equations (3) and (5), the mixing velocity is related to TKE as follows:
w m = C µ 1 4 k = C µ 1 4 D k u * e C k ξ
With the assumption D k = 1 / C µ 1 4 (based on log-law and local equilibrium assumption [39]) and with C k = 1 , the mixing velocity reverts to
w m ( ξ ) = u * e ξ ,
which shows that the mixing velocity decreases exponentially with distance from the bed and is in agreement with observations of turbulence intensity and TKE.

3.2. Damping Function for Free Surface

An additional correction is required in order to account for the damping effect of turbulence near the free surface [44,45,46]. In order to decrease turbulent viscosity near the free surface, Hosoda [47] proposed a damping function as
f ( ξ ) = 1 e B f ( 1 ξ )
where B f is a damping coefficient.
If we include the free surface damping function in Equation (3), the eddy viscosity reverts to
ν t = w m l m f ( ξ )
In the following sections, two eddy viscosity formulations will be presented.

3.3. First Formulation: Exponential-Type Profile of Eddy Viscosity

In the equilibrium region, where TKE production is balanced by dissipation, the velocity gradient is given by d U d y = C μ 1 / 4 k   l m . In the log-law layer, d U d y =   u * κ   y , and with a TKE given by k = C μ 1 / 4 u * f ( y ) , the mixing length should read as l m = κ y f ( y ) [38]. Since in the equilibrium region ( y + > 50 ), TKE is given by Equation (6), and the velocity profile is given by the log law for ( y + > 30 ). These two conditions that are given by d U d y k   l m =   u * κ   y show that the mixing length should be as l m = κ   y   e C k ( y / h ) . The 1st eddy viscosity is therefore given by
ν t ( y ) = α 1 κ u * y e C 1 ξ
where α 1 and C 1 are two coefficients, α 1 is related to C μ and D k , while C 1 is related to C k . Equation (11), i.e., the exponential-type profile of eddy viscosity is consistent with the exponentially decreasing mixing velocity (8). Equation (11) was proposed empirically and was used in the planetary boundary layer [48] and coastal engineering [49,50,51,52,53]. However, in order to allow accurate description for different flow conditions, Equation (11) was written in a Re * -dependent form as [35,36]
ν t ( y ) = u * y e y + + 0.34 Re * 11.5 0.46 Re * 5.98  
where in the wall units, y + = y u * / ν , and the friction Reynolds number Re * = h u * / ν . In other words,
ν t + ( y + ) = y + e y + + 0.34 Re * 11.5 0.46 Re * 5.98
where ν t + = ν t / ν . The link between Equations (11) and (12) is given by
C α = α 1 κ = e 0.34 Re * 11.5 0.46 Re * 5.98   and   C 1 = Re * 0.46 Re * 5.98

3.4. Second Formulation: Eddy Viscosity Formulation Based on Miwing Length Equation from Similarity Hypothesis

This formulation uses a mixing length, which was derived from an extension of von Karman’s similarity hypothesis, energy equilibrium assumption, and Nezu and Nakagawa’s (1993) TKE function as [37,38,39,54] (Appendix A)
l m ( ξ ) = κ h ( 1 e ξ )
Using Equations (8)–(10) and (13), the second proposed eddy viscosity is given by
ν t ( ξ ) = κ h u * e ξ ( 1 e ξ ) ( 1 e B f ( 1 ξ ) )

4. Results

The following ordinary differential equation for velocity distribution U in open-channel flows was obtained from analysis of the Reynolds-averaged Navier–Stokes equations [32,55].
d U d y = u τ 2 ν + ν t [ ( 1 y h ) α y h ]
where α is a parameter related to dip-phenomenon. For wide-open channels (ratio of channel width to flow depth > 5 ) α = 0 , and Equation (15) reverts in the outer region (in the wall units) to the following:
d U + d y + = 1 ν t + ( 1 y + Re * )
Mean streamwise velocities are obtained from the numerical resolution of Equation (16). To solve Equation (16), the eddy viscosity ν t + is calculated using the two proposed analytical eddy viscosity models given above, i.e., the first is given by Equation (12), while the second by Equation (14).
Both proposed models are validated by experimental data of velocities in open-channel flows for 923 < Re * < 6139 [18]. The measurements were carried out in a rectangular cross section and a hydraulically smooth wall open channel. The total length of the channel is 20 m with a cross-sectional size (60 cm wide × 65 cm deep). The width of the channel is sufficiently large to neglect the wall effect. A flow depth of 10 cm was kept constant with varying discharge to examine at various Froude and Reynold numbers. Moreover, measurement was carried out at a length of 18 m from the inlet position. Laser Doppler Anemometer (LDA) was used to carry out velocity measurements under different flow conditions. Table 1 summarizes experimental hydraulic conditions. The flow conditions are listed in Table 1 [18].

4.1. Velocity Profiles from the First Eddy Viscosity Formulation: Exponential-Type Profile

Computed mean streamwise velocity profiles are obtained from (16) with the first eddy viscosity given by Equation (12) and are validated by experimental data [18]. The following boundary condition for the velocity is applied at ξ = 0.2 (or y + = 0.2   R e * ): U + ( y + = 0.2 Re * ) = ( 1 / κ ) ln ( 0.2 Re * ) + B ; where κ = 0.41 and B = 5.29 [11].
Figure 1 and Figure 2 show comparisons between computed velocity profiles (red solid lines) and experimental data (symbols). The model allows the prediction of log law (black thin dashed lines) to be improved in the outer region.
In order to improve the results, the first eddy viscosity given by (12) is used with the condition of an eddy viscosity equal to zero at the free surface, which requires the use of the damping function given by Equation (9), Equation (12) reverts to
ν t + = y + e y + + 0.34 Re * 11.5 0.46 Re * 5.98 ( 1 e B f ( 1 y + Re * ) )
Results obtained with Equation (17) (magenta thick dashed curves) allow the prediction to be improved, particularly for high Reynolds numbers.

4.2. Velocity Profiles from the Second Eddy Viscosity Formulation Based on Von Karman’s Similarity Hypothesis

Computed mean streamwise velocity profiles are obtained from (16), and the second eddy viscosity given by (14) and are validated by experimental data [18] (Table 1). The same boundary condition is applied (velocity equal to the log-law value at y + = 0.2   R e * ). Figure 3 shows comparisons between computed velocity profiles (red solid lines) and experimental data (symbols). The model allows log-law profiles (black thin dashed lines) to be improved in the outer region. Figure 3 shows that computed velocity profiles (red solid lines) show good agreement and improve the prediction of log law (black dashed thin lines). In order to improve the results, the momentum equation is resolved with a second boundary condition at the lower limit of the logarithmic layer, i.e., y + = 30 , where the velocity is given by the log law as U + ( 30 ) = ( 1 / κ ) ln ( 30 ) + B . Results (magenta thick dashed curves) allow the prediction to be improved, particularly for high Reynolds numbers.

4.3. Eddy Viscosity Profiles

In order to compare the two proposed analytical eddy viscosity models (Equations (12), (14) and (17)) with the existing parabolic and log-wake-modified eddy viscosity profiles (Equations (1) and (2)), all equations are written in the same dimensionless form as
ν t   h   u * = κ y h ( 1 y h )
ν t   h   u * = κ y h ( 1 y h ) 1 + π Π y h sin ( π y h )
ν t   h   u * = ξ e ( ξ + 0.34 ) R e * 11.5 0.46 R e * 5.98
ν t   h   u * = ξ e ( ξ + 0.34 ) R e * 11.5 0.46 R e * 5.98 ( 1 e B f ( 1 ξ ) )
ν t   h   u * = κ e ξ ( 1 e ξ ) ( 1 e B f ( 1 ξ ) )
We notice that eddy viscosities from Equations (18), (19), and (22) are Re * -independent, while Equations (20) and (21) are Re * dependent.
Figure 4 and Figure 5 show comparisons of the vertical distribution of the different eddy viscosity models. The three R e * -independent Equations (18), (19), and (22) are first compared. Figure 4 shows that the shape of the eddy viscosity given by the second model (Equation (22)) (red solid line) is similar to the parabolic profile (green dashed line), where the maximum value is located at the same position, i.e., half water depth ( ξ = 0.5 ) . Even though the second model exhibits a similar shape, it predicts smaller values than the parabolic profile. The profile obtained from the second model is compared to log-wake-modified (dash-dotted lines) profiles. The magenta dash-dotted curve is from the log-wake-modified eddy viscosity given by Equation (19) with a Coles’ parameter Π = 0.2 . With a smaller value of Π , the blue dash-dotted curve (for Π = 0.15 ) is closer to the eddy viscosity given by the second model. This value ( Π = 0.15 ) is close to values found for open-channel flow experiments [18,30]. This provides an explanation about the accuracy of the computed velocity profiles obtained by the second model in the outer part of free surface flows.
For the Re * -dependent eddy viscosity given by the first model (Equations (20) and (21)), Figure 5 shows comparisons for the vertical distribution of the eddy viscosity for the four friction Reynolds number. The eddy viscosity curves (thick red solid lines) obtained from the first model (Equation (20)) show that the eddy viscosity increases from the bed—the maximum value is located at around the half water depth and then decreases until the free surface. However, the eddy viscosity does not vanish at the free surface. With the condition of an eddy viscosity equal to zero at the free surface, the profiles with the damping function (Equation (21)) predict everywhere smaller values than the parabolic eddy viscosity (thin magenta solid lines) and have a shape similar to the log-wake-modified eddy viscosity profile (blue dash-dotted lines). Interestingly, with a value of Coles’ parameter Π = 0.2 , eddy viscosities from both the first model (Equation (21)) and log-wake-modified profile (Equation (19)) are almost superimposed. The value Π = 0.2 of Coles’ parameter is the same as that proposed from open-channel flow experiments [18]. This provides an explanation about the accuracy of the computed velocity profiles obtained by the first model (Equation (21)) in the outer part of free surface flows.
The profiles obtained for the four R e * . numbers seem similar. For large values of Re * , the two coefficients C α = α 1 κ and C 1 . of the first model reach asymptotic values equal, respectively, to C α = α 1 κ = e ( 0.34 / 0.46 ) = 0.477 and C 1 = 1 / 0.46 = 2.17 (Figure 6). Equations (20) and (21) reverts to the following Re * -independent forms:
ν t   h   u * = C α ξ e C 1 ξ
ν t   h   u * = C α ξ ( 1 e B f ( 1 ξ ) ) e C 1 ξ
where C α = 0.477   , C 1 = 2.17 , and B f = 6 .
Figure 7 shows the comparison between the two proposed eddy viscosity models (24, 22), both with free surface damping function and parabolic and log-wake-modified eddy viscosity profiles. Figure 7 shows that with the free surface damping function both exhibit smaller values than the parabolic profile. Both models predict profiles similar to log-wake-modified eddy viscosity (Equation (19)). The first model (Equation (24)) (red solid line) is similar to (Equation (19))with a Coles’ parameter Π = 0.2 (red dash-dotted line), while the second model (Equation (22)) (blue solid line) is similar to (Equation (19))with a Coles’ parameter Π = 0.15 (blue dash-dotted line). Interestingly, the two coefficients of the first model are found to be equal, respectively, to C α = α 1 κ = 0.477 and C 1 = 2.17 . The value of the first coefficient C α results in α 1 = 1.16 (with κ = 0.41 ), which is close to α 1 = 1 (related to the assumption D k = 1 / C µ 1 4 ) (see also [49,50]). The value of C α is also between the two values 0.41 < 0.478 < 0.49 obtained, respectively, with the assumption D k = 1 / C µ 1 4 , D k C µ 1 4 κ = κ = 0.41 and the empirical value D k = 4.78 [56], D k C µ 1 4 κ = 0.49 . The value of the second coefficient allows the coefficient in TKE to be defined as C k = C 1 / 2 = 1.088 , which is close to the empirical value C k = 1 [55].
Figure 8 shows a comparison between the two proposed eddy viscosity models (24, 22) and experimental data of eddy viscosity [18,57]. In addition to the experimental data of Nezu and Rodi [18], data from experiments of Ueda et al. [57] are used, which seem to confirm the same behavior. Figure 8 shows that profiles obtained from both models show good agreement with experimental data of eddy viscosity.

5. Conclusions

The parabolic eddy viscosity is based on the log-law velocity profile and is valid only in the log layer. The improved log-wake-modified eddy viscosity was obtained from the log-wake law velocity profile. However, both were obtained from velocity profiles. The main shortcoming of the log-wake profile is the uncertainty in the value of the empirical Coles’ parameter.
In this study, the eddy viscosity is defined as a product between a velocity scale (related to the root-square of TKE, which is given by a semi-theoretical exponentially decreasing function) and a length scale (related to mixing length). From this definition, two analytical eddy viscosity models are proposed, namely, (1) the exponential-type profile of eddy viscosity and (2) an eddy viscosity based on an extension of von Karman’s similarity hypothesis with an additional correction in order to account for the damping effect of turbulence near the free surface. As for the parabolic and log-wake-modified profiles, the eddy viscosity from the second model is Re * independent, while the eddy viscosity from the first model is Re * dependent. The proposed analytical models are validated through computation of velocity profiles, obtained from the resolution of the momentum equation and comparisons to experimental data.
Mean streamwise velocity profiles were obtained by solving the momentum equation. For both the first and second proposed analytical eddy viscosity models, a boundary condition is applied at the lower limit of the outer region, i.e., y + = 0.2   Re * or ξ = 0.2 , where the velocity is given by the logarithmic law. Computed mean velocities are compared to experimental data of open-channel flows for 923 < Re * < 6139 . Computed velocity profiles show good agreement in the outer region.
In order to improve the results, the first eddy viscosity was used with the free surface damping function. Results allow the prediction to be improved, particularly for high Reynolds numbers. For the second analytical eddy viscosity model, a second boundary condition was used at the lower limit of the logarithmic layer, i.e., y + = 30 where the velocity is given by the log law. Results allow predicted velocity profiles to be improved, particularly for high Reynolds numbers. The results for velocity profiles show the ability of these analytical eddy viscosity models to predict accurately the velocities in the outer region from the momentum equation.
Finally, the vertical distribution of eddy viscosity from both proposed analytical models was analyzed. Both profiles from the first and second analytical eddy viscosity models with the free surface damping function are similar to the log-wake-modified profiles but with different values of Coles’ parameter. Π = 0.2 for the first model and Π = 0.15 for the second model. These values are close to values found from open-channel flow experiments. This provides an explanation of the accuracy of the computed velocity profiles in the outer part of free surface flows.
For large values of Re * ( Re * > 2000), the first model becomes Re * independent, and the two coefficients reach asymptotic values equal to C α = α 1 κ = 0.477 and C 1 = 2.17 . Interestingly, with a value of Coles’ parameter Π = 0.2 , eddy viscosities from both the first model and log-wake-modified profile are almost superimposed. The value Π = 0.2 of Coles’ parameter is the same as that proposed from open-channel flow experiments. The analysis of these two coefficients allowed the models’ assumptions to be verified and former empirical values to be found. The two proposed eddy viscosity models are validated by two experimental data. The comparison shows that profiles from both models show good agreement with experimental data of eddy viscosity.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author would like to thank Xiaofeng Liu from Penn State University for providing the experimental data of eddy viscosity of references [18,57] used in Figure 8.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A

The von Karman’s similarity hypothesis allows writing the mixing length as [58]
l m + = κ d U + / d y + d 2 U + / d y + 2
where U + is the streamwise mean velocity. With d U + / d y + k + / l m + , Equation (A1) becomes
l m + = κ k + / l m + d ( k + / l m + ) / d y +
Introducing the function f + = k + / l m + , (A2) becomes
k + = κ f + 2 d f + / d y +
We write (A3) in the following form:
d f + / d y + f + 2 = κ k +
The integration of the LHS term of Equation (A4) from A 0 to y + ¸ is given by
A 0 y + d f + / d y + f + 2 d y + = 1 f + ( y + ) 1 f + ( A 0 )
Integrating (A4) from A 0 to y + ¸ provides, therefore, the mixing length as
l m + ( y + ) = k + ( κ A 0 y + 1 k + d y + + l m + ( A 0 ) k + ( A 0 ) )
Using Equation (6) for TKE and taking the boundary condition l m + ( A 0 ) = κ A 0 yields
l m + ( y + ) = κ e C k y + / R e * ( R e * C k ( e C k y + / R e * e C k A 0 / R e * ) + A 0 e C k A 0 / R e * )
Rearranging the terms of (A7) allows writing the mixing length as [37,38,39,54]
l m + ( y + ) = α ( R e * ( R e * C k A 0 ) e C k ( y + A 0 ) / R e * )
where α = κ / C k . Since A 0 R e * we write (18) in a simplified form that does not depend on A 0 , as
l m + ( y + ) = α   R e * ( 1 e C k y + / R e * )
By taking C k = 1 , Equation (A9) reverts to Equation (13).

References

  1. Figuérez, J.; Galán, Á.; González, J. An Enhanced Treatment of Boundary Conditions for 2D RANS Streamwise Velocity Models in Open Channel Flow. Water 2021, 13, 1001. [Google Scholar] [CrossRef]
  2. Nazari-Sharabian, M.; Karakouzian, M.; Hayes, D. Flow Topology in the Confluence of an Open Channel with Lateral Drainage Pipe. Hydrology 2020, 7, 57. [Google Scholar] [CrossRef]
  3. Zhu, Z.; Yu, J.; Dou, J.; Peng, D. An Expression for Velocity Lag in Sediment-Laden Open-Channel Flows Based on Tsallis Entropy Together with the Principle of Maximum Entropy. Entropy 2019, 21, 522. [Google Scholar] [CrossRef] [Green Version]
  4. MacVicar, B.J.; Roy, A.G. Hydrodynamics of a forced riffle pool in a gravel bed river: 1. Mean velocity and turbulence intensity. Water Resour. Res. 2007, 43, W12401. [Google Scholar] [CrossRef] [Green Version]
  5. Liu, X.; Nayamatullah, M. Semianalytical Solutions for One-Dimensional Unsteady Nonequilibrium Suspended Sediment Transport in Channels with Arbitrary Eddy Viscosity Distributions and Realistic Boundary Conditions. J. Hydraul. Eng. 2014, 140, 04014011. [Google Scholar] [CrossRef]
  6. Johnson, E.D.; Cowen, E.A. Estimating bed shear stress from remotely measured surface turbulent dissipation fields in open channel flows. Water Resour. Res. 2017, 53, 1982–1996. [Google Scholar] [CrossRef]
  7. Macvicar, B.J.; Sukhodolov, A.N. Sampling strategies to improve scaling parameter estimates in rivers. J. Hydraul. Res. 2019, 57, 1–10. [Google Scholar] [CrossRef]
  8. Sukhodolov, A.N.; Krick, J.; Sukhodolova, T.A.; Cheng, Z.; Rhoads, B.L.; Constantinescu, G.S. Turbulent flow structure at a discordant river confluence: Asymmetric jet dynamics with implications for channel morphology. J. Geophys. Res. Earth Surf. 2017, 122, 1278–1293. [Google Scholar] [CrossRef]
  9. Sukhodolov, A.N. Field-based research in fluvial hydraulics: Potential, paradigms and challenges. J. Hydraul. Res. 2014, 53, 1–19. [Google Scholar] [CrossRef]
  10. Gualtieri, C.; Shao, D.; Angeloudis, A. Advances in Environmental Hydraulics. Water 2021, 13, 1192. [Google Scholar] [CrossRef]
  11. Nezu, I.; Nakagawa, H. Turbulence in Open Channel Flows; Balkema, A.A., Ed.; Balkema: Rotterdam, The Netherlands, 1993. [Google Scholar]
  12. Yokosi, S. The structure of river turbulence. Bull. Disaster Prev. Res. Inst. Kyoto Univ. 1967, 17, 1–29. [Google Scholar]
  13. McQuivey, R.S. Summary of turbulence data from rivers, conveyance channels and laboratory flumes. US Geol. Surv. Prof. 1973, N802B, B1–B66. [Google Scholar]
  14. Grinvald, D.I. Turbulence of Natural Flows (in Russian); Hydrometeoizdat: Leningrad, Russia, 1974. [Google Scholar]
  15. Iwasa, Y.; Asano, T. Characteristics of turbulence in rivers and conveyance channels. In Proceedings of the Third International Symposium on Stochastic Hydraulics, Delft, The Netherlands, 5–8 August 1980; pp. 565–576. [Google Scholar]
  16. Nikora, V.I.; Smart, G.M. Turbulence characteristics of New Zealand gravel-bed rivers. J. Hydraul. Eng. 1997, 123, 764–773. [Google Scholar] [CrossRef]
  17. Grass, A.J. Structural features of turbulent flow over smooth and rough boundaries. J. Fluid Mech. 1971, 50, 233–255. [Google Scholar] [CrossRef]
  18. Nezu, I.; Rodi, W. Open channel measurements with a laser Doppler anemometer. J. Hyd. Eng. ASCE 1986, 112, 335–355. [Google Scholar] [CrossRef]
  19. Absi, R. A simple eddy viscosity formulation for turbulent boundary layers near smooth walls. C. R. Mec. 2009, 337, 158–165. [Google Scholar] [CrossRef] [Green Version]
  20. Guo, J. Modified log-wake-law for smooth rectangular open channel flow. J. Hydraul. Res. 2014, 52, 121–128. [Google Scholar] [CrossRef]
  21. Sukhodolov, A.; Thiele, M.; Bungartz, H. Turbulence Structure in a River Reach with Sand Bed. Water Resour. Res. 1998, 34, 1317–1334. [Google Scholar] [CrossRef]
  22. Franca, M.J.; Lemmin, U. A field study of extremely rough, three-dimensional river flow. In Proceedings of the 4th IAHR International Symposium on Environmental Hydraulics, Hong Kong, China, 15–18 December 2004. [Google Scholar]
  23. Afzal, N.; Seena, A.; Bushra, A. Power Law Velocity Profile in Fully Developed Turbulent Pipe and Channel Flows. J. Hydraul. Eng. ASCE 2007, 133, 1080–1086. [Google Scholar] [CrossRef]
  24. Cheng, N.S. Power-law index for velocity profiles in open channel flows. Adv. Water Res. 2007, 30, 1775–1784. [Google Scholar] [CrossRef]
  25. Coles, D.E. The law of the wake in the turbulent boundary layer. J. Fluid Mech. 1956, 1, 191–226. [Google Scholar] [CrossRef] [Green Version]
  26. Hinze, J.O. Turbulence; McGraw-Hill: New York, NY, USA, 1975. [Google Scholar]
  27. Krug, D.; Philip, J.; Marusic, I. Revisiting the law of the wake in wall turbulence. J. Fluid Mech. 2017, 811, 421–435. [Google Scholar] [CrossRef] [Green Version]
  28. Cebeci, T.; Smith, M.O. Analysis of Turbulent Boundary Layers; Academic: San Diego, CA, USA, 1974. [Google Scholar]
  29. Cardoso, A.H.; Graf, W.H.; Gust, G. Uniform flow in a smooth open channel. J. Hydraul. Res. 1989, 27, 603–616. [Google Scholar] [CrossRef]
  30. Kirkgoz, M.S. The turbulent velocity profiles for smooth and rough channel flow. J. Hydraul. Eng. 1989, 115, 1543–1561. [Google Scholar] [CrossRef]
  31. Li, X.; Dong, Z.; Chen, C. Turbulent flows in smooth-wall open channels with different slope. J. Hydraul. Res. 1995, 33, 333–347. [Google Scholar]
  32. Absi, R. An ordinary differential equation for velocity distribution and dip-phenomenon in open channel flows. J. Hydraul. Res. 2011, 49, 82–89. [Google Scholar] [CrossRef]
  33. Absi, R. Rebuttal on A mathematical model on depth-averaged β-factor in open-channel turbulent flow. Environ. Earth Sci. 2020, 79, 113. [Google Scholar] [CrossRef]
  34. Gualtieri, C.; Angeloudis, A.; Bombardelli, F.; Jha, S.; Stoesser, T. On the Values for the Turbulent Schmidt Number in Environmental Flows. Fluids 2017, 2, 17. [Google Scholar] [CrossRef] [Green Version]
  35. Absi, R. Eddy viscosity and velocity profiles in fully-developed turbulent channel flows. Fluid Dyn. 2019, 54, 137–147. [Google Scholar] [CrossRef]
  36. Absi, R. Analytical eddy viscosity model for velocity profiles in the outer part of closed- and open-channel flows. Fluid Dyn. 2021, 56, 577–586. [Google Scholar] [CrossRef]
  37. Absi, R. Wave boundary layer instability near flow reversal. In Coastal Engineering 2002, Proceedings of the 28th International Conference Coastal Engineering, ASCE, Wales, UK, 19 November 2002; Smith, J.M., Ed.; World Scientific Publishing: Singapore, 2002; Volume 1, pp. 532–544. [Google Scholar]
  38. Absi, R. Comment on Turbulent diffusion of momentum and suspended particles: A finite-mixing-length theory. Phys. Fluids 2005, 17, 079101. [Google Scholar] [CrossRef] [Green Version]
  39. Absi, R. A roughness and time dependent mixing length equation. JSCE Doboku Gakkai Ronbunshuu B Jpn. Soc. Civ. Eng. 2006, 62, 437–446. [Google Scholar] [CrossRef] [Green Version]
  40. Yalin, M.S. Mechanics of Sediment Transport; Pergamon Press: Oxford, UK, 1977. [Google Scholar]
  41. Absi, R. Comments on Turbulent velocity profile in fully-developed open channel flows. Environ. Fluid Mech. 2008, 8, 389–394. [Google Scholar] [CrossRef]
  42. Absi, R. Comments on Modeling stratified suspension concentration distribution in turbulent flow using fractional advection–diffusion equation. Environ. Fluid Mech. 2021. Submitted. [Google Scholar]
  43. Absi, R. Analytical solutions for the modeled k-equation. ASME J. Appl. Mech. 2008, 75, 044501. [Google Scholar] [CrossRef]
  44. Nakagawa, H.; Nezu, I.; Ueda, H. Turbulence of open channel flow over smooth and rough beds. Jpn. Soc. Civ. Eng. 1975, 241, 155–168. [Google Scholar] [CrossRef]
  45. Langhi, M.; Hosoda, T.; Dey, S. Analytical Solution of k-ϵ Model for Nonuniform Flows. J. Hydraul. Eng. ASCE 2018, 144, 04018033. [Google Scholar] [CrossRef]
  46. Welderufael, M.; Absi, R.; Mélinge, Y. Assessment of velocity profile models for turbulent smooth wall open channel flows. ISH J. Hyd. Eng. 2021, in press. [Google Scholar]
  47. Hosoda, T. Turbulent Diffusion Mechanism in Open Channel Flows. Ph.D. Thesis, Kyoto University, Kyoto, Japan, 1990. [Google Scholar]
  48. Businger, J.A.; Arya, S.P.S. Heights of the mixed layer in the stable stratified planetary boundary layer. Adv. Geophys. 1974, 18A, 73–92. [Google Scholar]
  49. Hsu, T.W.; Jan, C.D. Calibration of Businger-Arya type of eddy viscosity model’s parameters. J. Waterw. Port Coastal. Ocean. Eng. ASCE 1998, 124, 281–284. [Google Scholar] [CrossRef]
  50. Absi, R. Discussion of Calibration of Businger-Arya type of eddy viscosity model’s parameters. J. Waterw. Port Coastal. Ocean. Eng. ASCE 2000, 126, 108–109. [Google Scholar] [CrossRef]
  51. Absi, R. Time-Dependent Eddy Viscosity Models for Wave Boundary Layers. In Coastal Engineering 2000, Proceedings of the 27th International Conference Coastal Engineering, Sydney, Australia, 16–21 July 2000; Edge, B.L., Ed.; ASCE: Reston, VA, USA; pp. 1268–1281.
  52. Absi, R. Concentration profiles for fine and coarse sediments suspended by waves over ripples: An analytical study with the 1-DV gradient diffusion model. Adv. Water Res. 2010, 33, 411–418. [Google Scholar] [CrossRef] [Green Version]
  53. Absi, R.; Tanaka, H.; Kerlidou, L.; André, A. Eddy viscosity profiles for wave boundary layers: Validation and calibration by a k-ω model. In Proceedings of the 33th International Conference Coastal Engineering, Santander, Spain, 1–6 July 2012. [Google Scholar]
  54. Absi, R. Stress length and momentum cascade in turbulent channel flows. Phys. Rev. Fluids 2021. Submitted. [Google Scholar]
  55. Yang, S.Q.; Tan, S.K.; Lim, S.Y. Velocity distribution and dip-phenomenon in smooth uniform open channel flows. J. Hydraul. Eng. 2004, 130, 1179–1186. [Google Scholar] [CrossRef]
  56. Nezu, I. Open-Channel Flow Turbulence and Its Research Prospect in the 21st Century. J. Hydraul. Eng. ASCE 2005, 131, 229–246. [Google Scholar] [CrossRef]
  57. Ueda, H.; Moller, R.; Komori, S.; Mizushina, T. Eddy diffusivity near the free surface of open channel flow. Int. J. Heat Mass Transf. 1977, 20, 1127–1136. [Google Scholar] [CrossRef]
  58. von Kármán, T. Mechanische ahnlichkeit und turbulenz. Math. -Phys. Kl. 1930, 5, 58–76. [Google Scholar]
Figure 1. Velocity profiles (red solid lines) obtained from (16) with the first eddy viscosity (12); symbols: experimental data [18]; dashed lines: log law (profiles shifted by 5 units).
Figure 1. Velocity profiles (red solid lines) obtained from (16) with the first eddy viscosity (12); symbols: experimental data [18]; dashed lines: log law (profiles shifted by 5 units).
Hydrology 08 00126 g001
Figure 2. Velocity profiles obtained from (16) with the first eddy viscosity (12) (red solid lines) and from (16) with (17) ( B f = 6 ) magenta thick dashed lines; symbols: experimental data [18]; dashed lines: log law (profiles shifted by 5 units).
Figure 2. Velocity profiles obtained from (16) with the first eddy viscosity (12) (red solid lines) and from (16) with (17) ( B f = 6 ) magenta thick dashed lines; symbols: experimental data [18]; dashed lines: log law (profiles shifted by 5 units).
Hydrology 08 00126 g002
Figure 3. Velocity profiles obtained from (16) with the second eddy viscosity (14) ( B f = 4 ) with the first boundary condition at ξ = 0.2 (red solid lines), and the second boundary condition at y + = 30 (magenta thick dashed lines); symbols: experimental data [18]; dashed lines: log law (profiles shifted by 5 units).
Figure 3. Velocity profiles obtained from (16) with the second eddy viscosity (14) ( B f = 4 ) with the first boundary condition at ξ = 0.2 (red solid lines), and the second boundary condition at y + = 30 (magenta thick dashed lines); symbols: experimental data [18]; dashed lines: log law (profiles shifted by 5 units).
Hydrology 08 00126 g003
Figure 4. Vertical distribution of eddy viscosity; green dashed line: parabolic eddy viscosity (Equation (18)); dash-dotted lines: log-wake-modified eddy viscosity (Equation (19)); solid line: second eddy viscosity (Equation (22)) with B f = 4 .
Figure 4. Vertical distribution of eddy viscosity; green dashed line: parabolic eddy viscosity (Equation (18)); dash-dotted lines: log-wake-modified eddy viscosity (Equation (19)); solid line: second eddy viscosity (Equation (22)) with B f = 4 .
Hydrology 08 00126 g004
Figure 5. Eddy viscosity profiles; thick red solid lines: first eddy viscosity (Equation (20)); thin magenta solid lines: first eddy viscosity with free surface damping function (Equation (21)); blue dash-dotted lines: log-wake-modified eddy viscosity (Equation (19)) with Π = 0.2 ; green dashed lines: parabolic eddy viscosity (Equation (18)).
Figure 5. Eddy viscosity profiles; thick red solid lines: first eddy viscosity (Equation (20)); thin magenta solid lines: first eddy viscosity with free surface damping function (Equation (21)); blue dash-dotted lines: log-wake-modified eddy viscosity (Equation (19)) with Π = 0.2 ; green dashed lines: parabolic eddy viscosity (Equation (18)).
Hydrology 08 00126 g005
Figure 6. Asymptotic behavior of the two coefficients of Equations (23) and (24): C α = α 1 κ = 0.477   and C 1 = 2.17 .
Figure 6. Asymptotic behavior of the two coefficients of Equations (23) and (24): C α = α 1 κ = 0.477   and C 1 = 2.17 .
Hydrology 08 00126 g006
Figure 7. Eddy viscosity profiles; red solid line: first eddy viscosity (Equation (24)); blue dashed lines: second eddy viscosity (Equation (22)) both with free surface damping function; red dash-dotted line: log-wake-modified (Equation (2b)) with Π = 0.2 ; blue dash-dotted line: Equation (2) with Π = 0.15 ; green dash-dotted line: parabolic eddy viscosity (Equation (18)).
Figure 7. Eddy viscosity profiles; red solid line: first eddy viscosity (Equation (24)); blue dashed lines: second eddy viscosity (Equation (22)) both with free surface damping function; red dash-dotted line: log-wake-modified (Equation (2b)) with Π = 0.2 ; blue dash-dotted line: Equation (2) with Π = 0.15 ; green dash-dotted line: parabolic eddy viscosity (Equation (18)).
Hydrology 08 00126 g007
Figure 8. Validation by experimental data; eddy viscosity profiles; red solid line: first eddy viscosity (Equation (24)); blue dashed lines: second eddy viscosity (Equation (22)) both with free surface damping function; symbols: experimental data of Nezu and Rodi (1986) and Ueda et al. (1976).
Figure 8. Validation by experimental data; eddy viscosity profiles; red solid line: first eddy viscosity (Equation (24)); blue dashed lines: second eddy viscosity (Equation (22)) both with free surface damping function; symbols: experimental data of Nezu and Rodi (1986) and Ueda et al. (1976).
Hydrology 08 00126 g008
Table 1. Flow conditions [18].
Table 1. Flow conditions [18].
CaseDepth of Flow, h (cm)Width to Depth Ratio Reynolds   Number   1 ,   R e = 4 R U m ν Froude   Number ,   F r = U m g h Friction   Reynolds   Number ,   R e * = h u * ν
P210.35.95.5 × 1040.189923
P310.06.014.3 × 1040.4882156
P410.06.021.0 × 1040.7043001
P510.55.744.0 × 1041.1706139
1  U m = mean bulk velocity, R = hydraulic radius.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Absi, R. Reinvestigating the Parabolic-Shaped Eddy Viscosity Profile for Free Surface Flows. Hydrology 2021, 8, 126. https://doi.org/10.3390/hydrology8030126

AMA Style

Absi R. Reinvestigating the Parabolic-Shaped Eddy Viscosity Profile for Free Surface Flows. Hydrology. 2021; 8(3):126. https://doi.org/10.3390/hydrology8030126

Chicago/Turabian Style

Absi, Rafik. 2021. "Reinvestigating the Parabolic-Shaped Eddy Viscosity Profile for Free Surface Flows" Hydrology 8, no. 3: 126. https://doi.org/10.3390/hydrology8030126

APA Style

Absi, R. (2021). Reinvestigating the Parabolic-Shaped Eddy Viscosity Profile for Free Surface Flows. Hydrology, 8(3), 126. https://doi.org/10.3390/hydrology8030126

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