DJE5329
DJE5329
DJE5329
Mémoire
Présenté par :
DJEBLI Ahmed
Soutenu le : 05 / 05 / 2009
i
REMERCIEMENTS
Je tiens à remercier toutes les personnes qui, par leur participation et leur
encouragement, m’ont permis de mener à bonne fin mon travail de mémoire.
Je tiens aussi à remercier mes amis Fatouh, Bachir et Mehdi pour leur aide
ii
RÉSUMÉ
iii
ABSTRACT
This study concerns a numerical simulation of a cooperating and opposing bi-diffusive natural
convection in a two-dimensional square enclosure having an aspect ratio (length / height)
equal to 4. The enclosure is filled with a homogeneous and isotropic porous media. The
horizontal walls of the enclosure are adiabatic and impermeable. The vertical walls are
maintained in different temperatures and concentrations constant. Fluid is air and the
substance diffused in the air of the enclosure is the vapor. Boundary Conditions of
temperatures and concentrations involve the natural convection bi-diffusive. The flow in the
porous media in the enclosure is modeled by Darcy-Brinkman-Forchheimer's model. The
numerical method of finite volumes is used for the resolution of the equations of continuity,
momentums, energy and mass. A case of bi-diffusive natural convection in a fluid medium
was study and considered as a reference. Results obtained allow the study of the influence of
certain parameter in the flow and the heat and mass transfer, in cooperating and opposing
cases, in the studied system. A comparison between Darcy model and Darcy-Brinkman-
Forchheimer model is realized to study the Forchheimer-Brinkman terms effects.
Key words: Porous medium, Heat and mass transfer, model of Darcy-Brinkman-Forchheimer,
Natural convection, finite volumes.
iv
HIـــــــــــــــــKLM
_ون وPPURk_ر اﻝPPbRcd^ اe_PPfg ^PPUIhi اﻝjPPklL^ ﻝS_PP`a رOPPQRS TPPLURM VPPWراY\[ اﻝPPه
_دةPkS ءOPLk_ق ﻡPif اﻝ.4 _ويP`_ع یPrMرd اsPLt لOPiن اﻝOPv یxPIlS yPSz_ق ﻡic ^} _دpRkاﻝ
_نPPPدیOkUاران اﻝYPPP اﻝ.\PPP}_c zPPPI_ن وPPPk
__ق آPPPifL_ن ﻝPPPI}اران اYPPP اﻝ.VPPP`c_R\ة ﻡPPP}_c
^P} ةzPbRfk_دة اﻝPk_ء و اﻝPk اﻝOP هye_k اﻝ.VRS_g VrLRK ﻡIآzM ارة وza _ن }^ درﺝ_تItOﺽOﻡ
zّIP`M PIآzRارة واﻝzPlرﺝ_ت اﻝYPود ﻝYPl اﻝYPft وطzPb اﻝ._ءPk_ر اﻝPKS ^P_ق هif اﻝj_ء داﺥkاﻝ
-^PWذج دارOPkfS \P_ق أﺥPif^ اﻝP} \ةP}_f_دة اﻝPk^ اﻝP} _نPیz اﻝ._رPbRcd^ اe_Pfg ^UIhi اﻝjklاﻝ
،VاریzkRPWdت اd_دPU ﻡjPl ﻝLkURPW اVPIRfkم اﻝOPlL ﻝVIkz اﻝVیzi اﻝ.zkIرﺵO}-_نkvczS
RPW دراPkM ye_P_ق ﻡPic ^P} _رPbRcd^ اe_Pfg ^PUIhi اﻝjPkl اﻝ._دةPk واﻝVP_i اﻝ،VPآzl اﻝVIkآ
_ن وPPیz اﻝsPPLt jPPاﻡOt ةYPPt zIgPPM VPPWراYS kPP`M _PPILt VLPPQlk اﻝe_PPRf اﻝ.yPPﺝzk\[ آPPوأﺥ
IPPS VPPc ﻡ_ر.فzPPUk_م اﻝPPf^ اﻝPP} VPPc_وURk_دة واﻝPPpRk اﻝIR_ﻝPPl^ اﻝPP} _دةPPkارة واﻝzPPl_ر اﻝPPbRcا
IیzPQfU اﻝzIgPM VPW دراjP اﺝP\ت ﻡPrc zkIPرﺵO}-_نkvczS-^Wذج دارOkc^ وWذج دارOkc
.zkIرﺵO} _ن وkvczS
،zkIPرﺵO} -_نkvczS -^Wذج دارOkc ،_ديkاري واﻝzl_ل اﻝRcd ا،\}_c W و: ت ا "ا ــ$%&آ
.VIRfkاﻝ مOl اﻝ،^UIhi اﻝjklاﻝ
v
TABLE DES MATIERES
DÉDICACE…………………………………………………………………………………… i
REMERCIEMENTS…………………………………………………………………………. ii
RÉSUMÉS…………………………………………………………………………………… iii
TABLES DE MATIERES…………………………………………………………………… vi
NOMENCLATURE…………………………………………………………………………. ix
1. INTRODUCTION ................................................................................................................. 1
1.1 GÉNÉRALITÉS .................................................................................................................. 1
1.2 CONTENU DU MÉMOIRE ............................................................................................... 2
1.3 REVUE BIBLIOGRAPHIQUE .......................................................................................... 3
2. MODÉLISATION PHYSIQUE ........................................................................................... 18
2.1 DESCRIPTION DU PROBLÈME…………………………………..…………………..18
2.2 HYPOTHÈSES SIMPLIFICATRICES ............................................................................ 19
2.3 MODÉLISATION MATHÉMATIQUE ........................................................................... 19
2.3.1 Equation de continuité..................................................................................................... 20
2.3.2 Equations de quantité de mouvement dans les milieux poreux ...................................... 20
2.3.2.1 L’équation de la quantité de mouvement suivant x .................................................... 20
2.3.2.2 L’équation de quantité de mouvement suivant y........................................................ 20
2.3.3 Equation de l’énergie ...................................................................................................... 20
2.3.4 Equation de transfert de matière ..................................................................................... 21
2.4 MISE EN FORME ADIMENSIONNELLE DES ÉQUATIONS ..................................... 21
2.4.1 Variables adimensionnelles ............................................................................................. 21
2.4.2 Equations adimensionnelles ............................................................................................ 22
2.4.2.1 Equation de continuité ................................................................................................ 22
2.4.2.2 Equation de la quantité de mouvement suivant x ....................................................... 22
2.4.2.3 Equation de quantité de mouvement suivant y........................................................... 22
2.4.2.4 Equation de l’énergie.................................................................................................. 22
2.4.2.5 Equation de transfert de matière ................................................................................. 22
2.5 CONDITIONS INITIALES ET AUX LIMITES .............................................................. 23
2.5.1 Conditions initiales ......................................................................................................... 23
2.5.2 Conditions aux limites .................................................................................................... 23
2.5.2.1 Cas coopérant ............................................................................................................. 23
vi
2.5.2.2 Cas opposant............................................................................................................... 23
2.6 TRANSFERTS THERMIQUE ET MASSIQUE .............................................................. 24
3. RÉSOLUTION NUMÉRIQUE ............................................................................................ 25
3.1 INTRODUCTION .............................................................................................................. 25
3.2 PRÉSENTATION GÉNÉRALE ........................................................................................ 25
3.2.1 Différences finies ............................................................................................................. 25
3.2.2 Eléments finis .................................................................................................................. 26
3.2.3 Volumes finis................................................................................................................... 26
3.3 MÉTHODE DES VOLUMES FINIS ................................................................................. 27
3.3.1 Introduction ..................................................................................................................... 27
3.3.2 Présentation générale de la méthode ............................................................................... 27
3.4 FORME GÉNÉRALE DE L’ÉQUATION DE TRANSPORT .......................................... 30
3.5 DISCRÉTISATION ........................................................................................................... 30
3.5.1 Intégration l’équation générale de transport .................................................................... 30
3.5.2 Différentes schémas de discrétisation ............................................................................. 32
3.5.3 Équation de continuité ..................................................................................................... 34
3.5.4 Équations de quantité de mouvement .............................................................................. 34
3.5.5 Équation de l’énergie ....................................................................................................... 36
3.5.6 Équation de transfert de matière ...................................................................................... 37
3.5.7 Résolution des systèmes d’équations de discrétisation UPu et VPv .................................. 37
3.6 MÉTHODE DE RÉSOLUTION ........................................................................................ 40
3.7 SOLUTION DES SYSTÈMES D’ÉQUATIONS DISCRÉTISÉES .................................. 42
3.8 CRITÈRE DE CONVERGENCE ...................................................................................... 42
3.9 DÉTAILS NUMÉRIQUES ................................................................................................ 43
4. RÉSULTATS ET DISCUSSIONS………………………………………………...……... 44
4.1 EFFET DE MAILLAGE ................................................................................................... 44
4.2 VALIDATION DU CODE DE CALCUL ........................................................................ 46
4.3 RÉSULTATS ET DISCUSSIONS .................................................................................... 49
4.3.1 Cas coopérant : milieu fluide avec Grc=105 .................................................................... 49
4.3.2 Modèle de Darcy-Brinkman-Forchheimer ...................................................................... 59
4.3.2.1 Première partie : cas coopérant ................................................................................. 59
4.3.2.1.1 Cas de Grc=105 ......................................................................................................... 59
4.3.2.1.1.1 Structure d’écoulement et distributions de la pression, de la température et de la
concentration ............................................................................................................................ 59
vii
4.3.2.1.1.2 Variation temporelle des nombres de Nusselt et de Sherwood ............................. 68
4.3.2.1.2 Cas de Grc=5.105 ...................................................................................................... 70
4.3.2.1.2.1 Structure d’écoulement et distributions de la pression, de la température et de la
concentration ............................................................................................................................ 70
4.3.2.1.2.2 La variation temporelle des nombres de Nusselt et de Sherwood ......................... 79
4.3.2.1.3 Cas de Grc=106 ......................................................................................................... 81
4.3.2.1.3.1 Structure d’écoulement et distributions de la pression, de la température et de la
concentration ............................................................................................................................ 81
4.3.2.1.3.2 La variation temporelle des nombres de Nusselt et de Sherwood ......................... 90
4.3.2.1.4 Cas de Da=5.10-7 et Grc=105 .................................................................................... 92
4.3.2.1.4.1 Structure d’écoulement et distributions de la pression, de la température et de la
concentration ............................................................................................................................ 92
4.3.2.1.4.2 La variation temporelle des nombres de Nusselt et de Sherwood ....................... 100
4.3.2.2 Deuxième partie : cas opposant, Da=10-2 ................................................................ 102
4.3.2.2.1 Cas de Grc=105 ....................................................................................................... 102
4.3.2.2.2 Cas de Grc=5.105 .................................................................................................... 107
4.3.3 Modèle de Darcy, cas coopérant ................................................................................... 113
4.3.3.1 Cas de Da=10-2 ......................................................................................................... 113
4.3.3.1.1 Cas de Grc =105 ...................................................................................................... 113
4.3.3.1.2 Cas de Grc=5.105 .................................................................................................... 122
4.3.3.2 Cas de Da=5.10-7 et Grc=105 .................................................................................... 129
4.3.4 Comparaisons des deux modèles Darcy-Brinkman-Forchheimer et Darcy .................. 137
CONCLUSION GÉNÉRALE ................................................................................................ 142
RÉFÉRENCES BIBLIOGRAPHIQUES ............................................................................... 144
viii
NOMENCLATURE
Symboles Grecs
ix
βc coefficient d’expansion solutal [g.mol.l-1]
βt coefficient d’expansion thermique [K-1]
µ viscosité dynamique [kg.m-1.s-1]
ν viscosité cinématiques [m2.s-1]
ρ masse volumique [kg.m-3]
variable générale dépendante du temps
coefficient de diffusion
∆ écart de la concentration
[g.mol.l-1]
∆ écart de la temperature
[K]
∆ incrément de temps adimensionnel
∆ longueur adimensionnelle d’un volume de contrôle dans la direction X
∆ longueur adimensionnelle d’un volume de contrôle dans la direction Y
Nombres adimensionnels
nombre de Darcy
Da
Sh nombre de Sherwood
nombre de Sherwood moyen
nombre de Lewis
!
"
Le
Nu nombre de Nusselt
&' nombre de Nusselt moyen
nombre de Prandtl
(
!
Pr
nombre de Schmidt
(
"
Sc
x
Exposant
n relatif au temps t
n+1 relatif au temps t + ∆t
Indices
xi
Chapitre 1 Introduction
1. INTRODUCTION
1.1 GÉNÉRALITÉS
Le choix de ce thème est justifié par l’importance fondamentale que représentent les
interactions des phénomènes de transfert de chaleur et de masse couplés, et par l'importance
pratique de ses applications. On peut citer comme exemple le transport de polluants dans les
sols, le stockage des déchets radioactifs, la migration de l'humidité dans des fibres isolantes,
la métallurgie, l’environnement, etc. Dans la plupart de ces applications, les fluides impliqués
sont non newtoniens. En effet, ce sont des fluides à structure moléculaire complexe, qui
n’obéissent pas à la loi classique de Newton. Des exemples de tels fluides sont les solutions
de polymères, les huiles, les pâtes à papier et les détergents liquides (Mamou 1998).
Ce thème ayant fait l’objet, dans le passé, de nombreuses recherches à cause des
applications potentielles possibles tant dans la nature que dans l’industrie. Une revue
exhaustive des publications disponibles sur ce sujet a été faite par Nield et Bejan (2006) pour
les milieux poreux.
1
Chapitre 1 Introduction
Dans le quatrième chapitre, une validation du code de calcul a été effectuée par une
comparaison avec les études disponibles dans la littérature. Les résultats obtenus concernent
en premier le cas d’un milieu fluide et ensuite celui d’un milieu poreux saturé d’un fluide
binaire. Pour ce dernier, les résultats sont divisés en trois parties. Le modèle de Darcy-
Brinkman-Forchheimer été employé dans la première partie où les forces de volumes sont
coopérantes et la deuxième partie où elles sont opposantes. La troisième a été destinée à
l’étude du modèle de Darcy pour déterminer l’effet des termes de Brinkman et Forchheimer.
Une comparaison des résultats obtenus dans ces trois cas sera faite à la fin de ce chapitre.
Le mémoire est terminé par une conclusion générale sur les résultats obtenus et
quelques recommandations pour les études futures.
2
Chapitre 1 Introduction
Une revue de littérature montre qu’il existe un grand nombre de travaux numériques
ou théoriques consacrés à l’étude de la convection naturelle à double diffusive dans une cavité
poreuse rectangulaire. De nombreuses références proposent des revues très documentées
(Nield et Bejan 2006) et nous nous limiterons à présenter les grandes catégories
d’investigations qui concernent ce sujet très vaste.
Un code précis en différences finies a été réalisé par Mota et Saatdjian (1997) et
employé pour résoudre les équations bidimensionnelles du modèle Darcy-Boussinesq pour un
annulaire excentrique, horizontal et rempli d'un milieu poreux saturé. Les surfaces des
cylindres sont imperméables et maintenues à des températures constantes uniformes. Des
solutions stationnaires ont été obtenues par la méthode ADI (alternating direction implicit
method). Une maille contenant approximativement 15103 points de maillage a été employée
pour capturer les cellules de contra-rotation secondaires dans le sommet de la couche poreuse.
Les résultats obtenus pour diverses valeurs de l’excentricité ε (distance entre les centre de
cylindre/rayon extérieur - rayon intérieur), variant de 0.01 à 0.9, montrent que la réduction du
rapport de rayon ou l'augmentation de l'excentricité a le même impact de géométrie dans la
partie supérieure de la couche où les effets convectifs sont plus prononcés. Pour ε<0.5, les
deux modèles d’écoulements bicellulaire et quadri-cellulaire restent stables pour la gamme du
nombre de Rayleigh étudiée. Pour ε ≥0.5. La transition d'un régime d’écoulement à un autre
se produit par une bifurcation parfaite, c'est-à-dire le changement d’écoulement arrive avec
une des branches de solution qui perd de stabilité. Les auteurs suggèrent que dans une
situation réelle, l'isolation soit plus efficace si ε est mise à la valeur maximale, pour laquelle le
régime de flux à quatre cellules est physiquement réalisable, plutôt qu'à la valeur qui réduit au
minimum le transfert de chaleur quand le modèle de flux est bicellulaire.
Beji et al. (1999) ont réalisé une étude numérique sur le phénomène de la convection
thermosolutale dans un annulaire concentrique verticale remplie d’une matière poreuse. Les
gradients de la température et de la concentration qui influent dans la dynamique du fluide
sont imposés horizontalement et les deux effets de buoyancy peuvent l’augmenter ou le
diminuer. La méthode des différences finis a été utilisée pour la discrétisation. Les résultats
obtenus donnent l’influence du nombre de Rayleigh thermique Rat, le rapport de buoyancy N,
le nombre de Lewis Le et le rapport des diamètres κ (diamètre extérieur/diamètre intérieur)
sur les nombres de Nusselt et de Sherwood.
3
Chapitre 1 Introduction
Une étude numérique détaillée de la convection forcée laminaire, dans une conduite
carrée tridimensionnelle rempli d'un milieu granulaire isotropique et saturé avec un liquide
newtonien, a été réalisée par Chen et Hadim (1999). Les résultats hydrodynamiques et ceux
du transfert de chaleur sont obtenus pour trois conditions thermiques aux limites différentes.
L’écoulement dans le milieu poreux a été modelé employant le modèle Darcy-Brinkman-
Forchheimer. Des modèles empiriques pour la porosité variable et la dispersion thermique ont
été déterminés basés sur des mesures expérimentales tridimensionnelles existantes. Des études
paramétriques ont été alors conduites pour examiner les effets de diamètre de particule, le
nombre de Reynolds, le nombre de Prandtl et le rapport de conductivité thermique. Il a été
trouvé que le nombre Reynolds et le diamètre des particules affectent principalement le
gradient des vitesses. L’augmentation du nombre de Reynolds ou la diminution du diamètre
des particules (c'est-à-dire, quand l'inertie et les effets de dispersion thermaux sont
augmentés), augmente le nombre de Nusselt. Les effets de la variation du rapport d'aspect de
la conduite, sur le facteur de friction, peuvent être négligés pour les petits diamètres de la
particule (Dp≤0.01), ou pour les nombres de Reynolds élevés (Re≥1000).
Mharzi et al. (1999) ont étudié la convection naturelle thermosolutale dans une cavité
allongée suivant l'axe horizontal, divisée par une couche verticale poreuse. Ils ont employé
l'approximation de Boussinesq, la fonction de courant, la vorticity et la formulation Darcy-
Brinkman pour décrire le transfert dans la couche poreuse. Les paramètres dimensionnels
choisis permettent d'obtenir des équations de conservation, valables dans les deux
compartiments de fluide et dans la couche poreuse. Ils ont résolu numériquement les
équations couplées avec l'approche des volumes de contrôle. Après le recouvrement de
quelques résultats de littérature, ils étudient, dans une cavité de section carrée qui est divisée
par une couche poreuse en deux compartiments égaux remplies de fluide. Les résultats
montrent que Le transfert de chaleur est sous l'influence du rapport de conductibilité
thermique Rk, tandis que le transport solutal est essentiellement sensible au rapport de la
diffusivité solutale Rd. L’augmentation de Rk (diminution de Rd) augmente l'échange
thermique par conduction dans la couche poreuse (l'échange diffusif solutal). En outre, leurs
effets globaux sur la convection thermosolutale sont inférieurs à ceux de Ra et Da.
L'augmentation du nombre de Lewis, qui caractérise le rapport de la diffusion thermique et la
diffusivité moléculaire, augmente successivement l'échange solutal moyen (Shmoy), tandis que
le nombre de Nusselt est exposé à une petite diminution. Cependant, quand les forces de
buoyancy thermiques et solutale se sont opposées (N < 0), l'augmentation de Le produit une
4
Chapitre 1 Introduction
Akrour et al. (2000) ont étudié numériquement la convection à double diffusion dans
une enceinte horizontale. La cavité soumise à un gradient de concentration verticale est
chauffée différentiellement au niveau de ses parois horizontales ou verticales. Les équations
couplées modélisantes sont discrétisées numériquement par la méthode des volumes finis. La
résolution a ensuite été effectuée par la procédure de double balayage (ADI) et le couplage
pression-vitesse par l’algorithme SIMPLER. Ils ont considéré l’effet, sur les transferts, du
nombre N qui caractérise le rapport des forces de volume d’origine thermique et solutale pour
Pr=7, Le=100 et A=2. Leurs résultats montrent que pour les faibles valeurs de N, des cellules
convectives dues aux forces thermiques se mettent en place et le transfert devient indépendant
de N. Pour des valeurs de N plus élevées, le transfert est diffusif et la stratification stable en
concentration supprime l’écoulement. Pour des valeurs intermédiaires de N, le transfert
diminue avec l’accroissement de N. Dans la zone de transition, la solution numérique dépend
de la condition initiale.
Une étude numérique a été effectuée sur la convection naturelle à double diffusive
dans un annulaire vertical circulaire poreux par Bennacer et al. (2000). Des mouvements
conduits par la constante de température extérieurement appliquée et des différences de
concentration imposées à travers les murs verticaux de la clôture. La formulation de Darcy-
Brinkman était utilisée. Le flux est supposé être laminaire et bidimensionnel. La variation de
la densité est traitée selon l'approximation de Boussinesq. L'approche des volumes de contrôle
est employée pour résoudre les équations gouvernantes. Pour accélérer la convergence
l'algorithme SIMPLER, à l'origine développé par Patankar, est couplé avec l'algorithme
SIMPLEC de Van Doormaal et Raithby. Des solutions pour les champs de flux, les
distributions de la concentration et des nombres de Nusselt et de Sherwood, sont obtenues
5
Chapitre 1 Introduction
dans les termes des paramètres gouvernants du problème. L'effet du nombre de Darcy Da et la
proportion de rayon κ sur Nui et Shi sont significatifs. Les résultats indiquent clairement que
pour l’écoulement coopérant (N > 0) Sh augmente significativement avec le rapport de
buoyancy N. D'autre part, pour un écoulement opposé (N < 0), les résultats sont beaucoup
plus complexes. Ainsi, les auteurs ont observé d'abord que Sh diminue avec N, et passe par
une valeur Shmin à une valeur donnée Nmin. Ensuite, Sh augmente avec N. On voit que Nmin,
dépend fortement du nombre de Darcy. Les résultats obtenus pour Da = 102 pour κ = 2 et N =
20, correspondent à un cas d’un liquide pur visqueux (Da→∞) qui a été examiné. Ce rapport
de buoyancy (N = 20) donne une structure d’écoulement multicellulaire.
Le début de convection double diffusive dans une cavité horizontale poreuse a été
étudié numériquement employant l'analyse de la stabilité linéaire par Mahidjiba et al. (2000).
Dans la formulation du problème, l'étude est faite selon le modèle Darcy avec l'approximation
de Boussinesq. Les conditions aux limites pour la température et la concentration sont
spécifiées sur les parois horizontales de la cavité tandis que les deux parois verticales sont
imperméables et adiabatique. Les méthodes de Galerkin et des éléments finis sont employées
pour résoudre les équations de perturbation. Le début de convection a été trouvé dépendant du
rapport d'aspect de la cavité A, la porosité φ, le nombre de Lewis Le, rapport des forces de
buoyancy thermique et solutale N et les conditions aux limites thermique et solutale. Pour une
6
Chapitre 1 Introduction
cavité limitée. Les auteurs ont montré l’existence d’un nombre de Rayleigh supercritique.
Pour le début, la convection supercritique est au dessus de la valeur critique de nombre de
Rayleigh pour que l’écoulement soit stable. En outre, les auteurs ont montré que le régime
d’écoulement stable est atteint lorsqu’on trouve un nombre de Rayleigh critique
correspondant à la transition du mode oscillatoire au mode de convection directe. Cependant,
pour une couche infinie les résultats indiquent l'absence d'un régime d’écoulement stable.
Khanafer et Vafai (2002) ont présenté une étude numérique de la convection mixte et
du transport de masse dans une enceinte rectangulaire avec couvert mouvant remplie d’une
matière poreuse saturée. Les deux parois horizontales sont constantes mais avec des
températures et des concentrations différentes. Par contre les deux parois verticales ont un
mouvement et la surface se déplace à une vitesse constante. L’équation de transport est traitée
avec la méthode des volumes finis avec l’ADI comme procédure de résolution. Les résultats
indiquent que le rapport de buoyancy, le nombre de Darcy, le nombre de Lewis et le nombre
de Richardson ont un grand effet sur le phénomène de double diffusion.
Bennacer et al. (2002) ont présenté une étude numérique et analytique concernant le
transfert combiné de chaleur et de masse dans un milieu poreux. Ce milieu est globalement
homogène et présente une anisotropie thermique. L'équation qui gouverne l’écoulement est
celle de Darcy-Brinkman. Le système des équations couplées est résolu par la méthode
classique des volumes finis. Le schéma hybride est utilisé pour la discrétisation spatiale des
équations. Le système des équations linéarisées est résolu par un double balayage (TDMA).
Le couplage pression-vitesse est assuré par une procédure itérative par l’accouplement de
deux algorithmes, SIMPLER et SIMPLEC. Dans le cas d'écoulements d'origine thermique,
l'analyse d'échelle est appliquée pour prédire analytiquement l'évolution du transfert de
chaleur et de masse en fonction de l'anisotropie thermique. Les simulations numériques sont
présentées pour une cavité carrée en faisant varier une large gamme de paramètres.
L’influence des différents paramètres caractérisant ce genre de problème a été explorée
numériquement en faisant varier le taux d’anisotropie thermique, le nombre de Darcy, le
nombre de Lewis, le rapport des forces de volume, le nombre de Rayleigh thermique et le
rapport des conductivités. Ces résultats numériques sont analysés en termes de transfert
moyen de chaleur et de masse sur les parois verticales de la cavité et montrent l’existence
d’un maximum de transfert de masse pour un rapport d'anisotropie critique. Cette situation
critique dépend de plusieurs paramètres dont le nombre de Lewis du fluide saturant. Une
corrélation globale a été établie permettant de prédire le transfert de masse pour les milieux
7
Chapitre 1 Introduction
thermiquement anisotropes en régime de Darcy. Par ailleurs les auteurs ont découplé dans
cette étude les effets de l'anisotropie thermique et hydraulique afin de pouvoir analyser plus
précisément l’influence de chaque anisotropie.
Younsi et al. (2002) ont étudié numériquement l’écoulement bidiffusive opposé dans
une cavité poreuse. Les gradients de température et de concentration sont imposés
horizontalement, et les deux effets de buoyancy peuvent augmenter la dynamique du fluide.
L’équation de Darcy avec les termes de Brinkman-Forchheimer, qui prend en considération
les effets d’inertie et les effets visqueux, a été utilisée pour déduire l’équation de mouvement.
L’algorithme SIMPLER est utilisé en se basant sur la méthode des volumes finis pour la
résolution du couplage de vitesse-pression. Les résultats montrent que le rôle du milieu
poreux est de diminuer la vitesse verticale dans la couche limite. Le transfert de chaleur et de
masse et le champ d’écoulement sont affectés profondément par les forces de buoyancy
8
Chapitre 1 Introduction
Bennacer et al. (2003) ont fait une analyse numérique sur la convection naturelle
thermosolutale dans une cavité remplie de liquide où est insérées dedans des couches
isotropiques ou anisotropiques poreuses. Des couches identiques poreuses sont attachées aux
parois verticales, lesquelles sont tenues à des températures et des concentrations constantes.
Les parois horizontales de la cavité sont adiabatiques et imperméables. Une approche des
volumes de contrôle finis est utilisée dans ce travail. L’algorithme SIMPLEC est employé
pour le couplage de vitesse-pression dans l'équation de mouvement. Le rapport d'aspect de la
cavité est égale à deux et le liquide de saturation est l’air (Pr = 0.71). L'analyse est exécutée
pour des nombres de Grashof thermique de 106 et 107 et de Schmidt de 7.1 pour différents
nombres de Darcy, d’épaisseur de couche poreuse et du rapport de perméabilité. Les résultats
sont présentés pour un flux thermique dominant, N = 0 et pour un flux solutale dominant, N =
10. L'effet d'anisotropie hydraulique sur le taux de chaleur et le transfert de masse est discuté.
Il a été montré que le taux des transferts de chaleur et de masse a une faible dépendance du
nombre de Darcy pour les régimes de perméabilité hauts et bas. Pour une certaine gamme de
paramètres, le taux de chaleur transféré diminue quand le flux pénètre dans la couche poreuse.
De là, il y a une valeur optimum (minimum) du nombre de Nusselt, qui est fonction du
paramètre anisotropie. Une corrélation pour la les transferts de chaleur et de masse a été
présentée.
Bahloul et al. (2004) ont effectué une étude analytique et numérique de la séparation
des composants d'un mélange binaire dans un milieu annulaire circulaire vertical poreux. Les
flux uniformes de chaleur sont appliqués sur les parois verticales de l’enceinte. La
formulation du problème est faite avec le modèle de Darcy. L'approximation Boussinesq est
appliquée à la variation de densité. Le problème est résolu analytiquement sur la base de
l’écoulement parallèle, valable dans la limite d’une cavité mince. La solution du champ
d’écoulement et de la distribution de la concentration dans l’enceinte a été trouvée employant
l’approche des volumes de contrôle. Les équations d'énergie et de concentration sont résolues
par la méthode implicite aux directions alternées (ADI). Le champ de fonction de courant est
obtenu par la méthode de surrelaxation (SOR). Les résultats illustrent clairement le fait que
pour une grande cavité allongée (A >> 1), le flux dans la région fondamentale de l’enceinte
est essentiellement unicellulaire et parallèle tandis que la température et la concentration dans
9
Chapitre 1 Introduction
10
Chapitre 1 Introduction
Barletta et Lazzari (2005) ont étudié numériquement la convection naturelle dans une
cavité verticale carrée remplie d'un liquide saturé comme un milieu poreux au moyen du
logiciel Femlab3.1i. La cavité possède des parois isothermes et elle est chauffée par une
limite interne concentrique circulaire soumise à un flux de chaleur uniforme. Le flux est
supposé être bidimensionnel, stable et laminaire. L’effet de dissipation visqueuse est négligé
et aucune génération de source interne n'est considérée. Le milieu poreux satisfait la loi Darcy
et la validité de l'approximation Boussinesq est utilisée. Le champ de vitesse et le champ de
Gobin et al. (2005a) ont étudié la convection bidiffusive dans un fluide binaire. La
géométrie utilisée dans ce travail est une cavité fermée rempli partiellement avec une couche
verticale poreuse. La description mathématique du problème est basée sur la formulation des
équations de conservations dans un seul domaine. Les résultats numériques montrent
l’influence de la couche poreuse sur la structure d’écoulement et le transfert dans la cavité.
Cette étude avait pour but l’analyse de l’influence des paramètres gouvernants qui
caractérisent la convection à double diffusion. Les transferts de masse et de chaleur sont
analysés en fonction de la perméabilité de la couche poreuse avec les forces de buoyancy, qui
donnent une amélioration spécifique de la structure d’écoulement et du transfert de chaleur
dans la cavité.
Gobin et al. (2005b) ont étudié la convection naturelle dans un milieu partiellement
poreux. Premièrement, ils ont traité les problèmes de modélisation d’écoulement des fluides et
de transfert de chaleur et de masse dans différents domaines, et dans divers modèles
alternatifs proposés dans la littérature. Ainsi, les auteurs ont discuté l’étude de l’écoulement
de fluide et des transferts de chaleur et de masse dus à la double diffusion dans une enceinte
verticale. L’influence de la couche poreuse et le problème de la stabilité linéaire de la
convection naturelle thermique dans la couche horizontale ont été analysés. Les résultats
montrent que la diminution du nombre de Darcy fait diminuer la perméabilité et l’écoulement
devient très difficile.
11
Chapitre 1 Introduction
12
Chapitre 1 Introduction
étude a été obtenu telle que la prévision de la structure d’écoulement et des transferts de
chaleur et de masse pour un grand choix des paramètres gouvernants soit précise. Des
expressions utiles approximatives ont été tirées pour décrire le régime des couches limites.
Des équations simplifiées pour Nu et Sh ont été obtenues pour les divers régimes
d’écoulement (chaleur dominant ou concentration dominante). Pour le cas de l’écoulement
opposé, l'existence des solutions multiples, pour des valeurs données des paramètres
gouvernants, a été démontrée analytiquement et numériquement.
13
Chapitre 1 Introduction
14
Chapitre 1 Introduction
15
Chapitre 1 Introduction
Récemment, Ouazaa (2008) a étudié le cas d’une cavité allongée soumise à des flux
uniformes de chaleur et de masse, il a utilisé l’hypothèse d’un écoulement parallèle pour
obtenir une solution analytique décrivant les champs de fonction de courant et de température
16
Chapitre 1 Introduction
dans la région centrale de la cavité. Il a démontré ainsi que des solutions multiples sont
possibles, certaines d’entre elles étant instables. La méthode des volumes finis a été employée
pour obtenir des solutions numériques à partir des équations gouvernantes complètes. Les
résultats obtenus ont montré un excellent accord entre la solution analytique et la simulation
numérique. Dans cette étude l'effet des paramètres de contrôle sur l'écoulement et sur les
transferts de chaleur et de masse dans le système a été déterminé.
Aussi, une autre étude numérique a été réalisée par Sivasankaran et al. (2008) pour
analyser l'effet combiné des gradients de température et de concentration sur l’écoulement de
la convection naturelle de l’eau froid près de sa densité extrême contenu dans une cavité
poreuse. Les équations gouvernantes sont discrétisées avec la méthode des volumes finis. Les
résultats sont discutés pour différentes valeurs de la porosité, du nombre de Darcy et du
nombre de Grashof. Les taux des transferts de chaleur et de masse calculés ne sont pas
linéaires avec la température de la paroi chaude. Ce comportement non linéaire des taux de
transfert de chaleur et de masse est dû à l'effet de la densité maximal. Le transfert de chaleur
et de masse augmente avec l'accroissement du nombre de Darcy et la porosité. Il a été observé
que la température à la densité maximale a des effets importants sur le transfert de chaleur et
de masse, dus à la formation d’une structure bicellulaire. Les résultats montrent aussi que les
plus hauts taux de transfert de chaleur et de masse peuvent être obtenus dans un milieu poreux
saturé avec l’eau en présence de porosité élevée.
17
Chapitre 2 Modélisation physique et mathématique
2. MODÉLISATION PHYSIQUE
ET MATHÉMATIQUE
Y, V Paroi adiabatique
et impermeable
Paroi active
T T0
H Milieu poreux
C0 C
Paroi active X, U
Paroi adiabatique
et imperméable
L
18
Chapitre 2 Modélisation physique et mathématique
1. Le milieu poreux est considéré homogène, isotrope et ne subit pas de variation spatiale
de la porosité dans l’enceinte.
2. Le fluide qui s’écoule est incompressible.
3. Le problème est bidimensionnel.
4. L’écoulement engendré est laminaire et transitoire.
5. Il n’y a ni réaction chimique ni source de chaleur ou de masse.
6. Les interactions entre les transferts de chaleur et de masse, connues sous le nom
d’effets Soret et Dufour respectivement, sont négligeables.
7. Le transfert de chaleur par rayonnement est négligeable.
8. Dans la convection naturelle bidiffusive, le couplage de transfert de chaleur et de
matière est nécessaire, parce que la densité de mélange du fluide, dépend de la
température et de la concentration. Les propriétés thermo-physiques du fluide sont
constantes et sont évaluées à la température et à la concentration de référence. On
considère approximativement que ρ = ρ0 est constante dans tous les termes, sauf celui
qui contient l’accélération de la pesanteur g (Fy = – ρf g) où ρ = f (T, C), (ρ est
fonction de la température et de la concentration).
On a alors : ρ = ρ0[1– βT (T–T0) – βc(C–C0)]
Dans notre cas : ρf (T, C) = ρf (T0, C0) [1–βT (T–T0) – βc(C–C0)].
19
Chapitre 2 Modélisation physique et mathématique
*
+,-
12 345 6
/0 12 7 0
+.
(2.1)
Fluide incompressible :
/ 0
+9 +:
+ +;
(2.2)
12 3V V7B
= 102G* IH / K 0J 6
12
K ?
-O,- Q6
12 Q6
12 / 4T
12 112 112 12
45 < / > 0 5S
2
= +? =
CDDDDDEDDDDDF
> +. A CDEDF ,CDEDF
> - > M CD
DED
P DF
(2.3)
c. (3) ,
K
12 : Terme visqueux de Brinkman.
0J 6
->
12 : Terme de Darcy.
d.(4) 6
K
e. (5)
- , -
O
12 Q 6
Q6 12 : Terme Forchheimer (forces inertielles microscopiques).
/ > / >
+
'
5 ' / > V+ [+ \ / +; [];\W
= +9 = +G99H = +G:9H +X K √9 Z: K + +9 + +9
> +. + +; √
(2.4)
/ /
^
5 ^ / > V+ [+ \ / +; [+;\W /
= +: = +G9:H = +G::H +X K √9 Z: K + +: + +:
> +. > + > +; +; √
45 S _` G
a H / ` G
a Hb (2.5)
G4cH
+
12 102 0
/ G4cH5 6 12 3d 102 7 / e
fff
+.
(2.6)
Où : G4cH G1
*HG4cHh / *G4cH5
20
Chapitre 2 Modélisation physique et mathématique
d G1
*Hdh / *d5
e
fff
G1
*Heifff / *e5fff
3,jk 7
On divise l’équation par 34cX 75 et on considère que : 3,jk 7
l
m 1.
-
Avec :
+
12 102 0
/6 12 3po 102 T7
+.
Donc on obtient : (2.7)
/ / rv [λα rv \ / rx [λα rx \
rT rGuTH rGwTH r rT r rT
rt rv rx
Ou : (2.8)
Pour mettre le système d'équations précédent sous une forme adimensionnelle, nous
avons choisi comme longueur de référence la hauteur H de l’enceinte, car c’est sur la hauteur
qu’agissent les forces motrices donnant naissance à l’écoulement convectif. Ainsi, nous avons
introduit les variables adimensionnelles suivantes :
21
Chapitre 2 Modélisation physique et mathématique
a. X = x / H
b. Y = y / H
d. U = u / (ν/H)
e. V = v / (ν/H)
g. Θ G
a H⁄∆ ; (∆
) ; pour la simplicité on prend Θ ≡ T.
h. Φ G
a H⁄∆ ; (∆
) ; pour la simplicité on prend Φ ≡ C
/ +; 0
+ +?
+
(2.11)
/ > / >
+
"
"- √
J / 6 J .
/ > V+ [+ \ / + [ + \W
= + = +GH = +G?H + = + + + +
> +. + +
(2.12)
/ > / >
+
"
"- √
J / 6 J . 6 / > V+ [+\ / + [+\W /
= +? = +G?H = +G??H + ? = + +? + +?
> +. + +
22
Chapitre 2 Modélisation physique et mathématique
• Gt
T ∆H
C ∆H
• G
• P
)
• S
D
K
• D
H
• A = L/H
wu u u uv
• φ
wu tt u u uv
• à t = 0 ; u = v = 0 et T = C = 0
0; 0;U=V=0
rT rC
rv rv
• Y = 0, X:
0; 0, U = V = 0
rT rC
rv rv
• Y= 1, X:
• X = 0, Y : T = 1 ; C = 0 ; U = V = 0.
• X = 4, Y : T = 0 ; C = 1 ; U = V = 0.
• X = 0, Y : T = 1 ; C = 1 ; U = V = 0.
• X = 4, Y : T = 0 ; C = 0 ; U = V = 0.
23
Chapitre 2 Modélisation physique et mathématique
¢&'|£.¤ ¥= [¦
+
\ §
a + £.¤
(2.16)
¢Q ¥a [c
= +
\ §
£.¤ + £.¤
(2.17)
24
Chapitre 3 Résolution numérique
3. RÉSOLUTION NUMÉRIQUE
3.1 INTRODUCTION
Depuis les 20 dernières années, nous vivons aussi une révolution qui a bouleversé
notre vie de tous les jours. Bien entendu, il s’agit d’une révolution informatique. Les
ordinateurs ont envahi les bureaux des chercheurs et aussi les foyers. Les ordinateurs sont,
chaque jour, de plus en plus rapides et, aujourd’hui, à la veille du troisième millénaire, nous
commençons à pouvoir résoudre des problèmes de transfert et de transport en 2-D et 3-D.
Ainsi, des codes commerciaux ont fait leur apparition sur le marché et ils sont
actuellement employés par les grands groupes industriels pour simuler le comportement réel
de leurs installations industrielles.
- minimisation de l’énergie ou principe des travaux virtuels pour les éléments finis,
- loi de conservation et calcul des flux pour la méthode des volumes finis.
La méthode des différences finies consiste à remplacer 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.
25
Chapitre 3 Résolution numérique
La méthode des volumes finis intègre, sur des volumes élémentaires de forme
simple, les équations écrites sous forme de loi de conservation. Elle fournit ainsi de manière
naturelle des approximations discrètes conservatives et elle est donc particulièrement bien
adaptées aux équations de la mécanique des fluides : équation de conservation de la masse,
équation de conservation de la quantité de mouvement et équation de conservation de
l’énergie. Sa mise en œuvre est simple si les volumes élémentaires sont des rectangles (ou des
parallélépipèdes rectangles en trois dimensions). Cependant, la méthode des volumes finis
permet d’utiliser des volumes élémentaires de forme quelconque, donc de traiter des
géométries complexes, ce qui est un avantage sur les différences finies.
Il existe une grande variété de méthodes selon le choix de la géométrie des volumes
élémentaires et des formules de calcul des flux. Par contre, on dispose de peu de résultats
théoriques de convergence.
Dans ce travail, nous avons choisi, d’employer la méthode des volumes finis intègre
pour ses avantages.
26
Chapitre 3 Résolution numérique
3.3.1 Introduction
Une grande majorité des codes commerciaux, actuellement disponible, résolvent les
équations de conservation en variable primitives en emploient la méthode des volumes finis.
Dans cette méthode, le domaine considéré est divisé en volumes finis, un seul volume entoure
un nœud du maillage. Les équations de conservation sont intégrées sur chaque volume fini. La
variation d’une grandeur φ entre le centre de deux volumes finis adjacents peut être linéaire,
quadratique, ou constante dans chaque volume fini. Cette méthode a quelques points de
similitude avec la méthode des différances finies. Tels que, les discrétisations, les méthodes
de résolution des systèmes linéaires tridiagonaux, les notions de stabilité et de consistance
ainsi que le traitement particulier des conditions aux limites
La méthode des volumes finis consiste à intégrer les équations de conservation des
différentes grandeurs physiques ayant la forme générale :
¨
1111111111111127 / «
127 §©^3Γ« S¦¬§Φ
G4ΦH / §©^34
¨
à l'intérieur de petits éléments de volumes, en veillant à ce que le flux d'une grandeur à travers
un élément de surface soit le même pour les deux volumes adjacents. Les équations sont
mises sous forme linéaire de façon à être résolues par une méthode matricielle.
Si on appelle P le point considéré, les points adjacents sont les points est E, ouest W,
sud S et nord N. Les faces du volume sont dénommées de la même façon e, w, s et n (fig. 3.1).
Une des particularités de cette méthode est d'utiliser des maillages décalés (staggered
grid) suivant les grandeurs considérées (fig. 3.2).
27
Chapitre 3 Résolution numérique
∆x
W E ∆y
w e
s
y
S
x
W P E W Ue
w U e Ue w Uw P e E
vs vs
s s
S
S
a- Maillage typique b- Maillage décalé dans la direction X
28
Chapitre 3 Résolution numérique
Y
dxw dxe
N ∆xN
dyn n
W w P e E
∆xP
dys s
S
∆xS
29
Chapitre 3 Résolution numérique
/ / [Γ + \ / + [Γ + \ / «
+« +G«H +GV«H + +« + +«
+. + + +
(3.1)
Equation Ф Г SФ
Continuité 1 0 0
1 ¨
5
Quantité de
J / 6 J .
* ¨ | |
mouvement suivant
U
X
Quantité de ¨ 6 5
1
J / 6 J . 6
mouvement suivant V ¨ | |
*
Y / . T / j
1
Energie T 0
¦
1
Concentration C 0
c
3.5 DISCRÉTISATION
30
Chapitre 3 Résolution numérique
¥. ¥h ¥ §§§ G°=
°a HΔXΔY
.Z∆.
¤̄ +®
+.
(3.2.1)
où °a et °= représentent les valeurs de la variable ° au point P aux temps t et t + ∆t,
respectivement.
¥. ¥h ¥ [Γ \ §§§ / ¥. ¥h ¥ [Γ \ §§§
.Z∆.
¤̄ + +® .Z∆.
¤̄ + +®
+ + + +
(3.2.3)
V[Γ \
[Γ + \ W ΔYΔ / V[Γ + \
[Γ + \ W ΔXΔ
+® +® +® +®
+ ¤ ¯
h
• Terme source :
¥. ¥h ¥ ® §§§ ® ΔXΔYΔ
.Z∆.
¤̄
(3.2.4)
ΔXΔY / _G
°H¤
G
°H¯ bΔY / _G6°H
G6°Hh bΔX V[Γ + \
[Γ + \ W ΔY /
¶
3®Ó µ®´ 7 +® +®
·. ¤ ¯
V[Γ + \
[Γ + \ W ΔX / ® ΔXΔY
+® +® (3.3)
h
Φ¤ ¹ ´ , Φ¯ ´ »
« Z« « Z«
¸ J J ¢ (3.3.1)
Φ
, Φh
«¼ Z«´ «´ Z«½
J J
¢
+«
¾
«¹ µ«´
, ¢+«¾ «´µ«»
¸+« ¤
+ G]H¿ + ¯ G]H À ¢
¢ ¾ ¢
(3.3.2)
, ¾
«¼ µ«´ +« «´ µ«½
+
G]HÁ + h G]HÂ
En remplacent les termes de l’équation (3.3) par leurs expressions (3.3.1) et (3.3.1) et
après réarrangement on obtient l’équation algébrique suivante :
31
Chapitre 3 Résolution numérique
AP ФP = AE ФE + AW ФW + AN ФN +AS ФS + b (3.4)
AP = AE + AW + AN + AS + AaP (3.4.1)
a
··
·.
(3.4.2)
°a est la valeur de la variable ° calculé au temps t et Å est le terme source.
Ǥ
¤ Δ, ǯ
¯ Δ ¢
Æ
Ç
Δ, Çh
h Δ (3.4.4)
|¤ GÉH , |¯ GÉH
¿ È · À È ·
À¢
¸ ¿
|
GÉH , |h GÉH
ÈÁ · ÈÂ ·
(3.4.5)
Á Â
¤ ¿ , ¯
Ê ÊÀ
"À ¢
¸
"¿ (3.4.6)
Á , h
Ê ÊÂ
"Á "Â
|¤ G|¤ |H / ϬÐG
Ǥ , 0H
Í Î
Ñ |¯ G|¯ |H / ϬÐGǯ , 0H¢ (3.4.7)
ÌÒ |
G|
|H / ϬÐG
Ç
, 0H
Ë i |h G|h |H / ϬÐGÇh , 0H
Tous les coefficients de l’équation algébrique (3.4) doivent avoir le même signe
(positif par convenance). Selon le schéma de discrétisation choisi, la fonction G||H peut
prendre différentes expressions. Cinq schémas sont été cités par Patankar (1980) et sont
résumés dans le tableau 3.2.
32
Chapitre 3 Résolution numérique
Î |¤ G1
0.5|¤ |H / ϬÐG
Ǥ , 0H
Ñ |¯ G1
0.5|¯ |H / ϬÐGǯ , 0H
G||H 1
0.5||
Différences
centrées Ò |
G1
0.5|
|H / ϬÐG
Ç
, 0H
i |h G1
0.5|h |H / ϬÐGÇh , 0H
Î |¤ / ϬÐG
Ǥ , 0H
Ñ |¯ / ϬÐGǯ , 0H
G||H 1
Ò |
/ ϬÐG
Ç
, 0H
Upwind
i |h / ϬÐGÇh , 0H
Ò Ï¬ÐG0, 1
0.5|
|H. |
/ ϬÐG
Ç
, 0H
Hybride • P < -2 et P > 2 : upwind
• -2 < P < +2 : différences centrées
i ϬÐG0, 1
0.5|h |H. |h / ϬÐGÇh , 0H
G||H maxG0.0, 1.0
0.5||H
|Ǥ |
Î / ϬÐG
Ǥ , 0H
×ÐI |¿|
1
|ǯ |
Ñ / ϬÐGǯ , 0H
|| ×ÐI |À|
1
G||H
×ÐI||
1 |Ç
|
Exponentiel
Ò / ϬÐG
Ç
, 0H
×ÐI |Á|
1
|Çh |
i / ϬÐGÇh , 0H
×ÐI |Â|
1
ΠϬÐ_0, G1
0.5|¤ |HR b. |¤ / ϬÐG
Ǥ , 0H
Ñ Ï¬Ð_0, G1
0.5|¯ |HR b. |¯ / ϬÐGǯ , 0H
G||H max G0.0, 1.0
0.5||R H
Loi de
puissance Ò Ï¬Ð_0, G1
0.5|
|HR b. |
/ ϬÐG
Ç
, 0H
i ϬÐ_0, G1
0.5|h |HR b. |h / ϬÐGÇh , 0H
33
Chapitre 3 Résolution numérique
/ 0
+9 +:
(3.5)
+ +;
(3.6)
Î9
a
/
9
a
9
a
/
Ñ9
a
6Îa / 6a 6Îi
a
/ 6ia
Ú Ü ∆
Ú Ü ∆ / Ú Ü §¤
Ú Ü §¤ 0
2 2 2 2
a) Suivant X :
= + = +() = +(?) + = + + + +
> +.
/ > +
/ > +
+
"
"- √
J / 6 J .
/ > V+ [+ \ / + [ + \W (3.7)
Après le remplacement des termes avec leurs valeurs d’intégrales dans l’équation
différentielle, on obtient une équation algébrique de la forme :
UNu
Nu
nu
Wu wu Pu eu Eu
UWu UPu UEu
su
USu
Su
34
Chapitre 3 Résolution numérique
VNv
Nv
nv
VWv VPv VEv
Pv
wv sv ev Ev
Wv
VSv
Sv
X
Figure 3.5 : Maillage décalé suivant y : le volume fini est décalé
vers le haut du volume fini typique.
Après le remplacement des termes avec leurs valeurs d’intégrales dans l’équation
différentielle, on obtient une équation algébrique de la forme :
¶
= = ¹Ù Zk¶ = ∆´ ∆. ¶
= = ´Ù ¶
Z»Ù = ∆´ ∆.
9
Î9 <
> J Ý J
Þ ∆ ∆ / > ∆¹
B /
Ñ9 V> J [ J
\ ∆ ∆ / > ∆´
W/
¶
= = ?¹ Z?k¶ = ]¿ ∆. ¶
= = ?¹½ Z?½¶ = ]¿ ∆.
Ò9 <
> J Ý J
Þ §¤ ∆ / > ]Á
B /
i9 V> J [ J
\ §¤ ∆ / > ]Â
W/
.
9
a
§¤ ∆ ∆
= - a J a J
a §¤ ∆ ∆ / (
Î )∆ ∆
ß
X9 / 69
(3.8.1)
> 9 √"
Î
> J [ \ ∆ / > ∆´
/ |¤ ; Ǥ > [ \ ∆ ; |¤ > ∆´ (3.8.2)
¶ ¶ ¶ ¶
= = ¹Ù Z´Ù = ∆ Ê¿ = ¹Ù Z´Ù = ∆
J ¹ J J ¹
Ò
> J [ \ §¤ /
/ |
; Ç
[ \ §¤ ; |
= = ?¹¶ Z?´¶ = ]¿ ÊÁ = ?¹¶ Z?´¶ = ]¿
(3.8.4)
J > ]Á J > J > ]Á
35
Chapitre 3 Résolution numérique
i [ \ §¤ /
/ |h ; Çh [ \ §¤ ; |h
¶
= = ?¹½ Z?½¶ = ]¿ Ê = ¶
?¹½ Z?½¶ = ]¿
> J J > ]Â J > J > ]Â
(3.8.5)
Î / Ñ / Ò / i / / §¤ ∆
= ]¿ ∆´ =
> ∆. "
(3.8.6)
9
9 / G
Î H∆X
ß
X9 / 69 .
9 §¤ ∆
a = ]¿ ∆´ - a J a J a
> ∆. √"
(3.8.7)
b) Suivant Y :
= +? = +(?) = +(??) + ? √ Z? K + +? + +?
/ /
5 6/ V [ \/ [ \W / (3.9)
> +. > + > + + " √ > + + + +
S_` (
= ) / ` (
a )b
Avec :
¶
Ê¿ = ¼ Z´¶ = ]
Î
J
/ |¤ ; Ǥ > [ J
\ §
; |¤ > ]Á (3.10.1)
Á
¶ ¶
ÊÀ = »¼ Z» = ]
Ñ / |¯ ; ǯ > [ \ §
; |¯ > ]Á (3.10.2)
J J À
¶ ¶
ÊÁ = ?¼à Z?´à = ∆´
Ò
J
/ |
; Ç
>
[ J
\ ∆ ; |
> ∆¼
(3.10.3)
¶ ¶
Ê = ?´à Z?½à = ∆
i / |h ; Çh > [ \ ∆ ; |h > ∆´ (3.10.4)
J J ´
= ∆´ ]Á =
Î / Ñ / Ò / i / > / ∆ §
(3.10.5)
∆. "
Ò /
/j Ý Þ ∆ §
2
(3.10.6)
/ / V [λ \/ [λ \W
rT rGUTH rGVTH = r rT r rT
rt rX rY P rX rX rY rY
(3.11)
Î Î / Ñ Ñ / Ò Ò / i i / (3.12)
36
Chapitre 3 Résolution numérique
/ |¤ ; Ǥ
¤ ∆ ; |¤ ]´
Ê¿ á ∆
Î
(3.12.1)
J ¿
Ñ / |¯ ; ǯ
¯ ∆ ; |¯
ÊÀ á ∆´
J ]À
(3.12.2)
Ò
/ |
; Ç
6
∆ ; |
ÊÁ á ∆´
J ]Á
(3.12.3)
i / |h ; Çh 6h ∆ ; |h
ÊÂ á ∆´
J ]Â
(3.12.4)
Î / Ñ / Ò / i /
∆´ ∆´
∆.
(3.12.5)
a
∆´ ∆´
∆.
(3.12.6)
+ = +() = +(?) = + + + +
/ / V [ \/ [ \W (3.13)
+. > + > + i + + +; +;
Î Î / Ñ Ñ / Ò Ò / i i / (3.14)
/ |¤ ; Ǥ
¤a ∆ ; |¤
Ê¿ = = ∆´
Î
J > ij ]¿
(3.14.1)
Ò
/ |
; Ç
6
a ∆ ; |
ÊÁ = = ∆´
J > ij ]Á
(3.14.3)
Î / Ñ / Ò / i /
∆´ ∆´
∆.
(3.14.5)
a
∆´ ∆´
∆.
(3.14.6)
¥ ¥ [+ / N \ §Ð§â 0
+ +
(3.15)
G
¤
¯ H∆ / G6
6h H∆ 0 (3.16)
37
Chapitre 3 Résolution numérique
Avec
9 ã9 / (
Î )∆X (3.17.1)
.
9
a
§¤ ∆
= ]¿ ∆´ - J J
ã9
9
a
ß
X9
a
/ 69
a
(3.17.2)
> ∆. √"
ä¤ / §¤ G
Î H ;
¤
§¤
∆´
å´
ä¯ / §¯ GÑ
H ;
¯
§¯
∆´
å´
(3.19)
6
6æ / §
G
Ò H ; §
∆´
å´
6h 6æh / §h Gi
H ; §h
∆´
å´
6æ
¹ ¹à » »à ¼ ¼à ½ ½à à
å ? Zå ? Zå ? Zå ? Zç
å´
(3.20.2)
On remplace les vitesses par leur expression on obtient la forme algébrique suivante :
Î Î / Ñ Ñ / Ò Ò / i i / (3.21)
Avec :
Î / / Ò / i (3.21.2)
ä¤
3
ä¯ 7∆ / 36æ
6æh 7∆ 0 (3.21.3)
38
Chapitre 3 Résolution numérique
Pour résoudre les équations des vitesses, il nous faut la pression (inconnue) et pour
résoudre l’équation de la pression nous avons besoin des vitesses (inconnues).
On définit :
è / f ; 6 6è / 6f ; è / f (3.22)
9
è
Î
Î9
è
/ Ñ
Ñ9
è
/ Ò
Ò9
è
/ i
i9
è
/ ã9 / Gè
Îè H∆ (3.23.1)
6:
è
Î 6Î:
è
/ Ñ 6Ñ:
è
/ Ò 6Ò:
è
/ i 6i:
è
/ ã: / Gè
Òè H∆ (3.23.2)
Ces équations montrent qu’avec une estimation de la pression, on peut obtenir une
estimation des vitesses.
9
′
Î
Î9
′
/ Ñ
Ñ9
′
/ Ò
Ò9
′
/ i
i9
′
/ ã9 / 3′
Î′ 7∆ (3.24.1)
6:
′
Î 6Î:
′
/ Ñ 6Ñ:
′
/ Ò 6Ò:
′
/ i 6i:
′
/ ã: / 3′
Ò′ 7∆ (3.24.2)
On constate que les corrections des vitesses à un point sont fonctions des corrections
de la pression et des vitesses aux points voisins. Donc si on a la correction de la pression, on
peut résoudre les systèmes d’équations des corrections des vitesses pour obtenir ces dernières.
On peut donc redéfinir les corrections des vitesses seulement en fonction de la correction de la
pression comme suit.
9
′
3′
Î′ 7∆ (3.25.1)
6:
′
3′
Ò′ 7∆ (3.25.2)
39
Chapitre 3 Résolution numérique
¯
è̄ /
¯′
è̄ / §¯ 3Ñ
′
′ 7 (3.26.2)
6
6
è / 6
′ 6
è / §
3′
Ò′ 7 (3.26.3)
ê
¤è / §¤ 3′
Î′ 7ë∆
ê
è̄ / §¯ 3Ñ
′
′ 7ë∆ / ê6
è / §
3′
Ò′ 7ë∆
ê6hè / §h 3i′
′ 7ë∆ 0
(3.27)
′ Î Î′ / Ñ Ñ
′
/ Ò Ò′ / i i′ / ′ (3.28)
Avec :
Î / Ñ / Ò / i (3.28.2)
′ G
¤è
è̄ H∆ / G6
è
6hè H∆ 0 (3.28.3)
_b_Φb _b
Où :
40
Chapitre 3 Résolution numérique
Où A est une matrice tridiagonale, alors l'algorithme de Gauss peut être simplifié, et cela en
tenant compte de la structure particulière de A.
b1 c1 x1 y1
a b2 c2 0 x y
2 2 2
a3 b3 c3 x3 y 3
... ... ... ... = ...
ai bi ci xi y i
... ... ... ... ...
0 an −1 bn −1 c n −1 x n −1 y n −1
an bn xn y n
Alors l'algorithme TDMA peut être résumé dans les étapes suivantes :
• Triangularisation de la matrice A :
c1
γ1 =
b1
ci
γi = i = 2, n − 1
bi − ai γ i −1
y1
β1 =
b1
y i − ai β i − 1
βi = i = 2, n − 1
bi − ai γ i −1
xn = β n
xi = βi − γ i xi −1 i = n − 1, n − 2, n − 3,......,1
41
Chapitre 3 Résolution numérique
Dans le cas où le régime est permanant la convergence est vérifiée après chaque
itération. La différence entre deux valeurs successive ente deux itérations consécutives aux
temps t et ∆t, des variables T, C, U et V, doit négligeable telle que :
«ïð∆ï ï
í,î µ«í,î
Ýì «ïð∆ï
ìÞ ñ 10µN (3.29)
í,î
42
Chapitre 3 Résolution numérique
En plus, les bilans thermique et massique doivent être satisfaits en comparant les
nombres moyens de Nusselt et de Sherwood près des parois verticales entre deux incréments
consécutives de temps de telle sorte qu’il y a conservation de chaleur et de masse :
|&'£=
&'£N | ñ 10µL
Q£=
£N Q ñ 10µL
(3.30)
Un code de calcul basé sur l’algorithme précédent en langage Fortran a été réalisé et
exécuté sur un micro-ordinateur personnel Pentium 4 de 3.4 GHz et 1.256 Go de RAM. Les
calculs numériques ont été effectués pour un maillage de 164×88 avec un pas de temps 10-4.
43
Chapitre 4 Résultats et discussions
4. RÉSULTATS ET DISCUSSIONS
Le tableau 4.1 montre la sensibilité des résultats numériques aux maillages. Nous
constatons que les maillages 1 jusqu’à 4 donnent des résultats ayant des écarts relativement
important par rapport aux maillages 5, 6, 7 et 8. Par conséquent, ils ne sont pas pris en
considération. Nous remarquons que les différences relatives des résultats obtenus avec le
maillage 6, par rapport au maillage 7, sont : 0.03%, 0.33%, 0.25% et 0.18% respectivement
pour Nu , Sh , |Umax| et |Vmax| sont largement admissibles. On remarque aussi d'après le
tableau 4.1 que le transfert de masse est plus sensible au maillage que le transfert de chaleur.
Pour étudier effet de l’espacement de maillage, le maillage 164x88 a été raffiné près
des parois verticales avec une progression géométrique de valeur C. Les résultats obtenus sont
avec les mêmes paramètres sont illustrés sur le tableau 4.2.
44
Chapitre 4 Résultats et discussions
Nous remarquons que les différences relatives des résultats obtenus avec C=1.03 par
rapport à C=1.01, sont : 0.2%, 0.95%, 0.3% et 0.5% respectivement pour Nu , Sh , |Umax| et
|Vmax| sont largement admissibles.
1.00
Maillage 122x66
Maillage 132x82
Maillage 164x88
0.60
0.40
0.20
0.00
45
Chapitre 4 Résultats et discussions
Ces résultats montrent que le maillage 164x88 est suffisant pour assurer une bonne
résolution spatiale du problème étudié. Nous avons choisi d’utiliser dans la suite de nous
calculs un maillage raffiné avec C=1.03, car il permet de mieux décrire l’effet de des forces
de buoyancy près des parois verticales actives. Donc tous les cas présentés dans cette étude
sont obtenus avec ce maillage raffiné de 164x88.
Dans la plupart des cas, nos constatons un bon accord de nos résultats avec ceux de
Younsi et al. (2002). En effet, nos résultats présentent une différence relative inferieur à 4.2 %
par rapport aux autres résultats.
Nusselt
Lauriat et Trevesion et Nithiarasu et Younsi,
Ra*=Ra.Da Présent
Prassad Bejan al Harkati et
travail
Kalache
10 1.07 / 1.08 1.06 1.078
50 / 2.02 1.958 1.936 1.976
100 3.09 3.27 3.02 2.98 3.089
500 / / 8.38 8.32 8.69
1000 13.41 18.38 12.514 12.32 12.7
46
Chapitre 4 Résultats et discussions
Nusselt
Younsi, Présent
Ra* Da Lauriat et Nithiarasu et
Harkati et
Prassad al travail
Kalache
10 10-6 1.07 1.08 1.06 1.078
100 10-6 3.06 3.00 2.98 3.089
10 10-2 1.02 1.02 0.99 1.003
100 10-2 1.7 1.71 1.68 1.24
(A1)
(A2)
(B1)
47
Chapitre 4 Résultats et discussions
(B2)
(C1)
(C2)
Figure 4.2 : Lignes de courant, isothermes et isoconcentrations
1- Younsi et al. (2002), 2- Présent travail
a- N=3 ; b- N=6 ; c- N=10.
48
Chapitre 4 Résultats et discussions
Dans ce cas, avec les conditions aux limites choisies, la température maximale se
situe près de la paroi verticale gauche de l’enceinte étudiée. Par contre, la concentration
maximale se trouve près de la paroi verticale droite.
Pour différentes valeurs du nombre Grashof thermique, à savoir Grt=104, 105, 106,
107 et 108, les résultats obtenus sont illustrés sur les figures 4.3, 4.4, 4.5 et 4.6,
respectivement.
49
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
50
Chapitre 4 Résultats et discussions
d- Isoconcentrations
51
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
52
Chapitre 4 Résultats et discussions
d- Isoconcentrations
53
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
54
Chapitre 4 Résultats et discussions
d- Isoconcentrations
Les figures 4.6 et 4.7 représentent la répartition de la pression, les contours de lignes
de courant, de la température et de la concentration obtenues, respectivement, pour Grt=108 et
Grc=105. Dans ce cas, on constate que les plus grandes valeurs de la pression se situent en
haut de l’enceinte près de la paroi horizontale supérieure et ceci est dû aux conditions aux
limites thermique et solutale et aux valeurs prises par les nombres de Grashof thermique et
solutal (fig. 4.6a et 4.7a). Cependant, l’écoulement du fluide obtenue a une structure
multicellulaire ; une cellule principale près des parois de l’enceinte et des cellules secondaires
au milieu de l’enceinte. Ces différentes cellules ont pour sens de rotation le sens horaire (fig.
4.6b et 4.7b). On obtient une distribution de la température faiblement variable suivant X mais
variable suivant Y de façon à satisfaire la stratification thermique stable ; la température croît
suivant Y (fig. 4.6c et 4.7c). Alors que, les gradients de concentration sont très fortement
concentrés près des parois verticales isothermes. Tandis qu’au centre de l’enceinte la
concentration a une valeur constante, C≈0.5 (fig. 4.6d et 4.7d). On conclut donc que dans ces
conditions, l’écoulement du fluide est totalement conduit par la force de buoyancy thermique.
55
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
56
Chapitre 4 Résultats et discussions
d- Isoconcentrations
a- Contour de pression
b- Lignes de courant
57
Chapitre 4 Résultats et discussions
c- Isothermes
d- Isoconcentrations
Le tableau 4.5 représente les taux moyens de transfert de chaleur et de masse et les
vitesses maximales longitudinale et transversale dans l’enceinte. Les forces de volume sont
coopérantes. L’accroissement du nombre Grashof thermique de Grt=104 à Grt=108 a entrainé
une augmentation considérable et des nombres moyens de Nusselt et de Sherwood et une
accélération très importante de l’écoulement. C’est ainsi que les valeurs de Nu , Sh , |Umax|
et |Vmax| sont passées de 0.594, 3.073, 67.396 et 83.958 pour Grt=104 à 13.662, 27.708,
3493.656 et 4522.051 pour Grt=108, respectivement.
58
Chapitre 4 Résultats et discussions
Nu Sh |Umax| |Vmax|
Dans ce cas, avec les conditions aux limites choisies, les températures maximales se
situent près de la paroi verticale gauche de l’enceinte étudiée. Par contre, les concentrations
maximales se trouvent près de la paroi verticale droite.
4.3.2.1.1 Cas de Grc=105
Pour différentes valeurs du nombre Grashof thermique, à savoir Grt=104, 105, 106,
107 et 108, les résultats obtenus sont illustrés sur les figures suivantes.
La figure 4.8 illustre typiquement la répartition de la pression, les contours de lignes
de courant, de la température et de la concentration obtenues, respectivement, pour Grt=104 et
Grc=105. On constate que les plus grandes valeurs de la pression se situent en bas du coté de
la paroi verticale droite (où la concentration est maximale), et ceci est dû aux conditions aux
limites thermique et solutale et aux valeurs prises par les nombres de Grashof thermique et
solutal (fig. 4.8a). Un courant rotatif monocellulaire dans le sens horaire est obtenu dans
l’enceinte (fig. 4.8b). La variation de la température dans l’enceinte est présentée par des
isothermes verticales légèrement déformés de gauche à droite en haut de l’enceinte et de
droite à gauche en bas de l’enceinte (fig. 4.8c). Tandis que, la variation de la concentration
dans l’enceinte est présentée par des isoconcentrations verticales largement déformés de
59
Chapitre 4 Résultats et discussions
droite à gauche en bas de l’enceinte et de gauche à droite en haut de l’enceinte (fig. 4.8d). On
conclut donc que, dans ces conditions, la force de buoyancy solutale conduit,
particulièrement, la circulation du fluide dans le milieu poreux.
a- Contour de pression
b- Lignes de courant
c- Isothermes
60
Chapitre 4 Résultats et discussions
d- Isoconcentrations
61
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
62
Chapitre 4 Résultats et discussions
d- Isoconcentrations
63
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
64
Chapitre 4 Résultats et discussions
d- Isoconcentrations
a- Contour de pression
b- Lignes de courant
65
Chapitre 4 Résultats et discussions
c- Isothermes
d- Isoconcentrations
66
Chapitre 4 Résultats et discussions
la concentration a une valeur constante, C=0.5 (fig. 4.12d). On conclut donc que dans ces
conditions, l’écoulement du fluide dans le milieu poreux est totalement conduit par la force de
buoyancy thermique.
a- Contour de pression
b- Lignes de courant
c- Isothermes
67
Chapitre 4 Résultats et discussions
d- Isoconcentrations
La variation temporelle du nombre de Nusselt est illustrée dans la figure 4.13, qui
montre qu’avec le temps le nombre de Nusselt diminue et se stabilise en fin de compte pour
les cinq cas étudiés. Les valeurs de Nusselt sont plus élevées pour les valeurs de Grashof
thermique élevée.
La figure 4.14 illustre la variation du nombre de Sherwood avec le temps. On
constate qu’après une légère fluctuation, le nombre de Sherwood diminue avec le temps et se
stabilise en fin de compte pour les cinq cas. Les valeurs de Sherwood sont plus élevées pour
les valeurs de Grashof thermique élevée.
68
Chapitre 4 Résultats et discussions
4
Grt=10
14 5
Grt=10
Nu
13
6
12 Grt=10
11 7
Grt=10
10
8
9 Grt=10
8
7
6
5
4
3
2
1
0
0 10000 20000 30000 40000 50000
t
4
Grt=10
5
Grt=10
6
Grt=10
24
7
Grt=10
22
Sh
8
Grt=10
20
18
16
14
12
10
8
6
4
2
0
0 10000 20000 30000 40000 50000
t
69
Chapitre 4 Résultats et discussions
Le tableau 4.6 représente les transferts moyens de chaleur et de masse et les vitesses
maximales longitudinale et transversale dans l’enceinte. Les forces de buoyancy sont
coopérantes, l’augmentation du nombre Grashof thermique a entrainé l’accroissement et des
nombres de Nusselt et de Sherwood et l’accélération de l’écoulement. C’est ainsi que les
valeurs de Nu , Sh , |Umax| et |Vmax| sont passées de 0.280, 1.052, 16.978 et 12.964 pour
Grt=104 à 6.843, 18.810, 508.175 et 1341.123 pour Grt=108, respectivement.
Nu Sh |Umax| |Vmax|
Pour différentes valeurs du nombre Grashof thermique, à savoir Grt=104, 105, 106,
107 et 108. Les résultats obtenus sont illustrés sur les figures suivantes.
70
Chapitre 4 Résultats et discussions
dans l’enceinte est présentée par des isoconcentrations verticales largement déformés de
droite à gauche en bas de l’enceinte et de gauche à droite en haut de l’enceinte (fig. 4.15d).
On conclut donc que, dans ces conditions, la force de buoyancy solutale conduit,
particulièrement, la circulation du fluide dans le milieu poreux.
a- Contour de pression
b- Isothermes
71
Chapitre 4 Résultats et discussions
c- Isoconcentrations
d- Lignes de courant
72
Chapitre 4 Résultats et discussions
gauche à droite en haut de l’enceinte (fig. 4.16d). On conclut donc que, dans ces conditions, la
force de buoyancy solutale conduit, particulièrement, l’écoulement dans le milieu poreux.
a- Contour de pression
b- Isothermes
c- Isoconcentrations
73
Chapitre 4 Résultats et discussions
d- Lignes de courant
Figure 4.16 : Contour de la pression, isothermes, isoconcentrations et structure de
l’écoulement pour : Grt=105, Grc=5.105.
74
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Isothermes
c- Isoconcentrations
75
Chapitre 4 Résultats et discussions
d- Lignes de courant
a- Contour de pression
b- Lignes de courant
76
Chapitre 4 Résultats et discussions
c- Isothermes
d- Isoconcentrations
77
Chapitre 4 Résultats et discussions
(fig. 4.19d). On conclut donc que dans ces conditions, l’écoulement du fluide dans le milieu
poreux est totalement conduit par la force de buoyancy thermique.
a- Contour de pression
b- Lignes de courant
c- Isothermes
78
Chapitre 4 Résultats et discussions
d- Isoconcentrations
La variation de nombre de Nusselt avec le temps est illustrée dans les Figure 4.20.
Cette figure montre qu’avec le temps le nombre de Nusselt diminue et se stabilise en fin de
compte pour les cinq nombres de Grashof thermique. Les valeurs du Nusselt augmentent avec
les valeurs du Grashof thermique.
4
Grt=10
14 5
Grt=10
13
Nu
6
12 Grt=10
11 7
Grt=10
10
8
9 Grt=10
8
7
6
5
4
3
2
1
0
0 10000 20000 30000 40000 50000
t
79
Chapitre 4 Résultats et discussions
4
Grt=10
5
Grt=10
25 6
Grt=10
7
Grt=10
Sh
8
Grt=10
20
15
10
Le tableau 4.7 représente les transferts moyens de chaleur et de masse et les vitesses
maximales longitudinale et transversale dans l’enceinte. Ce cas présente une similarité avec le
cas précédent de Grc=105. Les forces de buoyancy sont coopérantes, l’augmentation du
nombre Grashof thermique a entrainé l’accroissement et des nombres de Nusselt et de
Sherwood et l’accélération de l’écoulement. C’est ainsi que les valeurs de Nu , Sh , |Umax| et
|Vmax| sont passées de 0.386, 2.737, 35.430 et 53.368 pour Grt=104 à 6.848, 18.835, 508.224
et 1342.978 pour Grt=108, respectivement.
80
Chapitre 4 Résultats et discussions
Nu Sh |Umax| |Vmax|
Pour différentes valeurs du nombre Grashof thermique, à savoir Grt=104, 105, 106,
107 et 108 Les résultats obtenus sont illustrés sur les figures suivantes.
Les figures 4.22 et 4.23 illustrent la répartition de la pression, les contours de lignes
de courant, de la température et de la concentration obtenues, respectivement, pour Grt=104 et
105, Grc=106. On constate pour ces deux cas, que les plus grandes valeurs de la pression se
situent près de la paroi horizontale basse du côté droit de l’enceinte (fig. 4.22a et 4.23a).
L’écoulement obtenu dans l’enceinte est monocellulaire dans le sens horaire (fig. 4.22b et
4.23b). La variation de la température dans l’enceinte est présentée par des isothermes
verticales légèrement déformés de gauche à droite en haut de l’enceinte et de droite à gauche
en bas de l’enceinte (fig. 4.22c et 4.23c). Tandis que, la variation de la concentration dans
l’enceinte est présentée par des isoconcentrations largement déformés de droite à gauche en
bas de l’enceinte et de gauche à droite en haut de l’enceinte. La distribution de la
concentration obtenu et faiblement variable suivant X mais variable suivant Y de façon à
satisfaire la stratification solutale stable ; la concentration décroît suivant Y (fig. 4.22d et
4.23d). On conclut donc que, dans ces conditions, la force de buoyancy solutale conduit,
particulièrement, la circulation du fluide dans le milieu poreux.
81
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
82
Chapitre 4 Résultats et discussions
d- Isoconcentrations
a- Contour de pression
b- Lignes de courant
83
Chapitre 4 Résultats et discussions
c- Isothermes
d- Isoconcentrations
84
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
85
Chapitre 4 Résultats et discussions
d- Isoconcentrations
86
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
87
Chapitre 4 Résultats et discussions
d- Isoconcentrations
88
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
89
Chapitre 4 Résultats et discussions
d- Isoconcentrations
La variation temporelle du nombre de Nusselt est illustrée dans la figure 4.27, qui
montre qu’avec le temps le nombre de Nusselt diminue et se stabilise pour les cinq nombres
de Grashof thermique. Les valeurs de Nusselt augmentent avec les valeurs du Grashof
thermique.
4
Grt=10
5
14 Grt=10
Nu
6
Grt=10
12 7
Grt=10
10 8
Grt=10
8
0
0 10000 20000 30000 40000 50000
t
90
Chapitre 4 Résultats et discussions
4
Grt=10
5
Grt=10
25 6
Grt=10
7
Grt=10
Sh
8
Grt=10
20
15
10
D’après ces résultats, on peut conclure que les transferts de chaleur et de masse
augmentent avec les valeurs des nombres de Grashof thermique et solutal. Ainsi, les valeurs
des nombres de Nusselt et de Sherwood sont plus élevées dans le deuxième cas (Grc=5.105)
que dans le premier (Grc=105), de même, celles du troisième cas (Grc=106) sont supérieures à
celles du deuxième cas. Le tableau 4.8 représente les transferts moyens de chaleur et de masse
et les vitesses maximales longitudinale et transversale dans l’enceinte. Ce cas est aussi
similaire aux cas précédents (Grc=105 et Grc=5.105) avec une légère augmentation des valeurs
caractéristiques. Les forces de buoyancy sont coopérantes, l’augmentation du nombre Grashof
thermique a entrainé l’accroissement et des nombres de Nusselt et de Sherwood et
l’accélération de l’écoulement. C’est ainsi que les valeurs de Nu , Sh , |Umax| et |Vmax| sont
passées de 0.484, 3.936, 46.708 et 92.622 pour Grt=104 à 6.854, 18.867, 508.279 et 1345.294
pour Grt=108, respectivement.
91
Chapitre 4 Résultats et discussions
Nu Sh |Umax| |Vmax|
Pour différentes valeurs du nombre Grashof thermique, à savoir Grt=104, 105, 106,
107 et 108 Les résultats obtenus sont illustrés sur les figures 4.29, 4.30, 4.31, 4.32 et 4.33.
On constate que pour les valeurs de Grt=104, 105 et 106, les transferts de chaleur et de
masse sont pratiquement diffusifs. Les valeurs de la pression engendrée sont relativement
faibles (fig. 4.29a et 4.30a) à modérées (fig. 4.31a) par rapport au cas précédent de Da=10-2.
La structure de l’écoulement est monocellulaire et a pour sens de rotation le sens horaire (fig.
4.29b, 4.30b et 4.31b). Les transferts de chaleur et matière sont représentés, respectivement,
par des isothermes (fig. 4.29c, 4.30c et 4.31c) et des isoconcentrations (fig. 4.29d, 4.30d et
4.31d) verticales quasi-parallèles. Pour Grt=107 et 108, l’écoulement et les transferts de
chaleur et de masse sont conduit par convection et diffusion. On constate que les plus grandes
valeurs de la pression se situent en haut du coté de la paroi verticale gauche et sont
relativement assez élevées (fig. 4.32a et 4.33a). Un écoulement rotatif monocellulaire dans le
sens horaire est obtenu dans l’enceinte (fig. 4.32b et 4.33b). La variation de la température
dans l’enceinte est présentée par des isothermes verticales quasi-parallèles (fig. 4.32c et
4.33c). Tandis que, la variation de la concentration dans l’enceinte est présentée par des
isoconcentrations verticales légèrement déformés de droite à gauche en bas de l’enceinte et de
gauche à droite en haut de l’enceinte (fig. 4.32d et 4.33d). On conclut donc que, dans ces
conditions, le mode de transfert par diffusion est dominant.
92
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
93
Chapitre 4 Résultats et discussions
d- Isoconcentrations
a- Contour de pression
b- Lignes de courant
94
Chapitre 4 Résultats et discussions
c- Isothermes
d- Isoconcentrations
a- Contour de pression
95
Chapitre 4 Résultats et discussions
b- Lignes de courant
c- Isothermes
d- Isoconcentrations
96
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
97
Chapitre 4 Résultats et discussions
d- Isoconcentrations
Figure 4.32 : Contour de la pression, isothermes, isoconcentrations et structure de
l’écoulement pour : Grt=107, Grc=105.
a- Contour de pression
b- Lignes de courant
98
Chapitre 4 Résultats et discussions
c- Isothermes
d- Isoconcentrations
Le tableau 4.9 illustre les valeurs des caractéristiques de l’écoulement des transferts
thermique et solutale. On remarque que les vitesses maximales obtenues, |Umax| et |Vmax|,
sont très faibles pour les Grt=104, 105 et 106, faible pour Grt=107 et s’élèvent légèrement pour
Grt=108 à 6.058 et 4.567, respectivement. Tandis que les valeurs des nombres de Nusselt et
Sherwood, Nu et Sh , n’ont pas évoluées d’une manière significative avec Grt.
99
Chapitre 4 Résultats et discussions
Nu Sh |Umax| |Vmax|
Les figures 4.34 et 4.35 représente l’évolution temporelle des nombre de Nusselt et
de Sherwood pour les différentes valeurs de Grt. On remarque qu’avec le temps Nu et Sh
diminuent, rapidement, et se stabilise en fin de compte pour les cinq cas étudiés autour d’une
valeur moyenne. Pratiquement, l’augmentation de Grt n’a pas d’influence sur les transferts de
chaleur et de masse. C’est ainsi que les valeurs de Nu et de Sh sont passées tout simplement
de 0.250 et 0.250, pour Grt=104 à 0.254 et 0.385 pour Grt=108, respectivement.
100
Chapitre 4 Résultats et discussions
4
Grt=10
3,0
5
Grt=10
Nu
6
2,5 Grt=10
7
Grt=10
2,0
8
Grt=10
1,5
1,0
0,5
0,0
0 10000 20000 30000 40000 50000
t
4
4,5 Grt=10
5
4,0 Grt=10
6
3,5
Grt=10
Sh
7
Grt=10
3,0
8
Grt=10
2,5
2,0
1,5
1,0
0,5
0,0
101
Chapitre 4 Résultats et discussions
Dans ce cas, avec les conditions aux limites choisies, les températures et les
concentrations maximales se situent prés de la paroi verticale gauche de l’enceinte étudiée.
Pour différentes valeurs du nombre Grashof thermique, à savoir Grt=104, 105 et 106.
Les résultats obtenus sont illustrés sur les figures suivantes.
a- Contour de pression
102
Chapitre 4 Résultats et discussions
b- Lignes de courant
c- Isothermes
d- Isoconcentrations
103
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
104
Chapitre 4 Résultats et discussions
c- Isothermes
d- Isoconcentrations
105
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
106
Chapitre 4 Résultats et discussions
d- Isoconcentrations
Pour différentes valeurs du nombre Grashof thermique, à savoir Grt=104, 105 et 106.
Les résultats obtenus sont illustrés sur les figures suivantes.
Les figures 4.39 et 4.40 illustrent la répartition de la pression, les contours de lignes
de courant, de la température et de la concentration obtenues, respectivement, pour Grt=104 et
105, Grc=5.105. On constate que les plus grandes valeurs de la pression se situent près de la
paroi horizontale basse du coté de la paroi verticale gauche (fig. 4.39a et 4.40). L’écoulement
obtenu dans l’enceinte est monocellulaire et a pour sens de rotation le sens antihoraire (fig.
4.39b et 4.40b). La variation de la température dans l’enceinte est présentée par des
isothermes verticales légèrement déformés de droite à gauche en haut de l’enceinte et de
gauche à droite en bas de l’enceinte (fig. 4.39c et 4.40c). Tandis que, la variation de la
concentration dans l’enceinte est présentée par des isoconcentrations verticales largement
déformés de gauche à droite en bas de l’enceinte et de droite à gauche en haut de l’enceinte
(fig. 4.39d et 4.40d). On conclut donc que, dans ces conditions, la force de buoyancy solutale
conduit, particulièrement, l’écoulement dans le milieu poreux.
107
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
108
Chapitre 4 Résultats et discussions
d- Isoconcentrations
a- Contour de pression
b- Lignes de courant
109
Chapitre 4 Résultats et discussions
c- Isothermes
d- Isoconcentrations
110
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
111
Chapitre 4 Résultats et discussions
d- Isoconcentrations
Le tableau 4.10 représente les transferts moyens de chaleur et de masse et les vitesses
maximales longitudinale et transversale dans l’enceinte. Les forces de buoyancy sont
opposantes, on peut distinguer trois situations différentes. Dans un premier temps, la force de
buoyancy solutale dominante conduit l’écoulement et les transferts de chaleur et de masse
(Grt=104 avec Grc=105, 5.105 et Grt=105 avec Grc=5.105). Quand les forces de buoyancy
thermique et solutale deviennent comparables (Grc=105 et Grt=105), pratiquement leurs
actions s’annulent, les transferts de chaleur et de masse sont purement diffusifs et
l’écoulement est presque négligeable (Nu =0.249, Sh =0.250, |Umax|=0.017 et
|Vmax|=0.0059). Ensuite les forces de buoyancy thermique deviennent dominantes,
l’écoulement s’accélère et les valeurs des nombres de Nusselt et de Sherwood accroissent
avec le nombre de Grashof thermique (Grt=106, Grc=105 et 5.105).
112
Chapitre 4 Résultats et discussions
Nu Sh |Umax| |Vmax|
Dans ce cas, avec les conditions aux limites choisies, la température maximale se
situe prés de la paroi verticale gauche de l’enceinte étudiée. Par contre, la concentration
maximale se trouve prés de la paroi verticale droite.
4.3.3.1 Cas de Da=10-2
Pour différentes valeurs du nombre Grashof thermique, à savoir Grt=104, 105, 106,
107 et 108, les résultats obtenus sont illustrés sur les figures suivantes.
La figure 4.42 illustre la répartition de la pression, les contours de lignes de courant,
de la température et de la concentration obtenues, respectivement, pour Grt=104, Grc =105. Ce
cas est similaire au cas du modèle de Darcy-Brinkman-Forchheimer (§ 4.3.2.1.1.1 page 60).
Cependant, on remarque que les taux de transfert de chaleur et de masse sont plus élevées
(tableau 4.11).
113
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
114
Chapitre 4 Résultats et discussions
d- Isoconcentrations
a- Contour de pression
115
Chapitre 4 Résultats et discussions
b- Lignes de courant
c- Isothermes
d- Isoconcentrations
116
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
117
Chapitre 4 Résultats et discussions
c- Isothermes
d- Isoconcentrations
118
Chapitre 4 Résultats et discussions
isothermes. Tandis qu’au centre de l’enceinte la concentration a une valeur constante, C≈0.5
(fig. 4.45d et 4.46d). Mais, on remarque aussi que les taux de transfert de chaleur et de masse
sont plus élevées (tableau 4.11). On conclut donc que dans ces conditions, l’écoulement du
fluide dans le milieu poreux est totalement conduit par la force de buoyancy thermique.
a- Contour de pression
b- Lignes de courant
119
Chapitre 4 Résultats et discussions
c- Isothermes
d- Isoconcentrations
a- Contour de pression
120
Chapitre 4 Résultats et discussions
b- Lignes de courant
c- Isothermes
d- Isoconcentrations
121
Chapitre 4 Résultats et discussions
Les figures 4.47 et 4.48 illustrent la répartition de la pression, les contours de lignes
de courant, de la température et de la concentration obtenues, respectivement, pour Grt=104 et
Grt=105, Grc=5.105. On constate que ces deux cas sont similaires aux deux cas du modèle de
Darcy-Brinkman-Forchheimer (§ 4.3.2.1.2.1 pages 71 et 73). Néanmoins, les taux de transfert
de chaleur et de masse sont plus élevées (tableau..).
a- Contour de pression
b- Lignes de courant
122
Chapitre 4 Résultats et discussions
c- Isothermes
d- Isoconcentrations
a- Contour de pression
123
Chapitre 4 Résultats et discussions
b- Lignes de courant
c- Isothermes
d- Isoconcentrations
124
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
125
Chapitre 4 Résultats et discussions
c- Isothermes
d- Isoconcentrations
126
Chapitre 4 Résultats et discussions
qu’au centre de l’enceinte la concentration a une valeur constante, C≈0.5 (fig. 4.50d et 4.51d).
Mais, on remarque aussi que les taux de transfert de chaleur et de masse sont plus élevées
(tableau 4.11). On conclut donc que dans ces conditions, l’écoulement du fluide dans le milieu
poreux est totalement conduit par la force de buoyancy thermique.
a- Contour de pression
b- Lignes de courant
c- Isothermes
127
Chapitre 4 Résultats et discussions
d- Isoconcentrations
a- Contour de pression
b- Lignes de courant
128
Chapitre 4 Résultats et discussions
c- Isothermes
d- Isoconcentrations
Pour les différentes valeurs du nombre Grashof thermique, à savoir Grt=104, 105, 106,
107 et 108 Les résultats obtenus sont illustrés sur les figures 4.52, 4.53, 4.54, 4.55 et 4.56.
On constate que ce cas est pratiquement comparable au cas du modèle de Darcy-
Brinkman-Forchheimer (§ 4.3.2.1.4.1 page 93). Pour les valeurs de Grt=104, 105 et 106, les
transferts de chaleur et de masse sont pratiquement diffusifs. Les valeurs de la pression
engendrée sont relativement faibles (fig. 4.52a et 4.53a) à modérées (fig. 4.54a) par rapport au
cas précédent de Da=10-2. L’écoulement est monocellulaire tournant dans le sens horaire (fig.
4.52b, 4.53b et 4.54b). Les transferts de chaleur et de masse sont représentés, respectivement,
par des isothermes (fig. 4.52c, 4.53c et 4.54c) et des isoconcentrations (fig. 4.52d, 4.53d et
4.54d) verticales quasi-parallèles. Pour Grt=107 et 108, l’écoulement et les transferts de
chaleur et de masse sont conduit par convection et diffusion. On constate que les plus grandes
valeurs de la pression se situent en haut du coté de la paroi verticale gauche et sont
129
Chapitre 4 Résultats et discussions
relativement assez élevées (fig. 4.55a et 4.56a). Un écoulement rotatif monocellulaire dans le
sens horaire est obtenu dans l’enceinte (fig. 4.55b et 4.56b). La variation de la température
dans l’enceinte est présentée par des isothermes verticales quasi-parallèles (fig. 4.55c et
4.56c). Tandis que, la variation de la concentration dans l’enceinte est présentée par des
isoconcentrations verticales légèrement déformés de droite à gauche en bas de l’enceinte et de
gauche à droite en haut de l’enceinte (fig. 4.55d et 4.56d). On conclut donc que, dans ces
conditions, le mode de transfert par diffusion est dominant.
a- Contour de pression
b- Lignes de courant
130
Chapitre 4 Résultats et discussions
c- Isothermes
d- Isoconcentrations
a- Contour de pression
131
Chapitre 4 Résultats et discussions
b- Lignes de courant
c- Isothermes
d- Isoconcentrations
132
Chapitre 4 Résultats et discussions
a- Contour de pression
b- Lignes de courant
c- Isothermes
133
Chapitre 4 Résultats et discussions
d- Isoconcentrations
a- Contour de pression
b- Lignes de courant
134
Chapitre 4 Résultats et discussions
c- Isothermes
d- Isoconcentrations
Figure 4.55 : Contour de la pression, isothermes, isoconcentrations et structure de
l’écoulement pour : Grt=107, Grc=105.
a- Contour de pression
135
Chapitre 4 Résultats et discussions
b- Lignes de courant
c- Isothermes
d- Isoconcentrations
136
Chapitre 4 Résultats et discussions
Modèle de Darcy-Brinkman-Forchheimer
Modèle de Darcy
10
Nu
137
Chapitre 4 Résultats et discussions
Modèle de Darcy-Brinkman-Forchheimer
Modèle de Darcy
28
26
Sh
24
22
20
18
16
14
12
10
8
6
4
2
0
La différence entre les nombres Nusselt et de Sherwood entre les deux modèles est
d’autant plus élevée, que l’on augmente le nombre de Grashof solutal (tableau 4.11).
Les figures 4.61 et 4.62 illustrent les variations des nombres moyens du Nusselt et du
Sherwood en fonction du nombre de Grashof thermique pour les deux Modèles et avec un
Grc=5.105. On constate une parfaite superposition des deux modèles.
138
Chapitre 4 Résultats et discussions
Modèle de Darcy-Brinkman-Forchheimer
Modèle de Darcy
10
9
Nu
8
Grt
Modèle de Darcy-Brinkman-Forchheimer
28 Modèle de Darcy
26
Sh
24
22
20
18
16
14
12
10
8
6
4
2
Grt
139
Chapitre 4 Résultats et discussions
Modèle de Darcy-Brinkman-Forchheimer
Modèle de Darcy
0,254
Nu
0,253
0,252
0,251
0,250
Grt
Modèle de Darcy-Brinkman-Forchheimer
0,40 Modèle de Darcy
0,38
Sh
0,36
0,34
0,32
0,30
0,28
0,26
0,24
Figure 4.62 : Nombre de Sherwood en fonction en fonction de Grt pour Grc=105 et Da=5.10-7.
140
Chapitre 4 Résultats et discussions
141
Conclusion générale
CONCLUSION GÉNÉRALE
Les résultats obtenus nous ont permis d’établir une description de la convection
naturelle bidiffusive, de sa structure d’écoulement et des transferts de chaleur et de masse. Ils
peuvent être résumés comme suit.
Dans le cas d’un milieu fluide, quand les forces thermique et solutale sont
coopérantes, pour des valeurs de Grt inférieur à Grc, l’écoulement obtenu est monocellulaire
tournant dans le sens horaire et il est pratiquement conduit par la force de buoyancy solutale.
Alors que, l’écoulement est dominé par la force de buoyancy thermique pour des valeurs de
Grt supérieur à Grc et sa structure est multicellulaire pour Grt≥107.
142
Conclusion générale
dans le sens horaire. Les transferts de chaleur et de masse déterminés par les valeurs des
nombres moyens de Nusselt et Sherwood sont plus élevés dans le cas coopérant.
Dans cette étude, les modèles mathématiques utilisés sont basés sur de nombreuses
hypothèses simplificatrices. Il serait intéressant de modifier ces hypothèses de façon à tenir en
compte certains effets, tels que les effets tridimensionnels et la variation des propriétés
physico-chimiques du fluide. En outre, le milieu poreux était supposé homogène et isotrope
et les interactions entre les transferts de chaleur et de masse étaient négligées. Donc, l’étude
de l’influence des effets d’anisotropie et des effets de Soret et de Dufour, est aussi fortement
recommandée.
143
Références bibliographiques
RÉFÉRENCES BIBLIOGRAPHIQUES
1. ABDEDOU A., Transferts convectifs dans un échangeur plan contenant un milieu poreux
réactif, en absence d’équilibre thermique local. Mémoire de Magister, école militaire
polytechnique EMP, Borj El Bahri, Alger, 2005.
3. ALAM M., RAHMAN M., SAMAD M., Numerical study of the combined free-forced
convection and mass transfer flow past a vertical porous plate in a porous medium with
heat generation and thermal diffusion. Nonlinear Analysis: Modeling and control, vol. 11,
N° 4, pp. 331-343, 2006.
5. BAHLOUL A., YAHIAOUI M. A., VASSEUR P., BENNACER R., BEJI H., Natural
convection of a two-component fluid in porous media Bounded by tall concentric vertical
cylinders. Journal of Applied Mechanics, vol. 73, pp. 26-33, 2006.
6. BAHLOUL A., VASSEUR P., ROUBILLARD L., Convection of binary fluid saturating a
shallow porous cavity subjected to cross heat fluxes. J. Fluid Mech., vol. 574, pp. 317-
342, 2007.
7. BARLETTA A., LAZZARI S., 2D free convection in porous cavity heated by an internal
circular boundary. Excerpt from proceedings of the COMSOL Multiphysics User’s
Conference Stockholm, 2005.
8. BEJI H., BENNACER R., DUVAL R., VASSEUR P., Double-diffusive natural
convection in a vertical porous annulus. Numerical Heat transfer, Part A, vol. 36, pp. 153-
170, 1999.
144
Références bibliographiques
9. BENNACER R., BEJI H., DUVAL R., VASSEUR P., The Brinkman model for
thermosolutal convection in a vertical annular porous layer. International Journal of Heat
Mass Transfer, vol.27, N° 1, pp. 69-80, 2000.
10. BENNACER R., TOBBAL A., BEJI H., Convection naturelle thermosolutal dans une
cavité poreuse anisotrope : formulation de Darcy-Brinkman. Rev. Energ. Ren. vol.5, pp.
1-21, 2002.
11. BENNACER R., BEJI H., MOHAMED A. A., Double diffusive in a vertical enclosure
inserted with two saturated porous layers confining a fluid layer. International Journal of
thermal Sciences, vol. 42, pp.141-151, 2003.
12. BERA A., KHALILI P., Double-diffusive natural convection in an anisotropic porous
cavity with opposing buoyancy forces: multi-solutions and oscillations. Int. J. Heat Mass
Transfer, vol. 45, pp. 3205-3222, 2002.
13. BOURICH M., HASNAOUI M., AMAHMID A., A Scale analysis of thermosolutal
convection in a saturated porous enclosure submitted to vertical temperature and
horizontal concentration gradients. Energy Conversion and Management, vol. 45, pp.
2795-2811, 2004a.
14. BOURICH M., HASNAOUI M., AMAHMID A., Double-diffusive natural convection in
a porous enclosure partially heated from below and differentially salted. International
Journal of Heat and fluid flow, vol. 25, pp. 1034-1046, 2004b.
15. BOUSRI A., BOUHADEF K., Etude numérique du transfert de chaleur et de masse dans
un cylindre rempli d’un milieu poreux réactif. Eurotherm seminar N°81, Reactive Heat
Transfer in Porous Media, école des mines d’Albi, France, June 4-6, 2007.
16. CHEN G., HADIM H. A., Numerical study of three-dimensional non-Darcy forced
convection in a square porous duct. International Journal for Numerical Methods for Heat
and Fluid Flow, vol. 9, N° 2, pp. 151-169, 1999.
18. DINA S. R., effet d’un champ Magnétique sur la convection naturelle d’un fluide binaire
en milieu poreux confiné. Thèse de Maîtrise es sciences appliquées, école polytechnique,
université de Montréal, 2006.
145
Références bibliographiques
19. GOBIN D., GOYEAU B., NECULAE A., convective heat solute transfer in partially
porous cavities. Int. J. Heat Mass Transfer, vol. 48, pp. 1898-1908, 2005a.
20. GOBIN D., GOYEAU B., HIRATA S., Natural convection in partially porous media.
EUROTHERM 82, Numerical Heat transfer, september 13-16, 2005b.
21. HOOMAN K., GURGENCI H., A theoretical analysis of forced convection in porous-
saturated circular tube: Brinkman-Forchheimer model. Transp. Porous med., vol. 69, pp.
289-300, 2007.
22. HOOMAN K., GURGENCI H., Effects of dependent viscosity on Bénard convection in
porous medium using a non-Darcy model. Int. J. Heat Mass Transfer, vol. 51, pp. 1139-
1149, 2008.
23. KALLA L., Convection naturelle au sein d’une couche horizontale poreuse soumise à des
flux de chaleur uniformes. Thèse de Maîtrise es sciences appliquées, école polytechnique,
université de Montréal, 1999.
24. KALLA L., GOUJON-DURAN S., VASSEUR P., ROBILLARD L., Double-diffusive
convection in a horizontal porous layer. CIMASI, EHTP Casablanca Maroc, octobre 23-
25, 2000.
25. KHANAFER K., VAFAI K., Double-diffusive mixed convection in a lid-driven enclosure
filled with a fluid-saturated porous medium. Numerical Heat transfer, Part A, vol. 42, pp.
465-486, 2002.
26. KUZNESTSOV A. V., NIELD D. A., The effects of combined horizontal and vertical
heterogeneity on the onset of convection in porous medium: double diffusive case. Transp
Porous Med, vol. 72, pp. 157-170, 2008.
27. MAHIDJIBA A., MAMOU M., VASSEUR P., Onset of double-diffusive convection in a
rectangular porous cavity subject to mixed boundary conditions. Int. Comm. Heat Mass
Transfer, vol. 43, pp. 1505-1522, 2000.
28. MAHIDJIBA A., ROBILLARD L., VASSEUR P., Onset of penetrative convection of
cold water in a porous layer under mixed boundary conditions. Int. J. Heat Mass Transfer,
vol. 49, pp. 2820-2828, 2006.
29. MAMOU M., Convection thermosolutale dans des milieux poreux et fluides confinés.
Thèse de doctorat (Ph. D.), école polytechnique, université de Montréal,1998.
146
Références bibliographiques
31. MHAZRI M., DAGUENET M., DAOUDI S., Thermosolutal natural convection in a
vertically layered fluid-porous medium heated from the side. Energy Conversion and
Management, vol. 41, pp. 1065-1090, 1999.
32. MOTA J. P. B., SAATDJIAN E., On the reduction of natural convection heat transfer in
horizontal eccentric annuli containing saturated porous media. International Journal for
Numerical Methods for Heat and Fluid Flow Vol. 7 No. 4, pp. 401-416, 1997.
33. NAWAF H. S., Natural convection in a square porous cavity with an oscillating wall
temperature. The Arabian Journal for Science and Engineering volume 31, N° 1B, pp.
153-164, 2005.
34. NAWAF H. S., POP I., Periodic mixed convection in horizontal porous layer Heated from
below by isoflux heater. The Arabian Journal for Science and Engineering volume 31, N°
2B, 2006.
35. NIELD D. A., BEJAN A., Convection in porous media. Third edition, Springer, 2006.
36. OUAZAA N., Étude de la convection thermosolutale dans les milieux poreux confinés.
Mémoire de magister, université Mentouri – Constantine, 2008.
37. PAKDEE W., PHADUNGSAK R., Natural convection in porous enclosure caused by
partial heating or cooling. The 20th Conference of Mechanical Engineering Network of
Thailand, 2006.
38. PATANKAR S.V., Numerical Heat transfer and fluid flow. Hemisphere, Washington,
DC, 1980.
40. VASSEUR P., ROBILLARD L., SEN M., The Brinkman model for natural convection in
shallow porous cavity with uniform heat flow. Numerical Heat Transfer, vol. 15, pp. 221-
242, 1989.
147
Références bibliographiques
42. YOUNSI R., HARKATI A., KALACHE D., Numerical simulation of double-diffusive
natural convection in porous cavity: opposing flow. The Arabian Journal for Science and
Engineering, vol. 27, N° 1C, 2002.
43. YOUNSI R., Simulation numérique du transfert de chaleur et de masse en milieux fluides
et poreux. Thèse de doctorat, université des sciences et technologies Houari Boumediene,
Alger, 2002.
148