Universite Paris-Sud Annee 2010-11

Master 1 de Physique Appliquee et Mecanique

Examen de Turbulence
Jeudi 31 Mars 2011
Duree : 3 heures - sans document

Exercice 1 : Sillage dun objet propulse

U0

x
l(x)

Figure 1: Sillage bi-dimensionnel dun objet propulse, montrant un prol de vitesse non monotone en
aval de lobjet.

Un objet propulse, comme un avion ou un sous-marin, avance grace a laction dune source
denergie interieure, donnant lieu a une force appelee poussee. Lorsque cet objet propulse
avance a vitesse constante, la poussee est exactement compensee par la force de tranee. Le
sillage dun tel objet propulse diere notablement de celui dun objet passif, comme une balle
de tennis, qui ne possede pas de source denergie interieure. Dans le cas dun objet propulse,
on peut montrer en particulier que le prol de vitesse moyen est non monotone (gure 1).
Lobjectif de ce probleme est de caracteriser ce sillage.
On considere pour simplier un objet propulse bi-dimensionnel, invariant selon la direction
z, et on etudie son sillage turbulent dans le plan (x, y). Lobjet se deplace vers la gauche a
vitesse constante U0 . On se place dans le referentiel de lobjet, de sorte quil existe un
ecoulement uniforme U0 > 0 loin de lobjet. On note (u, v) la vitesse du uide selon (x, y),
avec x = 0 a larriere de lobjet. Le nombre de Reynolds est suppose tres grand, et lon se
place loin de lobjet (x grand), de sorte que leet de la viscosite peut etre neglige ici. On
note (x) la largeur caracteristique du sillage a la distance x, et lon suppose (x)  x. On
sinteresse au defaut de vitesse u(x, y) = U0 u(x, y) dans le sillage. On note Us (x) = u(x, 0)

1
le defaut de vitesse au centre, et lon suppose Us  U0 . Contrairement au sillage dun objet
passif, il ny a pas de debit associe au defaut de vitesse dans le cas de lobjet propulse, et lon
a donc : 
u(x, y)dy = 0.

1. A partir du prol de vitesse moyen de la gure 1, indiquer pour quelles valeurs de y la

contrainte visqueuse change de signe. Que pensez-vous de la contrainte turbulente en
y=0?

2. On rappelle lequation de Reynolds pour la vitesse moyenne ui :

ui 1 p   2 ui
ui + uj = ui uj + 2 . (1)
t xj xi xj xj

Ecrire les projections selon x et y de cette equation.

3. Expliquer, en quelques lignes de justication, pourquoi lequation de Reynolds selon x

peut se simplier sous la forme :

u u v 
U0 = . (2)
x y
On pourra admettre ce resultat pour la suite.

4. On utilise un modele de viscosite turbulente pour le tenseur de Reynolds :

u
u v  = T (x) . (3)
y

Que pensez-vous du choix T (x) independant de y ?

5. En multipliant lequation (2) par y 2 et en integrant, montrer que lintegrale


I0 = y 2 u(x, y) dy (4)

est independante de x (on supposera que u tend vers 0 plus vite que 1/|y|2 pour |y|
).

6. On admet dans la suite que le sillage est autosimilaire, et lon pose

u(x, y) = Us (x)f () (5)

avec = y/(x) la coordonnee reduite transverse. En exprimant lintegrale I0 a partir

de ce prol autosimilaire, en deduire que lon a

Us (x)(x)3 = J0 , (6)

ou lon exprimera J0 en fonction de I0 et dune integrale sans dimension.

2
 7. Justier physiquement que le coecient de viscosite turbulente puisse etre pris sous la
forme
T (x) = (x)Us (x),
ou est une constante de lordre de 1.

8. Exprimer les derivees partielles /x et /y en fonction de /, , (x) et d/dx.

An de determiner les variations de Us (x) et (x) en fonction de x, on introduit le prol
auto-similaire (5) dans lequation de Reynolds (2). Montrer, en vous aidant de (6), que
lon obtient alors
1 d
= cste. (7)
Us (x) dx

9. Les 2 lois de conservations mises en evidence permettent decrire Us (x) et (x) sous la
forme
Us (x) xn , (x) xm . (8)
En deduire les valeurs des exposants n et m. Dessiner lallure de levolution de la largeur
du sillage. En raisonnant sur le nombre de Reynolds turbulent ReT (x) = Us (x)(x)/,
que pensez-vous du sillage lointain ?

Exercice 2 : Ecoulement turbulent dans un tuyau cylindrique

On considere un ecoulement deau dans une portion de tube circulaire daxe Ox, de rayon
R = 5 cm et de longueur L = 20 m (gure 2). Une dierence de pression de p = 2 105 Pa
est imposee entre lentree et la sortie du tube, et lon cherche a calculer le debit dans le tube.
On neglige les eets detablissement du prol de vitesse a lentree du tube (L  longueur
dentree). On donne la densite de leau, = 1000 kg m3 , et la viscosite cinematique, =
106 m2 s1 .

R
U0
O x

y u(y)

Figure 2: Prol de vitesse moyen dun ecoulement turbulent dans un tube.

1. En faisant un bilan des forces sexercant sur les sections dentree et de sortie et sur la
paroi laterale du tube, calculer la contrainte parietale 0 .

2. Rappeler lexpression de la vitesse de friction u en fonction de 0 . Calculer u , et en

deduire le nombre de Reynolds turbulent Re = u R/, ainsi que lepaisseur de la
sous-couche visqueuse v (faire les applications numeriques).

3. Rappeler (en quelques lignes) la signication physique des notions de sous-couche vis-
queuse et de sous-couche inertielle, en particulier du point de vue des mecanismes de
transport de quantite de mouvement.

3
 Dans lapproximation dune sous-couche visqueuse tres ne comparee au rayon R du tube,
on peut approximer lecoulement pres de la paroi a celui dune couche limite au-dessus dune
plaque plane. On note y = R r la distance a la paroi et U0 la vitesse au centre du tube.
On rappelle les deux expressions du prol de vitesse moyen dans la sous-couche inertielle
(v  y  R) :

u(y) 1 y u(y) U0 1 y

= ln +C (a) et
= ln + C  (b). (9)
u v u R
Les valeurs experimentales des constantes sont 1/  2, 5 (constante de von Karman), C  5
et C   0.

4. En admettant que la loi (9b) reste approximativement valable sur toute la section du
tube, montrer que le debit volumique vaut :
 
2 3 u
Q = R U0 .
2

On pourra utiliser les integrales suivantes :

 1  1
ln x dx = 1 et x ln x dx = 1/4.
0 0

5. A partir des deux lois (9a) et (9b), exprimer la vitesse au centre U0 en fonction de u , du
nombre de Reynolds Re , ainsi que des constantes sans dimension. Calculer les valeurs
numeriques de U0 et de Q.

6. En deduire que le debit Q peut secrire sous la forme

Q = R2 U0 f (Re )

ou f (Re ) est une fonction sans dimension dependant de Re uniquement. Tracer lallure
de cette fonction, en discuter physiquement la limite Re . Que pensez-vous de 0
dans cette limite ?