Dimensionnement des échangeurs thermiques industriels : méthodes, calculs et bonnes pratiques

💡 L'essentiel à retenir :

• Le programme thermique (débits, températures d'entrée/sortie, puissance Q) constitue le point de départ de tout dimensionnement d'échangeur de chaleur industriel.
• Deux méthodes structurent le calcul de surface : LMTD avec facteur de correction F (températures de sortie connues) et ε-NTU (sorties inconnues ou cas de changement de phase).
• Le coefficient global U se calcule comme une résistance en série : convection interne, conduction de la paroi, convection externe et encrassement des deux côtés.
• L'encrassement dégrade U et augmente ΔP en service ; le dimensionnement intègre une marge via les facteurs Rf pour tenir compte de cette dérive.
• Le dimensionnement thermique et hydraulique se boucle par itération : augmenter la vitesse améliore h mais accroît ΔP, ce qui peut dépasser le budget disponible côté pompe.
• La validation finale couvre le pincement thermique, les risques de croisement de température, les contraintes mécaniques (dilatation, vibrations) et la nettoyabilité.

Le dimensionnement d'un échangeur thermique industriel va bien au-delà du simple calcul de surface. Cela suppose de verrouiller un programme thermique cohérent, d'estimer un coefficient global U traçable, d'intégrer l'encrassement dès la phase de conception, et de vérifier que les pertes de charge restent compatibles avec le budget disponible côté procédé. Qu'il s'agisse d'un échangeur thermique air, d'un échangeur de chaleur à plaques ou d'un échangeur tubulaire, un dimensionnement adapté influe sur la performance énergétique et les gains.

Un dimensionnement rigoureux d'un échangeur thermique commence par la collecte exhaustive des données du procédé. Une donnée manquante ou approximative se traduit invariablement par une surface sous-estimée ou une perte de charge non maîtrisée.
Le tableau suivant présente les grandeurs à consolider avant tout calcul.

Aucun calcul de puissance d'échangeur thermique fiable ne peut démarrer sans au minimum les débits, les Cp, les températures d'entrée et de sortie, et une limite sur ΔP. Les propriétés à la température de film (moyenne entre paroi et fluide) influencent directement les corrélations de h utilisées ultérieurement.

Le devoir thermique Q exprime la puissance à transférer entre les deux fluides. Pour un fluide monophasique, Q s'écrit : Q = ṁ · Cp · ΔT
avec ṁ le débit massique (kg/s), Cp la capacité calorifique (J/kg·K) et ΔT la différence de température (K) entre entrée et sortie du fluide considéré.
Le bilan se ferme en écrivant Q côté chaud = Q côté froid (en négligeant les pertes vers l'environnement, hypothèse usuelle en conception industrielle). Si les quatre températures (deux entrées, deux sorties) sont toutes connues, le système est surdéterminé ; il faut alors vérifier la cohérence des débits. Si une sortie est inconnue, elle se déduit du Q calculé sur l'autre fluide.
En présence d'un changement de phase (condensation d'un fluide par exemple), le Q se décompose : une partie sensible liée au refroidissement vers la température de saturation, et une partie latente Q = ṁ · L. Cette distinction est déterminante pour la méthode ε-NTU traitée plus loin. Le bilan doit également tenir compte du sous-refroidissement éventuel côté condensat.

La méthode LMTD (Logarithmic Mean Temperature Difference) s'applique lorsque les quatre températures terminales sont connues. Elle repose sur la relation centrale :
Q = U · A · ΔTlm
où U est le coefficient global d'échange (W/m²·K), A la surface d'échange (m²), et ΔTlm la différence de température logarithmique moyenne.

En contre-courant (configuration à privilégier, car elle maximise la force motrice thermique) :

• ΔT₁ = Tchaud,entrée − Tfroid,sortie
• ΔT₂ = Tchaud,sortie − Tfroid,entrée
• ΔTlm = (ΔT₁ − ΔT₂) / ln(ΔT₁ / ΔT₂)

En co-courant (moins efficace) :

• ΔT₁ = Tchaud,entrée − Tfroid,entrée
• ΔT₂ = Tchaud,sortie − Tfroid,sortie

La formule reste identique, mais ΔTlm est plus faible, ce qui conduit à une surface plus grande pour le même Q. La configuration co-courant se justifie uniquement lorsque l'une des sorties doit être limitée en température (risque de dégradation d'un fluide sensible, par exemple).
Une fois ΔTlm calculé, la surface théorique s'obtient par : A = Q / (U · ΔTlm)
> À retenir : La LMTD suppose un coefficient U constant sur toute la surface, ce qui est acceptable pour des fluides monophasiques à propriétés peu variables. En présence de forts gradients de viscosité (huiles lourdes) ou de changement de phase partiel, la méthode devient moins précise.
> Erreurs fréquentes : Confondre ΔT₁ et ΔT₂ selon le sens d'écoulement ; omettre la correction par facteur F pour les configurations multi-passes ; utiliser ΔTlm contre-courant sur un échangeur en réalité co-courant.

La majorité des échangeurs de chaleur industriels ne réalisent pas un pur contre-courant idéal. Un échangeur calandre-et-faisceau avec deux passes côté tubes, ou un échangeur à plaques en configuration asymétrique, présente des écoulements croisés qui réduisent l'efficacité thermique.
Le facteur de correction F (sans dimension, compris entre 0 et 1) permet d'ajuster la LMTD idéale contre-courant pour tenir compte de la configuration réelle :
ΔTlm,corrigé = F · ΔTlm,contre-courant

F dépend de deux paramètres adimensionnels P et R, définis à partir des quatre températures terminales, et qui caractérisent respectivement l'efficacité côté froid et le rapport des capacités calorifiques.
Seuil d'alerte : lorsque F descend en dessous de 0,75, l'échangeur thermique approche d'une zone d'instabilité où une faible variation de température entraîne une très grande variation de surface. En pratique, F < 0,75 constitue un signal pour reconsidérer la configuration (ajouter des calandres en série, modifier le nombre de passes, ou opter pour une autre technologie).
Le cas des températures croisées, où la sortie froide dépasse la sortie chaude, est physiquement impossible en une seule calandre un seul passage ; il impose impérativement un montage en série (plusieurs calandres).

La méthode ε-NTU (efficacité – Nombre d'Unités de Transfert) complète la LMTD dans les situations où une ou plusieurs températures de sortie sont inconnues a priori.
Les grandeurs clés sont les suivantes :

• Cmin et Cmax : capacités calorifiques des deux fluides (C = ṁ · Cp, en W/K) ; Cmin est la plus faible.
• Cr = Cmin / Cmax : rapport de capacités (sans dimension, entre 0 et 1).
• NTU = U · A / Cmin : nombre d'unités de transfert, mesure de la taille thermique de l'échangeur.
• ε = Q / Qmax = Q / (Cmin · (Tchaud,entrée − Tfroid,entrée)) : efficacité, fraction de la puissance maximale théoriquement transférable qui est effectivement transférée.

La relation ε = f(NTU, Cr) dépend du type d'écoulement (contre-courant, co-courant, croisé). Pour un contre-courant pur, par exemple : ε = [1 − exp(−NTU·(1−Cr))] / [1 − Cr·exp(−NTU·(1−Cr))].

Lorsqu'un fluide change de phase (vapeur qui se condense, liquide qui s'évapore), sa capacité calorifique apparente tend vers l'infini : la température de ce fluide reste constante tout au long de l'échange. Dans ce cas, Cmax → ∞ et Cr = Cmin / Cmax → 0.
Cette simplification rend les formules ε-NTU particulièrement directes. Pour tous les types d'écoulement avec Cr = 0 : ε = 1 − exp(−NTU)
ce qui signifie que l'efficacité dépend uniquement du nombre d'unités de transfert, indépendamment de la configuration d'écoulement. Cette propriété est très utile en avant-projet pour estimer rapidement la surface nécessaire.
Le pincement thermique mérite une attention particulière dans ces cas : avec une condensation à température constante, la différence de température minimale se produit à l'extrémité froide de l'échangeur de chaleur (sortie côté liquide condensé). Un pincement trop faible (inférieur à 5 à 10 K selon les applications) génère une surface d'échange disproportionnée. Ce point conditionne également le choix technologique : un échangeur thermique à plaques tolère des pincements plus faibles qu'un calandre-et-faisceau, en raison de son fort coefficient d'échange.

Le coefficient global U se calcule en additionnant toutes les résistances au transfert thermique, de la même manière que des résistances électriques en série :
1/U = 1/hᵢ + Rf,ᵢ + e/λ + Rf,o + 1/hₒ

Le coefficient global U se calcule en additionnant toutes les résistances au transfert thermique, de la même manière que des résistances électriques en série :

Le coefficient de convection h se calcule via le nombre de Nusselt : h = Nu · k / Dₕ, où k est la conductivité thermique du fluide et Dₕ le diamètre hydraulique.
Le choix de la corrélation dépend du régime d'écoulement, identifié par le nombre de Reynolds Re = ρ · v · Dₕ / μ :

• Régime turbulent (Re > 10 000) : la corrélation de Dittus-Boelter (Nu = 0,023 · Re⁰·⁸ · Pr^n, avec n = 0,4 pour chauffage et 0,3 pour refroidissement) est fréquemment utilisée pour une première estimation dans des tubes lisses. Elle est valide pour 0,7 < Pr < 160 et Re > 10 000. La corrélation de Gnielinski offre une meilleure précision, notamment en zone de transition.
• Régime laminaire (Re < 2 300) : Nu tend vers une valeur asymptotique (Nu ≈ 3,66 pour flux thermique constant en tube, ou 4,36 pour température de paroi constante). Le coefficient h est alors indépendant de la vitesse, ce qui souligne l'intérêt de maintenir un régime turbulent.
• Zone de transition (2 300 < Re < 10 000) : les corrélations classiques ne s'appliquent pas directement ; la corrélation de Gnielinski étendue reste une option.

Pour les échangeurs de chaleur à plaques corrugués, les corrélations de type Nu = C · Re^m · Pr^(1/3) sont spécifiques à la géométrie (angle d'ondulation des plaques). Les valeurs de C et m varient significativement selon le fournisseur et la géométrie ; il est préférable d'utiliser les données du fabricant.
> Erreurs fréquentes : Appliquer Dittus-Boelter hors de son domaine de validité (Pr très élevé, tube court) ; négliger la correction de viscosité à la paroi (facteur μ/μparoi)^0,14 pour les fluides visqueux ; oublier que les corrélations côté calandre (autour des tubes avec chicanes) diffèrent des corrélations intra-tube.

L'encrassement est l'accumulation de dépôts (calcaire, biofilm, particules, coke, oxyde) sur les surfaces d'échange. Il augmente la résistance thermique totale et, simultanément, réduit la section de passage et donc augmente les pertes de charge.
Les facteurs Rf représentent cette résistance additionnelle en m²·K/W. Des références industrielles (généralement issues de standards comme ceux du TEMA pour les échangeurs thermiques calandre-et-faisceau) proposent des valeurs indicatives selon la nature des fluides :

• Eau douce traitée côté tubes : Rf ≈ 0,000 09 m²·K/W
• Eau de mer : Rf ≈ 0,000 09 à 0,000 18 m²·K/W
• Huiles légères : Rf ≈ 0,000 18 m²·K/W
• Fluides chargés ou procédés lourds : Rf pouvant atteindre 0,001 m²·K/W

L'impact sur U est direct : si U propre vaut 1 000 W/m²·K et que la somme des Rf atteint 0,0002 m²·K/W, le U en service devient environ 833 W/m²·K, soit une chute de 17 %. Cela se traduit par une surface A plus grande dès la conception pour maintenir Q.
Les stratégies de maîtrise de l'encrassement comprennent :

• Maintenir des vitesses minimales (au moins 1 à 1,5 m/s côté tube pour l'eau) pour éviter la sédimentation.
• Traiter les fluides en amont (filtration, adoucissement, dosage d'inhibiteurs).
• Concevoir l'échangeur thermique pour permettre un nettoyage chimique (CIP) ou mécanique (accès aux têtes côté tubes pour ramonage).
• Choisir des matériaux à faible adhérence ou des surfaces traitées pour les fluides fortement encrassants.

> À retenir : Un échangeur thermique surdimensionné de 15 à 25 % en surface par rapport au propre initial constitue une pratique courante pour absorber la dégradation de U due à l'encrassement sur une période entre deux nettoyages.

Les pertes de charge dans un échangeur de chaleur se décomposent en pertes régulières (frottement le long des tubes ou des canaux) et en pertes singulières (entrée/sortie, chicanes, virages de passes). La relation générale fait intervenir la vitesse v au carré, la densité ρ du fluide, et des facteurs de friction f liés à Re et à la rugosité de surface.
Le compromis thermique-hydraulique est au cœur du dimensionnement itératif :

• Augmenter la vitesse d'écoulement améliore Re, augmente h et donc U, ce qui réduit A.
• Mais cette même augmentation de vitesse accroît ΔP selon une loi approximativement quadratique (ΔP ∝ v²).
• Si le ΔP résultant dépasse le budget disponible côté procédé ou côté pompe, il faut réduire la vitesse, ce qui dégrade h et impose une surface plus grande.

La démarche itérative typique s'articule ainsi :

1. Hypothèse initiale sur U et la géométrie (diamètre de tube, nombre de tubes, nombre de passes).
2. Calcul de A via LMTD ou ε-NTU.
3. Estimation des vitesses à partir de A et de la géométrie choisie.
4. Calcul de ΔP côté tube et côté calandre.
5. Vérification des ΔP par rapport aux limites du procédé (typiquement 0,3 à 1,5 bar pour les liquides en service courant, 0,05 à 0,3 bar pour les gaz et vapeurs basse pression).
6. Ajustement de la géométrie (nombre de passes, diamètre de tube, espacement des chicanes) et retour à l'étape 1 jusqu'à convergence.

Les limites à surveiller en plus du ΔP sont l'érosion (vitesses excessives, particulièrement en présence de particules ou pour les fluides diphasiques) et les vibrations du faisceau de tubes (liées à la vitesse côté calandre et à la géométrie des chicanes).

Choisir le bon type d’échangeur thermique dépend directement du service, des contraintes de pression et température, du niveau d’encrassement et des exigences de maintenance.