Coasta 2
Coasta 2
COASTAL BOTTOM
BOUNDARY LAYERS AND
SEDIMENT TRANSPORT
Peter Nielsen
World Scientific
ABOUT
THE BOOK
This book is intended as a handbook
models.
COASTAL BOTTOM
BOUNDARY LAYERS AND
SEDIMENT TRANSPORT
ADVANCED SERIES ON OCEAN ENGINEERING
Series Editor-in-Chief
Philip L-F Liu
Cornell University, USA
Forthcoming titles:
Water Waves Propagation Over Uneven Bottoms
by Maarten W Dingemans (Delft Hydraulics, The Netherlands)
Ocean Outfall Design
by lan R Wood (Univ. Canterbury, New Zealand)
Tsunami Run-up
by Philip L- F Liu (Cornell Univ.), Costas Synolakis (Univ. Southern California),
Harry Yeh (Univ. Washington) and Nobu Shuto (Tohoku Univ.)
Physical Modules and Laboratory Techniques in Coastal Engineering
by Steven A Hughes (Coastal Engineering Research Center, USA)
Advanced Series on Ocean Engineering — Volume 4
COASTAL BOTTOM
BOUNDARY LAYERS AND
SEDIMENT TRANSPORT
Peter Nielsen
The University of Queensland
Yo World Scientific
Singapore « New Jersey * London * Hong Kong
Published by
World Scientific Publishing Co. Pte. Ltd.
P O Box 128, Farrer Road, Singapore 9128
USA office: Suite 1B, 1060 Main Street, River Edge, NJ 07661
UK office: 73 Lynton Mead, Totteridge, London N20 8DH
All rights reserved. This book, or parts thereof, may not be reproduced in any form
or by any means, electronic or mechanical, including photocopying, recording orany
information storage and retrieval system now known or to be invented, without
written permission from the Publisher.
ISBN 981-02-0472-8
981-02-0473-6 (pbk)
a 4 Woe re
‘ | = gstied
mts ™AAdte FA Ww ; ’ iad 7 i :
is
s
a ; go 7 -
| 7 :
ain > pi
en —
~ i :
- } ie ~ - : a
tone” Galligan
gue iesealable;
= =
ai }
~ >= —
=) i
awiod
SS
Sitos yey:
id
, Take et
forge © Teaag Shutrinebindlie |
+ vipa somes’ ite bei,
© () 2 1)a,
iKsh. €2:.0)-M74
G8. 3S Qhey
PREFACE
vii
Preface
equations for turbulent, steady flow. The new equations, however, consider
explicitly the periodic velocity component as well as the time-average and the
random component. These equations, like the Reynolds equations for steady
turbulent flow, are very helpful in analysing the mechanisms of momentum
transfer in combined wave-current flows.
Section 1.2 presents a review of wave boundary layers. It is followed by a
discussion of oscillatory boundary layer models in Section 1.3.
Section 1.4 discusses the wave-generated boundary layer currents including
boundary layer streaming and the surf zone undertow.
Wave-current boundary layer interaction is the subject of Section 1.5.
Experimental data on wave energy dissipation and on wave boundary layer
structure in the presence of currents indicate that currents of typical, relative
strengths have very little effect on the wave boundary layer structure. On the other
hand, the current boundary layer structure is usually strongly influenced by the
presence of waves. A flexible, traditional modelling framework is suggested for
currents in the presence of waves. It is realised however, that the traditional eddy
viscosity based models are theoretically unsatisfactory. Therefore, a new type of
model is suggested in Section 1.5.9.
Chapter 2 considers the initiation of sediment motion and the quasi-steady
processes of sediment transport over flat sand beds under waves. Emphasis is
placed on discussing similarity and differences between these processes and steady
bed-load and sheet-flow sediment transport.
Small scale coastal bedforms are discussed in Chapter 3 together with the
hydraulic roughness of beds of loose sand under waves.
Chapter 4 deals with the motion of suspended sediment particles in coastal
flows. Emphasis is placed on the important mechanism of sediment trapping by
vortices. At the same time, it is shown that a pure wave motion has but negligible
effect on sediment settling.
Suspended sediment distributions is the topic of Chapter 5. Section 5.2
presents the nature of suspended sediment concentrations in coastal flows through
a review of the existing experimental data. The concept of pickup functions for
suspended sediment in unsteady flows is discussed in Section 5.3. Practical pickup
function formulae are derived, partly in analogy with pickup functions from steady
flow and, partly based on time-averaged suspended sediment concentrations in
unsteady flows. Section 5.4 discusses the mechanisms which distribute the
suspended sediment through the water column. It is shown that pure gradient
diffusion is an unsatisfactory model of suspended sediment distributions in most
natural flows. The reason is that the sediment distribution process often includes
mechanisms which have mixing lengths of the same order of magnitude as the
Viii
Preface
a Cp, Nabe
ix
‘a es la (ddwll weow yi MALL aie oy
esa on
a sake ee mae e |ai ie i
ubouink 2). deowton (Dewy tot 6. eae alm
» eons blind av Seta vaisingeee seer ae
oe antbauthibs hore} Los oaoyawihad deat to eisai abut saiest20207]
Pe apcny ; vro?Z pi trantsveds2 ‘abana ieu Nib a7
yeni
ior} view! f . met ‘hy Ww. es a lb-qcaagvpga pee ii
ia
| oO ans +403 body
be Mltnagen igre | aasiplaoep. 4, gk cle
phazs Klienbli Habouy.1g add 2o Abo ny
ee a tsb wen
‘iw’ sfaideoup. Lore wrtiaam alias 1Metalig hak
88 la 9
ous ienculdLee \ coiged SOO
tga
Yo
Bb Mh ”
9,ganeht
nett,
ce been
aegl sxor
Mslors Ly, I | dae mB ober ih ’ sg ne is a m
7 Ppa, u i |
oT Give Rs
iti ad
sol all oat,
“ Na
CONTENTS
NOTATION XVii
xil
Contents
xiii
Contents
XV
Contents
REFERENCES 299
AUTHOR INDEX 309
SUBJECT INDEX 313
xvi
NOTATION
xvii
Notation
XVill
Notation
XX
Notation
XXi
Notation
XXil
Notation
General operators
yt Time-average, steady component.
~
Periodic component, Eq (1.1.17).
Random component, first derivative or skin friction.
>
A
Peak value.
Subscripts
b Break point value.
c Pertaining to current in combined flow.
0 Deep water quantity.
rms Root mean square value.
s Pertaining to sediment.
w Pertaining to waves in combined flow.
XXili
oe hire xemy tnodle bod peated iiidde
ot Cored
Broo ve)
the ellenae (Dip ‘olgmooto tasemagré
‘ Ll q ire anes : ai
-
—
i a“ng 084
cant Mid meee
Ir o.tot ee
wv bwatit ie nate Meee take 3S
mugenSze on
Ti bena
a
eens nite
pouraqo leans
i ae
Ls
REAR 1 Saas UE NEE -
a ‘ Ut
~ ANEAD) eee ‘ z. )
JvOra “Kne
a _ aol au nige
A 10 iinet ta . 1986 ig
ot, Trig ria
‘usin 6 s
i. LT }
iv hornets Sissengttic
ve) | i = cots sty sre th 3oe
Wee ' at | b l 2 abfeae wives,
vin intl”? ; Tike overaged ulio@ smep f
fy} (Mi. T* supe Saleh PM.
4 Mere) Pos 6 otthea
ese te
; Vy iogires ag! mictoon x = i s WT,
ct _ average ob @y2d,
¢ fq ~~“ aa och walege ¢
1.1.1 Introduction
With respect to sediment transport modelling, the most important part of the
flow is the bottom boundary layer through which the main flow influences the bed.
The present chapter therefore describes the types of boundary layer flows
which occur in the coastal environment.
For simplicity, the bottom topography is taken for granted and considered to
be stationary. The effects of the boundary layer flow on mobile sand beds, and the
resulting topographical changes will be considered in the following chapters.
The bottom boundary layer is intuitively defined as the layer inside which
the flow is significantly influenced by the bed. There are various ways of defining
the thickness 6 of this layer in quantitative terms. However, in general terms the
boundary layer thickness obeys the formula
5 < WT (1.1.1)
where vz is the eddy viscosity and T is the flow period. Thus, for fixed eddy
viscosity, the boundary layer for tidal flow with a period of 12 hours will be
approximately sixty six times thicker than that of a ten second wave motion, and
while the tidal boundary layer thickness is often equal to the water depth the wave
boundary layer is generally only a small fraction of the depth. See Figure 1.1.1.
The ability of a certain flow component to transport sediment is mainly a
function of the shear stresses it induces at the bed. Therefore, since thinner
boundary layers mean larger shear stresses for a certain free stream velocity, the
Chapter 1: Bottom boundary layer flow
waves will tend to dominate over the tide with respect to sediment entrainment and
bedform formation. However, because of the “one step forwards and one step
backwards” nature of the wave motion, the currents may still be very important
transporters of wave entrained sediment unless the sediment concentrations vary
periodically in step with the wave motion.
wave induced
velocity
Figure 1.1.1: The tidal
boundary layer is much
thicker than the wave
boundary layer. So, even if
the tidal velocity Uide is
much larger than the
wave-induced velocity
amplitude U(z) near the
surface, the waves will
dominate the situation at
the bed.
taken for granted and generally expressed as a simple harmonic function of time.
N
Directi f
a ee WR ee
10
K
Bs Kik
—s (oe) K
— u(x,z,t)
u=u(z,t) n::
oc
®::
aQ::
U:-
ox
::
° O::
@-:-:
O::
OF:
Ore
r/A
roughness,
Relative
0.01
0.001
103 104 105 106 107
Reynolds number, (A2 w/v)
Figure 1.1.3: An overview of the ranges of the Reynolds number A’ w/v and the
relative roughness r/A under likely field conditions and in some previous laboratory
experiments. The range shown for observed relative roughness values over loose
sand beds corresponds to all of the laboratory data of Carstens et al (1969) and
Lofquist (1986). If data from artificially flat beds at very low flow velocities are
excluded, the lower limit of r/A is 0.08. Legend, K: Kemp & Simons, flume (1982,
1983); S: Sleath, flume and tunnel (1982,1987); V: van Doom, flume and tunnel
(1981,1982); J: Jonsson & Carlsen, tunnel (1976); D: Jensen, tunnel (1989).
u=u(x, z,t) is based on considerations of the data of Kemp and Simons (1982,
1983) and those of van Doorn (1981, 1982). The former show pronounced
differences between the velocity structure over the bedform crests and over the
troughs, while the latter show reasonably uniform behaviour through most of the
boundary layer.
The value of 0.5 is only meant as an indication. The choice of an upper limit
of r/A for the application of horizontally averaged models does of course, in the
end, depend on the amount of detail one needs to consider.
The observed range from laboratory experiments of the relative roughness of
natural sand beds, is derived from all of the dissipation data of Carstens et al
(1969) and the bed shear stress data of Lofquist (1986). Values of r/A are
obtained by applying Equation (1.2.22) "in reverse" to the authors’ friction factors.
Rippled sand beds generally gave r/A > 0.2, while values of r/A below
0.08 were found only for flat beds with very little sand movement, see Section 3.6.
The lower limit of observed r/A-values of 0.08 for natural (flat) sand beds is
surprisingly high. It corresponds to roughness values of the order a hundred grain
diameters, and there is some doubt regarding its interpretation in relation to
effective sediment transporting stresses, as discussed in Section 2.4. The roughness
of natural sand beds is discussed in detail in Section 3.6.
Very large relative roughness values may apply in nature over rocky areas
and over reef platforms.
The message from the presently available data for oscillatory flow over loose
sand is, in essence, that the range of Reynolds numbers and relative roughness
values likely to be found under field conditions are
corresponding to a limiting condition of, for example, (A, T) = (0.3m, 6s), and
most wave flumes and basins is fairly limited it can be difficult to obtain
adequately large values of the Reynolds number A’w/v for modelling boundary
layer phenomena.
Therefore, a different type of apparatus was suggested by Lundgren and
Soerensen (1956), and later applied in many studies, namely the oscillating water
tunnel. This is essentially a large U-tube where the flow is driven by a piston in
one of the vertical legs, see Figure 1.1.4.
The horizontal test section of such tunnels can be several metres long so that
Reynolds numbers in excess of 10° can be obtained. The orbital motion in the test
section of the U-tube differs from real wave-induced flow by being totally uniform
in the x-direction (OUeo/dx = 0) and by having no vertical orbital motion, but these
dissimilarities are often of no concern.
RSS |
| Test Section |
Ree Se ee eee
| 2- 10m |
Figure 1.1.4: An oscillating water tunnel can be used to model many characteristics
of the wave boundary layer, but the vertical velocities of a wave motion and the
horizontal variation of the free stream are not modelled.
Due to the considerable cost of large wave flumes and tunnels a somewhat
simpler and cheaper type of facility was introduced by Bagnold (1946) and later
applied with modifications by many others. Instead of oscillating the bulk of the
water, Bagnold oscillated the bed in it’s own plane through otherwise still water,
see Figure 1.1.5.
u=AW cos wt
Oscillating Plate
Figure 1.1.5: Oscillating plates or trays can be used for modelling oscillatory
boundary layers and sediment transport. However, distortion of the inertia/pressure
forces on sediment particles is of some concem.
The flow over such an oscillating bed is, for all practical purposes, similar to
the velocity defect
in a tunnel, but the two types are in some respects dissimilar as far as sediment
transport experiments are concerned.
The reason is that the forces resulting from fluid pressure gradients on a
resting sand particle are exaggerated on the oscillating bed. A grain at rest in the
reference frame which follows the tray with velocity —v..(t) will experience an
inertial force of magnitude
Fi = p(s+ Cw VE (1.1.4)
where p is the fluid density, s is the specific sediment density, V_ is the particle
volume and Cy is the added mass coefficient. However, the corresponding force
on a particle on a fixed bed in moving fluid is only
dice=
Fy = p(l+Cn)V— (1.1.5)
with typical values of Cy=0.5 and s = 26 the difference amounts to a factor
two. Thus, discrepancies can be expected for sediment transport phenomena
where the pressure force Fp is significant in comparison with the fluid drag force.
Further discussion of the forces on sediment particles is given in Section 2.1.
ou, ou, Ou 1 op fe a
= +i.
Qi i novon fw edz = =e
agent ro
where u and w are the velocities in the x and z directions respectively, p is the
fluid density, p isthe pressure, and v_ is the kinematic viscosity of the fluid.
The flow inside the boundary layer can often be considered to be essentially
horizontal (w = QO). The equation of motion can then be further simplified to
where L_ is the wave length. Hence the relative importance of the convective
term originating from uo is given by the factor 27A/L.
The second criterion for horizontal uniformity in the boundary layer is that
non-uniformities introduced by individual roughness elements should be restricted
to a layer which is considerably thinner than the boundary layer itself , see Figure
1.1.6.
Since the scale of the disturbances introduced by the individual roughness
elements is the bed roughness r, this may be expressed by 6/r>> 1 which
corresponds to A/r>>1, since 6/r_ isan increasing function of A/r.
More precise information about the boundary layer thickness relative to the
roughness height can be gained from Figure 28 of Sleath (1987). The Figure
shows, that over a bed of three dimensional roughness elements, the ratio of
boundary layer thickness to roughness size is given approximately by
u=u(z),w=0O
Pi eel treet bre, vbr Pte ' OUTER LAYER
Figure 1.1.6: Even when the free stream velocity is horizontally uniform, the flow
near the bed may be non uniform throughout an inner layer of thickness similar to the
boundary layer thickness.
Ques
apalT_ _ Op
mde (1.1.11)
OG
ae Bho
ee 0 (1.1.14)
The tilde denotes the periodic component, which is the phase average over
several (N) wave periods minus the time average
N
(zit) }= - Y' u(z, t+ JT) — u(z) C117)
fl
We note that with these definitions we have ¥=x'=x =O , and
S| 3! ?= xy’ = xy’ = xy’ = 0, while XY = Xy-—Xy. Inserting the expressions
—— wie —_—fos) into the “horizontal” Navier Stokes equation (1.1.6) gives
2ay Grit)
(ut+u+u’) Deri) ++2ootiew’)
++Fada 57 w+w+w’) = 0 1.1.20
(1.1.20)
t)
Li + SGI) + 5 Gm + 27 +5, +5,
2 ii) + Su? + UV) +5 ew’)
= -1
—+ op on ou, ou |
vj —[+—+— 94.121)
yes & ay? az”
By collecting all terms which represent momentum flux in the vertical
direction under the label T this gives, by analogy with Equation (1.1.15),
because V = 0.
If the flow is horizontally uniform we have w =0 by continuity, so that the
second term on the right hand side of Equations (1.1.22) and (1.1.23) disappear. In
order to derive the equivalent equation to Equation (1.1.21) for the periodic flow
component, we take phase averages (with zero mean in accordance with Equation
(1.1.17)) on both sides of Equation (1.1.18) and make use of the continuity
equation as above. Then the following equation results
ou a Set
+22 Gm + Harri + 2 GH+TM
Ete
Se +5. Gin)4+7 ie)
Cig
+ Sl cae +S es
+ SV) We) (1.1.24)
wr an 2~ 2 2~
= ES Bal
ous ouHou
p ox ax" dy? az
Again we are particularly interested in the terms which represent vertical flux
of horizontal momentum, so we extract all such terms and put them under the
collective label T
_ ou
T= py) —puw -puw —p (un) — 9 (u’w’) ¢1.1.25)
The velocity product terms in Equations (1.1.22) and (1.1.25) play the same
role for wave-current motion as the Reynolds stresses —p u’w’ do for turbulent
steady flows. These velocity product terms will be used to discuss the total shear
stresses and the eddy viscosity concept for oscillatory flows and for combined
wave-current motions in Sections 1.2.8-1.2.9, and Section 1.5.3.
1.2.1 Introduction
The wave boundary layer is intuitively defined as the layer close to the
bottom, where the wave-induced water motion is noticeably affected by the
boundary. This layer is normally very thin, generally a few millimetres over a
smooth, solid bed and a few centimetres over a flat bed of loose sand. Bedforms,
like ripples, may change the structure of the boundary layer by introducing strong
rhythmic vortices. Hence, the boundary layer over sharp crested ripples will extend
to a height of four or five ripple heights, or a total of about 50 centimetres under
field conditions.
Although the water motion induced by natural waves is not simple harmonic,
it is instructive and useful to study the simple harmonic, oscillatory boundary
layer, which corresponds to Uo = AW cos wt, and use it as an approximation to
natural wave boundary layers.
The total picture of the velocity variation in such a layer is at first sight very
complicated because both amplitude and relative phase change with the distances
from the bed, see Figure 1.2.1. However, in the following section we shall see that
a much clearer picture can be obtained by applying simple transformations to the
data.
Figure 1.2.1: Instantaneous velocities u(z,t) [cm/s] plotted against elevation from
top of the roughness elements. Numbers on the curves refer to the phase of the free
stream velocity Ue(z,t) = AW@cosM@t. The measurements were made in an
oscillating water tunnel by Jonsson and Carlsen (1976). Note that the velocity near
the bed tums before the free stream velocity and that the velocity amplitude is largest
in the range Scm<z<J0cm notat z3~-,
and leads ue by 45°, but in a turbulent flow the variation with time is much more
complicated. The deviation from the simple harmonic behaviour increases with the
ratio between the bed roughness, r and the semi-excursion A. Thus, for small r/A
the variation of t(t) is still quite smooth and rather like a simple harmonic. This is
the case for the measurements of Jonsson and Carlsen (1976) Test 1 where r/A
was only 0.008, see Figure 1.2.2. In this case the phase shift between u.. and the
bed shear stress is still fairly well defined and we see that it is somewhat smaller
than the 45° of smooth laminar flow.
Ug (t)
Figure 1.2.2: Time variations of the bed shear stress t(o0, t) for rough turbulent flow
over relatively small roughness elements r/A=0.008. After Jonsson and Carlsen
(1976).
For flow over fully developed sand ripples, the ratio r/A is of the order of
magnitude one, and the flow near the bed is dominated by the rhythmic formation
and release of strong vortices. Lofquist (1986) measured instantaneous values of
t(o, t) under such conditions. Figure 1.2.3 shows some of his results and we see
that the behaviour is completely different from that of a simple harmonic and also
from that of sin(@t—)lsin(@t-@)| which has been assumed in several
“theoretical” studies.
Another typical feature of oscillatory boundary layers is the “overshoot” near
the bed. That is, there are elevations where the velocity amplitude U exceeds
Aq@ , see Figure 1.2.4.
Figure 1.2.3: Shear stress variation over fully developed sand ripples. For these
experiments A was fixed at 0.24m and the ripple length at 0.32m. u.(t) varied as
sin wt.
The velocity overshoot occurs because the velocity defect w..(t) — u(z, t)
has the nature of a damped wave which alternately adds to and subtracts from the
free stream velocity u..(t) , see Section 1.2.4.
which says that the shear stress at the level z is equal to the fluid density times
the total acceleration defect above z.
Sleath (1987) discussed shear stresses and related quantities in great detail for
turbulent oscillatory flows. One of his most striking findings was that the total
shear stress calculated from Equation (1.2.1) was, for his experiments, about a
factor ten larger than the periodic, turbulent Reynolds stress defined by
TR(Z, 1) = —p(u'w’
where uv’ and w’ are horizontal and vertical, turbulent velocity fluctuations
respectively, see Figure 1.2.13.
Hence, as Sleath pointed out, the turbulent fluctuations u’ and w’ are “mere
spectators” to the oscillatory boundary layer processes. Their contribution to the
momentum transfer is totally overshadowed by the analogous contribution — p u w
from the periodic velocity components. This latter contribution in turn was found
to agree rather closely with the values of t calculated from the defect integral
ye 2 ig
For the shear stress amplitude |t (z, t)|_ the data of Sleath (1987), and of
several previous authors, show that it tends to decrease roughly exponentially with
increasing distance from the bed. The vertical decay scale is approximately the
boundary layer thickness 6 , ie,
a Stee re (1.2.2)
W rms
which is in good agreement with the data except within the above mentioned layer
of constant turbulence intensity immediately above the bed. The thickness and
importance of this constant-intensity-layer is generally small, but may be judged
from Figure 4 of Sleath (1987).
An example of the distribution of turbulence intensities in a turbulent,
oscillatory boundary layer is shown in Figure 1.5.7.
With respect to the time-dependent turbulence intensities ws and
(w’)", both Sleath (1987) and Jensen (1989) found that the amplitudes were
greatest near the bed and that the peak values occured later at higher elevations.
The general picture is one of parcels of turbulence propagating upwards from the
bed. On the basis of his data, Sleath derived the expression
@ S05
Wr = 221
(1.2.3)
for the speed of vertical convection of turbulence, where the boundary layer
thickness, 5,05 may be estimated from Equation (1.1.9).
Uo(t) = Aw ei (1.2.4)
The real part A@ cos wt represents the physical velocity, see Figure 1.2.5.
For laminar flow the shear stress is proportional to the local velocity gradient
and the fluid viscosity
F) du
Hit te) == VIG
=—(U-—Uo) Vz (1.2.6)
2G
in terms of which the equation of motion takes the form of the diffusion equation
Awcoswt
Figure 1.2.5: The complex velocity Aw e/®! has the constant modulus Aw and
moves around a circle with angular velocity @. The real part Aw cos wt which
represents the physical velocity oscillates between A@® and —AQ.
dD
Va =
a°D
> Ew (1.2.8)
This is easily solved by separation of variables and assuming that the velocity
defect function has the form
co
Since the velocity defect must vanish for z—> 2, Ay must be zero for all n, and
the boundary condition at the bed where the velocity itself vanishes
u(0,t) = Uso(t)-
Aw D(o, t) = Awe’ - AwY Bre” = 0 (1.2.12)
gives B, = 1 and By = 0 forn # 1. Hence, the complete solution is
The complex velocity defect function D(z) gives the velocity different
phases as well as different magnitudes at different elevations, see Figure 1.2.6.
The shear stress distribution is found by inserting the expression (1.2.14) for
the velocity field into Newton’s formula (1.2.5) which gives
This shows that the shear stress magnitude decreases exponentially away
from the bed with a decay length scale of V2v/w . The bed shear stress is
: cs I
1(0,t) = pvAa(1 +i) Va/2v & = pVav Awe*4) (1.2.15)
so we see that the bed shear stress in smooth, laminar oscillatory flow leads the
free stream velocity by m/4 radians or 45°.
The value 2 = p V@v Aw for the maximum bed shear stress shows that,
with the definition 7 = +P tw (Aw)? from Jonsson (1966), the wave friction
factor for smooth, laminar flow is
2
PRE aii
APSA
Ge (1.2.16)
Equation (1.2.15) for the bed shear stress corresponding to u(t) = AM Aes
is also valid for the individual harmonic components of an arbitrary free stream
velocity u(t), so the frequency response function for the bed shear stress (o,f)
from input u.0(f) is
fw = fw oe 4] (1.2.20)
0.1
0.01
Numbers refer to A/d 99
x Smooth
0.001
Figure 1.2.7; Measured values of f, from oscillatory flow over flat beds of fixed
sand grains, after Kamphuis, 1975. The Nikuradse roughness for these sand beds was
taken to be 2do0,
Somewhat simpler formulae are adequate for the fully developed, rough
turbulent regime where A’w/v — © while r/A_ is finite, and for smooth
conditions (r/A — QO).
In Figure 1.2.7 the crosses correspond to experiments with a smooth bed. It
can be seen that, for Reynolds numbers up to about 3- 10° the smooth bed wave
friction factor is well described by
D
(i= ot (1.2.16)
VA2w/ Vv
which is the theoretical result for smooth laminar flow. The range
310° < A’w/v < 610° isa transition zone where fw increases. It is followed
by the fully developed, smooth turbulent regime where fy decreases again,
although at a slower rate than in the laminar regime. Justesen (1988) suggests
Kamphuis’ rough bed data are arranged in groups with fixed relative
roughness in Figure 1.2.7, the numbers referring to the ratio doo/A. . Note that
data points towards the upper right hand corner tend to fall on horizontal lines
corresponding to constant fy for constant do9/A.. We say that these data, for
which fy = fw(r7A), are in the fully developed, Rough Turbulent Regime. Several
authors have proposed formulae for fw in this regime, of which the most
commonly used is
\ 0.194
Sw = exp [5.213 g = 597 7h) (1.2.22)
0.5
0.1
0.05
0.01
0.005
0.001 0.01 0.1 1
Figure 1.2.8: Observed wave friction factors for rough turbulent flow. + : Riedel
(1972) with r = 2d90, 0: Kemp and Simons (1982), Sleath (1987), AJensen (1989),
e: Jonsson and Carlsen (1976).
for oscillatory flow. However, he found that a single constant would not apply
throughout the whole range of relative roughness. Instead he suggested two values,
namely 200 for fairly rough beds (r/A 2 0.01) while the value 70 from steady
flow applies asymptotically for r/A — 0.
While the peak value * is reasonably described by the formulae above, the
time dependence of the bed shear stress is not so well understood in general.
Two examples of the variation of (o,f) with time were given in Figures
1.2.2 and 1.2.3 for two cases of turbulent flow with different relative roughness.
These show how (o,f) is fairly similar to the smooth laminar solution (1.2.16) for
small r/A. That is, the variation is almost sinusoidal but the phase lead, relative
tO Uco , is less than the “laminar” 45°. Jonsson and Carlsen’s Test 1 shows a
phase lead of about 25°.
Lofquist (1980, 1986) reported comprehensive measurements of t(0,f) over
natural sand beds. The sequence shown in Figure 1.2.3 corresponds to a constant
value of the velocity amplitude Aw with the peak acceleration Aw growing
from top to bottom. The curves show that larger Aw” leads to a greater number
of more pointed peaks in (0, t). Lofquist explained how the peaks are related to
the growth and release of lee vortices.
=.= 2 3
De = 3 Pfe (Aw) (1.2.26)
1.2.26
fe = 3k = 4f = aii fi (1.2.28)
Tu.
fe "ah pias
Figure 1.2.9: The friction factor fy plotted against the energy dissipation factor
fe from measurements by Lofquist (1986) over rippled sand beds.
Jonsson (1966) used a different type of definition. He defined the top of the
boundary layer as the minimum elevation where u(z, ft) equals u(t) when the
latter is maximum. Jonsson’s boundary layer is quite thin i e, his definition
corresponds to Oj = V2v/@ for smooth laminar flow or approximately half
of 805.
Kajiura (1968) worked with the displacement thickness defined as
1 co
This again is a fairly thin boundary layer since for smooth laminar flow it
corresponds to S¢=Vv/w andhence ID; (84)l = exp (- V2/2) = 0.49.
However, the displacement thickness has the advantage that it is related in a
simple way to the other important boundary layer parameter, namely, the wave
friction factor. Their interrelation stems from the fact that the above definition for
da_ is very similar to the integrated momentum Equation (1.2.1) which in turn
defines the friction factor through Equation (1.2.18). Combining these equations
leads to
This formula is exact for simple harmonic flows with the form
u(z, t) = [1—D,(z)] u(t) (Equation 1.2.13). However, it also provides a useful
estimate of 5g in general. This is important because the two major data sets of
oscillatory boundary layer flow over natural sand beds, Carstens et al (1969) and
Lofquist (1986) provide only fw (or fe) but no details of the velocity distribution.
Therefore, 5g (or any other vertical scale for the boundary layer) cannot be
determined directly from the data, only via Equation (1.2.30).
The practical limit for measuring boundary layer structures with present day
technology lies around the level where the velocity defect is one percent of the free
stream velocity amplitude. This level, 501 relates to the other measures of
boundary layer thickness as follows
Ss = tx /@ = Vfw/2A (1.2.32)
which, for flow of the form u(z,f) = [1—Dy (z)] u(t) , is related to the
displacement thickness by 5« = 64 Gi72y and since fy is generally of the
order 0.2 orless we see that 65x is fairly large compared to the other 5-values.
Aw
p2 (u-ux) ae (1.1.12)
contains two unknowns namely u and _ tT. It can therefore only be solved when the
relation between these two is known or assumed to have a certain form. We know
the relationship for laminar flows where it is given by Newton’s formula (1.2.5)
but for turbulent flows it is not well understood.
Many schemes (turbulence closure schemes) have been suggested for getting
around this problem and the simplest of these involves the use of the eddy
viscosity concept, which was first introduced by Boussinesq. It is defined by
analogy with Newton’s formula for laminar shear stress i e
tT = pv (122.33)
= a +V (1.2.34)
az
in terms of the Reynolds stress -puw and it is a general feature of steady
turbulent flows that the first term dominates over the molecular viscosity except
inside the laminar sublayer.
For a uniform, oscillatory flow with zero net flow (u = u+u’), the
analogous expression for the eddy viscosity becomes (in accordance with Equation
f12))
Vt
_= &
En =
_ =-(aw) on-'w) + V (1.2.35)
Paz oz
31
AND SEDIMENT TRANSPORT
Chapter 1: Bottom boundary layer flow
dt _ dio | A (v al (1.2.36)
are of are or
and it is the definition which has been used in all existing eddy viscosity based
models of oscillatory boundary layer flow.
Except for laminar flow, and for one particular class of turbulent flows,
which will be discussed in Section 1.2.9, the eddy viscosity for oscillatory
boundary layers is a function of the distance from the bed and it can generally not
be considered constant (at least not a real valued constant) over time for fixed z.
The reason is that the shear stress and the velocity gradient tend not to be zero at
the same phase see Figure 1.2.11.
Figure 1.2.11: Time dependence of shear stress and velocity gradient, both phase
averaged. From Jonsson & Carlsen (1976) Test 1, elevation above ripple crest: 4.5em.
This was realised previously by Horikawa & Watanabe (1968), Jonsson &
Carlsen (1976) and Sleath (1987). An example of the corresponding variation of
the eddy viscosity derived directly from the definition (1.2.33) is shown in Figure
1.2.12. This eddy viscosity takes on negative and even infinite values.
Sleath (1987) studied the turbulent stress component
°8
° S
)
Eddy
(cm?/sec
Viscosity
ani) (1.2.37)
as well as the total given by Equation (1.2.35). He found that the turbulent stress
component contributed only about ten percent of the total and it showed a very
different variation with phase, see Figure 1.2.13. In contrast, the
contribution —(u w) from the phase-averaged flow corresponded rather closely to
the total shear stress derived from the momentum equation (1.2.1).
Sleath (1987) also studied the time average of ze and of the total v;. He
found that Vzr was “firmly negative” over several roughness heights near the bed,
then went positive and continued to increase away from the bed.
For the time average v; of the total eddy viscosity he found that it was
generally positive and tended to increase with distance from the bed. However, the
general magnitude was significantly smaller than the equivalent for steady flow.
That is, in terms of the average friction velocity u* = (\q|/ p)°° he found
where K_ is analogous to von Karman’s constant « (= 0.4) for steady flow, but
is seen to be smaller by a factor of 3 to 4. These results were obtained for relative
roughness of the order r/A = 0.01.
0.016
~
Ae
1/2p(Aw)?
: ua!
0.008 4X: 1/2 (Aw)2
wt (degrees)
Figure 1.2.13; Reynolds stress -p(uw’) compared to the total shear stress T
calculated from Equation (1.2.1). After Sleath (1987), Test 4.
which, with the velocity gradient and the shear stress written in the complex forms
Thus, treating the oscillatory boundary layer flow as simple harmonic with a
correspondingly simple harmonic shear stress (or analysing only the fundamental
harmonic modes) leads to an eddy viscosity .(z) which is a function of z but
not of ¢.
However, making provision for the shear stress to be out of phase with the
velocity gradient requires v, to be complex with an argument @y equal to the
phase shift between 1 and du;/0z.
z/6q
We note that if one chooses to use real valued cosine functions instead of the
complex exponentials in (1.2.40) and (1.2.41) then, the result of the phase shift
between shear stress and velocity gradient is a real valued eddy viscosity which
varies strongly with time and takes both negative and infinite values as indicated
by Figure 1.2.12.
z/byq
Z (cm)
0.5 1.0
Figure 1.2.16: Variation of lvil/de for a range of flow conditions. The straight
line corresponds to Equation (1.2.43). x : Jensen (1989) test 10, (A, T, r, 505 ) =
(3.1m, 9.72s, 0.0mm, 60mm); + : Jensen (1989) test 12, (A, T, r, 5.05 ) = (1.58m,
9.72s, 0.84mm, 50mm); e: Jonsson & Carlsen test 2 (A, T, r, 505 ) = (1.79m, 7.2s,
63mm, 150mm) ; 0 : van Doorn (1982) M-series (A, T, r, 505 ) = (0.33m, 2.0s,
21mm, 26mm); «: van Doorn (1982) S-series (A, T, r, 5.05 ) = (0.10m, 2.0s, 21mm,
14.5mm).
was given by
v4
D(z) = exp[-(1 +) ] (1.2.45)
V2vV/@
for smooth laminar flow. Thus, the real and imaginary parts of the complex
logarithm In(D;) = In |ID,| + i Arg Dj, are identical:
Z
In ID\| ArgD, = 7a, (1.2. 46)
100
MILLIMETRES
IN
ELEVATION
7
om 0:2 0:4 1-0 2 4 10
Figure 1.2.17: ArgD, and In\D,i for Test SOORA from van Doom (1982).
A=10cm, r=2.1cm, T=2.0s. The identity of ArgDi and /n|Djl and their
shared proportionality with z shows that the boundary layer structure is analogous
to smooth laminar flow. A deviation from this "smooth laminar analogy" is visible
below z = 3mm, but at these elevations, the assumption of horizontal uniformity (u =
u(z,t)) is not valid anyway.
Zz —
v4
ArgD, = In ID\| — ia Z1 pe \2v,;/@ (1.2.47)
1.3.1 Introduction
The existing models for oscillatory boundary layers fall into two broad
physical categories, namely horizontally uniform models where u = u(z,f) and
models which take into account the horizontal variability of u(x,z,t) between
crests and troughs of the bed roughness elements. The latter group is by far the
smallest although realistic modelling of flow over the commonly observed sand
ripples obviously calls for models which can describe localised vortex formation.
Longuet-Higgins (1981) developed an essentially inviscid model based on
conformal mapping of a sharp crested ripple profile and the discrete vortex
method. The model yields a reasonable description of (o,f), but the results are
numerical only. Also, because of the assumed sharp ripple crest the vortex
shedding is continuous throughout each half cycle while observations of flow over
natural ripples tend to show that separation only occurs during flow deceleration.
A different numerical model was developed by Sleath (1982). His ripple
profiles had rounded crests and separation did not always occur. This model also
reproduced (o,f) quite well, but the fact that it is laminar calls for some caution
when extrapolating the results.
Recently, a different type of numerical solution was presented by Blondeaux
and Vittori (1990) which gives a good reproduction of the vortex shedding and the
resulting sediment clouds over vortex ripples.
Horizontally uniform models are much simpler. They can however, only be
literally valid at elevations which are well clear of the top of the roughness
elements, i e for z>>r. Hence, unless 8 >> r they are at most relevant as
descriptions of the horizontal average of the flow. Estimates of the ratio 5/r can
be found from Equation (1.1.9).
wt (degr.) =
0° 15° 30° 45° 60° 75° 90°
y ‘ Sv hy lob lewtermn %
+ —p
x “Pe ° r
1000 ° % ‘5 ape
100 §‘ 2 €
ott. A ° ; ” 4
10 / / 7 js 7 / /
ut
1
0: 5".0" 5-0 (e) () fe) 0
7 i
Figure 1.3.1: Dimensionless velocity profiles for smooth turbulent oscillatory flow at
different phases of the free stream velocity Awsin @t ; A = 3.1m, T 9.725,
yt =z. /v, ut = ii(z,t)/i(t). After Jensen (1989).
From this expression we note that In|D,| and Arg Dy, are identical
quantities and both proportional to the distance from the bed
Dw){e)
ELEVATION
MILLIMETERS
IN
2 4 6 10 20 40 60 100
Figure 1.3.2: The analogies between smooth, laminar and rough, turbulent
oscillatory flows become visible when real and imaginary parts of In D; are plotted
together against z. Data from Jonsson & Carlsen (1976) Test 2, (A,T,r) = (1.79m,
7.28,63mm).
D, = exp[-(1+i) ie (1.3.6)
where the two parameters z; and p are derived as shown in the figure above.
From presently available data it appears that Equations (1.2.13) and (1.3.6)
describe w(z,t) reasonably well for turbulent and transitional oscillatory
boundary layers with relative roughness greater than about 0.01. We may note
here, that all the available measurements of friction factors over beds of loose sand
by Carstens et al (1969), and Lofquist (1986) indicate relative roughness values
well within this range, see Section 3.6.
For smooth and almost smooth turbulent flows, like those investigated by
Jensen (1989), there is no longer an identity between Arg D; and In\DjI, see
Figure 1.3.3.
100
60
40
Figure 1.3.3: Arg D and In\D\\ derived from the measurements of Jensen (1989)
Test 13, A = 3.1m, T = 9.7s, r = 0.84mm. Even for these, almost smooth flow
conditions the formula (1.3.11) predicts the defect magnitude ID;! with great
accuracy. For thes experiments :0.09 Vr A = 0.0046m.
For very low relative roughness values Arg D; is smaller than /n |D| down
to 40%, but the two quantities do behave in a very similar fashion throughout the
boundary layer.
As mentioned in connection with Figure 1.2.17, there is a very close analogy
between u1(z,t) from fairly rough turbulent flows (r/A > 0.06) and the smooth
laminar solution. For such flows one finds p = / and hence
D; = exp[-(1 ay (1.3.7)
and as mentioned in Section 1.2.9 this flow structure corresponds to a constant
eddy viscosity vr = 0.5 zi.
107!
Arg D, <In |D,] ; Arg: In |D,|
-2 2
10 z,=0-09VrA
10-
10-4 103 102 107) 10
Figure 1.3.4: The simple formula (1.3.10) is valid over the full range where
horizontally uniform (u = u(z,t)) velocity models make sense. G : van Doorn (1981),
strip roughness, *: Sleath (1987), sand roughness, @ : Jonsson & Carlsen (1976) strip
roughness, 0: Jensen (1989), sand roughness.
Z1
ss= 21olited
ARNE (1.3.8)
>
Nid ~ CL 3.9)
>|> <|2
SS
for the complete range 0.01 < r/A <0.5 where Equation (1.3.6) is applicable,
see Figure 1.3.4.
In fact, this.formula can also be used for prediction of the magnitude of the
velocity defect for almost smooth turbulent flows where Arg D; and InID,I are
no longer identical. Thus, for such flows we still have
ie , the parameter z; maintains its "laminar value", the Stokes length, for smooth
oscillatory flows, even when they become fully turbulent.
With respect to the second parameter p the message from the experimental
data is not quite as clear, see Figure 1.3.5.
However, it can be seen that the power p varies from unity for very rough
flows to about //3 for smooth turbulent flows. It was noted by Nielsen (1985) that
the p-value of approximately 1/3 for smooth turbulent flows corresponds to
maximum energy dissipation for fixed 2}.
The fact that p=1 for r/A2 0.06 indicates similarity with laminar flow,
and hence that the assumption of a constant, real-valued eddy viscosity
(vi = 0.5 @ zi) provides a reasonable model in this range. See also page 40.
0)
NOs 1073 102 0"! 10
Figure 1.3.5: Best fit p-values (Equation 1.3.11) from experiments with different
roughness types. Legend as for Figure 1.3.4.
Outer Layer
d
Overlap Layer
r
2
|nner Layer 0-185 K Gyr Vi
Figure 1.3.6: The eddy viscosity distribution applied by Kajiura (1968). The overlap
- layer, which is analogous to the logarithmic part of a steady boundary layer, was only
expected to exist for almost smooth beds. Kajiura suggested that it would probably
disappear completely for r/A > 1/30. The friction velocity is defined by
te = Vfw/2 Aw, and the thickness d of the wall layer is given by d = 0.05 ix /0.
in. rs)
millimetres
Elevation
Z
The comparison in Figure 1.3.7 shows that all three models predict the
magnitude of the velocity defect, and hence In |DI quite well, but only Myrhaug’s
model performs well with respect to the phase (Arg D). The performance generally
improves with increasing model complexity.
The fact that Myrhaug’s model fares best among the ones tested in Figure
1.3.7 corresponds to the fact that his assumed eddy viscosity distribution
corresponds rather closely to the empirical values shown in Figure 1.2.16 for the
relevant Reynolds number. Sleath’s Test 4 has a similar Reynolds Number to
Jensen’s test 12. Myrhaug’s vydistribution may be less suited to model some of
the other cases shown in Figure 1.2.16.
It is not surprising that the agreement between these theories and
measurements is fairly unimpressive. They all assume v; to be a real-valued
function of z only, while most data show a general tendency for the local shear
Stresses to be out of phase with the local velocity gradients, see Figure 1.2.11.
These phase shifts make it necessary for v; to be either a complex function of z
or to be strongly time-dependent, see Sections 1.2.8 and 1.2.9.
40
- £n 1D,1
fe)
ELEVATION
IN
MILLIMETRES
jl
0.1 0.2 0.4 1.0 20 40 10.0
Figure 1.3.8: Comparison of the Bakker & van Doom mixing length model with data
from van Doom (1982), Test M10 RAL. We see that the model predicts fairly
reasonable values of InID,| while the predicted values of Arg Dj are clearly
different and further removed from the experimental values. The identity between
Arg D; and In|DjI is very clearly shown for these conditions, (T, A, r) = (2s, 0.33m,
0.021m).
The only case where it seems well justified to assume that v; is real-valued
and independent of time is for fairly rough beds, r/A > 0.06 where vy is
independent of z as well and given by
see Section 1.2.8 and Figure 1.2.16. It will be noted that this formula is not in
general agreement with Equation (1.2.44) which was suggested by Sleath (1991).
The reason is that the two formulae are not really predicting the same thing.
Sleath’s formula is aimed at the outer layer only, where v;= constant even for
very low values of r/A (see Figure 1.2.16). Equation (1.3.14) is meant to predict a
global value of v; for fairly rough conditions (7/A > 0.06). The two formulae give
identical results for r/A = 0.4.
1.4.1 Introduction
Wave motion is often thought of as purely oscillatory. However,
measurements of the velocity field under real waves show the existence of time-
averaged velocity components almost everywhere.
The magnitude of these steady flow components is generally much smaller
than that of the oscillatory components. However, because their effect is
cumulative, their contribution to the net sediment transport may well be
significant.
Ustokes = (1.4.1)
stokes 0
Figure 1.4.1: Eulerian (left) and Lagrangian (right) net flow velocities resulting from
superimposing a uniform drift on a sine wave in order to obtain zero flow of water.
These velocities are based solely on continuity considerations. That is, the forces
required to set up such a flow pattern in a real fluid are not considered.
onto the sine wave velocity field. This uniform, seaward drift velocity is generally
referred to as the Stokes drift. The considerations behind it are based on continuity
only, ie nothing is said about the forces which would be needed to drive it with a
real (viscous) fluid.
2
uy = oO cosh 2kz (1.4.2)
ninety degrees out of phase, as they would be in a perfectly inviscid wave motion. _
This gives rise to finite time-averaged stress terms of the form -—p uw
which are analogous to the familiar Reynolds stresses —p u’w’.
These stress terms grow from zero at the bed (where w = 0), towards an
asymptotic value of
Figure 1.4.2: Distribution of —p i ina wave boundary layer with constant eddy
viscosity. By comparing the expression (1.4.4) with the expression t= pVav Aw
for the peak bed shear stress due to a purely oscillatory flow, we see that the
= kA
asymptotic value of —p uw amounts to we
55
AND SEDIMENT TRANSPORT
Chapter 1: Bottom boundary layer flow
Iey tir
oe pe ene)
uw’) = ay DRC
+ =
a: uw (1.4.6)
ob
Mans op
uw =e
= ae a A5Ho) (1.4.7)
or
Vt ou
Be _1lo
p ax =~
z+uwt _1-
p T(0) (1.4.8)
a aft
* Vt
Ps
ox
iw +p4H) \ae (1.4.9)
ve
22
Up = oe [3+ et (—4cos§ + 2sinE + ve 2Esin§ + 2Ecos&)]
(1.4.10)
which is based on the assumptions of ge= 0 and 70) = -—p (iWon.
We see that according to this solution, the Eulerian mean velocity tends
towards an asymptotic value of magnitude 3(Aw)*/4c for z>> ds. The
corresponding mass transport velocity, i e the asymptotic time-averaged,
(0) = tw-pgD dy
ai
_ Sx = -P(UW)oo .
arte
If that is changed, the velocity distribution changes accordingly. See Figure 1.4.3.
The behaviour of —p uw and hence of uz, in turbulent boundary layers in
general, is at present neither well documented nor well understood. However, the
results in Section 1.3.3 concerning the structure of turbulent oscillatory boundary
layers indicate that the laminar solution (with v replaced by v;) may have fairly
wide applicability. It was found in Section 1.3.3 that turbulent oscillatory boundary
layers in the roughness range 0.06 < r/A < 0.5 seem to correspond to a constant
eddy viscosity with V2v;/® = 21 = 0.09 VrA. Hence, the stress distribution
(1.4.5) with 8 = 0.09 VrA may provide reasonable estimates for these flows.
ae
Figure 1.4.3: The bed shear stress must in general balance the sum of wind stress on the surface,
the radiation stress gradient and the pressure force due to mean surface slope.
57
AND SEDIMENT TRANSPORT
Chapter 1: Bottom boundary layer flow
Figure 1.4.4a: The wave maker is tured on and sand is entrained under the surface
nodes where the near-bed velocities are strongest.
Figure 1.4.4b; The suspended sediment is convected by the boundary layer drift
towards the surface antinode.
Figure 1.4.4c: Some sand is rising towards the surface under the antinode.
Figure 1.4.4d: Eventually, the sand which settles under the antinodes forms bars.
1.4.5 Undertow |
The shoreward boundary layer drift indicated by the solution (1.4.10) is
usually not observed in the surf zone. In the surf zone, the u-picture tends to be
dominated by the so-called undertow, a seaward mean velocity between the bed
and the wave trough level, see Figure 1.4.5.
The undertow is a gravity driven current related to the phenomenon of wave
setup as illustrated very lucidly by the experiments of Longuet-Higgins (1983). It
occurs because the radiation stress gradient dSx,/dx is not uniform over the depth
under breaking waves while the opposing pressure gradient dp/dx = pg dn/dx
from the wave setup is (nearly) uniform, and therefore dominates near the bed.
Quantitatively, the undertow is described by the time averaged equation of
motion (1.1.21) which for the two dimensional case can be written
undertow
Figure 1.4.5: The undertow occurs in the surf zone because the radiation stress
gradient is concentrated near the surface, while the opposing pressure gradient due to
the setup slope is essentially uniform over the depth.
2 0=
+
ay
earns
+ —.
5, uw (1.4.11)
a.
rm Vel dy
Ao _Si Ee:
MAM, + Io,
ay + IF
5H
— AG
+ ay Ae
+ 5, Hw (1.4.12)
ou
val _= oa
oles oO « Be
(Be 09 +, 972 + OEuo)
au * de
(1.4.13)
u
where the boundary condition vc oe | = T(o) has been used.
Oz z=0
The bed shear stress (0) may include a wind stress as well as the radiation
dS
stress contribution — aa and the setup contribution —pgD A See Figure 1.4.3.
To obtain the undertow distribution, Equation (1.4.13) is integrated once
more with the use of the non-slip boundary condition u(o) = 0. This leads to
Z Zz =
(1.4.14)
From this expression it is obvious that in order to model the undertow, one
must consider the complete system of turbulence from wave breaking (to get the
eddy viscosity Vci), the wave setup 1, and the process of wave transformation and
decay across the surf zone. For examples of recent models see Svendsen et al
(1987), Roelvink & Stive (1989) and Deigaard et al (1991).
61
AND SEDIMENT TRANSPORT
Chapter 1: Bottom boundary layer flow
1.5.1 Introduction
The hydraulic conditions which cause problems of siltation or erosion on the
coast are nearly always mixed in the sense that the velocity field has both a steady
component u and a periodic component u which may have very different
relative magnitudes and different directions.
We are, in the present context, mainly interested in the sediment transport
which results from these combined flows and hence we shall concentrate on the
flow structure near the bed where most of the sediment transport occurs. The
changes in wave height and direction which occur when waves travel through a
variable current field are considered outside the scope of this text. For a recent
review of these phenomena see Jonsson (1990).
ELEVATION
MM
IN
COMBINED
ALONE
ALONE
O
O 1@) 20 30
VELOCITY (cm/s)
Figure 1.5.1: Profile of a pure current u (+), the velocity amplitude lid in a pure
oscillatory flow (x), and both u and luli (e) from the combination of the two
flows. Measurements in an oscillating water tunnel by van Doorn (1982).
e RIPPLE TROUGH
+ RIPPLE CREST
ELEVATION
MILLIMETRES
IN
ELEMENTS
5 Le) I5 20 25 30
VELOCITY, cms!
Figure 1.5.2: Measurements of the steady component u_ and the amplitude U 1 of
the primary oscillatory component from a combined flow (following current) in a
wave flume. Data from van Doom (1981) Test V20 RA+RB.
There is the wave-dominated layer near the bed. Then there is a logarithmic
layer and finally an upper layer which covers sixty to eighty percent of the water
depth.
Since most of the sediment transport normally occurs in the two lower layers
ie, the wave-dominated layer and the logarithmic layer, the following sections are
focussed on those.
ou op at
essere? et
which becomes
ot
5 _gee
Op (1.5.2)
op
ay PS on
> (1.5.3)
then, it follows from Equation (1.5.2), that the shear stress gradient is also constant
over the depth, and we have: In a steady flow with hydrostatic pressure
distribution, the shear stress distribution must be linear. See Figure 1.5.3.
BE
a v=eh reno
P&SOM
>, (1.5.4)
The bed shear stress t(0) defines the typical turbulent velocity at the bed
ux = Vt(0)/p (1.5.5)
a = = da (1.5.6)
Ux K Zo
where X is called von Karman’s constant and has a value of approximately 0.4,
see Figure 1.5.4.
; th z
| ~ factor e
This logarithmic velocity law which is also called the “law of the wall” is not
a good approximation to the velocity distribution near the water surface, see for
example Figure 1.5.2, but it does provide a good and workable model where most
of the sediment transport, in steady flows, occurs.
The logarithmic velocity profile given by Equation (1.5.6) goes to zero at a
finite elevation z) above the bed. This is of course not physically realistic, but it
reflects the fact that Equation (1.5.6) is a horizontally uniform description, which
does not attempt to model the three dimensional flow details close to the roughness
elements. See Figure 1.1.5.
If the bed is perfectly smooth, a thin bottom layer exists within which the
flow is laminar. It is called the /aminar sublayer, and its thickness is generally
taken to be
1
Z = 30 d (1.5.8)
when the flow is fully turbulent, and this has been observed to be the case for
rux/V > 70. Based on this observation we define the equivalent Nikuradse
roughness of other surface geometries by
r = 3020 (1.5.9)
where Zo is determined experimentally by fitting a logarithmic curve of the form
(1.5.6) to measured velocities.
The choice of origin for the z-axis: the theoretical bed level, is not obvious,
but it is generally defined as the level of the origin (z=0) which leads to the best fit
by Equation (1.5.6) to measured current profiles.
The equivalent roughness of many different surface types has been
investigated by Schlichting (1979), and the corresponding, theoretical bed level has
been discussed in detail by Jackson (1981).
iii ge (1.5.10)
for shear stresses in laminar flow, we can define a turbulent eddy viscosity by
tT = pv; (1.5.11)
Then, from the logarithmic velocity distribution (1.5.6) and the
corresponding shear stress distribution
This is a simple and useful formula, but since it is based on the law of the
wall, Equation (1.5.6), it is not reliable near the free surface.
If the layer of interest is very thin compared to the total flow depth, we often
apply the so-called constant stress assumption. That is, the last terms in the
expressions for shear stress, Equation (1.5.12) and for eddy viscosity, Equation
(1.5.13), are left out. Thus, in the constant stress layer we have 7 = T(o) and
Vt = K Ux Z.
The velocity distribution in the constant stress layer is still logarithmic and
given by Equation (1.5.6), and the eddy viscosity corresponds to the straight line in
Figure 1.5.5.
Wile ce bY (1.2.34)
du du
dz dz
where the term —uw_ vanishes if the flow is horizontally uniform, because then
w =0 by continuity.
For oscillatory flows the concept of eddy viscosity is less trivial as discussed
in Section 1.2.8. However, an analogous expression can be derived for a turbulent,
purely oscillatory flow (u = 0)
Based on the expressions (1.1.22) and (1.1.25), for the steady and the
periodic components respectively of the total momentum transfer, we obtain the
following eddy viscosities
Vv w =
Wp
a
_ -uw- iw-iw-ww’
= ai= = Ley (1.5.15)
Oz Oz
where Vc is the eddy viscosity felt by the steady flow component u(z) and Vy is
the one felt by the periodic component 1(z,t). The interesting thing about these
two expressions is that they are not similar at all. Judging from the appearance of
these expressions, it would be quite surprising if vc and Vy turned out to be
identical or just similar.
From (1.5.14) it is clear that vc must be constant in time, but vw may well
be a function of time as well as of z. As for the pure wave case, which was
discussed in Sections 1.2.8 and 1.2.9, it is quite possible for the expression (1.5.15)
to take both negative and infinite values.
Coffey & Nielsen (1984, 1986) showed examples of corresponding values of
Vc and Vw from wave flume data (van Doorn 1981 Tests V10 and V20) and from
tunnel data (van Doorn 1982, Test S10). They found that in both cases vc was
typically three to four times greater than vy .
A similar data set derived from the velocity data of van Doorn (1982) Test
M10 is shown in Figure 1.5.6. The plotted values of vy are calculated as lv|I
from the first harmonic of 4(z,f) in accordance with Equation (1.2.42) from
Section 1.2.9.
The plotted values of Vc are derived from
45 (z)/p Ux (1 z/D)
Ve = Se ( lJs
3) 16 )
du du
Elevation,
cm
1 2 3 4 5 6 tj 8 9
Eddy viscosity, cm
Figure 1.5.6: The magnitudes of the eddy viscosities vc (0) and Vw (+), based on
Equations (1.5.14) and (1.5.15), may well be very different. Data from oscillating
water tunnel, van Doom (1982) Tests M20, (A, T, r, <u>) = (0.33m, 2s, 0.0021m,
0.2m/s). The line corresponds to Vc = 0.4 ux z, see Figure 1.5.5.
Despite the scatter in the values of Vc it is quite clear that the magnitudes
of Ve and vw are quite different. Typically, vc is three to four times greater than
Vw.
While this conclusion is perfectly in line with those previously drawn by
Coffey & Nielsen (1984, 1986), it is acknowledged, that the determination of ux
through a log-curve fit introduces an element of arbitrariness when the extent of
the "logarithmic layer" is poorly defined.
The definitions (1.5.14) and (1.5.15), of Vc and vw are the ones applied,
with more or less explicit acknowledgement, in previous eddy viscosity based
models of wave current boundary layer interaction. There are however, other
options which may turn out to be more appropriate. i
It may for example be reasonable to give the dominating term — uw of the
total momentum transfer expression (1.1.22) explicit consideration in the equation
of motion for u(z) instead of hiding it under the collective label of T(z). This
approach is outlined in Section 1.5.9 page 91, and in Section 7.2 page 293.
With respect to the eddy viscosity concept, the advantage of this approach is
that the expression for the current eddy viscosity vc then becomes analogous to
that from steady turbulent flow as the term — Zw is removed. Compare the
formulae (1.2.34) and (1.5.14). The only difference will be due to the difference in
turbulence intensity.
With the definition (1.5.14) however, vc is a function of the angle between
the current and the direction of wave propagation. The term —iw , which is
positive or zero (see Figure 1.4.2) makes a positive contribution to vc for a
following current and a negative contribution for an opposing current.
Sleath (1987) found, for purely oscillatoryflow in an oscillating water tunnel,
thatuw was typically ten times greater than u’w’. It is therefore natural to expect
— uw to be the dominant term in ve defined by Equation (1.5.14). The magnitude
of —uw will be greater under real progressive waves than in a tunnel where the
contribution described by Longuet-Higgins Equation (1.4.5) is absent. The
corresponding difference between u-profiles in wave flumes and in tunnels can be
seen by comparing Figures 1.5.1 and 1.5.2. For more details, see Section 7.2.
millimetres
in
Elevation
Figure 1.5.7: Root mean square vertical velocity fluctuations in five cases; pure
current 0.Jm/s (*) and 0.2m/s (0) depth average, a pure wave motion (A, T) =
(0.10m, 2.0s) (x) and the waves superimposed on each of the currents (®) and (0).
Data from oscillating water tunnel, van Doorn (1982).
For all cases with waves, w’rms peaks at the top of the roughness elements
(z = 2mm ) with a maximum value roughly equal to 0.5 Ux .
For z < zj there is no apparent increase in turbulence intensity due to the
addition of currents, but for z > zy the data do not contradict the simple
and remembering that Sleath (1987) found, for purely oscillatory flow, that the
term — iw totally over-shadowed the turbulence term -uv’w’. If the turbulent
73
AND SEDIMENT TRANSPORT
Chapter 1: Bottom boundary layer flow
© NO CURRENT
x U,/U,*.51
2 e =.
ROUGHNESS
ELEMENTS -£n 1D)
1 0.2 04 O6 08 10 20 40 60 80100
Figure 1.5.8: Measured values of the dimensionless velocity defect Dj(z) for tests
VOORA, VIORA, and V20RA from van Doom (1981). Even for superimposed
currents as strong as iix/ttx = 0.79 there is no evidence of change to the wave
boundary layer structure.
It must of course, be expected that with very strong currents, such as some rip
currents, the wave boundary layer structure can no longer be assumed unchanged.
However, the presently available data show no significant effect for the parameter
range 0 < Ux/ix < 0.79, see Figure 1.5.8.
The data in Figure 1.5.8 correspond to colinear waves and currents, but
recent measurements of the perpendicular case by Sleath (1990) also show little or
no change to the wave boundary layer structure due to superimposed currents of
moderate strengths, ux/A @< 0.08.
75
AND SEDIMENT TRANSPORT
Chapter 1: Bottom boundary layer flow
er (1.5.18)
where Ef is the wave energy flux per unit length of wave crest and Dz is rate of
wave energy dissipation per unit bed area, see Section 1.2.6. The energy flux is -
according to linear wave theory - given by
2kD
Op anB 4n
le anit Dear sinh 2D?
(1.5.20)
100
0.2
0.1
1 eZ 2 S 10 20 S16)
Figure 1.5.9: Wave energy dissipation factors derived, through Equation (1.5.29),
from the wave height decays measured by Simons et al (1988). The symbols *, 0, 0,
A correspond respectively to the following relative current strengths (ux/ttx S70)
1.3-2.1, 4.1-9.8, 5.5-10.5. The curve corresponds to Swart’s formula, Equation
(E222);
applied in the frame of reference which travels with speed U. That is, Ty is the
wave period in this moving frame of reference and the wave number must be
determined from the dispersion relation
2
kDtanhkD = —>D = alt = | (1.5.21)
where Tq is the wave period seen from the laboratory frame of reference.
For the wave energy dissipation we apply Jonsson’s (1966) definition
2
De = 3 pfe(A@g)° (1.5.22)
so we have
Gala 2
me Epg H* (Cor + »)hee 0 fe (A@a)’ (1.5.23)
where Cygr+U is a constant when the depth is fixed and the near bed velocity
amplitude is given by
7H
A@q = Tq sinh kD (1.5.24)
H
Hx) = ; i (1.5.26)
1
Lot Hy igreTS; ao)
3gTa (Cer + U) sinh” kD
In most previous studies, the observed wave height variation has been fitted
by an exponential which corresponds to viscous or laminar energy dissipation, i e
and the measurements have been discussed in terms of the dissipation factor
a= -—— (1.5.28)
For the real wave height variation under turbulent conditions, which is given
by Equation (1.5.26) there is no constant corresponding to @ but for small
energy dissipation rates we have
8"
Wo a 3ee fe (1.5.29)
3 g Ta(Cgr
+U) sinh® kD
This shows how, for fixed fe, the observed o-values will decrease for
increasing U and vice versa. This effect of the mean current on dH/dx was also
noted by Simons et al (1988).
The fe-values in Figure 1.5.9 which are based on measured o-values and on
the analysis above show that fe and presumably the general structure of the wave
boundary layer are practically unchanged by the superposition of currents. This
agrees with the friction factor measurements of Sleath (1990) and with the findings
related to Figure 1.5.8 about the wave velocity defect function.
Hence, the data indicate that, for practical purposes, the wave boundary layer
structure and the wave energy dissipation factor for a combined flow can be
calculated as if the current was not there. That is of course, after the near bed
velocity amplitude Aq@ _ has been calculated with due respect to the current by
using the dispersion relation (1.5.21). The conclusion above should not be applied
uncritically to conditions outside the range of conditions represented by the data in
Figures 1.5.8 and 1.5.9.
a) Assuming that the average bed shear stress and hence ux is kept
unchanged, what is the change to the current profile due to the
superposition of waves? See Figure 1.5.10, left-hand side.
b) What are the changes to the current profile when waves are added in
such a way that the reference current velocity uy = u(zr) is kept
unchanged? See Figure 1.5.10, right.
Inz Inz
The first approach which was taken by Lundgren (1972) is the most rational
because it leads to logically straight model building, while the latter approach leads
to iterative thinking, as applied in the model of Grant & Madsen (1979).
In the following we shall follow Lundgren’s line of thought where u« is
taken to be known and fixed, - for example determined by a fixed mean surface
slope
Wave induced
mixing
Ve
l= tx /@ (1.5.32)
However, the relationship between / and the various boundary layer thickness
values (Section 1.2.7) may depend on the relative roughness r/A , the relative
current strength ux/ ux and possibly on the angle @ between the current and the
direction of wave propagation
tt = Fi(y>
I Ux
,Q). (1.5.33)
Zo A@
The shape of the current profile for z</ seems to be as much like a linear
function of z as any other simple shape including the logarithmic shape suggested
by Grant & Madsen (1979), Fredsoe (1984), and by Christoffersen & Jonsson
(1985), at least for moderate values of r/A, see Figure 1.5.12.
Outside this layer, the wave-generated mixing is relatively weak, so the law
of the wall may be assumed to apply in which case the current distribution must be
logarithmic and can be written as
= Ux vE
u(z) = a Lape forez¢> (1.5.34)
AuOMT = —In—
hg
K Zo
(1.5.35)
can be expected to depend on the relative current strength, the relative roughness,
and on the angle between the current and the direction of wave propagation
4 ~ F(—,-,9) (1.5.36)
oO *
81
AND SEDIMENT TRANSPORT
Chapter 1: Bottom boundary layer flow
1000
100
10
1.0
0-16 5 10 15 20
Figure 1.5.12; Even for fairly small relative roughness (r/A), the current profile in the
inner, wave-dominated layer is not logarithmic. Data from Sleath (1990), A =
0.141m, T = 1.76s, r= 1.49mm, ux = 1.9cml/s (x) or 2.3cm/s (0).
z
—Zo = 140.06 ig
ifs (1.5.37)
Ux
100
40
20
10
However, Sleath (1990) found that this formula seemed inadequate for data
with smaller relative roughness ( r/A < 0.J ). Subsequently Sleath (1991)
developed a model which leads to the expression
ET: fe
*
ue (1.5.38)
Zo
which is shown together with the presently available data in Figure 1.5.13.
Several field studies, e g Cacchione & Drake (1982), Grant et al (1983),
Coffey (1987), Lambrakos et al (1988) and Slaattelid et al (1990), provide useful
83
AND SEDIMENT TRANSPORT
Chapter 1: Bottom boundary layer flow
Ve = 4Vw (1.5.39)
Figure 1.5.14: Some of the eddy viscosity distributions assumed by previous authors
for combined wave-current flow (shown qualitatively). G&M: Grant & Madsen
(1979), C&J: Christoffersen & Jonsson (1985), M&S: Myrhaug & Slaattelid (1989),
Sleath (1991). In all of these studies, it was assumed that the waves felt the same
eddy viscosity as the current.
strength Ux/ its up to 0.79. In addition, the wave energy dissipation data in
Figure 1.5.9, page 76, indicate no change to the wave friction factor for much
higher values of the relative current strength uU+/ be .
The essence of the problem is then to specify the distribution of the eddy
viscosity Vc which governs the distribution of u(z) .
Following the line of thought of Lundgren (1972) we assume that an outer
layer exists, where the eddy viscosity felt by the current u(z) is given by
just as in a pure current boundary layer which is thin compared with the flow depth
(a constant stress layer), see page 68.
Ku, 2@ U,/K
Figure 1.5.15, left: The assumed distribution of the eddy viscosity vc which is felt
by the steady component u(z) of a combined wave-current boundary layer flow.
Right: The resulting velocity distribution.
u(z) =
||
x
4 for z</ (1.5.42)
in the inner layer.
In the outer layer the current distribution is logarithmic and may be written
& a 4
u(z) = PsIn D5 On *2>"] (1.5.43)
ie a ee (1.5.44)
At the present state of the art, it does not seem possible to derive a general
expression for / which is valid for all flow conditions (see the discussion related
to Equation 1.5.36). As more experimental information becomes available, for
example in the form of Equation (1.5.36), it will be straightforward to make use of
it using the framework of Figure 1.5.15 with ] and zg determined from
Oa he gD ae)
unknown. However, for fairly large relative roughness, 0.06 < r/A < 0.5, where
the structure of the wave boundary layer corresponds to a constant, real-valued
eddy viscosity (see Section 1.2.9, page 40), a simple estimate of F can be
obtained if the influence of @ is neglected.
For such wave boundary layers, the eddy viscosity can be estimated by
Then, since the current eddy viscosity in the lower layer is written as
Vc = Kul (see Figure 1.5.15) , this gives the following expression for /
1 = 0.016 Aw , (1.5.47)
K Ux
100
40
20
10
1 2 4 10 20 40 100
Figure 1.5.16: Apparent roughness increase as function of the relative current
strength for the same data which was shown in Figure 13.13, The straight line
corresponds to Equation (1.5.48). <
Ux Zr Ux Zr
u(zr) = nae eee ore (1.5.49)
0 Ux > A’
where the unknown is the time-averaged friction velocity ux. The form of the
function F is assumed known, for example, given by Sleath’s (1991) expression
(1.5.38) or by Equation (1.5.48).
For the purpose of iterative solution, it will generally be convenient to rewrite
this Equation (1.5.49) in the form
(1.5.50)
89
AND SEDIMENT TRANSPORT
Chapter 1: Bottom boundary layer flow
When the current friction velocity u* has been found from the iteration, the
value of F can be calculated and so can Zqg = Zo F and / =e'z )F. All
components are then available for drawing the velocity distribution
mk for z<l
Kk /
u(Zz) = (1.5.51)
Ux Z
in for z>/
K Za
Example 1.5.2: Current profile from wave data and reference current
To illustrate the use of the model above for deriving details of the current
distribution u(z) from the current velocity at a single point plus the wave data and
the bed roughness, consider the numerical example [A,7,r,u(1.0m)] =
[0.9m, 85, 6cm, 0.45m/s].
Inserting these data into Equation (1.5.50) with F given by Equation
(1.5.48) gives
ea wistae
ie Sea.
le finn”
“ngEs' teeP
pet Ou 0 tire (M
: a
0.06/3
which converges rapidly to the friction velocity value u« = 0.0425m/s . This
corresponds to F = 0.44 “2=7.32, and hence za= zo F =0.0146m, and | = zq
2
SESE Ot
EOSONG 90 S ee ee Co
Mast 0.18
Ux =
1.0 0.9: =2n/8 [0.9
0.06 ) 6.22 -In(1+
0.52
= )
In 0.06/30 ~ In (1 +0.19
(1.5.53)
which equally rapidly converges to a friction velocity of ux = 0.048m/s. With this
value, the model of Sleath (1991) gives the velocity distribution which is shown by
the broken line in Figure 1.5.17.
Figure 1.5.17: Velocity profiles obtained from the data set [A,7,r, u(1.0m)] =
[0.9m, 8s, 6cm, 0.45m/s], with two different wave-current boundary layer models.
The full line corresponds to the model defined in Figure 1.5.15 with F = za/zo = 0.44
A@/ux , while the broken line corresponds to the model of Sleath (1991).
indication by the data that opposing currents (@= 180°) experience a larger,
apparent roughness increase Za/Zo than following currents, or in other words,
written
—wWww’
Val == V + aa
— (1.5.55)
ES:
az
(which is different from the definition used in the simplistic discussion above, see
Equation (1.5.14)), and integrated once to yield
ou _map lop,
Velcazs fig ai =,
ot 1)
nw’ (1.5.56)
For a hydrostatic mean pressure gradient due to the mean surface slope a
where the mean bed shear stress equals —pgD a. Equation (1.5.56) can be
rewritten as
du
Vala = ¢ (-D) 1
2 +iw = = -gp0-5a Tigo
+ UW (1.5.57)
which, for a thin bottom boundary layer where z<<D, can be conveniently
approximated by
ln. = eS & 9D a
The current distribution is then obtained by integration and use of the
non-slip boundary condition u(o) = 0
7 luslite + iw
Z ———
10
ZO
ea 2 4 6 8 10 12
Figure 1.5.18: Current distributions for a pure current and for a current with
following, opposing and perpendicular waves superimposed according to Equation
(1.5.59) with the eddy viscosity distribution (1.5.60). The relative "drift strength"
| (GW) oo/un | is, for both following and opposing waves, equal to 0.5. Legend: x:
pure current, *: following waves , +: perpendicular waves or tunnel, rectangle :
Opposing waves.
since the wave friction factor for constant eddy viscosity amounts to
2/V\Re = 2Vv/(A’@) = V2 8s/A, see Equations (1.2.16) and (1.4.4).
It seems rather natural to incorporate uw into new models of wave-current
boundary layer interaction, along the lines suggested above. For further incentive
see also Section 7.2, page 293. However, some questions need to be addressed
regarding the general nature of “w in turbulent boundary layers over rough beds.
Longuet-Higgins’ expression (1.4.5) was derived for laminar flow over a
plane bed and the only vertical velocities involved are due to the thickening and
thinning of the boundary layer under a progressive wave. Over a rough bed
however, there will be additional periodic, vertical velocities caused by the bed
geometry, and the significance of their contribution to “uw must be considered.
2.1.1 Introduction
For the purpose of sediment transport modelling, it is necessary to consider
three types of forces which govern the behaviour of cohesionless sediment
particles whether they are resting at the bed or moving around in a slurry or a thin
suspension. These are: the gravity force Fg = Mg ; intergranular forces related to
collisions or continuous contact; and the fluid forces which may be due to surface
drag or fluid pressures.
Hanes & Inman(1985b) suggested a typical value of 3] degrees for beach sand.
When a horizontal sand bed is exposed to a fast, steady flow, a finite top
layer of sand will start to move with the flow, partly as bed-load and partly in
suspension. The fact that the moving layer is of finite thickness is significant
although seemingly trivial, because it shows that the moving sand has increased
the strength of the sand below.
Since the shear stress is not decreasing downward, the top layer of immobile
sand is able to withstand the shear stress which eroded the top layers when the
flow was started. This is due to the fact that the moving sand is transferring at least
part of its weight to the bed as effective stresses and thereby increasing the
effective normal stress in the bed. The effective nomal stress transferred by the
moving sand is generally referred to as the dispersive stress.
Bagnold (1954,1956) studied the normal and tangential stresses in granular
flows and suggested that they be given as functions of the shear rate du/dz and of
the linear sediment concentration 2 which is the relative surface proximity
between sediment particles, see Figure 2.1.2. The linear sediment concentration is
related to the volumetric concentration c by
1
= oo 2 ii
(CmaxZ0, eed Ge)
where Cmax is the maximum concentration corresponding to grain contact. The
linear concentration 4 increases drastically as c approaches cmax, see Figure
aah
k—
d —
which has been termed the Bagnold number. A purely inertial regime is found for
B > 450, and a purely macro-viscous regime corresponds to B < 40.
In the inertial regime the stresses are proportional to the linear sediment
concentration and to the shear rate squared. Bagnold obtained the following semi-
empirical formula
rn
Syren Oe =idee
55 ps(Ad )
Fis for B > 450 (2.1.3)
T du
Oe ae Oe = 13(1+A)(1+A/2)
pv i for..8 < 40.2.1)
with a somewhat greater stress ratio than in the inertial regime, namely
tan Qd = T/6e = 0.75.
With typical values (s, d, 4, du/dz, v) = (2.65, 2:10%m, 1, 100 s!,10° m/s),
corresponding to B = 10.6, sheet flow conditions under waves will generally be in
the macro viscous regime.
Fp = ve
Er Mek
ax (2.1.5)
a dz paV
Pg
dp
Ox
Figure 2.1.3; The total pressure force is, by Green’s theorem, equal to the volume
times the pressure gradient when the pressure gradient can be considered constant
over the body.
If the shear stresses in the fluid are relatively small (and in particular for
inviscid fluids), the equation of motion (1.1.10) may be used to write the horizontal
pressure force in terms of the horizontal fluid acceleration
ye
Va _ pV du
oa (2.1.6)
When calculating the pressure force Fp ona body which is held fixed while
the fluid is accelerating past it, an extra mass oCywV_- must be added
corresponding to the volume of surrounding fluid which the body keeps from
accelerating, so for a fixed body in a horizontally accelerated fluid we get
Correspondingly, the force required to give a sediment particle with volume V and
du ' He
relative density s the acceleration ee through aresting fluid is
d
F = p(s+Cw) Vi (2.1.8)
where p(s+Cm) V_ is called the virtual mass of the body, and pCmV_ is called
the added hydrodynamic mass. The added mass coefficient for a sphere is 0.5, for
a long cylinder itis 1.0, see e g Lamb (1936).
For a particle which is fixed in a wave motion with the homogeneous velocity
field u=A@sin wt Equation (2.1.7) gives
Fp = p(1+Cmu)
VA cos wt (2.1.9)
Apart from the pressure forces described above, which can be evaluated on
the basis of inviscid flow theory, a particle exposed to a viscous or turbulent flow
will in addition feel drag forces. Drag forces occur in two varieties: skin friction
and form drag. Skin friction contributes most of the drag on slender, streamlined
bodies like kajaks, while form drag dominates for plump shapes like spheres and
most sediment grains. The drag force is normally given on the form
ti pAColulu (2.1.10)
where A (~1d’/4) is the cross sectional area facing the flow, and Cp is the
drag coefficient which depends on the sediment shape and on the Reynolds
number, dlul/v. When the flow becomes laminar the drag force is actually
proportional to the flow velocity and, for small (d lul/v < 1) spherical particles, the
drag force is given by Stokes’ law
Fp = 3npvdu Qainty
24
Co = dlul/v
(2.1.12)
Drag coefficients for spherical particles are plotted in Figure 4.2.1, page 164.
By comparing the expressions (2.1.9) and (2.1.10) for the pressure force and
the drag force respectively on a spherical particle in a wave motion, we see that the
force amplitude ratio
we nx (te
FF ) (2.1.15)
1 op
( = —K| .P8& % (2.1.16)
w
—FP1 op 44
pg dz
The water velocity u = (u, w) is in this definition the equivalent clear water
velocity, i e, the flow rate per unit area. The actual velocities of the water in the
pores are somewhat greater.
For the situation in Figure 2.1.4, where the seepage rate corresponds to the
vertical, equivalent clear water velocity w, the total vertical pressure gradient in the
pore water must be
OD. yiSevens
retin dbaWw.
pondiy (2.1.17)
2.2.1 Introduction
The initiation of sediment motion under steady flows and under waves has
attracted considerable interest in the past because it is a philosophically appealing
concept. In practical terms however, it is a very difficult concept to deal with.
Firstly, because “initiation of motion” is difficult to define - is it when one in a
thousand grains moves or when one in a hundred moves? Secondly, because the
complicating variables in a natural situation is very large. For example, the sand
bed is never left perfectly smooth from previous events. Relict bed forms will be
present and initiation of motion will occur sooner near the crest of these bedforms
due to local enhancement of the bed shear stresses.
In addition, biological activity will also complicate the micro-topography and
excretions from animals may tend to glue the sand particles together. Nevertheless,
we shall consider a few classical approaches to the description of incipient
sediment motion in the following.
_ Ao)
velocity %* are quantities which are easily measured, t(0) = p g DJ, where D is
the flow depth and / is the hydraulic gradient.
In connection with wave motion, the Shields parameter (corresponding to
total stress) is generally defined in terms of the peak bed shear stress
q _ WfwAoy _ 1
=usip (ele dol a Dédala wale aa i ow
where fw is the wave friction factor, defined on page 23, and this is the notation
which will be used in the following. |
,
= peDed
Tv’
228
is therefore frequently used to predict initiation of motion and the magnitude of
moving sediment concentrations.
Correspondingly, 1’ is often referred to as the effective stress in connection
with sediment transport.
A comprehensive discussion of the concept of effective stress in steady flows
is given by Engelund & Hansen (1972).
If the bed is flat, the form drag is absent so, from that point of view, 7 = T
and @’= 80.
There is, however, some indication that the total shear stresses on flat sand
beds under waves may not be totally effective with respect to transporting
sediment, see Section 2.4.4, page 121.
5 =
Yoened
fps p (Aw) 175 W (2.2.5)
where the special grain roughness friction factor, f 2.5 is based Swart’s (1974)
formula (1.2.22) and a roughness of 2.5d50
0.194
fos = exp [5.213 ee Sty AMA (2.2.6)
The relationship between 025 and 9 for both rippled and flat sand beds in
oscillatory flows is illustrated by the data of Carstens et al (1969) and of Lofquist
(1986) in Figure 2.2.1.
The data show that for rippled beds 9 is generally an order of magnitude
greater than 625 with no systematic trend between different grain sizes.
For flat beds, @ (= 6’) is also considerably larger than 025 , by a factor five
or so, when the activity level is high (625 2 0.3). For a flat bed at low activity
level, one would expect to find 8 = 0’ = 025.
10
0.1
0.01 0.1 1
Figure 2.2.1: Relationship between the "total Shields parameter" @ and the grain
roughness Shields parameter 62,5 for the presently available oscillatory flow data.
Legend bar : Lofquist 0.55mm sand, + : Lofquist 0.J8mm, * : Carstens et al
0.19mm, rectangle : Carstens et al 0.30mm, X : Carstens et al 0.59mm . All of the
above correspond to rippled beds, while the triangles correspond to flat beds,
Carstens et al 0.19mm and 0.30mm.
time. Wave friction factors and the corresponding hydraulic roughness of sand
beds under waves are discussed in detail in Section 3.6, page 145.
0.60
0.40
0.20
0.10
0.06
0.04
0.02
dv (s-1)gd /4V
0.01
0.4 1 2 4 10% 2820700840 100 200
Figure 2.2.2: The Madsen-Grant diagram for the initiation of sediment motion in
oscillatory flow and on oscillating trays. Data from Manohar (1955) and from
Carstens et al (1969).
@¢ = 0.05 (2.2.7)
for the critical Shields parameter seems justified in most practical cases.
The effective shear stress Te which corresponds to O¢ [Tc = P (S—1) gd Oc]
is called the critical shear stress for the particular sediment.
The Shields criterion ®¢ = @c(u«d/v) is but one of many criteria for the
initiation of sediment motion. Many others have been suggested. For a
comprehensive review, see Hallermeier (1980).
T+?
Glia 0.04 (2.2.8)
2.3.1 Introduction
The following section summarises current experimental facts about steady
bed-load transport, mainly in the light of the theory of Bagnold (1956), and with
emphasis on concepts and formulae which are transferable to oscillatory flows and
combined wave current flows.
Most of the steady bed-load formulae are of the form ® = ®(6’). That is,
they state direct relationships between effective bed shear stresses and
dimensionless transport rates without considering underlying details such as, the
amount of moving sediment, and the typical speed with which that sediment is
moving.
It is necessary, however, to know about these details in order to adapt steady
flow models to unsteady flow situations. Sections 2.3.3 and 2.3.5 therefore attempt
to extract information about the amount of bed-load and about the typical speed
with which the bed-load moves in steady flows.
The bed-load has been defined in different ways depending on the context. In
relation to measurements, it is often defined as that part of the total load which
travels below a certain level or (very pragmatically) as the part which gets caught
in bed-load traps. For modelling purposes however, it is more convenient to apply
the definition of Bagnold (1956).
Bagnold defines the bed-load as that part of the total load which is supported
by intergranular forces. The rest, i e, the suspended load and the wash load are
supported by fluid drag.
Obviously a given grain may well be supported partly by intergranular forces
and partly by fluid drag and hence contribute to both the suspended load and the
bed-load. This makes the bed-load practically unmeasureable in situations where
suspension is present as well, and this is of course of some concern. We shall
however stay with Bagnold’s definition in order to make use of the advantages it
gives with respect to rational discussion and modelling.
onto the top-most layer of the immobile bed; cg(z) is the volumetric concentration
of bed-load.
Assuming then, that the yield criterion for the top layer of immobile grains is
see Figure 2.3.1, we see that the amount of bed-load which is in equilibrium with
tT’ is given by
TT -T%
cp(z) dz
0
Pp (s-1) g tangs Se
Topmost
non-moving
layer
zh pona T=T.+O0,tanyp,
Figure 2.3.1: For the equilibrium bed-load transport rate, the dispersive stress Oe
must satisfy the yield criterion t = T+6e¢ tangs.
Lp = .. Jca(z) dz (2.3.4)
0
Introducing this expression for Lg into Equation (2.3.3) we see that the
vertical distribution scale measured in grain diameters is
Lp 0’ -@-
Te SE ee CRT YEREST OEM
d Cmax tan Qs ( :
Bagnold (1956) gave tangs = 0.63 as a typical value for fairly rounded
grains corresponding to a maximum concentration of the same value,i € Cmax =
0.63 [vol/vol], and he noted that the product cmax tan @s is fairly constant at
about 0.4 for different grain shapes.
Hence, as rules of thumb we have
and
is well understood because neither the sediment velocity distribution wUs(z) nor
the concentration distribution cg(z) through “the bed-load layer” are well
understood. We can, however, still predict Qs empirically with reasonable
confidence for steady flow because it has been measured directly in a large number
of experiments. Most of the data support transport formulae of the form
@, =i Qs
unease @.3:10)
®z = 8(0’-6,)!° (2.3.11)
of Meyer-Peter & Muller (1948), except that for high stress values (@ > 1), the
numerical constant 8 is too small and a value of about 72 seems more
appropriate. The slightly different formula
100
2 SYMBOL MATERIAL d(mm) §S SOURCE i009
Gravel 28.65 268
Sand 5.20 2.68
Lignite Breeze 5.20 1.25 [ Meyer-Peter
Baryta 5.20 4.22
Sand 0.79 2.68 Gilbert
Sand 0.69 2.67 Wilson :
100 ae 100
10
12 (6'-0.05)
Meyer-—Peter
0.1
0.01 0.01
0.01 0.1 1 10
Figure 2.3.2: Total sediment transport rate under steady flows over flat beds. Due to
the fact that the sediment was fairly coarse (0.7mm), suspended load contributed little
to the transport rates for the shown "high stress data" (0 > 0.5).
where the expression (2.3.7), page 112, has been inserted for the amount of
bed-load per unit area Cmax Lg . This reveals that the typical sediment velocity
defined by
1
oo
Oz
CmaxLe Ny~Cmax 25( 66) d poo spy
6"— Oa pag3.1al4
w ~ 750
Up = 4.8 ux (2.3.45)
nla
eecaok es (2.3.16)
That is, for a given friction velocity us, the velocity gradients in the
bed-load layer are inversely proportional to the weight of bed-load. This weight
being proportional to Lg.
From the form of Equation (2.3.16), it can be seen that, if the constitutive
relationship for the granular flow is written in terms of an effective viscosity, Ve
2 dus
T/p = Uk = Ve<—>
- ail
this viscosity must be proportional to ux Lg
Ve ~ U«xLp (2:3.17)
As mentioned above the result (2.3.15) for the typical sediment velocity and
the corresponding result (2.3.16) for the averaged velocity gradient disagree with
Bagnold’s constitutive relationships from Section 2.1.2. Thus, Bagnold’s formula
dus 1
92
epee
dz
which is in conflict with Equation (2.3.16), and with Lg ~ 0’ it leads to
d
Qa ~ CmaxLle Ug = Cmax Lg (La <*>) ~ 9” (2.3.19)
A model of this type was developed and discussed in detail by Hanes &
Bowen (1985).
If Bagnold’s formula (2.1.4) for the hyper-viscous regime is used instead of
(2.3.18) above, the result for the transport rate is
2.4.1 Introduction
The following section describes the processes of bed-load and sheet-flow
sediment transport under waves. Or in other words, the quasi-steady processes of
sediment transport which occur under waves over effectively flat sand beds.
Under sheet-flow conditions, the sediment transport rate is not entirely due to
bed-load transport in Bagnold’s sense, see Section 2.3.2. However, the assumption
of quasi-steadiness may be applied to the total transport rate under certain
conditions. That is, when the thickness of the sediment transporting layer is small
compared to the typical distance woT settled by a suspended sediment particle
during one wave period.
Such conditions can be observed at moderate flow intensities (025 S 1)
over artificially flattened beds before vortex ripples form, and for high flow
intensities where vortex ripples are naturally absent, see Section 3.4, page 135 ff.
The approach taken is essentially to adapt the steady flow sediment transport
models from the previous section to oscillatory flows and to combined
wave-current flows through quasi-steady considerations.
The quasi-steady approach is used because, it seems fairly reasonable to
consider the processes of bed-load and sheet-flow sediment transport as quasi
steady. There are however two main problems with this transfer of technology
from steady flows to oscillatory and combined flows.
Firstly, steady flow transport formulae are generally of the form ® = ®(6’),
i e, they assume knowledge of the effective (skin friction) stress 6’, and our
knowledge about the effective Shields parameter 9’ for waves over sand beds is
sparse. See Figure 2.2.1, page 106, and Section 3.6, page 145 ff.
The most commonly attempted way around this problem is to replace 6’ by
the more "calculable" 625 and we shall see that this approach is quite justifiable in
the sense that it leads to good predictions.
It should be kept in mind, however, that 625 and 6’ are conceptually
different as discussed in Sections 2.2.4 and 2.2.5.
The second main problem with transferring steady flow sediment transport
formulae to coastal conditions, is to deal with those effects of unsteadiness which
cannot be ignored.
This problem is addressed in detail in Section 2.4.4, page 121, and the
conclusion which is reached is in essence as follows.
It may well be possible to model bed-load transport over flat beds under
The remarkably close agreement with the directly adapted steady flow
formula is unexpected because the data in Figure 2.2.1, page 106, indicate that the
total Shields parameter on a flat sand bed is generally much greater than 625. If
the data in Figure 2.2.1 are to be trusted, in the sense that measured energy
dissipation rates are true indications of the effective sediment transporting stress,
this effective stress should be of the order 5625 when @25 = 1. Hence the
expected magnitude of Lgmax , based on Equation (2.3.7) and assuming "perfect
quasi-steadiness", would be Lgmax =~ 2.50’d = 2.5(5625)d = 12.5 5d.
That is however, about five times more than what was observed.
The writer can offer no clear explanation for this at present, but the observed
amounts can of course still be reasonably predicted by Equation (2.4.1), or maybe,
a fine-tuned version of it, using a power of about 3/4 rather than / as suggested
by Sawamoto and Yamashita.
Figure 2.4.1: Measured values of the peak amount of bedload on a flat bed under
"sine waves" in an oscillating water tunnel. The data are from Sawamoto &
Yamashita (1986) and include relative sediment densities in the range /.58<s<2.65,
and grain sizes in the range 0.2mm<d<1.6mm. The curve corresponds to Equation
(2.4.1).
Pr).
10(8.; — 0.05)1.5
Sleath (1978) studied the instantaneous bed-load transport rates under sine
waves as well as ®7% . He found that ®(t) varied approximately as sin*(ot + sp)
with a phase shift @g of the order ten to twenty degrees ahead of the free stream
velocity, ie
Assume further that the instantaneous sediment transport rate is given by the
adapted Meyer-Peter formula
The instantaneous bed-load transport rate will then have essentially the shape of
Sleath’s sin‘-expression (2.4.3), and since sin*t = 38 it leads to the expression
(2.4.2) for the averaged transport rate.
We note that since Equation (2.4.5) includes no phase lag between @7(t) and
@(t), we get sz = Qr. In reality however, it would be reasonable to expect some
time lag of @(t) behind 6’(¢) and hence, in general @g < @r.
-1 0 1 2 3 4 2, 6 ei
Figure 2.4.3: Sawtooth waves with zero mean flow and opposite asymmetry. These
waves would both generate zero net bed-load transport according to formulae of the
type (2.4.7).
While such formulae have the advantage of simplicity, they have the
disadvantage of not being able to model the net sediment transport rate which
results from certain types of wave asymmetry.
That is, asymmetry which are described essentially by the velocity moments
such as that of a 2nd order Stokes wave can be represented, while "sawtooth
asymmetry" is overlooked. Thus, the two "opposite" sawtooth waves in Figure
2.4.3 would lead to exactly the same transport rate according to formulae of the
type (2.4.7).
This is important to note, because several experiments, e g those of King
(1991), clearly show that the sawtooth wave with the steep front tends to generate
shoreward bed-load transport while the one with the steep rear generates seaward
bed-load transport. See Figure 2.4.4.
King’s observations show that sawtooth skewness matters, and hence that the
approach which leads to bed-load transport formulae of the type (2.4.7) is too
simplistic. The "next one up" from Equation (2.4.6) as a starting point is to assume
that the instantaneous transport rate is determined by the instantaneous effective
bed shear stress, ie
0.01 0.1 { 10
3s<T<l1ls.
where the findings in Sections 2.4.2 and 2.4.3, page 117 ff, seem to indicate that
the effective stress t’ can be calculated as the grain roughness stress 125.
Basing the transport rate on 1’(0,f) instead of u.(f) gives the potential of
capturing the effect of sawtooth asymmetry because the bed shear stress under a
‘sawtooth wave is likely to have the shape shown in Figure 2.4.5, with higher
(absolute) stress values under the rapidly accelerated half cycle.
Figure 2.4.5: Bed shear stress under a sawtooth wave for a boundary layer which is
laminar or has a constant (independent of z and ft) eddy viscosity. Although the
shoreward and seaward peak velocities are of equal magnitude, the shoreward peak
shear stress is considerably larger than its seaward counterpart. Qualitatively, the
reason is that the boundary layer has had less time to grow during the rapid shoreward
acceleration.
The influence of the steady flow component on the bed shear stress is
difficult to estimate from a time series of u(t) at a single level, and if u is
relatively large as in some rip currents or in a strong undertow, this calls for special
consideration, see e g Section 1.4.5 and Example 1.5.2.
If the steady flow component is weak however, in the sense that its influence
on the bed shear stress is small, it seems reasonable to try and derive t(0,t) from
Uco(t) by means of a simple transfer function based on our knowledge from simple
harmonic boundary layer flows.
For an oscillatory boundary layer with constant eddy viscosity (independent
of z and ¢ ), the bed shear stress corresponding to uo(t)=AM@e™
iot is given by
The phase shift @r is likely to depend on the Reynolds number A’w/v and
on the relative bed roughness _ r/A.
At present, it seems reasonable to approximate the instantaneous, effective
bed shear stress under waves of arbitrary shape (arbitrary uco[t]) by an adaptation
of Equation (2.4.12), for example
, ; d Uco
t(0,t) = eP f2.5 Arms ( COSOt Wp Uoo(t) + sings (2.4.13)
dt
Here, the friction factor f25 should be based on the representative semi-excursion
Anns = N21Urms/Wp , where Wp is the spectral peak angular frequency (2 7% Sp):
The bed shear stress time series corresponding to a given Ueo-record may then
1 5
based on the expression (2.4.13), in analogy with Equation (1.2.17), page 23.
As a simple digital filter to generate the time series 1’(0,tn), on a flat bed,
from Ueo(tn) one might correspondingly suggest
(2.4.14)
The corresponding time series for the instantaneous effective Shields parameter is
then given by
1
= fr5 Arms =
’ roe? ; Uoo(tn+1) Uco(tn—-1)
O'(tn) = eed:
ah ( COS@q Wp Uoo(tn) + singr Sate ity 5;age ge )
(2.4.15)
As indicated by the review of bed-load transport under sine waves in Section
2.4.3, page 117 ff, a respectable estimate of the instantaneous bed-load transport
rate is given by the Meyer-Peter formula with the effective Shields parameter
based on 62.5, see Equations (2.4.4) and (2.4.5), page 120.
In order to incorporate information about the changing direction of the bed
shear stress, Equation (2.4.5) must now be augmented to
For the half cycle average ©74, Equation (2.4.16) leads to results very similar
to Equation (2.4.2) which corresponds to the curve in Figure 2.4.2 with minor
deviations depending on the relative magnitude @25/@c of the peak effective
stress.
The expressions (2.4.15) and (2.4.16) are easily applied to any u..-record.
However, since the derivation of the model involves several assumptions and
simplifications, which cannot be checked completely at present, the model should
be calibrated as more data become available.
Application of the sediment transport model consisting, of Equations (2.4.15)
and (2.4.16) to skew waves over flat beds is illustrated in the following example.
Figure 2.4.6: Free stream velocity and calculated, dimensionless sediment transport
rates corresponding to the conditions in runs 473-485 and 478-520 of King (1991).
Legend: &: free stream velocity, * : ®(t) with gr = 20°, rectangle: ®(t) with
Qt = 30°, x: (ft) with or = 40°, -: O(t) withgr = 45°,
The calculations give the correct general magnitudes for ®74. Thus, for the
steep-front-halves (0<t<2.5s) the calculations give 7.09, 7.13, 1.18, and 1.19 in
order of increasing ©. The corresponding numbers for the last half of the waves
(2.5s<t<5s) are -1.04, -0.985, -0.855, -0.764. King’s measured values are /.09
and -0.60 respectively.
Inspecting the ratios between corresponding values reveals that the right
order of transport skewness is achieved by the model only when the phase shift @r
is set the maximum value of 45 °.
For a phase shift of 20° which corresponds roughly to the observations of
Sleath (1978) the predicted transport rate ratio is only /.05 as compared to the
measured value of /.7.
A complete explanation for this is not possible with the available
information, but it may be that the bed-load transport actually lags a bit behind the
bed shear stress, so that when Sleath observed a bed-load phase shift of @g = 20° ,
the bed shear stress might have been a bit more ahead of u.o(f) .
3.1 INTRODUCTION
The sea bed is very rarely flat. On the contrary, it tends to be covered by
sedimentary structures with a large range of sizes and of many different shapes,
and these structures: bars; dunes; ripples and animal mounds interact with the flow
in different ways.
The larger sedimentary structures such as bars modify the main flow pattern,
ie they make the incoming waves refract, diffract or break, and they reflect part of
the wave energy back into deep water. They may also determine the positions of
rip currents.
The small scale structures, such as ripples, have no immediate impact on the
main flow patterns, but they strongly influence the boundary layer structure and
the turbulence intensity near the bed. Hence, they have great influence on the
sediment transport.
The interaction between the large scale topography and the main flow will
not be considered here, but the dynamics and geometry of the bedforms (mainly
vortex ripples) will be discussed in detail together with their influence on the
boundary layer structure. Finally, an attempt is made to quantify the influence of a
given bed topography on the flow by a single linear measure, namely the
equivalent Nikuradse roughness.
If the flow is too weak to cause appreciable sediment motion (62.5 $ 0.05),
the bed topography will be dominated by relict bedforms from previous more
vigorous events and if no such events have occurred recently, the topography will
be dominated by bioturbation.
Under flows of intermediate strength (0.05 S @25S 1.0 ) the bed will be
active and will be covered with bedforms which are more or less in equilibrium
with the flow conditions. However, the shape of these bedforms will depend on
the detailed nature of the flow. If the flow is purely oscillatory and almost
symmetrical, the bed will be covered by regular, long crested vortex ripples, which
are so named because a vortex is formed twice every wave period in the lee of their
crests.
Shore
break
Reformed Swash
| Breaker zone waves
The regular, symmetrical flows required for the formation of vortex ripples
occur in nature seaward of the breaker zone and in bar troughs where the waves
have reformed after breaking on the bar, see Figure 3.2.1.
Vortex ripples are rarely seen in the breaker zone. Here, the bed tends to be
either flat or covered by megaripples. Megaripples or lunate megaripples (Clifton,
1976), are irregular features with typical lengths of the order one to two metres and
heights of the order ten to twenty centimetres, and with rounded crests. There is no
regular, rhythmic vortex shedding associated with megaripples, but occasionally a
large plume of turbulent, sediment laden water rises from a certain spot due to
complicated instabilities. These sand plumes can be of the order one metre high.
Vortex ripples are sometimes found superimposed on megaripples (Southard
et al 1990).
Under very vigorous flows (825 2 1.0) vortex ripples cannot exist even under
perfectly regular, oscillatory flows. However, megaripples have been found to
exist under such flow conditions (82.5 up to about 2.4) both in large wave flumes
and in oscillating water tunnels, see e g Ribberink & Al-Salem (1989).
Q(x) + a dx, and hence, the accumulated volume of solid sediment is — < dx
dp
nh _ _aQ
.e (3.3.1)
Note that this equation is only strictly valid when sediment storage in the
water column is insignificant. See Equation 5.3.1, page 223.
ak f py ’(x -ct)
n (-c) aie +
through cross sections where zp(x,t) = 0. The result (3.3.2) shows that if the
bedforms move downstream (c > 0) with constant shape, the sediment transport
rate must be maximum over the crest and minimum over the trough, ie Q varies
in step with zp, and the difference in transport rates between the crest and trough
sections is
Figure 3.3.2: If the peak in the sediment transport rate is shifted dx away from the
crest, the bedforms will be growing or diminishing at a rate determined by 6x/A.
Ole a
OG = 0.
where the shift 5; is the downstream distance between the bedform crest and the
point of maximum sediment flux, see Figure 3.3.2.
Under these conditions, the continuity Equation (3.3.1) gives
where we can cancel the common factor — ike *%—) and solve for the speed of
migration to get
Secondly, it is hard to predict the shift 5, between the bedform crest and the
point of maximum transport.
For a comprehensive discussion of the formation and growth of vortex ripples
see Lofquist (1978) and Sleath (1984).
3.4.1 Introduction
Vortex ripples are of special interest for coastal sediment transport studies
because their influence on the boundary layer structure and the sediment transport
mechanisms is very strong. That is, over vortex ripples, the suspended sediment
distribution will scale on the ripple height, while for other bedforms like
megaripples and bars, the suspension distribution will scale on the flat bed
boundary layer thickness which is much smaller than the height of those bedforms.
The shape and size of vortex ripples was first studied in detail by Bagnold
(1946), who also described the flow and the sediment transport mode above them.
Inman (1957) investigated their natural occurrence in a comprehensive field study,
and their development from a flat bed and their adaptation to new flow conditions
have been studied comprehensively by Sleath (1984) and Lofquist (1978)
respectively. All of these sources provide excellent illustrations of ripple shapes
and associated flow patterns. The occurrence of different overall ripple patterns, i e
neatly two-dimensional versus confused three-dimensional has been discussed by
Carstens et al (1969), Nielsen (1979) and Sleath (1984).
Figure 3.4.1; The size of vortex ripples is closely linked to the orbit length 2A of
the wave-induced fluid motion near the bed.
2.0
1.0
0.40
0.20
Figure 3.4.2: Field data of relative ripple length versus mobility number based on
significant wave heights for field data. The curve corresponds to Equation (3.4.3) .
Based on the field data of Inman (1957), Dingler (1974) and Miller & Komar
(1980), Nielsen (1981) suggested
y
ep |£23-=.0.37 In 8u fae
>> 1000 + 0.75 In’y
for ripples under field conditions, i e, under irregular waves.
Parabola
Figure 3.4.3: If the maximum slope on the ripple profile is assumed equal to tan 9,
the height to length ratio must for geometric reasons be O.Stan @ for a triangular
ripple, and 0.25 tan @ for a parabolic ripple.
It turns out that vortex ripples have a maximum height to length ratio (n/A)
which falls within the range corresponding to these two idealised geometries, see
Figure 3.4.4. The data indicate a maximum steepness, at vanishing flow speed, of
When the grain roughness Shields Parameter 625 exceeds 0.2, the flattening
effect of the contracted flow over the crest tends to increasingly overpower the
constructive effect of the lee vortices which build up the crests by scooping sand
towards them.
As a result, n/A becomes a decreasing function of @25, see Figure 3.4.5,
and eventually the ripples are flattened completely. This happens at a 025-value of
approximately /.0.
0.50
0.40
0.30
0.20
0.10
0.05 0.10 0.15 0.20 0.30
Figure 3.4.4: The maximum value of n/A for a given bed material lies in the range
0.25 tang <(N/A)max < 0.5 tan , where the two limiting values correspond to an
equilateral triangle and to a parabola respectively, both with a maximum local
steepness of tang, see Figure 3.4.3.
The trend of the laboratory data (regular waves) shown in Figure 3.4.5 is
reasonably described by
0.40
0.20
0.10
0.04
LEGEND T(sec) d(mm) s
+ DHI TUNNEL 7.7-13.9 0.082 2.65
on " 4.4-9.0 0.2+0.135 2.65
0.02 |4 CARSTENS et AL 3.55 0.19 2.65
A " " 3.55 0.30 2.65 | FLAT BED |
O YALIN & RUSSELL 1.82 0.48 1.19 | +0+++4 |
V KENNEDY & FALCON 1.8-2.8 1.00 1.035 | |
O NIELSEN (1979) 1.70 0.082 2.65 |
Figure 3.4.5: Under gentle flow conditions, the ripple steepness N/A is only limited
by the angle of repose, but for 62.52 0.2 it becomes a decreasing function of 02.5.
For field conditions (irregular waves), the ripple steepness tends to be smaller
than for regular laboratory waves. On the basis of the field observations of Inman
(1957) and of Dingler (1974), Nielsen (1981) suggested the formula
Dingler suggested, on the basis of his field observations, that the ripples
disappear when y reaches a value of approximately 240, and that n/A could, in
general, be described as a function of the mobility number y.
However, while this may be true, when only quartz sand is considered, it
does not hold in general. This can be seen by comparing Figures 3.4.6 and 3.4.5.
In Figure 3.4.6 1/A has been plotted as a function of w for the same data
as in Figure 3.4.5 and it is quite obvious, that ripples formed by sediments of
different densities form different trends in this plot.
0.40
0.04
T(sec) d(mm)
0.02 | + DHI TUNNEL 7.7-13.9 0.082
eon " 4.4-9.0 0.2 0.135
f CARSTENS et AL 3.55 0.19
4 " " 3.55 0.30 O
O YALIN & RUSSELL 1.82 0.48
0.01 | SO KENNEDY & FALCON 1.8-2.8 1.00
O NIELSEN (1979) 1.70 0.082
4 10 20 40 100 200
Figure 3.4.6: Ripple steepness as function of the mobility number for sediments of
different densities.
0.40
0.20
0.10
0.04
0.02
conditions however, where most sediments will have densities close to that of
quartz (s ~ 2.65), it is quite reasonable to describe n/A as a function of the
mobility number.
Consequently, since the relative ripple length X/A, is also reasonably well
described as a function of y (cf Figure 3.4.2, page 137) one might expect that the
relative ripple height 1/A can be described as a function of y. Figure 3.4.7
shows that this is indeed the case.
Nielsen (1981) suggested the formula
was suggested, where y and A should be based on the significant wave height.
0.01 On 1
From their field study (weak to moderate flows in 22m depth on the shelf off
Sable Island), Amos et al identified eight different bedform regimes and mapped
them in terms of the grain roughness Shields parameter 625 versus a current
Shields parameter given by
ae 14 p 0.003 uioo
p (s-l) gd
where Uioo is the mean current velocity measured one metre above the bed. See
Figure 3.5.1.
The bedforms in strong rip currents are not well researched for the good
reason that such rips are extremely dangerous work environments. However,
qualitative observations indicate that the bedforms tend to be current-dominated,
e, they have steeper offshore faces and migrate offshore.
In fairly deep bar troughs over which waves have reformed after breaking on
the bar, very regular wave ripples are often found. One such case was reported by
Nielsen (1983) and Wright et al (1986). The depth was 1.2m to 1.5m, the
significant height and period of the reformed waves were 0.5m and 7s
respectively, and the depth averaged longshore current was approximately 0.5m/s.
In shallower bar troughs or rip feeder channels however, the conditions are
often quite obviously current-dominated, with the bedforms resembling current
dunes. This may be due to relatively weaker wave activity, Put,it may also be due
to the shallower depth itself via the current Froude number Vv’/gD. The stability
analysis of Engelund (1970) indicates that current dunes are generally formed at
Froude numbers between 0.5 and 2.0, but are unlikely to form at Froude
numbers below 0.25.
With respect to the megaripples which are sometimes found in oscillatory
flows of high intensity (@ 2 1, see Section 3.4.2) , Arnott and Southard (1990)
found that they generally disappeared if a moderate current (<u> > 0.05m/s) was
superimposed. Their experiments were carried out in an oscillating tunnel with
0.09mm sand and a fixed wave period of 8.5s.
Coffey (1987) observed bedforms, current profiles and sediment suspension
profiles in depths between one and two metres in the mouth of the Port Hacking
estuary south of Sydney and found ripples superimposed on the shoals. These
ripples were in the relative height and length ranges : 50 < n/d < 200 and 300 <
A/d < 1540, which corresponds to the expected ranges for current ripples. The
relative current strengths ux/ ux for these observations were in the range [0.15;
0.72).
3.6.1 Introduction
In order to formulate simple models of natural flows, it is generally necessary
to apply a simplified description of the bed geometry, and in the extreme, we often
try to summarise the bed geometry in terms of a single length. Most commonly, the
chosen length is the equivalent Nikuradse roughness or briefly the hydraulic
roughness, r.
The only bed geometry for which the definition of the roughness is obvious is
a layer of densely packed spheres for which the roughness equals the grain
diameter, r=d. For all other geometries, the definition is indirect.
145
AND SEDIMENT TRANSPORT
Chapter 3: Bedforms and hydraulic roughness
That is, the roughness is determined from the structure of a purely steady
flow above the bed, see Equation (1.5.9), page 66.
The following section deals with the concept of the equivalent Nikuradse
roughness of natural sand beds, particularly those exposed to oscillatory flows. It is
found from the available friction and energy dissipation data that the roughness of
such beds are generally one to two orders of magnitude larger than that of sand
paper with the same sand size.
This great roughness is in many cases obviously due to bedforms which
generate roughness of the order of their height (r= 1 ).
However, also flat, mobile sand beds dissipate wave energy at a high rate and
thus, in this sense, appear very rough, particularly at high flow intensities where a
substantial layer of sediment is in motion. Thus, flat beds in oscillatory sheet-flow
(252 1.0) generally exhibit roughness values (based on total friction or on total
energy dissipation rates) of the order 100 to 200 grain diameters.
This is somewhat surprising because the corresponding roughness in steady
sheet-flows are generally one order of magnitude smaller. It therefore prompts the
question whether a substantial part (more than half) of the energy dissipation
measured by Carstens et al (1969) and Lofquist (1986) on flat sand beds, could
have been due to other mechanisms than skin friction. If that is the case, the
obvious candidate for an additional dissipation mechanism under waves is
percolation, see e g Sleath (1984) or Dean & Dalrymple (1991).
The suspicion above is reinforced by the observations of sediment transport
rates over flat beds under waves, which indicate that the effective stress, with
respect to sediment transport, is considerably smaller than the bed shear stress
measured either directly or via the rate of energy dissipation.
Thus, the bed-load data analysed in Sections 2.4.2 and 2.4.3, page 117 ff,
indicated that the effective sediment transporting stress corresponds to a roughness
of only about 2.5d50, while energy dissipation measurements under similar
conditions indicate roughness values of the order 100d50.
This difference between steady and oscillatory flows remains unexplained,
except for the possibility that it may be due to percolation.
0.1
0.01
0.01 0.1 1
Figure 3.6.1: Measured wave energy dissipation factors from Carstens et al (1969)
and from Lofquist (1986). Legend: bar: Lofquist 0.55mm sand, +: Lofquist 0.18mm
sand, *: Carstens et al 0.19mm, rectangle: Carstens et al 0.30mm, x: Carstens et al
0.59mm. All of the above correspond to equilibrium ripples while the triangles
correspond to artificial flat beds where measurements were taken before ripples had
time to form (Carstens et al ).
(fe) which are plotted in Figure 3.6.1. fe is expected to be smaller than fw, but the
differences are generally small compared to the scatter see Figure (1.2.9).
0.1
0.01
0.01 0.1 1 10
Figure 3.6.2; Steady flow friction factors over natural sand beds with median sand
sizes in the range [0./9mm; 0.93mm] and depths in the range [0./4m, 0.33m]. Data
from flume experiments by Guy et al (1966). Legend: filled rectangle: plane bed and
transition., + : rippled beds, *: dunes, open rectangle: antidunes, x : standing
waves, triangle: chute and pool flow.
Roughly speaking, these friction factors for sand beds in oscillatory flows are
an order of magnitude larger than those found in steady flows; namely between
0.04 and 0.4, compared to a range of about [0.005; 0.04] for steady flows, see
Figures 3.6.1. and 3.6.2.
When data from artificially flattened beds at low flow intensities (02.5 < 0.25)
are excluded, the range of observed energy dissipation factors for sand beds under
oscillatory flows is
t== 5 Pf <u>
2 3 =
t=pgDI
was fitted to steady flows over beds of densely packed sand grains, at sufficiently
large grain Reynolds numbers, then the zero intercept level zo was approximately
one thirtieth of the sand size (cf Schlichting 1979)
i, == SUZ (3.6.3)
Figure 3.6.3: Hydraulic roughness derived through Equation (3.6.5) for sand beds
with different grain sizes in a variety of steady flows. Same legend as for Figure
S102.
If not the whole velocity distribution, but only the depth-averaged velocity
<u> together with the depth D and the hydraulic gradient J has been
measured, a crude estimate of the roughness can be obtained based on a total-flow
relationship of the form
r= 50d (3.6.6)
1000
10
0.01 0.1 1
Figure 3.6.4: Hydraulic roughness corresponding to the fe-data in Figure 3.6.1. For
rippled beds the roughness is generally in the range [100ds0; 1000dso0], and if the
trend of the limited flat bed data can be extrapolated beyond 62,5 = 1, the bed
roughness under oscillatory sheet-flow is also of the order J00ds0 or more. Same
legend as in Figure 3.6.1, page 147.
2
r=Settee ti:A
0.0081 (3.6.7)
If only the peak bed shear stress f or the time-averaged wave energy
dissipation D ee has been measured, the hydraulic roughness may be inferred via
the wave friction factor fy or the energy dissipation factor fe, Thatis, fw or fe
is found from Equation (1.2.18) or Equation (1.2.26) respectively, and the
roughness is then found by applying a wave friction factor formula like Equation
(1.2.22) in reverse ie
Sls)
Inf + 5.977
r= a[ 5.213 (3.6.8)
The roughness values derived with this formula from the energy dissipation
factor data in Figure 3.6.1 are shown in Figure 3.6.4.
1000
0.1
AND SEDIMENT TRANSPORT 153
Chapter 3: Bedforms and hydraulic roughness
We see that, except for (artificially) flat beds at low flow intensities
(25S 0.20), this roughness is of the order 100d50 to 1000d50.
This is interesting because, the steady flow data in Figure 3.6.3 indicate
roughness values as low as 2ds5q occurring in the upper (steady) flow regimes,
and Wilson’s (1989) formula (3.6.6) predicts only r= 5d __ for Shields
parameters of the order unity, see Figure 3.6.5.
This is of the order ten times more than predicted by Wilson’s (1989) formula
for steady flows.
No explanation has been given so far to this large difference between the
roughness values of sand beds under steady and oscillatory sheet-flows, except
that part of the energy dissipation or total drag under waves may be due to
percolation.
In terms of the grain roughness Shields parameter 025 the roughness,
corresponding to total drag on flat sand beds under oscillatory sheet-flow, may be
expressed as
If this is a valid indication of the boundary layer structure in general over flat
beds of loose sand, it has important implications for the relevance to flow over
natural sand beds of many of the existing models of oscillatory boundary layers.
The reason is that many of these have been developed to match experimental
evidence from flows with a relative bed roughness of less than 0.03, and are not
suited for flows with r/A>0.1, see Section 1.3, pp 40-52.
0
9
8
7
6
5
4
3 Figure 3.0.75 The
hydraulic roughness of
rippled beds is between
2 one and three _ ripple
heights as long as the flow
1 is not too strong
(62.55 0.5). Same legend
as in Figure 3.6.1.
0
0 0.2 0.4 0.6 0.8
Indeed , r/1_ is seen to vary very little up to 625= 0.5 even though the
ripple steepness generally starts to decrease from about 825=0.25. See Figures
3.4.5 and 3.4.7.
The momentum transfer between a sediment laden flow and the bed can be
seen as consisting of two components: The usual drag contribution, which would
also be there for a fixed bed with the same geometry, and a component due to the
transfer of momentum by accelerating and subsequently crashing sediment
particles.
Correspondingly the roughness may be seen as consisting of a fixed bed part
rf and of a moving sediment part rs:
r = Mft+rs (3.6.12)
The magnitude of the fixed bed contribution rg for rippled beds may be
0.01 0.1 1
Figure 3.6.8: Relative ripple roughness as function of the grain roughness Shields
parameter for oscillatory flow data. Same legend as in Figure 3.6.1.
inferred from some of the experiments which have been carried out with ripple-like
fixed beds.
It has generally been assumed that the roughness of a rigid ripple profile
would be proportional not only to the ripple height, but also to the steepness and
hence relations of the form r= const-n 1/A or rA/n* = const have been
sought.
Thus, Bagnold’s (1946) sharp crested, parabolic ripples correspond to r Wr
=20.3, while the triangular concrete ripples of Jonsson & Carlsen (1976)
correspond to r WA =10.9 and 7.4 for TestI and Test II respectively.
Figure 3.6.8, where r Wn? is plotted versus 25 for the available movable
bed data, indicates that a value for r Wn? of about 8 is appropriate for vortex
ripples at low flow intensity.
In order to account for the roughness contribution (momentum transfer) from
the moving sand over ripples, we may add a term corresponding to the expression
(3.6.10), page 154, for moveable flat bed roughness and get the following estimate
Estimated
(given
ripples)
[cm]
r
10
Observed r [cm]
Figure 3.6.9: Predicted versus measured roughness values corresponding to the
Se-data of Carstens et al (1969) and of Lofquist (1986). The predictions are based on
Equation (3.6.13) with measured ripple data inserted.
where the ripple geometry may be measured or estimated from Equations (3.4.5,
3.4.7) or (3.4.6, 3.4.8).
Figure 3.6.9 shows a comparison between this formula and the roughness
values derived via Equation (3.6.8), page 153, from the fe-data of Carstens et al
(1969) and Lofquist (1986). The figure shows reasonable agreement.
We note that the nature of the first term makes it impossible to estimate the
roughness by a simple formula of the form r/d = F(@25) or r/d = F(w) as
suggested by Raudkivi (1988). The reason is that the ratio 17/Ad_ varies strongly
with the wave period for a fixed value of 25 or of wy, see the discussion of
Nielsen et al (1990).
While, according to Figure 3.6.9, energy dissipation rates based on Equation
(3.6.13) are in good agreement with experiments, there is reason to suspect that the
roughness determined from (3.6.13) and (3.6.10) for flat beds are too large for the
estimation of the vertical scale of suspended sediment profiles, see Section 5.4.9,
page 255.
Also, the discussion of bed-load transport rates and corresponding effective
bed shear stresses for flat beds under waves in Section 2.3 indicates that bed shear
stresses determined from a roughness given by (3.6.13) or (3.6.10) for flat beds are
too large.
The optimal choice of roughness for both bed-load transport rates and
suspension distributions over flat sand beds under waves lies closer to Wilson’s
steady flow formula (3.6.6), page 151, with 6 replaced by 025.
Thus, for purposes other than estimation of frictional dissipation, it might be
appropriate to replace Equation (3.6.13) by an equivalent formula with the
"moving grain contribution" given by Wilson’s expression for example
4.1 INTRODUCTION
Transport of suspended sand is of great practical importance in coastal and
fluvial engineering. Therefore, the mechanics of suspended sediment transport
have been the object of substantial research efforts through several decades.
However, there is still much to be learned about the concentration magnitude and
about the distribution of suspended sediment in waves, in currents, and in
combined wave-current flows.
In order to describe suspended sediment distributions we must first seek to
understand the behaviour of suspended sand grains in different types of flows.
That is the object of the following sections.
The most important flow structure in connection with suspended sediment is
that of a vortex with horizontal axis because such vortices are able to trap sand
grains and carry them along over considerable distances, see Figure 4.1.1 (right).
In contrast, a wave motion where the fluid orbits are also circles or ellipses,
has no trapping capability, see Figure 4.1.1 (left). The mechanism of sand trapping
in vortices was first pointed out by Tooby et al (1977), but we shall discuss it in
detail in Section 4.6, p 181, and refer to some empirical evidence of its importance.
In coastal sediment transport the nearly homogeneous oscillating flow field
induced by surface waves is, of course, important, and there have been many
attempts to derive a suspension-maintaining mechanism from this flow structure.
For example, the delay distance idea of Hattori (1969) and the time lag hypothesis
of Bhattacharya (1971).
The former of these has little physical substance but the latter will be proven
and quantified in Section 4.5.
A third mechanism, which might influence the sediment settling rate under
waves, is the possible settling velocity reduction due to nonlinearity of the drag
force. We shall see in Section 4.5.5, through an analytical approach to this
mechanism, that the changes in settling rate due to the wave motion itself are either
purely oscillatory, and therefore without net effect, or they are so small that they
have no practical importance.
Wave Vortex
sediment
water ath
particle
U, = Wo + Wo (€7) u,=Oo+w,O(e)
Figure 4.1.1; While vortices and wave motions are similar in so far as both have
circular or elliptical fluid orbits, they differ greatly with respect to their influence on
suspended sediment particles. Sand will settle almost unhindered through the wave
motion, but it will tend to get trapped in vortices.
The second step takes terms of magnitude ewe into account, where € is the
relative magnitude of the fluid accelerations
1 du | (4.1.1)
metieruceulelt
Finally, as a third step, effects of the order Wo are considered. Hence the
sediment particle velocity vector us is built up as a series of the form
Us = U+Wottritunt...
= ut vol
(SJreuneet +. (4.1.2)
us = U+Wo (4.1.3)
Secondly, the fluid accelerations induced by waves are generally much smaller
than the acceleration of gravity so that the perturbation parameter € is
conveniently small.
Us
Wo = V A(s-1) gd (4.2.2)
3CpD
STOKES
Cp=24/R
Fp = 3npvdwo (4.2.4)
_ (s-lgd Ps
am 18v
settle faster than the typical settling velocity of the individual particles. This effect
may be the reason for the unexpectedly large, measured settling velocities for the
smaller ( d<0.04 cm.) particles in Figure 4.2.2.
These small, angular sand grains may have been settling in clouds.
Alternatively they may have obtained abnormally large averaged settling velocities
by hitch-hiking during parts of their trip in the wakes of larger particles.
These effects must be considered when settling velocity measurements from
settling tubes, where many sand sizes are dropped together, are interpreted.
w(cm/s)
Figure 4.2.2; Settling velocities (w0, ws0, w90) determined with a settling tube for
narrow sieve fractions of a sediment sample from the swash zone at Palm Beach,
Sydney, Australia. This sample contained a large fraction of shell hash, so the
measured settling velocities lie considerably below the curve given by Gibbs et al for
quartz spheres, Equation (4.2.6). The relatively higher, measured values for the
smaller grain sizes are probably due to the particles settling as a cloud rather than as
individual particles.
The following section briefly discusses the various terms in the equation of
motion for suspended sediment particles. Subsequently, the equation is brought
into the most convenient form for deriving a perturbation solution in the form of
Equation (4.1.2), page 163.
Consider a quasi-spherical particle with diameter d and relative density s
which moves under the action of gravity and pressure gradients and drag from the
surrounding fluid.
We neglect the history term of Basset (1888) which accounts for changes in
the fluid drag due to changes in the flow structure around the particle, and neglect
any effects of the particles rotation such as the Magnus effect (Magnus 1853). The
reader is referred to, for example, Clift et al (1978) for a more detailed discussion.
With those simplifications applied, the equation of motion reads
~y
Co _ bie; (4.3.2)
Wo
The last term in Equation (4.3.1) corresponds to fluid pressure on the added
hydrodynamic mass, and the coefficient Cy may be assumed to have a value
close to 0.5 which is the theoretical Cy-value for a sphere.
By introducing the pressure gradient vector in the form
Vp = pe - pm (4.3.3)
and rearranging, the equation of motion becomes
dus _ dus
i ae
= £ (usw) + (utur)-V(u+ur)
a ou
= ar + =e + u-Vu + (ut+uy)-Vuy + urVu
_BA du,asdale
duy om Vu (4.3.5)
dur Our
“aks = 5 + (u + ur)-Vur (4.3.6)
d ,
(4.3.7)
and with
we get
dur + Cp
NEARLY MOIan (4.3.9)
diy arpdcd(sveiGi) hintvo Gesoawn
It is convenient to eliminate the drag coefficient Oy from this equation. To
this end, we apply the definition of C,, , Equation (4.2.1), page 164, which gives
I-¥
du u
ae «a(t Wied Oey 5 me pugLE (4.3.11)
6)
This general differential equation for the relative sediment particle velocity
ur is the basis for the analysis of the motion of suspended particles in the
following sections.
The nature of the solutions is found to depend in a very important way on the
third term u;-Vu , which, in essence, represents the effect of the small scale flow
structure.
This term is negligible in a pure wave motion (for woT/L << 1), and this is
the reason why a pure wave motion has very little effect on the sediment settling
rate, see Figure 4.1.1, page 162, and Section 4.5.5, page 178.
On the other hand, for a sediment particle inside a vortex with diameter of the
order of ten centimetres, the term u,-Vu is very significant, and it provides the
vortex trapping mechanism which is discussed in Section 4.6, pp 181-189.
Equation (4.3.11) is linear for y=1, and exact analytical solution is
therefore possible, for small particles for which y= 1, see Figure 4.2.1.
For larger sediment particles, which have y-values in the interval [0;/],
perturbation solutions of the form
Us =~ U+uUr = U+wWotn
t+ unt.
are provided in Section 4.5 for homogeneous flows and in Section 4.6 for particles
in vortex flows.
d
et
OO 2 =e 2 (4.4.1)
(a)
which results from Equation (4.3.11) when the fluid is at rest (or moving with a
steady, uniform velocity), and the particle is so small that the flow around it is
laminar and Stokes law applies (y= /).
This equation is equivalent to
i (ur—-Wo)
—, = eR
wa (ur ae— Wo) (4.4.2)
where the vector g has been written as gwo/Wo, and this equation has the solution
From this we see that deviations from the terminal settling velocity wo decay
exponentially with a time scale of wo/ag which, for beach sand with relative
density 2.65 and settling velocity 0.02m/s, is approximately 4 - / Os.
This means that, starting from rest, such a sediment particle will be moving at
99% of its terminal settling velocity after only 0.002 seconds.
In other words, this means that deviations from the first approximation
us = U + Wo (4.1.3)
4.5.1 Introduction
In the previous section we saw that sediment particles in a resting fluid will
attain a relative velocity approximately equal to the terminal settling velocity after
only about one hundredth of a second.
However, if the fluid motion is accelerated, deviations of the order ew, will
remain where
€ = 1 du | (4.1.1)
g at
The nature of these deviations is the subject of the following section.
Ur = WortUri (4.5.1)
dur uri du
dt og—
A Wo + rl un:Vu = —-Aa—
— WoVu
i) (4.5.2 )
dt
This equation is particularly simple to solve if the flow field may be assumed
homogeneous so that the gradient matrix is zero, Vu = 0.
As an example, a solution (Equation 4.5.27, page 176) is given in Section
4.5.4, for settling through homogeneous, oscillatory flow.
Grom.
For sediment particles moving in a homogeneous, oscillatory flow, it can then
be shown that all effects of order € Wo are analogous to the results for small
particles derived from the linear Equation (4.5.2). These effects are all purely
oscillatory, so there is no settling velocity reduction of the order € Wo. The settling
velocity reduction due to non-linear drag is of the order Wo .
We introduce the following dimensionless variables
U = u/R®
Uy = ur/Wo
T = ot (4.5.3)
where R is the typical scale of the water motion and @ its radian frequency.
Then Equation (4.3.11) for particles in a uniform flow becomes
Wr 08 yi-ry, _ 28 _,,Rw dU
aT WoW WoW a Wo® dT See,
or
dU.
eT en " a eg
dT Wo® WoW WoW dT oe
Wh . : ° :
where €=Ro'/g. We seek solutions to this equation in the form of a
perturbation series
ur = Wotunt+unt... (4.5.6)
Up = Up eu te Up F
Ur ib= “)
(0 Un
+ (wn) + 2(Un
| oe (4.5.7)
and for the linearisation of the drag term in Equation (4.5.5) we find
. i 7 waiEy) (4.5.11)
WoW —] WoW | -8
UM,) ,, Sasi
dF \Wn ried nN WwW
(4.5.12)
dT Wo® (2-y)Wr1 WoW dT
a(U si y ‘ (-yUn Wn
d (Un no.ge( Urs E\patiocons. Te) ud
aT | Wo lop] WoW = U2, + Peas Ww,
(4.5.14)
+ a 0
u= ut) = rs cea (4.5.16)
To find the first order solution wri, we use the vertical part of Equation
(4.5.13) which becomes, in this case
Ro Wo, : =i
Wri = sin (wf — tan Bz) (4.5.21)
Sent ok oP
+B?
where
Wo ®
Bz = ws
Si
(Q-) agSo (4.5.22 )
We note that for sand size particles in water the value of Bz is of the order
of magnitude 10-*. Hence, for practical purposes, the expression (4.5.21)may be
approximated by
Wo du
Us A u + Wo oe eS
Q-ve dt (4.5.26 )
Us = u(t-d:) + Wo (4.5.27)
5; i eee (4.5.28)
(2-y)g
or approximately wo/g, which is of the order 0.00]/s to 0.0]s for beach sand.
We note that the solution (4.5.21) has the time average zero so there is no
settling velocity reduction of magnitude ew.
There is however a net reduction of the settling velocity of order €°Wo when
the drag force is non-linear. We find its magnitude by inserting the first order
solution (4.5.21) into the second order equation (4.5.15).
In order to simplify the algebra however, let us write the first order solution
(4.5.21) in the form
da
te Zot: 8 02) we
a
res 2-3y+7 rt sin? w t (4.5.30)
dt oO we 2
or
da
Tt Ww
=
t
= mg 2+7 21 4 oso @ 2!) (4.5.31)
Wo Wo D 2
a
wr2 =
ey
aes
2
rj (4.5.32 )
Wr,
/Wo
-
]
Figure 45.1: Settling
velocity reduction for a
Imm brass sphere in
vertically oscillating water
———— Ho’s Numerical Solution Wo = 0.39mls, @ =
Equation (4°5:33) 22.7rad/s. |Measurements
and numerical solution by
Ho (1964).
This shows that the settling velocity reduction due to non-linear drag has the
order of magnitude €’wo. The solution (4.5.33) is compared to measurements and
numerical results from Ho (1964) in Figure 4.5.1.
Note that even though the derivation was based on the formal assumption of
€ << 1, the agreement with Ho’s experiments is good up to € = 3.
Us = U + Wo (4.1.3)
are either purely oscillatory, and therefore without a time-averaged effect, or they
are for all practical purposes negligibly small.
In most practical cases the water motion above the bottom boundary layer
under waves can be considered quasi-uniform as far as the motion of suspended
sediment is concerned.
More precisely, this is the case when the distance settled by a sediment
particle over one wave period is small compared to the vertical scale on which the
wave induced velocities change. That is when
uote he
WR; cos wt
4.5.35
(4.5.35)
We now use Equation (4.5.13), page 174, to find the horizontal and vertical
components of the first order solution uy1, [us(t) = u(t)+wotun(t)+...]. In
analogy with the solution (4.5.21) for purely vertical motion, we find
ur) =
g M+ (4.5.36)
sete
Ro Wo_
ie sin (wt
— tan 'B:Bz)
where
oe _ (4.5.37)
and 8, is given by Equation (4.5.22), page 175.
Note that the time average uy; =0 . Therefore, as indicated by Figure 4.1.1
page 162, there is no net drift induced by the waves at order ewo.
The solution is simplified considerably when the particle is so small that
y= 1 in which case we will also have Bx, Bz << 1. Then we get
2
Wo COS Wt
ur) = : ee Wog ales
du
dt 4.5.38
(4.5.38)
Wo Sin @t
and hence
Wo du
Hes ul) Wes = eins
= (4.5.39)
This result proves and quantifies Bhattacharya’s (1971) time lag hypothesis
for small particles in the two-dimensional case. For larger particles (y=0) the
time lag for the vertical component is smaller, namely wo/(2-)g instead of
Wo/, g.
ud
i gg Figure 4.5.2: In a
homogeneous flow field
where the water particles
follow elliptical orbits, the
component u,i of the
sediment velocity is
directed away from the
orbit centre and _ thus
represents a _ centrifugal
effect.
For larger sediment particles (y~0) there is a net delay of settling wy due
to the non-linearity of the drag force. It can be evaluated by inserting the first order
solution (4.5.36) into the second order Equation (4.5.15), page 174, and taking the
time average.
The solution is analogous to Equation (4.5.33). Hence, the settling velocity
reduction Wr is ofthe order €’wo.
In most practical cases the reduction is insignificant. For example, for a
typical wave motion with vertical semi-excursion Rz = 0.1m and @ = 1s! we
find from Equation (4.5.33).
wn = 6:10 wo (4.5.41)
4.6.1 Introduction
One of the most important mechanisms for entraining and suspending
sediment is the trapping of sediment by vortices with horizontal axes. See Figures
4.1.1, page 162, and Figure 4.6.1.
This mechanism was first illustrated experimentally by Tooby et al (1977).
The reader is referred to Tooby et al’s magnificent photographic illustration. It
shows a heavy particle as well as air bubbles trapped on circular paths inside a
water filled cylinder which rotates around its horizontal axis.
In the following we shall develop an analytical description of the
phenomenon along the same lines as used in the previous sections. That is, the
sediment particle velocity is written in the form
Us = U+WotUnt+urt.... (4.1.2)
= u t+ wo]
[Steuart |
where uw is the fluid velocity vector and wo _ is the still water settling velocity. The
perturbation parameter € is the ratio between fluid accelerations and gravity
er =t—-r4 (4.1.1)
us = U+Wo (4.6.1)
which is that of a rigid body rotating around the origin with angular velocity
w, see Figure 4.6.1. Applying Equation (4.6.2), we find the approximate sediment
particle velocity
-Z 0 4
us = U+ Wo = of) + ee = (sre) (4.6.4)
Figure 4.6.1: Sediment particle with settling velocity wo moving in a forced vortex
with velocity field u given by Equation (4.6.3). The particle paths corresponding to
u+Wo are circles and thus the time-averaged velocity of a sediment particle will (at
this level of approximation) be zero.
Hence, the sediment paths are analogous to the fluid paths except for a
horizontal shift of wo/@. That is, since all circles around (0, 0) are fluid particle
paths, any circle around (wo/@, 0) is a possible sediment path, see Figure 4.6.1.
Note that the sediment paths in a pure, deep water wave motion are not
closed, although the water particles move in circles like in the vortex, see Figure
4.1.1, page 162, and Section 4.5.5, page 178 ff.
It is important to understand the difference between wave motion and vortex
motion with respect to sediment suspension. It stems from the fact that the wave
motion is essentially homogeneous while the vortex motion is not.
Quantitatively, in terms of the equation of motion (4.3.11), page 169, the
difference is that the term u;-Vu is approximately 0 for the wave motion while
it is significant for sediment particles inside vortices.
Tooby et al (1977) showed experimentally that small sediment particles in a
vortex do, in fact, follow circular paths very closely and only very slowly spiral
away from them. This spiralling process is discussed in Section 4.6.3, page 186.
The fact that the settling is strongly delayed in a vortex flow field has been
noticed by Reizes (1977), who found the phenomenon in a numerical study. He
concluded that the sediment particle must tend to spend more time in the upward
moving parts of the flow than in the downward moving parts.
This is indeed the case. The sediment path shown in Figure 4.6.1 lies entirely
in the “upward moving” part of the vortex.
For buoyant particles or air bubbles the centre of the “sediment path” would
lie on the negative x-axis. Thus, they would spend the majority of the time within
downward moving fluid.
The described trapping mechanism will work for all sand grains with settling
velocity smaller than the maximum velocity in the vortex.
The next question to be asked is whether the trapping is a feature of the rather
unnatural, forced vortex only.
The answer is no for the following reason. In general, the velocity field of a
circular vortex with angular velocity @ can be written
For the components us and ws of this velocity field, we therefore have the
following symmetries
and
This means that any particle path which crosses the x-axis twice must be
closed. Since, as a result of the symmetry, a particle that has travelled along the
curve P;P2 in Figure 4.6.2 must travel back to P; via the mirror image of PiP2.
Figure 4.6.2: Ina vortex, where the flow field has the form given by Equation (4.6.5),
any particle path given by us = 4+ Wo must be closed if it crosses the x-axis twice.
Therefore, the trapping capability of a vortex is not conditional on the flow being
rotational.
A fair model of many natural vortices is the Rankine vortex in which the
fluid velocity field is given by
ie Saal Mamet A
BO) = GRE ORY ie cael)
It is characteristic for this vortex model that the core rotates as a rigid body,
with the velocity being proportional to the distance from the centre, while
u_ becomes inversely proportional to this distance farther away.
In this vortex a sand grain with settling velocity wo can theoretically be at
rest at the two singular points in Figure 4.6.3. At these two points we have
u+ws = 0, and their coordinates are (R’w/2w, + V(R? @/ 2wo)” =e 0) The
circle which is shown, is given by the equation
z/R
ZZ 2 2
fe
= ae ra | = ey = (4.6.11)
R 2Wo R 2Wo
and it is the locus of all points with ws= 0. Sand grains will move
upward (ws >0) in the interior and downward (ws <0) outside the circle.
Figure 4.6.3 also shows the sediment particle paths corresponding to
Us = u+ws. Some of these are closed and could thus keep sand grains trapped.
Trapping is only possible if the maximum upward water velocity in the
vortex exceeds the sediment settling velocity which, for the Rankine vortex, is the
case when Wo < WR/2.
The equation of the sediment path through a point (Xo, 0) is
Now it could be argued that, since the trapping paths are closed, sediment
particles are no more likely to get onto them than to get off. Hence the trapping
mechanism might not be effective. However, the situation in practice is that the
sediment particles get into the vortex during its formation. This process is easy to
observe with the vortices behind ripples and dunes.
Us = uU + Wo (4.6.1)
This equation is only strictly valid however, if the flow is steady and uniform,
and that is never the case for vortex flow. Vortex flow is always accelerated, and
the accelerations will cause deviations from the closed sediment particle paths
described above. In the following we shall describe these acceleration effects and
show that they consist of a spiralling effect and a steady horizontal drift.
Qualitatively the spiralling effect is familiar and easy to understand. The
vortex is a centrifuge. Heavy particles will spiral outwards and light particles or air
bubbles will spiral inwards towards its orbit centre.
The steady horizontal drift is best understood by considering a sediment
particle at rest at the point (wo/q@, 0) in a forced vortex with angular velocity
@, see Figure 4.6.1, page 182.
Resting at this position is possible for the sediment particle under the
assumption (4.6.1). However, a fluid particle whose orbit passes through this point
is travelling on a circular path with radius ws/@ and radian frequency «. It is
therefore accelerated towards the vortex centre with acceleration @ wo.
The horizontal pressure gradient which provides the fluid particle with this
acceleration will tend to push the “resting” sediment particle towards the centre of
the vortex. This describes the driving mechanism for the steady horizontal drift
which is not only present at the point (Wo/@ ,0) but everywhere in the vortex.
In order to obtain a quantitative description of the drift, consider the sediment
particle in Figure 4.6.4 which moves around the circle
Ros ot + Wo/®@
= 4.6.1
rs) R sin ot ) kAi6:12)
in a forced vortex in accordance with the zero order solution us = u+Wo.
A water particle on the fluid orbit, which passes through the sediment
particles position, would be travelling on a circular path with radius rs(t) and
radian frequency @ and hence would have the acceleration
d
= = -w'rs (t) (4.6.13)
dulg i 2 » {cos ot ie Wo @
6.14
dt oF Sra 0 ) a )
(1)- of3)
Wey =z
an
and with this and Equation (4.6.14) inserted, the equation of motion (4.5.2) for a
small particle reads
du
<a + Cea — @wr = awRcos@t + (a-1)Wo @
dt Wo
dw. ag ;
—T + Wri + Our = aorR sin wt (4.6.16)
dt Wo
(4.6.17)
Uri) _ @WR_
2
(cos at 1
bes strat Die © (ises — (-a) v6 (0| (4.6.18)
or
2
ur) = a Wor (1) = ad = Q)Wo B (| (4.6.19)
which shows that the first term represents a centrifugal effect with a spiralling time
scale
=i
ip = Eal => (4.6.20)
while the second term represents a horizontal, steady drift towards the z axis.
4.7.1 Introduction
In order to understand sediment suspension in turbulent flows it is first
necessary to analyse the influence of turbulence on a single particle which is
settling through it.
Experiments by Murray (1970) indicate that the turbulence has two main
effects on settling particles, which may be quantitatively described in the following
way.
Consider an ensemble of many identical particles which are dropped into a
water filled jar where a certain level of turbulence is maintained by some kind of
stirring, see Figure 4.7.1.
The distances L; , which the particles settle during the interval ¢, are
measured and the corresponding time averaged velocities wj=Lj/t are
calculated. It may be observed then that the ensemble average wj;_ is less than the
still water settling velocity
Wi < Wo (4.7.1)
and that the variance Var{wj} is roughly proportional to: the turbulence intensity
w as to the Lagrangian integral scale of the turbulence Ty ; and to //t. That is
: 0
Wot
STILL WATER
z
Figure 4.7.1: Sediment particles settling through turbulent water will, on the average,
settle more slowly than they would through still water.
X
Figure 4.7.2: A sediment particle which settles through this arrangement of vortices
will experience a considerable delay compared with settling through still water
because it spends more time in the upward moving parts of the fluid than in the
downward moving parts. Note that the spatial average of the particle velocity is wo .
To quantify the effect, consider the special flow pattern in Figure 4.7.2. For
simplicity, let the flow structure be such, that the particle velocity varies with z as
where D is the vortex diameter and A <1. Then the resulting settling velocity is
Ws = = (4.7.4)
where T2p_is the time it takes to settle past two complete vortices
2D
r dz
4.7.5
Ce JWs(Z)
0
2»
dz
RPS JWo (1+A cos 1 z/D)
0
2D
Typ2D = ——a
wo dA
Hence, the settling velocity is reduced from its still water value wo to
Ws = Wo(1-A’)°> (4.7.6)
This expression becomes meaningless for A > 1 but in compliance with the
physical “reality” of Figure 4.7.2, we may extend it to
This corresponds to the fact that when the maximum upward water velocity
along the symmetry line in Figure 4.7.2 exceeds wo there will be positions where
u+Wo = 0 at which the particle will stop and could be trapped forever.
In practice, the time limit for trapping is set by the typical lifetime of local
flow structures Tg which is called the Eulerian time scale for the turbulence.
ad 4
£7?he = 2zw 4.7.
(4.7.8)
we get
t t
E{2(0} = E{2) w(t) |w(e) dr dr} (4.7.9)
oO oO
Pww(Ar) = +Ad)}/ow
E{w(to) w(to (4.7.10)
TL Th
Figure 4.7.3: Typical shape of the Lagrangian autocorrelation coefficient.
Measurements of Pww(Ay are quite scarce, but they tend to confirm the
general picture presented above, see e g Snyder & Lumley (1971) or Sato &
Yamamoto (1987). Also, the integral
Tr = Jpw at (4.7.12)
0
. =) \2/pw©) (4.7.13)
In terms of these parameters, two asymptotic approximations for the spread
given by Equation (4.7.11) can be derived, which are valid for relatively short and
relatively long times respectively
Since, for t1< tz, the inner integral in Equation (4.7.11) is approximately equal
to t, , while for large t; , itis approximately 7, .
Considering a heavy particle with settling velocity wo we may, as a first
approximation, use the approximation ws = w—wo for its vertical velocity. The
accuracy of this approximation was discussed in Section 4.5.
Using this simple superposition law, Taylor’s dispersion model can be used
to estimate the variance of the settling velocities of the particles settling through
the turbulence in Figure 4.7.1. If the settling time is large compared to the Integral
time scale, we can use the approximation (4.7.15) and get
Elsie Ege
Or} 8 = Tot (4.7.18)
In order to capture the loitering effect we must consider the velocity derivative in
its Eulerian form
(4.7.20)
for the unconditional expected value of the absolute velocity increment.
Wishing to model the loitering effect, the essence of which is that small
velocities correspond to slow changes of velocity, see Figure 4.7.2, we consider the
conditional expected value of [dy
(4.7.21)
The value of this expression will depend on the possible correlations between the
Me ow dw ow
velocity components (u,v,w) and the partial derivatives wer
ax oy’ Oz
However, knowing nothing about these correlations it seems coe ire as the
simplest option, to assume zero correlation so that
Ef \8wi jui-i} =
(4.7.22)
E(\SA}
ow
= >
ns Ow
(4.7.23)
and
Ow TE
AE = iF (4.7.26)
Based on Equation (4.7.25), we see that the replacement for Taylor’s constant
correlation
p = exp[-d/Tz) (4.7.17)
For a non-neutral particle with still water settling velocity wo which moves
in accordance with ws = w - Wo, the analogous expression is
The last term in this expression generates the loitering effect because it will be
small if wo is nearly balanced by w;-1 and hence j-; will then be closer to
unity.
The loitering effect expressed by Equation (4.7.28) is much weaker for three
dimensional flow than it would be for purely vertical flow. This is because of the
first two terms inside the square bracket whose randomness will smooth things out.
However, the loitering effect may still be significant.
Figure 4.7.4 shows results of a numerical simulation for a particle with still
water settling velocity wo settling through three-dimensional turbulence with
u’ =v" =w’ =o". The model is described by
0.8
0.6
0.4
® Wo = 1 cm/s
0.2 | X Wo= 2cm/s
© Wo = 3cm/s
A Wo = 4cm/s
0.1 02 04 1.0 2.0 4.0 10
F igure 4.7.4: Measured settling velocities from grid turbulence by Murray (1970)
and simulated values based on the random walk model described by Equations
(4.7.28), (4.7.29), and (4.7.30) with Ag = J.
Vi = Piven + OVI-pr Gy
and
The actual value of this correlation coefficient for various types of natural
turbulence is presently unknown.
This means that the validity of the assumption (4.7.22) is basically unknown.
If the correlation turned out to be zero for some flow, it would mean that there was
no settling delay due to loitering in this flow.
If the correlation turned out to be negative, a random turbulence field would
tend to enhance the settling velocity.
SEDIMENT SUSPENSIONS
5.1 INTRODUCTION
balance the settling rate cws. The vertical sediment velocity ws is practically
equal to the local, vertical fluid velocity minus the still water settling velocity,
Ws ~ W— Wo, See Section 4.5.5, page 178.
Stationary sediment suspensions do exist. Hence, the upward fluxes which
balance settling must be present. They occur in turbulent flows because water
parcels which travel upward through a given plane generally contain larger
sediment concentrations than water parcels which travel downward through the
same plane. In other words, the upward sediment flux required to balance the
settling rate in a steady situation, is generated by vertical mixing. This type of
process is called gradient diffusion.
Random vertical mixing can generate an upward sediment flux provided the
average sediment concentration decreases with height above the bed.
To obtain a quantitative description of gradient diffusion, consider the simple
mixing process which is illustrated by Figure 5.1.1.
VM MMI
NS BSS Ss
SUL Sd hoe
Lm process of gradient
diffusion may be
quantified in terms of
dz { IN QS \ ekJae exchange
Assume that the sediment concentration varies linearly with z and that
mixing is provided by exchanging layers of thickness dz and separation Jm once
every ¢ seconds.
Since the concentration difference between such a pair of exchanging layers
; dc
is lm aE each exchange corresponds to transport of the volume — dx dy dz lm dc
dz
upward through an area dx dy of an intermediate plane, say z = zo . The minus is
elie mo 92 ,
due to the fact that, if a positive, the exchange results in a net downward
: dew: ;
transport of sediment. In most natural cases, a negative and the exchange
process provides the upward transport required to balance settling.
Summation of the contributions from all exchanges which intersect the plane
Z = Zo, gives a total transport of — dx dy Im = through the area dx dy.
With the exchanges occurring at the rate of one every ft seconds, the vertical
transport rate per unit area is therefore
lm de
Epes rae 51.2)
dc
qz = —&5—
es (5.1.3)
el
corresponding to
2
Es = bm (5.1.4)
Thus, sediment diffusivity has the dimension of Lee like viscosity, and the
SI units are m/s.
The equation which expresses a time-averaged balance between settling and
the diffusive flux, given by Equation (5.1.3), reads
(5.1.6)
Gradient diffusion can be used to describe the upward sediment flux if the
mixing length Jm is small compared to the overall height of the concentration
profile. However, if the two lengths are of the same order of magnitude, a different
model must be applied. See the discussion in Section 5.4.2, page 233 for details.
Nevertheless, gradient diffusion has been used exclusively for the modelling
of suspended sediment distributions since the Nineteen Thirties without critical
assessment of its validity as a description of the actual physical process. Worrying
experimental data like those of Coleman (1970) and of Nielsen (1983) have
generally been bypassed by ad hoc modifications to the sediment diffusivity, see
Figure 5.1.2.
Coleman’s data show that, in order to model the distribution of suspended
sediment in steady, open channel flows as pure gradient diffusion, the sediment
diffusivity €s must generally be of a very different magnitude than the eddy
viscosity Vr, (Vt =~ Kuxz(1-z/p) for z<D/2). Furthermore, the magnitude of
€s must be a strongly increasing function of the relative settling velocity wo/ux.
/U.D
Es
° 5
e a
@
tc) oO
® B
4 8
A +
v x
Figure 5.1.2: Sediment diffusivities derived through Equation (5.1.6) from measured
concentration profiles c(z) under the assumption of pure gradient diffusion.
Different sediment sizes (different values of the relative settling velocity wo/u)
give very different values for the diffusivity €s. Data from Coleman (1970).
The data of Nielsen (1983) and McFetridge & Nielsen (1985) from
oscillatory flow over ripples (Figures 5.2.11 and 5.2.12, p 222-223) give further
evidence that pure gradient diffusion is inadequate as a model of suspended
sediment distributions. In order to explain their measurements, the magnitude of €s
would have to be different for different grain sizes, as for Coleman’s
measurements. More importantly, however, the distribution of €s(z) would have
to be significantly different for different sand sizes in the same flow.
This would be very unsatisfactory. Using gradient diffusion as a model
makes sense only if the same diffusivity can be used for all particles and if this
diffusivity is closely related to the diffusivity of momentum, i e, the eddy viscosity.
In recognition of the increasing weight of experimental evidence against pure
gradient diffusion as the sole distribution mechanism, Section 5.4 is devoted to a
new modelling framework which acknowledges the presence of large scale or
convective mixing mechanisms.
The new, combined convection-diffusion model is discussed qualitatively in
general terms and developed quantitatively for the special case of pure,
non-breaking waves over rippled beds. For that case, the new model leads to
improved agreement with observations, but the details of the model outlined in
Section 5.4.8 should not be seen as definitive. As more detailed data become
available, the details of the model should be adjusted as appropriate.
The main point at present is, that a combined convection-diffusion model
with the general features outlined in Section 5.4.5 has the potential of vastly
improved modelling of sediment suspensions, including aspects which could not
possibly be accounted for by pure gradient diffusion.
It is hoped that the combined convection-diffusion approach will make it
possible to reconcile the concepts of sediment diffusivity and eddy viscosity, so
that €s ~ v;. This is not possible under the assumption of pure gradient diffusion,
as illustrated by Coleman’s data in Figure 5.1.2.
Before the development of the new model, an illustrated introduction to
suspended sediment distributions in coastal areas is given in Section 5.2, page 206.
For flows that are both steady and uniform, it is reasonable to assume
instantaneous equilibrium between near-bed sediment concentrations and the flow
conditions (represented for example by the Shields parameter). However, this is
not possible if the flow is unsteady or non-uniform. In unsteady or non-uniform
flows, the near-bed concentrations are often considerable when the instantaneous
bed shear stress is zero because sand is arriving from above.
The alternative approach, which will be discussed in detail in Section 5.3
(page 222), is to consider entrainment and deposition separately and to try and
specify the pickup rate p(t) in terms of the instantaneous, local flow parameters. It
is acknowledged, however, that a complete description of the sediment pickup
process under waves may not become available for some time.
Bagnold’s (1956) definitions, may be partly bed-load and partly suspended load
near the bed. Thus, after measuring the total concentrations near the bed, the
problem of separating bed-load and suspended load remains.
As an example, this separation will be attempted with the total concentration
measurements by Horikawa et al (1982), using Bagnold’s definitions.
Bagnold’s definition of bed-load can be written as
co
where Ge(z) is the dispersive stress which supports the bed-load that is travelling
above the level z.
Solving this equation with respect to the bed-load concentration gives
-1 doe
(5.2.2)
sihealed <6 oe
Consider Case 1-1 of Horikawa et al (1982) at the phase 90 degrees which
corresponds to maximum free stream velocity. The general experimental
conditions are described by (A, 7, d,s) = (0.72m, 3.6s, 0.2mm, 2.65).
The linear sediment concentration A, derived from the total, volumetric
concentration in accordance with the definition (2.1.1) page 96, is of the order
unity, and the velocity gradient is of the order / 00s close to the bed. Hence, the
Bagnold number, defined by Equation (2.1.2), page 97, is about 20. Therefore, the
relevant formula for the dispersive stress according to Bagnold (1954) is
Between the outer, "exponential region" of pure suspension and the level of
the undisturbed bed, the concentration gradients are considerably greater. Through
this layer, which for these experiments has a thickness of five millimetres, the total
concentrations drop by a factor 100.
e Measured total
x Estimated bed-load
O Suspended
x
x
UNDISTURBED
In the moving layer below the level of the undisturbed bed, the
concentrations decrease rather slowly with elevation, from cmax in the stationary
bed to between 50% and 80% of cmax at the level of the undisturbed bed.
advantage of being simple and reliable and their output is easy to interpret. The
concentrations, determined from drying and weighing the captured samples, are
independent of such factors as grain size, shape, colour and grading. These factors
influence the outputs of other devices to varying extents.
When using suction samplers, consideration must be given to nozzle
orientation and intake velocity because these may influence the amount of sand
that gets captured. If the suction speed is low compared to the ambient flow and to
the sediment settling velocity, and if the grains have to turn sharp corners to get in,
the larger grains will tend to escape.
For steady flow the criterion is simple. The intake velocity should match the
flow velocity with respect to both direction and magnitude. This is of course not
possible in a wave motion so compromises must be made and corrections may be
required.
Bosman et al (1987) investigated the efficiency of suction samplers in
oscillatory flows. They found that for nozzles perpendicular to the main flow
oscillation the capture rate increased with increasing intake velocity up to
Uintake ~ A @. For higher intake velocities, the capture rate was fairly constant at
around 0.8.
In order to capture instantaneous suspended sediment concentrations, many
other types of concentration gauges have been developed over recent years. Most
of these rely on the absorption or back-scatter of sound (Jansen, 1978), or
electromagnetic radiation (Sleath 1982, Downing et al 1981). The output of these
types of instruments depends strongly on such sediment characteristics as grain
size, shape, grading and colour which vary with both location and height. Precise,
quantitative information is therefore difficult to derive. Nevertheless, even
qualitative information such as correlations between velocities and concentrations
is valuable at the present state of the art.
The complications mentioned above are less severe for devices which are
based on X-ray or y-ray absorption. The absorption of these types of radiation
depends almost exclusively on the density of the absorbing medium.
Also the dependence of the electric conductivity on sediment concentration
has been used to measure sediment concentrations, see e g Horikawa et al 1982.
detailed descriptions.
Theoretical considerations often deal with the horizontally averaged quantity
c(z,t). It is, however, important to acknowledge the difference between this
theoretical quantity and conveniently measurable quantities such as the average
concentration along a sensor beam or over a small measuring volume. While the
time series of the real horizontal average c(z,t) may be fairly well behaved, the
time series of the measured quantities may look very noisy due to the "spotted
carpet effect". That is, since natural sediment concentrations are horizontally
non-uniform, the "carpet of concentrations" at a certain level c(x,y,zo,t) is spotted
and these spots are convected past. the sensor. This may result in much high
frequency variation which is irrelevant to the modelling of the horizontal average
c(z,t).
Figure 5.2.2 shows suspended sediment concentrations under irregular waves
measured with y-absorption along a 30cm shore parallel path. The bed was
covered by rounded megaripples which do not shed vortices regularly as vortex
ripples do. The measured concentrations seem random to a great extent although
some correlation obviously exists with the measured near-bed velocities.
u (0.20,t) [m/s]
c (0.04,t) BY VOLUME
0 1 2 3 4 5 6an
MINUTES
LER ENNEN eee ceemeeeeeeee eS eS
0 102030 SECONDS
Sometimes, a certain area of a sand bed with megaripples will generate large
suspended sediment concentrations for a given set of flow conditions. At other
times, it may not.
Thus, individual entrainment events are not entirely predictable in terms of
the free stream velocity history. They result from complicated flow instabilities
which also depend on subtle details of the micro-topography.
Nadaoka et al (1988) described an important mechanism for the creation of
large clouds of suspended sediment in the outer surf zone, where the bed is
commonly covered by megaripples. They attributed the formation of these clouds
to obliquely descending vortices which are able to lift considerable amounts of
sand into suspension after they reach the bed. Some of the large concentration
peaks in Figure 5.2.2 may be associated with this type of sediment clouds.
The concentration variations in Figure 5.2.2 result not only from clouds of
sand rising from the bed or settling from above, but also from clouds being
convected horizontally past the probe.
More details about concentration time series under irregular waves are given
by, for example, Hanes & Huntley (1986) and Vincent & Green (1990).
The picture is somewhat less confused if the wave motion is regular,
especially under sheet-flow where the conditions are essentially constant in a
— @—
C+C —— _z=1.5mm
<<a -- z=14mm
—-— z=26mm
2x10°
This is because the sand clouds remain identifiable for a considerable part of
a wave period. Hence, other sand clouds than the latest one generated at the nearest
ripple are also making an organised imprint. In some cases a sensor may also see
the same sediment cloud twice.
When the sand bed has high-profile features like vortex ripples, the
suspended sediment concentrations cannot be expected to be horizontally uniform
and the velocity field will likewise be non-uniform in a horizontal plane.
Consequently, the time-averaged vertical flux of sediment wsc_ will vary
greatly in a horizontal plane. In fact, the upward transport of sediment over vortex
ripples occurs through quite narrow spatial corridors which originate near the
ripple crest and follow the typical upward paths of the released lee vortices, see
Figure 5.2.5.
The existence of special corridors for upward sediment movement is very
obvious from the observations of Nakato et al (1977). They found that the time-
averaged, vertical sediment transport rate cCws was clearly different from zero
almost everywhere in the two vertical sections which they studied. One above the
& L
<<
<—
~<—
<—
<<
>
«ee
<—
Figure 5.2.5: The upward transport of sand over vortex ripples is very organised.
The sand is carried upwards by the lee vortices which are released at the free stream
velocity reversals. The arrows in this figure indicate the local, time-averaged, vertical
sediment flux c Ws.
ripple crest and one above the trough. In general, they found that wsc was
negative everywhere in their two vertical sections except very close to the ripple
crest, see Figure 5.2.5. At the same time, the stationarity of the situation obviously
requires zero net vertical transport on the average (averaged over time and over
any horizontal plane).
Nakato et al (1977) estimated the average upward sediment flux as
(w-—wo)¢+wc in accordance with the assumption that the sediment velocity
equals the fluid velocity plus the still water settling velocity of the sediment
Us =U +Wo. The limitations of this assumption are insignificant for the present
purpose, see Sections 4.5 and 4.6 in general and Equation (4.5.26) page 176 for a
simple estimate.
ie)
nm
ELEVATION
CENTIMETRES
IN
4 6 10-5 2 4 6
Figure 5.2.6: Vertical distributions of the time-averaged concentration ¢ and the
root-mean-square of the periodic Crms. The values are obtained by averaging the
values given by Nakato et al (1977) for two vertical sections, one over the ripple crest
and one over the trough.
10
ae:
gv
aot
3 7(3)
E
S6
5
5
3
bea
5
ue
jam
2
1
iene Kt ee (5.2.5)
35
30
25
20
15
elevation
Local
[cm]
8.01 0.1 1 .
Averaged concentration [g/l]
The main reason for the difference in concentration profiles is the difference
in ripple size. Thus, for the two cases shown in Figure 5.2.9, the regular flow
formed ripples with an average height of /5cm, while the ripples in the irregular
flow had an average height of only 4cm.
As a consequence of the bed roughness (~ ripple height) being smaller, the
shape parameter wo T/NArr will be larger for the irregular waves and the profile
shape therefore tends to be more upward concave.
20
*)
160.°2 4- 10,529 “4 10°2 4 10°2 4
c(m*/m?)
Figure 5.2.10: Time-averaged concentration profiles under non-breaking (x) waves
and under spilling breakers (0) of the same height. Wave flume data from Nielsen
(1979).
@ d=0.42mm - 0.60mm
+ d=0.06mm - 0.11mm
ELEVATION
ABOVE
(cm)
RIPPLE
CREST,
z
fe)
10° 10° 1007
CONCENTRATION BY VOLUME
Figure 5.2.1]: Concentration profiles for two sieve fractions of sand suspended in
the same flow (waves over natural ripples). The profiles are not similar. The trend
for the fine sand is convex upward while the trend for the coarse sand is concave
upward. (h, T, H, dso, A, nN) = (0.3m, 1.51s, 0.13m, 0.19mm, 0.078m, 0.01 1m).
Z(cm)
SQ
-42-0.6 0.3-0.42 0.18-0.3 0.15-0.18 0.11-0.15 0.06-0.11
Different grain sizes are not picked up in direct proportion to their abundance
in the undisturbed bed. Small particles have a greater chance of becoming
suspended than large particles. Details of this phenomenon are discussed in Section
5.3.7, page 230.
c (x,y,z,t)
> —*
q, q,+ al dx
d (x,y,t) p (x,y,t)
ae a Sn a ee eRe ah ie cer
eee eee
CB OM YS
OeSeerie:
Sree a am Se) eer isoo
6) eecayenne
fel ew a tae UG:
eer Meer
Meee ee eee ater
OS CC6 feria Gretta:
RAC SEIS oie at he PAC OTee OSES ao
Sere
Oo
eenee ete
Shel te
SY ehOLaKaTe te Wh 6 ale te gere wiinlets OE bk Le
Shh 75 5
Sire ‘ea etas
(rsSko
STO gue
BawahaS pele ve ter gtae
tererscatc
Ouin terol lame,
Figure 5.3.1: The imbalances between local pickup rate p(x,y,t) and deposition rate
d(x,y,t) = Wo C(x,y,zr,t) through the near bed reference level z, generate the changes
in bed level.
)
Sy bot d-p
hee. no. (abhlG of0ofBen+ Gay
oq
Ra I(+ y |e (5.3.1)
a)
COO :
SSCS nar ae sity tte ®
Figure 5.3.2: Van Rijn (1984) studied pickup functions in steady flow via the scour
rates and the build up of the concentration profile over a sand bed following a fixed
concrete bed.
P = woCo (5.3.3)
A considerable amount of information is available about Co [=c(o)] under
waves, see e g Nielsen (1986).
Alternatively, we can try to derive an expression for p(t) from van Rijn’s
formula (5.3.2) (or a similar steady-flow formula) if we assume that this
steady-flow formula can be applied instantaneously to an unsteady flow. With a
slight rearrangement van Rijn’s formula reads:
In order to apply this formula it is, of course, necessary to know how the
effective Shields parameter 6’(t) varies with time. Attempts to estimate p(t) for
three special cases are presented in the following sections.
and to apply the general value @,~0.05 for the critical Shields parameter. This
model for 6’(f) was found to work well for estimation of bed-load transport rates
over flat beds in Section 2.4.4, page 121.
Hence, we shall try and insert these expressions for 9’(f) and Qc into van
Rijn’s pickup formula (5.3.4) and compare the result with experimental data.
For the purpose of comparing van Rijn’s formula (5.3.4) with p-values
determined from concentration measurements via Equation (5.3.3), we note that
the last factor in Equation (5.3.4) is almost proportional to the settling velocity
within the relevant range of quartz grain sizes. Thus, the expression
na (ae ge3
Wo = yo
is within twenty percent of Gibbs et al’s (1971) expression (4.2.6) for quartz
spheres in the diameter range 0.3mm < d < 1.5mm.
Inserting these approximations for 6’(f) and wo into van Rijn’s formula
(5.3.4) we find
DP ~
uv 0.007 wo (025-0.05)!° (5.3.7)
Combined with Equation (5.3.3), this formula corresponds to
O'(tn) =
V4Coed
fr.s Arms (COSt Uco(tn+1) — Uco(tn—1)
Wp Ucoo(tn) + sinGr as ay gaat
(2.4.15)
on the basis of measured or simulated values of the free stream velocity Ueo(t). Wp
is the spectral peak angular frequency for the velocity record and
Arms = V2 (Uco)rms/Wp. The grain roughness friction factor 25 is calculated
from Equation (2.2.6), page 105, and @r is the assumed phase shift of the bed
shear stress ahead of the free stream velocity at the peak frequency. Details are
given in Section 2.4.4, page 121.
Co = 0.005 62 (5.3.10)
for the near-bed reference concentration, which through Equation (5.3.3)
corresponds to
D = 0.005 wo @ (5.3.11)
These formulae for Co and p apply reasonably well for both rippled and flat
beds, see Nielsen 1986. The modified effective Shields parameter 6, is given by
es 025
= ——-_ 53ah2
(l—- n/a)? ee
where the grain roughness Shields parameter 625 is calculated as described in
Section 2.2.5, page 105. The correction factor (1 —m N/A)’ is the square of the
velocity correction suggested by Du Toit & Sleath (1981) for the flow
enhancement near the crest of vortex ripples.
For irregular waves, Equations (5.3.10) and (5.3.11) may be used with the
regular wave height H_ replaced by Hrms_ or, correspondingly with Aw
replaced by V2 (ucc)rms in the calculation of 025.
With respect to the time dependence of p(t) over rippled beds, Nielsen
(1988) noted that most of the resulting sediment transport over vortex ripples
occurs above the ripple crest level and that the input of sand into this domain
happens virtually as instantaneous puffs at the time of free stream velocity
reversals. Thus, the pickup function must vary with time as shown qualitatively in
Figure 5.3.4. It has two distinct peaks at the times (7 and ¢“ where the free stream
velocity changes direction. These are the times when the lee vortices with their
clouds of sand move upwards into the main flow.
Quantitatively, these pickup functions may be described in terms of functions
of the form
u(t)
p(t)
E
= ~
Figure 5.3.4: The sediment pickup function p(t) under regular waves has two peaks
per wave period. For rippled beds the peaks are very sharp and correspond to the
release of the lee vortices. For flat beds (sheet-flow) the peaks are flatter and occur
closer to the extreme velocity in either direction, see Figures 5.2.3 and 5.2.4.
The periodic delta function 67(t-to) has the dimension T!. It is periodic with
period T and is further specified by
The coefficients V4 and V in Equations (5.3.13) and (5.3.14) are then seen
to represent the total amount of sand (per unit area) which is picked up in each
pickup event. They have the dimension of length.
The sum of V@ and V“ must therefore, if the process is stationary, balance
the total amount of sand which settles out in one wave period , ie
Vt neuawveiCo Tt (5.3.17)
Their relative magnitude V/V" must reflect the relative ability of the forward and
the backward velocities to entrain sand.
6) 1 2 3 4 5 6
e : w,lcom/s],x :dx10f[mm]
Figure 5.3.5: Variation of the average grain size and the average settling velocity
with elevation above the ripple crest. The abrupt decrease in both quantities from z
=0 to z=0.01m indicates that most of the sorting happens in the pickup process.
After that, the size distribution remains fairly constant. Data from Nielsen (1983).
Noting that the Shields parameter, see Equation (2.2.2), p 103, is proportional
tod” fora given bed shear stress. One might expect the relative pickup rate for
sand with size d to be proportional to dso/d. It might also be argued that, the
pickup rate of a sand fraction with settling velocity w should be proportional to
u,/w . Hence, the pickup rate relative to that of the median size should be
proportional to ws0/w. On the other hand, small particles are sometimes sheltered
against the eroding effect of the flow by larger grains (the armouring effect), and
this makes things more complicated.
There is precious little information available about the selective entrainment
of different sand sizes under waves. However, the data in Figure 5.3.6 indicate that
the simple formula
fraction in bed d
fraction in near—bed suspension ~ dso (5.3.18)
can be used as a rule of thumb.
6 % in bed
in nearbed suspension
Equation (5.3.18) corresponds to the time-averaged pickup rate for sand with
parameters (d,w) being given by
= d "
D(dw) = i * fraction in bed * p(dso) (5.3.19)
where p(dso) is given by a suitable formula for the bulk pickup rate. For example
Equation (5.3.11), page 228.
5.4.1 Introduction
In the previous section the processes of sediment pickup and deposition were
considered, i e the processes by which sediment goes from a state of rest to a state
of movement and vice versa. The present section deals with the processes which
take the sand upwards from the immediate vicinity of the bed.
These processes are under one referred to as distribution processes but two
distinct categories aie defined and aiialysed in detail. These are convective
processes and gradient diffusion.
A quantitative framework is suggested for describing convective entrainment.
The characteristic behaviour of sediment concentrations c(z,t) resulting from this
model are compared to that resulting from the gradient diffusion description.
It is shown that natural entrainment processes, in general, contain elements of
both convection and diffusion. Therefore, a combined convection-diffusion
description for the distribution of suspended sediment is developed in Sections
5.4.5 through 5.4.9.
CONVECTION DIFFUSION
the overall scale, the process may be described as diffusive, see Figure 5.4.1.
For both types of processes, and indeed for combined convection-diffusion
processes, the conventional means of describing concentration profiles of
suspended sediment is the conservation equation for the volume of sediment which
may be written as
which expresses that a divergence of the sediment flux field g = c us must result
in a change of the local sediment concentration.
For the sake of simplicity, we shall in the following, consider a horizontally
uniform sediment concentration field c = c(z,t) and a correspondingly uniform
sediment velocity field us = us(z,t). The conservation equation (5.4.1) is then
simplified to
Oc Ne Nias Cl pea
Be 1Ta inigyll WS) tara) RAD
where qz denotes the total sediment flux (transport rate per unit area) in the
z-direction.
It is generally considered too complicated to describe the vertical sediment
velocity ws in detail in a turbulent flow. So instead, a broader approach is
generally applied to the sediment flux. In a horizontally averaged description, the
vertical sediment flux gz is considered to consist of a downward component -woc
due to gravitational settling, and an upward flux which can be of convective
(subscript C) or diffusive (subscript D) nature or a combination of both. The total,
vertical sediment flux is thus written as
q, = ~Wolt+Gnt+ |e (5.4.3)
dc _ yy 9C _ Op _ 94e (5.4.4)
see Figure 5.4.2.
‘c i
fe}°
Qo = —€sVec (5.4.5)
ie, the diffusive flux vector is directed from larger towards smaller concentrations.
It is proportional to the concentration gradient vector Vc and to the sediment
diffusivity ¢€s. The diffusivity has dimension length squared per unit time (or
velocity times length), the same as kinematic viscosity and eddy viscosity which
can both be considered as diffusivities of momentum. The diffusivity may be
expressed in terms of a typical mixing length as outlined in Section 5.1.2, p 201.
The concept of gradient diffusion was developed in connection with the
kinetic theory of gasses and statistical mechanics where macroscopic distances are
covered in many random steps between which the velocity changes due to
collisions. In that case the diffusivity has a clear physical meaning. It is simply the
typical step length (the mixing length) squared divided by the typical time between
collisions. Thus, diffusion is a useful concept when the motion considered is of this
random walk nature.
However, gradient diffusion cannot describe details of a mixing process on a
smaller scale than the mixing length. Processes where the mixing length is of the
same magnitude as the overall scale of the sediment distribution can therefore not
be described in terms of gradient diffusion.
Examples of such convective processes are the entrainment of sediment
from rippled sand beds under waves by travelling vortices, and the lifting of sand
straight from the bed to the surface by the rising plumes generated behind plunging
breakers. In these processes, the sediment flux is obviously not necessarily related
to the concentration gradient as expressed by the diffusion equation (5.4.5).
For the purpose of modelling natural sediment suspension processes, which
generally contain elements of both diffusion and convection, we must first develop
a quantitative description of the convective sediment flux. To this end, the vertical
convective sediment flux is written in the form
to
Figure 5.4.3: In the simple convection model considered here, a sediment particle
which is picked up from the bed at time ¢o will travel upwards with speed we until it
reaches its entrainment level ze at time to + ze/wc. After that, it is assumed to settle
with its still water settling velocity wo.
5
We = Wt = ee G23)
where the minus is due to the ">" in the definition (5.4.7) of F(z), which is
opposite to the standard definition of distribution functions in probability texts.
a, SE (5.4.9)
and with the convective flux described in accordance with Equation (5.4.6) and
Figure 5.4.3. With the expression (5.4.6) inserted for g * Equation (5.4.9) becomes
3 = Mog ty oc ey
P't-| FO - pt-F FO
? a | ae ’
(5.4.10)
: d
where the brief notation p’ = oh and F’ = has been used for the derivatives.
dc se
Wo = PF) (5.4.11)
and is easily integrated to yield
The constant of integration is found by considering this equation at the bed level,
Z=0
where woCo is the time-averaged rate_of deposition , see Figure 5.3.1, page 222.
Hence, for a stationary situation with dzp/dt = 0, we get
Wo C(0) = p (5.3.3)
That is, the average pickup rate p must balance the average settling rate at the
bed.
For such a stationary situation, Equation (5.4.12) takes the form
(5.4.19)
This equation is linear with constant coefficients. It can therefore be solved
analytically for most realistic functions F(z), but it has particularly simple
solutions if the distribution function F(z) is an exponential
F@) = e7 (5.4.20)
In this case the solution is
Bn = Pues (5.4.23)
We
The fact that Cyn is complex with a small negative argument shows that the
concentration cy(o,f) at the bed lags behind the pickup function py, . The
imaginary part of By, makes this phase lag grow with distance from the bed. In
terms of real-valued functions, the solution (5.4.21)-(5.4.23) can be written
This solution, which corresponds to the pure convection process, differs from
the pure gradient diffusion solution, which will be discussed below, in two
important respects. Firstly, the magnitude of all harmonic concentration
components decay as e © since Re{Bn} = J for all m. Secondly, the time lag
zlwe-Arg{Cn}/n @® grows at the same rate with z for all frequencies.
The expression (5.4.24) shows that, when the distribution process is
convective, a concentration peak will travel upwards with the convection speed we
ac
ae
_ Mee
dc _ 9p
5 (5.4.25)
Ip = ~8555 (5.4.26)
we get
Cer
or =
Ee 27.
Wo a7 ap -e (Es a
OC (5.4.27)
Wo Dcuinds fice.
am SF ay (Es 7) = 0 (5.4.28)
since the constant of integration is zero when there is no sediment flux at infinity.
The time-averaged Equation (5.4.29) can also be written as
feine (5.4.30)
dz Es
Bea
C(z) = Clo)e”?Sés (5.4.31)
The constant c(o) ie, the concentration at the bed, may be given in terms of
the time-averaged bed shear stress if the flow is uniform.__
Alternatively, if the bed-level is stationary dep _ 0), c(o) can be
dt —
expressed in terms of the time averaged pickup rate p , since the deposition rate
Wo C(o) must be balanced by p. See Figure 5.3.1. This leads to the expression
Z
= dz
Cz) = a e-wo] (5.4.32)
dc
—Es 7 pm p(t) for z=0 (5.4.33)
lee tre
— 3_
Es Es
in = 0 (5.4.36)
which corresponds to Equation (5.4.19), page 240, for the pure convection
problem. In order to get an impression of similarities and differences between the
two descriptions, consider the special case of constant diffusivity. For constant
diffusivity, Equation (5.4.36) is reduced to
With the boundary condition (5.4.33) and the simple harmonic pickup function
Pn —OnZWo/Es pinwt
cn(z,b) ———————
Woh é tena 2 (5.4.38)
where
cniéceugibl
2 4
wlan.
INWE,
Wo
(5.4.39)
(5.4.40)
Im (An)
ao,
ASYMPTOTE
Figure 5.4.5: The variation of the complex coefficient o» of the pure diffusion
solution. It differs from the corresponding B, (Equation (5.4.23)) of the pure
convection solution in that Re{a,} is a function of @ and, Im{an} is not
proportional to n®. This means that the different harmonic components of the
diffusion solution decay at different rates with z, and their time lags grow at
different rates with z.
from which we see that, in contrast with the convection solution, the rate of decay
with elevation of individual Fourier components cp(z,t) increases with frequency.
In the pure convection solution (5.4.24), all Fourier components of c(z,t) decay at
the same rate.
Furthermore, the time lags of cp(z,t) relative to the pickup function grow at
different rates with z. In the pure convection solution, all the time lags grow at the
same rate, namely, z/wc. This means that in a convection-dominated process,
concentration peaks will show up with very similar shape in records from different
elevations. In a diffusion-dominated process, the peaks will become blurred more
rapidly.
The length scale L of the convection solution is replaced by €s/wo in the
pure gradient diffusion solution.
scale turbulence is present. The role of the small scale turbulence is to smooth
concentration differences by gradient diffusion.
The combined nature of the entrainment process is also evident from
measured time-averaged concentration profiles of different grain sizes in the same
flow. As discussed in Section 5.2.3, page 209, the distributions of different sand
sizes are not always similar, as they should be in any one of the pure processes, see
Equation (5.4.15) page 239 and Figure 5.4.4 page 242.
For a combined convection-diffusion process, the behaviour of the
horizontally averaged concentrations c(z,t) can be described by the conservation
equation (5.4.4) page 235 with both convective and diffusive fluxes included and
given by Equations (5.4.6) page 238, and Equation (5.4.26) page 241 respectively
dc Or ee S + oa “a vEt 0G
Finlay iis we? errr z era ah dz (Esey
(5.4.41)
By taking time-averages, this equation is reduced to
WOaide)
Pe ileal ee cies (5.4.42)
which describes the time-averaged suspended sediment concentrations in a
horizontally uniform, combined convection-diffusion process. This equation can
immediately be integrated once with respect to z
au = dc
WoC — pF(z) + es on = const (5.4.43)
0z
Wo (0,1) — p(t) Fo) = nSa
from the pure convection case, and
Wo C(0,t) + Es oC =n oe
dc
Wo C(0,t) — p(t) F(o) + esa, =p =
2b
LEW
The second and the third terms on the left-hand side are minus the convective
and diffusive fluxes upward from the bed. Together, these terms amount to the
total upward sediment flux from the bed, i e , to the pickup rate
C]
p(t) F(o) - €s Ee = p(t) (5.4.45)
or
~ e541, = [1-F
Thus, for a given pickup rate p(t), the concentration solution depends on the
chosen value F(0) of the convective distribution function at the bed. In the
following we shall mostly be using F(0)= 1 which through Equation (5.4.46)
corresponds
~€5 a es (5.4.47)
That is, the diffusive sediment flux is assumed to vanish at the bed. If the
diffusivity €s tends smoothly towards zero at the bed, then the diffusive flux
vanishes smoothly. However, it vanishes abruptly if €s is not tending towards zero
at the bed. Thus, the general validity of the choice of F(0) = 0 is perhaps
debatable. It is a convenient choice however, and it leads to reasonable agreement
with measurements, see Example 5.4.2, page 256.
For the time-averaged concentration we always have the simple bottom
boundary condition
WoC(o)
-p =n Ai
see Figure 5.3.1, page 222 which for a stationary situation with no net erosion or
deposition becomes
Wo C(0) - p = 0 (5.3.3)
yak
Ga = feWooe (5.4.50)
as long as the same function is used both inside and outside the brackets.
The first term inside the brackets of Equation (5.4.49) grows very rapidly
with the settling velocity wo. Thus, while the pure diffusion solution, which
corresponds to the last term, may dominate for the finer sand fractions, the first
term will become dominant for the coarser sand fractions. The first term of the
solution accounts for the convective entrainment, and for vanishing diffusivity, it
tends towards the pure convection solution
This may not be immediately obvious from the form of the solution (5.4.49).
It will however become evident from the Examples 5.4.1 and 5.4.3.
The fact that the shape of the convective distribution function F(z) can be
expected to be displayed by the coarsest sand through Equation (5.4.15) is
important. It shows that information about F(z) can be extracted directly from
measured concentration distributions of the coarsest sand fractions.
Fa) = et
and where the diffusivity €s is a constant. In this case, the integral (5.4.50) may be
taken as
G@ = aa
2rye D 1
gee A = sae liu = 1 ned —Wo2/E
Ea (5.4.51)
tee ESSE
WoL WoL
te F(z) (5.4.15)
Co
for vanishing relative diffusivity, €s/wol — 0
elevation
relative
(z/L)
.01 0.1 4
relative concentration
Figure 5.4.6: Relative time-averaged sediment concentrations corresponding to the
combined convection-diffusion solution (5.4.49) for different values of the relative
diffusivity €s/woL. For small values of €s/wol the solution approaches the pure
convection solution (5.4.15).
Bn = ieee (5.4.23)
We
(5.4.57)
The factor 4 corresponds to the factor of approximately 4 between the
diffusivities of time-averaged momentum (u) and of oscillatory momentum (7)
in a wave-dominated boundary layer. See Equation (1.2.38) and Figure 1.5.6, page
70. The eddy viscosity expression applied in Equation (5.4.57) is Equation
(1.3.14), page 52. This eddy viscosity estimate should be a reasonable choice for
the large relative roughness values (0.085<r/A<1.15) which are generally
exhibited by sand beds under waves. See Sections 1.3.3, p 42, and 3.6.4, p 152.
The vertical speed of sediment convection we should be closely related to
the speed w; with which boundary layer turbulence propagates upwards in
oscillatory boundary layers. Sleath (1987) found, from purely oscillatory flow that
w 5.05
Wr = 2919)
(1.2.3)
where 605 is the boundary layer thickness which corresponds to a dimensionless
velocity defect IDI of 0.05 at the top of the boundary layer.
Assuming the simple boundary layer structure which corresponds to constant
eddy viscosity, the defect magnitude |DI decays exponentially in accordance with
Equation (1.3.7). Then Equation (1.2.3), for the convection velocity can be written
edre ] (5.4.59)
c(o) Es
Hence, considering the ¢-distribution for the coarsest sand in Figure 5.2.12
we see that, in this semi-logarithmic presentation, F(z) has an upward concave
shape which is typical of functions of the form (J +z/L)" where L is an
appropriate vertical scale. Furthermore, F(z) is seen to decrease by roughly a
factor 10° through nine centimetres of height for those experimental conditions.
It seems natural to chose L = z1, where z; is the equivalent of the laminar
Stokes length in an oscillatory boundary layer with constant eddy viscosity, see
Equation (1.3.7), page 46. For the experiments reported in Figure 5.2.12, page 221,
Z, was approximately three millimetres. Hence, the elevation of 9cm corresponds
to z/L = z/z, ~ 30 anda power of -2 _ is thus required to give the observed
concentration decrease of the order 10°? over 9cm. Hence we may suggest
=n)C(z) = P_
ae -Gle) |©
| _Wo_ F(C)e"’ Go) do + 1 (5.4.49)
o dz
see page 248, where G(z) = lee and the reference level z; is generally
S
chosen to be the bed level, ie, zy = 0. The elevation z is measured from the
ripple crest level if the bed is rippled. The time-averaged pickup function is
predicted by
F(z) = 5.4.60
oe (1+2/z1) is
see page 254. The optimal choice of the bed roughness _r seems to be
or, with the dimensionless elevation § = C/z; and the dimensionless settling
velocity w* = WoZ1/E€s
+ 1) (5.4.62)
d(mm) w(ms)
—-— 0-125-0-18
seseseees O18 -0-25
0-25 -0-35
0-35 -0-50
IO LOO A
—-— 0:71-1-0 0-106
Figure 5.4.7; Time-averaged concentration profiles for different sand size fractions
over vortex ripples under non-breaking waves. Field data from Nielsen (1983).
elevation
Relative
(Z/Z1)
0
0.001 0.01
Relative concentration
Figure 5.4.8: Relative sediment concentrations calculated from the combined
convection-diffusion model, Equation (5.4.61). The numbers on the curves
correspond to the dimensionless settling velocity w* = woz1/€s.
It can be seen, in qualitative agreement with the data in Figure 5.4.7, that the
concentration profiles for the finer sand are upward convex while those for the
coarser sand are upward concave.
In detail, the calculations behind the c(z)-profiles in Figure 5.4.8 are as
r = 0.12m
C(Z,Wo) = Co e Wo2/Es
which would plot as a straight line in Figure 5.4.8. The deviations from this
straight line, very close to the bed, are due to the form of the bottom boundary
condition (5.4.47) applied in the combined convection-diffusion model.
In order to obtain upward convex concentration distributions for the fine sand
fractions, it is necessary to assume a diffusivity distribution for which the pure
diffusion solution (5.4.32) is upward convex. This is obtained when €,(z) is a
decreasing function of z. There are good physical reasons for suggesting that the
sediment diffusivity is a decreasing function of z over ripples under non-breaking
waves, but this fine tuning of the model will not be attempted at the present stage.
At present, the main point is, that the combined convection-diffusion model
is capable of explaining the qualitative difference between the distributions of
different grain sizes suspended in the same flow. That is not possible within the
framework of either pure gradient diffusion or pure convection.
The distribution of total concentrations (the "sum" of the size fraction profiles
discussed above) over a bed of fairly well graded sediment is not necessarily itself
given by Equation (5.4.61) with any particular value of w*. Nevertheless, it is
often attempted to predict the distribution of total concentrations with a formula
like Equation (5.4.61) and using the mean (or median) settling velocity of the bed
sediment.
For the case corresponding to the data in Figure 5.4.7, that leads to the total
c(z)-profile which is shown in Figure 5.4.9. The reference concentration
Co = c(O) = p/wo is calculated from Equation (5.3.10), page 228.
0.4
© DS)
flow and sediment parameters were (T, AQ, dso, Wo) = (85, 1.50m/s, 0.21mm,
0.033m/s ) corresponding to a grain roughness Shields parameter of 025 = 2.43
>> 1. Hence, the bedforms which were present (n, A) = (0.17m, 2.5m) would
have been rounded megaripples rather than sharp crested vortex ripples.
The boundary layer parameters and the concentration profile therefore are
calculated as for a flat bed. Then we find the hydraulic roughness r = 0.0026m
from Equation (3.6.14), page 159, leading to z; = 0.09VrA = 0.0063m, and the
diffusivity es; = 2@zi = 8.2:10°m’/s, see Equation (5.4.57), p 252. The
reference concentration is found from Equation (5.3.10), p 228
0.3
0 : ; ;
1E-05 0.0001 0.001 0.01 0.1
time averaged concentration
Figure 6.1.1: The modelling of coastal sediment transport processes still presents many
interesting, unresolved problems. The challenges of the swash zone are particularly daunting.
6.1 INTRODUCTION
models. The sediment transport rate, through a unit width of a vertical plane
perpendicular to the x,u-direction, can be calculated as
D
where c(z,t) is the local, instantaneous sediment concentration and us(z,t) is the
instantaneous, horizontal sediment velocity. In most cases we are mainly interested
in the time-averaged transport rate, Q which, for a periodic process can be
calculated as
t+TD
Q= =) Jcz) u(z,t) dz dt (6.2.2)
tz=o
where T is the wave period. In terms of the steady, periodic and random
components, this may be written as
O = |\(cu+cut+cw)az (6.2.3)
Z=0
see Equation (1.1.17), page 11.
The last term c’w’ in this:expression is generally expected to be small. The
first two terms may be of similar magnitude or either of them may be dominant,
when transport in the wave direction is considered.
Transport perpendicular to the wave direction will, of course, be dominated
by cv_ which is the equivalent of the first term in Equation (6.2.3).
By far the majority of the existing sediment transport models are based on
this "cu-integral approach". However, the alternative models, using the "particle
trajectory approach” will, in some cases, lead to simpler and more accurate models.
The basic idea of the particle trajectory approach is illustrated in Figure
6.2.1. It can be seen that the time-averaged transport rate Q can be expressed
very simply in terms of the time-averaged pickup rate p , and the average distance
I, travelled by sediment particles.
In a short time interval 6; around f! , an amount p(t’) 5: of sediment is
picked up per unit area. During the time interval from this pickup to the settling of
the last grain, these grains will make a total transport contribution of J p(t) 5r.
This corresponds to the time-averaged transport rate
ror:
2 ES ES GCE JRaCR SS PEDO BOPEORe) wo Gina TRO ia SSR ORe ans ORIORaAT icmkaealn Kaas
JSS TAT
eiran ei e ESE
ees nae oh figeMirRPO
a es Sig Mi PRT
MgriAl vwERT sg,Sr aan
PTR ROB
g YS iaiea 18g
Mel PT RT
eo Ye
os” IPT
>
eiSig,IP Pe PS > IWael
Oni e recog! TaharMeee
RP PeSig Bieta
Figure 6.2.1: Ina short time interval 8; the amount of sand picked up per unit area
is p. If the average distance travelled by the moving sediment is J/,, the
corresponding sediment transport through one unit width of the plane A is seen to be
lx,Pp &: ‘
Q = kp (6.2.4)
Dp = 8(0’-6)'> (2.3.11)
adAte
ae = ) (Ee
= 00003ieee z | 32
(5.3.2)
Oz = 8(0’-6)'!>
Vs-l gd d
where the definition (2.3.10), page 112, for ®g has been used, and
el et Pie
p = 000039") ee | V(s-l) gd
pres seope na
p
Thus, for sand which is being picked up at the rate given by van Rijn’s
pickup formula to yield the bed-load transport rate given by the Meyer-Peter
formula, the grains must, on the average, jump a distance equivalent to 44/ grain
diameters in each jump. Interestingly, this result is independent of the Shields
parameter @’.
predominant type of bedforms in laboratory flumes. The ripples themselves and the
sediment transport above them have therefore been studied in considerable detail.
Sediment transport over ripples is interesting because it is at first sight
counter intuitive. That is, the net sediment transport tends to be in the opposite
direction to the strongest flow velocities, see Figure 6.3.1. Thus, the net sediment
transport over vortex ripples under Stokes waves without boundary layer drift
(Uoo(t) = Uicos mt + U2 cos 2t ) is always in the offshore direction, see the
tunnel data of Sato (1986). Correspondingly, Inman & Bowen (1963) found that
when a current was superimposed on fairly sinusoidal waves over rippled beds the
@
~
E
=
Oo
2
®
> R 4 5 6 f§ 8 9 10 11
vm 6
~
N
E 2
iS
= -2
Sa -6
= .
- -10 A
oa Test 12
@®
=-14 — — Test24
"O Spiny ohha YS eT:
Figure 6.3.1: Measured net sediment transport rates over rippled beds for three
different bed sediments under essentially the same flow conditions. Data from
Schepers (1978).
D
DT+t
0 = Ji
s§
@aeae + FJ Jiten cen arae
Ont
and, with the expressions (5.2.4), page 215, and (5.4.38), page 244 for the
sediment concentrations inserted, leads to the diffusion based transport formula
a)
Gy
9 a —_
aps JWZ)
aps ay —2z/Ls
Taz
6.3.3 A particle trajectory model for shore normal transport over ripples
The second model, the convective transport model, discussed by Nielsen
(1988a+b) was of the "particle trajectory type" and based on the assumption that
the sediment was distributed vertically by a simple convection process similar to
the one described in Section 5.4.3, page 238. For simplicity however, the vertical
convection velocity we was assumed very large, so that all sand was assumed to
reach its entrainment level ze instantaneously when the lee vortex was released. |
In this description the sand which is entrained at time f° and which totals Vv
per unit area, contributes to the total sediment distribution by
—| = FO xii(z) dz - |F'(ze)
& Jixcartie dt dze
D D Wo
(6.3.6)
zZ=0 Ze=O tet
where the first term has been transformed using integration by parts.
While this may be a fairly complicated expression, its evaluation is straight
forward compared to solving the time-dependent diffusion equation with even a
moderately complicated €;-distribution.
To illustrate the basic principles of transport modelling with models of the
particle trajectory type, we shall evaluate Equation (6.3.6) on the basis of a
particularly convenient expression for the velocity,
D D Ze/Wo
D D
1
mete
Bits Fe =a U1
u(z) dz +Birt
an. Je—Ze/Ls (1lL-— cos wo
eo dz e
Z=0 Ze=0
; ols
= w
f= +frei@at
Wo
“1,
Wo
—~— Sed,
(6.3.8)
=O Lots)
0)
; als
~ V4"
2. == pri FOIu
—— +. emi
v4 EL = Woaa | 6.3.9 OO
2=0 1 + OE.
(4)
D
oLs
—_
0, = Co]
F@)ue) dz + oti L sia
ageaatlt
oat 4020)
z=0 . 1+ Cyan
This result is very similar to the result (6.3.1) for the diffusion based model,
although the model philosophies are very different. The actual transport rates are
practically identical, see Figures 6.3.2 and further examples given by Nielsen
(1988a). It would therefore seem logical to use the convection model rather than
the diffusion model since Equation (6.3.10) is much simpler than Equation (6.3.1).
The particle trajectory approach, which is used in the convective transport
model, has a further advantage however, namely that the expression (6.3.6) can be
15
r
®
we 5
1=
S 0) Sates ~:,Gee i ee | fe | a
5
: +5
3 -10
-15
grab and dump
-20
3 4 5 6 v 8 9 10 11 12
distance from wavemaker, [m]
evaluated for almost any convective distribution function F(z) while an analytical
solution for the diffusion model is impossible for all but a few €s-distributions.
Furthermore, the convective transport model can, with minor modifications,
be used on a wave by wave basis for irregular waves, while the analytical solution
for the diffusion based transport model relies on the assumption of periodicity.
6.3.4 The grab and dump model for shore normal transport over ripples
With the present level of accuracy of coastal bottom boundary layer
calculations, it is difficult to justify very complicated sediment transport formulae.
Hence, it is worthwhile checking if further simplification is possible, and indeed it
is. It is, in fact, possible to model the shore normal sediment transport over rippled
ae 1
VT = 5a a es (6.3.11)
OF =“
for the grab and dump model. This is a remarkably simple expression which,
according to Nielsen’s (1988a) comparisons with Schepers’ (1978) laboratory
experiments, turns out to be more accurate than Equations (6.3.1) and (6.3.10)
from the more complicated models above.
The result (6.3.11) does not include a transport contribution carried by the
steady flow component u, but such a contribution may be included by adding the
usual Jem u(z) dz provided that u(z) and c(z) are known.
sediment
net
rate
transport
3 4 5 6 7 8 9 10 11 12
distance from wavemaker, [m]
Figure 6.3.3: Estimated and measured sediment transport rates as function of distance from the
wavemaker in a wave flume. For this coarse sand (ds50 = 0.465mm) only the simple grab and
dump model predicts the correct magnitude of the transport rates. Data from Schepers (1978) Test
36. For further details see Figure 6.3.1, page 267.
To understand this, consider the three formulae (6.3.1), (6.3.10) (both with the
c u-contribution ignored) and (6.3.11) in slightly rewritten and simplified forms. If
only the contribution from the first harmonic of u(z,t) is included in the diffusion
based transport formula (6.3.1), the three formulae can be written as follows
asia
nO) obs
a8 Ig SF ols - (6.3.12)
Pompe
oe ols
ete (6.3.13)
Se atercerigli (6.3.14)
pe pa, UiS
6.4.1 Introduction
Shore normal sediment transport over flat beds may occur on the beach face,
ie, in the swash zone, throughout or in parts of the surf zone and outside the surf
zone if the waves are large and the sand is fine.
normal sediment transport through the swash zone in order to model beach change,
it is beyond the present state of the art to model swash zone sediment transport.
Too little is known, at present, about the boundary layer flow and the
corresponding shear stresses in the swash zone to even attempt a description of the
basic sediment transport mechanisms in this area. See also Figure 6.1.1, page 262.
Furthermore, swash zone sediment transport includes the extra complication
of significant flows of water perpendicular to the sand surface. Experience has
shown, that increasing the inflow by pumping from the ground water in the beach,
tends to increase the rate of beach accretion.
Thus, experiments show that the flow perpendicular to the sand surface in the
swash zone affects the rate of beach accretion. The details of the mechanisms are,
however, not understood.
== 3
H.
Beaches erode if —° < 0.00070 te e )
(6.4.1)
Beaches accrete —*
Lo
where Ho and Lo are the deep water wave height and wave length respectively.
___ They found the best correlation by using the mean offshore wave height
(Ho). For a Rayleigh distribution, the mean wave height is related to the root mean
square and significant heights by H = 0.886 Hrms = 0.626 Hs.
QO = 0.00018ua (6.4.2)
where the total net transport Q is measured in m/s and ut) in m/s.
This formula was, however, only designed for 0.2mm quartz sand, and it
may be under predicting Q systematically for shorter wave periods than the ones
used by the authors (6.5s and 9./s). This was mentioned by the authors in relation
to the data of Horikawa et al (1982), and is not unexpected since shorter period
means larger d59/A and hence a larger friction factor for fixed uyms, see Section
1.2.5, page 23. It should also be remembered that formulae of this type fail to
predict the net sediment transport due to "sawtooth wave asymmetry" as discussed
in Section 2.4.4, page 121.
Not all of these limitations apply to the alternative sheet-flow sediment
transport formula (2.4.16), which was developed in Section 2.4.4, page 125. It
should be equally valid for all (non-cohesive) sediments and for waves of all
shapes and periods. However, neither of the two formulae applies, without major
modifications, to flows with large steady flow components such as rip currents or
in surf zones with a strong undertow.
In such flows, the contribution to the bed shear stress from 4u(z) must be
given special consideration. As a first, rough approximation, this problem might be
approached as if there was a linear transfer function between near-bed velocities
and the bed shear stress. The procedure would then be to subtract the time-average
from the near-bed velocities and apply Equation (2.4.15), page 125, to the residual
u(zr,t) = u(Zr,t) — Uu(Zr) ie
1
5 S25 Arms as
at Tire as ; U(Zr,tn+i1) — U(Zr,tn—1)
0’[u(tn)] = ea ad; (cos@r Wp H(zr,tn) + sings ao iea wt )
(6.4.3)
To this expression, a contribution due to the time-averaged bed shear stress (0)
must then be added in order to get the total, instantaneous effective Shields
parameter
"n)
O(n) OIF
= Olli] + gateed
UO)Fins (6.4.4)
This total value is then inserted into the sheet-flow sediment transport formula
of the momentum transfer term p (iw)oo which was discussed in Section 1.5.9,
page 91. For an example of shore normal sediment transport calculations over a
flat bed see Example 2.4.1, page 126.
6.5.1 Introduction
Sediment transport in the shore parallel direction or more precisely, the
direction perpendicular to the wave motion, is simpler to model than the transport
in the wave direction because V = 0 and the equivalent formula to Equation
(6.2.3), page 264, for the time-averaged sediment transport rate therefore becomes
Oy = \(ev+ev)d (6.5.1)
Z=0
where it seems reasonable to assume that the last term c’v’ is insignificant so that
the time-averaged sediment transport rate is given simply by
D
While this formula is valid in principle, both inside the surf zone and outside,
the details of the calculations are quite different for the two cases. This is because
of the strong influence on both c and v_ from the extra turbulence and from
additional convective mechanisms which are associated with wave breaking.
The presence of steady flow components (u(z), v(z) ) superimposed on the
wave motion i..(t) will increase the effective bed shear stresses and hence the
magnitude of c(z) compared with a pure wave motion. However, direct
experimental data which quantify this effect are not yet available.
As an educated guess however, it might be adequate to base the estimation of
the reference concentration Cg on a pure-wave formula like Equation (5.3.10),
page 228, with 025 replaced by an augmented effective Shields parameter
analogous to Equation (6.4.4). For example
Za Awr ,t
= —.
— +—
Zo Nee ey
determined in agreement with data from perpendicular waves and currents, and
from tunnel data where iw = 0, see Figures 1.5.13 and 1.5.16, pp 67 and 69.
AQ@
Ls = 0.075 pres be 0.021m (5.25)
oO
a
ig oS KV
eee (1.5.50)
eo = oe = 42
Zo i eae?
With the apparent roughness increase given by the simplistic formula (page 89)
AMO r ,v AO
ie)
ary 5 eg e
0.44 25:
(1.5.48)
and «=0.4 and with other values inserted Equation (1.5.50) becomes
2a=
F 79=0.015m, and | = e'zq = 0.041m.
With these values inserted, the current velocity distribution given by
Equations (1.5.42) and (1.5.43), page 87, becomes
Wz) = Se (6.5.5)
Vx Z Z
7 In goose 0.073 In 0015 [m/s] for z>0.04lm
This current distribution is shown in Figure 6.5.1 together with the sediment
concentration profile (Equation 6.5.4), and the distribution of the product c v.
elevation
[m]
Figure 6.5.1: The sediment concentration distribution (6.5.4) and the current
distribution (6.5.5) together with the distribution of the product Cv ( not to scale).
The estimated ripple height is 0.056m.
2
Ve gol9(Ls) 0.79
Oy = —K Co Ls {1 _ (6.5.6)
For //Ls 2 2, as in the present example, the second term can be neglected
completely so, we find the simple formula
Oy =
Ve
x ©? ;
B for I/Ls = 2 (6.5.7)
With the values from above inserted, this gives a net sediment transport rate
of Oy =2.7- 107 'm’/s.
This transport rate seems to be too small since the Meyer-Peter formula
(Equation 2.3.11, page 112), with the effective Shields parameter 0’ based on the
current friction velocity alone, gives the higher sediment transport rate
— att 1.5
Qy = ee - 005| Vs-l) gd d = 13-10 °m’/s (6.5.8)
This ought to be a lower estimate, since the additional bed shear stress contribution
of the waves is neglected.
Augmentation of the effective bed shear stress in accordance with Equation
(6.5.3), page 280, and keeping the bedform shape unchanged leads to an increase
of about 36% for the expression (6.5.7). But that is not enough to make it realistic.
Nevertheless, Rasmussen & Fredsoe (1981) used a very similar method to
that which lead to Equation (6.5.7) above, and they found good agreement with
laboratory experiments. The difference is, that their sand was fine (d59 = 0.18mm)
compared to the sand in the example above (d59 = 0.80mm).
Thus, we get a similar picture to that for the shore normal sediment transport
over rippled beds discussed in Section 6.3, p 266 ff. That is: traditional cu-integral
formulae like Equation (6.5.7) seem to work for fairly fine sand. For coarse sand,
however, the process is a very organised "grab and dump" process and is probably
better described by a "grab and dump" model.
A "grab and dump" model for longshore transport over ripples under non-
breaking waves may be constructed as follows.
Sand is picked up twice every wave period and the total amount picked up is
V = CowoT. The average distance ly travelled in the direction along the ripple
crests by that sand would be about 7/2 times the velocity at a certain near-bed
level, say one ripple height. Making use of the velocity distribution (6.5.5), that
leads to
L sioPiaxAVe
and hence, (see Figure 6.2.1, page 265) to the transport contribution
Ve
Vly = Cowol — abT
from sand entrained in one wave period, or an average grab and dump sediment
transport rate of
taist? ve 1)
Oy, = CowoT x 21 (6.5.9)
:
6)... = Co F(z). = Commman OSa (5.4.15)
(1+ 2/2)
with Co = 0.072 and z; = 0.0063m. The hydraulic roughness was found from
Equation (3.6.14) tobe 0.0026m corresponding to zo = r/30 = 0.000085m.
As in the previous example and Example 1.5.2, page 90, the velocity
distribution is derived from the information above by first finding the current
friction velocity from
be K V(z
7. = (2r) (1.5.50)
inze — nH 42 © 42
Zo vx’ A’ 2
which, with the simplistic formula for the apparent roughness increase (page 89)
Z
cay | ( | 5|= 0.44 es (1.5.48)
za = 04452 7 = 0.0028m
*
and
1 = e!zq = 0.0076m
From these parameters, the velocity distribution is found by inserting into
Equations (1.5.42) and (1.5.43), page 87
The net shore parallel sediment transport rate can now be found from
D
Oy = Javea
4)
which becomes
my if A D Stn s
—_— Vx Z. a
= Co—s |ea + |— az
2 Big J
3)
(142/21) | (142/21)
_ o D/z\ : =
Oy = Cony
— Vx Z1
2 tntt- “7
1
1+ 2
14+//7
Z1
+ J—™a +
Vz i+
nx
*)
Za
+e
(6.5.11)
which for D/z; > °° and within the range 0./< I/z; < 80 may be approximated
by
In=
ne Bs
0.52 + 0.38 cos (0.8In=> ra = Sea:
=e Vx
Oye aCe
Z1
(6.5.12)
see Figure 6.5.2. The contents of the brackets in Equation (6.5.12) generally have
magnitude about / so, for order of magnitude estimates it is reasonable to use
a v*
Oy = on 21
Goues 2 (6.5.13)
for the shore parallel sediment transport rate under non-breaking waves over a flat
bed.
0. it H oe SO
0.1 1 10
100
Figure 6.5.2: The first two terms (x) inside the brackets of Equation (6.5.11) can be
approximated by the corresponding terms in Equation (6.5.12), (the curve).
Qy = 2.9. 10°m’/s
If the bed shear stress contribution from the current is included, using
Equation (6.5.3),page 280, in the calculation of Co from Equation (5.3.10), one
finds Co = 0.074 _ instead of 0.072 or an increase of 2.8% which is
insignificant.
On the other hand, the effect of the current on the concentration distribution
through increased diffusivity may be significant. Thus, if the diffusivity €s(z) in
the combined convection-diffusion solution (5.4.49), p 248, is assumed equal to
the eddy viscosity on which the current velocity distribution (6.5.10) is based, ie
then, the concentration profile is stretched vertically, compared to the pure wave
solution, as shown in Figure 6.5.3. The net shore parallel sediment transport rate,
0.3
[m]
elevation
0.1 edocesenesascasnsaceacaccnesessees ee eee eee ne Wyereenenee et eeaccceweneccenccencccnsccescscese! sucwsesceutesvandabaec=ccbsadeunscl
pure convection
TRV ABA
Figure 6.5.4: Sediment suspension in large plunging waves at Whale Beach, Sydney,
Australia. Large amounts of sand are brought to the surface near "the plunge point"
and spread horizontally along the surface.
where the subscript "b” refers to the break point, y is the wave height to depth
ratio at the break point, a» the breaker angle, and = "ri" stands for the runup
limit. The quantity J is called the immersed weight longshore sediment transport
rate.
The value of the constant K for use in connection with the root mean square
breaker height for field data is about 0.77 (Shore Protection Manual 1984).
Recommended values of K and its scatter have recently been discussed by Bodge
& Kraus (1991) in relation to the presently available data.
One of the interesting aspects of the CERC formula is that it seems to work
quite well without considering the size of the sand - "The CERC Formula
Paradox" (Nielsen 1988b).
The weak grain size dependence was confirmed by Kamphuis (1990) who,
based on dimensional analysis and consideration of both field and laboratory data,
recommended
rl
J, te * ~1.25 fs 0.25
=
Hbrms/Tp
26-109 ||
Lop
m8 |EL
dso
sin™(208)
(6.5.16)
where mp is the beach slope at the break point, Lop is the deep water wave
length corresponding to the peak wave period Tp, and Hb,rms is the root mean
square breaker height. ©
This weak grain size dependence of the total longshore sediment transport
rate is interesting in relation to the observations of shore normal transport over
rippled beds discussed in Section 6.3.1, page 266, and to the corresponding shore
parallel transport rates discussed in Example 6.5.1, page 282. In both cases it was
found that sediment transport rates over rippled beds seems to be less sensitive to
grain size than classical transport models predict.
Classical transport models based on gradient diffusion and the assumption of
a flat bed, e g Deigaard et al (1986), predict a very strong grain size dependence,
for a given topography.
The most likely reason for the observed weak grain size dependence for the
total longshore sediment transport rate seems to lie in the topographical difference
between beaches of coarser and finer sand. Beaches of fairly coarse sand
(dso2 0.4mm) tend to develop topographies which are more conducive to
sediment transport. An example is shown in Figure 6.5.4. If the sand on that beach
had been finer, the beach profile would have been flatter and the waves would
have broken as spilling breakers with much less sediment entrainment capability.
7.1 Introduction
It has been the aim of the previous chapters to point out some of the
remaining, unanswered questions in coastal sediment transport modelling as well
as to summarise our present knowledge. Identifying loose ends or gaps in our
knowledge is the first step to any serious research project. Therefore, the following
sections summarise some of the presently unresolved problems and indicate areas
in which new insights are urgently needed.
du
arerage (7.1)
7.1
for the current component. However, the eddy viscosity vc which is defined by
this equation has some rather unfortunate characteristics. This can be seen by
considering the expression
ve + Vv (7.2)
The first deterministic term, (— uw), in this expression may often be ignored
since w_must, by continuity, be zero for horizontally uniform flows. The second
term, (—i w ), however, will often be dominant and this gives the eddy viscosity Vc
a problematic nature.
Hence, vc is not a "turbulent eddy viscosity", because the leading
contribution (iw) is deterministic.
It is not isotropic, but depends strongly on the direction of the current relative
to the direction of wave propagation.
___ The eddy viscosity given by Equation (7.2) may in fact become negative if
—uw is dominant and of the form, derived by Longuet-Higgins (1956), (see pages
55 and 56), while a is negative. This corresponds to the commonly observed
situation which is shown in Figure 7.1.
___ The bed shear stress is positive but smaller than the asymptotic value of
—u w. See also Figure 1.4.2, page 55.
The steady flow velocity is positive close to the bed, but starts to decrease at
the elevation where
Above this point, the eddy viscosity vc, is negative because T and du/dz
have opposite signs.
A negative eddy viscosity is unsatisfactory in itself and even more so in
relation to sediment diffusivity. It cannot be avoided however, if the time-averaged
equation of motion is written in the commonly used form (7.1).
These conceptual problems with the eddy viscosity are, however, avoided if
_—Oo
(oe)
relative
z/5s
elevation
0
=1.9 -1 -0.5 0 0.5 1 1.5
Figure 7.1: In situations where the time-averaged bed shear stress is positive but
smaller than the asymptotic value of —iuw , the shear stress T(z) and the current
gradient a have opposite signs in the upper part of the flow. This corresponds to
negative values of the eddy viscosity defined by Equation (7.1).
current becomes
du UZ 7
Vel ‘de = asta] + Uw (7.4)
With this equation it is easier to explain the situation shown in Figure 7.1.
The new eddy viscosity, Vc1, which is given by (for uw = 0)
ee ee (7.5)
is, unlike vc, basically a turbulent eddy viscosity and it is unlikely to be strongly
anisotropic. It is also more probable that vc1 can be meaningfully related to the
sediment diffusivity. The latter is important with respect to sediment transport
modelling.
With the time-averaged equation of motion written in the form (7.4) it is also
possible to explain the decrease in current velocity near the suface which can be
seen in Figure 1.5.2, page 63. The situation in Figure 1.5.2 corresponds to a
positive bed shear stress (0) of the order —1.5p (UW). According to
Equation (7.4), with
under waves, which is not directly related to the effective bed shear stress. It
would, however, be very interesting to see some direct measurements of the flow
structure over flat beds of loose sand under oscillatory flows.
This equation, which is the usual type of boundary condition for steady,
uniform flow, expresses equilibrium between suspended sediment concentrations
at the bed and the effective Shields parameter 9’. Such an equilibrium does,
however, not exist in non-uniform or unsteady flows, see Section 5.3.1 (pp
222-223). A different boundary condition has therefore been suggested by Nielsen
et al (1978) and van Rijn (1984).
Their approach is to treat deposition and entrainment separately, and then try
to relate the entrainment rate, or pickup rate p(t), directly to the effective,
instantaneous Shields parameter 9’(t).
The pickup function model is, however, still a simplified, formal description
Allen, J R L (1982): Sedimentary structures: Their character and physical basis. Elsevier,
Amsterdam.
Amos, C L & MB Collins (1978): The combined effect of wave motion and tidal currents
on the morphology of intertidal ripple marks: the Wash. U K J Sedimentary Petrology
Vol 48, No 3 , pp 849-856.
Amos, CL, A J Bowen, D A Huntley & C F M Lewis (1988): Ripple generation under the
combined influences of waves and currents on the Canadian continental shelf.
Continental Shelf Res, Vol 8, No 10, pp 1129-1153.
Amott, R W & J B Southard (1990): Exploratory flow duct experiments on combined flow
bed configurations and some implications for interpreting storm event stratification. J
Sedimentary Petrology, Vol 60, No 2, pp 211-219.
Asano, T, M Nagakawa & Y Iwagaki (1986): Changes in current profiles due to wave
superimposition. Proc 20th Int Conf Coastal Engineering, Taipei, pp 925-940.
Bagnold, R A, (1946): Motion of waves in shallow water: Interaction of waves and sand
bottoms, Proc Roy Soc Lond, A 187, pp 1-15.
Bagnold, R A (1954): Experiments on a gravity-free dispersion of large solid spheres in a
newtonian fluid under shear. Proc Roy Soc Lond A 225, pp 49-63.
Bagnold, R A, (1956): The flow of cohesionless grains in fluids. Phil Trans Roy Soc Lond,
No 964, Vol 249, pp 235-297.
Bailard, J A (1981): An energetics total load sediment transport model for a plane sloping
beach. J Geophys Res, Vol 86, No C11, pp 10938-10954.
Bakker, W T, & T van Doorn, (1978): Near bottom velocities in waves with a current. Proc
16th International Conference on Coastal Engineering, ASCE, Hamburg, pp 1394-1413.
Basset, A B, (1888): On the motion of a sphere in a viscous fluid. Phil Trans Roy Soc Lond,
Ser A, Vol 179, pp 43-63.
Bhattacharya, P K, (1971): Sediment suspension in shoaling waves, Ph D thesis, University
of Iowa, Iowa City, Iowa.
Bijker, E W (1967): Some considerations about scales for coastal models with movable
beds. Delft Hydraulics Lab Pub No 50.
Blondeaux, P & G Vittori (1990): Oscillatory flow and sediment motion over a rippled bed.
Proc 22nd Int Conf Coastal Engineering, Delft, pp 2186-2199.
Boczar-Karakiewicz, B & R D L Davidson-Arnott (1987): Nearshore bar formation by
non-linear wave processes - A comparison of model results and field data. Marine
Geology, Vol 77, pp 287-304.
Boczar-Karakiewicz, B & L A Jackson (1991): Beach dynamics and protection measures in
the Gold Coast area, Australia. Proc 10th Australasian Conf on Coastal and Ocean
Engineering, Auckland, pp 411-415.
Bodge, K R & N C Kraus (1991): Critical examination of longshore sediment transport rate
magnitude. Proc " Coastal Sediments ’91", AS CE, pp 139 - 155.
Bosman, J J (1982): Concentration measurements under oscillatory water motion. Delft
Hydraulics report M1695 part 2.
Bosman, J J, ET J M van der Velden & C H Hulsbergen (1987): Sediment concentration
measurements by transverse suction , Coastal Engineering, Vol 11, pp 353-370.
Bretschneider , C L (1954): Field investigation of wave energy loss of shallow water ocean
waves, Beach Erosion Board, Tech Memo 46.
Brevik, I, (1981): Oscillatory rough turbulent boundary layers, J Waterway Port Coastal
Ocean Div. ASCE, 107 , pp 175-188.
Bruun , P (1962):Sea level rise as a cause of shore erosion. Proc A SC E,, Vol 88, WW1, pp
117-130.
Bruun, P (1983): Review of conditions for uses of the Bruun Rule of erosion. Coastal
Engineering , Vol 7, pp 77-89.
Bruun, P (1990): Port Engineering, IV edition, Vol 2, Gulf Publishing Company.
Cacchione, D A & D E Drake (1982): Measurements of storm generated bottom stresses on
the continental shelf. J Geophys Res, Vol 87, No C3, pp 1952-1960.
Carstens, M R, FM Neilson & H D Altinbilek (1969): Bedforms generated in the laboratory
under an oscillatory flow, CE RC Tech Memo 28.
Carter T G, P L-F Liu, and C C Mei (1973): Mass transport by waves and offshore
bedforms. J Waterways Harbours & Coastal Division, A S C E, Vol 99, WW2, pp
165-184.
Christoffersen, J B & I G Jonsson (1985): Bed friction and dissipation in combined current
and wave motion, Ocean Eng, Vol 12, No 5, 387-423.
Clift, R, R J Grace, & M E Weber, (1978): Bubbles, drops and particles. Academic Press,
New York
Clifton, H E (1976): Wave-formed sedimentary structures - A conceptual model. In Davis &
Ethington ed: Beach and nearshore sedimentation. S EP M special publication 24, pp
126-148.
Coffey, F C (1987): Current profiles in the presence of waves and the hydraulic roughness
of natural sand beds. Ph D Thesis, Univ of Sydney, 252 pp.
Coffey, F C & P Nielsen (1984): Aspects of wave-current boundary layer flows. Proc 19th
Int Conf Coastal Engineering, Houston Texas, pp 2232-2245.
Coffey, F C & P Nielsen (1986): The influence of waves on current profiles. Proc 20th Int
Conf Coastal Eng, Taipei, pp 82-96.
Coleman, N L, (1970): Flume studies of the sediment transfer coefficient, Water Resource
Res, 6(3), pp 801-809.
Dean, RG & R A Dalrymple (1991): Water wave mechanics for engineers and scientists.
World Scientific, Singapore.
Deigaard, R, J Fredsoe & I B Hedegaard (1986a): Suspended sediment in the surf zone. J
Waterway, Port, Coastal and Ocean Eng, ASCE, Vol 112, No 1, pp 115-128
Deigaard, R, J Fredsoe & I B Hedegaard (1986b): A mathematical model for littoral drift. J
Waterway, Port, Coastal and Ocean Eng, ASCE, Vol 112, No 3, pp 351-369.
Deigaard, R, P Justesen & J Fredsoe (1991): Modelling undertow by a one-equation
turbulence model. Coastal Engineering, Vol 15, pp 431-458.
Delft Hydraulics (1989): Bedforms, near-bed sediment concentrations and sediment
transport in simulated regular wave conditions. Report No H840.
Dick, JE & J FA Sleath (1991): Velocities and concentrations in oscillatory flow over beds
of sediment. J Fluid Mech, Vol 233, pp 165-196.
Dingler, J R (1974): Wave formed ripples in nearshore sands. Ph D thesis Univ of
California, San Diego, 136 pp.
Downing, J P, R W Sternberg & C R B Lister (1981): New instrumentation for the
investigation of sediment suspension processes in shallow marine environments. Marine
Geology, Vol 42, pp 14-34.
Du Toit, CG, & JF A Sleath, (1981): Velocity measurements close to rippled beds in
oscillatory flow, J Fluid Mechanics 112, pp 71-96
Engelund,F (1970): Instability of erodible beds. J Fluid Mech, Vol 42, pp225-244.
Engelund, F (1981): Transport of bed load at high shear stress. Progr Rep 53, Inst
Hydrodynamic and Hydraulic Eng, Tech Univ Denmark, 31-35.
Engelund, F & E Hansen (1972): A monograph on sediment transport in alluvial streams.
Teknisk Forlag, Copenhagen.
Fernandez-Luque, R (1974): Erosion and transport of bed-load sediment. Dissertation,
KRIPS Repro BV, Meppel, The Netherlands.
Fredsoe, J (1984): Turbulent boundary layer in wave-current motion, J Hydr Eng, ASC E,
Vol 110, HY8, 1103-1120.
Gibbs, R J, M D Mathews & D A Link (1971): The relationship between sphere size and
settling velocity. J Sed Petrology, Vol 41, pp 7-18.
Gilbert, G K (1914): The transportation of debris in running water. U S Geol Survey, Prof
Paper 86.
Grant, W D, & O S Madsen, (1979): Combined wave and current interaction with a rough
bottom, JGeophysical Res, 84, 1808
Grant, W D, & O S Madsen, (1982): Movable bed roughness in unsteady oscillatory flow. J
Geophys Res, Vol 87, pp 469-481.
Grant, W D, J A Williams, S M Glen, D A Cacchione & D E Drake (1983): High frequency
bottom stress variability and its prediction in the Code Region. Woods Hole Oceanogr
Inst, Tech Rep 83-19.
Guy, H P, DB Simons, & E V Richardson (1966): Summary of alluvial channel data from
flume experiments, 1956-1961. U S Geological Survey, Prof paper 462-I, Washington
DC.
Hallermeier R J (1980): Sand motion initiation by water waves: Two asymptotes. J
Waterway, port,coastal and ocean Div, A SC E, Vol 106, 299-318.
Hallermeier R J (1981): A profile zonation for seasonal sand beaches from wave climate.
Coastal Engineering, Vol 4, pp 253-277.
Hallermeier, R J (1981b): Seaward limit of signficant sand transport by waves. CETA 81-2,
CERC, Ft Belvoir, Va.
Hammond, T M & MB Collins (1979): On the threshold of transport of sand-sized sediment
under the combined influence of unidirectional and oscillatory flows. Sedimentology,
Vol 26, pp 795-812.
Hanes, D M (1990): The structure of events of intermittent suspension of sand due to
shoaling waves, in The Sea, Vol 9 Part B, John Wiley & Sons, pp 941-954.
Hanes, D M & A J Bowen (1985): A granular fluid model for steady intense bed-load
transport. J Geophys Res, Vol 90, No CS, pp 9149-9158.
Hanes, D M & D A Huntley (1986): Continuous measurements of ‘suspended sand
concentration in wave dominated nearshore environment. Continental Shelf Res, Vol 6,
No 4, pp 585-596.
Hanes, D M & D L Inman (1985a): Observations of rapidly flowing granular-fluid
materials. J Fluid Mech, Vol 150, 357-380.
Hanes, D M & D L Inman (1985b): A dynamic yield criterion for granular-fluid flows. J
Geophys Res, Vol 90, No BS, 3670-3674.
Hansen, J B & I A Svendsen (1986): Experimental investigation of the wave and current
motion over a longshore bar. Proc 20th Int Conf Coastal Eng, Taipei, pp 1166-1179.
Hardistry, J & J P Lowe (1991): Experiments on particle acceleration in stationary and
oscillating flows. J Hydraulic Engineering, ASC E, Vol 117.
Hattori, M, (1969): The mechanics of suspended sediment due to standing waves, Coastal
Engineering Japan, 12, 69-81
Hedegaard, I B, R Deigaard & J Fredsoe (1991): Onshore/offshore sediment transport and
morpho- logical modelling of coastal profiles. Proc Coastal Sediments ’91, Seattle, pp
643-657.
Ho, H W, (1964): Fall velocity of a sphere in an oscillating fluid. PhD thesis, University
Iowa, Iowa City, Iowa
Homma, M, K Horikawa & R Kajima (1965): A study of suspended sediment due to wave
action. Coastal engineering in Japan, Vol 8, pp 85-103.
Horikawa, K, & A Watanabe, (1968): Laboratory study on oscillatory boundary layer flow,
Coastal Engineering Japan 11, 13-28
Horikawa, K, A Watanabe & S Katori (1982): Sediment transport under sheet flow
condition. Proc 18th Int Conf on Coastal Eng, Capetown, pp 1335-1352.
Inman, D L (1957): Wave generated ripples in nearshore sands. Beach Erosion Board, U S
Army Corps of engineers, Tech Memo 100.
Inman, D L & A J Bowen (1963): Flume experiments on sand transport by waves and
currents. Proc 8th Int Conf Coastal Engineering, A S C E, pp137-150.
Iwagaki Y & T Kakinuma, (1967): On the bottom friction factors off five Japanese coasts,
Coastal Eng Japan, Vol 10, 13-22.
Jackson, P S (1981): On the displacement height in the logarithmic velocity profile. J Fluid
Mech Vol 111, pp 15-25.
Murray, S P (1970): Settling velocities and vertical diffusion of particles in turbulent water,
J Geophysics Res. Vol 75, No 9, 1647-1654.
Myrhaug, D (1982): On a theoretical model of rough turbulent wave boundary layers.
Ocean Eng, Vol 9, No 6, pp 547-565.
Myrhaug, D & O H Slaattelid (1989): Combined wave and current boundary layer model for
fixed, rough seabeds. Ocean Engineering Vol 16, No 2, pp 119-142.
Nadaoka, K, S Ueno & T Igarashi (1988): Sediment suspension due to large eddies in the
surf zone. Proc 22nd Int Coastal Eng Conf, Torremolinos, pp 1646-1660.
Nakato, T, F A Locher, J R Glover, and J F Kennedy (1977): Wave entrainment of sediment
from rippled beds, Proc. ASCE, 103 (WW1), 83-100.
Natarajan, P (1969): Sand movement by combined action of waves and currents. Ph D thesis,
University of London.
Nielsen, P, (1979): Some basic concepts of wave sediment transport, Ser. Paper 20, Inst
Hydrodyn Hydraul Eng, Tech Univ Denmark, 160 pp.
Nielsen, P (1981): Dynamics and geometry of wave generated ripples. J Geophys Res, Vol
86, No C7, pp 6467-6472.
Nielsen, P (1983): Entrainment and distribution of different sand sizes under water waves. J
Sedimentary Petrology, Vol 53, No 2, pp 423-428.
Nielsen, P (1983): Analytical determination of nearshore wave height variation due to
refraction, shoaling and friction. Coastal Engineering, Vol 7, pp 233-251.
Nielsen, P (1984a): On the motion of suspended sand particles. J Geophys Res, Vol 89, No
C1, pp 616-626.
Nielsen, P (1984b): Field measurements of time-averaged suspended sediment
concentrations under waves. Coastal Engineering, Vol 8, pp 51-72.
Nielsen, P, (1985): On the structure of oscillatory boundary layers, Coastal Engineering 9,
261-276.
Nielsen, P (1986): Suspended sediment concentrations under waves. Coastal Engineering,
Vol 10, pp 23-31.
Nielsen , P (1988a): Three simple models of wave sediment transport. Coastal Engineering,
Vol 12, pp 43-62.
Nielsen, P (1988b): Towards modelling coastal sediment transport. Proc 21st Int Conf
Coastal Eng, Torremolinos, pp 1952-1958.
Nielsen, P (1990): Coastal bottom boundary layers and sediment transport. In P Bruun ed
Port Engineering (4th edition), Vol 2, pp 550-585.
Nielsen, P, I A Svendsen & C Staub (1978): Onshore-offshore sediment transport on a
beach. Proc 16th Int Conf Coastal Eng, Hamburg, pp 1475-1492.
Nielsen, P, N R Sena & ZJ You (1990): The roughness height under waves. J Hydraulic
Res, Vol 28, No 5, pp645-647.
Nikuradse, J (1933): Stromungsgesetze in glatten und rauhen rohren. V D J Forschungsheft
361, Berlin.
Owen, P R (1964): Saltation of uniform grains in air. J Fluid Mech, Vol 20, No 2, pp
225-242.
Peregrine, J H & I A Svendsen (1978): Spilling breakers, bores and hydraulic jumps. Proc
16th Int Conf Coastal Eng, Hamburg, pp 540-551.
Rasmussen, P & J Fredsoe (1981): Measurements of sediment transport in combined waves
and current. Inst Hydrdyn & Hydraul Eng (ISVA), Tech Univ Denmark, Progress report
53, pp 27-30.
Raudkivi, A J (1988): The roughness height under waves. J Hydraulic Res, Vol 26, No 5,
pp569-584.
Reizes, J A, (1977) A numerical study of the suspension of particles in a horizontally
flowing fluid, paper presented at 6th Australasian Hydraulics and Fluid Mechanics
Conference, Adelaide .
Ribberink, J S and A Al-Salem (1989): Bed forms, near-bed sediment concentrations and
sediment transport in simulated regular wave conditions. Delft Hydraulics Report H840,
part 3.
Ribberink, J S and A Al-Salem (1990): Bed forms, sediment concentrations and sediment
transport in simulated wave conditions. Proc 22nd Int Conf Coastal Engineering, Delft,
pp 2318-2331.
Riedel, H P (1972): Direct measurement of bed shear stress under waves. Ph D Thesis, Dept
Civ Eng, Queen’s University, Kingston, Ontario.
Roelvink, J A & M J F Stive (1989): Bar generating cross shore flow mechanisms on a
beach. J Geophys Res , Vol 94, No C4, pp 4785 - 4800.
Sato, S (1986): Oscillatory boundary layer flow and sand movement over ripples. Ph D
thesis, Dept of Civ Eng, Univ of Tokyo.
Sato, Y & K Yamamoto (1987): Lagrangian measurement of fluid particle motion in an
isotropic turbulent field, J Fluid Mech, Vol 175, 183-199.
Sawamoto, M & T Yamashita (1986): Sediment transport rate due to wave action. J of
Hydroscience and Hydraulic Eng, Vol 4, No 1, pp1-15.
Schepers , J D (1978): Zandtransport onder invloed van golven en een eenparige stroom bij
varierende korreldiameter. M Eng Thesis, Delft Univ of Technology.
Schlichting, H (1979): Boundary layer theory, 7th ed, McGraw-Hill, New York.
Shields, A (1936): Anwendung der Aehnlichkeitsmechanik und Turbulenzforchung auf die
Geschiebebewegung. Mitt Preuss Versuchsanstalt fur Wasserbau und Schiffbau, No 26,
Berlin.
Shore Protection Manual (1984), U S Army Coastal Engineering Research Centre,
Vicksburg Mississippi.
Simons, R R, A J Grass & A Kyriacou (1988): The influence of currents on wave height
attenuation. Proc 21st Int Conf Coastal Engineering, Malaga.
Slaattelid, OH, D Myrhaug & K F Lambrakos (1990): North Sea bottom steady boundary
layer measurements. J Waterway, Port, Coastal and Ocean Engineering, Voll16, No 5,
pp 614-633.
Sleath, J F A (1970): Measurements close to the bed in a wave tank. J Fluid Mech, Vol 42,
pp 111-123.
Sleath, J F A (1978): Measurements of bed-load in oscillatory flow. Proc A SC E, Vol 104,
No WW4, 291-307.
Sleath, J F A (1982): The suspension of sand by waves. J Hydraulic Res, Vol 20, No 5, pp
439-452.
Sleath, J F A (1984): Sea Bed Mechanics, Wiley Interscience.
Sleath, J F A (1985): Energy dissipation in oscillatory flow over rippled beds, Coastal Eng,
Vol 9, 159-170.
Sleath, J F A (1987): Turbulent oscillatory flow over rough beds, J Fluid Mech, Vol 182,
369-409.
Sleath, J F A (1990): Bed friction and velocity distributions in combined steady and
oscillatory flow. Proc 22nd Int Conf Coastal Eng, Delft, pp 450-463.
Sleath, J F A (1991): Velocities and shear stresses in wave-current flows. J Geophys Res,
Vol 96, No C8, pp 15237-14244.
Snyder, W H & J L Lumley (1971): Some measurements of particle velocity autocorrelation
functions in a turbulent flow, J Fluid Mech, Vol 48, 41-71.
Southard, J B, JM Lambie, D C Federico, H T Pile, & C R Weidman (1990): Experiments
on bed configurations in fine sands under bidirectional purely oscillatory flow, and the
origin of hummocky cross-stratification. J Sedimentary Petrology, Vol 60, No 1, pp 1-17.
Spalart, PR & B S Baldwin (1987): Direct simulation of a turbulent oscillating boundary
layer, NASA Tech Memo 89460, Ames Res Centre, Moffett Field, Ca.
Staub, C, IG Jonsson & I A Svendsen (1984): Variation of sediment suspension in
oscillatory flow. Proc 19th Int Conf Coastal Eng, Houston, pp 2310-2321.
Stive, M J F (1988): Cross-shore flow in waves breaking on a beach. Dissertation, Delft
University of Technology.
Svendsen, I A & P A Madsen (1984): A turbulent bore on a beach. J Fluid Mech, Vol 148,
pp73-96.
Svendsen, I A, H A Schaffer & J Buhr Hansen (1987): The interaction between the
undertow and the boundary layer flow on a beach. J Geophys Res, Vol 92, pp
11845-11856.
Swart , D H (1974): Offshore sediment transport and equilibrium beach profiles. Delft Hydr
Lab Publ No 131.
Taylor, G I (1921) Diffusion by continuous movement, Proc Lond Math Soc, Vol 20,
196-211.
Tooby, P F, G L Wick, & J D Isacs, (1977): The motion of a small sphere in a rotating
velocity field: A possible mechanism for suspending particles in turbulence, J
Geophysical Res, 82 (15), 2096-2100.
Trowbridge, J & O S Madsen (1984): Turbulent wave boundary layers I. Model formulation
and first order solution. J Geophys Res, Vol 89, No C5, pp 7989-7997.
van Doorn, T, (1981): Experimental investigation of near-bottom velocities in water waves
without and with a current. TOW Report M 1423 part 1, Delft Hydraulics Laboratory .
A\-Salem, A 121, 131, 136, 278 35, 36, 37, 41, 44, 46, 48
Carter, TG 59
Allen, JRL 151
Carstens, MR 4, 5, 27, 28, 29, 45,
Altinbilek, H D 4, 5, 27, 28, 29, 45,
105, 106, 135, 139, 140, 141, 142,
105, 106, 107, 135, 139, 140, 141,
146, 147, 152, 153, 156, 157, 158,
142, 146, 147, 152, 153, 156, 157,
159, 296
158, 159, 296
CERC 291, 292
Amos, CL 108, 143, 144
Christoffersen, J B 73, 75, 81, 85
Arnott, R W 143, 145
Clift, R 167
Asano, T 75, 83, 88
Clifton, HE 131, 136
Coffey, FC 69, 70, 74, 82, 83, 143,
Bagnold, R A 6, 25, 27, 28, 96, 97, 145
98, 109, 110, 111, 114, 115, 116, Coleman, N L 204
135, 157, 207 Collins, M B 108, 143
Bailard, JA 121
Bakker, W T 51, 52
Dairymple, R A 28, 146
Baldwin, BS 52
Davidson-Arnott ,R DL 277
Basset, A B 167, 170
Dean, R G 28, 146
Bhattacharya, PK 161, 176
Deigaard, R 61, 277, 291, 292
Bijker, E W 84, 143
Delft Hydraulics 140, 141, 142,
Blondeaux, P 41
216, 218, 259
Boczar-Karakiewicz, B 277
Dick, JE 298
Bodge, KR 292
Dingler, JR 137, 140, 141, 142
Bosman, J J 209, 216
Downing, J P 209
Bowen, AJ 108, 115, 143, 144, 267
Drake, DE 83
Bretschneider, C L 28
Du Toit, CG 228
Brevik, I 38, 48, 49, 50
Bruun, P 109
Engelund, F 104, 105, 112, 134,
Cacchione, DA 83
145, 149
Carlsen, N A 4, 14, 15, 26, 27, 32,
Convective distribution model 205, Diffusion equation 19, 203, 236, 241
238-242, 244-245, 249-251, Diffusive distribution model
253, 260-261 241-245, 248, 257-258
Convective sediment flux 234-236, Diffusive sediment flux 202-203,
238-239, 241, 246-247 234-236, 241, 243, 246-247
Convective transport model 270-276 Diffusivity 201-205, 236, 241, 243-
Current 245, 252-255, 258, 269, 289, 298
boundary layer concepts 64-68 Dispersion 190, 193-200
eddy viscosity 36, 67-68, 84-94, Dispersive stress 96, 111, 207
293-296 Displacement thickness 29-30
Eulerian drift 53-54, 56-57 Dissipation factor, see Energy
influence on wave boundary layer dissipation factor
73-75 Distribution function 237, 252-253
Lagrangian drift 54-60 Distribution models 233-261, 298
mass transport 54-60 Drag coefficient 100, 165
profile 62-64, 78-94, 293-296 Drag force 7, 98-104, 107, 110, 147,
shear stress distribution 64 162, 164, 167, 169, 171-172, 180
undertow 60-61 Drift velocity
wave-current boundary layer Eulerian 53-54, 56-57
62-94 Lagrangian 54-60
wave-generated 52-61 Dunes 186
Critical Shields parameter 107-108,
144-145
Eddy viscosity
Critical stress 108
based models 48-52, 70, 84-87,
293
Darcy’s law 102 combined flow 68-71, 92,
Decay, wave height 77 293-296
Defect function 19-22, 29-30, oscillatory flow 31-40, 57, 68, 71
38-40, 43-47,49-52, 74 relation with sediment diffusivity
Defect velocity 7, 19, 29-30, 152 204-205
Deposition 131, 222-223, 239 steady flow 36, 67-68
Depth of closure 109 turbulent component 32-34
Diffusion 190, 193-205, 219, 221, wave-current flow 68-71, 92,
233-236, 241-261 293-296
Diffusion-convection distribution Effective normal stress 95-98,
model 205, 245-261 110-111
seston. 471% _—
= - a
ABOUT THE AUTHOR