1 / 28

15. Physics of Sediment Transport William Wilcock (based in part on lectures by Jeff Parsons)

This lecture/lab covers the characterization of sediment, the definitions of viscosity, sediment settling, and bottom boundary layers in sedimentary systems.

Download Presentation

15. Physics of Sediment Transport William Wilcock (based in part on lectures by Jeff Parsons)

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. OCEAN/ESS 410 15. Physics of Sediment TransportWilliam Wilcock (based in part on lectures by Jeff Parsons)

  2. Lecture/Lab Learning Goals • Know how sediments are characterized (size and shape) • Know the definitions of kinematic and dynamic viscosity, eddy viscosity, and specific gravity • Understand Stokes settling and its limitation in real sedimentary systems. • Understand the structure of bottom boundary layers and the equations that describe them • Be able to interpret observations of current velocity in the bottom boundary layer in terms of whether sediments move and if they move as bottom or suspended loads – LAB

  3. There are number of ways to describe the size of sediment. One of the most popular is the Φ scale. φ= -log2(D) D = diameter in millimeters. To get D from φ D = 2-φ Sediment Characterization

  4. Sediment grain smoothness Sediment Characterization Sediment grain shape - spherical, elongated, or flattened Sediment sorting % Finer Grain size

  5. Sediment Transport Two important concepts Gravitational forces - sediment settling out of suspension Current-generated bottom shear stresses - sediment transport in suspension (suspended load) or along the bottom (bedload) Shields stress - brings these concepts together empirically to tell us when and how sediment transport occurs

  6. Definitions

  7. 1. Dynamic and Kinematic Viscosity The Dynamic Viscosity μ is a measure of how much a fluid resists shear. It has units of kg m-1 s-1 The Kinematic viscosity ν is defined where ρf is the density of the fluid. ν has units of m2 s-1, the units of a diffusion coefficient. It measures how quickly velocity perturbations diffuse through the fluid.

  8. 2. Molecular and Eddy Viscosities Molecular kinematic viscosity: property of FLUID The molecular kinematic viscosity (usually referred to just as the ‘kinematic viscosity’), ν is an intrinsic property of the fluid and is the appropriate property when the flow is laminar. It quantifies the diffusion of velocity through the collision of molecules. (It is what makes molasses viscous). The Eddy Kinematic Viscosity, νe is a property of the flow and is the appropriate viscosity when the flow is turbulent flow. It quantities the diffusion of velocity by the mixing of “packets” of fluid that occurs perpendicular to the mean flow when the flow is turbulent Eddy kinematic viscosity: property of FLOW In flows in nature (ocean), eddy viscosity is MUCH MORE IMPORTANT! About 104 times more important

  9. 3. Submerged Specific Gravity, R Typical values: Quartz = Kaolinite = 1.6 Magnetite = 4.1 Coal, Flocs < 1 f

  10. Sediment Settling

  11. Settling Velocity: Stokes settling Settling velocity (ws) from the balance of two forces - gravitational (Fg) and drag forces (Fd)

  12. Balance of Forces Write balance using relationships on last slide k is a constant Settling Speed Use definitions of specific gravity, R and kinematic viscosity ν k turns out to be 1/18

  13. Limits of Stokes Settling Equation Assumes smooth, small, spherical particles - rough particles settle more slowly Grain-grain interference - dense concentrations settle more slowly Flocculation - joining of small particles (especially clays) as a result of chemical and/or biological processes - bigger diameter increases settling rate Assumes laminar settling (ignores turbulence) Settling velocity for larger particles determined empirically

  14. Boundary Layers

  15. Bottom Boundary Layers The layer (of thickness δ) in which velocities change from zero at the boundary to a velocity that is unaffected by the boundary • Inner region is dominated by wall roughness and viscosity • Intermediate layer is both far from outer edge and wall (log layer) • Outer region is affected by the outer flow (or free surface) δ is likely the water depth for river flow. δ is a few tens of meters for currents at the seafloor

  16. Shear stress in a fluid Shear stresses at the seabed lead to sediment transport force rate of change of momentum τ = shear stress = = area area

  17. The inner region (viscous sublayer) • Only ~ 1-5 mm thick • In this layer the flow is laminar so the molecular kinematic viscosity must be used Unfortunately the inner layer it is too thin for practical field measurements to determine τ directly

  18. The log (turbulent intermediate) layer • Generally from about 1-5 mm to 0.1δ (a few meters) above bed • Dominated by turbulent eddies • Can be represented by: where νe is “turbulent eddy viscosity” This layer is thick enough to make measurements and fortunately the balance of forces requires that the shear stresses are the same in this layer as in the inner region

  19. Shear velocity u* Sediment dynamicists define a quantity known as the characteristic shear velocity, u* The simplest model for the eddy viscosity is Prandtl’s model which states that Turbulent motions (and therefore νe) are constrained to be proportional to the distance to the bed z, with the constant, κ, the von Karman constant which has a value of 0.4

  20. Velocity distribution of natural (rough) boundary layers z0 is a constant of integration. It is sometimes called the roughness length because it is often proportional to the particles that generate roughness of the bed (a value of z0 ≈ 30D is sometimes assumed but it is quite variable and it is best determined from flow measurements) From the equations on the previous slide we get Integrating this yields

  21. What the log-layer actually looks like Slope = κ/u* = 0.4/u* lnz0 Plot ln(z) against the mean velocity u to estimate u* and then estimate the shear stress from Z0

  22. Shields Stress When will transport occur and by what mechanism?

  23. Hjulström Diagram

  24. Shields stress and the critical shear stress • The Shields stress, or Shields parameter, is: • Shields (1936) first proposed an empirical relationship to find θc, the critical Shields shear stress to induce motion, as a function of the particle Reynolds number, Rep = u*D/ν

  25. Shields curve (after Miller et al., 1977) - Based on empirical observations Transitional Sediment Transport Transitional No Transport

  26. Initiation of Suspension If u* > ws, (i.e., shear velocity > settling velocity) then material will be suspended. Suspension Transitional transport mechanism. Compare u* and ws Bedload No Transport

More Related