Faculté des Sciences et des Sciences Appliquées, BOUIRA

Département Génie Mécanique : 2eme Année Master énergétique

Module: Echangeur de chaleur

Solution série TD N° 01 : Le coefficient global de transfert de chaleur

Exercice 1
Pour un échangeur mince (épaisseur négligeable), le coefficient global d’échange est
déterminé par la relation suivante :
1
U
1 1

hi he
- Calcul de hi :
On détermine le régime de l’écoulement dans le tube interne
La vitesse de l’écoulement de l’eau:
m m 0,5
v    1,61 m / s
S 1 2 1 2
   d  990   0,02 
4  4 
v  dh 1,61  0,02
Re    53488,37
 0, 602 106
Re >10  Régime est turbulent
4

Donc le nombre de Nusselt est calculé à partir de la corrélation suivante :

Nu  0, 023Re0.8 Pr n  0, 023(53488,37)0.8 (3,91)0,4  240,55

Alors :
hi  dh Nu  k 240,55  0,637
Nu   hi    7661,52 W / m2 C
k dh 0,02
- Calcul de he :
On détermine le régime de l’écoulement dans l’espace annulaire (entre les deux tubes)
La vitesse de l’écoulement de l’huile:
m m 0,8
v    2,39 m / s
S 1  1 
   ( D2  d 2 
852   (0,032  0,022 
 4   4 
v  dh v  (D-d) 2 ,39  0,01
Re     637,33
  37,5 106
Re <2000  Régime est laminaire
Les nombres de Nusselt pour un écoulement laminaire complètement développé dans un
espace annulaire avec une surface isotherme et l'autre adiabatique (Kays et Perkins) sont
donnés par le tableau suivant :

1
 Di/De Nui Nue
0 - 3.66
0.05 17.46 4.06
0.1 11.56 4.11
0.25 7.37 4.23
0.5 5.74 4.43
1 4.86 4.86

Par interpolation, le nombre du Nusselt sur la surface extérieure du tube interne,

correspondant à d/D= 0,67 est :
Nu  5, 45
Alors :
Nu  k 5,45  0,138
he    75, 21
Dd 0,01

Le coefficient d’échange global devient :

1 1
U   74, 49 W / m2 C
1 1 1 1
 
hi he 7661,52 75, 21

Pour améliorer le coefficient d’échange, on peut mettre des ailettes dans le côté où circule
l’huile

Exercice 2
1- la résistance thermique de l'échangeur de chaleur par unité de longueur.
Les surfaces d’échanges intérieure et extérieure :
Ai    di  L
Ae    de  L
On prend L= 1 m
1 R fi ln(d e / di ) 1 R fe
R    
hi Ai Ai 2 kL he Ae Ae
1 0,0005 ln(1,6 / 1, 2) 1 0,0002
    
700  (  0,012  1)   0,012  1 2    380  1 1400  (  0,016  1)   0,016  1
 0,0695 C / W

2- les coefficients globaux de transfert thermique Ui et Ue

1 1
Ui    381, 67 W / m2 C
RAi 0, 0695  (  0, 012 1)

1 1
Ue    286, 25 W / m2 C
RAe 0, 0695  (  0, 016 1)

2
Exercice 3
1- Le coefficient d’échange global de cet échangeur de chaleur.
On calcul he :
La vitesse de l’écoulement de l’eau dans l’espace annulaire:
m m 0,3
v    0,73 m / s
S 1  1 
   ( D2  d 2 
998   (0,0252  0,012 
 4   4 
v  dh 0,73  0,015
Re    10906,37
 1, 004 106
Re >104  Régime est turbulent
Donc le nombre de Nusselt est calculé à partir de la corrélation suivante :

Nu  0, 023Re0.8 Pr 0,4  0, 023(10906,37)0.8 (7, 01)0,4  85,145

Alors :
Nu  k 85,145  0,598
he    3354,58 W / m2 C
dh 0,015
Donc le coefficient global d’échange est égale à :
1 1
U   2007,63 W / m2 C
1 1 1 1
 
hi he 5000 3354,58
2-

1 1
U   1001,9 W / m2 C
1 1 1 1
  R fe   0,0005
hi he 5000 3354,58

Ce résultat démontre l’influence de l’encrassement sur la performance d’un échangeur de

chaleur, où on constate que l’encrassement diminue de la moitié le coefficient d’échange
global.

Exercice 4
1- Le coefficient d’échange global côté de l’eau froide Ui
1 1
Ui  
RAi 1  di ln(d e / di )  di
hi 2k he d e
1
  248,53W / m 2 C
1 0, 02 ln(26 / 20) 20
 
8000 2  50 200  26

2- Ailettes sur surface extérieure du tube

Efficacité global des ailettes :

A
o  1  f (1   f )
A
3
Dimension d’une ailette : (l, L, e) = (15mm, 1m, 2mm)
Surface des ailettes :
Af  N (2  l  L  e  L)  16  (2  0,015 1  0,002 1)  0,512m2
Surface d’échange totale (surface des ailettes + surface sans ailettes) :
A  Af  ( de L  N  e  L)  0,512  ((  0,026 1)  (16  0,002 1))  0,562m2
Af 0,512
o  1  (1   f )  1 
(1  0, 78)  0,80
A 0,562
1 1
Ui  
RAi 1  di ln(d e / di )   di L
hi 2k o he A
1
  1141, 23W / m 2 C
1 0, 02 ln(26 / 20)   0, 02 1
 
8000 2  50 0,8  200  0,562

En ajoutant des ailettes sur la surface extérieure (côté du faible coefficient de convection), on
constate une augmentation considérable du coefficient d’échange global U.

Exercice 5

1- La chaleur échangée entre la vapeur et l’eau circulant dans un seul tube.

Le coefficient de transfert thermique moyen pour la condensation en film sur un tube

extérieur horizontal est déterminé par la corrélation suivante:
0.25
 g l ( l  v ) Lvkl 3 
h  0.729  
 l (Tsat  Tp ) D 
5 5
Lv  Lv   Cpl  (Tsat  Tp )  2407 103   4178  (40  30)  2433112,5J / kg
8 8
0.25
 9,81 994  (994  0, 05)  2433112,5  0, 6233 
h  0.729  3   9292,35 W / m2 C
 0, 720 10 (40  30)  0, 03 
Q  hA(Tsat  Tp )  h DL(Tsat  Tp )  9292,35    0,03 1 (40  30)  8757,83 W

2- La chaleur échangée entre la vapeur et l’eau circulant dans l’ensemble des tubes.
0.25
 g  (   v ) Lvkl 3  1
htubes  0.729  l l    htube
 l (Tsat  Tp ) D  N  N 0,25
Où : N est le nombre de rangée de tubes en vertical

1
htubes  0,25
 9292,35  7060, 66 W / m2 C
3

Q  htubes Ntotal DL(Tsat  Tp )  7060,66 12    0,03 1 (40  30)  79854,18 W

4
 3- Le débit de condensation de la vapeur.

Q 79854,18
mcondensation    0,033 kg / s
Lv 2407 103

Exercice 6

le régime d’ébullition est déterminé à partir de la courbe d’ébullition en fonction de T ,

dans notre cas T  Tp  Tsat  350  100  250C , donc le régime correspondant est celui de
l’ébullition en film.
Pour l’ébullition en film sur un cylindre horizontal la densité du flux thermique est donné
par :

0.25
 gkv 3 v ( l  v )  Lv  0.4  Cpv  (Tp  Tsat )  
q film  0, 62   (Tp  Tsat )
 v D(Tp  Tsat ) 
0.25
 9,81 0, 03573  0, 441(957,9  0, 441)  2257  103  0.4 1977  (250)  
 0, 62    (250)  5,93 104 W / m2
5 3
 1, 7310  510 (250) 

La densité du flux par rayonnement est calculée par cette relation:

qray      (Tp 4  Tsat 4 )  0,05  5,67 108  (250  273)4  (100  273)4   157W / m2

La densité du flux total est calculée par la formule :

3 3
qtotal  q film  qray  5,93 104  157  5,94 104W / m2
4 4

Donc le flux de chaleur échangé :

Qtotal  A  qtotal    D  L  qtotal    5 103 1, 2  5,94 104  1119,66W