Department de Physique

Master
Modélisation et Simulation en Mecanique (MSM)
&
Mécanique et ingénierie

Effets de l'angle d'inclinaison sur la

convection naturelle dans des enceintes
remplies de nanofluide à base de cuivre-eau

Réalisé par : Demandé par :

MARBAH SANAA PR. R. SEHAQUI
SAIH LAMYAA
EL FADILI EL MAHDI

2023/2024
 Table des matières
Chapitre 1 : Etude bibliographique .......................................................................................... 9
1) Convection naturelle .................................................................................................... 9
2) Convection de Rayleigh-Bénard ................................................................................... 9
3) Généralité sur les Nanofluide ......................................................................................10
4) Définition Nanofluide .................................................................................................10
Chapitre 2 : Formulation mathématique de La convection naturelle dans les enceintes remplies
avec nanofluide Cu – Eau ...................................................................................................... 11
1) Configuration de problème .......................................................................................... 11
2) Equation de vorticité ...................................................................................................12
3) Equation de l’énergie ..................................................................................................12
5) Analyse adimensionnelle .............................................................................................13
5-1) Equation de vorticité sans dimension ................................................................13
5-2) Equation de l’énergie sans dimension ...............................................................13
5-3) Equation cinématique sans dimension ...............................................................14
6) Les Conditions aux limites adimensionnelles ..............................................................14
7) Changement de variable ..............................................................................................14
Chapitre III : Résolution Numérique ......................................................................................15
1) Définition de la méthode de volume fini ......................................................................15
2) Principe de la méthode des volumes finis ....................................................................15
3) Maillage ......................................................................................................................16
4) Définition de la méthode de volume fini ......................................................................16
5) Les Etapes de la méthode volume fini .........................................................................16
6) La procédure numérique ..............................................................................................16
7) Discrétisation de notre système ...................................................................................17
8) Enceinte ondulée .........................................................................................................17
9) Discrétisation ..............................................................................................................17
10) Condition sur les murs inférieure et supérieure ........................................................19
11) Le nombre moyen de Nusselt ...................................................................................19
12) Validation de code ...................................................................................................19
Chapitre IV : Résultats et discussion ......................................................................................21
1) Contexte de l'étude : ....................................................................................................21
2) Paramètres de l'étude : .................................................................................................21

2
3) Effets du Nombre de Rayleigh sur les Lignes de Courant et les Isothermes pour un
Nanofluide Cu-Eau à φ = 0,1..............................................................................................21
4) Influence des Nanoparticules sur l'Épaisseur de la Couche Limite Thermique et le
Nombre de Nusselt .........................................................................................................22
5) Effets des Angles d'Inclinaison sur les Lignes de Courant et les Isothermes pour Ra
= 𝟏𝟎𝟓 et φ = 0,1 ............................................................................................................23
6) Effets de la Fraction Volumique sur les Champs d'Écoulement et la Distribution de
Température pour Ra =𝟏𝟎𝟓 et ϕ = 30°................................................................................25
7) Effets des Angles d'Inclinaison et des Concentrations de Nanoparticules sur la Variation
du Nombre de Nusselt par Rapport au Nombre de Rayleigh ...............................................26
8) Variation du Nombre de Nusselt Moyen en Fonction du Nombre de Rayleigh pour
Différents Angles d'Inclinaison ...........................................................................................27
9) Variation du Nombre de Nusselt Local avec un Angle d'Inclinaison de 90° et l'Impact des
Nanoparticules ...................................................................................................................28
10) Variation du Nombre de Nusselt Local en Fonction des Angles d'Inclinaison pour φ =
0,1 et l'Amélioration du Transfert de Chaleur .....................................................................29
11) Amélioration du Transfert de Chaleur en Fonction de l'Angle d'Inclinaison et du
Nombre de Rayleigh pour φ = 0,1 ......................................................................................30
12) Conclusions .............................................................................................................30
13) Références ...............................................................................................................31

3
 Abstract

Les effets de l'angle d'inclinaison sur le transfert de chaleur par convection naturelle et
l'écoulement de fluide dans une enceinte bidimensionnelle remplie de nanofluide à base de
cuivre ont été analysés numériquement. La performance des nanofluides est testée à l'intérieur
d'une enceinte en tenant compte de la dispersion des particules solides. L'angle d'inclinaison
est utilisé comme paramètre de contrôle pour l'écoulement et le transfert de chaleur. Il a été
varié de 0° à 120°.

Les équations régissant le problème sont résolues à l'aide de la technique des volumes
finis pour une gamme de nombres de Rayleigh entre 103 ≤ 𝑅𝑎 ≤ 105. Il est constaté que
l'effet de la concentration de nanoparticules sur le nombre de Nusselt est plus prononcé à faible
fraction volumique qu'à haute fraction volumique. L'angle d'inclinaison peut être un paramètre
de contrôle pour une enceinte remplie de nanofluide. Le pourcentage d'amélioration du
transfert de chaleur en utilisant des nanoparticules diminue pour des nombres de Rayleigh plus
élevés.

4
 Nomenclature
Cp Chaleur spécifique à pression constante (kJ 𝑘𝑔−1 𝑘 −1 )
g Accélération gravitationnelle (m/𝑠 2 )
H Hauteur de l'enceinte (m)
h Coefficient de transfert thermique local (W. 𝑚−2 . 𝑘 −1 )
h Longueur sans dimension du chauffage partiel, h’/H
h’ Longueur du radiateur (m)
Gr Nombre de Grash
k Conductivité thermique (W . 𝑚−1 . 𝑘 −1 )
Nu Nombre de Nusselt, Nu = hH/k
𝑁𝑢𝑎𝑣𝑔 Nombre moyen Nusselt
Pr Nombre de Prandtl
𝑞𝑤 Flux thermique (W. 𝑚−2 )
Ra Nombre de Rayleigh
T Température (K)
u, v Composantes de la vitesse (m. 𝑠 −1 )
U, V Composantes sans dimension de la vitesse
W Longueur de l'enceinte (m)
x, y Coordonnées sans dimension
𝑦𝑝 Centre de chauffage sans dimension
𝑦𝑝′ Centre d'emplacement du radiateur (m)

Symboles grecs
𝛼 Diffusivité thermique en (𝑚2/𝑠)
𝛽 Coefficient de dilatation thermique du fluide en (1/𝐾)
ℇ Tolérance numérique
φ Fraction volumétrique de nanoparticules
ϕ Quantité transportée
𝜗 Viscosité cinématique en (𝑚2/𝑠)
𝜃 Température sans dimension

5
 Ψ Fonction de flux sans dimension
𝜓 Fonction de flux avec dimension (𝑚2/𝑠)
Ω Fonction de vorticité adimensionnelle
𝜔 Tourbillon dimensionnel (𝑠 −1 )
𝜌 Densité en (𝑘𝑔/𝑚3 )
𝜇 Viscosité dynamique en (N.s. 𝑚-2)
Indices
Avg moyenne
nf nanofluide
f fluide
H chaud
L froid
s solide
W paroi
P particule
Nombres sans dimensions
𝑅𝑒 = 𝑢𝐿/𝜗 nombre de Reynolds
𝑅𝑎 = 𝑔H3(𝑇𝐻 − 𝑇𝐿 )𝛽/𝜗𝛼 nombre de Rayleigh
𝑃𝑟 = 𝜗 /𝛼 nombre de Prandtl
𝑁𝑢 = hH/k nombre de Nusselt
Gr= gβ𝐻3 (𝑇𝐻 − 𝑇𝐿 )/ ϑ2 Grashof

6
 Liste de figures
Figure 1 schéma d’une enceinte inclinée
Figure 2 Volume de contrôle bidimensionnel
Figure 3 Schéma représentant le maillage décalé
Figure 4 (a) Nombre de Nusselt par rapport au nombre Ra et comparaison avec d'autres
publications publiées travaux, et (b) comparaison entre les travaux actuels et d’autres données
publiées pour la répartition de la température sur la paroi gauche (Ra = 105 , Pr = 0,7, ϕ = 0°).
Figure 5 Streamlines (à gauche) et isothermes (à droite) pour ϕ = 0° et φ = 0,1, (a) Ra = 103
(𝜓𝑚𝑖𝑛 = -1,156), (b) Ra = 104 (𝜓𝑚𝑖𝑛 = -5,228) et (c) Ra = 105 (𝜓𝑚𝑖𝑛 = 11,22), (fluide pur
(lignes pointillées —), lignes (nanofluides)).
Figure 6 Streamlines (à gauche) et isothermes (à droite) pour Ra =105 et φ = 0,1, (a) ϕ = 30°,
(𝜓𝑚𝑖𝑛 = -17,42 (fluide pur)), 𝜓𝑚𝑖𝑛 = -22,55 (nanofluide)), (b) ϕ = 60° ,( 𝜓𝑚𝑖𝑛 = -24,08 (fluide
pur), 𝜓𝑚𝑖𝑛 = -31,44 (nanofluide)), (c) ϕ = 90° , (𝜓𝑚𝑖𝑛 = -16,13 (nanofluide)), (𝜓𝑚𝑎𝑥 = 20,88
(nanofluide)), (d) ϕ = 120° , (𝜓𝑚𝑎𝑥 = 26 (fluide pur), 𝜓𝑚𝑎𝑥 = 31,44 (nanofluide)), (fluide pur
(lignes pointillées —), lignes (nanofluides)).
Figure 7 Rationalises (à gauche) et isothermes (à droite) pour Ra = 105 et ϕ = 30°, (a) φ = 0,
(𝜓𝑚𝑖𝑛 = -17,42), b) φ = 0,05, (𝜓𝑚𝑖𝑛 = -18,49).
Figure 8 Variation des nombres de Nusselt moyens avec les nombres de Rayleigh à différentes
valeurs de fraction volumique pour différents angles d'inclinaison (a) ϕ = 0° , (b) ϕ = 30° , (c)
ϕ = 60° , (d) ϕ = 90° ,(e) ϕ = 120°
Figure 9 Variation des nombres de Nusselt moyens avec des angles d'inclinaison pour φ = 0,1.
Figure 10 Profils de vitesse au milieu de la paroi chaude pour différents angles d'inclinaison
pour φ = 0,1 et Ra =105 .
Figure 11 Variation du nombre de Nusselt le long de la paroi chauffée pour ϕ =90°.
Figure 12 Variation du nombre de Nusselt local le long de la paroi chauffée pour φ = 0,1 à
différents angles d'inclinaison.
Figure 13 Rapport d'amélioration du transfert de chaleur dû à l'ajout de nanoparticules.

7
 Introduction générale
L'angle d'inclinaison joue un rôle crucial dans la dynamique de la convection naturelle au
sein des enceintes remplies de nanofluide CU - eau. Cette introduction générale vise à fournir
un aperçu des effets de cet angle sur ce phénomène thermique complexe.
La convection naturelle, un processus de transfert de chaleur fondamental, est grandement
influencée par les variations de température et de densité d'un fluide. L'ajout de nanoparticules
dans un fluide de base, comme l'eau, crée un nanofluide avec des propriétés thermiques et
fluidiques uniques. Le cuivre (CU), en raison de sa haute conductivité thermique, est souvent
utilisé dans la fabrication de ces nanoparticules, ce qui les rend particulièrement intéressantes
pour l'amélioration des performances thermiques.
Lorsque ces nanofluides sont confinés dans des enceintes et soumis à un angle
d'inclinaison, les forces gravitationnelles entrent en jeu, modifiant la distribution de la
température et du mouvement du fluide. Comprendre comment l'angle d'inclinaison affecte la
convection naturelle dans ces conditions est essentiel pour optimiser le transfert de chaleur
dans diverses applications, telles que le refroidissement électronique, les systèmes de
chauffage solaire et les réacteurs chimiques.
Cette étude explore les variations de la convection naturelle en fonction de l'angle
d'inclinaison dans des enceintes remplies de nanofluide CU - eau. En examinant les profils de
température, les flux de chaleur et les schémas de circulation du fluide, elle vise à identifier les
tendances et les comportements spécifiques à différents angles d'inclinaison. Les résultats
obtenus offriront des perspectives précieuses pour concevoir des systèmes de refroidissement
et de chauffage plus efficaces, en exploitant pleinement les avantages des nanofluides et en
tenant compte des conditions d'inclinaison spécifiques aux applications réelles.

8
 Chapitre 1 : Etude bibliographique

1) Convection naturelle
La convection naturelle, souvent appelée convection libre, est un processus de transfert de
chaleur qui se produit dans les fluides, tels que les liquides et les gaz, en raison de différences
de température. Contrairement à la convection forcée, qui est induite par des forces externes
telles que des pompes ou des ventilateurs, la convection naturelle se produit spontanément en
raison des variations de densité causées par le chauffage ou le refroidissement d'une substance.
La convection naturelle joue un rôle important dans de nombreux phénomènes naturels et
applications industrielles. Par exemple, elle est responsable du mouvement de l'air chaud près
d'un radiateur, des courants océaniques et atmosphériques, ainsi que du refroidissement de
composants électroniques et des systèmes de chauffage solaire.
2) Convection de Rayleigh-Bénard
La convection de Rayleigh-Bénard est un type de convection naturelle qui se produit dans une
couche horizontale plane de fluide chauffée par dessous, dans laquelle le fluide développe un
motif régulier de cellules de convection dite cellules de Bénard. La convection de Rayleigh –
Bénard est l'un des phénomènes de convection les plus étudiés en raison de son accessibilité
analytique et expérimentale.
Le phénomène de Rayleigh-Bénard correspond à l’état instable dans lequel se trouve une
couche de fluide dilatable, couche confinée entre deux plans rigides horizontaux, d’épaisseur d
et soumise à un écart de température ΔT. La structure du fluide est encore stable et au repos, si
cet écart dépasse une valeur critique Tc, la structure devient instable et des mouvements
naissent à l’intérieur du fluide. Ces mouvements augmentent l’échange de la chaleur bien au-
delà de l’échange thermique par simple conduction. C’est-à-dire qu’à la conduction vient
s’ajouter un flux thermique convectif : les particules du fluide en mouvement transportent de
la chaleur en l’évacuant autour d’elles vers d’autres particules plus froides, puis vers
l’extérieur du fluide. L’exemple le plus simple de ce phénomène est celui de l’eau que l’on fait
bouillir dans une casserole sur une plaque chauffante : au voisinage de celle-ci, le fluide se
réchauffe, devient plus léger et se met à monter sous l’effet des forces d’Archimède, tandis
qu’à la surface, au contact de l’air, il se refroidit, devient plus lourd et se met à descendre.
Les instabilités de Rayleigh-Bénard sont contrôlées par une concurrence entre plusieurs forces,
la poussée d’Archimède (la force motrice), la dissipation de quantité de mouvement (le terme
visqueux) et la diffusivité thermique. Ces deux dernières représentent les forces de freinage.

9
3) Généralité sur les Nanofluide
Le développement considérable récent des recherches traitant des nanofluides est dû au fait
qu’il est possible pour certaines applications L'amélioration du transfert de chaleur dans ces
systèmes est un sujet essentiel du point de vue des économies d'énergie en introduisant dans un
fluide pur une faible concentration de nanoparticules. Cette nouvelle génération de fluides
porte le nom de ″nanofluides″. Ce terme a été introduit par Choi au laboratoire d’Argonne en
U.S.A 1995 et reste couramment utilisé pour caractériser ce type de suspension colloïdale.
Récemment, une nouvelle classe de nanofluides, les « nanofluides hybrides », est utilisée pour
améliorer encore le taux de transfert de chaleur obtenus en dispersant, dans un fluide base,
deux sortes de nanoparticules.
4) Définition Nanofluide
L’ajout des particules nanométriques dans un fluide de base "appelé nanofluide" Les
nanofluides sont des suspensions colloïdales. Le nanofluide est un fluide contenant des
particules de taille nanométrique (inférieure à 100 nm) appelées nanoparticules, l’ajout des
nanoparticules aux fluides de base avec des proportions adéquates (ne dépassant pas 10%)
pour éviter l’agrégation. Dans le cas des fluides caloporteurs, un des premiers paramètres à
prendre en compte afin d'évaluer le potentiel d'échange de chaleur est la conductivité
thermique. Les types de nanoparticules sont :
− Métalliques : le cuivre (𝐶𝑢), l’aluminium (𝐴𝑙), l’or (𝐴𝑢), l’argent (𝐴𝑔).
− Oxydes : l’oxyde de cuivre (𝐶𝑢𝑂), l’oxyde de l’aluminium (𝐴𝑙2𝑂3), le de titane (𝑇𝑖𝑂2).
− Nanotubes de carbone
Le choix des liquides de base est primordial pour un nanofluide, il permit d’assurer la stabilité
de la suspension dans le temps et d’éviter tout phénomène d’agrégations, la sélection d’un tel
fluide se fera en fonction de la nature de nanoparticules, les solvants les plus utilisés sont :
− L’eau.
− L’éthylène glycol, EG.
− Les huiles.
− Le Toluène.

10
 Chapitre 2 : Formulation mathématique de La convection
naturelle dans les enceintes remplies avec nanofluide Cu –
Eau
1) Configuration de problème
Le tracé du modèle considéré est représenté sur la figure 1 avec les coordonnées. Il s'agit d'une
enceinte carrée bidimensionnelle avec une inclinaison angle. Sa paroi verticale a une
température isotherme différente lorsque l'angle d'inclinaison est nul. Les murs restants sont
adiabatiques. L'enceinte est remplie de nanofluide newtonien, incompressible et laminaire.

Figure 1 schéma d’une enceinte inclinée

Tableau 1 :
Propriétés thermophysiques du fluide pur et des nanoparticules.

La figure 1 montre un diagramme schématique de l'enceinte partiellement chauffée. Le fluide

contenu dans l'enceinte est un nanofluide à base d'eau contenant des nanoparticules de cuivre.
On suppose que le fluide de base (c.-à-d’eau) et les nanoparticules sont en équilibre thermique
et non un glissement se produit entre eux. Les propriétés thermophysiques du nanofluides sont
donnés dans le tableau 1. La paroi gauche est maintenue à une température constante (TH)
supérieure à la paroi droite (TL). Les propriétés thermophysiques du nanofluide sont
supposées être constante sauf pour la variation de densité, qui est approximée par le modèle
Boussinesq.
Les équations qui régissent la convection naturelle laminaire et stationnaire en termes de la
fonction courant-vorticité.

11
2) Equation de vorticité
L’équation de vorticité pour la convection naturelle laminaire et stationnaire exprime comment
la vorticité évolue en raison des forces de flottabilité (provoquées par les gradients de
température) et de la diffusion visqueuse.
La forme typique de l'équation de vorticité pour un fluide incompressible et en régime de
convection naturelle peut être écrite comme suit :
𝜕𝜔 𝜕𝜔 µ𝑛𝑓 2 𝜑𝜌𝑠 𝛽𝑠 + (1 − 𝜑)𝜌𝑓 𝛽𝑓 𝜕𝑇 𝜕𝑇
𝑢 +𝑣 = ∇ 𝜔+ 𝑔( ′ cos 𝜙 − ′ sin 𝜙)
𝜕𝑥′ 𝜕𝑦′ 𝜌𝑛𝑓 𝜌𝑛𝑓 𝜕𝑥 𝜕𝑦
𝜕𝜓 𝜕𝜓
En injectant ces expressions 𝑢 = 𝜕𝑦′ , 𝑣 = − 𝜕𝑥′ (u la vitesse radiale, v la vitesse tangentielle),
dans l’équation, on obtient l'équation de vorticité avec dimension.
𝜕 𝜕𝜓 𝜕 𝜕𝜓 µ𝑛𝑓 𝜕2𝜔 𝜕2 𝜔 𝜑𝜌𝑠 𝛽𝑠 +(1−𝜑)𝜌𝑓 𝛽𝑓 𝜕𝑇 𝜕𝑇
𝜕𝑥 ′
(𝜔
𝜕𝑦
)
′ − 𝜕𝑦 ′
(𝜔 𝜕𝑥 ′ ) = 𝜌 (𝜕𝑥′2 + 𝜕𝑦′2 ) + 𝜌𝑛𝑓
𝑔(𝜕𝑥 ′ cos 𝜙 − 𝜕𝑦′ sin 𝜙) (1)
𝑛𝑓

Avec,
 µ𝑛𝑓 : La viscosité du nanofluide peut être approximée à celle d'un fluide de base
contenant une suspension diluée de fines particules sphériques et est donnée par
µ𝑓
Brinkman (1952). Donnée par µ𝑛𝑓 = (1−𝜑) 2,5

 𝜌𝑛𝑓 : La densité effective du nanofluide est donnée par 𝜌𝑛𝑓 = (1 − 𝜑)𝜌𝑓 + 𝜑𝜌𝑠

3) Equation de l’énergie
L'équation de l'énergie dans la convection naturelle est une expression mathématique qui décrit
le transfert de chaleur dans un fluide en mouvement sous l'effet de la gravité. Cette équation
tient compte des variations de température, de densité et de mouvement du fluide. Elle est
basée sur la conservation de l'énergie thermique et intègre les effets de la conduction, de la
convection et parfois de la dissipation visqueuse.
En général, l'équation de l'énergie pour la convection naturelle dans un fluide incompressible
est donnée par :
𝜕𝑇 𝜕𝑇 𝐾 𝜕𝜓 𝜕𝜓
𝑢 𝜕𝑥′ + 𝑣 𝜕𝑦′ = 𝛼𝑛𝑓 ∆𝑇 Avec 𝛼𝑛𝑓 = (𝜌𝐶𝑒𝑓𝑓 si l’on injecte ces exprissions 𝑢 = 𝜕𝑦′ , 𝑣 = − 𝜕𝑥′
) 𝑝 𝑛𝑓

l’équation devient
𝜕 𝜕𝜓 𝜕 𝜕𝜓 𝜕 𝜕𝑇 𝜕 𝜕𝑇
𝜕𝑥 ′
(𝑇
𝜕𝑦 ′
)−
𝜕𝑦 ′
(𝑇
𝜕𝑥 ′
)=
𝜕𝑥 ′
(𝛼𝑛𝑓
𝜕 𝑥′
)+
𝜕𝑦 ′
(𝛼𝑛𝑓
𝜕 𝑦′
) (2)

Avec,
 (𝜌𝐶𝑝 ) : La capacité thermique du nanofluide est exprimée selon Abu-Nada (2008) et
𝑛𝑓
Khanafer et al. (2003) : (𝜌𝐶𝑝 ) = (1 − 𝜑)(𝜌𝐶𝑝 )𝑓 + 𝜑(𝜌𝐶𝑝 )𝑠
𝑛𝑓
 La conductivité thermique effective du nanofluide est approximée par le modèle de
Maxwell-Garnetts.
𝑘𝑛𝑓 𝑘𝑠 +2𝑘𝑓 −2𝜑 (𝑘𝑓 − 𝑘𝑠 )
=
𝑘𝑓 𝑘𝑠 + 2𝑘𝑓 + 𝜑 (𝑘𝑓− 𝑘𝑠 )

12
4) Equation cinématique
L’équation cinématique avec dimension s’écrit sous forme
𝜕2 𝜓 𝜕2 𝜓
𝜕𝑥′2
+ 𝜕𝑦′2 = −𝜔 (3)

5) Analyse adimensionnelle
L'analyse dimensionnelle permet de simplifier les équations de base, de mettre en évidence les
paramètres les plus influents et de fournir des expressions générales qui peuvent être utilisées
pour prédire le comportement de différents systèmes.
Les équations précédentes peuvent s’écrire d’une façon adimensionnelle en introduire les
𝑥′ 𝑦′ 𝜓 𝜔𝐻 2 𝜗𝐻 𝑢𝐻
paramètres adimensionnels suivantes : 𝑥 = , 𝑦= , Ψ=𝛼 , Ω= , 𝑉= , 𝑈= ,
𝐻 𝐻 𝑓 𝛼𝑓 𝛼𝑓 𝛼𝑓
𝑇−𝑇𝐿
𝜃=𝑇
𝐻 −𝑇𝐿

5-1) Equation de vorticité sans dimension

En remplaçant les paramètres adimensionnels dans l'équation (1), on trouve
𝜕 𝜕Ψ 𝜕 𝜕Ψ 𝜇𝑛𝑓 𝜕2 Ω 𝜕2 Ω 𝜑𝜌𝑠 𝛽𝑠 +(1−𝜑)𝜌𝑓 𝛽𝑓 𝜕𝜃 𝜕𝜃
𝜕𝑥
(Ω
𝜕𝑦
)−
𝜕𝑦
(Ω
𝜕𝑥
)=
𝜌𝑛𝑓 𝛼𝑓
(
𝜕𝑥 2
+ 𝜕𝑦2 ) + 𝜌𝑛𝑓 𝛼𝑓2
𝐻3 (𝑇𝐻 − 𝑇𝐿 )𝑔(𝜕𝑥 cos 𝜙 − 𝜕𝑦 sin 𝜙)

Et par une simplification, on trouve l'expression finale de l'équation de vorticité sans

dimension

𝜕 𝜕Ψ 𝜕 𝜕Ψ 𝑃𝑟 𝜕2 Ω 𝜕2 Ω 1 𝛽𝑠
(Ω )− (Ω ) = [(1−𝜑)0.25 𝜌 ]( + 𝜕𝑦2 ) + 𝐺𝑟𝑃𝑟 2 [(1−𝜑)𝜌𝑓 +
𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑥 ((1−𝜑)+𝜑 𝑠 ) 𝜕𝑥 2 +1 𝛽𝑓
𝜌𝑓 𝜑 𝜌𝑠

1 𝜕𝜃 𝜕𝜃
𝜑 𝜌𝑓
]( cos 𝜙 − 𝜕𝑦 sin 𝜙) (4)
+1 𝜕𝑥
(1−𝜑) 𝜌𝑠

5-2) Equation de l’énergie sans dimension

En injectant les paramètres adimensionnels dans l’équation (2), on trouve
𝜕 𝜕Ψ 𝜕 𝜕Ψ 𝛼𝑛𝑓 𝜕 2 𝜃 𝜕 2 𝜃
(𝜃 ) − (𝜃 )= ( + )
𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑥 𝛼𝑓 𝜕𝑥 2 𝜕𝑦 2
En simplifiant, on obtient
𝜕 𝜕Ψ 𝜕 𝜕Ψ 𝜕 𝜕𝜃 𝜕 𝜕𝜃
𝜕𝑥
(𝜃
𝜕𝑦
)−
𝜕𝑦
(𝜃
𝜕𝑥
)=
𝜕𝑥
(𝜆
𝜕𝑥
)+
𝜕𝑦
(𝜆
𝜕𝑦
) (5)
𝑘𝑛𝑓
𝑘𝑓
Avec 𝜆= (𝜌𝐶𝑝)
(1−𝜑)+𝜑 𝑠
(𝜌𝐶𝑃 )𝑓

13
 5-3) Equation cinématique sans dimension
En injectant les paramètres adimensionnels dans l’équation (3), on parvient à
𝛼𝑓 𝜕 2 Ψ 𝛼𝑓 𝜕 2 Ψ 𝛼𝑓 𝜕2 Ψ 𝜕2 Ψ
𝐻2 𝜕𝑥 2
+ 𝐻 2 𝜕𝑦2 = −Ω 𝐻 2 Après simplification, l'équation devient 𝜕𝑥 2
+ 𝜕𝑦2 = −Ω (6)

6) Les Conditions aux limites adimensionnelles

Les conditions aux limites adimensionnelles sont écrites comme suit :
 Sur le mur gauche : x = 0 ; Ψ = 0 ; Ω = - ∂² Ψ /∂x² ; 𝜃 = 1
 Sur le mur droit : x = 1 ; Ψ = 0 ; Ω = - ∂² Ψ /∂x² ; 𝜃 = 0
 Sur les murs supérieur et inférieur : Ψ = 0 ; Ω = - ∂² Ψ /∂y² ; ∂𝜃/∂y = 0
7) Changement de variable
Le changement de variable utilisé est défini comme suit :
𝑋 = 𝑋 (, ɳ) , 𝑌 = 𝑌(, ɳ)
𝑑𝑈 𝑑𝑢
𝑑𝑈 = 𝑑𝑋 𝑑𝑋 + 𝑑𝑦 𝑑𝑦

Finalement on obtient une nouvelle équation de transport en fonction de  et  :

𝜕 𝜕 𝜕 𝑇 𝜕 𝜕 𝜕 𝑇 𝜕 𝜕
( 𝑈 ∗ ) + (𝑉 ∗ ) = ( (𝑎 𝜕 − 𝛽 𝜕ɳ ) + 𝜕ɳ ( 𝐽 (𝑎 𝜕 − 𝛾 𝜕ɳ ) (7)
𝜕 𝜕ɳ 𝜕 𝐽

Avec

𝑑𝑦 𝑑𝑥 𝜕𝑥 𝑑𝑦 𝜕𝑥 2 𝜕𝑦 2
𝑈 ∗ = (( 𝑑ɳ ) 𝑈 − ( 𝑑ɳ ) 𝑉) , 𝑉 ∗ = ((𝜕𝜉 ) − ( 𝑑 ) 𝑈) , 𝑎𝜙 = ((𝜕𝜂) + (𝜕𝜂 ) )

𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑥 2 𝜕𝑦 2
𝛽𝜙 =(𝜕𝜂) (𝜕𝜉 ) + (𝜕𝜉 ) (𝜕𝜂 ) , 𝛾 = ((𝜕𝜉 ) + ( 𝜕𝜉 ) )

14
 Chapitre III : Résolution Numérique
1) Définition de la méthode de volume fini
La méthode des volumes finis est une technique utilisée pour résoudre numériquement
des équations aux dérivées partielles (EDP). Elle est comparable à la méthode des différences
finies et à celle des éléments finis. Cependant, contrairement à la méthode des différences
finies, qui approxime les dérivées, la méthode des volumes finis exploite directement la forme
dite « forte » de l’équation à résoudre.
2) Principe de la méthode des volumes finis
La méthode des volumes finis est une approche numérique couramment utilisée pour résoudre
des équations aux dérivées partielles (EDP) dans des domaines complexes. Voici une
reformulation de votre description :
Dans la simulation par la méthode des volumes finis, le domaine de calcul est discrétisé en un
ensemble fini de sous-domaines élémentaires appelés “volumes de contrôle”. La méthode des
volumes finis consiste à intégrer les équations aux dérivées partielles (EDP), décrites
précédemment, sur chaque volume de contrôle. Chaque volume de contrôle contient un nœud
principal. Par exemple, pour un nœud principal (noté ‘P’), les points E (Est) et W (Ouest) sont
des voisins dans la direction ξ, tandis que les points N (Nord) et S (Sud) sont des voisins dans
la direction η. Le volume de contrôle entourant le nœud principal (‘P’) est délimité par des
lignes discontinues. Les faces du volume de contrôle sont localisées aux points (e) et (w) dans
la direction ξ, et aux points (n) et (s) dans la direction η.

 

  
(η)s

(ξ)w
 (ξ)e

S
Figure 2 Volume de contrôle bidimensionnel.

15
3) Maillage
Le maillage est la subdivision du domaine d’étude en grilles longitudinales et transversales
dont l’intersection représente un nœud.
La discrétisation du domaine est obtenue par un maillage constitué d’un réseau de points
(nœuds). Ainsi un élément de volume (volume de contrôle) est défini autour de chaque nœud.
Les grandeurs scalaires sont stockées dans le nœud (P’) du maillage, tandis que les grandeurs
vectorielles sont stockées aux milieux des segments reliant les nœuds. L’équation générale de
transport est intégrée sur le volume de contrôle associé aux variables scalaires et les équations
de quantité de mouvement sont intégrées sur le volume de contrôle associé aux composantes
de la vitesse.
Le volume de contrôle de la composante longitudinale (U) est décalé suivant la direction (ξ)
par rapport au volume de contrôle principal, celui de la composante transversale décalé suivant
la direction (η). Ce type de maillage dit : « maillage décalé » permet une bonne approximation
des flux convectifs et une stabilisation numérique de la solution.
La construction des volumes de contrôle et le maillage décalé sont montrés dans le schéma
suivant

Figure 3 schéma représentant le maillage décalé

16
4) Définition de la méthode de volume fini
La méthode des volumes finis est une technique utilisée pour résoudre numériquement
des équations aux dérivées partielles (EDP). Elle est comparable à la méthode des différences
finies et à celle des éléments finis. Cependant, contrairement à la méthode des différences
finies, qui approxime les dérivées, la méthode des volumes finis exploite directement la forme
dite « forte » de l’équation à résoudre.
5) Les Etapes de la méthode volume fini
Voici les étapes générales impliquées dans la méthode de volume fini :
 Définir l’équation du phénomène :
Tout d’abord, identifiez l’équation aux dérivées partielles (EDP) que vous souhaitez
résoudre. Cela peut être une équation de convection, de diffusion, ou une combinaison des
deux.
 Discrétisation du domaine :
Divisez le domaine de calcul en volumes finis. Ces volumes finis sont de petits volumes
disjoints en 3D (ou des surfaces en 2D, des segments en 1D) qui couvrent l’ensemble du
domaine d’étude.
 Intégration de l’équation sur chaque volume élémentaire :
Utilisez le théorème de Gauss (ou le théorème d’Ostrogradski) pour transformer les intégrales
de volume d’un terme de divergence en intégrales de surface. Évaluez ensuite les termes de
flux aux interfaces entre les volumes finis. Pour élaborer une approximation des flux aux
interfaces, utilisez une fonction de flux numérique.

 Traitement des conditions aux limites :

Prenez en compte les conditions aux limites spécifiques au problème que vous résolvez. Ces
conditions peuvent être des conditions de Dirichlet (valeurs imposées), des conditions de
Neumann (flux imposés), ou d’autres types de conditions.
6) La procédure numérique
L'équation d'énergie Eq. (5) est initialement résolue dans le temps en utilisant le schéma de
l'implicit Direction Alternée (ADI). Les résultats de température obtenus sont ensuite utilisés
dans l'équation de vorticité Eq. (1). L'équation de transport de la vorticité instationnaire dans
le temps est également résolue en utilisant le schéma de l'implicit Direction Alternée. Un pas
de temps constant sans dimension de 0,001 est utilisé lors de l'analyse numérique. Des critères
de convergence de 10−6 sont utilisés lors des calculs de la vitesse et de la fonction de courant,
comme indiqué dans l'Eq suivante
𝑗=𝑀
∑𝑗=1 ∑𝑖=𝑁
𝑖=1 |𝜙
𝑛+1
+𝜙𝑛 |
ℇ= 𝑗=𝑀 𝑛+1
(8)
∑𝑗=1 ∑𝑖=𝑁
𝑖=1 | 𝜙 |

Où ∈ représente la fonction de courant (𝜓) ou les composantes de la vitesse.

16
7) Discrétisation de notre système
Les équations de conservation présentées au chapitre précédent peuvent être écrites sous une
forme commune pour chaque enceinte étudiée (carrée et ondulée). Cette formulation permet
de ne pas répéter le travail de discrétisation pour chaque équation
8) Enceinte ondulée
Si on note  la variable étudiée, chacune des équations peut être réduite à une seule équation
générale, en coordonnées cartésiennes selon la forme :

𝜕 𝜕 𝜕 𝑇 𝜕 𝜕 𝜕 𝑇 𝜕 𝜕
(𝑈 ∗ ) + (𝑉 ∗ ) = ( 𝐽 (𝑎 𝜕 − 𝛽 𝜕ɳ ) + 𝜕ɳ ( 𝐽 (𝑎 𝜕 − 𝛾 𝜕ɳ ) (9)
𝜕 𝜕ɳ 𝜕

Avec
𝜕 𝜕
(𝑈 ∗ ) + (𝑉 ∗ )= Terme convectif
𝜕 𝜕ɳ

𝜕 𝑇 𝜕 𝜕 𝜕 𝑇 𝜕 𝜕
( 𝐽 (𝑎 𝜕 − 𝛽 𝜕ɳ ) + 𝜕ɳ ( 𝐽 (𝑎 𝜕 − 𝛾 𝜕ɳ )= Terme diffusif.
𝜕

T = Coefficient de diffusion
L’évaluation des différentes intégrales de l’équation de transport ci-dessous demande que l’on
fasse un certain nombre d’approximations de façon à former un système d’équation linéaires
reliant entre elles les valeurs de la fonction  aux points du maillage(𝑃(𝑖,𝑗)
′
).

𝑇 𝜕 𝜕 𝑇 𝜕 𝜕
(𝑉 ∗ )) 𝑑𝑣 = ∫ ( (  (𝑎
𝜕 𝜕 𝜕 𝜕
∫ ( (𝑈 ∗ ) + − 𝛽 )+ ( (𝑎 − 𝛾 ))𝑑𝑣
𝜕 𝜕ɳ 𝜕 𝐽 𝜕 𝜕ɳ 𝜕ɳ 𝐽 𝜕 𝜕ɳ
𝑣𝑐 𝑣𝑐

Avec
𝑑𝑣 = 𝑑. 𝑑ɳ. 1
9) Discrétisation
 Discrétisation différents termes de l’équation de transport
Terme convectif suivant 
𝜕
∫ (𝑈 ∗ ) 𝑑𝑣 = ((U ∗ )𝒆 − (U ∗ )𝑤 )∆ɳ = 𝐹𝑒 𝑒 − 𝐹𝑤 𝑤 (10)
𝜕
𝒗𝒄

Avec 𝐹𝑒 = 𝑈𝑒∗ ∆ɳ , 𝐹𝑤 = 𝑈𝑤∗ ∆ɳ

Terme convectif suivant η
𝜕
∫ (𝑉 ∗ ) =( (𝑉 ∗ )𝑛 − (𝑉 ∗ )𝑠 )∆ = 𝐹𝑛 𝑠 − 𝐹𝑠 𝑠 (11)
𝜕ɳ
𝑣𝑐

Avec 𝐹𝑛 = 𝑉𝑛∗ ∆ , 𝐹𝑠 = 𝑉𝑠∗ ∆

17
Terme diffusif suivant 
On a fait l’intégrale de terme diffusif entre e et w on trouve expression suivante
𝑇 𝑇 𝜕 𝜕 𝑇 𝜕 𝜕
∫ ( (  (𝑎 − 𝛽 ) 𝑑 𝑑ɳ=(( (  (𝑎 − 𝛽 ) − ( (  (𝑎 − 𝛽 ) )∆ɳ
𝜕 𝜕 𝜕 𝜕 𝜕
𝜕 𝐽 𝜕 𝜕ɳ 𝜕 𝐽 𝜕 𝜕ɳ 𝜕 𝐽 𝜕 𝜕ɳ 𝑒 𝑤
𝒗𝒄

Finalement on obtient
𝜕 𝑇 𝜕 𝜕
∫ ( (  (𝑎 − 𝛽 ) 𝑑 𝑑ɳ = (D𝑛𝑒 (𝑒 − 𝑝 )- 𝐷𝑤 (𝑝 − 𝑤 )+
𝜕 𝐽 𝜕 𝜕ɳ
𝑣𝑐

𝐷𝑛𝑒 (𝑁 − 𝑠 + 𝑁𝐸 − 𝑆𝑊 )-𝐷𝑛𝑤 (𝑁 −𝑆 +𝑁𝑊 − 𝑆𝑊 ) (12)
𝑇𝑒 𝑎 ∆ɳ 𝑇𝑤 𝑎 ∆ɳ 𝑇𝑒 𝛽 𝑇𝑤𝛽
Avec 𝐷𝑒 = , 𝐷𝑤 = , 𝐷𝑛𝑒 = −( ) , 𝐷𝑛𝑤 = −( )
𝐽𝑝𝑒 𝐽𝑤𝑝 4𝐽 4𝐽

Terme diffusif suivant η

On a fait l’intégrale de terme diffusif entre n et s on trouve expression suivante
𝜕 𝑇 𝜕 𝜕 𝑇 𝜕 𝜕 𝑇 𝜕 𝜕
∫ ( (𝑎 𝜕 − 𝛾 𝜕ɳ ) 𝑑𝜉 ⅆ𝜂=[( 𝐽 (𝑎 𝜕 − 𝛾 𝜕ɳ ))𝑛 − ( 𝐽 (𝑎 𝜕 − 𝛾 𝜕ɳ ))𝑠 ]∆
𝜕ɳ 𝐽
𝑽𝒄

Finalement on obtient
𝜕 𝑇 𝜕 𝜕
∫ ( 𝐽 (𝑎 𝜕 − 𝛾 𝜕ɳ ) 𝑑𝜉 ⅆ𝜂= (𝐷𝑛 (𝑁 − 𝑃 ) − 𝐷𝑠 (𝑃 − 𝑠 )+𝐷𝑠𝑒 (𝐸 −𝑤 + 𝑁𝐸 −
𝜕ɳ
𝑣𝑐
𝑁𝑊 ) − 𝐷𝑁𝑊 (𝐸 − 𝑤 + 𝑆𝐸 − 𝑤 ) (13)
𝑇𝑛 𝛾 ∆ 𝑇𝑠𝛾 ∆ 𝑇𝑛𝛽 𝑇𝑠𝛽
Avec 𝐷𝑛 = , 𝐷𝑠 = , 𝐷𝑠𝑒 = −( ) , 𝐷𝑠𝑤 = −( )
𝐽ɳ𝑃𝑁 𝐽ɳ𝑃𝑆 4𝐽 4𝐽

Après une substitution dans l’équation de transport on obtient la forme discrétisée suivante :

(𝐹𝑒 𝑒 − 𝐹𝑤 𝑤 )+( 𝐹𝑛 𝑠 − 𝐹𝑠 𝑠 )= (D𝑛𝑒 (𝑒 − 𝑝 )- 𝐷𝑤 (𝑝 − 𝑤 )+ 𝐷𝑛𝑒 (𝑁 − 𝑠 + 𝑁𝐸 −

𝑆𝑊 )-𝐷𝑛𝑤 (𝑁 −𝑆 +𝑁𝑊 − 𝑆𝑊 )) + (𝐷𝑛 (𝑁 − 𝑃 ) − 𝐷𝑠 (𝑃 − 𝑠 )+𝐷𝑠𝑒 (𝐸 −𝑤 + 𝑁𝐸 −
𝑁𝑊 ) − 𝐷𝑁𝑊 (𝐸 − 𝑤 + 𝑆𝐸 − 𝑤 )) (14)
Résolution des systèmes d’équations linéaires (méthode de Jacobi)
Afin de discrétiser l’équation de transport et de Poisson, on a la relation générale ayant la
forme suivante
′
𝑎𝑤 𝜙𝑖−1,𝑗 + 𝑎𝑒′ 𝜙𝑖+1,𝑗 + 𝑎𝑠′ 𝜙𝑖,𝑗−1 + 𝑎𝑛′ 𝜙𝑖,𝑗+1 + 𝐷𝑛𝑒
′ ′
𝜙𝑖+1,𝑗+1 + 𝐷𝑛𝑤 𝜙𝑖−1,𝑗+1
𝑎𝑃 𝑖,𝑗 = { ′ ′ } (15)
𝐷𝑠𝑒 𝜙𝑖+1,𝑗−1 + 𝐷𝑠𝑤 𝜙𝑖−1,𝑗−1
𝑘
La méthode de Jacobi consiste à estimer une solution (𝜙𝑖,𝑗 ) dans l’itération (k) pour calculer
(𝑖,𝑗 ) comme suit :
𝑘+1

𝑖−1,𝑗 + 𝑎𝑒′ 𝑘𝑖+1,𝑗 + 𝑎𝑠′ 𝑘𝑖,𝑗−1 + 𝑎𝑛′ 𝑘𝑖,𝑗+1 + 𝐷𝑛𝑒

′ 𝑘
𝑎𝑤 ′
𝑘𝑖+1,𝑗+1 + 𝐷𝑛𝑤
′
𝑘𝑖−1,𝑗+1
𝑘+1 ={ } (16)
𝑖,𝑗 ′ 𝑘
+𝐷𝑠𝑒 𝑖+1,𝑗−1 + 𝐷𝑠𝑤
′
𝑘𝑖−1,𝑗−1
La condition de convergence dans (8)

18
Si la condition (8) n’est pas encore satisfaite, on pose que : 𝑘𝑖,𝑗 = 𝑘+1
𝑖,𝑗
et on recalcule la
même formule jusqu’à ce qu’elle vérifie la condition de convergence.
 Discrétisation les vitesses

1 (𝛹𝑖+1,j) 𝜕𝑦 (𝛹𝑖,𝑗+1) 𝜕𝑦
𝑉𝑖,𝑗 = 𝑗 [ (𝜕𝜂 ) − ( 𝜕𝜉 ) ]
2𝛥𝜉 𝑖,𝑗 2𝛥𝜂 𝑖,𝑗

1 (𝛹𝑖,𝑗+1 ) (𝛹𝑖,𝑗+1 ) 𝜕𝑥
𝑈𝒊,𝒋 = [ − ( ) ]
𝑗 2𝛥𝜉 2𝛥𝜂 𝜕𝜂 𝑖,𝑗

 Discrétisation la fonction de courant

𝛹𝑖,𝑗 = 𝑃∗ (𝛹𝑖+1,𝑗 + 𝛹𝑖−1,𝑗 ) + 𝑄 ∗ (𝛹𝑖,𝑗+1 + 𝛹𝑖,𝑗−1 ) − 𝑅
𝛥𝑦 2 𝛥𝑥 2
Avec 𝑃∗ = (𝛥𝑥 2 + 𝛥𝑦 2 ) , 𝑄∗= (𝛥𝑥 2 + 𝛥𝑦 2 ), 𝑅 = 2. 𝑃∗ 𝛥𝑥 2
2 2

10) Condition sur les murs inférieure et supérieure

Une représentation précise de la vorticité à la surface est la plus étape critique dans la
formulation de la fonction du flux-vorticité. Une formule précise de second ordre est utilisée
pour la condition limite de vorticité. La vorticité au mur inférieur s’exprime par :
(8𝛹𝑖,𝑗−𝛹𝑖+1,𝑗 )
𝛺(𝑖,𝑗) = 2(𝛥𝑦)

11) Le nombre moyen de Nusselt

𝟏
𝑵𝒖𝒗𝒂𝒈 = ∫𝟎 𝑵𝒖(𝒚)𝒅𝒚

12) Validation de code

Pour garantir une solution indépendante du maillage, nous avons effectué une procédure de
test de maillage approfondie. Nous avons utilisé sept combinaisons de maillage différentes
pour le cas de Ra = 105, Pr = 0,7 et un angle d’inclinaison nul. En calculant le nombre de
Nusselt moyen sur la paroi gauche, nous avons déterminé qu’une taille de maillage de 51 × 51
assure une solution indépendante du maillage. La valeur convergée (Nu = 4,644) a été
comparée à d’autres valeurs connues rapportées par d’autres chercheurs. Notre solution
numérique est en bon accord avec d’autres travaux rapportés dans la littérature. La solution
numérique actuelle est validée en comparant les résultats du code actuel pour Ra = 105 et Pr =
0,70 (avec un angle d’inclinaison nul) avec les expériences menées par Krane et Jessee (1983)
ainsi qu’avec les simulations numériques réalisées par Khanafer et al. (2003). Il est clair que
le code actuel est en bon accord avec d’autres travaux rapportés dans la littérature, comme le
montre la figure 4b.

19
Figure 4 (a) Nombre de Nusselt par rapport au nombre Ra et comparaison avec d'autres
publications publiées travaux, et (b) comparaison entre les travaux actuels et d’autres données
publiées pour la répartition de la température sur la paroi gauche (Ra = 105 , Pr = 0,7, ϕ = 0°).

20
 Chapitre IV : Résultats et discussion
1) Contexte de l'étude :
Une analyse numérique a été menée pour étudier les effets de l'angle d'inclinaison sur un
nanofluide Cu-eau remplissant une enceinte carrée inclinée bidimensionnelle. L'enceinte est
chauffée différemment sur deux parois opposées, tandis que les deux autres parois opposées
sont adiabatiques.
2) Paramètres de l'étude :
Les calculs ont été effectués pour diverses valeurs de fraction volumique de nanoparticules (0
≤ φ ≤ 0,1), d'angles d'inclinaison (0° ≤ ϕ ≤ 120°) et de nombres de Rayleigh (103 ≤ Ra ≤
105 ). Les propriétés thermophysiques du fluide pur et des nanoparticules de cuivre sont
répertoriées dans le Tableau 1.
3) Effets du Nombre de Rayleigh sur les Lignes de Courant et les
Isothermes pour un Nanofluide Cu-Eau à φ = 0,1

Figure 5 Lignes de courant (à gauche) et isothermes (à droite) pour ϕ = 0° et φ = 0,1, (a) Ra =

10³ (𝜓𝑚𝑖𝑛 = -1,156), (b) Ra = 10⁴ (wmin = -5,228), et (c) Ra = 10⁵ (𝜓𝑚𝑖𝑛 = -11,22), (fluide
pur (lignes pointillées —), lignes (nanofluides)).

21
  Figure 5a-c :
La figure 5a–c montre les effets du nombre de Rayleigh sur les lignes de courant (à gauche) et
les isothermes (à droite) pour φ = 0,1. La figure présente une comparaison du champ de flux
et de la distribution de température de nanofluide (lignes continues) et fluide pur (lignes en
pointillés). De plus, cette figure est présentée pour comparer les résultats obtenus avec la
littérature (Khanafer et al., 2003). Une cellule singulière se forme à bas nombre de Rayleigh
et la forme de la cellule principale est circulaire. Comme le montre la figure, la forme de la
cellule devient ellipsoïdale avec l'augmentation du nombre de Rayleigh. L'ajout de
nanoparticules affecte le diamètre de la cellule. Dans tous les cas, l'écoulement tourne dans le
sens des aiguilles d'une montre.

4) Influence des Nanoparticules sur l'Épaisseur de la Couche Limite

Thermique et le Nombre de Nusselt
 Figure 5b :
De plus, la Fig. 5 illustre comment l'épaisseur de la couche limite thermique près de la paroi
chauffée est influencée par l'ajout de nanoparticules. Cette sensibilité de l'épaisseur de la
couche limite thermique à la fraction volumique de nanoparticules est liée à l'augmentation de
la conductivité thermique du nanofluide. En effet, des valeurs plus élevées de conductivité
thermique s'accompagnent de valeurs plus élevées de diffusivité thermique. La valeur élevée
de la diffusivité thermique entraîne une diminution des gradients de température et, par
conséquent, augmente l'épaisseur de la couche limite, comme le montre la Fig. 5b. Cette
augmentation de l'épaisseur de la couche limite thermique réduit le nombre de Nusselt,
𝑘 𝜕𝑇
cependant, selon l'équation (𝑁𝑢 = − ( 𝑘𝑛𝑓 ) 𝜕𝑥 ), le nombre de Nusselt est le produit du gradient
𝑓
de température et du rapport de conductivité thermique (conductivité du nanofluide sur la
conductivité du fluide de base). Étant donné que la réduction du gradient de température due à
la présence de nanoparticules est beaucoup plus faible que le rapport de conductivité
thermique, une amélioration du Nusselt se produit en augmentant la fraction volumique de
nanoparticules. L'augmentation du nombre de Nusselt sera abordée plus en détail dans cette
section de discussion.

22
5) Effets des Angles d'Inclinaison sur les Lignes de Courant et les
Isothermes pour Ra = 𝟏𝟎𝟓 et φ = 0,1

23
Figure 6 Streamlines (à gauche) et isothermes (à droite) pour Ra =105 et φ = 0,1, (a) ϕ = 30°,
(𝜓𝑚𝑖𝑛 = -17,42 (fluide pur)), 𝜓𝑚𝑖𝑛 = -22,55 (nanofluide)), (b) ϕ = 60° ,( 𝜓𝑚𝑖𝑛 = -24,08 (fluide
pur), 𝜓𝑚𝑖𝑛 = -31,44 (nanofluide)), (c) ϕ = 90° , (𝜓𝑚𝑖𝑛 = -16,13 (nanofluide)), (𝜓𝑚𝑎𝑥 = 20,88
(nanofluide)), (d) ϕ = 120° , (𝜓𝑚𝑎𝑥 = 26 (fluide pur), 𝜓𝑚𝑎𝑥 = 31,44 (nanofluide)), (fluide pur
(lignes pointillées —), lignes (nanofluides)).

 Fig. 6a–d :
La Fig. 6a–d illustre les lignes de courant et les isothermes à différents angles d'inclinaison
pour Ra = 105 et φ = 0,1. Encore une fois, cette figure compare le nanofluide Cu-eau et le
fluide pur. On observe que la forme de la cellule principale est sensible à l'angle d'inclinaison
et à l'ajout de nanoparticules. De plus, les isothermes indiquent que l'ajout de nanoparticules
devient plus efficace pour ϕ = 90°. Dans ce cas, un écoulement de type Rayleigh–Bénard se
forme à l'intérieur de l'enceinte en raison de l'existence d'un gradient thermique vertical dans
la cavité. L'apparition du panache est supposée être la raison de l'amélioration plus importante
du transfert de chaleur pour ϕ = 90°. Plus de discussions sur ce comportement seront abordées
ultérieurement lors de l'étude de la Fig. 10. La direction de rotation de l'écoulement change
complètement pour ϕ = 120°, passant de dans le sens des aiguilles d'une montre à dans le sens
contraire des aiguilles d'une montre.

24
6) Effets de la Fraction Volumique sur les Champs d'Écoulement et la
Distribution de Température pour Ra =𝟏𝟎𝟓 et ϕ = 30°

Figure 7 Rationalises (à gauche) et isothermes (à droite) pour Ra = 105 et ϕ = 30°, (a) = 0,

(𝜓𝑚𝑖𝑛 = -17,42), b) = 0,05, (𝜓𝑚𝑖𝑛 = -18,49).

 Fig. 7a et b :
Les effets de la fraction volumique sur les champs d'écoulement et la distribution de
température sont affichés dans la Fig. 7a et b pour Ra = 105 et ϕ = 30°. Comme indiqué dans
la littérature (Khanafer et al., 2003), les composantes de vitesse du nanofluide augmentent en
raison d'une augmentation du transport d'énergie dans le fluide avec l'augmentation de la
fraction volumique. Ainsi, la dimension de la cellule principale devient plus petite avec
l'augmentation de la fraction volumique du nanofluide. Les valeurs absolues des fonctions de
courant indiquent que la force de l'écoulement augmente avec l'augmentation de la fraction
volumique du nanofluide.

25
7) Effets des Angles d'Inclinaison et des Concentrations de Nanoparticules
sur la Variation du Nombre de Nusselt par Rapport au Nombre de
Rayleigh

Figure 8 Variation des nombres de Nusselt moyens avec les nombres de Rayleigh à
différentes valeurs de fraction volumique pour différents angles d'inclinaison (a) ϕ = 0° ,
(b) ϕ = 30° , (c) ϕ = 60° , (d) ϕ = 90° ,(e) ϕ = 120°

26
  Fig. 8a-e :
La Fig. 8a–e présente la variation du nombre de Nusselt par rapport au nombre de Rayleigh
pour différents angles d'inclinaison. Trois concentrations différentes de particules sont
utilisées, φ = 0, 0,05 et 0,1. La figure montre que l'effet de la concentration de nanoparticules
sur le nombre de Nusselt est plus prononcé à bas nombre de Rayleigh qu'à haut nombre de
Rayleigh. Cela est lié au mécanisme dominé par la conduction pour le transfert de chaleur à
bas nombre de Rayleigh par rapport au mécanisme de convection à nombre de Rayleigh élevé.
Par conséquent, l'effet des nanoparticules hautement conductrices sur le transfert de chaleur
devient plus significatif à bas nombre de Rayleigh. Cependant, les forces de flottabilité
augmentent et surpassent les forces visqueuses et le transfert de chaleur est dominé par la
convection à haut nombre de Rayleigh. De plus, il est intéressant de noter qu'à une
concentration élevée de nanoparticules, le modèle de viscosité de Brinkman rend le fluide
plus visqueux, ce qui explique la différence relativement plus faible entre le nombre de
Nusselt pour les cas de Ra =1 × 104 et Ra = 1 × 105 par rapport à la différence entre Ra =
1 × 104 et Ra = 1 × 103 . L'angle d'inclinaison affecte le transfert de chaleur. Globalement, le
transfert de chaleur diminue avec l'angle d'inclinaison et Le transfert de chaleur est formé
pour ϕ = 90°.

8) Variation du Nombre de Nusselt Moyen en Fonction du Nombre de

Rayleigh pour Différents Angles d'Inclinaison

Figure 9 Variation des nombres de Nusselt moyens avec des angles d'inclinaison pour φ= 0,1.

 Figure 9 :
La figure 9 montre la variation du nombre de Nusselt moyen en fonction du nombre de
Rayleigh pour différents angles d'inclinaison. Le transfert de chaleur augmente avec
l'augmentation du nombre de Rayleigh de manière monotone. Il est démontré que le cas de
l'angle d'inclinaison de 90° au plus bas nombre de Nusselt en raison de la faible vitesse
rencontrée à cet angle ; voir la figure 10, où la vitesse verticale est tracée au milieu du mur
chauffé (y = 0,5) pour φ = 0,1 et Ra = 105 . De plus, la figure 9 montre que le nombre de
Nusselt maximal se produit pour le cas de l'angle d'inclinaison de 30° en raison de
l'augmentation de la vitesse d'écoulement à cet angle comme le montre la figure 10. De plus,
la figure 9 montre que le régime de conduction prévalait pour Ra = 1000, puis le nombre de
Nusselt moyen a commencé à s'écarter de la valeur de Nu ≈ 1.

27
  Figure 10 :

Figure 10 Profils de vitesse au milieu de la paroi chaude pour différents angles d'inclinaison
pour φ = 0,1 et Ra =105 .

9) Variation du Nombre de Nusselt Local avec un Angle d'Inclinaison de 90°

et l'Impact des Nanoparticules
 Figure 11 :

Figure 11 Variation du nombre de Nusselt le long de la paroi chauffée pour ϕ =90°.

La figure 11 est présentée pour montrer la variation du nombre de Nusselt local le long du
mur chauffé pour le cas de l'angle d'inclinaison de 90°. La figure 11 montre que la valeur
minimale du nombre de Nusselt se situe autour de l'axe central de la plume. Le gradient de
température au niveau du mur chauffé est minimal dans la région de la plume, ce qui entraîne
une valeur minimale du nombre de Nusselt dans cette région. Le gradient de température est
minimal dans la région de la plume car l'épaisseur de la couche limite thermique dans la

28
région de la plume devient très grande. Cela s'accompagne d'une réduction des gradients de
température. De plus, la figure 11 montre que l'ajout de nanoparticules entraîne une
augmentation du nombre de Nusselt.

10) Variation du Nombre de Nusselt Local en Fonction des Angles

d'Inclinaison pour φ = 0,1 et l'Amélioration du Transfert de Chaleur
 Figure 12 :

Figure 12 Variation du nombre de Nusselt local le long de la paroi chauffée pour φ = 0,1
à différents angles d'inclinaison.

La figure 12 illustre la variation du nombre de Nusselt local le long du mur chauffé à

différents angles d'inclinaison pour φ = 0,1. Comme on peut le voir sur la figure, le nombre de
Nusselt local diminue le long de la moitié inférieure du mur chaud pour les valeurs de l'angle
d'inclinaison de ϕ = 0° à ϕ = 60°. Cependant, le nombre de Nusselt augmente sur la partie
supérieure du mur. Une tendance ondulée est observée pour ϕ = 90° en raison de la formation
d'un écoulement de type Rayleigh-Bénard à l'intérieur de l'enceinte. Ainsi, deux points
maximaux ont été formés à cet angle mais leurs valeurs sont inférieures à celles des autres
valeurs des angles d'inclinaison en raison de la valeur réduite de la vitesse comme le montre la
figure 10. Les discussions précédentes indiquent que généralement le transfert de chaleur est
amélioré avec l'ajout de nanoparticules. Pour estimer l'amélioration du transfert de chaleur
entre le cas de φ = 0,1 et le cas de fluide pur (fluide de base), l'amélioration est définie comme
𝑁𝑢(φ = 0.1 )−𝑁𝑢(𝑏𝑎𝑠𝑒 𝑓𝑙𝑢𝑖𝑑𝑒)
𝐸= × 100%
𝑁𝑢(𝑏𝑎𝑠𝑒 𝑓𝑙𝑢𝑖𝑑𝑒)

29
11) Amélioration du Transfert de Chaleur en Fonction de l'Angle
d'Inclinaison et du Nombre de Rayleigh pour φ = 0,1

 Figure 13 :

Figure 13 Rapport d'amélioration du transfert de chaleur dû à l'ajout de nanoparticules.

L'amélioration du transfert de chaleur (φ = 0,1) est tracée en fonction de l'angle d'inclinaison à

différents nombres de Rayleigh dans la figure 11. Pour l'ensemble de la plage de nombres de
Rayleigh, la figure illustre que l'amélioration du transfert de chaleur est presque constante
jusqu'à ϕ = 90°. De plus, le pourcentage d'amélioration du transfert de chaleur diminue avec
l'augmentation du nombre de Rayleigh pour tous les angles d'inclinaison. Il est intéressant de
noter que l'amélioration du transfert de chaleur pour Ra =1 × 104 est la même pour tous les
angles d'inclinaison. De plus, pour Ra = 105 , l'amélioration du transfert de chaleur est réduite
en augmentant l'angle d'inclinaison.

12) Conclusions
Dans cet article, l'influence de l'angle d'inclinaison et de la fraction volumique de
nanoparticules a été étudiée pour une enceinte carrée. Le cuivre-eau a été utilisé comme
nanofluide. Les résultats ont clairement indiqué que l'ajout de nanoparticules de cuivre a
produit une amélioration remarquable du transfert de chaleur par rapport à celui du fluide pur.
Le transfert de chaleur s'améliore presque linéairement avec l'augmentation du nombre de
Rayleigh, mais l'effet de la concentration de nanoparticules sur le nombre de Nusselt est plus
prononcé à faible nombre de Rayleigh qu'à nombre de Rayleigh élevé. L'angle d'inclinaison
de l'enceinte est proposé comme paramètre de contrôle pour l'écoulement du fluide et le
transfert de chaleur. Il a été constaté qu'un transfert de chaleur plus faible se produit pour un
angle de 90 degrés. Cependant, des valeurs plus élevées de fraction volumique deviennent
insignifiantes du point de vue de l'écoulement du fluide à cet angle d'inclinaison. Les effets de
l'angle d'inclinaison sur le pourcentage d'amélioration du transfert de chaleur deviennent
insignifiants à faible nombre de Rayleigh, mais il diminue l'amélioration du transfert de
chaleur avec le nanofluide. Enfin, l'angle d'inclinaison est un bon paramètre de contrôle pour
les enceintes remplies de fluide pur et de nanofluide.

30
13) Références
Abu-Nada, E., 2008. Application des nanofluides pour l'amélioration du transfert de chaleur
des écoulements séparés rencontrés dans une marche arrière. Int. J. Heat Fluid Flow 29, 242–
249.
Akbarinia, A., Behzadmehr, A., 2007. Étude numérique de la convection mixte laminaire d'un
nanofluide dans des tubes courbés horizontaux. Appl. Therm. Eng. 27, 1327–1337.
Aounallah, M., Addad, Y., Benhamadouche, S., Imine, O., Adjlout, L., Laurence, D., 2007.
Investigation numérique de la convection naturelle turbulente dans une cavité carrée inclinée
avec une paroi chaude ondulée. Int. J. Heat Mass Transfer 50, 1683–1693.
Aydin, O., Unal, A., Ayhan, T., 1999. Étude numérique de l'écoulement dû à la flottabilité
dans une enceinte carrée inclinée chauffée et refroidie sur des parois adjacentes. Numer. Heat
Transfer Part A 36, 585–589.
Bairi, A., Laraqi, N., Garcia de Maria, J.M., 2007. Étude numérique et expérimentale de la
convection naturelle dans des cavités parallélépipédiques inclinées pour de grands nombres de
Rayleigh. Exp. Therm. Fluid Sci. 31, 309–324.
Bejan, A., 1995. Convective Heat Transfer, deuxième éd. Wiley, New York.
Brinkman, H.C., 1952. La viscosité des suspensions et solutions concentrées. J. Chem. Phys.
20, 571–581.
Catton, I., 1978. Convection naturelle dans les enceintes. In: Proceedings of the Sixth
International Heat Transfer Conference, 6.
Choi, U.S., 1995. Amélioration de la conductivité thermique des fluides avec des
nanoparticules. In: Siginer, D.A., Wang, H.P. (Éds.), Développements et Applications des
Écoulements Non-Newtoniens, FED, vol. 231, 66, pp. 99–105.
Cianfrini, C., Corcione, M., Dell’Omo, P.P., 2005. Convection naturelle dans des cavités
carrées inclinées avec des parois opposées chauffées différemment. Int. J. Therm. Sci. 44,
441–451.
Daungthongsuk, W., Wongwises, S., 2007. Une revue critique du transfert de chaleur
convectif des nanofluides. Renew. Sust. Eng. Rev. 11, 797–817.
De Vahl Davis, G., Jones, I.P., 1983. Convection naturelle dans une cavité carrée : une
solution numérique de référence. Int. J. Numer. Meth. Fluid 3, 227–248.
Duangthongsuk, W., Wongwises, S., 2008. Effet des modèles de propriétés thermophysiques
sur la prédiction du coefficient de transfert de chaleur convectif pour un nanofluide à faible
concentration. Int. Commun. Heat Mass Transfer 35, 1320–1326.
Elsherbiny, S.M., 1996. Convection libre dans des couches d'air inclinées chauffées par le
dessus. Int. J. Heat Mass Transfer 39, 3925–3930.
Hwang, K.S., Ji-Hwan, L., Jang, S.P., 2007. Transfert de chaleur par flottabilité des
nanofluides à base d'eau-Al2O3 dans une cavité rectangulaire. Int. J. Heat Mass Transfer 50,
4003–4010.

31
Jang, S.P., Choi, S.U.S., 2004. Convection libre dans une cavité rectangulaire (convection de
Bénard) avec des nanofluides. In: Proceedings of the IMECE, Anaheim, Californie, États-
Unis.
Jou, R.Y., Tzeng, S.C., 2006. Recherche numérique sur l'amélioration du transfert de chaleur
convectif naturel rempli de nanofluides dans des enceintes rectangulaires. Int. Commun. Heat
Mass Transfer 33, 727–736.
Kang, H.U., Kim, S.H., Oh, J.M., 2006. Estimation de la conductivité thermique d'un
nanofluide en utilisant un volume efficace de particules expérimental. Exp. Heat Transfer 19,
181–191.
Khalifa Bdul-Jabbar, N., 2001. Coefficient de transfert de chaleur par convection naturelle –
une revue, II. Surfaces verticales et horizontales isolées. Energy Convers. Manage. 42, 505–
517.
Khanafer, K., Vafai, K., Lightstone, M., 2003. Amélioration du transfert de chaleur par
flottabilité dans une enceinte bidimensionnelle utilisant des nanofluides. Int. J. Heat Mass
Transfer 46, 3639–3653.
Krane, R., Jessee, J., 1983. Quelques mesures de champ détaillées pour un écoulement de
convection naturelle dans une enceinte carrée verticale. In: Proceedings of the First ASME–
JSME Thermal Engineering Joint Conference, 1, pp. 323–329.
Lee, T., Lin, T.F., 1995. Convection naturelle tridimensionnelle de l'air dans une cavité
cubique inclinée. Numer. Heat Transfer Part A 27, 681–703.
Lo, D.C., Young, D.L., Murugesan, K., Tsai, C.C., Gou, M.H., 2007. Formulation vitesse-
vorticité pour la convection naturelle en 3D dans une cavité inclinée par la méthode DQ. Int.
J. Heat Mass Transfer 50, 479–491.
Maiga, S.E.B., Palm, S.J., Nguyen, C.T., Roy, G., Galanis, N., 2005. Amélioration du transfert
de chaleur en utilisant des nanofluides dans les écoulements en convection forcée. Int. J. Heat
Fluid Flow 26, 530–546.
Ostrach, S., 1988. Convection naturelle dans des enceintes. J. Heat Transfer 110, 1175–1190.
Oztop, H.F., Abu-Nada, E., 2008. Étude numérique de la convection naturelle dans une
enceinte rectangulaire partiellement chauffée remplie de nanofluides. Int. J. Heat Fluid Flow
29 (5), 1326–1336.
Patankar, S.V., 1980. Transfert de chaleur numérique et écoulement de fluide. Hemisphere
Publishing Corporation, Taylor and Francis Group, New York.
Polidori, G., Fohanno, S., Nguyen, C.T., 2007. Note sur la modélisation du transfert de
chaleur des nanofluides newtoniens en convection libre laminaire. Int. J. Therm. Sci.46, 739–
744.
Putra, N., Roetzel, W., Das, S.K., 2003. Convection naturelle des nanofluides. Heat Mass
Transfer 39, 775–784.
Soong, C.Y., Tzeng, P.Y., Chiang, D.C., Sheu, T.S., 1996. Étude numérique sur la transition de
mode de la convection naturelle dans des enceintes inclinées chauffées différemment. Int. J.
Heat Mass Transfer 39, 2869–2882.

32
Trisaksri, V., Wongwises, S., 2007. Revue critique des caractéristiques de transfert de chaleur
des nanofluides. Renew. Sust. Energy Rev. 11, 512–523.
Varol, Y., Oztop, H.F., Koca, A., Ozgen, F., 2008. Convection naturelle et écoulement de
fluide dans une enceinte inclinée avec un chauffage d'angle. Appl. Therm. Eng. doi:
10.1016/j.applthermaleng.2008.02.033.
Versteeg, H.K., Malalasekera, W., 1995. Introduction à la dynamique des fluides
computationnelle : La méthode des volumes finis. John Wiley & Sons Inc, New York.
Wang, X.-Q., Mujumdar, A.S., 2007. Caractéristiques de transfert de chaleur des nanofluides :
une revue. Int. J. Therm. Sci. 46, 1–19.
Wang, X.-Q., Mujumdar, A.S., Yap, C., 2006. Transfert de chaleur par convection libre dans
des cavités rectangulaires horizontales et verticales remplies de nanofluides. In : Conférence
internationale sur le transfert de chaleur IHTC-13, Sydney, Australie.
Xuan, Y., Li, Q., 2000. Amélioration du transfert de chaleur des nanofluides. Int. J. Heat Fluid
Flow 21, 58–64.

33