Chapitre 

Exercices complémentaires

O. Thual,  octobre 

Sommaire
 Exercice sur le chapitre  . . . . . . . . . . . . . . . . 
Milieu très profond . . . . . . . . . . . . . . . . . . . 
Cisaillement de deux couches infinies . . . . . . . . . 
Tension superficielle et vagues . . . . . . . . . . . . . 
Stabilité de l’écoulement de Poiseuille . . . . . . . . . 
 Exercice sur le chapitre  . . . . . . . . . . . . . . . . 
Rouleaux et relation de dispersion . . . . . . . . . . . 
Milieu poreux en rotation . . . . . . . . . . . . . . . 
Convection thermohaline . . . . . . . . . . . . . . . . 
 Exercice sur le chapitre  . . . . . . . . . . . . . . . . 
Bifurcation fourche imparfaite . . . . . . . . . . . . . 
Bifurcation avec hysteresis . . . . . . . . . . . . . . . 
Système dynamique et polynôme de degré trois . . . 
Équations de Saint-Venant . . . . . . . . . . . . . . . 


 Chapitre . Exercices complémentaires

 Exercice sur le chapitre 

EXERCICE . Milieu très profond

On considère une couche fluide infiniment profonde dont la surface libre, d’équa-
tion z = η(x, t), est en contact avec l’atmosphère de pression constante pa . On
suppose que le mouvement du fluide est décrit par les équations d’Euler incom-
pressibles D
∂U 1
div U = 0, + U · grad U = − grad p − g ez , (.)
∂t ρ
où ρ est la masse volumique, U (x, z, t) = u ex + w ez est le champ de vitesse et
p(x, z, t) est le champ de pression.

z pa
η
x
U

Figure . – Couche infiniment profonde à surface libre d’équation z = η(x, t).
∂η ∂η
) Interpréter les conditions aux limites ∂t +u ∂x = w et p = pa en z = η(x, t)
et limz→−∞ w = 0.
La condition aux limites cinématique de surface énonce que la vitesse normale à la
surface libre, d’équation F (x, z, t) = z − η(x, t) = 0, est égale à la vitesse normale
de cette surface. Cette condition s’écrit dF dt = ∂t + U · grad F = 0 ou encore
∂F
∂η ∂η
∂t + u ∂x = w. La condition aux limites dynamique exprime l’égalité des forces de
contacts qui sont des forces de pression. La condition cinématique w = 0 sur un fond
plat est ici reportée à l’infini z → −∞.
) Indiquer les étapes qui permettent de modéliser les petites oscillations
irrotationnelles de la surface libre à l’aide du système d’équations U =
grad φ, ∆φ = 0 et p = pa − ρ g z − ρ ∂φ ∂t avec les conditions aux limites
∂η ∂φ ∂φ ∂φ
∂t = ∂z et ∂t = −g η en z = 0 et limz→−∞ ∂z = 0.
L’hypothèse irrotationnelle permet d’écrire le champ de vitesse sous la forme U =
grad φ où φ(x, z, t) est le potentiel des vitesses. La condition div U = 0 s’écrit alors
∆φ = 0. L’équation de quantité de mouvement linéarisée ∂U ∂t = − ρ grad p − g ez
1
∂φ
s’intègre en ∂t +p/ρ+g z = pa /ρ, la constante d’intégration pa /ρ pouvant être choisie
ainsi dans la mesure où φ est défini à une constante C(t) près. En développant les
valeurs f [x, η(x, t), t] = f (x, 0, t)[1 + O(η)] des champs en surface pour les petites
oscillations, on peut appliquer les conditions aux limites en z = 0 plutôt qu’en
z = η(x, t). On linéarise alors ces conditions aux limites en utilisant w = ∂φ ∂z et
∂φ
p = −ρ ∂t + pa − ρ g z.
) Justifier la recherche de solutions de la forme φ = Φ(z) ei kx x−i ω t et
η = ηm ei kx x−i ω t où Φ(z) et ηm sont complexes. Montrer que Φ(z) =
Φm exp(k z) où k = |kx |.
Les équations étant linéaires à coefficients réels, la partie réelle d’une solution com-
plexe est aussi solution. Comme les coefficients sont constants, les solutions sont
des exponentielles. On vérifie a posteriori que le taux de croissance 00
temporel de ces
ondes est nul. En reportant dans l’équation ∆φ = 0, on obtient Φ (z) + k 2 Φ(z) = 0.
En utilisant la condition aux limites en −∞, on voit que Φ(z) = Φm exp(k z) où
Φm est une amplitude complexe arbitraire.
√
) En déduire la relation de dispersion ω = ± g k des ondes de surfaces
en eaux profondes. Montrer que l’on peut choisir ηm réel sans perte de
Exercice sur le chapitre  

généralité et exprimer la solution réelle sous la forme η = ηm cos(kx x−ω t)

et φ = g ηωm ek z sin(kx x − ω t). Montrer que les trajectoires de ces petites
oscillations sont des cercles de centres (x0 , z0 ) et de rayon a(z0 ) que l’on
calculera. Quel est l’ordre de grandeur de la profondeur au-dessous de
laquelle une houle de longueur d’onde Lx devient négligeable ?
Les conditions aux limites s’écrivent −i ω ηm = k Φm et −i ω Φm = −g ηm . On en
déduit que ω 2 = g k et Φm = −i g ηm /ω. On peut se ramener à ηm réel en changeant
l’origine des x ou des t. La partie réelle de la solution complexe est alors de la forme
indiquée. Comme u = ∂φ ∂φ
∂x et w = ∂z , la trajectoire de centre (x0 , y0 ) est décrite
par x(x0 , z0 ; t) = x0 − a sin(kx x − ω t) et z(x0 , z0 ; t) = z0 + a cos(kx x − ω t) où
a(z0 ) = ηm exp(k z0 )/ω est le rayon du cercle. À la profondeur Lx , l’amplitude du
mouvement a diminué d’un facteur exp(−2 π) ∼ 2 10−3 .

EXERCICE . Cisaillement de deux couches infinies

On considère la superposition de deux couches fluides incompressibles de

masses volumiques respectives ρ1 et ρ2 avec ρ1 ≥ ρ2 . À l’équilibre, la couche la

U z
ρ2 U2 η(x, t)
x
ρ1 O
U1

Figure . – Cisaillement de deux couches infinies avec la vitesse U ex .

plus lourde est au repos et occupe le demi-espace z ≤ 0, tandis que la couche la

plus légère est animée d’une vitesse constante U ex dans le demi-espace z ≥ 0
(Figure .). On s’intéresse aux petites perturbations modélisées par les équa-
tions d’Euler D :
∂U 2 1
div U 2 = 0 et + U 2 · grad U 2 = − grad p2 − g ez ,
∂t ρ2
∂U 1 1
div U 1 = 0 et + U 1 · grad U 1 = − grad p1 − g ez , (.)
∂t ρ1
où x = x ex + z ez est le vecteur des coordonnées, U 1 (x, t) = u1 ex + w1 ez
et U 2 (x, t) = u2 ex + w2 ez sont les champs de vitesses respectifs des deux
couches et p1 (x, t) et p2 (x, t) leurs champs de pression. Le vecteur unitaire ez
est vertical et g est l’intensité de la gravité.
) D’où proviennent les conditions aux limites ∂η ∂η ∂η ∂η
∂t + u1 ∂x = w1 , ∂t + u2 ∂x =
w2 et p1 = p2 en z = η(x, t) ? Quelles sont les conditions aux limites en
z = ±∞ ?
) Indiquer brièvement les étapes qui permettent de modéliser les petites os-
cillations irrotationnelles de la surface libre à l’aide du système d’équa-
tions U 1 = grad φ1 , ∆φ1 = 0, U 2 = grad (U x + φ2 ), ∆φ2 = 0,
 Chapitre . Exercices complémentaires

p1 = pr1 − ρ1 g z − ρ1 ∂φ ∂φ2 ∂φ2

∂t , p2 = pr2 − ρ2 g z − ρ2 ( ∂t + U ∂x ). Comparer
1

les constantes d’intégration pr1 et pr2 entre elles en considérant l’équilibre.

) Justifier
  aux limites en z = 0 suivantes :
les conditions

ρ1 ∂t + g η = ρ2 ∂φ
∂φ1 ∂φ2 ∂η ∂φ1 ∂η ∂η
∂t + U ∂x + g η , ∂t = ∂z et ∂t + U ∂x = ∂z .
2 ∂φ2

) Justifier la recherche de solutions de la forme φ1 = Φ1 (z) ei kx x+s t , φ2 =

Φ2 (z) ei kx x+s t et η = ηm ei kx x+s t où Φ1 (z), Φ2 (z) et ηm sont complexes
ainsi que s = σ − i ω. Montrer que l’on peut écrire Φ1 (z) = Φ1m F1 (k z)
et Φ2 (z) = Φ2m F2 (k z) avec k = |kx | où F1 et F2 sont deux fonctions que
l’on explicitera.
) En déduire la relation de dispersion généralisée s’écrit ρ1 [γ + s2 ] = ρ2 [γ −
(s + i kx U )2 ] où γ est une constante que l’on déterminera.
) En déduire que les racines sont de la forme
 q 
2 2 2
s = β −i ρ2 kx U ± ρ1 ρ2 k U − (ρ1 − ρ2 ) g k 2

où β est une constante que l’on déterminera.

) Exprimer le taux de croissance σ pour une perturbation de longueur d’onde
Lx lorsqu’elle conduit à une instabilité.
) On suppose que ρ2 = 1040 kg/m3 est la salinité de la Mer Rouge dont
une couche de vitesse U = 1 m/s est mise en contact avec une couche
dont la salinité ρ1 = 1440 kg/m3 est celle de la Mer Morte. Estime r la
taille caractéristiques Lc des plus gros tourbillons que l’on voit se former
lorsque les deux couches fluides se rencontrent. Calculer leur vitesse de
propagation. On donne g = 10 m/s2 .
) On suppose que l’on observe, en plus des tourbillons de taille L ≤ Lc , des
structures d’extension horizontale plus grande que Lc . Quelle peut être la
nature de ces structures ? À quelles vitesses se propagent-elles ?

Corrigé . Cisaillement de deux couches infinies

)Les conditions aux limites cinématiques de surface énoncent que la vitesse

normale à la surface libre, d’équation F (x, z, t) = z − η(x, t) = 0, est égale à la
vitesse normale de cette surface. Cette condition s’écrit dF ∂F
dt = ∂t +U ·grad F =
∂η ∂η ∂η ∂η
0 ou encore ∂t + u1 ∂x = w1 pour le fluide du bas et ∂t + u2 ∂x = w2 pour le
fluide du haut. La condition aux limites dynamique exprime l’égalité des forces
de contacts qui sont des forces de pression. On a donc p1 = p2 . La condition
cinématique w = 0 sur un plan est ici reportée à l’infini. Elle s’écrit w1 = 0
pour z → −∞. et w2 = 0 pour z → ∞.
)L’hypothèse irrotationnelle permet d’écrire le champ de vitesse de la couche
du bas sous la forme U 1 = grad φ1 où φ1 (x, z, t) est le potentiel des vitesses.
La condition div U 1 = 0 s’écrit alors ∆φ1 = 0. Pour la couche du haut, on
peut choisir U 2 = grad (U x + φ2 ) de sorte que φ2 (x, z, t) soit le potentiel des
écarts de vitesse au champ constant U ex . L’équation de quantité de mouvement
linéarisée ∂U 1 ∂φ1
∂t = − ρ grad p − g ez s’intègre, pour la couche du bas, en ρ1 ∂t +
p1 + ρ1 g z = pr1 , la constante d’intégration pr1 pouvant être choisie de manière
arbitraire dans la mesure où φ1 est défini à une constante C1 (t) près. Pour la
∂φ1 ∂φ1
couche du haut, on obtient de même ρ2 ∂t + U ∂x + p2 + ρ2 g z = pr2 .
Exercice sur le chapitre  

À l’équilibre, lorsque φ1 = φ2 , la continuité du champ de pression impose

pr1 = pr2 .
)En développant les valeurs f [x, η(x, t), t] = f (x, 0, t)[1 + O(η)] des champs
en surface pour les petites oscillations, on peut appliquer les conditions aux
limites en z = 0 plutôt qu’en z = η(x, t). On linéarise alors les conditions aux
limites cinématique à l’interface en utilisant w1 = ∂φ ∂φ2
∂z et w2 = ∂z . On écrit
1

ensuite la condition aux limites dynamique p1 = p2 en z = 0 en ne négligeant

pas η devant φ1 et φ2 .
)Les équations étant linéaires à coefficients réels, la partie réelle d’une solution
complexe est aussi solution. Comme les coefficients sont constants, les solutions
sont des exponentielles. En reportant dans les équations ∆φ1 = 0 et ∆φ2 = 0,
00 00
on obtient Φ1 (z) − k 2 Φ1 (z) = 0 et Φ2 (z) − k 2 Φ2 (z) = 0. En utilisant les
conditions aux limites à l’infini, on voit que Φ1 (z) = Φ1m exp(k z) et Φ2 (z) =
Φ2m exp(−k z) où Φ1m et Φ2m sont des amplitudes complexes arbitraires. On
a donc F1 (Z) = eZ et F2 (Z) = e−Z .
)En utilisant Φ01 (0) = k Φ1m et Φ02 (0) = −k Φ2m , les conditions aux li-
mites cinématiques en z = 0 s’écrivent s ηm = k Φ1m et (s + i kx U ) ηm =
−k Φ2m . La condition aux limites dynamique s’écrit ρ1 [g k ηm + s Φ1m ] =
ρ2 [g k ηm + (s + i kx U )Φ1m ]. En éliminant Φ1m et Φ2m
 et en écartant la  solu-
2 2

tion triviale ηm = 0, on en déduit ρ1 g k + s = ρ2 g k − (s + i kx U ) . On
a donc γ = g k.
)Le discriminant réduit de l’équation (ρ1 +ρ2 ) s2 +2 i ρ2 kx U s+(ρ1 −ρ2 ) g k−
ρ2 k 2 U = 0 s’écrit ∆0 = −ρ2 k 2 U 2 − (ρ21 − ρ22 ) g k + ρ2 (ρ1 + ρ2 ) k 2 U =
ρ1 ρ2 k 2 U − (ρ21 − ρ22 ) g k. Les racines ont bien la forme indiquée avec β =
1/(ρ1 + ρ2 ).
p
)Le taux de croissance s’écrit σ = ρ1 ρ2 k 2 U 2 − (ρ21 − ρ22 ) g k/(ρ1 + ρ2 ) avec
k = 2 π/Lx .
)Seules les perturbations de longueur d’onde L ≤ Lc avec Lc = 2gπ ρρ21−ρ
ρ2 2
2 U ∼
1 2
0, 95 m sont instables. Les plus gros tourbillons, de taille voisine de Lc , se
propagent à la vitesse de phase cϕ = ρ1ρ+ρ
2
2
U ∼ 0, 4 m/s.
)Les structures
h de tailles L ≥ Lc sont des ondes
i de gravité de vitesse de phase
p
cϕ = β −ρ2 U ± −ρ1 ρ2 U 2 + (ρ21 − ρ22 ) g/k . Il faudrait également calculer
la vitesse de groupe pour rendre compte de l’observation de la propagation de
paquets d’ondes.

EXERCICE . Tension superficielle et vagues

Si R est le rayon de courbure de l’interface z = η(x, t) délimitant deux fluides

non miscibles, la tension superficielle induit une différence de pression p2 −p1 =
T /R où T est une constante, p2 = p|z=η+ et p1 = p|z=η− (figure .). Lorsque
η est petit, cette condition s’écrit
∂2η
p2 − p1 = T en z = η(x, t) . (.)
∂x2

On suppose ici que le deux fluides sont incompressible de masses volumiques

 Chapitre . Exercices complémentaires

z
ρ2 U2 R p η(x, t)
2

ρ1 O p1 x
U1

Figure . – Interface de deux couches infinies en prenant en compte la tension

superficielle.

respectives ρ1 et ρ2 avec ρ1 ≥ ρ2 . À l’équilibre, la couche la plus lourde est au

repos et occupe le demi-espace z ≤ 0. On s’intéresse aux petites perturbations
modélisées par les équations d’Euler D :
∂U 2 1
div U 2 = 0 et + U 2 · grad U 2 = − grad p2 − g ez ,
∂t ρ2
∂U 1 1
div U 1 = 0 et + U 1 · grad U 1 = − grad p1 − g ez , (.)
∂t ρ1
où x = x ex + z ez est le vecteur des coordonnées, U 1 (x, t) = u1 ex + w1 ez
et U 2 (x, t) = u2 ex + w2 ez sont les champs de vitesses respectifs des deux
couches et p1 (x, t) et p2 (x, t) leurs champs de pression. Le vecteur unitaire ez
est vertical et g est l’intensité de la gravité.
) D’où proviennent les conditions aux limites ∂η ∂η ∂η ∂η
∂t + u1 ∂x = w1 , ∂t + u2 ∂x =
2
∂ η
w2 et p1 = p2 − T ∂x 2 en z = η(x, t) ? Justifier les conditions aux limites
w1 = 0 en z = −∞ et w2 = 0 en z = ∞.
) Indiquer brièvement les étapes qui permettent de modéliser les petites os-
cillations irrotationnelles de la surface libre à l’aide du système d’équations
U 1 = grad φ1 , ∆φ1 = 0, U 2 = grad φ2 , ∆φ2 = 0, p1 = pr1 − ρ1 g z − ρ1 ∂φ 1
∂t ,
p2 = pr2 − ρ2 g z − ρ2 ∂φ∂t . Comparer les constantes d’intégration pr1 et pr2
2

entre elles.    
∂2η
) Justifier les conditions aux limites ρ1 ∂φ ∂t
1
+ g η = ρ 2
∂φ2
∂t + g η +T ∂x 2,
∂η
∂t = ∂φ ∂η ∂φ2
∂z et ∂t = ∂z en z = 0.
1

) Justifier la recherche de solutions de la forme φ1 = Φ1 (z) ei kx x+s t , φ2 =

Φ2 (z) ei kx x+s t et η = ηm ei kx x+s t où Φ1 (z), Φ2 (z) et ηm sont complexes.
Montrer que l’on peut écrire Φ1 (z) = Φ1m F1 (k z) et Φ2 (z) = Φ2m F2 (k z)
avec k = |kx | où F1 et F2 sont deux fonctions que l’on explicitera.
) En déduire la relation de dispersion s’écrit ρ1 (γ + s2 ) = ρ2 (γ − s2 ) − T k 3
où γ est une constante que l’on déterminera.
) En déduire que les racines sont de la forme s = ±i ω avec
r
ρ1 − ρ2
ω = ge(k) k
ρ1 + ρ2
où ge(k) est une fonction de k que l’on exprimera.
) Exprimer la relation de dispersion des ondes de surface dans le cas ρ1  ρ2 .
Exercice sur le chapitre  

) Lorsque la couche supérieure est animée d’une vitesse moyenne U2 , la

relation de dispersion dans le cas T = 0 (pas de tension superficielle) est

s2
 
= ρ2 g k − (s + i kx U2 )2 .
 
ρ1 gk+
tanh(k hr )

Comment cette relation se généralise-t-elle au cas T 6= 0 (prise en compte

de la tension superficielle).
) Dans le cas de l’interface air-eau, on mesure T = 7 10−2 N/m. En dé-
duire l’ordre de grandeur des longueurs d’onde à partir desquelles on peut
négliger la tension superficielle. On pourra prendre ρ1 = 1000 kg/m3 et
g = 10 m/s2 .

Corrigé . Tension superficielle et vagues

)Les conditions aux limites cinématiques de surface énoncent que la vitesse

normale à la surface libre, d’équation F (x, z, t) = z − η(x, t) = 0, est égale à la
vitesse normale de cette surface. Cette condition s’écrit dF ∂F
dt = ∂t +U ·grad F =
0 ou encore ∂η ∂η ∂η ∂η
∂t + u1 ∂x = w1 pour le fluide du bas et ∂t + u2 ∂x = w2 pour le
fluide du haut. La condition aux limites dynamique exprime le saut de pression
dû à la tension superficielle. La condition cinématique w = 0 sur un plan est ici
reportée à l’infini. Elle s’écrit w1 = 0 pour z → −∞. et w2 = 0 pour z → ∞.
)L’hypothèse irrotationnelle permet d’écrire le champ de vitesse de la couche
du bas sous la forme U 1 = grad φ1 où φ1 (x, z, t) est le potentiel des vitesses.
La condition div U 1 = 0 s’écrit alors ∆φ1 = 0. Pour la couche du haut, on a
de même U 2 = grad φ2 . L’équation de quantité de mouvement linéarisée ∂U ∂t =
− ρ1 grad p − g ez s’intègre, pour la couche du bas, en ρ1 ∂φ∂t
1
+ p 1 + ρ 1 g z = pr1 ,
la constante d’intégration pr1 pouvant être choisie de manière arbitraire dans
la mesure où φ1 est défini à une constante C1 (t) près. Pour la couche du haut,
on obtient de même ρ2 ∂φ ∂t + p2 + ρ2 g z = pr2 . À l’équilibre, lorsque φ1 = φ2 ,
1

le champ de pression est continu (car η = 0) et l’on a donc pr1 = pr2 .

)En développant les valeurs f [x, η(x, t), t] = f (x, 0, t)[1 + O(η)] des champs
en surface pour les petites oscillations, on peut appliquer les conditions aux
limites en z = 0 plutôt qu’en z = η(x, t). On linéarise alors les conditions aux
limites cinématiques à l’interface en utilisant w1 = ∂φ ∂φ2
∂z et w2 = ∂z . On écrit
1

ensuite la condition aux limites dynamique p1 = p2 en z = 0 en ne négligeant

pas η devant φ1 et φ2 .
)Les équations étant linéaires à coefficients réels, la partie réelle d’une solution
complexe est aussi solution. Comme les coefficients sont constants, les solutions
sont des exponentielles. En reportant dans les équations ∆φ1 = 0 et ∆φ2 = 0,
00 00
on obtient Φ1 (z) − k 2 Φ1 (z) = 0 et Φ2 (z) − k 2 Φ2 (z) = 0. En utilisant les
conditions aux limites à l’infini, on voit que Φ1 (z) = Φ1m exp(k z) et Φ2 (z) =
Φ2m exp(−k z) où Φ1m et Φ2m sont des amplitudes complexes arbitraires. On
a donc F1 (Z) = eZ et F2 (Z) = e−Z .
)En utilisant Φ01 (0) = k Φ1m et Φ02 (0) = −k Φ2m , les conditions aux limites
cinématiques en z = 0 s’écrivent s ηm = k Φ1m et s ηm = −k Φ2m . La condi-
tion aux limites dynamique s’écrit ρ1 [g k ηm + s Φ1m ] = ρ2 [g k ηm + s Φ1m ] −
 Chapitre . Exercices complémentaires

T kx2 ηm . En éliminant Φ1m

et Φ2m et en
écartant la solution triviale ηm = 0,
on en déduit ρ1 g k + s2 = ρ2 g k − s2 − T k 3 . On a donc γ = g k.
T
)On en déduit la relation de dispersion exprimée avec ge(k) = g + ρ1 −ρ2 k2 .

)Dans le ρ1  ρ2 , on a ω = ge(k) k avec ge(k) = g + ρT1 k 2 .

p

)En reprenant les étapes du calcul, on voit qu’il suffit de remplacer g par
T
ge(k) = g + ρ1 −ρ2
k 2 dans la relation de dispersion.
T
)On peut négliger la tension superficielle lorsque g  k 2 , c’est-à-dire pour
ρ1q
des longueurs d’ondes l = 2 π/k telles que l  lc avec lc = 2 π gTρ1 ∼ 1.7 cm.

EXERCICE . Stabilité de l’écoulement de Poiseuille

On considère un fluide pesant de masse volummique ρ compris entre deux

plaques planes horizontales délimitant les frontières z = −h et z = h. On note
g = −g ez le champ de gravité. On suppose que le mouvement de ce fluide est
régi par le modèle des équations de Navier-Stokes incompressibles : div U = 0,
dU /dt = −grad Π − g ez + ν ∆U , où g est la gravité, p = ρ Π est le champ
de pression et ν la viscosité cinématique. On suppose que le mouvement est
D avec U (x, z, t) = u ex + w ez . Les conditions aux limites sur les deux parois
s’écrivent u(x, h, t) = w(x, h, t) = 0 et u(x, −h, t) = w(x, −h, t) = 0.
) On considère le profil de vitesse U e = U (z) ex avec U (z) = U0 [1 − (z/h)2 ].
Caculer la dérivée particulaire dU e /dt et tracer le profil U (z).
) Calculer le champ Πe (x, z) tel que (U e , Πe ) soit un équilibre du modèle en
supposant que Πe (0, 0) = U02 .
) On définit les variables adimensionnées par x = h x∗ , U = U0 U ∗ ,
t = (h/U0 ) t∗ et Π = U02 Π∗ . Écrire les équations de Navier-Stokes adi-
mensionnées en faisant apparaitre le √ nombre de Reynolds Re = U0 h/ν
et le nombre de Froude F r = U0 / g h. On pourra omettre les ∗ pour
présenter le résultat final.
) Écrire les conditions aux limites ainsi que U e = U (z) ex et Πe (x, z) après
adimensionnalisation, les ∗ étant omises désormais dans toutes les expres-
sions.
) Pour ces équations adimensionnées, on note u = U (z) + u e(x, z, t), w =
w(x,
e z, t) et Π = Πe (x, z) + Π(x,
e z, t) où u
e, we et Πe sont des petites pertur-
bations. Écrire les équations linéarisées du modèle en ne conservant que
les termes dominants.
) Justifier l’existence d’une fonction ψ(x, z, t) telle que ue = − ∂ψ e = ∂ψ
∂z et z ∂x .
Comment se nomme-t-elle ?
) Montrer que l’élimination de la pression conduit à l’équation d’Orr-Som-
merfeld
 
∂ ∂ ∂ψ 00
+ a1 U (z) ∆ψ + a2 U (z) = a3 Re−1 ∆2 ψ
∂t ∂x ∂x

où a1 , a2 et a3 sont des constantes dont on indiquera la valeur.

) Expliciter les conditions aux limites à l’aide de ψ.
Exercice sur le chapitre  

0, 02
<latexit sha1_base64="vRZYayh4gC2zefUbrPaA2hUHH+4=">AAAC7XicjVLLSsNAFD2Nr1pfVZdugkVwIWXSjV0W3LisYB9QiyTptB2aJiGZCKX4C+5EF+LWb/IP1K/wznQKPvAxIcmZc+85M/fOeHEgUsnYc85aWFxaXsmvFtbWNza3its7zTTKEp83/CiIkrbnpjwQIW9IIQPejhPujr2At7zRiYq3rniSiig8l5OYd8fuIBR94btSUeyIVS6LJVZmetjfgWNACWbUo+IbLtBDBB8ZxuAIIQkHcJHS04EDhpi4LqbEJYSEjnNco0DajLI4ZbjEjug7oFnHsCHNlWeq1T6tEtCbkNLGAWkiyksIq9VsHc+0s2J/8p5qT7W3Cf094zUmVmJI7F+6eeZ/daoWiT6qugZBNcWaUdX5xiXTXVE7tz9UJckhJk7hHsUTwr5Wzvtsa02qa1e9dXX8RWcqVs19k5vh9dfq1A5mvaZzoWvgfD3076BZKTus7JyxUq1qLkQee9jHIZ36MWo4RR0NWmWIG9zh3oqsW+vBepylWjmj2cWnYT29AwIHmD0=</latexit>
sha1_base64="WSGCsVgenbSkP9TR+v7TCK2HFas=">AAAC7XicjVK7TsMwFD0Nr1JehY4sEVUlBlQlLHSsxMJYJPqQSoUS122jpnGUOEhVxcbMhmBArHwTGyPwFVy7qQRUPBwlOT73nmPfa7uh78XSsp4zxsLi0vJKdjW3tr6xuZXf3mnEIokYrzPhi6jlOjH3vYDXpSd93goj7oxcnzfd4bGKNy95FHsiOJPjkHdGTj/weh5zpKKsA+vwIl+0ypYe5jywU1CsFkrXLwBqIv+Oc3QhwJBgBI4AkrAPBzE9bdiwEBLXwYS4iJCn4xxXyJE2oSxOGQ6xQ/r2adZO2YDmyjPWakar+PRGpDRRIo2gvIiwWs3U8UQ7K/Yn74n2VHsb099NvUbESgyI/Us3y/yvTtUi0UNF1+BRTaFmVHUsdUl0V9TOzU9VSXIIiVO4S/GIMNPKWZ9NrYl17aq3jo6/6kzFqjlLcxO8/Vqd2sG013QudA3s74c+DxqHZdsq26d0HyqYjix2sYd9OvUjVHGCGuq0ygA3uMO9IYxb48F4nKYamVRTwJdhPH0AY6KaGQ==</latexit>
sha1_base64="mubtl/C6fCldwesg4XDt2dmlJDQ=">AAAC7XicjVLLSsNAFD2Nr1pf1S7dBEvBhZSkG7ssuHFZwT6gFknSaTs0TUIyEUpx59qd6ELc+if+g36AqGs/wDvTFNTiY0KSM+fec2bunbEDl0fCMB5T2tz8wuJSejmzsrq2vpHd3KpHfhw6rOb4rh82bStiLvdYTXDhsmYQMmtou6xhDw5kvHHGwoj73rEYBaw9tHoe73LHEpIy9ozSaTZvFA019FlgJiBfyRUu3h+en6p+9g0n6MCHgxhDMHgQhF1YiOhpwYSBgLg2xsSFhLiKM5wjQ9qYshhlWMQO6NujWSthPZpLz0ipHVrFpTckpY4CaXzKCwnL1XQVj5WzZH/yHitPubcR/e3Ea0isQJ/Yv3TTzP/qZC0CXZRVDZxqChQjq3MSl1h1Re5c/1SVIIeAOIk7FA8JO0o57bOuNJGqXfbWUvEXlSlZOXeS3Bivv1YndzDpNZ0LXQPz+6HPgnqpaBpF84juQxmTkcY2drBLp76PCg5RRY1W6eMS17jRfO1Ku9XuJqlaKtHk8GVo9x99AJyz</latexit>

Re = 9772
<latexit sha1_base64="TIGRBu4bfV3pmIzzM0qS3Tm7svg=">AAAC8HicjVLLSsNAFD2Nr1pfVZdugkVwVRI3rQuh4MZlFfuAWiRJp3VoXkwmQil+hDvRhbj1j/wD9Su8M01BLT4mJDlz7j1n5t4ZN/Z5Ii3rJWfMzS8sLuWXCyura+sbxc2tZhKlwmMNL/Ij0XadhPk8ZA3Jpc/asWBO4Pqs5Q6PVbx1zUTCo/BcjmLWDZxByPvccyRRrTN2dFipHFwWS1bZ0sOcBXYGSshGPSq+4wI9RPCQIgBDCEnYh4OEng5sWIiJ62JMnCDEdZzhBgXSppTFKMMhdkjfAc06GRvSXHkmWu3RKj69gpQm9kgTUZ4grFYzdTzVzor9yXusPdXeRvR3M6+AWIkrYv/STTP/q1O1SPRR1TVwqinWjKrOy1xS3RW1c/NTVZIcYuIU7lFcEPa0ctpnU2sSXbvqraPjrzpTsWruZbkp3n6tTu1g0ms6F7oG9vdDnwXNg7Jtle1Tq1SrZhcijx3sYp9OvYIaTlBHQ69yi3s8GMK4Mx6Np0mqkcs02/gyjOcPJOmZag==</latexit>
sha1_base64="UzCaRJ8V1sjzxOMovtQA9vOFOA4=">AAAC8HicjVLLSsNAFD2Nr/quunQTLIKrknTTuhALblzWYh9QiyTptIbmxWQilOJHuBNdiFu/wJX/4R+oX+GdaQpq8XFD2jPnnnNn7p3YkefGwjBeMtrM7Nz8QnZxaXlldW09t7HZiMOEO6zuhF7IW7YVM88NWF24wmOtiDPLtz3WtAdHMt+8ZDx2w+BUDCPW8a1+4PZcxxJENWvsYL9UKp7n8kbBUKFPAzMF+cPnpxooqmHuHWfoIoSDBD4YAgjCHizE9LRhwkBEXAcj4jghV+UZrrBE3oRUjBQWsQP67dOqnbIBrWXNWLkd2sWjl5NTxy55QtJxwnI3XeUTVVmyP9UeqZrybEP6t9NaPrECF8T+5Zso/+uTvQj0UFY9uNRTpBjZnZNWSdRU5Mn1T10JqhARJ3GX8pywo5yTOevKE6ve5WwtlX9VSsnKtZNqE7z92p08wXjWdC/0GZjfL30aNIoF0yiYJ0a+UsY4stjGDvbo1kuo4BhV1NUu17jFnca1G+1eexhLtUzq2cKX0B4/AO6Tm3E=</latexit>
sha1_base64="TIbWKoP1kWGNeeylxiEfUiwnIPI=">AAAC8HicjVLLSsNAFD2Nr1pfVZeCBIvgqiRuWhdiwY3LVuwDapEkTjWYF5OJUIpLP8Cd6ELc+gVd+R/+gY+f8M40BR/4mJDkzLn3nJl7Z+zIc2NhGE8ZbWx8YnIqO52bmZ2bX8gvLjXiMOEOqzuhF/KWbcXMcwNWF67wWCvizPJtjzXts10Zb54zHrthcCB6Eev41kngdl3HEkQ199n2Vqm0eZQvGEVDDf07MFNQ2Hkc1N4uVwfVMP+KQxwjhIMEPhgCCMIeLMT0tGHCQERcB33iOCFXxRkukCNtQlmMMixiz+h7QrN2ygY0l56xUju0ikcvJ6WOddKElMcJy9V0FU+Us2R/8u4rT7m3Hv3t1MsnVuCU2L90o8z/6mQtAl2UVQ0u1RQpRlbnpC6J6orcuf6hKkEOEXESH1OcE3aUctRnXWliVbvsraXizypTsnLupLkJXn6tTu5g2Gs6F7oG5tdD/w4am0XTKJo1o1ApYziyWMEaNujUS6hgD1XU1SpXuMGtxrVr7U67H6ZqmVSzjE9De3gHQkWdtg==</latexit>

<latexit sha1_base64="s9QHnEmm0myZUMo2MA5EV8Cv3rI=">AAAC73icjVLJTsMwFJyGrZStwJFLRIXEqUq4wBGJC8ci0UVqK+S4bjFkI3GQqop/4IbggLjySfwB8BU8u6kEVCyOkoznvRn7PduLfZkqx3kpWDOzc/MLxcXS0vLK6lp5faORRlnCRZ1HfpS0PJYKX4airqTyRStOBAs8XzS9yyMdb16LJJVReKqGsegGbBDKvuRMEdXopHIQsLNyxak6ZtjTwM1BBfmoReV3dNBDBI4MAQRCKMI+GFJ62nDhICauixFxCSFp4gI3KJE2oyxBGYzYS/oOaNbO2ZDm2jM1ak6r+PQmpLSxQ5qI8hLCejXbxDPjrNmfvEfGU+9tSH8v9wqIVTgn9i/dJPO/Ol2LQh8HpgZJNcWG0dXx3CUzXdE7tz9VpcghJk7jHsUTwtwoJ322jSY1teveMhN/NZma1XOe52Z4+7U6vYNxr+lc6Bq43w99GjT2qq5TdU+cyuFBfiGK2MI2dunU93GIY9RQp1UucIt7PFhX1p31aD2NU61CrtnEl2E9fwCSrZoA</latexit>
sha1_base64="ll5ixb/zbA96H+AFqyS1ABFhSr4=">AAAC73icjVK5TsNAEH0xV7gDlDQWERJVZNOQjkg0lAGRQ0oiZJtNWOKL9RopivgHOgQFouUPqPgP/gD4CmY3RgIijrGcvH3z3uzOrN3Y54m0rOecMTE5NT2Tn52bX1hcWi6srNaTKBUeq3mRH4mm6yTM5yGrSS591owFcwLXZw23v6fyjQsmEh6FR3IQs07g9ELe5Z4jiaq3E94LnONC0SpZOsxxYGeguPv0eAiKalR4QxsniOAhRQCGEJKwDwcJPS3YsBAT18GQOEGI6zzDJebIm5KKkcIhtk+/PVq1MjaktaqZaLdHu/j0CnKa2CRPRDpBWO1m6nyqKyv2p9pDXVOdbUD/blYrIFbilNi/fB/K//pULxJdlHUPnHqKNaO687IqqZ6KOrn5qStJFWLiFD6hvCDsaefHnE3tSXTvaraOzr9opWLV2su0KV5/7U6dYDRruhf6DOzvlz4O6tsl2yrZB1axUsYo8ljHBrbo1ndQwT6qqNEuZ7jCDW6Nc+PauDPuR1Ijl3nW8CWMh3dcZpwH</latexit>
sha1_base64="2iSPbGBlJ6ok04PDJEIwed41MtA=">AAAC73icjVLLSsNAFD2Nr1pfVZeCBIvgqiRu7M6CG5ct2Ae0RZI41tG8TCZCKS7duxNdiFv/oCv/wz/w8RPemaagFh8Tkpw5954zc++MHbo8FobxnNEmJqemZ7Kzubn5hcWl/PJKPQ6SyGE1J3CDqGlbMXO5z2qCC5c1w4hZnu2yhn22J+ONCxbFPPAPRC9kHc/q+vyYO5Ygqt6OedezDvMFo2iooY8DMwWF3adB9f1qfVAJ8m9o4wgBHCTwwOBDEHZhIaanBRMGQuI66BMXEeIqznCJHGkTymKUYRF7Rt8uzVop69NcesZK7dAqLr0RKXVskiagvIiwXE1X8UQ5S/Yn777ylHvr0d9OvTxiBU6I/Us3yvyvTtYicIySqoFTTaFiZHVO6pKorsid65+qEuQQEifxEcUjwo5SjvqsK02sape9tVT8RWVKVs6dNDfB66/VyR0Me03nQtfA/H7o46C+XTSNolk1CuUShiOLNWxgi059B2Xso4IarXKKa9ziTjvXbrR77WGYqmVSzSq+DO3xA7AJnkw=</latexit>

0
<latexit sha1_base64="TA0lMnNje99cfY+WHlB5XFn1UVg=">AAAC5nicjVLLTsMwEJyGVymvwpVLRIXEqUq40GMlLhyLRB9SqVDiusWqm0SJg1RVfAE3BAeu/BV/AHwFa9eVgIqHoyTj2Z2xd+0wkSJTnvdScJaWV1bXiuuljc2t7Z1yabeVxXnKeJPFMk47YZBxKSLeVEJJ3klSHoxDydvh6FTH2zc8zUQcXahJwnvjYBiJgWCBIurcuypXvKpnhrsIfAsqsKMRl99xiT5iMOQYgyOCIiwRIKOnCx8eEuJ6mBKXEhImznGLEmlzyuKUERA7ou+QZl3LRjTXnplRM1pF0puS0sUhaWLKSwnr1VwTz42zZn/ynhpPvbcJ/UPrNSZW4ZrYv3TzzP/qdC0KA9RMDYJqSgyjq2PWJTdd0Tt3P1WlyCEhTuM+xVPCzCjnfXaNJjO1694GJv5qMjWr58zm5nj7tTq9g1mv6VzoGvjfD30RtI6rvlf1K/WavQ5F7OMAR3TmJ6jjDA00zfp3eMCjI51752mW6BSsYg9fhvP8AYvOlmg=</latexit>
<latexit sha1_base64="tqBnW07EXs5m+1LbvTdMtTHMzJM=">AAAC6nicjVK5TsNAEH0xVwhXgJLGIkKiimwa0hGJhjJB5JBChOzNJlisD9lrpCjiC+gQFIiWj6DiP/gD4CuY3TgSEHGMZfvtm3mzM7PrRsJLpGW95IyZ2bn5hfxiYWl5ZXWtuL7RTMI0ZrzBQhHGbddJuPAC3pCeFLwdxdzxXcFb7sWh8rcueZx4YXAihxHv+s4g8PoecyRRdeusWLLKljZzGtgZKB08Px2DrBYW33GKHkIwpPDBEUASFnCQ0NOBDQsRcV2MiIsJedrPcYUCaVOK4hThEHtB3wGtOhkb0FrlTLSa0S6C3piUJnZIE1JcTFjtZmp/qjMr9qfcI51T1Takv5vl8omVOCf2L90k8r861YtEHxXdg0c9RZpR3bEsS6qnoio3P3UlKUNEnMI98seEmVZO5mxqTaJ7V7N1tP9VRypWrVkWm+Lt1+5UBeNZ07nQNbC/H/o0aO6Vbats161StYKx5bGFbezSqe+jiiPU0NAVXOMWd4Ywbox742EcauQyzSa+mPH4AfJvmZg=</latexit>
sha1_base64="6xJbv2DxrZjeqtgd2qBj7kPD0MI=">AAAC5nicjVLLSsNAFD2Nr1pf1a0gwSK4Kokbu7PgxmUL9gG1SDKd1tBpEpKJUIpLV+5EF279iq78D//Ax094Z5qCWnxMSHLm3HvOzL0zbii8WFrWc8aYm19YXMou51ZW19Y38rnNehwkEeM1FoggarpOzIXn85r0pODNMOLOwBW84faPVbxxyaPYC/xTOQx5e+D0fK/rMUcSVbXO8wWraOlhzgI7BYWjp3H1/XpnXAnybzhDBwEYEgzA4UMSFnAQ09OCDQshcW2MiIsIeTrOcYUcaRPK4pThENunb49mrZT1aa48Y61mtIqgNyKliT3SBJQXEVarmTqeaGfF/uQ90p5qb0P6u6nXgFiJC2L/0k0z/6tTtUh0UdI1eFRTqBlVHUtdEt0VtXPzU1WSHELiFO5QPCLMtHLaZ1NrYl276q2j4y86U7FqztLcBK+/Vqd2MOk1nQtdA/v7oc+C+kHRtop2oVzCZGSxjV3s05kfoowTVFDT69/gDveGMG6Nh0mikUkVW/gyjMcPl/qatA==</latexit>

<latexit sha1_base64="vnd5zgD+4aU/cUnrpyqT36+9OpU=">AAAC7nicjVLLSsNAFD2Nr1pfVZdugkVwoWXSjV0W3LisYNpCLZKk0zqYF8lEKMVvcCe6ELf+kn+gfoV3pimoxceEJGfOvefM3Dvjxr5IJWMvBWNufmFxqbhcWlldW98ob2610ihLPG57kR8lHddJuS9Cbkshfd6JE+4Ers/b7tWxireveZKKKDyTo5j3AmcYioHwHEmUfcgOWO2iXGFVpoc5C6wcVJCPZlR+xzn6iOAhQwCOEJKwDwcpPV1YYIiJ62FMXEJI6DjHDUqkzSiLU4ZD7BV9hzTr5mxIc+WZarVHq/j0JqQ0sUeaiPISwmo1U8cz7azYn7zH2lPtbUR/N/cKiJW4JPYv3TTzvzpVi8QAdV2DoJpizajqvNwl011ROzc/VSXJISZO4T7FE8KeVk77bGpNqmtXvXV0/FVnKlbNvTw3w9uv1akdTHpN50LXwPp+6LOgVatarGqdskqjnl+IInawi3069SM0cIImbFpF4Bb3eDBi4854NJ4mqUYh12zjyzCePwCaNZh0</latexit>
sha1_base64="6RQC7ta5t/mP6368xyUZkK3bIQ0=">AAAC7nicjVK7TsMwFD0Nr1JehY4sEagSA1RJFzpWYmEsEmkrlQolqVus5qXEQaoqNnY2BANi5ZfYGIGv4NpNJaDi4SjJ8bn3HPte24k8ngjDeM5pc/MLi0v55cLK6tr6RnFzq5mEaewyyw29MG47dsI8HjBLcOGxdhQz23c81nKGRzLeumRxwsPgVIwi1vXtQcD73LUFUdaBsW9Uz4u7RsVQQ58FZgZ266Xy9QuARlh8xxl6COEihQ+GAIKwBxsJPR2YMBAR18WYuJgQV3GGKxRIm1IWowyb2CF9BzTrZGxAc+mZKLVLq3j0xqTUUSZNSHkxYbmaruKpcpbsT95j5Sn3NqK/k3n5xApcEPuXbpr5X52sRaCPmqqBU02RYmR1buaSqq7IneufqhLkEBEncY/iMWFXKad91pUmUbXL3toq/qoyJSvnbpab4u3X6uQOJr2mc6FrYH4/9FnQrFZMo2Ke0H2oYTLy2MYO9ujUD1HHMRqwaBWOG9zhXou0W+1Be5ykarlMU8KXoT19APvQmlA=</latexit>
sha1_base64="87/yCGeXZMHdSiAJ8YbqCFxA2eE=">AAAC7nicjVLLSsNAFD2Nr1pf1S7dBEvBhZakG7ssuHFZwT6gFknSaR2aJiGZCKW4c+9OdCFu/RL/QT9A1LUf4J1pCmrxMSHJmXPvOTP3ztiByyNhGI8pbWZ2bn4hvZhZWl5ZXcuub9QjPw4dVnN81w+bthUxl3usJrhwWTMImTWwXdaw+/sy3jhjYcR970gMA9YeWD2Pd7ljCaJqu8aOUTrJ5o2ioYY+DcwE5Cu5wsX7w/NT1c++4Rgd+HAQYwAGD4KwCwsRPS2YMBAQ18aIuJAQV3GGc2RIG1MWowyL2D59ezRrJaxHc+kZKbVDq7j0hqTUUSCNT3khYbmaruKxcpbsT94j5Sn3NqS/nXgNiBU4JfYv3STzvzpZi0AXZVUDp5oCxcjqnMQlVl2RO9c/VSXIISBO4g7FQ8KOUk76rCtNpGqXvbVU/EVlSlbOnSQ3xuuv1ckdjHtN50LXwPx+6NOgXiqaRtE8pPtQxniksYktbNOp76GCA1RRo1U4LnGNGy3QrrRb7W6cqqUSTQ5fhnb/ARU9nOo=</latexit>
0, 02
Re = 5772
<latexit sha1_base64="BtxfXxAwyfQbCyXYdFWCSwmqb2c=">AAAC8HicjVLLSsNAFD2Nr1pfVZdugkVwVZKCtBuh4MZlFfuAKpKk0xqaF5OJUIof4U50IW79I/9A/QrvTKegFh8Tkpw5954zc++MmwR+KizrJWfMzS8sLuWXCyura+sbxc2tVhpn3GNNLw5i3nGdlAV+xJrCFwHrJJw5oRuwtjs8kvH2NeOpH0dnYpSwi9AZRH7f9xxBVPuUHR5Uq5XLYskqW2qYs8DWoAQ9GnHxHefoIYaHDCEYIgjCARyk9HRhw0JC3AXGxHFCvooz3KBA2oyyGGU4xA7pO6BZV7MRzaVnqtQerRLQy0lpYo80MeVxwnI1U8Uz5SzZn7zHylPubUR/V3uFxApcEfuXbpr5X52sRaCPmqrBp5oSxcjqPO2Sqa7InZufqhLkkBAncY/inLCnlNM+m0qTqtplbx0Vf1WZkpVzT+dmePu1OrmDSa/pXOga2N8PfRa0KmXbKtsnVqle0xcijx3sYp9OvYo6jtFAU61yi3s8GNy4Mx6Np0mqkdOabXwZxvMHGsWZZg==</latexit>
sha1_base64="Kr9nXiBtOPNBDxhU1mOB5mUHXmw=">AAAC8HicjVLLSsNAFD2Nr1pfVZdugkVwVZKCtBux4MZlLfYBVSRJpzU0LyYToRQ/wp3oQtz6Ba78D/9A/QrvTFNQi48b0p4595w7c+/Ejjw3FobxktFmZufmF7KLuaXlldW1/PpGMw4T7rCGE3ohb9tWzDw3YA3hCo+1I84s3/ZYyx4cynzrkvHYDYMTMYzYmW/1A7fnOpYgqlVn+3vlcuk8XzCKhgp9GpgpKBw8P9VBUQvz7zhFFyEcJPDBEEAQ9mAhpqcDEwYi4s4wIo4TclWe4Qo58iakYqSwiB3Qb59WnZQNaC1rxsrt0C4evZycOnbIE5KOE5a76SqfqMqS/an2SNWUZxvSv53W8okVuCD2L99E+V+f7EWgh4rqwaWeIsXI7py0SqKmIk+uf+pKUIWIOIm7lOeEHeWczFlXnlj1LmdrqfyrUkpWrp1Um+Dt1+7kCcazpnuhz8D8funToFkqmkbRPDYK1QrGkcUWtrFLt15GFUeooaF2ucYt7jSu3Wj32sNYqmVSzya+hPb4AeRvm20=</latexit>
sha1_base64="2bYqW1HrY3ZuRkyntWJqfS5qarE=">AAAC8HicjVLLSsNAFD2Nr1pfVZeCBIvgqiSCtBux4MZlK/YBtUgSpxrMi8lEKMWlH+BOdCFu/YKu/A//wMdPeGeagg98TEhy5tx7zsy9M3bkubEwjKeMNjY+MTmVnc7NzM7NL+QXlxpxmHCH1Z3QC3nLtmLmuQGrC1d4rBVxZvm2x5r22a6MN88Zj90wOBC9iHV86yRwu65jCaKa+2x7q1TaPMoXjKKhhv4dmCko7DwOam+Xq4NqmH/FIY4RwkECHwwBBGEPFmJ62jBhICKugz5xnJCr4gwXyJE2oSxGGRaxZ/Q9oVk7ZQOaS89YqR1axaOXk1LHOmlCyuOE5Wq6iifKWbI/efeVp9xbj/526uUTK3BK7F+6UeZ/dbIWgS7KqgaXaooUI6tzUpdEdUXuXP9QlSCHiDiJjynOCTtKOeqzrjSxql321lLxZ5UpWTl30twEL79WJ3cw7DWdC10D8+uhfweNzaJpFM2aUaiUMRxZrGANG3TqJVSwhyrqapUr3OBW49q1dqfdD1O1TKpZxqehPbwDOCGdsg==</latexit>

0, 04
<latexit sha1_base64="P107o87R9bwg7MYsBDe6An8Rci8=">AAAC7nicjVLLSsNAFD2Nr1pfVZdugkVwoWUigl0W3LisYNpCFUnSaR2aF8lEKMVvcCe6ELf+kn+gfoV3pimoxceEJGfOvefM3Dvjxr5IJWMvBWNmdm5+obhYWlpeWV0rr2800yhLPG57kR8lbddJuS9Cbkshfd6OE+4Ers9b7uBYxVvXPElFFJ7JYcwvAqcfip7wHEmUvc/22OFlucKqTA9zGlg5qCAfjaj8jnN0EcFDhgAcISRhHw5SejqwwBATd4ERcQkhoeMcNyiRNqMsThkOsQP69mnWydmQ5soz1WqPVvHpTUhpYoc0EeUlhNVqpo5n2lmxP3mPtKfa25D+bu4VECtxRexfuknmf3WqFokearoGQTXFmlHVeblLpruidm5+qkqSQ0ycwl2KJ4Q9rZz02dSaVNeueuvo+KvOVKyae3luhrdfq1M7GPeazoWugfX90KdB86Bqsap1yir1Wn4hitjCNnbp1I9QxwkasGkVgVvc48GIjTvj0XgapxqFXLOJL8N4/gCfQZh2</latexit>
sha1_base64="hCIuMJzmceNp4TbMne4dvGYE91Y=">AAAC7nicjVLLSsNAFD2Nr1pf1S7dBKXgQksigl0W3LisYNpCLZKk0zo0L5KJUIo79+5EF+LWX3LnUv0K70xTUIuPCUnOnHvPmbl3xok8ngjDeM5pM7Nz8wv5xcLS8srqWnF9o5GEaewyyw29MG45dsI8HjBLcOGxVhQz23c81nQGRzLevGRxwsPgVAwj1vHtfsB73LUFUdaesWscnBe3jYqhhj4NzAxs10rl6xcA9bD4jjN0EcJFCh8MAQRhDzYSetowYSAiroMRcTEhruIMVyiQNqUsRhk2sQP69mnWztiA5tIzUWqXVvHojUmpo0yakPJiwnI1XcVT5SzZn7xHylPubUh/J/PyiRW4IPYv3STzvzpZi0APVVUDp5oixcjq3MwlVV2RO9c/VSXIISJO4i7FY8KuUk76rCtNomqXvbVV/FVlSlbO3Sw3xduv1ckdjHtN50LXwPx+6NOgsV8xjYp5QvehivHIYxNb2KFTP0QNx6jDolU4bnCHey3SbrUH7XGcquUyTQlfhvb0AQDrmlI=</latexit>
sha1_base64="G2HugOK/JnqA04rJoiZvnrzc4g0=">AAAC7nicjVLLSsNAFD2Nr1pf1S7dBEvBhZZEBLssuHFZwT6gFknSaR2aJiGZCKW4c+9OdCFu/RL/QT9A1LUf4J1pCmrxMSHJmXPvOTP3ztiByyNhGI8pbWp6ZnYuPZ9ZWFxaXsmurtUiPw4dVnV81w8bthUxl3usKrhwWSMImdW3XVa3e/syXj9jYcR970gMAtbqW12Pd7hjCaKq28aWsXuSzRtFQw19EpgJyJdzhYv3h+enip99wzHa8OEgRh8MHgRhFxYiepowYSAgroUhcSEhruIM58iQNqYsRhkWsT36dmnWTFiP5tIzUmqHVnHpDUmpo0Aan/JCwnI1XcVj5SzZn7yHylPubUB/O/HqEytwSuxfunHmf3WyFoEOSqoGTjUFipHVOYlLrLoid65/qkqQQ0CcxG2Kh4QdpRz3WVeaSNUue2up+IvKlKycO0lujNdfq5M7GPWazoWugfn90CdBbadoGkXzkO5DCaORxjo2sEmnvocyDlBBlVbhuMQ1brRAu9JutbtRqpZKNDl8Gdr9BxpJnOw=</latexit>

Re = 1772
<latexit sha1_base64="aEk0JqtZn0n2SsDqn4eLd7Nkgxg=">AAAC8HicjVLLSsNAFD2Nr1pfVZdugkVwVZJu2o1QcOOyin1ALZKk0xqaF5OJUIof4U50IW79I/9A/QrvTKegFh8Tkpw5954zc++MmwR+KizrJWcsLC4tr+RXC2vrG5tbxe2dVhpn3GNNLw5i3nGdlAV+xJrCFwHrJJw5oRuwtjs6lvH2NeOpH0fnYpywXugMI3/ge44gqn3GjuxqtXJZLFllSw1zHtgalKBHIy6+4wJ9xPCQIQRDBEE4gIOUni5sWEiI62FCHCfkqzjDDQqkzSiLUYZD7Ii+Q5p1NRvRXHqmSu3RKgG9nJQmDkgTUx4nLFczVTxTzpL9yXuiPOXexvR3tVdIrMAVsX/pZpn/1claBAaoqRp8qilRjKzO0y6Z6orcufmpKkEOCXES9ynOCXtKOeuzqTSpql321lHxV5UpWTn3dG6Gt1+rkzuY9prOha6B/f3Q50GrUratsn1qleo1fSHy2MM+DunUq6jjBA001Sq3uMeDwY0749F4mqYaOa3ZxZdhPH8AEKGZYg==</latexit>
sha1_base64="0w2MSd+FEARMwWcA1Ua+fyMQCSs=">AAAC8HicjVLLTsJAFD3UF+ILdemmkZi4Ii0b2BhJ3LhEIo8EiWnLgA19ZTo1IcSPcGd0Ydz6Ba78D/9A/QrvDCVRiY/bFM6ce86duXdqR54bC8N4yWhz8wuLS9nl3Mrq2vpGfnOrGYcJd1jDCb2Qt20rZp4bsIZwhcfaEWeWb3usZQ+PZL51yXjshsGpGEWs61uDwO27jiWIatXZgVkul87zBaNoqNBngZmCwuHzUx0UtTD/jjP0EMJBAh8MAQRhDxZiejowYSAirosxcZyQq/IMV8iRNyEVI4VF7JB+B7TqpGxAa1kzVm6HdvHo5eTUsUeekHScsNxNV/lEVZbsT7XHqqY824j+7bSWT6zABbF/+abK//pkLwJ9VFQPLvUUKUZ256RVEjUVeXL9U1eCKkTESdyjPCfsKOd0zrryxKp3OVtL5V+VUrJy7aTaBG+/didPMJk13Qt9Bub3S58FzVLRNIrmiVGoVjCJLHawi3269TKqOEYNDbXLNW5xp3HtRrvXHiZSLZN6tvEltMcP2kubaQ==</latexit>
sha1_base64="fN2MnOAIh/6GzPJlBmmvuYxkaP4=">AAAC8HicjVLLSsNAFD2Nr/quuhQkWARXJemm3YgFNy5V7ANqkSSd1tC8mEyEUlz6Ae5EF+LWL3Dlf/gHPn7CO9MU1OJjQpIz595zZu6dsSPPjYVhPGe0icmp6Zns7Nz8wuLScm5ltRaHCXdY1Qm9kDdsK2aeG7CqcIXHGhFnlm97rG739mS8fs547IbBsehHrOVb3cDtuI4liKofsR2zVCqe5vJGwVBDHwdmCvK7T4+H75cbjwdh7g0naCOEgwQ+GAIIwh4sxPQ0YcJARFwLA+I4IVfFGS4wR9qEshhlWMT26NulWTNlA5pLz1ipHVrFo5eTUscWaULK44TlarqKJ8pZsj95D5Sn3Fuf/nbq5RMrcEbsX7pR5n91shaBDsqqBpdqihQjq3NSl0R1Re5c/1SVIIeIOInbFOeEHaUc9VlXmljVLntrqfiLypSsnDtpboLXX6uTOxj2ms6FroH5/dDHQa1YMI2CeWjkK2UMRxbr2MQ2nXoJFezjAFW1yhVucKtx7Vq70+6HqVom1azhy9AePgAt/Z2u</latexit>

a) 0, 06 k
<latexit sha1_base64="IXFE5t28m/IQk3BKhVIEqqvaw/w=">AAAC6nicjVLLSsNAFD2Nr/quunQTLIKrkrixy4Ibly3YB9QiSTqtQycPJhOhFL/AnehC3PpT/oH6Fd6ZpqAWHxOSnDn3njv3MX4ieKoc56VgLSwuLa8UV9fWNza3tks7u600zmTAmkEsYtnxvZQJHrGm4kqwTiKZF/qCtf3Rqba3r5lMeRydq3HCeqE3jPiAB54iqjG6LJWdimOWPQ/cHJSRr3pcescF+ogRIEMIhgiKsICHlJ4uXDhIiOthQpwkxI2d4QZrpM3Ii5GHR+yIvkPadXM2or2OmRp1QKcIeiUpbRySJiY/SVifZht7ZiJr9qfYExNT5zamv5/HColVuCL2L93M8786XYvCAFVTA6eaEsPo6oI8Sma6ojO3P1WlKEJCnMZ9skvCgVHO+mwbTWpq1731jP3VeGpW74PcN8Pbr9XpDKa9prnQNXC/D30etI4rrlNxG065Vs0vRBH7OMARTf0ENZyhjqbJ4Bb3eLCEdWc9Wk9TV6uQa/bwZVnPH72nl8w=</latexit>
<latexit sha1_base64="xkRW35A/XieE00d3e7jnwUiDqf4=">AAAC6nicjVK5TsNAEH0xVwhXgJLGIkKiimwa0hGJhjJB5JBChOzNBqysD9lrpCjiC+gQFIiWj6DiP/gD4CuY3TgSEHGMZfvtm3mzM7PrRsJLpGW95IyZ2bn5hfxiYWl5ZXWtuL7RTMI0ZrzBQhHGbddJuPAC3pCeFLwdxdzxXcFb7uBQ+VuXPE68MDiRw4h3fec88PoecyRR9cFZsWSVLW3mNLAzUDp4fjoGWS0svuMUPYRgSOGDI4AkLOAgoacDGxYi4roYERcT8rSf4woF0qYUxSnCIXZA33NadTI2oLXKmWg1o10EvTEpTeyQJqS4mLDazdT+VGdW7E+5Rzqnqm1IfzfL5RMrcUHsX7pJ5H91qheJPiq6B496ijSjumNZllRPRVVufupKUoaIOIV75I8JM62czNnUmkT3rmbraP+rjlSsWrMsNsXbr92pCsazpnOha2B/P/Rp0Nwr21bZrlulagVjy2ML29ilU99HFUeooaEruMYt7gxh3Bj3xsM41Mhlmk18MePxA4dgmdM=</latexit>
sha1_base64="Svt6pl9lIwlraG5BJAsRk5cwB2s=">AAAC6nicjVLLSsNAFD2Nr1pfVZeCBIvgqiRu7M6CG5ct2AfUIkk6raGTB8lEKMWlK3eiC3HrR3Tlf/gHPn7CO9MU1OJjQpIz595zc8+d2CF3Y2EYzxltZnZufiG7mFtaXlldy69v1OMgiRxWcwIeRE3bihl3fVYTruCsGUbM8mzOGnb/SMYbFyyK3cA/EYOQtT2r57td17EEUdX+Wb5gFA219GlgpqBw+DSqvl9tjypB/g2n6CCAgwQeGHwIwhwWYrpaMGEgJK6NIXERIVfFGS6RI21CWYwyLGL79OzRrpWyPu1lzVipHfoKpzsipY5d0gSUFxGWX9NVPFGVJftT7aGqKXsb0NtOa3nECpwT+5dukvlfnfQi0EVJeXDJU6gY6c5JqyRqKrJz/ZMrQRVC4iTuUDwi7CjlZM660sTKu5ytpeIvKlOycu+kuQlef3UnOxjPms6FfgPz+6FPg/p+0TSKZtUolEsYryy2sIM9OvUDlHGMCmqqg2vc4k7j2o12rz2MU7VMqtnEl6U9fgDbA5wY</latexit>

0 0, 5 1 1, 5
<latexit sha1_base64="6uXk0oY2HYwRtq1mYanQgJxunyA=">AAAC7nicjVLLSsNAFD2Nr1pfVZdugkVwoWXiQrssuHFZwbSFKpKk0zo0L5KJUIrf4E50IW79Jf9A/QrvTFNQi48JSc6ce8+ZuXfGjX2RSsZeCsbM7Nz8QnGxtLS8srpWXt9oplGWeNz2Ij9K2q6Tcl+E3JZC+rwdJ9wJXJ+33MGxireueZKKKDyTw5hfBE4/FD3hOZIoe5/tscPLcoVVmR7mNLByUEE+GlH5HefoIoKHDAE4QkjCPhyk9HRggSEm7gIj4hJCQsc5blAibUZZnDIcYgf07dOsk7MhzZVnqtUereLTm5DSxA5pIspLCKvVTB3PtLNif/IeaU+1tyH93dwrIFbiiti/dJPM/+pULRI91HQNgmqKNaOq83KXTHdF7dz8VJUkh5g4hbsUTwh7Wjnps6k1qa5d9dbR8VedqVg19/LcDG+/Vqd2MO41nQtdA+v7oU+D5kHVYlXrlFXqtfxCFLGFbezSqR+hjhM0YNMqAre4x4MRG3fGo/E0TjUKuWYTX4bx/AGkTZh4</latexit>
sha1_base64="XSjLL24WOpsAsdyFZD9QTvV/CTY=">AAAC7nicjVK7TsMwFD0Nr1JehY4sEagSA1QJA3SsxMJYJNJWKhVKUrdYzUuJg1RVbOxsCAbEyi+xMQJfwbWbSkDFw1GS43PvOfa9thN5PBGG8ZzTZmbn5hfyi4Wl5ZXVteL6RiMJ09hllht6Ydxy7IR5PGCW4MJjrShmtu94rOkMjmS8ecnihIfBqRhGrOPb/YD3uGsLoqw9Y9c4OC9uGxVDDX0amBnYrpXK1y8A6mHxHWfoIoSLFD4YAgjCHmwk9LRhwkBEXAcj4mJCXMUZrlAgbUpZjDJsYgf07dOsnbEBzaVnotQureLRG5NSR5k0IeXFhOVquoqnylmyP3mPlKfc25D+TublEytwQexfuknmf3WyFoEeqqoGTjVFipHVuZlLqroid65/qkqQQ0ScxF2Kx4RdpZz0WVeaRNUue2ur+KvKlKycu1luirdfq5M7GPeazoWugfn90KdBY79iGhXzhO5DFeORxya2sEOnfogajlGHRatw3OAO91qk3WoP2uM4VctlmhK+DO3pAwX3mlQ=</latexit>
sha1_base64="B6ttGm3SCXo9xEpxyN4QoYkGFpc=">AAAC7nicjVLLSsNAFD2Nr1pf1S7dBEvBhZbEhXZZcOOygn1ALZKk0zo0TUIyEUpx596d6ELc+iX+g36AqGs/wDvTFNTiY0KSM+fec2bunbEDl0fCMB5T2tT0zOxcej6zsLi0vJJdXatFfhw6rOr4rh82bCtiLvdYVXDhskYQMqtvu6xu9/ZlvH7Gwoj73pEYBKzVt7oe73DHEkRVt40tY/ckmzeKhhr6JDATkC/nChfvD89PFT/7hmO04cNBjD4YPAjCLixE9DRhwkBAXAtD4kJCXMUZzpEhbUxZjDIsYnv07dKsmbAezaVnpNQOreLSG5JSR4E0PuWFhOVquorHylmyP3kPlafc24D+duLVJ1bglNi/dOPM/+pkLQIdlFQNnGoKFCOrcxKXWHVF7lz/VJUgh4A4idsUDwk7Sjnus640kapd9tZS8ReVKVk5d5LcGK+/Vid3MOo1nQtdA/P7oU+C2k7RNIrmId2HEkYjjXVsYJNOfQ9lHKCCKq3CcYlr3GiBdqXdanejVC2VaHL4MrT7Dx9VnO4=</latexit>

<latexit sha1_base64="TA0lMnNje99cfY+WHlB5XFn1UVg=">AAAC5nicjVLLTsMwEJyGVymvwpVLRIXEqUq40GMlLhyLRB9SqVDiusWqm0SJg1RVfAE3BAeu/BV/AHwFa9eVgIqHoyTj2Z2xd+0wkSJTnvdScJaWV1bXiuuljc2t7Z1yabeVxXnKeJPFMk47YZBxKSLeVEJJ3klSHoxDydvh6FTH2zc8zUQcXahJwnvjYBiJgWCBIurcuypXvKpnhrsIfAsqsKMRl99xiT5iMOQYgyOCIiwRIKOnCx8eEuJ6mBKXEhImznGLEmlzyuKUERA7ou+QZl3LRjTXnplRM1pF0puS0sUhaWLKSwnr1VwTz42zZn/ynhpPvbcJ/UPrNSZW4ZrYv3TzzP/qdC0KA9RMDYJqSgyjq2PWJTdd0Tt3P1WlyCEhTuM+xVPCzCjnfXaNJjO1694GJv5qMjWr58zm5nj7tTq9g1mv6VzoGvjfD30RtI6rvlf1K/WavQ5F7OMAR3TmJ6jjDA00zfp3eMCjI51752mW6BSsYg9fhvP8AYvOlmg=</latexit>
<latexit sha1_base64="tqBnW07EXs5m+1LbvTdMtTHMzJM=">AAAC6nicjVK5TsNAEH0xVwhXgJLGIkKiimwa0hGJhjJB5JBChOzNJlisD9lrpCjiC+gQFIiWj6DiP/gD4CuY3TgSEHGMZfvtm3mzM7PrRsJLpGW95IyZ2bn5hfxiYWl5ZXWtuL7RTMI0ZrzBQhHGbddJuPAC3pCeFLwdxdzxXcFb7sWh8rcueZx4YXAihxHv+s4g8PoecyRRdeusWLLKljZzGtgZKB08Px2DrBYW33GKHkIwpPDBEUASFnCQ0NOBDQsRcV2MiIsJedrPcYUCaVOK4hThEHtB3wGtOhkb0FrlTLSa0S6C3piUJnZIE1JcTFjtZmp/qjMr9qfcI51T1Takv5vl8omVOCf2L90k8r861YtEHxXdg0c9RZpR3bEsS6qnoio3P3UlKUNEnMI98seEmVZO5mxqTaJ7V7N1tP9VRypWrVkWm+Lt1+5UBeNZ07nQNbC/H/o0aO6Vbats161StYKx5bGFbezSqe+jiiPU0NAVXOMWd4Ywbox742EcauQyzSa+mPH4AfJvmZg=</latexit>
sha1_base64="6xJbv2DxrZjeqtgd2qBj7kPD0MI=">AAAC5nicjVLLSsNAFD2Nr1pf1a0gwSK4Kokbu7PgxmUL9gG1SDKd1tBpEpKJUIpLV+5EF279iq78D//Ax094Z5qCWnxMSHLm3HvOzL0zbii8WFrWc8aYm19YXMou51ZW19Y38rnNehwkEeM1FoggarpOzIXn85r0pODNMOLOwBW84faPVbxxyaPYC/xTOQx5e+D0fK/rMUcSVbXO8wWraOlhzgI7BYWjp3H1/XpnXAnybzhDBwEYEgzA4UMSFnAQ09OCDQshcW2MiIsIeTrOcYUcaRPK4pThENunb49mrZT1aa48Y61mtIqgNyKliT3SBJQXEVarmTqeaGfF/uQ90p5qb0P6u6nXgFiJC2L/0k0z/6tTtUh0UdI1eFRTqBlVHUtdEt0VtXPzU1WSHELiFO5QPCLMtHLaZ1NrYl276q2j4y86U7FqztLcBK+/Vqd2MOk1nQtdA/v7oc+C+kHRtop2oVzCZGSxjV3s05kfoowTVFDT69/gDveGMG6Nh0mikUkVW/gyjMcPl/qatA==</latexit>

<latexit sha1_base64="9FS8I7JF1eXwRedl9+pWkYxP9D8=">AAAC6nicjVLLSsNAFD2Nr/quunQTLIKrkrixy4Ibly3YB9QiSTqtQycPkolQil/gTnQhbv0p/0D9Cu9Mp6AWHxOSnDn3njv3MX4ieCYd56VgLSwuLa8UV9fWNza3tks7u60sztOANYNYxGnH9zImeMSakkvBOknKvNAXrO2PTpW9fc3SjMfRuRwnrBd6w4gPeOBJohruZansVBy97HngGlCGWfW49I4L9BEjQI4QDBEkYQEPGT1duHCQENfDhLiUENd2hhuskTYnL0YeHrEj+g5p1zVsRHsVM9PqgE4R9KaktHFImpj8UsLqNFvbcx1ZsT/FnuiYKrcx/X0TKyRW4orYv3Qzz//qVC0SA1R1DZxqSjSjqgtMlFx3RWVuf6pKUoSEOIX7ZE8JB1o567OtNZmuXfXW0/ZX7alYtQ+Mb463X6tTGUx7TXOha+B+H/o8aB1XXKfiNpxyrWouRBH7OMARTf0ENZyhjqbO4Bb3eLCEdWc9Wk9TV6tgNHv4sqznDytLl5I=</latexit>
<latexit sha1_base64="8mvYjjA4+aS01/ys0uxYPCVLe0k=">AAAC6nicjVK5TsNAEH0xVwhXgJLGIkKiimwa0hGJhjJB5JBChOzNJlisD9lrpCjiC+gQFIiWj6DiP/gD4CuY3TgSEHGMZfvtm3mzM7PrRsJLpGW95IyZ2bn5hfxiYWl5ZXWtuL7RTMI0ZrzBQhHGbddJuPAC3pCeFLwdxdzxXcFb7sWh8rcueZx4YXAihxHv+s4g8PoecyRRdfusWLLKljZzGtgZKB08Px2DrBYW33GKHkIwpPDBEUASFnCQ0NOBDQsRcV2MiIsJedrPcYUCaVOK4hThEHtB3wGtOhkb0FrlTLSa0S6C3piUJnZIE1JcTFjtZmp/qjMr9qfcI51T1Takv5vl8omVOCf2L90k8r861YtEHxXdg0c9RZpR3bEsS6qnoio3P3UlKUNEnMI98seEmVZO5mxqTaJ7V7N1tP9VRypWrVkWm+Lt1+5UBeNZ07nQNbC/H/o0aO6Vbats161StYKx5bGFbezSqe+jiiPU0NAVXOMWd4Ywbox742EcauQyzSa+mPH4AfT1mZk=</latexit>
sha1_base64="7g2YK/N7y2oDo16IsjNGlzesNzY=">AAAC6nicjVLLSsNAFD2Nr1pfVZeCBIvgqiRu7M6CG5ct2AfUIkk6raF5MZkIpbh05U50IW79iK78D//Ax094Z5qCWnxMSHLm3Htu7rkTO/LcWBjGc0abmZ2bX8gu5paWV1bX8usb9ThMuMNqTuiFvGlbMfPcgNWEKzzWjDizfNtjDbt/JOONC8ZjNwxOxCBibd/qBW7XdSxBVNU8yxeMoqGWPg3MFBQOn0bV96vtUSXMv+EUHYRwkMAHQwBB2IOFmK4WTBiIiGtjSBwn5Ko4wyVypE0oi1GGRWyfnj3atVI2oL2sGSu1Q1/x6Oak1LFLmpDyOGH5NV3FE1VZsj/VHqqasrcBve20lk+swDmxf+kmmf/VSS8CXZSUB5c8RYqR7py0SqKmIjvXP7kSVCEiTuIOxTlhRyknc9aVJlbe5WwtFX9RmZKVeyfNTfD6qzvZwXjWdC70G5jfD30a1PeLplE0q0ahXMJ4ZbGFHezRqR+gjGNUUFMdXOMWd5qn3Wj32sM4Vcukmk18WdrjB0inm94=</latexit>

<latexit sha1_base64="Uo+NQaiOOnlvcn0QuGxY7Q4cA5c=">AAAC7HicjVLLSsNAFD2Nr1pfVZdugkVwIWUiiF0W3LisaB9QRZJ0WkPTJEwmQil+gjvRhbj1n/wD9Su8M52CWnxMSHLm3HvOzL0zXhIGqWTsJWfNzM7NL+QXC0vLK6trxfWNRhpnwud1Pw5j0fLclIdBxOsykCFvJYK7Ay/kTa9/pOLNay7SII7O5DDhFwO3FwXdwHclUads7+CyWGJlpoc9DRwDSjCjFhffcY4OYvjIMABHBEk4hIuUnjYcMCTEXWBEnCAU6DjHDQqkzSiLU4ZLbJ++PZq1DRvRXHmmWu3TKiG9gpQ2dkgTU54grFazdTzTzor9yXukPdXehvT3jNeAWIkrYv/STTL/q1O1SHRR0TUEVFOiGVWdb1wy3RW1c/tTVZIcEuIU7lBcEPa1ctJnW2tSXbvqravjrzpTsWrum9wMb79Wp3Yw7jWdC10D5/uhT4PGftlhZeeElaoVcyHy2MI2dunUD1HFMWqo0yo93OIeD1Zk3VmP1tM41coZzSa+DOv5A2n0mAY=</latexit>
sha1_base64="w/fDVciiDMJOzrWLuPk/1pk+8U0=">AAAC7HicjVLLSsNAFD2Nr1pf1S7dBKXgQkoiiF0W3LisaB9QiyTptA5NkzCZCKW4c+tOdCFu/Sd3LtWv8M40BbX4mJDkzLn3nJl7Z9zI57G0rOeMMTM7N7+QXcwtLa+sruXXN+pxmAiP1bzQD0XTdWLm84DVJJc+a0aCOQPXZw23f6jijUsmYh4Gp3IYsfbA6QW8yz1HEnVi7e6f57etkqWHOQ3sFGxXCsXrFwDVMP+OM3QQwkOCARgCSMI+HMT0tGDDQkRcGyPiBCGu4wxXyJE2oSxGGQ6xffr2aNZK2YDmyjPWao9W8ekVpDRRJE1IeYKwWs3U8UQ7K/Yn75H2VHsb0t9NvQbESlwQ+5dukvlfnapFoouyroFTTZFmVHVe6pLorqidm5+qkuQQEadwh+KCsKeVkz6bWhPr2lVvHR1/1ZmKVXMvzU3w9mt1agfjXtO50DWwvx/6NKjvlWyrZB/TfShjPLLYxBZ26NQPUMERqqjRKj3c4A73RmDcGg/G4zjVyKSaAr4M4+kDy4+Z4g==</latexit>
sha1_base64="fN35LE3bsxmoo7xyt2HrMs5cFUk=">AAAC7HicjVLLSsNAFD2Nr1pf1S7dBEvBhZREELssuHFZ0T6gFknSaQ1NkzCZCKW4c+tOdCFu/RT/QT9A1LUf4J1pCmrxMSHJmXPvOTP3ztih50bCMB5T2tT0zOxcej6zsLi0vJJdXatFQcwdVnUCL+AN24qY5/qsKlzhsUbImdW3PVa3e3syXj9jPHID/0gMQtbqW13f7biOJYg6NLZ2TrJ5o2iooU8CMwH5cq5w8f7w/FQJsm84RhsBHMTog8GHIOzBQkRPEyYMhMS1MCSOE3JVnOEcGdLGlMUowyK2R98uzZoJ69NcekZK7dAqHr2clDoKpAkojxOWq+kqHitnyf7kPVSecm8D+tuJV59YgVNi/9KNM/+rk7UIdFBSNbhUU6gYWZ2TuMSqK3Ln+qeqBDmExEncpjgn7CjluM+60kSqdtlbS8VfVKZk5dxJcmO8/lqd3MGo13QudA3M74c+CWrbRdMomgd0H0oYjTTWsYFNOvVdlLGPCqq0SheXuMaN5mtX2q12N0rVUokmhy9Du/8A5O2cfA==</latexit> <latexit sha1_base64="E3zwC36lHCIyS9XsPipapBidcmk=">AAAC7HicjVLLSsNAFD2Nr1pfVZdugkVwISURxC4LblxWtA+oIsl0WoemSZhMhFL8BHeiC3HrP/kH6ld4Z5qCWnxMSHLm3HvOzL0zfhyIRDnOS86amZ2bX8gvFpaWV1bXiusbjSRKJeN1FgWRbPlewgMR8roSKuCtWHJv4Ae86fePdLx5zWUiovBMDWN+MfB6oegK5imiTt29g8tiySk7ZtjTwM1ACdmoRcV3nKODCAwpBuAIoQgH8JDQ04YLBzFxFxgRJwkJE+e4QYG0KWVxyvCI7dO3R7N2xoY0156JUTNaJaBXktLGDmkiypOE9Wq2iafGWbM/eY+Mp97bkP5+5jUgVuGK2L90k8z/6nQtCl1UTA2CaooNo6tjmUtquqJ3bn+qSpFDTJzGHYpLwswoJ322jSYxteveeib+ajI1q+csy03x9mt1egfjXtO50DVwvx/6NGjsl12n7J44pWoluxB5bGEbu3Tqh6jiGDXUaZUebnGPByu07qxH62mcauUyzSa+DOv5A2x8mAc=</latexit>
sha1_base64="65ICjlvNolUkCwhllc/UR89cuFA=">AAAC7HicjVLLSsNAFD2Nr1pf1S7dBKXgQkoiiF0W3LisaB9QiyTptA5NkzCZCKW4c+tOdCFu/Sd3LtWv8M40BbX4mJDkzLn3nJl7Z9zI57G0rOeMMTM7N7+QXcwtLa+sruXXN+pxmAiP1bzQD0XTdWLm84DVJJc+a0aCOQPXZw23f6jijUsmYh4Gp3IYsfbA6QW8yz1HEnVi7+6f57etkqWHOQ3sFGxXCsXrFwDVMP+OM3QQwkOCARgCSMI+HMT0tGDDQkRcGyPiBCGu4wxXyJE2oSxGGQ6xffr2aNZK2YDmyjPWao9W8ekVpDRRJE1IeYKwWs3U8UQ7K/Yn75H2VHsb0t9NvQbESlwQ+5dukvlfnapFoouyroFTTZFmVHVe6pLorqidm5+qkuQQEadwh+KCsKeVkz6bWhPr2lVvHR1/1ZmKVXMvzU3w9mt1agfjXtO50DWwvx/6NKjvlWyrZB/TfShjPLLYxBZ26NQPUMERqqjRKj3c4A73RmDcGg/G4zjVyKSaAr4M4+kDzheZ4w==</latexit>
sha1_base64="KVxl4/gl/mg9NfgzFBroTO/yVa4=">AAAC7HicjVLLSsNAFD2Nr1pf1S7dBEvBhZREELssuHFZ0T6gFknSaQ1NkzCZCKW4c+tOdCFu/RT/QT9A1LUf4J1pCmrxMSHJmXPvOTP3ztih50bCMB5T2tT0zOxcej6zsLi0vJJdXatFQcwdVnUCL+AN24qY5/qsKlzhsUbImdW3PVa3e3syXj9jPHID/0gMQtbqW13f7biOJYg6NLd2TrJ5o2iooU8CMwH5cq5w8f7w/FQJsm84RhsBHMTog8GHIOzBQkRPEyYMhMS1MCSOE3JVnOEcGdLGlMUowyK2R98uzZoJ69NcekZK7dAqHr2clDoKpAkojxOWq+kqHitnyf7kPVSecm8D+tuJV59YgVNi/9KNM/+rk7UIdFBSNbhUU6gYWZ2TuMSqK3Ln+qeqBDmExEncpjgn7CjluM+60kSqdtlbS8VfVKZk5dxJcmO8/lqd3MGo13QudA3M74c+CWrbRdMomgd0H0oYjTTWsYFNOvVdlLGPCqq0SheXuMaN5mtX2q12N0rVUokmhy9Du/8A53WcfQ==</latexit>

/k =
0
<latexit sha1_base64="s9QHnEmm0myZUMo2MA5EV8Cv3rI=">AAAC73icjVLJTsMwFJyGrZStwJFLRIXEqUq4wBGJC8ci0UVqK+S4bjFkI3GQqop/4IbggLjySfwB8BU8u6kEVCyOkoznvRn7PduLfZkqx3kpWDOzc/MLxcXS0vLK6lp5faORRlnCRZ1HfpS0PJYKX4airqTyRStOBAs8XzS9yyMdb16LJJVReKqGsegGbBDKvuRMEdXopHIQsLNyxak6ZtjTwM1BBfmoReV3dNBDBI4MAQRCKMI+GFJ62nDhICauixFxCSFp4gI3KJE2oyxBGYzYS/oOaNbO2ZDm2jM1ak6r+PQmpLSxQ5qI8hLCejXbxDPjrNmfvEfGU+9tSH8v9wqIVTgn9i/dJPO/Ol2LQh8HpgZJNcWG0dXx3CUzXdE7tz9VpcghJk7jHsUTwtwoJ322jSY1teveMhN/NZma1XOe52Z4+7U6vYNxr+lc6Bq43w99GjT2qq5TdU+cyuFBfiGK2MI2dunU93GIY9RQp1UucIt7PFhX1p31aD2NU61CrtnEl2E9fwCSrZoA</latexit>
sha1_base64="ll5ixb/zbA96H+AFqyS1ABFhSr4=">AAAC73icjVK5TsNAEH0xV7gDlDQWERJVZNOQjkg0lAGRQ0oiZJtNWOKL9RopivgHOgQFouUPqPgP/gD4CmY3RgIijrGcvH3z3uzOrN3Y54m0rOecMTE5NT2Tn52bX1hcWi6srNaTKBUeq3mRH4mm6yTM5yGrSS591owFcwLXZw23v6fyjQsmEh6FR3IQs07g9ELe5Z4jiaq3E94LnONC0SpZOsxxYGeguPv0eAiKalR4QxsniOAhRQCGEJKwDwcJPS3YsBAT18GQOEGI6zzDJebIm5KKkcIhtk+/PVq1MjaktaqZaLdHu/j0CnKa2CRPRDpBWO1m6nyqKyv2p9pDXVOdbUD/blYrIFbilNi/fB/K//pULxJdlHUPnHqKNaO687IqqZ6KOrn5qStJFWLiFD6hvCDsaefHnE3tSXTvaraOzr9opWLV2su0KV5/7U6dYDRruhf6DOzvlz4O6tsl2yrZB1axUsYo8ljHBrbo1ndQwT6qqNEuZ7jCDW6Nc+PauDPuR1Ijl3nW8CWMh3dcZpwH</latexit>
sha1_base64="2iSPbGBlJ6ok04PDJEIwed41MtA=">AAAC73icjVLLSsNAFD2Nr1pfVZeCBIvgqiRu7M6CG5ct2Ae0RZI41tG8TCZCKS7duxNdiFv/oCv/wz/w8RPemaagFh8Tkpw5954zc++MHbo8FobxnNEmJqemZ7Kzubn5hcWl/PJKPQ6SyGE1J3CDqGlbMXO5z2qCC5c1w4hZnu2yhn22J+ONCxbFPPAPRC9kHc/q+vyYO5Ygqt6OedezDvMFo2iooY8DMwWF3adB9f1qfVAJ8m9o4wgBHCTwwOBDEHZhIaanBRMGQuI66BMXEeIqznCJHGkTymKUYRF7Rt8uzVop69NcesZK7dAqLr0RKXVskiagvIiwXE1X8UQ5S/Yn777ylHvr0d9OvTxiBU6I/Us3yvyvTtYicIySqoFTTaFiZHVO6pKorsid65+qEuQQEifxEcUjwo5SjvqsK02sape9tVT8RWVKVs6dNDfB66/VyR0Me03nQtfA/H7o46C+XTSNolk1CuUShiOLNWxgi059B2Xso4IarXKKa9ziTjvXbrR77WGYqmVSzSq+DO3xA7AJnkw=</latexit>

kc = 1, 02
<latexit
sha1_base64="/zNlodeAICcl2KdvO06cOjc8ZKk=">AAAC83icjVLJSsRAEH3GfR/1IngJiuBpTLw4F0Hw4lHBUVEHSWLP2Ew2ko4win/hTfQgXv0gwQ9Qv8LXPRFccKmQ5PWretVdVe2nocyV4zz2WL19/QODQ8Mjo2PjE5OVqendPCmyQNSDJEyyfd/LRShjUVdShWI/zYQX+aHY89sb2r93JrJcJvGO6qSiEXmtWDZl4ClSB0e5bEXecnvNOa4sOFXHmP0duCVYWJ89fwJtK6m84ggnSBCgQASBGIo4hIeczyFcOEjJNXBBLiOSxi9wiRFqC0YJRnhk2/y2uDos2ZhrnTM36oC7hHwzKm0sUpMwLiPWu9nGX5jMmv0p94XJqc/W4d8vc0VkFU7J/qV7j/yvTtei0ETN1CBZU2oYXV1QZilMV/TJ7Q9VKWZIyWl8Qn9GHBjle59to8lN7bq3nvE/m0jN6nVQxhZ4+bU6fYJurzkXXgP369C/g92VqutU3W3ehxq6NoQ5zGOJU1/FOjaxhTp3iXGFG9xahXVt3Vn33VCrp9TM4JNZD2/EJJzd</latexit>
sha1_base64="Mft0SXKNharoFiSk/uN3toDEnNI=">AAAC83icjVLLSsNAFD2Nr7a+qm4EN0ERXNXEjd0IRTcuK9gH2iJJOq1D8yKZCLUIfoQ70YW4dePfCH6A+hXemabgAx8Tkpw5954zc++MHbo8FobxlNHGxicmp7K5/PTM7Nx8YWGxFgdJ5LCqE7hB1LCtmLncZ1XBhcsaYcQsz3ZZ3e7tyXj9jEUxD/xD0Q9Zy7O6Pu9wxxJEHTVj3vWszd6OcVJYM4qGGvp3YKZgrbx8/py7fNytBIU3NNFGAAcJPDD4EIRdWIjpOYYJAyFxLQyIiwhxFWe4QJ60CWUxyrCI7dG3S7PjlPVpLj1jpXZoFZfeiJQ61kkTUF5EWK6mq3iinCX7k/dAecq99elvp14esQKnxP6lG2X+VydrEeigpGrgVFOoGFmdk7okqity5/qHqgQ5hMRJ3KZ4RNhRylGfdaWJVe2yt5aKv6hMycq5k+YmeP21OrmDYa/pXOgamF8P/TuobRVNo2ge0H0oYTiyWMEqNujUt1HGPiqo0io+rnCDWy3RrrU77X6YqmVSzRI+De3hHUGonlo=</latexit>
sha1_base64="f7P01K3ZMOXE4IQEedDTmlvSXKQ=">AAAC83icjVLLSsNAFD2Nr1pfVZdugkVwVRM3diMU3LisYB/YFknSaR2aF5mJUIp/4U50IW79IP9A/QrvTFNQi48JSc6ce8+ZuXfGjX0upGW95Iy5+YXFpfxyYWV1bX2juLnVEFGaeKzuRX6UtFxHMJ+HrC659FkrTpgTuD5rusMTFW9es0TwKDyXo5h1A2cQ8j73HEnURUfwQeAcDI+ty2LJKlt6mLPAzkAJ2ahFxXd00EMEDykCMISQhH04EPS0YcNCTFwXY+ISQlzHGW5QIG1KWYwyHGKH9B3QrJ2xIc2Vp9Bqj1bx6U1IaWKPNBHlJYTVaqaOp9pZsT95j7Wn2tuI/m7mFRArcUXsX7pp5n91qhaJPiq6Bk41xZpR1XmZS6q7onZufqpKkkNMnMI9iieEPa2c9tnUGqFrV711dPxVZypWzb0sN8Xbr9WpHUx6TedC18D+fuizoHFYtq2yfWaVqpXsQuSxg13s06kfoYpT1FCnVULc4h4PRmrcGY/G0yTVyGWabXwZxvMHzaSbLw==</latexit>

/k = 0
, 004
k
<latexit
sha1_base64="4+OcFM3VnAXyaSh3bpA+sQROj9o=">AAAC93icjVLLattAFD1W+nDdR9xkU+hGJBS6KM4oBOpNwZBNly7UiSGPIikTd7BejEYF1+Q/sgvJonTbzwn0A5J+Rc+MZUgTmuQKSWfOvefO3HsnKhJVGiHOG97Cg4ePHjeftJ4+e/5isf1yaavMKx3LQZwnuR5GYSkTlcmBUSaRw0LLMI0SuR2NN61/+5vUpcqzz2ZSyL00HGXqUMWhIbW/W6pRGq6NP4h3Qmx8aa+KjnDm3wRBDVZ7r77/Bq2ft/9gFwfIEaNCCokMhjhBiJLPDgIIFOT2MCWniZTzSxyhRW3FKMmIkOyY3xFXOzWbcW1zlk4dc5eEr6bSxxtqcsZpYrub7/yVy2zZ/+Weupz2bBP+ozpXStbgK9m7dPPI++psLQaH6LoaFGsqHGOri+ssleuKPbl/pSrDDAU5iw/o18SxU8777DtN6Wq3vQ2d/8JFWtau4zq2wuWt1dkTzHrNufAaBNeHfhNsrXcC0Qk+8T50MbMmXmMFbzn19+jhI/oYcBeNY5zizJt4J94P7+cs1GvUmmX8Y96vv0+hncU=</latexit>
sha1_base64="5xqcIZ35fjydvHn6O9mq2gszF/k=">AAAC93icjVLLSsQwFD3W9/gadSO4KYrgQsZUBGcjiG5cKjjOgC/aGscwfZGmwjgIfoY70YW41b8R/AD1K7zJdMAHPlLanpx7z0nuTbwkEKli7KnL6u7p7esfGCwMDY+MjhXHJ3bTOJM+r/hxEMua56Y8EBGvKKECXkskd0Mv4FWvsaHj1TMuUxFHO6qZ8IPQrUfiRPiuIupwPxX10F1srLIFxpaPirOsxMywvwMnB7NrU+fPg5eP61tx8Q37OEYMHxlCcERQhAO4SOnZgwOGhLgDtIiThISJc1ygQNqMsjhluMQ26Fun2V7ORjTXnqlR+7RKQK8kpY050sSUJwnr1WwTz4yzZn/ybhlPvbcm/b3cKyRW4ZTYv3SdzP/qdC0KJyibGgTVlBhGV+fnLpnpit65/aEqRQ4JcRofU1wS9o2y02fbaFJTu+6ta+IvJlOzeu7nuRlef61O76DdazoXugbO10P/DnaXSg4rOdt0H8pojwFMYwbzdOorWMMmtlChVSSucINbq2ldW3fWfTvV6so1k/g0rId3zRafQg==</latexit>
sha1_base64="pshhDoDzYAH1ThLM3wSHKF+nvRs=">AAAC93icjVLLSgMxFD0dX7W+qi7dDBbBhdSMCLoRCm5cVrBVqA9mxlhD50WSEYr4H+5EF+LWz/EP1K/wJk7BBz4yzMzJufec5N4kyCKhNGNPJWdoeGR0rDxemZicmp6pzs61VZrLkLfCNErlQeArHomEt7TQET/IJPfjIOL7QW/bxPcvuFQiTfZ0P+NHsd9NxJkIfU3U8aES3dhf7W2xFcbWT6o1Vmd2uN+BV4AaitFMq684xClShMgRgyOBJhzBh6KnAw8MGXFHuCROEhI2znGFCmlzyuKU4RPbo2+XZp2CTWhuPJVVh7RKRK8kpYsl0qSUJwmb1Vwbz62zYX/yvrSeZm99+geFV0ysxjmxf+kGmf/VmVo0zrBpaxBUU2YZU11YuOS2K2bn7oeqNDlkxBl8SnFJOLTKQZ9dq1G2dtNb38afbaZhzTwscnO8/Fqd2cF7r+lc6Bp4Xw/9O2iv1T1W93ZZrbFZXIgyFrCIZTr1DTSwgyZatIrENW5x5/SdG+feeXhPdUqFZh6fhvP4BlkhnBc=</latexit>

<latexit sha1_base64="IXFE5t28m/IQk3BKhVIEqqvaw/w=">AAAC6nicjVLLSsNAFD2Nr/quunQTLIKrkrixy4Ibly3YB9QiSTqtQycPJhOhFL/AnehC3PpT/oH6Fd6ZpqAWHxOSnDn3njv3MX4ieKoc56VgLSwuLa8UV9fWNza3tks7u600zmTAmkEsYtnxvZQJHrGm4kqwTiKZF/qCtf3Rqba3r5lMeRydq3HCeqE3jPiAB54iqjG6LJWdimOWPQ/cHJSRr3pcescF+ogRIEMIhgiKsICHlJ4uXDhIiOthQpwkxI2d4QZrpM3Ii5GHR+yIvkPadXM2or2OmRp1QKcIeiUpbRySJiY/SVifZht7ZiJr9qfYExNT5zamv5/HColVuCL2L93M8786XYvCAFVTA6eaEsPo6oI8Sma6ojO3P1WlKEJCnMZ9skvCgVHO+mwbTWpq1731jP3VeGpW74PcN8Pbr9XpDKa9prnQNXC/D30etI4rrlNxG065Vs0vRBH7OMARTf0ENZyhjqbJ4Bb3eLCEdWc9Wk9TV6uQa/bwZVnPH72nl8w=</latexit>
<latexit sha1_base64="xkRW35A/XieE00d3e7jnwUiDqf4=">AAAC6nicjVK5TsNAEH0xVwhXgJLGIkKiimwa0hGJhjJB5JBChOzNBqysD9lrpCjiC+gQFIiWj6DiP/gD4CuY3TgSEHGMZfvtm3mzM7PrRsJLpGW95IyZ2bn5hfxiYWl5ZXWtuL7RTMI0ZrzBQhHGbddJuPAC3pCeFLwdxdzxXcFb7uBQ+VuXPE68MDiRw4h3fec88PoecyRR9cFZsWSVLW3mNLAzUDp4fjoGWS0svuMUPYRgSOGDI4AkLOAgoacDGxYi4roYERcT8rSf4woF0qYUxSnCIXZA33NadTI2oLXKmWg1o10EvTEpTeyQJqS4mLDazdT+VGdW7E+5Rzqnqm1IfzfL5RMrcUHsX7pJ5H91qheJPiq6B496ijSjumNZllRPRVVufupKUoaIOIV75I8JM62czNnUmkT3rmbraP+rjlSsWrMsNsXbr92pCsazpnOha2B/P/Rp0Nwr21bZrlulagVjy2ML29ilU99HFUeooaEruMYt7gxh3Bj3xsM41Mhlmk18MePxA4dgmdM=</latexit>
sha1_base64="Svt6pl9lIwlraG5BJAsRk5cwB2s=">AAAC6nicjVLLSsNAFD2Nr1pfVZeCBIvgqiRu7M6CG5ct2AfUIkk6raGTB8lEKMWlK3eiC3HrR3Tlf/gHPn7CO9MU1OJjQpIz595zc8+d2CF3Y2EYzxltZnZufiG7mFtaXlldy69v1OMgiRxWcwIeRE3bihl3fVYTruCsGUbM8mzOGnb/SMYbFyyK3cA/EYOQtT2r57td17EEUdX+Wb5gFA219GlgpqBw+DSqvl9tjypB/g2n6CCAgwQeGHwIwhwWYrpaMGEgJK6NIXERIVfFGS6RI21CWYwyLGL79OzRrpWyPu1lzVipHfoKpzsipY5d0gSUFxGWX9NVPFGVJftT7aGqKXsb0NtOa3nECpwT+5dukvlfnfQi0EVJeXDJU6gY6c5JqyRqKrJz/ZMrQRVC4iTuUDwi7CjlZM660sTKu5ytpeIvKlOycu+kuQlef3UnOxjPms6FfgPz+6FPg/p+0TSKZtUolEsYryy2sIM9OvUDlHGMCmqqg2vc4k7j2o12rz2MU7VMqtnEl6U9fgDbA5wY</latexit>

/k = 0, 00
Rec = 5772

<latexit
sha1_base64="ZP/SqFem4F760/5EeS1px5OMP5M=">AAAC93icjVJNb9NAEH11obSFtim9VOJiESFxQOm6F3KpVIkLx1QiaaS0INvZhFX8pfUaKUT9H70hOFRc+TmV+gNofwVvN44EVLSMZfvtm3mzOzMbFYkqjRCXS97yg4crj1bX1h8/2djcamw/7ZV5pWPZjfMk1/0oLGWiMtk1yiSyX2gZplEij6PJG+s//iR1qfLsnZkW8jQNx5kaqTg0pN6flGqchnuTA/FKiPaHRlO0hDP/Nghq0Dzc/XwFWidv3OAEQ+SIUSGFRAZDnCBEyWeAAAIFuVPMyGki5fwSZ1intmKUZERIdsLvmKtBzWZc25ylU8fcJeGrqfTxgpqccZrY7uY7f+UyW/ZfuWcupz3blP+ozpWSNfhI9j7dIvJ/dbYWgxHargbFmgrH2OriOkvlumJP7v9WlWGGgpzFQ/o1ceyUiz77TlO62m1vQ+f/6SIta9dxHVvh+s7q7AnmveZceA2Cv4d+G/T2W4FoBUe8D23MbRXP8BwvOfXXOMRbdNDlLhrn+Ipv3tT74l143+eh3lKt2cEf5v34BVm5nck=</latexit>
sha1_base64="K8MeMX3tgGdWyGO1KX4Vr/4+RV0=">AAAC93icjVLLSsNAFD2Nz/qsuhHcBEVwIXXixm4E0Y3LCrYK9UESxzo0LyYToYrgZ7gTXYhb/RvBD1C/wjvTFHzgY0KSM+fec2bunfGSQKSKsaeC1dPb1z8wWBwaHhkdGy9NTNbTOJM+r/lxEMtdz015ICJeU0IFfDeR3A29gO94rQ0d3znlMhVxtK3aCd8P3WYkjoXvKqIO9lLRDN2l1ipbZKxyWJpjZWaG/R04OZhbmz57Ll4+rlfj0hv2cIQYPjKE4IigCAdwkdLTgAOGhLh9nBMnCQkT57jAEGkzyuKU4RLbom+TZo2cjWiuPVOj9mmVgF5JShvzpIkpTxLWq9kmnhlnzf7kfW489d7a9Pdyr5BYhRNi/9J1M/+r07UoHKNiahBUU2IYXZ2fu2SmK3rn9oeqFDkkxGl8RHFJ2DfKbp9to0lN7bq3rom/mEzN6rmf52Z4/bU6vYNOr+lc6Bo4Xw/9O6gvlx1WdrboPlTQGYOYwSwW6NRXsIZNVFGjVSSucINbq21dW3fWfSfVKuSaKXwa1sM71y6fRg==</latexit>
sha1_base64="tfbwWQVnZIdjVW6BkW9iBkdhFVE=">AAAC93icjVLLSsNAFD2Nr1pfVZdugkVwIXXixm4EwY3LCvYBbZUkjnVoXkwmQhH/w53oQtz6Of6B+hXemabgAx8Tkpw5954zc++MlwQiVYw9F6yJyanpmeJsaW5+YXGpvLzSTONM+rzhx0Es256b8kBEvKGECng7kdwNvYC3vMGBjrcuuUxFHB2rYcJ7oduPxLnwXUXUSTcV/dDdHuyxLcZqp+UKqzIz7O/AyUEF+ajH5Td0cYYYPjKE4IigCAdwkdLTgQOGhLgeroiThISJc1yjRNqMsjhluMQO6NunWSdnI5prz9SofVoloFeS0sYGaWLKk4T1araJZ8ZZsz95XxlPvbch/b3cKyRW4YLYv3TjzP/qdC0K56iZGgTVlBhGV+fnLpnpit65/aEqRQ4JcRqfUVwS9o1y3GfbaFJTu+6ta+IvJlOzeu7nuRlef61O72DUazoXugbO10P/Dpo7VYdVnSNW2a/lF6KINaxjk059F/s4RB0NWkXiBne4t4bWrfVgPY5SrUKuWcWnYT29A2M5nBs=</latexit>

/k <
0
b) Re
sha1_base64="808KUZxxT9CMXt2H8OFYcXPvz0M=">AAAC8nicjVLLSsNAFD3GV61Wqy7dBEvBVUkK0m6EgiAuVawPtJRkOq2heTGZCKX4Fe5EFyK48iv8C/9A/QrvTCP4wMeEJGfOvefM3Dvjxr6XSMt6GjPGJyanpnMz+dm5wvxCcXHpIIlSwXiTRX4kjlwn4b4X8qb0pM+PYsGdwPX5odvfVPHDcy4SLwr35SDmrcDphV7XY44k6niPt9nGeq1WbRdLVsXSw/wO7AyUGoXHtLyVv9+Jiq84RQcRGFIE4AghCftwkNBzAhsWYuJaGBInCHk6znGBPGlTyuKU4RDbp2+PZicZG9JceSZazWgVn15BShNl0kSUJwir1UwdT7WzYn/yHmpPtbcB/d3MKyBW4ozYv3Tvmf/VqVokuqjrGjyqKdaMqo5lLqnuitq5+aEqSQ4xcQp3KC4IM61877OpNYmuXfXW0fFnnalYNWdZboqXX6tTOxj1ms6FroH99dC/g4NqxbYq9i7dhzpGI4cVrGKNTr2GBraxgyatEuAS17gxpHFl3Bp3o1RjLNMs49MwHt4AgQCc2w==</latexit>
sha1_base64="+y1N5itdRvXDZsk4wJcnKLjNQRs=">AAAC8nicjVLLSsNAFD2Nr/qoVl26CRbBVUkK0m6EghuXVexDainJdFpD8yKZCKX4Fe5EF+LWr/Av/AP1K7wzTUEtPiYkOXPuPWfm3hk7dJ1YGMZLRpubX1hcyi6vrK7l1jfym1uNOEgixusscIOoZVsxdx2f14UjXN4KI255tsub9vBIxptXPIqdwD8To5B3PGvgO32HWYKo81PeZYcH5XKpmy8YRUMNfRaYKShUc8+JDqAW5N9xgR4CMCTwwOFDEHZhIaanDRMGQuI6GBMXEXJUnOMaK6RNKItThkXskL4DmrVT1qe59IyVmtEqLr0RKXXskSagvIiwXE1X8UQ5S/Yn77HylHsb0d9OvTxiBS6J/Us3zfyvTtYi0EdF1eBQTaFiZHUsdUlUV+TO9U9VCXIIiZO4R/GIMFPKaZ91pYlV7bK3loq/qkzJyjlLcxO8/Vqd3MGk13QudA3M74c+CxqlomkUzRO6DxVMRhY72MU+nXoZVRyjhjqt4uEGd7jXhHarPWiPk1Qtk2q28WVoTx808Jvg</latexit>
<latexit sha1_base64="MxkfomHOWj+vI1G9OtRRMZa+QRo=">AAAC8nicjVLLSsNAFD2Nr1pfVZdugkVwVRJB2o1QcOOyin1IlZJMpzWYF5OJUIpf4U50IW79If9A/QrvTFNQi48JSc6ce8+ZuXfGjX0vkZb1kjNmZufmF/KLhaXlldW14vpGM4lSwXiDRX4k2q6TcN8LeU

<latexit sha1_base64="Jzo0uLdNXPwfMHfVRPUiswG97p8=">AAAC7HicjVLJTsMwFJyGrZStwJFLRIXEqUq40GMlLhzL0kUqFUpct1hNk8hxkKqKT+CG4IC48k/8AfAVPLupBFQsjpKM570Z+z3bjwORKMd5yVlz8wuLS/nlwsrq2vpGcXOrkUSpZLzOoiCSLd9LeCBCXldCBbwVS+4N/YA3/cGRjjevuUxEFJ6rUcw7Q68fip5gniLq7JTbl8WSU3bMsGeBm4ESslGLiu+4QBcRGFIMwRFCEQ7gIaGnDRcOYuI6GBMnCQkT57hBgbQpZXHK8Igd0LdPs3bGhjTXnolRM1oloFeS0sYeaSLKk4T1araJp8ZZsz95j42n3tuI/n7mNSRW4YrYv3TTzP/qdC0KPVRMDYJqig2jq2OZS2q6onduf6pKkUNMnMZdikvCzCinfbaNJjG16956Jv5qMjWr5yzLTfH2a3V6B5Ne07nQNXC/H/osaByUXafsnjilaiW7EHnsYBf7dOqHqOIYNdRplT5ucY8HK7TurEfraZJq5TLNNr4M6/kDGySYTA==</latexit>
sha1_base64="0AYvdRIbqKuW+sXQRsdXjAJQDcc=">AAAC7HicjVLLTsJAFD3UF+ILdemmkZi4Iq0b2UnixiWiPBIkpi0DTihtM52aEOInuDO6MG79Blf+h3+gfoV3hpKoxMdtCmfOPefO3Dt1I5/H0rJeMsbM7Nz8QnYxt7S8srqWX9+ox2EiPFbzQj8UTdeJmc8DVpNc+qwZCeYMXJ813P6hyjcumYh5GJzKYcTaA6cX8C73HEnUSZWZ5/mCVbR0mNPATkHh4PmpCopKmH/HGToI4SHBAAwBJGEfDmJ6WrBhISKujRFxghDXeYYr5MibkIqRwiG2T789WrVSNqC1qhlrt0e7+PQKcprYIU9IOkFY7WbqfKIrK/an2iNdU51tSP9uWmtArMQFsX/5Jsr/+lQvEl2UdA+ceoo0o7rz0iqJnoo6ufmpK0kVIuIU7lBeEPa0czJnU3ti3buaraPzr1qpWLX2Um2Ct1+7UycYz5ruhT4D+/ulT4P6XtG2ivaxVSiXMI4strCNXbr1fZRxhApqtEsP17jFnREYN8a98TCWGpnUs4kvYTx+AOTOmlM=</latexit>
sha1_base64="Qj2T/DIemAm0d4EjgOo1S8XFCfI=">AAAC7HicjVLLSsNAFD2Nr1pfVZeCBIvgqiRu7M6CG5et2gfUIkk6rUPTJCQToRSXLt2JLsSt39CV/+Ef+PgJ70xTUIuPCUnOnHvPmbl3xg5cHgnDeE5pU9Mzs3Pp+czC4tLySnZ1rRr5ceiwiuO7fli3rYi53GMVwYXL6kHIrJ7tsprdPZDx2gULI+57J6IfsGbP6ni8zR1LEHV8xPSzbM7IG2rok8BMQG7/aVh+v9oclvzsG07Rgg8HMXpg8CAIu7AQ0dOACQMBcU0MiAsJcRVnuESGtDFlMcqwiO3St0OzRsJ6NJeekVI7tIpLb0hKHduk8SkvJCxX01U8Vs6S/cl7oDzl3vr0txOvHrEC58T+pRtn/lcnaxFoo6Bq4FRToBhZnZO4xKorcuf6p6oEOQTESdyieEjYUcpxn3WliVTtsreWir+oTMnKuZPkxnj9tTq5g1Gv6VzoGpjfD30SVHfzppE3y0auWMBopLGBLezQqe+hiEOUUKFVOrjGLe40T7vR7rWHUaqWSjTr+DK0xw84gJyY</latexit>

Figure . – a) Taux de croissance σ en fonction du nombre d’onde k. b)

Iso-contours positifs de σ/k dans le plan (Re, k).

) Justifier la recherche de solutions de l’équation d’Orr-Sommerfeld sous la

forme ψ(x, z, t) = Ψ(z) exp(i k x + s t) où s = σ − i ω ∈ C
I.
) Montrer que Ψ(z) vérifie l’équation du quatrième ordre
hσ  ω i 2 1
2
+ i 1 − z2 − (D − k 2 )Ψ(z) − i a4 Ψ(z) = D2 − k 2 Ψ(z) ,
k k k Re

où D = d/dz et a4 est une constante que l’on précisera. Écrire les condi
tions aux limites vérifiées par Ψ(z) en z = ±1.
) Les taux de croissance σ = Σ(k, Re) obtenus en résolvant numériquement
cette équation avec ses conditions aux limites sont représentés sur la figure
.. Discuter la stabilité de l’écoulement de Poiseuille plan en fonction de
la valeur du nombre de Reynolds.
 Chapitre . Exercices complémentaires

) Au voisinage du seuil de l’instabilité, les modes instables sont appelés

“ondes de Tollmien-Schlichting”. Exprimer leur longueur d’onde lc en fonc-
tion de la longueur h.
) La déstabilisation de l’écoulement de Poiseuille est-elle une bifurcation
super-critique ou sous-critique ?

Corrigé . Stabilité de l’écoulement de Poiseuille

dU
)On a dte = 0 et le tracé de U (z) est un parabole. )On a Πe (x, z) = Πr −
2 ν U0 x/h2 − g z. )En omettant les ∗, on obtient div U = 0, ∂U ∂t + U ·grad U =
−grad Π − F r−2 ez + Re−1 ∆U . )Les conditions aux limites adimensionnées
s’écrivent u(x, 1, t) = w(x, 1, t) = 0 et u(x, −1, t) = w(x, −1, t) = 0. L’équilibre
s’écrit U e = U (z) ex avec U (z) = (1 − z 2 ) et Πe = 1 − 2 Re−1 x − F r−2 z.
)En remarquant que les termes d’ordre zéro s’annulent et en négligeant les
e U 0 (z) =
termes d’ordre deux, on obtient ∂∂xue + ∂∂zwe = 0, ∂∂tue + U (z) ∂∂xue + w
2 2 2 2
− ∂∂x
Π
− Re−1 ( ∂x
∂ ∂
e et ∂∂twe + U (z) ∂∂x
2 + ∂z 2 ) u
w
= − ∂∂zΠ − Re−1 ( ∂x
∂ ∂
2 + ∂z 2 ) w.
e e e
e
)L’incompressiblité justifie le choix d’une fonction de courant. )En dérivant
∂ ∂
par − ∂z l’équation pour u e, par ∂x l’équation pour w
e puis en les sommant, on
obtient l’équation indiquée pour a1 = a2 = a3 = 1. )On a ∂ψ ∂ψ
∂x = ∂z = 0
pour z = ±1. )On cherche des solutions sous cette forme dans la mesure òu le
système est invariant par translations en x et en t, mais pas en z. )Comme
U 00 = 2, on a a4 = 1. Les conditions aux limites sont Ψ = DΨ = 0 en z = ±1.
)L’écoulement est stable jusqu’au seuil Rec = 5772. )Comme kc = 1, 02,
on a lc = 2 π h/kc ∼ 6, 16 h. )Comme pour tous les écoulements parallèles,
la bifurcation est sous-critique.

 Exercice sur le chapitre 

EXERCICE . Rouleaux et relation de dispersion

On considère un champ de vitesse D défini par la fonction de courant ψ(x, y, t)

à l’aide des relations u = − ∂ψ ∂ψ
∂z et w = ∂x . On suppose que ψ est solution
d’une équation aux dérivées partielles linéarisée autour d’un équilibre dont les
coefficients sont réels et constants en espace et en temps.
) En l’absence de conditions aux limites, c’est-à-dire dans un milieu in-
fini, justifier la recherche de solutions bornées en espace de la forme
ψ = Ψm ei kx x ei kz z es t où Ψm est une amplitude complexe.
) On suppose que p la relation de dispersion généralisée est de la forme s =
Σ(k) où k = kx2 + kz2 est le module du vecteur d’onde k = kx ex + kz ez
et Σ(k) une fonction réelle. Quelles conditions doit vérifier la fonction Σ
pour que l’équilibre soit stable ?
) On suppose maintenant que ψ doit vérifier les conditions aux limites
ψ(x, 0) = 0 et ψ(x, d) = 0 pour tout x ∈ R. Montrer que l’on peut trouver
une solution réelle du problème avec conditions aux limites de la forme

ψ(x, z, t) = A exp(s t) cos(kx x) sin(n π z/d)

Exercice sur le chapitre  

où A ∈ IR et n est un entier quelconque avec s = Σ(k). Détailler le

raisonnement permettant d’affirmer que ψ est bien une solution.
) Tracer les lignes de champ du champ de vitesse ainsi obtenues dans le cas
n = 1.
) Tracer, pour ce même cas, les vecteurs vitesses sur le même graphe en
supposant A > 0.

Corrigé . Rouleaux et relation de dispersion

)Commme le système d’équations aux dérivées partielles est invariant par les
translations d’espace et de temps, on sait que les modes propres sont de la
forme indiquée. Comme les coefficients sont réels et que le système est linéaire,
la partie réelle d’une solution complexe est solutions.
)Pour que l’équilibre étudié soit stable, il faut que la partie réelle de Σ soit
strictement négative pour tous les modes, c’est-à-dire pour toutes les valeurs
de k.
)Les fonctions φ+ = 21i A ei kx x ei π z eS(k) t et φ− = − 21i A ei kx x e−i π z eS(k) t
sont des solutions complexes élémentaires dans la mesure où les vecteurs
d’ondes k + = kx ex + kz ez et k − = kx ex − kz ez ont le même module. Leur
somme est égale à la fonction ψ = A eS(k) t cos(kx x) sin(π z) qui est donc
solution du problème linéaire. On peut donc choisir A réel et A(t) = A eS(k) t .
)Les isolignes sont celles des rouleaux de convection : des cercles au centre,
des carrés aux frontières et une transformation de cercles en carrés au milieux
(voir figure ??). La période en x est Lx = 2 π/kx . La hauteur en z est de .
Les maxima et minima de ψ alternent tous les Lx /2 dans la direction x.
)En supposant A > 0, le champ de vitesse U = −ey ∧ grad ψ tourne dans le
sens trigonométrique autour des maxima de ψ et dans le sens horaire autour
de ses minima.

PROBLÈME . Milieu poreux en rotation

On considère un fluide pesant s’écoulant dans un milieu poreux d’extension

infini (éloigné de toute frontière). On note g = −g ez le champ de gravité.
L’ensemble du dispositif est soumis à une rotation uniforme d’axe vertical
Ω0 = Ω0 ez . On suppose que la température du fluide est constante et que le
mouvement de ce fluide, dans le repère tournant, est régi par le modèle :

div U = 0
dU 1
+ 2 Ω0 ∧ U = − grad p − g ez − β U , (.)
dt ρr

où β est une constante positive rendant compte du frottement du milieu poreux

d ∂
(loi de Darcy). On note U (x, y, z, t) = u ex + v ey + w ez et dt = ∂t + U · grad
(dérivée particulaire).
) Montrer le repos U = 0 est une solution du problème. Expliciter, pour cet
équilibre, le champ de pression pc (z) associé en supposant que pc (0) = pr .
 Chapitre . Exercices complémentaires

) On pose p(x, t) = pc (z) + ρr Π(x, t). Montrer que les équations du modèle
∂
linéarisées autour du repos peuvent s’écrire div U = 0, ∂t U + 2 Ω0 ez ∧ U =
−grad Π − β U .
) On choisit d’adimensionnaliser les équations avec les unités [L] = g/(4 Ω20 ),
[τ ] = 1/(2 Ω0 ), [U ] = [L]/[τ ] et [Π] = [U ]2 . Écrire les équations du mo-
dèle linéarisé ainsi adimensionnées en mettant en évidence le nombre sans
dimension E = β/(2 Ω0 ).
) On note désormais les grandeurs adimensionnées avec les mêmes notations
que celles des grandeurs physiques. On s’intéresse aux petits mouvements
quasi-D définis par

U = u(x, z, t) ex + v(x, z, t) ey + w(x, z, t) ez Π = Π(x, z, t) .

et
(.)
Écrire les équations adimensionnées pour les composantes u, v et w.
) Justifier, sans démonstration, l’existence d’une fonction de courant
ψ(x, z, t) telle que

∂ψ ∂ψ
u=− et w= . (.)
∂z ∂x
∂ ∂
) En déduire que le système peut s’écrire sous la forme ∂t ∆ψ = −A1 ∂z v−
∂ ∂
A2 ∆ψ et ∂t v = A3 ∂z ψ − A4 v où A1 , A2 , A3 et A4 sont des constantes
que l’on déterminera.
) On cherche des solutions complexes de la forme

[ψ(x, z, t), v(x, z, t)] = [Ψm , Vm ]ei kx x+i kz z+s t , (.)

où Ψm et Vm sont des constantes complexes. Justifier la recherche de solu-

tions sous cette forme.
) En déduire la relation de dispersion généralisée des perturbations de l’état
de repos. Montrer que l’on peut écrire s = −A5 ± i Ω(kx , kz ) où A5 est un
constante que l’on précisera.
p note θ(kx , kz ) l’angle défini par kx = k cos θ, kz = k sin θ et k =
) On
kx2 + kz2 . Exprimer Ω(kx , kz ) en fonction de θ(kx , kz ).
) Le repos est-il stable ou instable vis-à-vis des petites perturbations ?
) Tracer très schématiquement une composante de vitesse en fonction du
temps pour un point quelconque dans le cas kx = kz et E = 1.

Corrigé . Milieu poreux en rotation

)On a pc (z) = p1 − ρr g z. )En remplaçant dans les équations, les termes

−grad pc /ρr − g ez = 0 disparaissent. Il suffit ensuite de négliger le terme
U · grad U . )En utilisant les mêmes notations que pour les grandeurs avec
dimensions, les équations adimensionnalisées s’écrivent div U = 0, ∂U
∂t + ez ∧
U = −grad Π − E U . )Dans le cas quasi-D, le modèle s’écrit ∂u∂x + ∂w
∂z = 0,
∂u ∂Π ∂v ∂Π ∂w ∂Π
∂t − v = − ∂x − E u, ∂t + u = − ∂y − E v, ∂t = − ∂z − E w. )Les champs
de vitesse D incompressibles s’écrivent à l’aide d’une fonction de courant.
)En multipliant l’équation en u = − ∂ψ ∂ ∂ψ
∂z par − ∂z et l’équation en w = ∂x par
∂
∂x puis en les sommant, on obtient A1 = 1 et A2 = E. De plus, A3 = 1 et
Exercice sur le chapitre  

A4 = E. )Comme les équations sont linéaires à coefficients réelles, on peut

chercher des solutions complexes et prendre ensuite leur partie réelle. La forme
des solutions recherchées est justifiées par l’invariance du problème par les
translations en temps et en espace, dans toutes les directions. )En reportant
dans les équations, on obtient (s + E) (−k 2 ) Ψm = −i kz Vm et (s + E) Vm =
i kz Ψm avec k 2 = kx2 + kz2 . On en déduit la relation de dispersion généralisée
(s +pE)2 = −kz2 /k 2 et donc s = −E ± i Ω(kx , kz ) avec Ω(kx , kz ) = kz /k =
kz / kx2 + kz2 . )On a Ω(kx , kz ) = sin θ(kx , kz ). )Le repos est toujours stable
car Re(s) = −E < √ 0. )La vitesse en un point est proportionnelle à la fonction
−t
f (t) = e cos(t/ 2+ϕ) on ϕ est une phase. Il s’agit d’une sinusoïde qui décroit
fortement au bout d’une période.

EXERCICE . Convection thermohaline

La convection thermohaline est le mouvement induit par les changements de

masse volumique dus aux variations de température et de concentration de sel
dans l’eau. On décrit ici ces mouvements de convection par le modèle
dU 1 ρ
div U = 0 , = − grad p − g e + ν ∆U ,
dt ρr ρr z
ρ dT dS
= 1 − α (T − Tr ) + β (S − Sr ) , = κT ∆T , = κS ∆S , (.)
ρr dt dt
où U est le champ de vitesse, T le champ de température, S le champ de salinité
et p le champ de pression. Les coefficients de diffusivité de la température et
de la salinité sont respectivement κT et κS . Les coefficients de dilatation dus à
la température et à la salinité sont respectivement α > 0 et β > 0. On a noté
ρr , Tr et Sr des constantes de référence. La gravité g est dans la direction Oz.
On s’intéresse aux rouleaux de convection qui se développent dans une couche
fluide comprise entre deux plaques planes d’équations respectives z = 0 et
z = d où d est l’épaisseur de la couche. Le forçage du mouvement est induit
par les conditions aux limites (T, S) = (T1 , S1 ) en z = 0 et (T, S) = (T2 , S2 ) en
z = d avec T1 > T2 et S1 > S2 . On suppose que les champs sont indépendants
de y et que la vitesse v dans cette direction est nulle. On suppose des conditions
aux limites sans glissement pour la vitesse, c’est-à-dire w = 0 et ∂u ∂z = 0

) Calculer les profils conductifs Tc (z) et Sc (z) obtenus dans le cas U = 0.

En déduire
T1 − T2 2
pc = pr − ρr g z + ρr α g (T1 − Tr ) z − ρr α g z
2d
S1 − S2 2
− ρr β g (S1 − Sr ) z + ρr β g z . (.)
2d

On a Tc = T1 − T1 −T
d
2
z et Sc = S1 − S1 −S
d
2
z. Comme ρc = 1−α(Tc −Tr )+β(Sc −Sr )
et pc (z) = −ρr g ρc (z), on en déduit l’expression de pc (z) indiquée.
0

) On pose T = Tc (z) + θ(x, z, t) et S = Sc (z) + Σ(x, z, t) et p = pc (z) +

ρr Π(x, t). Exprimer, à l’aide de ces nouvelles variables, les équations du
modèle avec ses conditions aux limites.

Les équations du modèle se mettent sous la forme div U = 0, dU dt = −grad Π +

T1 −T2 S1 −S2
α g θ ez − β g Σ ez + ν ∆U , dθ
dt = d w + κT ∆θ, dΣ
dt = d w + κS ∆Σ avec les
conditions aux limites θ = 0, Σ = 0, w = 0 et dudz = 0 en z = 0 et z = d.
 Chapitre . Exercices complémentaires

) On choisit d comme unité de longueur, d2 /κT comme unité de temps,

T1 − T2 comme unité de température et S1 − S2 comme unité de salinité.
Écrire les équations adimensionnées à l’aide des nombres sans dimension

α g d3 (T1 − T2 ) β g d3 (S1 − S2 )
RT = , RS = ,
ν κT ν κS
ν κT
P = et L=
κT κS
respectivement appelés “Rayleigh thermique”, “Rayleigh salin”, “Prandtl”
et “Lewis”.

On pose x = d x∗ , z = d z∗ , t = (d2 /κT ) t∗ , u = (κT /d) u∗ , w = (κT /d) w∗ , θ =

(T1 − T2 ) θ∗ , Σ = (S1 − S2 ) Σ∗ et Π = (κ2T /d2 ) Π∗ . En remplaçant dans les équations
du modèle puis en omettant les indices ∗, on obtient div U = 0, dU dt = −grad Π +
RT P θ ez − RS P L Σ ez + P ∆U , dt = w + ∆θ, et dt = w + L ∆Σ.
−1 dθ dΣ −1

) On s’intéresse aux petites perturbations de l’équilibre. Montrer que l’on

peut mettre le système sous la forme
∂ ∂θ ∂Σ
∆ψ = RT P − RS P L−1 + P ∆2 ψ
∂t ∂x ∂x
∂θ ∂ψ ∂Σ ∂ψ
= + ∆θ , = + L−1 ∆Σ .
∂t ∂x ∂t ∂x

On peut tout d’abord négliger les termes non linéaires en U · grad . On définit une
fonction courant ψ par les relations u = − ∂ψ ∂ψ
∂z et w = ∂x . En prenant le rotationnel
de l’équation de quantité de mouvement, on obtient l’équation pour ∆ψ indiquée.
) On cherche des solutions complexes du modèle linéarisé sous la forme

(ψ, θ, Σ) = (ψm , θm , Σm ) ei(kx x+kz z)+s t (.)

où (ψm , θm , Σm ) sont trois amplitudes complexes, k = kx ex + kz ez un

vecteur d’onde réel et s = σ − i ω une valeur propre complexe. Montrer
que s doit être la racine de l’équation

s k 2 P −1 = A kx2 (s + k 2 )−1 − B kx2 (s L + k 2 )−1 − k 4 , (.)

où A et B sont des constantes que l’on explicitera. En déduire que les

valeurs propres s sont les racines d’un polynôme d’ordre  (qu’il n’est pas
nécessaire de développer).

Le système linéaire peut se mettre sous la forme −s k 2 ψm = RT P i kx θm −

RS P L−1 i kx Σm + P k 4 ψm , s θm = i kx ψm − k 2 θm = 0 et s Σm = i kx ψm −
L−1 k 2 Σm = 0. En annulant le déterminant de ce système, on obtient s k 2 P −1 =
RT kx2 (s + k 2 )−1 − RS kx2 (s L + k 2 )−1 − k 4 . Donc A = RT et B = RS .
) Montrer que s = 0 est une racine de ce polynôme lorsque RT − RS =
F (kx , kz ) où F est une fonction que l’on explicitera. Dans le cas où les
deux autres racines ont des parties réelles strictement négatives, discuter
la stabilité de l’état conductif lorsque le paramètre R = RT −RS est proche
de F (kx , kz ).

En supposant que s = 0 pour le polynôme d’ordre  qui exprime la relation

de dispersion généralisée, on obtient la relation RT − RS = k 6 /kx2 . On a donc
F (kx , kz ) = (kx2 + kz2 )3 /kx2 . Si RT − RS > F (kx , kz ), la valeur propre s = 0 devient
positive et l’équilibre devient instable. Si RT −RS < F (kx , kz ), l’équilibre est stable.
Exercice sur le chapitre  

) Montrer que les choix kz = n π et kz = n π avec n entier permettent de

construire des solutions satisfaisant les conditions aux limites. En déduire
la taille des rouleaux de convection lorsque l’état conductif bifurque vers
une solution stationnaire.

La superposition de deux solutions ayant le même kx , des kz opposés avec kz = n π

et des amplitudes opposées, est une nouvelle solution. Elle satisfait les√conditions
aux limites. Le minimum de (kx2 + n2 π 2 )3 /kx2 est obtenu pour kc = π/ 2 et vaut
Rc = 27 π 4 /4. Le seuil d’instabilité est obtenu pour RT − RS = Rc et k = kc . La
taille des rouleaux
√ de convection au seuil de cette instabilité stationnaire est alors
Lc = 2 π/kc = 2 2.

 Exercice sur le chapitre 

EXERCICE . Bifurcation fourche imparfaite

On considère le système dynamique ẋ = µ x − x3 +  avec x(t) ∈ IR et µ ∈ IR

et  ∈ IR+ .
) Définir l’espace des phases et l’espace de contrôle de ce système. Dans
le cas particulier  = 0, calculer les équilibres et tracer le diagramme de
bifurcation en fonction de µ en dessinant les trajectoires. Comme se nomme
la bifurcation obtenue ?

L’espace des phases des états x est IR et l’espace des paramètres de contrôle (µ, )
est IR × IR+ . Le diagramme de bifurcation pour  = 0 est celui d’une bifurcation
fourche supercritique.
) Sans faire de calcul, mais en traçant éventuellement l’allure de fonctions
utiles, tracer le diagramme de bifurcation en fonction de µ dans le cas
0 <   1. Quelle est ou quelles sont les bifurcations obtenues dans ce cas.

Le tracé des fonctions f (µ, ; x) = µ x − x3 +  pour différentes valeurs de µ permet

de déterminer le nombre d’équilibres et leur stabilité et de tracer le diagramme de
bifurcation (figure .). La bifurcation fourche est remplacée par une bifurcation
noeud-col à cause de la brisure de symétrie induite par le (petit) terme .

f f f x
! x ! x ! x
0 0 µ
0 0
µ<0 µ∼0 µ>0

Figure . – Bifurcation fourche imparfaite

 Chapitre . Exercices complémentaires

EXERCICE . Bifurcations avec hysteresis

On considère le système dynamique réel Ẋ = µ + 3 X − X 3 où X(t) ∈ IR et

µ ∈ IR. On cherche à déterminer et interpréter son diagramme de bifurcation.
) Tracer la fonction G(X) = 3 X − X 3 en spécifiant les coordonnées de ses
extrema.
) Tracer la fonction F (X) = µ + G(X) pour toutes les valeurs significatives
de µ. Indiquer à chaque fois son nombre de racines.
) Tracer dans le plan (µ, X) les équilibres Xe (µ) du système dynamique.
) Déterminer la stabilité de chacun de ces équilibres et tracer en traits pleins
les équilibres stables et en traits pointillés les équilibres instables.
) Tracer, dans le plan (µ, X), les trajectoires du système dynamique.
) Décrire le phénomène d’hysteresis observé en augmentant µ d’une valeur
fortement négative vers une valeur fortement positive puis en le diminuant
de nouveau.

Corrigé . Bifurcations avec hysteresis

F X
2
µ
0 2
0 1
X

a) b)

Figure . – a) Fonction F (X) pour plusieurs valeurs de µ. b) Diagramme

de bifurcation.

)La fonction G(X) est tracée en traits gras sur la figure .a. Ses extrema
sont atteints en X = ±1 et valent respectivement ±2.
)Les fonctions F (X) sont tracées sur la figure .a. Elles ont une racine pour
|µ| > 2, deux (dont une double) pour |µ| = 2 et trois pour |µ| < 2.
)Le tracé des équilibres est la courbe µ = 3 X − X 3 dans le plan (µ, X)
(figure .b).
)La stabilité des équilibres Xe s’obtient en lisant sur la figure .a le signe de
la dérivée F 0 (Xe ). On obtient le graphe de la figure .b.
)Les trajectoires sont des points, des segments de droites ou des demi-droites
représentées sur la figure .b.
)En augmentant µ, on suit l’équilibre stable positif jusqu’à µ = 2 où l’on
transite brusquement vers l’équilibre stable négatif. Si l’on diminue alors µ, on
reste sur cet équilibre jusqu’à µ = −2 où l’on transite de nouveau.
Exercice sur le chapitre  

EXERCICE . Polynôme de degré trois

On considère le système dynamique réel défini par Ẋ = F (X) avec X(t) ∈ IR

, µ ∈ IR et F (X) = µ X − X 2 + X 3 .
) Déterminer les équilibres Xe en fonction des valeurs de µ.
) Tracer l’ensemble des ces équilibres dans un diagramme (µ, X).
) Indiquer la valeur de µ∗ au-delà de laquelle il n’existe plus qu’un seul
équilibre.
) Tracer F (X) pour les cinq cas suivants :
µ ∈] − ∞, 0[, µ = 0, µ ∈]0, µ∗ [, µ = µ∗ et µ ∈]µ∗ , ∞[.
) Discuter la stabilité de tous les équilibres.
) Tracer le diagramme de bifurcation dans le plan (µ, X) en représentant les
trajectoires à l’aide de flêches.

Corrigé . Polynôme de degré trois

)L’équilibre X = 0 existe pour toutes les valeurs de µ. Pour µ < µ∗ = 1/4, il

√
existe deux équilibres supplémentaires Xe± = (1± 1 − 4 µ)/2, qui ne sont plus
que Xe∗ = 1/2 pour µ = µ∗ . )Le tracé des équilibres est constitué de la droite
X = 0 et de la parabole µ = X −X 2 de sommet (µ, X) = (µ∗ , Xe∗ ) = (1/4, 1/2).
)Pour µ > µ∗ = 1/4, il n’y a plus qu’un équilibre. )Le tracé de Y = F (X)
est représenté sur la figure . pour les cinq cas considérés. )Un équilibre
Xe est stable si sa dérivée F 0 (Xe ) est négative. Comme F 0 (0) = µ, l’équilibre
Xe = 0 est stable pour µ < 0, marginal pour µ = 0 et instable pour µ > 0.
Comme F 0 (Xe ) = Xe (2 Xe − 1) pour Xe = Xe± , on en déduit que le seul cas
de stabilité pour ces équibres est pour Xe = Xe− et µ ≥ 0. )Le diagramme de
bifurcation avec ses trajectoires est tracé sur la figure ..

a)
<latexit sha1_base64="XXZfoyKNPcbam78K3eaf8VJ25ag=">AAAC63icjVLLSsNAFD2Nr7a+qm4EN0ERdFPSurDLohuXVewDapFkOq2heZFMhFoEP8Cd6ELc+gP+jeAHqF/hnWkKavExIcmZc+85M/fOWIFjR8IwnlPaxOTU9Ew6k52dm19YzC0t1yI/DhmvMt/xw4ZlRtyxPV4VtnB4Iwi56VoOr1u9fRmvn/Mwsn3vWPQD3nLNrmd3bGYKoo7M7dPchpE31NDHQSEBG+XVi5fM1dNexc+94wRt+GCI4YLDgyDswERETxMFGAiIa2FAXEjIVnGOS2RJG1MWpwyT2B59uzRrJqxHc+kZKTWjVRx6Q1Lq2CSNT3khYbmaruKxcpbsT94D5Sn31qe/lXi5xAqcEfuXbpT5X52sRaCDkqrBppoCxcjqWOISq67IneufqhLkEBAncZviIWGmlKM+60oTqdplb00Vf1WZkpVzluTGePu1OrmDYa/pXOgaFL4f+jioFfOFnXzxkO5DCcORxhrWsUWnvosyDlBBlVbp4Bq3uNNc7Ua71x6GqVoq0azgy9AePwCnepsk</latexit>
F
<latexit sha1_base64="BJ5hOQW/IxQg9ADJPonmMpyuL6I=">AAAC6nicjVLLSsNAFD2Nr1pfVZeCBIvgqiS6sDsLgrhswT6gFknSaR2aJiGZCKW4dOVOdCFu/Yiu/A//wMdPeGeaglp8TEhy5tx7bu65EztweSQM4zmlTU3PzM6l5zMLi0vLK9nVtWrkx6HDKo7v+mHdtiLmco9VBBcuqwchs3q2y2p291DGaxcsjLjvnYh+wJo9q+PxNncsQVT56CybM/KGWvokMBOQO3galt+vNoclP/uGU7Tgw0GMHhg8CMIuLER0NWDCQEBcEwPiQkJcxRkukSFtTFmMMixiu/Ts0K6RsB7tZc1IqR36ikt3SEod26TxKS8kLL+mq3isKkv2p9oDVVP21qe3ndTqEStwTuxfunHmf3XSi0AbBeWBk6dAMdKdk1SJ1VRk5/onV4IqBMRJ3KJ4SNhRyvGcdaWJlHc5W0vFX1SmZOXeSXJjvP7qTnYwmjWdC/0G5vdDnwTV3by5l98tG7liAaOVxga2sEOnvo8ijlFCRXVwjVvcaa52o91rD6NULZVo1vFlaY8ffvOb9w==</latexit>
F
<latexit sha1_base64="BJ5hOQW/IxQg9ADJPonmMpyuL6I=">AAAC6nicjVLLSsNAFD2Nr1pfVZeCBIvgqiS6sDsLgrhswT6gFknSaR2aJiGZCKW4dOVOdCFu/Yiu/A//wMdPeGeaglp8TEhy5tx7bu65EztweSQM4zmlTU3PzM6l5zMLi0vLK9nVtWrkx6HDKo7v+mHdtiLmco9VBBcuqwchs3q2y2p291DGaxcsjLjvnYh+wJo9q+PxNncsQVT56CybM/KGWvokMBOQO3galt+vNoclP/uGU7Tgw0GMHhg8CMIuLER0NWDCQEBcEwPiQkJcxRkukSFtTFmMMixiu/Ts0K6RsB7tZc1IqR36ikt3SEod26TxKS8kLL+mq3isKkv2p9oDVVP21qe3ndTqEStwTuxfunHmf3XSi0AbBeWBk6dAMdKdk1SJ1VRk5/onV4IqBMRJ3KJ4SNhRyvGcdaWJlHc5W0vFX1SmZOXeSXJjvP7qTnYwmjWdC/0G5vdDnwTV3by5l98tG7liAaOVxga2sEOnvo8ijlFCRXVwjVvcaa52o91rD6NULZVo1vFlaY8ffvOb9w==</latexit>
F
<latexit sha1_base64="BJ5hOQW/IxQg9ADJPonmMpyuL6I=">AAAC6nicjVLLSsNAFD2Nr1pfVZeCBIvgqiS6sDsLgrhswT6gFknSaR2aJiGZCKW4dOVOdCFu/Yiu/A//wMdPeGeaglp8TEhy5tx7bu65EztweSQM4zmlTU3PzM6l5zMLi0vLK9nVtWrkx6HDKo7v+mHdtiLmco9VBBcuqwchs3q2y2p291DGaxcsjLjvnYh+wJo9q+PxNncsQVT56CybM/KGWvokMBOQO3galt+vNoclP/uGU7Tgw0GMHhg8CMIuLER0NWDCQEBcEwPiQkJcxRkukSFtTFmMMixiu/Ts0K6RsB7tZc1IqR36ikt3SEod26TxKS8kLL+mq3isKkv2p9oDVVP21qe3ndTqEStwTuxfunHmf3XSi0AbBeWBk6dAMdKdk1SJ1VRk5/onV4IqBMRJ3KJ4SNhRyvGcdaWJlHc5W0vFX1SmZOXeSXJjvP7qTnYwmjWdC/0G5vdDnwTV3by5l98tG7liAaOVxga2sEOnvo8ijlFCRXVwjVvcaa52o91rD6NULZVo1vFlaY8ffvOb9w==</latexit>
F
<latexit sha1_base64="BJ5hOQW/IxQg9ADJPonmMpyuL6I=">AAAC6nicjVLLSsNAFD2Nr1pfVZeCBIvgqiS6sDsLgrhswT6gFknSaR2aJiGZCKW4dOVOdCFu/Yiu/A//wMdPeGeaglp8TEhy5tx7bu65EztweSQM4zmlTU3PzM6l5zMLi0vLK9nVtWrkx6HDKo7v+mHdtiLmco9VBBcuqwchs3q2y2p291DGaxcsjLjvnYh+wJo9q+PxNncsQVT56CybM/KGWvokMBOQO3galt+vNoclP/uGU7Tgw0GMHhg8CMIuLER0NWDCQEBcEwPiQkJcxRkukSFtTFmMMixiu/Ts0K6RsB7tZc1IqR36ikt3SEod26TxKS8kLL+mq3isKkv2p9oDVVP21qe3ndTqEStwTuxfunHmf3XSi0AbBeWBk6dAMdKdk1SJ1VRk5/onV4IqBMRJ3KJ4SNhRyvGcdaWJlHc5W0vFX1SmZOXeSXJjvP7qTnYwmjWdC/0G5vdDnwTV3by5l98tG7liAaOVxga2sEOnvo8ijlFCRXVwjVvcaa52o91rD6NULZVo1vFlaY8ffvOb9w==</latexit>
F
<latexit sha1_base64="BJ5hOQW/IxQg9ADJPonmMpyuL6I=">AAAC6nicjVLLSsNAFD2Nr1pfVZeCBIvgqiS6sDsLgrhswT6gFknSaR2aJiGZCKW4dOVOdCFu/Yiu/A//wMdPeGeaglp8TEhy5tx7bu65EztweSQM4zmlTU3PzM6l5zMLi0vLK9nVtWrkx6HDKo7v+mHdtiLmco9VBBcuqwchs3q2y2p291DGaxcsjLjvnYh+wJo9q+PxNncsQVT56CybM/KGWvokMBOQO3galt+vNoclP/uGU7Tgw0GMHhg8CMIuLER0NWDCQEBcEwPiQkJcxRkukSFtTFmMMixiu/Ts0K6RsB7tZc1IqR36ikt3SEod26TxKS8kLL+mq3isKkv2p9oDVVP21qe3ndTqEStwTuxfunHmf3XSi0AbBeWBk6dAMdKdk1SJ1VRk5/onV4IqBMRJ3KJ4SNhRyvGcdaWJlHc5W0vFX1SmZOXeSXJjvP7qTnYwmjWdC/0G5vdDnwTV3by5l98tG7liAaOVxga2sEOnvo8ijlFCRXVwjVvcaa52o91rD6NULZVo1vFlaY8ffvOb9w==</latexit>

X
<latexit sha1_base64="p1j1zoOBRDWaYar4rxVus3IOhfg=">AAAC6nicjVLLSsNAFD2Nr1pfVZeCBIvgqqR1YXcW3LhswT6gFknSaQ3Ni8lEKMWlK3eiC3HrR3Tlf/gHPn7CO9MU1OJjQpIz595zc8+dWKHrRMIwnlPazOzc/EJ6MbO0vLK6ll3fqEdBzG1WswM34E3LjJjr+KwmHOGyZsiZ6Vkua1j9IxlvXDAeOYF/IgYha3tmz3e6jm0KoqrNs2zOyBtq6dOgkIDc4dOo+n61PaoE2TecooMANmJ4YPAhCLswEdHVQgEGQuLaGBLHCTkqznCJDGljymKUYRLbp2ePdq2E9Wkva0ZKbdNXXLo5KXXskiagPE5Yfk1X8VhVluxPtYeqpuxtQG8rqeURK3BO7F+6SeZ/ddKLQBcl5cEhT6FipDs7qRKrqcjO9U+uBFUIiZO4Q3FO2FbKyZx1pYmUdzlbU8VfVKZk5d5OcmO8/upOdjCeNZ0L/QaF74c+DerFfGE/X6wauXIJ45XGFnawR6d+gDKOUUFNdXCNW9xprnaj3WsP41QtlWg28WVpjx+sX5wJ</latexit>
X
<latexit sha1_base64="p1j1zoOBRDWaYar4rxVus3IOhfg=">AAAC6nicjVLLSsNAFD2Nr1pfVZeCBIvgqqR1YXcW3LhswT6gFknSaQ3Ni8lEKMWlK3eiC3HrR3Tlf/gHPn7CO9MU1OJjQpIz595zc8+dWKHrRMIwnlPazOzc/EJ6MbO0vLK6ll3fqEdBzG1WswM34E3LjJjr+KwmHOGyZsiZ6Vkua1j9IxlvXDAeOYF/IgYha3tmz3e6jm0KoqrNs2zOyBtq6dOgkIDc4dOo+n61PaoE2TecooMANmJ4YPAhCLswEdHVQgEGQuLaGBLHCTkqznCJDGljymKUYRLbp2ePdq2E9Wkva0ZKbdNXXLo5KXXskiagPE5Yfk1X8VhVluxPtYeqpuxtQG8rqeURK3BO7F+6SeZ/ddKLQBcl5cEhT6FipDs7qRKrqcjO9U+uBFUIiZO4Q3FO2FbKyZx1pYmUdzlbU8VfVKZk5d5OcmO8/upOdjCeNZ0L/QaF74c+DerFfGE/X6wauXIJ45XGFnawR6d+gDKOUUFNdXCNW9xprnaj3WsP41QtlWg28WVpjx+sX5wJ</latexit>
X
<latexit sha1_base64="p1j1zoOBRDWaYar4rxVus3IOhfg=">AAAC6nicjVLLSsNAFD2Nr1pfVZeCBIvgqqR1YXcW3LhswT6gFknSaQ3Ni8lEKMWlK3eiC3HrR3Tlf/gHPn7CO9MU1OJjQpIz595zc8+dWKHrRMIwnlPazOzc/EJ6MbO0vLK6ll3fqEdBzG1WswM34E3LjJjr+KwmHOGyZsiZ6Vkua1j9IxlvXDAeOYF/IgYha3tmz3e6jm0KoqrNs2zOyBtq6dOgkIDc4dOo+n61PaoE2TecooMANmJ4YPAhCLswEdHVQgEGQuLaGBLHCTkqznCJDGljymKUYRLbp2ePdq2E9Wkva0ZKbdNXXLo5KXXskiagPE5Yfk1X8VhVluxPtYeqpuxtQG8rqeURK3BO7F+6SeZ/ddKLQBcl5cEhT6FipDs7qRKrqcjO9U+uBFUIiZO4Q3FO2FbKyZx1pYmUdzlbU8VfVKZk5d5OcmO8/upOdjCeNZ0L/QaF74c+DerFfGE/X6wauXIJ45XGFnawR6d+gDKOUUFNdXCNW9xprnaj3WsP41QtlWg28WVpjx+sX5wJ</latexit>
X
<latexit sha1_base64="p1j1zoOBRDWaYar4rxVus3IOhfg=">AAAC6nicjVLLSsNAFD2Nr1pfVZeCBIvgqqR1YXcW3LhswT6gFknSaQ3Ni8lEKMWlK3eiC3HrR3Tlf/gHPn7CO9MU1OJjQpIz595zc8+dWKHrRMIwnlPazOzc/EJ6MbO0vLK6ll3fqEdBzG1WswM34E3LjJjr+KwmHOGyZsiZ6Vkua1j9IxlvXDAeOYF/IgYha3tmz3e6jm0KoqrNs2zOyBtq6dOgkIDc4dOo+n61PaoE2TecooMANmJ4YPAhCLswEdHVQgEGQuLaGBLHCTkqznCJDGljymKUYRLbp2ePdq2E9Wkva0ZKbdNXXLo5KXXskiagPE5Yfk1X8VhVluxPtYeqpuxtQG8rqeURK3BO7F+6SeZ/ddKLQBcl5cEhT6FipDs7qRKrqcjO9U+uBFUIiZO4Q3FO2FbKyZx1pYmUdzlbU8VfVKZk5d5OcmO8/upOdjCeNZ0L/QaF74c+DerFfGE/X6wauXIJ45XGFnawR6d+gDKOUUFNdXCNW9xprnaj3WsP41QtlWg28WVpjx+sX5wJ</latexit>
X
<latexit sha1_base64="p1j1zoOBRDWaYar4rxVus3IOhfg=">AAAC6nicjVLLSsNAFD2Nr1pfVZeCBIvgqqR1YXcW3LhswT6gFknSaQ3Ni8lEKMWlK3eiC3HrR3Tlf/gHPn7CO9MU1OJjQpIz595zc8+dWKHrRMIwnlPazOzc/EJ6MbO0vLK6ll3fqEdBzG1WswM34E3LjJjr+KwmHOGyZsiZ6Vkua1j9IxlvXDAeOYF/IgYha3tmz3e6jm0KoqrNs2zOyBtq6dOgkIDc4dOo+n61PaoE2TecooMANmJ4YPAhCLswEdHVQgEGQuLaGBLHCTkqznCJDGljymKUYRLbp2ePdq2E9Wkva0ZKbdNXXLo5KXXskiagPE5Yfk1X8VhVluxPtYeqpuxtQG8rqeURK3BO7F+6SeZ/ddKLQBcl5cEhT6FipDs7qRKrqcjO9U+uBFUIiZO4Q3FO2FbKyZx1pYmUdzlbU8VfVKZk5d5OcmO8/upOdjCeNZ0L/QaF74c+DerFfGE/X6wauXIJ45XGFnawR6d+gDKOUUFNdXCNW9xprnaj3WsP41QtlWg28WVpjx+sX5wJ</latexit>

b)
<latexit sha1_base64="QzEjC1xKBjZv1Tdc8Z8xI4hREoU=">AAAC63icjVLLSsNAFD2Nr7a+qm4EN0ERdFPSurDLohuXVewDapFkOq2heZFMhFoEP8Cd6ELc+gP+jeAHqF/hnWkKavExIcmZc+85M/fOWIFjR8IwnlPaxOTU9Ew6k52dm19YzC0t1yI/DhmvMt/xw4ZlRtyxPV4VtnB4Iwi56VoOr1u9fRmvn/Mwsn3vWPQD3nLNrmd3bGYKoo6s7dPchpE31NDHQSEBG+XVi5fM1dNexc+94wRt+GCI4YLDgyDswERETxMFGAiIa2FAXEjIVnGOS2RJG1MWpwyT2B59uzRrJqxHc+kZKTWjVRx6Q1Lq2CSNT3khYbmaruKxcpbsT94D5Sn31qe/lXi5xAqcEfuXbpT5X52sRaCDkqrBppoCxcjqWOISq67IneufqhLkEBAncZviIWGmlKM+60oTqdplb00Vf1WZkpVzluTGePu1OrmDYa/pXOgaFL4f+jioFfOFnXzxkO5DCcORxhrWsUWnvosyDlBBlVbp4Bq3uNNc7Ua71x6GqVoq0azgy9AePwCqAZsl</latexit>

X
<latexit sha1_base64="p1j1zoOBRDWaYar4rxVus3IOhfg=">AAAC6nicjVLLSsNAFD2Nr1pfVZeCBIvgqqR1YXcW3LhswT6gFknSaQ3Ni8lEKMWlK3eiC3HrR3Tlf/gHPn7CO9MU1OJjQpIz595zc8+dWKHrRMIwnlPazOzc/EJ6MbO0vLK6ll3fqEdBzG1WswM34E3LjJjr+KwmHOGyZsiZ6Vkua1j9IxlvXDAeOYF/IgYha3tmz3e6jm0KoqrNs2zOyBtq6dOgkIDc4dOo+n61PaoE2TecooMANmJ4YPAhCLswEdHVQgEGQuLaGBLHCTkqznCJDGljymKUYRLbp2ePdq2E9Wkva0ZKbdNXXLo5KXXskiagPE5Yfk1X8VhVluxPtYeqpuxtQG8rqeURK3BO7F+6SeZ/ddKLQBcl5cEhT6FipDs7qRKrqcjO9U+uBFUIiZO4Q3FO2FbKyZx1pYmUdzlbU8VfVKZk5d5OcmO8/upOdjCeNZ0L/QaF74c+DerFfGE/X6wauXIJ45XGFnawR6d+gDKOUUFNdXCNW9xprnaj3WsP41QtlWg28WVpjx+sX5wJ</latexit>

µ
0 µ⇤ <latexit sha1_base64="MemXgHv8JLvWE4CUAdBykVLay8Q=">AAAC7HicjVLLTsJAFD3UF+ILZemmkZC4IgUXsiRx4xKjPBIwpi0DNvSV6dSEEHdu3RldGLd+iv+gH2DUtR/gnaEkKvExTdsz595zZu6dsULXiYRhPKa0mdm5+YX0YmZpeWV1Lbu+0YiCmNusbgduwFuWGTHX8VldOMJlrZAz07Nc1rQGezLePGM8cgL/SAxDduyZfd/pObYpiDrsePFJNm8UDTX0aVBKQL6aK1y8Pzw/1YLsGzroIoCNGB4YfAjCLkxE9LRRgoGQuGOMiOOEHBVnOEeGtDFlMcowiR3Qt0+zdsL6NJeekVLbtIpLLyeljgJpAsrjhOVquorHylmyP3mPlKfc25D+VuLlEStwSuxfuknmf3WyFoEeKqoGh2oKFSOrsxOXWHVF7lz/VJUgh5A4ibsU54RtpZz0WVeaSNUue2uq+IvKlKyc20lujNdfq5M7GPeazoWuQen7oU+DRrlY2imWD+g+VDAeaWxiC9t06ruoYh811GmVPi5xjRvN1660W+1unKqlEk0OX4Z2/wGbgJ0t</latexit>

Figure . – a) Fonction F (X) dans ses cinq configurations. b) Diagramme

de bifurcation.
 Chapitre . Exercices complémentaires

EXERCICE . Équations de Saint-Venant

On considère le système des équations de Saint-Venant

∂h ∂(U h) ∂U ∂U ∂h g |U | U
+ =0, +U = −g cos α +g sin α− 2 4/3 , (.)
∂t ∂x ∂t ∂x ∂x K h
où U (x, t) et h(x, t) sont respectivement la vitesse et la hauteur d’une couche de
fluide peu profonde à surface libre coulant sur un plan incliné faisant un angle
α avec l’horizontale. La gravité est noté g et K est le “nombre de Strickler” qui
caractérise le frottement du fluide sur le fond. On suppose que U ≥ 0.
) Montrer que les solutions homogènes et stationnaires (hn , Un ) vérifient
2/3 √
Un = K hn I où I est une constante que l’on déterminera. Quelle est
la dimension du nombre de Strickler ?

2/3 √
On a Un = K hn I avec I = sin α. L’unité de K est en m1/3 /s.
) On choisit hn comme unité de longueur, Un comme unité de vitesse et
hn /Un comme unité de temps. Écrire les équations √ ainsi adimensionnées
en faisant apparaître le nombre de Froude F r = Un / g cos α hn ainsi que
−1/3
le nombre sans dimension M = g K −2 hn .

En notant les champs adimensionnés de la même façon que les champs

 avec dimen-
∂(U h) |U | U
sions, on obtient ∂h
∂t + ∂x = 0 et ∂U
∂t + U ∂U
∂x = − 1 ∂h
2
F ∂x + M 1 − h4/3
.

) Calculer la relation de dispersion des ondes dans le cas M = 0 et F = O(1).

Exprimer les vitesses de phases en revenant aux grandeurs dimensionnelles.
Ces ondes sont-elles dispersives ?

Dans le cas M = 0, la linéarisation U = 1 + U e et 1 + eh des équations conduit à

∂U
+ ∂U
= − 1 ∂e
h
et ∂e
h
+ ∂e
h
+ ∂U
= 0. La recherche de solution de la forme
e e e
∂t ∂x 2
F ∂x ∂t ∂x ∂x
(U h) = (A, B) exp(i k x + s t) conduit aux équations (s + i k)A = − F12 i k B et
e, e
(s + i k)B + i k A = 0. En annulant le déterminant de ce système linéaire on obtient
(s + i k)2 + k 2 /F 2 = 0. On en déduit s = −i ω avec ω = k (1 ± 1/F ). En grandeurs
dimensionnées, la vitesse de phase de ces ondes non dispersives est cϕ = Un ±
√
g cos α hn .
) Calculer la relation de dispersion des ondes dans le cas 1/M = 0 et
F = O(1). Exprimer les vitesses de phases en revenant aux grandeurs
dimensionnelles. Ces ondes sont-elles dispersives ?

Dans le cas M infini, on a U = h2/3 en grandeur adminensionnées et donc ∂h ∂t +

∂
∂x (h 5/3
) = 0. Les petites perturbations de l’équilibre (U, h) = (1, 1) vérifient ∂e
h
∂t +

3 ∂x = 0. La relation de dispersion des ondes est alors ω = 3 k. La vitesse de phase

5 ∂he 5

de ces ondes non dispersives est cϕ = 3 Un en grandeur dimensionnée.

5

) Dans le cas M = O(1) et F = O(1), écrire les équations linéarisées autour

de l’équilibre.

Dans le cas général, les équations linéarisées sont ∂U

+ ∂∂x
U
= − F12 ∂e
h
∂x −2 M
e+4 M e
U h
e e
∂t 3
et ∂∂th + ∂∂xh + ∂∂x
U
= 0.
e e e

) Montrer que la stabilité de l’équilibre dépend des racines complexes d’un

polynôme de degré deux en s = σ − i ω dont on explicitera les coefficients
en fonction des données du problème.
Exercice sur le chapitre  

La recherche de solution de la forme (U, h) = (A, B) exp(i k x + s t) conduit aux

équations (s + i k)A = − F12 i k B − 2 M A + 34 M B et (s + i k)B + i k A = 0. En
annulant le déterminant de ce système linéaire, on obtient l’équation (s + i k +
2 M )(s + i k) − i k i k F12 − 43 M = 0.

) Montrer que l’écoulement de base se déstabilise pour un nombre de Froude

critique Fc que l’on calculera. On pourra poser σ = 0 et calculer les parties
réelles et imaginaires du polynôme.

La courbe marginale dans l’espace des paramètres de contrôle est déterminée par
la relation σ = 0. Le polynôme
 caractéristique
 en s = −i ω s’écrit alors [i(−ω +
k) + 2 M ] i (−ω + k) − i k i k F 2 − 3 M = 0. En annulant les parties réelles et
1 4
c
imaginaires, on obtient (−ω + k) = − 23 k et (−ω + k)2 = k 2 /Fc2 d’où ω = 35 k et
Fc = 3/2.