Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

DJE5329

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 160

RÉPUBLIQUE ALGÉRIENNE DÉMOCRATIQUE ET POPULAIRE

MINISTÈRE DE L'ENSEIGNEMENT SUPÉRIEUR ET DE LA RECHERCHE SCIENTIFIQUE

UNIVERSITÉ DE MENTOURI CONSTANTINE


FACULTÉ DES SCIENCES DE L'INGENIEUR
DÉPARTEMENT DE GÉNIE MÉCANIQUE

N° d’ordre : …... / … / 2009


Série : ……. / GM / 2009

Mémoire

Présenté pour obtenir le diplôme de magister en


Génie mécanique
Option : thermo fluides

Convection naturelle bidiffusive dans une enceinte


rectangulaire remplie d’une matière poreuse

Présenté par :
DJEBLI Ahmed

Soutenu le : 05 / 05 / 2009

Devant le jury composé de :

Président : M. BESSAIH Rachid Prof. Université Mentouri Constantine


Rapporteur : M. BENISSAAD Smail M.C. Université Mentouri Constantine
Examinateurs : M. BELLAOUAR Ahmed Prof. Université Mentouri Constantine
M. TALBI Kamel M.C. Université Mentouri Constantine
A mes Parents.
Pour leur dévouement, leurs sacrifices, leurs endurances et leurs souffrances pour moi.
A mes frères et à ma sœur.
En leur souhaitant beaucoup de bonheur.

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.

Plus particulièrement, je remercie mon encadreur Monsieur BENISSAAD Smail pour


m’avoir suivi avec patience et intérêt et pour la confiance qu’il a placé en moi tout au long de
ce travail. Ses conseils précieux et ses encouragements m’ont été d’une aide très précieuse
dans la réalisation de ce projet.

Je tiens à exprimer mes vifs remerciements à Monsieur le professeur BESSAIH


Rachid pour avoir accepté de présider le jury d’examen de ce mémoire. Je remercie également
Monsieur le professeur BELLAOUAR Ahmed et Monsieur TALBI Kamel pour l’honneur qui
m’ont fait d’avoir accepté de faire partie du jury et d’avoir consacré de leur temps à la lecture
de ce mémoire.

Je tiens aussi à remercier mes amis Fatouh, Bachir et Mehdi pour leur aide

ii
RÉSUMÉ

Cette étude concerne la simulation numérique de la convection naturelle bidiffusive,


coopérante et opposante, dans une enceinte rectangulaire bidimensionnelle ayant un rapport
d’aspect (longueur/ hauteur) égal à 4. L’enceinte est remplie d’une matière poreuse homogène
et isotrope. Les parois horizontales de l’enceinte sont adiabatiques imperméables. Les parois
verticales sont maintenues à des températures et des concentrations constantes et différentes.
Le fluide est l’air et la substance diffusée dans l’air de l’enceinte est la vapeur d’eau. Les
conditions aux limites de températures et de concentrations entraînent la convection naturelle
bi-diffusive. L’écoulement dans le milieu poreux de l’enceinte est modélisé par le modèle de
Darcy-Brinkman-Forchheimer. La méthode numérique des volumes finis est utilisée pour la
résolution des équations de continuité, de quantité de mouvement, de l’énergie et de la masse.
Un cas de la convection naturelle bidiffusive dans un milieu fluide a été étudié et considérés
comme référence. Les résultats obtenus ont permis l’étude de l’influence de certains
paramètres sur l’écoulement et sur les transferts de chaleur et de masse, dans les deux cas
coopérant et opposant, du système étudié. Une comparaison entre le modèle de Darcy et le
modèle de Darcy-Brinkman-Forchheimer est réalisée pour déterminer les effets des termes de
Brinkman et Forchheimer.

Mots clés : Milieu poreux, Transfert de chaleur et de masse, Modèle de Darcy-Brinkman-


Forchheimer, Convection naturelle, Volume finis.

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‫ ا‬VPI”Rfk‫م اﻝ‬OPlL‫ ﻝ‬VIk–z‫ اﻝ‬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 zIgšPPM VPPW‫را‬YS ™kPP`M _PP”ILt VLPPQlk‫˜ اﻝ‬e_PPRf‫ اﻝ‬.yPP‫ﺝ‬zk‫\[ آ‬PP‫وأﺥ‬
›IPPS VPPc‫ ﻡƒ_ر‬.‫ف‬zPPUk‫_م اﻝ‬PPf‫^ اﻝ‬PP} VPPc‫_و‬URk‫_دة واﻝ‬PPpRk‫› اﻝ‬IR‫_ﻝ‬PPl‫^ اﻝ‬PP} ‫_دة‬PPk‫ارة واﻝ‬zPPl‫_ر اﻝ‬PPbRc‫ا‬
›I‫ی‬zPQfU‫ اﻝ‬zIgšPM VPW‫ درا‬jP‫› اﺝ‬P‫\ت ﻡ‬Prc zkIP‫رﺵ‬O}-‫_ن‬kvczS-^W‫ذج دار‬Okc‫^ و‬W‫ذج دار‬Okc
.zkI‫رﺵ‬O} ‫_ن و‬kvczS

،zkIP‫رﺵ‬O} -‫_ن‬kvczS -^W‫ذج دار‬Okc ،‫_دي‬k‫اري واﻝ‬zl‫ƒ_ل اﻝ‬Rcd‫ ا‬،\}_c žW‫ و‬: ‫ت ا "ا ــ‬$%&‫آ‬
.VI”Rfk‫اﻝ‬ ‫م‬Ol‫ اﻝ‬،^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

A rapport d'aspect géométrique = L/H


Coefficient dans l’équation de discrétisation
|| fonction du schéma numérique en fonction du nombre de Peclet
C concentration adimensionnelle.
Cf coefficient du terme de Forchheimer (Cf = 0.55)
D coefficient de diffusion [m2.s-1]
dX distance adimensionnelle entre deux centres de volume consécutifs dans la direction X
dY distance adimensionnelle entre deux centres de volume consécutifs dans la direction Y
g accélération de pesanteur [m.s-2]
H hauteur de l’enceinte [m]
K perméabilité du milieu poreux [m/s]
k conductivité thermique [W·m-1·K-1]
L longueur de l’enceinte [m]
P pression adimensionnelle
S terme source
t temps adimensionnel
T température adimensionnelle
u vitesse horizontale [m/s]
U vitesse horizontale adimensionnelle
v vitesse verticale [m/s]
V vitesse verticale adimensionnelle
x coordonnée horizontale dimensionnelle [m]
X coordonnée horizontale adimensionnelle
y coordonnée verticale dimensionnelle [m]
Y coordonnée verticale adimensionnelle

Symboles Grecs

φ Porosité du milieu poreux


α Coefficient de diffusion thermique [m2.s-1]
λ rapport des conductivités thermiques

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

nombre de Graschof solutal 


C ∆H

Grc

nombre de Graschof thermique 


T ∆H

Grt

Sh nombre de Sherwood
 nombre de Sherwood moyen

nombre de Lewis 
!
"
Le

rapport des forces de volume solutale et thermique 


#$ ∆
#% ∆
N

Nu nombre de Nusselt
&' nombre de Nusselt moyen

nombre de Prandtl 
(
!
Pr

nombre de Rayleigh thermique 


T ∆H
)
Rat

nombre de Schmidt 
(
"
Sc

x
Exposant

n relatif au temps t
n+1 relatif au temps t + ∆t

Indices

P point au centre du volume fini


E point à l’est du point P
W point à l’ouest du point P
N point au nord du point P
S point au sud du point P
référence au terme source S
f référence au fluide.
m référence au mélange.
s référence au solide.
max minimum des valeurs à X=0 et X=4
min maximum des valeurs à X=0 et X=4
0 état de référence
i indice d’un point du maillage suivant x
j indice d’un point du maillage suivant y
référence à la variable générale dépendante

xi
Chapitre 1 Introduction

1. INTRODUCTION

1.1 GÉNÉRALITÉS

La présence de gradients de température et de concentration dans un milieu poreux


saturé entraîne l’apparition d’un écoulement de convection naturelle avec transport de chaleur
et de masse. En effet, les différences de température et de concentration provoquent une
distribution non uniforme de la densité du milieu, ce qui donne naissance au mouvement du
fluide sous l’effet de la gravité. Ce phénomène est appelé convection thermosolutale.

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.

L’étude de ce mode de convection en milieu poreux s’est particulièrement


développée ces dernières années vu son intérêt pour la résolution de certains problèmes
technologiques importants tel que le changement de phase où la zone diphasique, dite zone
pâteuse, est représentée par un milieu poreux saturé par un fluide binaire. Il est donc pertinent,
pour définir la distribution des transferts dans de tel système, de caractériser les écoulements
de convection thermosolutale dans le milieu poreux (Danquigny 2003).

Le travail présenté dans ce mémoire a pour but l’étude des mouvements de


convection naturelle à double diffusive se développant dans une cavité horizontale, remplie
par un milieu poreux saturé par un fluide. Les gradients de température et de concentration
sont engendrés par des conditions limites thermique et solutale, imposées sur les parois

1
Chapitre 1 Introduction

actives du système. Bien entendu, la configuration géométrique de l’enceinte et les conditions


limites imposées à ses frontières ont une influence majeure sur la structure de l’écoulement et
sur les taux de transfert de chaleur et de matière résultants.

L’investigation porte essentiellement sur les effets de la variation des paramètres de


contrôle du problème sur la structure de l’écoulement et sur les transferts de chaleur et de
matière au sien du système.

1.2 CONTENU DU MÉMOIRE

Le premier chapitre a été consacré à la présentation du sujet et à une recherche


bibliographique des divers travaux réalisés sur la convection bidiffusive dans les milieux
poreux en général. Les différentes études théoriques, numériques et expérimentales qui ont été
passées en revue, nous ont permis d’assimiler les concepts fondamentaux du problème étudié.

La définition de la géométrie du problème étudié, sa modélisation mathématique,


basée sur le modèle de Darcy-Brinkman-Forchheimer, la mise sous forme adimensionnelle
des équations gouvernantes et des conditions aux limites et la mise en évidence des grandeurs
caractéristiques qui servent de paramètres de contrôle du problème sont abordés dans le
deuxième chapitre.

Le troisième chapitre est consacré à l’exposition de la méthode des volumes finis


utilisée pour la discrétisation des équations gouvernantes et l’algorithme SIMPLER employé
pour la résolution du système d’équations résultant de cette discrétisation.

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

1.3 REVUE BIBLIOGRAPHIQUE

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

réduction importante de l'échange de chaleur moyen (Numoy). Le rapport des forces de


buoyancy N, est le principal paramètre dirigeant de la convection thermosolutale par son
ampleur et son signe. En particulier, quand N= –1 et Le = Rk = Rd = 1, l’écoulement
convective dans la cavité est complètement réprimée. Par conséquent, dans ce cas particulier,
la moyenne du nombre de Nusselt (ou de Sherwood) en fonction N est une courbe en forme
de V symétrique par rapport à la valeur N = –1 (Numoy (│N │) = Numoy (–│N│) = Shmoy
(│N│) = Shmoy (–│N│)). De plus, pour Ra=104, quand –3 ≤ N ≤ – 0.5 et si au moins un des
trois paramètres, Rk, Rd et Le diffère de l'unité, un écoulement multicellulaire peut être
observé dans la cavité. En outre, la structure et l'intensité de ce fluide de convection
multicellulaire dépendent des valeurs d'autres paramètres gouvernants.

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 modèle Darcy avec l'approximation Boussinesq a été employé pour étudier la


convection naturelle à double diffusive dans une cavité peu profonde poreuse par Kalla et al.
(2000). Les parois horizontales sont soumises à des flux uniformes de chaleur et de masse,
tandis que les parois verticales sont exposées à un flux de chaleur. La méthode de différence
finie a été employée pour obtenir la solution numérique des équations gouvernantes. Les
équations de l'énergie et de concentration sont résolues avec la méthode ADI. Le champ de la
fonction de courant est obtenu avec la méthode SOR (over relaxation method). Un bon accord
est observé entre une solution approximative analytique et les simulations numériques. Les
résultats analytiques indiquent que pour la convection supercritique ou les deux contributions
solutale et thermique se déstabilisent, il y a trois solutions différentes. Une de ces solution est
instable pour une gamme de donnés des paramètres gouvernants, à condition que Rat soit plus
grand que la valeur critique qui dépend de Ras et Le. D'autre part, pour la convection
souscritique, qui se produit quand les deux forces de buoyancy sont opposées, il y a cinq
solutions différentes sont possibles (dont deux sont instables). L'existence des solutions
multiples quand Rat = 0 a aussi été démontrée.

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.

La convection naturelle avec le transfert de chaleur et de masse et avec des gradients


de chaleur et de concentration horizontaux et opposé a été examinée par Bera et Khalili
(2002) dans une cavité anisotrope poreuse utilisant le modèle Darcy. Le milieu poreux est
supposé être hydrodynamiquement et thermiquement anisotrope. Les directions principales
du tenseur de perméabilité sont prises obliques au vecteur de gravité, tandis que ceux des
diffusivités thermique et solutale coïncident avec les axes de coordonnée horizontaux et
verticaux. Les auteurs ont donné une grande attention pour comprendre l'effet de paramètres
anisotropes sur l'existence d'oscillations instables permanentes et des solutions multiples
stationnaires. De l'étude de solutions analytiques, qui peuvent être considérées comme une
vérification des résultats numériques, simultanément, ils ont trouvé qu'il existe un intervalle
de rapport de buoyancy INM, selon les valeurs paramétriques, dans lesquelles des solutions
multiples existent. Pour le cas instable un intervalle semblable, INO, pour le rapport de
buoyancy a été observé numériquement, dans lequel des oscillations permanentes existent. La
périodicité de l'oscillation change brutalement en changeant la perméabilité du milieu. Les
résultats indiquent que le maximum INM et INO des intervalles est atteint à un angle
d'orientation de θ=45°. La direction locale d’écoulement change à cause de la variation des
couches thermique et de concentration et des paramètres anisotropes. Il a été observé qu'une
petite rotation du tenseur de perméabilité fait un changement significatif de la force
d’écoulement aussi bien que les profiles de la température et de la concentration. En général,
la chaleur et la masse suivent le modèle complexe selon l'interaction entre les coefficients de
diffusion et le rapport de buoyancy.

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

thermique et solutale. Le transfert de chaleur et de masse diminue avec la diminution de la


perméabilité.

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

le cœur sont linéairement stratifiées dans la direction verticale. La solution analytique


développée confirme ces observations.

L’étude de Bourich et al. (2004a) présente une analyse d’échelle de la convection


naturelle à double diffusive bidimensionnelle dans une cavité carrée poreuse chauffée plus
tard. La surface supérieure est froide par contre les limites verticales sont soumises à un
gradient solutal horizontale. Dans la formulation du problème, le modèle de Darcy est
considéré avec l’approximation de Boussinesq. L’analyse est établie en se basant sur
l’approche nommée « scaling law » pour prédire l’ordre de l’importance des caractéristiques
des transferts de masse et de chaleur dans des cas limites qui correspondent aux écoulements
dirigés par les forces thermique ou solutale. Les équations de conservation sont résolues
numériquement avec la méthode des différences finies. La procédure itérative utilisée est
l’ADI et les valeurs de la fonction de courant sont calculées par la méthode SOR. Les résultats
de l'analyse d'échelle ont été validés numériquement et les auteurs ont proposé les corrélations
utiles de Nu et Sh dans ces cas limites.

Bourich, et al. (2004b) ont étudié la convection naturelle bidimensionnelle


bidiffusive dans une cavité carrée poreuse, partiellement chauffée dont la surface supérieure
est refroidie à une température constante. Les parois verticales de la matrice poreuse sont
soumises à un gradient de concentration horizontal. Les équations gouvernantes sont
discrétisées selon le schéma des différences centrales finies. La procédure itérative est
exécutée employant la méthode ADI. Le champ de la fonction de courant est calculé avec une
la méthode SOR. Les paramètres gouvernants le problème sont le nombre de Rayleigh
thermique (Ra = 100 et 200), le nombre de Lewis (Le = 0.1 ; 1 et 10), le rapport de buoyancy
(–10≤N≤10) et la position relative de l'élément chauffant avec la ligne médiane verticale de la
cavité (δ = 0 et 0.5). L'effet des paramètres gouvernants sur les caractéristiques liquides a été
analysé. Les résultats indiquent que l'effet de N sur le comportement dynamique du fluide et
les transferts de chaleur et de masse a été trouvé dépendant du type de solution, des nombres
de Lewis et de Rayleigh et de la position de la source de chauffage isotherme. L'existence
d’écoulement bicellulaire asymétrique est prouvée dans le cas de δ = 0.5. Les transferts de
chaleur et de masse dans la cavité poreuse dépendent de la solution considérée, quand des
états multiples stables sont possibles. Les solutions bicellulaires réduisent considérablement le
transfert de chaleur comparé à aux solutions monocellulaires. La multiplicité des solutions
disparaît en présence de gradient solutal horizontal quand des conditions critiques sont
atteintes.

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

température dépendent du nombre de Rayleigh modifié Ra qui est le produit du nombre de


Rayleigh Ra et du nombre de Darcy Da. La moyenne du nombre de Nusselt et la moyenne de
l'énergie cinétique Ekin ont été évaluées. Les résultats obtenus montrent que pour les valeurs
données du rapport d’aspect A et l’angle d’inclination de cavité α, Nu et Ekin augmentent avec
Ra. De plus, pour des valeurs données de A et de Ra, Nu et Ekin augmentent aussi avec α.
D’autre part, pour les valeurs données de α et de Ra, Nu diminue avec la diminution de A
tandis que Ekin augmente.

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

Une étude numérique bidimensionnelle de la convection naturelle dans une cavité


carrée avec un milieu poreux et avec une paroi chauffée a été présentée par Nawaf (2005). Le
modèle de Darcy est utilisé et le fluide est supposé un standard Boussinesq. La température de
la paroi chauffée est oscillante avec le temps et avec une valeur constante. Par contre, la paroi
froide est mise à une température constante. Les parois horizontales sont adiabatiques. Les
résultats montrent la variation temporelle des lignes de flux, les isothermes et le nombre de
Nusselt. Quand la température de la paroi chauffé oscille avec des grandes amplitudes et
fréquences, le nombre de Nusselt devient négative et ceci est due à l’insuffisance du temps
donnée pour le transfert de paroi chauffé jusqu’à la paroi froide.

Le problème de combinaison de la convection libre et forcée et du transfert de masse


dans une plaque verticale poreuse semi infinie, en présence de génération de chaleur et de
diffusion thermique, a été étudié numériquement par Alam et al. (2006). Les équations
différentielles non linéaires et leurs conditions aux limites, qui caractérisent le problème, sont
transformées à un système d’équations différentielles ordinaires. Ce système est résolu
numériquement par l’utilisation de la technique d’itération de Nachtsheim-Swigert avec le
schéma d’intégration Runge-Kutta d’ordre six. Les résultats indiquent que les profiles de la
vitesse et de la température augmentent et le profile de la concentration diminue avec
l’augmentation du paramètre de génération de chaleur. Aussi, les profiles de la vitesse et de la
concentration augmentent avec l’augmentation du nombre de Soret. Dans le régime de
convection mixte le nombre de Sherwood augmente. Par contre le nombre de Nusselt diminue
avec la croissance du paramètre de génération de chaleur et du nombre de Soret.

Bahloul et al. (2006) ont étudié analytiquement et numériquement le comportement


d'un mélange binaire saturant un milieu poreux dans un annulaire vertical. Des flux de chaleur
uniformes ont été appliqués aux parois verticales, tandis que les parois horizontales sont
imperméables et adiabatiques. Des gradients solutaux constants ont été appliqués sur les
parois verticales pour obtenir la convection à double diffusion avec ou sans l’effet de Soret.
Les paramètres gouvernants sont le nombre de Rayleigh thermique Rat, le rapport de
buoyancy N, le nombre de Lewis Le, le rapport d'aspect A et la constante a. Dans la gamme
des paramètres gouvernants considérés dans cette étude, un bon accord a été trouvé entre les
prédictions analytiques et les résultats numériques obtenus. Dans les conditions aux limites
thermique et solutale imposées aux parois verticales de l’enceinte, il a été démontré
numériquement que l’écoulement dans la cavité est presque parallèle à condition que le
rapport d'aspect A soit assez grande. Le modèle d’écoulement parallèle proposé dans cette

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.

Mahidjiba et al. (2006) ont étudié numériquement et analytiquement le début de la


convection dans une couche horizontale bidimensionnelle poreuse saturée avec de l'eau
froide. La méthode des éléments finis est employée pour résoudre les équations linéaires de
perturbation. Le début de convection est dépendant du rapport d'aspect de la cavité A et du
paramètre d'inversion γ (défini par 2(Tmax-T)/∆T). L'effet de ces paramètres de contrôle est
étudié pour des cas de couches infinies et finies. L'effet des diverses conditions aux limites
thermiques imposées à la couche horizontale a aussi été examiné dans cette étude. Des
nombres de Rayleigh critiques définissant le seuil de convection ont été établis en fonction de
γ dans le cas de la couche infinie. Ces résultats indiquent aussi qu'une situation asymptotique
est atteinte quand γ≤1. Il a aussi été trouvé, que pour γ>2, les solutions atteignent
asymptotiquement le cas classique d'un liquide avec une relation linéaire entre la température
et la densité.

L'étude numérique de la convection transitoire mixte dans une couche


bidimensionnelle horizontale poreuse chauffée par une source de flux de chaleur constante a
été effectuée par Nawaf et Pop (2006). Les équations soumises aux conditions aux limites
initiales sont intégrées numériquement en utilisant la méthode des volumes finis. Le champ
thermique transitoire, le champ d’écoulement et la moyenne du nombre de Nusselt sont
présentés en fonction du nombre de Péclet Pe. Pour le cas particulier de nombre de Rayleigh
Ra = 102 et le rapport de longueur d'appareil de chauffage à l'épaisseur de la couche poreuse
A, il a été trouvé que pour A = 3 et A = 5 avec des petites valeurs du nombre de Péclet, le
mode de convection naturelle domine, tandis que pour de grandes valeurs du nombre de
Péclet, le mode de convection forcée domine. Cependant, pour des valeurs modérées la
convection oscillante mixte est observée et une variation périodique de la moyenne du nombre
de Nusselt est obtenue. Quand la longueur d'appareil de chauffage est égale à l'épaisseur de
couche poreuse (A=1) des résultats stationnaires sont obtenus pour la gamme de Pe = 0.1 - 10.

13
Chapitre 1 Introduction

Pakdee et Rattanadecho (2006) ont étudié numériquement l’écoulement convectif


naturel transitoire dans un milieu poreux saturé avec un liquide dans une cavité rectangulaire.
La cavité est isolée sauf la paroi supérieure qui est partiellement exposé à un extérieur
ambiant. La surface exposée permet le transport convectif par le milieu poreux, produisant
une stratification thermique et des circulations de flux. Le système d'équations différentielles
est résolu avec des conditions aux limites et initiales appropriées employant la méthode des
différences finies. Les résultats obtenus montrent qu’un gradient latéral de température près
de la paroi supérieure pousse la force de buoyancy dans une condition instable. L’écoulement
bidimensionnel est caractérisé principalement par deux tourbillons symétriques qui sont
produits par l'effet de buoyancy. Les directions de rotation des tourbillons produits dans deux
conditions différentes sont opposées. Il a été trouvé aussi que, le coefficient de transfert de
chaleur, le nombre de Rayleigh et le nombre de Darcy influencent considérablement sur les
caractéristiques des mécanismes de transfert de chaleur et d’écoulement. Les grandes valeurs
du nombre de Rayleigh augmentent les intensités des lignes de courants, augmentant ainsi la
pénétration d’écoulement en bas. Les petites valeurs du nombre de Darcy gênent les
circulations d’écoulement, donc le transfert de chaleur par convection est considérablement
réduit.

Bahloul et al. (2007) ont étudié analytiquement et numériquement la convection


naturelle dans un mélange binaire chauffé différentiellement. Le fluide est contenu dans une
cavité peu profonde rectangulaire poreuse et soumise à l’effet de Soret. Toutes les quatre
faces sont exposées à des flux de chaleur uniformes, les faces opposées étant chauffées et
froidies, respectivement. Des solutions analytiques pour les champs de fonction de courant, de
température et de concentration sont obtenues employant une supposition d’écoulement
parallèle dans la région principale de la cavité et une forme intégrale de l'énergie et des
équations constitutives. La confirmation numérique des prédictions analytiques est aussi
obtenue. L’approche de volume de contrôle est utilisée pour la discrétisation des équations
gouvernantes. L’équation de la fonction de courant a été résolue par la méthode ADI et les
équations de température et de concentration on été résolue par la méthode. Les résultats sont
présentés d'abord en présence d'un gradient vertical de température. Le cas spécial où les
forces de buoyancy sont opposées et d'intensité égale (N = –1) a été aussi discuté. Pour cette
situation particulière, le nombre de Rayleigh supercritique pour le début de convection a été
prévu sur la base d'une analyse de stabilité linéaire. Des états multiples stables près du seuil
de convection on été trouvés.

14
Chapitre 1 Introduction

Une étude numérique du transfert de chaleur et de masse dans un cylindre rempli


d’un milieu poreux réactif a été réalisée par Bousri et Bouhadef (2007). Leur analyse était
basée sur le modèle de Darcy-Brinkman et les équations de conservations, qui décrivent
l’évolution spatiale et temporelle des phénomènes couplés des réactions chimiques de
diffusion et de transfert thermique. L’approche est faite à partir d’un modèle 2D. Le système
des équations gouvernantes obtenu est discrétisé et traité par la méthode des volumes finis.
Les résultats obtenus indiquent que la réduction de perméabilité accentue la filtration.
L’augmentation du nombre de Reynolds accroît l’inertie du fluide et réduit le temps de
transition. L’augmentation du terme source lié à la réaction chimique augmente la température
et diminue la concentration à la sortie.

La convection forcée entièrement développée à l'intérieur d'un tube circulaire rempli


d’un milieu poreux saturé soumis à un flux de chaleur uniforme sur la paroi, a été examinée
sur la base du modèle Brinkman-Forchheimer par Hooman et Gurgenci (2007). Pour les
grands et petits nombres de Darcy, la solution pour l'équation de mouvement Brinkman-
Forchheimer a été trouvée en termes d'une méthode asymptotique. Une fois que la distribution
de vitesse est déterminée, l'équation d'énergie est résolue employant la même technique
asymptotique. Les résultats obtenus montrent que le profil de vitesse dépend fortement du
paramètre s (s = [M.Da]-1/2). Avec l’augmentation de ce paramètre, la région centrale
contenant une distribution relativement uniforme de la vitesse s'étend plus loin vers la paroi.
Pour les grandes valeurs de s, le profil de vitesse est limité à une couche très mince adjacente
à la paroi. Pour les grandes valeurs de s, l'influence du coefficient de Forchheimer Cf et du
coefficient M (M = µeff/µ) sur les résultats est négligeable. Cependant, pour les petites valeurs
de s, le nombre de Nusselt croît avec l'augmentation de Cf et diminue avec l’augmentation de
M.

Marcelo et al. (2007) Ont étudié analytiquement les équations macroscopiques de


transport de chaleur et de masse pour l’écoulement turbulent dans des structures perméables.
Le mécanisme de la convection naturelle à double diffusive est examiné pour la phase liquide
dans le régime turbulent. Les équations sont basées sur le concept de double décomposition,
qui considère, les fluctuations de temps et les déviations spatiales des valeurs moyennes. Ce
travail avait l'intention de démontrer que des mécanismes complémentaires de transport sont
mathématiquement tirés si la température, la concentration et la vitesse présentent
simultanément des fluctuations de temps et des déviations spatiales dans le domaine
d'analyse. L'analyse de stabilité des mélanges, pour des composants plus lourds sous les

15
Chapitre 1 Introduction

gradients de la température et de la concentration, a été discutée. Les résultats montrent que


des écoulements hydrodynamiquement stables sont obtenus sous certaines distributions de
température et de concentration qui refroidissent la turbulence et mènent finalement à un
processus de relaminarisation. Des situations inconditionnellement instables ont été aussi
étudiées.

La température dépendant de l'effet de la variation de la viscosité dans la convection


Bénard, d'un gaz ou un liquide dans une enceinte remplie d'un milieu poreux, a été étudiée
numériquement par Hooman et Gurgenci (2008). Cette étude est basée sur le modèle général
de transfert de mouvement dans un milieu poreux. Des solutions numériques des équations
gouvernantes sont obtenues par la méthode de différences finies, employant la technique de
Gauss-Seidel avec la méthode SOR. Le modèle Arrhenius qu’est appliqué, propose une forme
exponentielle de la relation entre la température et la viscosité : constant (µ=µC), diminuant
(jusqu'à 0.13µC) et augmentant (jusqu'à 7.39 µC). Les effets de la variation de viscosité du
liquide sur les isothermes, les lignes de courant et le nombre de Nusselt ont été étudiés. La
définition d'une température de référence, qui ne change pas avec le nombre de Rayleigh,
mais augmente avec le nombre de Darcy, a été trouvée comme une variable pour représenter
la dépendance de la variation de la viscosité en fonction de la température. Il a été trouvé
aussi que la température de référence, à laquelle les propriétés des liquides doivent être
évaluées, est approximativement indépendante des autres paramètres considérés.

Les effets d'hétérogénéité hydrodynamique horizontale et verticale, thermique et


solutal, sur le début de convection dans une couche horizontale d'un milieu poreux saturé
uniformément chauffé par le bas, ont été étudiés analytiquement par Kuznestsov et Nield
(2008) employant la théorie de la stabilité linéaire, pour le cas d'hétérogénéité faible. Une
température et un flux de concentration uniformes sont imposés aux parois supérieure et
inférieure. Les parois latérales verticales sont supposées isolées et imperméables. Le modèle
Brinkman est employé. Il a été trouvé que l'effet d'une telle hétérogénéité sur la valeur critique
du Nombre de Rayleigh Ra basé sur des propriétés moyennes est du deuxième ordre si les
propriétés varient d’une manière linéaire. Les effets d'hétérogénéité horizontale et verticale
sont alors comparables, une fois que le rapport d'aspect est pris en compte et à une première
approximation, sont indépendant.

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

2.1 DESCRIPTION DU PROBLÈME

Nous allons considérer la simulation numérique de la convection naturelle à double


diffusion dans milieu poreux contenu dans une enceinte rectangulaire bidimensionnelle. La
géométrie du problème considéré est illustrée sur la figure 2.1. Il s’agit d’une enceinte
rectangulaire de longueur L et de hauteur H, ayant un rapport d’aspect égal à 4. L’enceinte est
remplie d’une matière poreuse homogène et isotrope. Les parois verticales de l’enceinte sont
supposées soumises à des températures et concentrations constantes mais différentes, tandis
que les parois horizontales sont considérées comme adiabatiques et imperméables. Le fluide
est l’air et la substance diffusée dans l’air de l’enceinte est le vapeur d’eau. Les gradients de
température et de concentration résultants des conditions aux limites produisent la convection
naturelle à double diffusion.

Les forces de buoyancy, thermique et solutale, sont considérées soient coopérantes


ou opposées en fonctions des conditions aux limites imposées à l’enceinte étudiée.

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

Figure 2.1 : Description de la géométrie du problème.

18
Chapitre 2 Modélisation physique et mathématique

2.2 HYPOTHÈSES SIMPLIFICATRICES

Pour une formulation simple du modèle mathématique, nous considérons les


approximations, les plus couramment employées dans type de problème, entre autres les
approximations de Boussinesq. Nous supposons alors que :

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)].

2.3 MODÉLISATION MATHÉMATIQUE

En tenant compte de toutes les hypothèses considérées, le système d’équations


différentielles aux dérivées partielles décrivant ce problème est composé de l’équation de
continuité, les équations de quantité de mouvement et les équations du transfert de chaleur et
de matière dans un milieu poreux. Pour tenir compte de tous les effets, on associe à l'équation
de Darcy, les termes de Brinkman et de Forchheimer. Le modèle ainsi formé est connu sous le
nom du modèle de Darcy-Brinkman-Forchheimer ou loi de Darcy étendue. Les équations
modélisantes à travers un milieu poreux s'écrivent alors sous les formes suivantes (Bejan et
Nield 2006).

19
Chapitre 2 Modélisation physique et mathématique

2.3.1 Equation de continuité

*
+,-
12 345 6
/0 12 7  0
+.
(2.1)

Fluide incompressible :

/ 0
+9 +:
+ +;
(2.2)

2.3.2 Equations de quantité de mouvement dans les milieux poreux

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

GJH GNH GRH GUH


G=H GLH

(2.3)

a. (1) 45 V> +. / > 102


12 127
12?
W : Forces inertielles macroscopiques.
= +? = 3?
>

12 G*IH : Gradient de pression.


b.(2) > 6
=

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).


f. (6) 45 S2 : Forces de volume.

2.3.2.1 L’équation de la quantité de mouvement suivant x

/ > / > 
+
 '
5 ' / > V+ [+ \ / +; [];\W
= +9 = +G99H = +G:9H +X K √9 Z: K + +9 + +9
> +. + +; √
(2.4)

2.3.2.2 L’équation de quantité de mouvement suivant y

/ / 

 ^
5 ^ / > V+ [+ \ / +; [+;\W /
= +: = +G9:H = +G::H +X K √9 Z:  K + +: + +:
> +. > + > +; +; √

45 S _` G
a H / ` G
a Hb (2.5)

2.3.3 Equation de l’énergie

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

Dans notre cas il n’y a pas de génération de chaleur donc : e


fff
 0.

3,jk 7
On divise l’équation par 34cX 75 et on considère que : 3,jk 7
l
m 1.
-

Avec :

i. G4cH : capacité calorifique du mélange.


ii. G4cHh : capacité calorifique de la matrice solide.
iii. G4cH5 : capacité calorifique du fluide.
iv. e
fff
: source d’énergie.
v. 34cX 7 : capacité calorifique à pression constante.

 o : coefficient de diffusion thermique.


n-
3,jk 7
vi.
-

 p : rapport des conductivités thermiques.


nl
n-
vii.

+
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)

2.3.4 Equation de transfert de matière


+
12 {12  0J G| H
* +. /6 (2.9)

Dans notre cas on peut écrire :

/> />  V+ [| + \ / +; [| +;\W


+ = +G9H = +G:H + + + +
+. + +;
(2.10)

(Avec Dm= * D : diffusivité massique du milieu poreux)

2.4 MISE EN FORME ADIMENSIONNELLE DES ÉQUATIONS


2.4.1 Variables adimensionnelles

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

c. t’ = t / (H2/ν) ; pour la simplicité on prend t’≡ t.

d. U = u / (ν/H)

e. V = v / (ν/H)

f. P = (p – p0) / [ρf / (ν/H)]

g. Θ  G
a H⁄∆ ; (∆  
 ) ; pour la simplicité on prend Θ ≡ T.

h. Φ  G
a H⁄∆ ; (∆  
 ) ; pour la simplicité on prend Φ ≡ C

2.4.2 Equations adimensionnelles


2.4.2.1 Equation de continuité

/ +;  0
+€ +?
+
(2.11)

2.4.2.2 Equation de la quantité de mouvement suivant x

/ > / > 

"
"- √… J / 6 J . … / > V+ [+ \ / +‚ [ +‚ \W
= +€ = +G€€H = +G?€H +ƒ €  = + +€ + +€
> +. + +‚
(2.12)
„ „

2.4.2.3 Equation de quantité de mouvement suivant y

/ > / > 

"
"- √… J / 6 J . 6 / > V+ [+\ / +‚ [+‚\W /
= +? = +G?€H = +G??H +ƒ ?  = + +? + +?
> +. + +‚ „ „

‡ˆ.  / ‡ˆj  (2.13)

2.4.2.4 Equation de l’énergie

/ /  P VrX [λ rX\ / rY [λ rY\W


rT rGUTH rGVTH = r rT r rT
rt rX rY
(2.14)


2.4.2.5 Equation de transfert de matière

/> />  i V+ [| +\ / +‚ [| +‚\W


+ = +G€H = +G?H = + + + +
+. + +‚
(2.15)
Ž

Cette adimensionalisation nous a permis la mise en évidence de sept paramètres de


contrôle qui caractérise le problème étudié à savoir, respectivement, les nombres de Grashof
thermique Grt et solutal Grc, le nombre de Prandtl Pr, le nombre de Schmidt Sc, le nombre de
Darcy Da, le rapport d’aspect A et la porosité .

22
Chapitre 2 Modélisation physique et mathématique

• Gt 
T ∆H

C ∆H
• G‘ 


• P 
)

• S‘ 
D
K
• D– 
H
• A = L/H
w˜™uš› œ› ž˜› œu šŸ™Ÿ›u ž˜›uv
• φ
w˜™uš› t˜t–™› œu šŸ™Ÿ›u ž˜›uv

2.5 CONDITIONS INITIALES ET AUX LIMITES

La solution du problème non dimensionnel dépend des valeurs des paramètres de


contrôle et des conditions initiales et aux limites adimensionnelles nécessaires pour la
résolution des équations de ce système. Ces conditions sont les suivantes.

2.5.1 Conditions initiales

• à t = 0 ; u = v = 0 et T = C = 0

2.5.2 Conditions aux limites

0;  0;U=V=0
rT rC
rv rv
• Y = 0, X:

0;  0, U = V = 0
rT rC
rv rv
• Y= 1, X:

2.5.2.1 Cas coopérant

• X = 0, Y : T = 1 ; C = 0 ; U = V = 0.

• X = 4, Y : T = 0 ; C = 1 ; U = V = 0.

2.5.2.2 Cas opposant

• 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

2.6 TRANSFERTS THERMIQUE ET MASSIQUE

Les flux moyens de transfert de chaleur et de matière sont exprimés en termes


adimensionnels par les nombres moyens de Nusselt &' et de Sherwood  définis
respectivement par :

¢&'|£.¤  ¥= [¦…
+
\ §
a + £.¤
(2.16)

¢Q  ¥a [c…

= +
\ §
£.¤ + £.¤
(2.17)

24
Chapitre 3 Résolution numérique

3. RÉSOLUTION NUMÉRIQUE

3.1 INTRODUCTION

Au cours de ces dernières années, nous assistons à une augmentation fulgurante de la


performance des ordinateurs, qui a permis l’utilisation de méthodes numériques de plus en
plus performantes et efficaces pour la résolution de systèmes d’équations différentielles non
linéaires.

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.

3.2 PRÉSENTATION GÉNÉRALE

En vue du passage d’un problème exact (continu) au problème approché (discret), on


dispose de plusieurs techniques concurrentes et complémentaires : les différences finies, les
éléments finis et les volumes finis. Chacune de ces trois méthodes correspond à une
formulation différente des équations de la physique :

- équilibre des forces en chaque point pour les différences finies,

- 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.

3.2.1 Différences finies

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

• Avantages : grande simplicité d’écriture et faible coût de calcul.

• Inconvénients : limitation de la géométrie des domaines de calcul, difficultés de prise


en compte des conditions aux limites portant sur les dérivées ou les gradients de
l’inconnue et en général absence de résultats de majoration d’erreurs.

3.2.2 Eléments finis

La méthode des éléments finis consiste à approcher, dans un sous-espace de


dimension finie, un problème écrit sous forme variationnelle (comme minimisation de
l’énergie, en général) dans un espace de dimension infinie. La solution approchée est dans ce
cas une fonction déterminée par un nombre fini de paramètres comme, par exemple, ses
valeurs en certains points (les nœuds du maillage).

• Avantages : traitement possible de géométries complexes, détermination plus naturelle


des conditions aux limites, possibilité de démonstrations mathématiques de
convergence et de majoration d’erreurs.

• Inconvénients : complexité de mise en œuvre et coût en temps de calcul.

3.2.3 Volumes finis

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 MÉTHODE DES VOLUMES FINIS

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 de volumes finis, les algorithmes de résolution des équations de


transport en variables primitives, le schéma upwind pour la discrétisation des termes
d’advection sont quelques idées qui ont été popularisées par Patankar (1980), qui est une
référence incontournable dans ce domaine. Dans ce chapitre, la méthode des volumes finis et
son utilisation pour résoudre des problèmes de transport est expliquée ainsi que l’algorithme
SIMPLER utilisé pour résoudre les équations en variables primitives.

3.3.2 Présentation générale de la méthode

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

Figure 3.1 : Volume fini typique

Volume de contrôle Volume de


centré pour P, C et T contrôle décalé
N N
vn vn
n n

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

Figure 3.2 : Volumes de contrôle typique et décalé.

28
Chapitre 3 Résolution numérique

La pression et la température sont calculées en considérant un volume de contrôle


centré autour du point P, donc entre les faces e, w, s et n en bidimensionnel (fig. 3.3). Les
vitesses sont évaluées dans un volume centré sur une face du volume précédent, par exemple
pour la vitesse horizontale U sur la face est entre deux faces passant par les points P et E et
entre les faces sud et nord.

Physiquement, la température et la pression sont des grandeurs scalaires plus


facilement calculables en un point (supposé représentatif du volume), tandis que la vitesse est
plutôt liée à un flux (que l'on considérera constant) à travers une surface. Vues les hypothèses
effectuées (grandeurs et flux constants), il conviendra de prendre des maillages suffisamment
petits de façon à ce que les moyennes effectuées soient représentatives.

Y
dxw dxe

N ∆xN

dyn n

W w P e E
∆xP

dys s

S
∆xS

∆xW ∆xP ∆xE X

Figure 3.3 : Maillage typique uniforme suivant X et suivant Y.

29
Chapitre 3 Résolution numérique

3.4 FORME GÉNÉRALE DE L’ÉQUATION DE TRANSPORT

Les équations différentielles aux dérivées partielles adimensionnelles (2.11 à 2.15)


présentées dans le chapitre 2, peuvent se mettre sous la forme générale d’une équation de
transport (3.1) dont les expressions des grandeurs physiques, correspondant à chacune des
équations (2.11 à 2.15), sont données dans le tableau 3.1.

/ /  [Γ + \ / +‚ [Γ +‚ \ / «
+« +G€«H +GV«H + +« + +«
+. + +‚ +
(3.1)

Tableau 3.1 : Variable Ф, coefficient de diffusion Г et terme source SФ


pour les équations de notre problème.

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

3.5.1 Intégration l’équation générale de transport

On intègre l’équation générale de transport sur un volume de contrôle et en fonction


du le temps, on trouve alors :

30
Chapitre 3 Résolution numérique

¥. ¥h ¥ §§§ / ¥. ¥h ¥ §§§ / ¥. ¥h ¥ §§§ 


.Z∆. ¤̄ +® .Z∆. ¤̄ +€® .Z∆. ¤̄ +?®
+. + +‚ (3.2)
¥. ¥h ¥ [Γ \ §§§ / ¥. ¥h ¥ [Γ \ §§§ / ¥. ¥h ¥ ® §§§
.Z∆. ¤̄ + +® .Z∆. ¤̄ + +® .Z∆. ¤̄
+ + +‚ +‚

L’intégration des différents termes s’effectue comme suit :

• Terme de la variation temporelle :

¥. ¥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.

• Termes de transport convectif :

¥. ¥h ¥ §§§ / ¥. ¥h ¥ §§§  _G…°H¤


G…°H¯ bΔYΔ /
.Z∆. ¤̄ +€® .Z∆. ¤̄ +?®
+ +‚ (3.2.2)
_G6°H
G6°Hh bΔXΔ

• Termes de transport diffusif :

¥. ¥h ¥ [Γ \ §§§ / ¥. ¥h ¥ [Γ \ §§§ 
.Z∆. ¤̄ + +® .Z∆. ¤̄ + +®
+ + +‚ +‚
(3.2.3)
V[Γ \
[Γ + \ W ΔYΔ / V[Γ +‚ \
[Γ +‚ \ W ΔXΔ
+® +® +® +®
+ ¤ ¯ h

• Terme source :

¥. ¥h ¥ ® §§§  ® ΔXΔYΔ
.Z∆. ¤̄
(3.2.4)

Après l’intégration et la division par Δ l’équation générale de transport devient :

Δ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)

• Le terme source est :

Ä  ƒÅ ΔΔ / aƒ °ƒa (3.4.3)

Où : AaP est le coefficient de la variable ° au temps t.

°ƒa est la valeur de la variable ° calculé au temps t et ƒÅ est le terme source.

• Les flux convectifs sont :

Ǥ  …¤ Δ, ǯ  …¯ Δ ¢
Æ
Ç  … Δ, Çh  …h Δ (3.4.4)

• Les flux diffusifs sont :

|¤  GɁH , |¯  GɁH
¿ È ·‚ À È ·‚
À¢
¸ ¿

|  GɁH , |h  GɁH
ÈÁ ·‚ È ·‚
(3.4.5)
Á Â

• Les nombres de Peclet sont :

¤  ¿ , ¯ 
Ê ÊÀ
"À ¢
¸
"¿ (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

3.5.2 Différentes schémas de discrétisation

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.

Dans notre étude, pour discrétiser les équations différentielles adimensionnelles,


nous avons utilisé le schéma hybride et une discrétisation temporelle d’ordre un.

32
Chapitre 3 Résolution numérique

Tableau 3.2 : Les différents schémas de discrétisation (Patankar 1980).

Schéma Formule de G||H Expressions des différents coefficients

Î  |¤ 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

Schéma couplé de deux schémas


Î  ϬÐG0, 1
0.5|¤ |H. |¤ / ϬÐG
Ǥ , 0H
des différences centrées et upwind
Ñ  ϬÐG0, 1
0.5|¯ |H. |¯ / ϬÐGǯ , 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

3.5.3 Équation de continuité

/ 0
+9 +:
(3.5)
+ +;

L’intégration de l’équation de continuité sur un volume de contrôle décalé (fig. 3.4 et


3.4) donne :

¥¯ ¥h [+ / +‚\ §§  ¢…ƒaÙ Q¯ ∆X / ¢6ƒaÙ Qh ∆¤  (…¤9


¤ +€ +? ¤Ù a )∆ hÙ

…¯9 ƒ / (6 9
6h9 )§¤
Ù Ù a a a
Ù Ù Ù Ù

(3.6)
…Î9
a
/ …ƒ9
a
…ƒ9
a
/ …Ñ9
a
6Îa / 6ƒa 6Îi
a
/ 6ia
Ú Ü ∆ƒ
Ú Ü ∆ƒ / Ú Ü §¤
Ú Ü §¤  0
2 2 2 2

3.5.4 Équations de quantité de mouvement

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

Figure 3.4 : Maillage décalé suivant X : le volume fini est décalé X


vers la droite du volume fini typique.

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 :

ƒ …ƒ9  Î …Î9 / Ñ …Ñ9 / Ò …Ò9 / i …i9 / 9 (3.8)


= = €¹Ù Z€k¶ = ∆‚´ ∆. ¶
= = €´Ù ¶
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 / 6ƒ9
(3.8.1)
> ƒ9 √"

On obtient les coefficients suivants :

Î 
> J [ \ ∆ƒ / > ∆´ 
/ |¤ ; Ǥ  > [ \ ∆ƒ ; |¤  > ∆´ (3.8.2)
¶ ¶ ¶ ¶
= = €¹Ù Z€´Ù = ∆‚ Ê¿ = €¹Ù Z€´Ù = ∆‚
J ¹ J J ¹

Ñ  > J [ \ ∆ƒ / > ∆´  / |¯ ; ǯ  > [ \ ∆ƒ ; |¯  > ∆´ (3.8.3)


¶ ¶ ¶ ¶
= = €´Ù 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 / 6ƒ9 . …ƒ9 §¤ ∆ƒ
a = ]¿ ∆‚´ - a J a J a
> ∆. √"
(3.8.7)

b) Suivant Y :

= +? = +(€?) = +(??) +ƒ ? √€  Z?  K + +? + +?
/ / 


5 6/ V [ \/ [ \W / (3.9)
> +. > + > +‚ +‚ " √ > + + +‚ +‚

S_` (
= ) / ` (
a )b

De la même manière que pour U on obtient l’équation algébrique suivante :

ƒ 6ƒ:  Î 6Î: / Ñ 6Ñ: / Ò 6Ò: / i 6i: / : (3.10)

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)
∆. "

. 6ƒ: ∆ƒ § / ‡ˆ. [ ¼ ´ \ ∆ƒ §


- a J a J a  Z
=
:  > 6ƒ:
a ∆´ ]‚Á
/ (ƒ
Ò )∆X
߅ƒ: / 6ƒ:
∆. √" J

Ò / ƒ
/‡ˆj Ý Þ ∆ƒ §
2
(3.10.6)

3.5.5 Équation de l’énergie

/ /  V [λ \/ [λ \W
rT rGUTH rGVTH = r rT r rT
rt rX rY P rX rX rY rY
(3.11)

Après l’intégration des différents termes sur un volume de contrôle typique, on


obtient une équation algébrique de la forme :

ƒ ƒ  Î Î / Ñ Ñ / Ò Ò / 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)

3.5.6 Équation de transfert de matière

+ = +(€) = +(?) = + + + +
/ /  V [ \/ [ \W (3.13)
+. > + > +‚ iŽ + + +; +;

Similairement à l’équation de transfert de l’énergie, on a une équation de la forme :

ƒ ƒ  Î Î / Ñ Ñ / Ò Ò / i i /  (3.14)

/ |¤ ; Ǥ  …¤a ∆ƒ ; |¤ 
Ê¿ = = ∆‚´
Î 
J > ij ]¿
(3.14.1)

Ñ  / |¯ ; ǯ  > … ā ∆ƒ ; |  ij ]´


ÊÀ = = ∆‚
J
(3.14.2)
À

Ò 
/ | ; Ç  6 a ∆ƒ ; | 
ÊÁ = = ∆´
J > ij ]‚Á
(3.14.3)

i  / |h ; Çh  > 6ha ∆ƒ ; |h  ij


Ê = = ∆´
J ]‚Â
(3.14.4)

ƒ  Î / Ñ / Ò / i /
∆´ ∆‚´
∆.
(3.14.5)

  ƒa
∆´ ∆‚´
∆.
(3.14.6)

3.5.7 Résolution des systèmes d’équations de discrétisation UPu et VPv

Il est nécessaire de connaître les valeurs de la pression aux points du maillage


typique. Donc, il nous faut une équation de discrétisation de la pression à chaque point
intérieur du maillage typique. Pour cela, on utilise l’équation de continuité discrétisée aux
nœuds du maillage typique.

¥ ¥ [+ / N‚ \ §Ð§â  0
+€ +€
(3.15)

G…¤
…¯ H∆ƒ / G6
6h H∆ƒ  0 (3.16)

37
Chapitre 3 Résolution numérique

Or, on sait déjà que suivant la direction X (équation 3.8) :

ƒ …ƒ9  Î …Î9 / Ñ …Ñ9 / Ò …Ò9 / i …i9 / 9

Qu’on peut mettre sous la forme suivante :

ƒ …ƒ9  Î …Î9 / Ñ …Ñ9 / Ò …Ò9 / i …i9 / ã9 / (ƒ


Î )∆X (3.17)

Avec

9  ã9 / (ƒ
Î )∆X (3.17.1)

. …ƒ9
a
§¤ ∆ƒ
= ]¿ ∆‚´ - J J
ã9  …ƒ9
a

߅X9
a
/ 6ƒ9
a
(3.17.2)
> ∆. √"

On remarque que la position Pu du maillage décalé est confondue avec la position e


du maillage typique donc : UPu = Ue.

ƒ …¤  Î …Î9 / Ñ …Ñ9 / Ò …Ò9 / i …i9 / ã9 / Gƒ


Î H∆X (3.18)

ä¤ / §¤ Gƒ
Î H ;
…¤  … §¤ 
∆‚´
å´

ä¯ / §¯ GÑ
ƒ H ;
…¯  … §¯ 
∆‚´
å´
(3.19)

6  6æ / § Gƒ
Ò H ; § 
∆´
å´

6h  6æh / §h Gi
ƒ H ; §h 
∆´
å´

ä¤  幀¹Ù Zå» €»ÙZå¼ €¼Ù Z彀½Ù ZçÙ


…
å´
(3.20.1)

D’une manière similaire on peut obtenir :

6æ  ¹ ¹à » »à ¼ ¼à ½ ½à à
å ? Zå ? Zå ? Zå ? Zç
å´
(3.20.2)

ä¤ et 6æ sont appelés les pseudo-vitesses suivant X et Y, respectivement.


…

On remplace les vitesses par leur expression on obtient la forme algébrique suivante :

ƒ ƒ  Î Î / Ñ Ñ / Ò Ò / i i / ƒ (3.21)

Avec :

Î  §¤ ∆ƒ ; Ñ  §¯ ∆ƒ ; Ò  § ∆ƒ ; i  §h ∆ƒ (3.21.1)

ƒ  Î /  / Ò / 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).

Ce couplage des équations des vitesses et de la pression nécessite une solution


simultanée et itérative. La solution commence par des estimations qui doivent être corrigées.

On définit :

…  …è / …f ; 6  6è / 6f ;   è / f (3.22)

• U, V et P sont les valeurs exactes.

• U*, V* et P* sont les valeurs estimées.

• …é , 6é et é sont les valeurs des corrections.

D’où on peut écrire :

ƒ …ƒ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.

D’autre part, on peut aussi écrire :

ƒ …ƒ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)

Donc, avec ces définitions si on a la correction de la pression, on calcule directement


les corrections des vitesses.

Alors on peut écrire :

39
Chapitre 3 Résolution numérique

…¤  …¤è / …¤′  …¤è / §¤ 3ƒ′


Î′ 7 (3.26.1)

…¯  … è̄ / …¯′  … è̄ / §¯ 3Ñ


ƒ′ 7 (3.26.2)

6  6 è / 6 ′  6 è / § 3ƒ′
Ò′ 7 (3.26.3)

6h  6hè / 6h′  6hè / §h 3i′


ƒ′ 7 (3.26.4)

On utilise ensuite ces expressions de Ue, Uw, Vn et Vs dans l’équation de continuité


discrétisé, on obtient alors :

ꅤè / §¤ 3ƒ′
Î′ 7ë∆ƒ
ê… è̄ / §¯ 3Ñ


ƒ′ 7ë∆ƒ / ê6 è / § 3ƒ′
Ò′ 7ë∆ƒ

ê6hè / §h 3i′
ƒ′ 7ë∆ƒ  0
(3.27)

Qui est de la forme :

ƒ ƒ′  Î Î′ / Ñ Ñ

/ Ò Ò′ / i i′ / ƒ′ (3.28)

Avec :

Î  §¤ ∆ƒ ; Ñ  §¯ ∆ƒ ; Ò  § ∆ƒ ; i  §h ∆ƒ (3.28.1)

ƒ  Î / Ñ / Ò / i (3.28.2)

ƒ′  G…¤è
… è̄ H∆ƒ / G6 è
6hè H∆ƒ  0 (3.28.3)

Donc l’équation de discrétisation de la correction de la pression est similaire à celle


de l’équation de la pression, la différence est que la source S’P contient les estimations des
vitesses alors que la source SP, contient les pseudo-vitesses.

3.6 MÉTHODE DE RÉSOLUTION

Après la discrétisation de l’équation de transport avec la méthode des volumes finis,


nous obtenons un système d’équation algébrique non linéaire. Pour résoudre ce système
d’équation, on les considère sous une forme matricielle.

_b_Φb  _b

Où :

_b est une matrice de (il-2 x jl-2) éléments.

_Φb Vecteurs des inconnues.

40
Chapitre 3 Résolution numérique

• Méthode TDMA (Algorithme de Thomas)

On résolue le système d’équations itérativement par double balayage en utilisant


l’algorithme de Thomas ou la méthode TDMA (Tri-Diagonal-Matrix-Algorithme).

Considérons le système linéaire suivant : Ð  ã

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.

En adoptant la notation suivante pour les éléments de A, x et B :

 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

• Résolution du système à matrice bidiagonale :

xn = β n

xi = βi − γ i xi −1 i = n − 1, n − 2, n − 3,......,1

41
Chapitre 3 Résolution numérique

3.7 SOLUTION DES SYSTÈMES D’ÉQUATIONS DISCRÉTISÉES

On suit l’algorithme SIMPLER (Semi-Implicit-Pressure-Linked-Equation-Revised).

1. On commence par une initialisation du champ des vitesses (Estimer un champ de


vitesses).
ä et 6
2. Avec ce champ des vitesses initiales, ou calcule les pseudo-vitesses … ä.

3. Avec les pseudo-vitesses, on peut résoudre (par la méthode itérative de balayage) le


système d’équations de discrétisation de la pression et obtenir une estimation de cette
dernière P*.
4. L’estimation de la pression est utilisée dans les équations de discrétisation des
vitesses. Alors les systèmes d’équations de discrétisation des vitesses sont résolus (par
la méthode itérative de balayage) et on obtient une estimation des vitesses U* et V*.
5. Avec les estimations des vitesses, on peut résoudre (par la méthode itérative de
balayage) le système d’équations de discrétisation de la correction de la pression et
obtenir une estimation de cette dernière P’.
6. Avec la correction de la pression, on calcule les corrections des vitesses et on corrige
le champ de vitesse.
7. On résout le système d’équations de discrétisation de la température (par la méthode
itérative de balayage) et on obtient le champ de température.
8. On résout le système d’équations de discrétisation de la concentration (par la méthode
itérative de balayage) et on obtient le champ de concentration.
9. On vérifie l’atteinte du régime permanent (stationnaire), s’il est atteint, on arrête le
calcul, sinon on augmente le temps d’un pas et avec le champ des vitesses corrigées
comme une nouvelle initialisation, on retourne à l’étape (2).

3.8 CRITÈRE DE CONVERGENCE

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)

3.9 DÉTAILS NUMÉRIQUES

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

4.1 EFFET DE MAILLAGE

Afin d’illustrer l’influence du maillage, sur l’écoulement et le transfert de chaleur et


de masse, représentés respectivement par les nombres moyens de Nusselt et de Sherwood aux
parois actives, nous avons exécuté le code pour différents maillages uniformes, dans le cas
coopérant, avec : A=4, Pr=0.71, Sc=0.65, Grt= 105, Grc=105 Cf=0.55, Da=10-2, φ=0.7.

Tableau 4.1: Effet du maillage sur les caractéristiques de l’écoulement.

N Maillage òó ôõ |ö÷øù| |ú÷øù|

1 42x22 0.3892608 2.735322 32.655810 44.032120

2 82x44 0.3834970 2.604439 32.434700 43.058940

3 100x50 0.3825613 2.573911 32.367570 43.740600

4 122x66 0.3816253 2.547732 32.289580 42.996600

5 132x82 0.3354265 1.731527 28.961690 20.737450

6 164x88 0.3349295 1.721992 28.882630 20.972090

7 192x92 0.3348125 1.716262 28.809710 21.011330

8 262x122 0.3339869 1.704769 28.643840 21.162850

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

Tableau 4.2 : Effet de l’espacement sur les caractéristiques de l’écoulement.

N Espacement Nu Sh |Umax| |Vmax|


1 uniforme 0.3349295 1.721992 28.882630 20.972090
2 C=1.01 0.344696 1.709807 28.824840 21.073890
3 C=1.03 0.337537 1.693545 28.737910 21.180290

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.

La figure 4.1 représente les différents profils de la vitesse horizontale au milieu de


l’enceinte (X=2) pour différents maillages. On remarque qu’à l’exception du premier maillage
de 122×66, tous les autres maillages ont des profils de vitesse identiques.

1.00
Maillage 122x66
Maillage 132x82
Maillage 164x88

0.80 Maillage 192x92


Maillage 262x122
Maillage 164x88, C=1.03

0.60

0.40

0.20

0.00

-40.00 -20.00 0.00 20.00 40.00

Figure 4.1 : Profils de la vitesse horizontale U à X=2 pour différents maillages.

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.

4.2 VALIDATION DU CODE DE CALCUL

Pour valider notre code de calcul et vérifier l'exactitude numérique de la présente


étude, nous avons procédé à diverses comparaisons avec différents résultats disponibles dans
la littérature et repris par Younsi et al. (2002). Les tableaux 4.3 et 4.4 donnent une
comparaison du transfert de chaleur obtenu par notre code et les résultats de différents
travaux, dans le cas de la convection thermique en cavité carrée dans un milieu poreux. Les
modèles de Darcy et de Darcy-Brinkman on été utilisés, respectivement. Cette validation a été
exécutée employant un maillage uniforme 81×81, avec les paramètres suivants : Ra= 3.106,
Pr=1, Le=100, A=1, φ=0.7, Cf=0.25.

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.

Concernant la convection double diffusive dans un milieu poreux, une comparaison


de la structure d’écoulement et des distributions de chaleur et de matière a été effectuée avec
cette référence. Qualitativement, la figure 4.3 montre une concordance totale entre nos
résultats et ceux de Younsi et al. (2002), obtenus avec un maillage 81×81 et les paramètres
suivant : Ra= 3.106, Da=10-4, Pr=0.149, Le=10, A=1, φ=0.7, Cf=0.25.

Tableau 4.3 : Modèle de Darcy (transfert de chaleur pure : N=0)

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

Tableau 4.4 : Modèle de Darcy-Brinkman (transfert de chaleur pure : N=0).

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

4.3 RÉSULTATS ET DISCUSSIONS

Les paramètres qui contrôlent l’écoulement du fluide et les transferts de chaleur et de


masse en convection naturelle bidiffusive sont : A, Grt, Grc, Pr, Cf, Sc, Da, φ. Tous ces
paramètres sont maintenues fixes : A=4, Pr=0.71, Sc=0.65, Cf=0.55, Da=10-2 et 5.10-7, φ=0.7,
à l’exception des nombres de Grashof thermique et solutal qui seront variés de 104 à 108 et de
105 à 106, respectivement. Pour la première et la deuxième partie, on a utilisé le modèle de
Darcy-Brinkman-Forchheimer. La troisième partie est consacrée à la comparaison entre les
modèles de Darcy et Darcy-Brinkman-Forchheimer.

4.3.1 Cas coopérant : milieu fluide avec Grc=105

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.

La figure 4.3 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,
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.3a). Un courant rotatif monocellulaire dans le sens horaire est obtenu dans
l’enceinte (fig. 4.3b). La variation de la température dans l’enceinte est présentée par des
isothermes verticales déformées de gauche à droite en haut de l’enceinte et de droite à gauche
en bas de l’enceinte (fig. 4.3c). Tandis que, la variation de la concentration 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.3d). On conclut donc que, dans
ces conditions, la force de buoyancy solutale conduit, particulièrement, la circulation dans le
milieu fluide.

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

Figure 4.3 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=104, Grc=105

La figure 4.4 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=105 et Grc=105, ce
qui correspond à N=1. On remarque que les plus grandes valeurs de la pression se situent en
haut du coté de la paroi verticale gauche (où la température est maximale) et en bas du coté de
la paroi verticale droite (où la concentration est maximale), tandis que dans le reste de
l’enceinte la pression a une valeur quasi-constante selon une direction oblique. Ceci est dû
aux conditions aux limites thermique et solutale et aux valeurs prises par les nombres de
Grashof thermique et solutal (fig. 4.4a). Similairement au cas précédent, un courant rotatif
monocellulaire dans le sens horaire est obtenu dans l’enceinte (fig. 4.4b). Mais les Lignes de
courant sont plus confinées, ce qui signifie que l’écoulement obtenu est très important. Les
valeurs numériques des vitesses obtenues confirment cette analyse (tableau 4.5). Aussi, la
variation de la température dans l’enceinte est présentée par des isothermes verticales
largement déformés de gauche à droite en haut de l’enceinte et de droite à gauche en bas de
l’enceinte (fig. 4.4c). Tandis que, la variation de la concentration dans l’enceinte est présentée
par des isoconcentrations fortement concentrées près des parois verticales isothermes. Tandis
qu’au centre de l’enceinte la concentration a une valeur constante, C≈0.5 (fig. 4.4d). On
conclut donc que dans ces conditions, l’écoulement dans le milieu fluide est conduit
conjointement par les forces de buoyancy thermique et solutale.

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

Figure 4.4 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=105, Grc=105

La figure 4.5 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=106 et
Grc=105. Dans ce cas, on remarque que les plus grandes valeurs de la pression se situent en
haut de l’enceinte du coté de la paroi verticale gauche (où la température 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.5a). L’écoulement du fluide a une structure
monocellulaire tournant dans le sens horaire (fig. 4.5b). Les Lignes de courant sont plus
confinées, ce qui signifie que l’écoulement obtenu est encore plus important que le cas
précédent (Grt=105). Les valeurs numériques des vitesses obtenues confirment cette analyse
(tableau 4.5). Les isothermes sont fortement concentrées près des parois verticales, loin des
parois isothermes et des parois horizontales (vers le milieu de l’enceinte), 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.5c). 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.5d). On conclut donc que dans ces conditions, l’écoulement dans le milieu fluide est
conduit particulièrement par la force de buoyancy thermique.

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

Figure 4.5 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=106, Grc=105

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

Figure 4.6 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=107, Grc=105.

a- Contour de pression

b- Lignes de courant

57
Chapitre 4 Résultats et discussions

c- Isothermes

d- Isoconcentrations

Figure 4.7 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=108, Grc=105.

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

Tableau 4.5 : Caractéristiques de l’écoulement, milieu fluide, cas coopérant : Grc=105.

Nu Sh |Umax| |Vmax|

Grt=104 0.5944415 3.073611 67.396520 83.958370

Grt=105 1.213076 4.036448 136.717500 133.954700

Grc=105 Grt=106 3.092439 7.161746 425.380700 429.689600

Grt=107 6.789007 14.045480 1273.842000 1464.307000

Grt=108 13.662350 27.708770 3493.656000 4522.051000

4.3.2 Modèle de Darcy-Brinkman-Forchheimer

4.3.2.1 Première partie : cas coopérant

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

4.3.2.1.1.1 Structure d’écoulement et distributions de la pression, de la température et de


la concentration

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

Figure 4.8 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=104, Grc=105.

La figure 4.9 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=105 et
Grc=105, ce qui correspond à N=1. On remarque que les plus grandes valeurs de la pression se
situent en haut du coté de la paroi verticale gauche (où la température est maximale) et 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.9a). Similairement au cas précédent, un courant rotatif
monocellulaire dans le sens horaire est obtenu dans l’enceinte (fig. 4.9b). Mais les Lignes de
courant sont plus confinées, ce qui signifie que l’écoulement obtenu est très important. Les
valeurs numériques des vitesses obtenues confirment cette analyse (tableau 4.6). Aussi, 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.9c). Tandis que, la variation de la concentration dans l’enceinte est présentée
par des isoconcentrations verticales fortement déformés de droite à gauche en bas de
l’enceinte et de gauche à droite en haut de l’enceinte (fig. 4.9d). On conclut donc que dans ces
conditions, l’écoulement du fluide dans le milieu poreux est conduit conjointement par les
forces de buoyancy thermique et solutale.

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

Figure 4.9 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=105, Grc=105.

Les figures 4.10 et 4.11 illustrent 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=106, 107 et Grc=105. Dans ce cas, on remarque que les plus grandes
valeurs de la pression se situent en haut de l’enceinte du coté de la paroi verticale gauche (où
la température 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.10a et 4.11a).
L’écoulement du fluide a une structure monocellulaire tournant dans le sens horaire (fig.
4.10b et 4.11b). Les Lignes de courant sont plus confinées, ce qui signifie que l’écoulement
obtenu est plus important. Les valeurs numériques des vitesses obtenues confirment cette
analyse (tableau 4.6). Les isothermes sont largement déformées de gauche à droite en haut de
l’enceinte et de droite à gauche en bas de l’enceinte (fig. 4.10c) dans le cas de Grt=106. Avec
Grt=107, loin des parois isothermes et des parois horizontales (vers le milieu de l’enceinte), 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.11c). 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.10d et 4.11d). On conclut donc que dans ces conditions, l’écoulement
du fluide dans le milieu poreux est conduit particulièrement par la force de buoyancy
thermique.

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

Figure 4.10 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=106, Grc=105.

a- Contour de pression

b- Lignes de courant

65
Chapitre 4 Résultats et discussions

c- Isothermes

d- Isoconcentrations

Figure 4.11 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=107, Grc=105.

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, Grc=105 sont présentés sur la
figure 4.12. Dans ce cas, on remarque 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.12a). L’écoulement du fluide a une structure monocellulaire tournant dans le
sens horaire (fig. 4.12b). Les Lignes de courant sont plus confinées, ce qui signifie une
accélération de l’écoulement. Les valeurs numériques des vitesses obtenues confirment cette
analyse (tableau 4.6). 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.12c). 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

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

Figure 4.12 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=108, Grc=105.

4.3.2.1.1.2 Variation temporelle des nombres de Nusselt et de Sherwood

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

Figure 4.13 : Variation temporelle du nombre de Nusselt.

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

Figure 4.14 : Variation temporelle du Nombre de Sherwood.

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.

Tableau 4.6 : Caractéristiques de l’écoulement, milieu poreux : cas coopérant, Grc=105.

Nu Sh |Umax| |Vmax|

Grt=104 0.2809787 1.052308 16.978080 12.964030

Grt=105 0.3337525 1.693558 28.737910 21.180290

Grc=105 Grt=106 0.9193805 4.279505 86.903780 93.295250

Grt=107 1.967479 6.544524 191.009800 178.738600

Grt=108 6.843239 18.810100 508.17500 1341.123000

4.3.2.1.2 Cas de Grc=5.105

4.3.2.1.2.1 Structure d’écoulement et distributions de la pression, de la température et de


la concentration

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.15 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 et Grc=5.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.15a). Un écoulement rotatif monocellulaire dans le sens horaire est obtenu dans
l’enceinte (fig. 4.15b). 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.15c). Tandis que, la variation de la concentration

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

Figure 4.15 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=104, Grc=5.105.

La répartition de la pression, les contours de lignes de courant, de la température et


de la concentration obtenues, respectivement, pour Grt=105 et Grc=5.105 sont illustrés par la
figure 4.16. Similairement au cas précédent, 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.16a). Un écoulement rotatif
monocellulaire dans le sens horaire est obtenu dans l’enceinte (fig. 4.16b). 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.16c).
Tandis que, la variation de la concentration 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

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.

Les figures 4.17 et 4.18 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=106 et 107, Grc=5.105. Dans ce cas, on remarque aussi que les plus grandes valeurs de la
pression se situent en haut de l’enceinte du coté de la paroi verticale gauche (où la
température 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.17a et 4.18a).
L’écoulement du fluide a une structure monocellulaire tournant dans le sens horaire (fig.
4.17b et 4.18b). Les Lignes de courant sont plus confinées, ce qui signifie que l’écoulement
obtenu est plus important. Les valeurs numériques des vitesses obtenues confirment cette
analyse (tableau 4.7). Les isothermes sont largement déformées de gauche à droite en haut de
l’enceinte et de droite à gauche en bas de l’enceinte (fig. 4.17c) dans le cas de Grt = 106. Avec
Grt = 107, loin des parois isothermes et des parois horizontales (vers le milieu de l’enceinte),
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.18c). 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.17d et 4.18d). On conclut donc que dans ces conditions,
l’écoulement du fluide dans le milieu poreux est conduit particulièrement par la force de
buoyancy thermique.

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

Figure 4.17 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=106, Grc=5.105.

a- Contour de pression

b- Lignes de courant

76
Chapitre 4 Résultats et discussions

c- Isothermes

d- Isoconcentrations

Figure 4.18 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=107, Grc=5.105.

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, Grc=5.105 sont présentés sur la
figure 4.19. Similairement au cas de Grc=105, on remarque 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.19a). L’écoulement du fluide a une structure
monocellulaire tournant dans le sens horaire (fig. 4.19b). Les Lignes de courant sont plus
confinées, ce qui signifie une accélération de l’écoulement. Les valeurs numériques des
vitesses obtenues confirment cette analyse (tableau 4.7). 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.19c). 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

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

Figure 4.19 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=108, Grc=5.105.

4.3.2.1.2.2 La variation temporelle des nombres de Nusselt et de Sherwood

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

Figure 4.20 : Variation temporelle du nombre de Nusselt.

79
Chapitre 4 Résultats et discussions

La figure 4.21 illustre La variation du nombre de Sherwood avec le temps. On


remarque 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.

4
Grt=10
5
Grt=10
25 6
Grt=10
7
Grt=10
Sh

8
Grt=10
20

15

10

0 10000 20000 30000 40000 50000


t

Figure 4.21 : Variation temporelle du nombre de Sherwood.

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

Tableau 4.7 : Caractéristiques de l’écoulement, milieu poreux : cas coopérant, Grc=5.105.

Nu Sh |Umax| |Vmax|

Grt=104 0.3869724 2.737755 35.430850 53.368250

Grt=105 0.4326397 3.087615 41.145760 59.657010

Grc=5.105 Grt=106 0.9661006 4.982965 88.110660 119.430100

Grt=107 1.984040 6.654623 191.373200 183.929500

Grt=108 6.848080 18.835520 508.224000 1342.97800

4.3.2.1.3 Cas de Grc=106

4.3.2.1.3.1 Structure d’écoulement et distributions de la pression, de la température et de


la concentration

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

Figure 4.22 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=104, Grc=106.

a- Contour de pression

b- Lignes de courant

83
Chapitre 4 Résultats et discussions

c- Isothermes

d- Isoconcentrations

Figure 4.23 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=105, Grc=106.

La Figure 4.24 représente la répartition de la pression, les contours de lignes de


courant, de la température et de la concentration obtenues, respectivement, pour Grt=106 et
Grc=106, ce qui correspond à N=1. On remarque que les plus grandes valeurs de la pression se
situent en haut du coté de la paroi verticale gauche (où la température est maximale) et en bas
du coté de la paroi verticale droite (où la concentration est maximale) (fig. 4.24a).
L’écoulement obtenu monocellulaire dans le sens horaire (fig. 4.24b). Aussi, la variation de la
température dans l’enceinte est présentée par des isothermes verticales déformées de gauche à
droite en haut de l’enceinte et de droite à gauche en bas de l’enceinte (fig. 4.24c). Tandis que,
la variation de la concentration dans l’enceinte est présentée par des gradients fortement
concentrés près des parois verticales. La concentration au centre de l’enceinte est
pratiquement constante C≈0.5 (fig. 4.24d). On conclut donc que dans ces conditions,
l’écoulement du fluide dans le milieu poreux est conduit conjointement par les forces de
buoyancy thermique et solutale.

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

Figure 4.24 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=106, Grc=106.

La répartition de la pression, les contours de lignes de courant, de la température et


de la concentration obtenues, respectivement, pour Grt=107 et Grc=106 sont représentés par la
figure 4.25. Dans ce cas, on remarque que les plus grandes valeurs de la pression se situent en
haut de l’enceinte du coté de la paroi verticale gauche (fig. 4.25a). L’écoulement du fluide a
une structure monocellulaire tournant dans le sens horaire (fig. 4.25b). Les Lignes de courant
sont plus confinées, ce qui signifie que l’écoulement obtenu est plus important. Les valeurs
numériques des vitesses obtenues confirment cette analyse (tableau 4.8). Loin des parois
isothermes et des parois horizontales (vers le milieu de l’enceinte), 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.25c). 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.25d). On conclut donc que dans ces conditions, l’écoulement du fluide dans le milieu
poreux est conduit particulièrement par la force de buoyancy thermique.

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

Figure 4.25 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=107, Grc=106.

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=106 sont illustrés sur la
figure 4.26. Dans ce cas, on remarque 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.26a). L’écoulement du fluide a une structure monocellulaire tournant dans le
sens horaire (fig. 4.26b). Les Lignes de courant sont plus confinées, ce qui signifie une
accélération de l’écoulement. Les valeurs numériques des vitesses obtenues confirment cette
analyse (tableau 4.8). On obtient une stratification thermique parfaitement stable, la
distribution de la température et faiblement variable suivant X mais variable suivant Y (fig.
4.26c). 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.26d). On conclut donc que dans ces conditions, l’écoulement du
fluide dans le milieu poreux est totalement conduit par la force de buoyancy thermique.

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

Figure 4.26 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=108, Grc=106.

4.3.2.1.3.2 La variation temporelle des nombres de Nusselt et de Sherwood

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

Figure 4.27 : Variation temporelle du nombre de Nusselt.

90
Chapitre 4 Résultats et discussions

La figure 4.28 représente 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 pour les cinq cas étudiés. Les valeurs de Sherwood sont augmentent les valeurs de
Grashof thermique élevée.

4
Grt=10
5
Grt=10
25 6
Grt=10
7
Grt=10
Sh

8
Grt=10
20

15

10

0 10000 20000 30000 40000 50000


t

Figure 4.28 : Variation temporelle du nombre de Sherwood.

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

Tableau 4.8 : Caractéristiques de l’écoulement, milieu poreux : cas coopérant, Grc=106.

Nu Sh |Umax| |Vmax|

Grt=104 0.4843733 3.936440 46.708630 92.622730

Grt=105 0.5740218 4.177529 50.945190 97.951780

Grc=106 Grt=106 1.104095 5.715921 88.156890 150.260600

Grt=107 3.640586 11.117050 250.842800 558.130700

Grt=108 6.854128 18.867330 508.279700 1345.294000

4.3.2.1.4 Cas de Da=5.10-7 et Grc=105

4.3.2.1.4.1 Structure d’écoulement et distributions de la pression, de la température et de


la concentration

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

Figure 4.29: Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=104, Grc=105.

a- Contour de pression

b- Lignes de courant

94
Chapitre 4 Résultats et discussions

c- Isothermes

d- Isoconcentrations

Figure 4.30 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=105, Grc=105.

a- Contour de pression

95
Chapitre 4 Résultats et discussions

b- Lignes de courant

c- Isothermes

d- Isoconcentrations

Figure 4.31 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=106, Grc=105.

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

Figure 4.33 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=108, Grc=105.

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

Tableau 4.9 : Caractéristiques de l’écoulement, milieu poreux : cas coopérant, Da=5.10-7 et


Grc=106.

Nu Sh |Umax| |Vmax|

Grt=104 0.2500296 0.2505172 6.710425.10-3 4.968148.10-3

Grt=105 0.2500301 0.2505074 1.220169.10-2 9.036320.10-3


Da=5.10-7
5
Grt=106 0.2500380 0.2505665 6.711295.10-2 4.971865.10-2
Grc=10
Grt=107 0.2500812 0.2520148 0.6160218 0.2520148

Grt=108 0.2540514 0.3857225 6.058145 4.567605

4.3.2.1.4.2 La variation temporelle des nombres de Nusselt et de Sherwood

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

Figure 4.34 : Variation temporelle du nombre de Nusselt.

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

0 10000 20000 30000 40000 50000


t

Figure 4.35 : Variation temporelle du nombre de Sherwood.

101
Chapitre 4 Résultats et discussions

4.3.2.2 Deuxième partie : cas opposant, Da=10-2

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.

4.3.2.2.1 Cas de Grc=105

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.

La figure 4.36 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 gauche, 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.36a). L’écoulement
obtenu dans l’enceinte est monocellulaire et a pour sens de rotation le sens antihoraire (fig.
4.36b). 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.36c). 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.36d). On conclut donc que,
dans ces conditions, la force de buoyancy solutale conduit, en particulier, l’écoulement dans
le milieu poreux.

a- Contour de pression

102
Chapitre 4 Résultats et discussions

b- Lignes de courant

c- Isothermes

d- Isoconcentrations

Figure 4.36 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=104, Grc=105.

103
Chapitre 4 Résultats et discussions

La répartition de la pression, les contours de lignes de courant, de la température et


de la concentration obtenues, respectivement, pour Grt=105 et Grc=105 (N= –1), sont illustrés
par la figure 4.37. En raison de la compétitivité entre les forces de volume thermique et
solutale de valeur égale et de sens opposé, on constate que la différence entre les valeurs
maximale et minimale de la pression est relativement négligeable (fig. 4.37a). On obtient un
écoulement bicellulaire contrarotatif (fig. 4.37b). Les variations de la température et de la
concentration dans l’enceinte sont présentées par des isothermes et des isoconcentrations
verticales quasi-parallèles (fig. 4.37c et 4.37d). On conclut donc que, dans ces conditions, les
transferts thermique et massique sont purement diffusifs.

a- Contour de pression

b- Lignes de courant

104
Chapitre 4 Résultats et discussions

c- Isothermes

d- Isoconcentrations

Figure 4.37 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=105, Grc=105.

La figure 4.38 représente la répartition de la pression, la structure de l’écoulement et


les distributions de la température et de la concentration obtenues, respectivement, pour
Grt=106 et Grc=105. On constate que les plus grandes valeurs de la pression se situent en haut
du coté de la paroi verticale gauche, 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.38a). On
obtient un écoulement monocellulaire tournant dans le sens horaire (fig. 4.38b). 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.38c). 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 haut de l’enceinte et de
droite à gauche en bas de l’enceinte. Les gradients de concentration sont concentrés près des
parois verticales (fig. 4.38d). On conclut donc que, dans ces conditions, la force de buoyancy
thermique conduit, en particulier, l’écoulement dans le milieu poreux.

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

Figure 4.38 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=106, Grc=105.

4.3.2.2.2 Cas de Grc=5.105

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

Figure 4.39 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=104, Grc=5.105.

a- Contour de pression

b- Lignes de courant

109
Chapitre 4 Résultats et discussions

c- Isothermes

d- Isoconcentrations

Figure 4.40 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=105, Grc=5.105.

La figure 4.41 représente la répartition de la pression, la structure de l’écoulement et


les distributions de la température et de la concentration obtenues, respectivement, pour
Grt=106 et Grc=5.105. On constate dans ce cas que les plus grandes valeurs de la pression se
situent en haut du coté de la paroi verticale gauche (fig. 4.41a). L’écoulement obtenu est
monocellulaire et son sens de rotation devient le sens horaire (fig. 4.41b). La variation de la
température dans l’enceinte est présentée par des isothermes verticales déformées de gauche à
droite en haut de l’enceinte et de droite à gauche en bas de l’enceinte (fig. 4.41c). Alors 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 haut de l’enceinte et de droite à gauche
en bas de l’enceinte et très concentrés près des parois verticales. Tandis que la concentration a
une valeur pratiquement constante dans le centre de l’enceinte C≈0.5 (fig. 4.41d). On conclut
donc que, dans ces conditions, la force de buoyancy thermique conduit, en particulier,
l’écoulement dans le milieu poreux.

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

Figure 4.41 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=106, Grc=5.105.

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

Tableau 4.10 : Caractéristiques de l’écoulement, milieu poreux : cas opposant, Da=10-2,


Grc=5.105 et 5.105.

Nu Sh |Umax| |Vmax|

Grt=104 0.2726482 0.8850822 14.176030 10.955580

Grc=105 Grt=105 0.2499751 0.2503316 0.01775341 0.005912934

Grt=106 0.9618942 3.901042 84.147990 80.958320

Grt=104 0.4025944 2.652867 33.667740 51.943340

Grc=5.105 Grt=105 0.3578133 2.239670 28.453250 45.312160

Grt=106 0.8365300 2.830520 81.048230 56.048230

4.3.3 Modèle de Darcy, cas coopérant

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

4.3.3.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.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

Figure 4.42 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=104, Grc=105.

La figure 4.43 représente la répartition de la pression, les contours de lignes de


courant, de la température et de la concentration obtenues, respectivement, pour Grt=105 et
Grc=105. Ce cas est aussi similaire au cas du modèle de Darcy-Brinkman-Forchheimer déjà
traité (§ 4.3.2.1.1.1 page 62). Mais, on remarque que les valeurs maximales de la pression
sont relativement faibles. Cependant, les taux de transfert de chaleur et de masse sont plus
élevées (tableau 4.11).

a- Contour de pression

115
Chapitre 4 Résultats et discussions

b- Lignes de courant

c- Isothermes

d- Isoconcentrations

Figure 4.43 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=105, Grc =105.

116
Chapitre 4 Résultats et discussions

La figure 4.44 représente la répartition de la pression, les contours de lignes de


courant, de la température et de la concentration obtenues, respectivement, pour Grt=106,
Grc=105. Ce cas est aussi similaire au cas du modèle de Darcy-Brinkman-Forchheimer déjà
cité (§ 4.3.2.1.1.1 page 64). Mais, on remarque que les taux de transfert de chaleur et de
masse sont encore plus élevées (tableau 4.11).

a- Contour de pression

b- Lignes de courant

117
Chapitre 4 Résultats et discussions

c- Isothermes

d- Isoconcentrations

Figure 4.44 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=106, Grc =105.

Les figures 4.45 et 4.46 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, Grc=105. Ces deux cas présentent certaines différences avec les cas du modèle de
Darcy-Brinkman-Forchheimer (§ 4.3.2.1.1.1 pages 64 et 67). On remarque que la pression
varie suivant Y, les plus grandes valeurs 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.45a et
4.46a). 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
tournants dans le sens horaire (fig. 4.45b et 4.46b). 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.45c et 4.46c). Alors que,
les gradients de concentration sont très fortement concentrés près des parois verticales

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

Figure 4.45 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=107, Grc=105.

a- Contour de pression

120
Chapitre 4 Résultats et discussions

b- Lignes de courant

c- Isothermes

d- Isoconcentrations

Figure 4.46 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=108, Grc=105.

121
Chapitre 4 Résultats et discussions

4.3.3.1.2 Cas de Grc=5.105

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

Figure 4.47 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=104, Grc=5.105.

a- Contour de pression

123
Chapitre 4 Résultats et discussions

b- Lignes de courant

c- Isothermes

d- Isoconcentrations

Figure 4.48 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=105, Grc=5.105.

La figure 4.49 représente la répartition de la pression, les contours de lignes de


courant, de la température et de la concentration obtenues, respectivement, pour Grt=106,

124
Chapitre 4 Résultats et discussions

Grc=5.105. De même, ce cas est aussi similaire au cas du modèle de Darcy-Brinkman-


Forchheimer (§ 4.3.2.1.2.1 page 75). Mais, on remarque que les valeurs minimales de la
pression sont relativement faibles. Cependant, les taux de transfert de chaleur et de masse sont
plus élevées (tableau 4.11).

a- Contour de pression

b- Lignes de courant

125
Chapitre 4 Résultats et discussions

c- Isothermes

d- Isoconcentrations

Figure 4.49 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=106, Grc=5.105.

Les figures 4.50 et 4.51 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. Ces deux cas présentent aussi certaines différences avec les cas du
modèle de Darcy-Brinkman-Forchheimer (§ 4.3.2.1.2.1 pages 75 et 78). On remarque que la
pression varie suivant Y, les plus grandes valeurs 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.50a et
4.51a). 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 tournants dans
le sens horaire (fig. 4.50b et 4.51b). 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.50c et 4.51.c). Alors que, les gradients de
concentration sont très fortement concentrés près des parois verticales isothermes. Tandis

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

Figure 4.50 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=107, Grc=5.105.

a- Contour de pression

b- Lignes de courant

128
Chapitre 4 Résultats et discussions

c- Isothermes

d- Isoconcentrations

Figure 4.51 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=108, Grc=5.105.

4.3.3.2 Cas de Da=5.10-7 et Grc=105

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

Figure 4.52 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=104, Grc=105.

a- Contour de pression

131
Chapitre 4 Résultats et discussions

b- Lignes de courant

c- Isothermes

d- Isoconcentrations

Figure 4.53 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=105, Grc=105.

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

Figure 4.54 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=106, Grc=105.

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

Figure 4.56 : Contour de la pression, isothermes, isoconcentrations et structure de


l’écoulement pour : Grt=108, Grc=105.

136
Chapitre 4 Résultats et discussions

4.3.4 Comparaisons des deux modèles Darcy-Brinkman-Forchheimer et Darcy

La figure 4.57 illustre La variation de nombre de Nusselt avec le Grashof thermique


pour les deux Modèles avec Grc=105. Il est montré que pour des petites valeurs de Grashof
thermique la différence entre les nombres de Nusselt pour les deux Modèles est petite et elle
augmente avec l’augmentation des valeurs de Grashof thermique.

Modèle de Darcy-Brinkman-Forchheimer
Modèle de Darcy

10
Nu

0,00E+000 2,00E+007 4,00E+007 6,00E+007 8,00E+007 1,00E+008


Grt

Figure 4.57 : Nombre de Nusselt en fonction en fonction de Grt pour Grc=105.

La figure 4.58 illustre la variation de nombre de Sherwood avec le Grashof


thermique pour les deux Modèles avec Grc=105. Il est montré que pour des petites valeurs de
Grashof thermique la différence entre les nombres de Sherwood pour les deux Modèles est
petite et elle augmente avec l’augmentation des valeurs de Grashof thermique.

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

0,00E+000 2,00E+007 4,00E+007 6,00E+007 8,00E+007 1,00E+008


Grt

Figure 4.58 : Nombre de Sherwood en fonction en fonction de Grt pour Grc=105.

La figure 4.59 illustre la variation de nombre de Nusselt avec le Grashof thermique


pour les deux Modèles et avec un Grashof solutal égale à 5.105. Il est montré que pour des
petites valeurs de Grashof thermique la différence entre les nombres de Nusselt pour les deux
Modèles est petite et elle augmente avec l’augmentation des valeurs de Grashof thermique.

La figure 4.60 illustre la variation de nombre de Sherwood avec le Grashof


thermique pour les deux Modèles et avec un Grashof solutale égale a 5.105. Il est montré que
pour des petites valeurs de Grashof thermique la différence entre les nombres de Sherwood
pour les deux Modèles est petite et elle augmente avec l’augmentation des valeurs de Grt.

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.

On conclut donc que le modèle de Darcy-Brinkman-Forchheimer peut être


approximé par le modèle de Darcy dans le cas de φ<0.8, Re<1 et Da<10-6, Re est le nombre
de Reynolds. Ce résultat est conforme avec les travaux de Nield et Bejan (2006) et Vasseur et
al. (1989). Ces observations sont aussi illustrées sur le tableau 4.11.

138
Chapitre 4 Résultats et discussions

Modèle de Darcy-Brinkman-Forchheimer
Modèle de Darcy

10

9
Nu
8

0,00E+000 2,00E+007 4,00E+007 6,00E+007 8,00E+007 1,00E+008

Grt

Figure 4.59 : Nombre de Nusselt en fonction en fonction de Grt pour Grc=5.105.

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

0,00E+000 2,00E+007 4,00E+007 6,00E+007 8,00E+007 1,00E+008

Grt

Figure 4.60 : Nombre de Sherwood en fonction du Grashof thermique pour Grc=5.105.

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

0,00E+000 2,00E+007 4,00E+007 6,00E+007 8,00E+007 1,00E+008

Grt

Figure 4.61 : Nombre de Nusselt en fonction de Grt pour Grc=105 et Da=5.10-7.

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

0,00E+000 2,00E+007 4,00E+007 6,00E+007 8,00E+007 1,00E+008


Grt

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

Tableau 4.11 : Caractéristiques de l’écoulement, modèle de Darcy, cas coopérant.

Grt Nu Sh |Umax| |Vmax|

Grt=104 0.3347596 1.635408 37.680940 22.096320

Grt=105 0.5038062 2.703228 71.484860 39.043220

Grc=105 Grt=106 2.101811 7.427006 251.450000 244.297100

Grt=107 2.639083 7.983262 573.168000 295.656000

Grt=108 9.532064 26.081860 1708.551000 2498.12900


-2
Da=10
Grt=104 0.5990508 4.430360 91.851190 119.560800

Grt=105 0.7290959 5.042470 112.827500 135.026500

Grc=5.105 Grt=106 2.218518 8.126999 272.241100 315.60700

Grt=107 2.664422 8.617319 580.238700 305.963400

Grt=108 9.539660 26.122680 1710.352000 2503.65500

Grt=104 0.2501644 0.2501398 0.006773777 0.00506373

Grt=105 0.2500308 0.2504906 0.01231496 0.00918984

Da=5.10-7 Grc=105 Grt=106 0.2500398 0.2505076 0.06773916 0.00505695

Grt=107 0.2500786 0.2520482 0.6219234 0.4644038

Grt=108 0.2540710 0.3862819 6.128740 4.651964

141
Conclusion générale

CONCLUSION GÉNÉRALE

Une étude de la convection naturelle bidiffusive coopérante et opposante dans une


enceinte rectangulaire complètement rempli par un milieu poreux saturé par un fluide binaire
a été entreprise. Le fluide est l’air et la substance diffusée dans l’air de l’enceinte est la vapeur
d’eau. Les parois verticales de l’enceinte sont soumises à des gradients de températures et de
concentrations. Les parois horizontales sont adiabatiques et imperméables. Le milieu poreux
est supposé homogène et isotrope. L’écoulement dans l’enceinte est modélisé par le modèle
de Darcy-Brinkman-Forchheimer. L’écoulement et le transfert thermique et massique sont
modélisés par les équations de continuité, de quantité de mouvement, d’énergie et de masse.
Ces équations ont été discrétisées avec la méthode des volumes et résolues par la méthode
itérative de double balayage selon l’algorithme de Thomas (TDMA).

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.

Avec le modèle de Darcy-Brinkman-Forchheimer, quand les forces de volumes sont


coopérantes, pour des valeurs de Grt inférieur à Grc, l’écoulement obtenu est monocellulaire, a
pour sens de rotation le sens horaire et il est pratiquement conduit par la force de buoyancy
solutale. Tandis que, pour des valeurs de Grt supérieur à Grc, la force de buoyancy thermique
domine et conduit l’écoulement. Lorsque les forces de volumes sont opposantes, l’écoulement
obtenu est monocellulaire et a pour sens de rotation le sens antihoraire et il est conduit par la
force de volume solutale. Pour Grt égale à Grc, un faible écoulement bicellulaire contrarotatif
apparait et le transfert de chaleur et de masse devient diffusif, à cause de la compétition entre
les deux forces de même valeur et de sens opposés. Avec des valeurs de Grt supérieur à Grc, la
force de buoyancy thermique domine et on obtient un écoulement monocellulaire tournant

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.

Une comparaison du modèle de Darcy avec le modèle de Darcy-Brinkman-


Forchheimer dans le cas coopérant a été faite en employant différentes valeurs de Grt et Grc,
pour deux nombres de Darcy 10-2 et 5.10-7. On a constaté que pour Da=10-2, il y a des
différences notables entre les deux modèles de point de vue structure d’écoulement et transfert
de chaleur et de masse. Dans le cas de Da=5.10-7, les résultats obtenues avec les deux modèles
sont parfaitement identiques. On conclu que le modèle de Darcy-Brinkman-Forchheimer ne
peut être approximé par le modèle de Darcy que sous certaines conditions : φ≤0.8, Re<1 et
Da<10-6.

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.

2. AKROUR D., BENNACER R., BETROUNI M., Effet d’hystérésis en convection


thermosolutal avec des gradients croisés (Etang de stockage). Rev. Energ. Ren., vol. 3 pp.
39-47, 2000.

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.

4. BAHLOUL A., YAHIAOUI M. A., VASSEUR P., ROBILLARD L., Thermogravitional


separation in a vertical annular porous layer. Int. Comm. Heat Mass Transfer, vol. 31, N°
6, pp 783-794, 2004.

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.

17. DANQUIGNY C., Etude expérimentale du transfert de masse en milieu poreux


hétérogène. Thèses de doctorat, université de Louis Pasteur, Strasbourg, 2003.

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

30. MARCLO J. S. L., LUZIA A. T., Modeling of double-diffusive turbulent natural


convection in porous media. Int. J. Heat Mass Transfer, vol. 47, pp. 4233-4241, 2007.

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.

39. SIVASANKARAN S., KANDASWAMY P., NG C. O., Double-diffusive convection of


anomalous density fluids in a porous cavity. Transp. Porous Med., vol. 71, pp. 133-145,
2008.

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.

41. VERSTEEG H.K., MALALASEKERA W., An introduction to computational fluid


dynamics, the finite volume method. Longman Group Ltd, London, 1995.

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

Vous aimerez peut-être aussi