PROBLMES CORRIGS

Il nest pas de problme quune absence de

solution ne finisse par rsoudre.
Aphorisme attribu Henri QUEUILLE

PROBLME N 1 : Coefficient dchange

nonc
Le transfert de chaleur entre deux fluides seffectue travers un tube dacier de
diamtres intrieur/extrieur 18 / 21 mm.
On donne :
- ct intrieur : h1 = 1000 W / m 2 K ; temprature moyenne de mlange T1 = 10 C
- ct extrieur : h2 = 2000 W / m 2 K ; temprature T2 = 25 C
- acier : = 46 W / m .K
1. Calculer le coefficient global dchange k.
2. Aprs un an de fonctionnement, on estime avoir une rsistance dencrassement
Re = 4.10 4 W 1 m 2 K . Dterminer le nouveau coefficient dchange global.
3. En attribuant une efficacit de 1 au tube neuf, que devient cette efficacit au bout
dun an ?
4. Quel est alors le flux chang dans un tube de longueur L = 1 m ?

Solution
1.
La paroi du tube a pour paisseur :
D d 21 18
e=
=
= 1, 5 mm = 1, 5.10 3 m
2
2
Comme elle est mince par rapport aux diamtres, en ngligeant sa courbure on peut
calculer k partir de la formule (6.2a) relative une paroi plane, avec Re = 0 :
1, 5.10 3
1
1
e
1
1
1
=
+ +
=
+
+
k h1 h2 1000
46
2000

1
= 10 3 + 0 ,03.10 3 + 0 , 5.10 3 = 1, 53.10 3
k
On constate que la rsistance thermique e / de la paroi ne reprsente ici que 2% de
la rsistance totale. Enfin :

k = 653 W / m 2 K
2.
En prsence dune rsistance dencrassement, on applique maintenant la formule (6.2a)
complte :
1
1
e
1
=
+ Re + +
k h1
h2

= ( 1 + 0 ,4 + 0 ,03 + 0 ,5 ) 10 3 = 1,93.10 3
k = 518 W / m 2 K

rel

3.
Lefficacit dont il est question dans lnonc doit tre comprise comme un rapport
/ max (dfinitions 3.5 et 4.37), soit ici :
E=

1 an
neuf

k1 an
k neuf

518
653

E = 0 ,793

4.
La surface dchange nest pas la mme des deux cts. Suite la question 1, on

calcule une valeur approche de par (6.10) ( 6.2.3) :

D+d
21 + 18
=
L=
10 3 1
2
2
= 61.10 3 m 2
Le flux chang est donn par (6.2c) qui scrit avec les notations de lnonc :
= k ( T2 T1 ) = 518 61.10 3 15

= 475 W
Commentaires
Cet exercice trs lmentaire a surtout pour objet de matrialiser des ordres de
grandeur.
Pour les tudiants, sa principale difficult rside dans le calcul de lpaisseur de la
paroi, qui nest pas D d comme on le lit souvent !!

PROBLME N 2 : Mthode NUT

nonc
Un changeur contre-courant fonctionne dans les conditions suivantes :
Tce = 350 C
T fe = 120 C
Tcs = 200 C

T fs = 290 C

qt min = qtf

Puissance = 415 kW

1. Quelle est la puissance change si on fait travailler lchangeur en mode cocourant, avec les mmes tempratures dentre et les mmes dbits ? (Utiliser la mthode
NUT).
2. Quelles sont les nouvelles tempratures de sortie ?

Solution
1.
En contre-courant, avec qt min = qtf , on a pour efficacit (cf. 3.11) :
E=
R=

T fs T fe
Tce T fe

290 120
= 0 ,74
350 120

Tce Tcs
350 200
=
= 0 , 882
T fs T fe
290 120

et, partir du tableau 3.1 :

1 RE
1
NUT =
Ln
= 2 , 45
1R
1E
En circulation co-courant, les dbits ntant pas modifis, les coefficients dchange ne
le sont pas non plus. On garde donc le mme NUT (vu que NUT = k / qt min ). Par contre,
la nouvelle efficacit E scrit (tableau 3.1) :
1
E' =
{1 exp [1 ( 1 + R ) NUT ]}
1+ R
Il vient, aprs remplacement de R et NUT par leur valeur :
E' = 0 , 526
Puisque les conditions dentre sont identiques dans les deux cas, la nouvelle
puissance ' est telle que :
' E'
=

E
ceci daprs la relation (3.12). Alors :
0 , 526
' = 415
0 ,74

' = 295 kW
cest--dire 70% de la puissance en contre-courant.

2.
La nouvelle temprature de sortie froide T 'fs sobtient partir de la nouvelle efficacit :
E' =

T 'fs T fe
Tce T fe

T 'fs = T fe + E' ( Tce T fe )

= 120 + 0 , 526 ( 350 120 )
T 'fs = 241 C

et la nouvelle temprature de sortie chaude Tcs' partir de R :

R=

Tce Tcs'
T 'fs T fe

= 0 , 882 inchang

Tcs' = Tce R ( T 'fs T fe )

= 350 0 , 882 ( 241 120 )

Tcs' = 243 , 3 C

Commentaires
Il y a dautres faons de rsoudre ce petit exercice, par exemple en utilisant les
rsultats du chapitre 2. Mais la dmarche prconise permet de se roder la mthode NUT
dans un cas simple.
Avec cette valeur de R, on est dj dans la zone asymptotique dun changeur cocourant (fig. 3.1), et les tempratures de sortie sont donc trs voisines.

PROBLME N 3 : changeur bitube

nonc
Pour refroidir un dbit de 9,4 kg/h dair de 616 C 178 C, on le fait passer dans le
tube central dun changeur bitube contre-courant de 1,5 m de long, de 2 cm de diamtre et
de faible paisseur.
1. Calculer la puissance calorifique vacuer. On donne pour lair :
C pc = 1060 J / kg K .
2. Le fluide de refroidissement est de leau, qui pntre dans la section annulaire la
temprature de 16 C avec un dbit de 0,6 l/mn. Calculer la temprature de cette eau la
sortie de lchangeur. On prendra C pf = 4180 J / kg K .
3. Calculer le coefficient dchange hc ct fluide chaud (on ne tiendra pas compte
dune ventuelle correction en / p ).
4. Dterminer lefficacit de cet changeur, puis son NUT. En dduire le coefficient
dchange global, puis le coefficient dchange h f ct fluide froid.
5. La paroi extrieure de lchangeur est isole. Quelle est approximativement
lpaisseur b de lespace annulaire qui permettrait dobtenir cette valeur de h f ? (On
admettra dabord lcoulement laminaire, et on vrifiera ensuite cette proprit).

Solution
1.
Le flux total peut se calculer ct chaud :
= qtc ( Tce Tcs )
Daprs les donnes, le dbit thermique unitaire chaud est :
9 ,4
qtc = q mc C pc =
1060
3600
qtc = 2 ,77 W / K
et alors :
= 2 ,77 ( 616 178 )

= 1213 W
2.
Le calcul du flux total ct froid va maintenant nous donner T fs .

= qtf ( T fs T fe ) avec ici :

qtf = q mf C pf =
qtf = 41, 8 W / K

do :
T fs = T fe +

qtf

0 ,6
4180
60

= 16 +

( q mf = 0 , 3 kg / mn )

1213
16 + 29
41, 8

T fs 45 C

3.
Il faut dabord connatre le rgime dcoulement de lair, donc le Reynolds ct chaud.
La temprature moyenne de lair est approximativement ( 6.2.1) :
T + Tcs 616 + 178
< Tc > = ce
=
= 397 C = 670 K
2
2
A cette temprature, les tables donnent :
c = 0 ,525 kg / m 3 ; c = 6 , 20.10 5 m 2 / s
d 2 ( 2.10 2 )2
=
= 3 ,14.10 4 m 2
4
4
On en dduit la vitesse dbitante :
q
9 ,4
1
Vc = mc =
c S c 3600 0 ,525 3,14.10 4

La section du tube est : S c =

Vc = 15 , 8 m / s
do le nombre de Reynolds :
V d 15 ,8 2.10 2
Rec = c =
5100
c
6 , 20.10 5
Il sagit dun rgime de transition. On peut donc utiliser la formule (4.26), en notant
que le rapport d / L = 2 / 150 est ngligeable (lnonc nous demande galement dignorer la
correction en / p ):

0 ,116
Rec2 / 3 125 Prc 2 / 3
Rec
Pour calculer hc , il est un peu plus rapide ici de passer par le nombre de Nusselt
(4.10d, 4.1.5) :
St c =

Nu c = St c Rec Prc = 0 ,116 Rec2 / 3 125 Prc1 / 3

A 670 K, le nombre de Prandtl de lair est : Prc = 0 ,68 . On trouve :

Nu c = 17 , 25
Toujours 670 K, la conductivit de lair est : c = 0 ,0505 W / m K .
h d
17 , 25 0 ,0505
Nu c = c
do hc =
c
0 ,02
hc = 43 , 5 W / m 2 K

4.
On constate que :
qt min = qtc = 2 ,77 W / K , do lefficacit (formule 3.10) :
E=

Tce Tcs 616 178 438

=
=
Tce T fe
616 16
600

E = 0 ,73

Daprs le tableau 3.1, pour un changeur contre-courant :

1 R E
1
NUT =
Ln
1 R
1 E
et dans le cas prsent :
qt min
2 ,77
R=
=
= 0 ,066
qt max
41, 8
donc :
NUT =

1 0 ,066 0 ,73
1
Ln
= 1,07 Ln 3 , 52
1 0 ,066
1 0 ,73

NUT = 1, 35

De la dfinition du NUT (3.14a) on tire alors :

NUT qt min
k =

Puisque lpaisseur du tube central est faible, on ne fait pas la distinction entre surface
dchange ct chaud et ct froid, et on nglige la rsistance thermique de la paroi. Donc :
= surface latrale du tube = d L

= 2.10 2 1,5 = 0 ,094 m 2

k=

1, 35 2 ,77
0 ,094

k = 39 ,8 W / m 2 K
Le coefficient dchange global sexprime aussi partir de (6.2a) (e et Re tant
ngligs) :
1
1
1
1
1
1
=
+
=
+
=
k hc
hf
43 , 5 h f
39 , 8
On en dduit :
h f 500 W / m 2 K

5.
Dans un change eau-air, la temprature de paroi est proche de celle de leau. Dans le
cas prsent, celle-ci varie peu. On peut donc admettre la condition T p cte . La paroi
extrieure (concave) tant isole, et lcoulement suppos laminaire, la formule (4.44)
sapplique. Elle donne Nu en fonction de R2 / R1 .
Dautre part, Nu = h f Dh / f , avec ici Dh = 2 ( R2 R1 ) = 2 b (formule 4.42). En
partant dune valeur arbitraire mais raisonnable de b, on peut par approximations successives
ajuster Nu avec la valeur donne par (4.44).
La temprature moyenne approche de leau (fluide froid) est :
T fe + T fs 16 + 45
< Tf > =
=
= 30 , 5 C
2
2

La conductivit correspondante (voir tables) est f = 0 ,612 W / m .K .

Essayons avec R2 R1 = b = 3 mm (ce qui fait

R2 10 + 3
=
= 1,3 ).
R1
10

500 2 3.10 3
= 4 , 92
0 ,612
Daprs (4.44), pour obtenir cette valeur de Nu, il faudrait un rapport R2 / R1 = 1,05
environ, soit b = 0 , 5 mm , ce qui est trop faible. Essayons avec b un peu plus lev, pour
augmenter Nu.
Par exemple, avec b = 3 , 2 mm (soit R2 / R1 = 1,32 ) , Nu = 5 , 23 ce qui correspond
peu prs un rapport R2 / R1 = 1,35 dans (4.44). On admettra donc comme valeur
approche :
Nu =

b = 3 , 2 mm

On vrifie enfin le Reynolds :

V f Dh
q mf 2 b
Re f =
=
f
f Sf f
avec

f = 0 ,083.10 5

< T f > = 30 , 5 C

et,

en

premire

approximation

S f = b ( d + 2 ( b / 2 )) surface dun rectangle de hauteur b et de longueur gale la

circonfrence moyenne de lannulaire, soit ( d + b ) (N.B. le diamtre moyen est d + b , et

le diamtre extrieur d + 2 b , voir Problme 1).
0 ,6
Dautre part, q mf = 0 ,6 kg / mn =
kg / s
60
q mf
2b
0 ,6
1
2
1
Re f =
=
3

3
f b ( d + b ) f
60 10 ( 20 + 3 , 2 ) 10
0 ,083.10 5
Re f = 331
Lcoulement est bien laminaire.

Commentaires
Dans cet exercice, on doit en particulier chercher une caractristique gomtrique de
lchangeur permettant de respecter les conditions thermiques imposes.
On a en plus loccasion daborder rgime de transition et coulement annulaire.

PROBLME N 4 : Chemine
nonc
On veut estimer la chute de temprature des fumes dans une chemine, en
considrant le conduit comme un changeur dont les fumes constituent le fluide chaud, et
lair ambiant le fluide froid. On admet que la temprature Ta de lair est constante le long de
la paroi extrieure de la chemine. On dsigne par k le coefficient global dchange travers
la paroi.
1. En adaptant le calcul dun changeur co-courant au cas particulier ci-dessus
( T f = Ta = cte ), montrer que la temprature des fumes dans la chemine obit la loi :

Tcs

Tc Ta
= exp
S
Tce Ta
qtc
2. Le conduit est cylindrique, de diamtre D et de longueur L. crire la temprature
de sortie des fumes.
3. Calculer Tcs avec les valeurs suivantes : L = 20 m ; D = 30 cm ; Tce = 320 C ;

Ta = 10 C ; k = 20 W / m 2 K . Pour les fumes : q mc = 0 , 5 kg / s ; C pc = 1050 J / kg K .

Solution
1.
A travers un lment de paroi dS de la chemine schappe un flux de chaleur :
d = k ( Tc T f ) dS = k ( Tc Ta ) dS
d est aussi la chaleur perdue par le fluide chaud :
d = qtc dTc

On regroupe les deux quations :

dTc
k
=
dS
Tc Ta
qtc
Nous avons ici Ta = cte , de sorte que :
dTc d ( Tc Ta )
do en intgrant :
k

Tc Ta = cte exp
S
qtc
A lentre, S = 0 et Tc = Tce , donc cte = Tce Ta , et

Tc Ta
= exp
S
Tce Ta
qtc

C.Q.F.D.

2.
A la sortie, Tc = Tcs et S = (surface latrale totale de la chemine) :

 = DL
La formule de la question 1 devient :
k

Tcs = Ta + ( Tce Ta ) exp

D L
qtc

3.
Les calculs numriques donnent :
qtc = q mc C pc = 0 , 5 1050 = 525 W / K
20

Tcs = 10 + ( 320 10 ) exp

0 ,3 20
525

Tcs = 161 C

La chute de temprature est donc de 160 C environ dans la chemine.

Commentaires
1. La formule donnant Tcs (question 2) est recommande par le Centre Scientifique et
Technique du Btiment (CSTB) pour le calcul des chemines.
2. Dans la question 1, on ne peut pas crire le bilan local sur le fluide froid
( d = qtf dT f ) du fait que dT f = dTa = 0 ( Ta = cte ), ce qui est quivalent qtf . On
se retrouve dans la mme situation quavec un vaporateur ou un condenseur ( 2.4 et 3.5).
3. Lhypothse Ta = cte est acceptable ( Ta nest pas ici une temprature de mlange
mais la temprature du fluide extrieur loin de la paroi) ; elle conduit une formule simple
qui donne des estimations numriques correctes.
4. Dans une revue technique, la formule de la question 2 tait accompagne du
commentaire suivant : est la section de la chemine !! Ce qui donnerait ici une chute de
temprature infrieure 1 C. De limportance de la rigueur dans le vocabulaire !!

PROBLME N 5 : changeur changement de phase

nonc
De leau froide circule dans un tube de chaudire condensation. Sa temprature
dentre est T fe = 18 C et son dbit q mf = 400 kg / h . Le rchauffage est assur par
condensation de vapeur deau lextrieur du tube, la temprature Tc = 104 C .
On donne : diamtre intrieur d = 12 , 5 mm ; diamtre extrieur D = 16 mm ;
longueur L = 2 , 4 m ; conductivit de la paroi p = 46 W / m K . Pour leau, on admettra
dans la gamme de temprature considre : f = 0 , 7.10 6 m 2 / s ; Pr f = 5 , 5 ;
C pf = 4180 J / kg K .
1. Calculer le coefficient dchange h f lintrieur du tube.
2. On donne le coefficient h ct vapeur : hc = 8000 W / m 2 K . Calculer le coefficient
global dchange k.
3. Calculer le NUT et lefficacit de lappareil.
4. Dterminer la temprature de sortie deau T fs , puis la quantit de chaleur Q
rcupre annuellement grce au dispositif condensation, si lon considre que la saison de
chauffe dure 150 jours et que la chaudire fonctionne 5 heures par jour.

Solution
1.
Il faut dj savoir quel est le rgime dcoulement. Pour cela, calculons Re :
4 qv
V d qv d
Re =
=
=

S d
Ici, q v = 400 l / h =

400.10 3
m3 / s
3600
1

0 ,4

3600 12 , 5 10 3 0 ,7 10 6
Re = 16170
Re = 4

Avec cette valeur de Re, le rgime est turbulent et hydrauliquement lisse. Donc, St est
donn par (4.21). On a bien Pr > 0 ,66 . Dautre part, L / d = 2400 / 12 , 5 > 60 . Les
conditions de validit de la formule sont donc satisfaites.
Pour le fluide froid :
St f = 0 ,023 Re 0 , 2 Pr 0 ,6 = 0 ,023 16170 0 ,2 5 , 5 0 ,6
St f = 1,19.10 3

et aussi

St f =

hf

CpV

do :

h f = St f C p V = St f C p

4 qv

d2

On a = 10 3 kg / m 3 pour leau et d = 12 , 5 mm , ce qui donne :

h f = 4500 W / m 2 K

2.
Lpaisseur de la paroi est :
D d 16 12 , 5
e=
=
= 1,75 mm
2
2
La question ne prcise pas de quel ct on demande k. Cela sous-entend que lon
raisonne comme avec une paroi plane ( 6.2.1. et 6.2.3), et k scrit :
1,75.10 3
1
1
e
1
1
1
=
+ +
=
+
+
= 3 , 85.10 4
hf
k hc
4500
46
8000

k 3000 W / m 2 K
3.
Lchangeur est un appareil fluide isotherme, pour lequel R = 0 ( 3.5.2), avec ici
qt min = qtf . Ainsi :
NUT = k / qtf
o est la surface moyenne dchange ( 6.2.2) :
D+d
16 + 12 , 5
=
L=
10 3 2 , 4 = 0 ,107 m 2
2
2
Le dbit thermique unitaire de leau a pour valeur :
400
4180 = 465 W / K
qtf = q mf C pf =
3600
et le NUT vaut :
NUT = 0 ,69

Lefficacit se trouve sur le tableau 3.1. Pour R = 0 , quelque soit le modle

dchangeur :
E = 1 exp ( NUT ) = 1 exp ( 0 ,69 )
E 0 ,5

4.
Toujours pour la raison que qt min = qtf , E est lefficacit ct froid (formule 3.11),
avec Tce = Tc = cte :
E=

T fs T fe
Tc T fe

On en tire T fs :
T fs = T fe + E ( Tc T fe ) = 18 + 0 , 5 ( 104 18 )
T fs = 61 C

La puissance thermique du systme sexprime par le bilan sur le fluide froid :

= qtf ( T fs T fe ) = 465 ( 61 18 ) = 20000 W
= 20 kW
Sur une saison de chauffe qui compte 150 5 heures, on peut donc rcuprer par
condensation au maximum :
Q = 20 150 5 kWh
Q 15000 kWh

Commentaires
Il sagit encore dun problme allg o on ne fait mention ni des tempratures de
rfrence, ni de la distinction entre coefficient dchange global ct froid et ct chaud, et o
le calcul de h ct vapeur nest pas demand.
La temprature de condensation de 104 C sous-entend que la vapeur est en lgre
surpression par rapport la pression atmosphrique.
Les questions 3 et 4 peuvent aussi tre traites avec les rsultats relatifs aux
changeurs fluide isotherme ( 2.4).

PROBLME N 6 : changeur faisceau de tubes et calandre

nonc
Dans la sous-station de chauffage collectif dun immeuble on dsire installer un
changeur faisceau tubulaire et calandre, destin porter de 40 60 C un dbit deau de
20000 kg/h. Le fluide primaire qui circule dans les tubes est de leau surchauffe arrivant
180 C, raison de 10000 kg/h. Les tubes ont un diamtre intrieur d = 20 mm ; la vitesse
dcoulement adopte est telle que Re = 10000. Le coefficient dchange global k est estim
450 W / m 2 K . On admet pour leau surchauffe les caractristiques thermophysiques
suivantes : C pc = 4315 J / kg K ; = 920 kg / m 3 ; = 19.10 5 kg / m . s .
1. Calculer la puissance change et la temprature de sortie du fluide chaud.
2. Lchangeur est contre-courant, avec une seule passe sur chaque fluide, les tubes
tant monts en parallle (modle de la fig. 1.5 sans les chicanes). Dterminer :
- la surface dchange ncessaire
- la vitesse dans les tubes
- la section totale des tubes
- le nombre de tubes et la longueur du faisceau

Solution
1.
Calculons tout dabord les dbits thermiques unitaires.
- Sur le fluide froid :
qtf = q mf C pf
On prend gnralement pour leau C pc = 4180 J / kg K (voir tables en fin de volume).
20000
4180
3600
qtf = 23200 W / K
- Sur le fluide chaud :
10000
qtc = q mc C pc =
4315
3600
qtc = 12000 W / K
qtf =

La puissance thermique sobtient partir des donnes relatives au fluide froid :

= qtf ( T fs T fe ) = 23200 ( 60 40 ) = 46400 W

= 464 kW
On tire maintenant Tcs du bilan sur le fluide chaud :
= qtc ( Tce Tcs )
464000 = 12000 ( 180 Tcs )
Tcs = 180 38 ,7
Tcs = 141, 3 C

2.

Les calculs prcdents montrent que nous avons ici :

qt min = qtc

et lefficacit scrit :
T Tcs
38 ,7
E = ce
=
Tce T fe
140
E = 0 , 276
Quant au facteur de dsquilibre :
q
12000
R = tc =
qtf
23200
R = 0 , 517

Du tableau 3.1 on tire le NUT :

1 RE
1 0 , 517 0 , 276
1
1
NUT =
Ln
=
Ln
1R
1 E
1 0 , 517
1 0 , 276
NUT = 0 , 35
et aussi :
NUT = k / qtc
ce qui nous donne la surface dchange :
0 , 35 12000
=
450

= 9 , 33 m 2

Le nombre de Reynolds est impos :

Vd
Re = 10000 =

do la vitesse V :
19.10 5 10 4
V =
920 20.10 3
V = 0 ,103 m / s

La vitesse V est la mme dans tous les tubes. Si q v dsigne le dbit-volume total,
la section totale des tubes est donc :
q
q
1 10000 1
1
S = v = mc
=
V
V
3600 920 0 ,103

S = 2 ,9.10 2 m 2

Le nombre de tubes ncessaires pour assurer le dbit demand est gal :

N =

d2 / 4

2 ,9.10 2 4

( 2.10 2 ) 2

= 92 , 3

Comme il faut bien avoir un nombre entier de tubes, on prendra lentier

immdiatement suprieur :
N = 93 tubes

Si la longueur du faisceau est L (cest donc aussi la longueur de chaque tube), la

surface totale dchange a pour valeur :
= N dL
On en dduit :
9 , 33

L=
=
N d 93 0 ,02
L = 1, 60 m
soit une longueur totale de tubes de 149 m environ.
Commentaires
1. Lexercice aborde dune manire encore lmentaire le dimensionnement dun
changeur faisceau de tubes.
2. Les tudiants qui nont jamais vu un changeur tubulaire nimaginent pas toujours
trs bien la structure de lensemble et les proprits qui dcoulent du montage des tubes en
parallle (FEMM, 7.6).
3. Le NUT de lchangeur est faible. Le modle de lappareil na donc pas beaucoup
dimportance ( 3.4.3 et fig. 3.1).
4. Un calcul approch de pourrait tre effectu partir du bilan global en
linarisant la distribution de temprature de leau ( 6.2.1) :
T fe + T fs
= k ( Tc < T f > ) avec < T f > =
2
Cest plus rapide et on trouve ici pratiquement le mme rsultat (parce que NUT est
petit). Mais ce nest pas toujours comme a. Alors prudence !!
5. Les paramtres k et se rfrent implicitement la surface intrieure des tubes
puisquon ne donne pas le diamtre extrieur dans lnonc.

PROBLME N 7 : changeur plaques

nonc
De lair chaud 80C produit par une installation industrielle est vacu lextrieur
raison de 5000 kg/h. Il traverse pralablement un changeur servant prchauffer 20C
lair neuf admis dans le btiment une temprature dentre de 5C.
Lchangeur est un appareil plaques planes et courants croiss, sans brassage des
fluides, dont la surface dchange est = 50 m 2 . Dans la plage de fonctionnement prvue,
son coefficient global dchange k est estim par le constructeur :
k = 3 ,6 q mf + 65
avec q mf en kg/s et k en W / m 2 K .
1. Estimer a priori sur quel circuit le dbit thermique unitaire est minimal.
2 En procdant par approximations successives, dterminer quel dbit dair neuf cet
changeur peut porter 20C.
3. Calculer la temprature de sortie du fluide chaud et la puissance de lchangeur.
P.J. Abaque E = f(NUT), fig. P.7.

FIG. P.7. changeur courants croiss, fluides non brasss

Solution
1.
On a T = ( 20 5 ) = 15 C sur le fluide froid. Le fluide chaud entre Tce = 80 C ,
ce qui doit permettre un T nettement plus important. Les fluides tant les mmes des deux
cts, q mf est certainement suprieur q mc , do qtf > qtc car la chaleur massique C p
dpend trs peu de la temprature. Par consquent, a priori :

qt min = qtc
puisque le fluide qui accuse le plus grand T est celui qui a le plus petit dbit thermique
unitaire ( 3.1).
2.
Les donnes ne portent que sur trois tempratures aux bornes et un dbit. Le dbit
inconnu q mf ne peut pas tre calcul directement ; lefficacit non plus puisque ici E = E c et
on ne connat pas Tcs . On va donc procder par approximations successives partir dune
estimation empirique de q mf , en calculant NUT, puis E au moyen de labaque ( !! dernier
rappel : abaque est un nom masculin, on dit un abaque, fermez la parenthse) et enfin T fs que
lon va comparer la valeur impose 20 C.

Choisissons comme valeur de dpart :

q mf = 15000 kg / h = 4 ,166 kg / s
NUT =

k
k
=
qt min
qtc

avec

= 50 m 2

q mc = 5000 / 3600 = 1, 39 kg / s

Daprs les tables numriques, C pc 1006 J / kg K trs peu dpendant de T. Soit :

qtc = q mc C pc = 1400 W / K

Dans ces conditions :

k = 3 ,6 q mf + 65 = 3 ,6 4 ,166 + 65 = 80 W / m 2 K
80 50
= 2 , 86
1400
q
q
puisque
R = tc = mc
qtf
q mf
NUT =

C pc C pf

R = 1, 39 / 4 ,166 = 0 , 333

En reportant sur labaque, on trouve pour lefficacit :

E 0 , 85

Puisque qt min = qtc , nous avons (formules 3.10) :

E = Ec =

1
1 T fs T fe
Ef =
R
R Tce T fe

do lon tire :
T fs = T fe + R E ( Tce T fe )

= 5 + 0 ,333 0 ,85 ( 80 5 )
soit :

q mf

T fs = 26 , 25 C
Cette valeur est suprieure la temprature impose (20 C) ; on a donc pris un dbit
trop faible.

Recommenons lopration avec :

q mf = 25000 kg / h = 6 , 945 kg / s
La trame du raisonnement est identique ; qtc est inchange = 1400 W / K
k = 3 ,6 6 , 945 + 65 = 90 W / m 2 K ;

NUT =

90 50
= 3, 2
1400

R = 1, 39 / 6 ,945 = 0 , 2
et labaque donne : E 0 , 9
T fs = 5 + 0 , 9 0 , 2 ( 80 5 )
T fs = 18 , 5 C

Nous nous sommes nettement rapprochs de 20 C, mais maintenant T fs est un peu

trop faible ; donc q mf est un peu trop grand.

Faisons une troisime tentative avec q mf = 22500 kg / h = 6 , 25 kg / s

k = 3 ,6 6 , 25 + 65 = 87 ,5 W / m 2 K
87 , 5 50
NUT =
= 3 ,125
1400
R = 1, 39 / 6 , 25 = 0 , 222
Sur labaque, on observe que E est pratiquement inchange 0 , 9
T fs = 5 + 0 , 9 0 , 222 ( 80 5 )
T fs 20 C :

La convergence est atteinte. La valeur demande est donc :

q mf = 22500 kg / h = 6 , 25 kg / s

3.
Nous pouvons par exemple dduire Tcs du bilan global de lchangeur (3.1) :
qtc ( Tce Tcs ) = qtf ( T fs T fe )
do :
Tcs = Tce

qtf

( T fs T fe ) = 80

qtc
puisque qtf / qtc = q mf / q mc . Finalement :

6 , 25
( 20 5 )
1, 39

Tcs = 12 , 5 C

Le bilan sur le fluide chaud (ou froid) donne enfin la puissance de lchangeur :
= qtc ( Tce Tcs ) = 1400 ( 80 12 ,5 )

= 94500 W

Commentaires
Cet exercice illustre lun des rares cas o lon doive procder par itration avec la
mthode NUT .
Pour le calcul de T fs (question 2), on peut aussi calculer pralablement Tcs partir de
E puis dduire T fs du bilan enthalpique de lchangeur.

PROBLME N 8 : changeur tubulaire courants croiss

nonc
Une chaudire mazout produit des fumes 360 C qui sont utilises pour
prchauffer de 20 120 C lair ncessaire la combustion.
Linstallation consomme 1500 kg de mazout par heure. La combustion dun kilog de
mazout demande 17,95 kg dair et produit 18,95 kg de fumes.
Lappareil utilis est un changeur tubulaire courants croiss, deux passes ct
tubes ; lair circule dans les tubes, verticalement, et les fumes autour des tubes,
horizontalement (cette disposition limite les obligations de ramonage).
La temprature de paroi T p des tubes doit tre partout suprieure 160 C pour
viter la corrosion due la condensation des composs soufrs.
Les tubes forment un faisceau en ligne, pas carr, de pas relatif e + = 1, 4 . Ils ont
pour diamtre extrieur D = 55 mm et pour diamtre intrieur d = 50 mm .
Pour chaque coulement, le nombre de Reynolds est fix : Re f = 25000 dans les tubes
(fluide froid) et Rec = 5000 en calandre (fluide chaud, vitesse de rfrence = vitesse
dbitante en calandre vide).
Dans le domaine de temprature envisag, on prendra comme caractristiques
moyennes des fumes :

= 3.10 5 kg / m . s

C p = 1212 J / kg K

= 0 ,033 W / m K

= 0 ,585 kg / m 3

1. Calculer la temprature de sortie des fumes Tcs .

2. Lchangeur peut fonctionner soit avec des entres opposes, soit avec des entres
du mme ct (voir figure P8). Quel sens de fonctionnement doit tre choisi pour viter le
risque de condensation des composs soufrs ?
3. Calculer quelles vitesses de rfrence Vc0 et V f0 correspondent les nombres de
Reynolds imposs.
4. Dterminer le coefficient global dchange k c ct fluide chaud (ngliger la
rsistance thermique de la paroi).
5. Calculer la puissance de lchangeur, son efficacit E et son NUT. Dterminer
la surface dchange c ct chaud, le nombre de tubes ncessaires et la longueur L de
chaque tube.

FIG. P8 changeur tubulaire courants croiss.

Fluide chaud : fumes ; fluide froid : air
A gauche : entres opposes. A droite : entres du mme ct

Solution
1.
La temprature inconnue Tcs sera obtenue partir du bilan enthalpique de
lchangeur :
qtc ( Tce Tcs ) = qtf ( T fs T fe )
Il faut dabord calculer les dbits thermiques unitaires.
- Pour le fluide froid (air), le dbit-masse est :
1500 17 ,95
= 7 ,5 kg / s
qmf =
3600
La chaleur massique de lair dpend trs peu de la temprature. On larrondit
C pf 1010 J / kg K , soit :
qtf = q mf C pf = 7 , 5 1010 = 7 , 57.10 3 W / K
- Pour le fluide chaud (fumes) :
1500 18 , 95
q mc =
= 7 , 9 kg / s
3600
qtc = 7 , 9 1212 = 9 , 57.10 3 W / K
On observe au passage que :

qt min = qtf = 7 , 57.10 3 W / K

Du bilan enthalpique on tire :
qtf
7 , 57
Tce Tcs =
( T fs T fe ) =
( 120 20 ) = 79 C
qtc
9 , 57
Tcs = 281 C

2.
Air et fumes ont des proprits voisines, et les coefficients dchange lextrieur et
lintrieur des tubes doivent tre du mme ordre de grandeur. Comme temprature locale de
paroi, on admettra donc en premire approximation :
T p ( Tair + T fumes ) / 2
Avec entres opposes, on a :
* ct (a) : air 20 C ; fumes 281 C, soit :
20 + 281
T pa
= 150 , 5 C < 160 C
do risque de condensation
2
* ct (b) : air 120 C ; fumes 360 C, soit :
360 + 120
T pb
= 240 C
pas de problme
2
Cette solution est proscrire : il y a un fort risque de condensation vers la sortie des
fumes.
Avec entres du mme ct :
* ct (a) : air 120 C ; fumes 281 C. Donc :
120 + 281
200 C > 160 C
cest bon
T pa
2
* ct (b) : air 20 C ; fumes 360 C :
20 + 360
T pb
= 190C > 160C
cest encore bon
2
Conclusion : il vaut mieux placer les entres du mme ct

3.
Vitesse de rfrence V f0 ct fluide froid :
Il sagit ici de lcoulement de lair dans les tubes.
Re f =

V f0 d

= 25000

donn, avec d = 50 mm

f doit tre value la temprature moyenne de mlange de lair ( 4.1.5 ), soit :

Tmf =

20 + 120
= 70 C
2

f = 2 ,02.10 5 m 2 / s
On en dduit :
25000 2 ,02.10 5
0
Vf =
50.10 3
V f0 = 10 ,1 m / s
Vitesse de rfrence Vc0 ct fluide chaud
Il sagit maintenant des fumes, sachant que :

Rec =

Vc0 D

= 5000 avec c =

c 3.10 5
=
= 5 ,13.10 5 m 2 / s
c
0 , 585

on en tire :
Vc0

Rec c
5000 5 ,13.10 5
=
=
D
55.10 3

Vc0 = 4 ,66 m / s

4.
Coefficient dchange hc ct chaud (en calandre)
Le nombre de Reynolds est donn : Rec = 5000 , avec une vitesse de rfrence en
calandre vide, donc compatible avec les formules du parag. 4.3.1.
Le faisceau est en ligne ; de ce fait il rsulte que : eT+ = e L+ = e + = 1,4 . En admettant
que le nombre de nappes est au moins gal 10 ( 4.3.1.3 ), le nombre de Stanton est donc
prendre dans la formule (4.28) avec = 0 ,7 pour le fluide chaud :
0 ,6

1, 4 + 0 , 90
0 ,2

St c = 0 ,023 1 + 6 , 2

1
,
4
5000 0 ,32 Pr 0 ,7

1
,
4

0
,
98

Le nombre de Prandtl est calcul la temprature Tmc , ce qui donne daprs les
valeurs de lnonc :
C p 3.10 5 1212
Prc = =
=
= 1,10
a

0 ,033
et lon obtient :
St c = 0 ,024

Le coefficient dchange hc se dduit de St :

hc = St c c C pc Vc0 = 0 ,024 0 , 585 1212 4 ,66

(la vitesse de rfrence V 0 tant la mme que pour Re, cf. 4.3.1.2).
hc = 79 , 3 W / m 2 K
Coefficient dchange h f ct fluide froid (dans les tubes)

On sait que Re f = 25000 ; en labsence dinformation sur la rugosit de la paroi, il est

raisonnable dadmettre que lon est en rgime turbulent lisse . On utilisera donc la formule
de Dittus et Boelter (4.21) avec = 0 ,6 (fluide froid), sous rserve que lon ait L / d > 60 ,
ce qui devra tre vrifi en fin de calcul :
St f = 0 ,023 Re f 0 ,2 Pr f 0 ,6

Pour lair Tmf = 70 C , Pr f = 0 ,7 .

St f = 0 ,023 25000 0 ,2 0 ,7 0 ,6
St f = 0 ,00376

De l nous dduisons h f :
h f = St f f C pf V f0

A la temprature Tmf = 70 C , les tables donnent pour lair : f = 1,02 kg / m 3 et

C pf 1015 J / kg K . Ainsi :
h f = 0 ,00376 1,02 1015 10 ,1
h f = 39 , 3 W / m 2 K
Coefficient global ct fluide chaud
Le coefficient global k c ct fluide chaud est ensuite calcul avec (6.6) :
1 1
1
kc =
+
c c hc f h f

1
1
+ c
hc
f hf

Le rapport des surfaces dchange est gal au rapport des diamtres D / d = 55 / 50 .

1
kc =
1
55 1
+
79 , 3 50 39 , 3
k c = 24 ,6 W / m 2 K

5.
Puissance
On peut calculer la puissance par exemple ct froid :

= qtf ( T fs T fe ) = 7 ,57.10 3 ( 120 20 )

= 757 kW
Efficacit
Daprs (3.12) :
E=

qt min ( Tce T fe )

757.10 3
7 , 57.10 3 ( 360 20 )

E = 0 , 295
NUT

Lchangeur est courants croiss, un fluide brass ( 1.2), avec une seule passe sur
le fluide brass. Le dbit thermique maximum est celui des fumes, cest--dire du fluide
brass. Le tableau 3.1 donne dans ce cas :
1

NUT = Ln 1 +
Ln ( 1 ER )
R

R=

qt min
qt max

7 , 57
9 , 57

soit

R = 0 ,79

et finalement :
NUT 0 , 41
Surface dchange
De NUT on dduit la valeur de la surface dchange c ct fluide chaud :
NUT =

kc c
24 ,6 c
=
= 0 , 41
qt min 7 , 57.10 3

c 126 m 2
Faisceau de tubes
La surface dchange c ct fluide chaud correspond au diamtre extrieur des

tubes : D = 55.10 3 m . Si Lt est la longueur totale des tubes, on a donc :

126
Lt =
c = D Lt soit
55.10 3
Lt 729 m
Appelons n le nombre de tubes ncessaires dans chaque passe pour assurer le dbit
demand, et S t la section totale des tubes :
St = n

d2

4
Pour avoir n, il faut calculer pralablement S t partir de lexpression du dbit dair :
q mf = f S t V f

Dans un calcul de dbit, la masse volumique f doit tre value la temprature

moyenne de mlange, ici Tmf ( 4.1.4 et 4.1.5 ), ce qui donne f = 1,02 kg / m 3 . En outre,
V f = V f0 = 10 ,1 m / s , et q mf = 7 , 5 kg / s . Alors :
7 ,5
= 0 ,728 m 2
1,02 10 ,1
4 St
n=
avec d = 50.10 3 m
2
d
n = 370 tubes par passe
St =

soit au total :
N = 740 tubes

La longueur de chaque tube (cest--dire la longueur du faisceau) est donne par :

L
729
L= t =
m
n
740
L1m

Vrification
Il y a lieu maintenant de sassurer que la longueur de chaque tube rpond bien la
condition L / d > 60 (question 4, calcul de h f ). Or ce nest pas le cas, puisque nous avons ici
L / d 20 : daprs (4.24), St f doit donc tre multipli par un facteur 1,12. Cela remontera
de quelques pour cent le coefficient dchange global, et il faudrait normalement procder
une itration avec cette nouvelle valeur. Mais on est excusable si on renonce recommencer
les calculs.
En ce qui concerne h f , vu le nombre total de tubes, la condition nombre de nappes
suprieur 10 est raliste.

Commentaires
- La question 2 illustre le genre de raisonnement sur des donnes incertaines que
lingnieur doit frquemment pratiquer (ou devrait tre apte pratiquer .). Bien quelles
soient par essence nuances, les conclusions auxquelles on aboutit sont appeles tayer des
choix.
- Dans le cas prsent, la suite du calcul montre que hc > h f , ce qui a pour effet de
relever la temprature de paroi estime, et donc de faire reculer le risque de condensation
ct (a) dans la formule entres opposes . Malgr tout, cette solution apparat moins sre
que la disposition entres du mme ct .
- Les caractristiques thermophysiques des fumes sont fournies dans lnonc pour
raccourcir les calculs. On doit normalement les estimer partir de la composition des dites
fumes, avec les tables relatives O2 , N 2 et CO2 .

PROBLME N 9 : changeur tubulaire, faisceau en quinconce

nonc
On tudie un changeur air-eau faisceau de tubes et calandre. Leau circule dans
les tubes et lair autour des tubes, perpendiculairement eux. Les tubes sont disposs
verticalement, chaque nappe transversale constituant une passe pour la circulation de leau.
La calandre est rectangulaire, de hauteur H c = 2 ,615 m et de largeur l c dterminer.
Le faisceau est en quinconce, de pas triangulaire quilatral. Le pas transversal entre
tubes est fix ST = 95 mm . Dans chaque nappe transversale, les tubes sont monts en
parallle. Ils sont en acier. Leurs diamtres intrieur/extrieur sont d = 31 mm et
D = 35 mm . Une distance de 0 ,5 ST est impose entre laxe des tubes dextrmit et la
paroi de la calandre. On se fixe une vitesse dbitante maximale pour leau (fluide froid)
Vf = 1 m / s .
Les donnes du problme sont :
Fluide
Eau
Air
Temprature dentre
T fe = 20 C
Tce = 250 C
Temprature de sortie

T fs = 80 C

Dbit

q mf = 12 kg / s

q mc = 50 kg / s

1. Calculer la puissance thermique change et la temprature de sortie de lair.

2. a) Dterminer le nombre N L de tubes ncessaire dans chaque nappe transversale.
(prciser quel est dans ce cas le nombre de nappes longitudinales du faisceau de tubes).
b) En dduire la largeur l c de la calandre et la vitesse Vc0 de lair (fluide chaud)
dans la calandre vide.
3. Calculer le coefficient global dchange k f ct fluide froid. Commenter le
rsultat.
4. a) Dterminer lefficacit de lchangeur, son NUT et la surface dchange f
ncessaire (ct fluide froid).
b) Calculer la longueur totale L des tubes, et le nombre N T de nappes
transversales.
5. La rugosit intrieure des tubes est = 10 2 mm . Pour chaque demi-tour du
fluide (passage dune nappe la suivante), le coefficient de perte de charge singulire est
valu = 0 , 4 . On comptera e = 0 , 5 pour lentre et s = 1 pour la sortie du circuit.
Calculer la perte de charge totale et la puissance utile Pu de la pompe ncessaire pour
alimenter lchangeur en eau.

Solution
1.
La puissance de lappareil se dduit du bilan thermique sur leau :
= qtf ( T fs T fe )
Avec C pf = 4180 J / kg K environ pour leau, le dbit thermique unitaire vaut :
qtf = q mf C pf = 12 4180

qtf = 50160 W / K
do la puissance :
= 50160 ( 80 20 ) W
= 3010 kW

Le bilan thermique sur lair scrit :

= qtc ( Tce Tcs ) = q mc C pc ( Tce Tcs )
En toute rigueur, pour connatre C pc , il nous faut la temprature moyenne Tmc de lair,
qui est inconnue. Mais les chaleurs massiques varient assez peu avec la temprature. En se
basant sur une estimation approximative de 210 C , on a C pc 1025 J / kg K . Do :
qtc = 50 1025
qtc = 51100 W / K
On en dduit la temprature de sortie de lair :

3010.10 3
Tcs = Tce
= 250
qtc
51100
Tcs = 191 C

La temprature moyenne Tmc est donc de 220 C , et lincidence sur la valeur de C pc

est ngligeable. On en reste l.
2.
2a). Nombre de tubes
Le dbit deau est donn :
q mf = f S f V f = 12 kg / s
Les tubes tant en parallle, leur section totale dans une nappe transversale est :
q mf
d2
S f = NL
=
4
f Vf
Avec V f = 1 m / s et f = 10 3 kg / m 3 , il faut donc dans chaque nappe :

NL =

4 12

( 31.10 3 ) 2 10 3 1

= 15 ,9 tubes

La vitesse de 1 m/s tant un maximum, on doit arrondir lentier suprieur :

N L = 16 tubes
Le rajustement qui en rsulte sur la vitesse est ngligeable, et on peut conserver dans
la suite des calculs V f = 1 m / s .
Le nombre de tubes dune nappe transversale est gal au nombre de nappes
longitudinales, soit ici 16 tubes ( 4.3.1.1).
2b). Largeur de la calandre et vitesse de lair
Une range de 16 tubes au pas transversal de 95 mm ncessite un espace de 15 entreaxes = 15 ST , plus 0 ,5 ST gauche et droite entre les tubes et la paroi de la calandre
(condition pose par lnonc).

Mais le faisceau tant en quinconce, la seconde range est dcale de 0 ,5 ST vers la

droite ou vers la gauche, ce qui fait en tout :
l c = ( 15 + 2 0 , 5 + 0 , 5 ) S T = 16 , 5 95.10 3
l c = 1, 57 m

La calandre vide a une section S c = H c l c = 2 ,615 1, 57

S c = 4 ,10 m 2

A la temprature moyenne Tmc = 220 C (question 1), les tables donnent pour la
masse volumique c = 0 ,715 kg / m 3 , do la vitesse dbitante dans la calandre vide :
q
50
Vc0 = mc =
c S c 0 ,715 4 ,10
Vc0 = 17 m / s

3.
Pour calculer le coefficient dchange global, il faut dabord connatre les coefficients
dchange de chaque ct des tubes.
3a). Coefficient h f lintrieur des tubes
Le nombre de Reynolds dans les tubes est : Re f =

Vf d

A la temprature moyenne de mlange de leau :

20 + 80
Tmf =
= 50 C
2
la viscosit cinmatique vaut (cf. tables) : f = 0 ,057.10 5 m 2 / s , do :

Re f =

1 31.10 3

0 ,057.10 5
Re f = 54400
Avec ce nombre de Reynolds ( < 10 5 ), on doit pouvoir utiliser la formule de DittusBoelter (4.21), condition davoir un allongement des tubes L / d > 60 . Ici la longueur des
tubes est la hauteur de la calandre, cest--dire H c . On a effectivement :
Hc
2 615
=
= 84 ,3
d
31
Le compte est bon, et alors :
St f = 0 ,023 Re f 0 ,2 Pr f 0 ,6
A la temprature Tmf = 50 C , les tables donnent Pr f = 3 ,67 .
St f = 0 ,023 54400 0 ,2 3 ,67 0 ,6 = 0 ,023 0 ,113 0 , 456
St f = 1,185.10 3
Daprs la dfinition du nombre de Stanton, avec les valeurs numriques dj connues:
h f = St f f C pf V f = 1,185.10 3 10 3 4180 1

h f = 4950 W / m 2 K
3b). Coefficient h f autour des tubes
Le faisceau est en quinconce, pas triangulaire quilatral. On a donc :

- pour le paramtre gomtrique adimensionnel ( 4.3.1.3 ) : = 1 / 2

S
95
- et pour le pas transversal relatif : eT+ = T =
D
35
+
eT = 2 ,7
Le nombre de Stanton est donn par (4.29b) (dans le cas < 1 et avec
= 0 ,7 pour le fluide chaud) :
0 ,66

29 eT+

0 ,346
St c = 0 ,023 1 + +
1
Prc 0 ,7
Rec
e 1

A la temprature moyenne de lair Tmc = 220 C (question 1), Prc = 0 ,68 et

c = 3,74.10 5 m 2 / s . La vitesse de rfrence pour le nombre de Reynolds est ici la vitesse

dans la calandre vide Vc0 = 17 m / s (question 2), et la longueur de rfrence est le diamtre
extrieur D des tubes do :
Rec =

Vc0 D

17 35.10 3
3 ,74.10 5

Rec = 15900
0 ,66

29 2 ,7

0 ,346
St c = 0 ,023 1 +
1
0 ,68 0 ,7
15900
2 ,7 1

St c = 0 ,0132
On en dduit le coefficient dchange moyen :

hc = St c c C pc Vc0 = 0 ,0132 0 ,715 1025 17

hc = 164 , 5 W / m 2 K
3c). Coefficient global k f
Le coefficient global dchange ct froid est donn par (6.6d) :
f 1
1
e f
1
=
+
+
kf
c hc m h f

avec ici ( 6.2.3) f / c = d / D et f / m = 2 d /( d + D ) , et aussi e = ( D d ) / 2

= 2 mm . Pour la conductivit thermique de lacier = 46 W / m . K .
1
31 1
2.10 3 2 31
1
=
+
+
kf
35 164 , 5
46
31 + 35 4950
1
= 5 , 38.10 3 + 0 ,04.10 3 + 0 , 21.10 3
kf

k f = 177 , 5 W / m 2 K

Dans ce total, le terme dominant est le coefficient dchange ct air. Le coefficient

dchange ct eau compte pour moins de 4%. La rsistance thermique de la paroi est
ngligeable.
4.
4a). Efficacit, NUT, surface dchange
Daprs la 1re question, qt min = qtf . Lefficacit de lchangeur est donc exprime
par (3.11) :
E=

T fs T fe

Tce T fe

80 20
250 20

E = 0 , 26

Lchangeur est de type 1-N. Son NUT est donn dans le tableau 3.1 :
NUT =

2 E (1 + R 1 + R2 )

Ln
1 + R2
2 E (1 + R + 1 + R2 )
Le facteur de dsquilibre vaut :
qt min
50160
R=
=
= 0 , 98
qt max 51100

Tous calculs faits :

NUT = 0 , 36

Nous savons que NUT =

K
qt min

kf f
qt min

La connaissance du NUT permet donc de dterminer la surface dchange f :

f = NUT

qt min
kf

= 0 , 36

50160
177 , 5

f = 101,7 m 2
4b). Nombre de nappes transversales
Si L est la longueur totale des tubes, la surface dchange f (ct fluide froid, cest-dire lintrieur des tubes) scrit :
f =dL
do :
101,7
L=
31.10 3
L = 1045 m

Dans une nappe transversale, il y a N L tubes, dont la longueur totale est N L H c . Le

nombre de nappes transversales N T est donc tel que :
L
1045
NT =
=
= 24 , 97
N L H c 16 2 ,615
La valeur trouve est trs proche de lentier 25. On prendra donc, sans avoir besoin de
faire une correction :

N T = 25 nappes
5.
5a). Perte de charge
Lorsquun circuit possde des branches en parallle, la perte de charge X est la
mme sur chacune des branches et sur lensemble (FEMM, 7.6.1).
Dans une passe, o chaque tube a une longueur H c , daprs FEMM (6.140b) la perte
V f2
Hc
de charge en ligne vaut : 4 C f

, et la perte de charge sur chaque singularit slve

d
2
V f2
:
2
De lentre la sortie de leau, le circuit comporte N T passes ; chaque particule fluide
suit un itinraire qui comporte donc N T tubes, N T 1 coudes de demi-tour, une entre et une
sortie, avec une vitesse uniforme V f , do une perte de charge totale (FEMM. 6.177 et
6.178) :
2
NT H c

Vf
X = 4 C f
+ ( NT 1 ) + e + s
d
2

On a vrifi la question 3a) que lallongement H c / d des tubes est suprieur 60.
Vu le nombre de Reynolds Re f = 54400 , le coefficient de frottement C f est donn par la
formule de Colebrook (4.23 et aussi FEMM 7.10) :

0 , 883
1

= 2 , 5 Ln 0 , 285
+

D Re C f / 2
Cf / 2

Le calcul est itratif. Dans la pratique, les valeurs des coefficients de frottement vont
peu prs de 10 1 10 4 . Comme point de dpart, on peut prendre une valeur moyenne
C f = 2.10 3 dans le membre de droite :

0 , 883
10 2

= 2 , 5 Ln 0 , 285
+

3
31
Cf / 2
54400 2.10
/2

et on obtient comme premire estimation :

C f = 5 , 83.10 3
On reprend le calcul prcdent en remplaant la valeur de dpart 2.10 3 par cette
nouvelle valeur, ce qui conduit une seconde estimation :

C f = 5 , 20.10 3
Rinjectons cette nouvelle valeur dans le membre de droite. Le rsultat est :
C f = 5 , 26.10 3
On est prs de la convergence. Une ultime itration donne :
C f = 5 , 25.10 3

La perte de charge est alors :

25 2 ,615

3 12
3
X = ( 4 5 , 25.10

) + ( 24 0 , 4 ) + 0 , 5 + 1 10
2
31.10 3

10 3
2
X = 27700 Pa = 0 , 277 bar
= ( 44 , 28 + 11,1 )

5b). Pompe
La puissance utile de la pompe doit tre gale la puissance P dissipe dans
lcoulement (FEMM, 7.4.2 et 7.4.4) :
Pu = P = q v X

Cette formule est encore valable pour un circuit avec branches en parallle, q v tant le
dbit total (FEMM, 7.6.1) :
Pu = 12.10 3 27700
Pu = 332 , 4 W

Commentaires
Ce problme prend en compte de nouveaux lments : changeur multipasses,
montage des tubes en parallle, gomtrie du faisceau, perte de charge .
- Un calcul approch de la surface dchange pourrait tre conduit partir du bilan
dnergie sur le fluide froid, en appliquant (6.5), soit : = k f f ( < Tc > < T f > ) . On a
< Tc > = Tmc = 220 C , < T f > = Tmf = 50 C , soit f = 3010.10 3 / [177 ,5 ( 220 50 )] ,

ce qui donne une surface dchange de 99 ,75 m 2 , au lieu de 101,7 m 2 par la mthode NUT
(question 4a). Lcart nest que de 2%. Doit-on en dduire que cette mthode est un
raffinement inutile ? Certainement pas, dabord parce que lcart nest pas toujours aussi
modeste (ici, le NUT est petit, cf. Problme 6), ensuite parce que, en accumulant les
approximations, on finit par largir sensiblement la marge dincertitude.
- Le calcul itratif du coefficient de frottement converge toujours rapidement, mme si
lestimation de dpart est un peu loin du rsultat final.
- On peut tre surpris par la faible puissance de la pompe, due en particulier au fait que
la vitesse de circulation nest pas trs leve (1 m/s). Mais il faut prendre aussi en compte son
rendement (rarement suprieur 0,7) et tout le circuit dans lequel lchangeur est insr. Au
bout du compte, on arriverait facilement une puissance relle de 1 2 kW.

PROBLME N 10 : Condenseur
nonc
On tudie un condenseur de butane, destin produire 15 kg/s de butane liquide
partir dune vapeur saturante la temprature Tv = 50 C et la pression p = 5,08 bars.
Lappareil utilis est un changeur faisceau de tubes horizontal, quatre passes
ct tubes. Ceux-ci ont un diamtre intrieur d = 15 ,75 mm et un diamtre extrieur
D = 19 ,05 mm .
On dispose comme fluide froid dune eau de rivire dont la temprature maximale est
de 15 C en t, et qui ne doit pas subir un chauffement suprieur 10 C. Elle circule
lintrieur des tubes.
Les caractristiques physiques du butane p = 5,08 bars sont donnes ci-dessous :

Temprature (C)
Chaleur latente de vaporisation LV (J/kg)
Conductivit thermique phase liquide (W/m K)
Viscosit dynamique phase liquide (kg/m s)

20

0 ,17.10 3

50
335.10 3
0,125
0 ,13.10 3

Masse volumique phase liquide ( kg / m 3 )

579

542

0,131

1. Calculer le dbit de fluide froid q mf ncessaire dans lchangeur.

2. Quel est le dbit-masse deau q' mf dans chaque tube, le nombre de tubes N
ncessaire pour chaque passe et le nombre total N t de tubes ? (On se fixe une vitesse
dbitante V = 2 m/s dans les tubes).
3. Dterminer le coefficient dchange h f ct fluide froid.
4. On se base a priori sur une temprature moyenne de paroi ct fluide chaud
T pc 30 C . Dterminer quel est alors le coefficient dchange moyen hc lextrieur des
tubes (raisonner ici comme pour un tube seul).
5. En tenant compte de la disposition des tubes en nappes verticales (qui diminue hc )
et de la rsistance thermique de la paroi, on admet finalement un coefficient global dchange
ct butane : k c = 730 W / m 2 K .
Calculer la surface dchange ncessaire. En dduire la longueur du faisceau de
tubes.
6. La rugosit moyenne lintrieur des tubes est = 0 ,01 mm . Les coefficients de
pertes de charge singulires sont donns : entre = sortie = 1, 2 pour chaque tube, compte
tenu de la forme des ttes de calandre.
Calculer la perte de charge dans le circuit deau et la puissance de la pompe qui doit
assurer la circulation (en admettant pour celle-ci un rendement global = 0 ,6 ).

Solution
1. Dbit deau
La puissance thermique de lchangeur se calcule en crivant le bilan sur le fluide
chaud :

= q mc LV = 15 335 10 3

5000.10 3 W
Sachant que lon doit avoir T fs T fe 10 C , le bilan sur le fluide froid donne, dans
la condition la plus dfavorable :
= qtf ( T fs T fe ) = 10 qtf
qtf = / 10 = 500.10 3 W / K = q mf C pf

Dans la gamme de temprature considre pour leau, C pf 4180 J / kg K . Do le

dbit minimal ncessaire :
q mf = 500.10 3 / 4180
q mf 120 kg / s

2. Nombre de tubes
Le dbit-masse dans un tube est :
d2
q' mf = f S V = f
V
4

V = 2 m/s

( 15 ,75.10 3 ) 2 2
4
= 0 , 39 kg / s

q' mf = 10 3
q' mf

avec

Dans chaque passe, les N tubes sont videmment monts en parallle, et le dbit total
scrit :
q mf = N q' mf
do : N = 120 / 0 , 39
N = 308 tubes par passe
et avec quatre passes il y a au total :
N t = 1232 tubes

3. Coefficient dchange ct froid

Pour calculer h on passe par lintermdiaire du nombre de Stanton, qui sexprime en
fonction de Re et Pr.
Conformment au 4.1.5, le nombre de Reynolds dans les tubes sera valu la
temprature moyenne de mlange de leau :
1
Tmf = ( T fs + T fe ) avec T fe = 15 C et T fs = ( 15 + 10 ) C
2
Tmf = 20 C

Daprs (4.13) : Re f =

Vd

Les tables donnent la viscosit cinmatique de leau 20 C :

f = 0 ,101.10 5 m 2 / s
Re f =

2 15 ,75 10 3
5

= 31200 : lcoulement est turbulent lisse ( Re < 10 5 ) ;

0 ,101 10
on calculera donc St f avec la formule de Dittus - Boelter (4.21) pour le fluide froid :

St f = 0 ,023 Re 0 ,2 Pr 0 ,6
Avant dentamer le calcul, assurons-nous quand mme que les conditions de validit
de la formule sont respectes :
- nombre de Prandtl : pour leau 20 C, Pr 7 , donc Pr > 0 ,7 . Cest correct.
- allongement des tubes L / d > 60 ? On ne sait pas. Admettons-le provisoirement. Nous
verrons la fin.
Donc :
St f = 0 ,023 ( 31200 ) 0 ,2 7 0 ,6
St f 9.10 4

et aussi en revenant la dfinition (4.6) : St f =

hf

CpV

Sagissant de leau : = 10 3 kg / m 3 et C p = 4180 J / kgK

h f = 9.10 4 10 3 4180 2
h f = 7525 W / m 2

4. Coefficient dchange ct chaud (butane)

La condensation se produit lextrieur dun tube horizontal ; le coefficient moyen
nous est donn par (5.1) :
0 ,25

g l 3l LV

hc = 0 ,725
l D ( Tv T p )

Nous avons comme temprature moyenne de film :

Tv + T p 50 + 30
TF =
=
= 40 C
2
2
A cette temprature, les caractristiques thermophysiques du butane sont obtenues en
interpolant les valeurs donnes dans lnonc ; on trouve :

l = 554 kg / m 3
l = 0 ,127 W / m K

l = l = 0 , 26.10 6 m 2 / s
l
En outre :
LV = 345.10 3 J / kg
Tv = 50 C ; T p = 30 C (valeur admise)

D = 19 ,05.10 3 m
et le calcul donne :
hc 1800 W / m 2 K

5. Dimensions de lchangeur
La surface dchange correspondant la valeur donne de k est la surface ct fluide
chaud, cest--dire lextrieur des tubes. On a donc, daprs (6.5) et avec le rsultat de la
question 1 :

= k c c ( Tv Tmf ) = 730 c 30 = 5000.10 3

c = 228 ,3 m 2
c est la surface extrieure des 1232 tubes, soit en appelant L la longueur dun tube
(qui est aussi la longueur du faisceau) :

c = N t D L et de ce fait

L=

228 ,3
1232 19 ,05.10 3

L = 3 ,10 m

Nous pouvons maintenant vrifier que la formule donnant St a t correctement

employe dans la question 3 : nous avons ici L / d 200 , donc trs suprieur 60.

6. Perte de charge
Dans une passe, les tubes sont monts en parallle : la perte de charge X p est donc
la mme pour un tube et pour la passe (FEMM, 7.6.1 et formule 6.178) soit :
2
L

V
X p = 4 C f
+
d
2

o est la somme des coefficients de pertes de charge singulires. Dans ce problme,

= 1, 2 + 1, 2 = 2 ,4 .
Compte tenu de la valeur de Re, le coefficient de frottement est donn par la formule
de Blasius (FEMM . 7.8a) :
C f = 2 0 ,0395 Re 0 ,25 = 0 ,006
On peut donc se passer de la rugosit . Cest une donne surabondante.
Pour les quatre passes, la perte de charge totale stablit :
2

3 ,10
10 3 2
X = 4 4 0 ,006
+
2
,
4

2
15 ,75.10 3

X 57000 Pa = 0 ,57 bar

Daprs FEMM (6.167a) et (7.73b), la puissance dissipe dans le circuit est :
P = q v X
En admettant que lentre et la sortie du circuit sont au mme niveau, la puissance
utile fournie par la pompe doit tre gale la puissance dissipe dans le circuit. On a donc
partir de FEMM ( 7.4.2, 7.5.2.1, et relation 7.55) :

P=

Pu

q v X

avec

qv =

q mf

= 0 ,12 m 3 / s

et

= 0 ,6

P = 11400 W

Commentaires
1. Lnonc contient certaines informations qui ont pour but de faciliter la rsolution
du problme, mais qui doivent normalement tre estimes par lingnieur : il sagit du
coefficient moyen dchange k c et de la temprature moyenne de paroi T p .
En fait, ces deux grandeurs sont interdpendantes puisque la formule qui donne k c
contient T p . L comme dans dautres exemples cits, cest une procdure itrative qui est
recommander :
- premire estimation de T p
- calcul de hc , k c et du nombre de tubes
- choix de la disposition du faisceau et en particulier du nombre de tubes dans une nappe
verticale, et rvaluation de k c par (5.2)
- nouvelle valuation de T p et ritration du calcul. Souvent, la convergence est obtenue ds
cette seconde tape.
2. Pour le calcul de C f , la formule de Colebrook (FEMM, 7.10) est plus sre, mais
Blasius est plus rapide et donne ici un rsultat trs voisin.
3. La valeur de h relative la condensation du butane est nettement plus faible que ce
que lon obtiendrait avec de leau.

PROBLME N 11 : changeur compact

nonc
On aborde ltude davant-projet dun changeur compact courants croiss, avec un
fluide brass.
Lchangeur est constitu dune batterie de tubes parallles, relis par des ailettes
perpendiculaires aux tubes et jouant le rle dentretoises (fig. 1: vue de face, et fig. 2: coupe
longitudinale dun tube). On admettra que les tubes sont de section rectangulaire (fig. 3 :
coupe transversale dun tube), disposs en quinconce (fig. 4 : vue de dessus).
Il sagit dun changeur air-eau. Leau circule dans les tubes, lair autour des tubes,
entre les ailettes. Les dimensions du caisson sont : 1 m (largeur) 1 m (hauteur) 0,5 m
(profondeur). Les autres dimensions sont portes sur les figures, en mm. Lappareil comporte
1786 tubes et 500 ailettes. Les caractristiques aux entres sont les suivantes :
- pour leau : T fe = 12 C ; q mf = 28 kg / s
- pour lair : Tce = 125 C ; q mc = 14 kg / s

FIG. P11.1

FIG. P11.3

FIG. P11.2

FIG. P11.4

1. Dterminer le coefficient dchange h f ct eau.

2.
2.a) Calculer le diamtre hydraulique pour lcoulement de lair (la section prendre en
compte est la section frontale de passage de lair, voir fig. 1).
2.b) La section totale de passage de lair est S c = 0 ,782 m 2 . Evaluer le nombre de Reynolds
Rec ct air. Quelle est linfluence de la pression de lair ?
2.c)

Labaque de la figure 5 donne St Pr 2 / 3 en fonction de Rec . La temprature de

rfrence est la temprature moyenne de mlange Tmc . Calculer le coefficient dchange hc

ct air. Comment est-il influenc par la pression de lcoulement ?

FIG. P11.5
3. La surface dchange c ct fluide chaud est de 489 m 2 , dont 418 m 2 pour les
ailettes. On admet pour ces dernires une efficacit = 0 , 86 . Dterminer la conductance
globale K de lchangeur (on nglige la rsistance thermique des parois).
4. Calculer lefficacit de lchangeur
5. Calculer les tempratures de sortie pour lair et leau. Vrifier si les tempratures
moyennes adoptes en cours de calcul sont suffisamment proches des tempratures calcules.

Solution
1. Coefficient dchange ct eau
La premire chose faire est destimer approximativement les tempratures de
rfrence pour valuer ensuite les caractristiques thermophysiques des fluides. Sil savre
la fin du problme que cette estimation est trop imparfaite, on recommencera en se basant
cette fois sur les rsultats du premier calcul.
Ecrivons dj le bilan enthalpique global :
q mc C pc ( Tce Tcs ) = q mf C pf ( T fs T fe )

Les chaleurs massiques varient peu avec T. Prenons C pc = 1007 J / kg K pour lair et
C pf = 4180 J / kg K pour leau, soit compte tenu des valeurs imposes aux dbits :
14 1007 Tc = 28 4180 T f

ce qui donne

Tc
8
T f

Un ordre de grandeur plausible serait Tc 90 C et T f 11 C , ce qui respecte

peu prs le rapport de 8. Cela correspondrait des tempratures de sortie Tcs 35 C et
T fs 23 C , do des tempratures moyennes de mlange ( 4.1.5) :
125 + 35
12 + 23
= 80 C ; Tmf
= 17 , 5 C
2
2
Nous partirons sur ces bases-l.
Tmc

La section du tube est rectangulaire ( 4.2.1.3, 2 b = 2 mm , l = 18 mm ) et Re f est

donn par (4.17) :
V f Dhf
2 2 b l 2 2 18
Re f =
o Dhf =
=
= 3 ,6 mm
f
l + 2b
18 + 2
La section intrieure dun tube est : S = 18 2 = 36 mm 2 . Elle est parcourue par un
dbit :
q' vf =

q vf

28.10 3
= 15 ,68.10 6 m 3 / s
1786

nb. de tubes
do une vitesse :
15 ,68.10 6
Vf =
= 0 , 435 m / s
36.10 6
Re f est calcul la temprature moyenne de mlange. A Tmf = 17 , 5 C , les tables
numriques donnent pour la viscosit de leau :

f = 1,06.10 6 m 2 / s
On en dduit : Re f =

3,6.10 3 0 ,435
1,06.10 6

Re f 1480
Lcoulement est donc laminaire.

La formule utiliser pour calculer St est en principe (4.18) ou (4.19). Pour trancher, il
faut comparer L / Dhf et 0 ,014 Pe . Daprs (4.17) :
Pe f =

V f Dhf
af

avec a f 0 ,142.10 6 m 2 / s pour leau Tmf = 17 , 5 C , ce qui donne :

Pe f 11.10 3

Ici, L = 1 m , soit :

L
1000
=
0 ,025 > 0 ,014
Dhf Pe f
3 ,6 11000

Cest donc la formule tubes longs (4.19) qui sapplique, sous rserve que les
conditions sur Re, Pr et / p soient bien vrifies, ce qui est le cas :
0 ,14

En ce qui concerne / p , on doit sattendre une temprature moyenne de paroi T p

8 , 23
St f =
Pe f

assez proche de Tmf puisque le coefficient dchange ct eau sera en principe beaucoup plus
lev que ct air. En admettant T p = Tmf + 10 C , on trouve / p = 1,07 / 0 , 88 do une
correction ( / p ) 0 ,14 1,03 que lon pourrait la rigueur ngliger.

On obtient finalement :
8 , 23
St f =
1,03
11.10 3
St f 0 , 77.10 3 , et alors :
h f = St C p V = 0 ,77.10 3 10 3 4180 0 , 435
h f 1400 W / m 2 K

2.
2.a) Le nombre de Reynolds se calcule dans un canal lmentaire, dont la section est
11, 5 1, 9 mm (fig. 1). Le diamtre hydraulique est donn par (4.17) :

2 2 b l 2 1, 9 11,5.10 6
Dhc =
=
l + 2b
( 11,5 + 1,9 ) 10 3
Dhc = 3, 26 mm
2.b) Nous avons ici : Rec =

Vc Dhc

et c =

c
c

En admettant une bonne quirpartition des dbits dans les canaux, la vitesse est la
mme partout ; elle peut tre calcule directement partir du dbit et de la section totale :
q
Vc = mc , soit :
c Sc
q D
Rec = mc hc
c Sc
On ne nous prcise pas la pression moyenne de lair, mais cela est sans importance car
la viscosit dynamique dpend trs peu de la pression (proprit indique dans les tables
numriques). Comme cest le dbit-masse q mc qui est fix, Rec est donc indpendant de la
pression.

2.c) Le calcul doit tre fait la temprature moyenne de mlange, soit Tmc 80 C . Alors :

c = 2 ,08.10 5 kg / m . s
Rec =

14 3 , 26.10 3

2 ,08.10 5 0 ,782
Rec 2800

Labaque donne approximativement :

St c Pr 2 / 3 = 6 , 5.10 3
Pour lair, le nombre de Prandtl est peu dpendant de T : Pr 0 ,7
St c = 6 , 5.10 3 0 ,7 2 / 3
St c = 8 , 24.10 3

On a aussi :
hc
St c =
c C pc Vc

avec

c Vc =

q mc
= Gc vitesse massique (4.10)
Sc

De ce fait nous navons pas besoin de c et le coefficient hc ne dpend pas de la

pression :
q
hc = St c C pc mc avec C pc = 1007 J / kg (question 1)
Sc
Adoncque :
8 , 24.10 3 1007 14
hc =
0 ,782
hc = 149 W / m 2 K

3.
Les ailettes sont situes ct fluide chaud. La conductance globale est donne par (6.8)
et (6.6a) :

1
1
K =
+
hc ( L + a ) f h f

avec a = 418 m 2 ; L = c a (surface latrale des tubes) = 71 m 2 ; = 0 ,86 .

Ct eau, la surface dchange est la surface intrieure des 1786 tubes :
f = primtre 1 m 1786 = 2 ( 18 + 2 ) 10 3 1786

f = 71,45 m 2
Revenons K :

1
1

K =
+
149 ( 71 + 0 ,86 418 ) 71,45 1400
K = 39370 W / K

4.
En revenant la question 1, on trouve pour les dbits thermiques unitaires :
qtc = 14 1007 = 14100 W / K
qtf = 28 4180 = 117.10 3 W / K

Donc, qt min = qtc = 14100 W / K

Lair chaud qui circule autour des tubes est brass. Si lon se reporte au tableau 3.1, on
est la rubrique changeur courants croiss, un fluide brass qt min :
E = 1 exp ( / R ) avec
qt min
14100
R=
=
= 0 ,12
qt max 117000

= 1 exp ( R NUT )

NUT est alors calculable par (3.14b) :

K
39370
=
NUT =
qt min 14100
NUT = 2 ,79
= 1 exp ( 0 ,12 2.79 ) = 0 , 285

0 , 285

E = 1 exp
0 ,12
E = 0 , 907
5.
5.a) Temprature de sortie dair
Puisque qtc = qt min , lefficacit sexprime selon (3.10) :
E=

Tce Tcs
Tce T fe

On en tire :
Tcs = Tce E ( Tce T fe ) = 125 0 , 907 ( 125 12 )
Tcs = 22 , 5 C
5.b) Temprature de sortie deau :
Le bilan enthalpique nous dit que :
qtc ( Tce Tcs ) = qtf ( T fs T fe )
14700
( 125 22 , 5 )
117000
T fs = 24 , 5 C

T fs = 12 +

5.c) Vrifications
- Pour la temprature moyenne de leau, nous avions choisi 17,5 C. Avec les valeurs
calcules elle est de ( 12 + 24 , 5 ) / 2 = 18 , 2 C . La diffrence est insignifiante.

- Pour la temprature moyenne de lair, lestimation avait donn 80 C. Nous sommes

( 125 + 22 , 5 ) / 2 74 C . Lcart na pas de rpercussion significative sur les constantes
thermophysiques.
- Quant lestimation de T p (question 1), son ordre de grandeur est compatible avec
les valeurs calcules de hc et h f . Sa rvision ne changerait pas grand-chose la correction en

/ p qui reste faible de toute faon.

Commentaires
Le problme sinspire dune application prsente dans louvrage de J.F. Sacadura
[Initiation aux transferts thermiques].
Lappareil est assez performant. Son efficacit est leve, et on observe que la
temprature de sortie du fluide froid (24,5 C) est suprieure la temprature de sortie du
fluide chaud (22,5 C), situation trs favorable (cette ventualit a t signale dans le cas des
changeurs contre-courant, 2.3.5). En pratique, cet excellent rsultat est probablement
nuancer un peu, du fait notamment que les dbits ne sont jamais parfaitement quirpartis
(solution question 2.b) dans les tubes ou dans les diffrents chemins de passage de lair.
Concernant la section rectangulaire des tubes, on notera que malgr un aplatissement
marqu ( 2 / 18 mm ) lapproximation 2 b ` l est viter : elle donnerait un diamtre
hydraulique de 4 b = 4 mm alors que lon trouve Dhf = 3 ,6 mm avec la formule complte
(question 1). Mme remarque pour la section de passage de lair (question 2a). On peut
signaler galement que, en ralit, les extrmits des tubes seront arrondies pour amliorer
lcoulement et diminuer les pertes de charge ct air.

PROBLME N 12 : Assemblage de deux changeurs

nonc
On dispose de deux changeurs identiques et de deux sources deau (une source
chaude et une source froide) dont les tempratures Tce , T fe et les dbits q mc , q mf sont
donns.
Lobjet du problme est de choisir, entre plusieurs dispositions des changeurs, celle
qui permet dobtenir la valeur de T fs la plus leve.
1. Avec un changeur seul, montrer que quelque soit son modle, si qt min = qtc , on a :
Tcs T fe
Tce T fs

1 E
1 RE

2. Les deux changeurs sont placs en srie comme lindique la figure P12.1.
2.1. En raisonnant sans faire de calculs, montrer que les deux appareils ont la mme
efficacit (dsigne par E). On admet que linfluence de la temprature sur les constantes
thermophysiques est ngligeable.
2.2. Soit Et lefficacit de lchangeur quivalent lassemblage ralis. En admettant
qt min = qtc , et en utilisant le rsultat de la question 1, exprimer Et en fonction de R et de E.
Dterminer T fs en fonction de Tce , T fe , Et et R.

FIG. P12.1

3. Les changeurs sont maintenant disposs en srie sur le fluide chaud et en

parallle sur le fluide froid (fig. P12.2). On ne prcise pas ici quel est le qt minimum.
3.1. Montrer sans calculs que les deux appareils ont la mme efficacit ( E' c ct chaud ou
E' f ct froid).
3.2. En raisonnant sur le circuit srie, cest--dire ct fluide chaud, exprimer Tcs puis T fs
en fonction de E' c et des donnes.
4. Ceci nest pas une question.
En permutant les circuits (fig. P12.3) on obtient :
T fs = Tce E' f ( 2 E' f ) + T fe ( 1 E' f )2

FIG. P12.2

5. Les deux appareils sont des changeurs courants croiss, avec fluide chaud
brass une seule passe, de surface dchange = 170 m 2 . On a dtermin
exprimentalement le coefficient global dchange k. Pour q mc > 4 kg / s et q mf < 10 kg / s :
k = 52 ( q mf + 0 , 8 q mc ) 0 ,5 + 65

a) q mf

en W / m 2 K

On donne : Tce = 90 C ; T fe = 10 C
Dterminer lefficacit dun changeur seul dans les trois cas suivants :
= 9 ,6 kg / s ; q mc = 8 , 2 kg / s

b) q mf = 4 , 8 kg / s ; q mc = 8 , 2 kg / s
c) q mf = 9 ,6 kg / s ; q mc = 4 ,1 kg / s
Calculer T fs dans le cas a).

FIG. P12.3

6. On rappelle les deux sources deau dont on dispose :

- eau chaude : Tce = 90 C ; q mc = 8 , 2 kg / s
- eau froide : T fe = 10 C ; q mf = 9 ,6 kg / s
Calculer T fs dans les trois dispositions suivantes (on dsignera avec le symbole les
caractristiques relatives un seul changeur) :
I - en mettant les changeurs en srie
II - en les mettant en srie sur leau chaude et en parallle sur leau froide (le fluide brass
est le fluide chaud)
III - en les mettant en srie sur leau froide et en parallle sur leau chaude (le fluide brass
tant toujours le fluide chaud).
Conclure

Solution
1.
Sachant que qt min = qtc , daprs (3.10) lefficacit et le facteur de dsquilibre
scrivent :
E=

Tce Tcs
Tce T fe

R=

qtc T fs T fe
=
qtf
Tce Tcs

On vrifie alors aisment que :

T Tcs
1 ce
Tce T fe
Tcs T fe
1 E
=
=
T fs T fe Tce Tcs
1 RE
Tce T fs
1
Tce Tcs Tce T fe

CQFD

2.
2.1) Les deux appareils sont identiques : ils ont en particulier la mme surface dchange .
Etant traverss par les mmes dbits chauds et froids, ils ont aussi le mme coefficient global
dchange k, et par consquent le mme R et le mme NUT. Leurs efficacits sont donc
gales.
2.2) Daprs la question 1, avec les notations de la figure 1 :
Tc1 T f 2
1 E'
=
Tce T fs
1 R E'
Tcs T fe
Tc1 T f 2

1 E'
1 R E'

Effectuons le produit de ces quations membre membre :

Tcs T fe
( 1 E' )2
=
Tce T fs ( 1 R E' )2
Si lon considre lchangeur quivalent lensemble, pour lequel qt min est toujours
gal qtc , on a aussi :

Tcs T fe
Tce T fs

1 Et
1 R Et

Des deux relations prcdentes on tire :

Et =

( 1 R E' ) 2 ( 1 E' ) 2
( 1 R E' )2 R ( 1 E' ) 2

Puisque qt min = qtc dans le rseau, lefficacit de lchangeur quivalent scrit :

Et =

Tce Tcs
Tce T fe

do :

Tcs = Tce Et ( Tce T fe )

Nous disposons galement du bilan global du rseau :

qtf ( T fs T fe ) = qtc ( Tce Tcs )
do lon peut extraire T fs :
T fs = T fe +

qtc
( Tce Tcs )
qtf

En remplaant Tcs par son expression ci-dessus, on obtient la temprature de sortie du

fluide froid :
T fs = T fe + R Et ( Tce T fe )

3.
3.1) Le raisonnement est analogue celui de la question 2 : les changeurs sont identiques et
travaillent dans les mmes conditions (dbits, coefficients dchange etc). Par consquent ils
ont le mme R et le mme NUT, donc la mme efficacit.

3.2) Lefficacit ct fluide chaud scrit, selon lchangeur considr :

T Tc1 Tc1 Tcs
E' c = ce
=
Tce T fe Tc1 T fe

Ces deux quations constituent un systme deux inconnues Tc1 et Tcs . Tirons Tc1 de
lune delles :
Tcs E' c T fe
Tc1 = E' c ( Tc1 T fe ) + Tcs
soit
Tc1 =
1 E' c
et reportons dans lautre :
Tcs E' c T fe
E' c ( Tce T fe ) = Tce Tc1 = Tce
1 E' c
On obtient ainsi la temprature de sortie chaude :
Tcs = Tce E' c ( 2 E' c ) ( Tce T fe )

Le bilan enthalpique global du rseau donne ensuite la temprature de sortie froide :

qtf ( T fs T fe ) = qtc ( Tce Tcs )
En remplaant Tcs il vient :
T fs = T fe +

qtc
E' c ( 2 E' c ) ( Tce T fe )
qtf

4.
La question nest pas pose mais le calcul serait analogue en raisonnant encore sur le
circuit srie, cest--dire cette fois-ci sur le circuit froid.
5.
5.a) La chaleur massique de leau dpend trs peu de la temprature ; de ce fait les C p sont
pratiquement identiques sur les deux circuits, et lon voit que :
qt min = qt du fluide brass = qtc
Il y a une seule passe sur le fluide brass. Du tableau 3.1 on tire alors :
E = 1 exp ( / R )
avec = 1 exp ( R NUT )
et toujours :
NUT = k / qt min
Les donnes permettent de calculer k :

k = 52 ( 9 ,6 + 0 ,8 8 , 2 ) 0 ,5 + 65
k = 274 W / m 2 K
Dans la plage de temprature indique on prendra pour leau :
C p = 4190 J / kg K
qt min = q mc C p = 8 , 2 4190 = 34360 W / K
On en dduit NUT et R :
274 170
NUT =
= 1, 35
34360
q
q
8,2
R = tc = mc =
= 0 , 854
qtf
q mf
9 ,6
De l on tire :
= 1 exp ( 0 ,854 1,35 ) = 0 ,684
E = 1 exp ( 0 ,684 / 0 , 854 )
E 0 , 55

En se rfrant (3.10) on obtient enfin :

1
1 T fs T fe
do :
E = Ec =
Ef =
R
R Tce T fe

T fs = T fe + R E ( Tce T fe ) = 10 + 0 , 854 0 , 55 80
T fs = 47 ,6 C
5.b) Maintenant, qt min = qt du fluide non brass = qtf . Le tableau 3.1 nous dit que :
1
[1 exp ( R )] avec
R
Nous avons donc cette fois :
E=

= 1 exp ( NUT )

k = 52 ( 4 ,8 + 0 ,8 8 , 2 ) 0 ,5 + 65 = 240 W / m 2 K
qt min = q mf C p = 4 ,8 4190 = 20110 W / K
NUT =
R=

k
240 170
=
= 2 ,03
qt min
20110

q mf

4 ,8
= 0 , 585
8 ,2

q mc
ce qui donne :
= 1 exp ( 2 ,03 ) = 0 ,869
1
E=
[1 exp ( 0 ,585 0 ,869 )]
0 , 585
E 0 ,68
5.c) La situation est la mme que dans 5.a) : qt min = qt du fluide brass = qtc .

La succession des calculs est analogue :

k = 52 ( 9 ,6 + 0 , 8 4 ,1 ) 0 ,5 + 65 = 251,5 W / m 2 K
251, 5 170
NUT =
= 2 ,5
4 ,1 4190
4 ,1
R=
= 0 , 43
9 ,6
= 1 exp ( 0 ,43 2 ,5 ) = 0 ,658
et enfin :
E = 1 exp ( 0 ,658 / 0 , 43 ) 0 ,78
E = 0 ,78

6.
6.I. Montage en srie
Avec les valeurs numriques donnes, cette disposition correspond au cas 5.a). On a
donc pour chaque changeur :
E' = 0 , 55 ; R' = 0 , 854
Daprs la question 2, lefficacit de lchangeur quivalent a pour valeur :
( 1 0 ,55 0 ,854 )2 ( 1 0 ,55 )2
Et =
( 1 0 ,55 0 ,854 )2 0 ,854 ( 1 0 , 55 )2

Et = 0 ,73
et toujours daprs 2. :
T fs = 10 + 0 , 854 0 ,73 80
T fs 60 C
6.II. Srie sur circuit chaud et parallle sur circuit froid
Les donnes correspondent cette fois au cas 5.b).
q mc = 8 , 2 kg / s
Le dbit froid dans chaque changeur est gal la moiti du dbit total :
9 ,6
1
q' mf = q mf =
= 4 , 8 kg / s
2
2
et par consquent :
q' t min = q' tf

Daprs 5.b) on a donc :

R' = 0 , 585
E' = 0 ,68
A partir de la question 3, on a :
q
T fs = T fe + tc E' c ( 2 E' c ) ( Tce T fe )
qtf
o qtc et qtf sont les dbits globaux dans le rseau.
Pour calculer T fs , nous avons besoin de E' c (efficacit dun changeur ct circuit
srie, cest--dire circuit chaud). Mais ici E' = E' f puisque qt min = q' tf et :
E' c = R' E' f = 0 , 585 0 ,68 = 0 , 4

donc :
T fs = 10 +

8,2 C p
9 ,6 C p

0 , 4 ( 2 0 , 4 ) 80

T fs = 53 ,7 C
6.III. Srie sur fluide froid et parallle sur fluide chaud
La situation est celle du cas 5.c). Pour chaque changeur :
1
q' mc = q mc = 4 ,1 kg / s
q' mf = 9 ,6 kg / s
2
E' = 0 ,78
La formule de la question 4 qui donne T fs fait intervenir E' f :
T fs = Tce E' f ( 2 E' f ) + T fe ( 1 E' f )2

Or, ici E' = E' c puisque qt min = q' tc . Donc E' f = R' E' c :
R' =

4 ,1 C p
q' tc
=
= 0 , 427
q' tf
9 ,6 C p

E' f = 0 , 427 0 ,78 = 0 , 333

et finalement :
T fs = 90 0 , 333 ( 2 0 , 333 ) + 10 ( 1 0 , 333 )2
T fs = 54 , 3 C

Conclusion :
- avec un changeur seul : T fs = 47 ,6 C
- avec deux changeurs en srie : T fs = 60 C
- les deux montages en srie-parallle sont presque quivalents : T fs = 53 ,7 ou 54 , 3 C .
Le montage en srie est ici nettement meilleur.

Commentaires
Le problme reprend sous une forme plus simple, avec deux changeurs seulement, les
raisonnements du chapitre 7 ( 7.2 et 7.3). Mais on ne doit pas, partir des rsultats obtenus,
dduire une rgle gnrale sur la supriorit de tel ou tel assemblage.
- Le fait de donner une loi phnomnologique k = f ( q ) permet de calculer k
rapidement et donc dabrger la solution.
- Les questions sans calcul (2.1. et 3.1.) droutent souvent les tudiants, mais elles
permettent de voir si les notions essentielles ont t comprises.
- Dans les questions 5.2. et 5.3., lerreur ne pas faire est de confondre le dbit total et
les dbits partiels sur les branches en parallle, ou de confondre E' c et E' f .
- On doit bien noter que sur le circuit en parallle, les tempratures de sortie des deux
changeurs ne sont pas les mmes ( Ts1 Ts 2 , 7.3). Pour calculer les efficacits en revenant
la dfinition / max il faudrait connatre ces tempratures. Cest pour cela quen srieparallle on raisonne sur le circuit srie avec les efficacits relatives.

PROBLME N 13 : changeur en pingle

nonc
On tudie un changeur tubulaire en pingle, selon le schma ci-dessous (un seul tube
a t reprsent en pointills). Ses caractristiques sont E, R et NUT. Le fluide chaud circule
lintrieur des tubes.

FIG. P13.1
Ce modle dchangeur peut tre considr comme quivalent lassemblage de deux
demi-changeurs de mmes dimensions, lun (N 1) fonctionnant co-courant, lautre (N 2)
contre-courant (fig. P13.2).

FIG. P13.2

Les donnes du problme sont :

- les dbits thermiques unitaires (avec la condition qtc <

1
qtf )
2

- les tempratures dentre Tce et T fe

- les diamtres d, D, et la longueur totale L de chaque tube ( L1 = L2 = L / 2 ). On admettra
que les tubes sont tous identiques, et on ngligera les problmes lis aux coudes qui assurent
le demi-tour. La calandre est isole sa surface externe.
1.

Justifier la relation suivante : T fs = ( T fs1 + T fs 2 ) / 2 . Exprimer R1 et R2 en

fonction de R , ainsi que NUT1 et NUT2 en fonction de NUT.

2. Calculer E1 et E 2 en fonction de R et de NUT. Calculer ensuite Tc 2 (temprature
de sortie chaude du demi-changeur N1) et Tcs , puis lefficacit E de lchangeur complet

en fonction de E1 , E 2 et R (comparer avec le rsultat de la thorie gnrale des rseaux,

ch.7) . En dduire T fs .
3. Lappareil est un changeur eau-eau. On donne :
Tce = 90 C ; T fe = 20 C ; qvf = 82 m 3 / h ; q vc = 19 , 80 m 3 / h
faisceau de 60 tubes en parallle ; L = 5 m ; d = 15 mm ; D = 18 mm
pas carr ST = S L = 28 mm ; diamtre de la calandre Dc = 0 , 40 m
Effectuer sans calculs une premire valuation des tempratures moyennes de
mlange Tmc et Tmf .
Calculer le diamtre hydraulique Dh et le nombre de Reynolds Re f de
lcoulement en calandre. Dterminer le coefficient dchange h f autour des tubes.
4. Calculer le coefficient dchange hc lintrieur des tubes. Dterminer le
coefficient global dchange k (on raisonnera comme si la paroi tait plane, et on ngligera
sa rsistance thermique).
5. Calculer le NUT et lefficacit E, puis les tempratures de sortie T fs et Tcs .

Solution
1.
1a). Le dbit froid dans chaque demi-changeur est la moiti du dbit total :
1
q mf 1 = q mf 2 = q mf
2
Donc la temprature de sortie T fs est la moyenne des deux tempratures de sortie T fs1
et T fs 2 .
1b). Dans le demi-changeur N1, on a : qtf 1 =

Mais on sait aussi que qtc <

1
qtf et qtc1 = qtc .
2

1
qtf . Donc qtc1 < qtf 1 . En consquence :
2

qt 1 min = qtc1 = qtc

R1 =
R1 =

qt 1 min
qt 1 max

qtc1
qtf 1

2 qtc
qtf

Pour lchangeur complet : qt min = qtc et R =

R1 = 2 R
Le NUT du demi-changeur N1 est :
k
NUT1 = 1 1
qt 1 min

qtc
.
qtf

Le coefficient dchange est indpendant du sens de circulation des fluides, donc

k 1 = k 2 = k coefficient dchange global de lchangeur complet.
La surface dchange 1 est la moiti de la surface totale .
k
NUT1 =
2 qtc
k
k
=
.
Le NUT de lchangeur complet vaut : NUT =
qt min
qtc
Do :
NUT1 =

1
NUT
2

Le raisonnement est identique pour le demi-changeur N2 : qt 2 min = qtc , et :

R2 = 2 R

NUT2 =

1
NUT
2

2.
2a). Efficacits
Le demi-changeur N1 fonctionne en co-courant. Daprs le tableau 3.1 :
1
E1 =
{1 exp [ ( 1 + R1 ) NUT1 ]}
1 + R1
E1 =

1
1 + 2R

1 + 2R

NUT
1 exp
2

Le demi-changeur N2 est contre-courant. Toujours daprs le tableau 3.1, son

efficacit est :
1 exp [ ( 1 R2 ) NUT2 ]
E2 =
1 R2 exp [ ( 1 R2 ) NUT2 ]
1 2R

1 exp
NUT
2

E2 =
1

2
R

1 2 R exp
NUT
2

2b). Tempratures de sorties chaudes

Dans le demi-changeur N1, qt 1 min = qtc . Alors, daprs (3.10) :
E1 =

Tce Tc 2
Tce T fe

Tc 2 = Tce E1 ( Tce T fe )

Dans le numro 2, cest pareil : qt 2 min = qtc . Mais ici, la temprature dentre chaude
est Tc 2 . De ce fait :

E2 =

Tc 2 Tcs
Tc 2 T fe

Tcs = Tc 2 E 2 ( Tc 2 T fe )

([

= Tce E1 ( Tce T fe ) E 2 Tce E1 ( Tce T fe ) T fe

Tcs = Tce ( Tce T fe ) ( E1 E1 E 2 + E 2 )

Dans lchangeur complet, qt min = qtc et lefficacit scrit :

E = Ec =

Tce Tcs
Tce T fe

E = E1 E1 E 2 + E 2
Les deux demi-changeurs sont assembls en srie-parallle, le circuit srie tant celui
du fluide chaud. Le rsultat est conforme aux expressions (7.29) 7.3.1, ou (7.32a) 7.3.3.
On a galement, toujours daprs (3.10) :
1
1 T fs T fe
E=
Ef =
R
R Tce T fe
do lon tire :
T fs = T fe + R E ( Tce T fe )

3.
3a). Tempratures moyennes
Vu le rapport des dbits, qui est voisin de 4, le T f de leau froide ( q max ) vaudra
environ du Tc de leau chaude. On peut admettre, pour commencer le calcul T f 10 C
et Tc 40 C , ce qui donnerait pour les moyennes :
Tmf 25 C ; Tmc 70 C
3b). Diamtre hydraulique de lcoulement en calandre
Le faisceau est longitudinal ( 4.3.2) ; Dh est donn par (4.31c) :

Dh = 4
Dh = 4

ST S L
D
D

(28.10 )

3 2

18.10

18.10 3 = 0 ,0554 0 ,018

Dh = 0 ,0374 m

3c). Nombre de Reynolds

V f Dh
avec V f = vitesse dbitante du fluide froid, quil faut calculer.
Re f =

La section de passage du fluide froid est :

S f = sec tion calandre sec tion tubes
Le faisceau comporte 60 tubes, mais comme ils sont en U il y a 120 tubes dans une
section droite de la calandre :
Dc2
D2
Sf =
120
= ( 0 ,40 2 120 0 ,018 2 )
4
4
4
2
S f = 0 ,095 m
La vitesse dbitante est dfinie par :
q vf
82
1
Vf =
=

Sf
3600 0 ,095
V f = 0 , 24 m / s

La viscosit de leau Tmf = 25 C est f = 0 , 92.10 6 m 2 / s .

Re f =

0 , 24 0 ,0374

0 , 92.10 6
Re f = 9760

3d). Coefficient dchange ct froid

Pour le nombre de Stanton externe des changeurs tubulaires courants parallles, on

dispose de la formule (4.32), valable si 5.10 3 < Re < 10 5 , soit avec le fluide froid :
St f = 0 ,026 Re f 0 ,18 Pr f 0 ,6

A 25 C, le nombre de Prandtl de leau est : Pr f = 6 , 4 .

St f = 0 ,026 9760 0 ,18 6 , 4 0 ,6
St f = 1,63.10 3

On calcule le coefficient dchange h f en revenant la dfinition de St :

h f = St f f C pf V f

et en prenant C pf = 4180 J / kg K pour leau :

h f = 1,63.10 3 10 3 4180 0 , 24
h f = 1635 W / m 2 K

4.
4a). Coefficient dchange interne
Dans chacun des 60 tubes, le dbit est :
q
19 ,80
qv0 = vc =
60
60 3600
q v0 = 0 ,0917.10 3 m 3 / s
do une vitesse de circulation :

Vc =

q v0

d2 / 4

4 0 ,0917.10 3

15.10 3

Vc = 0 , 52 m / s

Le nombre de Reynolds vaut :

V d
Rec = c

On a admis sur le circuit deau chaude une temprature moyenne Tmc = 70 C , do la

viscosit c = 0 , 42.10 6 m 2 / s .

Rec =

0 ,52 15.10 3

0 ,42.10 6
Rec = 18600
La formule de Dittus-Boelter (4.21) pour le fluide chaud sapplique si L / d > 60 , ce
qui est largement vrifi (5000/15) :
St c = 0 ,023 Rec 0 ,2 Prc 0 ,7
A Tmc = 70 C , Prc = 2 ,62
St c = 0 ,023 18600 0 ,2 2 ,62 0 ,7
St c = 1,64.10 3

Le coefficient dchange est :

hc = St c c C pc Vc
La chaleur massique varie peu avec T. On gardera celle du circuit froid.
hc = 1,64.10 3 10 3 4180 0 , 52
hc = 3565 W / m 2 K
4b). Coefficient dchange global
La formule 6.2a (paroi plane) scrit en labsence de rsistance :
1
1
1
1
1
=
+
=
+
= ( 0 , 28 + 0 ,612 ) 10 3
k hc
hf
3565 1635

k = 1120 W / m 2 K

5.
5a). NUT
Pour lchangeur complet, puisquon raisonne comme avec une paroi plane, on a
daprs (6.10 a et b) une conductance globale :
K = k m
avec
D+d
m =
L totale
2

Alors :
NUT =

K
qt min

k m
qtc

Pour les 60 tubes :

18 + 15

m =
10 3 5 60 = 15 , 56 m 2
2

Il faut encore le dbit thermique unitaire de leau chaude :

19 , 8
qtc = c q vc C pc = 10 3
4180
3600
qtc = 23.10 3 W / K
do :
NUT =

1120 15 , 56

23.10 3
NUT = 0 ,76

5b). Efficacit
Pour calculer E, il nous faut aussi le facteur de dsquilibre :
qt min
q
R=
= tc
qt max
qtf

Les deux fluides ont la mme masse volumique et la mme chaleur massique, donc :
q
19 , 8
R = vc =
qvf
82
R = 0 , 24

Les efficacits partielles E1 et E 2 ont t dtermines dans la question 1 :

1
1 + 2 0 , 24
E1 = 0 , 29

E1 =

1 + 2 0 , 24

0 ,76
1 exp
2

1 2 0 , 24

1 exp
0 ,76
2

E2 =
1 2 0 , 24

1 0 , 48 exp
0 ,76
2

E 2 = 0 , 297

On a montr que :
E = E1 E1 E 2 + E 2
E = 0 , 29 0 , 29 0 , 297 + 0 , 297
E = 0 , 50

5c). Tempratures de sortie

Toujours daprs la premire question :
T fs = T fe + R E ( Tce T fe ) = 20 + 0 , 24 0 , 50 70
T fs = 28 , 4 C
Tcs = Tce ( Tce T fe ) ( E1 E1 E 2 + E 2 ) = Tce E ( Tce T fe )
Tcs = 90 0 , 50 70
Tcs = 55 C

On doit enfin revenir sur lvaluation pralable des tempratures moyennes. Avec les
rsultats obtenus, on a :
20 + 28 , 4
Tmf =
= 24 ,2 C pour une estimation de 25 C : cest trs correct.
2
90 + 55
Tmc =
= 72 , 5 C pour une estimation de 70 C : cest assez prs.
2
Une itration nest donc pas indispensable.

Commentaires
Le problme combine un coulement longitudinal dans un faisceau de tubes et un
montage en rseau.
Les deux demi-changeurs ont des efficacits trs voisines, malgr leur diffrence de
catgorie (co-courant et contre-courant). Ceci est du au fait que le NUT est petit, ainsi que R
(voir linfluence de ces deux paramtres dans les exemples des figures 3.1 et 3.2). On observe
aussi que lefficacit totale est nettement suprieure chacune dentre elles.

PROBLME N 14 : Rseau 3 fluides

nonc
On considre un rseau ouvert mettant en uvre deux changeurs (fonctionnant
contre-courant) et trois fluides, avec un montage en srie sur le fluide froid (fig. P14.1).
Les tempratures dentre T fe , Tce1 , Tce 2 et les dbits thermiques unitaires des trois
fluides ( qtf , qtc1 , qtc 2 ) sont donns, ainsi que les efficacits E1 et E 2 des changeurs. On
impose galement qtf < qtc1 et qtf < qtc 2 .
1. Calculer les tempratures T f 1 , T fs , Tcs1 , Tcs 2 en fonction des donnes, puis le flux
total chang . A partir de lexpression de , prciser quelles conditions doivent satisfaire
les tempratures dentre pour que lensemble fonctionne de faon correcte.

FIG. P13.1
2. On se place dans le cas o qtc1 = qtc 2 . Les efficacits des deux changeurs sont
notes E A et E B .
Comparer les puissances globales changes AB et BA selon que lchangeur A
est plac en position 1 et lchangeur B en position 2, ou vice-versa (on tudiera le signe de
la diffrence AB BA ). En dduire les rgles dune disposition optimale pour le rseau
lorsque les tempratures dentre sont imposes.
3. Application pratique :
Le fluide froid est de lair, avec un dbit q vf = 800 m 3 / h . Les deux dbits chauds
sont gaux (cas de la question 2). On donne : T fe = 25 C ; Tce1 = 80 C ; Tce 2 = 65 C .
Les efficacits sont : E A = 0 ,38 ; E B = 0 , 25 .
Comparer les valeurs AB et BA selon la position relative des deux changeurs.
Calculer les tempratures de sortie froides T fs ( AB ) et T fs ( BA ) . Vrifier que les rsultats
sont en accord avec les conclusions de la question 2, et que la condition de fonctionnement
(question 1) est satisfaite.

Solution
1.
1a). Sorties de lchangeur N 1
Le dbit thermique unitaire minimum dans lchangeur est celui du fluide froid :
qtf < qtc1 , do lefficacit :
E1 = E f 1 =

T f 1 T fe
Tce1 T fe

T f 1 = T fe + E1 ( Tce1 T fe )
1
Ec , donc ici :
R
qt min
qtf
avec R1 =
=
qt max
qtc1

Dans ce cas, on a aussi daprs (3.11) : E f =

Ec1 = R1 E1 =

Tce1 Tcs1
Tce1 T fe

qtf

Tcs1 = Tce1

qtc1

E1 ( Tce1 T fe )

1b). Sorties de lchangeur N 2

Les conditions sont les mmes : qtf < qtc 2 . Le fluide froid entre T f 1 :
E2 = E f 2 =

T fs T f 1
Tce 2 T f 1

T fs = T f 1 + E 2 ( Tce 2 T f 1 )

En remplaant T f 1 par son expression obtenue en 1a) on obtient :

T fs = T fe + E 2 ( Tce 2 T fe ) + E1 ( 1 E 2 ) ( Tce1 T fe )

On a galement daprs (3.11) :

E c 2 = R2 E 2 =

Tce 2 Tcs 2
Tce 2 T f 1

Tcs 2 = Tce 2

qtf

Tcs 2 = Tce 2

qtf

avec

R2 =

qtf
qtc 2

E 2 ( Tce 2 T f 1 )
qtc 2
En remplaant T f 1 :
qtc 2

E 2 Tce 2 T fe E1 ( Tce1 T fe )

1c). Puissance
Le plus simple est de calculer la puissance du systme en faisant le bilan sur le fluide
froid (circuit srie) :
= qtf ( T fs T fe )

= qtf E 2 ( Tce 2 T fe ) + E1 ( 1 E 2 ) ( Tce1 T fe )

1d). Condition de fonctionnement

Pour un fonctionnement correct, dans chaque changeur, la temprature dentre
chaude doit tre suprieure la temprature dentre froide.
Tce1 > T fe

Tce 2 > T f 1
Ceci est vident pour le premier. Pour le second, on devra donc vrifier que :
Tce 2 > T fe + E1 ( Tce1 T fe )
Les conditions de fonctionnement de lchangeur N 2 doivent tre dfinies aprs
celles de lchangeur N 1.

2.
Montage avec A en 1 et B en 2 :
Daprs lexpression de obtenue dans la question 1 :
AB = qtf E B ( Tce 2 T fe ) + E A ( 1 E B ) ( Tce1 T fe )

Montage avec B en 1 et A en 2 :

BA = qtf E A ( Tce 2 T fe ) + E B ( 1 E A ) ( Tce1 T fe )

Aprs simplification, on obtient pour la diffrence :

AB BA = qtf ( E A E B ) ( Tce1 Tce 2 )

E A > E B et Tce1 > Tce 2

AB > BA si
E < E et T < T
B
ce1
ce 2
A
On en dduit une rgle doptimisation :

et vice-versa.

- Si Tce1 > Tce 2 : placer lchangeur le plus performant en position 1

- Si Tce1 < Tce 2 : placer lchangeur le plus performant en position 2
- Si Tce1 = Tce 2 ou E A = E B : AB = BA . La disposition est indiffrente.
3.
Il faut dabord dterminer le dbit thermique unitaire du fluide froid :
qtf = q mf C pf = f q vf C pf
La temprature moyenne de lair est inconnue. Une valeur de la masse volumique
f = 1,10 kg / m 3 (correspondant Tmf 37 C ) est plausible pour dmarrer le calcul.
Pour la chaleur massique, peu dpendante de la temprature, on prendra
C pf = 1006 J / kg K .
qtf = 1,10

800
1006
3600

qtf = 246 W / K
Disposition A-B
Daprs la question 2 :
AB = qtf E B ( Tce 2 T fe ) + E A ( 1 E B ) ( Tce1 T fe )

AB = 246 [0 , 25 ( 65 25 ) + 0 ,38 ( 1 0 , 25 ) ( 80 25 )]
AB = 6315 W

Du bilan sur le fluide froid :

= qtf ( T fs T fe )
on tire :
T fs ( AB ) = T fe +

AB
qtf

= 25 +

6315
246

T fs ( AB ) = 50 ,7 C

Disposition B-A
BA = qtf E A ( Tce 2 T fe ) + E B ( 1 E A ) ( Tce1 T fe )

BA = 246 [0 ,38 ( 65 25 ) + 0 , 25 ( 1 0 ,38 ) ( 80 25 )

BA = 5830 W
On calcule T fs de la mme faon :
T fs ( BA ) = T fe +

BA
qtf

= 25 +

5830
246

T fs ( BA ) = 48 ,7 C

On a ici Tce1 > Tce 2 : la meilleure disposition est bien celle o E1 > E 2 , donc A-B. Le
manque gagner avec B-A dans ce cas particulier est de 485 W, soit 7 ,7 % .
La condition de fonctionnement (question 1) est : Tce 2 > T f 1
T f 1 = T fe + E1 ( Tce1 T fe ) avec E1 = E A
T f 1 = 25 + 0 , 38 ( 80 25 )
T f 1 = 45 , 9 C

La condition est vrifie puisque Tce 2 = 65 C .

DONNES NUMRIQUES

Chaleur de vaporisation LV de diffrents corps la pression de 1 bar

Dsignation

LV

( kJ / kg )

Air
Alcool thylique
Alcool mthylique
Ammoniac
Butane
Eau
Ethane
Ethanol
Fluides frigorignes non chlors
Fuel domestique ( 100 C)
Mthane
Propane
Tolune

197
846
1101
1369
402
2256
490
846
120 240
260
511
448
356

Tension superficielle de leau ( Tsat = 373 K ) : = 38.10 3 N / m

REPRES BIBLIOGRAPHIQUES

Traits gnraux
CRABOL J. Transfert de chaleur (3 vol.), Masson, 1990.
FRAAS A.P. Heat exchanger design. John Wiley, 1989.
INCROPERA F.P., DE WITT D.P. Fundamentals of heat and mass transfer. John Wiley,
1985.
KAKA S., LIU H., PRAMUANJAROENKIJ A. Heat exchangers. Selection, rating and
thermal design. 3rd edition, CRC Press, 2012.
PADET J. Fluides en coulement ; mthodes et modles. 2 dition, Socit Franaise de
Thermique, www.sft.asso.fr. , 2011.
PADET J. Principes des transferts convectifs. 2 dition, Socit Franaise de Thermique,
www.sft.asso.fr. , 2010.
SACADURA J.F., Coord. Initiation aux transferts thermiques. 6 dition, Tec. et Doc.,
Lavoisier, 2000.
TAINE J., PETIT J.P. Transferts thermiques ; mcanique des fluides anisothermes. 2me
dition, Dunod, 2003.
ZUKAUSKAS A. High-performance single-phase heat exchangers. Hemisphere, 1989.

Ouvrages et publications portant sur des thmes plus prcis

Efficacit, mthode NUT, optimisation (chapitre 3)

FEIDT M. Thermodynamique et optimisation nergtique des systmes et procds. Tec. et

Doc., Lavoisier, 1987.
HEGGS P.J. Minimum temperature difference approach concept in heat exchangers
network. Congrs TEC 88 : Progrs rcents dans les changeurs thermiques, Grenoble, 1988.
PIERRE B. Dimensionnement des changeurs de chaleur. Revue Gnrale de Thermique,
260, p.587, 1983.
Calcul des coefficients dchange et des pertes de charge (chapitres 4 6)

BOISSIER A. et al. Les pertes de charge et le transfert thermique, ct gaz, dans les
changeurs de chaleur tubes lisses, circulations orthogonales. Bull. Direction Etudes et
Recherches, EDF, Srie A, N 2/3, 1971.
CHAI H.C. A simple pressure drop correlation equation for low finned tube crossflow heat
exchangers. Int. Comm. Heat Mass Transfer, 15, p. 95-101, 1988.
DE VRIENDT A.B. La transmission de la chaleur. Gatan Morin, 1982 (calcul des
ailettes).
GRETh, Coord. Fouling mechanisms. Ed. Eur. Thermique et Industrie, 1992.
GOSSE J. Guide technique de thermique. Dunod, 1981.
GOTH Y., FEIDT M., LAURO F., BAILLY A. Transferts de chaleur et pertes de charge
en coulement monophasique eau-eau dans les tubes corrugus. Revue Gnrale de
Thermique, N 294-295, 1986 (un mastic a permut les formules de Nu entre cts intrieur et
extrieur).

KAKA S., BERGLES A.E., FERNANDES E.O., Coord. Two-phase flow heat exchangers.
NATO ASI Series, Kluwer Acad. Pub., Dordrecht, 1988.
MAHFOUD M., ABDELAZIZ F., LEBOUCHE M. Evolution du coefficient dchange et
des pertes de charge dans un faisceau de tubes, en fonction de langle dattaque. Int. J. Heat
and Mass Transfer, 30, N12, p.2671, 1987.
MARNER W.J., BERGLES A.E., CHENOWETH J.M. On the presentation of performance
data for enhanced tubes used in shell-and-tubes heat exchangers. J. of Heat Transfer, 105,
p.358, 1983.
MARTINET J. Elments de thermocintique. Tec. et Doc., Lavoisier, 1989 (calcul des
ailettes).
Rseaux (chapitre 7)

BARRERE M. La thermoconomie. Revue Gnrale de Thermique, N 255, p.243, 1983.

BELKEBIR M., et al. Synthse dun rseau dchangeurs de chaleur. Chem. Eng. J., 42,
p.119, 1989.
PIERRE B., dj cit pour le chapitre 3.
Rgimes variables (chapitre 8)

AZILINON D., PIERSON P., PADET J. Constante de temps des changeurs thermiques.
Revue Gnrale de Thermique, 338, p.731, 1991.
BAGUI F., ABDELGHANI-IDRISSI M.A. Rgimes transitoires des changeurs thermiques
tubulaires. Congrs Franais de Thermique SFT 05, Actes vol. 2, p.189, Reims, 2005.
CHITOU N., MAI T.H., PADET J. tude de lefficacit dun changeur en rgime
variable. Entropie, N 220-221, p.87-91, 1999.
EL WAKIL N., REBAY M., PADET J. Numerical study of transient forced convection in
paralle-plate heat exchangers. ICHMT Int. Symposium on Transient Convective Heat Transfer,
Cesme, Turquie, Begell House, 1996.
HADIDI M., GUELLAL M., LACHI M., PADET J. Loi de rponse dun changeur
thermique soumis des chelons de temprature aux deux entres. Int. Comm. in Heat Mass
Transfer, 22, N1, p.145-153, 1995.
HENRION M., FEIDT M. Comportement en rgime transitoire de divers types
dchangeurs de chaleur; modlisation et consquences. Int. Comm. in Heat Mass Transfer,
18, p.731, 1991.
LACHI M., EL WAKIL N., PADET J. The time constant of double pipe and one pass
shell-and-tube heat exchangers in the case of varying flow rates. Int. J. Heat Mass Transfer,
40, N9, p.2067-2079, 1997.
MAI T.H., CHITOU N., PADET J. Heat exchanger effectiveness in unsteady state. Eur.
Phys. J. Applied Physics, 8, p.71-75, 1999.
PIERSON P., AZILINON D., PADET J. Simulation du fonctionnement des changeurs
thermiques soumis des conditions aux limites variables. Revue Phys. Appl. 24, p.93-107,
1989.
PIERSON P., PADET J. Time constant of solar collectors. Solar Energy, 44, N2, p.109114, 1989.

INDEX ALPHABTIQUE DES MATIRES

A
Ailettes, 4.3.4. ; 6.2.3 ; problme 11.
C
Chemine, problme 4.
Chicanes (gomtrie), 4.3.3.1.
Coefficient dchange, 4.1.4 ; 4.3.1.3 ; 4.3.4.3 ; 5.2 ;
5.3 ; 6.2 ; 8.1.
Coefficient de frottement, 4.2.2.3 ; problmes 9, 12.
Condenseurs, 2,4 ; 5.2 ; problmes 5, 10.
Conductance globale, 3.3. ; 6.2.
Constante de temps, ch. 8
Croisements de temprature, 3.7.2 ; 7.6 ; 8.8.3.

- laminaire, 4.2.1 ; 4.4.2 ; 5.2.2 ; 8.5 .

- turbulent, 4.2.2 ; 4.4.3 ; 8.4.2 ; 8.5.
Efficacit
- dun changeur, 3.2 ; 3.4.3 ; 7.2.3 ; 7.3.4.
- dune ailette, 4.3.4
- moyenne en rgime variable, 8.8.
- relative, 3.2 ; 7.2.2 ; 8.8.
Encrassement, 6.1 ; 6.2.3 ; 8.8.
vaporateurs, 2.4 ; 5.3.
F
Facteur de dsquilibre, 3.2 ; 7.2.1 ; 8.8.2.
Facteur de forme, 4.4.1.
Fluide brass, 1.2 ; 3.4.3 ; 3.5.3 ; problme 8.
Formule de Colebrook, 4.2.2.3 ; problme 9.

D
L
Diamtre hydraulique, 4.2.2.6 ; 4.3.1.5 ; 4.3.2 ;
4.4.1 ; 5.2.2 ; problme 11.
DTLM, 2.5.1.

Longueur caractristique, 4.3.1.5 ; 4.3.2 ; 4.5.

Longueur de rfrence, 4.1.5 ; 4.2.2.6 ; 4.3.1.5 ; 4.5.

bullition, 5.3.
changeurs
- chicanes, 1.1 ; 1.3 ; 3.7.2 ; 4.3.3 ; 7.2.4.2 ;
7.4.2.
- contre-courant, 2.3 ; 3.4.2 ; 3.7.1 ; 7.2.5 ;
problmes 2, 3, 6, 13.
- courants croiss, 1.2 ; 3.4.3 ; 3.5.3 ; 4.3.1 ;
problmes 7, 8, 11.
- courants parallles, 1.2 ; 4.3.2 ; problme 13.
- fluide isotherme, 2.4 ; 3.5.2 ; 4.2.1.2 ; 4.2.1.3 ;
problmes 4, 5, 10.
- modules, 1.3.1 ; 7.2.4.2.
- passes, 1.3.2 ; 3.7.1 ; 7.4.2 ;
problmes 8, 9, 10.
- plaques, 1.1.2 ; 1.3.2 ; 6.1 ; 7.4.3 ; problme 7.
- tubes ailets, 4.3.4 ; problme 11.
- bitubes, 4.4 ; 8.2 ; problme 3.
- co-courant, 2.2 ; 3.4.1 ; problme 13.
- compacts, 4.3.4.4 ; problme 11.
- en pingle ( faisceau en U ), 1.1 ; problme 13.
- en srie, 7.2 ; 7.4.1 ; problme 12.
- en srie-parallle, 7.3 ; problme 12.
- quivalents, 7.2.3 ; 7.2.5 ; 7.3.4 ; 8.6.1.
- P-N, 1.3.2 ; 3.4.3 ; 3.7.2 ; 3.7.3 ; 7.2.4.1.
- tubulaires, 1.1.1 ; 4.3.1 ; 4.3.3 ; 7.4.2 ;
problmes 6, 8, 9, 10, 13.
coulement
- autour de faisceaux de tubes, 4.1.4 ; 4.3.1.
- de transition, 4.2.3.

Mthode NUT, 3.6 ; 3.7 ; 7.2 ; problmes.

N
Nombres de Nusselt, Pclet, Prandtl, Reynolds,
Stanton, 4.1.5.
NUT, 3.3 ; 3.4.3 ; 6.2 ; 7.2.5 ; 7.3.4.
P
Pas dun faisceau, 4.3.1.1 ; problmes 9, 13.
Perte de charge, 4.2.2.3 ; 4.3.2 ; problmes 9, 10.
Pincement, 3.7.1 ; 3.7.3 ; 7.6.
Puissance dun changeur, 3.1 ; 6.2.
R
Rseau, problmes 12, 13, 14.
- fluide intermdiaire, 7.5.1.
- maill, 7.4.
- ouvert, 7.5.2.
Rugosit, 4.2.2.3.
T
Temprature
- de film, 4.1.3 ; 4.1.4 ; 5.1 ; problme 10.
- de mlange, 2.1 ; 2.2. ; 4.1.2 ; 4.1.4.
- de rfrence, 4.1 ; problmes 8, 9.
Temps de retard, 8.3.1 ; 8.6.3.

Thermodpendance, 4.2.1.1 ; 4.2.2.1.

Tubes
- corrugus, 4.2.2.5 ; 4.4.4.
- faisceaux, 4.3.1.1 ; 4.3.1.6 ; 4.3.3 ; 4.3.4.3 ;
problmes 9, 13.
- nappes, 4.3.1.1 ; problme 9.
- section annulaire, 4.4 ; problme 3.
- section circulaire, 4.2.1.2 ; 4.2.2.2.
- section rectangulaire, 4.2.1.3 ; problme 9.

V
Viscosit, 4.2.1.1.
Vitesse massique, 4.1.5 ; 4.3.1.2.
Vitesse de rfrence, 4.3.1.2 ; 4.3.4.3.