Ministère de l’Enseignement Supérieur et de la Recherche Scientifique
École Nationale Polytechnique de Constantine, Malek BENNABI
Présentation
SIMULATIONS NUMERIQUE D’UN ECHANGEUR DE CHALEUR EN
FORAGE PROFOND « BOREHOLE » DANS UN SYSTEME
GEOTHERMIQUE
Projets de fin d’étude d’ingénieur
Présenté par : Encadré par :
Membres de Jury
ENP Constantine Présidant
ENP Constantine Examinateur
ENP Constantine Rapporteur
En vue l’obtention de diplôme d’ingénieur d’état en Génie mécanique
Soutenue publiquement à Constantine le 13/09/2021
Plan de la présentation
1. Introduction
2. Étude théorique simplifiée
3. Étude de la géométrie et du maillage
4. Résultats
5. Conclusion
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 2
1. Introduction
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
La Géothermie
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 4
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Les types de la géothermie
Haute Moyenne Très basse
Haute énergie Basse énergie
énergie énergie énergie 10-
>150°C
>150°C
30-90°C
90-150°C 30°C
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 5
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
L’exploitation de la géothermie
Au monde
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 6
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
L’exploitation de la géothermie
En Algérie
Figure. 1 : Répartition de la capacité installée de géothermie à usage direct par pays en Afrique.
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 7
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Le forage
Figure. 2 : Schéma simplifié d'une installation de forage fond de trou (MFT) .
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 8
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Les pompes à chaleur géothermique
Figure. 3 : Les différents types de la géothermie superficielle .
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 9
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Problématique?
Quel est l’influence de la géométrie d’un échangeur de chaleur (forage)
sur la température de sortie?
Quel est l’influence des paramètres ci-dessous sur la température
de sortie d’un échangeur de chaleur au sol?
Influence de diamètre de l’hélicoïdale
Influence du débit de l’eau à l’entrée de tube hélicoïdale
Influence du pas « e » de l’hélicoïde
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 10
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Objectifs?
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 11
2. Étude théorique simplifié
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Formulation mathématique du problème
1.Équation de continuité : (1)
2.Équation de quantités de mouvement : (2)
3.Équation de l’énergie (3)
3.Taux d’extraction de chaleur : (4)
4.Capacité de performance thermique : (5)
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 13
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Modélisation de la turbulence
Modèle k –epsilon
L’équation de l’énergie cinétique turbulente k :
(6)
L’équation de taux de dissipation de l’énergie cinétique ɛ :
(7)
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 14
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Méthode numérique
Méthode des volumes finis
Figure. 4 : Volume de contrôle élémentaire. Figure. 5 : Subdivision du volume de contrôle.
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 15
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Méthode numérique
Les schémas de discrétisation
Variables Schémas de
discrétisation
Pression Second Order
Quantité de mouvement Second Order
Energie cinétique de turbulence, k Second order Upwind
Le taux de dissipation de l’énergie cinétique de turbulence, ɛ Second order upwind
Energie Second order upwind
Tableau.1:Les schémas de discrétisations.
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 16
3. Étude de la géométrie et du
maillage
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Géométrie et description du problème
Figure. 6 : Les dimensions du Borehole « tube en U ». Figure. 7 : Dimensionnement de tube hélicoïde .
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 18
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Choix et optimisation du maillage
Figure. 8 : Vus du maillage sur la surface supérieure de tube en U . Figure. 9 :Vus du maillage d’une coupe vertical de tube en U.
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 19
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Organigramme et procédure générale d’utilisation du logiciel Gambit
et ANSYS Fluent
Gambit ( création de Gambit (conditions aux
Gambit (maillage) Ansys Fluent
géométrie) limites)
Condition sur les Mise à jour des
Création des points Maillage des arrêtes
faces propriétés
Condition sur les
Création des arrêtes Maillage des faces Initialiser la solution
volumes
Création des Régler la solution du
Maillage es volumes
surfaces « monitors »
Création des Export du ficher
Lancer le calcul
volumes mesh
Convergence de
solution
Analyse et
visualisation
Figure. 10 : Organigramme de calcul.
Soutenance mémoire d'ingénieur, Constantine 13– 09 – 2021 20
4. Résultats
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Résultats
• Validation
• Influence de la géométrie de l’échangeur de chaleur (forage)
• Influence de diamètre de l’hélicoïde
• Influence du débit de l’eau à l’entrée de l’hélicoïde
• Influence du pas « e » de l’hélicoïde -100 kV
30 kV 35 kV
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 22
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Résultats
• Validation
• Influence de la géométrie de l’échangeur de chaleur (forage)
• Influence de diamètre de l’hélicoïde
• Influence du débit de l’eau à l’entrée de l’hélicoïde
• Influence du pas « e » de l’hélicoïde -100 kV
30 kV 35 kV
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 23
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Validation
316
314
312
310
Température (K)
308
25 kV 25 kV 306
304
302 Presente étude
300 Lyu et al [5]
298 -100 kV 45 kV
0.2 0.6 1
Débit massique (kg/s)
30 kV 35 kV
Figure. 11 : Comparaison des contours de la température de sortie des différents débits entre le travail présent et des résultats
numérique de Lyu et al.
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 24
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Résultats
• Validation
• Influence de la géométrie de l’échangeur de chaleur (forage)
• Influence de diamètre de l’hélicoïde
• Influence du débit de l’eau à l’entrée de l’hélicoïde
• Influence du pas « e » de l’hélicoïde -100 kV
30 kV 35 kV
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 25
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Influence de la géométrie de l’échangeur de chaleur
Figure. 13 : Contours de température de tube en U et l’hélicoïde.
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 26
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Influence de la géométrie de l’échangeur de chaleur
Figure. 14 : Contour de vitesse de tube en U et l’hélicoïde. Figure. 15 : contour de la pression de l’hélicoïde et tube en U.
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 27
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Résultats
• Validation
• Influence de la géométrie de l’échangeur de chaleur (forage)
• Influence de diamètre de l’hélicoïde
• Influence du débit de l’eau à l’entrée de l’hélicoïde
• Influence du pas « e » de l’hélicoïde -100 kV
30 kV 35 kV
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 28
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Influence de diamètre de l’hélicoïde
Figure. 16 : Contours de la température des différents diamètres Figure. 17 : Contours de la vitesse magnitude des différents
de l’hélice. diamètres de l’hélice.
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 29
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Influence de diamètre de l’hélicoïde
75 316.5
316
Taux d'extraction de chaleur (kW)
70
Température de sortie (K)
315.5
315
65
314.5
60
314
313.5
55 Taux d'extraction de chaleur
313
Temperature de sortie
50 312.5
0.1 0.15 0.2
Diamètre de l'hélice (m)
Figure. 18 : Taux d’extraction et la température de sortie en fonction des différents diamètres.
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 30
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Résultats
• Validation
• Influence de la géométrie de l’échangeur de chaleur (forage)
• Influence de diamètre de l’hélicoïde
• Influence du débit de l’eau à l’entrée de l’hélicoïde
• Influence du pas « e » de l’hélicoïde -100 kV
30 kV 35 kV
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 31
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Influence du débit pour un diamètre de l’hélice de 0.2 m
Figure. 19 :Contours de température statique pour les différents Figure. 20 :Contour de la vitesse magnitude pour les différents
débits. débits.
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 32
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Influence du débit pour un diamètre de l’hélice de 0.2 m
80 345
70 340
Taux d'extraction de chaleur (kW)
335
60
Température de sortie (K)
330
50
325
40
320
30
315
20
310
Taux d'extraction de chaleur Temperature de sortie
10 305
0 300
0.2 0.6 1
Débit massique (kg/s)
Figure. 21 : Le taux d’extraction et la température de sortie en fonction de débit massique .
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 33
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Résultats
• Validation
• Influence de la géométrie de l’échangeur de chaleur (forage)
• Influence de diamètre de l’hélicoïde
• Influence du débit de l’eau à l’entrée de l’hélicoïde
• Influence du pas « e » de l’hélicoïde -100 kV
30 kV 35 kV
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 34
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Influence du pas « e » pour un diamètre de l’hélice de 0.2 m et un débit 0.2Kg/s
Figure. 22 : Contours de la température statique pour les différents Figure. 23 : Contour de la vitesse magnitude pour les différents
pas « e ». pas « e ».
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 35
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Influence du pas « e » pour un diamètre de l’hélice de 0.2 m et un débit 0.2Kg
34.5 340
34
339
33.5
Taux d’extraction de chaleur (kW)
Température de sortie (K)
33
338
32.5
32 337
31.5
336
31
30.5
335
30 Taux d'extraction de chaleur Temperature de sortie
29.5 334
0.1 0.15 0.2
Pas de l'hélice (m)
Figure. 24 : La variation de taux d’extraction et la température de sortie en fonction du pas « e ».
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 36
5. Conclusion générale
Introduction Étude théorique simplifiée Étude de la géométrie et du maillage Résultats Conclusion
Conclusion générale
L’échangeur de chaleur du tube spiralé VS L’échangeur de chaleur du tube
en U.
La diminution du diamètre de l’hélicoïde une augmentation de
température de l’eau à la sortie.
La diminution du débit de l’eau à l’entrée du tube hélicoïdal une
augmentation de température à la sortie.
-100 kV
Les meilleures performances concernant la température à la sortie sont
30 kV
obtenues pour le cas du pas le plus petit e = 0.1m.
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 38
Fin de
présentation
d’Ingéniorat
‹#›
Projet de Master
SIMULATION NUMERIQUE D'UN ECHANGEUR DE CHALEUR
EN FORAGE PROFOND « BOREHOLE » DANS UN SYSTEME
GEOTHERMIQUE : APPLICATION DANS LE CONTEXTE
ALGERIEN
En vue l’obtention de diplôme de Master de Master en Energétique
Soutenue publiquement à Constantine le 13/09/2021
Introduction Technique expérimentale Résultats Conclusion
Plan de la présentation
1. Introduction
2. Étude de la géométrie
3. Résultats
4. Conclusion
Soutenance mémoire de Master, Constantine 13 – 09 – 2021 41
1. Introduction
Introduction Etude de la géométrie Résultats Conclusion
Donner sur la géothermie en Algérie
Figure. 1 : Carte de température des ressources géothermales.
Soutenance mémoire de Master, Constantine 13 – 09 – 2021 43
Introduction Etude de la géométrie Résultats Conclusion
Etude intérieur
Figure. 2 : Source thermale de Hammam Meskhoutine, Guelma, Algérie .
Soutenance mémoire de Master, Constantine 13 – 09 – 2021 44
Introduction Etude de la géométrie Résultats Conclusion
Quelques réservoirs géothermique de l’Algérie
Figure. 3 : Les différents reliefs rencontrés dans la région de Figure. 4 : Extension géographique et principaux points d’eau
Hammam Meskhoutine. captant la nappe albienne en Algérie.
Soutenance mémoire de Master, Constantine 13 – 09 – 2021 45
Introduction Etude de la géométrie Résultats Conclusion
Les applications possible de la géothermie en Algérie
Figure. 5 : Schéma d’un chauffage des habitats. Figure. 6 : Schéma d'une serre agricole.
Soutenance mémoire de Master, Constantine 13 – 09 – 2021 46
Introduction Etude de la géométrie Résultats Conclusion
Problématique?
Quel est l’influence des ailettes dans un échangeur de chaleur
(forage) sur la température de sortie de l’eau de ce dernier ?
Quel est l’influence du débit de l’eau à l’entrée de tube en U avec ailettes sur la
température de sortie de l’eau d’un échangeur de chaleur au sol?
Soutenance mémoire d'ingénieur, Constantine 13 – 09 – 2021 47
Introduction Étude théorique simplifiée Dispositifs expérimentaux Résultats Conclusion
Objectifs?
Soutenance mémoire de Master, Constantine 13– 09 – 2021 48
2. Etude de la géométrie
Introduction Etude de la géométrie Résultats Conclusion
Géométrie et description du problème
Figure. 7 : Les dimensions du Borehole « tube en U ».. Figure. 8 : Les dimensions de « tube en U avec ailettes» .
Soutenance mémoire de Master, Constantine 13 – 09 – 2021 50
3. Résultats
Introduction Etude de la géométrie Résultats Conclusion
Résultats
• Influence des ailettes dans un tube en U
• Influence de débit dans un tube en U avec ailette
Soutenance mémoire de Master, Constantine 13 – 09 – 2021 52
Introduction Etude de la géométrie Résultats Conclusion
Résultats
• Influence des ailettes dans un tube en U
• Influence de débit dans un tube en U avec ailette
Soutenance mémoire de Master, Constantine 13 – 09 – 2021 53
Introduction Etude de la géométrie Résultats Conclusion
• Influence des ailettes dans un tube en U
Figure. 11: Contours de la température statique de tube en U Figure. 12: Contours de la température statique de tube
lisse. avec ailettes.
Soutenance mémoire de Master, Constantine 13 – 09 – 2021 54
Introduction Etude de la géométrie Résultats Conclusion
• Influence des ailettes dans un tube en U
360
350
Température de sortie de l'eau (K)
340
330
320
310
Température de sortie de U sans ailettes
300
Température de sortie de U avec ailettes
290
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1
Débit massique (kg/s)
Figure. 13: La comparaison des résultats numérique de température de sortie en fonction du débit des deux tubes en U
lisse et avec ailettes.
Soutenance mémoire de Master, Constantine 13 – 09 – 2021 55
Introduction Etude de la géométrie Résultats Conclusion
Résultats
• Influence des ailettes dans un tube en U
• Influence de débit dans un tube en U avec ailette
Soutenance mémoire de Master, Constantine 13 – 09 – 2021 56
Introduction Etude de la géométrie Résultats Conclusion
• Influence de débit dans un tube en U avec ailette
Figure. 14: Contours de température des différents débits.
Soutenance mémoire de Master, Constantine 13 – 09 – 2021 57
Introduction Etude de la géométrie Résultats Conclusion
• Influence de débit dans un tube en U avec ailette
120 355
350
Taux d'extraction de chaleur (kW)
100
345
80 340
Température (K)
335
60
330
40 325
320
20 Taux d'extraction de chaleur
315
Température de sortie
0 310
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1
Débit massique (kg/s)
Figure. 15: Le taux d’extraction et la température de sortie en fonction du débit.
Soutenance mémoire de Master, Constantine 13 – 09 – 2021 58
4. Conclusion générale
Introduction Etude de la géométrie Résultats Conclusion
Conclusion
L’ajoutant des ailettes la surface d’échange de chaleur augmente ainsi que la
température de sortie.
La comparaison entre un tube lisse et un autre avec ailettes se fait avec la même vitesse
pour obtenir des résultats corrects.
La diminution du débit de l’eau à l’entrée du tube en U avec ailettes une
augmentation de température à la sortie.
Soutenance mémoire de Master, Constantine 13 – 09 – 2021 60
Fin de
présentation de
Master
‹#›
Fin de la
présentation
Merci Pour Votre attention
‹#›