Transcript Sediment characterization - University of Washington

OCEAN/ESS 410
15. Physics of Sediment
Transport
William Wilcock
(based in part on lectures by Jeff Parsons)
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
Sediment
Characterization
• 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-

Diameter,
D
Type of
material
-6
64 mm
Cobbles
-5
32 mm
Coarse Gravel
-4
16 mm
Gravel
-3
8 mm
Gravel
-2
4 mm
Pea Gravel
-1
2 mm
Coarse Sand
0
1 mm
1
0.5 mm
Medium Sand
2
0.25 mm
Fine Sand
3
125  m
Fine Sand
4
63 μm
Coarse Silt
5
32  m
Coarse Silt
6
16  m
Medium Silt
7
8 m
Fine Silt
8
4 m
Fine Silt
9
2 m
Clay
Coarse Sand
Sediment Characterization
Sediment grain smoothness
Sediment grain shape - spherical, elongated or
flattened
Sediment sorting
Sediment Transport
Two important concepts
•Gravitational forces - sediment settling out
of suspension
•Current-generated bottom shear stresses sediment transport in suspension or along
the bottom (bedload)
Shield stress - brings these concepts
together empirically to tell us when and how
sediment transport occurs
Definitions
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
m
n=
rf
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
2. Molecular and Eddy Kinematic Viscosities
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
3. Submerged Specific Gravity, R
rp - r f
R=
rf
rp
ra
f
Typical values:
Quartz = Kaolinite = 1.6
Magnetite = 4.1
Coal, Flocs < 1
Sediment Settling
Stokes settling
Settling velocity (ws) from the balance of two forces gravitational (Fg) and drag forces (Fd)
Fd µ ( Diameter ) ´ ( Settling Speed )
´ ( Molecular Dynamic Viscosity )
µ Dws m
Fg µ ( Excess Density ) ´ ( Volume )
(
´ ( Acceleration of Gravity )
)
(
)
µ r p - r f Vg µ r p - r f D 3g
µ means "proportional to"
Settling Speed
Fd = Fg
Balance of Forces
(
)
Dws m = k r p - r f D 3 g
ws
r
(
=k
p
ws
r
(
=k
p
k is a constant
)
- r f D2g
m
)
- rf rf
rf
1 RgD 2
ws =
18 n
m
Write balance using
relationships on last slide
D g
2
Use definitions of specific
gravity, R and kinematic
viscosity 
k turns out to be 1/18
Limits of Stokes Settling
Equation
1.
2.
3.
4.
Assumes smooth 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
and has a bigger effect than decrease in specific
gravity as a result of voids in floc.
Assumes laminar flow (ignores turbulence)
Shear Stresses
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
u
z
y
x
Outer region
z ~ O(d)
Intermediate layer
Inner region
d
is likely
the water
depth for
river flow.
is a few
tens of
meters for
currents
at the
seafloor
• 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)
Shear stress in a fluid
Shear stresses at the seabed lead to sediment transport
z
y
x
force
= shear stress = area
¶u
¶u
t = m = rfn
¶z
¶z
rate of change of momentum
=
area
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
¶u
¶u
t = m = rfn
¶z
¶z
Unfortunately the inner layer it is too thin for practical field
measurements to determine  directly
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:
¶u
t = rn e
¶z
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
Shear velocity u*
Sediment dynamicists define a quantity known as the
characteristic shear velocity, u*
¶u
u* = n e
¶z
¶u
2
t = rn e
= ru* = Constant
¶z
2
The simplest model for the eddy viscosity is Prandtl’s
model which states that
n e = ku* z
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
Velocity distribution of natural (rough)
boundary layers
From the equations on the previous slide we get
du
rku* z = ru*2
dz
Integrating this yields
u ( z) 1 z
= ln
u*
k z0
Þ ln z = ln z0 +
k
u*
u ( z)
z0 is a constant of integration. It is sometimes called the
roughness length because it is generally proportional to
the particles that generate roughness of the bed (usually
z0 = 30D)
What the log-layer actually looks like
z
0.1d
lnz
not applicable because
of free-surface/
outer-flow effects
z
0.1d
not applicable because
of free-surface/
outer-flow effects
lnz0
~30D
slope = u* /k
Slope =  /u*
= 04/u*
log
layer
log layer
~ 30D
U
Plot ln(z) against the mean velocity
u to estimate u* and then estimate
the shear stress from
viscous sublayer
U
Z0~ 30D
viscous sublayer
U
t = r f u*2
Shields Stress
When will transport occur and by
what mechanism
Shields stress and the critical shear stress
• The Shields stress, or Shields parameter, is:
t
qf =
( r p - r f ) gD
• 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/
Shields curve (after Miller et al., 1977)
- Based on empirical observations
Sediment Transport
No Transport
Initiation of Suspension
If u* > ws, (i.e., shear velocity > Stokes settling velocity)
then material will be suspended.
Suspension
No Transport
Transitional transport
mechanism. Compare
u* and ws
Bedload