Modelisation Des Ecoulements Turbulents A Bulles: These
Modelisation Des Ecoulements Turbulents A Bulles: These
Modelisation Des Ecoulements Turbulents A Bulles: These
THESE
et de
par
Bellakhal Ghazi
Les milieux diphasiques mettent en jeu des interactions aux interfaces qui modifient d’une
maniere significative la structure des champs moyens et fluctuants des ecoulements. La conception
de modeles a deux fluides adaptes aux ecoulements industriels necessite la prise en compte de
l’effet de ces interactions au niveau des relations de fermeture adoptees.
Le travail developpe dans cette these porte sur le developpement de modeles a deux fluides au
premier ordre deduits par reduction des fermetures au second ordre. La demarche adoptee, basee
sur le principe de decomposition du tenseur des contraintes de Reynolds en contributions
Mots cles :
The two-phase flows involve interfacial interactions which modify significantly the structure of
the mean and fluctuant flow fields. The design of the two-fluid models adapted to industrial flows
requires the taking into account of the effect of these interactions in the closure relations adopted.
The work developed in this thesis concerns the development of first order two-fluid models
deduced by reduction of second order closures. The adopted reasoning, based on the principle of
decomposition of the Reynolds stress tensor into two statistically independent contributions
turbulent and pseudo-turbulent parts, allows to preserve the physical contents of the second order
relations closure. Analysis of the turbulence structure in two basic flows : homogeneous bubbly
flows uniform and with a constant shear allows to deduce a formulation of the two-phase turbulent
viscosity involving the characteristic scales of bubbly turbulence, as well as an analytical
description of modification of the homogeneous turbulence structure induced by the bubbles
presence. The eulerian two-fluid model was then generalized with the case of the inhomogeneous
flows with low void fractions. The numerical results obtained by the application of this model
integrated in the computer code MELODIF in the case of free sheared turbulent bubbly flow of
wake showed a satisfactory agreement with the experimental data and made it possible to analyze
the modification of the characteristic scales of such flow by the interfacial interactions. The two-
fluid first order model is generalized finally with the case of high void fractions bubbly flows where
the hydrodynamic interactions between the bubbles are not negligible any more.
Key words :
Le travail presente dans ce memoire a ete mene dans le cadre d'une these en cotutelle au
Laboratoire d'Hydraulique de l'Ecole Nationale d'Ingenieurs de Tunis et a l'Institut National de
Mecanique des Fluides de Toulouse de l'Institut National Polytechnique de Toulouse.
Il m'est agreable au terme de ce travail d'exprimer ma gratitude et mes vifs remerciements envers
tous mes collegues et amis a l'ENIT et l'IMFT qui, par leur aide, leur soutien et leurs conseils, ont
contribue a l'aboutissement de ce travail.
Je remercie tout particulierement Monsieur Olivier Simonin, Directeur de l'IMFT, pour le grand
interet qu'il a toujours manifeste envers ce travail, pour son aide precieuse et les conseils qu'il m'a
prodigues. Je le remercie egalement pour l'honneur qu'il m'a accorde en acceptant de presider le
jury.
Je remercie vivement Monsieur Alain Line, Professeur a l'Institut National des Sciences Appliquees
de Toulouse, d'avoir accepte juger ce document en tant que rapporteur et je lui suis tres
reconnaissant pour le grand interet qu'il porte a ce travail.
Je tiens a adresser mes sinceres remerciements a Monsieur Taieb Lili, Professeur a la Faculte des
sciences de Tunis, d'avoir bien voulu porter un jugement pour cette these en tant que rapporteur.
Je me rappelle maintenant en redigeant ces lignes, et apres plusieurs annees, que c'est avec le cours
de la mecanique des milieux continus du Professeur Taieb Lili, que j'ai commence a decouvrir
l'aspect passionnant et fascinant de la mecanique des fluides. Qu'il trouve dans ces quelques lignes
la reconnaissance de l'un de ces anciens etudiants.
Que pourrais je dire pour remercier mes deux directeurs de recherche Monsieur Lucien Masbernat,
Professeur a l'IMFT, et Monsieur Jamel Chahed, Maitre de conferences a l'ENIT. Le travail que j'ai
accompli sous leur direction durant toutes ces annees a ete pour moi un veritable epanouissement
aussi bien sur le plan scientifique qu'humain. L'accueil chaleureux que me reserve a chaque fois
Monsieur Lucien Masbernat m'a permis de me sentir, au cours de mes sejours a Toulouse, un peu
comme chez moi. Son regard critique, son appui scientifique et moral ont ete pour moi d'une aide
precieuse. Je n’oublierai jamais les prodigieux debats que nous avons menees. Qu’il trouve ici
l’expression de toute ma reconnaissance. Les longues heures de discussion passees avec Monsieur
Jamel Chahed ont ete un veritable apprentissage du metier de chercheur. Sa fagon d'aborder et de
simplifier les difficultes que nous rencontrons, son devouement pour la recherche ajoutes a ses
qualites humaines n'ont pu que m'apporter le courage et la sollicitude dont j'avais besoin. Je tiens a
lui exprimer ma profonde gratitude.
Je tiens aussi a adresser mes chaleureux remerciements a mes amis : Monaem, Bechir, Foued,
Nizar, Moez, Sadok et Daouser pour leur precieux soutien. Je suis persuade que nos longs debats
scientifiques et divers ont ete pour moi une vraie source de richesse d'esprit. Je ne me dois pas de
manquer d'avoir une vive pensee a mes amis de Toulouse : Chouaeb Labiod et Abdel Majid
Lemaissi avec lesquels j'ai partage de chaleureux moments.
Enfin je tiens a remercier profondement ma famille pour l'enorme amour qu'elle m'eprouve. Mon
cher pere, qui n'est plus des notres, etait et le restera pour toujours mon honorable exemple dans la
vie. Ma mere, qui ne cesse de m'offrir tant d'affection et de devouement, a su m'apporter tout le
courage dont j'avais besoin. J'espere que ce travail sera pour elle la preuve de mon amour. Une
tendre pensee a ma sreur et mon frere pour leur soutien et leur affection, ils tiennent sans doute
une enorme place dans mon coeur. Un tres grand merci pour mon beau frere Mohsen qui n’a cesse
de me soutenir et surtout je ne me dois pas d’oublier mes deux petits neveux Myriam et Aziz qui
bien qu’ils bien qu’ils soient nes au cours de ce travail, ont contribue a son aboutissement par leurs
sourires et la joie de vie qu’ils repandent. J’adresse egalement ma sincere gratitude a ma belle-
famille. Les cordiaux moments que nous partageons ensemble sont toujours un vrai reconfort pour
moi. Et comme le veut bien la coutume, je finis par exprimer ma profonde reconnaissance a ma
fiancee Rim. Je reconnais qu’elle a supporte, et qu’elle continue de le faire, ma consecration pour
ce travail. J’espere que son aboutissement lui apporte une preuve d’amour.
Sommaire 1
Introduction 5
1.1 Introduction.................................................................................................................................. 9
1.4.1 Introduction....................................................................................................................... 37
1.6 Conclusion.................................................................................................................................. 47
1
Chapitre 2 : Modeles de turbulence en ecoulement a bulles 49
2.1 Introduction................................................................................................................................ 49
Chapitre 3 : Modeles de turbulence a trois equations pour les ecoulements a bulles 108
3.1 Introduction...............................................................................................................................108
3.3 Modeles de turbulence au premier ordre pour les ecoulements a bulles......................... 122
3.4 Test des modeles de turbulence a trois equations en turbulence homogene a bulles......128
2
Chapitre 4 : Couches cisaillees minces 2-D en ecoulement a bulles 141
4. 1 Introduction.............................................................................................................................. 141
4.2.2 Fermetures du modele a deux fluides : Modele ENIT avec fermeture au premier
ordre de la turbulence.....................................................................................................................143
4.5 Conclusion.................................................................................................................................176
5.3 Turbulence induites par les bulles enecoulement a fort taux de vide................................ 187
3
5.3.2 Resultats et discussions................................................................................................... 190
5.4 Un nouveau modele pour la turbulence induite par les bulles :......................................... 191
Conclusion 195
4
Introduction
Les systemes de fluides industriels et environnementaux mettent souvent en jeu des ecoulements
turbulents multiconstituants et/ou multiphasiques (genie petrolier, genie chimique, genie
energetique, geo-fluides etc.). D’une maniere generale la conception et le controle des systemes de
fluides multiphasiques se basent sur des methodes globales de bilans matieres dans les reacteurs
dont les fermetures s’appuient sur des modelisations conceptuelles ajustees a partir d’experiences
realisees en laboratoires ou pratiquees sur des pilotes industriels ou semi-industriels. Cette
approche est a l’origine des differents modeles industriels developpes pour representer de fagon
predictive le fonctionnement global des reacteurs utilises en genie des procedes.
Au plan scientifique, la modelisation des systemes gaz-liquide pose des questions importantes
sur lesquelles se penchent les communautes de chercheurs en mecanique des fluides appliquee et en
genie des procedes. A la base de ces questions se trouve la complexite propre des systemes mis en
reuvre et qui sont souvent congus pour mettre en contact differentes phases afin de favoriser des
phenomenes de transfert ou de transformation. Ces milieux multiphasiques sont complexes et il est
necessaire d’ameliorer la connaissance sur les interactions aux interfaces pour pourvoir representer
les phenomenes de transfert et les cinetiques de transformations dans ces milieux fluides
multiphasiques.
Le progres important realise au cours des dernieres decennies dans le domaine de la Mecanique
des Fluides Numerique permet aujourd’hui une approche plus phenomenologique des problemes de
transferts dans les systemes de fluides mono et multiphasiques et le developpement de modeles
locaux recouvre aujourd’hui des enjeux majeurs au plan des applications. C’est particulierement
vrai pour ce qui concerne les contacteurs gaz-liquide ou des travaux experimentaux et de
modelisation ont permis de progresser dans l’analyse a l’echelle locale des phenomenes physiques
associes aux processus de transfert notamment en ce qui concerne le transfert de masse.
5
elaborees capables de representer ces interactions et leurs effets sur les champs moyens, turbulents
et de taux de presence des phases ; donc sur les mecanismes de transfert et de transformation. Elles
renvoient egalement a la necessite de developper des experiences pour analyser les phenomenes de
transfert aux interfaces et pour developper et ajuster les modelisations.
Plusieurs experiences, realisees au cours de la derniere decade, ont etabli des descriptions locales
d’ecoulements gaz-liquide disperses dans des configurations de base (ecoulement homogenes,
ecoulements cisailles minces, ecoulements en presence de paroi). Ces resultats experimentaux
mettent en evidence l'important effet de la presence des phases sur la structure des champs moyens
et fluctuants. Ils montrent egalement que le champ moyen des forces interfaciales et la turbulence
du liquide controlent la distribution spatiale des phases mettant en evidence les forts couplages
entre champs moyens, champs fluctuant et champs de taux de presence des deux phases. Les
experiences recentes en ecoulements gaz-liquide a forts taux de vide montrent d'importantes
modifications de la structure des ecoulements diphasiques. Dans ces ecoulements, l’agitation
induite par les bulles devient preponderante et la structure de l’ecoulement diphasique est
essentiellement gouvernee par le mouvement relatif des bulles. Il devient donc difficile de bien
representer l’ecoulement diphasique et ses proprietes de transport turbulent sans une description
precise des mecanismes qui controlent la turbulence "classique" du liquide, la "pseudo-turbulence"
induite par les bulles et leurs eventuels couplages.
Parallelement a ces travaux experimentaux, des progres importants dans la modelisation des
ecoulements a bulles, auxquels ont contribue les travaux conduits au Laboratoire de Modelisation
en Hydrauliques et Environnement (LMHE) ont permis de developper des modeles euleriens a deux
6
fluides dans une approche statistique avec fermeture en un point au second ordre de la turbulence.
Les travaux menes au Laboratoire de Modelisation en Hydraulique et Environnement, dans le
domaine de la modelisation des ecoulements diphasiques gaz-liquide disperses visent, precisement,
l’amelioration des connaissances des interactions interfaciales. Les developpements proposes
precisent plus particulierement deux aspects qui sont ou ignores ou peu pris en consideration dans
les modeles actuels :
- en premier lieu, dans la formulation de la force moyenne interfaciale, on fait figurer, le tenseur de
dispersion issu du terme de masse ajoutee. Le role des correlations des fluctuations de vitesse du
gaz, exprimees en fonction du tenseur de Reynolds du liquide, s’avere determinant pour le calcul de
la distribution de taux de vide dans les ecoulements a bulles inhomogenes
Les modelisations de la turbulence et des interactions interfaciales proposees ont ete testees dans
des configurations de base d'ecoulements paralleles ou quasi-paralleles pour lesquelles on dispose
de donnees experimentales. Le modele a ainsi ete mis en reuvre pour simuler des ecoulements
turbulents homogenes a bulles et des ecoulements de couches cisaillees minces, en turbulence libre
et en presence de parois, en gravite normale et en micro gravite. Ces applications se sont averees
tres utiles pour preciser, analyser et comprendre les mecanismes physiques qui controlent les
interactions interfaciales et leurs effets sur les champs cinematiques moyens et fluctuants des deux
phases ainsi que sur leurs champs moyens de taux presence.
Cette presentation rapide des progres recents dans la modelisation des ecoulements diphasiques
conduit a degager un certain nombre d’orientations ou on est en mesure d’attendre des avancees
significatives dans le perfectionnement et la generalisation des modeles. C'est dans cette perspective
que s'inscrivent les travaux presentes dans cette these. En s'appuyant sur une approche
phenomenologique essentiellement fondee sur l'analyse d'experiences, sur l'analyse d'echelles et sur
l'analyse dimensionnelle on propose de nouvelles idees de modelisation des interactions
interfaciales. Il s’agit d’adapter d’orienter et d’anticiper les progres scientifiques dans l’analyse
physique et dans la modelisation des ecoulements diphasique pour pousser les limites de validite des
7
modeles afin prendre en compte les complexites propres aux systemes gaz-liquide. En particulier on
propose un ensemble de fermetures qui visent deux objectifs essentiels : le premier concerne le
developpement, l'analyse et la validation de modeles de turbulence au premier ordre pour les
ecoulements a bulles. Ces modeles sont obtenus par analyse et reduction des fermetures au second
ordre ou on prend en compte les effets interfaciaux dans la modelisation des termes de
redistribution et de diffusion avec des echelles caracteristiques de la turbulence et de la pseudo
turbulence. Le second objectif est relatif a ^amelioration de la modelisation de la pseudo-turbulence
en ecoulements a forts taux de vide pour elargir le domaine d'application des modeles aux systemes
industriels souvent congus pour provoquer des transferts souvent en milieu dense.
propose dans le chapitre 3 au cas des ecoulements inhomogenes. Ce modele sera implante dans le
code industriel MELODIF et applique a l'ecoulement a bulles cisaille mince de sillage. Nous
generalisons en seconde etape le modele au premier ordre developpe au cas des ecoulements a
forts taux de vide ou les interactions hydrodynamiques entre les bulles ne sont plus negligeables.
8
C%4fD%E 7 ELEMENTS' DE EI^EE DE MDDEETEX^D^ EDEEE/EE/^E DEE ECDDEEMEATE ^ EDEEEE
1.1 Introduction
Le travail presente dans cette these s'inscrit dans le cadre des recherches developpees au
Laboratoire de Modelisation en Hydraulique et Environnement a l'ENIT dans le domaine de la
modelisation eulerienne des ecoulements turbulents multiphasiques. Ces travaux conduits en etroite
collaboration avec l'Institut de Mecanique des Fluides de Toulouse visent l'amelioration des
fermetures des modeles euleriens a deux fluides pour les ecoulements gaz-liquide disperses.
Dans les ecoulements diphasiques, les grandeurs physiques subissent des sauts de discontinuite
au passage des interfaces. La description eulerienne de l'ecoulement n'est possible qu'a l'aide du
formalisme des distributions qui permet ainsi d'etendre le domaine de definition des grandeurs
physiques a tout le domaine de l'ecoulement. Les equations instantanees de bilans de masse et de
quantite de mouvement dans chaque phase sont ainsi ponderees par une fonction caracteristique de
presence des phases. Le traitement statistique de ces equations instantanees permet de decrire le
milieu diphasique, en termes de grandeurs moyennes phasiques definies et continues sur tout le
domaine de l’ecoulement.
Le premier chapitre de la these est consacre a la formulation eulerienne des equations de base du
modele a deux fluides et a leur modelisation. Les equations locales et instantanees sont presentees a
l'aide du formalisme des distributions. Ces equations sont ensuite moyennees dans chaque phase et
les bilans moyens qui en decoulent sont commentes vis-a-vis du probleme de fermeture. La
fermeture des termes interfaciaux s'appuie sur la formulation lagrangienne de l'expression de la
9
CHAP/TEE 7 EEEMEWTE DE EASE DE EA MODEETSX^O^ E^EEE/EE/^E DEE ECOEEEMEATE A E^EEEE
force instantanee exercee par le fluide porteur sur les bulles. Plusieurs approches ont ete proposees
pour permettre le passage de la description lagrangienne de la force a la formulation eulerienne de
l'echange interfacial.
Le developpement de fermetures pour le modele eulerien a deux fluides renvoie a la necessite de
les confronter aux donnees experimentales. L'ensemble des resultats experimentaux disponibles
pour les ecoulements a bulles est analyse dans la perspective de degager une classification
hierarchique permettant de tester les modelisations. Nous evoquons egalement l'apport de la
simulation numerique directe (DNS) et la simulation des grandes echelles (LES) dans la
comprehension des phenomenes physiques qui controlent les ecoulements a bulles.
Un ecoulement a bulles peut etre decrit a un instant donne par une repartition aleatoire des
inclusions dans une phase continue. Cette repartition varie au cours du temps et on rencontre en
chaque point du domaine de l'ecoulement diphasique successivement, sur des intervalles de temps
variables, la phase continue ou la phase dispersee. La description eulerienne instantanee d'un tel
ecoulement peut etre realisee en definissant une fonction ck (x, t) qui indique la presence ou non de
10
CKFP/TRE 7 ELEMENTS' DE EI^EE DE E^ MODEETEX^O^ EDEEE/ENNE DEE ECODEEMENTE ^ EDEEEE
On supposera que l'interface separant les deux phases est d’epaisseur nulle de sorte que l’on passe
instantanement en la traversant d’une phase a l’autre. La fonction caracteristique de presence (1-1)
verifie les regles de derivation suivantes (Ishii, 1975) :
O)
1
(1-2)
<
II
»k = - U'V*' (1-3)
Avec :
n Ik
U Ik
Consequence directe de l’hypothese de l’epaisseur nulle des interfaces, la relation (1-2) indique que
le gradient de la fonction caracteristique n’est different de zero qu’au passage des interfaces. La
relation (1-3) signifie que lorsqu’on se deplace a la vitesse uI de l’interface on reste dans la meme
phase.
11
C%4f/TRE 7 ELEMENTS' DE EI^EE DE E^ MODEETEX^O^ EDEEE/ENNE DEE ECODEEMENTE ^ EDEEEE
(1-4)
ii
+
<
o
s
O
—(Pk»k) + V.(pkuku k) = V.(a k) + pkg (1-5)
Ot
ou k = 1,2 pour distinguer la phase continue (k=1) et la phase dispersee (k=2), pk est la masse
volumique, uk est le vecteur vitesse, sk est le tenseur des contraintes et g est l’acceleration de
pesanteur.
Les equations de conservation pour chaque phase definies, au sens des distributions, en tout point
du domaine de l’ecoulement s’obtiennent en multipliant les equations (1-4) et (1-5) par la fonction
Bilan de masse:
O
%-(P k % k) + V-(P k u k % k) = mk (1-6)
Ot
O
%-(P k u k % k ) + V-(P k u k u k % k ) = V-(s k % k) + P k % k g + L k (1-7)
Ot
Les termes interfaciaux font intervenir la distribution de Dirac definie sur les interfaces. Ils
representent les densites volumiques de transfert de masse et de quantite de mouvement a
l'interface.
12
CHAP/TEE 7 EEEMENTE DE EASE DE EA MODEE/EA^ON EDEEETENNE DEE ECODEEMENTE A EDEEEE
Les termes interfaciaux mk et Lk assurent le couplage entre les equations de bilan des deux
phases. En considerant qu'il n'y a pas accumulation de masse aux interfaces et que ces dernieres
sont soumises a des tensions superficielles, on etablit alors les relations de saut suivantes :
m, + m2 = 0 (1-10)
L + L 2 = fsi 51 (1-11)
L'equation (1-10) traduit la conservation de la masse au niveau des interfaces, tandis que le terme
ou H est la courbure algebrique de l’interface selon nI, VI est l’operateur gradient pris par rapport
La resolution directe des equations (1-6) et (1-7) suppose que l’on connaisse a chaque instant la
position des interfaces et leurs vitesses de deplacement ce qui est quasi impossible. Le passage a
une description moyennee de l'ecoulement diphasique suppose que la prise de moyenne doit etre
effectuee dans chacune des phases. Ischii (1975) introduit a cet effet la notion de moyenne
eulerienne phasique definie par :
13
CDXETTEE 7 EEEMENTE DE EAEE DE EA MODEETEA^ON EDEEETENNE DEE ECODEEMENTE A EDEEEE
Notons que l’operateur de moyenne phasique ne commute pas en general avec tout operateur de
derivation D et on montre aisement que :
L'operateur de moyenne statistique verifie les conditions de Reynolds. Il est defini a partir de N
realisations jN dans des experiences repetitives d’une grandeur j en un point (x,t) selon :
1 N
5
< j >= ™(n § ji)
(1-19)
Cependant, l’operateur de moyenne statistique n’est pas adapte pour l'analyse experimental et on a
generalement recours aux operateurs de moyenne temporelle ou spatiale definies respectivement
par :
T dt
<j> (1-20)
<j> — fj dv (1-21)
qui verifient les conditions de Reynolds respectivement dans les ecoulements stationnaires et
ergotiques. Cependant lorsque ces trois operateurs verifient telles conditions, ils permettent
d’aboutir a des equations formellement equivalentes et il est ainsi inutile de preciser la nature de
l’operateur choisi.
14
CHAPTTRE 7 EEEMENTE DE EAEE DE EA MODEE/EAT/ON EDEER/ENNE DEE ECODEEMENTE A EDEEEE
Ck jk X k jk (1-22)
avec :
L’utilisation de l’operateur de moyenne phasique permet ainsi d'obtenir des grandeurs moyennes
phasiques definies et continues sur tout le domaine de l’ecoulement. Elle decrit ainsi le milieu
diphasique comme un melange de deux fluides coexistants et en interaction. Cette formulation
efface l’information sur la geometrie des interfaces et decrit les interactions interfaciales par des
termes definis et continus sur tout le domaine de l’ecoulement.
Les bilans moyens de masse et de quantite de mouvement sont obtenus en moyennant les
equations de conservation locales et instantanees (1-6) et (1-7). Ils s’ecrivent :
e ==
y-(Pka%) + V.(Pkakuk) = < mk > (1-24)
et
Dans l’equation (1-25) le tenseur moyen des contraintes est ecrit comme la somme des contraintes
sk = -pkI+tk (1-26)
t k=-p k uk uk (1-27)
15
CHAP/TRE 7 EEEMENTE DE EAEE DE EA MODEE/EA^ON EDEER/ENNE DEE ECODEEMENTE A EDEEEE
Compte tenu de (1-8) et (1-9), les termes d’echange interfaciaux de masse et de quantite de
mouvement s'ecrivent :
Les interfaces sont supposees sans masse et d’epaisseur nulles dans lesquelles il n’y a pas
accumulation de matiere. La conservation de masse aux interfaces s’ecrit alors :
D’autre part, les interfaces sont soumises a des tensions superficielles. Le bilan de quantite de
mouvement sur ces interfaces s’ecrit :
En l’absence de transfert de masse, les equations de bilans moyens formulees dans le paragraphe
precedent contiennent 14 inconnues :
16
CHAP/TRE 7 EEEMENTE DE EAEE DE EA MODEE/EA^ON EDEER/ENNE DEE ECODEEMENTE A EDEEEE
La fermeture du systeme d’equations du modele a deux fluides necessite que l’on ait autant
d’equations que d’inconnues. Avec les 4 equations disponibles (deux par phase), il est possible
d’expliciter les relations suivantes :
• Les taux de presence des deux phases sont relies par la relation suivante : a1 + a 2 = 1 pour des
• Le tenseur des contraintes visqueuses tk peut etre explicite par des lois de comportement
Le probleme de fermeture va alors concerner cinq inconnues. Si l’on souhaite obtenir un modele
permettant la description des champs de vitesses moyennes dans les deux phases et leur repartition
dans le domaine de l’ecoulement, il faut que l’on choisisse les vitesses moyennes dans chacune des
phases et le taux de presence de l’une des phases comme inconnues principales. Il n'est donc pas
possible de conserver les deux pressions phasiques comme inconnues principales et le probleme de
fermeture concerne donc :
• La pression dans l’une des phases : il s’agit d’exprimer la pression moyenne dans l’une des
phases en fonction de la pression dans l’autre.
La capacite des equations (1-24) et (1-25) a decrire le comportement physique d’un milieu
diphasique depend de la precision avec laquelle les fermetures vont restituer la perte d’information
relative a la prise de moyenne des equations instantanees.
Dans la suite, nous nous limitons a des ecoulements a bulles sans transfert de masse aux interfaces.
Les fluides consideres sont supposes newtoniens a masses volumiques constantes. Dans ces
conditions, le tenseur des contraintes visqueuses de la phase k s’ecrit :
17
CHAP/TRE 7 EEEMENTE DE EAEE DE EA MODEE/EA^ON EDEER/ENNE DEE ECODEEMENTE A EDEEEE
avec vk est la viscosite cinematique de la phases k. Nous negligeons egalement l’ecoulement dans
les bulles. Pour les ecoulements a bulles millimetriques considerees ici, la vitesse en tout point de la
bulle sera identifiee a la vitesse de son centre d’inertie.
Il n'existe pas de methodes parfaitement deductives pour elaborer des fermetures des modeles a
deux fluides. Comme pour les ecoulements monophasiques turbulents, la procedure de fermeture
consiste en une construction d'approximations coherentes des termes inconnus par des formulations
qui doivent obeir a des criteres dictes par des considerations mathematiques et des principes
generaux de la physique. Drew et Wallis (1992) developpent dans le cas diphasique une
methodologie similaire permettant d'aboutir a la representation la plus complete possible des
interactions interfaciales dont nous resumons les criteres generaux:
1.3.2.1 L'homogeneite
18
CHAP/TRE 7 EEEMENTE DE EAEE DE EA MODEE/EA^ON EDEER/ENNE DEE ECODEEMENTE A EDEEEE
Ce critere consiste au fait que le modele propose ne doit pas reproduire des resultats irrealistes
tels qu'a titre d'exemple une energie negative ou une valeur negative d'un coefficient de correlation
cense etre positif.
Ces principes methodologiques constituent ainsi des conditions necessaries pour tout schema de
fermeture. Leur verification permet d'eviter le developpement de comportements physiquement
inacceptable.
En l’absence de transfer! de masse, les termes interfaciaux dans les bilans moyens de quantite de
mouvement s’ecrivent :
Ils representent la densite volumique en un point donne de la force exercee par le liquide sur les
surfaces des bulles. La modelisation de ces deux termes pose deux difficultes principales :
• La premiere est relative a la formulation, dans le cas general, de la force exercee par le liquide
sur une inclusion isolee : l’expression exacte de cette force n'est connue que dans des cas
asymptotiques d’ecoulements de reference (ecoulement potentiel de fluide parfait (Auton et al
(1988)), ecoulement rampant en regime de Stockes (Gatignol (1983)). L’extrapolation de ces
resultats analytiques hors du domaine lineaire reste encore hypothetique. Mais on assiste
depuis plusieurs annees au developpement des travaux numeriques tel que ceux de Rivero et al
(1991), Magnaudet et al (1995) etc. qui ont permis des avancees importantes.
19
C%4f/TRE 7 ELEMENTS' DE EI^EE DE E^ MODEETEX^O^ EDEEE/ENNE DEE ECODEEMEATE ^ EDEEEE
presente. Ceci souleve des problemes theoriques qui n’ont ete abordes que dans le cadre
d’hypotheses restrictives d’homogeneite et de stationnarite du champ turbulent de la phase
continue (Batchelor (1972), Lance (1986), Simonin et Viollet (1989)).
Dans l’etat actuel des connaissances, la solution analytique donnant lieu a l’expression exacte de
la force subie par une particule spherique plongee dans un liquide en mouvement non uniforme
n’est connue que dans les situations asymptotiques de l’ecoulement rampant dit egalement regime
de Stockes (Gatignol 1983 ; Maxey Riley 1983) ou de l’ecoulement de fluide parfait (Voinov
1973 ; Auton 1988). L’ecoulement peut etre considere comme rampant lorsque :
R„ = <1 (1-34)
VL
Rer est le nombre de Reynolds relatif bati sur la vitesse relative et le rayon de la bulle.
L’ecoulement de fluide parfait correspond a un etat asymptotique ou Rer tend vers l’infini. les
solutions etablies dans ces situations limites se basent sur la decomposition du champ des
contraintes en champ non perturbe rencontre lorsque le volume de la particule est occupe par le
liquide et un champ perturbe engendre par la presence de la particule.
Pour les nombres de Reynolds intermediates, on ne dispose pas de solutions exactes et les
expressions adoptees sont generalement fondees sur l’extrapolation des resultats dans les situations
20
CH4E/TEE 7 ELEMENTS' DE EI^EE DE E4 MDDEETEX^D^ EDEEE/EA/^E DEE ECDDEEMEATE 4 EDEEEE
Nous presentons dans la suite quelques resultats caracteristiques concernant l’ecriture exacte de
la force exercee sur une particule et nous discutons ensuite l’extrapolation de ces resultats dans le
cas general.
Les travaux de Gatignol (1983) sont, a ce propos, les plus remarquables. Ils concernent l’etude
des forces exercees par un fluide en ecoulement instationnaire et non uniforme en regime de
Stockes sur une particule solide. La methode adoptee se base sur une decomposition de
l’ecoulement en ecoulement perturbe et non perturbe par la particule. Elle permet d’ecrire la force
exercee par le fluide sur la particule sous la forme d’une superposition de deux contributions : la
premiere est due au champ hydrodynamique non perturbe qui verifie les conditions aux limites en
l’absence de particule et le second est liee a la perturbation due a la presence de la particule.
Apres linearisation des systemes decrivant ces contributions, Gatignol (1983) etablit l’expression
generale suivante de la force exercee par le fluide sur la particule spherique rigide :
avec
du(v)
Fi = Pu( g) (1-37)
dt
Dans ces expressions les vitesses u(v) et u(s) designent respectivement les moyennes sur la surface
vitesse de la particule et g est l’acceleration de la pesanteur. Une etude des ordres de grandeur a
permis a Gatignol de negliger le terme convectif devant le terme temporel dans la derivee
particulaire et de choisir la derivee materielle suivant la particule :
21
CHAP/TEE 7 EEEMEWTE DE EASE DE EA MODEETSX^O^ E^EEE/EE/^E DEE ECOEEEMEATE A E^EEEE
d_ 8_
+ up • v (1-39)
dt dt
Cette expression fait apparaltre d’une maniere coherente l’origine des forces qui peuvent s’exercer
sur une inclusion. Le terme Fext represente l’action des forces exterieures qui dans notre cas
perturbe, il comprend la force de Tchen (premier terme) qui est une force d’inertie representant la
force que le fluide doit exercer sur une particule fluide equivalente a la particule solide pour qu’elle
du(v)
soit animee de l’acceleration cette force n’existe pas si le fluide est au repos ou en
dt
mouvement uniforme, et le poids du liquide (second terme) qui occupe le volume de la particule
(poussee d’Archimede). Le terme FII represente l’action du champ perturbe par la presence de la
• La force de trainee stationnaire de Stockes (premier terme), cette force est due a la non
uniformite du champ des contraintes autour de l’inclusion. Elle apparalt des lors qu’il y a
mouvement relatif meme s’il n’est pas accelere.
• La force de masse ajoutee (second terme). C’est une force d’inertie de resistance instationnaire.
Elle traduit la reaction du fluide sur la particule suite a une acceleration relative de celle-ci.
• La force de Basset (dernier terme) appelee egalement force d’histoire. C’est une force de
resistance qui apparalt dans le cas d’un mouvement relatif de la particule. L’interpretation
physique de ce terme est delicate, toutefois Landeau et Lifchitz (1959) montrent qu’il a comme
origine la diffusion visqueuse retardee a travers la couche limite de la quantite de mouvement.
Voinov (1973) a montre que la force dynamique exercee par un fluide parfait en ecoulement non
uniforme, instationnaire et irrotationnel sur une particule rigide s’ecrit sous la forme :
fluide et de la particule, u est la vitesse absolue du fluide en l’absence de la particule prise au centre
22
CHAETTEE 7 EEEMEW7S DE EASE DE EA MODEETSAITO^ E^EEETE^^E DEE ECOLEEMEWE A E^EEEE
et up est la vitesse absolue de celle ci. On voit apparaltre dans cette expression le terme de Tchen et
de la masse ajoutee qui s’ecrivent d’une maniere identique dans les cas extremes du nombre de
Reynolds relatif. Auton et al (1988) ont generalise ces resultats au cas d’ecoulement parfait, non
uniforme, instationnaire et rotationnel et ont considere les ecoulements autour d’un cylindre et
d’une sphere. Ils montrent que la contribution du rotationnel a la force exercee par le fluide peut
etre superposee a la contribution du champ irrotationnel et etablissent l’expression suivante de la
force :
Le dernier terme qui apparalt dans cette expression represente l’action du champ rotationnel, il
introduit la force de portance qui est exprimee en fonction de la vorticite de l’ecoulement du fluide.
Cette force est de meme origine que la trainee : la non uniformite du champ des contraintes induit
une force resultante dont la composante tangentielle a l’ecoulement traduit la trainee et la
composante normale traduit la portance.
Remarque :
La nature des derivees qui interviennent au niveau des differentes forces a fait l'objet de
nombreuses discussions. Il est alors interessant d’en rappeler les principaux elements :
• L’etude des ordres de grandeur a permis a Gatignol de negliger le terme convectif u .(Vu)
8u
devant la derivee temporelle —, ce qui lui permet de choisir la derivee materielle suivant la
• un raisonnement similaire a conduit Maxey et Riley (1983) a un resultat semblable. Ils trouvent
des expressions de meme forme que celles obtenues par Gatignol. Toutefois, au niveau du terme
de Tchen, ils aboutissent au choix de la derivee suivant le fluide :
D 8
— = —+ u. V (1-42)
Dt 8t
Quant au terme de masse ajoutee, ils font intervenir la derivee suivant la particule (equation
(1-39)).
23
CHAP/IRE 7 EEEMEWE DE EAEE DE EA MODEETSAZ/O^ ELEERTE^^E DEE ECOEEEMEATE A ELEEEE
• Minier a repris les travaux de Gatignol, il a montre que dans le cadre de l’hypothese des nombres
de Reynolds relatifs faibles, quelque soit le mouvement du fluide, on aura l’identite des deux
operateurs de derivee : -d et ^D. Le terme du champ non perturbe s'ecrit ainsi sous la forme :
Du
Fi Pu D -pug (1-43)
• Auton et al (1988) font apparaltre analytiquement dans le cas d'un ecoulement parfait les deux
Plusieurs travaux experimentaux ont ete realises pour tenter de mesurer la force subie par une
particule spherique solide submergee dans un liquide accelere : chute d’une sphere dans un liquide
au repos (Odar 1966), oscillation dans un liquide au repos (Odar et Hamilton 1964), oscillation du
fluide autour d’une sphere (Scohoneborn 1975). Cependant la difficulte de l’experience ne permet
pas d’evaluer d’une maniere rigoureuse l’effet de l’acceleration ainsi que la distinction entre les
differentes contributions telles que la tralnee stationnaire, la force d’inertie ou la force d’histoire.
Les travaux d’Odar Hamilton (1964) concernent l’etude de la force subie par une sphere plongee
dans un liquide au repos et animee d’un mouvement oscillatoire. La sphere est immergee dans une
cuve carree remplie d’huile de cote 9d (d est le diametre de la sphere), la frequence de son
mouvement est reglable ce qui a permis d’explorer la plage des nombres de Reynolds relatifs :
Ces experiences generalisent les resultats analytiques etablis pour le cas des faibles nombres de
Reynolds en introduisant des coefficients aux differents termes de l’expression de la resultante des
forces dynamiques exercees par le fluide sur la sphere rigide. Cette force s’ecrit sous la forme :
24
CHAP/IRE 7 EEEMEWE DE EAEE DE EA MODEEZEAITOA ELEERTEAAE DEE ECOLEEMEAIE A ELEEEE
Jt 0
Ac (1-45)
l
avec :
1.2
cx - cx0(1 + —) (1-46)
Ac
0.066
c 1.05 (1-47)
(Ac2 + 0.12)
3.12
cH 2.88 + (1-48)
(Ac2 +1)3
Plusieurs travaux ont ete developpes pour l’etude de la force de tralnee exercee sur une inclusion
de forme quelconque par les effets visqueux du fluide environnant. Cette force s’ecrit en fonction
du tenseur des contraintes visqueuses sous la forme exacte suivante :
Fd - |t .n ds (1-49)
S
Par analogie a la force de frottement elle sera exprimee a l’aide d’un coefficient de tralnee qui sera
donne dans le cas d’une particule a forme quelconque par la relation :
F
cd 1 (1-50)
2 P Ap u - u (u - up)
Plusieurs modelisations ont ete proposees pour le coefficient de tralnee (Ex: Hadamar et
Rybczinsky (1911) ; Moore (1963)). La modelisation la plus pratique est celle de Wallis (1969).
Elle a generalise la tralnee de Stockes pour des nombres de Reynolds plus grands dans le cas d’une
sphere rigide en introduisant le coefficient empirique CD donne par :
25
CHAPTIRE 7 EEEMEWE DE EAEE DE EA MODEETEAITOA ELEERTEAAE DEE ECOEEEMEAIE A ELEEEE
Bel Fdhila et Duineveld (1996) trouvent une bonne correspondance entre des valeurs
experimentales de la tralnee sur une bulle dans une eau contaminee et une autre expression du
coefficient de tralnee CD pour des spheres rigides a grands nombre de Reynolds (Clift et al 1978) :
Nous signalons enfin que dans des travaux experimentaux recents sur les ecoulements a bulles a
forts taux de vide homogene de Garnier et al (2002) et de zone de melange, les auteurs constatent
une attenuation significative de la vitesse de glissement des bulles qu'ils attribuent a une
modification du coefficient de tralnee CD pour les forts taux de vide. Nous revenons sur cette
question dans le chapitre 5.
La force de portance a egalement fait l'objet d'etudes dans plusieurs travaux. Auton et al (1988) ont
etabli analytiquement l'expression de cette force qui figure dans l'equation (1-41) en ecoulement
potentiel avec un coefficient egal a 0.5. En s'appuyant sur des travaux numeriques Legendre et al
(1998) ont mis en evidence une forte dependance du coefficient de portance en fonction du nombre
de Reynolds relatif. D'autre part, les travaux numeriques de Rivero (1991) ont precise la nature de
la derivee particuliere qui intervient dans le terme de la masse ajoutee. Ils ont egalement confirme la
valeur 0.5 qui reste valable pour une plage de nombre de Reynolds relatif allant de 0.1 a 300. Ces
travaux indiquent aussi que la force de Basset peut etre negligee pour la plage des nombres de
Reynolds relatifs etudiee.
A supposer qu'on dispose d'une formulation exacte de la force exercee par le liquide sur une
inclusion, le probleme est alors d'en deduire une modelisation des termes interfaciaux de transfert
de quantite de mouvement. Il s'agit de definir une methode theoriquement valable permettant de
passer de l'information lagrangienne instantanee a une description macro scopique locale et
moyennee. Nous nous proposons dans ce paragraphe de presenter un bref apergu sur les differentes
approches proposees.
26
CHAPTIRE 7 EEEMEWE DE EAEE DE EA MODEETEAITOA E^EERTEAAE DEE ECOEEEMEAIE A E^EEEE
Methode d'homogeneisation
particules identiques. Un etat instantane de ce systeme sera decrit par le vecteur (Cn,Cn) dans
l'espace des phases forme par l'ensemble des positions des barycentres des inclusions CN et de leurs
1 • •
—f(CN,CN,t) dCN dCN represente la probabilite de trouver a l'instant t les inclusions dans
Cette formulation permet de definir la moyenne d'ensemble diphasique en un point (x, t) d'une
1 r • • •
<A> N! j A(x,t,CN,CN)f(CN,CN,t) dCN dCl (1-54)
La moyenne d'ensemble relative a une grandeur AG(x0,t,CN,CN) liees a la phase dispersee est
1 • • •
27
CHAPTIRE 7 EEEMEWE DE EAEE DE EA MODEETEAITOA E^EERTEAAE DEE ECOEEEMEAIE A E^EEEE
Tandis que la moyenne d'ensemble relative a une grandeur AL(x0,t,CN,CN) dans le liquide est
definie en introduisant la fonction definie sur l'ensemble des inclusions qui prend la valeur 1 quand
• k-N
h(x0,t,CN,CN) -1 H(R x Xk ,t,CN,CN) (1-56)
k-1
ou H est la fonction de heaviside, R est le rayon de l'inclusion et xk les positions des barycentres
des inclusions dans la configuration donnee. Cette moyenne est donnee par :
1 r • • • •
< AL >-N! jAL(xct,CN,CN)h(x,t,CN,CN)f(CN,CN,t/xg)dCN dC N (1-57)
1
Njf(CN,CN,t/x0) dCN dCN -a (1-58)
1 r • • •
N j h(x, t,CN,CN)f(CN,CN,t) dCN d C N - 1 - a (1-59)
f(CN,CN,t/x0) qui ne sont connues que dans des cas tres simplifies (Batchelor 1972, Van
Cette methode qui s'inspire de l'hydrodynamique des fluides (Massignon 1957) a ete proposee
par Lance (1986). Elle considere l'ecoulement diphasique disperse comme etant un fluide dans un
volume V porteur de N inclusions. La description eulerienne macroscopique de l'ecoulement
diphasique sera deduite a partir d'un traitement des equations dynamiques des N particules qui
s'ecrivent :
H11(n)
pG J(n)—— - F(n) avec n - 1,...,N (1-60)
dt
28
CHAPTIRE 7 EEEMEWE DE EAEE DE EA MODEETEAITOA E^EERTEAAE DEE ECOEEEMEAIE A E^EEEE
ou F(n) est la force hydrodynamique exercee par le fluide sur la neme inclusion. Cette force depend
en general des positions et des vitesses des (N-1) autres inclusions ce qui implique une dependance
entre les equations de ce systeme.
La demarche introduite par Lance (1986) consiste a construire un champ de vitesse eulerien et
instantane a partir de ce systeme a N equations lagrangiennes. Pour cela, Lance (1986) definit la
Lance (1986) montre que la fonction p(n) verifie l'equation de conservation suivante
N
%oW) -^p (1-62)
n-1
L'equation (1-63) fournit ainsi une description eulerienne, au point (x,t), de revolution de la
geometrie de l'ecoulement. Lance (1986) etablit egalement l'equation de bilan de quantite de
mouvement qui s'ecrit :
J + V.(P) - — (1-64)
et PG
avec :
NN
P -£p(n)(x,t) uGn)uGn) et F -£p(n)(x,t) F(n) (1-65)
n -1 n -1
29
CHAPTIRE 7 EEEMEWE DE EAEE DE EA MODEETEAITOA E^EERTEAAE DEE ECOEEEMEAIE A E^EEEE
L'equation (1-65), ecrite en variables d'Euler, est equivalente aux N equations lagrangiennes (1-60)
du mouvement. Un traitement statistique des equations (1-63) et (1-64) permet alors de deduire
une description eulerienne en grandeurs moyennes du milieu diphasique disperse.
Cette methode a ete introduite par Simonin et Viollet (1989) et adoptee depuis par plusieurs
auteurs (Ex : Belfdhila et Simonin (1992), Roig (1993), Chahed (1999)). Elle se base sur la
description de l'action de la phase continue sur la phase dispersee en tout point de l'ecoulement par
une densite volumique de force. Cette formulation suppose cependant que la taille des bulles est
suffisamment faible pour que l'hypothese d'homogenei2te a l'echelle de la bulle reste valable.
Les equations euleriennes de bilans instantanes de quantite de mouvement dans les deux phases
s'ecrivent :
% L P L lUf - % LV -S L +% L P L g (1-66)
Du G
%GpG % G fP +% G P G g (1-67)
Dt
D
avec — + u L .V est la derivee particulaire en suivant le mouvement instantane de la phase
Dt
liquide et fP est la densite volumique de la force exercee par le liquide sur la phase dispersee. Cette
densite de force s'exprime en fonction du tenseur des contraintes dans la phase liquide par la
relation suivante :
1 j (V.s L ) dt
fP sL.n ds (1-68)
UB
Si l'on admet l'homogeneite spatiale des moments a l'ordre 1 et 2 a l'echelle du diametre d'une bulle,
nous pouvons identifier les densites volumiques de forces a leurs moyennes volumiques et etablir
par consequent l'egalite suivante :
30
CHAPTIRE 7 EEEMEWE DE EAEE DE EA MODEETEAITO^ E^EERTEAAE DEE ECOEEEMEAIE A E^EEEE
fP - V.s L (1-69)
La densite de force est exprimee comme la somme des contributions de l'ecoulement perturbe et de
l'ecoulement non perturbe par la presence de la phase dispersee, nous ecrivons ainsi :
Soient u|0) et sL0) les champs des vitesses et des contraintes de l'ecoulement non perturbe, les deux
La prise de moyenne des equations instantanees (1-66) et (1-67) conduit aux resultats suivants:
DuT
PLaL « L V.S L -P L V.(a L u Lu L) + a L P L g - M G (1-72)
Dt
Du g (0)
PGaG a 2(V. < s 1 >) -P G V.(a G UG UG) +a G P G g + M G (1-73)
Dt
avec :
qui s'exprime, compte tenu de la formulation lagrangienne de la force qui s'exerce par le liquide sur
une inclusion, comme suit :
3 CT D
M G -< % G V .°L ) > - < % G T P L D lU G - U L)|(U G - U L)) > - < % G P LCA(~ U G - ^ U L^) >
4 L d dt Dt
- < %G2PLCLW X (UG - UD > .(0)^
(1-75)
Nous revenons sur cette approche dans le chapitre 4 dans la presentation du modele a deux fluides.
31
CHAPTIRE 7 EEEMEWE DE EAEE DE EA MODEETEAITOA E^EERTEAAE DEE ECOEEEMEAIE A E^EEEE
Methode basee sur la definition d'une distribution des inclusions dans l'espace des phases
Cette methode decrite par Haynes (2004) est fondee sur la definition a un instant donne d'une
distribution des inclusions dans l'espace des phases. En effet, selon une approche similaire a celle
decrite pour la methode d'homogeneisation, le milieu diphasique disperse est congu comme un
systeme a np particules ou np est le nombre moyen d'inclusions par unite de volume. Une
configuration instantanee de l'ecoulement sera decrite par un vecteur dans l'espace des phases. On
introduit ainsi le nombre probable fp des centres d'inclusions localisees autour d'un point (x, cP) de
l'espace des phases. La valeur moyenne d'une grandeur physique FG associee a la phase dispersee
np -/fpdcp (1-77)
La presence des bulles dont le diametre est de quelques millimetres introduit en general une
echelle de longueur superieure a l’echelle de Kolmogorov et une vitesse relative de deplacement
des bulles superieure aux fluctuations turbulentes de vitesse dans la phase continue (Lance et al,
1996). Il en resulte que les mecanismes de production de transfert et de dissipation de l'energie
turbulente peuvent etre considerablement modifies. Lance et Bataille (1990) dans leur etude des
ecoulements homogenes a bulles, evoquent trois types de mecanismes :
32
CHAPTIRE 7 EEEMEWE DE EAEE DE EA MODEETEAITOA E^EERTEAAE DEE ECOEEEMEAIE A E^EEEE
1. Le premier concerne les fluctuations turbulentes supplementaires generees par les bulles dans
leur mouvement relatif. C'est typiquement le genre de fluctuations qu'on peut enregistrer dans
le cas d’ascension des bulles dans un liquide au repos (Mareuge et Lance, 1995).
2. Le second mecanisme est relatif a la production turbulente par cisaillement au niveau des
sillages des bulles.
3. Le dernier est attribue a la deformation des bulles sous Faction des structures turbulentes du
liquide. Ce mecanisme subsiste meme a faible taux de vide lorsque l’interaction entre les bulles
peut etre negligee. Bien qu’il soit assez peu etudie dans la litterature, ce mecanisme semble etre
important. Serizawa et Kataoka (1990) montrent que la deformation des bulles est susceptible
d’absorber l’energie turbulente du liquide et modifier la dissipation.
Ces mecanismes sont loin d’etre independants, ils interagissent avec le champ de l’ecoulement
moyen et egalement entre eux. Ils introduisent egalement des effets de couplage non-lineaires. La
modelisation des ecoulements turbulents a bulles doit donc tenir compte de ces aspects en tentant
de les introduire dans les modelisations des differents mecanismes de la turbulence.
Les equations de transport des tensions de Reynolds dans le liquide peuvent etre obtenues de
maniere exacte en ecoulement diphasique disperse. En premiere etape, l’equation de transport de la
fluctuation de vitesse est deduite par soustraction de l’equation de la vitesse moyenne de celle de la
vitesse instantanee. En seconde etape on multiplie cette equation par la fluctuation de vitesse et on
ajoute l’equation obtenue a sa transposee et on applique l’operateur de moyenne. On obtient
l’equation de transport des contraintes de Reynolds, Lance (1986).
33
CHAPTIRE 7 EEEMEWE DE EAEE DE EA MODEETEAITOA E^EERTEAAE DEE ECOEEEMEAIE A E^EEEE
D =TT,
Dt(aL ULUL) = "aL u Lu L.Vu l +(u Lu L.Vu l) + - 2na L Vu L.(Vu L)T +vV.(u Lu Ln idI)
(1-78)
Cette equation fait apparaitre des termes equivalents a ceux rencontres dans le cas monophasique et
qui sont ponderes par le taux de presence de la phase liquide. On peut distinguer dans le second
membre respectivement les termes de production par le gradient de vitesse moyenne, de dissipation
visqueuse, de redistribution par la pression, de diffusion et un terme specifique aux ecoulements
diphasiques traduisant le mecanisme de production interfaciale :
Ce terme interfacial peut etre exprime sous la forme de correlations entre les fluctuations de vitesse
et les fluctuations du tenseur des contraintes sur les interfaces :
avec
Cette ecriture permet d’interpreter le terme interfacial comme la puissance developpee par les
fluctuations des contraintes aux interfaces dans le mouvement fluctuant. Dans le cas ou la phase
dispersee est constituee de particules solides, Lance (1986) montre que les conditions
d’impermeabilite et de non glissement a l’interface permettent d’identifier les composantes
normales et tangentielles des vecteurs vitesse dans le liquide et sur la particule (indice p) de sorte
qu’il devient possible d’etablir :
34
C%4f/TRE 7 ELEMENTS' DE EI^EE DE E^ MODEETEX^O^ EDEEE/ENNE DEE ECODEEMENTE ^ EDEEEE
ou M est la contrainte exercee par le liquide sur la phase dispersee a l’interface. Dans le cas ou la
phase dispersee est gazeuse, les conditions a l’interface deviennent des conditions d’impermeabilite
et de nullite des contraintes tangentielles et le resultat (1-82) n’est plus valable que pour la trace du
tenseur.
avec P0 - P : I est la trace du tenseur P. Le terme de production interfaciale dans l’equation (1-83)
Ce resultat montre que la production interfaciale de l’energie turbulente correspond d’une part a la
puissance developpee par la force d’interaction interfaciale dans le mouvement relatif moyen et
d’autre part a la puissance des fluctuations des contraintes associees aux fluctuations de vitesse des
interfaces.
Plusieurs travaux de modelisation ont permis de developper des modeles euleriens a deux fluides
avec fermeture en un point au premier et au second ordre de la turbulence. Ces differents modeles
seront presentes dans le second chapitre et seront commentes vis-a-vis a leurs aptitudes a rendre
compte des observations experimentales dans les situations simples de turbulence homogene en
ecoulement uniforme et avec cisaillement constant. Cette analyse nous conduit a proposer une
forme de classification des modeles de turbulence proposes pour les ecoulements a bulles et d'en
deduire le niveau minimum de fermeture necessaire pour prendre en compte les effets interfaciaux
sur la structure de la turbulence en ecoulement a bulles.
35
CHAP/TEE 7 EEEMENTE DE EASE DE EA MODEE/EA^ON EDEEETENNE DEE ECODEEMENTE A EDEEEE
Nous avons evoque plus haut que la fermeture du modele a deux fluides impose que l'on
exprime la pression dans l'une des deux phases en fonction de la pression dans l'autre. Plusieurs
travaux (Drew et Wood 1985 ; Lopez de Bertodano et al 1990 ; Lahey et Lopez de Bertodano
1991 ; Antal et al 1991), ont propose d'exprimer la difference de pression entre les deux phases
sous la forme generale :
ou pkI est la pression interfaciale du cote de la phase k. La difference de pression (pG - pGI) est
pG - pGI (1-86)
La difference des deux pressions interfaciales est creee par les tensions superficielles et de
contraintes de cisaillement a l’interface. Antal et al (1991) expriment cette difference de pression
comme suit :
4a
(pGI pLI) -1 (1-87)
d
ou 11 est le tenseur de cisaillement a l’interface qui necessite une modelisation. On remarque que
pour respecter la coherence tensorielle de l’equation (1-87) il faut que les termes non diagonaux de
1 2
(pLI - pL) p uG uL (1-88)
4
Cette relation exprime un effet de Bernouilli dans la phase liquide caracterise par une decroissance
de la pression lorsqu’on se rapproche dans la phase liquide de l’interface due a l’augmentation de la
vitesse.
36
CDXETTEE 7 EEEMENTE DE EAEE DE EA MODEETEA^ON EDEEETENNE DEE ECODEEMENTE A EDEEEE
Il faut noter que toutes les relations proposees sont deduites a partir de l’etude de l’ecoulement
autour d’une inclusion spherique isolee. Il n’est pour autant pas assure que ces resultats
conviennent a la description physique de la difference des pressions phasiques dans les ecoulements
a bulles.
Plusieurs auteurs (Drew et Lahey 1982 ; Simonin et Viollet 1990 ; Troshko et Hassan 2001)
adoptent l’egalite des pressions moyennes phasiques dans les ecoulements gaz-liquide. Cela
suppose que la pression exercee en un point par un ecoulement disperse ne depend pas de la phase
ce qui est assez acceptable dans le cas des bulles millimetriques supposees indeformables.
1.4.1 Introduction
Plusieurs experiences, realisees au cours des deux dernieres decades, ont etabli des descriptions
locales d’ecoulements gaz-liquide disperses dans des configurations de base :
37
CHAPTTRE 7 EEEMENTE DE EAEE DE EA MODEE/EAT/ON EDEER/ENNE DEE ECODEEMENTE A EDEEEE
Ces resultats experimentaux mettent en evidence l'important effet de la presence des phases sur la
structure des champs moyens et fluctuants. Ils montrent egalement que le champ moyen des forces
interfaciales et la turbulence du liquide controlent la distribution spatiale des phases mettant en
evidence les forts couplages entre champs moyens, champs fluctuant et champs de taux de presence
des deux phases. L’etude de ces couplages est determinante tant pour la comprehension des
mecanismes d'interaction en milieux diphasiques qu’en regard de tout objectif de modelisation.
La figure (1- 2) ci-dessous illustre les couplages pouvant se produire dans un ecoulement
• Les couplages mis en jeu par les termes interfaciaux (1-29) qui apparaissent dans les equations
(1-25) de bilans de quantite de mouvement
• Les couplages entre le mouvement moyen et fluctuant dans chacune des phases mis en jeu par
la production des fluctuations turbulentes par le gradient de vitesse moyenne qui se manifestent
terme de production par cisaillement aT «Vu l + (uL uLVu l)T dans l'equation (1-78)
• Les couplages induits par l'effet de la paroi ainsi que des conditions aux limites.
38
CfMP/7%E 7 Eza/Ems" DE A4&E DE E4 v\4DDEE/E477D# EfDEE/EAWE DEE ECOf/EEtfEATE,4 «D4JEE
T ciines interiaciaux
On peut egalement distinguer des couplages que nous qualifions "d’implicites" et qui se traduisent
par Vinfluence de la presence de la phase dispersee sur les mecanismes de production, de diffusion
turbulente, de redistribution et de dissipation de la turbulence de la phase continue que nous
symbolisons sur la figure (1-3).
Production
A■ :■
: v
Presence de la phase
Redistribution Diffusion
dispersee
: A
A ■
Dissipation
Si l'on considere en particular les couplages explicites et implicites induits par la presence de la
phase dispersee, ils seraient considerablement simplifies si la phase dispersee est uniformement
repartie spatialement. Cette hypothese est realisee dans les experiences d'ecoulements homogenes a
bulles d'ou on deduit le premier critere de classification lie a l'homogeneite de l'ecoulement.
39
C%4fD%E 7 ELEMENTS' DE EI^EE DE MDDEETEX^D^ EDEEE/EE/^E DEE ECDDEEMEATE ^ EDEEEE
Dans les ecoulements inhomogenes la non uniformite (spatiale ou temporelle) de la repartition des
phases et des gradients de vitesses moyennes fait que tous les couplages implicites et explicites
representes dans les figures (1-2) et (1-3) interviennent. Cependant, l'etude des ecoulements
inhomogenes peut egalement s'envisager graduellement en introduisant pour cela un deuxieme
critere de classification lie aux hypotheses de simplification du tenseur de gradient de vitesse
moyenne de la phase continue (Ex : bidimensionnalite de l'ecoulement ; presence ou non de paroi).
40
CfMP/7%E 7 Eza/Ems" DE A4&E DE E4 vt/ODEE/E477D# E[DEE/EAWE DEE ECO[/EEt4EATE,4 RD/JEE
Teimes interfaciaux
Production
f
Presence de la phase
Redistribution Diffusion
disperse*
Y
Dissipation
L'experience de Mareuge et Lance (1995) correspondant a Fascension des bulles a faibles taux
de vide dans un liquide au repos permet de focaliser Fetude sur les aspects caracteristiques de la
production par les interactions interfaciales ainsi que sur la structure du champ fluctuant induit par
les bulles. L'etude de la decroissance de la turbulence a l'aval de la grille en ecoulement a bulles a
faible taux de vide (Lance et Bataille (1991)) permet d'etudier la superposition de la turbulence
induite par les bulles a celle generee par la grille ainsi que F influence des bulles sur le mecanisme de
dissipation. Les travaux de Lance et al (1991) en turbulence homogene a bulles uniformement
cisaille constitue une continuity logique de ce qui precede. En situation monophasique, Fadjonction
d’un gradient de vitesse moyenne du liquide introduit une production par cisaillement de la
41
C%4fD%E 7 ELEMENTS' DE EI^EE DE MDDEETEX^D^ EDEEE/EE/^E DEE ECDDEEMEATE ^ EDEEEE
Les ecoulements a bulles turbulents libres sont caracterises par la non uniformite des tenseurs de
gradient de vitesse moyenne dans les deux phases et de la repartition du taux de vide. Ils se
developpent selon une direction privilegiee sans etre soumis a aucune contrainte de paroi. Ces
ecoulements introduisent par rapport aux ecoulements homogenes les differents couplages entre les
mouvements moyens, turbulents et la repartition des phases montres par les figures (1-2) et (1-3)
tout en evitant les difficultes liees a la presence de paroi.
Les ecoulements a bulles de paroi peuvent egalement etre etudies graduellement selon un sens
croissant de complexite en partant de l'ecoulement de couche limite a bulles jusqu'aux ecoulements
a bulles a geometries complexes. La configuration d'ecoulement de conduite a bulles en
microgravite (Kamp et al 1995) est particulierement interessante dans la mesure ou elle elimine les
transferts interfaciaux de quantite de mouvement puisque la vitesse de glissement est nulle. Dans ce
cas les schemas des couplages explicites et implicites se reduisent ainsi :
42
CfMP/TRE 7 Eza/Ems" DE A4&E DE M Ef/IEMEAWE DES" ECOf/ZafEATS",4 mt/AAES"
Tonnes interfaciaux
Production
Presence de la phase
Redistribution Diffusion
►
diqiersee
<
*
■
Dissipation
Nous repertorions dans le tableau (1-1) les experiences realisees dans les diverses configurations
d'ecoulements turbulents inhomogenes a bulles.
43
CK4P/TRE 7 ELEMENTS' DE EI^EE DE ET MODEETET^O^ EDEEE/ENNE DEE ECODEEMENTE T EDEEEE
Moursali et al (1995).
Couche limite a bulles
Gabillet et al (1998)
Serizawa et al (1975)
Michiyoshi et Serizawa
Ecoulement a bulles (1986)
turbulent de paroi
Wang et al (1987)
Ecoulement a bulles en conduite
Liu et Bankoff (1987)
Serizawa et al (1992)
Grossetete (1995)
44
CK4P/TRE 7 ELEMENTS' DE EI^EE DE E4 MODEETE4^O^ EDEEE/ENNE DEE ECODEEMENTE 4 EDEEEE
D'un point de vue physique, les ecoulements turbulents diphasiques mettent en jeu des echelles
supplementaires associees aux caracteristiques de la phase dispersee. La pertinence des methodes
de simulation depend du niveau de description qu'elles permettent et donc des echelles
caracteristiques qu'elles reproduisent. Nous nous proposons dans ce paragraphe de donner un
aperqu sur les differentes classes des methodes de simulation numerique des ecoulements a bulles.
D'un point de vue theorique, la simulation numerique directe constitue la meilleure reproduction
des ecoulements turbulents a bulles puisqu'elle ne necessite aucune modelisation explicite de la
turbulence. Elle consiste a la resolution des equations instantanees de bilans de masse et de quantite
de mouvement (de Navier Stockes) dans chacune des phases avec la prise en compte des conditions
aux limites aux interfaces qui meuvent dans le domaine du calcul numerique de l'ecoulement (Hirt
et Nichols 1981; Tomiyama et al. 1994; Susuman et al. 1994; Unverdi et Triggvason 1992; Takagi
et Matsumoto 1997).
Cette approche reste cependant limitee aux faibles nombres de Reynolds et a un nombre reduit de
bulles a cause de son cout excessif sur le plan du calcul numerique. En effet, a grand nombre de
Reynolds la reproduction des echelles caracteristiques des fluctuations turbulentes au niveau des
sillages des bulles necessite une tres grande resolution spatiale.
Une deuxieme methode est egalement utilisee dans le cadre de la simulation numerique directe
des ecoulements a bulles selon laquelle on neglige l'ecoulement a l'interieur des bulles. Ces dernieres
sont supposees indeformables dont les trajectoires sont calculees a partir des equations du
mouvement qui s'ecrivent pour N bulles de volume JP et de densite pP (P = 1,..., N) comme suit :
dx P
ut (1-89)
"dT
du p
p p J p "dT Fp (1-90)
45
CHAP/TEE 7 EEEMENTE DE EASE DE EA MODEE/EA^ON EDEEETENNE DEE ECODEEMENTE A EDEEEE
exercee par le liquide sur celle-ci. Les bulles sont ainsi modelisees par des points materiels disperses
sources de quantite de mouvement. Leur effet sur l'ecoulement de la phase liquide est pris en
compte par l'introduction dans les equations instantanees monophasiques de termes sources sous la
forme d'une distribution representant la force exercee par les inclusions sur le liquide (Squires et
Eaton (1990); Elgobashi et Truesdel (1993); Boivin et al. (2000); Sandaram et Collins (1999);
Miller et Bellan (1999); Ahmed et Elgobashi (2000)). Ces equations s'ecrivent dans le cas d'un
fluide incompressible comme suit :
V.uL = 0 (1-91)
du
PL L + P L u LVu L VpL +mV-(VuL) + IL (1-92)
dt
viscosite dynamique. IL est un terme source sous forme de distribution de type Dirac introduisant
la force par unite de volume exercee par les inclusions sur le liquide. Il est formule comme suit :
N r
1 l (xt) = "Zjf (x p, t)H(x - y)dy (1-93)
p=1
Dans cette expression, f(xP, t) represente la force exercee sur l'inclusion P dont le centre est situe
Cette methode constitue un outil assez performant pour l'etude des interactions des inclusions
avec la turbulence de la phase continue (Fevrier et Simonin (2001)). Cependant, elle reste limitee
par la difficulte de sa mise en reuvre des que la configuration d'ecoulement envisagee devient
complexe ou le nombre des bulles augmente. Un autre inconvenient de cette methode concerne
l'expression de la force exercee par le liquide sur une bulle qui n'est pas etablie sous une forme
analytique dans le cas general. Cette difficulte a ete contournee dans certains travaux recents en
calculant a partir de l'ecoulement de la phase continue simule le tenseur des contraintes exerce a la
surface de l'inclusion (Fadlun et al. (2000); Hu et al. (2001); Burton et Eaton (2002)).
46
CDXETTEE 7 EEEMENTE DE EAEE DE EA MODEETEA^ON EDEEETENNE DEE ECODEEMENTE A EDEEEE
Dans cette approche, les equations instantanees de Navier Stockes sont filtrees spatialement afin
d'obtenir des equations de transport permettant de calculer les grandes echelles de la turbulence de
la phase continue. Tandis que les micro-echelles a caractere universel sont modelisees par un
modele dit de sous-maille. Les trajectoires des inclusions sont calculees de la meme maniere que
dans le cas de la methode de simulation directe.
Par comparaison a la DNS, la LES induit une perte d'information au niveau de la modelisation des
micro-structures de la turbulence, cependant elle conserve un niveau de description assez
satisfaisant tout en permettant de depasser la limite des faibles nombres de Reynolds rencontrees
dans la DNS. Cette methode a permis en particulier d'analyser le phenomene de dispersion des
inclusions dans une turbulence homogene et isotrope (Fevrier et Simonin, 2001) et constitue un
outil important pour tester la validite des relations de fermeture plus ou moins empiriques adoptees
dans les modeles a deux fluides. Nous pouvons citer a titre d'exemples le developpement
d'approches numeriques par la technique de LES a permis des progres significatifs en precisant
certains aspects des interactions interfaciales en ecoulements a bulles (Rivero et al (1991) ;
Magnaudet et al (1995) ; Sangani et Didwania (1992), Merle et al (2004)). Cependant la prise en
compte de l'effet des bulles dans les modeles de sous-maille reste encore assez approximative.
1.6 Conclusion
47
CHAPTTRE 7 EEEMENTE DE EAEE DE EA MODEE/EAT/ON EDEER/ENNE DEE ECODEEMENTE A EDEEEE
Nous avons vu que le probleme de fermeture concerne essentiellement les termes interfaciaux de
transfert de quantite de mouvement et du tenseur de Reynolds dans la phase liquide. La fermeture
des termes interfaciaux se heurte a deux difficultes principales : la premiere concerne la formulation
de la force qui s'exerce sur une inclusion isolee. La deuxieme difficulte est associee au passage de
l'expression de la force de nature lagrangienne a une formulation eulerienne locale de l'echange
interfacial. Nous avons egalement souligne que la fermeture des termes interfaciaux fait apparaitre
le terme turbulent de double correlation des fluctuations des vitesses dans le liquide. Ce resultat
montre la grande importance de la precision de la fermeture de la turbulence dans le liquide. Nous
nous proposons de developper la modelisation de la turbulence dans le liquide dans le chapitre
suivant.
48
CHAP/TRE 2 MODEEEE DE TDREDEENCE EN ECODEEMENTE A EDEEEE
2.1 Introduction
Les fermetures de la turbulence dans les ecoulements a bulles se sont largement inspirees de
l'experience acquise dans l'etude des ecoulements turbulents monophasiques. D'une fagon generale,
les methodes de fermeture en un point de la turbulence donnent lieu a toute une classe de modeles
qui sont largement utilises dans la predetermination des ecoulements turbulents. Issus de la
hierarchie des equations des moments en un point, ces modeles sont d'autant plus performants que
leur ordre est eleve ; en revanche les modeles a degres eleves sont les plus complexes.
• D’abord au niveau de la formulation des fermetures : l’analyse physique des experiences permet
la comprehension des phenomenes qui est souvent a la base des schemas conceptuels de
modelisation.
• Ensuite dans l'ajustement des modeles : les experiences sont utilisees pour evaluer et ajuster les
modelisations
Les resultats obtenus dans les experiences menees sur les ecoulements diphasiques disperses
montrent que la repartition des phases est fortement correlee aux champs fluctuants dans ces
ecoulements : On sait depuis plus de deux decennies que la turbulence dans le liquide joue un grand
role dans la distribution des phases dans les ecoulements a bulles (Drew et Lahey 1982) et des
experiences recentes montrent l'important effet des contributions turbulentes des termes
49
CHAP/TRE 2 MODEEEE DE TDREDEENCE EN ECODEEMENTE A EDEEEE
Les developpements proposes ont permis des progres significatifs dans Amelioration des
modelisations de la turbulence dans les ecoulements a bulles. Cependant, on ne trouve pas dans la
litterature specialisee d'analyse de classification ni d'evaluation systematique des modelisations.
Pourtant l'examen des modeles proposes indique qu'une classification est possible en considerant le
niveau de fermeture et le nombre d'echelles mises en jeu. L'analyse avec classification des modeles
de turbulence en ecoulement a bulles va s'averer tres utile dans la mesure ou elles met en evidence
l'evolution des methodes de modelisation qui, en quelque sorte, a accompagne l'evolution de la
50
CHAP/TRE 2 MODEEEE DE TDREDEENCE EN ECODEEMENTE A EDEEEE
Nous nous proposons dans ce chapitre d'analyser les modeles de turbulence developpes pour les
ecoulements a bulles en essayant a chaque fois d'evaluer les fermetures vis-a-vis des donnees
experimentales de la turbulence homogene a bulles en ecoulement uniforme, (lance et Bataille
1991) et de cisaillement pur, (Lance et al 1991). Ces experiences marquantes ont permis de preciser
certains aspects essentiels des interactions interfaciales et de leurs effets sur la structure de la
turbulence dans des situations d'ecoulements simples ou il est possible d'interpreter et d'evaluer les
fermetures proposees.
Nous avons evoque au chapitre 1, que l'ecoulement de turbulence homogene avec et sans
cisaillement sont des situations qui permettent un decouplage entre les mouvements moyen et
fluctuant dans le liquide et quelles permettent a ce titre d'interpreter la modification de la structure
de la turbulence en ecoulements a bulles par les seuls effets interfaciaux.
En l'absence de gradients des grandeurs fluctuantes, la diffusion est negligee et les seuls
mecanismes qui subsistent sont :
- La dissipation par les effets visqueux.
- Le transfert spectral de l'energie induit par les interactions tourbillonnaires.
- La production par les gradients moyens (uniforme en ecoulement uniformement cisaille et nulle
en ecoulement uniforme)
- La production interfaciale uniforme dans tout le domaine (taux de vide constant).
Nous nous proposons d'analyser dans un premier temps la turbulence homogene uniforme et on
examinera ensuite l'ecoulement uniformement cisaille. La comparaison des donnees en ecoulements
a bulles avec les donnees de l'ecoulement monophasique nous permet de preciser les mecanismes
interfaciaux et leurs effets sur la structure de la turbulence.
51
CHAP/TRE 2 MODEEEE DE TDREDEENCE EN ECODEEMENTE A EDEEEE
dk
= e0 (2-1)
dt
Le transfert spectral de l'energie turbulente associe a cette evolution peut etre decrit par l'equation
suivante (Chassaing 2000) :
dE
T - 2nk2E (2-2)
"dt7
Dans la zone proche de la grille ou les forces d'inertie liees aux fluctuations turbulentes sont
importantes le mecanisme de transfert est dominant de sorte que l'on puisse considerer que les
tourbillons dissipatifs sont dissocies des structures turbulentes energetiques. En revanche, assez loin
en aval de la grille (pour x / M > 250, x etant la position longitudinale et M la taille de la maille de
la grille), les fluctuations turbulentes sont reduites a un point tel que les forces d'inertie deviennent
negligeables devant les forces de viscosite, Batchelor et Townsend (1948). On aboutit ainsi a une
situation asymptotique ou le transfert spectral est domine par le mecanisme de la dissipation
visqueuse.
52
CfMP/TRE 2 EZ/ZZES
respectivement aux structures turbulentes energetiques caracterisees par une echelle de temps xE
et aux structures turbulentes dissipatives caracterisees par une echelle de temps xe, tres faible
devant xE . II en resulte que les structures dissipatives s'adaptent rapidement aux variations des
structures energetiques. Get etat d'equilibre est qualifie par Kolmogorov (1948) d'universel du fait
de son independance des caracteristiques macroscopiques de l'ecoulement. Dans le cas de la
turbulence homogene pleinement developpee, cette zone d'equilibre universel s'elargit suffisamment
pour faire apparaitre une zone dite "sous-domaine inertiel" sur laquelle le mecanisme de transfert
est preponderant par comparaison a la dissipation. II en resulte que le taux de dissipation
correspond dans cette zone au flux conservatif de la cascade d'energie vers les petites structures
turbulentes et que la densite spectrale E est independante de la viscosite cinematique v. Une
analyse dimensionnelle permet alors d'etablir la relation suivante dite loi en puissance - ^.
2 --
3 (2-4)
Ce resultat theorique confirme par differents travaux experimentaux, met en evidence un caractere
universel de la turbulence monophasique. II indique qu'il existe une plage de nombre d'ondes ou la
53
CfMP/TRE 2 .' EZ/ZZES
Les experiences de Lance et Bataille (1991) de la turbulence homogene diphasique montrent que la
structure du spectre est fondamentalement modifiee par le presence des bulles. Leur analyse
conduit a un certain nombre de conclusions essentielles pour comprendre les mecanismes
fondamentaux par lesquels la turbulence est alteree par la presence des bulles.
Les experiences de Lance et Bataille (1991) sont realisees avec des faibles taux de vide ne
depassant pas 3 % et des bulles de diametre de l'ordre de 5 mm. Ces experiences, decrites en detail
dans Lance et Bataille (1991), produisent, en monophasique, des resultats similaires a ceux obtenus
par Comte-Bellot et Corsin (1966) : on retrouve le caractere isotrope de la turbulence et la meme
decroissance de l'energie turbulente a l'aval de la grille.
Les resultats obtenus dans le cas diphasique montrent une modification considerable du spectre de
l'energie turbulente en ecoulement a bulles (Figure (2-2)) : en particular la decroissance spectrale
5 8
en puissance - — est remplacee par une decroissance en puissance - — ; resultat d'ailleurs continue
iN -I
i i I iftn I
54
CK4E/7EE 2 MODELEE DIE mREUZEACE EA ECO^LEMEATE ^ E^LLEE
Lance et Bataille (1991) proposent une interpretation simple mais pertinente des mecanismes a
l'origine de la modification du comportement spectral de la turbulence en ecoulement a bulles. Ils
ecrivent les equations regissant les phenomenes de transport de l'energie turbulente et de son
transfert spectral en ecoulement a bulles sous la forme :
f = Pi-e
(2-5)
— = n- T - 2 Vk 2e (2-6)
dt
Par comparaison aux equations monophasiques (2-1) et (2-2), l'equation (2-5) fait apparaitre un
terme supplementaire associe a la production interfaciale de l'energie turbulente. Dans l'equation (2
6) apparait egalement un terme de production de la turbulence au niveau des sillages de bulles.
visqueuse caracterises respectivement par les echelles de temps tSl, tE , te. En s'appuyant sur les
donnees experimentales, l'analyse des ordres de grandeur montre que les echelles de production-
dissipation sont tres faibles en comparaison avec l'echelle caracteristique du transfert spectral ; cette
comparaison permet d'ecrire :
Ce resultat important indique que les fluctuations turbulentes produites dans les sillages sont
immediatement dissipees avant meme que leur transfert spectral ne puisse avoir lieu. Dans ces
conditions, l'equation (2-6) se reduit a l'equilibre suivant :
0 »n- 2 nk 2E (2-8)
n»-^ (2-9)
k
55
CfMP/TRE 2 EZ/ZZES
E* k 3 (2-10)
g
qui exprime une decroissance en puissance -3 voisine de la valeur experimentale de - —. Ce
resultat confirme, du moins pour les ecoulements a bulles a faibles taux de vide, l'hypothese de
l'equilibre production-dissipation dans les sillages.
Les experiences de Lance et Bataille (1991) indiquent d'autre part que la turbulence en ecoulement
homogene a bulles s'exprime comme la somme d'une turbulence induite par les bulles (uniforme
dans tout le domaine) et d'une turbulence associee a la cascade energetique, figure (2-3). Dans ces
conditions, l'equilibre production-dissipation dans les sillages des bulles implique que la
decroissance de l'energie cinetique a l'aval de la grille se produit avec un temps caracteristique
comparable a celui de l'ecoulement monophasique.
X/M
Figure (2-3) : Decroissance de l’intensite turbulente Vrf^/U derriere une grille, M = 4cm : la
56
CH4E/ZRE 2 MODEEEE DE mRBDZEACE EA ECODLEMEATE 4 BDLLEE
Plusieurs travaux ont ete consacres a l'etude de la turbulence homogene monophasique soumise
a un cisaillement constant. Nous citons a titre d'exemples les travaux experimentaux de Rose (1966,
1970), de Champagne et al (1970) ou de Mulhern et Luxton (1975) ainsi que les travaux
numeriques de Rogallo (1981) dans lesquels cet ecoulement a ete traite par simulation numerique
directe. L'ensemble des resultats obtenus montre que revolution longitudinale des composantes du
tenseur de Reynolds le long de l'ecoulement fait apparaitre deux phases : une premiere phase
transitoire caracterisee par une decroissance des composantes normales du tenseur de Reynolds et
une augmentation du frottement turbulent, suivie par une seconde phase caracterisee par un
equilibre asymptotique ou les composantes du tenseur de Reynolds se stabilisent a un niveau qui
demeure sensiblement constant. Cet equilibre, qui correspond a un regime anisotrope, ne peut etre
entretenu que par la presence de mecanisme d'apport d'energie turbulente par l'ecoulement moyen
afin de compenser la dissipation visqueuse.
L'ecoulement de turbulence homogene avec cisaillement uniforme constitue une etape essentielle
dans l'analyse des mecanismes de la turbulence. Le fait que le champ fluctuant soit homogene nous
autorise a negliger les phenomenes de diffusion et de focaliser l'analyse sur l'equilibre production-
redistribution-dissipation. C'est d'autant plus utile en ecoulement diphasique, que cette situation
d'ecoulement permet d'ecarter le probleme de la distribution des phases (consideree uniforme).
La turbulence homogene avec cisaillement uniforme en ecoulement a bulles a ete etudiee par Lance
et al (1991) avec des taux de vide moderes : les resultats experimentaux obtenus ont permis de
preciser les effets des interfaces sur les mecanismes de la turbulence notamment en ce qui concerne
les effets de redistribution. La figure (2-4) montre l'evolution longitudinale des composantes du
tenseur de Reynolds u'2, v'2 et u v' en ecoulement monophasique. On retrouve bien les deux
57
CfMP/TRE 2 EZ/ZZES
phases de decroissance et d'equilibre de l'energie, observees dans des experiences anterieures (Rose
(1966, 1970); Champagne et al. (1970, 1977)).
model
a = 0
SX/Ul
monophasique homogene a cisaillement constant, S = 2.9 s ' (D'apres Lance et al, 1991)
Dans le tableau (2-1) on compare les valeurs des composantes bn, b22 et b12 du tenseur
u'u 2
b ij (2-11)
k 3
58
CHTE/TEE 2 MODELEE DIE mREUZEACE EA ECO^LEMEATE 4 E^LLEE
(s-1)
Ces resultats indiquent que l'ecart entre les composantes normales b11 et b22 du tenseur
d'anisotropie dans le cas monophasique depend du taux de cisaillement. Plus precisement, les
resultats de Rose (1966, 1970) et de Champagne et al. (1970, 1977) montrent que cet ecart
augmente avec le taux de cisaillement (Chassaing 2000). En revanche un resultat essentiel des
experiences de cisaillement uniforme est de montrer que le frottement turbulent norme par l'energie
cinetique turbulente demeure quasiment constant voisin de 0.3.
Dans les experiences diphasiques de turbulence homogene avec cisaillement uniforme, Lance et
al (1991), les bulles sont injectees d'une maniere uniforme a l'entree de l'ecoulement de telle sorte
qu'on realise une repartition uniforme du taux de vide dans tout le domaine de l'ecoulement. Cette
injection ne modifie pas le profil de la vitesse moyenne du liquide (figures (2-5) et (2-6)). Le fait
d'avoir le meme cisaillement dans les experiences monophasique et diphasiques est essentiel ; cela
permet une interpretation aisee des mesures des composantes fluctuantes en ecoulement diphasique
: avec le meme taux de cisaillement, la variation du frottement turbulent (equivalente a une
variation de viscosite turbulente) peut par exemple ^interpreter vis-a-vis des effets des interactions
interfaciales sur les mecanismes de la turbulence.
59
CfMP/TRE 2 EZ/ZZES
1.5 -
X M
■ 12.5
□ 30
(cm) (cm)
Les figures (2-7)-(2-9) montrent revolution des composantes du tenseur de Reynolds normees par
la vitesse moyenne du liquide au centre de la section pour les taux de vide de 1 %, 1.4 % et 2 %
avec un taux de cisaillement de 2.9 s'1, La figure (2-10) montre les memes grandeurs avec un taux
de vide de 1 % et un cisaillement de 2.4 s'1. Ces figures indiquent une tendance a l'isotropie en
de la composante v'2 ainsi qu'une reduction du frottement turbulent. Pour mettre en evidence ces
effets, on presente dans le tableau (2-2) les valeurs des composantes du tenseur d'anisotropie a l'etat
d'equilibre en ecoulement monophasique et en ecoulements diphasiques.
60
CfMP/TRE 2 EZ/ZZES
model
— Ul
SX/U L
diphasique homogene a cisaillement constant, S = 2.9 s 1, a = 0.01 (D'apres Lance et al., 1991)
FT10
model
SX/U L
diphasique homogene a cisaillement constant, S = 2.9 s_1, a = 0.014 (D'apres Lance et al., 1991)
61
CfMP/TRE 2 EZ/ZZES
Ul
model
SX/U L
diphasique homogene a cisaillement constant, S = 2.9 s 1, a = 0.02 (D'apres Lance et al., 1991)
model
" 0.5
SX/Ul
diphasique homogene a cisaillement constant, S = 2.4 s 1, a = 0.01 (D'apres Lance et al., 1991)
62
CK4E/ZRE 2 MODELEE DIE mREUZEACE EA ECO^LEMEATE 4 E^LLEE
Tableau (2-2) : Variation des composantes du tenseur d'anisotropie en fonction du taux de vide
(experience de Lance et al. (1991))
Composantes du
tenseur d'anisotropie bn b22 b12
Taux de vide
Pour interpreter ces resultats, revenons aux equations de transport des composantes du tenseur de
Reynolds ; celles-ci s'ecrivent dans le cas d'ecoulement uniformement cisaille sous la forme suivante
(voir paragraphe 1.3.4.1 ; Chapitre 1) :
= -u v S-ejL+Zl
(2-15)
dt 2 2
63
C%4f/TRE 2 MODELEE DE TDREDLEACE EA ECODLEMEATE ^ EDLLEE
la cascade energetique et les equations de transport des composantes du tenseur de Reynolds (2-
12)-(2-14) se reduisent a l'equilibre au systeme suivant :
0 = -2uV S + 2p'°^: - 2 e 0
(2-16)
p Ox 3
0 = 2p^-2«.
(2-17)
p °y 3
^ _ 1 _Ou' Ov\
0 = - v' 2S + - p'(— + —) (2-18)
p Oy Ox
0 = -u 'v' S -e 0 (2-19)
L'experience montre que la structure de la turbulence est considerablement modifiee par la presence
de la phase dispersee ; on note une augmentation de l'isotropie et une diminution du frottement
turbulent. L'isotropisation de l'ecoulement se manifeste par une augmentation de la composante
transversale v'2 ; cette augmentation ne peut s'interpreter que par l'effet des bulles sur les termes
de correlation pression-deformation. Il semble que le passage aleatoire des bulles induit un
64
CH&P/TRE 2 MODELEE DE TDEEDLEACE EA ECODLEMEATE ^ EDLLEE
On se propose d'examiner dans la suite les differents modeles de turbulence proposes pour les
ecoulements a bulles en vue de proposer une forme de classification des fermetures; classification
qui s'appuie essentiellement sur l'ordre de fermeture et sur le nombre des echelles qu'elles mettent
en jeu. Dans cette etude nous commengons d'abord par les modeles qui ne pro cedent pas au
decouplage entre les fluctuations produites par cisaillement et celles induites par les bulles. Dans
cette classe de modeles qu'on appellera dans la suite modeles "a une echelle de temps", les modeles
sont developpes au premier et au second ordre ; ils n'introduisent pas explicitement des echelles
propres a la turbulence induite par les bulles. La seconde classe de modeles appelee "modeles a
deux echelles de temps" represente les fermetures ou l'on introduit des echelles propres au
mouvement des bulles pour representer les effets interfaciaux sur la turbulence du liquide. Ces
effets interfaciaux sont, selon les modeles, introduits dans les fermetures au premier et au second
ordre.
65
CK4ELTRE 2 MODELEE DE TDRBDLEACE EA ECODLEMEATE ^ EDLLEE
— 2
ULUL = -2V tS + 3kLI (2-20)
Plusieurs modeles de turbulence a deux equations a une echelle de temps (modeles k-e
essentiellement) ont ete proposes pour les ecoulements a bulles (Ex : Lee et al (1989), Simonin
(1991), Wang et al (1994), Morel (1997)). Dans ces modeles, l'ensemble des fluctuations
turbulentes y compris celles induites par les bulles est represente par une seule echelle et les effets
interfaciaux sur la structure de la turbulence sont modelises par l'introduction de termes interfaciaux
dans les equations de transport des grandeurs turbulentes. Les equations de transport de l’energie
1
turbulente k = — u L. u L et de son taux de dissipation e s'ecrivent ainsi sous la forme generale
suivante :
D
(a L e) = V. a L(— + v)Ve + a (c e1PL c e2 e) + Se (2-22)
Dt kT
Les equations (2-21) et (2-22) sont analogues aux equations de transport en ecoulement
monophasique a ceci pres qu'elles sont ponderees par le taux de vide a et qu'elle contiennent des
termes supplementaires d'echange interfacial Se et Se. Les termes du second membre representent
Les modeles au premier ordre a une echelle de temps ont ete mis en reuvre pour predeterminer la
structure de la turbulence dans differentes situations experimentales d'ecoulements a bulles y
compris celles ou une attenuation de la turbulence a ete observee (Ex : Wang et al 1987, Liu et
Bankoff 1990). Ces modeles sont en effet en mesure de reproduire aussi bien l'augmentation que la
diminution de l’intensite turbulente en ecoulement a bulles. Cela dependrait des poids respectifs des
66
CK4ELTRE 2 MODELEE DE TURBULENCE EA ECODLEMEATE ^ BULLEE
termes source Sk et Se ; leurs effets sur le bilan global de l'energie cinetique turbulence peut etre
Sk
represente par le rapport ou est une norme appropriee.
S
l'ajustement des constantes du modele de turbulence. Nous presentons dans la suite differentes
fermetures adoptees dans des modeles au premier ordre a une echelle de temps.
Lee et al. (1989) expliquent l’augmentation de l’energie turbulente due a la presence des bulles
par la diminution de l’energie potentielle du liquide qui accompagne l'ascension des bulles. Celle ci
induit le deplacement d’un volume du liquide equivalent dans la direction des pressions croissantes,
et se traduit par une diminution de l’energie potentielle du liquide qui contribue a l'accroissement de
l'energie turbulente. Ils en deduisent une formulation pour le terme d’echange interfacial d’energie
turbulente Sk sous la forme simple :
ou c3 est une fonction du nombre de Reynolds du liquide ajustee a partir des experiences de
2000
S =Sk (2-25)
kL
ou c4 est une constante ; dans le cas de l’ascension des bulles dans un liquide au repos, la
67
CH4E/TRE 2 MODELEE DE TURBULEACE EA ECODLEMEATE 4 BULLEE
V t = c„ ^ (2-26)
e
Cette relation augmente la valeur de cm de sorte que les profils de vitesse dans les ecoulements a
0.2
sk 0.037 + (0.21.10-5)ReL +---------- 2 (2-28)
ReL - 65000
1 + exp( )
5000
Ce modele a ete mis en reuvre pour simuler l’ecoulement a bulles dans une conduite cylindrique.
Les resultats numeriques sont compares aux donnees experimentales de Serizawa et al (1974).
Cette comparaison indique une concordance satisfaisante entre resultats numeriques et donnees
experimentales.
Wang et al (1994) modelisent les termes d’echange interfaciaux Sk et Se dans les equations de
A A V = ==
Sk = 2kL (CT -1) +------- (uG -UL).VaG (2-29)
PL P L a L a G SI
ou oI est une constante, d’apres les simulations de l’ecoulement a bulles en conduite, sa valeur est
68
CK4E/TRE 2 MODELEE DE TURBULEACE EA ECODLEMEATE ^ BULLEE
3aGPL C
Ad = UG UL (2-31)
4 dn 1
24
CD =---- (1 + 0.15R0687) pour ReL < 1000 et CD = 0.44 pour ReL > 1000 (2-32)
ReL
3+b
CT = (2-33)
1+b+2
PL
b= (2-34)
p dB El dB Ret
avec :
3/2
R et 3/4 k (2-35)
et 1l = cm
e
Dans l’expression de Sk le second terme a peu d’effet par rapport au premier. En le negligeant, le
proportionnalite egal a la frequence caracteristique des structures turbulentes les plus energetiques
comme pour le modele de Lee et al (1989) :
Se=f Sk (2-36)
k
Le modele de Wang et al (1994) a ete mis en reuvre pour la simulation d’un ecoulement a bulles
(eau/air) dans une conduite cylindrique et dans un elargissement brusque axisymetrique. Il a ete
egalement mis en reuvre dans la simulation de la dispersion d’un scalaire passif dans une zone de
melange (Issa et Oliveira, 1996).
69
CK4E/7EE 2 MODELEE DIE mREUZEACE EA ECO^LEMEATE ^ E^LLEE
developpee dans le mouvement fluctuant par l'ensemble des forces qui s'exercent sur les interfaces.
e
Se = Ce3—Sk ou Ce3 est une constante de l'ordre de l'unite (2-38)
kL
UR = (UG - UL ) - Ud (2-40)
Le coefficient FD qui apparait dans (2-37) correspond au coefficient de trainee donne par :
Fd = 4 P LT'
U, (2-41)
La relation (2-37) fait apparaitre une correlation entre les fluctuations de la vitesse dans la phase
continue et de la vitesse relative exprimee d'apres Simonin (1991) par :
ou qLG represente l'energie turbulente due a l'interaction des phases s'exprimant par :
Ce terme est modelise en se basant sur la theorie de dispersion de Tchen par le modele suivant:
b +h r
qLG = 2kL (2-44)
1 +h r
70
C%4f/TRE 2 MODELEE DE TDREDLEACE EA ECODLEMEATE ^ EDLLEE
b = (1 + CJ ^ + Ca (2-45)
P2
et hr est un parametre lie aux effets d'inertie qui agissent sur la particule :
tt
hr =-f (2-46)
t12
Ce modele a ete mis en reuvre dans plusieurs configurations d'ecoulements tels que l'ecoulement a
bulles dans un elargissement brusque (Bel Fdhila et Simonin, 1992), zones de melange a bulles
(Roig, 1993).
Ce modele a ete implante dans le code indistriel CATHARE (Bestion 1994) et a ete mis en reuvre
dans la simulation des ecoulements diphasiques en conduites verticales correspondant aux
experiences de Liu et Bankoff (1990) et de Grossetete (1995 a, b). Les resultats numeriques
obtenus concernant les grandeurs turbulentes montrent qu'a forts taux de vide avec transition vers
les ecoulements a poches, les resultats des modeles s'ecartent de maniere significative des donnees
experimentales de Grossetete (1995 a, b).
71
CH&P/TRE 2 MODELEE DE TDEEDLEACE EA ECODLEMEATE ^ EDLLEE
Drew et Lahey (1982) ont montre que la distribution du taux de vide en ecoulement a bulles en
conduite est essentiellement controlee par la turbulence. Ils deduisent analytiquement la distribution
radiale du taux de vide dans un ecoulement a bulles en conduite circulaire, et montrent que cette
repartition depend de la structure de la turbulence donc de la repartition de l'energie turbulente
entre les composantes normales du tenseur de Reynolds. Ces resultats ont ete a l'origine de
nombreux travaux experimentaux visant l'analyse de la structure de la turbulence en ecoulement a
bulles. Parallelement des travaux de modelisation se sont interesses au developpement de modeles
au second ordre de la turbulence en ecoulements a bulles.
D'une maniere analogue a la modelisation au premier ordre, les modeles au second ordre se sont
largement inspires de l'experience acquise dans l'etude des ecoulements monophasiques. L'equation
modelisee de transport du tenseur de Reynolds s’ecrit d'une maniere symbolique sous la forme
generale suivante :
Du' u
aT ~ = a L(Pij + Di, +f ij -S 0 ) + Sij (2-48)
Dt
ou Pij est la production turbulente par cisaillement dans la phase liquide, Dij et fij representent les
La dissipation est consideree isotrope conformement a l'hypothese d'isotropie locale et est decrite
dans les modeles au second ordre par l'equation generale suivante :
D e
— (a L e) = D e +a L~ (c e1PL - c e2e) + Se (2-49)
Dt kL
ou De est le terme de diffusion modelise selon le modele de Launder et al (1975) et PL est donne
par la relation la demi trace du terme Pij . Les constantes qui figurent dans cette equation
72
CH4ELTRE 2 MODELEE DE TDRBDLEACE EA ECODLEMEATE 4 EDLLEE
Sij = P% (2-50)
3 C 3
P1 = Ci(%lT P L a L) U G - U L (2-51)
4 d
basant sur la solution theorique de Biesheuvel et Van Wijingaarden (1984) etablie en ecoulement
potentiel a bulles sous la forme :
4 0 0
5
0 3 0
T (2-52)
5
0 0 3
5
Le coefficient empirique CI est fixe a la valeur 0.02 pour tenir compte du fait qu'une faible portion
de l'energie turbulente produite par les bulles est transferee vers les structures tourbillonnaires dans
la phase continue.
Le terme interfacial dans l'equation de transport du taux de dissipation est modelise par :
Ce modele a ete applique dans le cas de l'ecoulement a bulles en conduite et confronte aux
experiences de Serizawa et al. (1974) et Wang et al. (1987).
Lopez de Bertodano et al. (1994) montrent que dans le cas general, la modelisation a une echelle
de temps (au premier ou au second ordre) est insuffisante pour decrire la turbulence en ecoulement
diphasique. En effet, les equations modelisees de transport de l'energie turbulente et du taux de
73
CH4E/ZRE 2 MODELEE DIE mREUZEACE EA ECO^LEMEATE 4 E^LLEE
dissipation (2-21) et (2-22) se reduisent dans le cas d'un ecoulement turbulent homogene
monophasique uniforme aux systeme suivant :
DkL
---- L = -e (2-54)
Dt
De e
D = -c • -t:e (2-55)
D2k De
(2-56)
Dt2 Dt
=c ^(DkL)2
Dt2 e2 kL ( Dt ) (2-57)
dont la resolution nous donne une expression analytique de l'evolution au cours du temps de
l'energie turbulente :
1
kT 1 ( )
[--------- ] ce 2 -1 (2-59)
kT e0 t
1 + (c e2 - 1)
k L0
Avec l'approximation ce2 = 1.92 » 2, l'expression (2-59) peut s'ecrire sous la forme simplifiee :
Cette approximation montre que dans le cas d'une turbulence monophasique homogene l'energie
turbulente s'attenue selon une loi exponentielle avec la constante de temps
74
C%4f/TRE 2 MODELEE DE TDREDLEACE EA ECODLEMEATE ^ EDLLEE
t„ = ^ (2-61)
e0
En ecoulement diphasique, les equations modelisees (2-21) et (2-22) se reduisent dans le cas de
l'ascension des bulles dans un liquide au repos a :
0 = —a L e + Sk (2-62)
e
0 = —a L c e 2 1 e + Se (2-63)
kL
ce qui permet d'etablir la relation suivante entre les deux termes sources Sk et Se :
e
Se = c e 2\ Sk (2-64)
kL
a L^Dt^ = Sk —a L e (2-65)
De e
a L” = Se —a L c e ---e (2-66)
Dt kL
De
a ^=c ^ —a L e] (2-67)
Dt kL
1 De 1 DkL
(2-68)
e Dt e2 kL Dt
e (^)‘e 2
(2-69)
e0 kL0
En divisant les deux membres de l'equation (2-62) par aLkL0, on obtient, compte tenue de (2-69),
75
CH&P/TRE 2 MODELEE DE TDEEDLEACE EA ECODLEMEATE ^ EDLLEE
1 DkT Sk 8 ( kL )ce2
(2-70)
kL0 Dt aLkL0 kL0 kL0
a2 = (2-71)
L
La resolution de l'equation (2-70) permet de deduire une solution analytique de l'energie turbulente
k
—a
k L0 = (1—^)exp( — ^^at)
(2-72)
k 1+a tn
+a
k
On pose :
1 —a
c= (2-73)
1 +a
k
et on suppose — a > 0, l'equation (2-72) permet d'etablir la relation
kL0
1 ( 2at
+ - exp(----- )
1
k c ^0
=a (2-74)
k L0 1 z 2at
1 — exp(------)
c Tn
k
qui permet d'estimer enfin le rapport
k L0
kT 2 2at
a 1 +—exp(------ ) (2-75)
kT c Tn
Cette expression montre que l'energie turbulente en ecoulement diphasique homogene uniforme
s'attenue selon une loi exponentielle de maniere analogue au cas monophasique, mais avec un temps
caracteristique td donne par :
I0
Td = (2-76)
2a
76
CfMP/TRE 2 EZ/ZLES
a= (2-77)
0
Les resultats (2-76) et (2-77), obtenus a partir des equations du modele k-s a une echelle de temps,
sont en contradiction avec les resultats experimentaux obtenus par Lance et al. (1991) concemant
recoupment a bulles homogene uniforme derriere une grille. En effet, l'expression (2-76) montre
rapport depend de la valeur initiale du taux de dissipation ; on devrait s'attendre a ce que cette
valeur depende uniquement de l'energie turbulente induite par la presence des bulles. Ce
raisonnement, qui reste valable egalement pour la modelisation au second ordre, indique que la
modelisation de la turbulence diphasique a une echelle de temps est incomplete ; elle s'avere
insuffisante pour reproduire correctement le temps caracteristique de la decroissance de la
turbulence homogene a bulles a l'aval d'une grille. Au fait la turbulence produite par la grille et la
turbulence induite par les bulles sont deux phenomenes gouvemes par des mecanismes
fondamentalement differents associes a des echelles differentes et il est necessaire de les identifier
pour en preciser les effets sur la structure de recoupment.
77
CK4E/7EE 2 MODEEEE DE mRBDZEACE EA ECODLEMEATE ^ BDLLEE
Le modele Sato et al. (1981) suppose que la vitesse instantanee du liquide uL se decompose en
trois composantes statiatiquement independantes : une valeur moyenne uL, une composante
UL = u L + u L + uL (2-78)
Par consequent, l’energie turbulente dans le liquide kL s’ecrit comme la somme de l’energie
turbulente monophasique kL0 et de l’energie turbulente induite par l’agitation des bulles kLB :
Sato et Sekoguchi (1975) modelisent la viscosite turbulente selon le meme principe en ecrivant :
V t =n tm +n B (2-80)
ou vtm est la viscosite turbulente monophasique et nB est une viscosite turbulente exprimee
algebriquement par analyse dimensionnelle en fonction des grandeurs caracteristiques des bulles :
d
VB = cb 2 aG uT (2-81)
La formulation de cette viscosite turbulente induite par les bulles revient a introduire une nouvelle
echelle de temps caracteristique du mouvement des bulles. En effet, en exprimant l'energie
turbulente induite par l'agitation des bulles kLb selon la solution analytique de Biesheuvel et Van
78
C%4f/TRE 2 MODELEE DE TDREDLEACE EA ECODLEMEATE 4 EDLLEE
2
Wijingaarden, (1984) en ecoulement potentiel homogene kLb = 1 a u l'expression (2-80) de la
Ou t1 et tb representent les deux echelles de temps donnees par les relations suivantes :
k d (2-83)
11 =
e et tb
u
Cependant, on peut montrer aisement que le modele de Sato et al (1981), maintenant tres utilise
dans les codes de calcul, est incomplet ; il n'est pas en mesure de reproduire les diverses situations
de turbulence en ecoulement a bulles. En effet les modeles de double superposition de la turbulence
ne peuvent qu'augmenter la production par les gradients moyens et sont donc incapables de
reproduire l'attenuation de la turbulence observee dans plusieurs experiences d'ecoulements a bulles
(Wang et al 1987 ; Lee et al 1989 ; Liu et Bankoff 1990 ; Serizawa et al 1992). En particulier le
modele de Sato et al (1981) produirait pour des ecoulements a bulles uniformement cisailles, des
frottements turbulents toujours plus importants que ceux obtenus en ecoulement monophasique
avec le meme gradient de vitesse. L'experience de Lance et al (1991) montre que ce n'est pas
toujours le cas, figures (2-4), (2-7), (2-8) et (2-9).
Lopez de Bertodano et al (1994) supposent que la turbulence en ecoulement a bulles peut etre
consideree comme la superposition d'une turbulence monophasique produite par cisaillement dans
la phase liquide et d'une turbulence induite par les bulles de sorte que pour les faibles taux de
presence de la phase dispersee, l'energie cinetique turbulente s'ecrive selon la relation :
S = S SI +e BI (2-85)
ainsi que la viscosite turbulente qu'ils ecrivent en se basant sur le modele de Sato et al. (1981):
79
CH&P/TRE 2 MODELEE DE TURBULENCE EA ECODLEMEATE ^ BULLEE
avec :
V,(SI)
t ^ = cm ; c , = 0.09 (2-87)
SI
L'energie turbulente induite par les bulles kBI est determinee a l'aide d'une equation de transport
= V'[alv,(VkBI)]+-L(kBI -kBI)
(2-89)
Dt -n
ou tB est un temps caracteristique associe au deplacement des bulles et kaBI est l'energie turbulente
induite par les bulles dans un ecoulement potentiel calculee a partir de la solution analytique
d'Arnold (1988) :
2
kBi = ;2cvma G
u (2-90)
1
Avec le coefficient de la masse ajoutee cvm = — on retrouve le modele de l'energie pseudo-
L'energie turbulente induite par cisaillement kSI est modelisee par une equation de transport
Le taux de dissipation eSI associe a la turbulence induite par cisaillement est modelise a l'aide d'une
80
CK4ELTRE 2 MODELEE DE TDRBDLEACE EA ECODLEMEATE ^ EDLLEE
En comparant la somme des equations (2-89) et (2-91) a l'equation generale (2-21) nous identifions
les egalites suivantes :
k
aL8 a L 8SI + (2-93)
iB
k BI
Sk = (2-94)
-B
La relation (2-93) est en accord avec la decomposition (2-85) du taux de dissipation. Elle montre
que la dissipation totale en ecoulement a bulles peut etre interpretee comme une superposition
d'une dissipation induite par cisaillement dans le liquide et d'une dissipation induite par les bulles qui
se produit au niveau des sillages et qui est controlee par l'echelle de temps caracteristique -B . Ce
k
8 BI =' (2-95)
iB
3 cT
Sk =a G4-^T u, (2-96)
Les equations (2-90), (2-94) et (2-96) permettent d'etablir l'expression du temps -B qui s'ecrit :
2 0vm. d
- B = (l^) (2-97)
3 cT u,
Applique a l'ecoulement de turbulence homogene uniforme derriere une grille, le modele de Lopez
de Bertodano et al (1994) s'ecrit :
Dksi
= -e (2-98)
Dt
Desi = ^81 e
(2-99)
Dt 82 kL SI
81
CK4ELTRE 2 MODELEE DE TURBULENCE EA ECODLEMEATE ^ BDLLEE
du du
-u v = v — avec — = S » constante (2-102)
' dy dy
Revenons sur le modele de Lopez de Bertodano et al. (1994) ; les equations du modele s'ecrivent
en ecoulement de cisaillement uniforme :
2
kBI = kBI = 2 cv u (2-105)
82
CK4E/TRE 2 MODELEE DE TDRBDLEACE EA ECODLEMEATE ^ BDLLEE
k2
v (SI)
t =c m (2-106)
SI
v = c mb 2aG u (2-107)
que le modele de Lopez de Bertodano et al. (1994) calcule la meme valeur de la viscosite
turbulente induite par cisaillement v(tSI) que celle que produirait un modele monophasique dans les
== ^=(SI) ^=(BI)
u'v' = u'v' + u'v' (2-108)
Le modele de Lopez de Bertadano et al. (1994), bati sur le principe de superposition des viscosites
turbulentes de Sato et al (1981), ne peut pas reproduire la reduction du frottement turbulent
observee en ecoulement homogene avec cisaillement constant.
Ce modele a ete implante dans le code industriel de CFD "PHOENICS" (Spalding et al 1988) avec
les equations de bilans moyens de masse et de quantite de mouvement dans le liquide et dans le gaz
modelisees selon le modele de Ishii (1975) et mis en reuvre dans la simulation de deux ecoulements
ascendants a bulles en conduites cylindriques verticales et confronte aux experiences de Serizawa et
al. (1986) et Wang et al. (1987). Les resultats numeriques obtenus indiquent que le modele
reproduit d'une maniere satisfaisante les resultats experimentaux de Serizawa et al. (1986) avec
l'ajustement du coefficient de la masse ajoutee pour lequel on retient la valeur cvm = 2. Le modele
83
CfMP/TRE 2 .' M9DE1ESDE rZ/REZ/Z^MTEEA/ECOZ/LaZEVrS/t EZ/ZZES
0.1S0-
~0.100"
Figure (2-11) : Comparaison des composantes normales du tenseur de Reynolds avec les donnees
experimentales de Serizawa et al. (Lopez de Bertodano et al. (1994))
0.006
- □ ,
Figure (2-12) : Comparaison du frottement turbulent avec les donnees experimentales de Serizawa
et al. (Lopez de Bertodano et al. (1994))
84
CfMP/TRE 2 .' EZ/ZZES
j = 1.36 m/s
0.30 -
0.20 -
0.10 -
Figure (2-13) : Production par cisaillement de la turbulence dans le cas monophasique et diphasique
(Lopez de Bertodano et al. (1994))
1.20 -
0.80 '
Figure (2-14) : Profil de la vitesse moyenne du liquide dans le cas mophasique et diphasique (Lopez
de Bertodano et al. (1994))
85
CH4E/ZRE 2 MODEEEE DE mRBLZEACE EA ECODLEMEATE 4 EDLLEE
Cependant, ces resultats plutot positifs doivent etre interpretes vis-a-vis de l'effet des gradients
moyens sur la structure de la turbulence du liquide : l'ajustement de l'echange interfacial se traduit
par un ajustement des gradients de vitesse de sorte que la production de la turbulence par
cisaillement s'attenue. Celle ci s'ecrit :
ou r est la distance radiale par rapport a l'axe de la conduite (figures (2-13) et (2-14)).
et diphasiques avec deux taux de vide differents) montre qu'avec la meme valeur cvm = 2, les
niveaux de turbulence produits par le modele et confrontes aux resultats experimentaux n'indiquent
pas une concordance aussi bonne que celle obtenue dans le cas de l'experience de Serizawa et al.
(1986), figures (2-15) et (2-16). La figure (2-17) representant le frottement turbulent montre
d'importants ecarts avec l'experience qui expliquent la mauvaise prediction du profil de vitesse
moyenne, figure (2-18). D'ailleurs, les auteurs remarquent que les resultats sont ameliores en
utilisant la valeur cvm = 1.2 au lieu de 2, ce qui revient a un nouvel ajustement de l'echange
interfacial.
86
C/MP/7%E 2 .' rZ/REZ/Z^MTEE/yECOZ/LaZEVr^yt EZ/ZZES
model
0.16 -
j = 1.08 m/s
ic = 0.10 m/s
0.12 -
0.08 -
0.04 ■
Figure (2-15) : Comparaison des composantes normales du tenseur de Reynolds avec les donnees
experimentales de Wang et al. (Lopez de Bertodano et al. (1994))
model
0.16 ■
j = 1.08 m/s
j = 0.40 m/s
0.12 '
- 0.08 -
0.04 ■D
Figure (2-16) : Comparaison des composantes normales du tenseur de Reynolds avec les donnees
experimentales de Wang et al. (Lopez de Bertodano et al. (1994))
87
CfMP/TRE 2 .' EZ/ZZES
model
0.010 ■
0.005
0.000 -----$-
Figure (2-17) : Comparaison du frottement turbulent avec les donnees experimentales de Wang et
al. (Lopez de Bertodano et al. (1994))
1.50 -
1.00 -
0.50 ■
model
j = 1.08 m/s
Figure (2-18) : Profil de la vitesse moyenne du liquide dans le cas diphasique a deux taux de vide
differents (Lopez de Bertodano et al. (1994)
88
CH4E/ZRE 2 MODEEEE DE mRBDZEACE EA ECODLEMEATE 4 BDLLEE
Troshko et Hassan (2001) partent sur un principe de fermeture a deux echelles de temps pour
kL est modelisee selon (2-3). Le terme de production interfacial Sk est exprime par la puissance
developpee par la force interfaciale dans le mouvement relatif moyen essentiellement controle par la
ou PL est la production par le gradient de la vitesse moyenne du liquide et Sk est modelise par :
3
Sk a, u, (2-112)
4 d
Les deux echelles de temps sont introduites dans la modelisation de l'equation de transport du taux
de dissipation. Celle ci est modelisee selon l'equation (2-22) dans laquelle le mecanisme de
destruction du taux de dissipation est associe a des frequences propres a la turbulence et a la
pseudo-turbulence ; on considere que le phenomene de destruction est realise a deux echelles de
temps distinctes : une echelle controlee par les structures energetiques de la turbulence dans la
kT
phase liquide analogue a celle de l'ecoulement monophasique t0 et une echelle decrivant les
e
mecanismes de destruction au niveau des sillages des bulles pour laquelle on adopte l'expression (2
97) proposee par Lopez de Bertodano et al. (1994). L'equation modelisee de transport du taux de
dissipation est ainsi modelisee sous la forme :
D
(a L e) = V. a L(— + n)Ve + a L W0(c e1PL c e2 e) + W bc3Sk (2-113)
Dt ak
89
C%4f/TRE 2 MODELEE DE TDREDLEACE EA ECODLEMEATE ^ EDLLEE
e
wo = T 0 (2-114)
k
u,
-1 3 cT
wb = T b = (2-115)
2 c d
k2
nt = cm (2-116)
g
Les interactions interfaciales agissent par consequent d'une maniere implicite sur la viscosite
turbulente par la modification des valeurs de l'energie turbulente kL et de son taux de dissipation e
en relation avec l'effet des termes interfaciaux qui apparaissent dans les equations de transport.
Le modele de Troshko et Hassan (2001) ne se prete pas a une analyse simple dans le cas de la
turbulence homogene a cisaillement constant : l'expression de la viscosite turbulente, exprimee en
fonction de l'energie totale et de son taux de dissipation ne permet pas une separation des echelles.
Cependant, l'analyse de l'ecoulement de turbulence homogene de grille, permet une evaluation des
fermetures proposees. Dans ce cas, le modele de Troshko et Hassan (2001) se reduit a :
Dk L
e+Sk (2-117)
Dt
avec la relation :
>> rn (2-119)
Quand les bulles atteignent leur vitesse limite de glissement le terme Sk est constant, le systeme (2-
Dk
= —e = —e + S,, (2-120)
Dt
De*
= —c e 2 We + (c3Wb — c e 2 w)Sk (2-121)
Dt
90
CH&P/TRE 2 MODELEE DE TURBULENCE EA ECODLEMEATE ^ BULLEE
Rappelons que les resultats experimentaux de Lance et Bataille (1990) (Figure (2-3)) indiquent que
la decroissance de la turbulence en ecoulement diphasique avec le meme temps caracteristique que
dans le cas de l'ecoulement monophasique. Pour que ce modele reproduise la decroissance de la
turbulence en ecoulement diphasique a la frequence w de l'ecoulement monophasique equivalent,
c3wb — c e 2w = 0 (2-122)
qui est incompatible avec la relation (2-105). Remarquons egalement que la frequence w donnee
par la relation :
e e SI + e BI
w =— (2-123)
kL kSI + kBI
Le modele de Troshko et Hassan (2001) a ete implante dans le code industriel de CFD "CFX 4.2"
(AEA Technology, 1997) et applique a la simulation des experiences de l'ecoulement a bulles en
conduite verticale de Wang et al. (1987), Serizawa et al. (1975) et de Liu (1998). Les resultats
numeriques relatifs la simulation des experiences de Wang et al. (1987) sont compares a ceux
obtenus par le modele de Lopez de Bertodano et al. (1994), figures (2-19) - (2-24) ; la
comparaison montre que les deux modeles produisent des resultats comparables. Les donnees de
Wang et al. (1987) sont relatives a trois experiences d'ecoulements a bulles en conduite verticale
Tableau (2-3) : Parametres caracteristiques des experiences de Wang et al. (1987) considerees
W0 0.43 0 - Ascendant
91
CfMP/TRE 2 MODELED DE rZ/REZ/LEMTEE/ZECOZ/LEAZEAT^yt EZ/ZLES
Le modele est d'abord applique a l'ecoulement monophasique pour lequel il produit des resultats
parfaitement concordants avec les donnees experimentales. Le coefficient c3 a etc ajuste a la valeur
0.45 dans la simulation de l'experience W1 pour bien reproduire le pic du frottement turbulent au
voisinage de la paroi (figure (2-19)). Dans le meme temps, le modele reproduit convenablement
l'intensite turbulente longitudinale. Cependant avec cet ajustement, le modele a du mal a reproduire
la structure de la turbulence dans les deux autres ecoulements diphasiques : les intensites
turbulentes sont surestimees et le frottement turbulent est mal reproduit dans l'experience W2.
5.00E-03
. 4.50E-03 ♦ WO
4.00E-03
■ W1
3.50E-03
Lopez de Bertodtino et til., W1
—^ Troshko et Htisstin, W1
Troshko et Hassan, WO
-S 1.50E-03
tr 1.00E-03
5.00E-04
r/R
Figure (2-19) : Profils transversaux du frottement turbulent obtenus par le modele de Troshko et
Hassan en W0 et W1 et par le modele de Lopez de Bertodano et al. en W1 (D'apres Troshko et
Hassan, 2001)
92
CfMP/TRE 2 EZ/ZZES
O WO, ll'
—D— wo, v
—♦— Wl, ll
0.00 -
Figure (2-20) : Profils transversaux des intensites turbulentes obtenus par le modele de Troshko et
Hassan en WO et Wl et par le modele de Lopez de Bertodano en Wl (D'apres Troshko et Hassan,
2001)
■ W2
«— Troshko et Hassan, W2
Troshko et Hassan, WO
Figure (2-21) : Profils transversaux du frottement turbulent obtenus par le modele de Troshko et
Hassan en W0 et W2 (D'apres Troshko et Hassan, 2001)
93
CfMP/TRE 2 MODELED DE rZ/REZ/LEMTEE/ZECOZ/LEAZEjVDS'yt BUZZES
■£ 0.14 -
0.12 -
0.08
0.04 -
Figure (2-22) : Profil transversal de l'intensite turbulente definie par u' obtenu par le
♦ W3
— Troshko et Hassan, W3
-5 0.005
a 0.002
Figure (2-23) : Profils transversaux du frottement turbulent obtenus par les modeles de Troshko et
Hassan et Lopez de Bertodano et al en W3 (D'apres Troshko et Hassan, 2001)
94
CfMP/TRE 2 MODELED DE rZ/REZ/LEMTEE/ZECOZ/LElZEjVDS'yt EZ/ZLES
0 15 -
0.to *■
Figure (2-24) : Profils transversaux de l'intensite turbulente definie par u' obtenu par le
modele de Troshko et Hassan et des intensites turbulentes obtenues par le modele de Lopez de
Bertodano et al en W3 (D'apres Troshko et Hassan, 2001)
La figure (2-25) presente la viscosite turbulente calculee par le modele dans les differentes
simulations.
Figure (2-25) : Profils transversaux de la viscosite turbulente reproduits par le modele de Troshko
et Hassan dans les experiences WO, Wl, W2 et W3 (D'apres Troshko et Hassan, 2001)
95
CfMP/TRE 2 .' EZ/ZZES
Lance et al. (1991) ont propose un modele au second ordre pour les ecoulements turbulents
homogenes a bulles. Rappelons qu'en 1'absence de diffusion, l'equation de transport du tenseur de
Reynolds s'ecrit d'une maniere schematique sous la forme :
(2-124)
avec :
PL est le terme de production par cisaillement dans la phase liquide de la turbulence qui s'ecrit
96
CH4E/ZRE 2 MODEEEE DE mRBDZEACE EA ECODLEMEATE 4 BDLLEE
Dans le modele de Lance et al. (1991), le tenseur de Reynolds represente toutes les fluctuations de
vitesse du liquide ; celles produites par cisaillement et celles induites par les interactions
interfaciales. Cependant la dissipation visqueuse est decomposee en deux contributions
dans le liquide resultant de la cascade d'energie (supposee isotrope) et une dissipation eS dans les
Compte tenu des resultats experimentaux obtenus en turbulence homogene (paragraphe 2.2), Lance
et al (1991) supposent qu'il y a equilibre production - dissipation dans les sillages. Dans ces
conditions, la seule dissipation qui subsiste dans les equations de transport des tensions turbulence
e L = 3 e 01 (2-127)
ou e0 qui en ecoulement de turbulence homogene avec cisaillement constant, est modelisee comme
L'apport original dans le modele de Lance et al. (1991) est relatif a la modelisation du terme de
redistribution F de l'equation (2-124). Ce terme est decompose comme en monophasique en partie
lineaire (appelee egalement terme rapide) associee au couplage lineaire entre le mouvement moyen
et fluctuant et partie non lineaire associee au couplage non lineaire entre les fluctuations de vitesse
tel que :
Les auteurs retiennent la modelisation monophasique de Launder et al (1975) pour la partie lineaire
et proposent pour la partie non lineaire, d'introduire un effet supplementaire pour prendre en
97
CfMP/TRE 2 EZ/ZZES
compte la distorsion aleatoire des structures turbulentes provoquee par le mouvement turbulent des
bulles. Lance et al (1991) proposent done de modeliser la partie non-lineaire sous la forme:
uT
o(NL) “(Cl7+C3
a )(U>L (2-130)
k
Les resultats numeriques de 1'evolution longitudinale des composantes du tenseur de Reynolds sont
compares aux donnees experimentales de Lance et al. (1991), figures (2-4), (2-7), (2-8), (2-9) et
(2-10). La confrontation des resultats numeriques aux donnees experimentales montre une
concordance satisfaisante. En particular, le modele reussit a produire l'accentuation de l'isotropie et
la reduction du frottement turbulent induites par la presence des bulles.
Cependant, ce modele de retour a l'isotropie, qui comprend l'ensemble des fluctuations conduit, en
l'absence de contraintes exterieures imposees par le champ moyen, a une solution isotrope de la
turbulence totale. Cela est incompatible avec les resultats experimentaux obtenus en ecoulement ou
la turbulence induite par les bulles est importante. En effet, les resultats experimentaux de Mareuge
et Lance (1995) de l'ascension des bulles dans un liquide au repos indiquent que la turbulence
produite par les bulles est a caractere fortement anisotrope, figure (2-26). On remarque egalement
que la turbulence en ecoulement a bulles a l'aval de la grille est anisotrope, Lance (1986).
(■IQ2) o u" in s a v* in s
Figure (2-26) : Turbulence induite par des bulles en ascension libre dans un liquide au repos
(D'apres Mareuge et Lance, 1995)
98
CH4E/ZRE 2 MODEEEE DE mRBDZEACE EA ECODLEMEATE 4 EDLLEE
Chahed (1999) s'appuie sur les travaux de Lance et al. (1991) pour proposer un modele au
second ordre ou le tenseur de Reynolds est decompose en deux contributions statistiquement
independantes : une composante "pseudo-turbulente" associee aux fluctuations irrotationnelles non
dissipatives d'origine purement interfaciale et une composante associee aux fluctuations
"turbulentes" rotationnelles dissipatives qui comprend egalement les fluctuations dans les sillages
des bulles:
— (S) = = (T)
« =« +« (2-131)
Chacune des composantes est decrite par une equation de transport specifique. L'equation de
transport de la composante pseudo-turbulente est modelisee en s'appuyant sur la solution theorique
de Biesheuvel et Van Wijingaarden (1984) etablie en ecoulement potentiel homogene a bulles. Dans
cet ecoulement le tenseur des fluctuations turbulentes induites par le mouvement des bulles s'ecrit :
r(S)
3 I 2j 1
u'Lu L H =—a u 1 +------- a UjU (2-132)
20 " R 20
RR
(S) =(S)
D -(S)
Dt
uLuL = Diff(uLuL )+Dt(uLuL h ) (2-133)
Le dernier terme de l'equation (2-133) exprime l'ecart par rapport a la solution potentielle (2-132)
en ecoulement homogene. Il represente une diffusion due precisement au caractere inhomogene de
l'ecoulement. Le dernier terme de (2-133) est interprete comme la contribution, en situation
inhomogene, de la production interfaciale qui represente dans ce cas la puissance developpee par la
force de la masse ajoutee, et de la redistribution interfaciale par la correlation pression-deformation.
99
C%4f/TRE 2 MODELEE DE TDREDLEACE EA ECODLEMEATE ^ EDLLEE
D (T) (T) 2
— (uLuL ) = Diff(uLuL ) + Pl +F - 3 e 01 (2-134)
ou e0 represente la dissipation visqueuse isotrope qui resulte de la cascade energetique ; elle est
De e
—0 = Diff(e 0) + -- [c e1 Pl :I - c e2 e 0 ] (2-135)
Dt k
(T) ^k0l)
F(NL) c1(t-1 +at b1)(uLuL
(2-136)
kn 1 (T)
d
1 b = cR (2-138)
UR
La partie lineaire du terme de redistribution est modelisee d'une maniere analogue a la modelisation
monophasique de Launder et al. (1975).
Les termes de diffusion qui apparaissent dans les equations de transport (2-133), (2-134) et (2-135)
sont modelises en generalisant le modele de Launder et al. (1975). Ainsi, la diffusivite est exprimee
comme la somme des deux contributions une contribution turbulente associee a l'echelle de temps
100
CH&P/TRE 2 MODELEE DE TURBULENCE EA ECODLEMEATE ^ BULLEE
C = (T) (S)
Diff(j) (1 -a)(Tt uLuL +1 b uLuL )-vj (2-139)
(1 -a)
Ce modele a ete mis en reuvre pour la simulation de plusieurs ecoulements turbulents a bulles:
homogene uniforme derriere une grille, homogene a cisaillement constant et cisaille libre de sillage
correspondant respectivement aux experiences de Lance et Bataille (1990), Lance et al. (1991) et
Roig et al. (1998). Les equations du modele sont discretisees par la methode des differences finies
et integrees en mode parabolique. Les resultats numeriques obtenus montrent une concordance
satisfaisante avec les donnees experimentales dans les differentes configurations d'ecoulements
comme le prouvent les figures (2-27)-(2-30).
seul coefficient (le coefficient cR) qui apparait dans cette modelisation en plus des coefficients
cR = 0.67.
101
CK4E/7EE 2 MODELEE DIE mREUZEACE EA ECO^LEMEATE X E^LLEE
0,12
— a =0. ♦ data
0,1 - .... a =0.005
— a =0.01
0,08 -
__ a =0.02
0,06
-----•_
0,04
0,02 -
20 40 60
X/M
Figure (2-27) : Simulation par le modele de Chahed (1999) de l'intensite turbulente en turbulence
homogene de grille monophasique et a bulles. Comparaison avec les resultats experimentaux de
Lance et Bataille (1990)
Figure (2-28) : Simulation par le modele de Chahed (1999) des composantes du tenseur de
Reynolds en turbulence homogene monophasique soumise a un cisaillement constant. Comparaison
avec les resultats experimentaux de Lance et al (1990)
102
CK4E/7EE 2 MODELEE DIE mREUZEACE EA ECO^LEMEATE X E^LLEE
Figure (2-29) : Simulation par le modele de Chahed (1999) des composantes du tenseur de
Reynolds en turbulence homogene diphasique (a = 1 %) soumise a un cisaillement constant.
X=0.2 m
0,08 -
0,02 AAA’
Figure (2-30) : Simulation par le modele de Chahed (1999) de l'intensite turbulente longitudinale en
ecoulement de sillage monophasique et a bulles. Comparaison avec les resultats experimentaux de
Roig et al (1998)
103
CK4P/TRE 2 MODELEE DE TDREDLEACE EA ECODLEMEATE 4 EDLLEE
L'application du modele dans le cas de l'ecoulement du sillage a bulles, montre que le modele
reproduit bien les grandeurs turbulentes dans des situations ou la pseudo-turbulence est importante:
il predit bien l'augmentation importante de l'energie cinetique de l'ecoulement a bulles ainsi que
l'elargissement des profils de vitesse moyenne.
Examinons maintenant, comme pour les autres modeles, le comportement du modele dans les
situations de l'ecoulement homogene. Sans contribution turbulente (cas de l'ascension des bulles
dans un liquide au repos), on peut aisement verifier que le modele degenere vers une solution
asymptotique compatible avec la solution de Biesheuvel et Van Wijingaarden (1984) :
(T)
« =0 (2-140)
En ecoulement turbulent homogene derriere une grille, le modele de Chahed et al (2002) se reduit a
une simple superposition de la turbulence et de la pseudo-turbulence qui explique que la
decroissance de l'intensite turbulente a l'aval de la grille est comparable a celle de l'ecoulement
monophasique.
En turbulence homogene avec cisaillement uniforme, les equations de transport des tensions de
Reynolds donnent :
uv = uV avec u v = 0 (2-142)
—u 'v ' = -v '2S - C1^^ + at -1)u 'v ' + (g1v'2-g 2k + g 3u '2)S (2-145)
Dt kn
D ^=(T) == G =—(T) 2 2 == 2
—u '2 =-2u'v 'S - C1^^ + at -1)(u '2 - ^k0) + ^(2g1 -g 3)u 'v 'S - 3 G 0 (2-146)
Dt 1 k b
104
CH&P/TRE 2 MODELEE DE TDEEDLEACE EA ECODLEMEATE ^ EDLLEE
^D=(T) (T) 2 2
~v '2 = -C1 (” + at b1)(v '2 k0) -T" (g1 - 2g 3)u'v 'S -~s c (2-147)
Dt k 3 3
D ^2(S) 4 D =2
u au (2-148)
Dt 20 Dt
D ^2(S) 3 D =2
v au (2-149)
Dt 20 Dt
k = k0+kS (2-150)
avec :
tb etant le temps caracteristique du deplacement des bulles donne par (2-127) et le taux de
2
De 0 8
C1^-^u 'v 'S- C2G 0 (2-152)
Dt k0 k
diphasique. L'echelle de temps tb caracteristique de l'agitation des bulles a ete ajustee par rapport a
ce cas d'ecoulement et la valeur retenue pour le coefficient CR est egale a 0.67. Les autres
constantes dans cette modelisation sont egales aux valeurs conventionnelles adoptees dans la
modelisation au second ordre monophasique (Launder et al. 1975). Elles sont recapitulees dans le
tableau (2-4) suivant :
105
CK4E/TRE 2 MODELEE DE TURBULENCE EA ECOULEMEATE ^ BULLEE
Lorsque l'ecoulement atteint un etat d'equilibre asymptotique caracterise par des valeurs
sensiblement constantes des composantes du tenseur de Reynolds, l'equation (2-136) se reduit a un
equilibre production-redistribution qui permet d'exprimer algebriquement le frottement turbulent
sous la forme :
2.4 Conclusion
Cette analyse a montre que les modeles de turbulence au premier ordre ne sont pas capables de
reproduire la structure de la turbulence en ecoulement a bulles meme si on procede a la
decomposition des echelles mises en jeux. A ce niveau de fermeture, les modeles sont, dans leurs
formulations classiques, inadaptes pour representer les mecanismes de redistribution en particulier
lorsque ces mecanismes sont alteres par la presence des interfaces. L'analyse indique egalement que
la fermeture au second ordre de la turbulence apparait comme un niveau de fermeture ou il devient
106
CK4E/TRE 2 MODELEE DE TURBULEACE EA ECOULEMEATE ^ BULLEE
possible de rendre compte des effets interfaciaux a condition de les preciser au travers d'echelles
caracteristiques de la turbulence et de la turbulence induite par les bulles. Cela conduit a modeliser
des equations de transport specifiques pour la partie turbulence et pour la partie pseudo-turbulente
(induite par les bulles) du tenseur de Reynolds. Rappelons finalement que la reduction des
fermetures de la turbulence au second ordre developpees par Chahed (1999) permet d'exprimer le
frottement turbulent en ecoulement a bulles selon le concept de viscosite turbulente.
107
CK4E/TRE 3 MODELEE DE TURBULEACE ^ TRO/E EgU^T/OAE ROUR LEE ECOULEMEATE ^ BULLEE
3.1 Introduction
L'analyse des modeles de turbulence developpee dans le chapitre 2 indique qu'il faut remonter au
second ordre pour arriver a un niveau de fermeture capable de prendre en compte les effets
interfaciaux sur la structure de la turbulence en ecoulement a bulles. Elle indique egalement qu'il est
necessaire de proceder a une decomposition de la turbulence pour lui associer des echelles
differentes sensees representer la turbulence propre du liquide et la turbulence induite par les bulles.
Les relations de fermeture sont alors basees sur la decomposition du tenseur de Reynolds de la
phase continue en une partie « pseudo-turbulente » produite par le mouvement relatif des inclusions
et une partie « turbulente » produite par cisaillement dans la phase continue et dans les sillages ;
chaque partie est predeterminee separement par une equation de transport.
Nous nous proposons de developper un modele au premier ordre pour les ecoulements
turbulents a bulles deduit par reduction des fermetures developpees au second ordre. En effet, selon
108
CK4R/TRE 3 MODELEE DE TURBULEACE ^ TRO/E EgU^T/OAE ROUR LEE ECOULEMEATE ^ BULLEE
la formulation des fermetures des equations de transport des tensions de Reynolds, le bilan simplifie
production-redistribution dans l'equation du frottement turbulent en ecoulement de cisaillement pur
permet d'etablir une expression de la viscosite turbulente qui met en jeu des echelles
caracteristiques de la turbulence et de la pseudo-turbulence, Chahed (1999). Cette expression de la
viscosite turbulente est a la base de la formulation de modeles au premier ordre adaptes aux
ecoulements a bulles a faibles taux de vide que nous nous proposons de developper dans ce
chapitre.
Apres avoir analyse les implications des fermetures au second ordre sur la predetermination de la
turbulence en turbulence homogene avec et sans cisaillement, nous proposons des reductions de ces
modelisations. Ces reductions conduisent a une formulation analytique des composantes du tenseur
de Reynolds decompose en parties turbulente et pseudo-turbulente. Cette analyse est interessante
dans la mesure ou elle permet, dans le cadre de l'hypothese d'homogeneite, d'exprimer les parties
turbulentes et pseudo-turbulentes du tenseur de Reynolds ; on a en effet besoin de specifier ces
contributions pour definir les conditions aux limites des modeles ou on procede a la modelisation
d'equations de transport pour chacune des contributions. Les reductions des fermetures au second
ordre permettent egalement de developper des modeles au premier ordre (modele k-e et modele k-
ro) que nous appliquerons ensuite a deux ecoulements a bulles verticaux: il s'agit de l’ecoulement
homogene en aval d'une grille (Lance et Bataille 1991), et de l’ecoulement a cisaillement constant
(Lance et al 1991).
On utilise dans la suite les fermetures au second ordre de la turbulence pour interpreter les
resultats experimentaux de l'experience de l'ecoulement de turbulence homogene a bulles avec
109
CH4R/TRE 3 MODELEE DE TURBULEACE 4 TRO/E EQU4T/OAE ROUR LEE ECOULEMEATE 4 BULLEE
cisaillement constant. Cette analyse permet de mettre en evidence les echelles caracteristiques
relatives aux effets interfaciaux sur la structure de la turbulence. Le modele de Chahed (1999) est
formule ci-dessous dans le cas d'un ecoulement bidimensionnel turbulent homogene a bulles soumis
a un cisaillement constant dans la situation asymptotique de l'equilibre energetique. Les equations
(2-145)-(2-149) de ce modele de turbulence explicitees au chapitre 2 se reduisent sous la forme :
0 = -v '2 S - + at-1)u 'v ' + (g1 v '2 - g2k + g3u '2)S (3-1)
kn
2. . 2
0 — -2u'v'S — Cj (----+ atb )(u ' — k0) +— (2g^ — g3)u'v r S —e 0 (3-2)
k0 33 3
^ . 2
0 — -C1("----+ at b )(v ' -~ k0) — T" (g1 — 2g 3)u'v 'S - “e 0 (3-3)
k0 33 3
^=(S) 4 ^=2
u — — a uR (3-4)
20 R
=(S) 3 =2
v f2 J au
__ (3-5)
20
d
ou S — est le taux de cisaillement et tb est le temps caracteristique du deplacement
5y b u,
des bulles. Le bilan energetique se reduit pour la partie turbulente a un equilibre production-
dissipation (equation (3-6)) et la partie pseudo-turbulente de l'energie demeure constante et
uniforme (equation (3-7)) :
1 --- 2
kS — 4 auR (3-7)
110
CK4R/TRE 3 MODELEE DE TURBULEACE ^ TRO/E EgU^T/OAE ROUR LEE ECOULEMEATE ^ BULLEE
L'equation (3-8) est analogue a celle deduite dans le cas monophasique a ceci pres qu'elle fait
apparaitre au denominateur deux echelles de temps caracteristiques de la turbulence et de la
pseudo-turbulence. La decomposition des echelles au second ordre va s'averer essentielle dans
l'interpretation des resultats avec des modeles d'ordre inferieur.
Compte tenu de l'equilibre production-dissipation (equation (3-6)), les equations (3-2) et (3-3)
permettent de deduire les expressions suivantes des parties turbulentes des composantes normales
du tenseur de Reynolds :
2
(2 - 2g1 +g3)e0
u ,2 — - kn +, 3 (3-9)
3 qA + at-1)
k
2
= (T) ^=(T) 2 (g1 - 2g3 - 1)e 0
v '2 — w '2 ——2k. + ^ (3-10)
C1( + at -1)
k„
Il est egalement possible de deduire les expressions des contributions pseudo-turbulentes des
composantes normales du tenseur de Reynolds, a partir des equations (3-4), (3-5) et (3-7) :
^=(S) 4 ^2 4
u'2 — 20a uR — 5ks (3-11)
2
(2 - 2g1 +g3)e0
u2 — 2k„ + Ik +1
3 0 + 5= + (3-13)
€1^-^ + at -1)
k0
2
3 (g1 - 2g3 - 1)e 0
— 2kn + 3k +1
v —w — — k0 +— ks + - (3-14)
3 0 5 s C1(^ + at -1)
k0
111
CK4R/TRE 3 MODELEE DE TURBULEACE ^ TRO/E EgU^T/OAE ROUR LEE ECOULEMEATE ^ BULLEE
La relation (3-8) du frottement turbulent en ecoulement diphasique traduit un equilibre entre les
mecanismes de production par cisaillement du frottement turbulent et celui de la correlation
pression deformation. Elle met en jeu deux echelles de temps qui resultent de la decomposition des
mecanismes d'etirement tourbillonnaire modelise sous forme d'une superposition des effets
turbulents et d'effets induits par les bulles dans leur mouvement aleatoire. Notons egalement que les
composantes normales du tenseur qui apparaissent au numerateur comprennent les contributions
turbulentes et pseudo-turbulentes. En decomposant l'energie des fluctuations en contributions
turbulente et pseudo-turbulente, le numerateur de l'equation (3-8) peut se mettre sous la forme
suivante :
^(T) ^(T)
v2 u '2
C1C m — (1 -g1^ g 3^ + g2 (3-16)
kr kn
:5(1 -Y1) 4
C1C mb 5g3 +g 2 (3-17)
Compte tenue de la decomposition (3-15), la relation (3-8) permet d'etablir la formulation suivante
de la viscosite turbulente diphasique :
(1 + Cmbks)
n t — Cm f ^
(3-18)
(1 + ^^)
b
structure de la turbulence est modifiee par la presence des bulles, on peut s'attendre a ce que le
C
coefficient C m, et par consequent le rapport , soit dependant du taux de vide ainsi que des
112
CK4R/TRE 3 MODELEE DE TURBULEACE ^ TRO/E EgU^T/OAE ROUR LEE ECOULEMEATE ^ BULLEE
caracteristiques de la phase dispersee (diametre des bulles, vitesse relative). Les resultats
C
numeriques obtenus par le modele au second ordre indiquent que le rapport reste voisin de
l'unite pour la gamme des faibles taux de vide simules. La valeur du coefficient Cm est choisie egale
(1 + jk1)
kn
I
(1+ k1)
k0
vt — C„ — v t0 avec C m — 0.09 (3-19)
e0 (1 + a^^) T
(1 + a-L)
Tb Tb
Cette formulation met en evidence la contribution des interactions interfaciales dans la modification
de la structure de la turbulence par comparaison au cas monophasique. L’agitation des bulles,
caracterisee par l’energie pseudo-turbulente kS, induit une augmentation des fluctuations
turbulentes dans la phase continue. Par contre, retirement tourbillonnaire provoque par le
mouvement relatif des bulles engendre une attenuation du frottement turbulent et par suite une
attenuation de la production par cisaillement de la turbulence dans la phase continue. En
ecoulement a bulles a fort montant de l'energie pseudo-turbulente, la modelisation (3-19) reproduit
une augmentation de la viscosite turbulente diphasique. En revanche, lorsque le montant de
l'energie pseudo-turbulente est faible par rapport au montant de l'energie turbulente, la formulation
(3-19) reproduit dans certains cas une attenuation de la viscosite turbulente par comparaison au cas
monophasique equivalent.
Pour mieux preciser les mecanismes qui entrent en jeu dans la modification de la structure de la
turbulence en ecoulement a bulles, on s'appuie sur la formulation (3-19) de la viscosite turbulente
diphasique pour interpreter l'alteration de la structure de la turbulence, caracterisee par le tenseur
113
CK4R/TRE 3 MODELEE DE TURBULEACE ^ TRO/E EgU^T/OAE ROUR LEE ECOULEMEATE ^ BULLEE
u'u 2
b ij 5 ij (3-20)
k 3
2,
b(jS) +-5;
(1 + ks 3
)
k0 b(T) 2
2 Ij k^ - 2 «.
b,j — (b(jT) + 3 5 ,j)- (3-21)
(1+ A1)
k„
avec :
(T) ^=(S)
- 35 ^ b(1 — - 25
b (t) — u‘uj et (3-22)
bij — k
0
tenseur d'anisotropie de l'ecoulement monophasique bi(jm) . L'analyse est en particulier focalisee sur
les composantes du tenseur b11, b22 et b12 qui peuvent etre deduites a partir des resultats
experimentaux de Lance et al. (1991). Les relations algebriques (3-8), (3-9) et (3-10), permettent
d'ecrire ces composantes sous la forme :
1 +^
T k
(m) (3-23)
h(m) Tt Tt
12
1 +a
Ib
b 11(T) £ 1
(m) ,(m) (3-24)
b 11 22 1 +a—
Tb
114
CK4R/TRE 3 MODELEE DE TURBULEACE ^ TRO/E EgU^T/OAE ROUR LEE ECOULEMEATE ^ BULLEE
ou ttm) est l'echelle de temps de l'ecoulement monophasique donnee par la relation (3-30). D'un
b(1)—15 b™—-15
b(2) — 0 (3-25)
Les relations (3-21), (3-23), (3-24) et (3-25) permettent ainsi de calculer les composantes du
tenseur d'anisotropie en turbulence homogene a bulles uniformement cisaillee en fonction du
tenseur equivalent monophasique a condition de connaitre les rapports entre les echelles de temps
Tm-) et —, ainsi que le rapport entre les parties turbulente et pseudo-turbulente de l'energie
Tt Tb
k
d'agitation ^.
Nous ecrivons le rapport des contributions turbulente et pseudo-turbulente dans l'energie des
fluctuations turbulentes sous la forme suivante :
2
ks 1 u — 1 aCC R» Tt
k0
——
^
a- k0 , R T5 (0) T
(3-26)
4
et Lb
la composante non diagonale b12 du tenseur d'anisotropie se mettent alors sous les formes
suivantes :
1 +1 a uR 1 +1 aC,,C Reb Tt
V_ — 4 kn 4 R r (0) t
et b
(3-27)
Vtn
t0 ,1 +
. a Tt 1 + aTi
Tb Tb
115
CK4R/TRE 3 MODELEE DE TURBULEACE ^ TRO/E EgU^T/OAE ROUR LEE ECOULEMEATE ^ BULLEE
1 + 4 aC , C ]R^ ^
r
b(T^ Tt
et b (3-28)
b(m) T(m)
U12 Lt
1+a
Tt
ib
R4
La relation (3-27) montre que selon que la difference (—(0----------- ) est positive ou negative, le
V C,CR
k2
v
rapport —- est superieur ou inferieur a l'unite. Il faut cependant noter que vt0 — C, — ne
vt0 e0
represente pas la viscosite turbulente de l'ecoulement monophasique equivalent, mais plutot la
viscosite turbulente qu'aurait l'ecoulement monophasique soumis au meme cisaillement et ayant le
montant d'energie turbulente k0 et le meme taux de dissipation. Nous soulignons egalement que
(T)
b 12
l'expression (3-28) ne permet pas de deduire directement le sens de variation du rapport x(m)
12
tres faibles taux de vide pour exprimer les echelles qui entrent en jeu dans les relations (3-23), (3
24), (3-27) et (3-28). Dans ce cas, nous pouvons identifier le champ associe a la partie turbulente
du tenseur de Reynolds (avec l'indice (0)) au champ monophasique. Cela signifie que la presence de
la phase dispersee a tres faibles taux de vide engendre une modification de la structure de la
turbulence du liquide tout en preservant les echelles caracteristiques de la turbulence cisaillee
T(m) — (3-30)
t
116
CK4R/TRE 3 MODELEE DE TURBULEACE ^ TRO/E EgU^T/OAE ROUR LEE ECOULEMEATE ^ BULLEE
Tt Tt 1 S,
(3-31)
Tb Tb S
ou Sb — —- est un taux de cisaillement caracteristique de la turbulence dans les sillages des bulles.
d
Dans ces conditions, l'expression (3-27) de la viscosite turbulente ainsi que les relations (3-23) et
(3-24) correspondant au tenseur d'anisotropie peuvent s'ecrire sous les formes suivantes :
1 + ^ a R=A
4 R^ S
,(m) (3-32)
v^ ‘ a^
1 +---- W a Sb
cjc, S
1 + ^ a -Ri^5b
+ ’ ^ a^eb U’b
b(T) y 4 R^:') S
(3-33)
,(m) 1
b12 1 +-----W a^
S
b lT' — 1
(3-34)
b im' b 2^) 1 S
1 + ------
C^/C,
:=a —S
En particulier, la relation (3-33) indique que la presence des bulles a faibles taux de vide peut selon
les cas, provoquer l'attenuation ou l'augmentation de la composante non diagonale b12 du tenseur
facteur dont depend la composante b12 pour un taux de cisaillement et un taux de vide donnes est
117
CK4R/TRE 3 MODELEE DE TURBULEACE ^ TRO/E EgU^T/OAE ROUR LEE ECOULEMEATE ^ BULLEE
R eb
le rapport Selon que ce rapport est inferieur ou superieur a une valeur seuil qui vaut
R (m)
et
4
67 , l'equation indique respectivement une attenuation ou une augmentation de b12.
c^c
Pour verifier ces resultats on se propose de calculer cette composante dans le cas de l'experience de
Lance et al. (1991). Rappelons que dans ces experiences une attenuation du frottement turbulent a
ete mesuree en comparaison avec l'ecoulement monophasique. En particulier, pour l'ecoulement a
bulles avec le taux de cisaillement de 2.9 s-1, le nombre de Reynolds de l'ecoulement equivalent
monophasique R^ est de l'ordre de 150, cette valeur est deduite des simulations numeriques. Si
on considere que les diametres des bulles sont comprises dans l'intervalle 3-5 mm avec une vitesse
relative comprise entre 0.2 et 0.4 m/s, le nombre de Reynolds associe au mouvement des bulles sera
R eb
compris entre 600 et 2000 ce qui correspond a un rapport compris entre 4 et 14 qui reste
R (m)
et
4
assez faible par rapport a la valeur 67 .
c^c
Pour les tres faibles taux de vide, le developpement limite au premier ordre de la relation (3-32)
conduit a la simplification suivante de l'expression de la viscosite turbulente :
v t »v +a ^
1
4
Reb -
' CR^
R (m)
et (3-35)
t t S
Cette approximation a tres faibles taux de vide est tres interessante : elle permet une interpretation
aisee de l'effet des interactions interfaciales sur la structure de la turbulence dans cette situation
asymptotique. En effet, elle indique que pour un taux de vide et un taux de cisaillement donnes, la
deviation de la viscosite turbulente en presence des bulles par rapport a celle de l'ecoulement
monophasique equivalent s'exprime en fonction de la difference entre les nombres de Reynolds
Reb et Rem . Cette relation indique que l'effet de la presence des bulles sur la turbulence peut etre
interprete comme une forme de substitution a l'echelle des bulles de la turbulence monophasique
par la turbulence induite par les bulles. Tout se passe comme si la bulle dans son mouvement relatif
entraine une turbulence dans son sillage qui vient se substituer a la turbulence de l'ecoulement
118
CK4R/TRE 3 MODELEE DE TURBULEACE ^ TRO/E EgU^T/OAE ROUR LEE ECOULEMEATE ^ BULLEE
"monophasique". Cette turbulence modifie la viscosite turbulente qui selon les cas va etre plus
importante ou plus faible que la viscosite turbulente "monophasique".
La relation (3-35) montre, qu'a faibles taux de vide, la deviation de la viscosite turbulente par
S
rapport a la situation monophasique depend lineairement du rapport . Rappelons que ce rapport
Notons enfin que pour des ecoulements bulles gouvernes par la pseudo-turbulence de telle sorte
que nous pouvons negliger R^ devant Reb, la relation (3-49) se simplifie encore sous la forme :
c S
vt — v +v b avec v„ — v)-S^a Reb (3-36)
des fluctuations turbulentes produites par cisaillement dans l'ecoulement diphasique. En identifiant
ces deux echelles pour les tres faibles taux de vide, nous avons pu etablir des relations simples entre
les composantes du tenseur d'anisotropie de l'ecoulement diphasique et de l'ecoulement
119
C%4f/IRE 3 MODEEEE DE IDREDEEECE ^ IRO/E EgD^I/OEE EODR EES' ECODEEMEEIE ^ EDEEEE
relations entre les echelles caracteristiques de l'ecoulement monophasique t(m) et l'echelle de temps
t1 caracteristique des fluctuations turbulentes produites par cisaillement dans la phase continue de
l'ecoulement diphasique
^ t2
(1 + 7^) (—)2 -a() -1 = 0
(3-37)
t (m)2
k
tb tb
0 t
k^ t,
a + a2 + 4(1 + -^)
_t_L =__ ^ k0
2
(3-38)
tb 2(1 + ^) t
k 0 t (m)2
t
Nous verifions que la relation (3-38) admet comme valeur asymptotique lorsque le taux de vide
tend vers zero 11 = ttm). Les relations (3-21), (3-23), (3-24) et (3-25) permettent ainsi, a l'aide des
d'agitation —. Ce rapport est calcule a partir des donnees experimentales : l'energie turbulente est
k„
120
CHTR/IRE 3 MODEZEE DE IURgUEEECE ^ IRO/E EgU^I/OEE ROUR EEE ECOUEEMEEIE ^ gUEEEE
Les resultats de cette analyse sont recapitules dans le tableau (3-1) ou les composantes du tenseur
d'anisotropie sont calculees dans le cas d'un ecoulement a bulles uniformement cisaille pour deux
taux de vide (Bellakhal et al 2004 a). Ces resultats sont compares aux valeurs des composantes du
tenseur d'anisotropie deduites des resultats experimentaux de Lance et al. (1991) obtenus en
situation d'equilibre (voir paragraphe 2.3 au chapitre 2).
Les resultats presentes dans le tableau (3-1) indiquent un accord satisfaisant entre les resultas
analytiques developpes a partie des modelisations de la turbulence et les resultats experimentaux.
L'interpretation des modelisations au second ordre en terme d'echelles ramenees a des grandeurs
primaires de l'ecoulement permet de reproduire la modification de la structure de la turbulence
induite par la presence des bulles en ecoulement de turbulence homogene uniformement cisaillee.
Ces relations peuvent etre utilisees pour specifier les conditions aux limites des modeles de
turbulence ou il est necessaire de specifier les contributions turbulente et pseudo-turbulente du
tenseur de Reynolds.
121
CH4P/7RE 3 MODEZEE DE IURgUZEECE ^ IRO/E EgU4I/OEE ROUR ZEE ECOUZEMEE7E ^ gUZZEE
2
, 1
kS = — a u, (3-40)
S 4
122
CK4R/IRE 3 MODEZEE DE IURgUZEECE ^ IRO/E EgU^I/OEE ROUR ZEE ECOUZEMEEIE ^ gUZZEE
Avec :
(1+Ik1) (1+Ik1)
kn k0
nt = Cm = n t0
e0 (1 + a^^) T
(1 + a—L)
Tb Tb
^=(T)
- la diffusivite turbulente tt uV exprimee suivant le modele de Launder et al (1975) avec le
^=(T)
. , A A trace(u'u' )
temps de retournement tt =------ --------- -
t 2e0
- la diffusivite associee aux mouvements des bulles que nous proposons de modeliser sous la forme
=(S)
tensorielle tb u'u' avec le temps caracteristique Tb associe au mouvement relatif des bulles,
d
defini par : tb = CR
uT
Le modele des termes de transport diffusif qui generalise ainsi le modele de Launder et al (1975)
aux ecoulements a bulles s’exprime sous la forme:
C -(T) T(S)
SY
Diff(Y) = (1 -a)(T t u u + Tb u u )VY (3-42)
1 -a
123
CK4R//RE 3 MODEZEE DE IURgUZEECE ^ IRO/E EgU^I/O^E ROUR ZEE ECOUZEMEmE ^ gUZZEE
n
ou CsfTtk0 = —5— est une diffusivite turbulente equivalente a celle de l'ecoulement monophasique
sf
et CsfT bks =~^ ^ad u est une diffusivite induite par le deplacement des bulles. En
sf —
2
exprimant l'energie cinetique induite par les bulles sous la forme kS = 1 a
u la diffusivite
induite par les bulles apparait comme une generalisation du modele de Sato et al (1981)
Avec cette modelisation des termes de diffusion, les equations de transport de l'energie des
fluctuations (parties turbulente et pseudo-turbulente) et du taux de dissipation du modele de
turbulence s'ecrivent sous la forme suivante, elles constituent les equations du modele k - e a trois
2
C sk d dk« 1D
kS [(1 -a)(T tk0 +T bkS^^"^] + "%:- a (3-46)
Dt (1 -a) dx: dx, 4 Dt
D du'(du-+dj - e„s 0
—e 0 =------------- [(1 -a)(T tk0 +T bkS^"^] + 7"°" C en, (3-47)
Dt (1 - a) dx j dxj k0 dxj dxj dx
Dans le cas monophasique, le modele de turbulence k-ro developpe par Wilcocx (1988) utilise
pour seconde grandeur caracteristique de la turbulence une frequence caracteristique m calculee a
partir d'une equation de transport. La viscosite turbulente est definie par la relation suivante :
124
CH4P/7RE 3 MODEZEE DE IURgUZEECE ^ IRO/E EgU4I/O^E ROUR ZEE ECOUZEMEAIE ^ gUZZEE
k
nt (3-48)
w
e = wk (3-49)
L'equation (3-49) implique que la frequence w peut etre interpretee comme une frequence
caracteristique de la decroissance de la turbulence.
modele propose par Chahed (1999) pour exprimer la partie non lineaire du terme de redistribution
s'ecrit :
r(T 2 r(T 2
F ijNL) =-C1(t t1 +aT b1)(u iu j -3k0d ij) = F (jNL) =-C1(w 0 + aw b)(u iu j -3k0d ij) (3-50)
w = w0 +aw b (3-51)
Cette decomposition permet de donner une signification physique pour la seconde echelle de temps
introduite dans les fermetures au second ordre. Rappelons que l'echelle de temps introduite dans la
d
Tb C (3-52)
u
ou CR est l'unique coefficient supplementaire introduit dans le modele de turbulence pour les
ecoulements a bulles. Il a ete ajuste a partir des experiences de l'ecoulement homogene a bulles
125
CH4R/7RE 3 MODEZEE DE IURgUZEECE 4 IRO/E EgU4I/O^E ROUR ZEE ECOUZEMEAIE 4 gUZZEE
caracteristique de la turbulence induite par les bulles est interpretee par Lopez de Bertodano et al
(1994) comme le rapport de la dissipation de la turbulence dans les sillages des bulles par la
puissance developpee par la force de trainee, ce qui permet d'exprimer wb par l'expression
suivante :
3 C
wb (3-53)
2 CA d
Il est interessant de remarquer qu'avec les valeurs des coefficients CA = 0.5 (valeur communement
utilisee (Rivero 1991)) et que pour des bulles millimetriques CD = 0.44 (Wallis et al 1964), la
frequence wb s'ecrit :
w b = 1.32 b (3-54)
d
On retrouve ainsi une valeur tres proche de la valeur du coefficient CR ajuste a partir des donnees
experimentales.
1 + k^
k = k0 k0
vt - - (3-55)
w w 01 + a^
w0
Deux principaux resultats peuvent etre deduits a partir de cette formulation de la viscosite
turbulente : d'une part, nous retrouvons un resultat similaire a celui obtenu par reduction des
fermetures au second ordre mais d'une maniere plus "naturelle" : il n'est pas necessaire de
considerer la situation particuliere de la turbulence homogene a bulles uniformement cisaillee pour
126
CH4R/7RE 3 MODEZEE DE IURgUZE^CE 4 IRO/E EgU4I/O^E ROUR ZEE ECOUZEMERIE 4 gUZZEE
e 0 =w0kg (3-56,)
et de son taux de dissipation (3-41) permettent, compte tenu de (3-56), d'etablir l'equation de
Dw0 6 u (3ui + 6u
[C,, -1] ")-(C e2 - w0 (3-57)
Dt 6^: 6^: 6x
Dans le cas du modele de turbulence k-s diphasique, les equations de transport de l'energie
turbulente k0 et de son taux de dissipation e0 (equations 3-45 et 3-47) ont les memes formes que
Les equations de transport du modele de turbulence k-ro diphasique que nous proposons s'ecrivent
ainsi :
2
C sk k
6 (1 -a)(-^ k )^]
6k + -—
1D
kS + -^ a u (3-59)
Dt (1 -a) 6x. w 0 w b 6xj 4 Dt
127
CK4R/IRE 3 MODEZEE DE IUREUZE^CE ^ IRO/E EgU^I/O^E ROUR ZEE ECOUZEMERIE ^ gUZZEE
avec la formulation (3-55) de la viscosite turbulente et l'expression (3-53) de la frequence wb. Les
Les modeles de turbulence a trois equations (k-e et k-w) developpes dans ce chapitre sont
formules ci-dessous pour des ecoulements homogenes a bulles, uniforme et a cisaillement constant,
dans un plan vertical oxy (Bellakhal et al 2004 a, b). On note u et v les composantes longitudinale
(direction x) et horizontale (direction y) de la vitesse moyenne du liquide et u' et v' les fluctuations
correspondantes. Dans ces deux ecoulements la composante transversale de la vitesse du liquide est
composante longitudinale est a glissement constant par rapport a phase continue egale a la vitesse
limite.
Le tenseur de Reynolds comprend une partie turbulente et une partie pseudo-turbulente soit avec
les hypotheses de l’ecoulement homogene et conformement aux principes de modelisation on a :
= (T)
2
= (S) (T)
2
(S) ^=(S) = = (T) -(T)
u'2 = u'2 +u v'2 = v'2 +v ; u'v' = 0 et u' =v (3-61)
Selon les fermetures proposees, les equations des modeles de turbulence a trois equations pour les
ecoulements de turbulence homogenes a bulles s'ecrivent :
128
CK4P/IRE 3 MODEEEE DE IDREDEEECE ^ IRO/E EQD4I/OEE EODR EES' ECODEEMEEIE ^ EDEEEE
Modele k-e :
_D
k0 vtS2 -e 0 (3-62)
Dt
_D 1_D_ (a| fiR|2)
(3-63)
kS
Dt 4 Dt
Modele k-rn :
D
kr = v ,S2-p'k„ w 0 (3-65)
Dt
4D (a t.i2)
_D (3-66)
kS
Dt
D = i s2-pw 2
—w 0 (3-67)
Dt
Avec S le taux de cisaillement nul dans le cas de l’ecoulement uniforme de grille ; il est constant
dans le cas de l’ecoulement a cisaillement constant. Les composantes du tenseur de Reynolds sont
deduites a partir des relations de fermeture (3-13) et (3-14).
Les equations des modeles sont discretisees par la methode des differences finies et integrees en
mode parabolique : Le domaine continu (x,y) est remplace par un maillage discret dont nous avons
choisi qu’il soit irregulier dans la direction de l’ecoulement (le maillage est resserre dans la premiere
zone de l’ecoulement). Les resultats numeriques sont calcules de proche en proche en allant de
l’amont vers l’aval.
Les deux modeles a trois equations (k,e et k,w) sont equivalents en ecoulement de turbulence
homogene de grille, ils reproduisent les memes resultats numeriques. Dans les simulations
d’ecoulement de grille, les bulles sont injectees avec un glissement nul au niveau de la grille ou la
turbulence est supposee etre celle de l’ecoulement monophasique. La turbulence au droit de la
grille a ete ajustee, par le calage des valeurs de l’energie turbulente et du taux de dissipation (de la
129
CKFP/IRE 3 MODEEEE DE IDREDEEECE 4 IRO/E EQD4I/OEE EODR EES' ECODEEMEEIE 4 EDEEEE
frequence w0 pour le modele k-ro) a l’entree de l’ecoulement, afin de reproduire les resultats des
experiences realisees en monophasique. Sur la figure (3-1) on observe une bonne validation du
calcul de l’intensite turbulente en aval de la grille dans le cas monophasique et trois ecoulements a
taux de vide differents (vitesse moyenne 0.6 m/s, taille de la maille M = 0.04 m) : le glissement
genere par le bilan moyen de quantite de mouvement dans le gaz permet alors de reproduire la
turbulence interfaciale, de sorte que le niveau de la turbulence totale soit correctement predit, ainsi
que sa decroissance pour les differents taux de vide
0,1
------ alpha=0 n experience
■ ■ ■ alpha=0,005 a experience
0,08 ------alpha=0,01 ° experience
------ alpha=0,02 ♦ experience
0,06
0,04
0,02
20 40 60
X/M
Figure (3-1) : Decroissance de l’intensite turbulent^vu'2 /U derriere une grille. Comparaison des
resultats des deux modeles avec les donnees experimentales de Lance et Bataille (1991),
(M = 4cm).
L’intensite turbulente mesuree a une distance donnee de la grille depend fortement du taux de vide
et de la vitesse de l’ecoulement comme le montre la figure (3-2).
130
CKFP/IRE 3 MODEEEE DE IDREDEEECE 4 IRO/E EQD4I/OEE EODR EES' ECODEEMEEIE 4 EDEEEE
0,075
0,05
— modele » experience
0,025
Figure (3-2) : Intensite turbulente Vu'2 /U en fonction du taux de vide a pour differentes valeurs
de la vitesse moyenne. Comparaison des resultats des deux modeles avec les donnees
experimentales a X/M =36.4 de Lance et Bataille (1991).
Ces experiences suggerent que l’effet des bulles sur la turbulence du liquide est d’autant plus grand
que la pseudo-turbulence est importante vis-a-vis de la turbulence engendree par la grille.
La difference entre l’ecart type des fluctuations de vitesse transversale et verticale evaluee
experimentalement par Lance et Bataille (1991) reste de meme ordre de grandeur que l’incertitude
de mesure (de l’ ordre de 10%) et les composantes non diagonales du tenseur de Reynolds sont
negligeables, ce qui a mene ces auteurs a deduire que la presence des bulles n’altere pas l’isotropie
initiale de l’ecoulement monophasique. Ce n'est que dans des experiences ulterieures (Mareuge et
Lance, 1995) qu'on a mis en evidence le caractere anisotrope de la pseudo-turbulence, dans
lesquelles on mesure les fluctuations turbulentes induites par la montee de bulles dans un liquide au
repos (voir figure (2-26) dans le paragraphe 2.3.2.2 au chapitre 2). Les resultats obtenus indiquent
que le rapport entre l’ecart type des fluctuations de vitesse transversale et verticale est de l’ordre de
0.8
La figure (3-3) montre que les resultats numeriques indiquent une anisotropie en ecoulement
diphasique, anisotropie due au modele de pseudo-turbulence et qui tend a augmenter avec le taux
de vide.
131
CKFP/IRE 3 MODEEEE DE IDREDEEECE 4 IRO/E EQD4I/OEE EODR EES' ECODEEMEEIE 4 EDEEEE
Figure (3-3) : Evolution du rapport des composantes normales du tenseur de Reynolds en fonction
du taux de vide pour differentes valeurs de la vitesse moyenne. Comparaison des resultats des deux
modeles a X/M =36.4 avec les donnees experimentales de Lance et Bataille (1991).
Pour les deux modeles, on procede comme pour l’ecoulement homogene de grille par la
simulation du cas monophasique en ajustant des valeurs de l’energie turbulente et, d'une part du
taux de dissipation pour le modele k-e, d'autre part de la frequence caracteristique w0 pour le
modele k - w0, a l’entree de l’ecoulement afin de reproduire les mesures experimentales relatives a
ce dernier. Les bulles sont egalement injectees au droit de la grille sans vitesse de glissement et la
turbulence a cet endroit est supposee etre celle du cas monophasique.
La figure (3-4) montre le meilleur ajustement qui reproduit l'ecoulement monophasique. Les
figures (3-5), (3-6) et (3-7) temoignent aussi d’un accord satisfaisant entre simulations et
observations dans trois cas d’ecoulements a bulles uniformement cisailles. Le modele rend bien
compte des differences observees par rapport au cisaillement uniforme monophasique. Il predit bien
notamment l’augmentation de l’isotropie dans les ecoulements a bulles qui se traduit d’une part par
132
CKFP/IRE 3 MODEEEE DE IDREDEEECE 4 IRO/E EQD4I/OEE EODR EES' ECODEEMEEIE 4 EDEEEE
0,1
— u'/U ■ experience
0,08
— - - v'/U & experience
— -sqrt(u'v')/U • experience
0,04
0,02
SX/U
0,1
— u'/U ■ experience
0,08
■ - - v'/U a experience
.....-sqrt(u'v')/U • experience
c 0,04
0,02
SX/U
133
CKFP/IRE 3 MODEEEE DE IDREDEEECE 4 IRO/E EQD4I/OEE EODR EES' ECODEEMEEIE 4 EDEEEE
0,04
0,02
SX/U
— u'/U experience
0,08
- - v'/U experience
experience
0,02 -
SX/U
Ce comportement est bien mis en evidence en comparant la figure (3-8), relative a l’ecoulement
monophasique, qui presente l’evolution le long de la conduite des composantes du tenseur de
134
CKFP/IRE 3 MODEEEE DE IDREDEEECE ^ IRO/E EgD^I/OEE EODR EES' ECODEEMEEIE ^ EDEEEE
Reynolds normalisees par l’energie turbulente, aux figures (3-9), (3-10) et (3-11) relatives au cas
diphasique pour differents taux de vide.
Figure (3-8) : Evolution longitudinale des composantes du tenseur de Reynolds normalisees par
l’energie turbulente k dans le cas monophasique. Comparaison avec les donnees experimentales de
Lance et al (1991) (U= 1. m/s, a = 0 %, S = 2.9 s-1).
SX/U
Figure (3-9) : Evolution longitudinale des composantes du tenseur de Reynolds normalisees par
l’energie turbulente k dans le cas d’ecoulement a bulles. Comparaison avec les mesures
experimentales de Lance et al (1991) (U= 1. m/s, a = 1 %, S = 2.9 s-1).
135
CKFP/IRE 3 MODEEEE DE IDREDEEECE 4 IRO/E EQD4I/OEE EODR EES' ECODEEMEEIE 4 EDEEEE
0 2 SX/U 4 6
Figure (3-10) : Evolution longitudinale des composantes du tenseur de Reynolds normalises par
l’energie turbulente k dans le cas d’ecoulement a bulles. Comparaison avec les mesures
experimental de Lance et al (1991) (U= 1. m/s, a = 1.4 %, S = 2.9 s-1).
1,2
------ u'u'/k (modele) ♦ experience
0,6
0,3
SX/U
Figure (3-11) : Evolution longitudinale des composantes du tenseur de Reynolds normalisees par
l’energie turbulente k dans le cas d’ecoulement a bulles. Comparaison avec les mesures
experimentales de Lance et al [4] (U= 1. m/s, a = 2 %, S = 2.9 s-1).
136
CKFP/IRE 3 MODEEEE DE IDREDEEECE ^ IRO/E EgD^I/OEE EODR EES' ECODEEMEEIE ^ EDEEEE
On retrouve des resultats similaires avec le modele k - w (Bellalhal et al 2004 b). La figure (3
— u'/U ■ experience
0,08
■ - v'/U & experience
-sqrt(u'v')/U • experience
0,04 -
0,02 -
SX/U
Les figures (3-13) - (3-15), semblables aux figures (3-5) - (3-7), presentent revolution
longitudinale des composantes du tenseur de Reynolds reproduite par le modele k - w pour les
trois taux de vide de 1 %, 1.4 % et 2 %.
— u'/U ■ experience
0,08
- - - v'/U A experience
.....-sqrt(u'v')/U • experience
0,04 -
0,02 -
SX/U
137
CKFP/IRE 3 MODEEEE DE IDREDEEECE 4 IRO/E EQD4I/OEE EODR EES' ECODEEMEEIE 4 EDEEEE
0,1
— u'/U experience
0,08
experience
experience
0,04
0,02
SX/U
0,1
-sqrt(u'v')/U experience
* " ff —
■ ■ ■?-----
0,02 -
SX/U
138
C%4f/IRE 3 MODEEEE DE IDRBDEEECE ^ IRO/E EgD^I/OEE RODR EEE ECODEEMEEIE ^ RDEEEE
3.5 Conclusions
equations de transport (k0 -kS -e) et (k0 -kS - w) ont ete mis en reuvre dans le cas de deux
139
C%4f//RE 3 MODEEEE DE IDRRDEEECE ^ IRO/E EgD^I/OEE RODR EEE ECODEEMEEIE ^ RDEEEE
140
CH&P/IRE 4 CODCEEE C/EX/EEEEEM/ECEE 2-D EE ECODEEMEEI ^ EDEEEE
4.1 Introduction
Nous avons developpe, dans le chapitre 3 des modeles de turbulence au premier ordre (k - e et
k - w) pour les ecoulements homogenes a bulles a faibles taux de vide. Ces modeles ont ete
deduits par reduction des fermetures au second ordre ; en particulier cette reduction a permis
d'obtenir une formulation originale de la viscosite turbulente diphasique qui conserve, a travers les
echelles quelle met en jeu, l'aptitude de cette formulation a rendre compte des effets de la
turbulence et de la pseudo-turbulence sur le frottement turbulent. L'application de ces modeles de
turbulence dans deux cas d'ecoulements homogenes a bulles (uniforme et a cisaillement constant) a
montre une concordance satisfaisante avec les observations experimentales. Les modeles ont
reproduit, en particulier dans le cas de la turbulence homogene a cisaillement constant, la reduction
du frottement turbulent et de l'anisotropie de l'ecoulement observee experimentalement.
L'etude de la turbulence homogene permet d'eviter les difficultes liees d'une part a la repartition
des phases puisque le taux de repartition des phases est uniformement distribue et d'autre part a la
modification du champ de la vitesse moyenne dans le liquide par retroaction des fluctuations
turbulentes qui comprend, en ecoulement a bulles, celles induites par l'effet des interactions
interfaciales. Nous nous proposons dans ce chapitre d'appliquer le modele a deux fluides avec
fermeture au premier ordre de la turbulence aux ecoulements inhomogenes a faibles taux de vide.
Le modele a deux fluides sera mis en reuvre dans la simulation d'un ecoulement diphasique
bidimensionnel turbulent de couche cisaillee mince (le sillage vertical a bulles) pour lequel on
dispose de donnees experimentales completes (Roig 1993). Cet ecoulement a par ailleurs fait l'objet
d'analyse a l'aide d'un modele a deux fluides avec fermeture au second ordre de la turbulence
(Chahed et al 2002).
141
CH&P/IRE 4 CODCEEE C/E4/EEEEE M/ECEE 2-D EE ECODEEMEEI ^ EDEEEE
Apres avoir rappele les equations dynamiques, les fermetures du terme de transfert interfacial de
quantite de mouvement et les equations du modele de turbulence, nous presentons les applications
numeriques et les resultats des simulations du sillage a bulles. Ces resultats seront commentes vis-a
vis du comportement auto-preservatif du sillage a bulles et des proprietes qui en decoulent.
Les equations de bilans moyens de masse et de quantite de mouvement dans les deux modeles
euleriens s'ecrivent :
5a1 5 . —.
~5T + ^^(a1uu) = 0 (4-1)
Du 5 5 :
p1a^:L = a^™s ij - p^^[a1uLiuLj] + pl(1 -a)gi + M1i (4-3)
Dt 5x 5x■
d^2 5 = 5
P 2a 2 = a 2^" Sij -P 2—[a 2 u 2iu 2j] + P 2 a 2gi + M2i (4-4)
dt 5x^ 5x::
On designe par l'indice "1" la phase continue et par l'indice "2" la phase dispersee. Les derivees
D 5=5 d 5=5
— = —+ u,j et — 5t + u2i 5xi designent les operateurs de derivees en
Dt 5t 5xi dt
tenseur des contraintes dans le liquide. En adoptant l'approche fondee sur la notion de densite
locale de force (voir paragraphe 1.3.3.3 au chapitre 1), les termes interfaciaux de transfert de
quantite de mouvement qui figurent dans ces bilans sont donnes par la relation :
142
CK4P/IRE 4 CODCEEE C/E4/EEEEE M/ECEE 2-D EE ECODEEMEET ^ EDEEEE
ou s[j( ) est la composante fluctuante du tenseur des contraintes exercees par l'ecoulement non
Nous designons dans ce paragraphe la phase gaz dispersee correspondant aux bulles par l'indice
G et la phase liquide par l'indice L.
Le terme interfacial adopte dans ce modele s'ecrit MGl =< c Gf® > ou f]P’) est donnee par la
relation suivante :
uB
(»„ ,. = 5s j
fPi’ =uL IIs (("n,ds
(4-6)
B 9un 5x,
avec sj est le tenseur des contraintes associe a l'ecoulement perturbe, uB et 5uB designent
uB
respectivement le volume et la surface d'une bulle et B est un operateur de moyenne volumique
sur le volume d'une bulle. La densite fp1) comprend les contributions des forces de trainee, de
masse ajoutee et de portance. Son expression lagrangienne retenue s'ecrit comme suit (Rivero et al
1991):
fPi) =- 4 P LT .(0)
(uGi - uLi)) - CApL(dtuGi - Dt uLi)) - 2pLCLwLij) X (uGj - uLLj)) (4-7)
143
CfMP/7%E -/ 2-D EVECOr/LEta'AT^ E&/Z2EE
avec les derivees particulaires — = — + uT -^— et — = —+ u_-^-. Dans cette expression, les
Dt a ^ 6Xj dt a ^ 6Xj
champs de vitesse uRf et de vorticite to,jj' correspondent a l'ecoulement non perturb e, d est le
diametre des bulles et (CD ; CL ; CA) sont respectivement les coefficients de trainee, de masse
ajoutee et de portance.
Nous retenons dans ce travail, en suivant en cela Chahed et al (2002), la formulation suivante du
terme interfacial de transfert de quantite de mouvement qui inclue les contributions des champs
moyens dans les forces de trainee et de portance ainsi que les contributions des champs moyen et
turbulent dans la force de masse ajoutee:
CT d = D =
Mg, = -ap. ut u Rl 2ap L C L co LlJ xu rj ~ aPLCA(—uGl uLl)
dt Dt
(4-8)
5
-PLCA—[a(uGluGj -<,<,)]
‘Li^Lj,
Ac
avec a est le taux de presence de la phase gaz (en omettant l'indice G) et uRl est la vitesse
relative donnee par la relation uRl = uGl - uLl . Nous remarquons que dans cette modelisation
144
C%4f/7EE 4 CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECOEEEME^T ^ E^EEEE
trois equations
La turbulence de la phase continue est decrite par le modele de turbulence au premier ordre
(k0 -kS -e) a trois equations developpe dans le chapitre 3. Nous rappelons ici ses equations qui
D C.v t ^ + ^) - C-e 0
—e o =------------- [(1 -a)(t tko + T bks^"^] + r0" (4-11)
Dt (1 - a) 5x j 5Xj ko OXj OXj OX;
-.= C, kL-^
(1+k4
k0
= V t0 avec C ,, = 0.09 (4-12)
e0 (1 + a^) T
(1 + a-L)
Tb Tb
Nous notons enfin que le tenseur de Reynolds est calcule dans ce modele a l'aide de la fermeture de
Boussinesq qui s'ecrit dans ce cas :
145
C%4f/ZRE 4 CO^C^ES C/SX/LLEES MWCES 2-D E^ ECOLLEME^T X EUZEES
Tchen en turbulence homogene (Hinze 1975) en prenant notamment en compte les inhomogeneites
et les effets de croisement des trajectoires (Csanady 1963). Si on note u et v les composantes
longitudinale (selon la direction verticale x) et transversale (direction y) de la vitesse, les tensions
de Reynolds du gaz sont eXprimees en fonction de celles de la phase continue sous la forme :
Les coefficients Cn, C22 et C12 sont formules conformement a la theorie de dispersion de Tchen-
^ + CA
b2 + t r b=-1+Ca. 1 + CA
p
C11 = C22 = C T avec i = ^- et
Tp = (4-15)
1+ T + C
C. Tp ilf' 'R
P
CT est une constante introduite pour tenir compte des pertes de correlation susceptible de se
uG vG ^G
(4-16)
Vt
uLvL
ou vtG est la viscosite turbulente du gaz. Csanady et al (1963) eXprime cette viscosite en fonction
de celle du liquide en tenant compte de l'effet de croisement des trajectoires selon le modele
suivant :
UR
V tG = C Cv avec x
V12 et C b = 0.45 (4-17)
Vt V1+C2
146
CfMP/7%E -/ .' C7&47EEEEEv\EVCEE 2-D EVECOr/LEta'VT.I E&/EEEE
Le code MELODIF utilise une combinaison des methodes de differences finies - volumes finis
basee sur un processus de type pas fractionnaire (Thai Van et al 1994). Cette methode consiste a la
resolution de 1'equation devolution en plusieurs etapes separees correspondant chacune a une
equation elementaire. Des differentes methodes numeriques adaptees sont utilisees pour chacune
des etapes en function de la nature des equations a resoudre. Les differentes etapes de calcul sont
resumees dans la figure (4-1) :
147
CKLP/ZRE 4 COUC^EE C/EX/EEEEEMZMCEE 2-D EV ECOTEEMEAT H EUEEEE
CONDITION LIMITES
PROPRIETES PHYSIQUES
ETAPE DE CONVECTION
ETAPE DE PRESSION
TURBULENCE
SORTIE IMPRESSION
La discretisation des equations est de type volumes finis, elle delink une grille dite de vitesse ou
sont calcules les composantes des vitesses moyennes et les grandeurs turbulentes. Afin d'assurer la
validite de la discretisation des equations de bilan de masse dans les deux phases, la pression et le
taux de vide sont calcules au centre de chaque element de volume. Une nouvelle grille appelee
"grille de pression" decalee par rapport a la premiere est ainsi definie (voir figure (4-2)) :
148
CfMP/7%E -/ .' 2-D EVECOr/LEta'AT^ E&/Z2EE
Le premier terme du second membre represente une correlation entre la fonction de presence de
la phase dispersee et les fluctuations du tenseur des contraintes de l'ecoulement non perturbe. Ce
terme qui participe a la dispersion des bulles est modelise par Simonin et Viollet (1990) comme
suit :
da,
<X2^-^> = PlU>2 (4-18)
OX; dx;
Quant au second terme associe a la contribution de l'ecoulement perturbe, il est modelise par :
La vitesse relative qui intervient dans cette modelisation est donnee par la relation suivante :
(*2URi=(u2j-uJ-<%2u;^> (4-20)
ou le dernier terme du second membre represente la vitesse de drift qui traduit d'un point de vue
physique une correlation entre la fonction de presence de la phase dispersee et les fluctuations du
champ non perturbe. En effet, cette expression est deduite par la prise de la moyenne de la vitesse
relative instantanee qui s'ecrit :
%2URi=%2(u2j-u^) (4-21)
Ce terme de drift est modelise par Simonin (1990) scion la formulation suivante :
149
C%4f/7EE 4 CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECOEEEME^T ^ E^EEEE
avec :
D|2 =T1t2u;iu2i et TjT2 est un temps caracteristique du transit d'une particule dans un tourbillon
Le modele eulerien avec les termes interfaciaux de transfert de quantite de mouvement ainsi
formula necessite une fermeture du tenseur de Reynolds dans le liquide et dans le gaz ainsi que du
tenseur de double correlation binaire fluide-inclusions u1iu,2j . Cette fermeture est realisee a l'aide
t , duk^, + ^ + vt du^l
=-n k[ -+- ] + - [kk +n k-—]5 y pour k = 1 ; 2 (4-23)
dX: 3 dx
ou kk designe les energies turbulentes dans la phase continue (k = 1) et dispersee (k = 2). Les
viscosites turbulentes sont donnees en fonction des echelles caracteristiques de l'ecoulement par les
relations suivantes :
2.
v1 = T k1t
(4-25)
tF 2 b
V2 = (1 + b)v12 -bv1 + (1 -b)-^(jk2 -^3k12) (4-26)
150
CK4P/ZRE 4 CO^C^ES C/SX/LLEES MWCES 2-D E^ ECOLLEME^T X EUZEES
"■ = 2c "7"
(4-28)
L1 = T^| 3 kl (4-29)
t _ L1
T12 = (4-30)
T12 =— [” + C A ] (4-31)
FD P 2
tt2 = 1 (4-32)
V1 + c bV 2
OU V r = est un parametre lie aux effets d'inertie qui agissent sur la particule.
k,
3
La turbulence de la phase continue est predite dans le code MELODIF a l'aide du modele de
turbulence de Simonin (1991). C'est un modele au premier ordre k-s diphasique appartenant a la
classe des modeles a une echelle de temps (voir paragraphe 2.3.1.1 au chapitre 2). Tandis que
k12 = u1iu2i sont deduites a partir des caracteristiques de la turbulence de la phase continue et des
echelles caracteristiques des interactions entre mouvements fluctuants en se basant sur la theorie de
dispersion de Tchen (1947). Ces deux energies sont donnees par :
151
CHXE/TEE 4 COLC^ES C/SX/LLEES MACES 2-D E^ ECOLLEME^T X ELZLES
b2 +hr
k2 = k1 (4-33)
1 + hr
b2 +h r
k12 = 2k1 (4-34)
1 + hr
1 + C,
avec hr = ~^ et b =
12 P 2 +C
P
152
CHXE/TRE 4 COLC^ES C/SX/LLEES MACES 2-D EA ECOLLEMEAT X ELZLES
Les tensions de reynolds dans le gaz seront exprimees en fonction du tenseur de Reynolds dans le
liquide conformement aux relations (4-14) - (4-17).
Les deux versions du code MELODIF presentees plus haut (version initiale et version modifiee)
sont appliquees a la simulation du sillage a bulles et les resultats des simulations sont confrontes
aux resultats experimentaux de Roig (1993). Rappelons a cet effet que la version modifiee
correspond en fait a une implantation des fermetures du modele a deux fluides version ENIT dans
le code MELODIF. Rappelons aussi que le code parabolique avec fermeture au premier ordre
(modele de turbulence a trois equations) est une forme reduite du modele a deux fluides avec
fermeture au second ordre ou l'on modelise des equations de transport distinctes pour les
composantes turbulente et pseudo-turbulente du tenseur de Reynolds. L'application de ce modele
en ecoulement de turbulence homogene avec et sans cisaillement (voir chapitre 3) indiquent que le
modele de turbulence a trois equations fournit des resultats sensiblement identiques a ceux obtenus
a l'aide du modele au second ordre. Dans cette section on veut focaliser l'analyse sur l'aptitude des
fermetures au premier ordre a predeterminer la turbulence dans des situations inhomogenes en
l'appliquant a l'ecoulement de sillage a bulles. A cet effet on adopte la formulation du terme
interfacial et l'ajustement de son expression proposees par Chahed et al (2002) : cet ajustement a
permis de preciser les effets du champ moyen et du champ fluctuant sur l'echange interfacial de
quantite de mouvement et par voie de consequence sur la distribution des phases. Applique a
Dans la suite on presente les resultats numeriques obtenus avec le code MELODIF dans ses
deux versions en ecoulement de sillage a bulles avec un ajustement initial aux donnees
experimentales de Roig (1993). Ensuite, le code MELODIF est mis en reuvre pour produire des
experiences numeriques typiques notamment pour etudier des ecoulements de sillages a bulles ou
153
CHXP/TRE 4 COLCAES C/SX/LLEES M/ACES 2-D EA ECOLLEMEAT X ELLLES
l’on assure des conditions pour que se realise un comportement auto-preservatif en ecoulement a
bulles. Dans toutes les simulations des ecoulements diphasiques, le modele est d'abord teste dans le
cas des ecoulements monophasiques equivalents. Ces resultats seront interpretes vis-a-vis des effets
des bulles sur les proprietes des sillages : en s'appuyant sur la nouvelle formulation de la viscosite
turbulente en ecoulement a bulles, on cherchera a exprimer les lois d'evolution de l'ecoulement
diphasique en fonction de celles de l'ecoulement monophasique equivalent.
D'une maniere analogue a la procedure adoptee dans la mise en reuvre du modele dans le cas des
ecoulements homogenes, nous avons commence par la simulation de l'ecoulement du sillage
monophasique.
4.3.2.1 Maillage
Nous avons teste la sensibilite des resultats numeriques au raffinage du maillage dans les
directions verticale et horizontale dans le cas monophasique en choisissant comme variable de
reference l'energie turbulente. Le maillage retenu du domaine de l'ecoulement pour les simulations
est regulier dans la direction verticale y avec une longueur de maille de 1 cm, et est irregulier dans
la direction horizontale x avec des mailles resserrees au niveau de l'axe du sillage selon la loi
suivante :
154
C%4f/ZRE 4 CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECOEEEME^T ^ E^EEEE
Dx(1) = 1.1 mm
Dx(i) = Dx(i -1)1'1 (4-35)
si Dx(i) > 1.2 alors Dx(i) = Dx(i -1)
Dx(i -1)
Les conditions aux limites sont ajustees a partir des donnees experimentales acquises a 1 cm
avant le debut du sillage. Le profil de la vitesse moyenne du liquide est fixe par interpolation
lineaire du profil experimental. Le profil de l'energie turbulente est deduit a partir des donnees
experimentales relatives a l'intensite turbulente longitudinale du liquide en utilisant une loi deduite
Le profil du taux de dissipation a l'entree est deduit par la relation (4-39) obtenue a partir de
l'hypothese d'equilibre production-dissipation:
3 3
c 4 k2
e = ce —----- ou ce est une constante d'ajustement (4-39)
%x
155
CKFP/ZRE 4 CO^C^ES C/SX/LLEES MWCES 2-D E^ ECOLLEME^T X EUZEES
conditions de symetrie pour les limites verticales exterieures du domaine de l'ecoulement et des
conditions de pression constante a la sortie supposee libre.
Signalons enfin que les profils de la vitesse moyenne et de l'intensite turbulente a la section d'entree
ne sont pas parfaitement symetriques. Cette legere dissymetrie de l'ecoulement monophasique va
s'accentuer en ecoulement diphasique ; on pense qu'elle est due aux conditions d'injection.
Cependant, l'adoption des donnees experimentales dans les conditions aux limites a l'entree permet
implicitement d'introduire cette dissymetrie et nous autorise a confronter dans la suite les resultats
numeriques aux donnees experimentales.
Les figures (4-3) - (4-6) presentent les profils transversaux de la vitesse moyenne et de l'energie
turbulente dans les deux sections X = 20 cm et X = 30 cm. Ces figures temoignent d'une bonne
concordance entre les resultats numeriques et experimentaux.1
0,8
^ 0,6
----------------------■—■—i-----■ ■ " " 1 1 1 i i
g 0,4
— Vitesse monophasique a X = 20 cm ■ experience monophasique
0,2
0
0 5 10 15 20 25 30 35 40
Y (cm)
156
C%4f/7EE 4 CO^C^EE C/EX/EEEEEMZXCEE 2-D E^ ECOEEEME^T ^ E^EEEE
0,005
0,004
— ke monophasique a X = 20 cm
0,003
0,002
0,001
10 15 20 25 30 35 40
Y (cm)
Figure (4-4) : Profil transversal de l'energie turbulente a la section X = 20 cm. Comparaison des
resultats du modele avec les donnees experimentales de Roig (1993).
0,8
^0,6
---------------- ■—■—■ ■ ■ —>./.......................... ■ ■
I
| 0,4
0,2
0
0 5 10 15 20 25 30 35 40
Y (cm)
0,004
— ke monophasique a X = 30 cm
0,003
0,002
0,001
Y (cm)
Figure (4-6) : Profil transversal de l'energie turbulente a la section X = 30 cm. Comparaison des
resultats du modele avec les donnees experimentales de Roig (1993).
157
C%4f/ZRE 4 CO^C^ES C/SX/LLEES MWCES 2-D E^ ECOLLEME^T X EUZEES
Dans les simulations du sillage diphasique, les conditions a l'entree de l'ecoulement sont
egalement fixees a partir des donnees experimentales. Les profils de la vitesse du liquide et du taux
de presence sont interpoles sur le maillage a partir des mesures a la section d'entree tandis que le
profil de la vitesse du gaz est calcule comme la somme de la vitesse du liquide et la vitesse limite de
glissement des bulles. Pour ce qui concerne les grandeurs turbulentes, on identifie la partie
turbulente de l'energie a l'energie turbulente monophasique dans la section d'entree. La partie
pseudo-turbulente de la composante longitudinale du tenseur de Reynolds dans cette section est
calculee comme la difference entre la composante totale mesuree dans le cas du sillage a bulles et
celle mesuree dans l'ecoulement monophasique. Le profil de l'energie pseudo-turbulente a l'entree
est ensuite deduit a l'aide de la solution analytique de Biesheuvel et Van Wijngaarden (1984) etablie
en ecoulement potentiel. Bien que ces approximations des deux contributions turbulente et pseudo-
turbulente a l'entree restent assez intuitives, leur validite a ete prouvee par la concordance des
resultats numeriques obtenus avec les observations experimentales.
La figure (4-7) presente les resultats numeriques et experimentaux de Roig (1993) concernant le
profil transversal du taux de vide a la section X = 0.2 m de l'entree. Cette figure montre une
prediction satisfaisante de la distribution des phases ce qui prouve la bonne evaluation par le
modele du role de la turbulence et des termes interfaciaux dans les bilans de quantite de mouvement
des deux phases.
158
C%4f/7EE 4 COLC^EE C/EX/EEEEEMWCEE 2-D E^ ECOEEEME^T ^ ELEEEE
0,04 -t
0,02
10 15 25 30 35 40
Figure (4-7) : Profil transversal du taux de vide a X=20 cm. Comparaison des resultats du modele
avec les donnees experimentales de Roig (1993)
Dans les figures (4-8) et (4-9) nous confrontons les resultats numeriques et experimentaux de
Roig (1993) concernant les profils transversaux de l'energie turbulente et de la vitesse moyenne du
liquide dans les cas monophasique et a bulles a la meme section X = 0.2 m de l'entree. Nous
presentons egalement dans les figures (4-10) et (4-11) les profils numeriques et experimentaux de
ces grandeurs a la section X = 0.3 m de l'entree. Ces figures montrent une concordance satisfaisante
avec les resultats experimentaux. La figure (4-8) indique une augmentation importante de l'energie
turbulente dans le cas diphasique par rapport a l'ecoulement monophasique equivalent. Ceci traduit
un montant important de l'energie (pseudo-turbulente) induit par la presence des bulles. Dans ces
simulations, la vitesse relative des bulles a l'entree de l'ecoulement est prise egale a sa valeur limite,
elle reste ainsi constante dans tout le domaine de l'ecoulement ; dans ces conditions, le terme de
production dans l'equation de transport de la pseudo-turbulence est nul. Cette partie non
dissipative de l'energie est donc transportee par convection et diffusion turbulente.
Les figures (4-8) - (4-11) indiquent une bonne concordance entre les resultats du modele et
l'experience : ces resultats sont comparables a ceux obtenus avec le code parabolique avec
fermeture au second ordre de la turbulence, Chahed et al (2002). Cela signifie que le modele au
premier (a trois equations de transport) ordre deduit des fermetures au second ordre conserve le
contenu physique des modeles initiaux et reussit a l'aide des nouvelles echelles introduites dans
l'expression de la viscosite turbulente et dans la modelisation du terme de diffusion de rendre
159
C%4f/7EE 4 CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECOEEEME^T ^ E^EEEE
compte des effets des bulles sur la structure de la turbulence dans un ecoulement ou la pseudo
turbulence est importante.
0,014 -r
0,01
I
B 0,008
| 0,006
.22
P
5 0,004
0,002 -
0 -F
0 5 10 15 Y(cm) 20 25 30 35 40
Figure (4-8) : Profils transversaux de l'energie turbulente a X=20 cm dans le cas du sillage
monophasique et a bulles. Comparaison des resultats du modele avec les donnees experimentales de
Roig (1993).
10 15 20 25 30 35 40
Y(cm)
Figure (4-9) : Profils transversaux de la vitesse moyenne du liquide a X=20 cm dans le cas du
sillage monophasique et a bulles. Comparaison des resultats du modele avec les donnees
experimentales de Roig (1993).
160
C%4f/7EE 4 CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECOEEEME^T ^ E^EEEE
0,02
„ 0,015
— kt (X=30 cm) ■ mesures
Y (cm)
Figure (4-10) : Profil transversal de l'energie turbulente totale a X=30 cm dans le cas du sillage a
bulles. Comparaison des resultats du modele avec les donnees experimentales de Roig (1993)
10 15 20 25 30 35 40
Y (cm)
Les figures (4-9) et (4-11) montrent que le modele reussit a reproduire l'augmentation significative
de la largeur du sillage dans le cas diphasique. Cette augmentation de l'ouverture du sillage en
ecoulement a bulles s'interprete par l'accroissement important de la diffusion de quantite de
mouvement en relation avec l'augmentation du frottement turbulent. La nouvelle formulation
161
C%4f/7EE 4 CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECOEEEME^T ^ E^EEEE
proposee pour exprimer la viscosite turbulente en ecoulement a bulles indique en effet, qu'en
ecoulement ou la pseudo-turbulence est importante, l'agitation supplementaire induite par les bulles
se traduit par une augmentation de la viscosite turbulente. Elle rend bien compte de cet effet en
ecoulement de sillage a bulles ou le montant de la pseudo-turbulence est important. Pour illustrer
cela on porte sur la figure (4-12) le profil transversal du rapport des viscosites turbulentes en
ecoulement a bulles et en ecoulement monophasique equivalent a X = 0.2 m. Cette figure indique
que dans la zone ou se developpe le sillage a l'aval de la plaque, la viscosite turbulente en
ecoulement a bulles est 2 a 5 fois sa valeur en ecoulement monophasique equivalent. Avec le meme
taux de cisaillement le frottement turbulent est multiplie d'autant ; ce qui explique l'important
elargissement du sillage en ecoulement a bulles en comparaison avec le profil monophasique.
— nut(diph)/nut(monoph)
5
0
10 15 20 25 30 35 40
Y (cm)
Le modele de turbulence a trois equations permet de determiner deux echelles de temps : une
echelle relative a la turbulence et une echelle specifique des effets induits par les bulles.
L'expression de la viscosite turbulente qui met en jeux ces deux echelles de temps a permis de
rendre compte de la diminution du frottement turbulent en turbulence homogene avec cisaillement
uniforme (voir chapitre 3). Elle permet en ecoulement de sillage a bulle de reproduite l'importante
augmentation du frottement turbulent dans un ecoulement a bulles ou la pseudo-turbulence est
importante.
Dans la version initiale de MELODIF, la viscosite turbulente est generee par un modele k-e a
deux equations et son expression, analogue a celle du modele monophasique ne met en jeu qu'une
162
C%4f/7EE 4 CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECOEEEME^T ^ E^EEEE
seule echelle de temps. Pour comparer les deux versions du code, nous avons egalement calcule le
sillage a bulles a l'aide de la version initiale et les resultats dans les deux sections X = 0.2 m et X =
0.3 m des deux versions sont compares dans les figures (4-13) - (4-16). Notons que nous avons
procede de la meme maniere qu'en monophasique pour fixer les conditions aux limites.
0,8
0,4
o experience diphasique — VL diphasique a X=20 cm
Figure (4-13) : Profils transversaux de la vitesse moyenne du liquide a X=20 cm dans le cas du
sillage monophasique et a bulles. Comparaison des resultats numeriques de la version initiale de
MELODIF et du modele avec les donnees experimentales de Roig (1993).
0,8
UO 0,6
0,4
— VL (Y=30 cm)
0,2 ■ mesures
—VL (MELODIF) a Y=30 cm
10 15 20 25 30 35 40
Y (cm)
Figure (4-14) : Comparaison des profils transversaux numeriques de la vitesse moyenne du liquide
a X=30 cm dans le cas du sillage a bulles obtenus par la version initiale de MELODIF et la nouvelle
version. Confrontation aux donnees experimentales de Roig (1993).
163
C%4f/7EE 4 CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECOEEEME^T ^ E^EEEE
0,014
2 0,008
0,006
0,004
LU
0,002
0
0 5 10 15 Y(cm) 20 25 30 35 40
Figure (4-15) : Profils transversaux de l'energie turbulente a X=20 cm dans le cas du sillage
monophasique et a bulles. Comparaison des resultats numeriques de la version initiale de
MELODIF et du modele avec les donnees experimentales de Roig (1993).
0,02
— kt (X=30 cm)
■ mesures
„ 0,015
— kt (MELODIF) a X = 30 cm
Y (cm)
Figure (4-16) : Comparaison des profils transversaux numeriques de l'energie turbulente du liquide
a X=30 cm dans le cas du sillage a bulles obtenus par la version initiale de MELODIF et la nouvelle
version. Confrontation aux donnees experimentales de Roig (1993).
Les figures (4-13) et (4-14) montrent que la version initiale du code MELODIF produit une
surestimation de l'augmentation de l'ouverture du sillage diphasique en ecoulement a bulle. Ce
resultat peut s'interpreter par une surevaluation de la viscosite turbulente. Correlativement la
164
C%4f/ZRE 4 CO^C^ES C/SX/LLEES MWCES 2-D E^ ECOLLEME^T X EUZEES
diffusivite turbulente est surestimee et le profil de l'energie cinetique turbulente est beaucoup plus
etale que le profil experimental, figures (4-15) et (4-16). Ces figures indiquent d'autre part que le
niveau de l'energie turbulente calcule par la version initiale est inferieur a celui calcule par la version
modifiee. Ce resultat peut s'interpreter par le fait que sans distinction entre la partie turbulente et la
partie induite par les bulles, l'energie totale est supposee participer a la cascade energetique avec la
meme echelle de temps ; elle se trouve par consequent soumise au meme taux de dissipation.
Rappelons enfin que des applications similaires de la version initiale du code ont ete realisees
dans des ecoulements couches de melange a bulles, Roig (1993). Les resultats obtenus ont conduit
a exprimer la viscosite turbulente selon la formulation de Sato et al (1981), c'est a dire comme la
somme de la viscosite turbulente monophasique et d'une turbulence induite par les bulles. En
ameliorant, les resultats du modele, cette demarche produit en revanche une sous estimation de
l'ouverture des couches cisaillees en ecoulement a bulles. En utilisant des fermetures au second
ordre, Chahed (1999) indique en effet que la prise en compte d'une viscosite turbulente induite par
les bulles selon le modele de Sato et al (1981) n'engendre qu'une faible augmentation du frottement
turbulent "monophasique" qui n'explique pas l'importante augmentation de l'ouverture des couches
cisaillees.
165
CfMP/7%E -/ C7&47EEEEEv\EVCEE 2-D EVECOr/LEta'AT^ E&/EEEE
En particulier, a taux de vide constant, le systeme d'equations dynamiques du modele a deux fluides
ecrites dans l'hypothese de l'ecoulement quai-parallele se reduit a une forme similaire au systeme
d'equation de l'ecoulement monophasique qui, par consequent, peut se preter au meme traitement
analytique pour obtenir des solutions affines. La difference entre les deux systemes reside dans
l'expression du frottement turbulent qui en ecoulement a bulle va etre exprimee avec une viscosite
turbulente modifiee par les effets interfaciaux.
166
CK4P/ZRE 4 CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECOEEEME^T 4 E^EEEE
Nous rappelons que les gradients de vitesse moyenne et la turbulence sont confines en
ecoulement de sillage dans une couche mince d'epaisseur l(m) (echelle diffusive) faible vis-a-vis de
l'exterieur du sillage et du deficit de vitesse a l'axe du sillage U(m) = U0 - uaxedu sillage monoph.. Tandis
que les fluctuations turbulentes sont decrites par une echelle des fluctuations u0. On definit le
Une analyse des ordres des grandeurs des differents termes dans les bilans de masse et de quantite
de mouvement montre lorsque les echelles verifient les relations suivantes :
l(m)
® 0 ; Rel ® ¥ (4-40)
L
u i(m)
u0 O(1) et O( (4-41)
)
Um) U L
6u 6v ^
—+— = 0 (4-42)
6x 6y
6u 6u 'v'
U0 (4-43)
6x 6y
_L + v2 = PW
(4-44)
PL P
167
C%4f/ZRE 4 CO^C^ES C/SX/LLEES MWCES 2-D E^ ECOLLEME^T X EUZEES
+^ _ c
M0 = PIu(U - u)dy = ^2 1 pu2 d (4-45)
—¥
place dans un ecoulement uniforme de vitesse U0. On definit l'epaisseur de quantite de mouvement
M0 =pU2 5 (4-46)
La relation (4-45) montre que l'epaisseur de deficit de quantite de mouvement est proportionnelle
au coefficient de trainee :
5 = Cd| (4-47)
On peut interpreter d'un point de vue physique le comportement auto-preservatif du sillage par
la perte de memoire, a partir d'une certaine distance de l'obstacle, des conditions d'entree.
L'ecoulement reproduit une structure identique du champ de vitesse moyenne et de la turbulence
dans les differentes sections transversales dont revolution longitudinale sera decrite par la variation
des echelles caracteristiques.
On s'attend dans ce cas que l'equation (4-43) admette des solutions de la forme :
U0 - u - u'v'
U+ F(h) et G(h) (4-48)
Um)
avec h1 = ~yr
l (m)
: est une coordonnee adimensionnelle transversale.
L'injection des expressions (4-48) dans l'equation (4-43) permet d'etablir que la condition d'affinite
impose les relations suivantes des echelles :
168
CHXE/TEE 4 COLC^ES C/SX/LLEES MACES 2-D E^ ECOLLEME^T X ELZLES
condition suivante :
Cette expression definit le champ de viscosite turbulente comme une propriete de l'ecoulement
La resolution de l'equation (4-43) compte tenu des conditions (4-48), (4-49) et (4-40) permet
d'etablir le profil affine adimensionnel de la vitesse qui s'ecrit :
deficit de quantite de mouvement 5 et du nombre de Reynolds RET par les relations suivantes :
169
C%4f/7EE 4 CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECOEEEME^T ^ E^EEEE
uSm) _ i
rEtN W 2 (4-55)
"uT _ 2V
— F(x)=exp(0,5*x2)
U+ en monophasique
U+ en diphasique
-15 -10 -5 10 15
170
C%4f/ZRE 4 CO^C^ES C/SX/LLEES M/XCES 2-D E^ ECOLLEME^T X EUZEES
Les equations de bilans moyens en ecoulement a bulles (4-1) - (4-4) s'ecrivent dans le cas d'un
ecoulement stationnaire bidimensionnel sous les formes generales suivantes :
Dans le liquide :
= d v,1 ^vv 1
= d dp d d
1 1 u^T + v1 dy
p a _ -a1 ~-p1 T~(a1ul vi ) -p^T(a1vi' ) - My (4-56)
dx dy dx dy
= d u,1 = du dp d d
1 1 u^T + v1
p a _ -a1 ”-p^7(a1u1 )-p^7(a1u1v1)-p1a1g-Mx (4-57)
dx dy dx dx dy
Dans le gaz :
En negligeant l'acceleration et le poids des bulles vis-a-vis des forces qu'exerce sur elles le liquide,
les equations (4-58) et (4-59) se reduisent a :
171
CHXE/TEE 4 COLC^ES C/SX/LLEES MACES 2-D E^ ECOLLEME^T X ELZLES
0 _ -a2 ^ + My
(4-60)
2 dy y
0 _-a 2 + Mx (4-61)
2 dx x
= d v,11
= dv dp d d
1 1 u1^- + v1
p a _ -d"-p1 d~ (a1u1v1)-p1 (a1v12) (4-62)
dx dy dy dx dy
= du,1 ,,,
=du^^1 dp d d
1 1 u^T + v1
p a _ -d~-p1 ”(a1u1 )-p^7(a1u1v1)-p1a1g (4-63)
dx dy dx dx dy
Dans le cadre de l'hypothese d'une repartition uniforme du taux de vide, les equations de bilans dans
le liquide s'ecrivent :
— d v1 — d v1 dP d d
p1a1 u "d^+v1 "dy" _—™—p1~(a1u1v1)-p^~(a1v1 ) (4-64)
dy dx dy
— du1 —du1 dP d d
p1a1 u, —- + v, _ -~-p1 — (a1uT) -p1 d"(a1u1v1) (4-65)
1 dx 1 dy dx dx dy
Avec P _ p + p1a1g est la pression motrice, ces equations sont equivalentes a celles
Si on suppose que les sillages monophasique et a bulles possedent le meme profil de vitesse a
l'entree. Ceci s'interprete d'un point de vue physique par le meme deficit de quantite de mouvement
induit par l'obstacle dans les deux cas. Nous assistons alors dans le sillage a bulles a une
modification des echelles caracteristiques l(d) et u S(d) par rapport a celles de l'ecoulement
172
CfMP/7%E -/ 2-D EVECOr/LEta'VT^ E&/Z2EE
Les simulations de l'ecoulement de sillage a bulles a taux de vide uniforme sont effectuees avec
le meme profil de vitesse a l'entree que pour le sillage monophasique. Les resultats obtenus dans ce
cas montrent que les profils de vitesse ont un comportement auto-preservatif (figure (4-18)) avec, a
differentes sections, les echelles caracteristiques presentees dans le tableau (4-2) :
L'examen des tableaux (4-1) et (4-2) indique que les calculs produisent des valeurs voisines des
produits l(m)Ugm) (en ecoulement monophasique) et l(d)Ugd) (en ecoulement diphasique) ; cela
signifie que le deficit de quantite de mouvement induit par la plaque separatrice est le meme dans
les deux ecoulements. Ce resultat est cependant obtenu avec des valeurs d'ouverture et de vitesse
En sillage diphasique, la propriete d'affinite produite par les resultats numeriques peut s'interpreter
par l'introduction d'une viscosite turbulente diphasique :
y(d) _ ^(m)p
ou E est un coefficient attribue a l'effet des bulles sur la diffusivite de quantite de mouvement. II
s'en suit alors la definition d'un nombre de Reynolds diphasique par la relation :
R(m)
lC=-^ (4-68)
173
C%4f/7EE 4 CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECOEEEME^T ^ E^EEEE
u (d) 1 2 1
Re^ AfWr' (4-70)
- V
u
Compte tenu des equations (4-54), (4-55), (4-56) et (4-70) on deduit les relations :
W=VG et uk 1
Vr
(4-71)
qu'aurait le sillage monophasique ayant les echelles caracteristiques l(d) et USd). La fonction
j(x, y) est une fonction adimensionnelle qui tient compte de la competition entre les mecanismes
k t
de la turbulence et la pseudo-turbulence a travers les rapports — et —1-. Lorsque le comportement
k0 tb
affine est etabli, l'expression de ces rapports (supposes constants a travers la section de
coefficient r .
~du —(0) gu
e0 = -uv -u v — (4-72)
gy gy
avec :
—(0) = n gu
-u v =n (4-73)
gy
174
C%4f/ZRE 4 CO^C^ES C/SX/LLEES MWCES 2-D E^ ECOLLEME^T X EUZEES
1 l(d)
tt = l— ,(d) (4-74)
cm us
d
t b = CR (4-75)
u,
1 = 1 —2
ks = 4auR (4-76)
l(d)U(d)
S
nm — C m kn T + — (4-77)
R ET
ko (4-78)
VCTR ET
Compte tenu des expressions (4-74), (4-75), (4-76) et (4-78) le coefficient G s'ecrit :
c = \
1 + (lV rE«) a u
(d)
4c Us 0
r— (4-79)
a l(d) uR
1+
cTT TuSdy
(d)
u
Sachant que les rapports —(Ry et —— (d)
augmentent avec l'abscisse x ; en se plagant suffisamment
us us
loin du debut du sillage, on peut exprimer le coefficient G, par l'approximation suivante :
C^RET' duR
r (4-80)
4CR l(d)USd)
175
CHXE/TEE 4 COLC^ES C/SX/LLEES MXCES 2-D E^ ECOLLEME^T X ELLLES
Le nombre de Reynolds RET) est ajuste par rapport a l'experience a la valeur 12.5 (Lumley 1973).
En se referant aux experiences de Roig (1993) pour des bulles de diametre moyen de l'ordre de
3 mm et une vitesse relative de l'ordre de 0.4 m/s, on calcule une valeur du coefficient G » 1.5.
l (d) U (m)
S
Le tableau (4-3) presente les rapports (m) et (d) deduits a partir des resultats numeriques dans
l US
l (d) US"'
Sections
l(m) USd)
On deduit une valeur de G telle que VG »1.6, soit alors G » 2.6. Nous remarquons que les
4.5 Conclusion
La mise en reuvre du modele eulerien a deux fluides avec fermeture au premier ordre de la
turbulence (modele a trois equations) dans la simulation d'un ecoulement inhomogene a frontiere
libre (sillage a bulles) a permis de valider les modelisations et de tester la capacite de ce dernier a
reproduire les fermetures dans une situation d'ecoulement a bulles ou la pseudo-turbulence est
importante. Avec des ajustements appropries des conditions d'entree vis a vis des donnees
176
CHXE/TRE 4 COLC^ES C/SX/LLEES MACES 2-D EA ECOLLEMEAT X ELLLES
experimentales, les resultats obtenus dans cette simulation sont comparables a ceux obtenus avec le
code parabolique avec fermeture au second ordre de la turbulence. L'ensemble de ces resultats
conduit aux conclusions suivantes :
• L'accord entre les resultats numeriques et les resultats experimentaux montre que le modele
genere correctement l'augmentation du frottement turbulent induite par l'agitation des bulles via
la formulation de la viscosite turbulente adoptee.
• Le profil du taux de vide montre que le modele reproduit convenablement la repartition des
phases. Cela indique que les termes interfaciaux dans les equations dynamiques peuvent etre
correctement representees en s'appuyant sur les fermetures au premier ordre proposees.
177
CHXP/TRE J . COLCAES C/SX/LLEES M/ACES 2-D EA ECOLLEMEATS X ELLLES
5.1 Introduction
Le modele eulerien a deux fluides avec des fermetures au premier et au second ordre de la
turbulence presente et mis en reuvre dans les chapitres 3 et 4 est developpe pour les ecoulements
turbulents a bulles spheriques a faibles taux de vide ou les interactions hydrodynamiques entre les
bulles sont supposees faibles. Les limitations du modele se situent en premier lieu au niveau de la
formulation des termes interfaciaux dans les equations des champs des vitesses moyennes : la prise
de moyenne du terme interfacial, suppose que le diametre des bulles est relativement petit vis-a-vis
des tailles des structures turbulentes les plus energetiques. Ceci est vrai pour les petites bulles mais
cette hypothese constitue une serieuse limitation pour les bulles de plus grandes tailles. En effet,
l'expression meme de la force instantanee n'est valable que pour des bulles sans deformation
notable. L'autre limitation du modele vis-a-vis du diametre des bulles concerne la modelisation de la
turbulence de la phase continue et de la phase dispersee : il s'agit notamment de la decomposition
lineaire du tenseur de Reynolds et de l'hypothese d'equilibre "production-dissipation" dans les
sillage.
Les forts couplages entre la phase continue et la phase dispersee dans les ecoulements
relativement denses n’autorisent pas les simplifications adoptees et renvoie a la necessite de
developper des modelisations plus generales qui prennent en compte l'effet des interactions
hydrodynamiques dans les lois de fermeture du transfert interfacial de quantite de mouvement et
dans la modelisation de la turbulence.
178
CHXP/TRE J . COLCAES C/SX/LLEES M/ACES 2-D EA ECOLLEMEATS X ELLLES
Les experiences recentes en ecoulements gaz-liquide a forts taux de vide montrent d'importantes
modifications de la structure des ecoulements diphasiques. Dans les ecoulements a forts taux de
presence du gaz (milieu dense), les interactions hydrodynamiques sont fortes et elles ont des effets
importants sur la structure des champs moyens, fluctuant et de distribution des phases dans les
ecoulements diphasiques. En particulier, les experiences realisees en ecoulements a bulles a fort
taux de vide (Zuber et Ischii 1979, Larue de Tournemine 2001, Garnier et al 2002) indiquent que la
vitesse relative du gaz decroit avec l'augmentation du taux de vide. Les experiences de Larue de
Tournemine (2001) et de Garnier et al (2002) montrent d'autre part que l'energie cinetique
turbulente induite par les bulles, qui s'exprime selon la solution potentielle proportionnellement a
aUR, ne suit pas la diminution de la vitesse relative a fort taux de vide. Les mesures de Larue de
Tournemine (2001) semblent indiquer que l'energie turbulente induite par les bulles devrait etre
normee par la vitesse relative determinee par rapport au champs dit "lointain" dont la valeur est
plus importante que la vitesse relative determinee localement (par rapport au champs dit "proche").
Ces resultats indiquent que le modele de l'energie pseudo-turbulente base sur la solution potentielle
et qui attribue la variation de l'energie turbulente induite par les bulles aux seuls effets de la force de
la masse ajoutee est incomplet et il est necessaire de tenir compte de l'effet des interactions
hydrodynamiques ainsi que de l'energie produite au niveau des sillages pour les ecoulements a taux
de vide relativement eleves. Les experiences recentes en ecoulement homogene a forts taux de vide
de Larue de Tournemine (2001) et de Garnier et al (2002) fournissent des donnees remarquables
pour tester de nouvelles idees de modelisation et pour developper des modeles plus generaux.
On se propose dans la suite de batir un modele qui permet une interpretation simple de ces
resultats experimentaux.
179
CHXP/TRE J / COLC^ES C/SX/LLEES M/ACES 2-D EA ECOLLEMEATS X ELLLES
On se place dans le cas d’un reseau de bulles homogene et nous considerons le volume de
control e L3 contenant une bulle (Chahed et al 2004). Le taux de vide est donne par :
pd3
a = -L- (5-1)
1
1
P~| 3
L = hd a 3 avec h= (5-2)
6
Cette estimation place, d’emblee ce developpement dans le cadre des ecoulements a bulles tres
dilues.
180
CH4E/TEE J / CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECO^EEME^TE 4 E^EEEE
!r
La figure (5-1) montre le domaine de reference sur lequel va porter l’analyse. Ce domaine est
cubique, le sillage cree par la bulle est un sillage axisymetrique qui met en jeux comme pour le
sillage plan (voir paragraphe 4.5.2 au chapitre 4) :
- Deux echelles de longueur caracteristiques : l'echelle diffusive l tres faible devant l'echelle
advective L.
Us — Uo u axe du sillage .
du 1 d(yv)
-+- —0 (5-4)
dx y dy
U du — 1 d(y u 'v')
(5-5)
0 dx y dy
1 dp 1 d(y v'2)
— 0 (5-6)
r L dy y dy
181
C%4f/T%E J / CO^C^ES C/SX/EEEES MACES 2-D EA ECO^EEMEATS ^ E^EEES
On remarque que pour les ecoulements a tres faibles taux de vide on identifie la vitesse relative a
celle d'une bulle evoluant en milieu infini U0 — UR. En revanche quand la bulle evolue en milieu
confine comme c'est le cas en ecoulement a bulles a taux de vide plus important, on a U0 > UR du
fait de l’entrainement du liquide associe au confinement. Comme dans le sillage libre U0 est definie
Le flux de quantite de mouvement a travers la section transversale qui se conserve s'ecrit dans le
cas du sillage axisymetrique sous la forme suivante :
2 1
M0 — pL j[U(U0 -U)]2pydy — -CDAU2 (5-7)
0 2
M0 —p LU2 9 2 (5-8)
9 —, CdA
(5-9)
2
D'une maniere analogue au cas du sillage plan, la recherche d'une solution affine pour l'equation
(5-5) de la forme :
—exp(-V (5-11)
182
CK4E/TEE J / COCCAES C/SX/EEEES MACES 2-D EA ECOCEEMEATS ^ ECEEES
Us(x)
2
xN
— 2.8( 9) 3 (5-12)
Un
M — 0.6(—)3
(5-13)
9 9
2
u'v' 1 „ z X
2 GXp( ^^) (5-14)
US 14.1
2KS)
k — -3.75U (5-15)
14.1
mouvement reste inchangee puisque ce bilan est controle par la poussee d'Archimede, la meme dans
les deux cas. Dans ces conditions la vitesse moyenne de glissement peut etre obtenue a partir de la
conservation du debit liquide a travers la section du volume de controle. Pour des bulles
ellipsoi'dales, la conservation du debit liquide impose :
183
C%4f/TEE J COCC^ES C/SX/EEEES MACES 2-D EA ECOCEEMEATS ^ ECEEES
perpendiculairement au mouvement relatif des bulles). En notant d le diametre equivalent des bulles
1
on a d — g3dp ou g est le coefficient d'aplatissement defini par le rapport de la longueur du petit
4g 3
p 3
En remarquant que L — hd a 3 avec h — l'equation (5-18) s’ecrit
6
2 2
(U0 - UR ) p p
2
a3 Soit UR — U0(1 - 2
a3) (5-19)
3^2
4g3 h2 4g3 h
2
La puissance en de a signifie que l’intensite du reflux est inversement proportionnelle au taux
3
de presence surfacique du liquide. Tout se passe comme si la bulle ajuste sa vitesse relative pour
compenser l’augmentation du cisaillement due a l’entrainement du liquide du fait du confinement
lateral du sillage. Rappelons que le debit entraine dans le sillage est independant de x. Garnier et al
(2002) proposent une correlation en puissance 1/3 de a supposant que la modification de la vitesse
relative est due a des interactions hydrodynamiques en relation avec la distance inter-bulles.
importants, Larue de Tournemine (2001). On adopte cette demarche pour commenter le resultat
(5-18). La conservation de quantite de mouvement impose :
2
CpA 4 3\2 CA
C^U2
UL(1 2
a3) (5-20)
2 2
pg3 h2
184
C%4f/TEE J COCC^ES C/SX/EEEES MACES 2-D EA ECOCEEMEATS ^ ECEEES
ou CD0 est le coefficient de trainee sur une bulle isolee. Le coefficient de trainee s’exprime alors
par :
C
CD ^D 0 (5-21)
4 2 0
(1-------- — a 3)2
3 2
pg h
Les experiences de Garnier et al (2002) ont ete menees dans une colonne cylindrique de
diametre 60 mm. Les bulles sont injectees uniformement dans la section. Les auteurs ont procede
de telle sorte qu'ils conservent l'homogeneite de l'ecoulement dans le noyau central de la colonne.
La gamme des valeurs de taux de vide etudiee s'etale de 0.01 a 0.4 avec des bulles de diametres
equivalents compris entre 1.5 mm et 3 mm. D'autre part les auteurs ont observe que les bulles ont
une forme ellipsoidale avec un coefficient d'aplatissement estime de l'ordre de 0.5.
2
UR — U0(1 - 1.92a3) (5-22)
et
C
^D0
CD 2 (5-23)
(1 - 1.92a 3)2
La relation (5-22) indique que la vitesse UR admet la vitesse U0 comme valeur asymptotique
lorsque le taux de vide tend vers zero. Nous comparons dans les figures (5-3) et (5-4) le rapport
des vitesses calcule a partir de la relation (5-22) aux donnees experimentales de Garnier et al
2
Les figures (5-2) et (5-3) presentent la variation du rapport en fonction de a et a3 (Chahed
Un
et al 2004). Ces figures montrent que l’expression de la vitesse relative reproduit bien le glissement
moyen mesure par Garnier et al (2002). La figure (5-3) a l’air d’indiquer que meme si les
185
CfMP/7%E J C7&47EEEEEviaVCEE 2-D EVECOC/IEta'ATE^ E&/EEEE
interactions hydrodynamiques, a forts taux de vides sont importantes la loi du glissement en 2/3 du
taux de vide est maintenue jusqu'a des taux de vide atteignant 15%. Au dela de cette valeur le
confinement devient tel que les interactions hydrodynamiques deviennent importantes et le schema
de calcul se trouve insuffisant, Lance (1986). La loi de variation en puissance 1/3 de a proposee
par Gamier et al (2002) semble mieux approximer les resultats experimentaux pour les forts taux de
vide. On peut interpreter ce resultat par l'effet important des interactions hydrodynamiques dans ce
cas qui sont en relation avec la distance inter-bulles.
JL-0 0U
vitesse relative
Figure (5-2): Variation de la vitesse relative en ecoulement a bulles avec le taux de vide.
Comparaison aux resultats experimentaux de Gamier et al (2002)
♦ JL=0.044
Vitesse relative
■ JL=0.033
JL=0.016
JL=0.062
-------Calcul
alpha2/3
186
CH4E/TEE J / CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECO^EEME^TE 4 E^EEEE
L'equation (5-14) indique que le profil adimensionnel (sillage lointain) du frottement turbulent a
uV
une valeur maximale pour X = 1 ou on a —j— = 0.043 et il s'annule sur l'axe et aux frontieres du
sillage. D'autre part, on sait par l'experience que dans les couches cisaillees on a :
on a alors en premiere estimation pour l'ecoulement externe lorsque X > 1 la relation suivante :
1
k = -3.75U2u 'v'»- uv (5-25)
Le resultat (5-24) est bien connu pour les ecoulements proches d'un equilibre production-
dissipation. On doit donc avoir en premiere estimation :
2 f '(X)
k = -3.75U S (5-26)
R ET
X2
avec f (X) = exp( ■) et Ret 12.5 (5-27)
2
On remarque que cette solution n'est pas valable dans la zone centrale du sillage (X < 1). Les
experiences montrent que l'energie turbulente n'est pas nulle sur l'axe du sillage alors que u'v' = 0
et f' (X) = 0 a X = 0. On pourra etablir une estimation acceptable de k dans la zone X < 1 en
187
CfMP/7%E J Ci&LEEEEEvtgVCEE 2-D EVECO[/IEt7EATE4 E&/EEEE
Dans la zone d'affinite du sillage plan, 1 < c, < V3 , Tennekes et Lumly (1973) montrent que la
On retrouve des resultats similaires pour le sillage axisymetrique (figure (5-4)) ou le modele de
l'equilibre production-dissipation donne par l'expression (5-26) produit un profil (trait continue) qui
sous-estime l'energie turbulente a l'axe du sillage ou l'apport par convection montre que le niveau
est en realite plus important en comparaison avec les mesures (trait discontinu) de Uberoi et
Freymuth (1970). A l'exterieur du sillage l'expression (5-26) donne des niveaux d'energie plus
importants et le modele de l'equilibre production-dissipation est insuffisant dans cette zone ou le
bilan de l'energie cinetique est plutot domine par l'equilibre advection diffusion.
II n'en reste pas moins que le modele de l'equilibre production-dissipation semble bien representer le
bilan global de l'energie cinetique turbulente a travers toute la section du sillage.
----- k/kmax
188
C%4f/TRE J / CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECO^EEME^TE ^ E^EEEE
Autrement dit, si on s'interesse au montant global de l'energie turbulente dans le sillage on peut se
contenter d'integrer l'expression (5-17) pour representer le montant global.
Considerons que l’origine fictive du sillage est ax = 0. On considere que la largeur de reference a
Pour les ecoulements dilues (L>>d) on considere que l’abscisse finale xF est tel que
x0 1 x
(xF - x0) » L. On pose — =m=
0-------
3_2 ^
en notant X = — on a :
d 8b md
et XF = =hal+1 (5-27)
X=md=1
md m
On obtient :
2
l(x) = dX3 et Us(x) = 4 %%X (5-28)
2
XF 2 1
ks = V2P%2md33R754x,^T^V
2U0 jX-3
3 dX = 22md
mj33 3^4%
3-75 2U0(XF -1) (5-29)
^ET Rr
En remarquant que le volume du liquide dans le volume de controle est VL = q3d3a '(1 -a),
3 75 (— +1)3 -1 (+1)3 -1
= W2P 4c4m3^75—m
ks = 4c4md3RET"U2 hJdm-a-'(1 -a) a -u2 (5-30)
R, h (1 -a)
2
P
avec Ur = U,(1 - 2 a3)
h2
4g,33 h
189
CH4E/TEE J / CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECO^EEME^TE 4 E^EEEE
L’energie cinetique induite par les bulles est ecrite comme la somme de l'energie produite dans le
sillage et de la pseudo-turbulence gouvernee par les effets de masse ajoutee :
ou kSP est l’energie cinetique produite par la masse ajoutee dans le mouvement relatif moyen et qui
On compare dans la suite les expressions de l’energie cinetique induite par les bulles avec les
resultats experimentaux de Garnier (2001). A cet effet on pose :
,2 2 == 4
uS — kSP et u Sp =— kSP conformement a la solution potentielle.
La figure (5-5) montre une comparaison des resultats experimentaux avec la solution exprimant
l’energie cinetique totale induite par les bulles en fonction de U0. La figure (5-6) presente la meme
comparaison mais l’energie cinetique induite par les bulles est exprimee en fonction de UR
□ JL=0.044m/s
Bubble Induced Turbulence
A JL=0.033m/s
X JL=0.016m/s
JL=0.062m/s
— ksp
0,08
- - ks
kBI
0,06
0,04
0,02
0,1 0,2
alpha
Figure (5-5) : Energie totale induite par les bulles et energie deduite de la solution potentielle
normees par U0. Comparaison avec les resultats experimentaux de Garnier et al (2002).
190
CH4E/7EE J 2-D EVECOr/IEta'ATE,! E&/Z2EE
□ JL=0.044m/s
Bubble induced turbulence
A JL=0.033m/s
X JL=0.016m/s
+ JL=0.062m/s
------- ksp
■ ■ ■ ■ ks
P 0,15
X A
Figure (5-6) : Energie totals induite par les bulles et energie deduite de la solution potentielle
normees par UR . Comparaison avec les resultats experimentaux de Gamier et al (2002).
La figure (5-5) indique que Venergie cinetique normee par U0 est relativement bien reproduce
C’est un resultat qui est en accord avec les experiences d’Amaury (2001) qui indiquent que c’est en
effet la vitesse du champ lointain qui semble bien normer Fenergie cinetique turbulente induite par
les bulles.
ou ksp est F energie cinetique produite par la masse ajoutee dans le mouvement relatif moyen et qui
191
C%4f/TRE J / CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECO^EEME^TE ^ E^EEEE
et ks est l’energie cinetique produite dans le sillage qui en ecoulement homogene s’ecrit :
qa 3 1
+1)3i 3 -1 p
3.75 m a
ks = V2P 4c4m -U2 avec UR = U0 (1 - - 2 a3) (5-35)
h (1 -a) 3^2
4g3 h
En milieu inhomogene, on ne dispose pas de resultats sur la generation du champ fluctuant par le
mouvement relatif des bulles et nous retenons l’ equation suivante de transport de la
turbulente induite par les bulles:
L’equation (5-36) ainsi modelisee, exprime que l’ecart de la turbulente induite par les bulles par
rapport a la solution en ecoulement homogene represente une diffusion due precisement au
caractere inhomogene de l’ecoulement.
ainsi que des niveaux satisfaisant de l'energie turbulente pour les differents taux de vide.
192
CfMP/7%E J 2-D EVECOr/IEta'ATE^ E&/Z2EE
0,1
—. monophasique * data (alpha = 0%)
■ ■ ■ diphasique (alpha=1 %) * data (alpha = 1 %)
0,08
1 diphasique (alpha=0,5%) ° data (alpha = 0,5%)
diphasique (alpha=2%) ■> data (alpha = 2%)
0,06
=i
0,02 -
20 WM 40 60
5.5 Conclusions
Les resultats experimentaux obtenus en ecoulements a bulles a forts taux de vide montrent un
resultat essentiel qui consiste a la decroissance de la vitesse relative des bulles lorsque le taux de
vide augmente. Par consequences, le modele de l'energie pseudo-turbulente base sur la solution
potentielle qui attribue la production de la pseudo-turbulence a l'effet de la force de masse ajoutee
s'avere incomplet.
Nous avons propose dans ce chapitre un modele de turbulence base sur une modelisation de la
vitesse relative des bulles en function du taux de vide et une decomposition de l'energie turbulente
en milieu dense en une contribution kSI produite par cisaillement dans le liquide et une contribution
kBI induite par les bulles comprenant l'energie produite au niveau des sillages des bulles ainsi que
La mise en oeuvre du modele "de la turbulence induite par les bulles" pour la simulation des
experiences de Gamier et al (2002) montre une reproduction satisfaisante des niveaux de l'energie
193
CH4E/TEE J / CO^C^EE C/EX/EEEEEMWCEE 2-D E^ ECO^EEME^TE 4 E^EEEE
turbulente. Les resultats numeriques obtenus montrent egalement que la contribution essentielle a
cette energie correspond a l'energie produite dans les sillages des bulles.
194
Conclusion
Le travail developpe dans cette these porte sur le developpement de modeles euleriens a deux
fluides au premier ordre deduits par reduction des fermetures au second ordre, et ceci dans
l'objectif d'elaborer des modeles a deux fluides adaptes aux applications industrielles qui preservent
le mieux le contenu physique des fermetures au second ordre.
second ordre. Cette demarche nous a permis de deduire deux modeles (k0 -kS - e et k0 -kS - w)
dont la mise en reuvre en ecoulement homogene a prouve qu'ils realisent un niveau de description
de la structure de la turbulence aussi satisfaisant que celui obtenu par la modelisation au second
ordre. Le modele k0 -kS -e a ete generalise dans le chapitre 4 au cas des ecoulements
inhomogenes. Son implantation dans le code numerique industriel "MELODIF" nous a permis de
l'appliquer dans le cas de l'ecoulement a bulles cisaille libre de sillage. L'analyse des resultats
numeriques du point de vue du comportement affine nous a permis d'interpreter la modification des
echelles caracteristiques de l'ecoulement a bulles par comparaison au cas monophasique
correspondant. Nous avons considere dans le chapitre 5 le cas des milieux denses caracterises par
des forts taux de vide, le modele developpe dans ce cas tient compte de la modification de la vitesse
relative lorsque le taux de vide augmente ainsi que de l'energie turbulente produite dans les sillages
des bulles.
Les developpements de modelisations plus generales porteront sur la prise en compte de l'effet
des interactions hydrodynamiques dans les lois de fermeture du transfert interfacial de quantite de
195
mouvement, sur l'approfondissement de la modelisation de la turbulence et sur la modelisation de la
distribution des phases
Les phenomenes de transfert sont en relation directe avec la mise en contact des phases, donc
avec le taux volumique de presence de phases mais aussi avec le taux surfacique de presence de
l'aire interfaciale. Se precise alors toute l’importance de la maitrise de la distribution des bulles et de
leurs tailles dans le phenomene de transfert dans les systemes diphasiques gaz-liquide. Il s'agit, du
point de vue de la modelisation, de progresser dans l’ amelioration du contenu physique des
modeles pour mieux representer les differentes interactions en milieu diphasique : interaction
interfaciales, interactions hydrodynamiques en milieu dense, interaction et echange interfacial en
ecoulement de paroi, transport de l'aire interfaciale, coefficient de transfert de masse en ecoulement
turbulent diphasique.
Les experiences recentes en ecoulements gaz-liquide a forts taux de vide montrent d'importantes
modifications de la structure des l’ecoulement diphasique. Dans ces ecoulements, l’agitation induite
par les bulles devient preponderante et la structure de l’ecoulement diphasique est essentiellement
gouvernee par le mouvement relatif des bulles. Il devient donc difficile de bien representer
l’ecoulement diphasique et ses proprietes de transport turbulent sans une description precise des
mecanismes qui controlent la turbulence "classique" du liquide, la "pseudo-turbulence" induite par
les bulles et leurs eventuels couplages.
Des travaux en cours sur les transferts en milieu gaz-liquide turbulent, non rapportes dans ce
memoire de these mais auxquels nous contribuons (Lemaissi et al 2003), montrent que les echelles
caracteristiques des transferts dependent des echelles de la turbulence en general mais dependent
egalement des echelles propres qui caracterisent la bulle et son mouvement. Les echelles qui
semblent determiner les transferts aux interfaces sont celles qui sont inferieures a la taille de la bulle
(les echelles vues par la bulle); les echelles plus grandes ne sont pas sans effets sur le transfert
196
puisqu'elles determinent les mouvements des grandes structures donc l’ hydrodynamique des
systemes gaz-liquide. Le probleme est d’autant plus complexe que les champs moyens et fluctuants
des deux phases sont perturbes par les echanges interfaciaux: c’est donc dans une perspective
globale de modelisation des ecoulements diphasiques qu’il faudrait placer le probleme de l’etude
des transferts dans les systemes gaz-liquide.
D'autre part et compte tenu de l'importance de la distribution des tailles des bulles dans les
phenomenes de transfert en milieu gaz-liquide, il nous parait indispensable de prendre en compte la
distribution des diametres pour mieux preciser la formulation des forces interfaciales et aboutir a
une meilleure representation de la distribution des bulles. Pour y parvenir il va falloir construire des
modeles a N+1 fluides (N dispersions de tailles differentes et la phase continue). Une telle
demarche mecaniste, bien que theoriquement envisageable, va poser d’enormes difficultes de
modelisation notamment en ce qui concerne la modelisation des interactions entre la phase continue
et les differentes classes de la dispersion (interactions avec le champ moyen et fluctuant) et la
description des echanges entre les N dispersions (interactions hydrodynamiques, coalescence,
rupture).
197
References bibliographiques
Ahmed A. M. et Elgobashi S., (2000), "On the mechanisms of modifying the structure of turbulent
homogeneous shear flows with dispersed particles", Physics Fluids, vol. 12, pp. 2906-2930
Antal S.P, Lahey JR, Flaherty J.E, (1991), "Analysis of phase distribution in fully developped
laminar bubbly two-phase flow", Int. J. Multiphase Flow, vol. 5, pp. 635-652
Arnold G. S., Drew D. A., Lahey R. T., (1988), "Derivation of constitutive equation for interfacial
force and Reynolds stress for a suspension of sphere using ensemble averaging", Journal of
Chemical Engineering Communications, vol. 86, pp. 43-54
Auton T. R, Hunt J. C. R, Prud'Homme M., (1988), "The force exerted on a body in inviscid
unsteady non-uniform rotational flow", J. Fluid Mech., vol. 197, pp. 241-257.
Batchelor G. K. et Townsend A. A., (1948), "Decay of turbulence in the final period", Proc. Roy.
Soc., A (London), vol. 194(1039), pp. 527-543
Batchelor G. K., (1972), "Sedimentation in a dilute dispersion of spheres", J. Fluid Mech., vol. 52,
pp. 245-268
Bel Fdhila R., (1991), "Analyse experimental et modelisation d'un ecoulement vertical a bulles
dans un elargissement brusque", These de Doctorat de l’Institut National Polytechnique de
Toulouse
Bel Fdhila R., Simonin O. (1992). "Eulerian prediction of a turbulent bubbly flow downstream of a
sudden pipe expansion". Workshop on Two-phase flow predictions, 30 march-2 April, Erlangen
Bel Fdhila R. et Duineveld P.C., (1996), "The effect of surfactant on the rise of a spherical bubble
at high Reynolds and Peclet numbers", Phys. Fluids, vol. 8, pp. 310-321
198
.REFERENCES E/gE/CGRARN^gCES
Bellakhal G., Chahed J., Masbernat L., (2004) a, "Analysis of the turbulence structure in
homogeneous shear bubbly Flow using a turbulent viscosity model", Journal of Turbulence, vol.
5,2004,036
Bellakhal G., Chahed J., Masbernat L., (2004) b, "k-Omega turbulence model for bubbly flows",
5th International Conference on Multiphase Flow, ICMF’05 Yokohama, Japan
Biesheuvel A., Van Wijngaarden L., (1984), "Two-phase flow equation for a dilute dispersion of
gas bubbles in liquid", J. Fluid Mech., vol. 148, pp. 301-318
Boivin M., Simonin O., Squires K. D., (2000), "on the prediction of gas-solid flows with two-way
coupling using Large Eddy Simulation", Physic Fluids, vol. 12, pp. 2080-2090
Boussinesq J., (1885), "Sur la resistance qu'oppose un fluide indefini au repos, sans pesanteur, au
mouvement varie d'une sphere solide qu'il mouille sur toute la surface, quand les vitesses restent
bien continues et assez faibles pour que leur carres et produits soient negligeables, C.R.A.S. Paris,
100, pp. 935-937
Burton T. et Eaton J. K., (2002), "Analysis of a fractional-step method on overset grid", J. Comp.
Physics, vol. 177, pp. 1-29
Chahed J., (1999), "Forces interfaciale et turbulence dans les ecoulements a bulles", These de
Doctorat de l’Ecole Nationale d'Ingenieurs de Tunis
Chahed J., Colin C., Masbernat L., (2002), "Turbulence and phase distribution in bubbly pipe flow
under micro-gravity condition", J. Fluid Eng., vol. 124(4), pp. 951-956.
Chahed J., Roig V., Masbernat L., (2003), "Eulerian-Eulerian two-fluid model for turbulent gaz-
liquid bubbly flows", Int J. of Multiphase Flow, vol. 29-1, pp. 23-49
Chahed J., Bellakhal G., Masbernat L. (2004), "Turbulence and pseudo-turbulence modelling in
low and high void fraction bubbly flows', 5th International Conference on Multiphase Flow,
Clift R., Grace J. R., Weber M. E., (1978), "Bubbles, Drops and Particles", Academic Press, New
York.
199
Champagne F. H., Harris V. G., Corsin S., (1970), "Experiments on nearly homogeneous turbulent
shear flow", J. Fluid Mech., vol. 41(1), pp. 81-139
Comte-Bello G., Corsin S., (1966), "The use of a contraction to improve the isotropy of grid
generated turbulence", J. Fluid Mech., vol. 25(4), pp. 657-682
Csanady G.T., (1963), "Turbulent diffusion of heavy particles in the atmosphere" J. Atm. Sc vol 20
pp. 201-208
Drew D. A., Wallis G. B., (1964), "Fundamentals of two-phase modelling", Third International
Workshop on Two-phase Flow Fundamentals, Imperial College, London, June 15-19.
Drew D.A., Lahey R.T., (1982), "Phase distribution mechanisms in turbulent low-quality two-
phase flow in circular pipe", J. Fluid Mech., vol. 117, pp. 91-106.
Drew D. A., (1983), "Mathematical modelling of two-phase flow", Ann. Rev. Fluid Mech., vol.
15, pp. 261-291.
Drew D.A., Wood R.T., (1985), "Overview and taxonomy of models and methods for workshop
on two-phase flow fundamentals", National bureau of standards, Gaitherburg, Maryland,
September 22-27
Elgobashi S. et Truesdell G. C., (1993), "On the two-way interaction between homogeneous
turbulence and dispersed solid particles. I. Turbulence modification", Phys. Fluids, vol. 5, pp.
1790-1801
Fadlun E. A., Verzicco R., Orlandi P., Cabot W. H., (2000), "Combined immersed boundary
finite-difference methods for three-dimensional complex flow simulations", J. Comp. Fhysics, vol.
161, pp. 35-60
Fevrier P. et Simonin O., (2001), "On the spatial distribution of heavy particle velocities in
turbulent flow : from continuous field to particulate chao", Turbulence and shear flow phenomena,
second Intnl. Symp, KTH, Stockholm
Gabillet C., Colin C., Fabre J., Larrauri D., Briere E., (1998), "Experimental study in a turbulent
boundary layer", Third Int. Conf. On Multiphase flow, Lyon, June 8-12
200
Garnier C., Lance M., Marie J.L., Mareuge I., (2002), "Experiments on buoyancy-driven bubbly
flows at low and high void fractions". Proceedings 10th Workshop on Two-Phase Flow
Gatignol R., (1983), "The faxen formulae for a rigid particle in an unsteady non-uniform stokes
flow", Jounal de mecanique theorique et appliquee. 9, 143-160.
Grossetete C., (1995) b, "Experimental investigation and preliminary numerical simulations of void
profile development in a vertical cylindrical pipe", 2nd Int. Conference on Multiphase Flow, Kyoto,
Japan.
Hadamar J. S., (1911), Comptes Rendus a l’Academie des Sciences de Paris, C.R.A.S., Paris,
vol. 152, pp 1735-1738.
Harris V. G., Graham J. A., Corsin S., (1977), "Further experiments in nearly homogeneous
turbulent shear flow", J. Fluid Mech., vol. 81(4), pp. 657-687
Hirt C. W. et Nichol B. D., (1981), "Volume of fluid method for the dynamics of free surfaces
boundaries", J. Comp. Phys., vol. 39, pp. 201
Hu H. H., Patankar N. A., Zhu M. Y., (2001), "Direct numerical simulations of fluid-solid systems
using Arbitrary-Lagrangian-Eulerian technique", J. Comp. Physics, vol. 169, pp. 427-462
Ishii M., (1975), "Thermo-fluid dynamic theory of two-phase flow", Eyrolls, Collection de la
Direction des Etudes et Recherche d'Electricite de France.
Issa R. I. et Oliveira P. J., (1996), "Validation of two-fluid model in shear-free mixing layers.
Fluids Engeneering Division Conference, 236
Kamp A. Colin C. Fabre J., (1995), "The local structure of a turbulent bubbly pipe flow under
different gravity conditions" 2nd International Conference on Multiphase Flow, Kyoto, Japan.
Kistler A. L., O'brien V., Corsin S., (1954), "Preliminary measurements of turbulence and
temperature fluctuations behind a heated grid", Reas. Mem. N.A.C.A., 54D19
201
Kolmogorov A. N., (1941), "On the generation of isotropic turbulence in an incompressible
viscous fluid", Dokl Akad. Nauk SSSR, Physics, vol. 31, pp. 538-540
Lahey R.T, Lopez de Bertodano M., (1991), "The prediction of phase distribution using two-fluid
models", ASME/JSME Termal Engineering proceeding, vol. 2, pp. 193-200"
Lance M., (1986), "Etude de la turbulence dans les ecoulements diphasiques disperses", These de
doctorat es Sciences, Universite Claude Bernard, Lyon, France.
Lance M. et Bataille J., (1991), "Turbulence in the liquid phase of a uniform bubbly air water
flow", J. Fluid Mech., vol. 222, pp. 95-118.
Lance M., Marie J.L., Bataille J., (1991), "Homogeneous turbulence in bubbly flows". Journal of
Fluids Engineering, vol. 113, pp. 295-300.
Lance M., Marie J. L., Moursali E., Bataille J., Suzanne C., Roig V., Bel Fdhila R., Masbernat L.,
(1996), "Experimental study of turbulent bubbly shear flows", Chem. Eng. Commun., vol. 141, pp.
51-70
Landau L. D., Lifshitz E. V., (1959), "Course of theoretical physics", Fluid Mechanics, vol. 6,
Launder B.E., Reece G.J., Rodi W., (1975), "Progress in the development of a Reynolds stress
turbulence closure", J. Fluid mech., vol. 68, part 3, pp 537-566.
Larue de Tournemine A., (2001), "Etude experimental de l'effet du taux de vide en ecoulements
diphasiques a bulles". These de Doctorat a l'Institut de Mecanique des Fluides de Toulouse
Lee S.J, Lahey Jr R.T, Jones Jr O.C., (1989), "The prediction of two phase turbulence and phase
distribution phenomena using k-e model", Japanese Journal of Multiphase Flow, vol. 3, pp 335
368
Legendre D., Boree J., Magnaudet J., (1998), "thermal and dynamic aspects of the evolution of a
spherical bubble moving steadily in a superheated or subcooled liquid", Physics of fluids, vol.
10(6), pp. 1256
Ling S. C., Huang T.T., (1970), "Deacay of weak turbulence", Phys. Fluids, vol. 13(12), pp. 2912
2924
202
Liu T.J., Bankoff S.G., (1990), "Structure of air-water bubbly flow in a vertical pipe: I- Liquid
mean velocity and turbulence measurements", Int. J. Heat and Mass Transfer., vol. 36 (4),
pp 1049-1060.
Lopez de Bertodano M., Lee S.J., Lahey R.T., Drew D.A., (1990), "The prediction of two-phase
turbulence and phase distribution using a Reynolds stress model", Journal of Fluids Engineering,
vol. 112, pp. 107-113.
Lopez de Bertonado M., Lee S.J., Lahey R.T., Jones. O. C., (1994), "Development of a k-e
model for bubbly two-phase flow", Journal of Fluids Engineering, vol. 116, pp. 128-134
Magnaudet J., Rivero M., Fabre J. (1995), "Accelerated flows past a rigid sphere or a spherical
bubble. Part 1. Steady staining flow", J. Fluid. Mech., vol.284, pp. 97-135
Mareuge I., Lance M. (1996), "Bubble-induced dispersion of a passive scalar in bubbly flows",
Proceeding of the 2nd International Conference on Multiphase flow, Kyoto, April 3-7
Maxey R., Riley J., (1983), "Equation of motion for a small rigid sphere in non-uniform flow",
Phys. of Fluids, vol. 26 (4), pp. 883-889.
Miller R. S. et Bellan J., (1999), "Direct numerical simulation of a confined three-dimensional gas
mixing layer with one evaporating hydrocarbon-droplet-laden stream", vol. 384, pp.293-338.
Mills R. R., Kistler A. L., O'brien V., Corsin S., (1958), "Turbulence and temperature fluctuations
behind a heated grid", Tech. Note N.A.C.A., 4288
Minier J. P., (1992), "Construction de la trajectoire d'une particule dans un ecoulement turbulent",
Electricite de France, Direction des Etudes et de Recherches.
Moursali E., Marie J.L., Bataille J., (1995), "An upward turbulent bubbly layer along a vertical flat
plate", Int. J. Multiphase Flow, vol. 21, pp. 107-117.
203
Mulhearn P. J., Luxton R. E., (1975), "The development of turbulence structure in a uniform shear
flow, homogeneous turbulent shear flow", J. Fluid Mech., vol. 68(3), pp. 577-590
Odar F., Hamilton W. S., (1964), "Forces on a sphere accelerating in a viscous fluid", J. Fluid
Mech., vol. 18 Part 2, pp. 229-240.
Odar F., (1966), "Verification of the proposed equation for the calculation of the forces on a
sphere accelerating in a viscous fluid", J. Fluid Mech., vol. 25, pp591-592.
Reeks M. W., (1991), "On a kinetic equation for the transport of particles in turbulent flows",
Physics of Fluids, vol. 3, pp. 446-456
Rivero M., (1991), "Etude par simulation numerique des forces exercees sur une inclusion
spherique par un mouvement accelere", These de Doctorat de l’Institut National Polytechnique de
Toulouse
Roig V., (1993), "Zones de melange d'ecoulements diphasiques a bulles", These de Doctorat de
l’Institut National Polytechnique de Toulouse.
Roig V., Suzanne C., Masbernat L., (1998), "Experimental investigation of a turbulent bubbly
mixing layer", Int. J. Multiphase Flow, vol. 24, pp. 35-54
Rose W. G., (1966), "Results of an attempt to generate a homogeneous turbulent shear flow", J.
Fluid Mech., vol. 25(1), pp. 97-120
Rose W. G., (1970), "Interaction of grid turbulence with a uniform mean shear", J. Fluid Mech.,
vol. 44(4), pp. 767-779
Sangani A. S.,Didwania A. K., (1993), " Dynamic simulations of flows of bubbly liquids at large
Reynolds numbers", J. Fluid Mech., vol. 250, pp. 307-337
Sato Y., Sadatomi L., Sekoguchi K. (1981), "Momentum and heat transfer in two phase bubbly
flow", Int. J. Multiphase Flow, vol. 7, pp. 167-190
Sato Y., Sekoguchi K. (1975), "liquid velocity distribution in two-phase bubble flow", Int. J.
Multiphase Flow, vol. 2, pp. 79-95
204
Schoneborn P. R., (1975), "The interaction between a single partical and an oscillating flow, Int. J.
on Multiphase flow, vol. 2, pp307-317.
Serizawa A., Kataoka I., Michiyoshi I., (1975), "Turbulence structure of air-water bubbly flow : I-
Measuring techniques", Int. J. Multiphase Flow, vol. 2, pp. 221-233
Serizawa A., Kataoka I., Michiyoshi I., (1975), "Turbulence structure of air-water bubbly flow :
II- Local properties", Int. J. Multiphase Flow, vol. 2, pp. 235-246
Serizawa A., Kataoka I., Michiyoshi I., (1975), "Turbulence structure of air-water bubbly flow :
III- Transport properties", Int. J. Multiphase Flow, vol. 2, pp. 247-259
Serizawa A. et Kataoka I., (1990), "Turbulence suppression in bubbly two-phase flow", Nuc. Eng.
Des., vol. 122, pp. 1-16
Serizawa A., Kataoka I., Michiyoshi I., (1986), "phase distribution in bubbly flow. Data set N°24",
Proceeding of the Second International Workshop on Two-phase Flow Fundamentals.
Serizawa A., Kataoka I., Michiyoshi I., (1992), "Phase distribution in bubbly flow", Multiphase
Science and technology, vol. 6, Hewitt, G.F., Delhay J.M., Zuber N., Eds, Hemisphere Publ. Corp,
pp. 257-301
Simonin O., Viollet P.L., (1989), "Numerical study on phase dispersion mechanism in turbulent
bubbly flows" Int. Conf. on Mechanics of two-phase Flows, June 12-15, Taipei, Taiwan.
Simonin O., Viollet P.L., (1990), "Prediction of an oxygen droplets pulverization in a compressible
subsonic coflowing hydrogen flow", Symposium on Numerical Methods for Multiphase Flows,
June 4-7, Toronto, Canada
Simonin O., (1990), "Eulerian formulation for particle dispersion in turbulent two-phase flows",
Proc 5 th Workshop on Two-Phase Flow Predictions, March 19-22, 1990, Erlangen, FRG
Simonin O. (1991), "Prediction of the dispersed phase turbulence in particle-laden jets", 4th
Squires K. D. et Eaton J. K., (1990), "Particle response and turbulence modification in isotropic
turbulence by suspensed particles", Phys Fluids, vol. 2, pp. 1191
Stockes G. C., 1851, "On the effect of the internal friction of fluids on the motion of pendulums,
Trans. Cambridge Phil. Soc., vol. 9, pp. 8-27.
205
Stuhmiller (1977), "The influence of interfacial pressure forces on the character of two-phase flow
model equations", Int. J. Multiphase Flow, vol. 3, pp. 551-560
Sun T.Y., Faeth G.M., (1986), "Structure of turbulent bubbly jets-I. Methods and centerline
properties", Int. J. Multiphase Flow, vol. 12, N°1, pp. 99-114
Sun T.Y., Faeth G.M., (1986), "Structure of turbulent bubbly jets-II. Methods and centerline
properties", Int. J. Multiphase Flow, vol. 12, N°1, pp. 115-126
Sundaram S., et Collins L. R., (1999), "A numerical study of the modulation of isotropic
turbulence by suspensed particles", J. Fluid Mech., vol. 379, Part 2, pp. 105-143.
Sussman M., Smerka P., Other S., (1994), "A level set approach for computing solution to
incompressible two-phase flow", J. Comp. Phys., vol. 114, pp. 146-159
Thai Van D., Minier J.P., Simonin O., Freydier P., Olive J., (1994), "multidimensional two-fluid
model computation of turbulent dispersed two-phase flows", Numerical Method of Multiphase
Flows, ASME, vol. 185
Temkin S., Mehta H. K., (1982), "Drpolet drag in an accelerating and decelerating flow", J. Fluid
Mech., vol. 116, pp. 297-313.
Tennekes H., Lumley J.L., (1973), "A first course in turbulence", MIT press.
Tomiyama A., Zun I., Sou A., Sakaguchi T., (1993), "Numerical analysis of bubble motion with
the VOF method", Nucl. Eng. Des., vol. 141, pp. 69-82
Townsend A. A., (1956), "The structure of turbulent shear flow", University press, Cambridge
Troshko A. A. et Hassen Y. A., (2001), "A two-equation turbulence model of turbulent bubbly
flows", Int. J. Mult. Flow, vol. 27, pp. 1965-2000
Unverdi S. O. et Trygvasson G., (1992), "A front tracking method for viscous incompressible
multi-fluid flows", J. Comp. Phys., vol. 100, pp. 25-37
Van Wijingaarden L., (1976), "Hydrodynamic interaction between gaz bubbles un liquid", J. Fluid
Mech., vol. 77(1), pp. 27-44.
206
Voinov O. V., (1973), "Force acting on a sphere in an homogenous flow of an ideal fluid", J. of
Applied Mechanics and Technical Physics, vol. 14
Wang S.K, Lahey Jr R.T, Jones Jr O.C., (1987), "Three dimensional turbulence structure and
phase distribution measurements in bubbly two phase flows", Int. J. Multiphase Flow., vol. 13,
pp. 327-343.
Wang D. M., Issa R. I., Gosman A. D., (1994),"Numerical prediction of dispersed bubbly flow in a
sudden enlargement", Num. Methods in Multiph. Flow, vol. 185
Wilcox, D. C., (1988), "Reassessment of the scale-determination equation for advanced turbulence
models", AIAA Jl., vol. 26 (11), pp. 1299-1310.
Zuber N., Ishii M., (1979), "Drag coefficient and relative velocity in bubbly, droplet or particulate
flows, AIChE Journal", vol. 25, pp. 843-855.
207