Next Article in Journal
Synoptic Time Scale Variability in Precipitation and Streamflows for River Basins over Northern South America
Next Article in Special Issue
Supercritical Flow over a Submerged Vertical Negative Step
Previous Article in Journal
Flood Early Warning and Risk Modelling
Previous Article in Special Issue
Impacts of Climate Change and Variability on Precipitation and Maximum Flows in Devil’s Creek, Tacna, Peru
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determination of Skin Friction Factor in Gravel Bed Rivers: Considering the Effect of Large-Scale Topographic Forms in Non-Uniform Flows

by
Masoud Kazem
1,
Hossein Afzalimehr
1,
Mohammad Nazari-Sharabian
2 and
Moses Karakouzian
2,*
1
Faculty of Civil Engineering, Iran University of Science and Technology, Tehran 13114-16846, Iran
2
Department of Civil and Environmental Engineering and Construction, University of Nevada Las Vegas, Las Vegas, NV 89154, USA
*
Author to whom correspondence should be addressed.
Hydrology 2022, 9(4), 58; https://doi.org/10.3390/hydrology9040058
Submission received: 5 March 2022 / Revised: 28 March 2022 / Accepted: 30 March 2022 / Published: 31 March 2022

Abstract

:
Determination of skin friction factor has been a controversial topic, particularly in gravel-bed rivers where total flow resistance is influenced by the existence of small-scale skin roughness and large-scale topographic forms. The accuracy of existing models predicting skin friction factors in conditions where small-scale skin roughness and large-scale topographic forms exist is very low. The objective of this study is to develop a modified model that improves the accuracy of the determination of skin friction factors in gravel-bed rivers. To this end, 100 velocity profile data obtained from eight gravel-bed rivers were utilized to develop an analytical method that considers the momentum thickness of the boundary layer and its deviation in large-scale topographic bedforms in a 1D force-balance model. The results show that the accuracy of the skin friction factors is enhanced when (1) the model is in the form of an exponential function of energy slope, and (2) the deviation of momentum thickness is considered in the model. The proposed model results in high accuracy of the predicted skin friction factors for energy slopes between 0.001 and 0.1, which exist in most gravel-bed rivers with different morphologies. Additionally, this study model was used to modify the classic Einstein–Strickler equation. The modified equation resulted in improved accuracy of the predicted skin friction factors in non-uniform flow conditions even when velocity profiles and energy slope were not available.

1. Introduction

Friction factor and its sources have always been a major research topic [1,2,3,4]. During the past decades, the decomposition of friction factor has been a prevalent approach in sedimentology. Usually, the total friction factor is assumed to be the sum of the skin friction factor and bedform friction factor [5,6]. Briefly, large-scale impediments, such as bedforms or submerged vegetation patches, are responsible for the formation of large-scale flow structures. These large-scale flow structures transfer large amounts of momentum in the upper and the middle region of flow [7,8]. On the other hand, skin friction is the key factor in the momentum transfer procedures in the near-bed region of rivers. To be brief, not only does the skin friction factor play a key role in the flow resistance, but also it is the main index of shear stress which governs sedimentation and erosion processes. While the determination of skin friction factor has been focused on in the field and experimental studies, the results can also be used in numerical modeling of sedimentation. Such models can be used to predict the changes in bed geometry under different flow scenarios [9].
In the main, velocity profile has been widely used to predict bed friction factor in natural waterways. However, in the presence of a highly rough bed, using of velocity profile may lead to some unexpected results. This condition arises from the multi-layer structure of the boundary layer above the coarse sediments. Indeed, the lower layer of flow (or roughness layer), which is in contact with the coarse sediments, may have different characteristics in closely similar bed materials of river beds with different geometries of bedforms [10]. For small values of h/D, where h is the water depth, and D is a geometry parameter representative of the bed roughness (generally a grain diameter), the streamwise velocity profile is not self-similar inside the roughness layer [11,12]. This characteristic can only be detected by highly sensitive and expensive velocimeter tools, such as acoustic Doppler velocity profiler (ADVP). On the other hand, in a wide range of engineering works, the accuracy and sensitivity of velocimetry are limited by common tools, such as rotary current meters with poor spatial resolution, particularly in the near-bed zone. Consequently, there is a fundamental need for the development of new approaches that be applicable in the absence of sensitive and expensive velocimeter tools.

2. Materials and Methods

Amongst the well-known friction and flow resistance indices, the Darcy–Weisbach friction factor is one of the most applied ones in hydraulic calculations. The total friction factor (f) can be written as the following [13,14]:
f = 8 u u m 2
where u m is the average streamwise velocity in a particular section of flow, and u is the shear velocity which can be calculated via different methods. Application of the boundary-layer parameters based on the ASCE Task Force recommendation has led to a series of new approaches, such as the boundary-layer characteristics method (BLCM) [15]. Considering the BLCM, the shear velocity can be calculated as:
u = δ θ u m a x c δ  
where C is a constant with a value of 4.4 [16], δ is the boundary layer displacement thickness, θ is the momentum thickness of the boundary layer, and, umax is the maximum velocity of a particular velocity profile.   δ and θ can be calculated as [16]:
δ = 0 h 1 u u m a x d y  
θ = 0 h u u m a x 1 u u m a x d y  
Gravel bed rivers have recently been investigated in a distinct topic in which the friction factor can be divided into two parts [13,17,18,19,20]:
f = f + f  
where f is the skin friction factor (emerged from grains) and f is the form friction factor (emerged from large-scale topographic forms of the river bed). For plane beds, bed shear stress is the representative stress of the skin friction factor. Figure 1 shows a typical velocity distribution in which bed shear stress can be defined via momentum thickness of the boundary layer:
τ b e d = 1 b d D d x = d d x ρ U 2 θ  
where D is the drag force emerged by the bed surface and b is the width of the channel. On the other hand, for a wide channel, the shear stress can be calculated based on the Darcy–Weisbach equation:
τ b e d   = γ h s f   = γ h 1 8 f   U ¯ 2 h g = γ f   U ¯ 2 8 g  
where γ is the specific gravity of water, h is the flow depth, g is the gravity acceleration, and u ¯ is the weighted average velocity. However, in the above-mentioned equations, the roughness of the bed surface is not considered directly. For a relatively plane bed and uniform flow, Einstein suggested a modified form of Stickler’s Equation to relate the bed roughness (nb) to the median grain size (d50) [21].
n b = d 50 1 / 6 24  
By considering Einstein’s Equation, the Darcy–Weisbach friction factor can be calculated as:
f = g 72 d 50 d 1 / 3  
where d is the flow depth in a wide channel. In the same approach to relative roughness, Keulegan (1938) [22] developed a relatively simple but practical equation:
f = 2.03 l o g 12.2 h k s 2  
In this equation, h is the depth of flow and ks is Nikuradse equivalent roughness size. Determination of ks is a controversial issue and there is a wide range of estimations from 1.23d35 to 3d90 or 6.6d50 [23,24,25]. In another approach, the Shields parameter can be applied in the Keulegan Equation via a series of calculations. Critical shields parameter (τ*cr) can be defined as:
τ c r = τ γ S G 1 ) d 50 = h S f S G 1 d 50  
where ϒ is the specific gravity of water, and SG is the ratio of sediments’ specific gravity to water specific gravity. This equation is valid for both uniform and non-uniform flow if S f presents the true gradient of the energy loss. On the other hand, τ*cr can be written in association with S f as the following [14,26]:
τ c r = 0.15 S f 0.25  
Assuming ks = d50 in Equation (5), and by comparing Equations (6) and (7),   f can be calculated via logarithm rules. The final equation can be written as [14]:
f = 0.9742 1.5225 l o g S f 2  
While Einstein’s Equation (Equations (12) and (13)) is based on field and empirical studies on uniform flow and relatively plane bed, derivation of f in Equation (13) is not restricted to those conditions. This fundamental difference will lead to considerably different results when the two approaches are applied to a unique case.

3. Motivation and Objective

Generally, skin friction factors calculated via classic approaches (e.g., [21]) are developed based on uniform flow conditions where plane bed morphology is assumed and where no emergent or submerged vegetation and obstacles exist. Accordingly, existing models for predicting the skin friction factors for riverbeds where large-scale topographic forms and vegetation exist yield inaccurate results. Determination of friction factor in non-uniform flow is a relatively complicated task that is very far from the classic approaches and needs advanced algorithms and methodologies [4,15,19].
The objective of this study is to develop a modified model that improves the accuracy of the determination of skin friction factors in gravel-bed rivers. Accordingly, prevalent equations have been modified to increase their accuracy in the presence of large bedforms and non-uniform flow.

4. Methodology and Technical Approach

Since many classic equations of roughness and friction factor were developed based on the uniform flow condition, considering the fundamental assumptions in the application of the uniform flow is essential. According to (1/7)th power velocity profile law of Prandtl, the velocity distribution of a uniform flow can be described as [27,28]:
u U = y δ 1 / 7  
where δ is the boundary layer which is defined as the region adjacent to a surface over which the velocity changes from zero to the free-stream velocity (0.99 U) [29]. The momentum thickness is also calculable as:
θ = 7 72 δ = 7 72 y U u 7  
By considering the definition of δ , the momentum thickness can be defined as a linear function of flow depth:
θ = 7 72 y U 0.99 U 7 0.1043 y  
where y is the vertical elevation in which u = 0.99 U and is typically equal to the total depth of flow. Inserting Equation (16) into Equation (6) results in a differential form of bed shear stress for uniform flow:
τ U n i f o r m = d d x 0.1043 ρ U 2 y  
Assuming a condition in which skin friction factor is calculated via a classic approach (such as Einstein–Strickler, which relies on uniform flow) and also via a boundary layer characteristic approach such as in Equation (13), the ratio of skin friction factors can be written as:
f b e d f U n i f o r m = τ b e d τ U n i f o r m  
where f b e d is the calculated value of skin friction factor via momentum thickness and f U n i f o r m is the skin friction factor calculated via Einstein’s Equation. By inserting Equations (6) and (17) into Equation (18), the ratio of skin friction factors can be obtained as:
f b e d f U n i f o r m = d d x ρ U 2 θ d d x 0.1043 ρ U 2 y 9.587 y 0.99 θ  
Subsequently, the skin friction in a non-uniform flow can be estimated as:
f = 9.587 y θ f U n i f o r m  
Nonetheless, there is a fundamental assumption in this equation: all the momentum thickness of the boundary layer is produced by skin friction. However, despite uniform and even quasi-uniform flow, it is a controversial assumption for non-uniform flow. Figure 2 shows a conceptual scheme of momentum thickness above a large-scale topographic form. Owing to Equation (5), the momentum transition is affected by skin friction and form friction in a complex process. Consequently, a portion of the momentum thickness, which is not restricted to one as the upper limit, must be considered in Equation (20). This phenomenon can be included in Equation (20) in the form of:
f = 9.587 y ( φ θ ) f U n i f o r m  
and
φ θ d = 0.1043   f f U n i f o r m    
where φ is an indicator representing the portion of momentum thickness which is shaped by the skin friction and d is the flow depth in relatively shallow rivers. Equation (22) and the value of φ θ d are investigated via data gathered in 100 velocimetry stations of eight gravel-bed rivers, including four reaches in Iran and four in Italy (Table 1). Total and skin friction factors were respectively calculated via BLCM (Equations (1)–(4)) and the Keulegan approach (With a wide range of Ks values). The equivalent skin friction of uniform flow was also calculated via the Einstein–Strickler Equation (Equations (8) and (9)) by applying the value of d50.
Consequently, considering φ θ d as a nonlinear multivariable function of a series of non-dimensional parameters, this parameter can be defined as:
φ θ d = m F r c 1 d 50 d c 2 S f c 3 R e g r a i n c 4  
Subsequently, the modified momentum thickness ( θ s k i n ) which is formed by the skin friction would be:
θ s k i n = φ θ = m F r c 1 d 50 d c 2 S f c 3 R e g r a i n c 4 d  
Using the conjugate gradient approach to minimize the root-mean-square deviation (RMSD) for the 100 velocity profile, “m”, “c1”, “c2”, “c3”, and “c4” were respectively calculated for different Ks values. The results are presented in Table 2. RMSD should always be a non-negative value including 0, which represents the perfect fit, and generally, the lower value of RMSD is considered the better result.
While none of the classic variants of the Keulegan Equation contain flow parameters, the Lamb–Shields-based form of the Keulegan Equation (Equation (13)) was considered as the basis for comparison (Figure 3). Normalized root-mean-square deviation (NRMSD) was calculated for this comparison as the following:
N R M S D = 1 θ s k i n ¯   K e u l e g a n S h i e l d s 1 N θ s k i n   K e u l e g a n S h i e l d s θ s k i n 2 N  
In order to describe the application of the developed equation in a practical situation, a field case study was established to evaluate the skin friction factor where bed topography was captured with high resolution.

5. Field Study and Data Collection

The Marbor River is located in the central Zagros Mountains in the Dena region of the Isfahan province in Iran. A relatively straight reach of a local branch was selected for this study. The bed materials consist of gravel size grains with sparse round cobbles and negligible sand-size grains in the banks. The reach width varies from 3 to 4.5 m and the length of the study area was 15.2 m. Grain size distribution and the flow velocity were captured in five sections through 13 stations. The grain size distribution was calculated by using the Wolman method in all stations [30]. Figure 4 shows the general grain size distribution of the reach. Water surface elevation and topographic map of the bed were produced via accurate surveying. The topography of the bed, river plan, and the location of the sections are shown in Figure 5. The topography and velocimetry results represent a completely non-uniform flow condition in the selected reach. Figure 6 shows the profiles of the riverbed and water surface and also velocity distributions in the centroid of the reach. Bed slope has been calculated via accurate surveying and the 2D gravity projections on global coordinates have been included in the water fluxes [31,32]. All in all, data from 13 stations were used for the evaluation of skin friction factor by different methods. The displacement and momentum thicknesses of the boundary layer were calculated for each station through the 2D velocity profile (Table 3).
Subsequently, the skin friction factor values were evaluated using the following four methods:
a: Applying the friction factor of the equivalent uniform flow and the measured momentum thickness (Equation (20)): In this approach, the momentum thickness must be calculated via velocity profile. The results are comparable to the exact values of the flow characteristics method with equivalent Lamb–Shields parameter (Figure 7a). While there is a general correlation between these values, the deviation is also clear because of the φ parameter.
b: Applying the friction factor of the equivalent uniform flow and the modified momentum thickness (Equation (21)): Using the φ parameter in this approach leads to an independent equation that only relies on the energy slope, and the results are very close to the exact values. This approach is independent of velocity profile which makes it a suitable method in engineering works (Figure 7b). The applied equation is shown in Table 2 ( f = 0.3 S f 0.33 ). The calculated values of φ (using Equation (24)) and the measured values of this parameter are also shown in Table 4. Investigation of the values of φ at the 37 stations of the two gravel-bed rivers (Marbor and Gamasyab) shows that for the majority of the velocimetry stations, φ < 1 . For the centerline of the explored reach of the Marbor River, higher values of φ ( φ 1 ) are observed in the local contractions (Sec. I and Sec III), while the bed profile seems to have a much lower impact on the values of φ (Figure 8). The higher values of φ in the contractions might be related to the morphological drag of walls. Essentially, wall drag is the main source of scattering of results, particularly in natural channels [33]. However, Lamb et al. (2008) [26] mentioned that the effect of wall drag increases in the high bed slops (S > 0.02). In this study, this phenomenon was observed for two sections, although more data sets are required to evaluate the value of φ for various ranges of contractions. Figure 9 shows the surface waves near the wall of a contraction zone, which clearly shows the development of a horizontal boundary layer.
c: Applying modified friction factor of the equivalent uniform flow: Modified Einstein—Strickler Equation was used (by considering Ks = D50. See Table 2). This method is also independent of the velocity profile. The results contain a higher level of errors in comparison with the two previous methods, although it is less scattered in comparison with the classic Einstein–Strickler Equation (Figure 7c).
d: Applying friction factor of the equivalent uniform flow: The classic Einstein–Strickler Equation according to Equation (9) (Figure 7d).

6. Discussion

The findings of this study indicate that the skin friction factor can be calculated via Equation (21) with high accuracy in comparison with other classic methods. Many studies have shown that the average velocity is not the only relevant velocity scale in determining bed friction, and the effect of flow characteristics needs to be considered [34,35,36,37,38]. Correspondingly, the presented approach does not rely on average velocity as a single variable. Actually, it contains the effects of momentum thickness and its deviations in topographic forms. This method calculates the skin friction factor as an exponential function of energy slope, which means the skin friction factor rises when the energy slope increases. Lamb et al. (2008) [26] explained this phenomenon as a result of backwater effects and an associated pressure differential, which increased the mobility of particles on steeper slopes. This method also eliminated the need for velocity profiles which play a significant role in many fluvial parameters. However, the concept of the φ parameter is not restricted to the measurements where energy slope is available. Applying this concept in the classic Einstein–Strickler method leads to a modified version of the method, evaluating skin friction factor without considering energy slope and velocity profiles.
In order to compare the accuracy of the methods which are developed based on the φ parameter, the relative error ((fBLCMf’)/fBLCM) was calculated for 100 velocimetry points (Figure 10). This was done for a wide range of energy slopes (0.001–0.1). The results showed that Equation (21) contains less than 5% of error in comparison to the boundary layer flow characteristics method for energy slopes between 0.0014 and 0.056, considering the normal range of friction factor in gravel-bed rivers (Figure 10a). For the high energy slope (Sf > 0.06), the error increases linearly. The results also showed that for the conditions where the energy slope is not available, the modified Einstein’s Equation (based on the φ parameter’s concept) can reduce the error up to 50% on average (Figure 10b). In addition to the Marbor river, for other rivers with different energy slopes, the accuracy of results obtained from Equation (21) can be compared with the results obtained from the Einstein equation (Figure 11). According to the Figure, the Einstein equation is not appropriately accurate in low and high energy slopes. One reason can be attributed to the assumptions taken into account by Einstein. The normal or semi-normal flow regime is not prevalent in rivers with very steep energy slopes because of the tendency of flow to accelerate in such conditions. On the other hand, very low energy slopes are prevalent in the presence of flow blockage or decelerating flow. Despite the Einstein equation, Equation (21) takes energy slope into account. Consequently, the accuracy of the estimated skin friction factor increased on extreme slopes. However, it must be considered that the findings of this research are limited by the number of rivers where flow data are available. As presented in Figure 11a, the accuracy of results decreased sparsely in higher energy slopes. However, the accuracy of the results is still acceptable in comparison with classic methods. Despite the energy slope, relative errors do not comply with a predictable pattern for average velocity and relative roughness (Figure 12). The relative error remains below 6% in most cases. However, large relative errors have been observed in the presence of higher relative roughness. Consequently, it is recommended that future research studies focus on the accuracy of the proposed approach in rivers with higher relative roughness.

7. Conclusions

The skin drag in open-channel flow is determinable through the modified momentum thickness of the boundary layer. In order to derive a general equation, a 1D force-balance model was applied based on the momentum thickness of the boundary layer. The values of the initial model were compared to standard uniform flow models of skin friction. The deviation of the predicted values (known as φ ) was formulated through dimensionless parameters and an optimization procedure that used data sets of eight gravel-bed rivers with 100 velocimetry stations. Applying the φ parameter in the Einstein–Strickler Equation, the skin friction factor reveals an exponential function of the energy slope. The comparison with previous studies (e.g., [26]) shows a similar relation for the total friction factor and bed slope. This method was also considered for a field case study and the results showed a high correlation between the new method and the Keulegan method. The high correlation rate is also expectable for an energy slope range between 0.001 and 0.1, which involves the majority of gravel-bed rivers. The application of φ parameters can lead to engineer-friendly equations which are applicable in engineering works where velocity profile is not available. The results showed that for the conditions where the energy slope is not available, the modified Einstein’s Equation (based on the φ parameter’s concept) can reduce the error up to 50% on average.

Author Contributions

Conceptualization, H.A.; methodology, H.A.; software, M.K. (Masoud Kazem); validation, M.K. (Masoud Kazem), H.A. and M.N.-S.; resources, H.A.; data curation, M.K. (Masoud Kazem); writing—original draft preparation, M.K. (Masoud Kazem) and M.N.-S.; writing—review and editing, M.N.-S.; visualization, M.K. (Masoud Kazem); supervision, H.A., M.N.-S. and M.K (Moses Karakouzian). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Colosimo, C.; Copertino, V.A.; Veltri, M. Friction factor evaluation in gravel-bed rivers. J. Hydraul. Eng. 1988, 114, 861–876. [Google Scholar] [CrossRef]
  2. Rennie, C.D.; Millar, R.G. Discussion of “Estimation of Gravel-Bed River Flow Resistance” by Colin, D. Rennie and Robert, G. Millar. J. Hydraul. Eng. 1999, 125, 1317–1319. [Google Scholar] [CrossRef]
  3. Julien, P.Y.; Klaassen, G.J.; Ten Brinke, W.B.M.; Wilbers, A.W.E. Case study: Bed resistance of Rhine River during 1998 flood. J. Hydraul. Eng. 2002, 128, 1042–1050. [Google Scholar] [CrossRef] [Green Version]
  4. Milukow, H.A.; Binns, A.D.; Adamowski, J.; Bonakdari, H.; Gharabaghi, B. Estimation of the Darcy–Weisbach friction factor for ungauged streams using Gene Expression Programming and Extreme Learning Machines. J. Hydrol. 2019, 568, 311–321. [Google Scholar] [CrossRef]
  5. Zhenlin Li, M. Direct skin friction measurements and stress partitioning over movable sand ripples. J. Geophys. Res. Ocean 1994, 99, 791–799. [Google Scholar] [CrossRef]
  6. Hey, R.D. Flow resistance in gravel-bed rivers. J. Hydraul. Div. 1979, 105, 365–379. [Google Scholar] [CrossRef]
  7. Kazem, M.; Afzalimehr, H.; Sui, J. Formation of Coherent Flow Structures beyond Vegetation Patches in Channel. Water 2021, 13, 2812. [Google Scholar] [CrossRef]
  8. Kazem, M.; Afzalimehr, H.; Sui, J. Characteristics of Turbulence in the Downstream Region of a Vegetation Patch. Water 2021, 13, 3468. [Google Scholar] [CrossRef]
  9. Juez, C.; Battisacco, E.; Schleiss, A.J.; Franca, M.J. Assessment of the performance of numerical modeling in reproducing a replenishment of sediments in a water-worked channel. Adv. Water Resour. 2016, 92, 10–22. [Google Scholar] [CrossRef]
  10. Franca, M.J.; Ferreira, R.M.; Lemmin, U. Parameterization of the logarithmic layer of double-averaged streamwise velocity profiles in gravel-bed river flows. Adv. Water Resour. 2008, 31, 915–925. [Google Scholar] [CrossRef] [Green Version]
  11. Nikora, V.I.; Smart, G.M. Turbulence characteristics of New Zealand gravel-bed rivers. J. Hydraul. Eng. 1997, 123, 764–773. [Google Scholar] [CrossRef]
  12. Nikora, V.; Goring, D.; McEwan, I.; Griffiths, G. Spatially averaged open-channel flow over rough bed. J. Hydraul. Eng. 2001, 127, 123–133. [Google Scholar] [CrossRef]
  13. Yalin, M.S. Mechanics of Sediment Transport; Pergamon Press: New York, NY, USA, 1972; Volume 5, pp. 99–100. [Google Scholar]
  14. Afzalimehr, H.; Singh, V.P.; Najafabadi, E.F. Determination of form friction factor. J. Hydrol. Eng. 2010, 15, 237–243. [Google Scholar] [CrossRef]
  15. Afzalimehr, H.; Anctil, F. Vitesse de frottement associée à un écoulement non uniforme et une rugosité relative intermédiaire. J. Hydraul. Res. 2001, 39, 181–186. [Google Scholar] [CrossRef]
  16. Schlichting, H.; Gersten, K. Fundamentals of boundary-layer theory. In Boundary-Layer Theory; Springer: Berlin/Heidelberg, Germany, 2000; pp. 29–49. [Google Scholar] [CrossRef]
  17. Yen, B.C. (Ed.) Channel Flow Resistance: Centennial of Manning’s Formula; Water Resources Publication, LLC.: Littleton, CO, USA, 1992. [Google Scholar]
  18. Bertin, S.; Groom, J.; Friedrich, H. Isolating roughness scales of gravel-bed patches. Water Resour. Res. 2017, 53, 6841–6856. [Google Scholar] [CrossRef]
  19. Mendicino, G.; Colosimo, F. Analysis of Flow Resistance Equations in Gravel-Bed Rivers with Intermittent Regimes: Calabrian fiumare Data Set. Water Resour. Res. 2019, 55, 7294–7319. [Google Scholar] [CrossRef]
  20. Luo, M.; Wang, X.; Yan, X.; Huang, E. Applying the mixing layer analogy for flow resistance evaluation in gravel-bed streams. J. Hydrol. 2020, 589, 125119. [Google Scholar] [CrossRef]
  21. Einstein, H.A. Flow on a movable bed. In Proceedings of the Second Hydraulics Conference, Iowa City, IA, USA, 1–4 June 1942; Howe, J.W., Rouse, H., Eds.; University of Iowa: Iowa City, IA, USA, 1943; pp. 332–341, Bulletin 27. [Google Scholar]
  22. Keulegan, G.H. Laws of turbulent flow in open channels. J. Res. Bur. Stand. 1938, 21, 707–741. [Google Scholar] [CrossRef]
  23. Ackers, P.; White, W.R. Sediment transport: New approach and analysis. J. Hydraul. Div. 1973, 99, 2041–2059. [Google Scholar] [CrossRef]
  24. van Rijn, L.C. Equivalent roughness of alluvial bed. J. Hydraul. Div. 1982, 108, 1215–1218. [Google Scholar] [CrossRef]
  25. Hammond, F.D.C.; Heathershaw, A.D.; Langhorne, D.N. A comparison between Shields’ threshold criterion and the movement of loosely packed gravel in a tidal channel. Sedimentology 1984, 31, 51–62. [Google Scholar] [CrossRef]
  26. Lamb, M.P.; Dietrich, W.E.; Venditti, J.G. Is the critical Shields stress for incipient sediment motion dependent on channel-bed slope? J. Geophys. Res. Earth Surf. 2008, 113, F02008. [Google Scholar] [CrossRef] [Green Version]
  27. Sarma, K.V.; Lakshminarayana, P.; Rao, N.L. Velocity distribution in smooth rectangular open channels. J. Hydraul. Eng. 1983, 109, 270–289. [Google Scholar] [CrossRef]
  28. Zagarola, M.V.; Perry, A.E.; Smits, A.J. Log laws or power laws: The scaling in the overlap region. Phys. Fluids 1997, 9, 2094–2100. [Google Scholar] [CrossRef] [Green Version]
  29. Gadbois, J.; Wilkerson, G. Uniform flow development length in a rough laboratory flume. In Proceedings of the World Environmental and Water Resources Congress 2014, Portland, OR, USA, 1–5 June 2014; pp. 1234–1242. [Google Scholar] [CrossRef]
  30. Wolman, M.G. A method of sampling coarse river-bed material. EOS Trans. Am. Geophys. Union 1954, 35, 951–956. [Google Scholar] [CrossRef]
  31. Juez, C.; Soares-Frazao, S.; Murillo, J.; García-Navarro, P. Experimental and numerical simulation of bed load transport over steep slopes. J. Hydraul. Res. 2017, 55, 455–469. [Google Scholar] [CrossRef]
  32. Francalanci, S.; Solari, L.; Toffolon, M. Local high-slope effects on sediment transport and fluvial bed form dynamics. Water Resour. Res. 2009, 45, W05426. [Google Scholar] [CrossRef] [Green Version]
  33. Parker, G.; Wilcock, P.R.; Paola, C.; Dietrich, W.E.; Pitlick, J. Physical basis for quasi-universal relations describing bankfull hydraulic geometry of single-thread gravel bed rivers. J. Geophys. Res. Earth Surf. 2007, 112. [Google Scholar] [CrossRef]
  34. Papanicolaou, A.N.; Diplas, P.; Evaggelopoulos, N.; Fotopoulos, S. Stochastic incipient motion criterion for spheres under various bed packing conditions. J. Hydraul. Eng. 2002, 128, 369–380. [Google Scholar] [CrossRef]
  35. Schmeeckle, M.W.; Nelson, J.M. Direct numerical simulation of bedload transport using a local, dynamic boundary condition. Sedimentology 2003, 50, 279–301. [Google Scholar] [CrossRef]
  36. Zanke, U.C.E. On the influence of turbulence on the initiation of sediment motion. Int. J. Sediment Res. 2003, 18, 17–31. [Google Scholar]
  37. Wu, F.C.; Yang, K.H. Entrainment probabilities of mixed-size sediment incorporating near-bed coherent flow structures. J. Hydraul. Eng. 2004, 130, 1187–1197. [Google Scholar] [CrossRef] [Green Version]
  38. Vollmer, S.; Kleinhans, M.G. Predicting incipient motion, including the effect of turbulent pressure fluctuations in the bed. Water Resour. Res. 2007, 43, W05410. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Definition of momentum thickness in a fully developed boundary layer.
Figure 1. Definition of momentum thickness in a fully developed boundary layer.
Hydrology 09 00058 g001
Figure 2. Schematic illustration of momentum thickness near a topographic bed form.
Figure 2. Schematic illustration of momentum thickness near a topographic bed form.
Hydrology 09 00058 g002
Figure 3. Comparison between values of modified momentum thickness of skin friction and their correlation with the Keulegan–Shields approach.
Figure 3. Comparison between values of modified momentum thickness of skin friction and their correlation with the Keulegan–Shields approach.
Hydrology 09 00058 g003
Figure 4. Grain size distribution in the selected reach (Wolman method).
Figure 4. Grain size distribution in the selected reach (Wolman method).
Hydrology 09 00058 g004
Figure 5. Topography of the bed, river plan, and the location of the sections in the selected reach (The lowest point has been set as the baseline of local elevation code)—Marbor River.
Figure 5. Topography of the bed, river plan, and the location of the sections in the selected reach (The lowest point has been set as the baseline of local elevation code)—Marbor River.
Hydrology 09 00058 g005
Figure 6. Water surface elevation (light blue dashed line), riverbed elevation (violet line), and velocity profiles (dark blue lines) along the center of mass’ line (centroid line)—selected reach of Marbor River.
Figure 6. Water surface elevation (light blue dashed line), riverbed elevation (violet line), and velocity profiles (dark blue lines) along the center of mass’ line (centroid line)—selected reach of Marbor River.
Hydrology 09 00058 g006
Figure 7. Skin friction factor calculated using Equation (20) (a), Equation (21) (b), modified Einstein–Strickler Equation (c), and classic Einstein–Strickler Equation (d), in the Marbor River. The vertical axis shows the skin friction factor calculated with the flow characteristic method. The ranges of deviations are shown with dashed blue lines.
Figure 7. Skin friction factor calculated using Equation (20) (a), Equation (21) (b), modified Einstein–Strickler Equation (c), and classic Einstein–Strickler Equation (d), in the Marbor River. The vertical axis shows the skin friction factor calculated with the flow characteristic method. The ranges of deviations are shown with dashed blue lines.
Hydrology 09 00058 g007
Figure 8. 3D illustration of the measured φ in the investigated reach.
Figure 8. 3D illustration of the measured φ in the investigated reach.
Hydrology 09 00058 g008
Figure 9. Surface waves and the contraction streamlines, which are visible near the wall of a contraction zone (Marbor River).
Figure 9. Surface waves and the contraction streamlines, which are visible near the wall of a contraction zone (Marbor River).
Hydrology 09 00058 g009
Figure 10. The relative errors of Equation (21) (a) and Einstein Equation (b) in comparison with BLCM for determining the skin friction factor. The horizontal axis represents the total energy slope.
Figure 10. The relative errors of Equation (21) (a) and Einstein Equation (b) in comparison with BLCM for determining the skin friction factor. The horizontal axis represents the total energy slope.
Hydrology 09 00058 g010
Figure 11. Comparison between the values of skin friction factor calculated via BLCM and Equation (21) (a) and Einstein equation (b).
Figure 11. Comparison between the values of skin friction factor calculated via BLCM and Equation (21) (a) and Einstein equation (b).
Hydrology 09 00058 g011
Figure 12. Relative error of skin friction factors calculated using Equation (24) and the BLCM method in comparison with average velocity (left) and relative roughness (right).
Figure 12. Relative error of skin friction factors calculated using Equation (24) and the BLCM method in comparison with average velocity (left) and relative roughness (right).
Hydrology 09 00058 g012
Table 1. The rivers and profiles used to derive Equation (24).
Table 1. The rivers and profiles used to derive Equation (24).
RiverLocationNo. of ProfilesAvg. Depth (cm)Avg. Umean (cm/s)Avg. d50 (mm)Avg. FrAvg. fAverage f’ (Keulegan)Avg. Sf
MelodariItaly1413.664.9200.5541.2040.1010.038
CerasiaItaly1519.556470.4160.780.0720.016
ValanidiItaly815.861.4350.5021.0890.0860.026
GallicoItaly1324.572.2520.4701.1820.0880.028
ZayanderudIran57377100.2910.0920.0440.001
KajIran82763.9100.3960.1530.0420.004
GamasyabIran243084.2190.2870.1070.0430.003
MarborIran132292.6170.4330.1160.0500.006
Table 2. Calculated parameters to determine modified momentum thickness emerged by skin friction for different equivalent roughness in the Keulegan Equation.
Table 2. Calculated parameters to determine modified momentum thickness emerged by skin friction for different equivalent roughness in the Keulegan Equation.
θ s k i n = φ θ = m * F r c 1 * d 50 d c 2 * S f c 3 * R e g r a i n c 4 * d
ksmc1c2c3c4Pearson Correlation CoefficientEquivalent Equation for Skin Friction Factor
D500.0130−0.8000.89 0.017 d 50 d 0.47
1.5 D500.0140−0.85000.93 0.018 d 50 d 0.52
2 D500.0150−0.9000.94 0.019 d 50 d 0.57
Lamb–Shields Method (Equation (13))0.2250−0.330.3300.80 0.3 S f 0.33
Table 3. Measured and calculated flow parameters in the Marbor River.
Table 3. Measured and calculated flow parameters in the Marbor River.
StationIAICIIAIIBIICIIIAIIIBIIICIVAIVBIVCVAVB
d50 (mm) 20222117231713162016191916
d (cm) 21232121172321172923113325
Umean (m/s) 1.261.220.941.110.581.161.080.940.940.760.460.880.71
u* (m/s)0.160.170.140.1620.0920.0980.0860.0740.0950.0940.0360.140.075
Fr0.880.810.650.770.450.770.750.730.560.510.440.490.45
δ (m) 0.0340.0430.0490.0480.0380.0240.0240.0140.0490.0330.0080.0820.043
θ (m) 0.0190.0230.0250.0270.0180.0160.0170.010.0280.0220.0060.0430.028
f (-) 0.1230.1530.1750.1680.1990.0560.0500.0490.0820.1230.0480.1950.090
f (-) 0.0660.0670.0610.0670.0500.0470.0450.0440.0400.0400.0340.0520.040
f (-) 0.0570.0860.1140.1010.1480.0090.0050.0050.0420.0830.0140.1430.049
f Uniform (-) 0.0620.0620.0630.0590.0700.0570.0540.0620.0560.0560.0760.0530.055
Table 4. Comparison between measured and calculated φ.
Table 4. Comparison between measured and calculated φ.
SectionIIIIIIIVV
StationLeftRightLeftCent.RightLeftCent.RightLeftCent.RightLeftCent.
Calculated   φ 1.211.140.850.950.721.231.071.300.850.900.810.820.68
Measured φ 1.181.110.830.930.711.221.071.300.780.770.870.810.69
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kazem, M.; Afzalimehr, H.; Nazari-Sharabian, M.; Karakouzian, M. Determination of Skin Friction Factor in Gravel Bed Rivers: Considering the Effect of Large-Scale Topographic Forms in Non-Uniform Flows. Hydrology 2022, 9, 58. https://doi.org/10.3390/hydrology9040058

AMA Style

Kazem M, Afzalimehr H, Nazari-Sharabian M, Karakouzian M. Determination of Skin Friction Factor in Gravel Bed Rivers: Considering the Effect of Large-Scale Topographic Forms in Non-Uniform Flows. Hydrology. 2022; 9(4):58. https://doi.org/10.3390/hydrology9040058

Chicago/Turabian Style

Kazem, Masoud, Hossein Afzalimehr, Mohammad Nazari-Sharabian, and Moses Karakouzian. 2022. "Determination of Skin Friction Factor in Gravel Bed Rivers: Considering the Effect of Large-Scale Topographic Forms in Non-Uniform Flows" Hydrology 9, no. 4: 58. https://doi.org/10.3390/hydrology9040058

APA Style

Kazem, M., Afzalimehr, H., Nazari-Sharabian, M., & Karakouzian, M. (2022). Determination of Skin Friction Factor in Gravel Bed Rivers: Considering the Effect of Large-Scale Topographic Forms in Non-Uniform Flows. Hydrology, 9(4), 58. https://doi.org/10.3390/hydrology9040058

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