1st National Conference on Applied Energetics.

NCAE 2020, 11 and 12 February 2020, NAAMA, Algeria.

ÉTUDE NUMÉRIQUE D’UN ÉCHANGEUR DE CHALEUR A PLAQUES ET

JOINT
K. Boukhadia1, H. Ameur2
1
K. Boukhadia, Centre Universitaire Salhi Amed CUN, Naama, , Algérie e-mail boukhadiakarima@yahoo.fr
2
H. Ameur, Centre Universitaire Salhi Amed CUN, Naama, , Algérie e-mail boukhadiakarima@yahoo.fr

RÉSUMÉ

Les échangeurs de chaleur à plaques permettent d’atteindre des compacités encore plus élevées que les
échangeurs multitubulaires. La conception des plaques optimise le transfert de chaleur en offrant une surface
totale large et compacte à travers laquelle la chaleur peut être transmise d'un fluide à un autre. On cherche à
améliorer le transfert de chaleur entre ces plaques tout en cherchant à diminuer toutefois la perte de charge au
cours de l’écoulement du fluide dans l’appareil. La simulation numérique permet de calculer le flux thermique
du fluide dans l'échangeur de chaleur ; la perte de charge et la performance de l’échangeur de chaleur tout en
variant les vitesses de l’écoulement ainsi que l’espace entre les plaques. Les résultats ont mené à dire qu’une
augmentation du flux thermique échangé entre les deux plaques pour différentes épaisseurs du joint est constaté
lorsque les plaques sont serrées, elles permettent un échange thermique plus large qu’en augmentant l’espace
inter plaques cette augmentation est accompagnée avec une augmentation de chute de pression.
Key Words: Heat Transfer, Finite Elements, Natural Convection. Échangeur de chaleur à plaques et joint,
épaisseur du joint, vitesse de l’écoulement entre plaque, flux thermique et chute de pression.

NOMENCLATURE

Symboles Nombres sans dimensions

A ou S surface d’échange, m² J facteur de friction
Cp : Chaleur spécifique à pression constante J/kg K f : facteur de friction
Dh Diamètre hydraulique moyen m Re Nonbre de Reynolds
ΔTLM Différence de température logarithmique moyenne °C Nu Nombre de Nusselt
ΔP Chute de pression Pa Pr Nombre de Prandt
E Efficacité Lettres grecs
f Coefficient de frottement Φ Puissance thermique kW
h Coefficient de transfert de chaleur W/ m2 K ρ density, kgm-3
L Longueur de la plaque m λ conductivité thermique, W·m-1·K-1
m.Cp Débit capacité calorifique du fluide kJ/h °C μ Viscosité dynamique kg /m /s
NUT Nombre dunité de transfert Indices
P Pression, Nm-2 1 temperature à l’éntre de l’échangeur
T Température, K 2 temperature à la sortie de l’échangeur
u, v w : Compostantes de la vitesse suivant les trois axes m/s
léchangeurnternal
 1. INTRODUCTION

Aujourd’hui avec une croissance permanente des prix de l’énergie, la maitrise de l’énergie est devenue un
enjeu majeur dans tous les domaines d’activités. Pour les professionnels de l’énergie, le premier enjeu est de
concevoir des systèmes et procédés énergétiques avec des meilleures efficacités.
Le transfert de chaleur représente l’un des modes les plus communs d’échange d’énergie. C’est un phénomène
que l’on trouve dans nombreux secteurs de l’industrie. Les ingénieurs et les techniciens se trouvent confronté à
ce genre de problème et essayent de maximiser ou de minimiser ce phénomène selon les besoins de l’industrie
et dans le souci d’économiser cette énergie qui revient chère. De ce fait, les transferts thermiques ont, aussi
bien dans le domaine des sciences pures que dans celui des applications technologiques, un rôle souvent
essentiel. Exemple : (échangeur, moteurs thermiques, calorifugeage, isolation thermique ….etc.). La
connaissance des lois physiques qui régissent ces modes de transfert thermique est une chose essentielle et très
importante, car elles nous permettent de maitriser la façon et la qualité de cet écoulement de chaleur.
Les échangeurs thermiques sont en effet le point de passage obligé de tous les flux d’énergie, et tout effort de
maîtrise de cette énergie doit s’appuyer sur des échangeurs de chaleur bien conçus et bien utilisés.
L'échangeur de chaleur est un appareil thermique de grande importance dans les installations thermiques et
énergétique. On rencontre au moins un échangeur de chaleur dans une installation thermique. Ceux-ci incluent
la production d'énergie ; les industries alimentaires, de produit chimique; dans le domaine d’électronique ;
technologie environnementale; industrie ; climatisation et réfrigération [1].
A la différence des autres appareils thermiques, l’échangeur de chaleur ne contient aucune pièce mécanique
mobile. Le calcul de cet appareil est très complexe, on doit connaître exactement : sa géométrie (surface
d’échange et section de passage des fluides), ses caractéristiques thermophysiques, les vitesses d’écoulements
des fluides, les températures d’entrées des fluides, etc

2.MODÈLE MATHÉMATIQUE

Lors des études générales relatives à des installations des échangeurs de chaleur apparait souvent la nécessité
d’évaluer [2]:
 Les coefficients d’échange thermique.
 Les surfaces d’échanges.
 Les pertes de pressions.
 Les écarts de températures entre les fluides caloporteurs.

2.1. Pertes de charges

En général on considère que les dimensions des tubulures sont suffisamment importantes pour que les pertes de
charges y soient négligeables devant celles des plaques. La perte de charge d’un fluide en écoulement
isotherme dans une plaque (ΔP) entre les collecteurs d’entrée et de sortie, dépend de nombreux paramètres tels
que la densité, la viscosité, la vitesse du fluide, ainsi que la géométrie des plaques.

ρU2 L
P = Pcollecteur + 4. f. 2 Dh
(1)

Le coefficient de frottement (f) s’exprime d’une manière générale comme suite

𝑘
𝑓=𝑅 (2)
e𝑥

K et x dépendent du régime d’écoulement et de la géométrie de la plaque. L’exposant x du nombre de

Reynolds varie de 0,1 à 0, 4 [3].

2.2. Nombre d’unité de transfert NUT

Le groupement sans dimension h.S/ Cmin représente ce qu’on appelle le nombre d’unité de transfert noté NUT.
H.A 1
NUT = (m.C) =R (3)
min eq .(m.C)min

2
2.3. Efficacité d’un échangeur de chaleur
Les relations suivantes donnent la relation reliant l’efficacité d’échange au nombre d’unité de transfert pour
deux modes de fonctionnement différents des échangeurs pour un fonctionnement à contre-courant:

1−e[−(1−Rf)NUT]
E= (4)
1−Rf e[−(1−Rf)NUT]

2.4. Détermination du coefficient thermique de convection h [4].

Le problème de la convection est en fait de déterminer en fonction des conditions d’écoulement du fluide, des
caractéristiques géométriques des parois et des éventuels changements d’état du fluide. On définit donc trois
nombres sans dimension :
ℎ.𝐷ℎ
 Le nombre de Nusselt donné sous la forme : 𝑁𝑈 = 
ρ.U.Dh
 le nombre de Reynolds donné sous la forme :R e =

 Cp
 le nombre de Prandt donné sous la forme :Pr =
λ

2.5. Flux thermiqu échangé entre les plaques

Dans l’étude thermique des différents types d’échangeurs de chaleur, on a souvent recours à l’équation de type
Le flux thermiqu échangé entre les plaques s’exprime comme [5].
 = h. A(T1 − T2 ) (5)
𝜙 = ℎ. 𝐴. Δ𝑇𝐿𝑀 (6)

3. RESULTS

La présente étude porte sur l’étude du transfert thermique à travers des plaques d’un échangeur de chaleur à
plaques et joint. L’écoulement est supposé laminaire tridimensionnel, compressible (eau), se refroidi entre les
plaque en changeant les épaisseurs (distance entre les plaques) suivant le tableau en dessous.

FIGURE 1. Schéma détaillé de la géométrie de la plaque de l’échangeur de chaleur.

3
3.1. Structure de l’écoulement entre deux plaques adjacentes d’un échangeur de chaleur
La conception des plaques optimise le transfert de chaleur en offrant une surface totale large mais compacte à
travers laquelle la chaleur peut être extraite d'un liquide ou d'un gaz à un autre.

La zone de transfert de chaleur des plaques est comprise entre deux plaques superposées avec des motifs en
chevrons opposés, l'écoulement en forme d'hélice et une turbulence élevée se traduisent par des coefficients de
transfert élevés et un auto-nettoyage efficace.

La zone de distribution de la plaque assure un écoulement uniforme du fluide sur toute la plaque afin de
maximiser la capacité de transfert de chaleur. Une répartition optimisée du débit réduit également
l’encrassement et les zones de température inégales, tout en maintenant les niveaux de performance à un niveau
élevé, sans pertes d’énergie inutiles, coûts de maintenance ou arrêts imprévus.

Les joints sont des composants clés dans la performance des échangeurs de chaleur. Le joint et la plaque sont
conçus de manière à assurer une étanchéité optimale. Chacun est adapté à la mission de l'échangeur de chaleur.
Le profil, la largeur, l'épaisseur, le type de polymère et le composé appropriés font toute la différence pour
éviter les risques de fuite prématurée ou de dommage aux joints ou aux plaques.

FIGURE 2. Structure de l’écoulement sur une plaque simple d’un échangeur de chaleur.

3.2. Conditions initiales

Le logiciel CFX 14.0 met à notre disposition un sous-programme dans lequel il est possible d’introduire une
distribution initiale pour chacune des variables de l’écoulement.

 Le fluide est de l’eau.

 La densité : ρ = 999.7 [kg.m-3].
 La viscosité dynamique : μ = 8.899*10-4[kg /m /s].
 La pression de référence : 0 [atm].
 Les caractéristiques physiques et de transports de l’eau sont définies par le code de calcul.
3.3. Conditions aux limites
Les conditions aux limites sont définies et doivent être appliquées à toutes les régions aux extrémités du
domaine de calcul à savoir :
 L’entrée du fluide (INLET) : La vitesse à l’entrée de la plaque est variée de 0.0006, 0.001,0.01, 0.05 et
0.1[m/s]

4
 
0
 La sortie du fluide (OUTLET) : x ; p = 0 [atm]
 Parois : Une condition de non adhérence est considérée pour la plaque avant et arrière ou le coefficient de
transfert thermique global est défini comme h=379.44 [W/m2 K]
u v  w0
 Paroi adiabatique : Les surfaces du joint sont définies comme paroi adiabatiques
3.4. Influence de la vitesse d’écoulement

FIGURE 3. Variation du profil de température entre deux plaques de 3 mm d’épaisseur pour

différentes vitesses.

L’évolution du profil de vitesse de l’écoulement entre deux plaques de 3 mm d’épaisseur varie de manière
significative en fonction de la vitesse en m/s. Cette évolution fait appel à un gradient de température variant
entre 343 K à l’entré jusqu’à une valeur de 326 K à la sortie ce qui permet un bon refroidissement du fluide
entre les plaques.

Cette constatation est traduite par une augmentation du flux thermique échangé entre les deux plaques pour
différentes épaisseurs du joint, lorsque les plaques sont serrées, elles permettent un échange thermique plus
large qu’en augmentant l’espace inter plaques.

350

300

250
Q [kg mm^2 s^-3 ]

200 e=3mm
e=4mm
150 e=6mm
e=8mm
e=10mm
100

50

0
0 500 1000 1500 2000 2500

Re

FIGURE 4. Variation du flux thermique en fonction de Reynolds pour différentes épaisseurs du joint.

L’échange thermique entre les plaques serrées est clairement enregistré sur les différents profils de
températures en variant l’épaisseur du joint pour la même vitesse. Le gradient de température varie de 343 K à
l’entrée de la plaque jusqu’à 323 K à la sortie pour une épaisseur de 3 mm du joint quant la différence est de
343 K à l’entré à 326 K à la sortie pour des épaisseurs variant de 6 à 10 mm.
5
 FIGURE 5. Variation du flux thermique en fonction de Re pour différentes épaisseurs du joint.

3.5. Influence de l’épaisseur du joint

La figure 6 montre que la chute de pression pour une large épaisseur de joint prend des valeurs maximales et
elle est légèrement inférieure pour des valeurs moindres de l’espace entre les plaques pour la même vitesse
d’écoulement. Ceci est principalement dû au cisaillement de l’écoulement entre les surfaces des plaques

FIGURE 6. Variation du profil de pression entre deux plaques de différentes épaisseurs du joint pour
une vitesse de 0.05 m/s.

Habituellement, l'augmentation du transfert de chaleur entre les plaques d’un échangeur de chaleur génère des
pertes de pression. L’évolution de la perte de charge le long de la plaque est présentée sur la figure 7.

60 e=3mm
e=4mm
Pressure Drop [ kg mm^-1 s^-2 ]

e=6mm
50
e=8mm
e=10mm
40

30

20

10

0
0 500 1000 1500 2000 2500

Re

FIGURE 7 . Variation de la chute de pression en fonction de Re pour différentes épaisseurs du joint.

6
La variation de la chute de pression tend à augmenter avec l’augmentation du nombre de Reynolds. A chaque
fois l’épaisseur de joint diminue, la valeur de chute de pression augmente.

3.6. Performance de l’échangeur de chaleur

Dans les terminologies de surfaces compactes, le coefficient de transfert de chaleur et la perte de charge sont
généralement exprimés en termes de facteurs sans dimension. Deux paramètres adimensionnels bien connus, le
facteur de Colburn, j, et le facteur de friction Fanning, f, ont été utilisés pour décrire les performances de
transfert de chaleur et les caractéristiques de perte de charge,
ℎ 2⁄
𝑗 = 𝜌𝑢𝐶 𝑃𝑟 3 (7)
𝑝

2∆𝑃 𝐷ℎ
𝑓 = 𝜌𝑢2 𝐿
(8)

où Dh= 2*e et e représente l’épaisseur du joint).

La variation de performance thermique pour le fluide circulant dans un échangeur de chaleur à plaques et joint
est montrée sur la figure 8. Sur cette figure, la performance de l’échangeur tend à augmenter en diminuant
l’épaisseur du joint pour différentes vitesse correspondant aux différents nombres de Re.

1,00E-012

e=3mm
8,00E-013
e=4mm
e=6mm
e=8mm
6,00E-013 e=10mm
Perf

4,00E-013

2,00E-013

0,00E+000
0 200 400 600 800

Re

FIGURE 8. Variation de la performance de l’échangeur de chaleur en fonction du nombre de Re pour

différentes épaisseurs.

4. CONCLUSIONS

L’échangeur thermique se définit en fonction du rôle joué (réchauffement, refroidissement) de la nature des
fluides en présence (vapeur ou liquide, risques de corrosion) et débits nécessaires.
On cherche à optimiser les frais d’investissement et de fonctionnement. Il faut donc étudier ensemble les
aspects purement thermiques (coefficient de transfert, surface d’échange) et les aspects liés à la mécanique des
fluides.
Dans cette étude, à l’inverse de la méthode classique de calcul et d’optimisation des échangeurs de chaleur
nous avons présenté des résultats obtenus grâce au logiciel CFX 14.0. Ces résultats interprètent d’une manière
simple les diverses variations que subit le fluide véhiculé à l’intérieur des plaques ainsi que son échange
thermique avec le deuxième fluide aux différentes vitesses et différentes épaisseurs du joint.
En effet les flux thermiques sur les plaques d’un échangeur de chaleur à plaque et joint sont d’autant plus forts
que les vitesses d’écoulement sont élevées. Donc si les fluides ont des débits élevés, il est possible de se

7
satisfaire d’une surface d’échange plus faible ce qui diminue l’investissement. Par contre les pertes de charge
seront plus élevées ce qui est nécessitera des pompes plus puissantes (augmentation du cout d’investissement)
et une dépense énergétique plus élevée (augmentation du cout de fonctionnement).
L’étude d’optimisation au cas par cas fournira la réponse pour déterminer le choix le plus économique.

RÉFÉRENCES

[1] Sacadura, J.-F., G. Charnay, and J.-p. Schon, Initiation aux transferts Thermiques; Technique et
Documentation. Centre d'Actualisation Scientifique et Technique, INSA de Lyon, France, 2012.
[2] Vidil, R., Les échangeurs à plaques et joints: description, éléments de dimensionnement 2014: Groupement
pour la recherche sur les échangeurs thermiques.
[3] Kreith, F. and B.a.-D. Huga, Transmission de la chaleur et thermodynamique:" Principles of heat transfer".
Traduction et adaptation par Kodja, Badr-el-Dine 2017: Masson et Cie.
[4] Benayad, N., Etude et optimisation des échanges de chaleur,Université Mentouri Constantine, 2008.
[5] Djossou, A.A., et al., Etude des variations des proprietes thermophysiques et de la durabilite, en fonction
de l’humidite des eco-materiaux de construction au Benin, 2013, EPAC/UAC.
[6] Bordjane, M., Modes de transfert thermiquen,Université des Sciences et de la Technologie d’Oran
«Mohamed Boudiaf» 2017.
[7] El Haroui, M., Modelisation et analyse de la convection thermosolutale dans un milieu poreux sature de
fluide binaire,l’Université Sidi Mohamed Ben Abdellah 2017.
[8] Nehal, S. and O. Meziani Ourrad, Contribution à l’étude des transferts thermiques dans les échangeurs,
2017, Université Abderrahmane Mira.
[9] Sotehi, N., Caractéristiques Thermiques des Parois des Bâtiments et Amélioration de L'isolation. Univ.
Mentouri-Constantine, 2010.
[10] Bianchi, A.-M., Y. Fautrelle, and J. Etay, Transferts thermiques 2004: PPUR presses polytechniques.