18ème Congrès Français de Mécanique Grenoble, 27-31 août 2007
Modélisation au second ordre d'une turbulence
homogène cisaillée en présence de la rotation
Besma Chebbi1, Mounir Bouzaiane1* & Taieb Lili1
Laboratoire de Mécanique des Fluides
Département de Physique
Faculté des Sciences de Tunis -Tunisie
*
Email mbouzaiane@[Link], [Link]@[Link]
Résumé :
Dans ce travail, nous étudions à l’aide de modèles au second ordre les effets de la rotation sur
l’évolution des champs cinématique et thermique d’une turbulence homogène cisaillée. Le modèle de
Launder –Reece Rodi (LRR) d’une part et le modèle de Sarkar Speziale et Gatski (SSG) d’autre part,sont
retenus pour les termes de corrélation pression déformation et l’équation d’évolution de la dissipation de
l’énergie cinétique alors que pour la corrélation pression température et pour l’équation d’évolution de
la dissipation thermique seulement le modèle (LRR) est retenu. La méthode de Runge-Kutta du
quatrième ordre est utilisée pour intégrer deux systèmes différentiels non-linéaires. L’intégration
numérique a été menée pour des différentes valeurs du nombre adimensionnel de rotation R. Les résultats
obtenus ont montré que la rotation affecte aussi bien le champ cinématique que le champ thermique. Le
modèle SSG est d’un apport considérable sur la modélisation du champ thermique.
Abstract
In this work, effects of rotation on kinematics and thermal field evolutions of homogeneous sheared
turbulence are studied using second order models of turbulence. The Speziale Sarkar Gatski models
(SSG) in one hand and the classic Launder-Reece-Rodi in the other hand are retained for the pressure
strain correlation and for the time evolution equation of the cinematic dissipation whereas only the
(LRR) model is retained for the pressure-temperature correlation and for the time evolution equation of
the dissipation of the variance of the temperature fluctuation. A fourth order Runge-Kutta method is used
for numerical integrating of two systems of nine non linear differential equations obtained after
modeling. The numerical integration has been conducted for different values of the non-dimensional
rotation number R and has shown that rotation not only affects the kinematic field but also the thermal
field. The SSG model shows a great superiority than the LRR one.
Mots-Clefs :
Turbulence cisaillée, effets de la rotation, modélisation au second ordre.
Key words :
Sheared turbulence, rotation effects, second order modeling.
1
18ème Congrès Français de Mécanique Grenoble, 27-31 août 2007
1 Introduction
L’étude d’une turbulence à l’aide des modèles au second ordre demeure l’une des approches
importantes menées par les auteurs en vue de comprendre et analyser les écoulements
turbulents. C’est à l’aide des modèles que second ordre que nous avons étudié des différentes
configurations d’une turbulence homogène cisaillée soumise à une stratification stable
Bouzaiane et al. 2003, 2004,. Le cas d’une turbulence homogène cisaillée associée à un champ
thermique est l’un des écoulements aux quels nous nous sommes intéressés les dernières
années Bouzaiane et al. 2000. La modélisation au second ordre d’une turbulence homogène
cisaillée, en présence de la rotation, associée à un champ thermique constituant une extension
de ce dernier travail au cas de la rotation est une première dans nos travaux et n’a pas fait, à
notre connaissance, l’objet d’une étude antérieure. Ce ci constitue une motivation essentielle de
ce travail en plus de la place considérable que prennent les liquides et les gaz dans les systèmes
en rotation. Les résultats récents des simulations numériques directes de G. Brethouewer 2005,
sont à notre connaissance les plus récents pour le type d’écoulement étudié et sont retenus dans
notre travail.
Nous commençons au paragraphe suivant par présenter les équations générales de l’écoulement.
Nous ensuite le problème de fermeture des équations et nous présentons le modèle au second
ordre retenu dans ce travail. Dans ce même paragraphe, des grandeurs adimensionnelles
cinématiques et thermiques sont introduites. L’intégration numérique des équations de
l’écoulement pour des différentes valeurs du nombre adimensionnel de rotation R et une
présentation des résultats obtenus par les différents modèles sera effectuée au paragraphe 3. Un
résumé et une perspective du travail présenté dans ce papier feront l’objet du paragraphe 4.
2 Equations générales de l’écoulement, modélisation au second ordre
Nous précisons ici que nous étudions une turbulence homogène cisaillée associée à un champ
thermique présentant un gradient moyen de température dans la même direction que le
cisaillement moyen x2. La rotation étant selon la direction x3. Ce cas de turbulence est décrit par
les équations d’évolution des tensions de Reynolds, des flux thermiques turbulents , de
l’énergie cinétique turbulente K, de la variance de la fluctuation de la température et des
dissipations cinématiques et thermiques . Les formes condensées de ces équations
d’évolution sont les suivantes :
(1)
(2)
(3)
Dans ces équations, les termes notés P sont des termes de production dus d’une part aux
gradients cinématique S=U1,2 et de scalaire moyen S =T,2 d’une part et à la rotation d’axe x3
d’autre part
(4)
(5)
(6)
2
18ème Congrès Français de Mécanique Grenoble, 27-31 août 2007
Les termes notés sont les termes de corrélation faisant intervenir la fluctuation de pression :
(7)
(8)
A partir de ces équations dévolution, nous déduisons les équations d’évolution des
tensions de Reynolds de l'énergie cinétique turbulente ainsi que les
équations d’évolution des flux thermiques turbulents :
(9)
(10)
(11)
(12)
(13)
R est ici le nombre adimensionnel de rotation .
2.1 Modélisation au second ordre
A ce niveau du travail, une fermeture au second ordre des équations déduites est abordée.
Cette fermeture consiste en la modélisation au second ordre des termes de corrélation faisant
intervenir les fluctuations de pression. Ces termes appelés respectivement termes de corrélation
pression-déformation et pression-gradient d'un scalaire sont décomposés en deux
contributions :
Les termes 1 sont dits termes de retour à l'isotropie et traduisent l'interaction de la turbulence sur
elle-même alors que les termes notés 2 sont les termes linéaires et représentent l'interaction
entre le mouvement moyen et le mouvement turbulent.
La modélisation des termes linéaires est basée sur l’écriture d’une équation de poisson pour la
fluctuation de pression. En présence de cisaillement et de la rotation, cette équation est écrite
sous la forme :
3
18ème Congrès Français de Mécanique Grenoble, 27-31 août 2007
La solution de cette équation est à l’origine des modèles des termes su-mentionnés qui
s’écriront :
Les tenseurs A et B sont respectivement des tenseurs du quatrième et du troisième ordre
vérifiant un groupe de relations dites de contraintes cinématiques (Ristorcelli et al. 1995). Alors
que désigne la composante de la rotation selon la direction .
Dans ce travail, le modèles de Sarkar Speziale Gatski (1989) (SSG) d’une part et le modèle
classique de Launder-Reece et Rodi (LRR) d’autre part sont retenus pour les corrélations
pression-déformation et pour l’équation d’évolution de la dissipation de l’énergie cinétique. Le
modèle (LRR) est à lui seul retenu pour les corrélations pression-température et équation
d’évolution de la dissipation de la variance de la fluctuation de la température.
Une forme adimensionnelle des équations d'évolution est ensuite obtenue en introduisant des
grandeurs adimensionnelles cinématiques et thermiques . Il s'agit des équations d'évolution des
composantes du tenseur anisotrope de Reynolds, du coefficient de corrélation
associé au flux thermique turbulent, du rapport des flux thermiques
turbulents et du rapport adimensionnel . L'intégration numérique du système
d'équations différentielles non linéaires et la présentation des résultats obtenus feront l'objet du
paragraphe suivant.
3 Résolutions numériques et résultats
La méthode de Runge–Kutta d'ordre 4 est utilisée ici pour l'intégration numérique des
équations différentielles non-linéaires obtenues après la modélisation au second ordre. Les
conditions initiales sont celles retenues par G. Brethouewer dans ses simulations numériques
directes. L'intégration numérique est menée pour des différentes valeurs du nombre de rotation
R comprises entre –1,5 et 0,5 avec un pas de 0,5.
Nous présentons aux tableaux 1 et 2 les valeurs d’équilibre des grandeurs adimensionnelles,
obtenues en poursuivant l’intégration numérique vers des intervalles de temps suffisamment
grands. Nous constatons d’après ces deux tableaux que les résultats obtenus par le (LRR) ne
prédisent des valeurs d’équilibre que pour les valeurs du nombre de rotation R=0,0 et R=0,5
alors que. Le modèle SSG couplé au LRR pour les phénomènes thermiques prédits quant à lui
des valeurs d’équilibre pour toutes les valeurs considérées du nombre adimensionnel de
rotation.
b11 b22 b12 1 2 3
R= 0,00 0,192 -0,090 -0,183 0,181 0,484 0,346 0,658
R= 0,50 -0,119 -0,310 -0,415 0,206 0,882 0,232 0,514
4
18ème Congrès Français de Mécanique Grenoble, 27-31 août 2007
R= 1,00 …. … … …. 4,52 -0,623 26,3
Tableau 1 Valeurs d’équilibre prédites par le modèle LRR
b11 b22 b12 1 2 3
R= 0,00 0,220 -0,147 -0,164 0,167 0,298 -1,11 2,44
R= 0,50 0,537 -0,201 -0,365 0,387 0,621 -1,46 2,27
R= 1,00 0,809 -0,235 -0,550 0,061 ….. -1,8 …..
R=-0,50 -0,040 -0,123 -0,166 0,172 0,461 -1,56 1,81
R=-1,00 -0,271 -0,123 -0,291 0,309 1,35 -2,63 0,93
R=-1,50 -0,516 -0,164 -0,402 0,426 ….. ….. ….
Tableau 2 Valeurs d’équilibre prédites par le modèle SSG+LRRR
L’anisotropie augmente en présence de la rotation et une forte anisotropie est observée pour le
nombre adimensionnel de rotation R égal à 1,0. La valeur maximale des composantes
d’anisotropie bij atteint 0,809. Une anisotropie moins forte est constatée pour les valeurs
négatives du nombre de rotation.
Sur ses mêmes tableaux sont aussi présentées les valeurs d’équilibre des grandeurs
adimensionnelles thermiques . La différence est nette entre les résultats des deux
tableaux 1 et 2. Le modèle SSG retenu pour le champ cinématique est d’un grand apport sur le
champ thermique. Ce modèle n’intervenant pas explicitement dans les équation d’évolution
décrivant le champ thermique a permis d’améliorer les résultats obtenus par le modèle de LRR
quand il est retenu à lui seul pour les champs cinématique et thermique et a permis d’obtenir des
valeurs d’équilibre pour 4 des 6 valeurs retenus pour le nombre de rotation R. La rotation
affecte clairement le champ thermique et c’est essentiellement le coefficient de corrélation
qui semble le plus touché par la rotation. La valeur d’équilibre de ce coefficient double en
passant de R=0,0 à R=0,50. Une aberration est aussi constatée pour la valeur du nombre de
rotation égale à –1,0. La valeur du coefficient de corrélation du flux thermique turbulent
dépasse l’unité. Cette aberration est expliquée si on rappelle que le modèle classique (LRR) de
la corrélation pression-température n’est pas soumis à des conditions strictes de réalisabilité
telles que celles introduites par Shih et Lumley . Ce modèle n’est soumis qu à des relations
simples dites contraintes cinématiques.
4 Conclusions
Dans ce travail, nous avons étudié les champs cinématique et thermique d’une turbulence
homogène soumise aux effets de la rotation à l'aide des modèles classiques au second ordre.
Nous nous sommes intéressés en particulier à la détermination des valeurs d’équilibre des
grandeurs adimensionnelles cinématiques et thermiques.
Une forme adimensionnelle des équations a été obtenue après l’introduction de grandeurs
adimensionnelles cinématiques et thermiques. L’intégration numérique des équations a été
effectuée par la méthode de Runge-Kutta d’ordre 4 et a montré que les évolutions prédites par
les modèles au second ordre de SSG confirme l’existence d’états asymptotiques d’équilibre
pour les grandeurs adimensionnelles cinématiques et pour les différentes valeurs considérées du
nombre de rotation R.
5
18ème Congrès Français de Mécanique Grenoble, 27-31 août 2007
Les valeurs d’équilibre des grandeurs adimensionnelles cinématiques et thermiques sont
fonctions du nombre adimensionnel de rotation R.
Le modèle (LRR) est incapable de décrire un champ thermique d’une turbulence homogène
cisaillée en présence de la rotation. Une modélisation des équation de l’écoulement à l’aide des
modèles sophistiqués de Craft et Launder d’une part, de ceux de Shih et Lumley d’autre part
dans le but d’obtenir une meilleure prédiction du champ thermique est en vue .
References
Bouzaiane M., Ben Abdallah H.& Lili T. 2003 A Study of the asymptotic behaviours of
dimensionneless parameters in a stably homogeneous sheared turbulence J. Turbulence, 4(1)
Bouzaiane M., Ben Abdallah H.& Lili T. 2004 A Second order modeling of a stably stratified
sheared turbulence submitted to a non-vertical shear J. Turbulence , 5(9)
Bouzaiane M. & Lili T. 2002 Etude des états d’équilibre d’une turbulence homogène associée à
un champ thermique I. J. Ther. Sc. 41(9) 871-882
Brethouewer G. 2005 The effect of rotation on rapidly sheared homogeneous turbulence and
passive scalar transport, J. Fluid Mech. 542 305-342.
Launder B.E. ,Reece G., & Rodi W. 1975 Progress in the development of a Reynolds stress
closure, J. Fluid Mech. 68 537-576.
Speziale C.G. , Sarkar S. & Gatski T.B. 1990a Modeling the pressure strain correlation of
turbulence-an invariant dynamical system approach, J. Fluid Mech. 227 , 245-272.
Shih, T. H. 1996 Constitutives relations and realisability of single-point turbulence closures,
in Turbulence transition and Modeling edited by Mr. Hallback D. S. Hennigson and A. V.
and Johansson and P. H. Alfredsson., Dordrecht