See discussions, stats, and author profiles for this publication at: [Link]

net/publication/29650053

Convection naturelle dans une cavité horizontale occupée par deux fluides
non miscibles : solution analytique et numérique

Article · August 2007

Source: OAI

CITATIONS READS

0 219

4 authors:

Fakhreddine S. Oueslati Rachid Bennacer

National Engineering School of Carthage, TUNISIA Ecole Normale Supérieure de Paris / Univ. Paris Saclay
47 PUBLICATIONS 729 CITATIONS 543 PUBLICATIONS 5,678 CITATIONS

SEE PROFILE SEE PROFILE

Habib Sammouda Abdelfettah Belghith

University of Sousse College of Computer and Information Sciences, King Saud University, Riyadh, KSA
118 PUBLICATIONS 1,007 CITATIONS 335 PUBLICATIONS 4,202 CITATIONS

SEE PROFILE SEE PROFILE

All content following this page was uploaded by Rachid Bennacer on 19 May 2014.

The user has requested enhancement of the downloaded file.

 Manuscrit auteur, publié dans "JITH 2007, Albi : France (2007)"

13èmes Journées Internationales de Thermique

CONVECTION NATURELLE DANS UNE CAVITE

HORIZONTALE OCCUPEE PAR DEUX FLUIDES NON
MISCIBLES: SOLUTION ANALYTIQUE ET NUMERIQUE
F. S. OUESLATI(1), R. BENNACER(2), H. SAMMOUDA(3) and A. BELGHITH(1)
(1) LETTM, Département de physique, Faculté des Sciences de Tunis, Tunisie.
(2) LEEVAM, Rue d’Eragny–Neuville sur Oise 95031 Cergy-Pontoise Cedex – France
(3) Faculté des Sciences de Monastir, Tunisie.
[Link]@[Link]

Résumé: La convection naturelle, dans le cas de deux fluides non miscibles superposés dans une
cavité horizontalement allongée, est étudiée numériquement et analytiquement. Les frontières
horizontales sont soumises à des flux thermiques croisés. La solution analytique, basée sur
l'approximation d’écoulements parallèles, s'avère en bon accord avec la solution numérique.
L'existence d’écoulements naturels et antinaturels est démontrée. Les solutions multiples d'état
hal-00158373, version 1 - 29 Aug 2007

d'équilibre sont possibles dans une gamme du rapport des viscosités.

Mots-clés: Convection naturelle, écoulements parallèles, fluides non miscibles, écoulement
naturel, écoulement antinaturel.

1. INTRODUCTION
La stabilité et la convection naturelle, dans un système composé de deux fluides superposés
non miscibles, sont étudiées numériquement et analytiquement. Nous considérons une cavité
horizontale allongée. Les frontières horizontales sont soumis à des flux thermiques uniformes q’ alors
que celles verticales sont exposés à un flux d’intensité aq’. Dans la formulation du problème, nous
employons les équations de Naviers-Sokes. Des résultats sont donnés pour une gamme des
valeurs de chacun des paramètres régissant. La solution analytique, basée sur l'approximation
d’écoulements parallèles, s'avère en bon accord avec la solution numérique des équations
régissantes. L'existence des écoulements naturels et antinaturels est démontrée. En effet, le
nombre de Rayleigh critique dépend de l'intensité latérale de chauffage. Dans ces conditions, des
solutions analytiques approximatives pour la convection unicellulaire dans la région centrale de
la cavité peuvent être obtenues en utilisant une prétention de flux parallèle. Des résultats sont
présentés pour des nombres de Rayleigh critiques comme l'effet de l’intensité de chauffage
latéral, du rapport des viscosités cinématique et d'autres paramètres sur l’écoulement et le
transfert de chaleur convectifs. Les solutions multiples d'état d'équilibre sont possibles dans une
gamme de la du rapport des viscosités cinématiques, qui dépend fortement de l’intensité de la
viscosité cinématique du fluide supérieur. Le fait qu'un des deux fluides est plus conducteur
affecte non seulement les transferts thermiques produits par l’un ou l'autre mais également le
nombre critique de Rayleigh.

2. POSITIONNEMENT DU PROBLEME
La géométrie du système et le problème physique sont montrés dans fig. 1. L'approximation
de Boussinesq est adoptée: la densité varie linéairement avec la température, comme:
ρ = ρ0 1 − βT (T − T '0 ) (1)

Albi, France du 28 au 30 Août 2007 1

 13èmes Journées Internationales de Thermique

Ψ 2 =∂Ψ 2 ∂z = 0 ∂T2 ∂z = −1 λr
q’
x = 0, A z =0
(z’, v’) Fluid2 6444447444448
64748 h2 H
Ψ1 = Ψ 2 = 0 aq’ (x’, u’) aq’ Ψ 1 = Ψ2 =0 ∂Ψ1 ∂z =∂Ψ 2 ∂z
∂T1 ∂x = −a
h1 L’ Fluid1 T1 = T2 = 0 ∂T1 ∂z == λr ∂T2 ∂z
∂ 2 Ψ1 ∂z 2 =µ r ∂ 2 Ψ 2 ∂z
2
∂T2 ∂x =− a q’
Ψ1 =∂Ψ1 ∂z = 0 ∂T1 ∂z = −1

Figure 1: configuration du problème et system coordonnée

Les équations adimensionnelles qui décrivent le problème, sont:

∂∇ 2 Ψ1 ∂T
+ℑ( Ψ1 ,∇ 2 Ψ1 ) = Pr1 ∇ 4 Ψ1 − Pr1 RaT 1 (2)
∂t ∂x
∂T1
+ℑ( Ψ1 ,T1 ) =∇ 2T1 (3)
hal-00158373, version 1 - 29 Aug 2007

∂t
∂∇ 2 Ψ 2 ∂T
+ℑ( Ψ 2 ,∇ 2 Ψ 2 ) =υr Pr1 ∇ 4 Ψ 2 − β r Pr1 RaT 2 (4)
∂t ∂x
∂T2
+ℑ( Ψ 2 ,T2 ) =α r ∇ 2T2 (5)
∂t
∂f ∂g ∂f ∂g
avec : ℑ ( f , g ) = − α r = α 2 α1 , υr = υ2 υ1 , β r = β 2 β1
∂z ∂x ∂x ∂z ,
Les fonctions de courant sont définies comme : u1,2 = ∂Ψ1,2 ∂z and w1,2 = − ∂Ψ1,2 ∂x
.
Le nombre de Rayleigh thermique est donné par: RaT = g β1q′H ′4 λ1α1υ1 .

3. SOLUTION NUMERIQUE
La méthode des volumes finis est utilisée. L’algorithme SIMPLER est employé pour
résoudre les équations. Les critères de convergence doivent conserver la masse et l'énergie
globalement et localement. Afin d'assurer que les résultats soient indépendants, différentes
mailles sont examinées. Des maillages non uniformes très fins près des frontières sont adoptées.
On assume que la solution est convergée quand l'erreur est moins de 10-7.

4. SOLUTION ANALYTIQUE
Dans la limite d'une cavité peu profonde ( A >> 1 ), les équations gouvernant le présent
problème peut être considérablement simplifié et résolu analytiquement en utilisant
l'approximation d’écoulement parallèle Ψ1,2 ( x, z ) ≈ Ψ1,2 ( x ) , T1,2 ( x, z ) ≈ CT x + θ1,2 ( x ) et dans la
partie centrale de la cavité. CT est la constante inconnue du gradient de la température,
respectivement dans la direction de x (voir par exemple la référence. (Mamou et al. 1998)). Avec
ces approximations, la forme des équations régissantes (2)–(5) peut être réduite à:
∂ 4 Ψ1 ∂z 4 = E (6)

Albi, France du 28 au 30 Août 2007 2

 13èmes Journées Internationales de Thermique

d 2θ1 dz 2 = CT d Ψ1 dz (7)
∂ Ψ 2 ∂z = ( β r υ r ) E
4 4
(8)
d θ 2 dz = ( CT α r ) d Ψ 2 dz
2 2
(9)
Où E = CT RaT . Les équations ci-dessus peuvent être résolues avec les conditions aux limites, pour
rapporter les distributions de fonctions de courant et de températures. Les bilans énergétiques à
chaque section verticale de la cavité donnent l’expression de CT (voir par ex. Kalla et al. 2001).

5. RESULATS ET DISCUSSION
La solution analytique prédise de l'approximation d’écoulements parallèles a été comparée
avec la solution numérique obtenue par un code aux volumes finis. Le problème actuel est régi
par six paramètres sans dimensions notamment, RaT, αr, βr, λr, υr et A. Figure2 (a) et (b) illustre
les profils verticaux de la vitesse, et de la température. On remarque que les résultats numériques
sont en bon accord avec la solution analytique. On observe que le mouvement des deux fluides
hal-00158373, version 1 - 29 Aug 2007

est dans le sens des aiguilles d’une montre et de même intensité et les gradients thermiques sont
de même importance, ce qui est illustré par les isolignes de la fonction de courant et de
température, Figure 2(c) et (d).
0,50 0,50
Num eric Numeric
Analytic Analytic

0,25
0,25
fluid2
fluid2
0,00
z 0,00
z
fluid1
-0,25

fluid1 -0,25

-0,50
-30 -20 -10 0 10 20 30
-0,50
U -0,2 -0,1 0,0 0,1 0,2

(a) Temperature

(b)

(c) (d)
Figure 2: résultats analytique et numérique pour la distribution de: (a) composante des vitesses
horizontales (b) température. Isolignes: (c) fonction de courant (d) température
pour A=8, υr = αr = βr = λr = Pr = 1 et RaT = 10+5

Le mouvement convectif considéré est induit par des forces de poussée thermiques
résultant des conditions aux frontières appliquées. La présence de deux fluides non miscibles
superposés dans la même cavité, engendre une structure d’écoulement composé de deux cellules,
chacune révèle l’intensité de l’écoulement de chaque fluide. Toutefois, en tenant compte des
propriétés thermophysque et à cause du cisaillement au niveau de l’interface entre les deux
fluides, chaque fluide a tendance à imposer le sens d’écoulement dans tout le domaine.
Figure3 montre l’effet du chauffage latéral (a) sur l’intensité de l’écoulement et le transfert
thermique. On note, la présence de diagramme bifurcation fourche et ce pour les deux fluides, ce
qui explique la multiplicité des solutions. En plus, l’existence de solution naturelle et solution

Albi, France du 28 au 30 Août 2007 3

 13èmes Journées Internationales de Thermique

anti-naturelle. On remarque, que l’accroissement de a n’affecte pas le Rayleigh critique (RaTcr)

mais plutôt affecte le Rayleigh supercritique (RaTsup) pour lequel apparition d’une troisième
cellule, pour assurer la liaison au niveau de l’interface, Fig. 3(a). En plus, les gradients
thermiques sont plus accentues ce qui révèle un transfert thermique plus important, Fig. 3(b).
4 0,50
a=
0.1
0.1
0.4
2 0,25
0.4
0.2
0.2 fluide2
z
Ψ1, Ψ2

0
0,00

a=
0.1 fluide1

-2
0.2
-0,25 0.4
0.6
0.8
-4
0 4 4 4 4 4 4 -0,50
1x10 2x10 3x10 4x10 5x10 6x10
-0,15 -0,10 -0,05 0,00 0,05 0,10 0,15
RaT
Temperature
(a) (b)
hal-00158373, version 1 - 29 Aug 2007

α β λ υ
Figure 3: effet de a : r= r= r= r =1 et A=8, (a) Diagramme de bifurcation (b) profiles de
température au plan médian vertical pour RaT = 10+5

Le paramètre le plus important dans ce problème est le rapport des viscosités (υr). On
expose sur la Figure4, l’effet de υr sur l’intensité de l’écoulement et le transfert thermique. On
note, en plus de la présence de diagramme de bifurcation fourche, que l’intensité de l’écoulement
décroît avec l’accroissement de la viscosité, inversement le mouvement du fluide2 est plus
intense lorsque sa viscosité est plus faible, ce qui est illustré par la figure 4(a). Les profils des
vitesses horizontales dans le plan médian vertical illustrent bien ce constat, Fig. 4(b). Également,
l’accroissement de υr affecte non seulement le (RaTcr) mais aussi le (RaTsup). En effet, le fait de
croître υr, fait translaté la valeur du (RaTsup), ce qui est conséquence de l’affaiblissement de
l’intensité de l’écoulement.
5,0 0,50
nu=
0.5
0.5
1
2,5 0,25
1
2
2 z fluid2
Ψ1, Ψ2

0,0
0,00
nu =
4 fluid1
2
-2,5
-0,25
1
0.5
0.25

-5,0
0 4 4 4 4 -0,50
1x10 2x10 3x10 4x10
-60 -40 -20 0 20 40 60
RaT
u
(a) (b)

fonctions de courant (υr =0.5) isothermes (υr =0.5)

Figure 4: effet de υr : a=0.1, β r= λr= υr =1 et A=8, (a) Diagramme de bifurcation (b) profiles
de vitesses horizontales dans le plan médian vertical pour RaT = 10+5

Albi, France du 28 au 30 Août 2007 4

 13èmes Journées Internationales de Thermique

6. CONCLUSION
Le problème de la convection naturelle dans un domaine horizontal, occupé par deux
fluides non miscible superposés, soumis aux flux croisés de chaleur a été étudié. Des solutions
analytiques pour la fonction de courant et la température, dans la région centrale de chaque
fluide, sont obtenues en utilisant une approximation d’écoulements parallèles. Les résultats
analytiques indiquent que pour la convection supercritique, quand les contributions de
cisaillement et thermiques déstabilisent, plus de trois solutions différentes, un de ces solutions
étant instables, sont possibles à un ensemble donné de paramètres régissant à condition que la
droite soit rendue plus haute qu'une valeur critique qui dépend de a, de υr et de RaT. L'existence
de solutions multiples à proximité de a= 0 a été également démontrée. On observe une bonne
concordance entre les prévisions analytiques et les simulations numériques.

Nomenclature λ1,2 thermal conductivity of the fluid1,2

A aspect ratio of the enclosure, L ′ H ′
λr ration of the thermal conductivity
g gravitational acceleration
hal-00158373, version 1 - 29 Aug 2007

υ1,2 kinematics viscosity of the fluid1,2

q′ constant heat flux, per unit area
dimensionless time, t ′ (σ H ′ α )
t 2 υr ration of the kinematics viscosity, ρ1,2
density of the fluid1,2
T dimensionless temperature
ρr ρ 2 ρ1
( x, z ) dimensionless coordinate system ration of the density,

(u 1,2 , w1,2 ) dimensionless velocity terms Superscript

′ dimensionless variable
Greek symbol
α1,2 thermal diffusivity (fluid1,2), Subscript
αr ratio of thermal diffusivity
0
T
refers to the center of the cavity
temperature
β1,2 thermal expansion coefficient fluid1,2 1,2 fluid 1 or 2
r ratio
βr ration of the thermal expansion coefficient,
β 2 β1

References
[1] Kalla, L., Vasseur, P., Bennacer, R., Beji, H. and Duval, R. (2001) “Double diffusive
convection in a horizontal porous layer salted from the bottom and heated horizontally”,
International communication of Heat and Mass Transfer, Vol 28, pp. 1-10.
[2] Lauriat, G. and Prasad, V. (1987) “Natural convection in a vertical porous cavity: a numerical
study for Brinkman-extended Darcy formulation”, Journal of Heat Transfer, Vol 109, pp.
688-696.
[3] Mohamad, A. A. and Bennacer, R. (2001) “Natural convection in a confined saturated porous
medium with horizontal temperature and vertical solutal gradients”, International Journal of
thermal Science, Vol 40, pp. 82-93.
[4] Mamou, M., Vasseur, P. and Bilgen, E. (1998) “Double-diffusive convection instability in a
vertical porous enclosure”, Journal of Fluid Mechanics, Vol 368, pp. 263-289.
[5] Vasseur, P., Robillard L. and Anochiravani, I. (1986) “Natural convection in a shallow cavity
heated from the side with a uniform heat flux”, Chemical Engineering Communication, Vol
46, pp. 129-146.

Albi, France du 28 au 30 Août 2007 5

View publication stats