TD Turbulence

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

Université Ibn Khaldoun de Tiaret-Faculté des Sciences Appliquées

Département de Génie Mécanique-Master 1 Energétique


MDF Approfondie - TD Turbulence – Année :2021 - 2022

Exercice 1 : Propriétés de la moyenne


En utilisant la relation suivante :
1 t0 T 2
  (t)dt ,
T t0 T 2
Montrer que :
d d 
a)  , b)   0 c)        , d) .  . ,e)  .    . '  0 f) .  . ,g) .  .  .  , h)  i)
ds ds
 1 
 ds   ds j)      
      

Exercice 2 : Équations de la turbulence :


En utilisant la décomposition de Reynolds dans les équations de base, retrouver les équations suivantes :
u j
a) 0,
x j
uj
b) 0,
x j

c)  u i  u j  u i   1  P    ij  Rij 
t x j  x i x j

avec Rij   ui uj le tenseur de Reynolds

Exercice 3 : Ecoulement de Couette turbulent 1 :


Soit l’écoulement de Couette turbulent entre deux plans infiniment longs (figure 1) qui se déplacent dans des directions
opposées. L’écoulement se fait à masse volumique constante, et sa vitesse moyenne u dépend uniquement de la composante y.
Les forces de volume sont négligeables et le gradient de pression est nul selon l’axe des x et toutes les quantités physiques
sont supposées statistiquement stationnaires et invariantes par translation selon x.
Les contraintes turbulentes de Reynolds sont déterminées en utilisant l’approche de la longueur de mélange de Prandtl donnée
par la relation :
l(y)  K(l 2  h2) (1)

Figure 1 : Ecoulement de Couette turbulent


1) Déterminer la constante K qui satisfait à la condition :
dl
   (2)
dy yh
où κ est la constante de Prandtl.
2) Déterminer l’équation de la contrainte de cisaillement turbulente  t pour la longueur de mélange donnée.
Rappel : le modèle de Prandtl de la longueur de mélange est :

1
Université Ibn Khaldoun de Tiaret-Faculté des Sciences Appliquées
Département de Génie Mécanique-Master 1 Energétique
MDF Approfondie - TD Turbulence – Année :2021 - 2022

du du
t   uv  l 2 (3)
dy dy

3) Comme la pression moyenne P dans l’écoulement de Couette est constante, la contrainte de cisaillement est aussi
constante. Comme résultat, nous avons :

du
μ   u v    u f2  C te (4)
dy

p
avec uf  la vitesse de frottement.

A l’extérieur de la sous limite visqueuse, la contrainte de cisaillement visqueuse (ou laminaire) μ du peut être négligée
dy
comparée à la contrainte de cisaillement turbulente.
Comme condition de symétrie, nous avons u ( y  0)  0 .
Déterminer le profil de vitesses u (y) .
4) Déterminer le profil de la vitesse moyenne en fonction de la distance à la paroi inférieure y   y  h .
5) Montrer que pour les petites valeurs de y  .
6) Déterminer le profil des vitesses pour le cas suivant :
2hU p
- Nombre de Reynolds : Re   34000 , - la constante de la loi logarithmique B5, - la constante de Karman :

  0, 4

Exercice 3 : Ecoulement de Couette turbulent

On considère un écoulement de Couette entre deux plaques planes parallèles infinies séparées d’une distance 2h , animées
d’une vitesse  2Up l’une par rapport à l’autre (figure 2). On se place dans le référentiel tel que la plaque supérieure ait une
vitesse  Up et la plaque inferieur une vitesse  U p . On note Re  Up h  le nombre de Reynolds de l’écoulement.

La direction transverse (perpendiculaire au plan de la figure) étant très grande devant h , on peut considérer le problème
comme purement bidimensionnel. Il n’existe pas de gradient de pression moyen selon x , et toutes les quantités physiques sont
supposées statistiquement stationnaires et invariantes par translation selon x .

Figure 2 : Géomètrie de l’écoulement de Couette plan, le profil de vitesse correspondant au cas laminaire.

1) On sait que, dans le cas laminaire, le profil de vitesse d’une telle géomètrie est donné par :
u x (y)  U p y h

2
Université Ibn Khaldoun de Tiaret-Faculté des Sciences Appliquées
Département de Génie Mécanique-Master 1 Energétique
MDF Approfondie - TD Turbulence – Année :2021 - 2022

On cherche à déterminer si, dans le cas turbulent, le profil de vitesse moyen u x (y) peut également s’écrire sous cette forme.
En introduisant la solution u x (y)  U p y h dans l’équation de Reynolds selon x , montrer que u x u y  0 dans tout l’écoulement.
Que concluez-vous ?
2) On va chercher à déterminer le profil moyen u x (y) dans le cas turbulent. On introduit la coordonnée réduite   y h , et
on pose :

 u x (y)  U p f ()

 (1)
2
 u x u y  U p g()
Montrer que l’on a :
 2f g
 Re (2)
2 

Préciser la parité et les conditions aux limites en   1 des fonctions f () et g() .
3) Pour déterminer la fonction f () , il faut une relation supplémentaire reliant f et g . Pour cela, on peut introduire un
coefficient de viscosité turbulente  T (y) tel que :
u x
u x u y  T (y) (3)
y
On suppose que les fluctuations responsables du transfert de quantité de mouvement sont de vitesse caractéristique Up et de
taille l (y) , et l’on pose donc :
T (y)  Up l (y)
où  est une constante sans dimension. Justifier que l’on peut choisir l  h  y .
Par symétrie, on pourra ne considérer que la moitie supérieure (y  0) , et ainsi oublier la valeur absolue. En déduire une
relation entre g et f  puis, en utilisant l’équation (2), obtenir l’équation différentielle suivante : (1  ( Re)1  )   où l’on
a pose ()  f ()

4) En introduisant le changement de variable Z   (1  ( Re)1 ) , intégrer cette équation différentielle (toujours pour   0
).

Montrer finalement que le profil de vitesse s’écrit :

ln(1   Re(1  ))


f ( )  1  (4)
ln(1   Re)

La courbe de cette fonction est donnée dans la figure (XXX) pour différentes valeurs du nombre de Reynolds. Vérifier (par
un développement limité) que l’on retrouve bien la solution laminaire dans la limite  Re  1 . Que pensez-vous de la limite
Re  

3
Université Ibn Khaldoun de Tiaret-Faculté des Sciences Appliquées
Département de Génie Mécanique-Master 1 Energétique
MDF Approfondie - TD Turbulence – Année :2021 - 2022

Figure 3: Profils de vitesse adimensionnés (4) pour différentes valeurs de  Re .

Exercice 4 : Dispersion d’un polluant :


Soit C(x, y, z, t) la concentration d’un polluant dans un écoulement turbulent. L’équation de diffusion-advection de la
concentration est :

C C  2C
 ui D (1)
t x i x i x i

Où D est la constante de diffusion du polluant et u i es le vecteur vitesse de l’écoulement. En introduisant les décompositions
de Reynolds pour le champ de vitesse et le champ de concentration, montrer que la concentration moyenne vérifie l’équation :

C C   C 
 ui   C u i  D
 
t x i x i  x i 

Parmi les deux contributions du membre de droite, laquelle sera dominante dans un écoulement atmosphérique ?

Vous aimerez peut-être aussi