Théorie de La Houle

ZORKANI Mohammed
Dpartement dHydraulique
Chapitre 1
Equations des ondes en bidimensionnelle
Caractristiques des Houles
1) Introduction :
Il existe dans la nature plusieurs types donde marine : les ondes dues
au vent, ondes (vagues) engendres par un navire, tsunamis [ondes
engendres par un tremblement de terre :secousse tellurique] et la
mare que le gnie maritime doit connatre.
battement de surface , ressac , seiches
nergie
relative
h
na
vi
ondes
capillaires
re
vagues
surf
beat
tsunamis
10
1
Vent
gravitation
M2
e
tempte
10 2
5mn
10 3
50mn 12h
10 4
transtidales
lune et soleil
mare
30s
0,1
activit sismique
S2
temptes
vents
24h
10 5
T en sec
Spctre nertique des ondes selon Munk

Non de londe
Capillarit
Ultragravit
Gravit
Infragravit
Priode longue
Transtidal
Intervalle de priode (sec)

de 0 0,1
de 0,1 1
de 1 30
de 30 300
de 300 86400
de 86400
Le spectre des ondes dues au vent stendent jusqu 30s. Au del,

les infra vagues sont dues des interactions non linaires entre
vagues de vent de priodes beaucoup plus courtes. Des seiches dues
des variations de la pression atmosphrique qui provoquent des
oscillations stationnaires libres de priode allant de 30s 5 mn [le lac de
Genve a des seiche de 63mn (uninodale) et 36mn (binodale)&dAmplitude ~50cm].
Les ondes marines sont donc complexes car dans un train donde le
spectre possde des composantes de diffrent dphasage dont les
priodes et les amplitudes obissent une loi de distribution alatoire,
on recourt pour dpasser cette difficult aux approches statistiques
depuis 1945. Dans ce chapitre on va dvelopper une thorie
monochromatique : priode constante et amplitude uniforme en
profondeur constante. On rencontre dans la littrature plusieurs
E.H.T.P.
1
Equations des ondes en bidimensionnelle et caractristiques des houles :ondes de surface
ZORKANI Mohammed
thories, la plus simple et la plus utilise est celle faible amplitude

[thorie linaire des ondes] propose par Airy en 1845 [Thorie dAiry].
2) Thorie linaire des ondes : faible amplitude (infinitsimales)
La thorie bidimensionnelle des ondes de gravit priodiques a faible
amplitude libre [pas en engendrement par le vent] est bas sur la
linarisation des conditions limites la surface libre et sur le fond ; on
admet que lcoulement est irrotationnel : lcoulement est donc
potentielle de vitesse. Les hypothses sont :
Fluide homogne, incompressible
Leffet de la tension superficielle ngligeable [longueur donde
grande devant 3cm].Coriolis ngligeable: T f
r
v r u
(2 u)
t
Ecoulement irrotationnel : pas de contraintes qui sexercent sur

la surface libre (pas de vent ) : on dit onde libre et pas de
frottement sur le fond (fluide parfait) : le fluide glisse librement
sur le fond et toutes autres surfaces solides.
r
Il existe alors un potentiel de vitesse : v = grad (x, z, t ) selon lquation
r
de continuit ( div v = 0 ) ce potentiel vrifie :
2
x2
2
z2
= 0 (1 1)
Le fond est horizontal fixe(pas de transport de sdiment), impermable.

La pression atmosphrique qui sexerce sur la surface libre est
constante donc pas de vent : la surface se dforme librement.
Lamplitude de londe est petite vis vis de la longueur donde
et de la profondeur deau. Puisque les vitesses des particules
fluides sont proportionnelles lamplitude alors que la clrit
(vitesse de phase) de londe varie avec la profondeur ; cette
condition se traduit par les particules fluides ont une vitesse
faible compares la clrit. Ce qui signifie que cette thorie
nest applicable que dans les eaux suffisamment profondes o
on a u c car vers les eaux peu profondes on a u ~ c est
londe dferle [breaking wave].
(x, t )
z
c
x
H
L
z = h
r
g
d ( z) = d + z
La position dune particule fluide dans le plan du mouvement est repre

sur son orbite par :
E.H.T.P.
2
ZORKANI Mohammed
: horizontalement
relativeme nt au centre de l' orbite
:
verticalem
ent
Par dfinition la clrit est de londe donne par :

( s)
T : la priode
L
c= =
(1 2) o
T k
L : la longueur d' onde (m)
Supposons que le profil en surface est donn par :
2
2
H
x t H
(1 3)
& =
(x, t ) = cos 2 = cos (k x t ) o k =
T
L
2
L T 2
Au fond pas de vitesse normale [fond tant par hypothse impermable] :
vn =
= v = 0 (1 4)
=
n au fond horizontal z z = h
Le thorme de Bernoulli Lagrange pour un coulement non
1 2
p
permanent et irrotationnel est :
u + v 2 + gz + +
= 0 (1 5)
2
t
On linarise cette condition dynamique la surface libre z = o on a
te
p = p a = C et en ngligeant les termes en vitesse car ils sont du second

ordre [quadratiques] ; et comme la pression atmosphrique constante
1
est prise comme origine des pressions alors =
(1 6)
g t z=
cette quation est applicable la surface libre et comme est faible par
1
la procdure de linarisation on peut crire : =
(1 6/)
g t z =0
si la pression en surface nest pas constante on a
p a (x, y, t )
et t z = 0 pour z = 0
g + t =
Voyons voir ce qui se passe si la pression atmosphrique varie

harmoniquement dans le temps dans une eau profondeur constante :
oscillations forces par latmosphre
p a (x, t ) = p(x ) sin t
p a (x, y, t )
1
2
On a rsoudre = 0 avec C L = t
en z = 0
t = z
Cherchons des solutions de la forme: = (x, z ) cos t + (x, z ) sin t

p 0 e z
Ainsi si p(x ) = p 0 sin x : x + (x, z ) =
sin x
g k
alors
E.H.T.P.
3
ZORKANI Mohammed
La solution gnrale est donc donne par :

p 0 e z
cos kx
kz cos kx
sin x + Aekz
cos
t
Ae
=
sin t
g k
sin kx
sin kx
On constate quil peut y avoir rsonance (amplification).

On signale que si on a pris pour la pression atmosphrique une onde
progressive : p a (x, t ) = p 0 sin(t x ) alors
p 0 e z
cos( t x )
g k
laquelle on peut ajouter nimporte quelle solution pour une pression
nulle en surface. On a rsonance si k cest dire sil existe dans
le spectre atmosphrique un mode proche ce celui propagatif dans leau.
Data : Note sur la technique mathmatique des perturbations
Cette mthode est base sur le dveloppement en srie de Taylor
autour dune position de rfrence 0 en plus de lintroduction dun
paramtre de petitesse disons ; pour cela on utilise :
2 2
(x, y, , t ) = (x, y,0, t ) +
+
+L
z z =0 2 z2
(x, z, t ) =
z =0
+ ( 0 )
=
x z
x z= x z=
2
+L
z = 0
2
2
=
+ 1 + 2 + L
x z=
x z
0
3
1
2
+ 1 + 2 + L
2
x z2
+L
z = 0
De mme pour les drives en y, z, et t :

t z = = t z = + 1 + 2 2 + L tz
0
z = 0
z = 0
= (1) + 2 (2 ) + 3 (3 ) + L
( )
()
( )
( )
= 0 + 1 + 2 2 + 3 3 + L
On pose :
D.L. en srie
(0 )
(1)
2 (2 )
3 (3 )
c = c + c + c + c + L
(0 )
(1)
2 (2 )
3 (3 )
p = p + p + p + p + L
1 2
p
Prenons par exemple la condition (1 5)
u + v 2 + gz + +
=0
2
t
la surface libre z = : p = p a 0 on a Bernoulli :
E.H.T.P.
4
ZORKANI Mohammed
2
2
= 0 Ainsi
+
+
g +
t
2
z
x
z =
z =
z =
0 = g (0 ) + (1) + 2 (2 ) + L + (t1) + 2 (t2 ) + L
(
)
+ (( ) + ( ) + L)( ( ) + ( ) + L)
1
+ (( ) + ( ) + L) ( ( ) + ( ) + L)
2
1
+ [( ( ) + ( ) + L) + (( ) + ( ) + L)( ( ) + ( ) + L) ]
2
g(( ) + ( ) + ( ) + L) + ( ) + ( ) + ( ) ( )
1
1
1
+ ( ( ) ) + ( ( ) ) + ( ( ) ) + O( ) = 0
2
2
2
1
1
tz
1
x
1
tzz
2
x
1 2
x
2
tz
2
tzz
1
t
1
xz
1 2
y
1 2
z
2
xz
+L
1
tz
(0 ) = 0 ordre zro
g(1) (x, t ) + (t1) (x,0, t ) = 0 premier ordre ( Condition de Poisson )

2
1 r
g(2 ) + (t2 ) + (1) (tz1) + (1) = 0 deuxime ordre
2
LLLLLLLLLLLLLLLLLLLLLLL
de mme pour les autre expressions. On donnera plus bas les rsultats
des calculs au 2imordre, 3imordre
Cependant la fonction potentielle sera cyclique en fonction de la position
horizontale x et du temps t et comme est la drive temporelle du
potentiel alors la technique de sparation des variables impose que :
= Z(z ) sin(k x t ) (1 7)
en la reportant dans (1 1) on obtient

Z(z ) = Ae
kz
+ Be
k z
d2 Z
dz
alors (x, z, t ) = Ae
k 2Z = 0
kz
+ Be
k z
] sin(k x t )
vrifie lquation de Laplace. Pour dterminer les 2 constantes

dintgration A et B il faut satisfaire les conditions limites(1 4) et (1 6/) :
kh
kh
2k h
v=
= k Ae
Be sin(k x t ) = 0 A = Be
z z = h
Alors
k h + z)
k (h + z )
+e
(x, z, t ) = Bekh e (
sin(k x t )
(x, z, t ) = 2Bekhch k (z + h) sin(k x t ) (1 8)
La condition (1 6/) implique que
1
H
cos (k x t )
g t z =0 2
E.H.T.P.
5
ZORKANI Mohammed
gH
Finalement on obtient pour une onde progressive :
2 chkh
g H ch k (z + h)
H
(x, z, t ) =
sin(k x t ) si cos (k x t ) (1 9)
2
ch kh
2
La vitesse verticale des particules fluides appartenant la surface libre
D
est donne par v =

qui en thorie linaire scrit :
=
+u
Dt t
x
2 B ekh =
1 2
1
or =
alors v =
v=
g t2
t
g t z =0

2
do
+g
=0
z
et comme v =
z =0
z = 0 . En y reportant notre solution (1 9 )
en
L
=
T k
gL
2h
gT 2
g T 2h
2h
th
alors c =
th
ou c =
th
ou L =
2
L
2
2 L
L
Comme en gnral donde se propage des eaux profondes vers les eaux
peu profondes sa longueur donde et sa clrit diminuent. Son profile
surfacique change, le profil de pression sur une verticale et le champ de
vitesses changent galement : Mais la priode de londe reste constante.
Classification des ondes de gravit en fonction de la profondeur
relative :
on obtient la relation de dispersion : 2 = g k th(kh ) avec c =
transitional
deep
shallow
2 h / L
1
th( 2 h / L )
0,5
n
h/L
0,05
0,5
Si la profondeur relative est plus grande que 0,5 alors th

ainsi : c 0 =
g L0
2
gT
2
et L 0 =
gT
2h
L
, Les orbites des particules sont
circulaires dont le diamtre exponentiellement dcroissant vers le

z
fond et sont proches de 0 pour 0,5

L
E.H.T.P.
6
ZORKANI Mohammed
[si par exemple T = 10s et une amplitude de 2m alors c=15,6 m/s et

L = 156 m et les particules en surface ont une vitesse
orbite circonfrentielle / priode = H / T = 0,63 m / s : on a alors
u
c
H / T H
H
=
3 4% selon la thorie de Stokes on a une
L/T
L
L 0
H
limite
L 0
max
1
7
14% donc pas de dferlement ].
Si la profondeur relative est infrieure 0,05 alors th
2h 2h
L
L
gT
et L = gh T cest la condition eau peu
2
profonde [shallow water] (exemple londe de mare) se sont des
ondes de translation [c d qui affectent uniformment la section
: onde longue] et cest galement le cas de la houle proche de la
cote [ pour T = 10s et H =2m on a L=44,3m et alors h/L<0,05 et
vitesse ~ 4,4m/s ] mais dans cette limite la thorie linaire cesse
dtre valable et il faut faire appel une thorie plus exigeante sur
les hypothses par exemple ne plus ngliger lacclration verticale
{Thorie de Boussinesq ou de Serre}.
ainsi c =
gh =
Pour 0,5 0,05 on a des conditions de transition : c et L

L
dcroissent quand h diminue cest dire vers le rivage car :

c
L
2h
=
= th
c0
L0
L
Remarque : On peut crire la relation de dispersion sous la forme
h
h h h
h
th 2 =
= f
qui facilite la rsolution numrique pour :
L
L L L0
L0
Mthode de rsolution : cest une mthode itrative
2h
L 2i + 1 = L 0 th
L 2i
2h itration
L = L 0 th
i = 0 , 1, 2 L
L
2 L 2 i + 1 + L 2i
L 2i + 2 =
3
Voici une classification des ondes (Swell) en priode :
Dsignation de londe
Intervalle de priode (sec)
Capillarit
Ultragravit
Gravit
Infragravit
Long priode
Transtidales
0 0,1
0,1 1
1 30
30 300 ~5mn
300 8,64 104 ~24h
T > 8,,64 104
E.H.T.P.
7
ZORKANI Mohammed
Rappel : Dveloppement limit de Taylor Young

1
1 2 //
//
//
+ k 2 f yy
f (a + h, b + k ) = f (a, b ) + h f x/ + k f y/ +
h f xx + 2 h k f xy
1!
2!
}](a, b) + ... +
1 p p
h f p + C1php 1k f pp 1 + ... + Cpqhp qk f pp q q + ... + k p f pp (a, b ) + Rp
x
y
x
y
y
p! x

1
+k
o Rp =
h
(p + 1)! x y
ainsi lordre 3 on a :
(p +1)
f (a + h, b + k ) avec 0 1
1
1
1

f (a + h, b + k ) = f (a, b ) + h
f (a, b ) + h
f (a, b ) +
+k
+k
1! x
2! x
y
y
3
1

h
k
f (a + h, b + k ) avec 0 1
+
3! x
y
Remarque : effet de la tension superficielle [ et les rides ]
On a tablit que dans tout le fluide en linaire que :
2
+ gz +
=0
t
1
R
1
1
2 Loi de Laplace
pint p ext = p =
+
R1 R 2
2
Or en z = on a selon la loi de Laplace : p = p a = p a
. Le
2
R
x
2
est positif cest dire la surface libre

x2
est concavit vers le haut on a une rduction de pression.
Alors si on prend : p a constante = 0 on peut crire en thorie linaire
signe moins signifie que si

pa 2
1
2
+ g +

= 0 =
+
x2
g t z = 0 g x 2
t z=
La vitesse verticale des particules fluides appartenant la surface libre
D
est donne par v =

qui en thorie linaire scrit :
=
+u
Dt t
x
1 2
3
et
comme
alors v =
+
v=
v
=
g t 2
g x 2 t
t
z
z =0
E.H.T.P.
8
ZORKANI Mohammed
alors
+g
=0
z z x2
en
z=0
t
reportons notre solution :
g H ch k (z + h)
H
(x, z, t ) =
sin(k x t ) pour cos (k x t )

2
ch kh
2
On obtient la relation de dispersion avec leffet de la tension surfacique :
3
k
2
th(kh ) la vitesse de phase est donne par c = L =
= g k +

T k
Cette relation de dispersion contient 2 termes, reprsentant les 2 sortes

de force de rappel agissantes sur le dplacement de la surface libre.
Le premier dpend de g mais pas de qui reprsente la tendance de
leau en crte revenir lquilibre sous leffet de la pesanteur (repos) ;
le deuxime terme reprsente leffet de la tension qui tend aplatir cette
surface et ainsi rduire sa courbure. La tension surfacique ne prend de
limportance et a un effet significatif que pour les courtes longueurs
donde nommes les rides. En effet les 2 termes sont du mme ordre de
k3
g
k2
grandeur quand : g k
L 2
g
Cette valeur critique de la longueur donde est de 17mm pour leau

200C
[ = 0,073 N / m , = 10 Kg / m ]
3
On peut maintenant classer les ondes :

Les rides : les ondes plus courtes que la valeur critique rendent la
surface plus courbe et la tension superficielle est alors dominante :
1
2
1
k
2
c = #
=
3
k 2
L
k

#
3 k
3
= c
#
c g =
k 2
2
On constate que c g c les rides (les trs courtes longueurs donde)
apparaissent se propager individuellement vers larrire relativement

au paquet donde.
Les ondes de gravit en eau profonde :
1
c = # = g L
k k
2
# gk
1 g 1
#
# c
c g =
k
2
k
2
E.H.T.P.
9
ZORKANI Mohammed
Les ondes de gravit en eau peu profonde :

g h
:
1 2
h
6
N.B. : les calculs prcdents rsultent du D. L. de la tangente

hyperbolique selon le cas :
alors
si kh 1
thkh 1
alors
1
2
3
(kh)5 L
si kh thkh = kh (kh ) +
2
3
15
Cinmatique de londe et les pressions : sans tension superficielle

On tablit que le potentiel de vitesse (en linaire), pour une onde
1
kh 1 2 # g h k 2 1 h2k 2 k k 3
3
progressive dont le profil de la surface libre est

donn par : (x, z, t ) =
gH
2
ch k (z + h )
ch kh
H
2
cos (k x t ) , est
sin(k x t ) on dduit alors le
champ de vitesse en thorie linaire :
H ch k (z + h)
cos(k x t )
=
u(x, z, t ) =
x T sh kh
v (x, z, t ) = = H sh k (z + h) sin(k x t )
y T sh kh
On observe que la vitesse rsulte de 3 termes :

H
Une vitesse surfacique des particules fluides
<<qui existe
T
en eau profonde>>. Proche du fond z = h la vitesse maximale
en thorie linaire hors de la couche limite est donne par
H
alors que lexcursion maximale des particules
u0 =
T sh kh
H
H
fluides proche du fond est donne par A = =
.
2
2 sh kh
Un terme hyperbolique dfinissant la dcroissance de la
vitesse en fonction de la profondeur relative.
Un terme de phase en fonction de la position x et du temps t.
Calculons maintenant les acclrations en thorie linaire :
u
u
u u 2 2 H ch k (z + h)
sin(k x t )
ax =
+u
+v
=
t
x
y t T 2 sh kh
ay =
v
v
v v 2 2 H
+u
+v
=
t
x
y t
T2
sh k (z + h)
cos(k x t )
sh
kh
E.H.T.P.
10
ZORKANI Mohammed
N. B. : On observe quentre lacclration et la vitesse existe un

dphasage : a x , a y / (u, v ) = 90 0 . Le mouvement des particules fluides
autour de leur position au repos, cest position moyenne, se calcule par :
D 2 2 H ch k (z + h)
(
)
u=
=
sin
k
x
t
D t t T 2 sh kh
D 2 H sh k (z + h)
(
)
=
v=
cos
k
x
t
Dt t
T 2 sh kh
H ch k (z + h)
sin(k x t )
= udt =
2 sh kh
= vdt =
H sh k (z + h)
cos(k x t )
2 sh kh
Lorbite des particules fluides est donne par :
2
2
2
2
=1
a
b
H ch k (z + h)
H sh k (z + h)
et b =
o on a pos : a =
2 sh kh
2 sh kh
H
Ce sont donc des ellipses. On constate que
est le rayon des orbites
2
des particules en surface en eau profonde. Quand londe des eaux
profondes vers les eaux peu profondes [en passant par les eaux
intermdiaires] les orbites subissent les transformations suivantes :
0,5 H
kz
L
2
M
eau profonde
eau faible profondeur
ch k (z + h) sh k (z + h)

ekz quand z
Eau profonde :
sh kh sh kh
ch k (z + h) 1
sh k (z + h)
z

Eau peu profonde :
1 +
et
sh kh kh
sh kh h
Calculons les vitesses en eau peu profonde :
H
z
H g
u(x, z, t ) =
cos(k x t ) & v (x, z, t ) =
1 + sin(k x t )
2 h
T h
On peut galement dterminer les dplacements comme avant.
E.H.T.P.
11
ZORKANI Mohammed
Reportons notre solution dans le thorme de Bernoulli Lagrange

on obtient alors le profil de pression :
gH ch k (h + z )
p = gz +
cos(kx t )
2
ch kh
pression statique
pression dynamique due lacclration verticale

(rgime non permanent)
S.W.L.
gHchk h + z
2chkh
gz
Surpressions et sous pressions sont donc

p
z=
= en surface = g
H
2
et p z = h = sur le fond =
gH
2ch kh
4 ) Energie de londe et Puissance :

Lnergie cintique par unit de largueur et pour une longueur donde est
gH2L
L 0 1
2
2
E c = 0 h dxdz u + v E c =
16
2
Si nous retranchons lnergie potentielle dune masse au repos de
lnergie potentielle totale du volume ondul (c d avec surface
dforme) {voir figure} on obtient lnergie potentielle due seulement
londe en propagation ; ainsi lnergie potentielle par unit de largueur et
h +
h
pour une longueur donde est : Ep = 0L g(h + )
dx gL
2
2
Ep =
gH L
16
2
v
dz
dx
E.H.T.P.
12
ZORKANI Mohammed
On constate quipartition dnergie : Ep =
gH2L
= E c . Lnergie totale
16
g H2 L
unit
est donc : E = E c + Ep =
(Joule / mtre de l argeur )

8
Prenons, par exemple, une onde qui se propage travers dune
structure poreuse et supposons que la profondeur est la mme des 2
cots de louvrage alors on aura la mme longueur donde de chaque
gT 2 2h
th
cot car L =
; par consquent on aura une rduction en
2 L
amplitude car la structure induit la rflexion dune partie de lnergie
incidente et une autre partie dissipe dans louvrage poreux : daprs le
principe de conservation dnergie on a :
g H2 L
g HR2 L
g H2T L
E=
=
+
+ E dissipe
8
8
8
g H2 L
g HR2 L
g H2T L
E dissipe =
8
8
8
2
2
2
g H L HR HT
E dissipe =
1 2 2
8
H H
or par dfinition les coefficients de rflexion et de transmission (dont le

carr sont les facteurs associs) sont respectivement :
alors
H H
H H
H
H
R R = R R et T T = T T
H
H
H H
H H
g H2 L
2
2
E dissipe = E 1 R T
o E =
8
On dfinit le coefficient dabsorption par : A = 1 R 2 T 2
de sorte
que E dissipe = E A 2 il en rsulte lgalit : R 2 + T 2 + A 2 = 1

Une rduction de ~50% de lnergie correspond une rduction de
lamplitude de ~71% car E H2 . On peut facilement atteindre le
E.H.T.P.
13
ZORKANI Mohammed
coefficient de rflexion R et de transmission T exprimentalement et on

en dduit par consquent le coefficient dabsorption A de louvrage .
Lnergie tant variable dun point un autre sur une longueur donde on
parle alors de lnergie moyenne par unit de surface :
E E c + Ep
g H2
unit
E= =
=
(Joule / surface horizontale )

L
L
8
qui reprsente la densit dnergie ou lnergie spcifique.
N.B. : rflexion donde par un absorbeur donde vertical
Par Madsen [J. WaterWays3(74) et Coastal Engineering7(83) ]
ar
ai
at
monticule de
gravats
x
w
Lquation qui gouverne lcoulement hors de la structure poreuse est :
ig
g
2
xx +
= xx + k 2 = 0 et U(x, t ) = x = m
h
gh
Les quations qui gouvernent lcoulement dans labsorbeur sont :
n t + hUx = 0 conservation de la masse
1
Ut + g x + U( + U ) = 0 conservation de la quantit de mouvement
n
o et tiennent compte respectivement de la perte par frottement en
rgime laminaire et turbulent et n la porosit de la structure. On pose :
( + U ) U f U linarisation du frottement (approximation)
n
On cherche des solutions priodiques de frquence de la forme :
= Re (x )e
it
] et U = Re[v(x)e ]
it
En les reportant dans les quations du mouvement et en liminant U on

2
(1 if ) = 0 et v = gn 1 x
obtient :
xx +
f +i
gh
La solution gnrale, donne par Madsen et White (1976) pour un
coulement dans une structure poreuse, est :
E.H.T.P.
14
ZORKANI Mohammed
= Re a1e ix + a 2 eix eit

:0 x w
1 if
g
= 1 if
eit : 0 x w
U = Re a1e ix a 2 eix
gh
h
A lextrieur de la structure poreuse SWW quations donnent :
i t kx )
i t + kx )
= Re ai e (
+ ar e (
:x 0

o k =
gh
g i(t kx )
i t + kx )
ai e
ar e (
:
x
0
U = Re

h
ai & ar sont les amplitudes de londe incidente et londe rflchie.
Nos inconnues sont les amplitudes complexes a1, a 2 & ar . Elles peuvent
tre dtermines par application des conditions aux limites : au niveau
de la face frontale de labsorbeur et sur le mur vertical en arrire.
Lamplitude a 2 peut tre limine en utilisant le faite que la vitesse est
nulle au mur (x = w ) alors :
a 2 = a1e
2 i x
Les amplitudes a1 & ar sont dtermines en admettant la continuit de la

pression (donc de llvation de la surface libre) et la continuit de la
masse ( donc de la vitesse) la face frontale de labsorbeur (x = 0 ) on
ai + ar = ar 1 + e 2ikw
ar
(
1 ) + (1 + ) e 2ikw
=
obtient ainsi :
2ikw
a
(1 + ) + (1 + ) e 2ikw
i
ai ar = a t 1 e
a
Le coefficient de rflexion est donn par : R = r = R(f , n , kw )
ai
R
1
f ~ 10
n ~ 0,5
0,5
n ~ 0,95
kw
1
ar 1
=
et pas doscillation de R.
kw ai
1+
Pour un absorbeur long on a : Limite
E.H.T.P.
15
ZORKANI Mohammed
On note que labsorbeur amortie fortement les ondes tel que kw est
grand c d les courtes longueurs donde. Un absorbeur de forte
permabilit n amortie mieux les ondes : faible rflexion. Notons que
pour un court absorbeur la rflexion est presque totale.
La vitesse dcoulement dans labsorbeur est obtenue en dterminant a t
par le systme et en reportant dans la solution gnrale ; soit :
i x 2 w )
2 e ix e (
eit : 0 x w
U = ai
Re
2ix
h
(1 + ) + (1 ) e
La dtermination des coefficients de frottement peut se faire par les

3
(
(1 n)
1 n)
formules empiriques dEngelund : = 0
&
0
n3 d
n 2 d2
o d la taille des grains, la viscosit cinmatique et 0 & 0 sont des
constantes qui tiennent de la forme des particules 0 ~ 1000 & 0 ~ 2,8
et qui augmentent avec lirrgularit des particules :
0 ~ 780 1500ou plus & 0 ~ 1,8 3,6 ou plus
La puissance de londe est lnergie de londe par unit de temps se
propageant dans la direction de londe. Cette puissance peut scrire
comme le produit de la force agissante sur un plan vertical normal la
direction de propagation par la vitesse des particules fluides traversant
1 0
ce plan. On a donc : P = 0T (p + gz ) udtdz selon BernoulliLagrange
T h
1
E 2x + 2y + 2z + gz dxdydz = (p + t ) dxdydz
D 2
D
do on a le flux dnergie cest dire la puissance or D
bouge appliquons le thorme de transport de Reynolds :
DE
= x xt + y yt + z zt dxdydz (p + t ) v n ds
P=
Dt
D (t )
= grad grad t dxdydz (p + t ) v n ds

D(t )
= t n (p + t ) v n ds
DE
P=
= F = [ t ( n v n ) pv n ] ds
Dt
S
On signale quon a utilis la formule de Green pour effectuer le calcul
r
r
r
r

pour nous
2
2
nds
avec 0
S
t
E.H.T.P.
16
ZORKANI Mohammed
Appliquons cela notre cas: v n = 0 (la surface gomtrique S est fixe)

= u(x, z, t )
dE
alors P =
= t n dS , t = (p + gz ) n
dt
dS = dz dx
Notons que la force mise en jeu est celle dynamique [p + gz ] ; en
2 kh E 1
2 kh
gH2L
1 +
= 1 +
effectuant le calcul alors : P =

16T sh 2kh T 2 sh 2kh
Posons :
n=
2k h
1
nE
1 +
on a alors P =
2 sh 2k h
T
g H2 L
(Joule / mtre de l arg eur )
avec E = E c + Ep =
8
Quand une onde se propage des eaux profondes [du large] vers les
eaux peu profondes [la cte] : lnergie par unit de temps (la puissance)
en un point le long de son chemin de propagation doit tre gale celle
cot e
te
nE nE
P=
=
=C
T 1 T 2
1
2
Le l arge
en un autre point plus loin augmente de lnergie rflchie et de celle

dissipe par unit de temps entre ces 2 points. Si on nglige la rflexion
et la dissipation on peut alors crire : P = Constante ; donc quand londe
se propage du large vers la cote son nergie E dcrot inversement
proportionnel n car la priode T reste constante.
Maintenant construisons les orthogonales aux lignes de crte, la
puissance contenue entre 2 orthogonales conscutives doit tre
constante si on nglige bien entendu la dissipation et la rflexion, on
dsignera par B leur espacement : T (la priode) est constante
On dfinit ainsi les coefficients suivant :
plage
breaking
B
te
nE nE
B
= B
=C
T 1 T 2
d' o en y reportant E on obtient :
H1
n2L 2 B 2
=
H2
n1L1 B1
E.H.T.P.
17
ZORKANI Mohammed
K S = coefficient de Shoaling (profondeur )[gonflement ] =

K R = coefficient de refraction =
B2
B1
n2L 2
n1L1
H1
= K S KR
H2
Si pas de rfraction alors K R = 1 cest le cas o les lignes de crtes sont

parallles aux linges bathymetriques dgale profondeur.
Si on veut tenir compte des pertes par : infiltration, frottement ou
dferlement il faut crire : H1 H21 = K S K R K F .
Vitesse de groupe :
(c,L)
(c + dc, L + dL )
x
SWL
addition
c + dc
des
L + dL
cg
2 ondes
train d' ondes : un groupe

Quand les ondes du train donde se propagent le long dun canal houle
les ondes au devant du groupe diminuent en amplitude et de nouveaux
ondes apparaissent en arrire, de ce fait le nombre donde augmente.
Ce qui signifie que le groupe se propage plus lentement que les
composantes individuelles du paquet.
Lexplication de ce phnomne se trouve dans la faon dont une onde
se propage et que seulement une fraction (n) de son nergie qui se
transport avec. Chaque onde laisse de lnergie derrire elle,
relativement au groupe, une nouvelle onde apparat chaque T secondes
et gagne progressivement de lnergie dans le temps. Comme lnergie
dans le groupe doit rester constante [en ngligeant la rflexion et la
dissipation] lamplitude moyenne du groupe (enveloppe) doit continuer
diminuer car le nombre donde augmente avec le temps.
Une consquence et une application pratique de la notion de vitesse de
groupe [ moins quon ne sintresse chaque onde caractrise par sa
clrit] est de prvoir le temps larriv du paquet donde en un point
donn en utilisant la vitesse du groupe. Pour dterminer la vitesse du
groupe prenons le cas de 2 ondes monochromatiques de priode
dL
,
diffrentes (voir figure) : pour parcourir dL , il faut dt avec dt =
dc
pendant dt le train donde avance de dx :
E.H.T.P.
18
ZORKANI Mohammed
(L + dL ) + L cdt L car x = c t x alors

(c + dc ) + c
dx =
dt
g
0
2
2
dx cdt L
L
cg =
=
=c
dc
dt
dt
dt
cg = c L
dL
dL
dt =
dc
On observe que cest du fait que c est variable avec L que le paquet
donde possde une clrit de groupe diffrente de la vitesse de phase
de chacune des composantes.
2k h
1
c
2k h
:
= n c avec n = 1 +
On obtient alors c g = 1 +
2 sh 2k h
2 sh 2k h
n = 0,5 en eau profonde
0,5 n 1 avec
n = 1 en eau peu profonde
Ainsi lnergie du groupe se propage vitesse du groupe.
Autre mthode : notion de vitesse de groupe
On peut comprendre plus clairement cette notion de vitesse de groupe
en superposant deux ondes monochromatiques de mme amplitude
mais de frquence trs voisine ( , k ) et ( + , k + k ) :
= a sin[kx t ] + a sin[(k + k )x ( + )t ]
1
1
1
= 2a sin k + k x + t cos k x t
2
2
2
2
4
=
1
k
k
2
Ainsi lenveloppe [lamplitude du groupe] se propage avec la clrit :

soit c g
=
limite
k
k 0 k k
t fixe :
2
2
~
k + k
k
c
kc
c
c L
= c L
=
=c +k
=c +k
L k
L
k
k
k
cest ce rsultat quon a tablit prcdemment.
On peut retrouver ce rsultat par une approche plus lgante en
calculons le flux dnergie F pendant un temps T qui traverse une
surface S fixe dans lespace :
t+T
n
F = t
d s d t
n
t S
Soit
cg =
E.H.T.P.
19
ZORKANI Mohammed
En supposant que le mouvement est la superposition de 2 ondes

harmoniques dans le temps, soit :
(x, y, z, t ) = 1(x, y, z ) cos t + 2 (x, y, z ) sin t
introduisons cette solution dans lexpression du flux dnergie qui
2
(T tant la priode doscillation) on obtient :
traverse S avec =
T
2
1
ds
1
F = 2
n
n
S
ainsi le flux dnergie est nul si le mouvement est stationnaire : ce qui
nest pas surprenons car le transfert dnergie naura lieu que si le
mouvement est progressif.
Maintenant si 1 et 2 sont harmonique et si la surface S est ferme fixe
dans le fluide qui enveloppe un domaine D, la formule de Green nous
donne :
2
1
d s = 2 2 1 2 2 1 dxdydz
1
F = 2
n
n
S
D
Si 1 et 2 nont pas de singularits dans D (sources o puits) alors le
flux dnergie est nul car 1 et 2 sont harmoniques. Calculons dans le
cas dune onde progressive la vitesse avec laquelle le plus dnergie se
propage ; prenons alors = Achk (z + h) cos(kx + t + ) alors :
t+
t+T
n
F = t
d s d t = A 2k h~ 0 ch2k (z + h)dz sin2 (kx + t ) dt
n
t S
t
ainsi le flux dnergie moyen par unit de temps est :

F A 2kh sh 2kh
1 +
Fav. = =
T
4
2
k
h
F A 22 2
comme on a = gk thkh et c = alors Fav. = =
ch kh . c g
k
T
2g
1
2k h
o c g a les dimensions dune vitesse donne par : c g = c 1 +

2 sh2kh
Calculons maintenant lnergie moyenne stocke dans leau due au
mouvement ondulatoire par rapport la direction de propagation. Nous
avons vu que lnergie stocke dans D est donne par :
1
E 2x + 2y + 2z + gz dxdydz
D 2
Calculons cette nergie pour une largeur unit sur une longueur donde
nimporte quel instant :
2
E.H.T.P.
20
ZORKANI Mohammed
1
E E 0 = k 2 h 0L A 2 sh2k (z + h) cos 2 (k x + t + ) +
2
1 2 2
A ch k (z + h) sin2 (k x + t + ) dxdz
2
+ h 0L gz dxdz
o E 0 est lnergie potentielle de leau de profondeur h quand elle est au
repos. En ngligeant les termes dordre lev en amplitude, nous
obtenons ainsi lnergie due uniquement la propagation donde
progressive entre 2 plans distants dune longueur donde L :
A 22
E E0 =
Lch 2kh
2g
ainsi lnergie moyenne E av. dans le fluide par unit de longueur dans la
direction de propagation x qui rsulte du mouvement ondulatoire de leau
A 22
ch2kh
est donne par : E av =
2g
Nous constatons alors que :
Fav flux d' nergie par unit de temps
travers un plan vertical
avec
Fav = E av c g
E av nergie moyenne par unit de temps
dans la direction de propagatio n

Ainsi sous lhypothse que pas dnergie cre ou dtruite dans le fluide,
celle ci est transmise dans la direction de propagation de londe la
vitesse de groupe c g .
Relation de dispersion des ondes linterface de 2 couches dun
fluide de masse volumique diffrente : ondes internes
z
air
h / et /
x
h et
On cherche des solutions pour les 2 couches de fluide sous la forme :
= A chk (z + h) cos(k x t )
/ = B chk z h / cos(k x t )
E.H.T.P.
21
ZORKANI Mohammed
Au niveau de linterface entre les 2 fluides la pression [dynamique] et la

vitesse doivent tre continues; il en rsulte daprs la continuit de la
pression que :
/
/ /
1

/
/
g +
= g +
=
/
t
t
=
t
g /
/ 2 /
2
t2
t 2
Alors que la continuit de la vitesse est :

/
D
=
mais comme par dfinition v =
on obtient :
=
z
z
Dt
t
z
g /
) z =
2 /
t2
t2
En y reportant nos 2 solutions on obtient un systme en A et B, qui a une
solution non triviale quand son dterminant est nul alors on a :
=
2
g k /
coth kh + / coth kh /
g
/

c = =
k coth kh + / coth kh /
k
2
Note : si nous avons tenu compte de la capillarit pour les ondes

internes courtes on aura comme relation de dispersion :
gk / + k 3
g
g k / + k

c = =
=
k coth kh + / coth kh /
coth kh + / coth kh /
k
N.B. :
Si kh 1 et kh / 1 cest dire en eau profonde on a alors :
2
= gk
2
g /
c =
k + /
2
+ /
tout se passe comme si le fluide suprieur rduit g g = g

/
/
+ /
car
g
en eau profonde.
k
On remarque que si / c d que le fluide suprieur a une
densit plus grande devient imaginaire : il y a donc instabilit.
Si kh 1 et kh / 1 cest dire en eau peu profonde on a alors :
pour un fluide homogne on a c 2 =
2 = k 2
g / k h/
h / + / h
E.H.T.P.
22
ZORKANI Mohammed
/
Si kh 1 et kh / 1 on a alors : 2 = k 2 g h /
Pour ce qui est de ltude dynamique et cinmatique des ondes internes

voir le chapitre sur locanographie physique.
5) Transport de masse et projection deau sur le rivage (wave setup):
La thorie de faible amplitude (linaire) prvoit pour les particules fluides
des trajectoires fermes, ainsi il en rsulte quil ny a pas de transport de
masse. Cependant dans la ralit et surtout en eau peu profonde on
observe un transport de masse qui rsulte du fait que les trajectoires des
particules fluides ne sont pas fermes donc les particules avance
pendant le parcourt de chaque orbite : la vitesse de transport augmente
quand la profondeur relative diminue. La vitesse de transport massique
est plus faible que la vitesse des particules mais significative pour induire
une remonte deau le long du rivage et contribuer par consquence
largement au transport des sdiments vers large proche du fond (par le
courant de retour) aprs leur mise en suspension par la turbulence car le
dferlement projette deau vers la cte dans une couche de surface et
S ww
vertical
Hb
SWL
zone de lavage
(Wash zone)
courant de retour
comme on a conservation de la masse il en rsulte un courant de retour.
En se basant sur des tudes exprimentales en laboratoire Saville
(1961) pour des ondes dferlantes sur une plage a tablit une quation
donnant la remonte deau sur le plage (wave setup at the shore) S ww
qui a t propose ultrieurement (1973) par U. S. Army Coastal
Engineering Research Center. Si Hb est lamplitude au dferlement
dans la zone du ressac (cest un retour violant des vagues)(surf zone) :
Hb
S ww = 0,19 1 2,82
Hb
gT 2
Pratiquement londe se projette en moyenne sur 15% de Hb .
On prsentera une modlisation mathmatique de ce genre de problme
pour dterminer la descente (wave setdown) et la remonte deau (wave
setup) sur un rivage et la courantologie (transport de masse induit par la
houle : circulation ctire lchelle de la houle), ce qui permettra une
meilleur reprsentation de la courantologie marine et une bonne
modlisation du transport de sdiments et une tude de la dynamique de
E.H.T.P.
23
ZORKANI Mohammed
la gomorphologie des cotes sableuses . A ce niveau on tient signaler

que les effets non linaires sont responsables de lexcitation
dharmoniques (de faibles priodes) qui modulent le profil de la surface
libre ainsi que le champ de vitesse qui en rsulte auquel sajoute la
rflexion induite au moins par la pente du fond : ce qui module
globalement la morphologie et installe une barre sous marine qui
constitue une protection naturelle de la cote (il se peut quil sexcite
galement des ondes de coin \\ adge waves // dont lamplitude dcrot
exponentiellement vers le large : ce qui modifie galement la
courantologie du littoral.
Des ondes dincidence normale la cte qui sont fortement rflchies
par la ligne du rivage sont instables la perturbation induite due aux
ondes de coin. Ces ondes de coin extraient leur nergie de londe
incidente par le billet dinteractions non linaires (Galvin 1965,Guza
et Inman 1975) le potentiel de vitesse de cette onde subharmonique est :
Ag k e x
e =
2
e
cos k e y sin(e t + ) avec :
e =
e
2 = g k tg
e
e
x
y
Les valeurs successives duprushes donnent une valeur approximative
R R1
A 2
de lamplitude :
tg
2
o R1 et R 2 sont les intrusions (horizontales) successives et maximales

de londe incidente sur le rivage (voir complment vers la fin de ce chapitre).
6) Rflexion donde : Clapotis
Quand une onde rencontre un changement dans les conditions limites
(comme un changement de profondeur deau, une convergence ou une
divergence dans un canal houle, un obstacle submerg ou flottant la
surface libre, un mur verticale ou en talus ) il en rsulte une rflexion
partielle ou totale de lnergie incidente.
On superpose 2 houles sinusodales progressives de mme
caractristiques mais qui se propagent en sens inverse, on obtient le
Clapotis qui est une onde stationnaire. Prenons le cas obstacle
vertical inlastique et lisse (pas de frottement) : londe incidente sera
alors compltement rflchie cest dire lamplitude rflchie est
gale celle incidente : la superposition de ces 2 ondes donne une onde
purement stationnaire avec bien entendu des nuds parfaits et des
E.H.T.P.
24
ZORKANI Mohammed
ventres maximaux. On dmontre quaux nuds lenveloppe a (1 R ) H

ER HR
=
est le coefficient de rflexion.
E
H
Les orbites en rflexion totale sont aplaties avec un retour sur la mme
trajectoire courbe. Les trajectoires sont verticales sous les ventres et
horizontales sous les nuds. Londulation de la surface libre rsulte de
la superposition des 2 ondes incidente et rflchie (on a le droit dutiliser
le thorme de superposition car on tudie les ondes linaires) :
H
H
(x, t ) = cos(k x t ) + cos(k x + t ) = H cos k x cos t
2
2
Lamplitude de londe stationnaire est le double de celle incidente : 2 H .
On verra que ceci se traduit par une surpression sur louvrage:
cest pour cela quon construit plus des ouvrages qui ont un faible
coefficient de rflexion (en talus) et qui dissipent mieux lnergie
incidente (poreux : utilisation des blocs pour revtir la carapace avec
distribution approprie cet effet ).
et aux ventres (1 + R ) H o R =
ventre
nud
Enveloppe de la surface : onde stationaire pure R =

ventre A (1 + R ) H
HR
=1
H
nud A (1 R ) H
Enveloppe de la surface : onde partiellem ent stationaire R =
HR
H
E.H.T.P.
25
ZORKANI Mohammed
Remarque : Formule de Healy technique exprimentale

La mesure par une sonde linaire cest dire dont la rponse est
proportionnelle loscillation de la surface libre nous permet de
dterminer exprimentalement le coefficient de rflexion R car :
A ventre = (1 + R ) H
1+ R
1 TOS
R=
T.O.S.(Taux Of Stationary wave ) =
1 R
1 + TOS
A nud = (1 R ) H
La mesure de TOS (en dterminant lenveloppe de londe) permet donc
de dterminer le coefficient de rflexion de lobstacle : Formule de Healy
De manire similaire on dtermine le potentiel de vitesse de londe
rsultante par la superposition du potentiel incident et rflchie :
R = 1
g H ch k (z + h)
(x, z, t ) =
cos(k x ) sin( t ) si
H
ch kh
(x, t ) = 2 cos(k x t )
La fonction de courant est : (x, z, t ) =
gH
sh k (z + h )
ch kh
sin(k x ) sin( t )
En dduit alors de lquation linarise de pression que :

ch kh(h + z )
p = gz + gH
cos(k x ) cos( t )
ch kh
pression statique
pression dynamique due lacclration verticale

(rgime non permanent)
On peut ainsi dterminer le profil de pression au mur parfaitement

rflchissant R = 1 (qui est un ventre en x = 0 ). Dans la raliser du faite
que le mur est rugueux on a une dissipation dnergie (mme faible) qui
donne un coefficient de rflexion plus petit que 1 : lamplitude de londe
rflchie est lgrement infrieure celle incidente.
Connaissant le potentiel de vitesse 2D on peut en dduire le champ des
r
r r
vitesses par : v = grad =
i+
k
x
z
Clapotis gaufre :
Une houle se propageant dans une direction faisant un angle ( + )
H
avec Ox sur un mur est : 1 = cos[ t + k (x cos + y sin )]
2
une rflexion sur la paroi Oy provoque une houle se propageant dans la
H
direction , soit : 2 = cos[ t k (x cos y sin )]
2
y
incidente
rflchie
E.H.T.P.
26
ZORKANI Mohammed
Leur somme est = 1 + 2 = H cos( t + k y sin ) cos(k x cos )

Il y a donc des maxima et des annulations damplitude sur des parallles
Oy : k x cos = n et k x cos = (2n + 1) respectivement. Il y a une

2
ligne de nuds et de ventres mais il ny a pas simultanment un
maximum ou une annulation ; ceux ci se propage paralllement Oy
c
la clrit
(vers Oy ngatif si est positif).
sin
Rflexion normale sur un talus :
h
b
Pour les houles de faible cambrure, il est constat que la rflexion tait
b 1
peu prs totale pour . Quand la pente du talus devient plus faible, il
L 4
se produit un dferlement qui dissipe une part apprciable de lnergie
incidente et lnergie rflchie est proportionnelle au carr de lamplitude
b 1
b 3
est voisin de 0,1 quand = alors que lorsque = lnergie rflchie
L 2
L 4
nest que de 0,04 0,05 celle incidente.
Pour les houles de cambrure plus grande, lnergie de londe rflchie
est sensiblement plus faible que celle des houles de cambrure plus faible
car les houles de forte cambrure dferle plus facilement.
LIngnieur espagnol Iribarren propose une formule qui donne la pente
minimale de lobstacle rflchissant une houle de priode T. Pour des
valeurs plus faibles la houle dferle :
2 H
=
Iribarren
T 2g
La cambrure limite que peut avoir la houle sans dferle sur la pente est
2 sin2
H
=
M. Miche

L max
ainsi plus une houle est longue plus elle se rflchie facilement sur une
pente donne : En particulier londe de marre se rflchie sur tous les
rivages, mme ceux pente trs douce. Ces rsultats sont pour un fond
lisse et monolithique :
la rflexion change avec la rugosit et la porosit du fond et de la pente :
E.H.T.P.
27
ZORKANI Mohammed
K
K
2
2 L b sin
H Hb
si
:R =K
o en radians
2
L Lb
Hb
K
K
= 1 asphalte, bton plat

= 0,9 bton rugueux
= 0,6 remblai
= 0,5 enrochement artificiel
7) Le dferlement des ondes :

Nous avons vue que les particules fluides appartenants la crte dune
onde ont une vitesse plus faible que la clrit de londe.
Comme en eau profonde la vitesse des particules en surface est
H
donc une augmentation
proportionnelle lamplitude de londe u0
T
de lamplitude correspond laugmentation de la vitesse de ces
gT
particules qui tendent vers la clrit de londe c 0 =
ainsi dans cette
2
2
2
L
1
gT
2 H
H
u
limite donde =
1 H ~
= 0 ~ = 31,8% devient
2
2
c 0 gT
L 0
2
instable et dferle. Miche en 1944 avait dtermin les conditions au
dferlement sur fond horizontal, il propose :
1 2h
h
H
= th
c = cambrure critique (ou limite ) =
pour : 0,11 0,36
L
L max 7 L
Soit pour 0,066 hL010,35 et si pas de rflexion ni dissipation dnergie

de londe on a : 01 = thkhb donc cest les houles de cambrure faible qui
gonflent le plus avant de dferler.
Ainsi cette formule admet les 2 limites :
1
gT 2
H
En eau profonde : = 14% ; L 0 =
2
L 0 7
1 2h
H
H
= 0,9
=

En eau peu profonde :
h max
L max 7 L
Miche en 1951 proposait galement pour la cambrure limite, en eau
profonde sur un fond en pente avec rflexion totale, lexpression :
2
H
=

L critique
0,5
sin2
= pente du fond
Des ondes de types Stockes en mer trs profonde : L = 0,318 th
2h
L
E.H.T.P.
28
ZORKANI Mohammed
On classe le dferlement sur une plage en 3 catgories :

Dversant : apparition dcume mi hauteur. La crte instable
scoule sur la face avant (le front de londe est mousseux).
Plongeant : formation de rouleaux. La crte se courbe et tombe ce
qui induit un emprisonnement dair et un bruit sourd.
Gonflant : crte presque intacte mais pas sa base et le front avant.
b ~ 1,11
SWL
b ~ 0,84
0,5 0 3,3
0 0,5
glissant (spilling)
ou Deversant
b ~ 1,25
plongeant (plunging )
c=
g (h + )
0 =
0 3,3
frontal (surging)
ou gonflant
tg
=
H0 L 0
1
tg
2 H0 gT 2
La zone du dferlement en eau peu profonde est caractrise par une

saturation en nergie de la houle. Les mesures tant en laboratoire quen
nature montre un contrle de lamplitude de la houle a par la profondeur
locale suivant une relation linaire :
+a
max min = H = 2 a = h
La constante , initialement introduite par Mac Cowan (1891 94) dans

ltude thorique de londe solitaire ( = 0,78 ) .
En eau peu profonde (kh 1) Miche (1944) en supposant que londe de
Stokes (sinusodale) dferle sur un fond horizontal quand la vitesse de la
particule fluide en crte est gale sa clrit, il trouva : ( = 0,88 ) .
Comme londe de Stokes est symtrique cette valeur reprsente le
dferlement Spilling.
Ultrieurement beaucoup dauteurs [Galvin et Eagleson (1965), Iverson
(1962), Sverdrup et Munk (1946) ] proposent la mme valeur de qui
est observe en laboratoire (houle monochromatique).
Mais dans lapproche statistique base sur lamplitude quadratique
moyenne A rms conduit une valeur de sensiblement plus faible :
0,3 rms 0,5 Nelson (1983).
E.H.T.P.
29
ZORKANI Mohammed
Sur un fond de pente on dispose du critre :

tg
= paramtre de similarit de Battjes =
2,3
H
L0
Hb eb
hs
at
hM
g T2
o L 0 =
2
niveau d' eau au repos
h t = cte du creux
niveau d' eau moyen
a t = distance entre le niveau

moyen et le creux
h M = cte du niveau moyen
ht
h s = cte du niveau d' eau au repos

e b = hM h s
fond
0,04 pour pente de 0,05 0,10

eb hM h s
=
= surlvation =
Galvin
Hb
Hb
0,08 pour pente de 0,20
h
h
On va dsigner par : hb = hM et b = bM = M = b
Hb Hb
Pour le dferlement sur une plage de pente m Weggel (1972) propose :
=
1
3
Hb
H0/
F(m)
H0/
1
3
H
F(m) b/
H
0 + G(m) o
(
)
+Gm =
1
Hb 3

L
L0
0
1
D1(1 + m) D 2 1,715 0,185 e 28m
(
)
D
=
0
,
01
+
0
,
5
m
3
G(m) =
1
o
D1 D 2
F(m) = D1[1 + m G(m)]
D1 = 0,01 0,01e 28m 3
qui est en accord avec les courbes dIversen (U. S. Army Coastal Eng.)
Lamplitude maximale de la houle au dferlement :
) (
a(m) = 23,75 1 e 19m en s 2 / m

H
Hb
o
= b(m) a(m) b2

1,56
sans dim.
hb
gT
b(m) =
19,5m
1
e
+
E.H.T.P.
30
ZORKANI Mohammed
Proposition :
Pour comprendre la physique du dferlement il est intressant de se
rfrer la thorie de Boussineq pour le calcul de la clrit cest 3 h2 2
gh
2
gh 1 +
+
c =
dire :
2
2
h
3
3 h2 2
x
1
2
2 h 3 x
quil faut comparer la vitesse des particules fluides en crte de donde
en thorie non linaire (Stokes 2e ordre) : critre cinmatique.
Nous conseillons pour prvoir les meilleurs conditions de dferlement
dutiliser les figures (7 1 et 7 2 ) qui sont bases sur des rsultats
exprimentaux. Selon Kana (1979) on a :
Spilling brea kers : = 0,55 0,65
Transition al brea kers : = 0,65 0,75

Plunging brea kers : = 0,75 0,90
En se basant sur des rsultats exprimentaux de laboratoire Weggel

H
1,56
43,8
(1972) propose : =
1 e 19th b
L0
2
1 + e 19,5th
En se donnant la profondeur deau et la pente de la plage : lamplitude

de londe au dferlement Hb se calcule par Fig(7 1) alors que la
profondeur au dferlement hb = db par la Fig(7 2). On signale que si
londe se rfracte sur la pente de la plage, lamplitude dune onde
hypothtique non rfracte est value par : H0/ = K R H0
Munk a montr en 1949 que les conditions de dferlement sont lies la
H
1
cambrure en eau profonde : b/ =
.
13
H0
H0/
3,3
L 0
Quelle est la forme limite en surface de la houle irrotationnelle en
profondeur deau infinie ? : Stokes 1880
vr
r
e
vr
r
er
r
g
E.H.T.P.
31
ZORKANI Mohammed
A mesure que lamplitude dune vague crot, les crtes deviennent de

plus en plus aigues. Donnons londe, une clrit c : londe progressive
est remplace par un mouvement densemble. Nous savons quune
Dp
particule de la surface libre ne la quitte pas car
0 . Elle se comporte
Dt
alors comme un mobile glissant sur un plan inclin. Au sommet de la
vague en O, la vitesse est nulle. La vitesse v r de la particule superficielle
est # r puisquil sagit dun glissement sans frottement sur un plan
inclin dans le champ de pesanteur : en effet
selon la loi de Newton on a :
d vr
1 te
te
te
m
= g cos = C v r = C t r = C t 2 t r
dt
2
Recherchons, au voisinage de la crte, un potentiel donn par :
2 2
+
=0
x 2 z2
vrifiant les conditions de symtrie imposes par lexistence de cette
crte.
Lquation de continuit permet lintroduction de la fonction de courant
telle que :

=
u
=
z
x
2 2
et la fonction de courant vrifie

+
=0
2
2
z
w =
=
x
z
te
te
On peut montrer que = C 1 et = C 2 sont orthogonaux.

Cherchons un potentiel de vitesse complexe du type :
= A r n cos n
n
n
n in
+ i = A (x + i z ) = A z = A r e
= A r n sin n
ce choix est justifi en raison des conditions de symtrie.
La surface libre est ligne de courant dont y est constante. Comme elle
sannule en r = 0 ; il en rsulte alors que : = 0 la surface libre
Au voisinage de la crte la vitesse est donne par :
1
v r = (2gr cos l )2 =
= nr n 1 sin n l
r
( v r 0 pour 0 et v r 0 pour 0 en raison du mouvement).
3
La valeur de lexposant se dduit de lgalit prcdente : n = .
2
E.H.T.P.
32
ZORKANI Mohammed
Ainsi prs de la crte et la surface :
3
Ar 2
3
cos l = 0
2
= 60 0
3
Langle limite (point anguleux) de la crte de la vague de Stokes est
donc de 120 0 alors que pour une onde solitaire l = 45 0 soit 90 0 .
La vitesse de la particule concidant avec la crte est nulle avec la
systme daxes considr. En eau profonde, la particule fluide concidant
avec le point anguleux est anim dune vitesse gale la clrit de
londe. Cest le critre de dferlement.
2a H
= = 14% . La houle
La cambrure limite pour la houle de Stokes est
L
L
progressive ne peut avoir une cambrure suprieure 14%, au del, la
lame dferle.
On va maintenant discuter leffet dun corps flottant sur la propagation
dune houle en eau peu profonde. Seulement le cas 2D sera prsent
(toutes les quantits hydrodynamiques sont indpendantes de la
coordonne transversale y). On va prsenter le cas dune plaque
mince dans une eau profondeur constante. On a donc rsoudre :
z
On en dduit que l =
x
h
xx =
1
pour
tt x a
gh
pour
t = h xx x a
x = a
x=a
Les conditions aux limites sont :
( x vitesse ) et ( t pression) sont continues en x = ma
On sintresse lefficacit de la planche en surface comme brise lame
des ondes venant du cot droit (x = + ). La solution gnrale de notre
quation diffrentielle a la forme: (x, t ) = F(x ct ) + G(x + ct ) o c = gh
Il est naturel de chercher des solutions harmoniques :
(x, t ) = (x )eit x a et (x, t ) = V (x )eit a x a
Nos quations deviennent alors
d 2 2
d2 i
0
x
a
et
+
+ V = 0 x a
dx 2 gh
dx 2 h
La premire quation a pour solution (x ) = Ae ikx + Beikx avec k =

gh
E.H.T.P.
33
ZORKANI Mohammed
Ainsi on a : (x, t ) = Ae i(kx t ) + Bei(kx + t ) cest la superposition de 2

ondes une progressive se propageant vers la droite Ae i(kx t ) et lautre
vers la gauche Bei(kx + t ) . On rappelle que dans ce cas londe incidente

provient de la droite, pour (x ) on peut crire (x ) = Be ikx + Reikx : x a
o B est lamplitude de londe incidente et R de celle rflchie : quil faut
dterminer. Sur la gauche on crit : (x ) = Teikx o T est lamplitude de
londe transmise ( dtermine galement). Si la planche est rdige et
fixe on a alors (x, t ) 0 il en rsulte alors que V (x ) = 0 donc sous la
planche on a : xx = 0 ainsi (x ) est une fonction linaire en x :
(x ) = x + . Puisque x (x, t ) est la vitesse horizontale de leau, il
rsulte alors que sous la plaque on a un courant donn par : eit
donc constant sous la plaque et sinusodal dans le temps.

Ecrivons maintenant les conditions de raccordement des solutions aux
Beika + Re ika = a +
Beika Re ika =
ik
discontinuits ; il en rsulte alors :
ika
= a +
Te
Te ika =
ik
On dispose donc de 4 quations pour 4 inconnues : R, T, et . Si on
dsire calculer la pression sous la plaque il suffit dappliquer Bernoulli
Lagrange : p(x, t ) = t = i (x ) eit . La solution de notre systme est
2a
2
o L =
par :
donne en fonction du paramtre adimentionnel =
L
k
ie 2i
R=B
i + 1
R
=
CR = cofficient de rflxion =
e 2i
B
T =B

1 + 2 2
i + 1
T
1

=
i e 2i C T = cofficient de transmissi on =
B
=B
1 + 2 2
a i + 1
= Be 2i
E.H.T.P.
34
ZORKANI Mohammed
Notons quon a : CR2 + C 2T = 1 qui exprime la conservation dnergie.

2a
= ka . Il est intressant dtudier comment varie la
Notons que =
L
pression sous la plaque : p(x, t ) = t = i (x ) eit = i (x + )eit
donc la partie relle, quil faut prendre, est donne par : si B est relle
p(x, t ) = p1(x ) cos t p 2 (x ) sin t avec
22 x
(
)
b
x
=
+1
1
p1 (x ) = B [b1 (x ) sin + b 2 (x ) cos ]
1 + 22 a
o
2 2
x
p 2 (x ) = B [b 2 (x ) sin b1(x ) cos ]
b (x ) =
2
1 + 22 a
8) Wave Run Up : Ascension des Lames

Le niveau auquel un ouvrage en mer (un mur, un revtement en pierre
comme un talus dune digue ) doit tre ras (nivel) est en fonction
principalement du runup lvation.
En laboratoire (1957) Saville propose la figure (8 1) pour dterminer le
runup Ru (hauteur dascension mesure verticalement des lames sur
une structure par rapport SWL) en fonction de la priode de londe, de
lamplitude de londe non rfracte et de la cotangente de la pente de
louvrage avec comme paramtre H0/ T 2 .
Ces courbes sont pour une paroi lisse et impermable avec une
profondeur deau de 1 3 fois H0/ . Ces courbes sont donnes par
U.S.Army Coastal Enginnering Research Center (1973) la figure (8 1)
rsume cela. Le tableau (8 2) propos par Battjes (1970) donne leffet
dune paroi non lisse et permable sur le runup . Le facteur r est le
rapport du runup donn par la figure (8 1) celui pour une paroi
permable et rugueuse.
Saville (1957) propose une procdure employer pour utiliser la figure
(8 1) aux parois composes de plusieurs pentes : une pente
hypothtique unique est construite partir du point de dferlement au
point dascension des lames sur louvrage compos de plusieurs
pentes
ainsi
on effectue le calcul comme si on a une paroi unique
run : up
= Ru
hb = db
SWL
pente hypothtique
E.H.T.P.
35
ZORKANI Mohammed
dont on connat la valeur avec laquelle on effectue les calculs comme

avant; ce calcul se fait par essais et ttonnements successifs garce la
figure (8 1) et on compare avec la valeur estime jusqu leur accord.
Si non on recommence jusqu laccord souhait.
Une formule empirique du run up dune onde dferlante sur un talus
lisse est propose par Hunt : [Proc. ASCE 85 WW3 Sept 1959 pp123-152]
Ru
H
==
tan g
H L0
pour : 0,1 2,3 &
hauteur Run down

hauteur Run up
Rd
Ru
= (1 0,4 )
On signale que si la hauteur de louvrage nest pas correctement

dimensionne il se produit un problme de franchissement de celui ci
qui engendre des inondations des quais et une agitation de lautre cot
labri ( dans ce cas un systme dvacuation est prvoir) : on traitera
ultrieurement le problme de franchissement : qui prsente une
grande importance pour le dimensionnement des ouvrages maritimes et
pour calibrer le rseau dassainissement portuaire (voir ch05).
N.B. : En eau profonde Mitchell (1893) propose par le critre
Hb
=
0
,
142
L
Hb
0,b max
2
cinmatique :
2 m / s = 0,267
T
L 0,b
gT 2
= 1,2 avec L 0 =
2
L0
On peut calculer le coefficient de rflexion R(module du rapport de

lamplitude rflchie par celle incidente) dun talus sur lequel londe
dferle par (selon Miche R=1 pour une onde nondferlante sur le talus)
tg
gT 2
= tg
o =
2H
H L0
R = 1 autrement
Rappel Mathmatique
Formules de transformation de Gauss
& Les formules intgrales (ou les identits) de Green
Si est un oprateur linaire on a lgalit :
r
r
d = (n) dS Thorme de Gauss
R 0,1 2 si R1
( )
N.B. : on pose
r r 3 Vi r
r
grad
;
V=
ei
r 3r
x
=
i
1
on associe
i
= ei

r
rr
r
3
i =1 x i
vr = rotvr ; V rr = V r
i
xi
i =1
( )
ainsi
E.H.T.P.
36
ZORKANI Mohammed
r
r r r r
r
r
r r
V = V = divV = (n) = n
r
r
r
r
r r
= = grad = (n) = n
r r r
rr
rr
r r
r
r
r = V r
= V (n) = Vn
r
r r r r
r
r
r r
V = V = rotV = () (n) = n ()
LLLLLLLLLLLLLLLLLLLLEtc.
On obtient alors les formules particulires du thorme de Gauss :
r r
r r
formule intgrale de la divergence
V
d
n
V
dS
(Gauss Ostrogradski)
D r
r
d = ndS formule intgrale du gradient
D
r r r
r r r
V r d = V n r dS
D
r
r r
rotV d = n V dS formule intgrale du rotationne l

D
r
r r
r
= r
r r
r r
r r
Si r = r r r
r
d = r n dS
D
(n) = r n
Applications :r Les identits
r
r de rGreen
r
Nous avons = + alors daprs le
thorme de la divergence on a :
r
r
r
r
r
dS
d = + d = ndS =
n
D
D
r
r
dS 1ire identit de Green
+ d =
n
D
r
r
Si on change on obtient + d =
dS et par
n
D
soustraction il en rsulte :
(
(
(
(
)
)
)
)
( )
( )
( )
( )
( )
( )
dS 2me identit de Green
( ) d =
n
n
D
Il sensuit que pour un coulement potentiel de vitesse cest dire

dS = 0 (1)

n
n
r 2
Si maintenant on prend alors

dS = + d ; si
n
D
en plus est le potentiel de vitesse dun champ dcoulement :
que = 0 lgalit :
( )
E.H.T.P.
37
ZORKANI Mohammed
( )
r 2
2 0 alors d =
dS
D
n
r r
Comme v : cette galit met en vidence lnergie cintique .
Les relations prcdentes mettent en jeu 2 fonctions et , do lide
dintroduire un potentiel lmentaire pour que nos relations ne sexprime
que par rapport une seule fonction ?
Par exemple on utilise les potentiels de vitesse singulier :
1
c
1
c
si D E 2 : =
= log
log rPQ
si D E 3 : =
=
2
rPQ
4rPQ
rPQ
i = 1,2,3
(3 dimensions ) : 3D
i = 1,2
(2 dimensions ) : 2D
d
te
22
r = (x z i ) remarquons que le flux = c car
=
pour rPQ = C
dr
n
i
o P et Q sont 2 points qui appartiennent au domaine fluide D ou sa
frontire D = D . Ainsi est une fonction dfinie pour P Q mais
singulire pour P = Q et elle satisfait lquation diffrentielle de Laplace
Dans ce cas la premire formule de Green donne :
r 1 r
1
Q (Q ) dQ (Q )
C (P ) = Q
r dS Q (2)
r
n
D
Q PQ
PQ
4 si P D
C = 2 si P
o
0 si P D
en
Si maintenant on change
alors la premire identit de Green
donne galement :
r
r
1
1
1 (Q )
(
)
(
)
Q
d
Q
d Q =
dS Q (3)
Q
Q
Q
r
r
r
n
PQ
D PQ
D
PQ
Q
Si on introduit la relation (1) dans (2) on obtient une reprsentation
intgrale du potentiel (tel que 0 ) nutilisant bien entendu que des
informations concernant le comportement de ces fonctions sur la surface
(enveloppe) du volume D liquide.
r
nQ
(D)
( )
E.H.T.P.
38
ZORKANI Mohammed
C (P ) =
1 (Q )
dS Q (Q )
rPQ nQ
nQ
r dS Q
PQ
b
b
simple couche
double couche
Cette formule trouve beaucoup dapplications dans la pratiqueSouvent
dite la 3me formule de Green.
Daprs la formule de Green (1), on pourrait introduire dans
reprsentation (3) un potentiel ayant dautres proprits :
1
(P, Q ) =
+ (Q ) tel que : = 0 dans Dint D ext
rPQ
D int
D ext
On obtient ainsi une nouvelle formule de reprsentation du potentiel :

(Q )
C (P ) = (P, Q )
dS Q (Q )
(P, Q ) dS Q
nQ
nQ
Les formules(1,2 et 3) contiennent 2 potentiels connus, on les utilisent

souvent en hydrodynamique navale :
Le potentiel de simple couche : engendr par une distribution
1
de sources 1(Q ) : = 1(Q )
dS Q
rPQ
Le potentiel de double couche : engendr par une distribution
1
de diples 2 (Q ) : = 2 (Q )
dS Q
nQ rPQ
La houle de Stokes irrotationnelle du 2eordre :

Pour les houles de grandes amplitude mais finie Stokes en 1880 pour
effectuer ses calculs dveloppe le potentiel de vitesse sous la forme :
= H(1) + H2 (2 ) + H3 (3 ) + H4 (4 ) + L
Le premier terme H(1) correspond la linarisation prsente

prcdemment (Thorie dAiry).
Houle du 2e ordre : qui est valable pour h 0,01 gT 2
Le potentiel de vitesse est donn par
3 H2 ch 2k (z + h)
HL ch k (z + h)
=
sin(kx t ) +
sin 2(kx t )
16 T
2T sh kh
sh 4kh
Le profil de la surface libre est donn par
E.H.T.P.
39
ZORKANI Mohammed
H2
H
3
1 +
coth kh cos 2(kx t )
= cos(kx t ) +
2
4L 2 sh2kh
Le paramtre du dveloppement est la cambrure H L1 .
On constate que la deuxime harmonique; excit par la non linarit
est : de priode 0,5T 2 et de longueur donde 0,5L 2k
Les profils sont symtriques par rapport des plans passants par les
crtes et les creux . Les crtes sont plus hautes et plus courbes alors
que les creux plus plats : ce rsultat nest pas mis en vidence par la
thorie linaire (Airy).
La hauteur des crtes est donne par :

H H2
3
1 +
coth kh
max = +
2
4L 2 sh2kh
H H2
3
1 +
coth kh
Les creux sont au niveau : min = +
2
4L 2 sh2kh
Ce qui montre que le terme introduit na dimportance que pour une
houle damplitude crte creux H importante. La vitesse orbitale sont :
3 H H ch 2k (z + h)
H ch k (z + h)
cos(kx t ) +
cos 2(kx t )
u = x = T
sh kh
4 T L
sh 4kh
w = = H sh k (z + h) sin(kx t ) + 3 H H sh 2k (z + h) sin 2(kx t )
4 T L
T
sh kh
z
sh 4kh
Les trajectoires des orbites sont obtenues par les intgrales :
= u dt = w dt
Il faut effectuer les calculs numriquement de proche en proche pour
D
D t = u + u x + w z
obtenir les trajectoires en utilisant les relations :

D = w + u + w
D t
z
x
e
Pour la houle de Stokes de 2 ordre, posant = (k x t ) on obtient :
E.H.T.P.
40
ZORKANI Mohammed
H ch k (z + h)
H2 3 ch 2k (z + h)
sin +
=
1 2 sh2kh sin 2
2
2
sh
kh
8
Lsh
kh
H2 ch 2k (z + h) 2t
+
4L
T
sh2kh
2
= H sh k (z + h) cos + 3H sh 2k (z + h) cos 2
2 sh kh
16L
sh 4kh
H2 sh2k (z + h)
+
8L
sh2kh
Les orbites (trajectoires) des particules fluides ne sont plus fermes;

ainsi le mouvement de londe progressive saccompagne dun
dplacement de matire ( entranement ) : cest un transport de masse
( un courant ) :
Dans le cas de la houle dAiry, cestdire linaire monochromatique,

on na pas de transport de masse. La valeur de la vitesse de ce transport
de matire dpend de la cambrure de la vague est donne par
2H2
ch 2 k (z + h) sh 2 k h
lexpression : Um (z ) =
2
2kh
2LTsh kh
Dont la valeur en eau peu profonde et au voisinage de la surface libre
2
est donne par :
h
h 1 L
H L
Um (0 ) = coth 2 coth 4
L
L 4 h
L T
2
H
Dans le cas dune profondeur infinie, ce courant vaut : Um (0 ) = c ,
L
o c est la clrit de londe. Du fait de ce courant dentranement les
trajectoires de type elliptique ne sont plus fermes car elles se dplacent
la vitesse Um (z ) . A chaque priode les orbites avancent de Um (z )T :
2
()
()
Um 0 T Um 0 T
Au fond le courant vaut :
Um ( h) =
()
Um 0 T
H
8hT
2
h
h
4
L
L
h
sh2 2
L
sh4
E.H.T.P.
41
ZORKANI Mohammed
Les trajectoires se rduisent un simple mouvement de va et vient,

dont la rsultante est en sens inverse de la propagation de la houle :
z
x
h
sens de propagation
Daprs cette thorie de Stokes au 2e ordre les vitesses maximales sous

la crte et le creux de londe sont respectivement donnes par :
H
3 k H2
Max
+
u ,crte = u,crte =
(
)
2
sh
kh
16 sh 4 (kh )
H
3 k H2
u
Max
,creux = u,creux = 2 sh (kh ) 16 sh 4 (kh )
Une comparaison des mesures des vitesses maximales montre un bon

accord avec la valeur sous la crte, cependant sous le creux les
mesures montrent que la vitesse mesure est plus petite que celle
thorique dans une eau de profondeur infrieur 3m. En se basant sur
ces mesures en eau peu profonde h 0,01 gT 2 Van Rijn propose :
~
u ,crte = u
H
o = 1 + 0,3 et u donne par la thorie linaire
~
h
u ,creux = (2 ) u
Au 3me ordre et plus la clrit de londe dpend de la cambrure :

c
2h H 14 + 4ch2 4h L
= th
1 +
c0
L L 16sh 4 2h L
on constate quon a une dispersion en frquence et en amplitude.
E.H.T.P.
42
ZORKANI Mohammed
Couche limite de la houle

La couche limite donde est une couche mince formant la transition entre
le fond la couche suprieur o lcoulement fluide est irrotationnel &
oscillatoire :
coulement
turbulent
vitesse maximale de
l' coulement irrotationne
coulement
laminaire
u
u
couche limite d' onde

Pour coulement laminaire on a :
Jonsson 1980 propose : = wave =
Manohar 1955 propose : =
4,6
2
= longueur de Stokes = (m)
T
= vis cos it cinmatiqu e m 2 / s
T = priode d' oscillation (s )
Dans le cas dun coulement turbulent Jonsson et Carlson (1976)

proposent :
H
30
30
H
= 1,2
log
pour 10 500 (a)

2Ks
2Ks
Ks
Ks
Cette quation peut galement tre reprsente par :
1
4
0,5 H
= 0,072
0,5 H
K
s
Lquation (a) est base sur des rsultats exprimentaux.

Les rsultats thoriques de Freds (1984) avec une erreur de 20%
sont approximativement donns par :
E.H.T.P.
43
ZORKANI Mohammed
1
4
0,5 H
= 0,15
0,5 H
K
s
La transition un coulement oscillatoire pleinement turbulent sur un

fond plat peut tre estimer par la formule :
(u,critique )
o d50
1
4
H
= 5770
2 d50
est le diamtre moyen des particules solides du fond.
Transport de masse par

des ondes non dferlantes
Ecoulement oscillatoire dun fluide parfait :
Stokes (en 1847) est le premier mettre en vidence que les particules
ne dcrivent pas exactement des orbites fermes dans une onde de
faible amplitude se propageant dans un fluide parfait (irrotationnel) en
coulement oscillatoire. Les particules fluides possdent une vitesse
Lagrangienne moyenne au seconde ordre (nomme : Stokes drift)
dans la direction de propagation de londe.
La vitesse orbitale horizontale augmente avec z au dessus du fond.
Par consquence, une particule au sommet dune orbite sous une crte
va plus vite dans la direction de propagation que si elle est sous un
creux donde. Par dfinition the Lagrangien Stokes drift ne peut pas
tre dtecter par des mesures en un point fixe.
La valeur instantane du <Lagrangien Stokes drift> horizontal Us dune
particule deau qui a une position moyenne (x1, z1 ) est Us (x1 + , z1 + )
o (, ) sont les coordonnes de la particule sur sa trajectoire. Une
approximation de Us est donne par :
U
U
Us (x1 + , z1 + ) = Us (x1, z1 ) +
+
z
x
En utilisant la thorie linaire et en prenant en suite la moyenne sur une
priode , on obtient ainsi la vitesse en moyenne temporelle (note ici par
1
ch 2 k (z h)
(qu. a2)
une barre au dessus) : Us (z ) = k H2
2
8
sh (kh )
o on a pris lorigine des z la surface libre du haut vers le bas :
z=0
z
z=h
E.H.T.P.
44
ZORKANI Mohammed
Au fond (z = h)
k H2
: Us =
A la surface (z = 0 ) : Us =
8 sh2 (kh )
k H2ch 2 kh
8 sh2 (kh )
Pour des ondes se propageant sur un fond horizontal dans un domaine
non limit le dbit volumique m 2 / s sur la profondeur deau h est :
Ms = h0Us (z )dz =
H2 sh(2kh )
g H2
H2
th(kh ) =
8
8c
16sh2 (kh )
gT
g
th (kh ) = th (kh )
o c = clrit de l' onde =

2
H2
8
Pour des ondes se propageant sur un fond horizontal dans un domaine
limit Il est logique dimpos un dbit volumique nul en chaque position
x, ce qui conduit :
k H2
sh (2k h)
Us (z ) =
ch 2 k (z h) 2 k h (qu. a1)
2
8 sh (kh )
Cette quation montre que le courant rsultant est la somme dun

<Lagrangien Stokes drift> dans le sens de propagation de londe et dun
courant de retour uniforme dans le sens oppos. Ainsi on a un dbit vers
dans le sens de propagation de londe proche de la surface (vers la cote)
et un dbit ngative proche du fond ( vers le large : dans le sens oppos
la propagation de londe) : ce mcanisme de transport de masse
ncessite la prsence dun gradient de pression horizontal (cisaillement
est absent car le fluide est parfait par hypothse) qui ne peut tre caus
que par une lvation de la surface libre vers la cote (wave set up).
Le dbit volumique m 2 / s en une position fixe (x) dans un fluide illimit
peut tre galement dterminer par une approche Eulrienne par :
1 T ( t )
Me 0 U(z, t ) dz dt
T h
o U est la vitesse horizontale instantane au niveau z, et est le
dplacement de la surface libre par rapport au niveau moyen deau MWL
g H2
La thorie linaire donne : Me =
. La mthode Lagrangienne et
8c
Eulrienne conduisent au mme rsultat de dbit volumique. Mais la
distribution verticale de la vitesse de transport massique moyenne est
diffrente pour les 2 approches.
Cette quation en eau profonde (kh ) 1 se rduit : Ms =
E.H.T.P.
45
ZORKANI Mohammed
Effet de la viscosit dans un coulement oscillatoire turbulent :

Longuet Higgens (1953) a mont quil existe, pour un fluide rel
visqueux , un transfert en moyenne temporelle de la quantit de
mouvement dans la direction de propagation dans la couche limite par la
U
induisant un courant Eulerien moyen Ue en
diffusion visqueuse
plus <Stokes drift> Us .
La vitesse de transport totale moyenne Um est dfinie par :

U
U
Um = Ue + Us = Us +
U dt +
V dt
z
x
Pour coulement dans une couche limite laminaire Longuet Higgens a
obtenue :
z
z
2
z
5 8 cos e + 3 e
Um =
2
16 sh (kh)
k H2
2
.
Lquation de Um possde une valeur maximale donne par :
o = paisseur de la couche lim ite la min aire =
Um = 1,376
k H2
4 sh2 (kh)
u2
= 1,376
o
c
u = valeur maximale de la vitesse hors C. L.
c = clrit de l' onde = k

En admettant un dbit nul sur toute la profondeur deau Longuet
Higgens a obtenu :
Um (z ) = Us (z ) + Ue (z ) =
k H2
z
F
(qu. a3)) o
2
8 sh (kh) h
z2
3 kh
z
z
F = ch2k (z h) + +
sh2k h 3 2 4 + 1
2 2
h
h
h
3 sh 2k h 3 z 2
+
+ 2 1
2 2 kh
2 h
Pour un coulement oscillatoire dans un fluide illimit la vitesse de

transport Eulrienne en moyenne temporelle (induite par les effets de
viscosit) peut tre dcrite par :
E.H.T.P.
46
ZORKANI Mohammed
3 k H2 1
+ k 2H2 (h z ) coth(kh) (qu. a4))
2
16 sh (kh) 2
Le dbit moyen sur la profondeur deau est :
0
g H 2 3 k h H2 1
M = Ue + Ue dz =
+
+ k 2h 2H2 (h z ) coth(k h)
2
8c
16 sh (k h) 4
h
Ue (z ) =
z
h
0,4
0,8
4Um
0,8
0,8
1,6
2,4
k H
< Stokes drift > dans fluide limit (qu (a1) )
< Stokes drift > dans fluide illimit (qu (a2) )
transport de masse Longuet - Higgens fluide il lim it (qu (a3) )

transport de masse de Craik fluide il lim it (qu (a4) )
Transport de masse
par
des ondes dferlantes
Quand une onde dferle elle gnre un courant parallle la ligne de
cote (longshore current) et un autre vers le large (undertow). On
prsentera par ailleurs le modle mathmatique de Longuet Higgens
pour dterminer la circulation marine lchelle de la houle (par
lintroduction du tenseur de radiation).
Au dessus du niveau du creux donde dferlante existe un dbit vers la
cote. En premire approximation, ce dbit peut tre estimer par :
g H2
M=
8c
En utilisant c = g h en eau peu profonde, il en rsulte que :
E.H.T.P.
47
ZORKANI Mohammed
1 g 2
H
8 h
En admettant quil y a pas un transport total net deau sur la profondeur
deau, la valeur moyenne du courant de retour sous le creux, est donn
par :
1 g 1 2
Um,off =
ht H
8 h
M=
En prenant h t = 0,8 h il rsulte donc :

0,75Hb
Um,on
0,25Hb
ht
Um,off = 0,15 g h
3
2
H2
SWL
Um,off
Complment
C1- Variation thorique de lamplitude et de la cambrure des vagues
par fond dcroissant en absence de la rfraction :
On a vu que lnergie transmise (flux dnergie = puissance) par unit de
temps et de longueur travers un plan vertical fixe par mtre linaire de
crte est :
gH2L
2 kh E 1
2kh n E gH2
1 +
= 1 +
=
=
P=
cg
16T sh 2kh T 2 sh 2kh
T
8
2k h
c
1 +
= nc
2 sh 2k h
Au large, cest dire en eau profonde kh 1 , on a :
o c g est la vitesse de groupe locale : c g =
n = 0,5 en eau profonde

gH02L 0 gH02
c g 0 car
=
P0 =
n = 1 en eau peu profonde
16T
8
Si le fond a une pente faible et quune nergie apprciable nest rflchie
non plus dissipe alors on aura conservation de lnergie transmise.
Si les lignes de niveau sont parallles aux lignes de crte (pas de
rfraction) une nergie qui franchie un mtre de crte au large est la
mme que celle qui franchie un mtre de crte prs de la cte :
Calcul de la variation de lamplitude : P0 = P
E.H.T.P.
48
ZORKANI Mohammed
c g0 c 0
H
1
=
=
=
H
c
cg
2
k
h
g
0
th kh 1 +
sh 2kh
Posant x = k h alors
2
H
1
1
=
=
dont la drive est :
x
H
2
x
0
thx + 2
th x 1 +
ch x
sh
2
x
1 x th 2x
1
2
sannule
pour
:
th
x
=
x = 1,198 alors
2
x
x
ch2 x thx + 2
ch x
H
h
h
0,191
=
= 0,158
= 0,9129 ; = 0,191 et
L
L
1
,
198
H
0
0 min imum
H
x
= 1 quand thx + 2 = 1 2x = 1 + e 2x
Cependant
ch x
H0
h
0,191
h
h
x = 0,639 = k h = 2 = 0,1016 et
=
= 0,057
L
L 0 1,198
L
cest le point isomtrique .
C2- Variation de la cambrure :
Les cambrures au large et en situation quelconque sont respectivement
H
0 = 0
L0
et
H
H
H
=
=
alors
L
L
H0
L
2 L2 2
= 2 2 = 2 th 2k h
0 L0 0
H0

1
H
1
=
=
2kh H0 thkh
0
th3 kh 1 +
sh 2kh
H
H
La drive nannule au environ de kh 2,17 soit 0,34 et
0,33
L0
L
2
0,97
do un minimum a peine marqu :
0,985 . On
0
0
remarque que ce sont les houles les moins cambres au large qui
sapprochent de la terre et subissent le plus forte augmentation de la
cambrure.
E.H.T.P.
49
ZORKANI Mohammed
C3 - Augmentation de la vitesse horizontale : en thorie linaire
H g
cos(k x t )
u(x, z, t ) =
2 h
On a tablit en eau peu profonde :

v (x, z, t ) = H 1 + z sin(k x t )
T h
Ainsi lamplitude de la vitesse horizontale est u ~
H g gH
=
et celle
2 h 2c
H
z gH
k (z + h) car c = g h .
1 + =
T h 2c
La composante horizontale maximale de la vitesse en eau profonde et
g H0
.
en surface est donne par : u0 =
2 c0
verticale est v ~
te
u
gH 2c 0
H c 0 H L0
car T = C .
=
=
=
=
u0 2c gH0 H0 c
H0 L
0
Le rapport des vitesses maximales est gale au rapport des cambrures.
Au dferlement, o atteint la valeur limite, on obtient alors la limite du
u
.
rapport des vitesses
u0
Ainsi le rapport est :
On peut aussi exprimer u en fonction de H0 , L 0 et h . Pour kh assez on a

H
H02
H02
g H
1
= g2
et u2 =
TH
2
16
2 kh
4c
u ~
1
1
3
2
avec T 2 = 2
L0
do
g
g 2
4 h 4
H
L
0 0
2
2(2 )4
Cette vitesse tend vers linfini si la profondeur deau tend vers zro
h 0 . Mais en ralit sa valeur est limit par le dferlement.
H
Au large (eau profonde) on a : u0 = 0
2
u
=
u0
1
g 2
1
1
(2)4
1

1
2 (
2g)2
2g
= 2g H0L 02 donc :
L0
L
0
3
4
0,18
L
0
3
4
E.H.T.P.
50
ZORKANI Mohammed
Superposition des mouvements Harmoniques

La superposition dondes de mme priodes est :
2t
2t
2t
= a cos
= cos
a cos + sin
a sin
T
T
T
r = ( a cos )2 + ( a sin )2
2t
= r cos
o
T
tg = ( a cos ) ( a sin )
Par exemple prenons 2 composantes :
2t
2t
2t
= a cos
+ a / cos
/ = r cos

T
r = a 2 + a / 2 + 2aa / cos /
o
a sin + a / sin /
tg =
a cos + a / cos /
Ainsi si les 2 composantes sont en phase = / alors :

2t
= a + a / cos
proche des ventres
T
Mais si la diffrence de phase diffre dune demi priode alors :
2t
= a a / cos
proche des nuds
T
si en plus a = a / = 0 oscillation nulle aux nuds (onde harmonique

stationnaire pur).
Si maintenant on superpose 2 ondes harmoniques de priodes
diffrentes mais voisines (battement) :
2t
2t
= a cos
+ a / cos / /
T
la rsultante ne peut pas tre reprsente par une onde harmonique.

Analytiquement si :
= a cos(2m t ) + a / cos 2n t / avec (m n) est petit alors
r 2 = a 2 + a / 2 + 2aa / cos 2(m n)t + /
= r cos(2m t ) o
a sin + a / sin 2(m n)t + /
tg =
a cos + a / cos 2(m n)t + /
On peut donc considrer la superposition de ces 2 ondes comme une

harmonique dont les lments sont r & , qui ne sont pas constants mais
variables lentement dans le temps ayant la frquence (m n) :
{
{
}
}
E.H.T.P.
51
ZORKANI Mohammed
{
}
} = 1 dont les valeurs
Lamplitude est maximale quand : cos 2(m n)t + / = +1 et
minimale quand cos 2(m n)t + /
correspondantes sont respectivement a + a / et a a / .
Edge waves on a sloping beach

by F. URSELL 1952 Proc. Roy. Soc. A214 pp79 97
On considr des ondes de gravit dans un canal houle dont :
(ii) la longueur est finie et de profondeur constante
(ii) la longueur est infinie et de profondeur constante
Dsignons par Ox laxe le long du canal et par Oy laxe vertical vers le
haut et Oz laxe transversal. On supposera avec Ursell que le canal a
une profondeur infinie (Cette hypothse ne changent pas le caractre du
spectre).
(ii)
Dans un canal de profondeur infinie limit verticalement par les

plans ( x = 0 , x = a ; z = 0 , z = b ) les potentiels de vitesse des
modes normaux (naturels) sont donns par :
m 2 n 2 imn t
mx
nz
imn t
(
)
mn x, y, z e
= Cmn cos
cos
exp y 2 + 2 e
a
b
b
a
o n et m sont 2 entiers et Cmn est une constante complexe. On a :
m2 n2
2
mn
= g 2 + 2 (1)
b
a
Le spectre de frquences donnes par (1) est discret et infini. On peut
2
2
facilement vrifier que : 0 m
+1,n m,n g a (2).
On prendra la partie relle de la solution. La solution du mouvement libre
est donne par : (x, y, z, t ) = mn (x, y, z ) eimn t (3)
m n
(ii) Supposons que la longueur du canal tend vers linfini a lquation

(2) nous suggre que le spectre devient continu. Les modes normaux
sont :
E.H.T.P.
52
ZORKANI Mohammed
n (x, y, z; k ) e
()
in k t
nz
= Cn cos(kx ) cos
exp y
b
2 n 2 2 in (k ) t
k +
e
2
n2 2
n2 (k ) = g k 2 + 2
b
et k est un nombre positif. Pour le mode bidimensionnel (n = 0 ) toutes les

valeurs relles de sont valeurs propres, pour les modes 3Dimensions
(n0) tout les rels sont gn b . Quand n est donn il existe une
limite infrieure (frquence de coupure) au dessous de laquelle
nexiste pas de modes normaux.
Le mouvement libre est maintenant de la forme :
2 n2 2 in t
nz
(x, y, z, t ) = cos
0 Cn (k ) cos(kx ) exp y k + 2 e dk (4)
b
n
b
Ondes de coin non visqueuses : Inviscid edge waves

On va maintenant tudier le cas dun spectre mixte. On considre le cas
de modes normaux excits par une plage de pente lextrmit dun
canal houle semi infini. Le potentiel de vitesse est dfini dans la
rgion 0 y x tg & 0 z b o est satisfaite :
o
+ g y = 0 en y = 0

tg en y = x tg
=
2 = 0 avec les conditions aux limites
y
x

= 0 en z = 0 et z = b
La dernire condition limite montre que le potentiel est de la forme :

mz
it
= e cos
fm (x, y )
b
m
1
Si est antisymtrique par rapport z = b ce que nous admettrons
2
z
it
alors la srie se rduit : = e cos(2r 1) f2r 1(x, y )
b
r
o f2r 1 vrifie lquation :
2
2
2
2
(2r 1) 2 f2r 1(x, y ) = 0
+
2
2
x
y
b
et les 2 premires conditions aux limites.

Il est remarquer que le cas r = 1 est particulier :
E.H.T.P.
53
ZORKANI Mohammed
f2r 1(x, y ) = f1[(2r 1)x, (2r 1)y ]

auquel il faut porter une intention particulire.
Maintenant considrons une solution de la forme :
b
Une solution est : = cos kz exp[ k (x cos + y sin )] ei t (5) et
= cos kz F(x, y ) ei t
k=
daprs la premire condition limite on a : 2 = gk sin

Le mode (5) est d Stokes (1846) il sera dsign par Stokes edge
wave.
gk sin
g sin
Puisque grad dxdydz est fini, la frquence

=
=
2
2
4b
4
est une frquence discrte du spectre.
Pour une pente du fond faible the Stokes adge wave nest pas le seul
mode discret mais le premier dune suite. En considrant le potentiel :
k (x cos + y sin ) n
k (x cos(2m 1) y sin(2m 1) )
+ A mn e
= e
m =1
k (x cos(2m + 1) + y sin(2m + 1) ) cos kz e i t

m =1
qui vrifie la 2me et la 3me conditions. La premire condition est
satisfaite si :
m m tg(n r + 1)
& 2 = gk sin(2n + 1)
A mn = ( 1)
r =1 tg(n + r )
La vitesse est finie dans le secteur quand x si (2n + 1) 2 ; le
mode de Stockes correspond n = 0 . Ceci a t exprimentalement
vrifi par F. Ursell au laboratoire.
Voir galement J. Fluid Mech (1995) vol301 par P. Blondeaux &
G. Vittori :
The nonlinear excitation of synchronous edge waves by
monochromatic wave normally approaching a plane beach.
n
+ A mn e
Thorie donde cnodale
Quand la profondeur deau relative dcrot vers peu prs h L = 0,1, la

thorie de Stokes cesse dtre valable. Une autre approche thorique est
necessaire : cest la thorie cnodale, qui tait initialement dveloppe
par Korteweg et De Vries en 1895 (Phil. Mag. 5 Ser. 39).
La clrit de londe cnodale est donne par :
H 1 1 ( )
H L2
2
o
c = gh 1 +
2
h3
h 2 ( )
E.H.T.P.
54
ZORKANI Mohammed
o est un paramtre elliptique ou module. ( ) est lintgrale elliptique

complte de de premire espce, ( ) est lintgrale elliptique
complte de de seconde espce. Dont on rappelle les dfinitions :
On continuera prsenter les rsultats de cette thorie aprs le rappel
Rappel mathmatique : Les intgrales & fonction elliptiques
En physique mathmatique applique apparaissent en non linaires
des solutions qui font appel aux fonctions et intgrales elliptiques qui
sont la gnralisation des fonctions sinusodales prisent sur le cercle
trigonomtrique alors que les premiers sur une ellipse. On va les
rsumer dans ce qui suit :
Les formes de Legendre :
La reprsentation de Legendre de ces intgrales elliptiques de premire
et deuxime espce est :
d
(, ) =
0 1
0 1 2 sin 2
ou
= sin , 0 ( 2)
(, ) = 1 2 sin2 d
o et sont dsigns respectivement par le module et lamplitude de

lintgrale elliptique en question. La quantit / = 1 2 est dsigne
par le module complmentaire.
Ces intgrales sont tabules pour les valeurs de = arcsin et entre
0 et 2 ; en se donnant le valeur de 2 dans une intgrale quon veut
valuer, on doit dabord prendre la racine carre pour avoir , puis on
cherche = arcsin par la fonctions trigonomtrique naturelle sinus ,en
suite trouvait (, ) ou (, ) par les tables des intgrales elliptiques.
Les intgrales elliptiques compltes :
Les intgrales elliptiques compltes de premire et deuxime espce
sont les valeurs de (, ) ou (, ) pour = 2 , on trouve :
d
2
ou ( ) = , =
2 0 1 2 sin2
2
ou ( ) = , = 1 2 sin2 d
2 0
Il existe des tables part pour ces intgrales elliptiques compltes.
N.B. : on peut valuer par une mthode numrique ces intgrales.
On a par dfinition des fonctions intgrales elliptiques compltes :
(, n ) = 2n (, )
(, n ) = 2n (, )
E.H.T.P.
55
ZORKANI Mohammed
ceci rsulte des proprits de la fonction priodique sinus.

Si la borne infrieure dintgration nest pas zro, on peut crire :
2
2
1
d
d
d
= (, 2 ) (, 1 )

=
2
2
2
2
2
2
1 1 sin
0
1 sin 0 1 sin
on a une formule similaire pour (, ) . Si lune des bornes dintgration
est ngative on peut utiliser le faite que (, ) est impaire en :
d
d
(,) =
=
= (, )
0 1 2 sin 2
0 1 2 sin 2
et (,) = (, ) .
N. B. : pour des valeurs petites de on valuer avec une approximation
bonne valuer les intgrales par un dveloppement en srie.
Les formes de Jacobi :
Si nous posons sin = x dans les formes de Legendre on obtient alors
les formes intgrale de Jacobi :
dx
dx
=
x = sin dx = cos d d =
cos
1 x2
correspond
= 2
x = 1
alors
(, ) =
1 2 sin2
dx
(1 x )(1 x )
2
(, ) = 1 sin d =
2
2 2
1 2x 2
1 x2
1
d
2
= ( ) = , =
=
2 0 1 2 sin2 0
dx
dx
(1 x )(1 x )
2
2 2
1 1 2 x2
2
2
2
dx
= ( ) = , = 1 sin d =
1 x2
2 0
0
Pourquoi dsigne t on ces intgrales par elliptique ? tout simplement
elles sont lies au calcul de la longueur dun arc le long dune ellipse,
analogue au calcul trigonomtrique sur un cercle, en effet :
x = a sin
Lquation dune ellipse sous forme paramtrique est :
on
y = b cos
prendra a b (dans le cas o a b utilise: x = a cos , y = b cos ), on a :
ds 2 = dx 2 + dy 2 = a 2 cos 2 + b 2 sin2 d 2
puisque a 2 b 2 0 on peut crire :
E.H.T.P.
56
ZORKANI Mohammed
ds = a a b sin d = a 1
2
a2 b2
2
sin2 d
a
cest une intgrale elliptique de second espce o 2 = a 2 b 2 a 2 = e 2
(e est lexcentricit de lellipse). Si on veut le primtre de lellipse doit
varier de 0 2 , le rsultat est : 4a(, 2) . Pour un petit arc on peut
introduire ses bornes dans lintgrale est obtenir ainsi (, 2 ) (, 1 ) .
Rappelons nous que :
x
dx
u=
dx = sin 1 x
0 1 x2
qui dfinit u en fonction de x , ou linverse ; ainsi x = sin u . Dune
x
dx
dfinit u en fonction de
manire similaire : (, ) =
2
2 2
0
1 x 1 x
)(
(ou bien fonction de x = sin ) [on admet que est constante]. On crit :
x
dx
u=
= sn 1x x = sn u (lire ess en de u)
0
1 x 2 1 2 x 2
)(
Comme = amp u est lamplitude de lintgrale elliptique u = (, ) et

x = sin , on a :
sn u = sin = sin (amp u)
sn u est une fonction elliptique. Il existe dautres fonctions elliptiques qui
ont une ressemblance avec les fonctions trigonomtriques. On dfinit :
cn u = cos = cos (amp u) = 1 sin2 (amp u)
= 1 sn2u = 1 x 2
d
1
dn u
=
= 1 2 sin2 = 1 2 sn2u = 1 2 x 2
du du
d
la valeur de du d se calcule partir des expressions donnant u = (, ) .
Il existe des formules, comme en trigonomtrie, qui relient ces fonctions :
comme les formules daddition, dintgrales, de drives, Etc.
Exemple :
d
(sn u) = d (sin ) = cos d = cn u dn u
du
du
du
Pour plus de dtail consulter louvrage :
DAbramowitz, Milton, and Irene A. Stegun :
[ Handbook of Mathematical Functions with Formulas,
Graphs, and Mathematical tables ]
E.H.T.P.
57
ZORKANI Mohammed
National Bureau of Standards, Applied Mathematics Series 55

U. S. Government Printing Office, Washington, D. C. 1964
H 1 1 ( )
La clrit c = gh 1 +

se situe entre 2 limites :
2 2
(
)
h
Celle des ondes sinusodales de cambrure H h faible et celle des ondes

solitaires pour une grande cambrure H h cest dire :
Pour les premiers (ondes sinusodales) :
2h2
H L2 tend vers
H
2 = 3
0 &
1 soit : c = gh 1
2
h
3
L
h
qui sapproche de la clrit donne en thorie linaire :

1
2h 2
1
2
2h
gL
c=
th
= gh1
2
2
L
3L

Pour les deuximes quand londe tend vers londe solitaire :
H L2 tend vers
H
H
2
= 3
1 &
0,73 soit : c = gh 1 +
2
h
h
h
qui sapproche 2% de la clrit de londe solitaire

H
c = g (h + H) = gh 1 +
h
on dmontre que dans ce cas T = .

H L2
Le paramtre 2 = 3 (dit paramtre dUrsell dont on a parl en page
h
52 : propagation ondes sur une plage) caractrise le passage des houles
de type Stokes celles de type cnodale.
La longueur donde L est donne par :
2
1
2
16 h
( )
L =
3
H
Llvation de la crte au dessus du SWL est donne par

3
a c 16 h L
=
( ) [ ( ) ( )]

H
3 H H
3
H
= 0,73
h
Ces rsultats sont au second ordre dapproximation.
Le dferlement selon cette thorie se produit pour :
E.H.T.P.
58

Théorie de La Houle

Transféré par

Droits d'auteur :

Formats disponibles

Théorie de La Houle

Transféré par

Informations du document

Description originale:

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Théorie de La Houle

Transféré par

Droits d'auteur :

Formats disponibles

ZORKANI Mohammed

Spctre nertique des ondes selon Munk

Intervalle de priode (sec)

Le spectre des ondes dues au vent stendent jusqu 30s. Au del,

thories, la plus simple et la plus utilise est celle faible amplitude

Ecoulement irrotationnel : pas de contraintes qui sexercent sur

Le fond est horizontal fixe(pas de transport de sdiment), impermable.

La position dune particule fluide dans le plan du mouvement est repre

Par dfinition la clrit est de londe donne par :

p = p a = C et en ngligeant les termes en vitesse car ils sont du second

Voyons voir ce qui se passe si la pression atmosphrique varie

Cherchons des solutions de la forme: = (x, z ) cos t + (x, z ) sin t

La solution gnrale est donc donne par :

On constate quil peut y avoir rsonance (amplification).

De mme pour les drives en y, z, et t :

g(1) (x, t ) + (t1) (x,0, t ) = 0 premier ordre ( Condition de Poisson )

en la reportant dans (1 1) on obtient

vrifie lquation de Laplace. Pour dterminer les 2 constantes

(x, z, t ) = 2Bekhch k (z + h) sin(k x t ) (1 8)

La condition (1 6/) implique que

est donne par v =

z = 0 . En y reportant notre solution (1 9 )

Si la profondeur relative est plus grande que 0,5 alors th

, Les orbites des particules sont

circulaires dont le diamtre exponentiellement dcroissant vers le

fond et sont proches de 0 pour 0,5

[si par exemple T = 10s et une amplitude de 2m alors c=15,6 m/s et

14% donc pas de dferlement ].

Si la profondeur relative est infrieure 0,05 alors th

Pour 0,5 0,05 on a des conditions de transition : c et L

dcroissent quand h diminue cest dire vers le rivage car :

Intervalle de priode (sec)

Rappel : Dveloppement limit de Taylor Young

est positif cest dire la surface libre

signe moins signifie que si

est donne par v =

sin(k x t ) pour cos (k x t )

Cette relation de dispersion contient 2 termes, reprsentant les 2 sortes

Cette valeur critique de la longueur donde est de 17mm pour leau

On peut maintenant classer les ondes :

On constate que c g c les rides (les trs courtes longueurs donde)

apparaissent se propager individuellement vers larrire relativement

Les ondes de gravit en eau peu profonde :

N.B. : les calculs prcdents rsultent du D. L. de la tangente

Cinmatique de londe et les pressions : sans tension superficielle

progressive dont le profil de la surface libre est

sin(k x t ) on dduit alors le

champ de vitesse en thorie linaire :

On observe que la vitesse rsulte de 3 termes :

N. B. : On observe quentre lacclration et la vitesse existe un

autour de leur position au repos, cest position moyenne, se calcule par :

Lorbite des particules fluides est donne par :

eau faible profondeur

Reportons notre solution dans le thorme de Bernoulli Lagrange

pression dynamique due lacclration verticale

Surpressions et sous pressions sont donc

4 ) Energie de londe et Puissance :

On constate quipartition dnergie : Ep =