Flux biphasé

Abdoul Basit Ahmad

1221ME02
Flux biphasé

IIT Patna
 Flux biphasé
• Cas particulier de l'écoulement multiphasique
• Un écoulement gaz-liquide est appelé
écoulement biphasique
• Systèmes électriques à grande échelle – cas le
plus couramment étudié
• Également en réfrigération, pompes à chaleur
 Flux biphasé
et climatisation
• La tension superficielle rend les problèmes
dynamiques non linéaires
• Des changements à grande échelle dans la
densité des supports entraînent des effets de
compressibilité
 Flux biphasé
• Le changement de phase n'est pas instantané
et le système liquide-vapeur peut donc ne pas
être en équilibre.
 Quantités importantes
• Vol. un V. fraction de gaz = fraction de vide = a _ A i a A
, e, Wq, e, ,
• Qualité Mmassique du flux X =- -— , W = débit massique = Aou
- Wg+Wf' •
• Qualité volumétrique , B = ——, Q = vol. débit ■ ' ' ' Qg+Qf"

Vitesse superficielle de la phase g, j g = A

Rapport de glissement, rapport des vitesses de la
phase gazeuse et liquide : ug _ w g P f A f _ r X X (P/) (
u
f w fPg A g \1 - x) \pg/ '
 Analyse du flux biphasé
• Trois méthodes d'analyse
• Modèle de flux homogène :
– Considérez que le flux est monophasé.
– Propriétés moyennes pondérées des phases.
• Modèle à flux séparé :
– Phases artificiellement séparées.
– Équations écrites par phase
– Besoin d'informations sur l'interaction entre les phases
– De phase et de mur
– Superficie occupée par chaque phase
• Modèle de modèle de flux :
– Le flux biphasé dans un canal prend certaines configurations spécifiques.
 (suite)
• Modèles de modèles de flux (suite)
– Flow prend l'une des 3 ou 4 configurations prescrites.
– Besoin de comprendre la configuration du flux.
– Les modèles de flux dépendent de divers facteurs, à savoir. géométrie, vitesse
d'écoulement, composition du flux, etc.
– Besoin d'identifier un modèle pour analyser le flux
– Dans les canaux transparents et à faible vitesse, on peut simplement voir et observer.
– Photographie au flash ou radiographie X sinon.
 Modèles de flux courants
• Flux pétillant
– Phase gazeuse distribuée sous forme de bulles discrètes dans une phase liquide
continue
– Est répandu à de faibles concentrations de vapeur
• Flux de limaces
– Les bulles de gaz ont approximativement le diamètre du tuyau
– Le nez de la bulle est typiquement sphérique
– Bulle et paroi séparées par un film de liquide
– Des bouchons de liquide séparent les bulles de gaz adjacentes.
• Flux de désabonnement
– À des vitesses plus élevées
– Formé par la rupture de grosses bulles de vapeur dans le flux de bouchons.
– Le flux de gaz est très chaotique
– Le flux a un caractère oscillatoire
 (suite)
• Flux annulaire vaporeux
– Film liquide épais au niveau des parois et entraînement liquide important dans la phase
gazeuse centrale
– Le film liquide est aéré par de petites bulles de gaz
– Le liquide entraîné apparaît sous forme de grosses gouttelettes, agglomérées en longs
filaments ou « mèches ».

• Flux annulaire
– Film liquide au mur
– Noyau de gaz avec une certaine quantité d'entraînement de gouttelettes
– Ces gouttelettes ne sont pas agglomérées mais finement dispersées
– Ondes cohérentes de grande amplitude
Fig 1 : Régimes en écoulement
multiphasique
 Lacunes
• Aucune méthode générale pour décider du
modèle de flux
– Le modèle de flux est subjectif
– Dépend de nombreux facteurs locaux
– Des paramètres moins bien définis affectent
également le modèle d'écoulement
• Départ de l'équilibre hydrodynamique local
• Présence de traces de contaminants
• par exemple injection d'air dans une paroi
poreuse
– Pétillant près du mur. Éloignez-vous-en
Modèles d'écoulement dans un tube
chauffé vertical
• Cas particulier important
• La présence d'un flux de chaleur modifie le modèle d'écoulement
• Raisons principales:
– Départ de l'équilibre thermodynamique
– Profil de température radial
– Départ de l’équilibre hydrodynamique local tout au long
• Couche limite thermique formée au niveau du mur
• La température du mur dépasse la température de saturation
• La nucléation a lieu aux sites préférés le long du mur
• Les bulles de vapeur se développent et forment un flux pétillant
• La population de bulles augmente le long du tube
 (suite)
• Les bulles finissent par fusionner, entraînant un
écoulement de limaces
• De plus, en raison des vitesses élevées, un
entraînement de liquide observé
• Également évaporation à l'interface liquide-vapeur
• La perte de liquide entraîne un écoulement annulaire
• Finalement, tout le liquide sèche, conduisant à un
écoulement de vapeur monophasé.
• De petites gouttelettes de liquide sont intercalées
 (suite)
• Ceux-ci s'évaporent lentement avec le temps
Fig 2 : Modèles d'écoulement dans un tube
chauffé vertical
 Cartes de modèles de flux
• Nécessité largement ressentie d'une
méthode simple pour prédire les modèles de
flux et les transitions entre eux
• Une façon d’y parvenir consiste à utiliser ce
que l’on appelle les « cartes de modèles de
flux ».
• À l'aide d'un ensemble donné de paramètres
de flux locaux, cartographiez le modèle de
 flux et les transitions entre eux

Développement de cartes
• Basé sur la phase superficielle réelle
vitesses, jg
,jf
• Variables secondaires impossibles à
représenter
• Il est peu probable que ces paramètres aient
le même
 influence dans tous les cas
• Examen de chaque transition
individuellement
en cours
• Attention : la carte n'est qu'un guide
approximatif
103

10
"

102

10

1
1

10

10
10*

+,
j'ai 10 ans 102 10’ 10* 10*
'31
r-Ir 101 je

Ir
1 Ib/s'n 10 102 105 10*

Je "1

10*

Fig 3 : Carte de configuration

d'écoulement pour un écoulement
vertical
• Obtenu par Hewitt et Roberts 1969
• À partir d'observations sur l'écoulement air-eau à basse
pression et vapeur-eau à haute pression dans des tubes
verticaux de petit diamètre (1-3 cm)
• Les axes représentent les flux de quantité de mouvement
superficiels de
phases liquide (PfGf)2) et vapeur (pg( jg ) 2 )
• Peut également exprimer les flux de quantité de mouvement
superficiels en termes de vitesse de masse et de qualité de
vapeur comme suit :
[G(1- X
Pf(f)2= )] 2
pf
 (j 2 [G(1- X
pg g) )] 2
= P. 9

Équations de base de l'écoulement

multiphasique
• Il existe de nombreuses formes sous lesquelles les équations ont été
écrites
• Peut soit être considéré en termes d'instantané local, soit en termes
d'espace ou de temps.
quantités moyennes
• La moyenne simplifie les équations et les rend plus faciles à gérer, au prix
d'une perte d'informations
• Considérant un écoulement en régime permanent unidimensionnel
• Considérons un système comme le montre la figure 4
• Ecoulement multiphasique stratifié dans un canal incliné
• Stratifié pour dériver des équations générales avec des phases en contact
les unes avec les autres et avec les murs
• Supposons des valeurs moyennes de vitesse et de densité de chaque
phase sur n'importe quel plan normal à l'écoulement
• La somme des surfaces occupées par toutes les phases dans n'importe
quel plan normal à l'écoulement est égale à la surface totale du canal
Fig 4 : Modèle simplifié pour un écoulement multiphasique dans un
élément du canal
 Conservation de la masse
• En l'absence d'ajout ou de retrait de masse
dans le système :
8,. un,. X _
— (AkPk) + — (AkPkug) = Tk
• T représente la génération de masse par
unité de longueur
2n=0
 Conservation de l'élan
• Bilan des forces agissant sur un élément
différentiel
• Équilibrer le taux d'afflux et le taux de
création d'impulsion par rapport à la somme
des forces et à la génération d'impulsion due
au transfert de masse
• Forces de pression, forces gravitationnelles,
cisaillement de paroi, forces de cisaillement
interfaciales et génération d'impulsion à partir
du transfert de masse
 Expression
• La conservation de la quantité de mouvement pour

—Un
dp
P gwd
b r
gw
T

un élément différentielz s'écrit :

g sin 9 [A f Pf + A g p g ]
= d ( WfUf + WgUg )
• Où, W k = Aa k p k u k pour la phase k c'est-à-dire le
 débit massique
Conservation d'énergie
• Bilan des différents gains et pertes d'énergie pour l'élément
de commande
• Taux d'augmentation de l'énergie totale (cinétique + interne)
de la phase k + énergie totale convectée dans l'élément de
contrôle
• Equilibré avec la somme de l'apport de chaleur à la phase k et
du travail effectué sur la phase k
• Avec la vitesse à laquelle l'énergie est transférée à travers
l'interface de l'élément de commande
 Expression
La conservation de l’énergie pour un élément différentiel s’écrit :
£ [wgig + W f je f \ + 4 “a(e2" + W/y2 + (w, + Wf)g sin e
= Qw
Où,
^wf^wf + PwgPwg
^ wk est le flux de chaleur du mur à la phase k/ longueur du
mur ik = uk +~ , l'enthalpie par unité de masse
Pk
 Évaluation de la chute de
pression
• Équation de base :
{dFg + dF^ = -T g W P g W dZ — T fw P fw dZ = -A (dF) dZ
• Chute de pression totale :
'^Pr} f d P (dp\
azF)+(aza)+(azz)

où:
2 jours (1-
-(2a)=l2(W,u,+W,u,)=G dzx2v
X)2Vf
Un

F
Un Afi r . 1
+
A P g A Pf] = g sin e V^Pg + (1- ^pf\
Modélisation d'un écoulement
diphasique
 Modèle de flux homogène
– Hypothèses du modèle de flux homogène
– Vitesses égales du liquide et de la vapeur
– Atteinte de l'équilibre thermodynamique entre
les phases
– Utilisation d'un facteur de frottement monophasé
convenablement bien défini pour un écoulement
biphasique, basé sur les propriétés des deux
phases
 Équations régissant
• Pour un écoulement homogène unidimensionnel
constant

• Continuité:
W = Apu
• Élan:
— Ad — dF — Apg sin 3 dz = W du
• Énergie:
3q — ^w = di + d («) + g sin 3dz
 (suite)
• La barre supérieure indique les quantités moyennes
• Aussi:
di = Sq + dE + vdp
dF est la force de cisaillement totale nette du mur
7=v le volume spécifique du fluide homogène
7=w= XVg+(1- X)Vrl = V, + XV f g] = a = 1
p est la densité du fluide homogène dF = r w Pdz
• Maintenant,
pu2\
2)
f TP est un facteur de frottement biphasé pratique qui est évalué à l'aide
d'une viscosité moyenne biphasée, évaluée empiriquement.
 Modèle à flux séparé
• Hypothèses du modèle à flux séparés
– Les vitesses des phases ne doivent pas
nécessairement être égales
– Atteinte de l’équilibre thermodynamique entre
phases
– Utilisation de corrélations empiriques ou de
concepts simplifiés pour relier le multiplicateur de
frottement biphasé (ϕ 2 et la fraction de vide ( a )
aux variables indépendantes de l'écoulement
• Trouvé comme étant le plus valable pour
l'écoulement annulaire
 Équation d'impulsion
• Les équations de base sont déjà données pour l'écoulement séparé
• Le gradient de pression de friction peut être exprimé en termes de
gradient de pression monophasé du flux total considéré comme un
liquide.
dp., . , ■
— Fl est la partie du gradient de
pression statique global nécessaire pour vaincre le frottement
• Nous avons:

—
2ffoG
Vf
• (Pfo) est le multiplicateur biphasé
 (suite)
• En termes de débit de liquide uniquement
2ffG2(1-x)2Vf
D (Pf)
• Donc relation entre ff et ffo

ffoJe- 1 X.
• De même pour (Pfo) : et (Pr)
(Wro)2 = (/)2n -x]175
• Une relation similaire pour ^Pf 0 ) et Qp g )
 Évaluation de (
^fo) et un
• Afin de calculer le gradient de pression net, nous devons évaluer les deux
quantités ci-dessus
• Ces grandeurs doivent être connues en termes de grandeurs de débit
indépendantes
• Ceci est réalisé en utilisant des corrélations empiriques
• Ces corrélations vont du simple au incroyablement complexe
• Aucune corrélation valable dans toutes les conditions d'écoulement
• Chaque corrélation a sa propre plage de validité
• Réalisé pour la première fois par Martinelli et al
• Les corrélations importantes sont les corrélations de Lockhart-Martinelli,
Martinelli-Nelson et Thom.
• Ne considérera que Lockhart-Martinelli, d'autres sont plutôt similaires
La corrélation Lockhart-Martinelli
• Attribuer une zone définie du flux à chaque phase
• Supposons que les équations conventionnelles de chute de pression de
frottement puissent être utilisées pour chaque phase
• Ignore l'interaction entre les phases – le rend quelque peu inexact
• Postulats de base :
– Quatre régimes d'écoulement, définis en fonction du comportement de l'écoulement
– Les chutes de pression en phase liquide et gazeuse sont considérées comme égales,
quels que soient les détails du modèle d'écoulement particulier.
– L'accélération et les composantes de la charge statique de la chute de pression sont
considérées comme négligeables.
– Encore une fois, la déclaration ci-dessus implique que les chutes de pression en phase
liquide et gazeuse sont égales.
Fig 5 : Corrélation de Lockhart-Martinelli
jhank VO)