LES ÉCHANGEURS DE CHALEUR

WEBER Christophe MANIFICAT André  Bases théoriques 29 Octobre 2015

Niveau de lecture 

Dans le contexte énergétique actuel, la hausse des prix de l’énergie et la raréfaction des sources d’énergies fossiles poussent la société industrielle contemporaine à devenir
de plus en plus performante aussi bien sur l’aspect financier et environnemental que sur l’aspect maitrise de l’énergie. Ce dernier point essentiel aujourd’hui, conduit les
industriels dans une recherche des meilleurs rendements de leurs processus et une baisse de leur consommation énergétique devenus indispensables. L’objectif de
performance énergétique par la maitrise et la rationalisation de l’énergie pour une meilleure efficacité passe en particulier par l’optimisation des différents organes process
(Lire : Les politiques d’efficacité énergétique). Cette nouvelle stratégie de performance énergétique n’est pas seulement judicieuse du point de vue de l’écologie ou de
l’énergie, elle est souvent très rentable pour de nombreux cas, voire même pour des installations de taille réduite. Cette rentabilité impose néanmoins une gestion minutieuse
dans les petits systèmes, davantage encore que dans les grandes installations. À partir d’une situation bien définie, la rentabilité ne peut être garantie que par le recours à
des solutions techniques optimales offrant le meilleur rendement possible par une gestion dynamique des systèmes. C’est là qu’intervient encore, et ceci depuis toujours,
l’échangeur de chaleur. Celui-ci est un élément essentiel de la stratégie de performance énergétique.

Sans conteste l’échangeur de chaleur (ou échangeur thermique) occupe une place primordiale et indispensable dans tous les systèmes thermiques, qu’il soit pour un usage
industriel (chimie, pétrochimie, sidérurgie, agroalimentaire –IAA-, production d’énergie), pour l’automobile, l’aéronautique ou le bâtiment résidentiel ou tertiaire. De manière
générale, on s’accorde à dire que plus de 90% de l’énergie thermique utilisée dans les procédés industriels transite au moins une fois par un échangeur de chaleur[1].

Quelles sont les principales fonctions, les principales technologies et les méthodes de dimensionnement d’un échangeur thermique ? Le  choix et le dimensionnement d’un
échangeur de chaleur, pour une application donnée, dépendent de nombreux paramètres tels que : les propriétés physiques des fluides, leur agressivité, les températures ainsi
que leur pressions de service ou les matériaux. Les contraintes d’encombrement et de maintenance doivent également être prises en compte, ainsi que les considérations
économiques. [2].

1. Principal but et fonctions d’un échangeur thermique

La fonction principale de l’échangeur thermique est de transférer de l’énergie thermique d’un fluide vers un autre à des niveaux de températures distincts. Les fonctionnalités
des échangeurs thermiques sont extrêmement diverses et variées : les principales en sont les suivantes[3][1]:

Préchauffeur ou refroidisseur d’un liquide ou d’un gaz (air par exemple) pour lesquels de très nombreux exemples pourraient être rappelés et qui se
caractérise par une fonction simple : le contrôle de la température du fluide en un point particulier du procédé.
Récupérateur thermique qui permet d’introduire la récente notion de valorisation de l’énergie thermique d’un procédé. La fonction de cet échangeur est
alors d’assurer le transfert d’une capacité thermique maximale afin de permettre une valorisation maximale du rejet thermique sur des critères
énergétiques et économiques.
Réfrigérants (et aéro-réfrigérant lorsque le vecteur de refroidissement est l’air, réfrigérant sur eau) qui assure la dissipation de l’énergie thermique non
valorisable d’un procédé vers le milieu extérieur – fonction essentielle dans de nombreux process.
Capteurs ou émetteurs thermiques dont les équipements ont pour fonction, respectivement, d’associer réception d’énergie thermique et transmission vers
l’usage. Les exemples les plus courants sont les capteurs solaires thermiques ainsi que les radiateurs domestiques qui assurent le confort thermique
dans un bâtiment.
Des humidificateurs ou condenseurs partiels qui assurent la condensation d’une vapeur en mélange avec un gaz incondensable pour obtenir, en fin
d’opération, un gaz appauvri en vapeur : l’exemple rencontré fréquemment est le déshumidificateur d’air humide qui permet d’assurer un contrôle de
l’humidité de l’air en sortie de centrale de traitement d’air (CTA).D’autres exemples industriels peuvent être cités comme la captation de vapeurs
organiques dans les cryo-condenseurs ou les condenseurs de buées fréquemment rencontrés dans les installations de génie chimique.
Évaporateurs qui assurent l’évaporation complète ou partielle d’un liquide dans différents procédés notamment de production d’énergie mécanique (cycle
moteur de Rankine, de Hirn) et de production frigorifique : cycle à compression (PAC), réfrigérateur. Condenseurs qui assurent la condensation complète ou
partielle d’un gaz (vapeur) là encore par exemple pour la production d’énergie mécanique et frigorifique ;
Équipements qui permettent la congélation et la fusion d’une phase liquide ou vapeur grâce à une paroi refroidie en dessous du point triple du fluide. Ces
dispositifs sont utilisés pour assurer la séparation de plusieurs corps, assurer le stockage d’énergie thermique (stockage de glace ou stockage par MCP),
produire une phase solide pour divers usages.
 Le caloduc (Heat Pipe), véritable système thermique diphasique, qui permet notamment d’assurer la dissipation de la chaleur générée par les éléments
électroniques (Microprocesseur ou électronique embarquée), la récupération d’énergie, le maintien en température stable et uniforme[4]. Pour plus
d’informations sur les caloducs vous pouvez vous référer à [8], [9] et [14].

On notera qu’un même échangeur peut assurer plusieurs fonctions dans des domaines applicatifs variés.

2. Technologie des échangeurs

Difficile d’être exhaustif sur la description technologique des échangeurs tant la diversité de ces appareils et de leurs variantes en fonction des constructeurs est importante.
La Figure 1 rappelle de façon synthétique les principales classes d’échangeurs suivant des critères strictement technologiques.

Les échangeurs tubulaires (tubes et calandre, coaxial et autres) sont les plus répandus et
représentaient, en 1998, quasiment 50% des parts de marché en Europe. Les échangeurs dits
compacts (échangeurs à plaques et joints, par exemple), qui ont connus une croissance importante
dans les années 1990, tendent à faire diminuer ce pourcentage mais la grande adaptabilité des
échangeurs tubulaires (tenue en pression et en température, diversité des matériaux) fait qu’il est
difficile d’assurer leurs remplacements par d’autres technologies (Figure 2).

 
Fig. 1 : Classification des échangeurs thermiques (à surface) en fonction de la

3. Fonctionnement général d’un échangeur thermique technologie de fabrication

Il existe de nombreuses technologies d’échangeurs de chaleur, mais tous fonctionnent  suivant les mêmes
processus physiques d’échanges, c’est-à-dire :

– la conduction qui représente les échanges à travers les parois (le plus souvent métalliques),

– la convection qui représente les échanges entre les fluides et les parois,

– le  rayonnement qui représente les échanges radiatifs entre les fluides et les parois (principalement infra
rouge) bien que ce dernier soit souvent négligeable (car pris en compte uniquement pour un fonctionnement à
haute température).
Fig. 2 : Répartition en pourcentage des recettes par type d’échangeur
On distingue également trois modes d’écoulement différents : sur le marché totale de l’échangeur en Europe (données 1998)

– celui à  co-courants : écoulements parallèles des fluides et dans la même direction ;

– celui à contre-courants : écoulements parallèles des fluides mais dans des directions inverses ;

– et celui à courants croisés : écoulements perpendiculaires entre les deux fluides.

Le mode de circulation influe également sur la qualité et l’efficacité du transfert. Un échangeur dit anti-méthodique (ou à co-courants) signifie que l’entrée des deux fluides
(chaud et froid) se situe du même coté de l’échangeur. La configuration dite méthodique est celle ou les fluides circulent à contre-courants. On l’appelle méthodique car c’est
cette configuration qui permet d’obtenir les meilleures performances d’échangeurs. En effet, pour une circulation à co-courant des fluides, la température de sortie du fluide
chaud ne pourra pas être plus basse que la température de sortie du fluide froid tandis qu’à contre courant, il est possible d’abaisser la température du fluide chaud à une
valeur qui tendra vers la température d’entrée du fluide froid. Des représentations graphiques illustreront plus loin ces distinctions.

3.1. Calcul de la puissance thermique d’un échangeur : le bilan thermique

L’équation de bilan d’énergie exprime, conformément au premier principe de la thermodynamique, l’égalité des flux de chaleur échangés par chacun des fluides circulant dans
l’échangeur thermique, en négligeant les pertes vers l’extérieur [1].

La puissance échangée entre les deux fluides, dépend essentiellement de l’écart de température que celui-ci subit lorsque l’échange se réalise par chaleur sensible, ou alors
essentiellement de la valeur de la chaleur latente (correspondant à l’enthalpie de changement d’état) lorsqu’il y a un changement de phase (fluide diphasique).

3.1.1. ECHANGE MONOPHASIQUE
Lorsque l’on parle de transfert monophasique, cela signifie que les fluides en présence restent sous une même et unique phase et ne changent pas d’état (il reste uniquement
à l’état gazeux ou liquide) : les échanges de chaleur se traduisent ainsi uniquement par une élévation ou une baisse de leur température. On parle alors d’échange sous forme
de chaleur sensible.

La puissance échangée (P en Watt) entre deux fluides à l’état monophasique, dépend essentiellement des trois critères suivants :
Le débit massique des fluides (qm en [kg/s]),
Les propriétés thermophysiques des fluides (masse volumique (ρ en [kg/m3]) et chaleur spécifique (Cp en [J/kg-K])),
L’écart de température que vont subir les fluides (ΔT en [K] ou en [°C]), l’échangeur étant représenté (configuration co-courant) sur la Figure 3.

La chaleur se dirigeant toujours du fluide chaud vers le fluide froid et en faisant l’hypothèse que les échanges vers l’extérieur sont nuls (cas d’un échangeur isolé), on
a logiquement:

$$mbox{T entree chaud} > mbox{T sortie chaud}$$

$$mbox{T entree froid} < mbox{T sortie froid}$$

De ce fait, on peut exprimer la puissance échangée par chacun des fluides avec les expressions suivantes qui
impliquent la notion de conservation de l’énergie :

$$P=qm_c times C p_c times Delta T_c = – q m_f times C p_f times Delta T_f$$

Avec :

$$Delta T_c = T_{e_c} – T_{s_c}$$

$$Delta T_f = T_{e_f} – T_{s_f}$$

On notera également que de manière générale, les valeurs des propriétés thermophysiques des fluides sont
retenues comme étant des constantes calculées à la valeur moyenne de température et de pression de chacun
des fluides.

Grâce à l’écriture de ces bilans thermiques, on peut ainsi facilement retrouver une grandeur (par exemple la Fig. 3 : Représentation simplifiée d’un échangeur de chaleur
température de sortie du fluide froid) si on connait l’ensemble des autres grandeurs.

Une représentation assez courante des échanges thermiques entre deux fluides est le diagramme Température-Enthalpie (ou Température-Puissance). Dans le cadre d’un
échange par chaleur sensible uniquement, la représentation est celle de la Figure 4.

Ici, on illustre une configuration à contre courant (comme pour les exemples qui suivent) où le fluide
chaud (en rouge) entre du coté opposé à celui du fluide froid (en bleu) : on rappelle que, dans cette
configuration, le fluide chaud peut être refroidit à un niveau plus bas que la température de sortie du
fluide froid  tandis que le fluide froid peut, lui aussi, être chauffé à une température plus haute que la
température de sortie du fluide chaud.

Fig. 4 : Représentation du diagramme Température-Puissance pour des échanges

monophasique/monophasique
3.1.2. ÉCHANGE DIPHASIQUE
Lorsque l’on parle de transfert diphasique, cela signifie que les fluides en présence vont changer d’état (ou de phase) entre l’entrée et la sortie de l’échangeur : les échanges
de chaleur se traduisent ainsi quasiment uniquement par un changement d’état sans élévation ni baisse de leur température (pour les corps pures à pression constante). On
parle alors d’échange sous forme de chaleur latente.

Les changements de phase les plus fréquemment rencontrés dans le domaine des échangeurs de chaleur sont l’évaporation (Liquide → Gaz) et la condensation (Gaz →
Liquide), notamment dans le domaine du génie climatique (cycle des pompes à chaleur pour la production de froid) et de la production d’énergie électrique (via un cycle de
Rankine pour la production d’énergie mécanique puis électrique via un alternateur).

Les autres changements de phase, tels que la fusion et la solidification, restent encore des domaines plus marginaux bien qu’appliqués de plus en plus par l’utilisation
grandissante des MCP (Matériaux à changement de phase). On citera également la sublimation (Solide→ Gaz) qui est utilisé, par exemple, dans les process de lyophilisation.

La puissance échangée (P en [W]) entre deux fluides qui vont changer de phase dépend essentiellement des trois critères suivants :

Le débit massique des fluides (qm en [kg/s]),

La chaleur latente ou enthalpie de changement de phase (liquide/vapeur le plus couramment) qui est une propriété thermophysique des fluides (noté Lv ou
Δhlv en [J/kg]).
Le titre massique (noté xv) d’une des phases (en général la phase vapeur pour un mélange liquide/vapeur). Ce titre massique est définit comme le rapport
du débit massique de la phase vapeur au débit massique des deux phases (c’est-à-dire le débit massique total de fluide), il est donc compris entre 0
(phase totalement liquide) et 1 (phase totalement vapeur) :

$$x_v=frac{qm_v}{(qm_v+qm_l)}=frac{qm_v}{qm_{total}}$$

Pour une évaporation (ou une condensation) totale, on peut donc dire que l’on passe d’un titre de vapeur de 0 à un titre de vapeur de 1 (ou inversement pour la condensation)
donc Δxv=1. Ainsi, l’expression de la puissance est la suivante :

$$P=qm_c times Delta h_{vl_c}=-qm_f times Delta h_{lv_f}$$

Dans le cadre de cette expression, le fluide chaud (indice c) va par exemple se condenser totalement et le fluide froid va s’évaporer totalement. On a ainsi les expressions
suivantes pour le calcul de l’enthalpie de changement de phase :

$$Delta h_{vl_c}=H_{v_c}-H_{l_c}$$

$$Delta h_{lv_f}=H_{l_c}-H_{v_c}$$

On notera que la détermination de l’enthalpie de changement de phase est dépendante de la pression ou de la température de saturation du fluide (en sachant que ces deux
paramètres sont dépendants l’un de l’autre).

Lorsque la condensation ou l’évaporation est partielle (i.e. 0 < xv < 1), on fait alors intervenir la différence du titre massique entre l’entrée et la sortie du fluide : La puissance
échangée s’exprime ainsi par :

$$P=qm_c times Delta h_{vl_c}times Delta x_v=-qm_f times Delta h_{lv_f} times Delta x_v$$

Pour un changement de phase sur le fluide froid (évaporation) avec un fluide chaud monophasique, on peut par exemple rencontrer l’évolution suivante dans le diagramme
Température-Puissance (Figure 5):
Ici, le fluide froid entre dans l’échangeur à l’état liquide saturé (par exemple eau liquide à 1 bar et
100°C donc xv=0) et s’évapore jusqu’à saturation vapeur (par exemple vapeur d’eau à 1 bar et 100°C
donc xv=1) par l’apport d’énergie thermique issue du fluide chaud qui par conséquent se refroidit.

Pour un changement de phase sur le fluide chaud (condensation) avec un fluide froid monophasique, on
peut avoir  l’évolution du diagramme Température-Puissance (Figure 6) :

Ici, nous somme dans la situation inverse : le fluide chaud entre dans l’échangeur à l’état vapeur
saturée et se condense jusqu’à saturation par l’extraction d’énergie thermique issue du réchauffement
du fluide froid.
Fig. 5 : Représentation du diagramme Température-Puissance pour un échange
de type Evaporation / Monophasique

Fig. 6 : Représentation du diagramme Température-Puissance pour un échange

de type Condensation / Monophasique
3.1.3. ÉCHANGE MONOPHASIQUE ET DIPHASIQUE
Lors d’un échange par changement de phase, il est fort probable que les fluides disponibles ne soit pas exactement à l’état saturé (pour l’eau, l’état de saturation à la
pression atmosphérique est atteint à environ 100°C). De ce fait, l’échangeur va devoir coupler les échanges monophasiques (sensibles) et diphasiques (changements de
phase=latent). Pour le calcul de la puissance, on devra alors additionner les deux expressions vues plus amont pour calculer la puissance totale échangée.

Par exemple, chauffer un fluide liquide jusqu’à saturation puis l’évaporer partiellement ou totalement (ou refroidir un fluide gazeux jusqu’à saturation pour le condenser
partiellement ou totalement), dans ce mode de fonctionnement, on pourra donc additionner les deux expressions vues en amont pour calculer la puissance totale échangée :

$$P=qm_ctimes Cp_c times Delta T_{c_{v_{sat}}}+qm_c times Delta h_{vl_c} times Delta x_{v_c}  =qm_c times (Cp_c times Delta T_{c_{v_{sat}}}+ Delta
h_{vl_c}times Delta x_{v_c})$$

$$-P=qm_ftimes Cp_f times Delta T_{f_{l_{sat}}}+qm_f times Delta h_{lv_f} times Delta x_{v_f}  =qm_f times (Cp_f times Delta T_{f_{l_{sat}}}+ Delta h_{lv_f}times
Delta x_{v_f})$$

Enfin, il arrive fréquemment de devoir surchauffer de la vapeur à la suite de l’évaporation totale ou bien de sous refroidir un liquide à la suite d’une condensation totale.
L’expression de la puissance devient alors la suivante :

$$P=qm_ctimes Cp_{c_v} times Delta T_{c_{v_{sat}}}+qm_c times Delta h_{vl_c}+qm_ctimes Cp_{c_l} times Delta T_{c_{l_{sat}}}$$

$$=qm_c times ( Cp_{c_v} times Delta T_{c_{v_{sat}}}+ Delta h_{vl_c}+Cp_{c_l} times Delta T_{c_{l_{sat}}})$$

$$-P=qm_ftimes Cp_{f_l} times Delta T_{f_{l_{sat}}}+qm_f times Delta h_{lv_f}+qm_ftimes Cp_{f_v} times Delta T_{f_{v_{sat}}}$$

$$=qm_f times (Cp_{f_l} times Delta T_{f_{l_{sat}}}+Delta h_{lv_f}+Cp_{f_v} times Delta T_{f_{v_{sat}}})$$

Avec :

$$Delta T_{c_{v_{sat}}}=T{e_c}-T_{sat(x=1)_c}$$ $$Delta T_{f_{v_{sat}}}=T{s_f}-

T_{sat(x=1)_f}$$
 
 
$$Delta T_{c_{l_{sat}}}=T_{sat(x=0)_c}-T{s_c}$$
$$Delta T_{f_{l_{sat}}}=T_{sat(x=0)_f}-T{e_c}$$

Pour un corps pur, on souligne que :

$$T_{c_{l_{sat(x=1)}}}=T_{c_{v_{sat(x=0)}}} $$T_{f_{l_{sat(x=1)}}}=T_{f_{v_{sat(x=0)}}}$
$$ $

Pour un fluide chauffé, puis évaporé puis surchauffé par un fluide chaud qui est refroidi, condensé puis sous refroidi, on peut avoir l’évolution du diagramme température-
puissance (Figure 7).

Ici, l’énergie thermique permettant le chauffage, l’évaporation et la surchauffe du fluide froid est fournie par le refroidissement, la condensation ainsi que le sous-
refroidissement du fluide chaud.

Nous avons vu l’expression de la puissance d’un échangeur par l’équation de la conservation de l’énergie à travers l’écriture des différentes expressions et en fonction des
modes de transferts thermiques (sensible, latent ou couplage sensible + latent).  Nous n’avons cependant pas abordé le fond du problème du dimensionnement de
l’échangeur. En effet, ces premiers calculs permettent de connaitre les principales grandeurs (puissance, débit, température) qui vont conditionner le dimensionnement, c’est-
à-dire les dimensions (la taille) de l’échangeur.

Pour cela, nous allons introduire les deux principales méthodes de dimensionnement qui nécessitent toutes deux le calcul du coefficient d’échange thermique global.
 Fig. 7 : Représentation du diagramme Température-Puissance pour un échange
de type Evaporation / Condensation
3.2. Méthodes de dimensionnement
Au côté de l’équation de bilan d’énergie, la résolution des problèmes d’échange thermique nécessite l’écriture d’une seconde équation, de nature tout à fait différente de la
première, qui pour rappel est une équation thermodynamique et traduit le principe de conservation de l’énergie. Cette seconde équation est l’équation de transfert de chaleur
qui traduit la relation entre le flux de chaleur échangé et la différence de températures disponibles entre les deux fluides qui circulent dans l’échangeur thermique.

Le dimensionnement d’un échangeur consiste à calculer la surface d’échange nécessaire pour obtenir les performances désirées (puissance, température de sortie). Il faudra
par la suite déterminer les pertes de charge (pertes de pression) que va engendrer la circulation des fluides (par la friction contre les parois du fait de la viscosité) afin de
vérifier la concordance avec le cahier des charges qui impose souvent une limite haute en terme de valeur d’une perte de charge. Cette valeur haute correspondant, la plupart
du temps, à une capacité de pompage existante et donc à ne pas dépasser[5].

Pour connaitre la surface d’échange de l’échangeur, il faut dans un premier temps connaitre les grandeurs dimensionnantes telles que la puissance, les températures, les
débits (voir plus haut la détermination de ces grandeurs par les bilans thermiques).

Ensuite, il faut calculer, et ce quelle que soit la méthode de dimensionnement utilisée, les coefficients d’échanges thermiques partiels et le coefficient d’échange thermique
global.

3.2.1. COEFFICIENT D’ÉCHANGE GLOBAL ET PARTIEL

Pour expliquer clairement la notion de coefficient d’échange, on considère une paroi métallique au travers de laquelle la chaleur est transférée par conduction, cette paroi
étant en contact avec un fluide chaud et sur son autre face avec un fluide froid : cette situation se traduit par un transfert de chaleur entre les deux fluides. On peut exprimer
par l’équation qui suit la relation entre le flux de chaleur transféré P, l’écart de température ΔT des deux fluides (le chaud (indice c) et le froid (indice f)) et l’aire S de la paroi
d’échange:

$$P=Utimes S times (T_c-T_f) = frac{(T_c-T_f)}{R_{global}}$$

Cette équation, plus connue sous le nom de loi de Newton, caractérise le transfert de chaleur entre les deux fluides et introduit ainsi de nouvelles grandeurs :

le facteur U dont l’unité est en [W/m²-K] qui exprime la qualité du transfert de chaleur entre les deux fluides s’écoulant de part et d’autre de la paroi ; ce
facteur est désigné comme le coefficient d’échange thermique global entre les deux fluides (CETG) ;
le facteur Rg, inverse du produit U.A, dont l’unité est en [K/W] appelé résistance thermique globale (RTG) dont le concept est issu directement de l’analogie
avec la résistance électrique de la loi d’Ohm dans le domaine du génie électrique ; cette analogie associe la différence de température à la différence de
potentiel et le flux de chaleur à l’intensité électrique. La loi de Newton est alors analogue à la loi d’Ohm.

Comme nous l’avons dit, une grande part du calcul d’un échangeur repose sur l’évaluation de ce coefficient d’échange global dont la valeur peut être très différente suivant le
type de paroi en contact avec les fluides, les fluides et leurs propriétés, l’état du fluide (simple phase gazeux ou liquide, double phase), les modes de transfert entre les
fluides et la paroi (réchauffement ou refroidissement, évaporation, condensation…). Bien que le calcul d’un échangeur nécessite l’évaluation précise du coefficient d’échange
global, il est indispensable de connaître les ordres de  grandeur des coefficients d’échanges globaux dans diverses conditions, notamment pour juger rapidement de la
justesse d’un calcul (Tableau 1).

Tableau 1 : Ordre de grandeur du coefficient d’échange global

Situation physique U en W/m²-K

condenseur de vapeur d’eau 1100-1600

réchauffeur d’eau 1100-5600

condenseur de fluide frigorigène de type HFC avec eau comme fluide secondaire 500-1000

échangeur eau-eau 850-1800

batteries à ailettes (tube : eau, ailette : air) 25-70

échangeur eau/huile 110-350

condenseur d’ammoniac avec eau comme fluide secondaire 850-1400

échangeur gaz/gaz 10-40

 évaporateur de HFC à tubes et calandre à détente directe 800-1600

Evaporateur noyé de HFC à tubes et calandre 500-1200

On notera que la valeur de ce coefficient est avant tout gouvernée par la résistance thermique convective (c’est-à-dire le transfert entre la paroi et le fluide), la résistance par
conduction étant très faible dans la majorité des cas (car les parois sont de faible épaisseur et de conductivité thermique élevée). Ce point souligne alors l’importance d’une
bonne évaluation des transferts thermiques par convection.

Le coefficient d’échange global U se calcul à partir des différentes résistances thermiques de convection (interne et externe), de conduction (paroi solide) et éventuellement
des résistances d’encrassement (dépôt biologique, tartre…etc.) :

1. La résistance de conduction due à la présence d’une paroi solide s’exprime grâce à l’équation de la chaleur, elle est fonction de la géométrie de la paroi (coordonnées plan,
cylindrique ou sphérique). Elle dépend aussi de la nature du matériau (conductivité thermique).

Si la conduction est monodirectionnelle comme cela est le cas dans une paroi plane, on utilise l’expression suivante :

$$R_{conduction}=frac{e}{lambda times S}$$

Si la conduction est réalisée dans une paroi cylindrique, on utilise l’expression qui suit :

$$R_{conduction}=frac{1}{2pi times lambda times L} times ln left(frac{D_{ext}}{D_{int}}right)$$

2. Les résistances thermiques de convection sont en fait l’inverse du coefficient d’échange convectif h (interne et externe). Chaque coefficient d’échange convectif est
calculé à partir de la valeur du nombre adimensionnel de Nusselt qui est évalué par des corrélations empiriques (par expérimentation) données en fonction de la géométrie
(diamètre par exemple), des régimes d’écoulement (laminaire, transitoire ou turbulent) et de nombres adimensionnels (Reynolds (Re), Prandtl (Pr)). De manière générale, les
expressions du nombre de Nusselt (Nu) apparaissent sous la forme suivante :

$$N u=htimes frac{D_h}{lambda}=Atimes R e^x times P r^y$$

Ces corrélations sont spécifiques aux produits de nombreux constructeurs et possèdent des domaines de validité bien précis (géométrie, intervalle de la valeur du nombre de
Reynolds et du nombre de Prandtl) qu’il est impératif de respecter pour minimiser les erreurs sur ces données. On trouve ces corrélations dans différents ouvrages notamment
au GRETh où une revue exhaustive des méthodes existantes est disponible[6].

À titre d’exemple, la plus courante de ces expressions – appelée corrélation de Colburn (nom de l’auteur) – permet de déterminer le coefficient d’échange convectif d’un fluide
s’écoulant dans une section circulaire (tube) en régime turbulent (régime d’écoulement le plus favorable aux échanges thermiques convectifs et intervenant pour un nombre
de Reynolds > 10 000 dans un tube de section circulaire):

$$N u=0.023times R e^{0.8} times P r^{1/3}$$

Comme : $$N u=htimes frac{D_h}{lambda}$$

On peut alors déterminer le coefficient d’échange convectif partiel h :

$$h=frac{lambda}{D_h}times 0.023times R e^{0.8} times P r^{1/3}$$

La résistance thermique convective s’exprime alors par :

$$R_{convection}=frac{1}{h times S}$$

3. Enfin, les valeurs des résistances d’encrassement. En général, l’encrassement des échangeurs de chaleur est défini comme le dépôt de matériaux ou de substances non
désirées sur les surfaces d’échanges. Ce dépôt peut être constitué de cristaux, de sédiments, de résidus biologiques, de produits d’une réaction chimique ou bien encore être
la combinaison de plusieurs de ces éléments. C’est un processus qui est dynamique et qui se déroule de façon continue ou ponctuelle, généralement jusqu’à atteindre un état
d’équilibre (valeur asymptotique) ou jusqu’à ce qu’il soit nécessaire de nettoyer l’échangeur – la période d’encrassement pouvant varier de quelques heures à plusieurs mois,
voire plusieurs années. Ces résistances thermiques d’encrassement sont en général données par les standards TEMA (Tubular Exchanger Manufacturer Association[7]) qui
font office de référence (sauf retour d’expérience propre à l’exploitant). Leur prise en compte dans le calcul du coefficient d’échange global permet de simuler le
fonctionnement de l’échangeur en conditions dégradées (avec une résistance d’encrassement suite à un nombre d’heures de fonctionnement) et permet de faire un
surdimensionnement justifié due à une potentielle future baisse de performance. En effet, l’ajout d’une ou plusieurs résistances d’encrassement a pour effet néfaste de
dégrader la valeur du coefficient d’échange thermique global, induisant une hausse de la surface d’échange nécessaire et donc un surplus en cout matière et conception
(figure 8)[8].

Pour plus d’informations sur l’encrassement, vous pouvez vous référer à [1], [3] ainsi qu’à [12] et [13] qui traitent le sujet de manière assez exhaustive.

La résistance globale (et en conséquence U) se déduit de la règle d’additivité des résistances

partielles :

$$R_{global}= frac{1}{U times S} = sum R_i$$

On notera que le coefficient d’échange global doit être défini par rapport à une surface d’échange de
référence. Dans le cas, d’une paroi plane, cette surface de référence est évidente puisqu’unique : c’est
la surface de la paroi. Dans le cas d’une paroi tubulaire (cylindrique), deux surfaces d’échange de
référence peuvent être retenues : la paroi interne ou externe. On peut alors choisir indistinctement
l’une ou l’autre à la condition de maintenir ce choix  tout au long du calcul.

Ainsi, le coefficient d’échange global U (défini par rapport à la surface d’échange interne – notée avec
un indice int – pour les échangeurs tubulaires) est déterminé par la relation suivante dans laquelle
interviennent la résistance de paroi, les résistances convectives et les résistances d’encrassement :
Fig. 8 : Schématisation des différentes résistances thermiques d’une paroi [9]
 $$R_{global}= R_{convection_{ext}}+R_{encrassement_{ext}}+R_{conduction}+R_{encrassement_{int}}+R_{convection_{int}}$$

Cette expression prend différentes formes suivant le type d’échangeur considéré (en négligeant les résistances d’encrassement) :

Pour un échangeur tubulaire (de diamètre externe Dext et interne Dint) :

$$frac{1}{U_{int}times S_{int}}=frac{1}{h_{ext}times S_{ext}}+frac{1}{2pi times lambda times L}times ln left(frac{D_{ext}}{D_{int}} right)+frac{1}{h_{int}times

S_{int}}$$

        En sachant que :

$$S_{int}=pi times D_{int} times L$$

        On a finalement :

$$frac{1}{U_{int}}=frac{1}{h_{ext}}times left(frac{S_{int}}{S_{ext}}right)+frac{D_{int}}{2lambda} ln left(frac{D_{ext}}{D_{int}} right) + frac{1}{h_{int}}$$

Pour un échangeur à plaques (d’épaisseur e) :

$$frac{1}{U}=frac{1}{h_{c}}+frac{e}{lambda} + frac{1}{h_{f}}$$

3.2.2. MÉTHODE DE L’ÉCART DE TEMPÉRATURE MOYEN LOGARITHMIQUE

La représentation des évolutions de températures moyennes des fluides en fonction d’une grandeur spatiale (par exemple, x, la distance d’un des fluides par rapport à l’entrée
de l’échangeur) suppose la résolution simultanée des deux équations :

équation de bilan énergétique ;

équation de transfert de chaleur.

L’écriture des équations d’énergie et de transfert, pour un échangeur à co-courant sur un élément de surface dS, est la suivante :

$$mbox{d}P=qm_c times C p_c times mbox{d} T_c = – q m_f times C p_f times mbox{d} T_f$$

$$mbox{d}P=Utimes mbox{d} S times Delta T$$

À partir de ces deux équations, nous déduisons l’équation différentielle :

$$mbox{d}(Delta T)=mbox{d}T_c-mbox{d} T_c = -mbox{d}P times left(frac{1}{qm_c times C p_c }+frac{1}{ q m_f times C p_f}right) times mbox{d} S$$

D’où :

$$frac{mbox{d}(Delta T)}{Delta T}=-U times left(frac{1}{qm_c times C p_c }+frac{1}{ q m_f times C p_f}right) times mbox{d} S$$

L’intégration de cette dernière équation donne l’évolution spatiale de l’écart de température entre les fluides à partir de celui existant entre les fluides à l’entrée de
l’échangeur (x=0) en supposant constant le coefficient d’échange U et les capacités thermiques Cp :

$$Delta T(x)=Delta T(x=0) times exp left(-U times S times frac{x}{L} times left(frac{1}{qm_c times C p_c }+frac{1}{ q m_f times C p_f}right)  right)$$

Pour un échangeur à contre-courant, l’expression précédente devient :

$$Delta T(x)=Delta T(x=0) times exp left(-U times S times frac{x}{L} times left(frac{1}{qm_c times C p_c }-frac{1}{ q m_f times C p_f}right)  right)$$

Cette méthode est la seule méthode analytique permettant une représentation spatiale exacte de la température moyenne des fluides. On aura noté toutefois qu’elle suppose
des hypothèses fortes telles que, entre autres, la constance du CETG tout le long de l’échangeur.

Les équations précédentes aboutissent, après intégration aux bornes de l’échangeur, à l’expression suivante qui relie la puissance échangée, le coefficient d’échange global et
la surface d’échange totale :

$$P=Utimes Stimes frac{Delta T_a – Delta T_b}{ln left(frac{Delta T_a}{Delta T_b}right)}=Utimes S times Delta T_{ln}$$

La différence de température moyenne logarithmique ΔTln est définie de la manière suivante :

$$Delta T_{ln}=frac{Delta T_a – Delta T_b}{ln left(frac{Delta T_a}{Delta T_b}right)}$$

Cette expression fait intervenir les paramètres ΔTa et ΔTb qui représentent les écarts de températures aux deux bornes de l’échangeur. Elle s’applique aussi bien aux
échangeurs à co-courants qu’à contre-courants (Figure 9).

Dans le cas où la condition de courant parallèle des fluides n’est pas respectée, c’est-à-dire pour une
configuration à courant croisés (comme c’est le cas pour des batteries à ailettes), un aménagement
de la méthode est réalisé en introduisant un facteur correctif F qui prend en compte la spécificité de
l’échangeur. Ainsi, la différence de température moyenne logarithmique correspond à la valeur en
contre-courants multipliée par le facteur de correction F (répertorié dans des abaques spécifiques aux
configurations à courants croisés possible (1 fluide brassé, 2 fluides brassés…etc.)). Ce facteur F à
donc une valeur comprise entre 0 et 1.
La méthode de dimensionnement par l’écart de température logarithmique se résume donc à une Fig. 9 : Représentation de ΔTa et ΔTb en fonction de la configuration
d’écoulement
succession d’opérations simples dont l’organigramme suivant synthétise les principales étapes (Figure
10)[10].

Cette méthode permet essentiellement la détermination de la surface S de l’échangeur et ne s’applique

que dans des conditions strictes[11]:

le coefficient d’échange global est supposé constant,

les échangeurs sont à circulation parallèle (co ou contre-courant) ou corrigé d’un facteur
F si ce n’est pas le cas,
les capacités thermiques massiques des fluides sont supposées constantes tout au long
du processus.

La méthode de l’écart logarithmique, si elle est commode pour évaluer une surface d’échange (les
quatre températures des fluides étant connues a priori) présente une limite d’utilisation lorsqu’il s’agit
d’évaluer la puissance thermique échangée pour un échangeur de surface connue. Une méthode – celle Fig. 10 : Démarche de dimensionnement par la méthode du ΔTln
de l’efficacité et du nombre d’unités de transfert (NUT) – est proposée pour éviter cette difficulté.
L’objet de cette dernière méthode de dimensionnement est d’évaluer la puissance thermique transférée dans un échangeur dont la géométrie est connue.

3.2.3. MÉTHODE DE L’EFFICACITÉ ET DU NOMBRE D’UNITÉ DE TRANSFERT

On définira l’efficacité ɛ d’un échangeur comme le rapport du flux de chaleur réellement échangé dans l’échangeur (d’arrangement géométrique quelconque) au flux de chaleur
maximum qu’il est possible d’échanger (c’est-à-dire le flux échangé avec une surface d’échange infinie) pour un échangeur à circulation à contre-courant fonctionnant avec
des débits et des températures d’entrée des fluides identiques à ceux du cas considéré.

Le flux de chaleur maximum qu’il est possible d’échanger est déterminé de la façon suivante : la variation de température maximale que peut subir le fluide de plus faible débit
de capacité thermique (le débit de capacité thermique noté $$/dot{C}$$ et exprimé en [W/K] est le produit du débit massique (qm) par la chaleur spécifique (Cp)) est égale
la différence entre les températures d’entrée des fluides.

Le fluide qui est  susceptible de subir cette variation maximale est le fluide qui présente le débit de capacité thermique massique minimale :

$$dot{C}_{min}=min (qm_c times Cp_c  ;qm_f×Cp_f )=min (dot{C}_c;dot{C}_f)$$

 Ce résultat serait obtenu avec une surface d’échange infinie et un échangeur à contre-courant. Le flux de chaleur maximum possible est donc :

$$P_{max}=dot{C}_{min}times (T_{e_c}-T_{e_f})$$

L’efficacité de l’échangeur thermique peut alors s’écrire :

si le coté chaud présente le produit qm*Cp minimum :

$$epsilon = frac{P_{reel}}{P_{max}}=frac{dot{C}_{c}times (T_{e_c}-T_{s_c})}{dot{C}_{c}times (T_{e_c}-T_{e_f})}=frac{ (T_{e_c}-T_{s_c})}{ (T_{e_c}-

T_{e_f})}=frac{Delta T_c}{Delta T_{max}}$$

si le coté froid présente le produit qm*Cp minimum :

$$epsilon = frac{P_{reel}}{P_{max}}=frac{dot{C}_{f}times (T_{s_f}-T_{e_f})}{dot{C}_{f}times (T_{e_c}-T_{e_f})}=frac{ (T_{s_f}-T_{e_f})}{ (T_{e_c}-

T_{e_f})}=frac{Delta T_f}{Delta T_{max}}$$

Ainsi, d’une manière générale, l’efficacité est égale à la variation de température sur le fluide de $$dot{C}$$ minimum sur l’écart de température maximale dans l’échangeur.

Lorsque l’on ne connait pas les températures de sortie (ce qui est souvent le cas dans le cadre d’une simulation), l’efficacité d’un échangeur ɛ est déterminée par le
coefficient d’échange global U, la surface d’échange S et la valeur du paramètre $$dot{C}_{min}$$ et $$dot{C}_{max}$$ via l’introduction du terme de rapport de débit de
capacité thermique noté $$C=frac{dot{C}_{min}}{dot{C}_{max}}$$.

Les expressions qui relient ces grandeurs se déduisent de l’intégration des équations de bilan et de transfert.

Ainsi, on peut déduire,

pour une circulation des fluides à contre-courants :

$$epsilon = 1- frac{exp (-NUT times (1+C))}{1+C}$$

pour une circulation des fluides à co-courants :

$$epsilon = frac{1-exp (-NUT times (1+C))}{1-C times exp (-NUT times (1+C))}$$

Le rapport $$U times frac{S}{dot{C}_{min}}$$ est désigné comme le  NUT (Nombre d’Unités de Transfert). Cette grandeur est égale, si l’on applique les équations de bilan et
de transfert à :

$$NUT = U times  frac{S}{dot{C}_{min}} = frac{P}{Delta T_{ln}} times frac{T_{e_{(C_{min})}} -T_{s_{(C_{min})}}}{P}=  frac{T_{e_{(C_{min})}} -T_{s_{(C_{min})}}}
{Delta T_{ln}}$$

Soit le rapport de la variation de température sur le fluide de $$dot{C}$$  le plus faible, à l’écart de température logarithmique.
La méthode dite de l’efficacité et du NUT (méthode ɛ-NUT) est une méthode aisée de calcul des performances d’un échangeur thermique (diagnostic). En effet  la formulation
proposée permet à partir du coefficient d’échange global U, la surface d’échange S, des débits massiques et chaleurs spécifiques :

en premier lieu d’évaluer le nombre d’unités de transfert NUT ;

en deuxième lieu, à partir des expressions précédentes ou de celles issues d’abaques pour des configurations plus complexes d’échangeurs, de calculer
l’efficacité de l’échangeur ;
l’efficacité connue, il est aisé de déterminer latempérature de sortie de l’un des fluides (celui présentant le minimum). Ensuite, la température de sortie de
l’autre fluide se déduit des bilans thermiques.

Les limitations de cette méthode sont identiques à celles de la méthode de l’écart de température moyen  logarithmique. La mise en œuvre de cette méthode nécessite la
connaissance de la relation entre l’efficacité et le nombre d’unités de transfert. Cette relation sous forme d’abaque ou d’expressions algébriques dépend :

du mode de circulation des fluides,

du type d’échangeur de chaleur,
du mélange ou non des filets fluides lorsque l’écoulement est à courant-croisé.

Dans le cas de condenseur ou d’évaporateur (avec  température constante du fluide à 2 phases), on considère Cmax tend vers une valeur infinie car la température de
changement de phase est constante correspondant au fluide frigorigène qui s’évapore ou se condense.

3.3. Remarques sur les méthodes de calcul des échangeurs de chaleur

Pour le dimensionnement (méthode du ΔTln) comme pour la simulation (méthode ɛ-NUT) du fonctionnement d’un échangeur, trois situations types peuvent être rencontrées :

1. Dans l’échangeur circulent des fluides présentant des propriétés thermophysiques constantes et le coefficient d’échange admet une valeur sensiblement
constante : la méthode de calcul dite point consiste à appliquer sur l’ensemble de l’échangeur l’une des méthodes présentée.
2.   Dans l’échangeur circulent des fluides dont l’un au moins présente une variation brutale de ses propriétés : cas d’un fluide liquide qui devient vapeur puis
d’une vapeur surchauffée, par exemple. On ne peut légitimement pas accepter l’hypothèse de constance des propriétés et du coefficient U. La méthode par
zones qui consiste à découper virtuellement l’échangeur en un nombre limité de zones (2 voire 3) caractérisées par la constance (ou quasi constance) des
propriétés et de U dans chaque zone, permet  d’utiliser les méthodes classiques précédentes pour chacune d’elles (Figure 11).
3. Dans l’échangeur circulent des fluides dont les propriétés, le coefficient U varient
continument de façon importante (par exemple pour un fluide à l’état supercritique très
proche du point critique) : les méthodes précédentes ne peuvent plus s’appliquer. Les
méthodes numériques à éléments ou volumes finis doivent s’appliquer en considérant
l’échangeur comme une succession de mailles élémentaires pour lesquelles on établit
les équations de bilan thermique, les équations locales de transfert de chaleur et les
évaluations de pertes de pression.Ces méthodes de complexité croissante sont toutes
utilisées. La dernière méthode (par éléments ou volumes finis), de par sa complexité et
la lourdeur des calculs qu’elle implique, nécessite des outils numériques lourds (voir
section 7).

Fig. 11 : Représentation du découpage par zone pour le dimensionnement dans

  un diagramme Température-Puissance pour un échange de type Evaporation /
Condensation
4. Critères de performance d’un échangeur

Il existe plusieurs critères de performance d’un échangeur :

L’efficacité d’un échangeur est un critère de performance important.

Le coefficient d’échange global U (ou K) est aussi un des plus utilisé. Plus ce dernier est important, plus l’échangeur est performant. Il faut néanmoins que
les pertes de charge soient les plus faibles possibles (mais suffisante afin d’assurer une bonne distribution des fluides) mais soient totalement
dépendantes du design de l’échangeur.
Enfin, le critère de compacité est aussi important car il est défini par le rapport de l’aire de la surface d’échange au volume de l’échangeur. R.K. Shah[12]
[4] propose qu’un échangeur soit considéré comme compact si sa compacité (rapport de la surface d’échange sur le volume de l’échangeur) est supérieure
à 700 m²/m3 ; cette valeur est susceptible de varier de 500 à 800 m²/m3.

5. Le choix technologique : Identi cation des caractéristiques d’usages nécessaires à une sélection adaptée d’un
échangeur thermique[13]

La sélection d’une technologie est une phase essentielle dans le projet d’un système ou d’un équipement thermique. Les critères techniques et économiques sont nombreux
alors que le choix des technologies d’échangeurs thermiques est extrêmement varié.

La démarche entreprise pour aider le lecteur à assurer une sélection raisonnée est :

d’identifier l’ensemble des critères techniques de choix ;

de quantifier dans la mesure du possible des critères qualitatifs tels que les caractéristiques encrassantes d’un fluide, les risques associés à des défauts
d’étanchéité…
Différentes caractéristiques sont nécessaires pour amorcer les étapes de choix technologique d’un échangeur thermique. Ces grandeurs sont en premier lieu les données du
procédé dans lequel l’échangeur doit être intégré, puis des caractéristiques associées aux contraintes spécifiques.

– Les fluides de transfert, moteur et frigorigène du procédé

En premier lieu, il faut pouvoir définir le plus précisément possible les différentes données du process dans lequel l’échangeur est prévu d’être intégré. Cela passe
inévitablement par la définition des fluides utilisés et de leurs conditions d’usage.  La définition des fluides du procédé est bien évidemment d’une importance capitale car il
permet de définir les différentes propriétés thermophysiques des fluides qui seront mis en jeu et qui servent aux calculs des différentes grandeurs dimensionnelles et
adimensionnelles.

– L’état des fluides de transfert (monophasique, diphasique) 

– Les conditions extrêmes d’usage de l’échangeur : les  températures maximales d’usage de l’échangeur ; la  pression maximale de service de l’échangeur sur l’un ou l’autre
des circuits ; le critère de pertes de charge maximales acceptables (ou de la perte de charge maximale).

– L’encrassement des échangeurs thermiques

Ce critère demande une connaissance du potentiel encrassant du fluide qui peut induire plusieurs effets (Bouchage, Colmatage) et qui se traduisent par des effets sur la
performance globale (Dégradation des performances thermiques, Augmentation des pertes de charge, Surdimensionnement). Ainsi, la question de l’encrassement est
primordiale. Lors du fonctionnement avec des fluides à potentiel encrassant, il se peut qu’au bout d’un certain temps de fonctionnement de l’échangeur, celui-ci ne soit plus
en mesure de fournir la puissance nécessaire. Il peut alors être judicieux d’écarter certaines typologies d’échangeurs qui sont très sensibles à ces phénomènes pour éviter,
soit une dégradation trop importante et/ou rapide des performances de l’échangeur, soit pour préconiser une technologie d’échangeur qui permette une inspection et un
nettoyage aisés pour assurer – et ceci tout au long de sa vie – les performances escomptées.

– Critère sur les exigences d’étanchéité de l’échangeur et les risques de fuite des fluides

La question des risques de fuite de fluides (et incidemment la qualité de l’étanchéité de l’échangeur) peut être un critère important notamment vis-à-vis de la sécurité de
l’installation, les risques sur l’environnement. Ce critère peut être très contraignant en matière de choix technologique et peut intervenir dans le dimensionnement de
l’échangeur lui-même par la nécessité d’une épaisseur de paroi importante (prenant en compte par exemple des scenarii de corrosion annuelle…), voir même des exigences de
double paroi.

– Critère associé à  la nécessité d’inspection des équipements techniques

Ce critère est d’une part attribué par des obligations d’ordre réglementaire de type DESP (Directive des équipements sous pression) mais est aussi intimement lié aux
critères précédents de l’encrassement et de l’étanchéité.

– Critère d’encombrement/compacité 

La compacité d’un échangeur permet ainsi de juger rapidement de la performance d’un échangeur thermique. Ce critère peut aussi être déterminant en termes de choix
technologique et de dimensionnement suivant la place (surface et/ou volume) que demande l’installation de l’échangeur dans le process.

– Critère de coût de l’échangeur

Ce critère reste le point sensible des relations fabricants/utilisateurs, surtout pour des applications bien connues, maitrisées et standardisées. Pour des applications plus
spécifiques ou la réussite du projet est fortement conditionnée par la conception et la fabrication de l’échangeur, le critère du prix peut prendre moins d’importance.

6. Logiciels de dimensionnement des échangeurs de chaleur

Les problèmes de dimensionnement de ce composant primordial et essentiel concernent avant tout le traitement des aspects thermo-hydrauliques. La complexité du
problème de calcul thermique des échangeurs vient de plusieurs facteurs et contraintes, notamment la diversité des technologies des appareils et la nature de la physique
des écoulements. Il faut donc pour cela[14] [3]:

caractériser un modèle d’écoulement, surtout lorsque celui-ci est complexe ;

choisir des corrélations empiriques adaptées : ce problème peut être précisément résolu par les fabricants en réalisant des campagnes d’essais qui
permettent d’établir des corrélations propres à chacune de ses gammes d’échangeur ; il est plus difficile pour les exploitants qui bénéficie de l’utilisation
de technologie issue de plusieurs fabricants : dans ce cas, ils peuvent alors se rapporter aux corrélations qui sont présentées dans la littérature[15] ;
adopter un algorithme de résolution : on peut en effet choisir des algorithmes de type dimensionnement ou de type simulation qui utilisent des méthodes
analytiques ou numériques.

On fournit ci-dessous une liste non exhaustive des principaux logiciels de dimensionnement d’échangeur thermique :

Xchanger suitede chez HTRI (Heat Transfer Research Incorporation) : logiciel de référence qui inclut le dimensionnement des principales technologies
d’échangeurs pour quasiment toutes les configurations standardisées (tubes et calandre, plaques et joints, plaques et ailettes, batterie à ailette, Kettle …
etc.) ;
EDR (Exchanger Design and Rating) de chez AspenTech (anciennement HTFS) : logiciel de référence qui inclut également le dimensionnement des
principales technologies d’échangeurs pour quasimenttoutes les configurations standardisées ;
Prosec de chez Prosim : logiciels très puissants pour la simulation des échangeurs à plaques et ondes multifluides (applications cryogéniques
(liquéfaction d’air) essentiellement) ;
CETUC du GRETh (Groupement pour la Recherche sur les Echangeurs Thermiques) : logiciel validé par expérimentation qui permet le dimensionnement
d’échangeurs à tubes et calandre uniquement ;
CEPAJ du GRETh : logiciel validé par expérimentation qui permet le dimensionnement d’échangeurs à plaques uniquement ;
EchTherm de NeoTherm Consulting et diffusé exclusivement par le GRETh : boite à outils du thermicien qui inclut une multitude d’outils numériques de
calculs de coefficients d’échanges et de pertes de charge pour la majeure partie des technologies d’échangeurs et pour quasiment tous les modes de
fonctionnement (monophasique et diphasique) ; il inclut également des modules simplifiée de dimensionnement (tube et calandre, plaques et joint,
batterie à ailette, caloduc…etc.) ainsi qu’un outil d’aide au choix technologique ;
 Les logiciels de CFD (Computational Fluid Dynamic) pour les méthodes par volumes ou éléments finis qui permettent la résolution couplée des équations
de transferts thermiques et de la mécanique des fluides (Navier-Stockes). On citera les principaux logiciels comme Fluentde chez ANSYS, Comsol
Multiphysics de chez COMSOL, SolidWorks Flow simulation de chez Dassault systèmes, STARCCM+ de CD-Adapco : ces logiciels pointus sont réservés à la
résolution de problèmes spécifiques liés au dimensionnement d’échangeur. Ils sont trop lourd (tarif et prise en main) pour la réalisation d’un
dimensionnement simple mais souvent indispensables pour la recherche de nouvelles géométries complexes permettant d’améliorer les échanges
thermiques, par exemple.

Cette liste n’est pas  exhaustive de tous les logiciels disponibles sur le marché mais elle permet d’apprécier les principaux et surtout ceux qui sont utilisés par les
industriels.  On voit que le nombre de logiciels reste réduit et que l’utilisateur n’a que peu de choix sur les outils de dimensionnement d’échangeur. Certains développements
ont permis de mettre en place des outils robustes, constamment mis à jour (cf. HTRI ; AspenTech) et qui sont parfois des outils imposés (code de calcul) dans certains Cahier
Des Charges (CDC). Sans conteste, le développement d’un logiciel à usage interne reste marginal et réservé aux fabricants d’échangeurs thermiques.

7. Conclusions

Le dimensionnement d’un échangeur, nous venons de le voir, est un processus complexe et souvent sujet à de nombreuses contraintes et hypothèses. Il suppose en effet de
préciser :

la modélisation de l’échangeur retenu ;

les propriétés thermophysiques des fluides ;
le choix technologique réalisé ;
les lois (corrélations) d’échange et de perte de pression spécifiques et adaptées ;
les contraintes d’intégration et de design (compacité, performance thermique et perte de charge acceptable) ;
les contraintes liées à l’utilisation de certains fluides (résistance d’encrassement, étanchéité, inspection);
le choix d’un logiciel de calculs appropriés aux besoins ;
le prix de l’échangeur (qui passe par sa conception, sa fabrication et sa mise en place).

À l’évidence, le fait de disposer d’un échangeur bien adapté, bien dimensionné, bien réalisé et bien utilisé permet un gain non négligeable de rendement et donc d’énergie dans
les process.

Notations et symboles

Symbole Unité Définition

A – Constante

C – Rapport des débits de capacités thermiques minimum sur maximum

($$C=frac{dot{C}_{min}}{dot{C}_{max}}$$)

$$dot{C}$$ W/K Débit de capacité thermique ($$dot{C}=qm. Cp$$

Cp J/kg-K Chaleur spécifique ou Chaleur massique ou Capacité thermique

D m Diamètre

Dh m Diamètre hydraulique ($$D_h=4.Sp/Pm$$(tube : Dh=D_int /

espace annulaire Dh=Dext-Dint)

e m Epaisseur d’une paroi

h W/m²-K Coefficient d’échange convectif partiel (parfois noté α)

H J/kg Enthalpie massique

L m Longueur

Lv J/kg Chaleur latente (ou enthalpie de changement de phase liquide-

vapeur)

Nu – Nombre de Nusselt

NUT – Nombre d’Unité de Transfert ($$NUT=U.S/C_{min}$$)

P W Puissance thermique
 Pm m Périmètre mouillé

Pr – Nombre de Prandtl ($$Pr=mu . C_p/lambda$$

qm kg/s Débit massique

Re – Nombre de Reynolds ($$Re=rho . V. D_h/mu$$

R K/W Resistance thermique

S m² Surface d’échange

Sp m² Section de passage

T K ou °C Température

U W/m²-K Coefficient d’échange global (parfois noté K)

x M Longueur

xv – Titre massique rapporté à la phase vapeur

Lettre grecque

Symbole Unité Définition

ɛ – Efficacité thermique

λ W/m-K Conductivité thermique

µ kg/m-s Viscosité dynamique

ρ kg/m3 Masse volumique

Δh J/kg Chaleur latente (ou enthalpie de changement de phase liquide-

vapeur)

ΔT K ou °C Différence de température

ΔTln K ou °C Ecart de température moyen logarithmique

Δxv – Différence du titre massique rapporté à la phase vapeur entre l’entrée

et la sortie du fluide

Indice :

Symbole Définition

c chaud

e entrée

ext externe

f froid

int interne

l liquide

min minimum

max maximum

s sortie

sat saturation (rapporté à l’état saturé d’un fluide)

 v vapeur

Notes et références

[1] BONTEMPS. A, GUARRIGUE.A, GOUBIER.C, HUETZ.J, MARVILLET.C, MERCIER.P, VIDIL.R. Echangeurs de chaleur, Techniques de l’ingénieur, traité Génie énergétique, Références
B2340, B2341, B2342, B2343, B2344.

[2] Dans cette notice, on se préoccupe des échangeurs de chaleur dit pariétal, c’est-à-dire ceux dont le fluide chaud et le fluide froid sont séparés par une paroi, généralement
métallique. L’autre catégorie d’échangeur dit à contact direct occupe une place importante dans certaines applications spécifiques et nous renvoyons le lecteur, pour plus
d’informations à ce sujet, aux références :

– BRICARD.A, TADRIST L (1996). Initiation aux échangeurs thermiques à contact direct. Collection « Echangeurs de chaleur », GRETh & IUSTI, Lavoisier TEC&DOC, Juin 1996.

–  BRICARD A, TADRIST  L( 1999). Echangeurs de chaleur à contact direct. Techniques de l’ingénieur, traité Génie énergétique, Référence BE9565, Octobre .

Tout au long de la lecture de cette notice, le lecteur pourra se reporter à la nomenclature et aux références en fin de notice.

[3] MARVILLET.C, MANIFICAT.A, BERRADA.N et WEBER.C (2014). E.Book TOME 1: Echangeurs Thermiques -: Technologies, Conception et dimensionnement, GRETh.

[4] Pour plus d’informations sur les caloducs, le lecteur peut se référer à :

–  CHAUDOURNE.S (1987).  Les échangeurs à caloducs.  Collection « Echangeurs de chaleur », GRETh, Lavoisier TEC&DOC, Septembre.

–  BRICARD.A, CHAUDOURNE.S  (1997).   Caloducs. Techniques de l’ingénieur, traité Génie énergétique, Référence b9545, Avril.

–  REAY.D.A, KEW. P.A, McGLEN.R.J (2014). Heat Pipes – Theory, design and Applications, 6th edition, Elsevier, ISBN: 978-0-08-098266-3.

[5] Pour la détermination des pertes de charge par le calcul du facteur de friction, le lecteur pourra se reporter aux ouvrages

– MARVILLET.C, MANIFICAT.A, BERRADA.N et WEBER.C  (2014). E.Book TOME 1, op. cit.

– Manuel technique du GRETH, édition GRETh.

– ROHSENOW W.M, HARTNETT J.P, CHO Y.I (1998).  Hanbook of Heat Transfer.   3rd  Edition – McGraw-Hill, 1998, ISBN 0-07-053555-8

– Verein Deutscher Ingenieure (VDI) – Gesellschaft Verfahrenstechnik und Chemieingenieurwesen (GVC), Editor VDI Heat ATLAS 2nd Edition, Springer, 2010, ISBN 978-3-540-
77876-9

et surtout l’ouvrage de référence :

– Idel’cik (1986). Mémento des pertes de charge, éditions Eyrolles.

[6] Les principaux ouvrages de référence sont :

– Manuel technique du GRETH, édition GRETh

– ROHSENOW W.M, HARTNETT J.P, CHO Y (1998).  Hanbook of Heat, op. cit.

– Verein Deutscher Ingenieure (VDI), op. cit.

[7] Standards of the Tubular Exchanger Manufacturers Association (TEMA), 9th Edition, 2007.

[8] Pour plus d’informations sur l’encrassement, le lecteur peut se référer à :

–  MARVILLET.C, MANIFICAT.A, BERRADA.N et WEBER.C  (2014). E.Book TOME 1, op. cit.

–  BONTEMPS. A, GUARRIGUE.A, GOUBIER.C, HUETZ.J, MARVILLET.C, MERCIER.P, VIDIL.R. Echangeurs, op.cit.

–  H. DEMASLES, P. MERCIER, P. TOCHON, B. THONON.  « Guide de l’encrassement des échangeurs de chaleur » (2007) Editions GRETH – ISBN : 2-9502555-5-8

–  T.R. BOTT (1995).  Fouling of heat exchangers. Elsevier Science & Technology Books; April, ISBN: 04444821864

qui traitent le sujet de manière assez exhaustive.

[9] MARVILLET.C, MANIFICAT.A, BERRADA.N et WEBER.C  (2014). E.Book TOME 1, op. cit.

[10] MARVILLET.C, MANIFICAT.A, BERRADA.N et WEBER.C  (2014). E.Book TOME 1, op. cit.

[11] MARVILLET.C, MANIFICAT.A, BERRADA.N et WEBER.C  (2014). E.Book TOME 1, op. cit.

[12] SHAH R.K. (2003). Fundamentals of Heat Exchanger Design.   John Wiley & Sons.

[13] MARVILLET.C, MANIFICAT.A, BERRADA.N et WEBER.C  (2014). E.Book TOME 1, op. cit.

[14] BONTEMPS. A, GUARRIGUE.A, GOUBIER.C, HUETZ.J, MARVILLET.C, MERCIER.P, VIDIL.R. Echangeurs, op. cit.
[15] Notamment :

– Manuel technique du GRETH, édition GRETh

– ROHSENOW W.M, HARTNETT J.P, CHO Y (1998).  Hanbook of Heat, op. cit.

– Verein Deutscher Ingenieure (VDI), op. cit.

[16]

Bibliographie complémentaire

Site internet fournisseurs d’échangeur de chaleur (non exhaustif)

www.ciat.com www.date.fr www.howden.com
www.alfalaval.com www.deltaairplus.fr www.sierra.it
www.tranter.com www.spx.com/en/apv www.faviersetrem.fr
www.spirec.fr www.balcke-duerr.com www.luvata.com
www.barriquand.com www.mersen.com www.guentner.com
www.klingenburg.de www.cepic- www.sanhuamc.com
www.fivesgroup.com/FivesCryog carbomeca.com www.serck-aviation.com
enie www.luve.it www.chartindustries.com
www.linde-engineering.com www.kochheattransfer.co www.wieland-
www.alpema.org m thermalsolutions.com
www.morgana.fr www.swep.net www.valtimet.com
www.basetec.net www.trianon- www.wlv.com
www.airec.se echangeur.com www.heatric.com
www.britheat.com www.btt-nantes.com www.apiheattransfer.com
www.hrs-heatexchangers.com www.gea.com www.funke.de
www.vahterus.com www.gea-erge.fr www.occocoolers.co.uk
www.calgavin.com www.motacoolingsystem. www.brownfintube.com
com
www.sglgroup.com
www.atherm.com

Site internet fournisseurs de logiciel (non exhaustif)

www.htri.net
htfs.aspentech.com/
www.aspentech.com/products/aspen-edr
www.prosim.net/en/software-prosec-11
www.greth.fr
www.nist.gov/el/building_environment/evapcond_software
Site internet Association (non exhaustif)
www.nist.gov
www.tema.org
www.alpema.org
www.greth.fr
 

L’Encyclopédie de l’Énergie est publiée par l’Association des Encyclopédies de l’Environnement et de l’Énergie (www.a3e.fr), contractuellement liée à l’université Grenoble
Alpes et à Grenoble INP, et parrainée par l’Académie des sciences.

Pour citer cet article, merci de mentionner le nom de l’auteur, le titre de l’article et son URL sur le site de l’Encyclopédie de l’Énergie.

Les articles de l’Encyclopédie de l’Énergie sont mis à disposition selon les termes de la licence Creative Commons Attribution – Pas d’Utilisation Commerciale – Pas de
Modification 4.0 International.