Formulae of Sediment Transport in Steady Flows (Part 1)

*Shu-Qing Yang and Ishraq AL-Fadhly*

### **Abstract**

This paper makes an attempt to answer why the observed critical shear stress for incipient sediment motion sometimes deviates from the Shields curve largely, and the influence of vertical velocity is analyzed as one of the reasons. The data with *d50* = 0.016 � 29.1 mm from natural streams and laboratory channels were analyzed. These measured data do not always agree with the Shields diagram's prediction. The reasons responsible for the deviation have been re-examined and it is found that, among many factors, the vertical motion of sediment particles plays a leading role for the invalidity of Shield's prediction. The positive/negative deviations are associated with the up/downward vertical velocity in decelerating/accelerating flows, and the Shields diagram is valid only when flow is uniform. A new theory for critical shear stress has been developed, a unified critical Shields stress for sediment transport has been established, which is valid to predict the critical shear stress of sediment with/without vertical motion.

**Keywords:** critical shear stress, non-uniform flows, shields diagram, vertical velocity, decelerating/accelerating flows

### **1. Introduction**

The incipient motion of sediment is one of the most important topics in sediment transport. Generally, two methods are available in the literature to express quantitatively it, the shear stress approach and velocity approach [1]. The latter assumes that if the mean velocity excesses its critical velocity, then the sediment motion can be observed. The former used by researchers represents the force acting on a particle. Shields [2] is the earliest one who used the shear stress approach, or Shields number τ/(ρs-ρ)gd50 versus the Reynolds number, and he obtained a famous Shields curve to express sediment initiation. Francalanci et al. [3] interpreted the Shields number as the ratio of streamwise/vertical forces using the following form:

$$\tau\_{\ast\varepsilon} = \frac{\tau\_{\varepsilon}}{(\rho\_{\varepsilon} - \rho)\text{g}d\_{\sharp 0}} = \frac{4}{3} \frac{\tau\_{\varepsilon}\frac{\pi}{2}\left(\frac{d\_{\sharp 0}}{2}\right)^2}{(\rho\_{\varepsilon} - \rho)\text{g}\frac{4\pi}{3}\left(\frac{d\_{\sharp 0}}{2}\right)^3} \tag{1}$$

where *τ<sup>c</sup>* (=*ρu*\*c 2 ) is the critical shear stress for the median grain size of sediment, *d*50; *g* is the gravitational acceleration; *u*<sup>∗</sup> *<sup>c</sup>* is critical shear velocity; *ρ<sup>s</sup>* and *ρ* are the sediment and fluid densities, respectively. The shear stress exerted by the fluid

must be higher than the critical shear stress *τ<sup>c</sup>* to initiate sediment motion at the bed. Based on available experimental data, Shields in 1936 found that the Shields number *τ* <sup>∗</sup> depends on the particle Reynolds number (*R*<sup>∗</sup> ), i.e.,

$$R\_\* = \frac{u\_{\*\mathcal{L}} d\_{\mathbb{S}0}}{\nu} \tag{2}$$

where *ν* is the kinematic viscosity of the fluid.

The original Shields diagram has been reproduced and modified by many researchers. A comprehensive review has been done by many researchers [4, 5], in which some significant deviations of the observed critical shear stress from the standard Shields curve were observed. This has attracted extensive research by notable investigators, and some factors leading to the data scatter have been identified and discussed.

Some researchers believe that the definition of the incipient motion may cause the invalidity of Shields diagram, as the incipient motion depends more or less on the experimental observers'subjective judgment. To address this, criteria like "individual initial motion", "several grains moving" and "weak movement" has been introduced to express the incipient motion [6]. Subsequently, an error band has been included in the modified Shields diagram [7].

Other researchers attribute the large discrepancy to the stochastic nature of turbulence and sediment shape, its orientation, or exposure, protrusion [8–11]. It is natural to expect that when sediment is non-uniform, the critical condition is very difficult to determine, as the larger particles could move relatively easily than the finer one that is sheltered [12].

Over the past eight decades, the incipient motion has been extensively studied again and again [4], because the Shields diagram has been found invalid to predict the critical shear stress of sediment transport in some circumstances. The invalidity is not fully explained, some researchers ascribe it to sediment's characteristics, the other believe these deviations are caused by the flow conditions i.e., non-uniformity of flow [13].

Iwagaki [14] firstly linked the wide scatter in Shields diagram with flow's nonuniformity based on his observation: for the same sediment by the same experimenter, the observed critical shear stress in non-uniform flows largely deviates from that in uniform flows. Afzalimhr et al. [13] confirmed Iwagaki's results, they found experimentally that in decelerating flows, the critical shear stress is considerably below the Shields' prediction, and their experimental data are in complete disagreement with the Shields diagram. Other experimental researchers [15, 16] obtain similar results as Afzalimhr et al.'s [13] who claimed " … there is no universal value for *τ* <sup>∗</sup> ". Likewise, Buffington and Montgomery [4] also agreed "less emphasis should be given on choosing a universal *τ* <sup>∗</sup> ".

Some researchers try to explain the large discrepancy between predicted and measured critical shear stress by considering channel's characteristics, such as the channel shape and channel slope [16–22]. "the well-known Shields criterion is insufficient for large slope" was observed by Graf and Suszka [23], while Lamb et al. [24] comprehensively re-visited and examined almost all published datasets, and concluded that the critical shear increases with channel slope, this is totally different from the common sense that predicts increased mobility with increasing channel slope due to the added gravitational force in the downstream direction. But Chiew and Parker's experiments [17] in very steep channels show that the critical shear stress is decreased, contrary to Lamb et al.'s [24] conclusion.

Therefore, the brief literature review shows that Shields diagram cannot predict the critical shear stress well and there are many different potential causes for the

deviation, among them, it is necessary to clarify how the channel-bed slope and non-uniformity of flow affect the critical shear stress for sediment motion. The primary objectives of the present study are to


## **2. Theoretical considerations of influence of vertical velocity on the critical shear stress**

The author has been systematically investigating the role of vertical velocity on the mass and momentum transfer and has obtained a series of important and interesting conclusions [25–27]. It is found that omission of vertical velocity in our existing theorem of sediment transport makes many phenomena unexplainable. For example the presence of vertical velocity in non-uniform flows leads to the deviation of measured Reynolds shear stress from the linear distribution from the free surface to the bottom, consequently the upward velocity causes the positive deviation of velocity from the log-law or the wake-law is needed to express the velocity distribution, and the downward velocity results in the dip-phenomenon, or the maximum velocity is submerged and does not occur at the free surface as the log-law predicts. As the momentum and mass transfers are closely related to each other, it is interesting to investigate how the vertical motion affects sediment transport.

As a continuous effort, this study investigates the influence of upward/downward velocity on sediment incipient motion and the validity of Shields' diagram. **Figure 1** shows how a river flow interchanges with groundwater and the Darcy law

**Figure 1.** *The upward and downward vertical velocity generating from seepage face injection seepage.*

tells that vertical velocity is proportional to the hydraulic gradient, i.e., the suctions and injections inside groundwater can be expected in flood/dry seasons alternatively. The upward flow or injection flow may increase the sediment particles' mobility, or the required critical shear stress is reduced due to the "buoyant effect", which reduces the net settling velocity, mathematically

$$
\alpha' = \alpha - V\_b \tag{3}
$$

where *ω* = particle's settling velocity in still water and *ω*<sup>0</sup> = the net settling velocity subject to the vertical velocity of groundwater, *Vb*. The submerged weight in Eq. 1 can be represented by a drag force with the falling velocity *ω* in still water (*Vb* ¼ 0) as:

$$C\_d \pi \frac{d^2}{4} \frac{\rho \alpha^2}{2} = \pi \frac{4}{3} \left(\frac{d}{2}\right)^3 \mathbf{g}(\rho\_\varepsilon - \rho) \tag{4}$$

where *d* is the particle diameter, *Cd* is the drag coefficient.

If the upward velocity *Vb* of groundwater is so high and *Vb* = *ω*, the net settling velocity of the particle becomes zero, thus the particle is neutrally suspended, i.e., liquefaction state. This is often observed during earthquakes. In such case, saturated soil loses its strength and stiffness, it is natural that the Shields diagram cannot predict the particle's critical shear stress. Similarly, if the groundwater in **Figure 2** is downward, then the net falling velocity *ω*<sup>0</sup> should be higher than *ω*, the threshold critical shear stress should be unpredictable using the existing Shields diagram.

The above discussion clearly demonstrates that velocity *Vb* in a sediment layer may cause the invalidity of Shields diagram, which is supported experimentally by many researchers [28] who conducted experiments by observing the critical shear stress subject to injection and suction flows. Lu et al. [29] has reviewed these experimental results comprehensively. The influence of vertical motion on the critical shear stress has been discussed by many researchers, the parameters used to express the vertical motion include (i) the hydraulic gradient of seepage, e.g., Cheng and Chiew [30]; (ii) the pressure variation in flows [3]. But, there is no research available to investigate the role of time-averaged vertical velocity on the incipient motion of sediment transport.

The introduction of apparent sediment density is similar to Francalanci et al.'s treatment [3]. Instead of modifying the sediment density, they modified the water's density to eliminate the effect of pressure variation over time and space (like pressure induced by waves or bridge piers) on sediment's critical shear stress. Their results show that higher pressure yields higher "apparent water density", and lower pressure corresponds to lower "apparent water density". They found that the

**Figure 2.** *Schematic diagrams showing interaction of streamwise and vertical motions after Ladson (2008), p99.*

*Formulae of Sediment Transport in Steady Flows (Part 1) DOI: http://dx.doi.org/10.5772/intechopen.96705*

Shields number shown in Eq. 1 is actually the ratio of friction force in the streamwise direction (i.e., *τc*π*d*<sup>2</sup> /2) to the net force in the vertical direction, i.e., downward gravitational force (=*ρsg*π*d*<sup>3</sup> /6) minus the upward buoyant force (=*ρg*π*d*<sup>3</sup> / 6), so that the effective density of the particle is reduced from *ρ<sup>s</sup>* to *ρs-ρ*. When a particle is experienced in the environment with upward velocity *Vb,* additional upward force will be generated, and the problem is how to determine the additional upward force *Fvb* induced by the upward velocity.

Eq. 4 shows that the denominator of Shields number, i.e., (*ρs-ρ*)*gd* can be replaced by 3*Cdρω* 2 /4, this is why the parameter of settling velocity *ω* can be used in the study of sediment incipient motion. Some researchers believe that ω is a parameter for suspended load and it should not be used to express the sediment initial motion. The first one who uses *ω* to discuss the critical velocity is Yang [1], and this treatment significantly simplifies the problem. But this treatment has been blamed by many researchers who believe that the initiation problem does not involve any settling process. Now, Eq. 4 clearly shows that it is logical to express the incipient motion with sediment settling velocity.

For the case shown in **Figure 2**, if the upward velocity is zero, this is a static problem and the net force balance in vertical direction is expressed in Eq. 4. When the upward velocity is non-zero, this becomes a dynamic problem where the lift force *Fvb* should be included, i.e., submerged weight minus *Fvb* must be balanced by the drag force with settling velocity *ω'*.

$$C\_d' \pi \frac{d^2}{4} \frac{\rho \omega'^2}{2} = \pi \frac{d^3}{6} \mathbf{g}(\rho\_s - \rho) - F\_{vb} \tag{5}$$

where *Fvb* = *Cd,liftπd<sup>2</sup> ρVb* 2 /8. Let *ρ*<sup>0</sup> *<sup>s</sup>* <sup>¼</sup> *<sup>ρ</sup><sup>s</sup>* � *Fvb<sup>=</sup> <sup>π</sup>d*<sup>3</sup> *g=*6 and inserting the apparent density of sediment into Eq. 5, then the net force in the vertical direction can be alternatively expressed as:

$$C\_d' \pi \frac{d^2}{4} \frac{\rho \alpha \rho^2}{2} = \pi \frac{d^3}{6} \mathbf{g} \left(\rho\_s' - \rho\right) \tag{6}$$

In this study, apparent sediment density is introduced, and it depends on the vertical velocity of groundwater. Therefore, it is expected to have a relationship between the apparent sediment density and the settling velocity, similar to Einstein's relativity theory that the length/time depends on velocity if the light's speed is assumed to be constant. Therefore, the effect of vertical motion caused by pressure variation or seepage on sediment transport is eliminated after the introduction of apparent density. In other words, real lightweight particles motion must be the same as those with reduced settling velocity *ω*<sup>0</sup> in terms of mobility when both have the same settling velocity. Hence, the apparent density can greatly simplify the mathematical treatment for the complex *Fvb* induced by vertical motions. From Eqs. (4) and (5), the relationship between the modified settling velocity and the apparent density can be expressed by:

$$\frac{\rho\_s'-\rho}{\rho\_s-\rho} = a\left(\frac{a-V\_b}{a}\right)^2\tag{7}$$

where *α* is a coefficient (= *Cd* 0 /*Cd*) and *α* = 1 is assumed in this study to simplify the mathematical treatment. Because the drag coefficient depends on Re, this assumption is approximately correct if the particle size is very coarse or the Re does not change significantly with/without *Vb*. Eq. 7 tells that if *Vb* is equal to zero, then *ρ*0 *<sup>s</sup>* is the same as the natural sediment; if *Vb* is positive or upwards then *ρ*<sup>0</sup> *<sup>s</sup>* is less than the density of natural sediment *ρs*, and the particles behave like "lightweight sand" or plastic sands; if *Vb* = *ω,* then *ρ*<sup>0</sup> *<sup>s</sup>* is same as the density of water or similar to neutrally buoyant milk; if *Vb* is negative or downward, the higher apparent density of sediment behaves like heavy metals.

The vertical velocity *Vb* in **Figure 2** has the similar effect for the particles' stability as the buoyancy effect, i.e., the submerged weight of the particles is no longer *ρ<sup>s</sup>* � *ρ*, but *ρs*<sup>0</sup> � *ρ*, one may give the general expression of Shields number

$$\pi\_\*' = \frac{\pi\_\varepsilon'}{(\rho\_s' - \rho)\mathbf{g}d\_{\mathfrak{S}0}} \tag{8}$$

Substituting Eq. 7 into Eq. 8, one obtains:

$$\pi\_\*^{\prime} = \frac{\mathsf{r}\_c^{\prime}}{(\rho\_s - \rho)\mathsf{g}d\_{50}} \left(\frac{\alpha}{\omega - V\_b}\right)^2 \tag{9}$$

or

$$\pi\_\*^{'} = \frac{\pi\_\varepsilon^{'}}{(\rho\_\varepsilon - \rho)\lg\_{50}} \left(\frac{1}{1 - Y}\right)^2 \tag{10}$$

where *Y* = *Vb/ω*.

Eq. 9 or 10 generally expresses the influence of vertical velocity *Vb* on the critical shear stress. It is clear that *Vb* can be induced by seepage in the sediment layer, it can be also inferred that *Vb* can be estimated by the Darcy Law using the hydraulic conductivity and hydraulic gradient. Obviously, *Y* = 1 means that the particles can be suspended in water, i.e., liquefaction. If *Y* > 1, it means that particles flow in the upward direction with a net velocity of *Vb –ω*, this may have a devastating impact on dikes in flood defense as it may cause piping failure. In the following section, the analysis shows that the vertical velocity, *Vb* is ubiquitous in open channel flows, which is induced by non-uniform flows.

## **3. Influence of non-uniform flow on the critical shields stress**

Ideal uniform flow is very rare in natural conditions, flow rate and water depth/ channel width keep always changing, i.e., non-uniform, as shown in **Figure 3**. It is interesting to discuss how accelerating or decelerating flows generate the vertical

#### **Figure 3.**

*Non-uniform flows in open channel and the variation of water depth, in which u and v are mean velocities in x and y direction, respectively.*

*Formulae of Sediment Transport in Steady Flows (Part 1) DOI: http://dx.doi.org/10.5772/intechopen.96705*

velocity. To simplify the discussion, it is assumed that the flow rate is constant, i.e., d*h/*d*x* (6¼ 0). The 2-D continuity equation is:

$$\frac{\partial \overline{u}}{\partial \mathbf{x}} + \frac{\partial \overline{v}}{\partial \mathbf{y}} = \mathbf{0} \tag{11}$$

where *u* and *v* are the mean local velocity at any point in *x* and *y* directions, respectively. By integrating Eq. 11, one has

$$
\overline{\boldsymbol{\nu}} = -\int\_0^\mathcal{V} \frac{\partial \overline{\boldsymbol{u}}}{\partial \boldsymbol{x}} d\boldsymbol{y} \tag{12}
$$

∂*u*/∂*x* > 0 means accelerating, thus Eq. 12 tells that accelerating flows yield a negative or downward vertical velocity; but decelerating flows generates the positive or upward vertical velocity, i.e., ∂*u*/∂*x* < 0. Hence, the vertical velocity can be generated by in non-uniform flows.

For a channel with a constant width, its discharge per unit width can be expressed by:

$$\mathbb{Q}/b = \mathbb{U}h \tag{13}$$

where *Q* = discharge; *U* is the depth-averaged velocity, *b* is the channel width and *h* is the water depth. If Q/b could be constant in x direction, one has:

$$\frac{d(Uh)}{d\mathfrak{x}} = \frac{d}{d\mathfrak{x}}\int\_0^h \overline{u}d\mathfrak{y} = \mathbf{0} \tag{14}$$

The vertical velocity *vh* at the free surface can be obtained from Eq. 12 using Leibniz's rule, i.e.

$$d\overline{v}\_h = -\int\_0^h \frac{\partial \overline{u}}{\partial \mathbf{x}} d\mathbf{y} = -\frac{d}{d\mathbf{x}} \int\_0^h \overline{u} d\mathbf{y} + \overline{u}\_h \frac{dh}{d\mathbf{x}} \tag{15}$$

where *uh* is the horizontal velocity at the surface in the *x* direction. By inserting Eq. 14 into 15, one obtains:

$$
\overline{v}\_h = \overline{u}\_h \frac{dh}{d\mathbf{x}}\tag{16}
$$

Eq. 16 shows that *dh=dx* > 0, i.e., a decelerating flow generates *vh* > 0, but *dh=dx* < 0 or an accelerating flow yields the negative *vh*. Therefore, the vertical velocity in the main flows can also generate vertical velocity. Its interaction with groundwater can be obtained as *Vb = Vg + Vs*, or the groundwater velocity at the bed is jointly caused by Darcy velocity *Vg* and the vertical velocity caused by the main flow, *Vs*. Even the velocity on the solid–liquid interface may be very small, its importance for sediment transport should not be underestimated [3], and Eq. 16 shows that the vertical velocity has similar amplitude like the secondary current, i.e., about 1% of mean velocity.

Julien [5] replaced the Reynolds number in Shields' diagram by dimensionless particle diameter:

$$d\_{\*} = \left[\frac{\rho\_s - \rho g d\_{50}^3}{\rho}\right]^{1/3} \tag{17}$$

Similarly, *d\** needs modification by introducing the apparent density with the following form:

$$d'\_\* = \left[\frac{\rho'\_s - \rho}{\rho} \frac{\text{gd}\_{50}^3}{\nu^2}\right]^{1/3} \tag{18}$$

Inserting Eq. 7 into Eq. 18, one has

$$d\_\*^{\prime} = \left[\frac{\rho\_\sharp - \rho}{\rho} (\mathbf{1} - Y)^2 \frac{\mathbf{g} d\_{50}^3}{\nu^2} \right]^{1/3} \tag{19}$$

or

$$\frac{d\_\*'}{d\_\*} = (1 - Y)^{2/3} \tag{20}$$

Therefore, the empirical equation of Shields curve by Yalin and Silva [31] can be modified with the following form:

$$\pi\_\*^{\prime} = 0.13 d\_\*^{\prime - 0.392} \exp\left(-0.015 d\_\*^{\prime 2}\right) + 0.045 \left[1 - \exp\left(-0.068 d\_\*^{\prime}\right)\right] \tag{21}$$

For the fall velocity, many empirical equations are available in the literature. Julien [5] related *cd* in Eq. 4 with the particle diameter *d\** and obtained the following empirical equation:

$$\frac{\rho d\_{50}}{\nu} = 8 \left[ \sqrt{1 + 0.0139d\_\*^3} - 1 \right] \tag{22}$$

The incipient motion in uniform flows has been extensively investigated, but no one investigates the influence of vertical velocity on incipient motion, probably because this vertical flow may not large enough to induce discernible seepage, thus it is useful to estimate *Vb* using some measured parameters. The depth-average vertical velocity can be determined by,

$$V = U \frac{dh}{d\mathbf{x}}\tag{23}$$

where *U*, the average streamwise velocity, and both *U* and d*h*/d*x* are measurable parameters, thus Eq. 23 is convenient to use.

The vertical velocity is jointly induced by either the groundwater or the surface variation, the joint effect can be assumed as the proportional *V* and the nominal seepage velocity *Vs*, i.e., *λV* + *λsVs*, or:

$$V\_b = \frac{\lambda V + \lambda\_s V\_s}{1 - \varepsilon\_0} \tag{24}$$

where *λ* and *λ<sup>s</sup>* are the coefficients to relate *Vb* with the mean vertical velocity *V* and nominal seepage velocity (*Vs*) defined by Darcy (*Vs* = *ki*, *k* = hydraulic conductivity, *i* = hydraulic gradient), *ε<sup>0</sup>* = porosity of granular materials.

Generally in laboratory flumes, the second term of Eq. 24 is negligible (i.e., *Vs* = 0), but in natural streams both the river flow and underground water flow can generate the velocity at the river bed, thus two terms co-exist in Eq. 24.

## **4. Re-analysis of data on the original shields diagram**

To verify whether Eq. 21 is applicable to non-uniform flows, 329 data points are comprehensively compiled [13, 17, 23, 29, 32–40]. The hydraulic conditions of the used data are summarized in **Table 1**, and the experimental conditions are briefly outlined as follows:

Neil conducted his experiments in a flume 0.9 m wide and 5 m long by using sands with different particle sizes and densities [32]. Among the data sets, 11 data points are obviously above the Shields curve. White collected his data from a recirculating flume 6 m long and 0.3 m wide, uniform sediment was used with diameter between (0.016–2.2) mm [33]. The experimental datasets by Everts included 35 runs with size *d*<sup>50</sup> from 0.127 to 1.79 mm and specific gravity of 2.65, and 11 runs having *d*<sup>50</sup> from 0.09 to 0.18 and specific gravity of 4.7 [34]. **Figure 4** shows that almost all his data points are located below the Shields' prediction. Carling's data [35] were collected from a narrow natural stream and in a broad stream. Graf and Suszka measured the critical shear stress in a flume 16.8 m long, 0.6 m wide and 0.8 m high, gravel sediment with uniform size was used [23]. Shvidchenko and Pender [36] used a flume to study the effect of relative depth on the incipient motion of coarse uniform sediments. Gaucher et al.'s [40] experiments were conducted in a horizontal, rectangular glass walled flume with dimensions of 6 m long, 0.5 m wide and 0.7 m deep, different types of non-cohesive materials were used ranged from *d*<sup>50</sup> = 0.91 to 4.36 mm. Cheng and Chiew [30] investigated the influence of upward seepage on the critical conditions of incipient motion, the experiments were conducted in a horizontal flume 7.6 m long, 0.21 m wide and 0.4 m deep, with particle sizes of *d*<sup>50</sup> ¼ 0.63, 1.02 and 1.95 mm, and the seepage velocity (injection) was measured with a range between (0–0.0138) m/s. They found that the upward seepage reduces significantly the critical shear stress required by Shields curve. Kavcar and Wright [38] conducted experiments in a 7.5 m long, 0.6 m wide flume with both injection and suction seepage using sediment particle of *d*<sup>50</sup> =0.16, 0.5 and 1.2 mm and the observed value of seepage velocity, i.e. *Vs* is range between (�0.0026–0.00223) m/s. Liu and Chiew [29] examined the critical shear stress for sediment with *d*<sup>50</sup> = 0.9 mm subject to downward seepage with velocity between (�0.00314–0) m/s. Their glass-sided flume was 30 m long, 0.7 m wide and 0.6 m deep, they observed that the upward seepage (injection) decreases the critical shear velocity while the downward seepage (suction) increases it.

Nineteen flume experiments from Sarker and Hossain [37] are also included in **Figure 4**. They investigated the initiation of sediment motion under non-uniform sediment mixtures. Afzalimhr et al. [13] conducted experiments to investigate the effect of non-uniformity of flow on the critical shear stress in a channel (14 m long, 0.6 m width and 0.5 m depth), the sediment size of *d*<sup>50</sup> ¼ 8 mm was used for their observation. Different from Lamb et al's [24] prediction, their experimental data reveal that the value of critical shear stress is smaller than Shields' prediction by at least 50%. Similarly, Emadzadeh et al. [39] conducted experiments in accelerating and decelerating flow conditions, his flume was 14 m long, 0.6 m wide and 0.6 m deep. The sediment size used were *d*<sup>50</sup> = 0.8 and 1.3, 1.8 mm for a total of 72 data sets. The decelerating/accelerating flows were obtained by adjusting negative and positive bed slope (�0.7%, �0.9%, �1.25% and � 1.5%). It is found that the critical shear stress and Shields parameter for incipient motion in accelerating flow are higher than those predicted by Shields in uniform flow while their values in decelerating flow are considerably lower than that in accelerating flow.

These data mentioned are plotted in **Figure 4**, where the observed critical shear stress highly deviates from the standard Shields curve. All has been noticed and


*Sediment Transport - Recent Advances*

#### **Table 1.**

*Summary of experimental conditions by previous researchers.*

#### **Figure 4.**

*Measured critical shear stress versus* d\* *and its comparison with shields curve or Eq. 21 at Y = 0 (i.e., uniform flow) and 100% error band.*

commented by many researchers [4, 24]. The consensus is that this discrepancy cannot be simply attributed to measurement errors or methodological bias. In **Figure 4**, the three lines are the Eq. 21 (*Y* ¼ 0) �100% error band.

## **5. Dependence of critical shields stress on channel slope**

Many researchers have noticed that high channel's slope can cause the deviation of data from the Shields curve. For example Chiew and Parker [17] proposed that

**Figure 5.** *Dependence of critical shear stress on the channel slope.*

where *ϕ* = angle of streamwise bed slope, *θ* = angle of repose. Eq. 25 shows that the Shields number decreases with the increase of channel slope.

However, the formula given by Lamb et al. (2008) shows that the steep channel has a higher Shields number with the following form:

$$\pi\_\* = \exp\left[0.0249X^4 + 0.107X^3 + 0.199X^2 + 0.476X - 3.57\right] \tag{26}$$

where *X* = 0.407*ln*(142 *S*), and the slope *S* is in the regime 10�4< *S* <0.5.

**Figure 5** demonstrates the comparison of the measured data from **Table 1** and Eqs. 25 and 26. Obviously these equations do not agree the data points well. The measured *τ\** could be largely different even the same type of sediment and channel slope are used. Therefore, the invalidity of Shields prediction cannot be simply explained by the dependence of channel slope, and there are some physics inside for the discrepancy.

### **6. Seepage on critical shields stress**

**Figure 5** demonstrates that for the same particle size in the same channel slope, the data points behave largely different, which cannot be explained by any existing theory. Beyond other factors, Eq. 24 shows that the scatter could be induced by either groundwater or the main flow's non-uniformity, or both of them. The effect of seepage on the critical shear stress is discussed first, the experimental data [29, 30, 38] are showed in **Figure 6**.

The modified Shields number in Eq. 8 (i.e., that with seepage) will be the same as that obtained from the Shields curve if one uses both the apparent sediment density and the apparent critical shear stress (that with seepage), i.e.

$$\pi\_\*^{\prime} = \frac{\pi\_{\mathfrak{c}}^{\prime}}{(\rho\_{\mathfrak{s}}^{\prime} - \rho)\mathfrak{g}d\_{\mathfrak{s}0}} = \frac{\pi\_{\mathfrak{c}}}{(\rho\_{\mathfrak{s}} - \rho)\mathfrak{g}d\_{\mathfrak{s}0}} \tag{27}$$

Using Eq. 7, one obtained the ratio of critical shear stresses with/without *Vb* in the following form:

**Figure 6.** *Comparison of measured and predicted critical shear stress subject to seepage.*

*Formulae of Sediment Transport in Steady Flows (Part 1) DOI: http://dx.doi.org/10.5772/intechopen.96705*

$$\frac{\tau\_c'}{\tau\_c} = (\mathbf{1} - \mathbf{Y})^2 \tag{28}$$

**Figure 6** shows the critical shear stress predicted by Eq. 28 and the empirical factor *λ<sup>s</sup>* is found to be 8.5, the experimenters determined the critical shear stress without seepage using the Shields curve. The good agreement between the measured and predicted critical shear stress indicates that the introduction of apparent sediment density is acceptable.

Similarly, local scour by large vortices (e.g., scour holes around bridge piers) is not caused by higher velocity or higher boundary shear stress, but the upward velocity Y. The mechanism is similar to the helicopter whose rotor blades generate the upward velocity by "vortices". Consequently, low water pressure induces the seepage or upward velocity, large particles like stones/helicopter can be lift. One can easily infer the relationship between the upward velocity and "vortices" in front of an electricity fan. Likewise, by observing how tornadoes damages large particles like cars, houses on surface, one can easily concluded that the upward velocity or lift force is the cause, by no means the shear force.

## **7. Effect of non-uniformity of flow on the critical shear stress**

**Figure 4** shows that the Shields' curve could be totally invalid sometimes, these noticeable deviations imply that the non-uniformity of flow could affect the predictability of Shields curve, for example, Afzalimhr et al.'s [13] data points locate below the curve when the flow was decelerating, Emadzadeh et al.'s data points [39] were far from the Shields' prediction, and his data points were obtained from both decelerating and accelerating flows. Hence, the large deviations from Shields curve shown in **Figure 4** can be used to verify Eq. 24, i.e., the vertical velocity induced by flow's non-uniformity is responsible for the invalidity of Shields curve.

To confirm whether the invalidity of Shields curve is caused by the nonuniformity of flow, the data without seepage in **Table 1** are used, and the water depth variation *dh=dx* is calculated using the following formula:

$$\frac{dh}{dx} = \frac{\text{S} - \text{S}\_f}{\text{1} - \text{U}^2/\text{g}h} \tag{29}$$

where *dh=dx* is the water depth's variation, *S* and *S <sup>f</sup>* are the bed and energy slopes, respectively. Manning coefficient (*n*) can be assessed using the Strickler's formula:

$$m = \frac{d\_{\rm 50}^{1/6}}{21.1} \tag{30}$$

The energy slope *S <sup>f</sup>* in Eq. 29 can be determined from the Manning equation using the hydraulic radius *R,* i.e.*,*

$$\mathbf{S}\_f = \frac{n^2 U^2}{R^{4/3}}\tag{31}$$

In **Table 1**, the calculated *dh=dx* could be either negative or positive and the data points in **Figure 4** are included and replotted in **Figure 7**, where the data point is represented by the sign "+" if the obtained d*h*/d*x* is positive, otherwise the data point is marked by "-" for all negative *dh=dx* cases. **Figure 7** clearly shows that nearly all

#### **Figure 7.**

*The variation of water depth dh=dx has different values based on the influence of vertical velocity on the initial motion, where (*�*0.024<* **dh***=***dx** *<0.0526) for all data sets from Figure 3.*

data points above the Shields curve have "-" signs, indicating the flows were accelerating, whilst almost all data points below the Shields curve have the sign of "+", or decelerating. Therefore, the non-uniformity of flow can play an important role for the deviation of measured critical shear stress from the Shields curve. **Figure 7** reveals that the presence of vertical velocity is one of the main causes responsible for the deviation of observed critical shear stress from the Shields curves for these data, the accelerating flow enhances particles'stability, and decelerating flow enables sediment's mobility. In **Figure 7**, the calculated positive *dh=dx* ranges from 0.000237 to 0.0526 *and* � 0.024 to �0.00073. It remains necessary to investigate whether the higher *dh=dx* has the higher deviation, and its analysis is shown below.

## **8. Modification of shields diagram**

To examine whether data points without seepage shown in **Figure 4** can be expressed by Eq. 21, we can analyze the datasets without artificial seepage or with negligible groundwater effects, only those data are analyzed in which *Vb* is caused by the non-uniformity of flow in the main flow. Therefore, Eqs. 23 and 24 can be simplified as follows:

$$Y = \frac{V\_b}{a\nu} = \frac{\lambda U}{(1 - \varepsilon\_0)a\nu} \frac{dh}{d\mathbf{x}}\tag{32}$$

Experiments [13, 34, 39, 40] are analyzed first. They reported that their measured critical shear stress is lower than Shields' prediction. Besides, the datasets [32, 39] are examined; they claimed that higher values of critical shear stress were observed.

In these studies, the experimental data sets from non-uniform flows are plotted in **Figure 8** where the empirical factor λ is found to be 8.5 for both decelerating and accelerating flows. The comparison of the predicted and measured critical shear stress in **Figure 8** shows that the agreement is reasonably good. Better agreement can be obtained if *λ* is calibrated as a function of sediment gradation and shapes, turbulence. Here, the assumption is that sediment particle size is uniform and can be represented by *d50*.

*Formulae of Sediment Transport in Steady Flows (Part 1) DOI: http://dx.doi.org/10.5772/intechopen.96705*

**Figure 8.**

*Comparison of experimental results on threshold condition without seepage with Eq. 32.*

**Figure 9.**

*Influence of vertical velocity on critical shear stress, the solid line is the original shields curve (or* Y *= 0) and other lines are calculated from Eq. 21 with different* Y*.*

**Figure 9** shows the comparison of measured and the predicted critical shear stress for the datasets [23, 33, 35–37]. Obviously, the observed critical shear stress largely deviates from the solid line, i.e., Shields curve (*Y* = *Vs/ω* = 0), all data points can be covered by Eq. 9 or 10 when the parameter *Y* is used. In other words, **Figure 9** suggests that the scatter might be explained by variation of *Y*.

## **9. Discussion on slope's influence**

As mentioned, some researchers have found the dependence of the critical shear stress on the channel slope, but it is still an open question about the validity of

Shields curve, especially when the bed slope is large, thus it is worthwhile to discuss this dependence.

This study reveals that the deviation from the Shields curve could be caused by the vertical velocity, the Shields curve is approximately valid only when the flow is uniform, when the vertical velocity is almost zero. As the true uniform is very rare in laboratory or nature, thus it is understandable why Shields curve is invalid to express most of observed critical shear stress. Hence, one needs to answer whether the dependence of *τ\** on the channel slope is also caused by the flow's acceleration.

Eqs. 23 and 24 show that in almost all cases, there always exists the vertical velocity caused groundwater and flow's non-uniformity. Therefore, the widely observed dependence by Lamb et al. [24] may be also caused by the parameter *Y* (6¼ 0). Obviously, they assumed that the data used in their analysis were collected from uniform flows, thus these data can be used to compare the data with the Shields curve, and conclusion of the slope-dependence can be drawn. It is useful to examine this assumption by checking whether Lamb et al's data [24] are observed from uniform flows. Their data are listed in **Table 2**, in which only the laboratory data are included as their field data were certainly collected from non-uniform conditions. The last column of **Table 2** shows the length of flumes, and from it one can see that almost half of the flumes were less than 10 m. Kirkgöz and Ardiçlioğlu [41] measured the minimum length to form a uniform flow in a flume and found that a channel should be longer than 10 m as there is a transition zone from non-uniform flow to uniform flow. Even for those data from flumes longer than 10 m, the flow still could be non-uniform also when the parameter *dh=dx* is calculated using Eq. 6. Paola and Mohrig [42] suggest that uniform flow can only be assumed when the channel length is longer than h/S, if the water depth is 0.1 m, and slope is 1%0, this means that the channel length should be longer than 100 m. Therefore, one can conclude that it is likely that the data were generated in non-uniform conditions,


#### **Table 2.**

*Previously reported data selected from lamb et al. (2008).*

*Formulae of Sediment Transport in Steady Flows (Part 1) DOI: http://dx.doi.org/10.5772/intechopen.96705*

this may lead to the different interpretation of the dependence of critical shear stress on the channel slope. **Figure 3** shows the accelerating flows in steep channels, thus it is likely that the downward velocity increases particles'stability.

Chiew and Parker's data [17] is used as an example, their observation is opposite to Lamb et al's prediction [24], they also found the dependence of critical shear stress on the channel slope based on their own data. In their experiments, the channel slope was specially adjusted from -10<sup>o</sup> to 31<sup>o</sup> , their channel lengths used were 4 m and 2 m only. Obviously, their experiments were conducted in the nonuniform flow conditions as the 2 3 m length is too short to form a uniform flow. In other words, both conclusions drawn by Lamb et al. and Chiew and Parker [17, 24] are not very convincing as they did not check the parameter of *dh=dx*, and the data they used may be generated from non-uniform conditions.

## **10. Conclusions**

This paper investigates why the observed critical shear stress widely deviates from the Shields curve, its discrepancy or validity could be caused by many factors like sediment shapes, gradation, measurement errors, turbulence and channel-bed slopes. However, this study reveals that the vertical motion also plays an important role, and the vertical velocity could be induced by non-uniformity of flow and seepage turbulence alike. After re-examining 329 data points from the literature, the following conclusions can be drawn:


All in all, high horizontal motion can make a plane (a big particle) to fly, high vertical velocity can also make the same particle called helicopter to fly.

Two mechanisms are totally different. It is wrong to ascribe all sediment transport phenomena to the horizontal motion only, without considering the vertical motion.

## **Notations**

*b* = channel width *Cd* = drag coefficient; *d50 =* median size of sediment particles; *Fvb* = force induced by the vertical velocity; *g* = gravitational acceleration; *h* = water depth; *i* = hydraulic gradient; *k* = hydraulic conductivity; *n* = Manning coefficient; *Q* = discharge; *R\** = Reynolds number; *Sf* = energy slope. *U* = mean velocity; *u*\* = shear velocity; *u\*c* = critical shear velocity (*τ<sup>c</sup>* = *ρu*\*c 2 ); *u* and *v* = time-averaged velocity in the streamwise and vertical directions; *uh*, *vh* = horizontal and vertical velocities at the surface; *V* = vertical velocity; *Vb* = vertical velocity at the bed; *Vs* = nominal seepage velocity at the bed; *X* = 0.407ln(142*S*); *y* = distance normal to the wall; *Y* = *Vb*/*ω ε<sup>0</sup>* = porosity of granular materials *θ* = angle of repose. *λ* and *λ<sup>s</sup>* = coefficients; *ν* = kinematic viscosity; *ρ* = fluid density; *ρ<sup>s</sup>* = sediment density; *ρs' =* apparent density of sediment; *τ* = boundary shear stress;

*τ<sup>c</sup>* = critical boundary shear stress;


*Formulae of Sediment Transport in Steady Flows (Part 1) DOI: http://dx.doi.org/10.5772/intechopen.96705*

## **Author details**

Shu-Qing Yang\* and Ishraq AL-Fadhly School of Civil, Mining and Environmental Engineering, University of Wollongong, NSW, Australia

\*Address all correspondence to: shuqing@uow.edu.au

© 2021 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## **References**

[1] Yang, C.T. (1996). Sediment transport: theory and practice. McGraw-Hill, Sydney.

[2] Shields, A. (1936). Application of similarity principles, and turbulence research to bed-load movement, California Institute of Technology, Pasadena (translate from German).

[3] Francalanci, S., Parker, G. and Solari, L (2008). Effect of seepage – induced nonhydrostatic pressure distribution on bed-load transport and bed morphodynamic. Journal of Hydraulic Engineering, vol. 134, no. 4, pp. 378–389.

[4] Buffington, J.M. and Montgomery, D.R. (1997). A systematic study of eight decades of incipient motion studies with special reference to gravel-bedded rivers. Water Resources Research, vol.33, no.8, pp. 1993–2029.

[5] Julien, P. Y. (1995). Erosion and Sedimentation. Cambridge University Press.

[6] Kramer, H. (1935), Sand mixtures and sand movement in fluvial models, ASCE Transactions, 100(1909), 798–838.

[7] Beheshti, A.A., B. Ataie-Ashtiani, (2008). Analysis of threshold and incipient conditions for sediment movement, Coastal Engineering, 55(5), Pages 423–430,

[8] Wiberg, P.L., and Smith, J.D. (1987). Calculations of the critical shear stress for motion of uniform and heterogeneous sediments, Water Resource Research, vol.23, pp.1471–1480.

[9] Wilcock, P.R. (1987). Bed-load transport in mixed-size sediment, PhD. Thesis, MIT, Cambridge.

[10] Johnston, C.E., Andrews, E.D., and Pitlick, J. (1998). In situ determination of particle friction angles of fluvial

gravels, Water Resource Research, vol.34, no. 8, pp. 2017–2030.

[11] Kirchner, J.W., Dietrich, W.E., Iseya, F. and Ikeda, H. (1990). The variability of critical shear stress, friction angle, and grain protrusion in water worked sediments. Sedimentology, vol.37, pp.647–672.

[12] Garde, R. J. and Ranga Raju, (1985). Mechanics of sediment transportation and alluvial stream problems. 2nd ed. (New York: Wiley).

[13] Afzalimhr, H., Dey, S. and Rasoulianfar, P. (2007). Influence of decelerating flow on incipient motion of a gravel-bed stream. Sadhana, vol.32, no.5, pp. 545–559.

[14] Iwagaki, Y. (1956). Fundamental study on critical tractive force. Trans. Jan. Soc. Civil Engineering, no. 41, pp. 1–21.

[15] Andrews, E.D. and Kuhnle, R.A. (1993). Incipient motion of sand-gravel sediment mixture. Journal of Hydraulic Engineering, vol. 119, pp. 1400–1415.

[16] Dey, S. and Raju, U. (2002), Incipient motion of gravel and coal beds. Sadhana, vol.27, pp.559–568.

[17] Chiew, Y. M. and Parker, G. (1994). Incipient sediment motion on nonhorizontal slopes. Journal of Hydraulic Research, vol. 32, pp. 649–660.

[18] Andrews, E. D. (1994). Marginal bed load transport in a gravel-bed stream, Sagehen Creek, California. Water Resource Research, vol.30, pp. 2241–2250.

[19] Church, M, Hassan, M. A., and Wolcott, J. F. (1998). Stabilizing selforganized structures in gravel-bed streams. Water Resource Research, vol.34, pp. 3169–3179.

*Formulae of Sediment Transport in Steady Flows (Part 1) DOI: http://dx.doi.org/10.5772/intechopen.96705*

[20] Patel, P. L., and Ranga Raju, K. G. (1999). Critical tractive stress of nonuniform sediments. Journal of Hydraulic Research, vol.37, pp.39–58.

[21] Dey, S., and Debnath, K. (2000). Influence of stream-wise bed slope on sediment threshold under stream flow. Journal of Irrigation Draining Engineering, vol.126, pp.255–263.

[22] Mueller E R, Pitlick J and Nelson JM (2005). Variation in the reference Shields stress for bed load transport in gravel-bed streams and rivers. Water Resource Research. vol.41, W04006, (doi: 101029/2004WR003692).

[23] Graf, W.H. and Suszka, L. (1987). Sediment Transport in Steep Channels. Journal of Hydroscience and Hydraulic Engineering, vol.5, no. 1, pp.11–26.

[24] Lamb, M.P., Dietrich, W.E. and Venditti, J.G. (2008). Is the critical Shields stress for incipient sediment motion dependent on channel-bed slope? Journal of Geophysical Research, vol.113, F02008

[25] Yang S.Q. and J. W. Lee (2007). "Reynolds shear stress distributions in a gradually varied flow". Journal of Hydraulic Research, IAHR. 45(4), 462– 471.

[26] Yang S.Q. and A. T. Chow (2008). "Turbulence structures in non-uniform flows", Advances in Water Resources, 31, 1344–1351.

[27] Yang S.Q. (2009). "Velocity distribution and wake-law in gradually decelerating flows". J. Hydr. Res., IAHR, 47(2), 177–184.

[28] Rao, A. R., and Nagaraj, S. (1999). Stability and mobility of sand-bed channel affected by seepage. Journal of Irrigation Draining Engineering, vol.125, no.6, pp.370–379.

[29] Liu, X.X. and Chiew, Y.M. (2012). Effect of seepage on initiation of

cohesionless sediment transport. Acta Geophysica, vol. 60, no. 6, pp. 1778– 1796, DOI: 10.2478/s11600-012-0043-7.

[30] Cheng, N-S and Chiew, Y-M (1999). Incipient sediment motion with upward seepage. Journal of Hydraulic Research, vol.37, no.5, pp665–681.

[31] Yalin, M.S. and Silva, A.M.F. (2001). Fluvial processes. IAHR, Delft, the Netherlands.

[32] Neill, C. (1967). Mean-Velocity Criterion for Scour of Coarse Uniform Bed-Material. International Association for Hydraulic Research, pp. 46–54.

[33] White, S J. (1970). Plane bed thresholds of fine grained sediment. Nature, vol.228, pp. 152–153.

[34] Everts, C.H. (1973). Particle overpassing on flat granular boundaries. Journal of the Waterways Harbors and Coastal Engineering Division, vol. 99, pp. 425–439.

[35] Carling, PA. (1983).Threshold of coarse sediment transports in broad and narrow natural streams. Erath Surface Processes and Landforms, vol.8, pp.1– 18.

[36] Shvidchenko, AB. and Pender, G. (2000). Flume study of the effect of relative depth on the incipient motion of coarse uniform sediments. Water Resources Research, vol. 36, no.2, pp. 619–628.

[37] Sarker, LK and Hossain, MM (2006). Shear stress for initiation of motion of median sized sediment of no uniform sediment mixtures. Journal of Civil Engineering (IEB), vol.34, no. 2, pp. 103–114.

[38] Kavcar, P.C., and Wright, S.J. (2009). Experimental results on the stability of non-cohesive sediment beds subject to vertical pore water flux. Proc. World Environmental and Water

Resources Congress 2009: Great Rivers, vol. 342, pp. 3562–3571, DOI: 10.1061/ 41036.

[39] Emadzadeh, A, Chiew, YM and Afzalimehr, H. (2010). Effect of accelerating and decelerating flows on incipient motion in sand bed streams. Advances in Water Resources, vol.33, pp. 1094–1104.

[40] Gaucher, J, Marche, C and Mahdi, T-F. (2010). Experimental investigation of the hydraulic erosion of non cohesive compacted soils. Journal of Hydraulic Engineering, vol.136, no. 11, pp. 901–913.

[41] Kirkgöz, M. S. and Ardiçlioğlu M. (1997). Velocity Profiles of Developing and Developed Open Channel Flow. Journal of Hydraulic Engineering, vol.123, no. 12, pp. 1099–1105.

[42] Paola, C., and D. Mohrig (1996), Paleohydraulics revisited: Paleoslope estimation in coarse-grained braided rivers, Basin Research, 8, 243–254.

## **Chapter 2**
