8 Handout
8 Handout
8 Handout
能流 )
1. Energy flux (能
Progressive water waves transmit energy. The rate at which the energy is transferred is
called the energy flux, F. The instantaneous rate (per unit width) at which work is
being done in the direction of wave propagation by the pressure due to a wave is
Z η Z η
F= p · udz = (p0 + pd ) · udz, (1)
−h −h
t+T η
cosh k(h + z)
Z Z
1
= (−ρgz) · ωη dzdt
T t −h sinh(kh)
, (3)
Z t+T Z η
ρgω 1
=− ηdt z cosh k(h + z) dz
sinh(kh) T t −h
Z t+T Z η
ρgω 1 H
=− cos(kx − ωt)dt z cosh k(h + z) dz = 0
sinh(kh) T t 2 −h
t+T η
cosh k(h + z) ih cosh k(h + z) i
Z Z
1 h
= ρgη ωη dzdt
T t −h cosh(kh) sinh(kh)
cosh2 k(h + z)
Z t+T Z 0
1 2 (4)
' ρgωη dzdt
T t −h cosh(kh) sinh(kh)
1
ESOE 5136 2021 May 10 - Handout Lo
Cg = C · n g (7)
is called the group velocity (波 波 群 速 度 ). The group velocity indicates how fast the
energy is transmitted by the wave.
From the expression for ng in (5), we see that
p
C g ' C = gh, in the shallow water limit (kh → 0)
. (8)
C ' 1 C = gT , in the deep water limit (kh → ∞)
g
2 4π
折射 )
2. Wave refraction (折
When the speed of a wave changes due to a change in water depth, its propagation direction
also changes. For a 2D bathymetry (海域地形), in which the water depth does not change
in the alongshore direction (沿岸方向), the refraction process obeys Snell’s law (司 司乃 耳
定律 ):
sin θ1 sin θ2
= , (9)
C1 C2
where θ1 is the incidence angle at location 1, C1 is the wave speed at location 1, θ2 is the
incidence angle at location 2, C2 is the wave speed at location 2. The refraction process is
sketched in Figure 1.
Water waves slow down as they get close to shore. Due to wave refraction, the wave
direction θ also decreases. As a result, waves tend to approach shore normally.
2
ESOE 5136 2021 May 10 - Handout Lo
波浪 淺 化 )
3. Wave shoaling (波
As water waves approach shore, the water depth gets smaller and the waves slow down.
Due to the conservation of energy, the wave height increases and the wavelength decreases,
but the wave period remains unchanged.
Consider two wave rays approaching shore as sketched in Figure 3. Under the assump-
tions that there are no reflection and energy changes, the conservation of energy (strictly
speaking, energy flux, since it’s per time) requires
E0 n0 C0 b0 = E1 n1 C1 b1 = E2 n2 C2 b2 . (10)
Using the expression for E from (6) and considering only location 1 and location 2 for
brevity, we then have
s r r r
H2 Cg1 b1 n1 C1 b1
= = , (11)
H1 Cg2 b2 n2 C2 b2
• For a 2D bathymetry that doesn’t vary pin the alongshore direction, b0 = l0 cos θ0 ,
b1 = l0 cos θ1 , and b2 = l0 cos θ2 . Then, b1 /b2 can be calculated as
r r
b1 cos θ1 1 − sin2 θ1 1/4
= =( ) < 1, (12)
b2 cos θ2 1 − sin2 θ2
which indicates that on a 2D bathymetry, the wave heights will decrease towards
shore due to refraction.
In
p the case of normal incidence, θ0 = θ1 = θ2 = 0, refraction does not occur and
b1 /b2 = 1.
• If the water depths at location 1 and location 2 are both shallow, (11) simplifies to
s√ r
H2 h b1
= √ 1 , in the shallow water limit. (13)
H1 h2 b2
3
ESOE 5136 2021 May 10 - Handout Lo
H2 h2
∝ ( )−1/4 , for shallow water waves, (14)
H1 h1
which states that the wave height increases as a shallow water wave propagates to-
wards shore.
In addition, in the shallow water limit,
r r r r s√ r
H2 C1 b1 L1 T2 b1 h1 b1
= = = √ . (15)
H1 C2 b2 T1 L2 b2 h2 b2
L2 h2
= ( )1/2 , for shallow water waves, (16)
L1 h1
which states that the wavelength decreases as a shallow water wave propagates to-
wards shore.
4
ESOE 5136 2021 May 10 - Handout Lo
r1
r2
wave crest
travel direction
In constant water depth, the wave speed is fixed. The wave spreads in all directions at the
same speed. As a result, for a radially symmetric wave shape spreading radially outwards,
the shape remains radially symmetric. Let us consider two instants, t = t1 and t = t2 . At
t = t1 , the distance between the wave crest and the origin is r1 – the wave forms a ring of
radius r1 . At t = t2 , the distance between the wave crest and the origin is r2 – the wave
forms a ring of radius r2 . If we assume the wavelength L remains unchanged, then the
areas covered by the wave at the two instants are
L L L L
A1 = π(r1 + )2 − π(r1 − )2 = 2πr1 L, A2 = π(r2 + )2 − π(r2 − )2 = 2πr2 L. (19)
2 2 2 2
Due to energy conservation, the total wave energy – wave energy density (which scales as
H 2 ) multiplied by wave area (which scales as r) – should remain constant. Thus, we write
r
2 2 r1
H1 r1 = H2 r2 ⇒ H2 = H1 , (20)
r2
which shows that the wave height scales as
H ∝ r−1/2 . (21)