Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

MFEhoule

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 38

Houle et Vagues.

Ecoulements en Milieux Naturels Cours MSF12, M1 UPMC


P.-Y. Lagree
CNRS & UPMC Univ Paris 06, UMR 7190,
Institut Jean Le Rond Alembert, Bote 162, F-75005 Paris, France
pierre-yves.lagree@upmc.fr ; www.lmm.jussieu.fr/lagree
9 fevrier 2016

Resume
Dans ce ce chapitre nous regardons la propagation de la Houle lineaire dite Houle de Airy : cest une perturbation non visqueuse infinitesimale de la surface de
leau, a profondeur quelconque. Nous etablirons la relation de dispersion (k) pour des ondes sinusodale de faible hauteur du type cos(t kx), nous parlerons de la
vitesse de phase et de la vitesse de groupe (vitesse de propagation de lenergie). En faible profondeur, nous etablirons que la vitesse est gh, on commence dailleurs
par un cas ou la hauteur deau est faible pour introduire quelques notions de base sur lequation donde.
Cours du module 4AF04 (ecoulements en milieux naturels) en salle 204 Tour 45/55 8h30-12h30 (5,/02/16 12/02/16, 19/02/16, 26/02/16 04/03/16)
planning : http://www.master.sdi.upmc.fr cours, notes.. http://www.lmm.jussieu.fr/lagree/COURS/MFEnv/index.html

1 Contexte
1.1 Importance pratique
La comprehension des ecoulements en milieu naturel est un enjeu important : enjeu scientifique (comprehension du monde qui nous entoure), industriel (navigation,
construction navale, energie) et enfin humain car la majeure partie des humains vivent le long des fleuves ou des cotes. (70 % des cotes sont en erosion, 80 % de la
population mondiale habite a basse altitude et plus de 20 % a proximite dun ocean ou dun estuaire. (source Wiki /C^ ote (geographie))). En France, 10% de la
population reside dans les communes pres des cotes, representant une densite de 2.5 par rapport a la moyenne metropolitaine. La longueur des cotes en France est
denviron 6000km dont le tiers est sableux.

Lobjet de ces notes de cours (le sujet est classique depuis Airy et Stokes et il existe des teraoctets de cours similaires sur le Ouaibe, sans oublier les kilometres de
rayonnages de livres a consulter) est de poser et de resoudre les equations regissant la propagation de perturbations de lelevation de la surface dun fluide dans un
champ de pesanteur. Cest a dire les vagues ! En faisant differentes hypotheses on simplifiera par lanalyse asymptotique les equations de Navier Stokes pour aboutir
a des systemes dequations aux derivees partielles plus simples. On rappelle dans un premier temps lequation donde, puis, apres avoir fait apparatre les effets de
propagation qui transportent la vague sans la deformer, on fera apparatre ceux de dispersion, qui la cassent en vaguelettes, enfin les effets des non linearites qui la
font deferler. Lorsque ces deux derniers effets se compensent des ondes qui se propagent sans changer de forme peuvent exister, ce sont les solitons. Les equations de
Saint Venant sont ensuite retrouvees dans le chapitre suivant dans une optique plus appliquee aux ecoulements dans les fleuves et les lagunes.

1
Houle et Vagues -2-

1.2 Modelisation
Attention, on va faire dans ce cours des simplifications importantes, cest ce que lon appelle la modelisation. Ces simplifications ne sont pas abusives, au contraire
elles permettent dextraire la quintessence du probleme en saffranchissant des epiphenomenes inutiles. Il faut trouver des equations modeles et des conditions aux
limites avec un niveau raisonnable de simplification... Cest toute la subtilite de la Mecanique et des sciences de lingenieur.

2 Introduction
Nous allons etablir les equations de propagation de la houle et des vagues. nous commencons par une presentation simple 1D. Nous presentons ensuite les equations
completes. Dans un premier temps, il nous faut definir ce quest une vague. La definition nest pas claire, tout comme la definition dun son et dun bruit nest pas
claire en acoustique (cest le cas pour plein dautres phenomenes de notre entourage).

Onde : perturbation qui se propage. ex : onde sonore (perturbation de la pres-


sion de lair), onde sur la corde vibrante (deformation dune corde tendue), ondes
electromagnetiques (perturbation des champs electrique et magnetique).
Vague : Perturbation de la surface de leau qui se propage.
Houle : Vague en profondeur infinie (loin de la cote et de la source).
Vitesse du fluide, vitesse de phase, vitesse de groupe : Il faut bien differencier la
vitesse de la particule de la vitesse de londe. Cette premiere est bien la vitesse
associee au deplacement dun petit volume deau, tandis que la seconde est la
vitesse a laquelle on a limpression que la vague avance. La troisieme est la
vitesse de deplacement du paquet dondes et de lenergie Figure 1 Une vague : une perturbation de la hauteur deau qui se deplace.
Deferlante : Cest une vague dont la crete se depace plus vite que le pied : elle se
casse ...
Mascaret : Cest en fait un simple ressaut hydrodynamique comme on en voit dans les eviers, mais il se deplace...
Soliton : Une vague non lineaire et dispersive qui se propage sans changer de forme...

2.1 modele simple en onde longue, ondes lineaires simples :


Nous commencons par un cas particulier qui permet de poser les concepts.
Dans cette premiere partie nous passons rapidement sur les problemes pour
aboutir a lequation dondes qui est la base de ce cours. Notre demarche
ici sera tres elementaire, fondee sur une modelisation 1D, dans les par-
ties suivantes nous preciserons les hypotheses faites. En fait, la presentation
suivante est peu rigoureuse (et meme fausse), mais nous ecrirons ensuite
toutes les equations de maniere rigoureuse. La maniere la plus simple
est de considerer une etendue deau initialement au repos et depaisseur
constante h0 . Supposons que nous perturbions cette surface, la hauteur devient
:

Cette elevation dune masse deau (x)x (par unite transverse au plan)
produit une variation de pression egale a g par nivellement barometrique. Figure 2 Le volume de controle pour etablir lequation de dAlembert a partir
La somme des forces agissant transversalement est egale a lacceleration de la dun modele simple (vitesse constante dans la section).
masse deau comprise dans le volume Sx, dou, puisque le fluide est pousse de
Houle et Vagues -3-

droite a gauche avec la force gS en x et de droite a gauche en x + x


g(x)S(x) g(x + x)S(x + x) = xSu/t.
On a suppose que la vitesse u est constante sur toute la section. En faisant tendre x vers 0, la conservation de la quantite de mouvement est :
u
= g (1)
t x

La conservation de la masse totale deau, appliquee au precedent petit volume donne par unite de longueur transverse, en introduisant la vitesse transverse v = t
(faible) secrit :
u(x + x)(h0 + (x + x)) + u(x)(h0 + (x)) vx = 0,
donc en faisant tendre x vers 0 :
(u(h0 + ))
= (2)
t x
puisque lamplitude est petite :
u
= h0 .
t x
Lelimination de u donne la fameuse equation donde (ou equation de dAlembert 1747) :
2 2
2 2 2 1 2
= c0 avec c0 = gh0 . Que lon note parfois a laide de loperateur appele dAlembertien :  = 0, avec  =
t2 x2 x2 c20 t2
Pour ecrire ces relations nous sommes alles un peu vite en besogne dans notre decompositions en tranches. Il est en effet plus correct de bien voir que le volume
que lon a choisi varie avec lecoulement : il se deforme. Il faudra aussi tenir compte des termes dinertie non lineaires. On pourra ainsi montrer que le raisonnement
plus haut permet de trouver les equations Shallow Water correctes telles que nous les etablirons plus loin (dans le chapitre suivant).

2.1.1 equation donde


Revenons a lequation de dAlembert :
2 2
2
= c20 2 avec c20 = gh0 .
t x
En posant = x c0 t, et = x + c0 t , on voit que comme

= + et = +
x x x t t t

= + et = c0 ( + )
x t
2

lequation devient alors par ce changement de variables 4c20 = 0.
Remarque, on peut ainsi ecrire
2 2
2

2
c0 2
= ( c0 )( + c0 )
t x t x t x
2

Les solutions de = 0 sont par integration en : 0
= cst, cette constante est constante par rapport a mais elle depend de , appelons la g (), ) puis = g()+cst,
mais cette constante est une fonction de , appelons la f (). Donc les solutions sont de la forme : f () + g() :
= f (x c0 t) + g(x + c0 t)...
La perturbation de forme f se deplace vers la droite sans changer de forme a la vitesse c0 (g se deplace vers la gauche).
Houle et Vagues -4-

2.2 Energie
Avant de compliquer on va ecrire la conservation de lenergie
u u 4
= g et = h0 .
t x t x 3.5
ce que lon ecrit 3
2.5
u u 2
uh0 = gh0 u et g = gh0
t x t x 1.5
1
on fait la somme des deux et on reconnait la derivee de 21 (h0 u2 + g 2 ) qui est
0.5
la somme dune energie cinetique et dune energie potentielle. En effet lenergie
potentielle dun element de -3 -2 -1 0 1 2 3 4
R masse dy a la hauteur y est gydy, lenergie du niveau
initial a la surface , est 0 gydy = 21 g 2 . On a :
Figure 3 Advection a celerite c0 dune perturbation a 4 temps differents . La
1 vague se deplace sans changer de forme, ici la solution f (x c0 t) est tracee ainsi
( (h0 u2 + g 2 )) = c20 (u)
t 2 x que quelques droites caracteristiques t = x/c0 + cst).
on a mis les equations sous forme conservative : variation denergie=flux. Pour
une onde allant de gauche a droite, = f (x c0 t), donc u = gf /c0 donc h0 u2 + g 2 = 2gf 2 et u = (g/c0 )f 2 , en substituant et en simplifiant par g la gravite,
il ne reste que :
2 2
(f )) = c0 (f 2 ) ou encore (f )) + c0 (f 2 ) = 0.
t x t x
Classiquement pour lequation des ondes non dispersives, on voit donc que
lenergie se deplace a la vitesse c0 , en effet on a une equation dadvection. Nous
verrons plus loin que lenergie ne se propage pas forcement a la meme vitesse que
le signal, lune est la vitesse de groupe, lautre de phase.

2.3 Premier exemple le Tsunami


Il sagit de l onde linearisee simple, la perturbation se deplace sur de tres
grandes distances en eau tres profonde, lequation precedente devient lequation
2
2
des ondes (le dAlembertien) : t2 c20 x 2 = 0 avec c = gh0 .
Le premier probleme est de determiner la source : cest en fait le fond qui bouge
tres vite et qui deplace la masse deau.
Or on voit que sur limage le tsunami a parcouru 32 degres en 5heures. 32 degres
correspondent a 32/90/26378 = 3500km divise par 5 heures, on trouve lordre
de grandeur : 720 km/h or (712/3.6)2 /9.81 = 4000m Cest bien la profondeur de
locean indien en moyenne. Figure 4 Gauche, calcul Saint Venant (voir chapitre suivant) effectue par le CEA
Remarque, un Mile est une minute dangle (60 minutes = 1 degres, le tour de de simulation du tsunami du 26/12/04 [click to launch the movie, Adobe Reader
la terre est 360 degres, et la circonference de la terre est 2Rterre ) : la longeur du required]. Droite, sachant que la profondeur moyenne est de 4km, la vitesse moyenne
Mille marin est 6378000 2 /(360 60) = 1852m la distance mesuree au bout est de 200m/s
de 5 heures est de 32 degres, donc 32*60=1920 Miles, soit en 5 heures : 1920/5=
384 noeuds, 1 noeud = 0.5m/s, on retrouve bien les 200m/s.
Pres des cotes, il faut rajouter des termes...
Houle et Vagues -5-

(0)
(0)
u(0)
u(0)
(0)
(0)
u(0)
x u(0)


(0)
x
u u(t)
u(0)
x

(t) (t)

Figure 5 Gauche : une perturbation initiale de vitesse nulle se casse en deux, une vague part a droite, lautre a gauche. Droite, une perturbation de vitesse et de
surface se deplace.

2.4 Une perturbation


Si on se donne une perturbation initiale = 0 (x) et u = 0, on voit que cette perturbation va se casser en deux : une qui va a droite, lautre qui va a gauche. on a
au bout dun temps petit t

u(t) = u(0) gt
x
donc, on cree une vitesse negative (x < 0) et ensuite positive. u
x a la forme de deux bosses de part et dautre de x = 0, ces deux bosses sont les corrections de (0)
donc la nouvelle hauteur deau presente 2 bosses :
u
(t) = (0) + (t)h0
x

En revanche, si on part dune perturbation se deplacant, gt x va corriger la position de la bosse de maniere a la faire avancer.
Houle et Vagues -6-

Figure 6 Mecanisme simple de la propagation dune vague. Gauche : examinons une vague a un instant t, par conservation du debit, la ou la vitesse augmente
la vague est etiree : elle samincit. En revanche la ou la vitesse diminue, il y a augmentation du niveau. La bosse de deplace sur la droite. Droite : Ensuite, la bosse
setant deplacee, la pression est maximale a son sommet et decrot en descendant a droite. La pression decrot, donc la vitesse augmente. En revanche dans la partie
gauche, la pression crot de gauche a droite (jusquau sommet) : le gradient de pression ralentit le fluide....
Houle et Vagues -7-

3 Resolution Numerique, un schema naf


Les equations 1 et 2 sont ecrites sans dimension, on se choisit une echelle de longueur L, et le
temps secrit L/c0
u
= et = (3)
t x t x
par Taylor, le plus simple est de faire un schema explicite :
u u
' (u(t + t, x) u(t, x))/t et ' (u(t, x + t) u(t, x x))/2/x
t x
in+1 = in t(ui+1 ui1 )/2/x et un+1
i = uni t(i+1 i1 )/2/x
ce qui se code
for(i=1;i<nx-1;i++) eta[i] = etao[i] - dt*(uo[i+1]-uo[i-1])/2/dx ;
for(i=1;i<nx-1;i++) u[i] = uo[i] - dt*(etao[i+1]-etao[i-1])/2/dx ;
puis dupdater

for(i=1;i<nx;i++) uo[i]=u[i];
for(i=0;i<nx;i++) etao[i]=eta[i];
Malheureusement cest instable ! Mais si on utilise le champ calcule dans la premiere equation
que lon injecte dans la seconde celle de u. Ce qui se code

for(i=1;i<nx-1;i++)eta[i] = etao[i] - dt*(uo[i+1]-uo[i-1])/2/dx ;


for(i=1;i<nx;i++) uo[i]=u[i];
for(i=1;i<nx-1;i++) u[i] = uo[i] - dt*(etao[i+1]-etao[i-1])/2/dx ; Figure 7 Haut cas instable t > x. Bas Cas stable, la bosse
for(i=0;i<nx;i++) etao[i]=eta[i]; se casse en deux ondes...

cela stabilise le schema tant que t < x. Cest un schema qui reste cependant tres naf. Nous
verrons des techniques plus elaborees ensuite.

Utiliser les codes onde1.c onde2.c et onde3.c pour reproduire la figure 7.


Houle et Vagues -8-

3.1 Ondes planes


Il est dusage de decomposer le signal en somme dexponentielles :

exp(ik(x c0 t)) = exp(ikx it)).

Cest en fait la transformation de Fourier du signal. est la pulsation ( en Hz, 1/s) k, la frequence spatiale (en 1/m), c0 = /k est la vitesse de phase, elle est constante.
Par la suite on va voir se qui se passe lorsque la vitesse de propagation du signal depend de la pulsation : c(k), nous intuitons que, puisque tout signal secrit sous
forme de serie de Fourier, si chaque mode se deplace a sa vitesse c(k), alors le signal se decompose. On dit que le milieu est dispersif.
Ce qui peut aussi se passer cest que c depende de la profondeur h, les equations sont alors non lineaires. On verra que le haut de la vague va plus vite que le bas.
Le signal se raidit, le haut de la vague depasse le pied dicelle et la vague deferle (dans la region du deferlement, nos equations ne sont plus valides, personne ne sait
de toutes facons resoudre les ENS dans ce cas excessivement severe).
Enfin, les frottements visqueux attenuent la vitesse, et donc la hauteur...

Figure 8 Deux deformations possibles du paquet donde. A gauche la dispersion, ce que nous allons voir ici. A droite, le raidissement non lineaire que nous verrons plus loin.
Houle et Vagues -9-

3.2 Reflexion sur un changement de fond


Ces equations de propagation des ondes suffisent a expliquer en premiere
approximation le deplacement du Tsunami. Nous allons voir quelle sont
les modifications pres des cotes. Ce qui suit est une premiere approxima-
tion.
ave
end of the second media dw
itte
plateau 2
x>0

nsm
tra Etudions la reflexion dun signal sur un changement de fond lorsque
la vague passe de locean au plateau continental en x = 0. On se
concentre sur ce point particulier et on ignore tout ce qui se passe au
junction between 1 and 2 loin ailleurs (ce qui est bien sur une approximation). Dans la mer, a
ave
gw gauche, h1 est la profondeur, sur le plateau a droite, la nouvelle profon-
c omin
i n refl deur est h2 , comme h2 est plus petit, la vitesse des ondes c2 = gh2
e cte
x<0

d wa est plus petite que la vitesse des ondes en pleine mer c1 = gh1 .
pleine mer 1

ve
Une vague damplitude a va etre a la fois reflechie et transmise au ni-
beginning of first media
veau du changement de profondeur. Ecrivons les relations de transfert/
time reflexion.

On va ecrire que londe reflechie (damplitude R a) se superpose a londe


Figure 9 La discontinuite de la pleine mer a un plateau, les vagues transmises sont
incidente en mer et que pour x > 0 il ny a que londe transmise damplitude
amplifiees.
T a

(x < 0, t) = a exp(i(k1 x k1 c1 t) + R a exp(i(k10 x k10 c1 t)) et (x > 0, t) = T a exp(i(k2 x k2 c2 t).


On ecrit que lamplitude de la vague est continue en x = 0 soit :
exp(i(k1 c1 t) + Rexp(i(k10 c1 t)) = T exp(i(k2 c2 t)
pour que ce soit toujours vrai, on voit que la frequence doit rester la meme = k1 c1 = k10 c1 = k2 c2 , ce qui veut dire que k10 = k1 , et k1 c1 = k2 c2 , donc 2 /1 = c2 /c1 ,
or comme la vitesse est plus faible quand la hauteur est plus petite c2 < c1 . La longueur donde diminue.
On en deduit
1 + R = T.
Lexpression de la vitesse, puisque k10 = k1 , la vitesse est obtenue par t u = gx
ga gR a T a
u(x < 0, t) = exp(i(k1 x k1 c1 t) + exp(i(k1 x k1 c1 t)) et u(x > 0, t) = exp(i(k2 x k2 c2 t).
c1 c1 c2
puis en ce meme point x = 0 on ecrit que le flux total deau est conserve, soit la hauteur multipliee par la vitesse donc
1R T
h1 = h2
c1 c2
On en deduit les coefficients de reflexion de de transmission dune vague au travers du changement de section
p
1 h2 /h1 c1 c2 2 2c1
R= p = , et T = p =
1 + h2 /h1 c1 + c2 1 + h2 /h1 c1 + c2
On voit donc quune vague est amplifiee lorsquelle franchit une remontee du fond. En effet 2 > R > 1 et 1 > T > 0 pour h2 /h1 < 1 et en revanche pour h2 /h1 > 1
(passage a un fond plus profond) on a 1 > R > 0 et 0 > T > 1. Cet effet est independant de la pulsation ou de la longueur donde.
Houle et Vagues -10-

2.0

1.5

1.0

R,T
0.5

0.0

0 1 2 3 4
h_2h_1

Figure 10 Coefficients de reflexion et de transmission R courbe du haut et T courbe du bas en fonction de h2 /h1 . Pour h2 /h1 < 1 on a 2 > R > 1 et 1 > T > 0 et en revanche
pour h2 /h1 > 1 (passage a un fond plus profond) on a 1 > R > 0 et 0 > T > 1.

3.3 Cas dun paquet dondes


Considerons un paquet dondes constitue dune somme de vagues sinusodales, par exemple

e16(t/2) = k=
k=1 ak sin((2k 1)t) =

= 0.197919 sin(t) 0.185926 sin(3.t) + 0.16407 sin(5.t) 0.135996 sin(7.t) + 0.105873 sin(9.t) + ...
cela nous permet de decomposer le signal et de tracer la reflexion/ transmission de la figure suivante.

1.2 1.2

1.0 1.0

0.8 0.8

0.6 0.6

0.4 0.4

0.2 0.2

0.0 0.0
-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0
x x

Figure 11 Une vague unique arrive a une discontinuite de fond (le fond remonte brusquement), elle est transmise et amplifiee, reflechie et attenuee. A gauche,
londe a passe la discontinuite, on voit londe transmise/ reflechie en noir, en rouge londe initiale telle quelle serait sil ny avait pas eu de discontinuite. A droite, la
meme chose en film [click to launch the movie, Adobe Reader required]
Houle et Vagues -11-

3.4 Un cas plus complique : remontee douce dun fond

0.4

0.2

J0x^12
0.0

x=L h0 0.2
x=0
0.4
0 200 400 600 800 1000
x

Figure 12 Elevation de leau a larrivee dune plage, londe vient de droite a gauche. Lamplitude augmente graduellement et la longueur donde diminue.

Si le fond remonte doucement (inspire de Billing & King p115), on peut en premiere approximation dire que h(x) est lentement variable par rapport a la longueur
donde et que les equations sont
u h(x)u
= g et = .
t x t x
la deuxieme est 2 non simplifee. en eliminant u des deux equations, (t h(x) = 0)

2
= g [h(x) ].
t x x
Si londe vient de la droite et rencontre une plage On se donne
x h0
h(x) = h0 pour 0 < x < L et h(x) = h0 pour L < x. Avec =  1.
L L
et donc
2 2
= gx 2 + g .
t x x
on se donne la perturbation initiale de la forme
(L, t) = acos(t + )
et cherchons des solutions de la forme = H(x)cos(t + ) en substituant, on voit que lamplitude verifie :

d2 H dH 2
xx 2
+ + H = 0.
dx dx g

posons x = 2s2 , dH/dx = (1/(4s))dH/ds et d2 H/dx2 = 1/(16s2 )d2 H/ds2 1/(16s3 )dH/ds

d2 H dH 8 2
2
+ + H = 0.
ds sds g
Houle et Vagues -12-

m
(1)
Or on sait que les equations de la forme y 00 + y 0 /x + (1 n2 /x2 )y = 0 sont les solutions des equations de Bessel de la forme Jn (x) = x 2m+n
m=0 m!(m+n+1) ( 2 ) On en
deduit la hauteur deau :
J0 ( 2
g
x1/2 )
(x, t) = a
J0 ( 2
g
L1/2 )
on constate donc que lamplitude de la vague augmente, que la la longueur donde diminue, tandis que la vitesse de londe diminue au fur et a mesure que lon se
rapproche de la cote (et ce de maniere continue, mais identiquement au cas abrupte pour ce qui est des conclusions). Bien entendu, ce modele nest plus valide pres
de la cote car les vagues finissent par deferler : nous verrons plus loin pourquoi. Neanmoins, ce raccourcissement des longueurs dondes et cette amplification semblent
etre observes en particulier pour les tsunamis.
Cet exercice est quelque peu inexact car londe est stationnaire, on a un developpement plus complet dans :
http://wikhydro.developpement-durable.gouv.fr/index.php/ANSWER - Propagation dune onde dans un estuaire a pente du fond inclinee

Figure 13 Elevation de leau a larrivee dune plage dans le canal a Houle de lENSTA (photo PYL), londe vient de la droite. On voit que lamplitude augmente legerement
vers la plage et que la longueur donde diminue aussi a cet endroit. On constate que le modele de la figure precedente tombe en defaut lorsque la profondeur est trop faible...

3.5 Cas en axi


Dans le cas axi, lequation de conservation de quantite de mouvement :
u
= g
t r
et celle de la masse change un peu car la masse se conserve circulairement
hru
=
t rr
on a donc
2
2
= gh [r ]
t rr r
Houle et Vagues -13-

ou
2 1 2
2
+ 2 2 = 0
r r r c0 t
pour une onde it
2 1 2
2
+ + 2=0
r r r c0
or J000 (x) + J00 (x)/x + J(x) = 0, avec
x2 x4 x6
J0 (x) = 1 + + ...
22 22 42 22 42 62

1.0 2.0

0.8
1.5

0.6
1.0
0.4
J0 HxL

J0 HxL
0.5
Out[103]= 0.2

0.0
0.0

-0.2 -0.5

-0.4 -1.0
0 5 10 15 20 0 2 4 6 8 10
x x

x2 x4 x6
p
Figure 14 fonction de bessel J0 et ses approximations 1 22
+ 22 42
22 42 62
et 2/(x) cos(x /4) ...

Reprenons maintenant ces notions en compliquant et en ecrivant enfin les equations completes.
Houle et Vagues -14-

4 Equation generales completes


4.1 Navier Stokes complet
Reprenons le probleme a zero, on a deux fluides lun au dessus de lautre (air et eau), dans un champ de gravite. Les equations a resoudre sont a priori pour le
probleme de lecoulement des fluides sont pour chaque phase : l incompressibilite

u=0

la conservation de la quantite de mouvement


(u/t + (u )u) = p gey + 2 u.
Ces systemes dequations sont a resoudre dans les deux milieux que sont lair et leau. Remarquer a ce propos lequation de conservation de la masse secrit bien :

/t + (u) = 0 <=> u = ((/t + (u )).

La densite est conservee le long des lignes de courant : leau reste de leau, et lair reste de lair. A la surface : cest delicat ! En toute rigueur, il faut exprimer la
continuite des contraintes normales et tangentielles ainsi que des vitesses des deux fluides. On ecrit la continuite de la contrainte normale : soit n la normale exterieure
a la surface libre :
(pe pa K)ni = ij nj aij nj
avec pour lair aij = 2a Dij , et pour leau ij = 2Dij , ou tension superficielle et K courbure de la surface. Comme la viscosite dynamique du gaz est faible par
rapport au liquide, on neglige le mouvement de lair reste :
(p p0 K)ni = ij nj
avec p0 pression atmospherique. Cette expression est a projeter sur la normale et la tangente. Nous necrivons pas pour linstant ces equations developpees. Nous
simplifions tout de suite au cas non visqueux. Le but du cours est de passer dequations 2D a des equations 1D plus simples (le cas general 3D donne des systemes 2D,
la generalisation va de soit).

4.2 equations dEuler :


En pratique la viscosite est negligeable. Les equations a resoudre sont :
i) lincompressibilite
u=0
sous forme developpee :
u v
+ = 0,
x y
ii) la conservation de la quantite de mouvement :
u
( + (u )u) = p gey
t
sous forme developpee :
u u u p v v v p
( +u +v )= et ( +u +v )= g
t x y x t x y y
iii) accessoirement irrotationnel (pas de creation de vorticite)

u v
u = 0, sous forme developpee : + = 0.
y x
Houle et Vagues -15-

et ce avec les conditions aux limites ad hoc : a la surface on a pour la pression :

p = p0 + n avec n = (1 + ()2 )1/2

la pression est egale a la pression atmospherique, au saut de pression pres du a la tension de surface. condition dynamique : le saut de pression est tel que :

p1 (x, (x, t), t) p2 (x, (x, t), t) = n12

la normale a la surface est :


n12 = (/x, 1)/(1 + (/x)2 )1/2
condition cinematique :
Les deux fluides ne decollent pas de la surface qui les separe et ne sinterpenetrent pas : intuitivement, les vitesses normales sont donc continues a la traversee de
linterface :
u1 n12 = u2 n12 en y (x, t) = 0.
On retrouve cette meme relation a partir de u = 0 integree sur une petite epaisseur (cf. electromagnetisme ! Voir la condition sur la composante normale du champ
magnetique B a une interface !).
Une autre maniere de voir les choses est de considerer linterface entre leau et lair. Leau reste eau et lair reste air, ce que nous dit la conservation de la masse
(voir plus haut)

u = ( + (u )).
t
dans chacun des deux fluides, lecoulement est incompressible donc u = 0, et il reste que d dt = 0 ; la densite est transportee dans lecoulement, sa derivee totale
est nulle. Nous allons exploiter cette idee en appelant F (x, y, t) = 0 la fonction qui est nulle a linterface entre les deux fluides. F (x, y, t) = 0 est la forme implicite de
lequation de linterface ecrite sous forme explicite y = (x, t). On a donc

F (x, y, t) = y (x, t) = 0
dF
cest lequation de linterface, on remarque egalement que linterface etant materielle, F est transporte le long des lignes de courant : dt = 0 ; (derivee totale) Par
definition de la derivee lagrangienne dy d
dt = v(x, , t) et dt = t + u x , donc dF/dt = 0 secrit :


v(x, , t) = (x, , t) + u (x, , t)
t x
Remarque
quon peut reecrire en faisant apparatre les vitesses normales

= u n(1 + (/x)2 )1/2
t
Pour la vitesse sur le fond on a (en y = h0 ) : on a le glissement : u n = 0. v = u/x.
En resume :
u = 0, (ut + (u )u) = p ey . et u = 0 a la surface : p = p0 + .n;
v() = t + u/x
sur le fond on a v = u/x.

Nous commencons par examiner quelques exemples simples de solutions de


ces equations.
Figure 15 la surface libre , et un fond .
Houle et Vagues -16-

5 Equation de la Houle
5.1 ordres de grandeur
5.1.1 Premieres hypotheses :
La Houle (swell en anglais) est ce train de vagues que lon observe en mer. Le
modele que nous presentons est du vraisemblablement a Airy 1845. Il a ete aussi
presente par Stokes qui a pousse le calcul aux ordres suivants.
Reprenons les hypotheses precedentes et examinons les (voir Landau Lifshitz
12, voir Lamb Chpat IX). Premierement, nous oublions la viscosite en supposant
que le nombre de Reynolds est grand. La surface dequilibre est plane. Letat de
base est donc un fluide immobile. Il ne reste que la stratification de pression :
u = v = 0 et 0 = p/y g. Le fond du recipient est en y = 0, la surface libre
en y = H ou il regne la pression atmospherique, donc :
p = p0 + g(H y).
(La pression double chaque fois que lon senfonce de 10m). La perturbation est
legere, soit a son amplitude, et soit la periode des oscillations. En effet si on
observe un fetu de paille a la surface de la mer, ou le bouchon de la ligne de canne
Figure 16 la surface libre , et un fond plat, des vagues de longueur donde . a a peche, on constate quil navance pas avec les vagues mais revient a la meme
amplitude des deplacements. H est la hauteur deau au repos. place. La vitesse est a/ . Soit la longueur donde, le terme non lineaire u/x
a pour ordre de grandeur par rapport au terme instationnaire /t :
(u/x)/(/t) ((a/ )/)( ) = a/,
donc on peut lineariser si a/  1. On a aussi a/H petit : on sinteresse a de petites perturbations, mais H/ quelconque.
On a donc a/ << 1 et a/H << 1 mais H/ quelconque (plus tard on posera = a/ et = H/ , ici est petit est quelconque).

5.1.2 relation de pression :


La pression dequilibre est :
p = p0 + g(H y)
On pose p0 + g(H y) la valeur hydrostatique (de repos), on etudie les variations autour de cette valeur : p0 + g(H y) + P La pression est donc ecrite sous la
forme :
p = p0 + g(H y) + P.
a linterface elle secrit exactement :
p(x, y = H + ) = p0 . donc p(x, y = H + ) = p0 = p0 + g(H H ) + P (x, y = H + )
On en deduit que, de maniere exacte :
P (x, y = H + ) = g.
Lordre de grandeur de P est donc bien ga, et P (x, y = H + ) = g. La partie active de la pression est bien reliee aux variations de linterface :
p P
= = g
x x x
On remarque que lon va lineariser cette relation en
P (x, y = H) = g
Houle et Vagues -17-

5.1.3 Relation a linterface


Reprenons cette estimation avec les equations Euleriennes : la vitesse de deplacement de linterface est :

v(x, , t) = (x, , t) + u(x, , t) (x, t), linearisee en v(x, 0, t) = (x, t)
t x t
donc (a/ ) est bien la jauge de v. En ordre de grandeur lincompressibilite : u/x + v/y = 0, donne v et u de meme jauge, a condition de mesurer toutes les
directions avec la meme longueur (ce qui est assez naturel).

5.1.4 equations du mouvement


Si letat de repos est perturbe, pour la vitesse transverse il ne reste que P , en effet lequation complete :
v v v p
( +u +v )= g
t x y y

devient par definition de P , et par linearisation (si a/ << 1) cela donne :

v P
=
t y

comparons en ordre de grandeur (v/t) a P/y, cela donne respectivement a/ 2 et ga/. Cela donne 1/ 2 = g/ . On remarque que lacceleration transverse
est negligeable en ordre de grandeur par rapport a lacceleration de la pesanteur : a/ 2 << g car 1/ 2 = g/ = (g/a)(a/) << (g/a) car (a/) << 1. Lordre de
grandeur de la variation de pression P etant estime a ga on doit aussi garder
u P
=
t x
qui est automatiquement verifie en ordre de grandeur : 1/ 2 = g/ Remplacons par sa valeur (/g)1/2 dans a/ :
p
a/ (a/)(a/H)c0 ou c20 = gH

or a/H petit et a/ << 1. La vitesse est plus petite que c0 qui est la vitesse naturelle (par analyse dimensionnelle) associee a un fluide au repos de profondeur H.
H aurait dailleurs pu etre considere comme une longueur naturelle. Ici on a choisi et le temps associe obtenu pendulaire (/g)1/2 , ce qui semble raisonnable.
le deplacement de linterface

v(x, y = (H/ + (a/)) = /t + (a/)u/x


va donc se lineariser (terme non lineaires negligeables), et autour de y = H on va aplatir la condition aux limites.
la condition de pression a linterface linearisee en (a/) va nous fournir une relation liant P et . En effet :

P (x, y = (H/ + (a/))) = P (x, y = (H/)) + (a/)P (x, y = (H/))/y + ... = P (x, y = (H/))

Donc au final :
P (x, y = (H/)) = .
Houle et Vagues -18-

5.1.5 resumons nous


mesure par a
v mesure par a/
u mesure par a/
x mesure par
y mesure par
P mesure par
gaH = gH(a/H) = c20 (a/H)
p = p0 + g(H y) + P = p0 + c20 (1 y/H + (a/H)P ) ou encore, p mesure
+ 1 y/H
avec gH : p = p0 /(gH) p + (a/H)P
t mesure par et = H/c0 = /g.
Par substitution directe :
u v
+ =0
x y
u P
=
t x
v P
=
t y

v = 0 au fond, et en y = H : v = en haut et P = g en haut. ce qui donne
t
par elimination de la vitesse par incompressibilite :

2P 2P
+ =0
x2 y 2
v
et = P
x avec v = t et P = g en y = H et
P
y = 0 en y = 0.
t

5.1.6 Resolution en modes normaux


On cherche des solutions en F (y)ei(kxt) (car les equadiffs sont invariantes
en t et x), les conditions aux limites portent sur y. La solution de

Figure 17 lobjet flotte entre deux eaux et tourne en rond. Photo issue de Billingham 2P 2P
&King 2
+ =0
x y 2

est donc en P = eky ei(kxt) . Substituons :

P = P+ eky ei(kxt) + P eky ei(kxt)

donne
k
v= (P+ eky ei(kxt) P eky ei(kxt) )
i
Houle et Vagues -19-

soit au fond v(x, 0) = 0 => P+ = P donc


kP+ ky 2kP+
v= (e eky )ei(kxt) = i (sinh(ky))ei(kxt)
i
par lincompressibilite on obtient u : u = 2kP
(cosh(ky))e
+ i(kxt)
On remarque que lecoulement est irrotationnel. En general, pour une houle lineaire sans courant, il
ne peut y avoir de rotationnel puisquil ny a pas de transport (on rappelle que est conserve le long des lignes de courant)...

Si on pose
= (2P+ //g) cosh(kH)ei(kxt) = 0 ei(kxt)
donc
kg kg
u= cosh(ky)/cosh(kH), v = i sinh(ky)/cosh(kH), P = g cosh(ky)/cosh(kH)

a la surface : v = /t = (P/t)/(g) donc on a :
i(k/) sinh(kH) = i cosh(kH)
Donc on obtient la fameuse relation de dispersion :
2 = gk tanh(kH)
Retenons donc la relation de dispersion des ondes lineaires de profondeur quelconque. en ordre de grandeur : 1/ 2 c20 //H est verifie comme nous lavions suggere.
On definit la vitesse de phase c = /k et donc
r r
2 g g g H
c = tanh(kH) c = tanh(kH) en terme de longueur donde c = tanh(2 )
k k 2
Sur le tableau ci dessous on donne la vitesse en fonction de la lon-
3 2
gueur donde (on bonne approximation pratique est dans [4]) = T
2
(metres et secondes) et V = 3T (en noeuds 0.5m/s et en se-
condes).

(m) 0.1 1 10 100 1000


c (m/s) 0.4 1.25 4 12.5 40
T (s) 0.25 0.8 2.5 8 25
Les premiers signes dune depression forte sur latlantique sont larrivee
sur les cotes de lOuest de lEurope dune houle de longueur donde tres
longue et de hauteur faible. Leur periode est typiquement de 30s et leur vi-
tesse de 47m/s. Elle parcourt environ 4000km/jour et donc elle arrive bien
avant la tempete qui va mettre plusieurs jours a se deplacer. Plus elle avance,
plus la longueur donde diminue car les vagues se sont deplacees plus lente-
ment.

5.1.7 2 cas limites Figure 18 Lorsque le fond diminue, la vitesse diminue Sur ce graphe suivant on
trace en fonction de H la profondeur a fixe. la vitesse est alors c = . Si on se
Si le fond devient de plus en plus profond : donne exemple une houle oceanique de periode 10s (iT = 2/=10s). On trace la courbe
implicite = (g/(2 2 )tanh(2H/)) Au 2
p loin, = (g/(2 )), on voit quelle est de plus
r r
2 g g
= gkth(kH) gk c = en terme de longueur donde c = 2
en plus courte pres de la cote car = (g/ )H. La vitesse decrot aussi. Ceci permet
k 2
dexpliquer pourquoi les vagues sont toujours paralleles a la plage !
Houle et Vagues -20-

H na plus de raison dintervenir, et effectivement, il disparat : 1/ 2 c20 //H


est en realite 1/ 2 g/. La vitesse c0 est un artifice de calcul. Sur le graphe
suivant on trace en fonction de H a fixe. Si on se donne exemple une houle
oceanique de periode 10s, on voit quelle est de plus en plus courte pres de la cote.
La vitesse 2 decrot aussi. Ceci permet dexpliquer pourquoi les vagues sont
toujours paralleles a la plage ! (la courbe suivante 18 montre la longueur donde fonction de H lorsque lon se donne , cette courbe serait a modifier pres de nul :
la tension de surface modifie alors la relation de dispersion : ce sont les ondes capillaires de longueur donde environ egale a 1.7cm, vitesse 23cm/s).

5.1.8 si le fond est peu profond


Si maintenant le fond est peu profond H/ petit

2 = gkth(k(H)) gk 2 H
p
c = /k = gH ne depend pas de la longueur donde

avec H/ petit, a k fixe devient de plus en plus petit. Le milieu nest plus dispersif. Les vitesses et pressions

kg g c0 kg g y
u= cosh(ky)/cosh(kH) = et v = i sinh(ky)/cosh(kH) kH(y/H) = , et P = g cosh(ky)/cosh(kH) g
c0 H c0 H t
La vitesse u et la pression ne dependent plus de la profondeur, examinons maintenant la vitesse transverse. On constate que la vitesse transversale devient plus
faible que la vitesse longitudinale, le rapport de ces deux vitesses (v/u) est dordre (H/) << 1. Le milieu nest plus dispersif... On retrouve le modele de la section 1
ou la vitesse ne varie pas sur lepaisseur.

= (Hu) et u = g
t x t x
Houle et Vagues -21-

5.1.9 trajectoires
Pour tracer les trajectoires, :

x = X0 + X1 (t) y = Y0 + Y1 (t)

ou X1 et Y1 sont de petits deplacements autour de X0 et Y0 , au premier ordre on


peut confondre la vitesse Eulerienne et Lagrangienne donc
kg
X10 = u = 0 cosh(kY0 )/ cosh(kH) cos(kX0 t)

kg
Y10 = v = 0 sinh(kY0 )/ cosh(kH) sin(kX0 t)

avec par exemple X1 (0) = Y1 (0),
kg
X1 = 0 cosh(kY0 )/ cosh(kH) sin(kX0 t)
2
kg
Y1 = 2 0 sinh(kY0 )/ cosh(kH) cos(kX0 t)

par elimination du temps :
g
(X1 /[cosh(kY0 )] 0 sin(kX0 ))2 +
kc2 cosh(kH)
g g
+(Y1 /[sinh(kY0 )] 0 cos(kX0 ))2 = (0 2 )2
kc2 cosh(kH) kc cosh(kH)
Les particules de fluide tournent dans le sens des aiguilles dune montre le long
dune petite ellipse avec la periode de londe. Figure 19 source http://fr.wikipedia.org/wiki/Fichier:2sec 3mdepth puv drift.png
Au fur et a mesure que lon sapproche du fond, la trajectoire devient de plus en plus plate.
Attention, on mesure la position a partir du fond, donc si le fond est profond, il est plus judicieux dutiliser la distance a la surface H z, donc Y0 + z = H
kg kg
X1 = 0 cosh(k(H z))/ cosh(kH) sin(kX0 t) et Y1 = 2 0 sinh(k(H z))/ cosh(kH) cos(kX0 t)
2
on voit mieux que si kH est grand, cosh kH ' sinh kH ' ekH donc le rayon est kg2 0 ekz
On trace ici le dessin des trajectoires... Que lon peut voir en animation sur : http://www.lmm.jussieu.fr/~ lagree/SIEF/SIEF97/HouleA.mov
En eau tres profonde les trajectoires sont des cercles de rayon decroissant rapidement. En eau moins profonde, ce sont des ellipses, le petit axe diminue jusqua
zero, le grand axe reste a peu pres constant.
Houle et Vagues -22-

Figure 20 En eau tres profonde les trajectoires sont des cercles de rayon decroissant rapidement. En eau moins profonde, ce sont des ellipses, le petit axe diminue jusqua zero,
le grand axe reste a peu pres constant. http://www.lmm.jussieu.fr/~ lagree/SIEF/SIEF97/HouleA.mov

Figure 21 Sous la houle, visualisation de trajectoires de bulles en pose longue, les trajectoires des ellipses, le petit axe diminue jusqua zero, le grand axe reste a peu pres
constant. photo issue du Van Dyke.

Surface libre de la houle d'Airy


X: position d'une particule

Figure 22 A gauche, lobjet flotte entre deux eaux et tourne en rond. Photo issue de Billingham &King. A droite Trajectoire circulaire a la surface de leau [click
to launch the movie, Adobe Reader required], le film est aussi sur http://www.lmm.jussieu.fr/~lagree/SIEF/SIEF97/HouleA.mov
Houle et Vagues -23-

5.2 Ecriture en potentiel


Le fait que la houle soit irrotationnelle fait que lon ecrit souvent les equations en fonction potentiel (notons que lon avait pas besoin de cette hypothese dans le
2 2
cas linearise precedent). On a u =
x et v = y . lincompressibilite donne x2 + y 2 = 0. La condition de surface libre secrit alors a laide de lequation de Bernoulli

1 2 2
+ ( + ) + g = 0
t 2 x y
cette equation se linearise en y = H a la surface live aplatie :

+ g = 0
t

La condition dynamique a la surface libre secrit y = t + x x . Elle se linearise en

y =

t . Ces conditions linearisees a la surface aplatie :
t + g = 0 et
y =
t
on les combine, en y = H :
2
2
+g =0
t y

et on a toujours en y = 0 le gisement y =0
Le systeme a resoudre est donc

2 2 2
2
+ 2 = 0, + g , en y = H et en y = 0 on a =0
x y t2 y y
Solution
= 0 ei(kxt) , avec 2 = gk tanh(kh)
ig cosh(ky) i cosh(ky)
= 0 = 0
cosh(kH) k sinh(kH)
on remarque la double ecriture de car /k = (g) tanh(kh)

5.3 Non linearites : houle de Stokes


5.3.1 Cas de loscillateur simple
prenons un oscillateur simple avec un terme non lineaire

t2 + 02 + 2 = 0
on sinteresse a une perturbation dode petit donc on pose = 0 + 2 1 + ..., par substitution et identification des puissances a lode 0 on a le probleme linearise
t2 0 + 02 0 = 0
de solution 0 = A0 cos 0 t, lordre suivant satisfait :
t2 1 + 02 1 + 20 = 0
or cos2 0 t = (1 + cos 20 t)/2, donc la partie variable est telle que
t2 1 + 02 1 + A20 cos(20 t)/2 = 0
dou 1 = (A0 /6) cos 20 t), le developpement de la solution finale est donc
= A0 cos 0 t + 2 (A0 /2 + (A0 /6) cos 20 t)...

5.3.2 Cas de la houle


Houle et Vagues -24-

Avec la presentation avec le potentiel , on a un avantage, cest que le probleme


secit :
2 2 1 2 2
+ = 0 A la surface la relation exacte est + ( + ) + g = 0
x2 y 2 t 2 x y
Notre presentation precedente, linearisait tout de suite les equations, puis on
obtenait un laplacien pour la pression. En fait la demarche en potentiel est plus
interessante, car la non linearite est rejetee sur la surface.
Ensuite, on a dit que cette equation se linearise en

+ g = 0
t
dont la solution sera en cos(kxt), mais dans les termes negliges, il y a des carres Figure 23 comparaison Code SWASH http://swash.sourceforge.net/ et solution
2 2
x , donc des termes en cos (kx t), qui seront en terme source (petits). Ce de la houle de Stokes
terme source du second ordre en cos2 (kx t) secrit aussi en cos(2(kx t)), car
2 cos2 1 = cos 2. La non linearite dans la surface libre va faire ainsi apparatre
un terme source de phase double qui secrit apres calcul :

0 3 tanh2 kH
= 0 (cos(kx t) + cos(2(kx t) + ...)
2 tanh3 kH
3tanh2 kH
leffet est daplatir le fond et de de rendre la crete plus fine. Le nombre 4 tanh3 kH
est appele le nombre dUrsell.

5.4 Retour sur les limites


Si on estime le terme relatif non lineaire :
uu uk
u/c = k0 /tanh(kH)
t u
en onde peu profonde, cest 0 /H
0 3tanh2 kH 0 3 0 2
Le terme non lineaire de Stokes 2 tanh3 kH est 4 k2 H 3 H3 ne doit pas etre trop grand, sil est trop grand on pousse a lordre suivant.

5.5 Dispersion
5.5.1 Capilarite
(la courbe suivante montre la longueur donde fonction de h lorsque lon se donne cette courbe serait a modifier pres de nul : la tension de surface modifie alors
la relation de dispersion : ce sont les ondes capillaires de longueur donde environ egale a 1.7cm, vitesse 23cm/s).

3
= (Hu) et u = g + 3
t x t x x
Houle et Vagues -25-

5.5.2 Effet faible de la profondeur


q
gk 2 H 13 gk 4 H 3 + ... pour les ondes longues, et
p
Nous venons de voir par letude de la houle que = gk tanh(kH), un developpement limite est =
3 1 2 3
= gHk(1 16 k 2 H 2 + ... Pour une onde eitikx , on identifie avec i t
et k avec i x
donc t
= c0 ( x 16 H 2 x
3 ) + ... ou encore t = c0 x + 6 H x2 t + ...,

au lieu decrire

+ (Hu) = 0, et u = g
t x t x
on va ecrire
1 2
H2 2 + (Hu) = 0 et u = g
t 3 x t x t x
2
ce systeme simplifie a pour relation de dispersion 2 = 1+gk1 kH2H2 , qui a le bon gout de redonner la relation de dispersion 2 = gk 2 H(1 13 k 2 H 2 + ...). Cest donc une
3
bonne approximation ainsi que la forme suivante
1 2
(Hu) = et u H 2 2 u = g
t x t t 3 x t x
Ce sont les equations de Boussinesq. Ces equations conduisent a la propagation du soliton et du mascaret (cf plus loin...)

5.5.3 Effet faible de la profondeur : solution lointaine


On a trouve des solutions en ondes telles que 2 = gkth(kh0 ) si on developpe dans le cas peu dispersif de leau peu profonde, a grande longueur donde kh 0,
on fait un developpement limite en kh0 de la tangente puis de la racite, (et c20 = gh0 ) :
(kh0 )3 (kh0 )2 (kh)2
2 = gk(kh0 + ....) = (gh0 )k 2 (1 + ....), = c0 k(1 + ....)
3 3 6
(h0 )2 2
on en deduit en prenant la valeur + des ondes qui se deplacent vers la droite. On la forme donde = 0 exp(i(t kx)), avec i = (ik)c0 (1 + 6 (ik) + ....) ce qui
veut dire que puisque t = i et que x = ik lequation de dispersion linearisee est celle correspondant au probleme suivant
c0 h20 3
+ c0 + = 0.
t x 6 x3
Par invariances par dilatations on obtient des solutions semblables... en effet la solution de

=
t x3
3 1/3
R
par similitude T = X et HX = 1 = t f ( t1/3 ) avec i = 1 dou f = f 0 /3 f /3 soit f 00 = f /3 donc la solution est
x 000

x
= 31/3 t1/3 Ai(31/3 1/3 )
t
Donc, la solution est :
(x c0 t) 1/3
(x, t) = t1/3 Ai[ (2) ]
t1/3
Houle et Vagues -26-

Figure 24 Onde dispersive

WATER WAVES 7395

STATIONARY- PHASE SOLUTION


L LEADING-WAVE
SOLUTION
O,I
ft EXPERIMENTAL
DATA
x%15
ft,X=2ft,
d=2.3ft
p <0 --/" --- x'/d=6.525,X/d=0.87
h=-O.3
ft
, / /f .,-,
- // k // k i/ ',? ',i / / ,, '--,

2<o
. _

TIME' f (seconds)
0 I 2 3 4 5 6 7 8 9 I0
I I I I i
I I I I, I
0 2 4 6 1 I0 12 It 16 IS 20 22 24 26 28 30 ;52 ;54 36
TIME
' t=t (dimensionless)
Fig. 6. Experimental and analytical amplitude-time relationships.

[REGION WHERE ASYMPTOTIC METHODS and right casesrespectively.Sinced -- 2.3 feet


Figure 25 Comparaisons issues de Noda JGR 76ARE VALID
(30) 1971, Water
wavesthen X/d -- 0.87
generated byanda Figure
local 2surface
indicates the
disturbance. il compare, des experiences, aux
Figure 2 shows a comparison of the i/4 same
solutions 00 1/2envelopecharacteristics, type C, as shown
deux approximations proposees, celle de la phase stationnaire i 2e (by t) F (k 0 )exp(i(k 0 x 0
0 the data with the first wave the largest. t)) que lon va voir plus loin a celle de la solution
1 1/3 xc0
obtainedby asymptoticapproximationsand by
1/3 The middle cases are for X -- 1.0 foot
= h2 c0 t/2 Ai(2 (h2 c0 t)1/3 ) que lon vient de numerical
direct voir integration.It is a surprising
result that the asymptotic solutionsprovide (,Vd - 0.435), h -- -----0.1foot and Figure 2
0 0

goodapproximationto (7) for quitelargevalues confirmsthat the envelopecharacteristicsare


of A,/d and small valuesof x*/d. of type A where the first wave is not the
Figure 2 also indicatesthe generalshapeof largest. The bottom casesare for X -- % foot
(X/d -- 0.145), h ---- -----0.2
foot and again, Fig-
Houle et Vagues -27-

6 Vitesse de groupe
6.1 Energie de la houle
Energie de la Houle dans le cas profondeur infinie

= Aei(kxt) , = Beky ei(kxt) , B = iAg/ 2 = gk

en repassant en reeles
= 0 cos(kx t)eky = 0 sin(kx t) 0 = (k/)0 .
et evaluons les moyennes, dans un premier temps la moyenne spatiale (sur une longueur donde) :
Z
< sin2 kx >=< cos2 kx >= 1 sin2 (2x/)dx = 1/2
0

La moyenne dun cos au carre est 1/2. En revanche on definit la moyenne sur la profondeur par
Z 0
< e2ky >= e2ky dy = 1/(2k)

Ayant fait ces remarques evaluons lenergie totale dune tranche de longueur de fluide :

* energie potentielle,
elle sobtient en remarquant que son accroissement est proportionnel a la variation de hauteur dy dune masse (y), dou :

Ep =< g 2 /2 >

En moyenne sur une longueur donde : g02 /4.

* energie cinetique :

Ec =< (u2 + v 2 ) > /2 = (1/4)k20 = (1/4)k((/k)0 )2 = g02 /4


Remarque : cas du clapotis, cest une onde stationnaire il faut diviser par 2. Il y a donc equipartition de lenergie.
1
E = Ec + Ep = g02 .
2
Pour une amplitude de 1m, et une longueur donde de 100m (periode de 8s), le flux denergie par metre lineque est de 3.2 104 W/m, ce qui donne pour environ
1000km de cotes (de Brest a la Bidassoa) soit 32 GW. La puissance totale est bien superieure.
Maintenant que nous avons defini lenergie, on veut savoir a quelle vitesse elle se deplace. On va voir que la vitesse de propagation de la houle est la moitie de la
vitesse de propagation des cretes.

6.2 Propagation de lEnergie


6.2.1 Profondeur faible

Nous avons deja vu plus haut, dans le cas kH  1 que lenergie et les cretes se propagent a la vitesse c0 = gH.
Houle et Vagues -28-

6.2.2 Profondeur infinie


houle de Stokes-Airy Observons se qui se passe maintenant dans ce cas dispersif, les equations devolution, apres multiplication par u pour la premiere et par v
pour la seconde et par integration par parties deviennent :
u P u u v P v v
u = +P & v = +P
t x x t y y

On additionne les deux et on ajoute t 2 /2 de chaque cote en remarquant que la condition de surface libre peut sy retrouver t
2 /2 =
t = v(x, 0)P (x, 0)
La somme des trois equations membre a membre fait se devoiler lenergie totale ainsi que la divergence de la vitesse qui est nulle :

u2 + v 2 + g 2 P u P v
( )= + v(x, 0)P,
t 2 x y
puis par integration en y de la surface au fond, il ne reste que

E =
t x
R0 2 2 2 R0
ou E = u( u +v2+g )dy est lenergie totale et = P udy est le flux dimpulsion. Il ne reste plus qua integrer en y (facteur 1/2k) puis a moyenner sur la
longueur donde (facteur 1/2) pour obtenir
g 2
E = 1/2g02 , puis de meme pour P u : < >= g02 /(2k)/2 = 2 0 ( 2k
) = ( 2k )E, mais comme 2 = gk, posons vg = /k = ( 2k
). On constate que :


E + vg E = 0
t x
Lenergie de tout une tranche moyennee en espace de la houle dispersive en eau infiniment profonde est transportee a la vitesse vg .

6.2.3 Profondeur quelconque


On peut faire le meme calcul dans le cas dune profondeur quelconque, mais cest plus penible ; le resultat est le meme : lenergie est transportee a la vitesse vg .

6.3 vitesse de groupe et vitesse de phase


Lelevation de la surface de leau est de la forme :
= A(x, t)exp(i(kx t)
ou = (kx t) est la phase, A est lamplitude que nous avons jusqua present pris constante, mais qui peut varier lentement avec t et x. On a donc pour :

k = /x, et = /t.

On remarque que lon a donc aussi : k/t + /x = 0 On peut faire apparatre une equation de continuite de la phase, lequation precedente peut aussi secrire
en effet :
k (kc)
+ =0
t x
(kc) k
Dans un milieu dispersif, (k) = kc(k), on a donc : x = k x

k k
+ vg = 0 avec vg =
t x k
Ce qui signifie que k est constant le long des courbes x, t verifiant dx/dt = vg , a k fixe, ce sont des droites.
Houle et Vagues -29-

Figure 26 Gauche Battement de deux ondes presentant un leger decallage. Droite paquet donde dans une porte de frequence donnant un signal module par sinus cardinal

6.4 Equations de pente douce


application equation Mild-slope equation ou equation de Berkhoff (1972)
(cvg ) + k 2 cvg = 0
un peu compliquee a etablir, mais si kH  1 elle secrit (en eliminant la vitesse cf debut du chapitre) :
2
(H(x)) + =0
g
ces equations sont comme loptique geometrique.

6.5 Modulation du paquet donde...


deux sinus purs On presente deux premieres approches simples... la premiere consiste a faire la somme de deux ondes sinusodale de frequence k0 et k0 + k0 /2 et
dobserver le battement... On fait la somme et on developpe grace aux identites trigonometriques :
cosp + cosq = 2cos(p + q)/2cos(p q)/2...
et comme on peut simplifier en developpant puisque lon suppose que lecart de frequences est petit :
k0 + k0 /2 k0 et (k0 + k0 ) (k0 ) + k0 /k + ...
donc la somme de deux ondes de frequences legerement differentes donne :
cos(k0 x (k0 )t) + cos((k0 + k0 )x (k0 + k0 )t) 2cos(k0 x (k0 )t)cos(k0 (x t/k)/2).
Lamplitude de londe somme 2cos(k0 x (k0 )t) est modulee parcos(k0 (x t/k)/2) lenveloppe se deplace a la vitesse /k. Cest ce que lon observe sur le
dessin suivant :
La porte :
La seconde etape classique lorsque lon introduit la notion de vitesse de groupe consiste a faire la somme dun paquet dondes de frequences piquees : une courbe en
cloche autour de k0 de largeur k0 . Une onde est consideree comme la somme de toutes ses frequences spatiales et alors on peut ecrire pour levolution dune onde
dont la forme au temps t = 0 a ete decomposee en modes F (k) : Z
= F (k)ei(kxt) dk

Houle et Vagues -30-

On privilegie pour commencer un cas simple : F (k) presente un pic en k0 et comme on peut developper :
k = k0 + (k k0 ) et = 0 + /k(k k0 ) + (1/2)000 (k k0 )/2 + ...
Posons = (k k0 ) et vg = /k et developpons le deplacement au voisinage de cette frequence :
Z
i(kxt)
=e F ()ei(xvg t) d

exemple : si F est reduite a la fonction porte F (k) = 1 pour k0 k0 /2 < k < k0 + k0 /2, et 0 sinon :
Z
F ()ei(xvg t) d = k0 sinc((x vg t)k0 /2)

la fonction sinus cardinal sinc(x) = sin(x)/x est une fonction piquee en 0. Lenveloppe se deplace a la vitesse de groupe vg , les petites ondes se deplacent plus vite
a la vitesse de phase c0 .

6.6 Plus complique : paquet au loin par la methode de la phase stationnaire


Le calcul precedent est un peu simpliste, reprenons le. On part de :
Z
= F (k)ei(kxt) dk

Cette integrale est en fait difficile a calculer car on la vu il y a beaucoup doscillations qui se compensent mutuellement. On va appliquer lidee de la phase stationnaire
(methode du col), pour cela on commence par privilegier un rayon x = V t, ce qui permet deliminer la dependance en x. Londe se developpe alors en :
Z
F (k)exp(i)dk ou la phase = V k kc(k)

La contribution principale de lintegrale est donc lorsque la phase varie peu avec k, cest a dire lorsque /k = 0. La derivee sannule justement pour un certain
k0 tel que :
V = vg = /k
On developpe au voisinage de ce k0 , on note 000 = 2 /k 2 :
1
= 0 + vg (k k0 ) + 000 (k k0 )2 + ... :
2
et on injecte dans lintegrant, Z p
= exp(i(kvg 0 )t)F (k)exp(000 ((k k0 )/( (2i))2 + ...))dk

(astuce i = 1/i). Or, seule la frequence k0 est selectionnee, en effet F (k) = F (k0 ) + (k k0 )F/k + ..., les termes autres que F (k0 ) ont une contribution negligeable,
il vient : Z p 2
exp(i(k0 vg 0 )t)F (k0 ) exp(000 ((k k0 )/( (2i)) ))dk

00
p 2
RLintegrale 0 ((k k0 )/(
de exp( (2i)) )))dk est reecrite comme lintegrale de (2i/(000 t))1/2 exp(s2 )ds par changement de variable en definitive puisque
2

exp(s )ds =

i 2ei/4 (000 t)1/2 F (k0 )exp(i(k0 x 0 t)), autour du rayon principal x/t = vg
1/2
Lamplitude de la surface surface decroit au loin en t , et ce le long du rayon vg , en dehors de ce rayon, les ondes sont inexistantes. Le paquet donde se deplace
bien a la vitesse vg , lamplitude decrot en 1/ t
Houle et Vagues -31-

6.7 Plus complique : paquet au loin par la methode de la phase stationnaire comparaison avec la fonction Airy

Malheureusement, si on part directement de = k gh0 (1 16 k 2 h20 + ...) il ny a pas de terme 000 a premiere vue. En fait, il y en a un si on suppose k petit et fini,
q
h2 c h2 c
on a alors vg = c0 02 0 k 2 donc 000 = h20 c0 k. dou k = (c0 t x)/( 02 0 t) en mettant dans lequation de la phase stationnaire ces deux racines on fait apparatre
le sinus et on obtient apres calcul lourd la forme de londe

2 (c0 t x)3/2
 

2 F (0) ((h20 c0 t/2)(c0 t x))1/4 sin +
3 (h20 c0 t/2)1/2 4

Sinon, si on ne veut pas faire ce calcul, on peut partir directement et en effectuant le developpement au voisinage de k = 0

ic0 H 3 k 3 t
Z
= F (0) exp(ik(x c0 t) + + ...)dk
6
en fait on retombe sur la definition integrale de lequation dAiry

1 sz 3
Z
Ai(z) = cos(sz + )ds
0 3

ce qui redonnerait par changement de variable


21/3 1/3 (x c0 t)
 
Ai 2
(h20 c0 t)
1/3 (h20 c0 t)1/3
cest la fameuse solution dAiry que lon a deja vue avec ds solutions semblables. Or la fonction dAiry a un comportement asymptotique
1 2 1 2 3/2
Ai(z) z 1/4 sin( |z|3/2 + /4) pour z < 0 ou pour z > 0 on a Ai(z) z 1/4 exp 3 z
3 2

et on fait la substitution
2 (c0 t x)3/2
 

2 F (0) ((h20 c0 t/2)(c0 t x))1/4 sin 2 1/2
+
3 (h0 c0 t/2) 4
on retrouve bien cette meme forme.
Au final, on voit quelle est la forme de la perturbation de surface libre : elle devient exponentiellement petite pour x > c0 t, elle est maximale en x = c0 t
(lapproximation avec (c0 t x)1/4 y diverge, il faut garder Airy), ensuite pour x < c0 t elle a un caractere ondulatoire, lamplitude decroit lentement en (c0 t x)1/4
au fur et a mesure que lon seloigne du front.
Il faut alors retourner aux comparaisons issues de Noda JGR 76 (30) 1971, Water waves generated by a local surface disturbance. Il compare, des experiences,
aux deux approximations proposees, Airy complet, la solution en phase stationnaire asymptotique donc avec la puissance -1/4 qui fait que londe diverge en x = c0 t
Houle et Vagues -32-

7 Viscosite
7.1 Amortissement visqueux en volume
Posons = 0 cos(kx t)eky = 0 sin(kx t) 0 = (k/)0 et definissons des moyennes en longueur
0
1 1 1
Z Z
< sin2 kx >=< cos2 (kx) >= sin2 (2x/)dx = et en prfondeur < e2ky >= e2ky dy =
0 2 (2k)

Lenergie perdue par la houle est due au frottement visqueux. La perte denergie est liee a la dissipation
dE
= < 2D : D >
dt
ou D est le tenseur des taux de deformation des vitesses, or lamplitude de la vague est 0 , et lamplitude la la vitesse u0 ils vont lentement varier en temps avec une
constante de temps bien superieure a la periode. De maniere tres simple, on peut estimer lamortissement des ondes gravitationnelles comme Landau (Landau p130
25), variation de lenergie pour tout un volume de controle :

de/dt = : D = pdV + 2D : D

Il nous faut evaluer < 2D : D >, Or par definition : Dij = (1/2)(i uj + j ui ) donc la dissipation :

D : D = (x v + y u)2 + 2(x u)2 + 2(y v)2 )

Or, = 0 cos(kx t)eky (surface libre en 0, profonduer en ) dans le cas de la houle de Airy-Stokes, 0 = (k/)0 lamplitude varie lentement sur une
periode... On calucle
y v + x u = 2k 2 0 sin(kx t)eky
donc (si on moyenne sur une periode pour x), et comme 0 varie lentement avec le temps on integre sur toute la profondeur... Donc la dissipation

< 2D : D >= [(4)(1/2) + 4(1/2)] < 20 > k 4 /(2k) = 220 k 3

varie comme le cube de k (cf Lamb 348) or cette dissipation correspond a la diminution denergie : dE/dt =< 2D : D > Lenergie mecanique totale E est la somme
de lenergie cinetique + potentielle, cest aussi : (2 *lenergie cinetique) = 2(< (u2 + v 2 ) > /2).
Donc la valeur moyenne de lenergie E est 2 < (u2 ) >= 2k 2 20 (1/2)1/(2k) sa variation par rapport au temps est donc telle que

d 2 d
(0 )(k/(2)) = 220 k 3 donc (20 ) = 4k 2 20
dt dt
dou une decroissance avec un temps caracteristique de 1/(4k 2 ) = g 2 /(4 4 ) or lamplitude de la vague est 0 et 0 = k0 , donc

d
(0 ) = 2k 2 0
dt
une decroissance avec un temps caracteristique de 1/(2k 2 ) = g 2 /(2 4 ) = 0.712 ( en cm) onde capillaire duree de vie courte (1cm dure 1seconde et parcourt 30
cm) Si =100m houle oceanique, duree de vie 900 jours et fait 30 fois le tour de la terre),
Houle et Vagues -33-

7.2 Amortissement visqueux au fond


Au fond, la vitesse de glissement est
kg0
ue = cos(kx t) = 0 /sinh(kH)cos(kx t)
coshkH
p
On peut dans la couche limite dans laquelle la viscosite va jouer, elle est depaisseur = /, trouver la solution :

t u = t ue + y2 u

ce qui donne  
y y
u = (0 /sinh(kH)) cos(kx t) exp( )cos(kx t +
2 2

pour une valeur de 1%, ca donne une epaisseur de couche limite de 0.28 T en cm (periode en secondes). Une houle de 10s de periode a une epaisseur de couche
limite laminaire de 1cm.
Il est ensuite interessant detudier le second ordre de cette equation. On constate que les termes en cos(kx t)2 sont present, ils produisent donc en moyenne un
leger courant de derive. La valeur de ce courant est fondamentale pour la connaissance du transport des particules en suspension.
Lenergie perdue par la houle est due au frottement visqueux. On peut montrer que si on calcule directement la dissipation :
dE
= < 2D : D >
dt
ou D est le tenseur des taux de deformation des vitesses, or lamplitude de la vague est 0 , et lamplitude la la vitesse u0 ils vont lentement varier en temps avec une
constante de temps bien superieure a la periode.
Z
< 2D : D >= ((x v + y u)2 + 2(x u)2 + 2(y v)2 )dV

l e terme dominant est


H
Z Z
(y u)2 dV = u20 car e2 ydy = 1/(2H)
2 2 0
dou la decroissance lente de lamplitude :
du20
r

= 2u20 , avec =
dt 8H 2
le temps associe f2 = H 2 /() le temps associe au frottement en volume est B = 1/(k 2 ), donc est tel que f2 = B2 (k 2 H 2 )(k 2 )/() et comme (k 2 )/()  1
(car B  1/) en prenant kH dordre 1, B  F  1/
Le frottement a la paroi est le plus fort compare au frottement dans le volume.
Houle et Vagues -34-

Gauche Tintin les Cigares du Pharaon, Herge, droits reserves, droite linspiration la grande vague de KanagawaKanagawa-oki nami-ura) de Katsushika Hokusai, noter
que linversion effectuee par Herge rend plus dramatique la grande vague que nous voyons a lenvers (sauf les specialistes du Manga) http://commons.wikimedia.org

Asterix, La Grande Traversee, Goscinny, droits reserves...


Houle et Vagues -35-

7.3 Generation de la Houle


Comment nait la houle ? Elle est generee par le vent. Un mecanisme communement invoque est celui de Kelvin-Helmholtz, nous verrons en cours de stabilite cette
theorie et ses limites.

7.4 Retour sur les vagues pres des cotes


Le deferlement

En pratique on definit un parametre appele parametre de similarite de surf :


p
= tan /H

avec la pente de la plage, de petites valeurs de la pente sont associees a un parametre petit et corresponde a une plage dissipative. Les vagues vont arriver comme
sur la figure 3.4 sur une pente douce, grandir en taille et se casser par deferlement. On parle de spilling breaking wave
Si ce parametre est tres grand, les vagues sont reflechies sur la cote, la plage est reflective. Sil est un peu moins grand, elle vont deferler. Nous verrons ce probleme
du deferlement au chapitre suivant avec un autre point de vue.

la langue anglaise distingue 4 types de vagues, spilling breaking (deferlement glissant ou progressif), plunging breaking (deferlement plongeant tres intense),
collapsing breaking (sans traduction en France) et surging breaking (deferlement frontal ou gonflant).
voir figure.
Houle et Vagues -36-

!!"#"$

(/0*123#$.#4$#5678#'%%4
!"#$$$%&'&$$%%#()*+,#--.
*.##C;=77=12#>+3*?=12#@*A3 >.##)7612=12#>+3*?=12#@*A3
%&'()*+,-)./0'+0/,12345

B.##C6+2=12#>+3*?=12#@*A3 9.##/:77*;<=12#>+3*?=12#@*A3

D=26+3#--&E&$F###G+3*?3+#,8;3<
Figure 27 dapres Richardson FloSci-TN44.pdf, vagues deferlant pres sur une pente douce
Houle et Vagues -37-

8 Conclusion
Au travers de ces pages nous avons presente la base mecanique des ecoulements de vagues linearisees en profondeur quelconque. Le probleme est resolu completement
pour une onde monochromatique. Retenons que le signal est dispersif, un signal quelconque va se deformer et se decomposer en un train dondes, les longueurs donde
grandes se deplacent plus vite que les courtes. Lenergie se deplace a la vitesse de groupe differente de la vitesse de phase lorsque le milieu est dispersif. Lorsque la
cambrure augmente, de nouveaux termes vont apparatre et complexifier la resolution (cest ce qua fait Stokes). Pres du rivage, une legere augmentation de lamplitude
apparat a cause de la remontee du fond.

la plage de Treach er Goured, Ile dHouat est convexe a son extremite (a droite, la photo (PYL) de gauche a ete prise sur la point de la fleche, carte GoogleMaps),
des vagues paralleles aux deux pans de plage se croisent !
Houle et Vagues -38-

References
[1] J. Billingham,A. C. King Wave motion sur googlebooks
[2] Guyon E., Hulin J.-P. & Petit L. (1991) Hydrodynamique Physique, InterEditions, Ed. CNRS 506p, (voir page 208 et suivantes)
EDP Science 2001 pour la deuxieme edition
[3] Guyon E., Hulin J.-P., Petit L. & Mitescu C.D. (2001) : Physical hydrodynamics Oxford University Press (version anglaise)
[4] Cours des Glenans.
[5] Paterson A.R. (1983) : A first course in fluid dynamics, Cambrige, 528p (voir chapitre XV) googlebooks
[6] A. Petitjean (1998) Recent Progresses in Dam-Break Modelling in France, http://kfki.baw.de/conferences/ICHE/1998-Cottbus/274.pdf Case Studies : Dam Break,
Dam Break Modelling
[7] O. Thual http://thual.perso.enseeiht.fr/xsee/2010/00main.pdf
[8] Mei C.C. (1989) the applied dynamics of ocean surface waves adv series on Ocean Eng. Vol I, World Scientific
http://fr.wikipedia.org/wiki/Vague
http://scs-ingenierie.pagesperso-orange.fr/fr/acceuil.html
http://apere.free.fr/vagues/propag.htm qui prennent mes images de la science en f^
ete sans me citer.

La vague dans la peinture et dans les arts, attention contenu sexuel explicite :
http://fr.wikipedia.org/wiki/Fichier:Guillaume Seignac La Vague.jpg/
la vague de Courbet
la vague de Hokusai
Nee de la Vague de Lucien Clergue

un film classique sur les vagues rtsp://18.39.0.144:554/ifluids/Waves in Fluids.rm?cloakport=8080,554 en Real Player


http://web.mit.edu/hml/ncfmf/14WF.pdf, le synopsis,
existe aussi sur youtube http://www.youtube.com/watch?v=MNyebpog i0.

La page ouebe de ce texte


http://www.lmm.jussieu.fr/lagree/COURS/MFEnv/

La derniere version (9 fevrier 2016)) :


http://www.lmm.jussieu.fr/lagree/COURS/MFEnv/MFEhoule.pdf

Raymond Subes Sans Titre 1961 (entree de Jussieu Quai Saint Bernard)

Vous aimerez peut-être aussi