UNIVERSITE MOULOUD MAMMERI TIZI-OUZOU
FACULTE DU GENIE DE LA CONSTRUCTION
DEPARTEMENT DE GENIE MECANIQUE
Exercices sur le calcul de l’efficacité par la méthode NUT
Exercice 1 :
De l’eau froide circule dans un tube de chaudière à condensation. Sa température d’entrée est
Tfe 18 C et son débit massique qmf 400 kg/h. Le réchauffage de cette eau est assuré par
condensation de la vapeur d’eau à l’extérieur du tube, à la température Tc 104 C.
On donne : diamètre intérieur d 12,5 mm; diamètre extérieur D 16 mm; longueur L 2,4
m; conductivité de la paroi p 46 W/m.K. Pour l’eau, on admettra dans la gamme de
température considérée :
f = 0,7 10-6 m2/s; Prf = 5,5 Cpf 4180 J / kg.K.
1. Calculer le coefficient d’échange hf à l’intérieur du tube. Pour cela,
a) Calculer le nombre de Reynolds et définir la nature du régime d’écoulement.
b) Calculer hf
2. On donne le coefficient h côté vapeur : hc =8000 W / m2 K
Calculer le coefficient global d’échange U (W/ m2.K).
3. Calculer le NUT et l’efficacité E de l’appareil.
4. Déterminer la température de sortie d’eau Tfs, puis la quantité de chaleur Q récupérée
annuellement grâce au dispositif à condensation, si l’on considère que la saison de chauffe
dure 150 jours et que la chaudière fonctionne 5 heures par jour.
Pour un régime turbulent avec Pr 0,66 ; Re 10 5 ; L / d 60
On utilise les corrélations suivantes (Dittus – Boelter)
Fluide froid : 𝑺𝒕𝒇 = 𝟎. 𝟎𝟐𝟑 𝑹𝒆−𝟎.𝟐 𝑷𝒓−𝟎.𝟔
Fluide chaud : 𝑺𝒕𝒄 = 𝟎. 𝟎𝟐𝟑 𝑹𝒆−𝟎.𝟐 𝑷𝒓−𝟎.𝟕
Re et Pr sont évalués à la température moyenne du film.
Le nombre de Stanton Reynolds, Prandtl sont donnés par les relations :
St h / Cp.V ; Re=V. d / Pr = .Cp/
Stanton peut s’écrire également en fonction de Nusselt et Peclet:
St Nu / Pe avec : Pe =Pr.Re et Nu h.d /
Exercice 2 :
Un échangeur de chaleur fonctionne à contre-courant dans les conditions suivantes :
Fluide chaud : Fluide chaud :
Entrée : Tce =350 C Entrée : Tfe =120 C
Sortie : Tcs 200 C Sortie :Tfs 290 C
Le débit thermique minimal étant le débit thermique du fluide froid (Cmin = Cf =mf .Cpf ) et la
puissance thermique échangée entre les deux fluides est estimée à : = 415 kW.
On demande :
1. Quelle est la puissance échangée si on fait travailler cet échangeur en mode co-courant,
avec les mêmes températures d’entrée et les mêmes débits ? (Utiliser la méthode NUT).
2. Quelles sont alors les nouvelles températures de sortie ?
Exercice 3 :
Un échangeur bitube à contre-courant de 1,5 m de long, de 2 cm de diamètre et de faible
épaisseur est utilisé pour refroidir un débit d’air de 9,4 kg/h de la température de 616 °C
jusqu’à une température de 178 °C. L’air circule dans le tube central et l’eau dans l’espace
annulaire.
1. Calculer la puissance calorifique à évacuer. On donne pour l’air : Cpc = 1,060 kJ/ kg.K
2. Le fluide de refroidissement est de l’eau, qui pénètre dans la section annulaire à la
température de 16 °C avec un débit de 0.6 l/mn. Calculer la température de l’eau à la sortie de
l’échangeur. On prendra Cpf = 4,180 kJ/ kg.K
3. Calculer le coefficient d’échange hc côté fluide chaud.
4. Déterminer l’efficacité E de cet échangeur, puis le nombre d’unité de transfert NUT. En
déduire le coefficient d’échange global U puis le coefficient d’échange hf côté fluide froid.
INDICATIONS
Nombre de Stanton : St h / Cp.V
Nombre de Reynolds : Re=V. d / V. d /
Nombre de Prandtl : Pr = .Cp/
Nombre de Péclet : Pe =Pr.Re
Nombre de Nusselt : Nu h d / Nu = St. Re. Pr
Stanton peut s’écrire également en fonction de Nusselt et Peclet: St Nu / Pe
Pour des tubes circulaires, on utilise la corrélation de Hausen :
2200 Re 10000 ; Pr 0,66
0,14
0,116 2
−
2 𝑑 2/3 𝜇
𝑆𝑡 = ( 𝑅𝑒3 −125) 𝑃𝑟 3 [1 + ( ) ] ( )
𝑅𝑒 𝐿 𝜇𝑝
𝜇 ∶ 𝑣𝑖𝑠𝑐𝑜𝑠𝑖𝑡é à 𝑙𝑎 𝑡𝑒𝑚𝑝é𝑟𝑎𝑡𝑢𝑟𝑒 𝑚𝑜𝑦𝑒𝑛𝑛𝑒 𝑑𝑢 𝑚é𝑙𝑎𝑛𝑔𝑒
𝜇𝑝 ∶ 𝑣𝑖𝑠𝑐𝑜𝑠𝑖𝑡é à 𝑙𝑎 𝑡𝑒𝑚𝑝é𝑟𝑎𝑡𝑢𝑟𝑒 𝑚𝑜𝑦𝑒𝑛𝑛𝑒 𝑑𝑒 𝑝𝑎𝑟𝑜𝑖
Dans le cas qui nous concerne on peut supposer que la viscosité reste constante, le
0,14
𝜇
rapport (𝜇 ) =1
𝑝
𝑑
Le rapport d/L = 2/150 est faible, le terme [1 + ( 𝐿 )2/3 ] 𝑝𝑒𝑢𝑡 ê𝑡𝑟𝑒 𝑝𝑟𝑖𝑠 é𝑔𝑎𝑙 à 1.
--------------------------------------------
L’échangeur évolue à fluide isotherme (le fluide chaud évolue à température constante :
changement de phase pour la condensation de la vapeur). On a donc dans ce cas :
𝐶
Ctmin = Ctf et 𝑅 = 𝐶 𝑡𝑚𝑖𝑛 = 0
𝑡𝑚𝑎𝑥
Type de Efficacité NUT Cas particulier Elim
circulation E ( NUT, R) NUT (E,R) pour NUTinfini
Contre- 1 − exp[−(1 − 𝑅)𝑁𝑈𝑇] 1 1 − 𝐸𝑅 R=0, E=1-exp[-NUT] 𝐸𝑙𝑖𝑚 = 1
𝐸= 𝑁𝑈𝑇 = ln( ) 1
courant 1 − 𝑅. exp[−(1 − 𝑅)𝑁𝑈𝑇] 1−𝑅 1−𝐸 NUT =ln( )
1−𝐸
Co-courant 1 1 R=0, E=1-exp[-NUT] 1
𝐸= [1 − exp[−(1 + 𝑅)𝑁𝑈𝑇]] 𝑁𝑈𝑇 = − ln[1 − 𝐸(1 + 𝑅)] 1 𝐸𝑙𝑖𝑚 =
1+𝑅 1+𝑅 NUT =ln( ) 1+𝑅
1−𝐸
Solution exercice 1
Il faut d’abord savoir quel est le régime d’écoulement dans les tubes. Pour cela, calculons le
𝑞
nombre de Reynolds Re : 𝑅𝑒 = 𝑚𝑓
𝜌𝑓 .𝑆
𝑽 = 𝟎. 𝟗 𝒎/𝒔
4 𝑞𝑣𝑓 4 𝑥 400 10−3
𝑉 = 𝜋𝑑2 = 3600 𝑥 3.14 𝑥 (12.5 10−3)2 = 0.9 𝑚/𝑠
𝑉.𝑑 𝑞 𝑑 𝑞𝑣𝑓 𝑑 4 𝑞𝑣𝑓 𝑑 4 𝑞𝑣𝑓
𝑅𝑒 = = 𝜌𝑚𝑓 = = =
𝜐 .𝑆 𝑓 𝜐 𝑆 𝜐 𝜋𝑑2 𝜐 𝜋𝑑𝜐
𝑅𝑒 = 16 176 𝑹𝒆 = 𝟏𝟔 𝟏𝟕𝟔
Avec cette valeur de Reynolds, pour L/d=2400/12.5 >60 et Prf>0.66, la corrélation qui donne
Stanton est :
𝑆𝑡𝑓 = 0.023 𝑅𝑒 −0,2 𝑃𝑟 −0,6 = 𝑆𝑡𝑓 = 0.023 𝑥 16176−0,2 5,5−0,6 =
0,023 𝑥 0,144 𝑥 0,359 = 1,19 10−3 𝑺𝒕 = 𝟏, 𝟏𝟗 𝟏𝟎−𝟑 𝒇
ℎ𝑓
Avec 𝑆𝑡𝑓 = 𝜌 ce qui permet de déduire hf :
𝑓 .𝐶𝑝𝑓.𝑉
4 𝑞𝑣 4 𝑞𝑣
ℎ𝑓 = 𝑆𝑡𝑓 . 𝜌𝑓 . 𝐶𝑝𝑓 . 𝑉 = 𝑆𝑡𝑓 . 𝜌𝑓 . 𝐶𝑝𝑓 2
= 0.023 𝑅𝑒 −0,2 𝑃𝑟 −0,6 𝜌𝑓 . 𝐶𝑝𝑓 .
𝜋𝑑 𝜋 𝑑2
−3
4 𝑥 400 10
= 1,19 10−3 𝑥 1000 𝑥 4180. = 4506 𝑊/𝑚2 . 𝐾
3600𝑥 3.14 𝑥 (12.5 10−3 )2
𝒉𝒇 = 𝟒𝟓𝟎𝟔 𝑾/𝒎𝟐 . 𝑲
Calcul du coefficient global d’échange U
L’épaisseur du tube est :
𝐷−𝑑 16−12,5
𝑒= 2 = = 1.25 𝑚𝑚
2
L’épaisseur étant faible devant le diamètre, on raison comme une paroi plane
La résistance thermique totale est :
1 𝑒 1 1 1,75 10−3 1
𝑅𝑡ℎ = ℎ + 𝜆 + ℎ = 8000 + + 4506 = 0,000125 + 0,000038 + 0,000222 =
𝑐 𝑓 46
3,85 10−4
Le coefficient de transfert global est : 𝑈 =
1 U = 2597 𝑾/𝒎𝟐 . 𝑲
𝑅𝑡ℎ
2. Calcul du nombre d’unité de transfert NUT et de l’efficacité E
L’échangeur est à fluide isotherme (le fluide chaud évolue est à température constante :
changement de phase)
Calcul de NUT :
𝑈 .𝐴
𝑁𝑈𝑇 = C A est la surface moyenne d’échange et Ctf = Cmin est le débit de capacité
tf
thermique minimal. 𝐴 = 0,107 𝒎𝟐
(16+12,5)10−3 𝑥 2.4
𝐴 = 3.14 𝑥 = 0,107 𝑚2
2
Le débit thermique de l’eau a pour valeur :
400
Ctf = 𝑚̇𝑓 𝐶𝑝𝑓 = 3600 4180 = 465 𝑊/𝐾 𝐂𝐭𝐟 = 𝟒𝟔𝟓 𝑾/𝑲
Et le nombre d’unité de transfert NUT vaut : 𝒎𝟐
𝑈 .𝐴 2597 𝑥 0,107
𝑁𝑈𝑇 = C = = 0,60 𝑵𝑼𝑻 = 𝟎, 𝟔𝟎
tf 465
L’efficacité de l’échangeur co-courant en changement de phase est donnée pour R=0 dans le
tableau ci-dessous :
𝐶𝑡𝑚𝑖𝑛
Ctmin = Ctf et 𝑅 = =0 𝑅=0
𝐶𝑡𝑚𝑎𝑥
R=0, E=1-exp[-NUT]
E= 1- exp (-0,6) = 1- 0,548 = 0,452 𝑬 = 𝟎, 𝟒𝟓𝟐
3. Calcul de la température de sortie du fluide froid
Pour le fluide froid, l’efficacité peut s’écrire en fonction des températures :
𝑇𝑓𝑠 −𝑇𝑓𝑒
𝐸= soit, 𝑇𝑓𝑠 = 𝐸(𝑇𝑐 − 𝑇𝑓𝑒 )+ 𝑇𝑓𝑒 = 0,452 x (104-18)+18= 38,9+18=56,9 =57°C
𝑇𝑐 −𝑇𝑓𝑒
𝑻𝒇𝒔 = 57 °C
La puissance thermique du système s’exprime par le bilan thermique du fluide froid :
Ctf (Tfs Tfe ) 465 (57 18 ) 18135 W
18 kW 18 kW
Sur une saison de chauffe qui compte 150 5 heures, on peut donc récupérer par
condensation au maximum :
Q 18 150 5 soit Q 13500 kWh Q = 13500 kWh
Commentaires
Il s’agit encore d’un problème allégé où on ne fait mention ni des températures de référence,
ni de la distinction entre coefficient d’échange global côté froid et côté chaud, et où le calcul
de h côté vapeur n’est pas demandé.
La température de condensation de 104 °C sous-entend que la vapeur est en légère surpression
par rapport à la pression atmosphérique.
Les questions 3 et 4 peuvent aussi être traitées avec les résultats relatifs aux échangeurs à
fluide isotherme.
Type de Efficacité NUT Cas particulier Elim
circulation E ( NUT, R) NUT (E,R) pour NUTinfini
Contre- 1 − exp[−(1 − 𝑅)𝑁𝑈𝑇] 1 1 − 𝐸𝑅 R=0, E=1-exp[-NUT] 𝐸𝑙𝑖𝑚 = 1
courant 𝐸= 𝑁𝑈𝑇 = ln( ) 1
)
1 − 𝑅. exp[−(1 − 𝑅)𝑁𝑈𝑇] 1−𝑅 1−𝐸 NUT =ln(
1−𝐸
Co-courant 1 1 R=0, E=1-exp[-NUT] 1
𝐸= [1 − exp[−(1 + 𝑅)𝑁𝑈𝑇]] 𝑁𝑈𝑇 = − ln[1 − 𝐸(1 + 𝑅)] 1 𝐸𝑙𝑖𝑚 =
1+𝑅 1+𝑅 NUT =ln( ) 1+𝑅
1−𝐸
Solution exo 2 :
1.
En contre-courant, avec Ctmin Ctf , on a pour efficacité :
𝑇𝑓𝑠 − 𝑇𝑓𝑒 290 − 120
𝐸= = = 0,739 = 0,74
𝑇𝑐𝑒 − 𝑇𝑓𝑒 350 − 120
Le rapport R des débits de capacité thermique est :
𝑇𝑐𝑒 − 𝑇𝑐𝑠 350 − 200
𝑅= = = 0,882
𝑇𝑓𝑠 − 𝑇𝑓𝑒 290 − 120
et, à partir de l’expression du nombre d’unité de transfert pour l’échange en contre-courant:
1 (1 − 𝑅. 𝐸) 1 1 − 0,882𝑥0,74
𝑁𝑈𝑇 = 𝑙𝑛 = 𝑙𝑛 = 8,475 ln(1,33585) = 2,44
1−𝑅 (1 − 𝐸) 1 − 0,882 1 − 0,74
En circulation co-courant, les débits n’étant pas modifiés, les coefficients d’échange ne le sont
pas non plus. On garde donc le même NUT (vu que NUT U.S/Ctmin ). Par contre, la nouvelle
efficacité E’ de cet échangeur en mode co-courant change et s’écrit:
𝐸 ′ = 1+𝑅 {1 − 𝑒𝑥𝑝[1 − (1 + 𝑅)𝑁𝑈𝑇]} = 1+0,882 (1 − 0,027027)
1 1
Il vient, après remplacement de R et NUT par leur valeur :
E' 0,517
Puisque les conditions d’entrée sont identiques dans les deux cas, la nouvelle
puissance ' est telle que :
∅′ ∅′
′
𝐸 = = 𝐶 (𝑇 − 𝑇𝑓𝑒 )
∅𝑚𝑎𝑥 𝐶𝑡𝑚𝑖𝑛 (𝑇𝑐𝑒 − 𝑇𝑓𝑒 ) 𝑡𝑚𝑖𝑛 𝑐𝑒
∅ ∅ 𝐸′ 𝐸′ ∅′
de plus : telle que : 𝐸 = ∅ =𝐶 et le rapport donne : =
𝑚𝑎𝑥 𝑡𝑚𝑖𝑛 (𝑇𝑐𝑒 −𝑇𝑓𝑒 ) 𝐸 𝐸 ∅
ce qui donne pour la nouvelle valeur du flux :
𝐸′ 0,517
∅′ = ∅ = 415 𝑥 = 290 𝑘𝑊
𝐸 0,74
' 290 kW
c’est-à-dire 70% de la puissance en contre-courant.
2.
′
La nouvelle température de sortie du fluide froid 𝑇𝑓𝑠 s’obtient à partir de la nouvelle
efficacité :
𝑇 ′ −𝑇𝑓𝑒
𝐸 ′ = 𝑇𝑓𝑠 −𝑇
𝑐𝑒 𝑓𝑒
𝑇𝑓𝑠
′
𝑇𝑓𝑒 + 𝐸 ′ (𝑇𝑐𝑒 − 𝑇𝑓𝑒 ) 120 +0,517x ( 350 -120 )
′
𝑇𝑓𝑠 C
et la nouvelle température de sortie chaude du fluide chaud 𝑇𝑐𝑠′ s’obtient à partir du rapport
𝑇 −𝑇 ′
𝑇𝑐𝑒 − 𝑇𝑐𝑠
𝑅 = 𝑇 𝑐𝑒−𝑇 𝑐𝑠 = ′ −𝑇 = 0,882 qui reste inchangé.
𝑓𝑠 𝑓𝑒 𝑇𝑓𝑠 𝑓𝑒
′
𝑇𝑐𝑠′ s’écrit ′
alors : 𝑇𝑐𝑠 = 𝑇𝑐𝑒 − 𝑅 (𝑇𝑓𝑠 − 𝑇𝑐𝑒 ) = 350 − 0,882 (240 − 120) = 244°𝐶
Commentaires
Il y a d’autres façons de résoudre ce petit exercice mais la démarche préconisée permet de se
roder à la méthode NUT dans un cas simple.
Avec cette valeur de R, on est déjà dans la zone asymptotique d’un échangeur co-courant, et
les températures de sortie sont donc très voisines.
Solution Exo 3
1. Le flux total peut se calculer côté chaud :
∅ = 𝐶𝑡𝑐 . (𝑇𝑐𝑒 − 𝑇𝑐𝑠 )
D’après les données, le débit thermique unitaire 𝐶𝑡𝑐 côté fluide chaud est :
9.4
𝐶𝑡𝑐 = 𝑚̇𝑐 . 𝐶𝑝𝑐 = 3600 . 1060 = 2.77 𝑊/𝐾
Ce qui donne pour le flux total :
∅ = 2.77 (616 − 178) = 1213 𝑊 ∅ = 𝟏𝟐𝟏𝟑 𝑾
2. Le calcul du flux total côté froid va nous permettre de calculer la température de sortie du
fluide froid 𝑇𝑓𝑠
∅ = 𝐶𝑡𝑓 . (𝑇𝑓𝑠 − 𝑇𝑓𝑒 )
0.6
avec : 𝐶𝑡𝑓 = 𝑚𝑓 . 𝐶𝑝𝑓 = . 4180 = 41.8 𝑊/𝐾
60
Ce qui donne pour la température 𝑇𝑓𝑠 :
∅ 1213 𝑇𝑓𝑠 = 𝟒𝟓 °𝑪
𝑇𝑓𝑠 = 𝑇𝑓𝑒 + 𝐶 = 16 + = 45 °𝐶
𝑡𝑓 41.8
3. Pour calculer le coefficient d’échange côté fluide chaud (air), il faut d’abord connaître le
régime d’écoulement de l’air, donc le nombre de Reynolds côté chaud.
On doit cependant calculer approximativement les caractéristiques moyennes de l’air, à
savoir , la température moyenne de l’air est approximativement :
𝑇 +𝑇 616+178
𝑇𝑐𝑚𝑜𝑦 = 𝑐𝑒 2 𝑐𝑠 = = 397 °𝐶 𝑠𝑜𝑖𝑡 670 𝐾
2
A cette température, les tables des caractéristiques de l’air donnent respectivement pour la
masse volumique, la viscosité cinématique et la conductivité thermique, les valeurs suivantes :
𝜌𝑐 = 0.525 𝑘𝑔/𝑚3, 𝜐𝑐 = 6.20 10−5 𝑚2 /𝑠 , 𝜆𝑐 = 0.0505 𝑊/𝑚. 𝐾
La section du tube est :
𝑑2 (2.10−2 )2
𝑠𝑐 = 𝜋 = 3.14 = 3.14 10−4 𝑚2 𝒔𝒄 = 𝟑. 𝟏𝟒 𝟏𝟎−𝟒 𝒎𝟐
4 4
On en déduit la vitesse débitante :
𝑞𝑣𝑐 𝑞𝑚𝑐 9.4
𝑉𝑐 = = = = 15.8 𝑚/𝑠
𝑠𝑐 𝜌𝑐 . 𝑠𝑐 3600 𝑥 0.525 𝑥 3.14 10−4
𝑽𝒄 = 𝟏𝟓. 𝟖 𝒎/𝒔
d’où le nombre de Reynolds :
𝑉𝑐 .𝑑 15.8 𝑥 2 10−2
𝑅𝑒 = = = 5097 𝑅𝑒𝒄 = 𝟓𝟎𝟗𝟕
𝜐𝑐 6.20 10−5
Cette valeur de Reynolds montre que le régime d’écoulement est transitoire (Re > 5000). Le
rapport d/L= 2/150 est négligeable, 0n peut donc utiliser la corrélation suivante :
0,116 2 2
𝑆𝑡𝑐 = ( 𝑅𝑒𝑐 3 −125) 𝑃𝑟𝑐 −3
𝑅𝑒𝑐
Pour le calcul de hc il est plus simple de le déduire de Nuc avec :
1 2
𝑁𝑢𝑐 = 𝑆𝑡𝑐 𝑅𝑒𝑐 𝑃𝑟𝑐 = 0,116 ( 𝑅𝑒𝑐 2/3 −125) 𝑃𝑟𝑐 3 = 0.116 (𝑅𝑒𝑐 3 − 125) 𝑅𝑒𝑐 1/3
A la température moyenne de 670 K la valeur de Prandtl suivante :
𝜌𝑐 𝑥 𝜐𝑐 𝑥 𝐶𝑝𝑐 0.525 𝑥 6.20 10−5 𝑥 1060
𝑃𝑟𝑐 = = = 0,68 𝑷𝒓𝒄 = 𝟎. 𝟔𝟖
𝜆𝑐 0,0505
D’où pour le nombre de Nusselt :
2
𝑁𝑢𝑐 = 0.116 (50973 − 125) 0.681/3 = 0.116 (297 − 125)𝑥0.88 = 17.47
𝑵𝒖𝒄 = 𝟏𝟕. 𝟒𝟕
Or le coefficient d’échange peut être déduit de la relation de Nusselt qui s’écrit:
ℎ 𝑥𝑑
𝑁𝑢𝑐 = 𝑐𝜆
𝑐
𝜆𝑐 𝑥 𝑁𝑢𝑐 0.0505 𝑥 17.47
soit pour ℎ𝑐 ∶ ℎ𝑐 = = = 44 𝑊/𝑚2 . 𝐾
𝑑 2.10−2
𝒉𝒄 = 𝟒𝟒 𝑾/𝒎𝟐 . 𝑲
4.
L’efficacité de cet échangeur se calcule par la connaissance des températures d’entrée et
sortie des deux fluides et du débit thermique unitaire minimal 𝐶𝑡𝑚𝑖𝑛 :
𝐶𝑡𝑐 = 𝑀𝑐 . 𝐶𝑝𝑐 = 2.77 𝑊/𝐾
𝐶𝑡𝑓 = 𝑚𝑓 . 𝐶𝑝𝑓 = 41.8 𝑊/𝐾
Soit 𝐶𝑡𝑚𝑖𝑛 = 𝐶𝑡𝑐 = 2.77 𝑊/𝐾
L’efficacité E s’écrit alors :
𝐹𝑙𝑢𝑥 𝑟é𝑒𝑙 𝐶𝑡𝑚𝑖𝑛 .(𝑇𝑐𝑒 −𝑇𝑐𝑠 ) (𝑇 −𝑇 )
𝐸= = = (𝑇𝑐𝑒−𝑇𝑐𝑠 )
𝑓𝑙𝑢𝑥 𝑚𝑎𝑥 𝐶𝑡𝑚𝑖𝑛 .(𝑇𝑐𝑒 −𝑇𝑓𝑒 ) 𝑐𝑒 𝑓𝑒
(616−178) 438
𝐸= = (600) = 0.73 𝑬 = 𝟎. 𝟕𝟑
(616−16)
Pour un échangeur contre courant, le nombre d’unité de transfert NUT est donné par la
relation suivante :
1 1−𝐸𝑅 𝐶 2.77
𝑁𝑈𝑇 = 1−𝑅 𝐿𝑛 1−𝐸 avec 𝑅 = 𝐶 𝑡𝑚𝑖𝑛 = 41.8 = 0.066 𝑹 = 𝟎. 𝟎𝟔𝟔
𝑡𝑚𝑎𝑥
Ce qui donne :
1 1−0.73𝑥0.066 0.9518
𝑁𝑈𝑇 = 1−0.066 𝐿𝑛 1−0.73 = 1.07𝑥𝐿𝑛 0.27 = 1.07𝑥𝐿𝑛 3.525 = 1.35
𝑵𝑼𝑻 = 𝟏. 𝟑𝟓
De la définition de NUT, on peut tirer la valeur du coefficient global k (W/m2.K) :
𝑁𝑈𝑇. 𝐶𝑡𝑚𝑖𝑛
𝑘= 𝑜ù 𝐴 𝑒𝑠𝑡 𝑙𝑎 𝑠𝑢𝑟𝑓𝑎𝑐𝑒 𝑑 ′ é𝑐ℎ𝑎𝑛𝑔𝑒
𝐴
Puisque l’épaisseur du tube central est faible, on ne fait pas la distinction entre surface
d’échange côté chaud et côté froid, et on néglige la résistance thermique de la paroi. Donc :
𝐴 = 𝜋. 𝑑. 𝐿 = 3.14 𝑥 2 10−2 𝑥 1.5 = 0.094 m2 𝑨 = 𝟎. 𝟎𝟗𝟒 𝐦𝟐
Et la valeur de U est alors :
1.35 𝑥 2.77
𝑈 = 0.094 = 39.78 𝑊. 𝑚−2 . 𝐾 −1 𝑼 = 𝟑𝟗. 𝟕𝟖 𝑾. 𝒎−𝟐 . 𝑲−𝟏
Le coefficient d’échange global s’exprime aussi à partir de l’expression de la résistance
thermique surfacique globale (en négligeant l’épaisseur e du tube central ainsi que la
résistance par conduction dans ce tube) :
1 1 1
= ℎ +ℎ
𝑈 𝑐 𝑓
Et ℎ𝑓 se déduit de l’expression ci-dessus :
1 1 1
ℎ𝑓 = 1 1 = 1 1 = (0.0251−0.0227) = 416 𝑊. 𝑚−2 . 𝐾 −1
( − ) ( − )
𝑈 ℎ𝑐 39.78 44
𝒉𝒇 = 𝟒𝟏𝟔 𝑾. 𝒎−𝟐 . 𝑲−𝟏