Next Article in Journal
A Review of the 21st Century Challenges in the Food-Energy-Water Security in the Middle East
Previous Article in Journal
Beta Diversity Partitioning and Drivers of Variations in Fish Assemblages in a Headwater Stream: Lijiang River, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fixed-Bed and Mobile-Bed Resistance of Channels with Steep Gradients in Mountainous Areas

1
Southwest Water Transport Engineering Research Institute, Chongqing Jiaotong University, Chongqing 400016, China
2
State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu 610065, China
*
Author to whom correspondence should be addressed.
Water 2019, 11(4), 681; https://doi.org/10.3390/w11040681
Submission received: 24 December 2018 / Revised: 27 March 2019 / Accepted: 28 March 2019 / Published: 2 April 2019
(This article belongs to the Section Hydraulics and Hydrodynamics)

Abstract

:
Flood discharge and sediment transport are closely linked to channel resistance in steep mountain streams. Previous research has mainly focused on the resistance of fixed-bed channels with steep gradients and mobile-bed channels in alluvial rivers. The present study performs an experiment and establishes a calculation method for the fixed-bed resistance of mountain channels. The basic expression of the mobile-bed resistance of steep mountain channels is derived by determining the controlling factors of the bed load movement on the riverbed resistance. The proposed formula can accurately predict the variation of the bed load resistance. The results of the present research improve the understanding of fluid dynamics and sediment transport in steep mountain channels.

1. Summary

The bed resistance of rivers is an important basis for flood control and engineering design. On the other hand, it is one of the basic challenges in hydraulics and river dynamics. Research on the resistance of mountainous rivers is even more challenging [1]. The most striking feature of rivers in the southwestern mountainous region of China is the narrow valleys with steep slopes. Bed material is mostly comprised of pebbles and gravel, which are coarse particles with wide gradation. Due to the unique topographical condition and bed surface composition, the flow characteristics are a low water depth, large flow velocity, large Froude number, and variable flow patterns. In order to understand the bed load transfer of mountainous rivers, explore the downward flow, and study the sediment transport law, it is necessary to understand the resistance of mountainous rivers and the important parameter of hydraulic factors.
In mountainous rivers, when the sediments are still or transported with low-intensity, the presence of sediment particles on the riverbed surface and their impacts on the water flow are regarded as the fixed bed resistance of the steep gradient channel. The resistance characteristics of the mountainous river, composed of coarse bed material, are quite different from those of alluvial rivers with fine bed particles. One of the important influencing factors is the relative roughness of the resistance of mountainous rivers [2,3,4,5]. Based on the relative roughness, researchers have divided the roughness degree into three levels. The division criteria of two studies are slightly different [5,6]. Other researchers have established resistance expressions [3,4,7,8,9]. In channels with a steep gradient in a mountainous area, the water flow is shallow, the bed material is coarse, and the intermediate scale roughness is more common, while large-scale roughness mostly occurs in a ladder-deep pool or in a pebble shoal. The flow resistance of a steep riprap, located in constructed channels, was investigated by previous researchers [10,11]. Flume experiments or field data show that the resistance of the steep gradient channels and alluvial rivers is related to the size, shape, and gradation of sediment particles and the geometric shape of rivers. The characteristics of the resistance of two conditions are different and the differences in the quantitative studies cannot be ignored. Most of the studies used the measured data of a certain river or river section to fit the expression [3,12,13,14,15]. Large error values existed in formulas due to the different choice of roughness size and the differences of the data. The present work continues to explore the fixed-bed resistance characteristics of steep gradient channels through a flume experiment. Moreover, it intends to analyze the fixed bed resistance of the non-uniform bed surface gradient, based on a large number of measured data. It is expected that the recommended formula can control the error within a good range.
In order to consider the sediment movement impact of the riverbed on the water flow structure, it is necessary to study the mobile bed resistance of the mountainous rivers. The mobile bed resistance for alluvial rivers is the resistance of sand waves and develops at various stages. However, the large amount of coarse grain mountainous pebbles in the riverbed makes it difficult to form sand waves. Therefore, the resistance of the mobile bed mainly comes from the resistance caused by the consumed flow energy to maintain the pebble bed movement. At present, researchers focus on the resistance of steep gradient channels with sediment-laden flow to derive semi-empirical and empirical formula, mostly based on flume experiments and field data. The formulas are mainly obtained through three methods.
First, a large amount of experimental and measured data is used to fit a logarithmic formula, similar to the one for the fixed bed resistance. For steep channels with high relative flow submersion, researchers [16,17] found that the flow resistance has a correlation with relative submersion and it was found that the correlation has a logarithmic form. The corresponding hydraulic roughness as a single characteristic roughness length scale is frequently determined based on the equivalent grain size of the riverbed. Researchers conducted a flume experiment study on the resistance law and derived the expression, as follows [15]:
8 f = 0.67 + 3.2 ln ( R D )
where f, R, and D denote the friction factor, hydraulic radius, and grain size, respectively.
The research [15] shows that when there is a bed load movement in the water flow, the resistance is 2.5 times higher than the fixed bed resistance. The river resistance is constant and the resistance value is between that for a stationary condition and initiation sediment motion. The resistance of sediment-laden flow increases as the sediment transport rate increases. There is a significant defect in such expressions for the resistance. The variables in the formula are the hydraulic radius and the sediment particle size, while the impact of the sediment transport rate is ignored.
The second method assumes that the resistance of the sediment flow is the superposition of the fixed bed resistance with the same roughness and the resistance, caused by the bed load movement. Under this circumstance, researchers mainly studied the law of resistance change of the bed load movement [18]. The law of the resistance change in the bed movement is
f = f c + f b
where f c is the fixed bed resistance. f b is the resistance, caused by bed load. There is a nonlinear correlation between the resistance expression, caused by the bed load and the sediment concentration:
f b = 0.00025 C v 1.23
where Cv is the dimensional mass concentration. It is the ratio of the volumetric flow of sediment ( Q s ) to the discharge flow which carries the sediment (Q):
C v = Q s / Q
Gao [18] believed that the transport concentration, relative roughness, and dimensionless particle size are the resistance affecting factors of the bed. They used the flume test and experiment data to establish an expression for the resistance of the bed load movement, based on dimensionless analysis [14]:
f b = 0.048 C v 0.25 d * 0.5 ( h D ) 0.75
where d * = ( r s r r g n 2 ) 1 / 3 D 50 . Moreover, ρ s , ρ , and ν are the sediment density, fluid density, and the coefficient of the fluid motion, respectively.
Based on the research results of Gao and Abrahams [18], a series of experiments are conducted to study the sediment resistance of the slope flow. The Froude number (Fr) of the water flow, as an effective factor in the flow resistance, is considered in accordance with the study of Gao and Abrahams [18] to establish an expression for calculating the resistance of the slope flow.
f b = 282.5 C v 0.579 d * 0.8 ( h D ) 0.25 F r 3.539 S 1.195
where S is the gradient of the water surface.
The third method is used to compare the resistance coefficient ratios of the mobile and fixed beds and fit the expression of the ratio and influencing factors. Researchers have conducted a mobile bed experiment [14]. They considered the dimensionless particle size and sediment transport concentration as the influence factors and derived the following expressions for the flow resistance under the condition of bed load movement:
f f c = ( 30.1 d * C v + 1 ) 0.92
h h c = ( f f c ) 1 / 3
where h and hc are the water depth during the bed load movement and the clear water depth under, respectively.
However, few researchers have focused on the flow resistance of mountainous channels. The general studies of moving bed resistance focus on the sand wave resistance of alluvial rivers at various stages. In the past, most of the experiments in the field of large-scale river channel mobile bed resistance were performed to study the resistance characteristics of water flow with sediment through the destruction of the bed surface. There are few studies focused on the influence of the sediment motion impact on river resistance of the stable bed surface. The present work carries out the flume experiment to analyze the variation of the resistance of the mobile bed in steep gradient channels.

2. Experiment Setup

The resistance of the bed with a uniform sediment experiment of steep gradient channels was carried out at the State Key Laboratory of Hydraulics and Mountain River Development and Protection, Sichuan University. The length, width, and height of the flume were 12.0 m, 0.4 m, and 0.4 m, respectively. Moreover, the slope was set to 10%, and was set in accordance with the mean gradients in mountain streams of southwest China. Figure 1 presents the configuration of the flume.
The sides of the flume were made up of 1 cm thick tempered glass, which were fixed to the aluminum alloy frame structure. The head of the flume had water control valves and static pools. Water flowed through the rectangular weir into the static pool. Then, the flow passed through the two-stage energy dissipation grille before it flowed into the flume. The flume end was provided with a flat gate and a sand grafting device to control the water level by adjusting the opening of the gate so that an approximately uniform flow was obtained. It should be indicated that the flow was considered to be uniform when the gradient of the water surface was close to the slope of the bed surface.
An automatic water level meter was mounted horizontally on the aluminum alloy rack at the top of the flume. In this experiment, the length of the test section was 3.75 m. Four automatic water level gauges were arranged on the test section, where each water level meter worked independently. The test data was recorded by connecting to the same computer. In this fixed bed with uniform sediment experiment, natural pebbles were laid in the base of whole flume, where the layer thickness was 3 cm. In order to ensure the reliability, representativeness, and comparability of the data, the bed material was made of two kinds of uniform sediments with particle sizes of 9.2 mm and 12.3 mm.
In the experiment of sediment laden resistance law in the steep gradient channel, the uniform sediment (non-uniform coefficient less than 1.3) with four particle sizes (3.8 mm, 6.7 mm, 9.2 mm, and 12.3 mm) was selected as the sediment entry channel. A funnel was arranged in the inlet section to control the sediment supplement rate through the opening of the funnel. Devices were installed at the end of the flume to measure the sediment transport rate. Four automatic water level meters simultaneously recorded the water depth variation of each section.
In this paper, the bed surface was assumed to be stable. Therefore, the pebble on the bed surface was fixed in the flume with a little cement slurry to simulate the stable bed surface. The relative discharge error was less than 3% when the flow passed through the rectangle thin-wall weir.
The first step of the sediment supplement was to put the prepared sediment into the flume through the sanding funnel. Then, different sediment transport rates were obtained through the bottom opening degree. The water level meter recorded the water level variations before and after the sanding funnel and recorded the bed surface elevation in advance. Based on the data from the spot measurement, a relatively accurate water depth variation could be obtained by averaging. The water level was controlled by the downstream baffle to make the water flow sufficiently stable with uniform turbulence. The added sediments neither accumulated nor exchanged with the sediment on the bed surface and were finally carried downstream by the flow.

3. Data Description and Method

3.1. Experimental Study on the Fixed Bed with Uniform Sediment Resistance of Steep Gradient Channnel

The coefficients of the river resistance are mainly expressed by the Chezy coefficient ( C ), the Darcy-Weisbach coefficient ( f ), and the Manning roughness coefficient (n). There is a mutual transformation correlation between these coefficients:
C g = R 1 / 6 g 1 n = 8 f = U U *
where U is the average flow velocity.
In order to adjust the test setup to large specific slope setting, the flow rate is set to higher than 1 m/s and the variation of the water depth is limited to 1.63 cm to 6.18 cm. Especially, the drag coefficient presents the effect of several phenomena, such as bed interaction and non-uniform velocity distribution [19]. Since the above expression of the drag coefficient is equal to the ratio of the average flow velocity to the friction velocity, the ratio can be obtained by integrating the flow velocity distribution. Therefore, the study of the resistance is actually a study of the flow velocity distribution.
According to the sidewall roughness and flow conditions, the fixed bed resistance can be divided into the resistance of smooth, transition, and rough regions [20]:
  • k s / δ > 10 ; Rough region
  • k s / δ < 0.25 ; Smooth region
  • 0.25 < k s / δ < 10 ; Transition region
Aguirre-Pe et al. [13] and Graf et al. [6] pointed out that the flow velocity above the bed surface in the vertical direction is basically in line with the logarithmic distribution, as follows [6,13]:
u u * = 1 κ ln y k s + B
where u * , κ , and k s are the friction velocity, Karman constant (value is 0.4), and bed roughness size, respectively. It should be indicated that the Karman constant is set to κ = 0.4. B is a constant. Montes [21] and Aguirre-Pe et al. [13] found that when the bed material composition in the steep gradient channel is uniform, the flow velocity on the surface of the river bed is constant within a certain range [11,21]. This range can be expressed as: z 0 = α k s , where α is a constant. The average flow velocity of the uniform sediment (see Figure 2) bed surface is obtained as follows:
U = u 0 z 0 h + z 0 h u d y
According to Equations (10) and (11):
8 f = 1 k ln h k s + B 1 k + z 0 k h + z 0 h ( u 0 u * 1 k ln z 0 k s B )
Since y = z 0 and u = u 0 , Equation (12) can be re-written as the following:
8 f = 1 κ ln h k s + B 1 κ + z 0 κ h
In the fixed bed resistance of the uniform sediment river bed, k s is the particle size of the sediment particles D . Moreover, according to the data of steep gradient channels, z 0 = 0.2 D [12,13,22,23]. Based on the experimental data in this paper, combined with the data of Song and Recking et al. (see Table 1) [15,24], the best value of parameter B is 6.5.
Graf et al. [6] deduced the following expression for the fixed bed resistance:
8 f = u u * = 1 κ ln R k s + B r
where B r is a constant as the following:
{ R / k S > 20 ,    B r = 6.5 4 < R / k S < 20 ,    3.25 < B r < 6.5 R / k S < 4 ,    B r = 3.25
Figure 3 compares the data of Song and Recking [15,24] and the experimental results of the present study, obtained from Equations (12) and (14). Table 2 shows the hydraulic conditions and results of each group.
The results show that all these points fall within the range of expressions proposed by Graf [6]. Equation (12) can better reflect the fixed bed resistance of the uniform sediment in the steep gradient channels.

3.2. Investigation of the Resistance Characteristics of the Non-Uniform Sediment River Bed

In natural mountainous rivers, the riverbed surface is mostly composed of uneven coarse and fine particles. It is of great importance to study the resistance of the non-uniform bed surface of the steep gradient rivers and the effects of concealment and exposure between the coarse and fine particles in the water flow [3,8,12,25].
For the fixed bed resistance of the non-uniform bed material in the cobble channel of the mountainous area, researchers have analyzed a huge amount of measured field data and flume experiments to obtain an applicable expression for the resistance in a certain river section or channel. For the fixed bed resistance in mountainous rivers, it is difficult to introduce a universally applicable empirical expression through the measured data or the flume experiment data. Therefore, the estimation error is 25–30%, which is acceptable for practical applications [3].
By collecting a large amount of field measured and experimental data, it is found that the error value of the expression proposed by Hey, Griffiths and Graf [8,12,25] for the large-scale pebble channel is greater than 50%.
The basic expressions of resistance of the mountainous rivers are basically consistent. However, the different data selected by the researchers in the process of derivation and verification caused minor differences in the k s value and coefficients. According to the analysis discussed above, the uniform expression of the fixed bed resistance can be written as
8 f = A log ( h k s ) + L
where A and L are constants. Different scholars have proposed different values for A. However, the proposed values mostly vary from 5.6 to 5.75. The average value (i.e., A = 5.70) is considered in this study. At present, most studies choose k s = α D 50 . Kikkawa et al [26] believed that the roughness size k s is independent of the particle size. However, the former value reflects the protruding height and geometry of the sand wave [26]. It is assumed that in mountainous rivers, the probability of sand waves occurrence is extremely small so the influence of sand waves is neglected in the present study. When the non-uniformity coefficient of the bed load is large and the gradation is wide, D84 is selected for the roughness size. For the case where the non-uniformity coefficient is small and the gradation is narrow, D50 is selected for the roughness size.
It is assumed that the fixed bed resistance is mainly affected by the bed material area on the riverbed plane, along the water depth. Therefore, the representative particle size D j is initially determined, and ks is calculated through k s = α D j . Based on the above assumptions, the proposed calculation method of the present study is
α = D 50 D j i n p i D i 4 D 50 4
where D i can be divided into the following cases (see Figure 4):
According to the choices of D i , Equation (17) can be re-written as
α = 0.4 D 20 4 + 0.2 D 50 4 + 0.1 D 65 4 + 0.3 D 84 4 D 50 4 D 50 D j
Based on Equation (18), several values of α can be obtained (see Table 3) [2,3,8,27,28].
It is found that L = 5.50 and Equations (16) and (17) can better calculate the fixed bed resistance than using the field measured data to verify the fitting curve. Figure 5 and Figure 6 show the comparison and error analysis of the corresponding values in Equation (16), respectively. Moreover, the results indicate that the selection of α is based on the gradation of the bed surface sediment so that there are different α values for different rivers. Equation (16) can basically reduce the error within a better range in calculating the fixed bed resistance of the river channel in the mountainous area.

3.3. Law of Sediment Laden Flow Resistance in the Steep Gradient Channels

Table 4 shows 58 groups of results with various test conditions.
Figure 7 compares the results from the fixed bed resistance equation, the ones from the large-scale flume experiments, and the results from the experiment carried out in the present work. It is observed that the sediment movement improves the resistance, so the resistance is larger than that for the fixed bed condition. There are several sets of tests in the experiments with different conditions of the sediment movement. Under this circumstance, the value of f / f c is less than or equal to 1.0, which indicates that the resistance of the mobile bed does not change in comparison with the fixed bed resistance. According to the resistance variation of the mobile bed (Equation (1)) [29,30] and Equations (3), (5), and (6), it is found that, when there is a mass movement in the water flow:
f b > 0
Then:
f / f c > 1
Flow resistance may decrease or remain unchanged, due to the movement of the bed load. Sediments need energy from water to move, when they enter the river, which increases the river resistance. However, for the river channel with a stable bed surface, the coarse particles in the bed load do not exchange with the moving sediment. When the finer sediment enters the channel, it moves in the form of rolling or jumping and fills the coarse grain skeleton of the bed load. With sufficient sediment supplements, the fine sediment covers the original bed surface to form a smoother bed surface. Therefore, although the bed load movement requires a part of the energy to be consumed from the water flow, the resistance decreases, when the fine-grained bed load enters the bed surface with a large bed-grain size. The movement of the fine particles makes the bed surface smoother, resulting in less energy loss of the water flow. If the energy loss is greater than the energy required for the movement of the bed load, then the river resistance may decrease.
However, the resistance may not decrease when small size sediments enter the channel. In the experiment carried out in the present work, when the D 50 sediment (i.e., particle parameter is 6.7 mm) enters the D 50 material (i.e., sediment diameter is 9.2 mm and 12.3.mm), the resistance still increases. Therefore, it is speculated that only when the particle size of the sediment carried by the water flow satisfies a certain range may the resistance decrease due to the movement of sediment in the water flow.
It is assumed that the flow velocity in the y-direction follows the expression below (see Figure 8):
u y = ( 1 c s y ) u c y + u s y c s y
where u c y and u s y are the flow velocity of the clean water and the flow velocity of the sediment in the y-direction, respectively. Moreover, c s y is the sediment concentration in the y-direction. Performing the integral on Equation (21) yields the average flow velocity as the form below:
U = 0 h ( ( 1 c s y ) u c y + u s y c s y ) d y h
Equation (22) is summarized as
U = ( 1 C s ) U c ¯ + U s C s
where U c and U s are the flow average velocity of the clean water and the flow average velocity of the sediment, respectively.
Bagnold [30] found that, when sediment particles move in the water stream, there is a certain difference between the velocity of the sediment particles ( u s ) and the water velocity ( u n ) at the position y n . This difference is
u r = u n u s
Moreover, the relative velocity between the sediment particles and the water flow movement is equivalent to the sedimentation speed of the sediment particles:
u r = w
where w = 1.72 ( γ s γ γ ) g D , R e > 10 3 .
According to Equation (25), it is assumed that there is a certain relationship between the average velocity of sediment particle motion and the sedimentation speed, as the following:
U s = U c ¯ η w
where α is a constant coefficient.
Combining Equations (26) and (23) yields the following:
U = U c ¯ η C s w
A series of experiments are carried out on the moving speed of the sediment particles (see Table 5) for the coefficient η. In the flume with a specific slope of 0.1, the velocity of the sediments with multiple sizes is measured at different water depths. The analysis of the experimental results shows that η is a function of the shield stress, in the following form:
η = m Θ
where m = 0.23 is a constant coefficient.
Θ = ρ u * 2 ( ρ s ρ ) g D
Figure 9 illustrates a comparison of the velocity of sediment particles, calculated by Equations (26) and (28) with the measured data. The results show that Equations (26) and (28) can accurately predict the velocity of sediment particles.
The resistance of the sediment laden flow and the clean water can be expressed as
8 f = U U *
8 f c = U c ¯ U * c
where U * c is the friction velocity in sediment-laden flow.
Combining Equations (30) and (31) yields the following:
f f c = ( U c ¯ U U * U * c ) 2
Considering Equation (27):
U U c ¯ = 1 η C s w U c ¯
U * U * c = g h J g h c J = h h c
Applying Equations (33) and (34), Equation (32) can be re-written as
f f c = ( 1 1 η C s w U c ¯ ) 2 h h c
Resistance and water depth has the following correlation:
h h c = ( f f c ) 1 / 3
Re-arranging Equations (35) and (36) results in the following equation:
f f c = ( 1 1 η C s w U c ¯ ) 3
The relationship between the sediment transport rate and the ultimate sediment transport rate is uncertain. Therefore, a constant coefficient β is added and Equation (37) takes the form below:
f f c = ( 1 1 η β C s w U c ¯ ) 3
Therefore, the resistance variation after water flow and sedimentation can be predicted by Equation (38). In order to verify the accuracy of the proposed equation, the results from Equation (38) are compared with data reported by Song et al. [14] and Gao and Abrahams [18] and the results from the experiment carried out in the present study (see Table 6). It is observed that the calculation error is still in the acceptable range and the error value is less than 10%. Figure 10 shows the comparison of the calculation result of Equation (38) and the measured value. It is found that, when f / f c is small, the calculation result is acceptable. On the other hand, when f / f c is too large, although the data points are slightly divergent, they are still in the acceptable range.

4. Conclusions

The law of water and sediment movement in the steep gradient channel is a vital basic theory in solving the problem of natural water and material destruction conditions of the mountainous rivers in Southwest China. Therefore, further in-depth exploration and research on this topic has significant theoretical and practical application value. This paper used the combination of theoretical analysis and experimental research to systematically explore and analyze the fixed and mobile bed resistance of mountainous rivers. The following conclusions have been obtained:
  • Based on the existing results, a calculation scheme is established to simulate the velocity distribution of the uniform sediment bed on the steep gradient channels. Based on the test data, the proposed expressions are verified to have a reasonable performance;
  • Through the existing research results, an expression is determined to calculate the resistance of the fixed bed with non-uniform sediment. Moreover, the selection of the roughness size is analyzed to put forward a new calculation method. The expression for the resistance of the non-uniform bed surface of the steep gradient channels is obtained through analyzing the measured field data;
  • When calculating the resistance of the mobile bed in mountainous rivers, the fixed bed resistance cannot be used. Experimental results show that, when the bed surface is stable and the bed load particle size is much smaller than that of the gravel, the bed surface is smooth and flat. At this point, there are several possibilities so the resistance may increase, may remain unchanged, or even may decrease;
  • By assuming the flow velocity distribution of clean water and sediment laden flows, the variation law for the resistance of sediment laden flow is derived. Moreover, the relationship between the velocity of the sediment movement and the average flow velocity is obtained. It is found that the proposed method can accurately predict the variation of the sediment carrying resistance.
Because of the various influencing factors and complicated water characteristics in the mountainous rivers, the research in this field is still not mature enough. There are still a lot of possibilities to be explored in the study of river resistance and sediment transport characteristics in steep gradient mountainous areas.

Author Contributions

Conceptualization, J.H.; Investigation, J.H. and Q.Z.; Methodology, C.Z. And J.H.; Validation, J.H and F.L.; Visualization, J.H.; Writing—original draft, J.H. and Z.Y.; Writing—review and editing, J.H. and D.W.

Funding

This project is supported by the Chongqing Research Program of Application Foundation and Advanced Technology [grant number cstc2016jcyjA1937/1935], science and technology project of Chongqing Education Committee [grant number KJ1600535/KJ1600514], and Open Fund Project of Key Laboratory for Inland Waterway Regulation Engineering of Ministry of communication in Chongqing Jiaotong University [grant number NHHD-201513/201505].

Acknowledgments

Many thanks to the Key Laboratory for Inland Waterway Regulation Engineering of Ministry of communication in Chongqing Jiaotong University.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Nitsche, M.; Rechenman, D.; Turoswski, J.M.; Badoux, A.; Kirchner, J.W. Evaluation of bedload transport predictions using flow resistance equations to account for macro-roughness instep-pool streams. Water Resour. Res. 2011, 47, W08513. [Google Scholar] [CrossRef]
  2. Bathurst, J.C. Flow resistance of large-scale roughness. J. Hydraul. Div. 1978, 104, 1587–1603. [Google Scholar]
  3. Bathurst, J.C. Flow resistance estimation in mountain rivers. J. Hydraul. Eng. 1985, 111, 625–643. [Google Scholar] [CrossRef]
  4. Maxwell, A.R.; Papanicolaou, A.N. Step-pool morphology in high-gradient streams. Int. J. Sediment Res. 2001, 16, 380–390. [Google Scholar]
  5. Yu, B.Y. Flow resistance of large-scale roughness in open channel. In Proceedings of the International Symposium on Hydraulic Research in Nature and Laboratory, Wuhan, China, 17–20 November 1992. [Google Scholar]
  6. Graf, W.H.; Altinakar, M. Hydraulique Fluviale; tome 1; Presses Polytechniques et Universitaires Romandes: Lausanne, Switzerland, 1993. [Google Scholar]
  7. Bathurst, J.C.; Thorne, C.R.; Hey, R.D. Secondary flow and shear stress at river bends. J. Hydraul. Div. 1979, 105, 1277–1295. [Google Scholar]
  8. Hey, R.D. Flow resistance in gravel-bed rivers. J. Hydraul. Div. 1979, 105, 365–379. [Google Scholar]
  9. Bathurst, J.C.; Li, R.M.; Simons, D.B. Resistance equation for large scale roughness. J. Hydraul. Div. 1981, 107, 1593–1613. [Google Scholar]
  10. Rice, C.E.; Kadavy, C.K.; Robinson, K.M. Roughness of loose rock riprap on steep slope. J. Hydraul. Eng. 1998, 124, 179–185. [Google Scholar] [CrossRef]
  11. Pagliara, S.; Chiavaccini, P. Flow resistance of rock chutes with protruding boulders. J. Hydraul. Eng. 2006, 132, 545–552. [Google Scholar] [CrossRef]
  12. Griffith, G.A. Flow resistance in coarse gravel bed rivers. J. Hydraul. Div. 1981, 107, 899–918. [Google Scholar]
  13. Aguirre-Pe, J.; Fuentes, R. Resistance to flow in steep rough streams. J. Hydraul. Eng. 1990, 116, 1374–1387. [Google Scholar] [CrossRef]
  14. Song, T.; Chiew, Y.M.; Chin, C.O. Effect of bed-load movement on flow friction factor. J. Hydraul. Eng. 1998, 124, 165–175. [Google Scholar] [CrossRef]
  15. Recking, A.; Frey, P.; Paquier, A.; Belleudy, P.; Champagne, J.Y. Feedback between bed load and flow resistance in gravel and cobble bed rivers. Water Resour. Res. 2008, 44, W05412. [Google Scholar]
  16. Ferguson, R. Flow resistance equations for gravel- and boulder-bed streams. Water Resour. Res. 2007, 43, W05427. [Google Scholar] [CrossRef]
  17. Aberle, J.; Smart, G.M. The influence of roughness structure on flow resistance on steep slopes. J. Hydraul. Res. 2003, 41, 259–269. [Google Scholar] [CrossRef]
  18. Gao, P.; Abrahams, A.D. Bedload transport resistance in rough open-channel flows. Earth Surf. Process. Landf. 2004, 29, 423–435. [Google Scholar] [CrossRef]
  19. Cassan, L.; Roux, H.; Garambois, P.A. A semi-analytical model for the hydraulic resistance due to macro-roughnesses of varying shapes and densities. Water 2017, 9, 637. [Google Scholar] [CrossRef]
  20. Chien, N.; Wan, C.H.; McNown, J.S. Mechanics of Sediment Transport; ASCE: Reston, VA, USA, 1999. [Google Scholar]
  21. Montes, S. A new theory about bidimensional flow over a rough boundary. In Proceedings of the 3rd Latin America Congress of the IAHR, La Plata, Agentia, 1968; pp. 100–118. (In Spanish). [Google Scholar]
  22. Bathurst, J.C.; Thorne, C.R.; Hey, R.D. Direct measurements of secondary currents in river bends. Nature 1977, 269, 504–506. [Google Scholar] [CrossRef]
  23. Thompson, S.M.; Campbell, P.L. Hydraulics of a large channel paved with boulders. J. Hydraul. Res. 1979, 17, 351–354. [Google Scholar] [CrossRef]
  24. Song, T.; Graf, W.H. Nonuniform open-channel flow over arough bed. J. Hydrosci. Hydr. Eng. 1994, 12, 1–25. [Google Scholar]
  25. Graf, W.H.; Cao, H.H.; Suszka, L. Hydraulics of steep, mobile bed channels. In Proceedings of the 20th Congress of the IAHR, Moscow, Russia, 5–9 September 1983; pp. 301–305. [Google Scholar]
  26. Kikkawa, H.; Hiroyoshi, S.; Fukuoka, S. On the effects of suspended sediments to the bed roughness. In Proceedings of the 12th Congress of IAHR, Fort Collins, CO, USA, 11–14 September 1967; pp. 49–56. [Google Scholar]
  27. Barnes, H.H., Jr. Roughness Characteristics of Natural Channels; U.S. Geological Survey Water Supply Paper 1849: Washington, DC, USA, 1967.
  28. Judd, H.E.; Peterson, D.F. Hydraulics of Large Bed Elements Channels; Report No. PRW17-6; Utah Water Research Laboratory, Utah State University: Logan, UT, USA, 1969. [Google Scholar]
  29. Abrahams, A.D.; Li, G. Effect of saltating sediment on flow resistance and bed roughness in overland flow. Earth Surf. Process. Landf. 1998, 23, 953–960. [Google Scholar] [CrossRef]
  30. Bagnold, R.A. The nature of saltation and of bed-load transport in water. Proc. R. Soc. A Math. Phys. Sci. 1973, 332, 473–504. [Google Scholar] [CrossRef]
Figure 1. Schematic of the flume layout.
Figure 1. Schematic of the flume layout.
Water 11 00681 g001
Figure 2. Schematic diagram of the flow velocity distribution for uniform sediment bed.
Figure 2. Schematic diagram of the flow velocity distribution for uniform sediment bed.
Water 11 00681 g002
Figure 3. Comparison of the test data with the calculated values from Equations (12) and (14).
Figure 3. Comparison of the test data with the calculated values from Equations (12) and (14).
Water 11 00681 g003
Figure 4. The distribution of D i .
Figure 4. The distribution of D i .
Water 11 00681 g004
Figure 5. The comparison and error analysis of Equation (16) with the measured data (a) D j = D 50 , α = 4.82 , (b) D j = D 50 , α = 4.51 , and (c) D j = D 50 , α = 5.64 .
Figure 5. The comparison and error analysis of Equation (16) with the measured data (a) D j = D 50 , α = 4.82 , (b) D j = D 50 , α = 4.51 , and (c) D j = D 50 , α = 5.64 .
Water 11 00681 g005
Figure 6. The comparison and error analysis of Equation (16) with the measured data (a) D j = D 84 , α = 2.25 , (b) D j = D 84 , α = 2.35 , and (c) D j = D 84 , α = 2.51
Figure 6. The comparison and error analysis of Equation (16) with the measured data (a) D j = D 84 , α = 2.25 , (b) D j = D 84 , α = 2.35 , and (c) D j = D 84 , α = 2.51
Water 11 00681 g006
Figure 7. Comparison of the results from the fixed bed resistance equation and those from other researchers.
Figure 7. Comparison of the results from the fixed bed resistance equation and those from other researchers.
Water 11 00681 g007
Figure 8. Comparison of the measured and calculated values of sediment particle motion velocity.
Figure 8. Comparison of the measured and calculated values of sediment particle motion velocity.
Water 11 00681 g008
Figure 9. Comparison of the measured and calculated values for the sediment particle motion velocity.
Figure 9. Comparison of the measured and calculated values for the sediment particle motion velocity.
Water 11 00681 g009
Figure 10. Comparison of the results from Equation (38) and those from other works.
Figure 10. Comparison of the results from Equation (38) and those from other works.
Water 11 00681 g010
Table 1. Data range of former researchers.
Table 1. Data range of former researchers.
TypeNumber of TestSD (cm)h (cm)
Song160.5–1.51.237.7–20.4
Recking621.0–9.00.23–0.91.07–8.33
Paper’s experiment1810.00.92, 1.231.63–6.18
Table 2. Uniform sediment bed resistance test data of the steep gradient channel.
Table 2. Uniform sediment bed resistance test data of the steep gradient channel.
No.Discharge Q (m3/s)Particle Size D (cm)Gradient RatioWater Depth h (cm)Flow Velocity U (m/s)R (cm)R/d(8/f)0.5
10.00510.920.11.630.7771.5071.6386.146
20.01000.920.12.281.0932.0472.2257.310
30.01240.920.12.581.2052.2852.4847.577
40.01350.920.12.681.2642.3632.5697.798
50.01630.920.12.981.3652.5942.8197.985
60.01790.920.13.181.4072.7442.9827.969
70.02440.920.13.881.5753.2503.5328.078
80.03020.920.14.281.7623.5263.8328.604
90.03120.920.14.381.7793.5933.9068.585
100.03730.920.14.881.9133.9234.2648.748
110.04280.920.15.332.0064.2084.5748.779
120.01001.230.12.680.9302.3631.9215.736
130.01321.230.13.081.0692.6692.1706.154
140.01961.230.13.831.2783.2142.6136.595
150.02681.230.14.531.4783.6933.0037.015
160.03581.230.15.281.6934.1773.3967.444
170.03841.230.15.531.7364.3323.5227.459
180.04791.230.16.181.9374.7213.8387.870
Table 3. Value of α in different equations.
Table 3. Value of α in different equations.
D j = D 84 α D j = D 50 α
Barnes (1967)2.35Barnes (1967)5.64
Judd & Peterson (1969)2.34Bray (1979)4.82
Bathurst (1978)2.52Griffiths (1978)4.51
Hey (1979)2.51
Bathurst (1985)2.25
Table 4. Test results for mobile bed resistance flume.
Table 4. Test results for mobile bed resistance flume.
No.Q (L/s)DB (mm)DT (mm)hc (cm)h (cm)SU (m/s)Cvffcf/fc
117.899.203.803.103.140.101.440.00910.11810.11671.0123
217.899.203.803.103.070.101.440.03020.11540.11670.9890
317.899.203.803.103.440.101.440.04590.12950.11671.0700
417.899.203.803.103.280.101.440.03370.12350.11671.0581
524.459.203.803.753.780.101.630.01670.11160.11071.0080
624.459.203.803.753.880.101.630.02760.11450.11071.0347
724.459.203.803.753.670.101.630.03430.10830.11070.9789
835.769.203.804.494.330.101.990.03090.08550.08880.9635
935.769.203.804.494.520.101.990.02160.08930.08881.0058
1042.229.203.805.105.200.102.070.00730.09520.09331.0196
1140.039.203.805.004.980.102.000.02670.09750.09790.9960
1241.679.203.805.085.030.102.050.02770.09380.09470.9902
1312.809.203.802.422.490.101.320.02580.11160.10851.0289
1412.809.203.802.422.600.101.320.03700.11660.10851.0744
1512.809.203.802.422.780.101.320.04420.12460.10851.0400
1617.899.206.703.103.860.101.440.03310.14530.11671.2452
1717.899.206.703.103.490.101.440.03070.13140.11671.1258
1817.899.206.703.103.310.101.440.00780.12460.11671.1123
1917.899.206.703.103.330.101.440.01660.12540.11671.0742
2024.459.206.703.754.230.101.630.02340.12480.11071.1280
2124.459.206.703.754.400.101.630.02800.12990.11071.1733
2224.459.206.703.754.170.101.630.02300.12310.11071.1120
2335.769.206.704.494.750.101.990.01790.09390.08881.0579
2435.769.206.704.494.850.101.990.01640.09590.08881.0802
2535.769.206.704.494.990.101.990.03100.09870.08881.1114
2641.679.206.705.085.440.102.050.01840.10140.09471.0987
2741.679.206.705.085.680.102.050.03730.10590.09471.2100
2840.039.206.705.005.450.102.000.02970.10670.09791.0900
2912.809.206.702.423.310.101.320.04930.14840.10851.3678
3012.809.206.702.422.740.101.320.01370.12280.10851.1322
3112.809.206.702.423.040.101.320.03980.13630.10851.2562
3212.809.206.702.422.750.101.320.01680.12330.10851.1100
3317.899.209.203.103.910.101.440.04170.14720.11671.3600
3417.899.209.203.103.680.101.440.01490.13850.11671.1600
3517.899.209.203.103.760.101.440.02420.14160.11671.1800
3624.459.209.203.754.320.101.630.02480.12750.11071.1520
3724.459.209.203.754.570.101.630.04250.13490.11071.3200
3835.769.209.204.495.100.101.990.02830.10090.08881.1359
3935.769.209.204.494.940.101.990.02020.09770.08881.1002
4035.769.209.204.494.950.101.990.01480.09790.08881.1024
4149.629.209.206.006.350.102.070.01470.11650.11001.0980
4249.049.209.205.906.500.102.080.02250.11800.10711.1017
4349.049.209.205.906.640.102.080.02550.12060.10711.1254
4412.809.209.202.423.530.101.320.03750.15830.10851.4587
4512.809.209.202.422.940.101.320.02200.13180.10851.2149
4612.809.209.202.422.840.101.320.01330.12730.10851.1736
4749.049.2012.305.906.290.102.080.00440.11420.10711.1600
4849.049.2012.305.906.210.102.080.00940.11280.10711.1900
4912.0712.303.803.003.080.101.010.02560.23870.23251.0267
5012.0712.303.803.003.140.101.010.04700.24330.23251.0467
5112.0712.306.703.003.470.101.010.04790.26890.23251.1567
5226.3112.306.704.304.520.101.530.02490.15150.14411.1100
5334.7212.306.705.005.160.101.740.02930.13420.13011.1230
5412.0712.309.203.003.750.101.010.04460.29060.23251.2500
5526.3112.309.204.304.970.101.530.03440.16650.14411.1558
5635.7612.309.205.105.640.101.750.02390.14390.13011.1059
5726.3112.3012.304.305.000.101.530.01140.16750.14411.1628
5833.1812.3012.304.905.410.101.690.01920.14800.13411.2000
Table 5. Velocity of the particle motion.
Table 5. Velocity of the particle motion.
No.h (m)D (m)Jw (u/s)Uc (u/s)Us (u/s) u *   ( u / s )
10.010.0050.10.4891.2601.0630.099
20.010.010.10.6921.2600.9410.099
30.010.020.10.9781.2600.7230.099
40.010.0250.11.0941.2600.6380.099
50.010.030.11.1981.2600.5810.099
60.010.040.11.3831.2600.5290.099
70.010.050.11.5471.2600.4390.099
80.020.010.10.6922.2232.1850.140
90.020.020.10.9782.2231.7690.140
100.020.0250.11.0942.2231.6250.140
110.020.030.11.1982.2231.5890.140
120.020.040.11.3832.2231.3810.140
130.020.050.11.5472.2231.2460.140
140.030.010.10.6922.7512.6440.171
150.030.020.10.9782.7512.4010.171
160.030.030.11.1982.7512.1960.171
170.030.040.11.3832.7512.0220.171
180.030.050.11.5472.7511.8520.171
190.040.010.10.6922.8142.7720.198
200.040.020.10.9782.8142.6240.198
210.040.030.11.1982.8142.5270.198
220.040.040.11.3832.8142.3810.198
230.040.050.11.5472.8142.1780.198
240.050.010.10.6922.8642.9120.221
250.050.020.10.9782.8642.7840.221
260.050.030.11.1982.8642.6220.221
270.050.040.11.3832.8642.5110.221
280.050.050.11.5472.8642.2690.221
Table 6. Comparison of the resistance for mobile beds of large gradient channels and those from the Song experiment.
Table 6. Comparison of the resistance for mobile beds of large gradient channels and those from the Song experiment.
No.h (cm)S (%)D (mm)Cvffcf/fc (Exp)f/fc (Auth)Difference (%)
120.100.5012.33.58 × 10−70.0660.0661.00001.00540.54
220.900.5012.31.15 × 10−60.0680.0651.04621.0381−0.77
321.300.5012.31.72 × 10−60.0660.0641.03131.04981.80
421.400.5012.31.76 × 10−60.0670.0641.04691.05060.36
521.800.5012.32.37 × 10−60.0650.0641.01561.05964.33
622.500.5012.32.91 × 10−60.0670.0631.06351.06660.30
722.700.5012.33.93 × 10−60.0630.0621.01611.07545.83
823.500.5012.35.74 × 10−60.0660.0611.08201.08760.52
924.000.5012.36.06 × 10−60.0650.0611.06561.09012.30
1024.500.5012.35.66 × 10−60.0650.0601.08331.08900.53
1125.000.5012.38.19 × 10−60.0640.0601.06671.10053.17
1225.300.5012.38.57 × 10−60.0630.0601.05001.10234.98
1316.500.7512.37.30 × 10−60.0740.0721.02781.06563.68
1417.100.7512.31.57 × 10−50.0750.0711.05631.08863.06
1517.700.7512.31.56 × 10−50.0760.0701.08571.09000.39
1618.300.7512.33.55 × 10−50.0770.0691.11591.1151−0.08
1719.100.7512.34.17 × 10−50.0800.0671.19401.1217−6.06
1819.400.7512.35.98 × 10−50.0770.0671.14931.1330−1.41
1919.900.7512.37.69 × 10−50.0770.0661.16671.1416−2.15
2020.300.7512.31.19 × 10−40.0760.0661.15151.15560.35
2121.100.7512.31.17 × 10−40.0740.0641.15631.15690.06
2213.000.9012.33.86 × 10−60.0890.0811.09881.0322−6.05
2313.600.9012.31.08 × 10−50.0880.0791.11391.0622−4.64
2414.100.9012.31.39 × 10−50.0870.0781.11541.0708−4.00
2514.600.9012.32.70 × 10−50.0860.0761.13161.0910−3.59
2615.300.9012.36.00 × 10−50.0890.0751.18671.1160−5.96
2715.700.9012.36.98 × 10−50.0860.0741.16221.1216−3.49
2816.000.9012.39.16 × 10−50.0820.0731.12331.13040.63
2916.500.9012.31.50 × 10−40.0820.0721.13891.14640.66
3016.900.9012.31.66 × 10−40.0800.0711.12681.15062.12
3117.300.9012.31.78 × 10−40.0780.0701.11431.15383.55
3217.800.9012.32.10 × 10−40.0790.0701.12861.16012.80
3318.700.9012.33.11 × 10−40.0840.0681.23531.1744−4.93
3412.201.0012.32.10 × 10−50.0940.0831.13251.0725−5.30
3512.801.0012.32.90 × 10−50.0950.0811.17281.0837−7.60
3613.401.0012.37.26 × 10−50.0950.0801.18751.1119−6.36
3713.901.0012.31.19 × 10−40.0940.0781.20511.1280−6.40
3814.301.0012.31.48 × 10−40.0900.0771.16881.1357−2.83
3914.901.0012.31.93 × 10−40.0920.0761.21051.1455−5.37
4015.401.0012.33.20 × 10−40.0910.0741.22971.1623−5.49
4116.201.0012.33.50 × 10−40.0870.0731.19181.1674−2.05
4210.301.2512.37.67 × 10−50.0980.0911.07691.09421.61
4311.101.2512.31.50 × 10−40.1040.0871.19541.1169−6.57
4411.701.2512.32.54 × 10−40.1050.0851.23531.1347−8.14
4512.201.2512.33.26 × 10−40.1030.0831.24101.1441−7.81
4612.601.2512.33.95 × 10−40.1000.0821.21951.1513−5.59
4713.801.2512.35.83 × 10−40.1030.0781.32051.2674−4.02
4814.601.2512.36.45 × 10−40.0980.0761.28951.1732−9.02
4915.701.2512.36.42 × 10−40.1000.0741.35141.2765−5.54
508.401.5012.31.73 × 10−40.1150.1011.13861.1023−3.19
519.101.5012.34.16 × 10−40.1190.0971.22681.1314−7.78
529.901.5012.35.88 × 10−40.1270.0931.36561.3455−1.47
5310.501.5012.38.05 × 10−40.1260.0901.40001.3577−3.02
5411.601.5012.39.16 × 10−40.1320.0851.55291.4620−5.86

Share and Cite

MDPI and ACS Style

Hou, J.; Zhang, C.; Wang, D.; Li, F.; Yu, Z.; Zhou, Q. Fixed-Bed and Mobile-Bed Resistance of Channels with Steep Gradients in Mountainous Areas. Water 2019, 11, 681. https://doi.org/10.3390/w11040681

AMA Style

Hou J, Zhang C, Wang D, Li F, Yu Z, Zhou Q. Fixed-Bed and Mobile-Bed Resistance of Channels with Steep Gradients in Mountainous Areas. Water. 2019; 11(4):681. https://doi.org/10.3390/w11040681

Chicago/Turabian Style

Hou, Ji, Chunze Zhang, Dan Wang, Feng Li, Zijian Yu, and Qin Zhou. 2019. "Fixed-Bed and Mobile-Bed Resistance of Channels with Steep Gradients in Mountainous Areas" Water 11, no. 4: 681. https://doi.org/10.3390/w11040681

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