Research Topic Updated on Oct. 9, 2014

Size: px

Start display at page:

Download "Research Topic Updated on Oct. 9, 2014"

Jane Bell
5 years ago
Views:

1 Research Topic Updated on Oct. 9, 014

2 Sediment Settling Velocity F d W s 3 ga 1 d W s s s d d d a C F 1/ 1 gd C a a s d s Submerged weight = drag force One derives where ω s is the settling velocity, d is the diameter, and C d is the drag coefficient. C d and ω s have general relations: n n n e d N R M C 1/ / 1 Author M N Rubey (1933) Zhang (1961) Zanke (1977) Raudkivi (1990) Julien (1995) Cheng (1997) n n s D M N Nd M / 3

Wu and Wang (006) plotted C d ~ Re relation shown in the right figure, and established empirical formulas: M 53.5e 0.65 S P N 5.65e.5 S P n 0.7 0.

3 Wu and Wang (006) plotted C d ~ Re relation shown in the right figure, and established empirical formulas: M 53.5e 0.65 S P N 5.65e.5 S P n S P Wu and Wang (006) compared more than ten sediment settling velocity formulas, and found that the formulas of Zhang (1961), Hallermeier (1981), Dietrich (198), Cheng (1997), Ahrens (000), Jimenez and Madsen (003), and Wu and Wang (006) have comparable and reasonable reliabilities for predicting the settling velocity of naturally worn sediment particles (with a Corey shape factor of about 0.7). The average errors normally are less than 9%. 3

Sediment Deposit Porosity Han et al. (1981) d formula for uniform 1 0.55 pm d 4 1 sediment deposit: 0.095d d 0 0.3 0.

4 Sediment Deposit Porosity Han et al. (1981) d formula for uniform pm d 4 1 sediment deposit: 0.095d d e N 1 pk 1 p 11 p m k mk 3 d 0 d 1 mm d 1 mm where d is the sediment size in mm; d 0 is a reference size of 1 mm; and δ 1 is the thickness of the water film attaching to sediment particles, given a value of about mm. For a nonuniform sediment deposit, fine particles may fill the voids among coarse particles. If the filling is negligible, the Colby (1963) method can be used: Wu and Wang (006) modified the Komura (1963) formula as 0.1 p m ( d ) 4

5 Movable Bed Roughness Wu and Wang (1999, JHE, No. 1) 5

6 Predicted vs. Measured Flow Depths (4,376 data points) Error Range % of Calculated Flow Depths in Error Range Li-Liu van Rijn Karim Wu-Wang ±10% ±0% ±30%

7 Interactions among Nonuniform Sediments on the Bed Coarse particles have higher chances of exposure to flow, while fine particles are more likely sheltered by coarse particles. It is important to consider the effect of this hiding and exposure mechanism on nonuniform sediment transport. The widely used approach is to introduce correction factors into the existing formulas of uniform sediment incipient motion and transport: f ( d d or d d ) k k m k 50 where d k = size of kth sediment fraction, d m = mean diameter of bed material; d 50 = median diameter of bed material. 7

8 Hiding and Exposure Correction Factor Egiazaroff (1965) Formula: ck c log19 log(19dk / d m ) with the Shields number ck ck /[( s ) dk ] Ashida and Michiue (1971) Formula: ck c log19 log(19d ck d m / d Paker et al. (198) Formula: k k c50 / d d d m k 50 ) m d d k k / d / d m m m is an empirical coefficient between

9 Hiding and Exposure (Wu et al., 000) Exposure height of a particle on the bed, Δ e, is defined as the difference between the apex elevations of it and the upstream particle. If Δ e >0, the particle is at an exposed state; if Δ e <0, it is at a hidden state. For a particle with diameter d k in the bed surface layer, the value of Δ e is in the range between d k and -d j. Here, d j is the diameter of the upstream particle. Δ e is assumed to have a uniform probability distribution function: f 1 /( dk d j ), d j e dk 0, otherwise 9

10 The hidden and exposed probabilities of particles d i on the bed are: p hk N j1 p bj d k d j d j p ek N j1 p bj d k dk d j The criterion for sediment incipient motion proposed by Shields (1936) is then modified as ck d s k c p p ek hk m where Θ c =0.03, and m=0.6, as calibrated. 10

11 Wu et al. (000) Bed Load Formula p 3/ q bk n b / 1 gd n ck 3 bk s k. 11

12 Wu et al. (000) Suspended Load Formula p q sk / 1 gd ck 3 bk s k U sk

13 Comparison of Sediment Transport Capacity Formulas Using Brownlie s (1981) Data of Uniform Bed Load Error Ranges Percentages (%) of Calculated Transport Rates in Error Ranges Van Rijn Engelund & Fredse Bagnold Meyer- Peter & Mueler Wu et al. 0.8r r r Note: r = calculation / measurement. 13

14 Comparison of Sediment Transport Capacity Formulas Using Brownlie s (1981) Data of Uniform Bed-Material Load Error Ranges Percentages (%) of Calculated Transport Rates in Error Ranges Ackers & White Yang Engelund & Hansen SEDTRA Wu et al. 0.8r r r Note: r = calculation / measurement. 14

15 Comparison of Sediment Transport Capacity Formulas Using Toffaleti s (1968) Data of Nonuniform Bed-Material Load Error Ranges Percentages (%) of Calculated Transport Rates in Error Ranges Modified Ackers & W. Modified Engelund &H. Karim Modified Zhang SEDTRA Wu et al 0.5r r r Note: r = calculation / measurement. 15

16 Verification Scores of Nonuniform Bed-Load Formulas (Ribberink et al., 00) Formula Transport Rate Mean Diameter Average Wu et al E&H A&W + Day Parker (surface) A&W + P&S van Rijn MP&M + Egiaz MP&M + A&M MP&M + Hunziker (Factor n over/underestimation gives a score of 1/n) 16

17 Reliable Transport Capacity Formulas Single-sized total load Ackers-White (1973) formula is good for coarse sediment, not for fine sediment Laursen (1958) formula is good for fine sand and silt, not for coarser sediment Yang s (1973, 1984) formula has two sets of coefficients for sand and gravel Wu et al. (000) and Engelund-Hansen (1967) are good for wider size ranges Sing-sized bed load Wu et al. (000) formula Meyer-Peter and Mueller (1948) formula Single-sized suspended load Zhang (1961) formula Multiple-sized total load Wu et al. (000) formula is the top choice *: Ultimately, calibration using measurements is the most reliable approach. 17

18 Nonuniform Sediment Transport under Currents and Waves (Wu and Lin, 014) 18

19 Nonuniform Sediments in Rivers View of the Ooi River, Japan, showing sorting of gravel and sand on bars. From Ikeda (001) Bed armoring, River Wharfe, UK. From Powell (1998) 19

Photograph of a mixed gravel-sand beach at La Jolla, California, USA.

20 Nonuniform Sediments on Coastlines Gravel-sand beach at Kachemak Bay, Alaska, USA at low tide showing uniformly mixed sand and gravel on the backshore and large cobbles with surficial mud in the lower foreshore. Photograph taken by Peter Ruggiero. Photograph of a mixed gravel-sand beach at La Jolla, California, USA. Here, gravel and sand material are sorted zonally, with sand at the foreshore and gravel/cobble material along the backshore. Cobbles result from cliff erosion at the back of the beach. Photograph taken by K. Todd Holland. 0

21 Sediments near Mississippi River Estuary 1

Sediment Size and Sorting on Baeksu Tidal Flat, Korea Mean grain size (A, B) and sorting (C, D) data published in Yang et al. (005) for sediments of the Baeksu tidal flat, South Korea.

22 Sediment Size and Sorting on Baeksu Tidal Flat, Korea Mean grain size (A, B) and sorting (C, D) data published in Yang et al. (005) for sediments of the Baeksu tidal flat, South Korea. Maps derived from grab samples. Data for (A) and (C) are for summer 1998 and (B) and (D) for winter These data indicate temporal dependence on heterogeneity as sediments become coarser and more poorly sorted in the winter due to increased incident wave energy. A cross-shore gradient in mean grain size N1ϕ persists during both seasons. (Holland and Elmore, 008)

23 Bed Shear Stress due to Current Bed shear stress (grain + form): gn b, c 1/3 U c where U c is current velocity, h is flow depth, n is the Manning roughness coefficient. n can be specified using reference values or determined by 1/6 n h 18log(1 h / k s ) with s s s h k k k Grain roughness k s 3d 90 Single-sized 1.5d 90 Multiple-sized Form roughness k A Bed shear stress due to grain roughness 1 fu s r r r b, c c c with d 0.53 r r mm Raudkivi (1998) n f c n 45d r 1.5 gn h 1/ mm 3

24 Bed Shear Stress due to Waves Total and grain bed shear stresses (Josson, 1966) 1 f U 4 b, wm w w with Soulsby (1997) formula 1 fu 4 b, wm w w 0.5 f 0.37 A / k w w s 0.5 f 0.37 A / k w w s Total roughness k k k s s s Form roughness (Soulsby and Whitehouse, 005) Aw r A w A exp d d w 1.5 d r 0.15 r 1.0 exp 5000 A A U T w w w / 50 w 3.5 4

25 Second-order Stokes waves with r w =u w,max /U w -1 Asymmetric Waves cos cos u t U t r t w w w Bed shear stress due to grain roughness (Camenen, 00) U sin a sin a 1 w 13 c 1 c b, wm. on f w 1 rw rw 6 ac 6 ac U sin a sin a 1 w 13 t 1 t b, wm. off f w 1 rw rw 6 at 6 at 5

26 Bed Shear Stress due to Currents and Waves Grain shear stress 1 b, on f cw Uc Uwm, on U cu wm, on cos 1 b, off f cw Uc Uwm, off U cu wm, off cos( ) with f X f (1 X ) f cw u c u w X U U 0.5U u c c w 6

27 Nonuniform Bed Load Bed-load transport rates in onshore and offshore half cycles q p gd 3 b, on bk, on bk ( s / 1) k 1 ck q p gd 3 b, off bk, off bk ( s / 1) k 1 ck.. Resultant fractional bed-load rate T T q q q wc wt bk bk, on bk, off Tw Tw where T wc = onshore half cycle T wt = offshore half cycle T w = wave period 7

28 Nonuniform Suspended Load Fractional suspended-load transport rate b Uc qsk pbk ( s / 1) gd k 1 cri, k sk where U c = current velocity; and bed shear stress cos b b, c b, wm b, c b, wm 8

29 Summary of Uniform Bed-Load Data under Current and Waves (Camenen and Larson, 007) Author(s) Exp. facility Cycle No. of runs s d 50 (mm) U c (m/s) U w (m/s) T w (s) Abou-Seida (1965) OT Half Ahilan and Sleath OWT Half (1987) Horikawa et al. (198) OWT Half Kalkanis (1964) OT Half King (1991) OWT Half Sawamoto and OWT Half Yamashita (1986) Sleath (1977) OT Half Dibajnia and OWT Full Watanabe (199) Watanabe and Isobe (1990) Ahmed and Sato (003) Ribberink and Chen (1993) Ribberink and Al Salem (1994) Dohmen-Janssen and Hanes (00) Dohmen-Janssen (1999) OWT Full , , , , 6.0 OWT Full OWT OWT Full Full LWF Full , 9.1 OWT Full Ramadan (1994) OWT Full Ribberink (1995) OWT Full Katopadi et al. (1994) OWT Full Jassen et al. (1996) OWT Full Van der Hout (1997) OWT Full , Cloin (1998) OWT Full Hassan et al. (1999) OWT Full

30 Uniform Bed-Load under Waves only 30

31 Uniform Bed-Load under Combined Currents and Waves 31

32 Summary of Nonuniform Bed-Load Data under Current and Waves Author(s) Exp. Facil. Cycle Ahmed (00) OWT Full Hassan and Ribberink (005) O Dononghue and Wright (004) De Meijer et al. (00) Inui et al. (1995) Dibajnia and Watanabe (000) OWT OFT OWT Full Full Full No. of runs No. of sizes s d 50 (mm) U c (m/s) U w (m/s) T w (s) , , , , OFT Full , 5.0 OFT Full

33 Fractional Bed-Load under Current and Waves 33

34 Summary of Uniform Suspended-Load Data under Current and Waves (Camenen and Larson, 007) Author(s) Nielsen (1984) Bosman (198) Roelvink (1987) Steetzel (1987) Nieuwjaar (1987) Havinga (199) Grasmeijer (1995) Sistermans (00) Location Australian beaches, Australia DHL, The Netherlands Delft Hydraulic, The Netherlands Delft Hydraulic, The Netherlands Delft University, The Netherlands Vinje Basin, Delft, The Netherlands DUT, The Netherlands DUT, The Netherlands Exp. facility No. of runs d 50 (mm) h (m) U c (m/s) H sig (m) T w (s) Field Wave flume Large Scale Flume Large Scale Flume Flume Basin Flume Flume

35 Uniform Suspended-Load under Currents and Waves 35

36 Summary of Nonuniform Suspended-Load Data under Current and Waves Author(s) Jacobs and Dekker (000) Sistermans (001) Location DUT, The Netherlands DUT, The Netherlands Exp. facility No. of runs No. of sizes Flume 3 7 Flume 3 13 d 50 (mm) h (m) U c (m/s) H sig (m) T w (s)

37 Fractional Suspended-Load under Current and Waves 37

38 Near-Bed Suspended-load Concentration Near-bed suspended-load concentration is related to bed-load transport rate: c bk = q bk δu bk Bed-load layer thickness: max.0 d,0.5,0.01h 50 r Bed-load velocity: / 1 s u bk b gd k cri, k 0.5 c bk = δ p bk d k τ b τ cri,k cos b b, c b, wm b, c b, wm 38

39 Suspended-load Discharge Calculated using Near-Bed Concentr. Suspended-load transport is determined by integrating the product of current velocity and suspended-load concentration over the flow depth: c k c z bk q sk h z c udz The current velocity is determined with van Rijn s two-layer log law, and the suspended-load concentration is by (Williams et al., 1999) k L L where U *wcr is the time-averaged bed-shear velocity for ripple-scale roughness; U *wg is the peak wave-only bed-shear velocity for grain-scale roughness; and L is coefficient defined by Nielsen (199). wc, k wc, k wc s U U wc, k * wcr * wg 39

40 Transport rate by Near-bed Concentration of Uniform Suspended-Load 40

41 Fractional Transport rate by Near-bed Concentration of Nonuniform Suspended-Load 41

42 Summary Methods have been developed or selected to determine the bed shear stress under current and waves, in order to apply the existing Wu et al. (000) formula for multiple-sized sediment transport in coastal water. The enhanced Wu et al. bed-load formula can calculate the onshore and offshore bed-load transport rates separately and then derive the net transport rate, whereas the enhanced suspended-load formula calculates only the net transport rate due to the limited available data. The near-bed suspended-load concentration at the reference level is related to bedload transport rate, velocity and layer thickness. The developed formulas have been tested using the single-sized sediment transport data sets compiled by Camenen and Larson (007) and several sets of nonuniform sediment transport data collected from literature. More than half of the test cases are predicted within a factor of of the measured values. This accuracy is generally acceptable for sediment transport, particularly under current and waves. It is found that the measurement data of multiple-sized sediment transport under current and waves are quite limited. Only six sets of experiments on graded bed load and two on graded suspended load have been collected from literature. More experiment studies on graded sediment transport under current and waves are needed. 4

43 Publications Related W. Wu and S. S.Y. Wang (1999). Movable bed roughness in alluvial rivers, J. Hydraulic Eng., ASCE, 15(1), W. Wu, S.S.Y. Wang, and Y. Jia (000). Nonuniform sediment transport in alluvial rivers, J. Hydraulic Research, IAHR, 38(6), W. Wu and S. S. Y. Wang (006). Formulas for sediment porosity and settling velocity, J. Hydraulic Eng., ASCE, 13(8), W. Wu (007), Computational River Dynamics, Taylor & Francis, UK, 494 p. W. Wu and Q. Lin (014). Nonuniform sediment transport under non-breaking waves and currents. Coastal Engineering, 90, 1 14,

A NEW MODEL FOR SEDIMENT TRANSPORT RATE UNDER OSCILLATORY SHEET FLOW CONDITIONS

A NEW MODEL FOR SEDIMENT TRANSPORT RATE UNDER OSCILLATORY SHEET FLOW CONDITIONS Proceedings of the 7 th International Conference on Asian and Pacific Coasts (APAC 2013) ali, Indonesia, September 24-26, 2013 A NW MODL FOR SDIMNT TRANSPORT RAT UNDR OSCILLATORY SHT FLOW CONDITIONS X.