République Algérienne Démocratique et Populaire

Ministère de l’Enseignement Supérieur

et de la Recherche Scientifique
Université de Tissemsilt

Faculté des Sciences et de la Technologie

Département des Sciences et de la Technologie

Mémoire de fin d’études pour l’obtention du diplôme

de Master académique en
Filière : Génie Mécanique
Spécialité : Installations énergétiques et turbomachines

Présentée par : CHEDNI Abdelmalek

Thème
___________________________________________________________

Transfert de chaleur dans un milieu poreux avec

et sans obstacle

Soutenu le, ………………………..

Devant le Jury :
BENLEKKAM Mohamed Amine Président M.C.A. Univ-Tissemsilt
KHERRIS Sahraoui Encadreur Prof. Univ-Tissemsilt
SAIDI Lilia Co-encadreur Dr. Univ-Tissemsilt
MEDKOUR Maihoub Examinateur M.C.A. Univ-Tissemsilt

Année universitaire : 2023-2024

 1

Dédicace
À l’aide de dieu ’’Allah’’ le tout puissant et miséricordieux
Qui m’a tracé le chemin de ma vie, j’ai pu réaliser ce travail.
Je dédie chaleureusement ce travail aux personnes les plus chères au monde.

Á ma très chère mère,

Vous m’avez comblé avec votre tendresse et affection tout au long de mon
parcours. Que Dieu

Vous donne santé, bonheur et longue vie afin que je puisse vous combler à
mon tour;

À mon cher père ;

Rien au monde ne vaut les efforts fournis jour et nuit pour mon éducation et
mon bien être. Ce travail est le
fruit de tes sacrifices qui tu as consentis pour mon éducation et ma formation,
et puisse dieu t'accorder santé
et longue vie.
À mes très chers frères:
Merci de me suivre et de m’encourager dans mes projets vous me faites
toujours sentir que vous croyez en
moi et votre support fait une grande différence dans ma vie.
À ma sœur : en gage de ma profonde estime pour l’aide que tu m’as
apporté. Tu m’as soutenu, réconforté et encouragé. Puissent nos liens
fraternels se consolider et se pérenniser
encore plus.
À mes chers amis Mahdi, Sifou, Zakit et spécialement à Dendane Ahmed Abdelbasset et
Oussama larkem qui m’ont toujours encouragé, et à qui je souhaite plus de succès

Abdelmalek
2
 3

Remerciements

Au nom d’Allah Clément et Miséricordieux et Louange à Allah,

Seigneur des univers, que la
prière et la paix soient sur notre prophète Mohammad, sa famille
et ses compagnons.
Louange à Allah qui nous a donné la santé, la volonté et la
sagesse de réaliser le présent travail.
Je tiens à remercier très sincèrement mon encadreurs
Monsieur Sahraoui
KHERRIS et Madame Lilia SAIDI, pour leurs suivis et leurs
soutiens
permanents, leurs grandes disponibilités et leurs nombreux
conseils durant toute la période de
l’élaboration du présent travail.
Nous voudrions remercier également les enseignants qui ont bien
voulu examiner le présent
travail et faire partie du jury.
Enfin et surtout, Nous tenons à remercier nos familles qui ont
toujours été présentes et qui nous
ont soutenues tout le long de notre travail.
 4

SOMMAIRE

Nomenclature
Liste des figures
Liste des tableaux
Introduction générale
Chapitre 1 : Etude bibliographique des milieux poreux...............................................................................9
Introduction.............................................................................................................................................10
1.1 Domaines d’application des milieux poreux................................................................................10
1.1.1 Géologie..............................................................................................................................10
1.1.2 Anatomie médicale..............................................................................................................10
1.1.3 Génie civil...........................................................................................................................10
1.1.4 Industrie...............................................................................................................................11
1.1.5 Environnement.....................................................................................................................11
1.2 Propriétés structurelles des milieux poreux.................................................................................11
1.2.1 Volume Élémentaire Représentatif (VER)...........................................................................11
1.2.2 Capillarité............................................................................................................................11
1.2.3 Porosité................................................................................................................................12
1.2.4 Degré de saturation..............................................................................................................12
1.2.5 Teneur en eau.......................................................................................................................12
1.2.6 Perméabilité.........................................................................................................................13
1.2.7 Tortuosité.............................................................................................................................13
1.3 Caractéristiques thermo-physiques des milieux poreux...............................................................14
1.3.1 Conductivité thermique........................................................................................................14
[Link] Modèle de la conductivité thermique effective................................................................15
[Link] Modèle de Maxwell.........................................................................................................16
1.3.2 Viscosité effective...............................................................................................................16
1.3.3 Capacité calorifique équivalente..........................................................................................17
1.4 Modèles d’écoulement dans les milieux poreux (5).....................................................................17
1.4.1 Model de Darcy 1856..........................................................................................................17
1.4.2 Modèle de Brinkman 1947...................................................................................................18
1.4.3 Modèle de Darcy-Forchheimer 1901...................................................................................18
1.4.4 Model général de darcy-brinkman- Forchheimer (dupuit)...................................................19
Conclusion...................................................................................................................................................20
Chapitre 2 : Etude bibliographique des transferts de chaleur dans les milieux poreux...............................21
Introduction.................................................................................................................................................22
2.5 Applications du transfert de chaleur dans les milieux poreux......................................................22
 5

2.5.1 Isolation thermique..............................................................................................................22

2.5.2 Géothermie..........................................................................................................................22
2.5.3 Séchage................................................................................................................................22
2.5.4 Stockage d’énergie thermique..............................................................................................22
2.6 Modes de transfert de chaleur dans les milieux poreux................................................................23
2.6.1 Transfert de chaleur par conduction.....................................................................................23
2.6.2 Transfert de chaleur par convection.....................................................................................24
[Link] Convection naturelle........................................................................................................24
[Link] Convection forcée............................................................................................................24
2.6.3 Transfert de chaleur par rayonnement..................................................................................24
2.7 Effet des obstacles sur le transfert de chaleur..............................................................................25
2.7.1 Transfert de chaleur dans les milieux poreux sans obstacles................................................25
2.7.2 Transfert de chaleur dans les milieux poreux avec obstacles (7)..........................................25
2.8 Revues bibliographiques sur les milieux poreux..........................................................................26
Conclusion...................................................................................................................................................27
Chapitre 3 : Formulation Mathématique....................................................................................................28
Introduction.................................................................................................................................................29
3.1 Modèles physiques......................................................................................................................29
3.1.1 Premier modèle physique....................................................................................................29
3.1.2 Deuxième modèle physique.................................................................................................29
3.2 Hypothèses simplificatrices.........................................................................................................30
3.3 Formulation mathématique..........................................................................................................30
3.4 Equations adimensionnelles du problème....................................................................................31
Chapitre 4 : Résultats.................................................................................................................................33
Introduction.............................................................................................................................................34
4.1 Maillage.......................................................................................................................................34
4.2 Validation némurique..................................................................................................................35
4.3 Résultats et discutions..................................................................................................................35
4.4 Analyse des Tests de Sensibilité..................................................................................................36
4.4.1 Porosité (0.9 à 0.5)...............................................................................................................36
4.4.2 Nombre de Reynolds (5 à 20)..............................................................................................37
Conclusion...................................................................................................................................................37
Conclusion générale..................................................................................................................................38
Bibliographie.............................................................................................................................................40
Bibliographie...........................................................................................................................................41
Résumé et mots clés...................................................................................................................................42
Abstract and keywoeds................................................................................................................................42
 ‫‪6‬‬

‫‪.....................................................................................................................42‬الملخص و الكلمات المفتاحبة‬

‫‪Nomenclature‬‬
 7

3
Volume totale du solide [ m ] Température chaude [K]

3
Volume de pore [m ] vitesse de référence [m/s]

3
Volume de solide [ m ] Composante de la vitesse [m/s]

Vitesse de filtration [m / s] Surface [m²]

Temps [s] Température froide [K]

Teneur en eau [%] Températures [K]

Poids de l’eau [N] Hauteur de la cavité [m]

Accélération de la pesanteur [m / s²] Tortuosité [-]

Coefficient de Forchheimer [-]

Paramètre géométrique de Forchheimer [-]

Pression [Pa]

Viscosité dynamique effective [kg / m s]

Viscosité dynamique du fluide [kg / m s]

Masse volumique [ kg / m3]

3
Masse volumique de fluide [kg / m ]

Conductivité hydraulique [m / s]

Longueur de l’échantillon, de la conduite ou de la cavité [m]

Longueur réelle du courant de fluide [m]

Coefficient d’échange thermique [W/m².K]

 8

Symboles grecques

Porosité

Perméabilité [m²]

Conductivité thermique effective

Conductivité thermique de fluide

Liste des figures

Fig. 1– Représentation schématique d'un VER

Fig.2– représentation géométrique de la notion de tortuosité

Fig.3– modèle en série

Fig.4– modèle en parallèle

Fig.5– cavité poreuse sans obstacle

Fig.6– cavité poreuse avec obstacle

Fig.7– maillage de modèle

Fig.8– variation du champ de température et vitesse pour une porosité de 0.9

Fig.9– variation du champ de température et vitesse pour une porosité de 0.5

Fig.10– profil de vitesse et température pour un nombre de Reynolds de 5

Fig.11–profil de vitesse et température pour un nombre de Reynolds de 20

Fig.12– profil de vitesse et température pour un nombre de Reynolds de 20

Liste des tableaux

Tableau.1– Résumé de quelques modèles de la viscosité effective

 9

Introduction
général
 10

Introduction général

L'étude des milieux poreux est un domaine de recherche passionnant et multidisciplinaire

qui s'intéresse aux propriétés et aux comportements des matériaux constitués de pores
interconnectés. Ces milieux poreux sont omniprésents dans notre environnement, que ce soit dans
le sol, les roches, les matériaux de construction, les membranes de filtration, ou encore les tissus
biologiques. Comprendre les propriétés et les comportements de ces matériaux est donc d'une
importance primordiale pour de nombreuses applications industrielles et environnementales.

Dans ce contexte, cette étude vise à approfondir la compréhension des propriétés thermo-
physiques des milieux poreux, ainsi que des mécanismes de transfert de chaleur et de masse qui
sont associés. Pour ce faire, nous avons exploré les différentes méthodes de simulation numérique
utilisées dans ce domaine, ainsi que les modèles d’écoulement et les données expérimentales les
plus récentes.

La problématique centrale de cette étude est donc de comprendre les mécanismes de

transfert de chaleur et de masse dans les milieux poreux, ainsi que leur impact sur les propriétés
thermo-physiques de ces matériaux. Nous avons donc posé plusieurs questions de recherche :
quelles sont les principales caractéristiques structurelles des milieux poreux, et comment
affectent-elles les transferts de chaleur et de masse dans ces matériaux ? Quels sont les modèles
d’écoulement les plus appropriés pour simuler les écoulements dans les milieux poreux, et
comment peuvent-ils être adaptés en fonction des conditions spécifiques d'une application
donnée ?

Pour répondre à ces questions, nous avons utilisé une méthodologie basée sur la simulation
numérique, en utilisant des outils de modélisation mathématique et de simulation informatique
pour simuler les transferts de chaleur et de masse dans des milieux poreux de différentes tailles et
formes. Nous avons également utilisé des données expérimentales pour valider nos simulations et
pour comparer les résultats obtenus avec les modèles existants dans la littérature.

Le plan de notre étude est le suivant : nous commencerons par une présentation générale des
milieux poreux, en définissant leur nature et leurs domaines d'application. Nous explorerons
ensuite les propriétés structurelles et thermo-physiques des milieux poreux, ainsi que les modèles
d'écoulement les plus couramment utilisés dans ce domaine. Nous nous concentrerons ensuite sur
les mécanismes de transfert de chaleur et de masse dans les milieux poreux, en étudiant les
principaux modes de transfert et en examinant les modèles mathématiques les plus appropriés
pour simuler ces phénomènes. Enfin, nous conclurons en résumant les principaux résultats de
 11

notre étude et en discutant les perspectives de recherche futures dans ce domaine passionnant et
en constante évolution.

Chapitre 1
Etude bibliographique

Des milieux poreux

 12

Introduction
Les milieux poreux sont des matériaux qui comportent des espaces vides, appelés pores,
peuvent être naturels (ex. sol, roches, sable) ou artificiels (ex. mousses, éponges, textiles
technique) à l’intérieur des quels l’eau, l’air.

Ce chapitre traite les milieux poreux et leurs diverses applications, en mettant l’accent sur
leurs propriétés structurelles et leurs caractéristiques thermo-physiques. Les milieux poreux sont
des matériaux qui présentent un réseau de pores interconnectés, offrant ainsi la possibilité
d’absorber, de stoker et de transporter différents fluides. Ces matériaux jouent un rôle crucial dans
de nombreux domaines, tels que la filtration, l’absorption, la catalyse, la géotechnique, l’énergie
thermique et la mécanique des fluides. Le modèle d’écoulement dans les milieux poreux est un
outil essentiel pour comprendre et prédire le comportement de ces matériaux dans diverses
applications.

Ainsi, une compréhension approfondie des milieux poreux et de leur caractéristique est
nécessaire pour exploiter pleinement leur potentiel dans différents domaines.

Domaines d’application des milieux poreux

L’importance des milieux poreux est multiple. En voici quelques exemples :

Géologie
Les milieux poreux sont essentielles dans nombreux domaines, tels que l’exploration et
l’exploitation des ressources naturelles (pétrole, gaz naturel, eau souterraine, etc.), elle comprend
également des aspects géo physiques et géo- mécaniques, tels que l’étude de la porosité (le
volume des pores dans une roche) et de la perméabilité (la capacité d’une roche a permettre le
passage des fluides).

Anatomie médicale
Dans l’anatomie médicale, il existe plusieurs milieux preux qui sont présents à différents
endroits du corps, tels que les poumons, les reins, les vaisseaux sanguins et les ganglions
lymphatiques. Ces milieux poreux dans l’anatomie médicale permettent des fonctions importantes
telles que l’échange gazeux dans les poumons, la filtration rénale, le transport sanguin et la
défense immunitaire.

Génie civil
Dans le domaine du génie civil, les milieux poreux jouent un rôle essentiel dans de
nombreuses application, à la fois sur le plan pratique et théorique. Ils sont couramment utilisés
pour la construction d’infrastructures de drainage, tel que les terrains de sport ou les routes et la
filtration de l’eau ou l’étude des phénomènes physique.
 13

Industrie
Les milieux poreux sont des éléments stratégiques dans l’industrie pour la filtration, la
catalyse, le stockage d’énergie, l’isolation thermique, séparation des mélanges miscibles,
développement de matériaux innovants.

Environnement
Gestion des écosystèmes et des ressources naturelles, traitement des déchets solides et
liquides, évaluation des impacts environnementaux.

Propriétés structurelles des milieux poreux

Les propriétés du milieu poreux varient en fonction de l’échelle de description utilisée, ce
qui a une signification physique macroscopique et des valeurs variables selon les mesures.
D'autre part, des études thermiques en été réalisées pour calculer des paramètres équivalents en
définissant un volume élémentaire représentatif.

Volume Élémentaire Représentatif (VER)

Un milieu poreux est un matériau qui contient un réseau de pores interconnectés, ces
pores peuvent être de différent taille et de forme, avec un diamètre qui vari généralement de 0.05
μm pour les nanopores à 0.5 mm pour les micropores. Dans un milieu poreux, le volume
élémentaire représentatif est un petit volume qui est suffisamment grand pour contenir une
représentation statistiquement significative de la structure et des propriétés du milieu poreux

Fig. 1– Représentation schématique d'un VER [1]

1.1.1 Capillarité
La capillarité dans un milieu poreux se réfère l’altitude d’un liquide de s’infiltrer et de se
déplacer en haut à travers les pores d’un matériau solide en raison de forces d’attraction entre les
 14

molécules du liquide et celles du matériau. Elle dépend de la tension superficielle du liquide, de la

taille des pores et peut avoir des implications importantes dans de nombreux domaines.

1.1.2 Porosité
La porosité est définie comme étant la capacité d’un matériau solide à retenir les fluides
dans ces pores, est aussi une grandeur physique définie comme le rapport entre le volume des

vides et le volume total d'un milieu poreux. La porosité est généralement donnée en
pourcentage (de 0 à 100%)

(1)

Ou : est la porosité, est le volume des pores, est le volume total du matériau c'est-à-
dire la somme du volume de solide et du volume des pores.

1.1.3 Degré de saturation

Le degré de saturation est une mesure qui indique la quantité d’eau présente dans les
espaces poreux d’un matériau par rapport à sa capacité maximale à retenir l’eau. Il est exprimé en
pourcentage et varie de 0% (aucune eau ne présente) à 100% (tous les pores sont remplis d’eau).

(2)

Ou : et le volume d’eau, et le volume des vides.

1.1.4 Teneur en eau

Le teneur en eau d’un matériau est le rapport du poids d’eau contenu dans ce matériau au
poids du même matériau sec.

(3)

Ou : E est le poids d’eau dans le matériau, est poids du matériau sec.

 15

1.1.5 Perméabilité
La perméabilité est un paramètre qui mesure la capacité d'un matériau à laisser passer des
fluides à travers ses pores [2]. Sous l'effet d'un gradient de pression ou d'un champ de gravité.
L’unité de la perméabilité et le Darcy (D), mais il est plus courant d’utiliser le millidarcy (mD)

dans la pratique. κ, ne dépend que de la porosité et de la géométrie de la matrice solide [3].

La perméabilité est estimée par :

Lois de Darcy [3] :

(4)

Ou : q est le débit volumique par unité de surface, ҝ est la perméabilité(en m2ou D)

μ est la viscosité dynamique du fluide, P est le gradient de pression. Cette loi est valide
pour les écoulements laminaires dans des milieux poreux à faible Re.

La corrélation de Kozeny-carmen [4]:

(5)

Ou : (Sv) est la surface spécifique (surface par unité de volume, en m2/m3), (ε) est porosité du
milieu poreux, (Cs) une constante dépendant de la forme des grains (généralement proche de 5
pour les matériaux sphériques). Cette corrélation s’applique principalement a des milieux
constitués de particules sphériques ou presque sphériques et a des écoulements laminaires.

La corrélation de Hazen [5] :

(6)

Ou : d10 est le diamètre effectif des grains (en cm), CH est une constante dépendante des
propriétés du sol, généralement entre 100 et 150 pour les sables propres. Cette corrélation est
simple à appliquer et souvent utilisée dans l’ingénierie géotechnique, bien qu’elle soit limités à
certains types de sols granulaires.

1.1.6 Tortuosité
La tortuosité (ou coefficient de tortuosité) est définie comme le rapport entre la longueur
effective (la longueur réellement suivi entre deux point) de pore Le, et la longueur du milieu
poreux selon l’axe principale de diffusion du fluide L, comme représenté sur la figure.
 16

(7)

Fig.2– représentation géométrique de la notion de tortuosité

1.2 Caractéristiques thermo-physiques des milieux poreux

Les caractéristiques thermo-physiques des milieux poreux décrivent comment ces matériaux
réagissent aux variations de température.

1.2.1 Conductivité thermique

La conductivité thermique est une grandeur intrinsèque du matériau qui caractérise sa
capacité à se laisser traverser par un flux de chaleur. La conductivité dépend de la proportion des
différentes phases mais aussi de la géométrie de leur arrangement. Différents méthodes et modèles
permettent de déterminer la conductivité thermique, parmi eux, le modèle de Maxwell, modèle de
la conductivité thermique effective et le modèle de la conductivité apparente.

Avant de mentionner les modèles de calcul de la conductivité thermiques, il convient de souligner

que quelle que soit la méthode utilisée, la conductivité thermique du matériau se situe toujours
entre deux valeurs extrêmes. La limite inférieure (modèle série, figure 3) correspond à un milieu
ou le vecteur densité de flux de chaleur est perpendiculaire aux strates. [5]

Fig.3– modèle en série

La conductivité thermique s’écrit :

 17

(8)

La limite supérieur (modèle parallèle, figure 4) correspond à une disposition des

Strates parallèlement à la direction de la densité du flux de chaleur.

Fig.4– modèle en parallèle

La conductivité thermique s’écrit :

(9)

[Link] Modèle de la conductivité thermique effective

La conductivité thermique effective dans les milieux poreux représente la capacité d’un
matériau à conduire la chaleur, en tenant compte à la fois du solide et des pores remplis de fluide.
Physiquement, elle quantifie comment la chaleur se propage à travers le matériau, intégrant les
effets de conduction dans le solide et de convection dans les pores.

Ce modèle utilise des relations empiriques ou des modèles physiques plus complexes pour
estimer la conductivité thermique effective d’un milieu poreux en fonction des propriétés
géométriques du réseau poreux, telles que la porosité, la tortuosité, la connectivité, etc.

La conductivité thermique effective s’écrit :

Modèle série

(10)
 18

Modèle parallèle

(4)

Avec :

Ou : ( λ d ) est la conductivité stagnante et ( λ T ) est la conductivité thermique de dispersion. [2]

[Link] Modèle de Maxwell

Ce modèle considère le milieu poreux comme un ensemble de sphères solides dispersées
dans un fluide. Il suppose que la conductivité thermique du milieu poreux est influencée par le
contact entre les particules solides et le fluide.

La conductivité thermique alors calculée à partir de l’addition des contributions de la

conductivité thermique du fluide et de celle des particules solides, pondérées par leur
concentration respective.

(5)

1.2.2 Viscosité effective

La viscosité effective est la mesure de la résistance d’un fluide à l’écoulement, qui tient
compte des effets de surface et de la taille des particules en suspension dans le fluide. Elle est
souvent utilisée pour caractériser la viscosité d’un mélange complexe, comme une suspension ou
une émulsion, ou les particules ont une influence sur la viscosité globale du fluide.

Brinkman a considéré la viscosité dynamique effective comme étant égale à la viscosité du

fluide considéré( μeff =μf ). Lundgren (1972) a montré que la viscosité effective était fonction de la
viscosité du fluide et de la porosité de milieu [6].

Le Tableau (1) résume quelques modèles des viscosités effectives trouvées dans la
littérature.

Tableau.1– Résumé de quelques modèles de la viscosité effective

Modèle Corrélation Limite de validité

 19

Einstein μ eff Cas de sphère

=1−1.25 ε
μ ε →0

Happel μ eff Modèle théorique pour

=exp ( 4.58 ε )
μ des particules sphériques
ε ≤ 0,3

Mooney μ eff
μ
=exp
[ 2.5 (1−ε )
1−K (1−ε ) ] 0,75 < K < 1,5

1.2.3 Capacité calorifique équivalente

La capacité calorifique équivalente (CCE) est la quantité de chaleur nécessaire pour élever
la température d’un certain volume ou masse du matériau d’un degré Celsius. Cela permet de
comparer la capacité calorifique de différentes substances en établissant une équivalence entre
elle, généralement en utilisant l’eau comme référence. La CCE est souvent utilisé en
thermochimie pour calculer la chaleur émise ou absorbée lors d’une réaction chimique, ou en
physique pour mesurer la quantité d’énergie nécessaire pour chauffer un objet.

La capacité calorifique équivalente s'exprime :

(6)

1.3 Modèles d’écoulement dans les milieux poreux [5]

Les modèles d’écoulements dans les milieux poreux sont utilisés pour décrire le mouvement
des fluides à travers des matériaux poreux tels que les sols, les roches ou les filtres. Ces modèles
permettent de comprendre et de prédire les phénomènes d’écoulements et de transport de
substances dans ces [Link] quelques modèles d’écoulement couramment utilisés :

1.3.1 Model de Darcy 1856

En 1856, l'ingénieur français Henry Darcy travaillait pour la ville de Dijon en France sur un
projet impliquant l'utilisation de sable pour filtrer l'approvisionnement en eau, où il a réalisé des
expériences en laboratoire pour examiner les facteurs qui régissent la vitesse d'écoulement de l'eau
dans le sable.
 20

Les résultats de ses expériences ont permis de définir les principes empiriques de base de
l'écoulement des eaux souterraines, qui se traduisent par une équation aujourd'hui connue sous le
nom de loi de darcy.

(7)

Ou : (V fil ) est la vitesse de filtration, ( μf ) est la viscosité dynamique de fluide, κ est la

perméabilité des milieux poreux (⃗
∇ P ¿ est le gradient de pression appliqué.

Remarque : La loi de darcy est valable pour les fluides incompressibles, les milieux poreux
homogènes et isotropes, et lorsque le fluide s’écoule de manière régulière.

1.3.2 Modèle de Brinkman 1947

Dans le cas où la porosité est importante, il faut tenir compte des effets de diffusion
visqueuse au niveau des parois pour cette situation Brinkman (1947) propose de modifier
l’équation de Darcy de la façon suivante:

(8)

Avec : (μeff) représente la viscosité effective du milieu poreux. En général (μeff) est considéré
équivalente à la viscosité dynamique (μ).

1.3.3 Modèle de Darcy-Forchheimer 1901

Lorsque la vitesse d'écoulement augmente, les forces d'inertie ne sont plus négligeables.
Dans ce cas, nous montrons expérimentalement que pour un gradient de pression fixe, le débit
mesuré est plus petit qu'il ne le serait avec la loi de Darcy Forchheimer a proposé une
modification empirique de la loi de Darcy en ajoutant une relation non linéaire du second ordre
entre le gradient de pression et le débit.
 21

Entre le gradient de pression et le taux de filtration. Cette formulation appelée de Darcy-

Forchheimer ;
−μf ρC (9)
⃗
P= V fil −F f f ‖⃗
⃗ V fil ‖V
⃗ fil + ρf ⃗g
k √k

f : est un paramètre caractéristique de la géométrie du milieu poreux, généralement appelé

paramètre de Forchheimer. Muskat propose (1946) une classification du domaine de validité des
modèles de Darcy et Darcy-Forchheimer.

1.3.4 Model général de darcy-brinkman- Forchheimer (dupuit)

On généralise les trois modèles en introduisant un terme transitoire et le terme convectif,

ce qui conduire à l’équation suivante:

[ ( )]
⃗ fil (10)
1 ∂ V fil 1 V ⃗ fil− μf ⃗ ρf C f ⃗ ⃗
ρf + V fil ∇ ⏟
=− ⃗
∇ p + ρ⏟g +⏟
f⃗
2
μeff ∇ V V fil −F ‖V fil‖V fil
⏟
ε ∂ t ε ε
b c d
⏟k ⏟ √ k
e f
a

Avec :

a) Effet d’accélération macroscopique.

b) Gradient de pression.
c) Force volumiques.
d) Terme de Brinkman.
e) Terme de Darcy.
f) Terme de Forchheimer.

Conclusion
Les milieux poreux sont des structures complexes qui jouent un rôle essentiel dans
nombreux domaines, tels que la géologie, l’ingénierie des sols, la filtration de l’eau, et la
production et le stockage d’énergie. Leurs propriétés de porosité et de perméabilité les rendent
adaptés à de multiples applications, en particulier en termes de circulation des fluides.
 22

De plus, leur comportement peut varier en fonction des facteurs environnementaux et des
interactions avec les substances présents dans le [Link] est donc essentiel de développer des
modèles et des approches expérimentales appropriés pour étudier les milieux poreux.

Chapitre 2
Étude bibliographique des
transferts de chaleur dans les
milieux poreux
 23

Introduction

De tous temps, la transmission d’énergie, particulièrement celle de la chaleur, a eu une

grande importance dans divers domaines tels que l'ingénierie, la science des matériaux et la
géologie,

Le transfert de chaleur est décrit le mouvement de la chaleur d’une source a une autre, en raison
d’une différence de température. Il existe trois modes de transfert de chaleur : la conduction, la
convection et le rayonnement.

Le transfert de chaleur dans un milieu poreux est un processus complexe qui dépend de nombreux
facteurs tels que la conductivité thermique du matériau, la porosité et la géométrie du milieu. Le
transfert de chaleur se produit principalement par conduction et convection.

Dans ce chapitre, nous allons discuter du transfert de chaleur dans les milieux poreux avec
et sans obstacles, qui est un sujet crucial dans de nombreux domaines.

Historique

L’historique de recherche sur les milieux poreux remonte au début du XXe siècle, avec des
contributions significatives dans les domaines de la géologie, de l’ingénierie et de la physique
[10][11}.

Début du XXe siècle Les travaux de darcy (1856) sur le flux dans les milieux poreux ont
posé les bases de l’hydrodynamique dans les sols.

Années 1950-1960 Les recherches se sont intensifiées avec l’application des principes de la
mécanique des fluides et de la thermodynamique pour modéliser les transferts de chaleur et de
masse.

Années 1980-1990 Développement des modèles numériques et des simulations

informatiques pour étudier les propriétés thermiques des milieux poreux, avec l’émergence de la
mécanique des milieux continus.

Années 2000 a aujourd’hui intérêt croissant pour les applications environnementales et

énergétique, incluant la géothermie et le stockage de CO2, avec une recherche multidisciplinaire
intégrant la chimie, la physique et l’ingénierie.

Application pratique
 24

Les transferts de chaleur dans les milieux poreux ont de nombreuses applications pratiques
dans divers domaines. Voici quelques exemples significatifs :

1. Géothermie [6] Utilisation de l’énergie géothermique pour le chauffage ou la

production d’électricité. Les modèles de transfert de chaleur permettent d’évaluer
l’efficacité des systèmes de chauffage géothermique.
2. Stockage Thermique [7] Systèmes de stockage d’énergie thermique dans des
matériaux poreux, comme le stockage de chaleur dans le béton ou les aquifères. Ces
systèmes optimisent l’utilisation de l’énergie dans les bâtiments.
3. Systèmes de Récupération de chaleur [8] Utilisation de milieux poreux pour
récupérer la chaleur des eaux usées dans des systèmes de traitement. Cela améliore
l’efficacité énergétique des stations d’épuration.
Concepts fondamentaux des transferts de chaleur [9][10]
 Conduction Transfert de chaleur par contact direct entre molécules, décrit par la loi
de Fourier.
 Convection Transfert de chaleur par le mouvement des fluides, divisé en convection
naturelle et forcée, souvent modélisé par les équations de Navier- Stocks.
 Rayonnement Transfert de chaleur par émission d’énergie électromagnétique, régi
par la loi de Stefan-Boltzmann.

1.4 Modes de transfert de chaleur dans les milieux poreux

La conduction thermique dans le solide peut être décrite par l’équation de la chaleur

 Conduction

Ou et la température du solide et est la diffusivité thermique du matériau.

Transfert par convection

Pour le fluide dans les pores, la convection est généralement modélisée par

 Equation de convection

Où est la température du fluide, est le vecteur vitesse du fluide, et est la

diffusivité thermique du fluide.
 25

 Equation globale du transfert de chaleur

L’équation générale peut être exprimée comme

Ou, est la température effective, est la conductivité thermique effective des

milieux poreux peut être estimée comme :

et Q représente les sources ou puits de chaleur.

1.5 Effet des obstacles sur le transfert de chaleur

La présence d'obstacles dans un milieu poreux peut avoir un effet significatif sur le transfert
de chaleur. Les obstacles peuvent perturber le mouvement des fluides, ce qui peut réduire la
convection et augmenter la résistance thermique du matériau.

Cependant, les obstacles peuvent également créer des zones de turbulence qui peuvent améliorer
le transfert de chaleur par convection. La forme, la taille et la distribution des obstacles sont donc
des facteurs importants à considérer lors de l'analyse du transfert de chaleur dans un milieu poreux
avec des obstacles.

1.5.1 Transfert de chaleur dans les milieux poreux sans obstacles

Dans les milieux poreux sans obstacles, le transfert de chaleur se produit principalement par
conduction thermique.

La conductivité thermique dépend de la structure poreuse du matériau et de ses propriétés

thermiques intrinsèques.

1.5.2 Transfert de chaleur dans les milieux poreux avec obstacles [11]

Lorsque des obstacles sont présents dans les milieux poreux, le transfert de chaleur peut se
produire par convection forcée ou naturelle

Les obstacles peuvent modifier la structure poreuse du matériau, affectant ainsi sa conductivité
thermique et son coefficient de transfert de chaleur.

Dans un milieu poreux avec des obstacles, tels que des particules ou des fibres, le transfert de
chaleur peut être plus complexe. Les obstacles peuvent perturber le flux de chaleur et créer des
zones de convection thermique.
 26

Pour modéliser le transfert de chaleur dans ces milieux, des équations aux dérivées partielles
couplées doivent être résolues numériquement. Des méthodes telles que la méthode des volumes
finis et la méthode des éléments finis peuvent être utilisées pour résoudre ces équations.

1.6 Revues bibliographiques sur les milieux poreux

 « Transport, reactions and scaling in porous media » par L. Bazant, publié dans la revue Reviews
of Modern Physics (2002).

Cette revue présente les derniers développements dans le domaine de la modélisation des
processus de transport et de réaction dans les milieux poreux, ainsi que les méthodes analyse et les
perspectives futures.

 « Characterization of Porous Media» par Arthur G. Palmer, publié dans la revue Reviews of
Geophysics (1991).

Cette revue se concentre sur les techniques de caractérisation des roches et des sols poreux,
en mettant l’accent sur l’imagerie par tomographie et sur les mesures de propriétés physiques.

 « Pore-scale modeling of multiphase flow in porous media » par D. A. DiCarlo, publié dans la
revue Current Opinion in Chemical Engineering (2017).

Cette revue présente les derniers développements en matière de modélisation numérique des
phénomènes de transport dans les milieux poreux, en mettant l’accent sur les aspects multiphases
et sur la simulation des écoulements à l’échelle des pores.

« Advances in the understanding of fluid flow in fractured and/or porous media » par H. Miller,
publié dans la revue Advances in Water Resources (2012).

Cette revue s’intéresse aux avancées récentes dans la compréhension de l’écoulement des fluides
dans les milieux poreux fracturés, en mettant l’accent sur les méthodes d’analyse des données de
terrain et sur les simulations numériques.

 « Porous Media From Microstructure to Transport Properties »par F. Nice et al, publié dans la
revue Annual Review of Chemical and Biomolecular Engineering (2013).

Cette revue examine les différentes approches pour modéliser les propriétés de transport
 27

Dans les milieux poreux, en mettant l’accent sur les relations entre la microstructure et les
propriétés macroscopiques.

Conclusion

Comprendre et modéliser le transfert de chaleur dans les milieux poreux est essentiel pour
diverses applications en ingénierie, en géologie, en science des matériaux et dans nombreux autres
domaines. Le transfert de chaleur dans les milieux poreux avec des obstacles est un domaine de
recherche important qui a de nombreuses applications pratiques. La compréhension des effets des
obstacles sur le transfert de chaleur peut aider à améliorer l'efficacité énergétique et à développer
de nouvelles technologies.
 28

Chapitre 3
Formulation
Mathématique
 29

Introduction
Dans le présent chapitre, nous décrivons et formulons le modèle physique, les hypothèses
simplificatrices et les conditions aux limites appropriées.

1.7 Modèles physiques

1.7.1 Premier modèle physique

Le modèle physique que nous allons étudier est une cavité poreuse rectangulaire de hauteur
H et de langueur L saturée par un fluide incompressible. L’une des parois verticales est chauffée,
tandis que l’autre est refroidie, les parois horizontales étant considérées comme adiabatiques.

Fig.5– cavité poreuse sans obstacle

1.7.2 Deuxième modèle physique

Le problème consiste un transfert de chaleur dans un milieu poreux avec un obstacle chaud

Fig.6– cavité poreuse avec obstacle

L’objectif de cette étude et de déterminé les transferts de chaleur et de masse induits par les
mouvements par convection naturel.
 30

1.8 Hypothèses simplificatrices

Les hypothèses simplificatrices utilisées dans notre étude des phénomènes de convection
naturelle sont les suivant :

 La matrice poreuse est supposée isotrope, homogène, perméable et en équilibre thermodynamique

avec le fluide.
 Le fluide newtonien.
 Le transfert thermique par rayonnement est négligeable.
 L’écoulement du fluide à l’intérieur de la cavité est laminaire et bidimensionnelle.
 L'air initial à l'extérieur de la cavité est considéré être stagnant.

1.9 Formulation mathématique

Comme on l’a déjà signalé, Le modèle physique considéré consiste en une cavité poreuse
horizontale, homogène, isotrope et saturé par un fluide incompressible. Les faces latérales sont
maintenues à des températures différentes alors que les faces horizontales sont adiabatiques. Dans
ce cas une convection naturelle laminaire bidimensionnelle s’établit dans l’enceinte poreuse.

L’ensemble des équations gouvernant la conservation de la masse, la quantité de

mouvement et d’énergie s’écrivent sous forme adimensionnelle comme suit :

Equation de continuité

(11)

Equation de quantité de mouvement suivant l’axe X

[
ρf ∂u
ε ∂t
+u
1 ∂u 1 ∂v ∂ p
ε ∂x
+v
ε ∂ y ∂x]
+ =μ f ⃗
2 μf
∇ u− u−F
ҝ
ρf C F
√ҝ
‖⃗
V fil ‖u
(12)

Equation de quantité de mouvement suivant l’axe Y

[
ρf ∂u
ε ∂t
+u
1 ∂u
ε ∂x
+v
1 ∂v ∂ p
ε ∂y ∂y ]
+ =μ f ⃗
2 μf
∇ v− v−F
ҝ
ρf CF
√ҝ
‖⃗
V fil ‖v +
(13)

ρ f [ 1−β f ( T −T f ) ] g

Equation de conservation d’énergie

 31

[ ]
∂T ∂T ∂T λ eff 2
∂T ∂ T
2 (14)
+u +v = +
∂t ∂x ∂ y ( ρC p )eff ∂ x 2 ∂ y 2

La vitesse V fil est la vitesse de filtration ; les grandeurs P et T sont respectivement les
grandeurs moyennées de la pression et de la température; ε est la porosité du milieu poreux ; ρ f la
masse volumique du fluide; µf et µeff sont respectivement la viscosité dynamique du fluide et du
milieu poreux; k eff la conductivité thermique,  le coefficient d'expansion volumique thermique
du fluide et g représente l’accélération de la pesanteur.

1.10 Equations adimensionnelles du problème

Pour rendre les équations précédentes adimensionnelles, elles seront modifiées en utilisant
les relations suivantes :

x y u v p T −T F
X= , Y= , U= , V= , P= 2 et θ=
L H U0 V0 ρU 0 T C −T F

Les équations adimensionnelles régissant le problème sont les suivantes :

∂U ∂V (15)
+ =0
∂ X ∂Y

( )
∂U ε ∂V ε ∂ ∂
2 2
∂ θ Ra (16)
+ = 2
+ 2 −ε ε +
∂X ∂V ∂ X ∂Y ∂ X Pr

( )
∂ Uθ ∂Vθ ∂ ∂ θ
2 2 (17)
+ = 2
+ 2
∂ X ∂Y ∂ X ∂ Y Pr

Le nombre de Nüsselt local est déterminé par :

Nuloc=−
∂θ
∂Y | y=0
(18)
 32

Le calcul du nombre de Nüsselt moyen est donné par la formule suivante :

1
(19)
NuMoy =∫ Nu ( X ) dX
0

Conclusion

Dans ce chapitre, nous avons vu le problème physique régissant le phénomène de

convection naturel dans une cavité poreuse avec et sans obstacles et les équations utilisées pour
simuler ce phénomène
 33

Chapitre 4
Résultats
et discutions
 34

Introduction
Dans ce chapitre, nous examinerons les résultats obtenus de la résolution numérique du
phénomène de convection naturelle dans une cavité poreuse avec et sans obstacle. En utilisant la
méthode des éléments finis (MEF) et le logiciel comsol-multiphysiques.

Parmi les paramètres régissant le système d’équations, on retrouve le nombre de Darcy (Da)
qui représente la perméabilité du milieu, le nombre de Rayleigh (Ra) qui caractérise la différence
de température et en fin le nombre de prandtl (Pr) qui décrit la nature de fluide.

1.11 Maillage
Le processus de maillage revêt une importance cruciale lors de l'analyse CFD, car il influe
directement sur la précision des résultats obtenus. Bien qu'il n'existe pas de méthode universelle,
quelques principes fondamentaux peuvent guider la création d'un maillage adapté, résumés par
l'acronyme QRLT :

Qualité : Garantir une qualité optimale des éléments pour une précision maximale.

Résolution : Assurer une résolution adéquate dans les régions présentant des gradients élevés pour
capturer les phénomènes locaux.

Lissage : Veiller à un bon lissage entre les zones de maillage fin et grossier afin de réduire au
minimum les erreurs.

Total : Minimiser le nombre total d'éléments pour des temps de calcul raisonnables.

Fig.7– maillage de modèle

De plus, nous avons examiné l'effet de la taille du maillage sur le modèle numérique de
convection naturelle dans une cavité carrée poreuse saturée d'air. Nous avons généré quatre
maillages uniformes de tailles différentes : 20x20, 40x40, 60x60 et 80x80.

Nous avons constaté que le maillage de 60x60 offre une résolution suffisante pour être
utilisé tout au long de la simulation.
 35

1.12 Validation némurique

La simulation numérique a été réalisée à l’aide du code COMSOL-MultiPhysics Pour la
validation numérique du code de calcul, on procède à la comparaison des résultats numériques de
la présente étude avec ceux obtenus par S. Chatti, C. Ghabi and A. Mhimid (8), en régime
permanant, l'écoulement en convection naturelle à travers un milieu poreux saturé de fluide dans
une enceinte carrée, en utilisant le modèle d'écoulement de Brinkman.

1.13 Résultats et discutions

Cette étude se concentre particulièrement sur l'influence de la porosité et du nombre de
Reynolds sur les phénomènes de transfert thermique et de dynamique des fluides dans des
configurations complexes.

Fig.8– variation du champ de température et vitesse pour une porosité de 0.9

Cette figure illustre clairement l'impact d'une porosité élevée sur la distribution de la
température, mettant en évidence une conduction thermique accrue à travers le milieu.
 36

Fig.9– variation du champ de température et vitesse pour une porosité de 0.5

La réduction de la porosité accentue les gradients thermiques, ce qui est crucial pour les
applications nécessitant une isolation accrue.

Fig.10– profil de vitesse et température pour un nombre de Reynolds de 5

Cette figure montre une faible turbulence dans le fluide, caractéristique des écoulements à
faible Reynolds dans un milieu poreux.

Fig.11–profil de vitesse et température pour un nombre de Reynolds de 20

 37

L'augmentation du nombre de Reynolds révèle des comportements de flux plus complexes

et des turbulences accrues, illustrant l'effet du régime inertiel.

1.14 Analyse des Tests de Sensibilité

1.14.1 Porosité (0.9 à 0.5)

L'analyse de la sensibilité à la porosité révèle que la variation de la porosité influence
significativement le transfert de chaleur et la dynamique des fluides dans un milieu poreux avec
un obstacle. Une porosité plus faible (0.5) entraîne une augmentation de la conductivité thermique
du milieu en raison de la proportion plus élevée de matériau solide, ce qui améliore le transfert de
chaleur par conduction. Cela provoque une augmentation de la température du fluide. En
revanche, une porosité plus élevée (près de 0.9) réduit la résistance au flux (augmentation de la
vitesse en haut et en bas du bloc), permettant une circulation plus facile du fluide et facilitant le
transfert de chaleur par convection. Par conséquent, la variation de la porosité modifie l'équilibre
entre la conduction et la convection, ce qui influence la distribution de la température et le
comportement des fluides autour de l'obstacle.

1.14.2 Nombre de Reynolds (5 à 20)

L'analyse de la sensibilité au nombre de Reynolds montre que l'augmentation du nombre de
Reynolds (de 5 à 20) intensifie la turbulence du flux, augmentant ainsi la vitesse du fluide et le
coefficient de convection. Un régime de flux plus turbulent améliore le taux de transfert de
chaleur par convection, ce qui peut entraîner une diminution de la température moyenne du milieu
en raison de l'évacuation plus rapide de la chaleur. À faible nombre de Reynolds, le flux demeure
laminaire, favorisant un transfert de chaleur plus prévisible mais potentiellement moins efficace.

L'interaction entre la porosité et le nombre de Reynolds est cruciale pour optimiser les
conditions thermiques et fluidiques, permettant de concevoir des systèmes adaptés à des
applications spécifiques nécessitant un contrôle précis de la température et du flux de fluide.

Exemple sans obstacle

 38

Fig.12– profil de vitesse et température pour un nombre de Reynolds de 20

Dans un milieu poreux sans obstacle, la convection naturelle serait moins entravée,
permettant une distribution plus uniforme de la température et des vitesses de fluide accrues. Ces
conditions favorisent un transfert de chaleur plus efficace, mais peuvent aussi nécessiter des
ajustements dans la conception des systèmes pour gérer l'augmentation potentielle de la
convection.

Conclusion
Dans ce chapitre, nous avons vu les résultats obtenus de la simulation numérique dans une
cavité poreuse avec et sans obstacle, en régime permanent. Nous avons confronté nos résultats à
ceux obtenus par d’autres chercheurs.
 39

Conclusion
générale
 40

Conclusion générale

L'étude des milieux poreux s'avère être un domaine de recherche captivant et

multidisciplinaire d'une importance capitale pour de nombreuses applications industrielles et
environnementales Dans cette étude, nous avons cherché à comprendre les mécanismes
complexes de transfert de chaleur et de masses dans les milieux poreux, ainsi que leur influence
sur les propriétés thermo-physiques de ces matériaux.

Il apparait clairement que la structure des milieux poreux- notamment la porosité, la

tortuosité, la taille et la distribution des pores- joue un rôle fondamental dans les phénomènes de
transport. Ces propretés influencent directement la manière dont la chaleur et les fluides se
déplacent à travers le matériau, et donc les performances globales des systèmes qui dépendent de
ces transferts, tels que les réacteurs, les échangeurs de chaleur et les systèmes de filtration.

L'analyse des résultats obtenus met en lumière l'influence déterminante de la porosité et du

nombre de Reynolds sur le transfert de chaleur et la dynamique des fluides dans les milieux
poreux. Une porosité plus élevée favorise la convection, tandis qu'une porosité plus faible accroît
la conduction. L'augmentation du nombre de Reynolds intensifie la turbulence du flux, améliorant
ainsi le transfert de chaleur par convection.

Ces résultats soulignent l'importance de considérer conjointement la porosité et le nombre

de Reynolds lors de la conception de systèmes impliquant des milieux poreux. Une
compréhension approfondie de ces paramètres permet d'optimiser les performances thermiques et
fluidiques, menant à des applications plus efficaces et innovantes.

L'analyse des résultats obtenus a révélé l'influence déterminante de deux paramètres clés sur le
comportement des milieux poreux : la porosité et le nombre de Reynolds.

 La porosité, définie comme la proportion de vide dans le matériau, joue un rôle crucial
dans la distribution de la chaleur. Une porosité élevée favorise la convection, permettant
une circulation aisée du fluide et un transfert de chaleur plus efficace par ce mode. À
l'inverse, une porosité plus faible accroît la conduction, le matériau solide conduisant la
chaleur de manière plus directe.
 Le nombre de Reynolds, caractérisant le régime d'écoulement du fluide, influence la
turbulence et le coefficient de convection. L'augmentation du nombre de Reynolds
intensifie la turbulence, améliorant ainsi le transfert de chaleur par convection.

Ces résultats soulignent l'importance de considérer conjointement la porosité et le nombre de

Reynolds lors de la conception de systèmes impliquant des milieux poreux. En optimisant ces
 41

paramètres, il est possible de contrôler les transferts de chaleur et de masse, menant à des
applications plus performantes et innovantes.
 42

Bibliographie
 43

ANNEXES
 44

Annexe A

Nombre de Reynolds dans un milieu poreux :

 Régime laminaire : <1

 Régime turbulent : >10

 45

Bibliographie
[1]. Melle SAIDI Lilia , 2023. Contribution à l’étude numérique de la convection des
écoulements de nanofluide dans les milieux poreux, Université Ibn Khaldoun de Tiaret, Algérie ;

[2]. Ramdane YOUNSI, 2002. SIMULATION NUMERIQUE DU TRANSFERT DE CHALEUR ET DE

MASSE EN MILIEUX FLUIDES ET POREUX, UNIVERSITE DES SCIENCES ET DE LA
TECHNOLOGIE HOUARI BOUMEDIENNE. Algérie ;

[3]. BOUABDALLAH, BOUJDERDA. Amélioration d'un système de chauffage passif par l'insertion d'un
matériau poreux . jijel : s.n., 12/09/2021 ;

[4]. FILALI, M. MOHAMED. Conductivité thermique apparente des milieux granulaire soumis à des
contraintes mécaniques : modélisation et mesures. TOULOUSE : s.n., 24 Février 2006 ;

[5]. CHEIKH, M. Kamel BEN. Etude du transfert de chaleur dans une conduite courbée annulaire
partiellement poreuse . 05/12/2006 ;

[6]. ZAOUIDI, DJELID. : Etude des transferts de chaleur et de masse dans un milieu poreux saturé par un
fluide ou mélange d’un fluide non-Newtonien. ;

[7]. From Rayleigh-bénard convection to porous media convection : how porosity affects heat transfert and
flow structure. cambridge university : s.n., 18-05-2020 ;

[8]. S. Chatti, C. Ghabi and A. Mhimid. 2016. Fluid Flow and Heat Transfer in Porous Media and Post
Heated Obstacle. Int. J .ht 34, 0392-8764 ;

[3]. Bear,J. (1972). Dynamics of fluids in Porous Media. Dover pulication.

[1]. Cheng, L., et al. (2017). Advances in Thermal Conductivity of Porous Media. Applied Energy,
194, 172-185.

[1]. Incropera, F. P.,& DeWitt, D.P. (2007). Fundamentals of Heat and Mass Transfer. Wiley.

[1]. Bergman, T. L, & Lavine, A.S. (2011). Introduction To Heat transfer. Wiley.

[1]. Lund, J. W., Freeston, D. H., & Boyd, T. L. (2010). Direct Utilisation of Geothermal Energy
2010 Worldwide Review. Geothermics, 39(3), 159-180.

[1]. Zhao, X., & Zhang, L. (2012). Thermal Properties of Porous Materials for Thermal Energy
Conversion and Management, 65, 9-15.

[1]. Zhao, X., & Zhang, L. (2015). Heat Transfer Performance of Porous Media in Wastewater
Heat Recovery Systems. Journal of Envirenmental Management, 160, 21-30.

[5] Hazen, A.(1892). Some physical propeties of sands and gravels, with special reference to their
use in filtration. Massachusetts State board of Health.
46
 47

Résumé et mots clés

Cette étude vise à développer et à valider une modélisation numérique pour étudier la
convection naturelle dans les cavités poreuse, avec et sans obstacle, en régime permanent. Pour
cela, une configuration expérimentale a été mise en place où l'une des parois verticales est
soumise à un rapport calorifique, tandis que l'autre est parfaitement isolée thermiquement. Les
parois horizontales sont maintenues à une température inférieure. Le fluide traverse la cavité, ce
qui permet d'étudier les écoulements dans les pores. La simulation numérique a été réalisée en
utilisant la Méthode des Eléments Finis sous le logiciel Comsol-MultiPhysics. Les résultats de
simulation ont été validés et interprétés pour comprendre les mécanismes de convection naturelle
dans les cavités poreuse, avec et sans obstacle.

Mots clés: cavité, convection naturelle, milieu poreux.

Abstract and keywoeds

This study aims to develop and validate a numerical model for studying natural convection
in porous cavities, with and without an obstacle, in a permanent regime. To achieve this, an
experimental setup was designed where one vertical wall was subjected to a thermal gradient,
while the other was perfectly thermally isolated. The horizontal walls were maintained at a lower
temperature. The fluid flows through the cavity, allowing the study of fluid flow in the pores.
Numerical simulation was performed using the Finite Element Method under the Comsol-
MultiPhysics software. The simulation results were validated and interpreted to understand the
mechanisms of natural convection in porous cavities, with and without an obstacle.

Keywords: Cavity, Natural convection, Porous medium

‫الملخص و الكلمات المفتاحبة‬

‫هذه الدراسة تهدف إلى تطوير و التحقق من نموذج رقمي لدراسة التبادل الطبيعي في‬
‫ تم إعداد تركيب تجريبي حيث‬،‫ لتحقيق ذلك‬.‫ في حالة دائمة‬،‫ مع وجود و بدون عائق‬،‫الخزانات المائية‬
‫ تم‬.‫ بينما تم عزل األخرى حراريًا بشكل كامل‬،‫تم تطبيق عائق حراري على أحد الجدران العمودية‬
‫ مما يسمح بالدراسة‬،‫ يمر المائع عبر الخزان‬.‫الحفاظ على الجدران األفقية على درجة حرارة اقل‬
‫تحت‬MEF) ( ‫ تم إجراء التمثيل باستخدام طريقة العناصر المحدودة‬.‫لحركة المائع في الفجوات‬
‫ مع وجود و بدون‬،‫ لالهتمام بالتجارب الطبيعية في الخزانات المائية‬Comsol- MultiPhysics ‫برنامج‬
.‫عائق‬

.‫ وسط مسامي‬،‫ حمل حراري طبيعي‬،‫ تجويف‬:‫الكلمات المفتاحية‬