Next Article in Journal
The River Influence Controls Water Quality and Spatio-Temporal Microalgal Distribution in Pacific Estuaries (Padre Ramos and Salinas Grandes) of Nicaragua
Next Article in Special Issue
Turbulent Flow through Random Vegetation on a Rough Bed
Previous Article in Journal
Current Progress on Marine Microplastics Pollution Research: A Review on Pollution Occurrence, Detection, and Environmental Effects
Previous Article in Special Issue
Hydrodynamic Characteristics of Flow in a Strongly Curved Channel with Gravel Beds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Study on the Friction Factor and Reynolds Number Relationship for Flow in Smooth and Rough Channels

Department of Civil and Environmental Engineering, Pusan National University, Busandaehak-ro, 63beon-gil, Geumjeong-gu, Busan 46241, Korea
*
Author to whom correspondence should be addressed.
Water 2021, 13(12), 1714; https://doi.org/10.3390/w13121714
Submission received: 1 June 2021 / Revised: 16 June 2021 / Accepted: 16 June 2021 / Published: 21 June 2021

Abstract

:
The shear velocity and friction coefficient for representing the resistance of flow are key factors to determine the flow characteristics of the open-channel flow. Various studies have been conducted in the open-channel flow, but many controversies remain over the form of equation and estimation methods. This is because the equations developed based on theory have not fully interpreted the friction characteristics in an open-channel flow. In this paper, a friction coefficient equation is proposed by using the entropy concept. The proposed equation is determined under the rectangular, the trapezoid, the parabolic round-bottomed triangle, and the parabolic-bottomed triangle open-channel flow conditions. To evaluate the proposed equation, the estimated results are compared with measured data in both the smooth and rough flow conditions. The evaluation results showed that R (correlation coefficient) is found to be above 0.96 in most cases, and the discrepancy ratio analysis results are very close to zero. The advantage of the developed equation is that the energy slope terms are not included, because the determination of the exact value is the most difficult in the open-channel flow. The developed equation uses only the mean velocity and entropy M to estimate the friction loss coefficient, which can be used for maximizing the design efficiency.

1. Introduction

Head loss, hf, is a very important physical parameter for both the experimental and the theoretical analyses of fluid phenomena. The mechanism of the head loss in the open-channel flow is very complex and is not clearly explained yet. Usually, friction losses can be calculated whenever the friction coefficient, f, is clearly defined by the Darcy-Weisbach equation [1].
The friction coefficient is a key factor to determine the flow velocity in channel flows, which is also important to ensure the optimum hydraulic design. Furthermore, most studies on the friction coefficient f are not performed in an open channel but in a circular pipe flow. However, theoretical research for the open-channel flow was performed in the case of a relatively slow flow with no secondary current and small distribution. The theoretical analysis of the pipe flow was relatively easy compared to that of the open-channel flow [2].
The experimental data from Bazin [3,4] in the late 1800s and Varwick [5] in the early 1900s showed relationships between the friction coefficient and Moody curves in a pipe flow. For the open-channel flow, similar results by Bazin and Varwick [5] were presented. These results showed that the flow characteristics in the laminar, the transition, and the turbulence flows were similar to those in the pipe flow.
The Bazin, Manning [6], and Ganguillet-Kutter [7] equations were developed by using experimental data performed by Bazin and Varwick [5]. Many researchers, including Chezy [8] and Manning [6], developed empirical equations by using observed data in the open channel. However, these equations do not adequately represent the flow characteristics in the open channel. In order to improve the accuracy for representing the flow characteristics in an open channel, a sufficient study on the adequate and reliable analyses is required.
The Manning’s equation [6] was less accurate, even in well-flowing controlled man-made waterways. Chow [9] also recommended the adjustment of Manning’s roughness coefficient as a function of depth. Studies on the calculation of the friction coefficient and the head losses in the past are still in use or are not actively underway because of difficulties in the physical solution. This is based on studies done about 100 years ago. Additionally, Choo [10] used the entropy concept to derive the mean velocity using Chiu’s velocity formula [11,12], which was also used in this paper to derive the friction coefficient.
For the safe design, operation, and management of an increasingly developed and complex water resource facility, research on the calculations of a more accurate friction coefficient should be based on this study.
Therefore, this paper proposes a new theoretical equation to reflect the probabilistic entropy concepts and hydraulic properties. This article proposes a theory using the two-dimensional velocity formula and the probabilistic entropy to get the equation of the friction coefficient calculations. Equations based on these theories can be expected to be much more reliable than the empirical equations. The relationship between smooth and rough waterways is shown by means of comparison of the measured and estimated values to verify the accuracy of the determined values.
The results from the developed equation are based on a theoretical background. The friction coefficient was calculated directly based on the guidance equation combined with the physical factors. We found that it can be used to calculate the friction coefficients with very high accuracy without using energy slopes.

2. Theoretical Background

2.1. Existing Friction Loss Coefficient

The Hagen-Poiseule equation is shown as Equation (1), which calculates the friction head loss and can be written as:
h L = 8 μ u ¯ ω 0 R H 4 l
where hL is the friction head loss, μ is the fluid viscosity coefficient, ω0 is the water unit weight, RH is the hydraulic radius of a pipe, and l is the pipe length.
In this case, the friction head loss coefficient takes form by rearranging Equation (1) with RH = d/4, generating Equation (2) (Darcy-Weisbach [1]):
h L = f l d u ¯ 2 2 g
The representative equation for the friction loss coefficient in uniform flow can be written as Equation (3):
f = 8 g R h S V 2
where f is the friction coefficient, V is the mean velocity, g is the acceleration of gravity, R h is the hydraulic radius, S = h f / L is an energy slope, h f is the friction head losses, and L is the length of a given section.
This equation can be applied in streams close to the uniform flow, because it is difficult to calculate the value accurately in a nonuniform or unsteady flow, because S is an energy slope.
In a smooth pipe flow, the relationship can be expressed by using the Blasius [13] equation, written as Equation (4):
f = 0.223 R e 0.25
where R e is the Reynold’s number.
The equation is limited to valid values of R e between 750 and 25,000. For higher values, von Karman [14] developed a general expression modified by Prandtl [15], which matches the data measured by Nikuradse more accurately [16]. The resulting Prandtl-von Karman formula [5,14,15] can be written as Equation (5):
1 f = 2 log ( R e f ) + 0.4
Therefore, it is possible to establish a relationship between f and R e by an experiment in an open-channel flow using the above equations. However, the relation factors to f are different between the open-channel flow and the pipe flow, because it is affected by multiple factors, such as free water surfaces in the waterways, hydraulic radius, and water surface slope. Figure 1 shows the relationship between smooth and rough channel flows by analyzing various overseas experimental data [17].
The shortcomings of the above studies are that it is difficult to calculate the energy slopes correctly in an open-channel flow. In addition, equations should be applied differently according to the scope of the R e . For the example, the Prandtl-von Karman’s equation has to consider an uncertainty when using Equation (3). This is because it is hard to obtain an accurate flow velocity at the bottom of the channel.

2.2. New Friction Coefficient Using Entropy

Shannon [18] first defined entropy by function H(x) and can be written as Equation (6):
H ( x ) = + p ( x ) ln p ( x ) d x
where p ( x ) is the probability density function, and ln p ( x ) is dimensionless but dx has dimension.
Equation (6) means maximizing the entropy that represents the uncertainty for x , given p ( x ) for the continuous state variable x . Applying this concept to the water velocity can be written as Equation (7):
H ( u ) = 0 u m a x p ( u ) ln p ( u ) d u
The following constraints are used, such as the average value and probability, which are available information about u , which is instantaneous (point) velocity and can be written as Equations (8) and (9):
0 u m a x p ( u ) d u = 1
0 u m a x u p ( u ) d u = u ¯
Arranging the independent constraint conditions can be given as Equation (10):
a b Φ i ( u , p ) d u    i = 1 ,   2
where a is the minimum value of u , b is the maximum value of u , i is the constraint number ( i = 1 is Equation (8) and i = 2 is Equation (9)).
Therefore, p ( u ) , which maximizes the entropy, can be obtained using the method of Lagrange as Equations (11)–(13):
I ( u , p ) p + i = 1 2 λ i ϕ i ( u , p ) p = 0
I ( u , p ) = p ( u ) l n p ( u )
where ϕ 1 ( u , p ) = p ( u ) ,   ϕ 2 ( u , p ) = u × p ( u ) .
ϕ 1 ( u , p ) p = 1 ,    ϕ 2 ( u , p ) p = u
where λ 1 and λ 2 are Lagrange multipliers.
Substituting Equations (12) and (13) into Equation (11) can be constructed as the following Equation (14):
1 l n p ( u ) + λ 1 + λ 2 u = 0
where λ 1 1 = a 1   and   λ 2 = a 2 are the Lagrange multipliers.
Differentiating Equation (14) with respect to p ( u ) results in the velocity as Equation (15):
p ( u ) = e a 1 + a 2 u   ( 0 u u m a x )
Equation (15) and M = a 2 u m a x (entropy coefficient) are substituted into Equation (8) to obtain Equation (16):
u m a x e a 1 = M ( e M 1 )
Then Equation (15) and M = a 2 u m a x are substituted into Equation (9) to obtain Equation (17) (This is the two-dimensional average velocity equation, which is Chiu’s velocity equation [11,12]):
u ¯ u m a x = ϕ ( M ) = ( e M ( e M 1 ) 1 M )
Equation (17) can be used to restructure Equation (16) to obtain Equation (18):
e a 1 = ϕ ( M ) M ( e M 1 ) u ¯
The shear stress is the product of the dynamic viscosity (kinematic viscosity is the dynamic viscosity divided by density) and the velocity gradient, which can be expressed as Equation (19):
τ = ρ ν ( d u h ¯ ξ d ξ )
where τ is the shear stress, ρ is the density of the fluid, ν is the kinematic viscosity of the fluid, h ¯ ξ is the mean value of h ξ , and h ξ is the scale factor, which has the length dimensions.
The shear stress at the channel boundary (bottom) is the shear stress when ξ is ξ 0 , as in Equation (20):
τ 0 = ρ ν ( d u h ¯ ξ d ξ ) ξ = ξ 0 = ρ g R h S f
where τ 0 is the waterway boundary shear stress, g is the gravitational acceleration, and S f is the energy gradient.
The velocity cumulative probability of u is suggested by Chiu [11,12] as Equation (21):
P ( u ) = 0 u p ( u ) d u = ξ ξ 0 ξ m a x ξ 0
where ξ is the spatial coordinates ( 0 ξ 1 ), u is the velocity at ξ , ξ 0 is the minimum value of ξ (occurring at the channel boundary where u = 0 ), and ξ m a x is the maximum value of ξ (where u is at its maximum (i.e., u m a x )) (see Figure 2).
The ξ - η coordinates are the isovel system, which was first developed by Chiu [11,12] to explain two-dimensional velocity distribution in the cross-section of an open channel.
In Equation (21), if ξ is ξ 0 at the bottom of the channel, u is 0 and ξ m a x ξ 0 is 1, and by differentiating the velocity gradient in the channel bed, where p ( u ) = e a 1 , Equation (22) can be obtained:
( d u h ¯ ξ d ξ ) ξ = ξ 0 = 1 h ¯ ξ e a 1
For a wide channel, ξ / ξ max can be simplified to y / D and, hence, h ξ = D . For the hydraulic radius, R h D . Therefore, substituting Equation (22) into Equation (20), Equation (23) is obtained:
e a 1 = ν g R h 2 S f
Equating Equations (18) and Equation (22) expresses the average water velocity in the open-channel flow as Equation (24):
u ¯ = g R h 2 S f ν F ( M )
where F ( M ) = ( e M 1 ) / ( ϕ ( M ) M ) .
For the friction velocity ( u * ), the relationship between the average water velocity and the friction velocity is shown as Equation (25):
u * = τ 0 ρ
To calculate the friction term in Equation (25), Choo’s mean velocity distribution [10] is used for Equation (26):
u = u ¯ K ( M ) ln [ 1 + ( e M 1 ) ( ξ ξ 0 ξ m a x ξ 0 ) ]
where K ( M ) = ϕ ( M ) M .
Choo’s mean velocity was used earlier for calculating the discharge. However, in this paper, it will be used for converting friction velocity, since it has already been modified for the average water velocity in the open-channel flow [10].
The water velocity slope d u / h ¯ ξ d ξ is differentiated from Equation (26), and ξ 0 0 and ξ m a x ξ 0 = 1 are applied as Equation (27):
( d u h ¯ ξ d ξ ) = u ¯ ( e M 1 ) R h K ( M ) ( 1 + ξ ( e M 1 ) )
where, because the bottom boundary layer ξ 0 = 0 , Equation (27) is equal to Equation (28):
( d u h ¯ ξ d ξ ) ξ = ξ 0 = u ¯ ( e M 1 ) R h K ( M )
Equation (28) is inserted into Equation (20), which is the shear stress at the channel boundary, to obtain Equation (29):
τ 0 = ρ ν ( d u h ¯ ξ d ξ ) ξ = ξ 0 = ρ ν u ¯ ( e M 1 ) R h K ( M )
The relationship between the average water velocity and the friction velocity of the friction loss coefficient of the pipe flow is shown in Equation (30):
u ¯ u * = 8 f
Equations (25) and (30) are substituted into Equation (29) to obtain Equation (31):
τ 0 = 8 ρ ν 2 f ( ( e M 1 ) R h K ( M ) ) 2
Therefore, if Equation (31) is substituted for Equations (24) ( F ( M ) = ( e M 1 ) / ( ϕ ( M ) M ) ) and (26) ( K ( M ) = ϕ ( M ) M ) and use R e = ( u ¯ d ) / υ , Equation (32) can be obtained:
f = 8 d F ( M ) R e R h
Equation (32) can be used to estimate the frictional loss coefficient (f) of the open-channel flow, which reflects its entropy. Equation (32) does not require the hydraulic factors used in the existing equations, such as shear velocity ( u * ) or energy gradient ( S f ). In addition, the friction loss coefficient ( f ) can be expressed with only the average water velocity and the entropy M , which are easy to obtain. In addition, there is also an advantage in that the energy gradient ( S f ) can be estimated by using Equation (32).
Therefore, in this study, we proposed the friction loss coefficient of Equation (32) in the open-channel flow by using the concept of entropy, which has been used in many fields recently. The data used to demonstrate the utility of the equation were obtained by Yuen [19] and Babaeyan-Koopaei [20] for each stream of water. It is shown in the Figure 1 that the estimated friction loss coefficient was compared with the measured friction loss coefficient.

3. Experimental Data

To evaluate the accuracy of the proposed equation, we calculated the friction coefficient based on the data measured at the rectangular channel. The estimated results were compared with the measured data, as shown in Figure 1. First, the data measured by Yuen [19] at the trapezoidal section were used. Then, the data were measured by Babaeyan-Koopaei [20] at the trapezoid, the parabolic round-bottomed triangle, and the parabolic-bottomed triangle trapezoidal channel.
Yuen obtained data in a fully developed turbulent flow of the smooth trapezoidal open-channel flow. The ranges of the data were 0.5 < F r < 3.5 , 1.9 · 10 4 < R e < 6.2 · 10 5 and 0.3 < 2 b / H < 15 (where 2 b / H was the aspect ratio). The subcritical flow was also studied for the compound trapezoidal channel, which ranged in depths of 0.05 < D r < 0.5 . Here, Dr is the relative depth ratio ( ( H h ) H , where H is the flow depth of the channel, and h is the depth of the lower main channel). Several series of experiments were undertaken by using the Preston tube technique. These experiments were performed in a 21.26-m-long tilting channel with a working cross-section of 0.615 m wide × 0.365 m deep. A total of three sets were measured under the equivalent conditions, varying the bed slope at 0.001, 0.004, 0.009, 0.015, and 0.023. In addition, the point velocities were measured across the whole cross-section for the selected flow depths. Particular attention was focused on understanding the Reynolds and Froude number effects on these distributions.
Babaeyan-Koopaei [20] measured the data in the trapezoidal, parabolic, round-bottomed triangle, and parabolic-bottomed triangle channel. For each section, the measured data were used with the changes in the flow velocity and water levels under three flow conditions: 1   m 3 / s , 10   m 3 / s , and 100   m 3 / s (see Babaeyan-Koopaei for more information).
The values of the measured Re data are shown in Table 1. It can be seen that the measured data in the rectangular section reflected the transition zone and the turbulence zone. In the trapezoidal, parabolic, round-bottomed triangle, and parabolic-bottomed triangle sections, the measured data reflected the full turbulence zone.

4. Estimation of the Entropy Parameter, M

An estimate of the entropy parameter M is needed to use Equation (29). For the estimation of entropy parameter M, most researches used Equation (15) to calculate the entropy parameters in which the equation required, essentially, the maximum flow velocity in an open channel.
However, the maximum velocity occurred at the center of the pipe flow, but the location of the maximum velocity was unclear at the open-channel flow. Additionally, a lot of manpower, time, and effort were required to measure the maximum velocity in the open-channel flow. Moramarco [21] calculated the M values by using Equation (15). For that, he used data obtained from the average and maximum velocities at the upper river basin. Moramarco [22] proposed an equation for calculating ϕ ( M ) by substituting Chiu’s theory and the Manning and Prandtl-von Karman equations. However, the disadvantage of these equations were that it was difficult to clearly identify the point where the maximum velocity occurred, y m a x , and imaginary distance, y 0 , where the velocity was zero in the riverbed.
This study determined the entropy parameter M by using the expression developed by Choo [23]. The advantage of this method was that the entropy parameters in the stream could be obtained at any time without using the uncertain maximum velocity.
The entropy parameters M and the R e were calculated by using the same characteristics as those shown in Figure 3 and Figure 4. As the entropy parameter, M was increased, and the R e was also increased. On the other hand, as the friction coefficient f increased, R e decreased. Based on the value of R e , the two flows were identified as turbulent flows.

5. Results Analysis

The entropy parameter M, defined in Section 4, was used in Equation (32) to calculate the coefficient of friction f in an open-channel flow. The relationship between f and R e is shown in Figure 5 and Figure 6.
The coefficient of friction f in Figure 5 and Figure 6 shows the same trend as in Figure 1, where f tends to decrease as Re increases. In addition, in Figure 7 and Figure 8, the friction coefficient, f , shows a tendency to decrease as the discharge increases. The discrepancy ratio is the error ratio between the measurement and calculated values, separated by a range. The proportions on the y-axis are the ratio of the total comparison quantity to the range of the x-axis. Figure 9, Figure 10, Figure 11 and Figure 12 show that the discrepancy ratio results of the proposed equation were all distributed near 0.
The above results are summarized as follows: The entropy parameter M is a linear function of L o g ( R e ) , which is an increasing function of R e . The friction coefficient f is a linear decreasing function of L o g ( R e ) , which is a decreasing function of R e . Figure 13 shows the above results calculated using the proposed Equation (32), along with the relationship between the friction coefficients f and R e . Comparing the scale with a previous empirical study resulted in rough flows; the coefficient of determination was observed to be 0.8 within the range of the flow of the rectangular channel and 0.75 within the range of the flow of the compound trapezoidal channel.
Figure 10 shows the results determined with Yuen’s data measured in the rectangular and trapezoidal channels. Particularly, the calculated f value was expressed as 1.7 × 10 5 < R e < 6.18 × 10 6 , which exceeded the existing R e value. This range means that the proposed equation can represent the actual phenomenon in the natural stream.
The picture on the right shows the results determined with Babaeyan-Koopaei’s data measured in the trapezoidal, parabolic, round-bottomed triangle, and parabolic-bottomed triangle channels. In this case, the calculated f value was expressed as 13.5 × 10 4 < R e < 46.84 × 10 5 , which was close to the maximum range of the previous R e value. The important thing is that even the compound trapezoidal sections that were similar to the natural section, which were not available in the existing graph, as shown in Figure 1, can easily be calculated for f and R e . In Table 2 and Table 3, the regression equations from the Yuen [19] and Babaeyan-Koopaei [20] data show the relationship between M , f , R e , and Q .
Comparing the scale with a previous empirical study results in rough flows; the coefficient of determination was observed to be 0.8 within the range of the flow of the rectangular channel and 0.75 within the range of the flow of the compound trapezoidal channel.
As can be seen in Figure 13, the values of the coefficient of friction determined from past experiences are properly correlated with the measured data. There is no bouncing value on the graph. For the rectangular channel, this seems expressed fairly well as Bazin no. 17. Other types of channels matched well with the extended lines of the Prandtl-von Karman equation. In other words, it is meaningful that the values of the friction coefficient determined from proposed equation can be accurately estimate based on a theoretical formula rather than on an empirical parameter under a number of conditions.

6. Conclusions

The results of a study conducted approximately 100 years ago are still used to estimate the friction coefficient in an open-channel flow. However, as with the pipe flow (perfusion), the friction coefficient must be correctly determined in order to interpret the correct flow. This paper proposes a new form of friction coefficient calculation by using the two-dimensional velocity formula of Chiu [11,12] and probabilistic entropy.
The advantage of this equation is that it eliminates the terms of energy slopes, which are difficult to measure or calculate in an open-channel flow, making their application simple and very accurate on a theoretical basis.
In uniform flow conditions, a channel bed gradient may be the same or almost the same as an energy slope or water surface gradient. The normal depth is maintained as long as the slope, cross-section, and the surface roughness of the channel remains unchanged; thus, the average flow velocity remains constant. However, in natural flow and human-made open channels, such as irrigation systems and sewer lines, there are mostly nonuniform or unsteady flows. Unlike uniform flow conditions, these varied flows do not share the same energy slope, bed gradient, and water surface gradient.
Based on the data measured in the rectangular section, the proposed equation was used to determine the entropy parameter M and the friction coefficient f . The induced entropy parameters were shown to be a linear function of L o g ( R e ) , and the friction coefficient was the decreasing function of L o g ( R e ) .
If this study is to be carried out continuously by hydraulic data measured in various channel shapes, laboratory channels, and natural streams, the friction coefficient value estimated from the proposed equation will be actively used in the flow analysis and the design of hydraulic structures.

Author Contributions

Y.-M.C. and S.-H.P. carried out the survey of the previous studies. Y.-M.C. wrote the manuscript. S.-H.P. conducted all the simulations. Y.-M.C., J.-G.K., and S.-H.P. conceived the original idea of the proposed method. All authors have read and agreed to the published manuscript.

Funding

This research was funded by the Institute of Industrial Technology, Pusan National University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

This research was supported by the Institute of Industrial Technology, Pusan National University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Darcy, H.; Bazin, H. Recherches Hydrauliques; Enterprises par M.H. D’Arcy; Imprimerie Nationale: Paris, France, 1865. [Google Scholar]
  2. Yu, K.K. Particle Tracking of Suspended-Sediment Velocities in Open-Channel Flow. Ph.D. Thesis, Civil and Environmental Engineering, University of Iowa, Iowa City, IA, USA, 2004. [Google Scholar]
  3. Bazin, H.E. Recherches Experimentales sur Lecoulement de Leau dans les Canaux Decouverts; Memoire Presentes par Divers Savants al Academie des Sciences: Paris, France, 1865. [Google Scholar]
  4. Bazin, H.E. Etude d’une nouvelle formule pour calculer le debit des canaux decouverts. In Annales des Ponts et Chaussées; P. Vicq-Dunod: Paris, France, 1897; Volume 14, pp. 20–70. [Google Scholar]
  5. Varwick, F. Zur Fließformel für offene Künstliche Gerinne; Dresden University: Dresden, Germany, 1945. (In German) [Google Scholar]
  6. Manning, R.; Griffith, J.P.; Pigot, T.F.; Vernon-Harcourt, L.F. On the Flow of Water in Open Channels and Pipes; Transaction of the Institution of Civil Engineers of Ireland: Dublin, Ireland, 1890; Volume 20. [Google Scholar]
  7. Ganguillet, E.; Kutter, W.R. An investigation to establish a new general formula for uniform flow of water in canals and rivers. Z. Oesterreichischen Ing. Archit. Ver. 1869, 21, 6–25. [Google Scholar]
  8. Chezy, A. Thesis on the Velocity of the Flow in a Given Ditch; des Ponts et Chaussees, Library in France: Paris, France, 1775; Volume 847. [Google Scholar]
  9. Chow, V.T. Open-Channel Hydraulics; Mcgrow-Hill Civil Engineering Series; McGraw-Hill: Tokyo, Japan, 1959. [Google Scholar]
  10. Choo, T.H.; Yoon, H.C.; Lee, S.J. An estimation of discharge using mean velocity derived through Chiu’s velocity equation. Environ. Earth Sci. 2013, 69, 247–256. [Google Scholar] [CrossRef]
  11. Chiu, C.L. Velocity distribution in open channel flow. J. Hydraul. Eng. 1988, 115, 576–594. [Google Scholar] [CrossRef]
  12. Chiu, C.L. Application of entropy concept in open-channel flow study. J. Hydraul. Eng. 1991, 117, 615–628. [Google Scholar] [CrossRef]
  13. Blasius, H. Das Ähnlichkeitsgesetz bei Reibungsvorgängen in Flüssigkeiten; Forschungsheft des Vereins Deutscher Ingenieure: Berlin, Germany, 1913; Volume 131. [Google Scholar]
  14. Kảrmản, T. Mechanische Ähnlichkeit und turbulenz. In Proceedings of the 3rd International Congress for Applied Mechanics, Stockholm, Sweden, 24–29 August 1930; Volume 1, pp. 85–93. [Google Scholar]
  15. Prandtl, L. >The Mechanics of Viscous Fluids; Durand, B., Ed.; Aerodynamic Theory Springer: Berlin, Germany, 1935; Volume III, division G. 142. [Google Scholar]
  16. Nikuradse, J. Gesetzmassigkeiten der turbulenten Stromung in glatten Rohren; Ver Deutsch. Ing. Forschungsheft: Berlin, Germany, 1932; Volume 356. [Google Scholar]
  17. Ryu, H.J.; Kim, J.H.; Lee, B.H.; Lee, W.H.; Jang, I.S. Choisin Hydrography; Donghwa, Korea, 2007. [Google Scholar]
  18. Shannon, C.E. A mathematical theory of communication. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef] [Green Version]
  19. Yuen, W.H. A Study of Boundary Shear Stress, Flow Resistance and Momentum Transfer in Open Channels with Simple and Compound Trapezoidal cross Section. Ph.D. Thesis, Department of Civil Engineering, University of Birmingham, Birmingham, UK, 1989. [Google Scholar]
  20. Babaeyan-Koopaei, K.; Valentine, E.M.; Swailes, D.C. Optimal design of parabolic-bottomed triangle canals. J. Irrig. Drain. Eng. 2000, 126, 408–411. [Google Scholar] [CrossRef]
  21. Moramarco, T.; Saltalippi, C.; Singh, V.P. Estimation of mean velocity in natural channels based on Chiu’s velocity distribution equation. J. Hydrol. Eng. 2004, 9, 442–450. [Google Scholar] [CrossRef]
  22. Moramarco, T.; Singh, V.P. Formulation of the entropy parameter based on hydraulic and geometric characteristics of river cross sections. J. Hydraul. Eng. 2010, 15, 852–858. [Google Scholar] [CrossRef]
  23. Choo, Y.M.; Yun, G.S.; Choo, T.H.; Kwon, Y.B.; Sim, S.Y. Study of shear stress in laminar pipe flow using entropy concept. Environ. Earth Sci. 2017, 76, 1–27. [Google Scholar] [CrossRef]
Figure 1. Relationship between f and Re for smooth channels (Flow conditions: Laminar, Transitional, and Turbulent [14]).
Figure 1. Relationship between f and Re for smooth channels (Flow conditions: Laminar, Transitional, and Turbulent [14]).
Water 13 01714 g001
Figure 2. ξ - η of an open-channel flow (Chiu [11,12]).
Figure 2. ξ - η of an open-channel flow (Chiu [11,12]).
Water 13 01714 g002aWater 13 01714 g002b
Figure 3. Relationship between M and R e calculated using Yuen’s data.
Figure 3. Relationship between M and R e calculated using Yuen’s data.
Water 13 01714 g003
Figure 4. Relationship between M and R e calculated using Babaeyan-Koopaei’s data.
Figure 4. Relationship between M and R e calculated using Babaeyan-Koopaei’s data.
Water 13 01714 g004
Figure 5. The relationship between f and R e calculated using Yuen’s data.
Figure 5. The relationship between f and R e calculated using Yuen’s data.
Water 13 01714 g005
Figure 6. The relationship of f and R e calculated using Babaeyan-Koopaei’s data.
Figure 6. The relationship of f and R e calculated using Babaeyan-Koopaei’s data.
Water 13 01714 g006
Figure 7. The relationship of f and Q calculated using Yuen’s data.
Figure 7. The relationship of f and Q calculated using Yuen’s data.
Water 13 01714 g007
Figure 8. The relationship of f and Q calculated using Babaeyan-Koopaei’s data.
Figure 8. The relationship of f and Q calculated using Babaeyan-Koopaei’s data.
Water 13 01714 g008
Figure 9. The discrepancy ratio for f and R e calculated using Yuen’s data.
Figure 9. The discrepancy ratio for f and R e calculated using Yuen’s data.
Water 13 01714 g009
Figure 10. The discrepancy ratio for f and R e calculated using Babaeyan-Koopaei’s data.
Figure 10. The discrepancy ratio for f and R e calculated using Babaeyan-Koopaei’s data.
Water 13 01714 g010
Figure 11. The discrepancy ratio for f and Q calculated using Yuen’s data.
Figure 11. The discrepancy ratio for f and Q calculated using Yuen’s data.
Water 13 01714 g011
Figure 12. The discrepancy ratio for f and Q calculated using Babaeyan-Koopaei’s data.
Figure 12. The discrepancy ratio for f and Q calculated using Babaeyan-Koopaei’s data.
Water 13 01714 g012
Figure 13. The comparison between the empirical and calculated friction coefficients using Yuen and Babaeyan-Koopaei (rough channel).
Figure 13. The comparison between the empirical and calculated friction coefficients using Yuen and Babaeyan-Koopaei (rough channel).
Water 13 01714 g013
Table 1. The range of Reynolds numbers with the cross-section shape and the channel slope.
Table 1. The range of Reynolds numbers with the cross-section shape and the channel slope.
DataCross-Section ShapeChannel SlopeReynolds Number Range
Yuen [19]Rectangular0.00116,920~156,400
0.00445,770~160,900
0.00971,450~358,000
0.015108,600~335,000
0.023124,400~618,300
Babaeyan-Koopaei [20]Trapezoidal0.001167,000~4,474,000
Parabolic0.001135,000~4,630,000
Round-bottomed triangle0.001167,000~4,684,000
Parabolic-bottomed triangle0.001167,000~4,630,000
Table 2. The results for the relationship between M , f , R e , and Q using Yuen’s data.
Table 2. The results for the relationship between M , f , R e , and Q using Yuen’s data.
DivisionRegression EquationR2
M and ReYuen [19] y = 0.9561 ln ( x ) 4.925 R 2 = 0.8644
f and ReYuen [19] y = 0.3798 x 0.262 R 2 = 0.8456
f and QYuen [19] y = 0.007 x 0.213 R 2 = 0.7536
Table 3. The results for the relationship between M , f , R e , and Q using Babaeyan-Koopaei’s data.
Table 3. The results for the relationship between M , f , R e , and Q using Babaeyan-Koopaei’s data.
DivisionRegression EquationR2
M and ReBabaeyan-Koopaei [20] y = 0.8846 ln ( x ) 0.379 R 2 = 0.9994
f and ReBabaeyan-Koopaei [20] y = 1.5096 x 0.194 R 2 = 0.9895
f and QBabaeyan-Koopaei [20] y = 0.012 ln ( x ) 0.1236 R 2 = 0.9533
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Choo, Y.-M.; Kim, J.-G.; Park, S.-H. A Study on the Friction Factor and Reynolds Number Relationship for Flow in Smooth and Rough Channels. Water 2021, 13, 1714. https://doi.org/10.3390/w13121714

AMA Style

Choo Y-M, Kim J-G, Park S-H. A Study on the Friction Factor and Reynolds Number Relationship for Flow in Smooth and Rough Channels. Water. 2021; 13(12):1714. https://doi.org/10.3390/w13121714

Chicago/Turabian Style

Choo, Yeon-Moon, Jong-Gu Kim, and Sang-Ho Park. 2021. "A Study on the Friction Factor and Reynolds Number Relationship for Flow in Smooth and Rough Channels" Water 13, no. 12: 1714. https://doi.org/10.3390/w13121714

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