N◦: / 2 0 2

Université Sultan Moulay Slimane

Faculté Polydisciplinaire
Département de Physique
- Béni Mellal --

Master : INGÉNIERIE DES ÉNERGIES RENOUVELABLES

& ÉFFICACITÉ ÉNERGÉTIQUE

Option : Procédés Thermique et Valorisation d’Energie

Projet de Fin d’Etudes

Etude numérique et analytique de la convection mixte dans une

cavité rectangulaire simplement entrainée et chauffée par le bas

Réalisé par : KCHIBEL Dalal

Encadré par : Pr. LAMSAADI Mohamed
Soutenu le 27 / 09 / 2022 devant le jury composé de :
Pr. Soufiane Belhouideg (FPBM), Président
Pr. Mohamed Lamsaadi (FPBM), Encadrant
Pr. Hassan El Harfi, Examinateur

Année universitaire 2021-2022

 Table des matières

Liste des figures : ..................................................................................................................... 5

Liste des tableaux : ................................................................................................................... 5

Nomenclature : ......................................................................................................................... 6

Remerciement........................................................................................................................... 7

Introduction générale................................................................................................................ 8

Chapitre I : Etude Bibliographique ........................................................................................ 9

I. Transfert de chaleur par convection : ............................................................................. 9

1. Convection naturelle : ............................................................................................. 9

2. Convection forcée : ................................................................................................. 9

3. convection mixte :................................................................................................. 10

II. Revue bibliographique : ........................................................................................... 10

Chapitre II : Géométries-modèles mathématiques .............................................................. 18

I. Introduction : ................................................................................................................ 18

II. Description physique et géométrique du problème .................................................. 18

1. Hypothèses simplificatrices : ................................................................................ 19

2. Hypothèse de Boussinesq : ................................................................................... 19

III. Formulation mathématique du problème ................................................................. 20

1. Equations gouvernantes ........................................................................................ 20

2. Les conditions initiales et aux limites :................................................................. 21

3. Formes adimensionnelles...................................................................................... 21

4. Equations gouvernantes en formulation secondaire adimensionnelle : ................ 23

5. Transferts thermique ............................................................................................. 24

IV. Conclusion :.............................................................................................................. 24

Chapitre III : Méthode numérique ....................................................................................... 25

I. Introduction .................................................................................................................. 25

II. Méthode des différences finis................................................................................... 25

Page | 2
 III. Discrétisation des équations gouvernantes ............................................................... 26

1. Termes transitoires................................................................................................ 27

2. Termes convectifs ................................................................................................. 27

3. Termes diffusifs et autres que diffusifs ................................................................ 27

IV. Résolution des équations gouvernantes .................................................................... 28

1. Méthode des directions alternées (ADI) ............................................................... 28

2. Equations de la vorticité, de l’énergie .................................................................. 29

3. Équation de Poisson pour la fonction de courant ................................................. 31

V. Conclusion ................................................................................................................ 33

Chapitre IV : Etude analytique ............................................................................................ 34

I. Introduction : ................................................................................................................ 34

II. Approximation d’écoulement parallèle : .................................................................. 34

III. Equations gouvernantes approximées ...................................................................... 35

IV. Solution analytique ................................................................................................... 35

1. Résolution de l’équation de conservation de quantité de mouvement : ............... 35

2. Résolution de l’équation de l’énergie : ................................................................. 36

3. Détermination de la constante C ........................................................................... 37

4. Détermination de nombre nusselt : ....................................................................... 37

Chapitre V : Résultat Et Discussions................................................................................... 38

I. Introduction .................................................................................................................. 38

II. Effet du Maillage ...................................................................................................... 38

III. Effet de rapport de forme : ....................................................................................... 38

IV. Profils des lignes de courant, isothermes ................................................................. 39

V. Profil de vitesse, de fonction de courant et la température le long de la section

verticale de la cavité :..................................................................................................................... 40

1. Profil de vitesse : .................................................................................................. 40

2. Température : ........................................................................................................ 40

3. Fonction de courant .............................................................................................. 41

Page | 3
 VI. Influence des paramètres de contrôles sur le transfert de chaleur ............................ 42

1. Effet du nombre de Peclet sur Nusselt :................................................................ 42

2. Effet du nombre de Rayleigh sur Nusselt : .......................................................... 42

VII. Influence des paramètres de contrôles sur l’intensité de l’écoulement au centre de

la cavité : 43

1. Effet de nombre de Peclet sur l’intensité de l’écoulement au centre de la cavité :

43

2. Effet du nombre de Rayleigh sur 𝚿c : .................................................................. 44

Conclusion générale : ............................................................................................................. 45

Bibliographie : ........................................................................................................................ 46

Page | 4
 Liste des figures :

Figure 1 : Modèle physique et système de coordonnées .................................................................... 18

Figure 2 : Maillage en (i,j) et repère cartésien en ( x,y ) .................................................................... 26
Figure 3: illustration du maillage de la méthode ADI ...................................................................... 28
Figure 4 : Influence du rapport de forme sur le taux de transfert de chaleur ..................................... 39
Figure 5 : Isothermes.......................................................................................................................... 39
Figure 6 : Lignes de courant .............................................................................................................. 40
Figure 7: profil de la vitesse U au long de la section verticale de la cavité avec Pe=50 et Ra varie . 40
Figure 8 : profil de température au long de la section de la cavité avec Pe = 50 et Ra varie ............ 41
Figure 9 : Fonction de courant au long de la section de la cavité avec Pe =50 et Ra varie .............. 41
Figure 10: nusselt en fonction de peclet............................................................................................. 42
Figure 11 : Nusselt en fonction de Rayleigh ...................................................................................... 43
Figure 12: c en fonction de peclet .................................................................................................. 44
Figure 13: c en fonction de Rayleigh ............................................................................................... 44

Liste des tableaux :

Tableau 1: les variables caractéristiques du problème ........................................................... 22

Tableau 2: effet du maillage sur les resulats nimériques ....................................................... 38

Page | 5
 Nomenclature :

H′ Hauteur de la cavité
L′ Longueur de la cavité
u′ Composante dimensionnelle longitudinale du
v′ vecteur vitesse
Composante dimensionnelle transversale du
p′ vecteur vitesse
Pression dimensionnelle
t′ Temps dimensionnel
T’ Température dimensionnelle
T0 ′ Température de référence
u Composante adimensionnelle longitudinale du
v vecteur vitesse
Composante adimensionnelle transversale du
p vecteur vitesse
Pression adimensionnelle
τ Temps adimensionnel
T Température adimensionnelle
g Pesanteur
A Rapport de forme
λ Conductivité thermique
Cp Capacité calorifique
υ Viscosité cinématique
α Diffusivité thermique
ρ Masse volumique
β Coefficient de dilatation thermique
ΔT Facteur d’échelle pour la température
 Vorticité adimensionnelle
Ψ Fonction de courant adimensionnelle
Ψc Fonction de courant adimensionnelle au centre
de la cavité
Pr Nombre de Prandtl
Ra Nombre de Rayleigh
Pe Nombre de Peclet
𝐍𝐮 Nombre de Nusselt moyen

Page | 6
 Remerciement

Au terme de ce modeste travail, je saisis cette merveilleuse occasion pour exprimer mes
profondes gratitudes et mes sincères remerciements à toute personne qui, de près ou de loin, m’est
aidée à surpasser toutes les difficultés rencontrées pendant la période de mon projet de fin d’étude.

Je tiens à remercier vivement mon encadrant Monsieur LAMSAADI Mohamed pour ces
précieux conseils afin de réussir mon projet de fin d’étude.

Enfin je remercie chaleureusement les membres du jury pour leur présence et leur
acceptation d’évaluer mon travail.

Page | 7
 Introduction générale

Cette mémoire est consacré à la résolution des problèmes de la convection mixte dans une
cavité de forme rectangulaire remplie d’un fluide newtonien et dans le paroi inférieure exposé a un
flux thermique et le paroi supérieure mobile d’une vitesse constante et les deux autre paroi droite et
gauche sont supposé adiabatique

Cette mémoire est structurée en cinq chapitres :

o Le premier chapitre est consacré à l’introduction au problème étudié et à la connaissance de

quelques travaux antérieurs concernant la convection mixte.
o Le deuxième chapitre est consacré à la mise en place du problème dans le cadre physique et
mathématique.
o Dans le troisième chapitre, on décrive la méthode numérique utilisée pour la résolution des
équations mathématiques adimensionnelles obtenues dans le deuxième chapitre.
o Dans le quatrième chapitre on présente la solution analytique basée sur l’approximation de
l’écoulement parallèle.
o Les résultats analytiques obtenus lors de cette étude sont présentés et discutés dans le
cinquième chapitre, où on a consacré une partie à l’analyse de ces résultats influencés par les
différents paramètres qui contrôlent le système.
o Enfin une conclusion générale vient résumer l'essentiel des résultats obtenus, quelques
perspectives sont présentées en dernier qui pourraient constituer une suite intéressante à
cette étude.

Page | 8
 Chapitre I : Etude Bibliographique

I. Transfert de chaleur par convection :

La convection est un mode de transport d’énergie par l’action combinée de la conduction, de

l’accumulation de l’énergie et du mouvement du milieu. C’est le mécanisme le plus important de
transfert d’énergie entre une surface solide et un liquide ou un gaz. Le transfert d’énergie par
convection d’une surface dont la température est supérieure à celle du fluide qui l’entoure s’effectue
en plusieurs étapes.

D’abord la chaleur s’écoule par conduction de la surface aux molécules du fluide adjacentes.
L’énergie ainsi transmise sert à augmenter la température et l’énergie interne de ces molécules du
fluide. Ensuite les molécules vont se mélanger avec d’autres molécules situées dans une région à
basse température et transférer une partie de leur énergie. Dans ce cas l’écoulement transporte le
fluide et l’énergie.

L’énergie étant emmagasinée dans les molécules du fluide, elle est transportée sous l’effet de leur
mouvement. La transmission de chaleur par convection est désignée, selon le mode d’écoulement
du fluide, par convection libre et convection forcée. Lorsqu’il se produit au sein du fluide des
courants dus simplement aux différences de température, on dit que la convection est naturelle. Par
contre si le mouvement du fluide est provoqué par une action externe, telle une pompe ou un
ventilateur, le processus est appelé convection forcée. Si les deux modes de convection sont
importants, la convection est dite mixte.

1. Convection naturelle :

C’est la convection ou l’écoulement se fait naturellement car le mouvement des molécules

provoqué par la différence de masse volumique du fluide en différents point du milieu causé par la
distribution de la température, engendrant ainsi son propre mouvement en créant des forces
d’Archimède. Donc la différence de température, avec la force de gravité, crée une force de
flottabilité, qui par conséquent crée une différence de quantité de mouvement

2. Convection forcée :

La convection forcée dans laquelle le mouvement de fluide est provoqué par l’action des forces
extérieures telle qu’un mécanisme comme une pompe ou un ventilateur…etc., qui lui imprime des
vitesses de déplacement assez importantes. C’est ce mode qui est généré lorsque l’on veut améliorer
c'est-à-dire augmenter l’échange thermique en conséquence, l’intensité du transfert thermique par
convection forcée sera en liaison directe avec le régime de mouvement du fluide.

Page | 9
 3. convection mixte :

La convection mixte correspond au couplage de deux phénomènes, à savoir la convection naturelle

et forcée qu’on peut qualifier de cas général de la convection, quand le fluide est soumis
simultanément à une augmentation de pression et à une différence de température, on peut citer les
exemples de l’écoulement atmosphérique et les courants marins.

II. Revue bibliographique :

L’étude des transferts de chaleur par convection mixte dans les cavités fermées suscite
l’intérêt de nombreux scientifiques et industriels. Les recherches menées dans ce domaine
s’étendent sur plusieurs années. La revue bibliographique que nous allons présenter montre
l’importance de ces études ainsi que l’existence d’un nombre considérable de travaux théoriques,
expérimentaux et numériques consacrés à l’étude des transferts de chaleur par convection mixte et
naturelle dans des milieux fluides contenus dans des cavités de différentes formes.

Travaux antérieurs sur la convection naturelle en milieux Newtoniens

dans différentes géométries :

Ganzarolli et Milanez [1] ont étudié numériquement la convection naturelle à l'intérieur

d'une enceinte rectangulaire avec un flux de chaleur constant à la paroi de fond et
symétriquement par rapport à froid à partir de ses parois latérales verticales tandis que la paroi
supérieure est maintenue à une condition limite adiabatique. Ils ont présenté les résultats sous
forme de nombre de Nusselt (Nu) en fonction du nombre de Rayleigh (Ra) dans un intervalle
de 103 à 107, tandis que le rapport d’aspect a varié de 1 à 9. Ils ont également signalé que, à un flux
de chaleur prescrit, et les isothermes rationalise occupées plus uniformément à l'intérieur de toute la
cavité à une faible valeur de Rayleigh.

Chu et Churchill [2] qui ont étudié la convection naturelle au sein d'une enceinte rectangulaire
bidimensionnelle munie d'une source de chaleur discrète. Avec des parois horizontales
thermiquement isolées, une paroi verticale froide isotherme et une source de chaleur installée sur
l'autre paroi verticale, ils ont obtenu numériquement la position de la source qui optimise le
transfert thermique, en fonction du nombre de Rayleigh. Cette position optimale se trouve près du
centre de la paroi verticale, mais glisse légèrement vers le bas à mesure que le nombre de Rayleigh
augmente.

Lamsaadi et al. [3,4], ont mené une étude analytique et numérique des transferts thermiques
par convection naturelle dans une cavité rectangulaire horizontale confinant des fluides non
Newtoniens de type Ostwald-De Waele et chauffée à l’aide de flux uniforme. Dans tous les cas,les

Page | 10
champs d’écoulement et de température et les transferts thermiques qui en découlent sont trouvés
plutôt sensibles au comportement non-Newtonien qu’au nombre de Prandtl, pour les grandes
valeurs de celui-ci ( 𝑃𝑟 ≥ 100). Ainsi, comparé au cas Newtonien, le comportement pseudo-
plastique tend à favoriser l’apport de la convection thermique alors que celui dilatant tend à réduire
cet apport. Aussi, pour le type de chauffage adopté et dans le cas de grands rapports de forme (𝐴 >>
1), les résultats analytiques semblent s’accorder parfaitement avec ceux numériques, ce qui
constitue une validation de l’hypothèse de l’écoulement parallèle et le code de calcul numérique
élaboré du moins dans la gamme de valeurs des nombres adimensionnels gouvernants.

En chauffant du bas [3], par comparaison au cas Newtonien, une diminution de l’indice
de comportement précipite l’apparition de la convection alors qu’une augmentation de ce même
paramètre produit l’effet contraire.

Avec un chauffage latéral [4], ils ont observé que la convection ne présente pas de seuil
d’apparition, contrairement au cas précédent qui nécessite, lui, un seuil critique.

F. ZOUIRI et al [5] Cette étude traite l’analyse numérique de la convection naturelle laminaire au
sein d'une cavité carrée dont les parois verticales sont maintenues à une température constante alors
que les parois horizontales sont isolées thermiquement, à l’exception d’une fraction occupant 20% à
80% de la surface inférieure de l’enceinte et centrée par rapport à celle-ci, qui est maintenue à une
température constante et uniforme, supérieure à celle des parois latérales grâce à une source de
chaleur placée en contact de cette paroi. La résolution des équations régissant l'écoulement et le
transfert thermique est approchée par la méthode des volumes finis et l'algorithme SIMPLER est
adopté pour traiter le couplage vitesse-pression. L'étude se focalise sur l'influence occasionnée par
les variations du nombre de Rayleigh et de la longueur de la fraction chauffée, sur la structure de
l'écoulement et du transfert thermique au sein de la cavité remplie entièrement d’un fluide
newtonien incompressible.

Les résultats obtenus montrent que le fait d'augmenter le nombre de Rayleigh favorise le
transfert thermique dans la cavité étant donné l’augmentation des courants de convection et
donc, la vitesse. En revanche, il y a diminution de la surface de la cavité exposée aux
températures élevées. De même, l’augmentation de la largeur de la source chauffée et donc,
lasurface d’échange, améliore le transfert thermique et fait augmenter la surface exposée aux
températures élevées.

El Alami SEMMA et al [6] ont fait une étude sur la convection naturelle dans une cavité
triangulaire semi-poreuse remplie d’air. Les résultats numériques sont présentés pour le champ de

Page | 11
vitesse, le champ de température et le transfert thermique à travers les parois isothermes. L’effet
d’une couche poreuse attachée à la paroi verticale est examiné. Ils montrent que les faibles
épaisseurs affectent nettement le transfert thermique et le mode d’écoulement. L’effet de la
perméabilité est analysé, autrement dit, Les résultats montrent qu’une forte dépendance du transfert
thermique à travers les parois isothermes et la perméabilité du milieu poreux. Une corrélation du
nombre de Nusselt en fonction de l’épaisseur de la couche poreuse et du nombre de Rayleigh est
identifiée.

R. Zarrit et al [7] étudient numériquement la convection naturelle instationnaire et

bidimensionnelle dans une cavité rectangulaire inclinée et chauffée différentiellement. Le
nombre de Rayleigh varie entre 103 ≤ Ra ≤ 104, le rapport de forme (1 ≤ A ≤ 6) et l’angle
d’inclinaison θ comprise entre θ=0° et 180°. Les calculs mettent en évidence l’influence de
l’angle d’inclinaison, le nombre de Rayleigh et le rapport de forme sur la structure de
l’écoulement et le transfert de chaleur. Les résultats trouvés révèlent l’existence d’un angle
critique cr dépendant du rapport de forme A et du nombre de Rayleigh Ra. La valeur maximale de
Nusselt moyen se trouve dans l’intervalle (0    cr) pour A  2 et dans l’intervalle (cr   
90) pour A  2. En général, au voisinage de ( 0), l’écoulement est multicellulaire (rouleaux
longitudinaux formant des cellules contrarotatives). Pour ( cr), l’écoulement est monocellulaire.
Dans le cas où ( 180), le fluide stagne et une stratification thermique est bien marquée.

Une étude analytique et numérique de la convection naturelle d'un fluide à double diffusion
contenu dans une fente rectangulaire soumise à des flux de chaleur et de masse uniformes le
long des côtés verticaux, a été présentée par M. Mamou et al [8]. Les paramètres directeurs du
problème à l'étude sont le nombre de Rayleigh thermique, RaT; rapport de flottabilité, N;
Numéro de Lewis, Le; Nombre de Prandtl, Pr et le rapport d'aspect de la cavité, A. Dans la
première partie de l'étude analytique, une analyse d'échelle est appliquée aux deux cas extrêmes
d'écoulements entraînés par transfert de chaleur et transfert de masse. Dans la deuxième partie, une
solution analytique, basée sur l'approximation d'écoulement parallèle, est rapportée pour les
enceintes hautes (A≫1). Des solutions pour les champs d'écoulement, les distributions de
température et de concentration et les nombres de Nusselt et Sherwood sont obtenues en termesde
paramètres régissant le problème. Dans les limites des flux thermiques et massique, un bon accord
est obtenu entre la prédiction de l'analyse à l'échelle et celles de la solution analytique.

Les solutions numériques sont basées sur les équations de gouvernance complètes pour les
écoulements bidimensionnels et couvrent la plage 1≤ RaT ≤107, 0≤N≤105,10-3≤Le≤103,

Page | 12
1≤A≤20 et Pr = 7. Un bon accord est trouvé entre les prédictions analytiques et la simulation
numérique.

Sarris et autres [9] qui ont effectué une étude numérique de la convection naturelle dans les
réservoirs rectangulaires chauffés localement au-dessous. Basé sur des prévisions numériques,
les effets du Ra et la géométrie de la source de chauffage et du réservoir sur les structures
d'écoulement et le transfert thermique ont été étudiés pour le nombre de Rayleigh dans la
gamme 102 à 107. Leur étude a indiqué que l'augmentation de l'intensité de circulation
d'écoulement augmente avec une augmentation de longueur de la source de chauffage.

Alloui et al. [10] ont utilisé Le modèle de Darcy avec l’approximation de Boussinesq,
pour étudier la convection naturelle dans un milieu poreux saturé par un fluide binaire. La
géométrie considérée est une cavité carrée dont la paroi inferieure est chauffée de manière
isothermique, la paroi supérieure est refroidie `a une température constante et toutes les autres
parois sont adiabatiques. La m´méthode des volumes finis a été employée pour la discrétisation des
´équations du système. Les paramètres qui régissent le problème sont le nombre de Rayleigh
thermique, Ra le nombre de Lewis, Le, le rapport de flottabilité N, la longueur sans dimension de la
plaque inferieure, le rapport de forme de la cavité A, la porosité normalisée du milieu poreux et la
position relative de l’´élément chauffant par rapport à l’axe vertical de la cavité. Deux principaux
modes convectifs sont étudies, à savoir la convection `a une et deux cellules, et leurs
caractéristiques sont d´écrites. L’existence possible d’´écoulements tri-cellulaires est également de
montrée. La fonction de flux maximale et le nombre global de Nusselt et Sherwood sont présentés
en fonction des paramètres externes. L’existence de trois solutions en régime permanent pour un
ensemble donné de paramètres est de montrée.

Sath et al. [11] ont examiné numériquement le transfert de chaleur par convection
naturelle thermique dans une cavité rectangulaire partiellement occupée par une couche
poreuse verticale. En utilisant le modèle de Darcy-Brinkman dans la couche poreuse et
l’équation de Navier-Stokes dans le milieu fluide, les auteurs ont montré que le terme de
Brinkman permet de satisfaire la condition d’un glissement aux parois solides et à l’interface.

Revnic et Grosan [12] ont étudié le phénomène de la convection naturelle dans une couche
poreuse à géométrie rectangulaire de dimension infinie, l’étude est basée sur le modèle de

Nield et Kuznetsov [13], la couche poreuse utilisée est constituée d’une phase solide et des pores
saturés par un fluide, la phase solide peut être considérée comme étant un autre milieu poreux avec

Page | 13
les mêmes propriétés physiques. Le nombre de Rayleigh pris dans cette investigation n’excède pas
103 et les surfaces latérales sont maintenues à des températures constantes et différentes.

Horton et Rogers [14] étaient parmi les premiers qui ont traité le phénomène de la convection
naturelle dans une couche poreuse isotrope à porosité uniforme, saturée par un fluide et chauffée par
le bas.

Prasad et Tuntomo [15] ont utilisé l’équation de Darcy-Forchheimer pour étudier la

convection naturelle dans une cavité poreuse rectangulaire pour une gamme étendue des
principaux paramètres. Ils ont analysé l’influence du paramètre inertiel sur l’écoulement et le
transfert de chaleur. Ils ont montré que le nombre de Nusselt diminue avec l’augmentation du
paramètre inertiel.
Travaux antérieurs sur la convection mixte en milieux Newtoniens dans des différentes
géométries :

M. Lamsaadi et al [16] ont traité La convection mixte dans une cavité rectangulaire allongée
de rapport d’aspect, A = 12 ,confinant un fluide newtonien de nombre de Prandtl, Pr = 7, est
étudiée analytiquement et numériquement, dans le cas où les parois horizontales, dont celle
supérieure est mobile, sont adiabatiques, et celles verticales sont soumises à une densité de flux de
chaleur uniforme. Les résultats obtenus montrent la validité de l’approximation de
l’écoulement parallèle et l'effet considérable des nombres de Richardson (Ri = Gr/Re2) (0,1 ≤
Ri ≤ 105) et de Reynolds (0,1 ≤ Re ≤ 10) sur les intensités de l’écoulement et du transfert
thermique

S. BOUABDALLAH et al [17] ont étudié numériquement la convection mixte laminaire

stationnaire dans une cavité carrée. La cavité est remplie par l’air (Pr = 0.71) et contenant deux
portes d’entrée et de sortie de l’air. L’une des parois de la cavité est soumise à une température
constante Tc (paroi inférieure), tandis que les autres parois sont maintenues adiabatiques. Les
équations régissent l’écoulement de convection mixte à l’intérieur de la cavité sont résolus par les
méthodes des volumes finis. Le nombre de Grashof fixé à Gr=105, et le nombre de Reynolds (Re)
varie de telle sorte que le nombre de Richardson prend les valeurs : Ri =0.1, 1, 2, 5, 10 et 20 (Ri =
Gr/Re2). Les résultats obtenus ont été validés avec l’expérimental et d’autres travaux numériques.
L'influence des paramètres de l’écoulement sur la structure d'écoulement et de champ thermique est
présentée. L’effet de variation des positions des portes d'entrée et de sortie à la cavité a été
également examiné. Nos résultats de simulation, permettent de contrôler l’intensité de la convection
mixte et les positions des ouvertures de la cavité, pour un bon refroidissement des montages
électriques.

Page | 14
BOUCHAIR et al [18], ont étudié numériquement le transfert de chaleur par convection mixte
dans d’une cavité carrée remplie par un fluide conducteur de l’électricité et soumise à un champ
magnétique externe. L’étude numérique est analysée par la méthode des volumes finis et
l’Algorithme de SIMPLER. On fin sont examiner l’influence du nombre de Richardson (Ri) et le
nombre de Hartman (Ha) sur l'intensité de l’écoulement et le transfert de chaleur à l'intérieur de la
cavité. . Ils ont montré que le nombre de Nusselt diminue avec l’augmentation du nombre de
Hartman et augmente avec une augmentation du nombre de Richadson.

Joo - Sik (1998) [19], a présenté une étude numérique de la convection mixte de l’écoulement
de l’air, entre deux cylindres concentriques maintenus à des températures constantes puis
différentes. L’écoulement forcé est induit par le cylindre extérieur froid qui tourne lentement
avec une vitesse angulaire constante. Les investigations ont été faites pour différentes
combinaisons de Ra, Re et du rapport du diamètre du cylindre intérieur sur l’écart des diamètres
intérieur et extérieur σ, pour Ra < 105, Re < 1500 et 0.5 <σ <5. Les résultats montrent qu’il y a trois
types d’écoulements : écoulements à deux tourbillons, écoulements à un tourbillon et à aucun
tourbillon. La circulation du fluide dans la direction de rotation des cylindres diminue en
augmentant le nombre de Rayleigh. Le transfert de chaleur global à la paroi est rapidement réduit
lorsqu’on atteint la valeur critique du nombre de Reynolds à la transition.

Chamkha et al [20], ont étudié le problème de l'instabilité, laminaire on convection mixte

dans une cavité carrée en présence de génération de chaleur interne dépendant de la température.
Les parois horizontales supérieure et inférieure de la cavité sont adiabatiques alors que les parois
verticales gauche et droit sont différentiellement chauffées. La paroi gauche verticale est mobile
avec une vitesse constante tandis que toutes les autres parois sont fixées. Un champ magnétique
uniforme est appliqué dans la direction horizontale normale à la paroi en mouvement. Les équations
avec conditions aux limites sont résolues numériquement par la méthode des volumes finis. Les
résultats obtenue présent un bon accord avec les travaux publiée.

Louaraychi et al. [21] ont étudié La convection mixte dans des cavités rectangulaires
horizontales à simple et double couvercle, remplies d'un fluide newtonien et soumises à un flux de
chaleur uniforme le long de leurs petits côté verticaux. L’investigation est traitée
numériquement et analytiquement. La méthode des volumes finis FVM est utilisée pour
résoudre les équations directrices complètes pour lesquelles l'approximation de Boussinesq est
adoptée. L'algorithme SIMPLER est adopté pour la résolution du problème de couplage
pression/vitesse. L'approche analytique repose sur l'hypothèse d'un flux parallèle et valable dans le
cas des enceintes peu profondes. Un bon accord entre les deux approches est observé dans les

Page | 15
plages explorées des nombres de Peclet et de Rayleigh. Les effets de ces paramètres sur les
caractéristiques d'écoulement et de transfert de chaleur sont analysés pour les deux types des
cavités. Les zones caractérisant la 16 prédominance des convections naturelles et forcées ainsi que
lorsque les deux phénomènes en concurrence (convection mixte) sont délimités. Il est constaté que
la transition d'un régime dominé à un autre dépend du rapport Ra/Pe3.

ZAHIRI et al [22], ont étudié numériquement et analytiquement la convection mixte dans

une cavité rectangulaire remplie d'un fluide newtonien ayant un nombre de Prandtl constant
Pr = 7, en présence d'un champ magnétique. Les parois horizontales sont mobiles et
adiabatiques, tandis que les parois verticales sont soumises à une densité de flux de chaleur
uniforme. Les résultats obtenus montrent que l’augmentation de nombre Hartman engendre une
diminution du transfert de chaleur.

Rahman et al. [23] qui ont étudié numériquement la convection mixte dans une cavité ventilée en
présence d'un générateur de chaleur sous forme de solide circulaire. Les parois sont considérées
adiabatiques. Ils ont constaté que le nombre de Nusselt moyen le long de la surface de la source de
chaleur augmente avec l’augmentation des nombres de Reynolds et de Prandtl.

ADNANI et al. [24] ont étudié numériquement la convection mixte dans une cavité carrée dont les
parois horizontales sont en mouvement et qui contient un bloc froid à l’intérieur. La modélisation
indique que le transfert de chaleur et le champ de vitesse augmentent avec
l’augmentation du Ri.

Chang [25] a examiné la convection mixte laminaire dans un canal vertical, dont la section est finie
et soumise à une variation linéaire de la température. Les variations de la vitesse locale, de la
température, des valeurs locales et moyennes de Nusselt sont tracées pour illustrer les influences
des paramètres étudiés, notamment le nombre de Reynolds, le nombre de Grashof et le degré de
variation de la température de la paroi. Les résultats suggèrent que le nombre de Nusselt moyen
augmente avec le nombre de Reynolds et le nombre de Grashof.

RADHAKRISHNAN et al. [26] ont analysé expérimentalement et numériquement de la

convection mixte dans une enceinte rectangulaire contenant des blocs chauffants de tailles
différentes. Les simulations prouvent que la distribution de la température donne une approche de la
gestion de l'alimentation entre les chauffages, de sorte que le réchauffeur le plus «froid» puisse être
chargé pour maximiser la dissipation de chaleur totale, pour un excès de température prescrit, pour
tous les réchauffeurs.

Page | 16
Basak et al. [27] ont étudié numériquement la convection mixte dans une cavité carrée avec
un échauffement uniforme et non uniforme de la paroi inferieure. Les deux parois verticales
sont maintenues à une température froide constante tandis que la paroi supérieure est considérée
adiabatique et se déplace avec une vitesse uniforme. Ils ont observé que l’intensité de la circulation
de l’écoulement augmente avec l'augmentation de la valeur du nombre de Grashof Gr, qui est due à
l’effet de la force d’Archimède (force de flottabilité). Les résultats montrent que la convection
naturelle ou forcée est basée sur les deux paramètres de contrôle qui sont le nombre de Prandtl Pr et
le nombre de Richardson Ri dont l’expression est Gr/Re2. Re représente le nombre de Reynolds.

Prasad et Das [28] ont étudié la convection mixte dans une cavité rectangulaire avec la
méthode des volumes finis (FVM). Les parois verticales sont maintenues à une température
froide tandis que celle supérieure est chaude et la paroi inférieure est froide pour un premier cas et
inversement pour le deuxième cas. L’étude numérique est faite pour des nombres
de Grashof allant de 0 jusqu’à 106. Le nombre de Reynolds est égal à 100 et le rapport de forme est
égal à 0,5. Les résultats sont présentés sous forme des isothermes, des lignes de courants, ainsi que
l’évolution du nombre de Nusselt.

Abu-Nada et Chamkha [29] ont traité la convection mixte dans une cavité carrée inclinée
remplie d’un nanofluide. Les parois gauche et droite verticales sont maintenues adiabatiques alors
que le bas et le haut mobile sont, respectivement, refroidi et chauffé à températures
constantes. Il a été trouvé que la présence de nanoparticules améliore sensiblement le transfert
de chaleur et que cela s’accentue avec l’inclinaison de l’enceinte a valeurs modérées et grandes du
nombre de Richardson.

Oztop et Dagtekin [30] ont étudié numériquement la convection mixte bidimensionnelle et

stationnaire dans une cavité carrée portant deux parois mobiles verticales et chauffées
différentiellement. Les parois horizontales sont maintenues adiabatiques. Trois cas ont été
considérés dépendants de la direction du mouvement des parois (Figure 1.4). Toutes les
simulations ont été réalisées avec des valeurs du nombre de Richardson appartenant à
l’intervalle 0.01 ≤ Ri ≤ 100. Le nombre de Prandtl est égal à Pr = 0.7. Les résultats montrent
que le nombre de Richardson influe sur l’écoulement du fluide et le transfert thermique pour
Ri > 1, et par conséquent, le transfert de chaleur est assez meilleur.

Page | 17
 Chapitre II : Géométries-modèles mathématiques

I. Introduction :

On a étudié dans ce chapitre le phénomène de la convection mixte dans une cavité

rectangulaire entrainée simplement et chauffée par le bas. On transforme ce problème physique à un
modèle mathématique , Ces équations mettent en évidence le comportement des fluides en
convection laminaire mixte. Elles sont fondées selon les lois physiques : de conservation de masse
(continuité), de quantité de mouvement (Navier-Stokes), et d’énergie. Et on détermine les
conduction au limites appropries et pour faciliter la résolution du problème, on fait appel à des
hypothèses simplificatrices , et on conclue le chapitre par un analyse adimensionnelle à ces
équations, ce qui fait savoir des grandeurs essentielles sont les nombres adimensionnels comme le
nombre de Rayleigh et le nombre de prendlt, dans un système de coordonnées cartésiennes compte
tenu de la configuration géométrique considérée.

Enfin, on définit le coefficient de transfert de chaleur

II. Description physique et géométrique du problème

La géométrie du problème étudié est une cavité de forme rectangulaire de langueur L et de

largeur H, dont les parois verticales sont immobiles et adiabatiques tandis que la paroi supérieure
est mobile avec une vitesse U0 et la paroi inférieure est immobile mais chauffée par une densité le
flux thermique.

Figure 1 : Modèle physique et système de coordonnées

Page | 18
 1. Hypothèses simplificatrices :

Pour élaborer un modèle mathématique simplifié ,on va adopter les hypothèses de

simplifications suivantes :

 La troisième dimension de la cavité est supposée assez large de sorte que le problème puisse
être ramené à deux dimensions
 Les vitesses du fluide sont relativement faibles pour considérer l’écoulement comme
laminaire.
 L’écoulement transitoire
 toute réaction chimique, ou source de chaleur ou de masse est absente.
 Le fluide est incompressible et Newtonien.
 Le transfert de chaleur par rayonnement est négligeable
 Les propriétés physiques de fluide ( ,  ,  et  ) sont supposées constantes
 L’enceinte se trouve dans le champ gravitationnel de la Terre
 On néglige le travail, induit par les forces visqueuses et de pression.
 Les propriétés thermo-physiques du fluide sont constantes et sont évaluées à la
température de référence, excepté la densité du fluide, qui varie linéairement avec la
température, dans le terme des forces de volume, selon l'approximation de Boussinesq.

2. Hypothèse de Boussinesq :

On a La masse volumique d’un fluide (liquide, gaz, mélange liquide ou gazeux) est
considérée comme étant constante dans tous les termes des équations de conservation, c’est-à-dire
(𝑇, 𝑃) = 𝜌0(𝑇0, 𝑃0) = 𝜌0 = 𝐶𝑠𝑡𝑒 (𝑇0 et 𝑃0 sont la température et la pression de référence) sauf dans le
terme générateur de la convection (c’est-à-dire dans le terme de poussée d’Archimède).

Dans ce terme 𝜌 est alors une fonction de la température et de la fraction massique :

 = 𝜌0𝑇0[1 – 𝛽(𝑇 – 𝑇0)]

𝜌0 : Masse volumique à la température et la concentration de référence (constante).

𝑇0 : Température de référence (constante) associée à 𝜌0..

1 𝜕𝜌
𝛽 : Coefficient d’expansion thermique (constant) égal :𝛽 = − 𝜌 (𝜕𝑇)
𝑝=𝑐𝑠𝑡

Page | 19
 Le terme 𝛽, est un coefficient possédant les dimensions de l'inverse d'une température, il
varie avec cette dernière d'une façon décroissante dans le cas des gaz. En général, 𝛽 est positif pour
tous les fluides, sauf l’eau en dessous de 4°C.

la masse volumique n'est considérée comme variable que dans le terme moteur ou bien
gravitationnel de l'équation de quantité de mouvement (D. D. Gray, A Giorcini, 1975), au
niveau de tous les autres termes, elle est constante et égale à 𝜌0.

L’approximation de Boussinesq est valable tant que les variations de la densité sont
inférieures de 10% de grandeur de la densité moyenne ((∆𝜌/𝜌) ≤ 0.1).

III. Formulation mathématique du problème

1. Equations gouvernantes

Dans la présente étude le phénomène de la convection mixte est gouvernée par les équations
mathématiques suivantes :

a. Equation de conservation de la masse

Il s’agit d’un principe de la physique, qui traduit la conservation de la masse au cours du

temps pour un élément de volume de contrôle, en absence de puits ou de source. Pour un fluide
saturant incompressible de vitesse moyenne V, il découle l’équation de continuité où la divergence
est nulle, telle que :

∂u′ ∂v′
+ ∂y′ = 0 (2.1)
∂x′

𝑢’ et v’ étant les composantes du champs de vitesse 𝑉⃗ (u’,v’) dans les directions (x’) et (y’)
respectivement.

b. Equation de conservation de quantité de mouvement

L’équation de quantité de mouvement traduit le principe fondamental de la dynamique établie par

Newton pour un fluide incompressible, en faisant introduire le mouvement, la nature du fluide et les
forces qui sont à la cause. Pour un élément de volume de contrôle, l’équation s’écrit :

suivant x :

∂u′ ∂u′ ∂u′ 1 ∂p′ ∂2 u′ ∂2 u′

+ u′ ∂x′ + v′ ∂y′ = − ρ ∂x′ + ν (∂x′2 + ∂y′2 ) (2.2)
∂t

Suivant y :

Page | 20
∂v′ ∂v′ ∂v′ 1 ∂p′ ∂2 v′ ∂2 v′
+ u′ ∂x′ + v′ ∂y′ = − ρ ∂y′ + ν (∂x′2 + ∂y′2 ) + gβ(T′ − T′0 ) (2.3)
∂t

c. Équation de conservation de l’énergie :

En négligeant les effets de la dissipation visqueuse, l’équation d’énergie établie à partir du

premier principe de la thermodynamique s’écrit comme suit :

∂T′ ∂T′ ∂T′ ∂2 T′ ∂2 T′

+ u′ ∂x′ + v′ ∂y′ = α ( ∂x′2 + ∂y′2 ) (2.4)
∂t
𝜆
Où 𝛼 = 𝜌 𝐶 est la diffusivité thermique du fluide et 𝜌 𝐶p et 𝜆 représentent respectivement la
𝑝

capacité calorifique et la conductivité thermique de fluide.

2. Les conditions initiales et aux limites :

a. Conditions initiales :

A un temps initial, le fluide confiné dans la cavité est considéré au repos et à une répartition
homogène de température et de concentration, donc on peut écrire :

A l’instant : t = 0 , on a : u’ = v’ = T’ = 0

b. Conditions aux limites

Dans cette étude la condition d’adhérence est appliquée sur toutes les parois de la cavité, cela
permet d’écrire les conditions aux limites comme suit :

𝜕𝑇′
Pour : x’ = 0 et y’  u’ = v’ = 0 et =0
𝜕𝑥′

𝜕𝑇′
Pour : x’ = L et y’  u’ = v’ = 0 et =0
𝜕𝑥′

𝜕𝑇′ 𝑞
Pour : y’ = 0 et x’  u’ = v’ = 0 et = −𝜆
𝜕𝑦′

𝜕𝑇′ 𝑞
Pour : y’ = H et x’  u’ – u0 = v’ = 0 et 𝜕𝑦′ = − 𝜆

Où q est le flux de chaleur appliqués sur le paroi inférieur de la cavité.

3. Formes adimensionnelles

Les équations précédentes sont mises sous une forme adimensionnelle en utilisant les variables
caractéristiques du problème. La forme adimensionnelles a pour avantage de faciliter la
simplification des équations et de généraliser les résultats. De façon à rendre les équations
précédentes adimensionnelles, la dimension H de la cavité a été choisie comme longueur de
référence. La diffusivité thermique est utilisée pour adimensionnaliser la vitesse, la température est

Page | 21
adimensionnalisée par rapport aux écarts de température ∆𝑇 entre les deux parois horizontales de la
cavité.

Tableau 1: les variables caractéristiques du problème

Coordonnées cartésiennes 𝑥′ 𝑦′
(𝑥; 𝑦) = ( ; )
𝐻 𝐻

Composantes de vitesse u′ v′ 𝛼
(u; v) = (u′′ ; u′′) avec 𝑢′′ = 𝐻
𝑢′′
Temps 𝜏= 𝑡
𝐻

Température 𝑇′ 𝑞𝐻
𝑇 = ∆𝑇 avec ∆𝑇 = 𝜆
p′
Pression P=
ρu′′2

a. Les équations gouvernantes adimensionnelles :

Le système d’équations scalaires gouvernant le problème, après avoir introduit les variables
adimensionnelles, s’écrit sous la forme suivante :

∂u ∂v
+ ∂y = 0 (2.5)
∂x

∂u ∂u ∂u ∂P ∂2 u ∂2 u
+ u ∂x + v ∂y = − ∂X + Pr ( ∂x2 + ∂x2 ) (2.6)
∂τ

∂v ∂v ∂v ∂P ∂2 v ∂2 v
+ u ∂x + v ∂y = − ∂y + Pr (∂x2 + ∂y2 ) + PrRaT (2.7)
∂τ

∂T ∂T ∂T ∂2 T ∂2 T
+ u ∂x + v ∂y = + ∂y2 (2.8)
∂t ∂x2

b. Les conditions aux limites adimensionnelles :

𝜕𝑇
Pour : x = 0 et y  u = v = 0 et =0
𝜕𝑥
𝜕𝑇
Pour : x = A et y  u = v = 0 et =0
𝜕𝑥
𝜕𝑇
Pour : y = 0 et x  u = v = 0 et = −1
𝜕𝑦
u0 H 𝜕𝑇
Pour : y = 1 et x  u = = Re Pr = Pe et v= 0 et 𝜕𝑦 = −1
α

c. Les nombres adimensionnels caractérisant le transfert de chaleur par

convection mixte :

Rapport de forme
Le rapport de forme de la cavité A, qui présente une caractéristique géométrique de la
cavité, est donnée par la relation :
𝐿
𝐴=
𝐻

Page | 22
 Nombre de Rayleigh
Le nombre de Rayleigh mesure l’importance relative de la poussée d’Archimède qui
provoque le mouvement du fluide par rapport aux forces visqueuses et à la conduction thermique
qui tendent à freiner ce mouvement,
gβqH 4
Ra =
ναλ
Nombre de Prandtl
Le nombre de Prandtl, dépend des caractéristiques du fluide, qui est le rapport de la viscosité
cinématique ν et la diffusivité thermique α, il caractérise l’importance relative des effets thermiques
et visqueux, il est défini par :
ν μ
Pr = =
α ρα
Nombre de Peclet
Nombre adimensionnel, Nombre adimensionnel, est le rapport du transfert par convection
forcée et par diffusion thermique, donné par :
u0 H
Pe = = Re Pr
α
4. Equations gouvernantes en formulation secondaire adimensionnelle :
La vorticité 𝜔 est le rotationnel du vecteur de vitesse du fluide, et la fonction de courant 𝜓
est une fonction qui représente les lignes de courant d’un fluide, correspondant aux trajectoires des
particules dans un écoulement stationnaire.
Les équations ont été formulées en termes de vorticité/fonction. Nous avons opté pour cette
formulation plutôt que la formulation en termes de variables primitives, pour les raisons suivantes :
 L’introduction de la vorticité et la fonction de courant, permet de satisfaire automatiquement
l’équation de continuité.
 La possibilité de l'élimination de la pression dans l’équation de conservation de mouvement.

Les relations mathématiques s’écrivent comme suit :

La fonction de courant et la vorticité sont définies par les relations suivantes :

∂v ∂u
Ω = ∂x − ∂y = −∇2 ψ (2.9)

∂ψ
u= (2.10)
∂y

∂ψ
v = − ∂x (2.11)

On introduit les relations ci-dessus dans les équations (2.6 - 2.8) et après plusieurs opérations
mathématiques elles deviennent :

Page | 23
 ∂Ω ∂(uΩ) ∂(vΩ) ∂2 Ω ∂2 Ω ∂T
+ + = Pr ( ∂x2 + ∂y2 ) + PrRa ∂x (2.12)
∂τ ∂x ∂y

∂T ∂T ∂T ∂2 T ∂2 T
+ u ∂x + v ∂y = + ∂y2 (2.13)
∂t ∂x2

Conditions aux limites adimensionnelles s’écrivent :

𝜕𝑇
Pour : x = 0 et y  u = v =  = 0 et =0
𝜕𝑥

𝜕𝑇
Pour : x = A et y  u = v =  = 0 et =0
𝜕𝑥

𝜕𝑇
Pour : y = 0 et x  u = v =  = 0 et = −1
𝜕𝑦

u0 H 𝜕𝑇
Pour : y = 1 et x  u = = Re Pr = Pe et v= 0 et  = 0 et 𝜕𝑦 = −1
α

5. Transferts thermique

Nombre de nusselt est un nombre adimensionnel utilisé dans les opérations de transfert thermique.
Il représente le rapport entre le transfert thermique total et le transfert par conduction.

Dans la présente étude on s’intéresse au calcul du taux de transfert de chaleur au sein du fluide à
une position x donnée, qui est exprimée à travers le nombres de Nusselt. Dans le cas où la paroi
inférieur de l’enceinte est soumis à un flux de chaleur constant ce nombre est défini par :

qH qH 1
Nu(x) = λ T′ = λ ΔT T = T(x,0)−T(x,1) (2.14)

T′ qH
Avec : T = ∆T et ∆T = , où ∆T= T (x, 0) - T (x, 1) représente la différence de température
λ

adimensionnelle locale entre les deux parois horizontales d’équations y=0 et y=1.

Le nombre de Nusselt moyen mesurant l’intensité du transfert de chaleur global dans la

direction verticale, s’écrit :

1 A
Nu = A ∫0 Nu(x)dx (2.15)

IV. Conclusion :

Après avoir développé les équations de conservation, de mouvement et de chaleur

ainsi que les conditions aux limites et initiale associées, la recherche des solutions, sera
conduite selon deux directions Une résolution par une voie numérique dans le chapitre suivant, qui
traitera les techniques utilisées pour simuler la convection et une approche analytique présentée au
quatrième chapitre, rendue possible en se basant sur l'approximation de l'écoulement parallèle.

Page | 24
 Chapitre III : Méthode numérique

I. Introduction

Dans le domaine de thermique, de la mécanique de fluides et de la combustion, les phénomènes

physiques sont décrits par des équations aux dérivées partielles (EDP) fortement couplées et non
linéaires. En général, Ces équations n’admettent pas des solutions analytiques sauf dans des cas très
simplifiés. Mais, une solution numérique peut être possible en transformant ces équations
différentielles en systèmes d'équations algébriques linéaires par une méthode de discrétisation avant
de résoudre ce système par des méthodes directes ou par itérations. Il existe plusieurs méthodes de
discrétisation sont utilisées actuellement on peut citer à titre d’exemple :

 La méthode des différences finis.

 La méthode des éléments finis.
 La méthodes des volumes finis

II. Méthode des différences finis

La méthode des différences finies est la plus ancienne des méthodes numériques. Elle a été
introduite au 18ème siècle par Euler, et s’avère la plus adaptée et la plus facile d’utilisation pour les
problèmes à géométrie simple, en plus d’assurer grande simplicité d’écriture et faible coût de calcul.

La variable (T, ψ, 𝜴) est décrite par plusieurs valeurs dans les points d’un maillage (Figure
2). Le développement de l’inconnue en séries de Taylor tronquées est utilisé pour approximer les
dérivées de l’inconnue, dans chaque point du maillage, par différences finies en utilisant les
inconnues voisines. En remplaçant les dérivées apparaissant dans le problème continu par des
différences divisées ou combinaisons de valeurs ponctuelles de la fonction en un nombre fini de
points discrets ou nœuds du maillage, on obtient un système d’équations algébriques pour les
valeurs de l’inconnue dans chaque point du maillage. L’un des principaux avantages de cette
méthode est d’assurer grande simplicité d’écriture et faible coût de calcul.

La méthode des différences finis, consiste à transformer un problème exact et continu en un

problème approché et discret. Comme chaque méthode numérique, elle comporte une phase de
maillage et une phase de discrétisation.

Page | 25
 Figure 2 : Maillage en (i,j) et repère cartésien en ( x,y )

La méthode des différences finies est donc, une technique de discrétisation pour la
résolution des équations qui s’écrivent sous la forme conservatrice comme suit :
Équation de fonction de courant :

∇2 Ψ = −𝛺 (3.1)

Équation de vorticité (2.12) :

∂Ω ∂T
⃗ Ω) = Pr∇2 Ω + PrRa ∂x
+ ∇(v (3.2)
∂τ

Équation d’énergie (2.13) :

∂T
⃗ T) = ∇2 T
+ ∇(v (3.3)
∂τ

III. Discrétisation des équations gouvernantes

Cette phase discrétise le problème sur le maillage, ce qui conduit à un problème approché
beaucoup moins difficile. Le principe de la méthode des différences finies, consiste qu’à l’intérieur
du domaine, les dérivées partielles spatiales, apparaissant dans les équations de base, sont
discrétisées selon un schéma centré aux différences finies. Ce choix est lié à la nature du problème
étudié, en effet les informations sont connues à l’amont et à l’aval (cavité fermée) et les termes
diffusifs se trouvent pondérés par un facteur supérieur ou égal à l’unité.

Page | 26
 1. Termes transitoires

En utilisant les différences finies, la dérivée dans le temps est approximée par une
discrétisation implicite du premier ordre. Ainsi, en considérant les valeurs de f entre les instants τ𝑛 =
𝑛∆τ et τ𝑛+1 = (𝑛 + 1)∆τ, on obtient :

∂f n+1 fn+1 n
i,j −fi,j
| = (3.4)
∂τ i,j ∆τ

2. Termes convectifs

Les termes convectifs des équations (2.12) - (2.13) écrites sous une forme conservative sont
approchés par :

∂uf ui+1,j fi+1,j −ui−1,j fi−1,j

| = (3.5)
∂x i,j 2∆x

ou f représente la température  ou la fonction de courant Ψ et ∆x, l’espacement entre deux

nœuds voisins.

3. Termes diffusifs et autres que diffusifs

A l’intérieur du domaine, les dérivées première, seconde, apparaissant dans les équations
(2.12) et (2.13), sont approximées par :

∂f fi+1,j −fi−1,j
| = (3.6)
∂x i,j 2∆x

∂2 f fi+1,j −2fi,j +fi−1,j

| = (3.7)
∂x2 i,j ∆x2

∂2 f
Une approximation de la dérivée croisée ∂x ∂y de la fonction de 2 variables f(x,y).La discrétisation

du domaine de calcul est bidimensionnelle et fait intervenir deux pas d’espace supposés constants
∆x et ∆y dans les directions x et y.

∂2 f fi+1,j+1 −fi+1,j−1 −fi−1,j+1 +fi−1,j−1

| = (3.8)
∂x ∂y i,j 4∆x∆y

La discrétisation des dérivées partielles exprimant les conditions aux limites se fait à l’aide d’un
schéma aux différences finies décentrées en amont :

∂f fi−2,j −4fi−1,j +3fi,j ∂f fi,j−2 −4fi,j−1 +3fi,j

| =
∂x i,j 2∆x
| =
∂y i,j 2∆y
(3.9)

Page | 27
 Ou en aval

∂f 4fi+1,j −fi+2,j −3fi,j ∂f 4fi,j+1 −fi,j+2 −3fi,j

| = | = (3.10)
∂x i,j 2∆x ∂y i,j 2∆y

IV. Résolution des équations gouvernantes

1. Méthode des directions alternées (ADI)

La méthode ADI (Alterning Direction Implicite) est applicable en deux étapes : pour
la première étape (prédicteur), elle est implicite dans la direction x et explicite dans la
direction y. Pour la seconde étape (correcteur), la méthode ADI est explicite dans la direction
x, et implicite dans la direction y.

Pendant la première étape les valeurs aux nœuds (i, j) sont connues à l’ instant (n) mais
inconnues à l’instant ( nt + 1/2 )( tnt+1/2 = tnt + t/2 ) , Les valeurs inconnues sont associées à la
direction x seulement (C’est à dire j fixe). Les valeurs `a l’instant inconnu (n +1/2) dans tous les
nœuds de maillage sont données par la résolution des systèmes d’équation pour i = 2, . . . , imax - 1 et
pour chaque j tel que j = 2, . . . , jmax - 1.

Dans la second étape les valeurs aux nœuds (i, j) sont inconnues à l’ instant (n+1), mais
connues à l’instant ( nt + 1/2 )( tnt+1 = tnt+1/2 + t/2 ), On obtient les valeurs de tous les nœuds à
l’instant inconnu (n + 1) par la résolution des systèmes d’équations pour j = 2, . . . ,jmax – 1 et pour
chaque i tel que i = 2, . . . , imax - 1 . Les valeurs des nœuds à l’ instant (n + 1) sont associées à la
direction y seulement (c’est-`a-dire i fixe).
Cette méthode, fréquemment utilisée pour les fluides newtoniens. Généralement, elle donne lieu à
des matrices tridiagonales dans les deux directions et le schéma de calcul est très efficace car il est
simple, stable et consistant.

Figure 3: illustration du maillage de la méthode ADI

Page | 28
 2. Equations de la vorticité, de l’énergie

La résolution numérique des équations de la verticité de l’énergie et de la concentration est

effectuée par la méthode implicite des directions alternées (méthode ADI). Puisque les équations de
de la verticité et de l’énergie ( , ) présentent des formes similaires, elles peuvent être écrites sous la
forme unique d’une équation de transport générale d’une quantité quelconque f dans un fluide en
mouvement :

D’autre part, puisque les équations (3.2) et (3.3) présentent des formes similaires,
elles peuvent être condensées sous la forme unique suivante :

∂f
⃗ f) = Λ∇2 f + Sf
+ ∇(V (3.11)
∂τ

f : La fonction à déterminer (𝜴 ou ) ;
𝜕f
: Le terme transitoire
𝜕𝜏

⃗ f) : Le terme convectif
∇(𝑉⃗
∇2 f : Le terme diffusif
Λ : Le coefficient de diffusion ;
𝑆𝑓 : Le terme source.

Tous ces termes sont représentés dans le tableau suivant :

Equation f Terme Terme Terme Terme Γ

transitoire convectif diffusif source 𝑆f
Vorticité 𝜴 𝜕Ω ∇(v
⃗ Ω) ∇2 Ω 𝑃𝑟𝑅𝑎 𝜕𝑥
𝜕𝑇 Pr
𝜕𝜏

Energie T 𝜕T ∇(v
⃗ T) ∇2 T 0 1
𝜕𝜏

a. Schéma implicite en x et explicite en y

On utilise une discrétisation implicite dans la direction x , en optant pour des différences
centrées au temps tn+1/2 pour les dérivées en x et au temps tn pour les dérivées en y, la discrétisation
des dérivées partielles est faite selon :

n+1/2
fi,j −fn
i,j ui+1,j fn+1/2 i+1,j −ui−1,j fn+1/2 i−1,j vi,j+1 fn i,j+1 −vi,j−1 fn i,j−1
+ + =
∆τ/2 2∆x 2∆y

fn+1/2 i+1,j −2fn+1/2 i,j +fn+1/2 i−1,j fn i,j+1 −2fn i,j +fn i,j−1
Λ i,j ( + ) + Sf ni,j (3.12)
∆x2 ∆y2

Page | 29
 𝑛 𝑛+1/2
Où 𝑓𝑖,𝑗 représente la valeur discrète de la fonction f à l’instant τ𝑛 = 𝑛∆τ et 𝑓𝑖,𝑗 la valeur discrète
𝑛+1
intermédiaire qui servira pour la détermination de 𝑓𝑖,𝑗 à l’instant τ𝑛+1 = τ𝑛+∆τ .Ainsi, après
réarrangement on obtient la forme discrète implicite en X sous la forme matricielle suivante :

n+1/2 n+1/2 n+1/2

Axi,j fi−1,j + Bi,j
x
fi,j x
+ Ci,j fi+1,j = Dxi,j (3.13)

Avec :

u Λ
Axi,j = ( 2∆x
i+1,j
− ∆xi,j2 ) (3.14)

x 2 2Λi,j
Bi,j =( + ) (3.15)
∆τ ∆x2

x i−1,j u Λ
Ci,j = − ( 2∆x + ∆xi,j2 ) (3.16)

vi,j+1 Λ 2 2Λi,j v Λ
Dxi,j = (− + ∆yi,j2 ) f n i,j+1 + (∆τ − i,j−1
) f n i,j + ( 2∆x + ∆yi,j2 ) f n i,j−1 + Sf ni,j (3.17)
2∆y ∆y2

b. Schéma implicite en y et explicite en x

On utilise une discrétisation implicite dans la direction y , en optant pour des différences
centrées au temps tn+1/2 pour les dérivées en x et au temps tn+1 pour les dérivées en y, la
discrétisation des dérivées partielles est faite selon :

n+1/2
fn+1
i,j −fi,j ui+1,j fn+1/2 i+1,j −ui−1,j fn+1/2 i−1,j vi,j+1 fn+1 i,j+1 −vi,j−1 fn+1 i,j−1
+ + =
∆τ/2 2∆x 2∆y

fn+1/2 i+1,j −2fn+1/2 i,j +fn+1/2 i−1,j fn+1 i,j+1 −2fn+1 i,j +fn+1 i,j−1
Λ i,j ( + ) + Sfni,j (3.18)
∆x2 ∆y2

n+1/2
Où fi,j et f n+1 i,j désignent les valeurs de la fonction f aux instants respectifs tn+1/2 et tn+1.
De la même manière que précédemment, l´équation ci-dessus peut se mettre sous la
forme suivante :

y n+1 y y y
Ai,j fi,j−1 + Bi,j fi,jn+1 + Ci,j fi,j+1
n+1
= Di,j (3.19)

Avec :

y i,j+1v Λ
Ai,j = ( 2∆y − ∆yi,j2 ) (3.20)

y 2 2Λi,j
Bi,j = (∆τ + ∆y2
) (3.21)

Page | 30
 y v
i,j−1 Λ
Ci,j = − ( 2∆y + ∆yi,j2 ) (3.22)

y ui+1,j Λ 2 2Λi,j u Λ
+ ∆xi,j2 ) f n+1/2 i+1,j + (∆τ − + ∆xi,j2 ) f n+1/2 i−1,j + Sf ni,j
i−1,j
Di,j = (− ) f n+1/2 i,j + ( 2∆x (3.23)
2∆x ∆x2

En tenant compte de l’ensemble des nœuds on obtient un système d’équation dont la solution
se ramène à L’inversion d’une matrice tridiagonale. Le champ de température obtenu pour
l’ensemble des nœuds du domaine sera utilisé pour déterminer le champ de fonction de courant.

3. Équation de Poisson pour la fonction de courant

L’équation de la fonction de courant, une fois discrétisée avec un schéma centré classique, est
résolue avec la méthode de sur-relaxation successive par point (PSOR) qui n’est autre que la
méthode de Gauss-Seidel relaxée. La valeur de la fonction de courant aux nœuds est alors
calculée à l’aide de la relation suivante :

n+1
ψn+1
i,j = (1 − K)ψni,j + Kψ′ i,j (3.24)

Avec :

n+1 1
ψ′ i,j = 2(1+η2 ) (ψni+1,j + ψni−1,j + η2 ψni,j+1 + η2 ψni,j−1 + Δx 2 Ωi,j (3.25)

Alors :

K
ψn+1
i,j = (1 − K)ψni,j + 2(1+η2 ) (ψni+1,j + ψni−1,j + η2 ψni,j+1 + η2 ψni,j−1 + Δx 2 Ωi,j (3.26)

∆x
D’où : η = ∆y et 𝐾 est le coefficient de sur-relaxation, dont la valeur optimale est donnée par la

relation de Frankel (1950) :

π π
|1−√1−ζ2 | cos( )+η 2 cos( )
mx ny
K opt = 2 et ζ= (3.27)
ζ2 ζ2

Dans l’expression (3.16), valable uniquement pour un maillage rectangulaire uniforme, 𝑚x et

ny et désignent respectivement les nombres d’intervalles dans les directions x et y. A noter ici que
pour K= 1 , la procédure est identique à celle de Gauss-Seidel.

a. Champs de vitesse

Les champs des vitesses sont évalués à partir de la définition de la fonction de courant, puis
on utilisant les différences centrées pour obtenir une expression discrétisée de ces composantes, ce
qui nous donne :

Page | 31
 ψn+1 n+1
i,j+1 −ψi,j−1
un+1
i,j = (3.28)
2∆y

n+1 ψn+1 n+1

i−1,j −ψi+1,j
vi,j = (3.29)
2∆x

b. Conditions aux limites

Les conditions aux limites sont discrétisées suivant un schéma aux différences finies
décentré avant ou arrière selon la paroi considérée.

- Paroi horizontale inférieure (y = 0) :

4Ti,2 −Ti,3 +2Δy

ui,1 = vi,1 = ψi,1 = Ti,1 − =0
3

- Paroi horizontale supérieure (y= 1) :

4Ti,n −Ti,n−1 +2Δy

ui,n+1 − Pe = vi,n+1 = ψi,n+1 = Ti,n+1 − =0
3

- Paroi verticale gauche (x = 0) :

4T2,j −T3,j
u1,i = v1,i = ψ1,i = T1,i − =0
3

- Paroi verticale droite ( x = A ) :

4Tm,j −Tm−1,j
um+1,i = vm+1,i = ψm+1,i = Tm+1,i − =0
3

c. Critère de convergence

A chaque pas de temps, Δt, qui varie dans cette étude entre 10-7 et 10-4 dépendamment des
paramètres gouvernants, la fonction de courant, ψ , estcalculée à partir du champ de vorticité et
comparée à celle obtenue à l’itérationprécédente. Ce processus est répété jusqu’à ce que le critère
de convergence

∑i,j|ψiter+1
i,j −ψiter
i,j |
< εC (3.30)
∑i,j|ψiter+1
i,j |

soit satisfait. A signaler que le paramètre 𝜀𝐶 , dont le choix n’est pas libre, doit être
suffisamment petit de sorte que de grandes erreurs de troncature ne puissent s’introduire, tout
en restant au-dessus des erreurs d’arrondis(10−5 < 𝜀𝐶 < 10−4 )

Page | 32
 d. Critère de stationnarité

Pour le problème traité, le régime permanent est atteint lorsque le critère

|fn+1 n
i,j −fi,j |
< εp
|fn+1
i,j |

(3.31) Où f représente les valeurs des quatre champs du système ( , 𝜓 et T ) prise en chaque point du
domaine à l’instant 𝑛∆t. 𝜀𝑃 est une constante (comprise entre 10-3 et 10-4)

V. Conclusion

Dans ce chapitre nous avons présenté la méthode numérique des différences finies pour la
résolution du modèle mathématique. Dans le prochain chapitre, nous présentons la méthode
analytique en se basant sur le concept de l’écoulement parallèle.

Page | 33
 Chapitre IV : Etude analytique

I. Introduction :

Les problèmes physiques rencontrés dans notre quotidien, sont décrits par des équations aux
dérivées partielles, fortement couplées et non linéaires. En général, Ces équations n’admettent pas
de solutions analytiques exactes, sauf dans des cas simplifiés.

L’écoulement convectif dans le système représenté sur la Figure 1 est aussi régi par les
équations de base (2.9)-(2.13), qui sont des équations aux dérivées partielles, non linéaires et
fortement couplées. La solution analytique n'est possible, que dans certains cas limites pour lesquels
les équations se simplifient considérablement. Ainsi, pour le cas des cavités élancées, A ≫1,
l'obtention d'une solution analYtique approximée est possible, en utilisant le concept de
l’écoulement parallèle. Cette méthode, permet une compréhension rapide des phénomènes
physiques, impliqués dans la convection naturelle, grâce à sa flexibilité et à sa manipulation facile

II. Approximation d’écoulement parallèle :

Plusieurs auteurs se sont intéressés au concept de l’écoulement parallèle, pour l’étude de la

convection naturelle dans des cavités élancées ou aplaties. Ils ont montré numériquement que, cette
approximation décrit d’une façon précise l’écoulement et le transfert de chaleur.

Le concept de l’écoulement parallèle consiste à supposer que dans une cavité de grande
extension, l’écoulement engendré peut être décomposé en trois zones. Ces dernières consistent en
deux régions à l’extrémité de la cavité dans lesquelles l’écoulement tourne à 180◦ et la région
centrale de la couche dans laquelle l’écoulement est parallèle relativement aux longues parois de la
cavité. ceci permet de négliger la composante de la vitesse perpendiculaire aux parois horizontales,
de telle sorte que :

u(x, y) = u(x) et v(x, y) = 0 (4.1)

Dans ce type de problèmes il est courant d’éliminer la pression en utilisant la formulation de la

fonction de courant,Ψ. En utilisant le concept de l’écoulement parallèle la fonction de courant Ψ
dépende uniquement de l’ordonnée y :

Page | 34
 ψ(x, y) = ψ(y) (4.2)

À partir des équations de base et des bilans de conservation de l’énergie dans le volume de
contrôle, les profils de la température est alors donné par la somme d’un terme définissant une
variation longitudinale linéaire et d’un autre terme donnant la distribution transversale :

A
T(x, y) = C (x − 2 ) + θ (4.3)

Où 𝐶𝑇 est une constante qui exprime le gradient de température selon la direction (x). f est une
fonction décrit la variation de température suivant la direction y.

III. Equations gouvernantes approximées

Comme on l’a vu précédemment, la fonction de courant ne dépend que de y, la température est

une fonction de x et y. de plus, on se place en régime permanent. en substituant les approximations
(4.1)-(4.3) dans les équations (2.9)-(2.13), on obtient l’équation suivante :

d3 u
= Ra C (4.4)
dy3

d2 θ
= Cu(y) (4.5)
dy2

A partir de l’expression de  ( x,y ) adoptées pour l’écoulement parallèles :

dT dT dθ dψ
=C ; = dy ;u= (4.6)
dx dy dy

d4 ψ
L’équation précédente devient : = Ra C (4.7)
dy4

En intégrant l’équation de l’énergie (2.13) par rapport à x, tout en respectant les conditions aux
limites associées à la fonction de courant et à la température ,on obtient :

dθ
= Cψ − 1 (4.8)
dy

IV. Solution analytique

Dans cette partie, nous présentons la solution analytique stationnaire du problème pour un
fluide newtonien.

1. Résolution de l’équation de conservation de quantité de mouvement :

On commence par la résolution de l’équation de vitesse u(y) qu’on a obtenue précédemment

après l’application de l’approximation de l’écoulement parallèle :

Page | 35
 d3 u
= Ra C (4.9) En
dy3

intégrant 3 fois excessives l’équation (4.9) suivant y , on a trouvé l’expression de vitesse u(y) sous
forme d’une polynôme :

1 1
u(y) = 6 Ra C y 3 + 2 c1 y 2 + c2 y (4.10) A

l’aide des conditions aux limites et la condition d’écoulement de retour (conservation de la

matière):

u (0) = 0 et u (1) = Pe (4.11)

1
∫0 u(y)dy = 0 (4.12)

1 1
Les constantes C1 et C2 sont : c2 = −2 Pe + 12 Ra C et c1 = 6 Pe − 2 Ra C

Donc la distribution de vitesse suivant y est :

Ra C
u(y) = (2 y 3 − 3 y 2 + y) + Pe(3 y 2 − 2y) (4.13)
12

Détermination de la fonction de courant  :

𝑑𝜓
On a : 𝑢(𝑦) = 𝑑𝑦

Apres l’intégration suivant y , tout en respectant les conditions aux limites de la fonction de
courant , on obtient :

Ra C
ψ(y) = ( y 4 − 2 y 3 + y 2 ) + Pe(y 3 − y 2 ) (4.14)
24

2. Résolution de l’équation de l’énergie :

L’équation de l’énergie (4.5), après le remplacement de u(y) par son expression trouvée, est :

d2 θ Ra C2
= (2 y 3 − 3 y 2 + y) + Pe C (3 y 2 − 2y) (4.15)
dy2 12

On a les conditions suivantes :

∂θ
) = −1 (conditions aux limites) (4.16)
∂y (0,1)

1
∫0 θ(y)dy = 0 (condition de l’écoulement de retour ) (4.17)

Apres l’intégral de l’équation d’énergie suivant y , tout en respectant les conditions aux limites et de
l’écoulement de retour associées à la température , on obtient :

Page | 36
 Ra C2 Pe C 1
θ(y) = (12 y 5 − 30 y 4 + 20 y 3 − 1) + (15 y 4 − 20 y 3 + 2) − y + (4.18)
1440 60 2

3. Détermination de la constante C

La constante C représente le gradient de température. Elle est déterminée en intégrant l’équation

(2.13) sur le volume de contrôle choisi arbitrairement dans la figure 1, en prenant en compte les
conditions aux limites (4.16) et (4.17) et en raccordant avec la région de l’écoulement parallèle :

1
C = ∫0 u(y)θ(y)dy (4.19) C =
1 dψ
∫0 θ(y)dy (4.20)
dy

Pour trouver l’expression de C, on intègre par partie la relation précédente :

1 dψ 1 dθ(y)
C = ∫0 θ(y)dy = [ψ(y)θ(y)]10 − ∫0 ψ(y)dy
dy dy

Compte tenus de la relation (4.8) et des conditions aux limites (4.16) et (4.17)

1 1 1
C = − ∫0 ψ(y)(Cψ(y) − 1)dy = −C ∫0 ψ2 (y)dy + ∫0 ψ(y)dy

On trouve l’expression de C sous la forme :

1
∫0 ψ(y)dy
C= 1 (4.21)
1+∫0 ψ2 (y)dy

Apres intégration de la fonction de courant et de son carré , on a obtenu une équation linéaire de
3eme ordre :

Ra2 Ra Pe Pe2 Ra Pe
(576∗630) C3 − (12∗280) C2 + (1 + 105 − 720) C + 12 = 0 (4.22)

4. Détermination de nombre nusselt :

On a le terme de nusselt que on a déjà définit est :

1
Nu(x) = T(x,0)−T(x,1)

A
Avec : T(x, y) = C (x − 2 ) + θ(y)

1
Alor : Nu = θ(0) − θ(1)

après le remplacement de l’expression de (y) , on obtient :

720
Nu = −Ra C2 +60 Pe C +720 (4.23)

Page | 37
 Chapitre V : Résultat Et Discussions

I. Introduction

Dans ce chapitre nous allons nous intéresser à l’étude analytique de la convection mixte dans
la géométrie confinant un fluide Newtonien ayant Pr = 7. Nous verrons en premier lieu, l’influence
des paramètres de contrôle à savoir.; le nombre de Peclet Pe, et le nombre de Rayleigh thermique
Ra, sur l’écoulement et sur le transfert de chaleur. Dans la limite d’une cavité allongée (A ≫ 1), la
solution basée sur l’approximation de l’écoulement parallèle, développée dans le chapitre IV.

II. Effet du Maillage

En vue d’aboutir à un maillage optimal assurant le bon accord entre les résultats numériques
et analytiques, plusieurs essais numériques sont effectués en utilisant un code de calcul basé sur la
méthode des différences finis élaboré en langage Fortran, pour différents maillages et une valeur de
Pe = 50 , 𝑅𝑎 = 104 .Le maillage (300×80) a été approprié pour modéliser le champ d’écoulement et
de température dans une cavité allongée , puis qu’il fournit des résultats quasiment en accord avec
ceux analytiques avec une erreur ne dépasse pas 0,5 % .

Les résultats comparatifs concernant l’effet du maillage sur les résultats numériques, sont
présentés dans le tableau suivant :

Tableau 2: effet du maillage sur les resulats nimériques

Maillage Nusselt moyen Nu c

Résultat Analytique - 3 ,704 -9 ,4784
Résultat Numérique 200 x 80 3,6414 -9 ,4657
240 x 80 3,66818 -9,47288
260 x 80 3,674 -9,4741
280 x 80 3,684 -9,4769
300x80 3,689427 -9,4783

III. Effet de rapport de forme :

Pour estimer la valeur minimale de rapport de forme A pour laquelle l’écoulement soit
parallèle dans une couche horizontale, des tests numériques ont été effectués. Des résultats typiques
sont rapportés dans la Figure 4, pour différentes valeurs de rapport de forme A, ou les résultats
numériques sont indiqués par des symboles en noires et la solution analytique approximée est
indiquée par une ligne, valide pour le cas d’une cavité verticale infini (A>>1). A partir ce figure, on

Page | 38
 peut estimer la valeur minimum du rapport de forme A, au-dessus de laquelle l’écoulement dans la
cavité peut être considéré comme étant parallèle. On constate donc qu’en augmentant le rapport de
forme de A=2 à A=24, les résultats indiquent que le nombre de Nusselt. Le fait d’augmenter le
rapport de forme A pour (𝑅𝑎 = 104, 𝑃𝑒 = 50), conduit donc à une situation asymptotique pour
laquelle les solutions sont indépendantes du facteur de forme A. Alors, ceci est en bon accord avec
la prédiction de concept de l’écoulement parallèle. Suit à plusieurs tests, nous pouvons donc
conclure que les résultats numériques peuvent être considérés indépendants du rapport de forme de
cavité lorsque A ≥ 22. Pour cette raison la majorité des résultats rapportés dans ce travail sont
obtenus pour une valeur de rapport de forme A=24.

̅̅̅̅
𝑁𝑢

Figure 4 : Influence du rapport de forme sur le taux de transfert de chaleur

IV. Profils des lignes de courant, isothermes

La figure (5 et 6) présente les lignes de courant et les isothermes.

A partir des résultats du logiciel teplot du traçage , . On remarque que les isothermes sont
stratifiés et les lignes de courants gardent un aspect parallèle et présentent une symétrie par rapport
au centre de la cavité ce qui prouve l’existence de la solution d’écoulement parallèle.

Figure 5 : Isothermes

Page | 39
 Figure 6 : Lignes de courant

V. Profil de vitesse, de fonction de courant et la température le long de la section verticale

de la cavité :

Les variations de la vitesse ,la fonction de courant Ψ et la température au centre de la cavité

en fonction de 𝑦 pour Pe=50 sont exposées dans les figure (7-8-9) ,l’analyse de ces figures montre
l’accord entre la solution numérique et analytique.

1. Profil de vitesse :

Les profils de vitesse présente simultanément deux extrêmes (maximal et nulle) près des
frontières horizontales.

On aperçoit une valeur maximale au voisinage de la paroi active ce qui est dû aux conditions
cinématiques imposées à la paroi horizontal supérieure et une valeur de zéro au voisinage de la
paroi immobile inférieure.

Nous constatons ses profils augmentent les points extrémales de u lorsque l’augmentation de
Ra c’est à dire la convection naturel est dominante inversement forcée domine lorsqu’on diminue le
nombre de Rayleigh.

Figure 7: profil de la vitesse U au long de la section verticale de la cavité avec Pe=50 et Ra varie

2. Température :

On aperçoit que la valeur maximale de la température est diminué lorsqu’on augmente la

fraction volumique. Et ceci du au flux thermique exercé sur la paroi inférieure.

Page | 40
 On constate que la température commence par diminuer en s’éloignant de la paroi inférieure
jusqu’à arriver à la température la plus froide de la paroi supérieure ,ceci est dû au fait que
l’écoulement de fluide dans le sens horaire permet de transporter la chaleur de la paroi inférieure
chaude et porter le fluide vers le haut après l’évacuation de la chaleur à travers la paroi supérieure.
,y )
A
2
T(

Figure 8 : profil de température au long de la section de la cavité avec Pe = 50 et Ra varie

3. Fonction de courant

On constate que l’augmentation de Ra résulte une augmentation progressive de la fonction

de courant vers des valeurs négatives.

Figure 9 : Fonction de courant au long de la section de la cavité avec Pe =50 et Ra varie

Page | 41
VI. Influence des paramètres de contrôles sur le transfert de chaleur

1. Effet du nombre de Peclet sur Nusselt :

La figure 10 présente l’´évolution de nusselt moyen en fonction du nombre de Peclet Pe

dans la plage de 0 ≤ Pe ≤ 103 pour Ra = 720 et Ra=5000 et Ra =10000.

On remarque l’existence de deux types de solutions. Des solutions stables et d’autre instable
tracée par des points et ces dernières se disparait dés Pe critique.

On aperçoit que pour des grandes valeurs de Peclet , les courbes se confond entre eux et
avec la courbe de convection forcé, donc on déduit que la convection forcée est dominante.

̅̅̅̅
Nu Ra =10000

Ra = 5000

Stable

.…… Instable
Ra = 720

convection forcée pure

Figure 10: nusselt en fonction de peclet

2. Effet du nombre de Rayleigh sur Nusselt :

L’évolution du nombre de nusselt Nu en fonction de Ra est illustrée dans la figure 11, on

remarque l’apparition de deux solutions après Ra critique , une stable augmente progressivement
avec Ra et l’autre instable tracée par des points est l’inverse, elle diminue quand Ra augment.

pour des grandes valeurs de Ra., on remarque que les courbes de convection mixte se confond entre
eux et avec la courbe de convection naturelle donc la convection naturelle qui va être dominante.

Page | 42
 ̅̅̅̅
Nu Pe = 50

Pe = 20

Stable

…… Instable
Pe = 20

Pe = 0
Pe = 50
convection Natuelle pure

Figure 11 : Nusselt en fonction de Rayleigh

VII. Influence des paramètres de contrôles sur l’intensité de l’écoulement au centre de la

cavité :

1. Effet de nombre de Peclet sur l’intensité de l’écoulement au centre de la cavité :

L’évolution de l’intensité d’écoulement | Ψc | en fonction de Pe est illustrée dans la figure

12, on aperçoit deux types des solutions , des solutions sont stables et d’autre tracées par des tirés
sont instables et ces derniers sont disparait au-dessus de Pe critique .

Plus le nombre de Peclet augment , la valeur absolu de Ψc s’approche de l’asymptote de la

convection forcée , ce qui augmente le % de la convection focée dans la convection mixte.

Page | 43
 c

Ra = 10000

Ra = 5000

Convection forcée pure

Ra = 720

Stable

…… Instable

Figure 12: c en fonction de peclet

2. Effet du nombre de Rayleigh sur 𝚿c :

La figure 13 présente l’´évolution de la valeur absolu de l’intensité de l’´écoulement au centre de la

cavité en fonction du nombre de Rayleigh Ra dans la plage de 0 ≤ Ra ≤ 108 pour Pe= 0 et Pe= 20.
On remarque l’existence de deux solutions à partir de Rac = 6485,3251 . Une, étant stable ou la
valeur absolu de l’intensité d’écoulement s’accroit avec Ra et l’autre instable représentée par des
points dans laquelle la valeur absolu de l’intensité d’écoulement diminue.

c Stable

……. Instable

Pe = 20

Pe = 20
Pe = 0

convection Natuelle pure

Figure 13: c en fonction de Rayleigh

Page | 44
 Conclusion générale :

Dans ce travail, nous avons réalisé une étude analytique et numérique de la convection mixte
dans une cavité rectangulaire allongée, remplie par un fluide Newtonien. Les parois verticales de
l’enceinte sont immobiles et adiabatiques tandis que les autres parois , on a le supérieure mobile et
l’inférieure chauffée par une densité de flux constante. L’approximation de Boussinesq est utilisée
pour établir le modèle mathématique décrivant le phénomène de la convection mixte au sein du
milieu fluide. Dans le but de résoudre analytiquement les équations adimensionnelles régissantes,
écrites en formulation de variables secondaires. Ce problème est caractérisé par deux paramètres de
contrôle, ce sont le nombre de Rayleigh Ra et le nombre Peclet Pe. Par ailleurs, dans la limite d’une
cavité allongée (A ≫ 1), une solution analytique, basée sur l’approximation de l’écoulement
parallèle, a été déterminée pour prédire qualitativement les caractéristiques de l’écoulement et du
transfert thermique. Dans l’étude némerique on a suivi la méthode de différences finis.

Les principales conclusions de la présente étude peuvent être résumées comme suit :

 la convection mixte dans une cavité allongée et remplit par un fluide newtonien est
essentiellement contrôlée par les paramètres Ra et Pe.
 Il y a un accord entre les résultats analytiques et numériques.
 l’augmentation des deux paramètres de contrôle Rayleigh Ra et Peclet Pe de favorise
et améliore le transfert thermique et l’intensité d’écoulement.
 Pour rester dans la zone de la convection mixte , il ne faut pas augmanté Ra ou bien
Pe à des grandes valeurs malgré les bonnes résultats de son accroissement. Car on a
chaque type de la convection ( naturelle , forcée , mixte ) a sa zone de définition qui
se caractérise par des valeurs de Ra et Pe bien définit et spécifique.

Dans la présente étude, la cavité rectangulaire disposée horizontalement, n’est qu’une

configuration géométrique parmi d’autres, donc pour approfondir les connaissances de nombreux
points nécessitent encore une étude plus approfondie, en l’occurrence :

 Reprendre la même étude pour un fluide non Newtonien.

 Etudié d’autres conditions aux limites.
 Examiner l’effet de la convection thermosolutale.
 Examiner l’effet de l’inclinaison de la cavité.

Ceci permettra de voir de près le rôle de l’effet de ces paramètres de contrôle combinés avec celui
des conditions aux limites, sur la convection mixte.

Page | 45
 Bibliographie :

[1]. M.M. Ganzarolli, L.F. Milanez. "Natural convection in rectangular enclosures heated from
below and symmetrically cooled from the sides". Int. J. Heat Mass Transfer, 38, 1063–1073.(1995).
[2]. Chu H.H.-S., Churchill S.W. et Patterson C.V.S. "The effects of heater size, location, aspect
ratio, and boundary conditions on two-dimensional, laminar, natural convection channels''. J. Heat
Transfer, Vol. 98, pp. 1194-201. (1976).
[3]. Lamsaadi. M, Naïmi. M & Hasnaoui. M, (2005) Natural Convection of Non-Newtonian
Power Law Fluids in a Shallow Horizontal Rectangular Cavity Uniformly Heated From Below,
Heat and Mass Transfer, Vol. 41, N° 3, pp. 239-249,.
[4]. Lamsaadi. M, Naïmi. M & Hasnaoui. M, (2005) Natural Convection of Non-Newtonian
Power Law Fluids in a Shallow Horizontal Rectangular Enclosure Uniformly Heated From the Side,
Soumise pour publication dans Journal of Energy Conversion and Management, sous presse.
[5]. F. ZOUIRI, M.A. OULMANE, N. LABSI, Y.K. BENKAHLA, A. BOUTRA, "Convection
naturelle au sein d’une cavité carrée munie d’une source chauffante placée sur sa paroi inférieure",
23ème Congrès Français de Mécanique, Lille Septembre 2017.
[6]. El Alami SEMMA, Abdelkhalek CHEDDADI, Antonio CAMPO, Hassan RIDOUANE,
"Convection Etude numérique de la convection naturelle dans une cavité triangulaire semi -
poreuse ",
[7]. R. Zarrit, M.S. Boumaza, S. Kherrour et B. Dadda, "Convection naturelle dans une cavité
rectangulaire inclinée de différents rapports de forme " Revue des Energies Renouvelables Vol. 19
N°1 (2016) 97 – 109.
[8] M. Mamou, P. Vasseur et E. Bigen, Analytical and numerical study of double diffusive
convection in a vertical enclosure, Heat and Mass Transfer, vol. 32, pp. 115-125, 1996.
[9]. I.E. Sarris, I. Lekakis, N.S. Vlachos. "Natural convection in rectangular tanks heated
locally from below". Int. J. Heat and Mass Transfer, 47, 3549–3563. (2004)

[10] ALLOUI Z., DUFAUA L., BEJI H., VASSEUR P., Multiple steady states in a porous
enclosure partially heated and fully salted from below Int. J. of Therm. Sci. vol. 48, pp.
521–534, 2009.
[11] S. B. Sathe et T. W. Tong, Measurements of natural convection in partially porous
rectangular enclosures of aspect ratio 5. Int. Comm. Heat Mass Transfer, Vol. 15, pp. 203-212
(1988).
[12] C. Revnic, T. Grosan, I. Pop et D.B. Ingham, (2009) «Free Convection in a Square Cavity
Filled with a Bidisperse Porous Medium » Int. Journal of Thermal Sciences, 48, 1876–1883
[13] D.A. Nield, A.V. Kuznetsov, (2008) « Natural Convection about a Vertical Plate

Page | 46
Embedded in a Bidisperse Porous Medium » Int. J. Heat Mass Transfer 51, 1658–1664
[14] C. W Horton et F. T Rogers, (1945) « Convection Curent in a Porous Medium» Journal of
Applied Physics, 16, 367-370
[15] V. Prassad et A. Tuntomo, Inertia effects on natural convection in a vertical porous cavity.
Num. Heat Transfer, Vol. 11, pp. 295-320 (1987).
[16] Abdellatif LOUARAYCHI, Mohamed LAMSAADI, Mohamed NAÏMI, Hassan EL
HARFI,Mourad KADDIRI , Étude de la convection mixte dans une cavité rectangulaire
horizontale chauffée latéralement. Nature & Technology Journal. Vol. A : Fundamentals and
Engineering Sciences, 17 (2017) 12-16.
[17]. S. BOUABDALLAH, A. MEDJELLED, A. BENCHATTI & A. HAMINI, "Convection
Mixte dans une Cavité Carrée Contenant des Portes d’entrée et de Sortie d’air " RSSI, Vol. 01, No.
02, Décembre (2011) 01-09.
[18] Abdelsalam Bouchair, Omar Mezghanni" étude de la convection à l'intérieur d'une cavité
remplie par un fluide conducteur et soumise à un champ magnétique externe", Journal of
Constructional Steel Research, vol. 83, pp. 51-61, 2013.
[19]. Joo – Sik, Y. (1998): Mixed convection of air between two horizontal concentric cylinders
with a cooled rotating outer cylinder, International Journal of Heat and Mass Transfer, vol. 41,
N°2, pp.293 – 302.
[20] Ali J. Chamkha, " hydromagnetic combined convection flow in a vertical lid-driven cavity with
internal heat generation or absorption" , Numerical Heat Transfer, Part A, 41, pp. 529-546,2002

[21]. Louaraychi, A., Lamsaadi, M., Naïmi, M., El Harfi, H., Kaddiri, M., Raji, A., & Hasnaoui, M.
(2019). Mixed convection heat transfer correlations in shallow rectangular cavities with single and
double-lid driven boundaries. International Journal of Heat and Mass Transfer, 132,394–406.
[22]A. ZAHIRI, M.LAMSAADI, H.EL HARFI, M.KADDIRI, M.NAÏMI : " L’effet de champ
magnétique sur la convection mixte dans une cavité rectangulaire doublement entrainée", 13ème
Congrès de Mécanique, 11 - 14 Avril 2017.
[23] M.M. Rahman , S. Parvin , N.A. Rahim , M.R. Islam, R. Saidur , M. Hasanuzzaman
“Effects of Reynolds and Prandtl number on mixed convection in a ventilated cavity with a
heat-generating solid circular block “, Applied Mathematical Modelling, 2012, 36 (5) : 2056-
2066.
[24] M. ADNANI, B. MEZIANI and O. OURRAD, ‟Convection mixte dans une cavité carrée
en présence d‟un bloc soumis à une température froide‟, International Journal of Scientific
Research & Engineering Technology, 2017, 5 : 17-21.
[25] T. S. Chang, "Effects of a finite section with linearly varying wall temperature on mixed

Page | 47
convection in a vertical channel", Int. J. Heat Mass Transfer, 2007, 50 : 2346 – 2354.
[26] T.V. RADHAKRISHNAN, C. BALAJI and S. P. VENKATESHAN,”Optimization of
multiple heaters in a vented enclosure - A combined numerical and experimental study”,
International Journal of Thermal Sciences, 2010, 49: 721-732.
[27] T. Basak, S. Roy, P. K. Sharma, I. Pop, Analysis of mixed convection flows within a square
cavity with linearly heated side walls, Int. J. Heat. Mass Transfer 52, pp.2224-2242, 2009
[28] Prasad, Y. S., & Das, M. K. (2007). Hopf bifurcation in mixed convection flow inside a
rectangular cavity. International Journal of Heat and Mass Transfer, 50(17-18),
[29] Abu-Nada, E. et Chamkha, A. J. (2010), “Mixed convection flow in a lid driven square
enclosure filled with a nanofluid,” European Journal of Mechanic B/Fluids, Vol. 29, No. 6,
pp.472–482.,
[30] Hakan F. Oztop, I. Dagtekin. «Mixed convection in two-sided lid-driven differentially
heated square cavity», International Journal of Heat and Mass Transfer Vol. 47, pp. 1761–
1769,(2004).

Page | 48