Le 3ème Séminaire International sur les Énergies Nouvelles et

Renouvelables
The 3nd International Seminar on New and Renewable
Energies
Unité de Recherche Appliquée en Énergies Renouvelables,
Ghardaïa – Algérie 13 et 14 Octobre 2014

Étude expérimentale et numérique d’un

échangeur de chaleur air/sol horizontal
M. E. Benhamza1; A. Brima2; S. Houda3; L.Serrir4
1
Université de Biskra département de génie mécanique
2
Université de Biskra département de génie mécanique
3
Université de Biskra département d’architecture
4
Centre de recherche appliquée en énergies renouvelables- Ghardaïa
1
[email protected]
2
[email protected]
3
[email protected]

énergétique.Cet article traite le comportement thermique d’un

échangeur de chaleur air/sol horizontalqui se compose d'une
Résumé :
canalisation tubulaire souterraine en métal ou enmatière
Les échangeurs de chaleur air/sol horizontaux sont des plastique, dans notre cas il a été utilisé un système de
dispositifs utiliséspour réduire la consommation d’énergie dans canalisation circulaire deØ110mmenterré horizontalement. Ce
les bâtiments, par l’échange de l’énergie thermique de la terre modèles permet l’échange de la chaleur de fluide avec le sol
afin de rafraîchir en été et chauffer en hiver. Dans cette étude
dans le butde rafraîchir en été et de chauffer en hiver. Le
nous avons procédé à une simulation par le biais de la mécanique
des fluides numérique (CFD), pour cela il a été nécessaire de dispositif expérimental est installé au niveau de l’université de
concevoir l’échangeur tout en respectant le design et les Biskra située dans zone a climat chaud, qui est employé pour
dimensions réelles du dispositif expérimental. Pour cet objectif fournir des données expérimentales a été installé au sud de
nous avons utilisé aussi bien pour la conception et le maillage le
logiciel mailleur «GAMBIT» tout en prêtant attention aux l’Algérie et bien précisément au niveau de la ville de Biskra,
conditions aux limites adéquates. Nous avons testé plusieurs pour lebutd’étudier son efficacité.
configurations de maillage avec différents raffinements pour en Dans la littérature, plusieurs travaux de recherche
choisir une qui fasse un compromis entre la précision et les
théoriques et expérimentaux sur les échangeurs de chaleur
ressources informatiques et aussi nous avons raffiné grandement
notre maillage à proximité de la paroi. Concernant le solveur, enterrée ont été menés dans plusieurs laboratoires de
nous avons en recours au logiciel commercial FLUENT 6.3. On a recherches.expliquent les divers types des échangeurs de
confronté les résultats de simulation avec les résultats du chaleur et étudient l’efficacité de chauffage et de
dispositif expérimental, on a constaté une bonne concordance
entre eux. Il est à noter que pour évaluer la température le long
refroidissement.Sehli et al. [1]ont développé un modèle
des deux tronçons (entrée et sortie) de l’échangeur, nous avons numérique régulier unidimensionnel pour estimer l'efficacité
opté une fonction (UDF) la fonction définir par l’utilisateur. des échangeurs de chaleur terre-à-air, installés au sud de
l’Algérie à des profondeursdifférentes. Les auteurs ont
Mots-Clés: Échangeurs de chaleur air/sol, Géothermie, conclus que le système d'EAHE ne peut pas maintenir le
Étudenumérique.
confort thermique à l'intérieur des locaux, mais il pourrait être
I.INTRODUCTION : employé pour réduire la demande énergétique dans le
bâtiment domestique dans la même région. Ce système peut
Les échangeurs de chaleur enterrés sont comme des méthodes être utilisé en même temps comme un dispositif de
passives réalisées pour réduire la consommation climatisation.Misraet al. [2] ontdéveloppé un modèle sous
solveurFLUENT 6.3,L'analyse basée sur CFD a été utilisée

1
 Le 3ème Séminaire International sur les Énergies Nouvelles et
Renouvelables
The 3nd International Seminar on New and Renewable
Energies
Unité de Recherche Appliquée en Énergies Renouvelables,
Ghardaïa – Algérie 13 et 14 Octobre 2014

pour résoudre le champ de la température autour de la

conduite enterrée horizontalement de EATHE, en utilisant
une grille non structurée, et pour étudier les effets de la durée
de temps de l'opération continue, de la conductivité
thermique de sol du diamètre de conduite et de la vitesse
d'écoulement sur l'exécution thermique d'EATHE sous la
coupure des conditions ont été analysées.Un modèle
numérique du transfert thermique dans la terre a été
développé parEsen et al [3] afin de déterminer la distribution
de la température à proximité de la conduite.
L'approximation de différence finie est employée pour
Fig.1Schéma du système d’échangeur de chaleur
l'analyse numérique.il est à noter les résultats numériques air/sol
Pour
obtenus sont en accord les mesures
différents types de maillage et précision (les espaces entre les
expérimentales.Pourréduire la chargederefroidissementdes
nœuds) et avec un maillage de la couche limite c.-à-d. un
bâtiments en été Bansal et al. [4]ont réalisés plusieurs travaux
maillage fin très proche de la paroi du tube, et définie les
sur les échangeurs de chaleur air-sol, l’un de ces travaux un
conditions de frontière (définie les parois comme un solide, et
model a été développé par la CFD avec le solveurFLUENT
l’entrée comme un ventilateur pour assurer la circulation de
(version 6,3), pour un climat chaud et secil s’agit de la ville
l’air,…).
d'Ajmer (Inde).Bansal et al. [5]ontdéveloppés un échangeur
de chaleur air-solet ont intégrés un refroidisseur évaporatif à
la sortie de l'EATHE. L'analyse horaire d’une année du
système intégré a été effectuée pour des conditions
climatiques principalement chaudes et sèches en utilisant le
modèle de la (CFD) modelé avec le solveurFLUENT
(version 6,3).

II. DESCRIPTION DU MODÈLE DE CFD :

Nous avons fait une étude numérique d’un échangeur

de chaleur air/sol.Lesolveur 6,3 FLUENTS de CFD a été
employé dans cette étude.Pour contrôler le changement de la
température tout au long des tubes de l’échangeur, on a
dessiné la géométrie de l’échangeur sur le logiciel de maillage
Gambit2.3 et on a respect toutes les dimensions (la
longueur,diamètredes pipes, la forme des coudes, l’épaisseur
Fig.3Différents types de maillage utilisé pour la
et l’espacement des tubes), pour les deux cas (une forme simulation
complète et autre symétrie).
La même chose avec une géométrie symétrique utilisée
A. La géométrie et le maillage : pour réduire le temps de convergence de la simulation.
Dans cette étape nous avons discrétisé la géométrie de
l’échangeur de chaleur avec des dimensions réelles (longueur,
diamètre…),en utilisant unegéométrietridimensionnelle
adaptée à la structure maillée.

2
 Le 3ème Séminaire International sur les Énergies Nouvelles et
Renouvelables
The 3nd International Seminar on New and Renewable
Energies
Unité de Recherche Appliquée en Énergies Renouvelables,
Ghardaïa – Algérie 13 et 14 Octobre 2014

dissipation, et un l'algorithme SIMPLE standard pour le

couplage vitesse- pression.En ci que concerne le critère de
convergence nous avons choisi 1e-6pour l’énergie et les
autres équations pour assurer la convergence.

C. Équations régissantes :

Les équations régissantes de transport dans les

coordonnées 3D cartésiennes pour un écoulement
incompressible et stationnaire, le transfert de la chaleur et de
la masse sont donnés ci-dessous équation de continuité:

Fig.4Géométrie symétrique de l’échangeur sous Gambit

Nous pouvons aussi controler notre maillagevolumique

qui touche toute la surface de l’échangeur de chaleur en [ ] [ ]
utilisant « examin mesh » et vérifier le maillage selon les axes
« x,y,et z ».

[ ] [ ]

[ ]

[ ]

Équation d’énergie, pour un écoulement stationnaire,

sans source interne et dissipation visqueuse.

Fig.5 Contrôle du maillage. [ ] [ ]

B. Le solveur Fluent 6.3: Les écoulements turbulents sont caractérisés par les
champs de fluctuations de vitesse. Ces fluctuations mélangent
Pour la simulation numérique on a exporté la géométrie vers
le code FLUENT et communiqué les données nécessaires et des quantités transportées telles que l'énergie, et la
définir les conditions aux limites (déclarer que la température concentration d'espèce, aussi bien ont flotté les quantités
de la paroi du tube horizontale de l’échangeur est égale à celle transportées Puisque ces fluctuations peuvent être de petite
du sol à une certaine profondeur), et nous avons utilisé le échelle et de haute fréquence. Cependant, les équations
modèle de l’énergie pour définir les paramètres liés à modifiées contiennent des variables inconnues additionnelles,
l'énergie ou au transfert thermique, et définir la conductivité et les modèles de turbulence sont nécessaires pour déterminer
thermique de la conduite.Les schémas de discrétisation
ces variables en termes de quantités connues:
utilisés dans le présent travail 1er ordre pour la quantité de
mouvement,énergie cinétique turbulente, et le taux de

3
 Le 3ème Séminaire International sur les Énergies Nouvelles et
Renouvelables
The 3nd International Seminar on New and Renewable
Energies
Unité de Recherche Appliquée en Énergies Renouvelables,
Ghardaïa – Algérie 13 et 14 Octobre 2014

En ci que concerne les modèles de turbulences on a Représente la contribution de la dilatation de

plusieurs modèles mathématiques sont présenté, le modèle fluctuation dans la turbulence compressible au taux global de
Spalart-Allmaras, le modèle (k- ε), le modèle (k- ω), le dissipation :
modèle RSM (Reynolds Stress Modèles), modèle (DES),
modèle (LES)…
, M nombre mach √
Notre choix dans cette étude c’est le Modèle standard
Sont des constants, sont les nombres
k-ε [Launder et spalding (1972)] qui se base sur le concept
bousinesq (1977), notre choix du modèle de turbulences est de Prandtl de la turbulence pour
dépond avec les ressources informatiques et la quantité de Sont les limites définies pour l'utilisateur de
temps disponible pour la simulation. source.

Modélisation de la viscosité turbulente :

Équations de transport Model : Standard :

L'énergie cinétique de turbulence et son taux de

Constantes Modèles :
dissipation, sont obtenus à partir des équations suivantes de
transport:
, , , ,
[( ) ]

En ce qui concerne les conditions aux limites variables le long

des tubes verticaux en fonction de la profondeur du sol, nous
avons opté pour une fonction qui nous avons définie (UDF),
[( ) ] permettant d’attribuer la fonction de la température des
tubesverticaux en fonction de la profondeur. La température
de la surface du sol peut être estimée en utilisant une fonction
sinusoïdale (Eq. (1)).
Ou
(11)

Représente la génération de l'énergie cinétique de La température de la profondeur du sol peut être

turbulence due aux gradients moyens de vitesse, calculée présentée comme suit:
comme :
( ⁄ ) ( ⁄ )(12)
̅̅̅̅̅̅̅

Représente la génération de l'énergie cinétique de √ ,

turbulence due à la flottabilité, calculée comme :
Avec la température du sol au temps (t)et à une
profondeur (z), est la température moyenne annuelle du
sol, est L'amplitude de la température sur la surface du sol
Est le nombre turbulent de Prandtl pour l'énergie, (z=0), et est la fréquence de la vague annuelle de la
est le composant du vecteur de la gravité dans la direction . température.Les caractéristiques physiques ainsi que les
Pour le standard et réalizable . paramètres de calcul sont illustrés respectivement sur les
tableaux(I et II).
Le coefficient de la dilatation thermique.

4
 Le 3ème Séminaire International sur les Énergies Nouvelles et
Renouvelables
The 3nd International Seminar on New and Renewable
Energies
Unité de Recherche Appliquée en Énergies Renouvelables,
Ghardaïa – Algérie 13 et 14 Octobre 2014

demande très élevée en matière d’énergie sur tous pour les

TABLEAUI : dernières années.
Propriétés physiques [1]

Propriété Valeur
Conductivité thermique du sol (W/m. °C) La Chaleur 2.01
spécifique du sol (J/kg. °C) 1380
Densité du sol (kg/m3) 2300
Diffusivité thermique du sol (m2/s) 6.33X10-
7
Conductivité thermique du PVC (W/m. °C) 1.17
Chaleur spécifique du PVC (J/kg. °C) 900
Densité du PVC (kg/m3) 1380

TABLEAUII :

Paramètres d'entrée pour la validation comparative.

Propriété Valeur
Fig.7Image de l’échangeur de chaleur horizontal pour différentes
Longueur de pipe 50 m
profondeurs.
Diamètre de pipe 0.11m
Profondeur de pipe 4m
Vitesse d'air 4.5 ; 3.5m/s
Température d'entrée 37,48°c
Température de sol (z=4m) en été 24 ; 28°c

III. ÉTUDE EXPÉRIMENTALE :

Cette partie est développée dont le but d’étudier

l’efficacité de l’échangeur de chaleur enterré dans les régions
trop chaudes et arides.Nous avons dit que l’échangeur a été
installer au site de Biskra voir les figures (7, 8,9 réalise par le
labo LGM enceinte de l’université de Biskra), qu’est situé au
sud de l’Algérie, possède un climat aride et très chaud en été.
Biskra parmi les régions de l’Algérie de haute consommation
d’énergie pour la climatisation des locaux, surtout avec les
périodes de température trop élevées qui peuvent atteindre Fig.8Emplacement des sondes tout au longue de l’échangeur de chaleur.
50°c. C’est pour cela le laboratoire de génie mécanique de
l’université de Biskra a lancé des recherches pour développer
un système de climatisation de basse consommation
d’énergie, afin de réduire la consommation d’énergie utilisée
pour la climatisation. Il est important de signaler qu’une

5
 Le 3ème Séminaire International sur les Énergies Nouvelles et
Renouvelables
The 3nd International Seminar on New and Renewable
Energies
Unité de Recherche Appliquée en Énergies Renouvelables,
Ghardaïa – Algérie 13 et 14 Octobre 2014

50
45

403

30

25
0 500 1000 1500 2000 2500 3000
Time of day in minute
La température de l’air d’entrée
La température de l’air de sortie

Fig.10 Comparaison entre la température d’entrée et de sortie de

Fig.9Emplacement détaillé des sondes de l’air à l’intérieur de la l’échangeur.
gaine.
IV. VALIDATION DU MODÈLE DE
SIMULATION:
Nous avons fait des essaies durant les mois de Juillet et
Aout, les deux mois les plus chaud de l’année, pour étudier Pour valider notre travail nous avons fait une étude
réellement le comportement de l’échangeur et son efficacité, expérimentale avec des sondes qui donnent la variation de la
la figure (10) présente l’écart de la température entre l’entrée température tout au long des tubes de l’échangeur de chaleur
et la sortie de l’échangeur de chaleur enterrée pour trois jours et qui nous permettons de faire la comparaison avec les
successifs pour les deux mois. données numériques.

A. Les iso-surfaces:
50
Les iso-surfaces sous FLUENTreprésente la méthode
45 de choisir les points de température tout au long des pipes, la
figure (8) présente la température de l’air (ligne rouge) et la
température de la paroi.
40

35

30

25
0 5 10 15 20 25 30 35 40 45
2
Time of day in minute *10 Fig.11Les iso-surfaces

6
 Le 3ème Séminaire International sur les Énergies Nouvelles et
Renouvelables
The 3nd International Seminar on New and Renewable
Energies
Unité de Recherche Appliquée en Énergies Renouvelables,
Ghardaïa – Algérie 13 et 14 Octobre 2014

TABLEAUIII:
Validation du modèle de simulation (température de l'air ambiant:
37°C).

Distance axiale de Données Résultat

l’admission expériment s
d'EATHE (m) ales (°C) numériq
ue (°C)
2.65 35.24 35,66
6.9 32.28 32,39
10.81 29.83 29,35
18 26.74 26,98
23.44 25.95 25,84 0 5 10 15 20 25 30 35 40 45
30.21 25.00 25,02 Fig.13Distribution de la température de l’air avec la longueur des
45.3 24.21 24,26 conduites.

48
46 35 40 45
44
42
40
38
36
34
32
30
28

0 5 10 15 20 25 30 35 40 45
Longueur de la conduite (m)
Longueur de la conduite (m)
La température simulée après 1 heure
La température simulée avec un maillage Hex/wedge de 1.5.
La température expérimentale après 1 heure La température simulée avec un maillage Tet/Hybrid de 1.5.
La température simulée après 4 heures La température simulée avec un maillage Hex/wedge de 1.
La température expérimentale après 4heures La température expérimentale.
35 40 45
Fig.14La variation de la température pour différent types de maillages.
Fig.12 Validation des résultats de CFD avec des résultats
expérimentaux.

7
 Le 3ème Séminaire International sur les Énergies Nouvelles et
Renouvelables
The 3nd International Seminar on New and Renewable
Energies
Unité de Recherche Appliquée en Énergies Renouvelables,
Ghardaïa – Algérie 13 et 14 Octobre 2014

Référence :

[1] Sehli, A., et al. (2012). "The Potential of Earth-

air Heat Exchangers for Low Energy Cooling of Buildings
in South Algeria." Energy Procedia18: 496-506.

[2] Misra, R., et al. (2013). "CFD analysis based

parametric study of derating factor for Earth Air Tunnel
Heat Exchanger." Applied Energy103: 266-277.

[3] Esen, H., et al. (2007). "Numerical and

experimental analysis of a horizontal ground-coupled heat
pump system." Building and Environment42(3): 1126-
1134.

Fig.15 Distribution de la température sur [4] Bansal, V., et al. (2010). "Performance analysis
la paroi des conduites. of earth–pipe–air heat exchanger for summer cooling."
Energy and Buildings42(5): 645-648.
V. CONCLUSION:
[5] Bansal, V., et al. (2012). "Performance analysis
Dans cet article n on peut conclure que l’utilisation of integrated earth–air-tunnel-evaporative cooling system in
hot and dry climate." Energy and Buildings47: 525-532.
de maillage tridimensionnel est efficace pour étudier la
variation de la température tout au long de l’échangeur. On [6] Khalajzadeh, V., et al. (2011). "Parameters
remarque quel’étude numérique est simple et efficace, optimization of a vertical ground heat exchanger based on
quidonne des résultats acceptables par rapport aux mesures response surface methodology." Energy and
expérimentales. Soit pour la conduction sur la paroi de Buildings43(6): 1288-1294.
l’échangeur ou pour la circulation de l’air au long de
[7] N. Moummi et al (2010). "Le rafraîchissement
l’échangeur, une concordance est remarquée. par la géothermie: étude théorique et expérimentale dans le
site de Biskra." Energies Renouvelables Vol. 13 N°3 399 –
Pour l’étude numérique, nous pouvons conclu que la 406.
géométrie symétrique de l’échangeur donne de bons
résultats avec un temps de simulation court par rapport à [8] P.M. Congedo et al (2012). "CFD simulations of
l’utilisation de géométrie complète. horizontal ground heat exchangers: A compariso namong

Nous avons opté pour un maillage suffisamment fin different configurations."Applied Thermal Engineering 33-
34 24e 32.
de sorte que les résultats des variables recherchés ne sont
pas influencés par un raffinement poussé.

Pour l’étude expérimentale, selon les résultats

obtenus en été, on a remarqué un écart très important de
température entre l’entrée et la sortie de l’échangeur qui
peut aller jusqu'à 20°C. Ces résultats donnent une preuve
que l’échangeur de chaleur enterré peut contribuer
fondement à l’économie de l’énergie aussi bien en été
qu’en hiver.

8
 Le 3ème Séminaire International sur les Énergies Nouvelles et
Renouvelables
The 3nd International Seminar on New and Renewable
Energies
Unité de Recherche Appliquée en Énergies Renouvelables,
Ghardaïa – Algérie 13 et 14 Octobre 2014

Nomenclature :

EATHE : échangeur de chaleur enterré Air-Sol.

: X vitesse, m/s.

: Y vitesse, m/s.

Z vitesse, m/s.

viscosité turbulente m2/s.

viscosité cinématique m2/s.

Masse volumique de l’air, kg / m3

Température de l’air, °C.

P : Pression de l’air.

Viscosité dynamique de l’air.

Taux de dissipation de turbulence.

L'énergie cinétique de turbulence.

Viscosité cinématique de l’air.

Diffusivité thermique moyenne.

m2 /s.

La fréquence de la vague annuelle de la

température

Profondeur des pipes, m.

Conductivité thermique du sol

: Capacité de chaleur spécifique du sol.