univdocs.

com
Convection - Échangeurs de chaleur

Exercice 1

Introduction aux critères adimensionnels. Formules liant Nusselt, Reynolds et Prandtl.

Exercice 2

De l'eau à 25°C, alimente un tube de 2,5 cm de diamètre, permettant de refroidir un réacteur nucléaire, avec un débit de 100 l/mn.
Déterminer le débit de chaleur transféré et la température de l'eau à la sortie pour un tube de 4,5 m de long, la température de la
surface interne du tube étant constante et égale à 150°C.

Viscosité μ = 10-3 Pa.s

Capacité calorifique C P = 4,18 kJ/kg.°C
Masse volumique ρ = 1000 kg/m3
Conductibilité thermique λ = 0,64 W/m.°C

Exercice 3

De l'eau à 20°C alimente avec un débit volumique total de 800 l/mn un échangeur constitué d'un faisceau tubulaire de 20 tubes
cylindriques en parallèle de 2 cm de diamètre intérieur.
A l'extérieur de ce faisceau se condense de la vapeur d'eau saturée à 100°C, ce qui permet de maintenir la température des parois
extérieures du faisceau à 100°C.
En admettant que leur épaisseur est négligeable, on demande de préciser la longueur de chacun des tubes du faisceau permettant
d'amener la température de l'eau de 20°C à 80°C.
Les propriétés physico-chimiques de l'eau seront supposées indépendantes de la température :

Viscosité μ = 8.10-4 Pa.s

Capacité calorifique C P = 4,18 kJ/kg.°C
Masse volumique ρ = 1000 kg/m3
Conductibilité thermique λ = 0,64 W/m.°C

Exercice 4

De l'eau à 10°C est préchauffée dans un serpentin cylindrique de 10 m de long et de 1 cm de diamètre intérieur. La température de
la paroi extérieure est maintenue constante et égale à 100°C.
1. En supposant que le coefficient de transfert global entre la paroi exté​rieure du serpentin cylindrique et le fluide intérieur est
constant tout au long de l'échangeur, quelle est la température de l'eau à la sortie du serpentin, le débit massique d'eau à
10°C étant égal à 3 kg/mn.
2. A la sortie du serpentin, l'eau chaude alimente un circuit de chauffage par l'intermédiaire d'un tube cylindrique de 4 cm de
diamètre recouvert par un manchon cylindrique d'isolant de 4 cm d'épaisseur et de conductibilité thermique λ = 0,07 W/m.
°C.
La température du milieu extérieur est constante et égale à 20°C et le coefficient de transfert convectif de chaleur avec le
milieu extérieur étant constant et égal à 11,6 W/m2 .°C, on demande de préciser la température de l'eau chaude à
l'extrémité de 25 m de conduite thermiquement isolée. Quel est le pourcentage d'énergie perdue au cours du transport par
rapport au débit de chaleur fourni dans le serpentin.
Les propriétés physico-chimiques de l'eau seront supposées indépendantes de la température :

ρ = 1000 kg/m3 μ = 10-3 Pa.s

C P = 4,18 kJ/kg.°C λ = 0,64 W/m.°C

Exercice 5

Un réacteur parfaitement agité est le siège d'une réaction exothermique. Il est alimenté par un débit massique de 45 t/h d'un

converted by W eb2PDFConvert.com
univdocs.com
mélange réactionnel à 15°C. On désire maintenir le fluide contenu dans le réacteur à une température uniforme de 45°C avec un
taux de conversion du réactif limite tel que le débit de chaleur produit par la réaction est de 106 W.
Le réacteur étant bien isolé extérieurement, on pense réaliser l'isothermicité en incorporant un serpentin interne de
refroidissement cylindrique de 2 cm de diamètre, alimenté par une saumure à -5°C, avec un débit massique de 5 t/h.
Déterminer la longueur du serpentin ainsi que la température de la saumure à la sortie du serpentin. Pour cela :
1. Écrire le bilan de chaleur dans un élément de longueur dx du serpentin, en supposant que la seule résistance au transfert de
chaleur est celle existant à l'intérieur du tube cylindrique constituant le serpentin.
2. En supposant constant le coefficient de transfert entre le serpentin et le fluide réactionnel, écrire le bilan global de chaleur
dans le réacteur.
3. Calculer le coefficient de transfert de chaleur .
4. Calculer la longueur du serpentin et la température de la saumure à la sortie du serpentin.

ρ Cp λ μ
(kg/m³) (kJ/kg°C) (W/m°C) mPa.s
Saumure 1230 4,1 0,780 1,15
Fluide réactionnel 980 2,26

Exercice 6

De l'eau chauffée à 80°C alimente, avec un débit massique de 360 kg/h, un échangeur à double tube cylindrique concentrique afin
de préchauffer un gaz de -5°C à +55°C.
Le gaz circule dans le tube intérieur avec un débit massique de 180 kg/h. Le tube intérieur, de diamètre 0,175 m, a une épaisseur
négligeable et le coefficient de transfert de chaleur convectif côté liquide est supposé constant et égal à 5225 W/m2 .°C. Le
coefficient de transfert convectif côté gaz sera évalué à l'aide de la relation de Petukhov-Gnielinski.
L'échangeur étant globalement adiabatique et le coefficient de transfert global étant supposé constant, on demande de préciser la
surface d'échange permettant de répondre à l'objectif proposé dans les deux cas suivants :
l'écoulement du gaz et de l'eau se faisant à co-courant;
l'écoulement du gaz et de l'eau se faisant à contre-courant.
Quelle est dans les deux dispositions la température de l'eau à la sortie de l'échangeur ?
(C P de l'eau : 4190 J/kg.K)

Les propriétés physico-chimiques du gaz, supposées indépendantes de la température, seront prises égales à:

ρ = 0,676 kg/m3 C P = 2238 J/kg.K λ = 32,81 mW/m.K μ = 10,28.10-6 Pa.s

On rappelle que dans un échangeur à double tube cylindrique concentrique, le débit de chaleur transféré entre le fluide chaud et
le fluide froid est égal à :
Φ = U.S.Δϴ
où
U est le coefficient de transfert de chaleur global;
S est la surface d'échange;
Δϴ est la moyenne logarithmique des températures.

Exercice 7

On veut refroidir 3000 kg/h de nitrobenzène de 80 °C à 30 °C en utilisant 1500 kg/h d'eau, la température initiale de l'eau est de
20 °C. On utilise un échangeur constitué par deux tubes concentriques (tube intérieur, diamètres 25/33; tube extérieur, diamètres
50/60). Le nitrobenzène passe dans le tube intérieur et l'eau dans l'espace annulaire.
a. Calculer la température de sortie de l'eau.
b. Peut-on opérer à contre- courant et à courants parallèles ?
c. Sachant que le coefficient h1 (nitrobenzène/paroi) est égal à 1410 W/m2 °C, que le coefficient d'échange h2 (paroi/eau) est
égal à 1935 W/m2 °C, calculer le coefficient d'échange global Uext rapporté au diamètre extérieur du tube intérieur.
d. En déduire la longueur de l'échangeur Le.
Données :

ρ Cp λ μ
(kg/m³) (kJ/kg°C) (W/m°C) mPa.s
eau 1000 4,18 0,630 0,7

converted by W eb2PDFConvert.com
univdocs.com
nitrobenzène 1200 1,38 0,159 1

conductivité thermique de la paroi λ = 46,4 W/m°C

Exercice 8

De l'acide sulfurique circule avec un débit de 4500 kg/h dans un circuit qui comprend deux réservoirs en série où il est en contact,
par agitation avec des serpentins de refroidissement. De l'eau circule dans les serpentins à contre-courant de l'acide. Sachant que
le premier réservoir traversé par l'acide a un coefficient global de transmission U1 = 1160 W/m².°C et le second un coefficient
U2 = 731 W/m².°C, calculer la surface totale de serpentins nécessaire pour le refroidissement. On négligera les déperditions
calorifiques dans le circuit. On commencera par calculer le débit d'eau de refroidissement et la température de l'eau entre 1 et 2.
Acide sulfurique : Cp = 1,50 kJ/kg.°C Eau : Cp = 4,18 kJ/kg°C
Les températures aux différents points des circuits sont indiquées sur le schéma ci-dessous :

Exercice 9

Un échangeur coaxial de 36 m de long est constitué par un tube cylindrique de 4,2 cm de diamètre extérieur à l'intérieur duquel
est situé un tube cylindrique coaxial de diamètre intérieur 2,5 cm. L'épaisseur du tube intérieur et du tube extérieur est de 2 mm
De l'eau chaude à 85 °C est introduite dans le tube intérieur avec un débit de 60 l/mn. On se propose de refroidir l'eau avec une
solution aqueuse à 5 °C alimentée dans l'espace annulaire avec un débit de 30 l/mn.
1. L'alimentation des deux fluides pouvant se faire à co-courant ou à contre-courant on demande d'évaluer théoriquement
dans les deux cas le débit transféré du fluide chaud au fluide froid avec les hypothèses suivantes :
le tube extérieur est adiabatique;
le coefficient de transfert global entre le fluide chaud et le fluide froid est constant tout au long de l'échangeur.
2. Les coefficients de transfert locaux supposés constants étant évalués dans tous les cas par la relation de Sieder et Tate. On
demande de préciser quantitativement pour les deux associations possibles :
les températures des deux fluides à la sortie de l'échangeur
le débit de chaleur transféré du fluide chaud au fluide froid
Les propriétés physico-chimiques des deux fluides seront supposées égales à celles de l'eau à une température moyenne de
45 °C, soit :
ρ = 990 kg/m³ Cp = 4,18 kJ/kg.°C λ = 0,602 W/m.°C μ = 0,65 mPa.s
La conductivité thermique du métal constituant les tubes est égale à λ = 45,2 W/m.°C
Quelle est l'association qui conduit au débit de chaleur transféré maximal ?

Exercice 10

On veut refroidir 4,5 t/h de benzène d'une température de 80°C à 30°C à l'aide d'un échangeur tubulaire simple à contre-courant.
L'eau de refroidissement circule en double enveloppe et entre dans l'échangeur à une température de 20°C.
Le tube interne, dans lequel circule le benzène, a un diamètre intérieur de 22 mm et un diamètre extérieur de 25 mm et est en
acier (λP = 45 W/m °C).
La double enveloppe est parfaitement isolée de l'extérieur.
Les coefficients d'échange convectif eau/paroi et benzène/paroi sont respectivement :
he = 850 W/m2 °C et hb = 1700 W/m2 °C

1. Quelle est la longueur totale Lt du tube nécessaire si on cherche à avoir le plus faible débit d'eau avec une température de
sortie de l'eau ne dépassant pas 50°C ?

converted by W eb2PDFConvert.com
univdocs.com
Cpbenzène = 1,900 kJ/kg °C et Cpeau = 4,185 kJ/kg °C
2. Pour des questions d'encombrement on désire remplacer l'échangeur tubulaire simple par un système à deux calandres -
quatre passes.
Calculer la nouvelle surface S2 nécessaire à l'échangeur (rapportée au diamètre extérieur du tube intérieur). En déduire la
longueur d'une passe.

Exercice 11

Un réservoir de vapeur cylindrique, horizontal, est rempli d'une masse M = 19,3 kg d'eau sous forme de vapeur saturante.
Ce réservoir est constitué d'un cylindre métallique de résistance thermique négligeable, de longueur L = 5 m et de diamètre
D = 50 cm.
Ce réservoir est parfaitement isolé sur les faces correspondantes aux sections droites, mais perd de la chaleur par sa face latérale
par convection naturelle et par rayonnement.
Ce réservoir est placé dans l'air ambiant supposé à la température constante de Ta = 15°C.
1. Calculer la température de surface du réservoir . On admettra que cette température est pratiquement égale à celle de la
vapeur saturante Tvs .
On pourra également utiliser avec une bonne approximation la relation de Bertin reliant la pression de vapeur saturante de
l'eau Pvs à son volume spécifique Vvs sous la forme :
Pvs Vvs = 2,038 (avec Pvs en bars et Vvs en m³/kg)
ainsi que la relation de Duperay donnant, pour l'eau, la température de la vapeur saturante Tvs en fonction de la pression
Pvs :
(avec Pvs en bars et Tvs en °C)
2. On désire maintenir la vapeur d'eau stockée dans l'état ci-dessus, malgré les pertes thermiques du réservoir. Pour cela, on
dispose une résistance électrique chauffante dans le réservoir. On va chercher à déterminer la puissance de chauffe de
cette résistance, par effet Joule.
a. Calculer le flux de déperdition Φ CV par convection libre du réservoir. On négligera les pertes aux extrémités du
cylindre.
b. Calculer le flux de déperdition du réservoir par rayonnement Φ R. On considérera la face latérale du cylindre
assimilable à un corps noir rayonnant vers l'extérieur sous 4π stéradians, en négligeant les faces extrêmes, les
surfaces réceptrices étant supposées également noires, de température constante et égale à la température de l'air
ambiant Ta = 15°C. On rappelle la valeur de la constante de Stefan σ = 5,67 10-8 W/m2 K4
c. En déduire la puissance de chauffe électrique P (en kW) nécessaire pour maintenir l'eau contenue dans le réservoir
sous forme de vapeur saturante.

converted by W eb2PDFConvert.com
univdocs.com
On rappelle qu'en convection libre et pour un cylindre horizontal (échelle caractéristique L = D diamètre) :

régime laminaire 104 < Ra < 109 C = 0,53 n = 1/4

régime turbulent 109 < Ra < 1012 C = 0,13 n = 1/3

Les propriétés physiques de l'air pourront être évaluées par :

T (K) ν (m2 /s) λ (W/m.K Pr
300 15,68 · 10 -6 0,02624 0,708
400 25,90 · 10-6 0,03365 0,689
500 37,90 · 10 -6 0,04038 0,680
600 51,34 · 10 -6 0,04659 0,680
700 68,20 · 10 -6 0,05132 0,680

Exercice 12

Un échangeur coaxial, parfaitement isolé extérieurement, est utilisé pour réchauffer un fluide de chaleur spécifique CP = 3,8
kJ/kg°C et de masse volumique ρ = 950 kg/m3 .
Ce fluide circule à l'intérieur du tube, tandis que de la vapeur saturante, à une température TVS = 130°C, circule dans la double
enveloppe, cédant sa chaleur par condensation externe aux tubes. On suppose que la vapeur est toujours saturante en sortie de
l'échangeur.
Le fluide entre dans l'échangeur à la température TE = 20°C et à un débit volumique Q V1 = 0,7 m3 /h. Sa température de sortie est
TS1 = 100°C.
On cherche à déterminer, dans ces conditions, le coefficient d'échange global U1 , ainsi que la surface totale d'échange A. Pour
cela, on augmente le débit du fluide jusqu'à une valeur Q V2 = 1,2 m3 /h, correspondant à une conductance globale U2 . La
température de sortie du fluide devient, dans ces conditions TS2 = 95°C, la température d'entrée restant fixée à TE = 20°C.
Les régimes d'écoulement dans les tubes étant supposés turbulents, on pourra supposer, a priori, une dépendance de la
conductance interne fluide/paroi hi , du type : hi = KV·Q V0,8 , où KV est une constante à déterminer.
La conductance externe en condensation, sera supposée rester constante et égale à he = 3 000 W/m2 °C. On pourra négliger la
résistance conductive des tubes.
1. Calculer les puissances Φ1 et Φ2 échangées (en W), lors de chacun des essais.
2. Calculer les moyennes logarithmiques des températures lors des deux essais, soient Δϴ L1 et Δϴ L2 .
3. En déduire le rapport U1 /U2 .
4. Déterminer la valeur de la constante KV (unités SI).
5. Calculer la conductance U1 , conductance globale de transfert correspondant au débit Q V1 (unités SI).
6. En déduire la surface totale d'échange A (m2 ) nécessaire.
Reprendre le problème en utilisant la méthode NUT - efficacité.

Rayonnement

Exercice 1

Présentation de la notion de facteur d'angle (facteur géométrique, facteur de forme, facteur de vue).
Cas d'un corps convexe dans un corps concave.

Exercice 2

Un four tubulaire vertical semi-ouvert de forme cylindrique, d'un diamètre D = 75 cm et de hauteur L = 150 cm, est ouvert à son
extrémité supérieure, vers l'extérieur considéré comme un corps noir à la température T∞ constante de 20°C .
Les parois intérieures du four, chauffées électriquement sont considérées comme des corps noirs maintenus par régulation aux

converted by W eb2PDFConvert.com
univdocs.com TF06_convection_01.xmcd

TF06 ‐Convecon‐Exercice1
Présentation des critères adimensionnels utilisés couramment, ainsi que les paramètres et leurs unités :
um × D ρ × um × D μ
ν= v
isc
osi
téc
iné
ma q
ue(
m²/
s)
Reynolds : Re = = ρ
ν μ
μ × CP λ
ν α= diffusivité thermique (m²/s)
Prandtl : Pr = = ρ × CP
α λ

hc × D ρ masse volumique (kg/m³)

Nusselt : Nu =
λ μ viscosité (Pa.s)

λ c
ond
ucv
itét
her
miq
ue(
W/m.
K)
hc coefficient d'échange thermique
conv
ecfo uc onduct
an ce(W/ m².K
) CP capacité calorifique (J/kg.K)

D diamètre de la conduite um vitesse moyenne d'écoulement (m/s)

ou vitesse débitante en fût vide

Remarque : en ce qui concerne la viscosité µ, de dimension M.L-1.T-1, la seule unité officielle du S.I. est le
Pas
cal.
seconde( Pa.s
).L ePoiseuilleaét épr opos éma i
sa bandonn éc a rn onu l isé.
Lemil
li
‐pasc al
.sec ond e(mP a.
s)e sté q ui
va lentàl ac enp oi
se(cP ).
N kg × m s kg Références : h p:
//www.bi
pm.fr
/et h p:/
/phy
s i
cs.
nis
t.
gov
/
Pa × s = ×s = × =
2 2 2 m× s BI
PM =Bur eauI nternaonaldesPoidsetMesures
m s m NIST=Na on alI
nst ut
eofStandar
dsa ndTechnol
ogy

Relations entre les différents critères adimensionnels :

Calc
ulpréalabledunombredeRey nol
ds.SiRe<2300l'éc
oulementestlaminair
e.ÀparedeRe>10000l 'éc
oul
ement
estt
urbul
ent .Ent
re2300et10000onpar l
edez onecr
iq ue(àé v
iters
ipossi
ble)
.
Àparrdesnombr esdeRe
ynoldsetPrandtl
,ondéter
mine,plusoumoinsfaci
lement,etav
ecplusoumoi nsde
préc
isi
on,lenombr edeNussel
t.C'
estl
eNus sel
tquinousdonneraalor
slecoeffic i
entdeconveconh c .

1. Écoulement laminaire Nu = 3.66 lorsque la température de paroi est constante pour une conduite
longue : le régime
Nu = 4.36 lorsque le débit calorifique à travers la paroi laminaire est établi
est constant

0.8 0.4 TP > Tm 4

2. Écoulement turbulent Nu = 0.023 × Re × Pr Re > 10
Dius‐Boelt
er 0.8 0.3 TP < Tm 0.6 < Pr < 160
Nu = 0.023 × Re × Pr
(1930)
0.14
4¸5 1¸3 æ μm ö Re > 10
4
Sieder-Tate Nu = 0.027 × Re × Pr ×ç ÷
(1936) è μP ø 0.7 < Pr < 16700

-2
Petukhov f = ( 0.79 × ln ( Re) - 1.64)
6
(1970)
( Re- 1000) × Pr 3000 < Re < 5 × 10
f
( Re - 1000) × Pr ×
8 × ( 0.79 × ln ( Re ) - 1.64)
Gnielinski 2
8
(1976) Nu = = 0.5 < Pr < 2000
( ) 12.7 × ( Pr - 1)
2¸3 f 2¸3
1 + 12.7 × Pr -1 × 1+
8
( 0.79 × ln ( Re) - 1.64) × 8

MH 20/05/2013
univdocs.com TF06_convection_02.xmcd

TF06 - Convection - Exercice 2 Données :

-3 -3 L
ρ := 1000× kg× m D := 2.5× cm μ := 10 × Pa× s V' := 100× débit volumique
min
W kJ
λ := 0.64 CP := 4.18 V' m
m× K kg× °C um := = 3.395
2 s
D
TE := 25°C Lt := 4.5× m TP := 150°C π×
4

1°) Bilan dans un élément d'épaisseur dx du fluide dans le tube :

hc× π× D
( )
dT
ρ× V' × CP× dT = hc× π× D× dx× TP - T = × dx
TP - T ρ× V' × CP
T Lt
ó S ó hc× π× D æ TP - TE ö hc× π× D
1
On intègre de l'entrée à la ô dT = ô dx ln ç ÷= × Lt
TP - T ô ρ× V' × CP
sortie, sur une longueur Lt ô
õ0 è TP - TS ø ρ× V'× CP
õT
E

2°) Calcul du coefficient de convection : Définition des critères adimentionnels :

ρ × u× D 4× ρ× V' 4× 1000× 0.1

Reynolds : Re := Re := = 84883 = 84883
μ π× D × μ -3
π× 0.025× 10 × 60

Prandtl : μ× C P 0.001× 4180

Pr := = 6.531 = 6.531 Régime turbulent
λ 0.64
hc× D
Nusselt : Nu =
λ
Formule de Petukhov
-2 -2 f -3 f -2
f := ( 0.790× ln ( Re) - 1.64) = 1.863 ´ 10 = 2.329 ´ 10 = 4.826 ´ 10
8 8
Formule de Gnielinski
f
× ( Re - 1000) × Pr -3
8 2.329× 10 × ( 84883 - 1000) × 6.531
Nu := = 504.6 = 504.6
1 + 12.7×
f
× Pr
2¸3
-1( ) 1 + 12.7× 4.826× 10
-2 (
× 6.531
2¸3
-1 )
8
Formule de Petukhov - Gnielinski
( Re- 1000) × Pr
8× ( 0.790× ln ( Re) - 1.64)
2
Nu× λ W
Nu := = 504.6 hc := = 12919 ×

1+
12.7× Pr( 2¸3
-1 ) D 2
m × °C
( 0.790× ln ( Re) - 1.64) × 8
æ hc× π× D ö
-ç × Lt÷
ρ× V' × CP
TS := TP - TP - TE × e è
( ) ø = 85.1× °C
3°) Calcul de la température de sortie :

150 - ( 150 - 25) × exp éê-æç

12919 × π× 0.025× 4.5× 60 öù
÷ú = 85.1
4°) Calcul du débit de chaleur : ë è 1000× 0.1× 4180 øû

(
Φ := ρ× V' × CP× TS - TE = 419× kW ) 1000×
0.1
60
× 4180× ( 85.1 - 25) = 419 ´ 10
3

MH 1/2 11/05/2012
univdocs.com TF06_convection_02.xmcd

Re = 84883 Pr = 6.531

4 1
Formule de Sieder & Tate
5 3
simplifiée Nu := 0.027× Re × Pr = 442.67

Nu× λ W 442.67× 0.64

hc := = 11332 × = 11332
D 2 0.025
m × °C

æ hc× π× D ö
3°) Calcul de la température de sortie : -ç × Lt÷
ρ× V' × CP
TS := TP - TP - TE × e è
( ) ø = 79.7× °C

150 - ( 150 - 25) × exp éê-æç

11332 × π× 0.025× 4.5× 60 öù
÷ú = 79.7
ë è 1000× 0.1× 4180 øû
4°) Calcul du débit de chaleur :

( )
Φ := ρ× V' × CP× TS - TE = 381× kW 1000×
0.1
60
× 4180× ( 79.7 - 25) = 381 ´ 10
3

MH 2/2 11/05/2012
univdocs.com TF06_convection_03.xmcd

TF06 - Convection - Exercice 3

3
-3 -4 800 L -4 m
ρ := 1000× kg× m D := 2× cm μ := 8× 10 × Pa× s V' := × = 6.667 ´ 10
20 min s
W kJ débit volumique
λ := 0.64 CP := 4.18
m× K kg× °C u :=
V'
= 2.122
m
2 s
D
TE := 20°C TS := 80°C TP := 100°C π×
4

1°) Bilan dans un élément d'épaisseur dx du fluide dans un tube :

hc× π× D
( )
dT
ρ× V' × CP× dT = hc× π× D× dx× TP - T = × dx
TP - T ρ× V' × CP Lt
T
ó S ó hc× π× D
1
ô dT = ô dx
On intègre de l'entrée à la sortie, sur une longueur L t ô TP - T ô ρ× V' × CP
õT õ0
E
æ TP - TE ö hc× π× D
ln ç ÷= ×L
TP - TS ρ× V' × CP t
è ø
2°) Calcul du coefficient de convection :
Définition des critères adimentionnels :
-3
ρ× u× D 4× ρ× V' 4× 1000× 800× 10
Reynolds : Re := Re := = 53052 = 53052
μ π× D × μ -4
20× 60× π× 0.02× 8× 10
μ× C P -4
8× 10 × 4180
Prandtl : Pr := = 5.225 = 5.225 Régime turbulent
λ 0.64
hc× D
Nusselt : Nu = Formule de Petukhov
λ
-2 -2 f -3 f -2
f := ( 0.790× ln ( Re) - 1.64) = 2.068 ´ 10 = 2.585 ´ 10 = 5.084 ´ 10
8 8
Formule de Gnielinski

f
× ( Re - 1000) × Pr -3
8 2.585× 10 × ( 53052 - 1000) × 5.225
Nu := = 305.8 = 305.9
( ) × ( 5.225 - 1)
2¸3 -2 2¸3
1 + 12.7×
f
× Pr -1 1 + 12.7× 5.084× 10
8
Formule de Petukhov - Gnielinski
Nu× λ W
( Re- 1000) × Pr hc := = 9786×
D 2
m × °C
8× ( 0.790× ln ( Re) - 1.64)
2
Nu := = 305.8
12.7× ( Pr - 1)
2¸3
305.8× 0.64 3
1+ = 9.786 ´ 10
( 0.790× ln ( Re) - 1.64) × 8 0.02

3°) Calcul de la longueur d'un tube :

ρ× V' × CP æ TP - TE ö -3
1000× 800× 10 × 4180 æ 100 - 20 ö
Lt := × ln ç ÷ = 6.28 m × ln ç ÷ = 6.283
hc× π× D
è TP - TS
ø 20× 60× 9786× π× 0.02 è 100 - 80 ø

MH 1/2 11/05/2012
univdocs.com TF06_convection_03.xmcd

Re = 53052 Pr = 5.225

1
Formule de Sieder & Tate 0.8 3 0.8 1¸3
simplifiée Nu := 0.027× Re × Pr = 282.2 0.027× 53052 × 5.225 = 282.2

Nu× λ W 282.153× 0.64 3

hc := = 9029× = 9.029 ´ 10
D 2 0.02
m × °C

3°) Calcul de la longueur d'un tube :

ρ× V' × CP æ TP - TE ö -3
1000× 800× 10 × 4180 æ 100 - 20 ö
Lt := × ln ç ÷ = 6.81 m × ln ç ÷ = 6.81
hc× π× D
è TP - TS
ø 20× 60× 9029× π× 0.02 è 100 - 80 ø

MH 2/2 11/05/2012
univdocs.com TF06_convection_04.xmcd

TF06 - Convection - Exercice 4

1°) serpentin
kg
ρ := 1000 × D := 1cm Ls := 10m μ := 1 × mPa × s Tp := 100°C
3
m

W kJ kg kg
λ := 0.64 CP := 4.18 m'e := 3 × = 0.05 T0 := 10°C
m × °C kg × °C min s

Définition des critères adimentionnels :

m'e
ρ × um × D ρ D 4 × m'e
Reynolds : Re = = ρ× × Re := = 6366 Re = 6366
μ 2 μ π×D×μ
D
π× Régime turbulent
4 4´ 3
= 6366
-3
60 ´ π ´ 0.01 ´ 10
μ × CP -3
10 ´ 4180
Prandtl : Pr := = 6.531 = 6.531 Pr = 6.531
λ 0.64

Nusselt :
-2 -2 f -3 f -2
Formule de Petukhov f := ( 0.790 × ln ( Re) - 1.64) = 3.588 ´ 10 = 4.485 ´ 10 = 6.697 ´ 10
8 8
f
Formule de Gnielinski × ( Re - 1000) × Pr
8
Nu := = 50.4
1 + 12.7 ×
f
× Pr (
2¸3
-1 )
8

hc × D 50.04 ´ 0.64 3 Nu × λ 3 W
Nu = = 3.203 ´ 10 hc := = 3.223 ´ 10 ×
λ 0.01 D 2
m × °C

Bilan thermique :

dT hc × π × D
(
m'e × CP × dT = hc × π × D × dx × Tp - T ) Tp - T
=
m'e × CP
× dx

æ hc × π × D ö÷
æ Tp - T0 ö hc × π × D -ç ×L
ln ç ç m'e × CP s ÷
÷= ×L
TS := Tp - ( Tp - T0) × e è ø = 99.3 × °C
Tp - TS m'e × CP
è ø

æ 3223 × π × .01 × 10ö

-ç
3
÷
ç × 4180 ÷
100 - ( 100 - 10) × e è ø = 99.29
60

( )
Φ0 := m'e × CP × TS - T0 = 18.662 × kW TS = 99.3 × °C

MH 1/2 11/05/2012
univdocs.com TF06_convection_04.xmcd

2°) tube 4cm sur 25 m Définition des paramètres :

W W
Di := 4 × cm ep := 4 × cm De := Di + 2 × ep he := 11.6 × λP := 0.07 ×
2
m × °C m × °C
De = 12 × cm

4 × m'e 4´ 3
Retube := = 1592 = 1592
π × Di × μ 60 ´ π ´ 0.04 ´ 10
-3

Régime laminaire
3.66 ´ 0.64 Nu' × λ W
Nu' := 3.66 = 58.56 hi := = 58.56 ×
0.04 Di 2
m × °C
Le := 0.03 × Retube × Di = 1.91 m

Bilan thermique : dΦ = -m'e × CP × dT

T - Tpi Tpi - Tpe Tpe - Tex T - Tex

dΦ = = = =
1 1 æ De ö 1 1 1 æ De ö 1
hi × π × Di × dx × ln ç ÷ he × π × De × dx + × ln ç ÷+
2 × π × λP × dx è Di ø hi × π × Di × dx 2 × π × λP × dx è Di ø he × π × De × dx

En prenant comme surface d'échange de référence la surface du tube interne (pDidx) :

(
π × Di × dx × T - Tex )
dΦ =
Di æ De ö Di
donne (
-m'e × CP × dT = U × π × Di × dx × T - Tex )
1
+ × ln ç ÷+
hi 2 × λ' è Di ø he × De dT æ U × π × Di ö
= -ç ÷ × dx
T - Tex
avec
è m'e × CP ø

1 1 W
U := = 2.78 U = 2.78 ×
1 Di Di æ De ö 1
+
.04
+
.04 æ 12 ö
× ln ç ÷
2
m × °C
+ + × ln ç ÷ 58.56 .12 ´ 11.6 2 ´ 0.07 è 4 ø
hi he × De 2λP
è Di ø

Lt := 25 × m la température d'entrée dans le tube T'E := TS T'E = 99.3 × °C Tex := 20°C

est la température de sortie du
serpentin

æ U × π × Di ö÷
-ç ×L
æ TS - Tex ö æ U × π × Di ö çè m'e × CP t÷ø
ln ç ÷ = -ç ÷ × Lt (
T'S := Tex + T'E - Tex × e ) = 96.0 × °C
è TE - Tex ø è m'e × CP ø
æ 2.78 × π × 0.04 × 25ö
-ç
3
÷
ç × 4180 ÷
20 + ( 99.29 - 20) × e è ø = 96.04
60

( )
Φ'0 := m'e × CP × T'S - T'E = -678.254 W

Φ'0
=
(
m'e × CP × T'S - T'E
=
T'S - T'E ) T'S - T'E
= -0.036
Φ'0
= -3.6 × %
Φ0 (
m'e × CP × TS - T0 TS - T0 ) TS - T0 Φ0

MH 2/2 11/05/2012
univdocs.com TF06_convection_05.xmcd

TF06 - Convection - Exercice 5

serpentin (saumure) ms

kg tonne ρs -1
ρs := 1.23× D := 2cm ms := 5× um := um = 3.594 m× s
L h 2
D
W kJ π×
λ := 0.78 CP := 4.1 débit massique 4 μ := 1.15× mPa× s
m× K kg× °C

kJ tonne
milieu réactionnel : CPM := 2.26× mM := 45× TsE := -5× °C
kg× °C h
6
Q := 10 × W T0 := 15× °C TM := 45× °C

1°) Bilan serpentin :

hc× π× D æ TM - TsE ö hc× π× D
( )
dT
ms× CP× dT = hc× π× D× dx× TM - T = × dx ln ç ÷= ×L
TM - T ms× CP TM - TsS ms× CP s
è ø

2°) Bilan réacteur :

( )
Q = mM× CPM× TM - T0 + ms× CP× TsS - TsE ( )
6 45000

TsS := TsE +
Q - mM× CPM× TM - T0 (
= 21.8 °C
) -5 +
10 -
3600
× 2260× ( 45 - 15)
= 21.78
ms× CP 5000
4100
3600
Définition des critères adimentionnels :

ρs× um× D 4× ms 4 ´ 5000

Reynolds : Re := Re := = 76886 Re = 76886
μ π× D × μ -3
3600× π× 0.02 ´ 1.15× 10
μ× C P -3
1.15× 10 × 4100
Prandtl : Pr := = 6.045 = 6.045 Régime turbulent
λ 0.78

hc× D
Nusselt : Nu = ( Re- 1000) × Pr
λ
8× ( 0.79× ln ( Re) - 1.64)
2
Nu := = 448.1
( )
Nusselt : Formule de Petukhov-Gnielinski
2¸3
( 76886- 1000) × 6.045
12.7× Pr -1
1+
( 0.79× ln ( Re) - 1.64) × 8
8× ( .79× ln ( 76886) - 1.64)
2
= 448.1

1+
(
12.7× Pr
2¸3
-1 )
( .79× ln ( 76886 ) - 1.64) × 8

Nu× λ kW 448.1 ´ 0.78 3

hc := = 17.47 = 17.48 ´ 10
D 2 0.02
m × °C

ms× CP æ TM - TsE ö 5000

4100
Ls := × ln ç ÷ = 4.0 m × ln æç
3600 45 + 5 ö
hc× π× D
èTM - TsS
ø ÷ = 3.98
17470 × π× 0.02 è 45 - 21.78 ø

MH 1/1 15/05/2012
univdocs.com TF06_convection_06.xmcd

TF06 ‐Convecon‐Exerci
ce6 Nouvelles valeurs pour le gaz (conformes au poly 2013)

Définiondespar
amèt
respourl
egaz(
flui
def
roi
dqui
ser
échauffe):

kg -6 mW J
ρ := 0.676 × μ := 10.28 × 10 × Pa × s λ := 32.81 CP := 2238 D := 0.175 × m
3 m× K kg × K
m
-1 kg ṁg
ṁ g := 180 × kg × h = 0.05
s ρ m
T GE := -5°C T GS := 55°C um := = 3.075
débit massique 2 s
D
π×
4
Paramètres de l'eau (fluide chaud qui se refroidit) :
kg kg J
ṁ e := 360 × = 0.1 CPe := 4190 × T eE := 80°C
h s kg × K

-1
Posons : (chaud) C1 := ṁ e × CPe = 419.0 × W × K T 1E = T eE = 80 × °C T 1S := T eS

-1
(froid) C2 := ṁ g × CP = 111.9 × W × K T 2E = T GE = -5 × °C T 2S = T GS = 55 × °C

( ) (
Φ = C1 × T 1E - T 1S = C2 × T 2S - T 2E ) ( )
Φ := C2 × T 2S - T 2E = 6.714 × kW

T
empé
rat
uredes
oredel
'
eau: C2
T 1S := T 1E -
C1
(
× T 2S - T 2E = 64.0 × °C )
Définiondesc
rit
èresadi
men onnel
spourl
ecal
cul
dehc
ôtégaz:

ρ × um × D 4 ×ṁg
Reynolds : Re := = 35387 Re := Re = 35387
μ π× μ× D
Régime turbulent
μ × CP
Prandtl : Pr := = 0.701
λ

hc × D
( Re- 1000) × Pr
Nusselt : Nu = 2
λ 8 × ( 0.79 × ln ( Re ) - 1.64)
Nu := = 79.9
Nusselt : Formule de 1+
(
12.7 × Pr
2¸3
-1 )
Petukhov-Gnielinski ( 0.79 × ln ( Re) - 1.64) × 8

donnée pour l'eau :

Nu × λ W
hc := hc = 14.97 × W
D 2 he := 5225 ×
m ×K 2
m ×K

1 1 1 1 W
Calcul du coefficient de transfert global = + U := = 14.9 ×
U he hc 1 1 2
+ m K
he hc

MH 1/3 27/05/2013
univdocs.com TF06_convection_06.xmcd

Débit de chaleur transféré ( )

Φ := C2 × T 2S - T 2E = 6.714 × kW 1=chaud
2 = froid
Φ = U × S × Δθ = U × S × DTML

ΔT 0 - ΔT L ΔT L - ΔT 0
Rappel : DTML = =
æ ΔT 0 ö æ ΔT L ö
ln ç ÷ ln ç ÷
è ΔT L ø è ΔT 0 ø

Écoulement à co-courant Écoulement à contre-courant

100 100

80 80

60 60

40 40

20 20

0 0

- 20 - 20
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1

Écoulement à co-courant Écoulement à contre-courant

DTMLcoc :=
( T1E - T2E) - ( T1S - T2S) = 33.8 × °C DTMLctr :=
( T1E - T2S) - ( T1S - T2E) = 43.3 × °C
æ T1E - T2E ö æ T1E - T2S ö
ln ç ÷ ln ç ÷
è T1S - T2S ø è T1S - T2E ø

( 80 + 5) - ( 64 - 55) = 33.8
( 64 + 5) - ( 80 - 55) = 43.3
æ 80 + 5 ö æ 64 + 5 ö
ln ç ÷ ln ç ÷
è 64 - 55 ø è 80 - 55 ø

Remarque : DTML est une différence de température, et peut donc s'exprimer en °C ou en K.

W
Rappel : U = 14.9 ×
2
m °C

Φ 2 Φ 2
Scoc := = 13.3 m Sctr := = 10.4 m
U × DTMLcoc U × DTMLctr

6714
= 13.3 6714
14.9 ´ 33.8 = 10.4
14.9 ´ 43.3

L'écoulement à contre-courant est (toujours) plus performant !

MH 2/3 27/05/2013
univdocs.com TF06_convection_06.xmcd

Onpeutégal ementu l i
serl
amé
thodeduNUT(
Nombr
edUni
tésdeT
rans
fer
t)
On dit NTU en anglais.

On commence par définir les débits de capacité

-1
thermique (appelés aussi débits calorifiques) : (chaud) C1 = 419 × W × K
débit massique × capacité thermique massique -1
C=m∙ c (froid) C2 = 111.9 × W × K
P =ρ ∙
V∙cP

La puissance échangée s'écrit alors : ( ) (

Φ = C1 × T 1E - T 1S = C2 × T 2S - T 2E )
( )
3
Φ := C2 × T 2S - T 2E = 6.714 ´ 10 W

Ondét
erminelepl uspe tetl epl us Cmin := min ( C1 , C2) = 1 × C2 ( )
Cmax := max C1 , C2 = 1 × C1
grand de ces paramètres

Cmin
Ainsi que le rapport des débits de CR := = 0.267
Cmax
capacité thermique

( )
3
Onév
aluel
edébi
téner
gé quemax imal Φmax := Cmin × T 1E - T 2E = 9.511 ´ 10 W
échangé

Φ
ε := = 70.6 × %
Calcul de l'efficacité Φmax

Calcul du NUT.
A en onl
esr
ela onsentr
eεetNUTdépendentdel
adi
spos
ion
co-courant contre-courant

NUTcoc :=
ë (
-ln é1 - ε × 1 + CR ù
û) = 1.774 NUTctc :=
1 æ 1 - ε × CR ö
× ln ç ÷ = 1.385
1 + CR 1 - CR è 1-ε ø

U× S
NUT = W
Cmin Rappel : U = 14.9 ×
2
m °C

NUTcoc × Cmin 2 NUTctc × Cmin 2

Scoc := = 13.3 m Sctc := = 10.4 m
U U

On retrouve bien les valeurs obtenues avec la DTML

εcoc =
(
1 - exp éë-NUT × 1 + CR ùû ) εctc =
(
1 - exp éë-NUT × 1 - CR ùû )
1 + CR (
1 - CR × exp éë-NUT × 1 - CR ùû )

MH 3/3 27/05/2013
univdocs.com TF06_convection_07.xmcd

TF06 - Convection - Exercice 7

D4 ‚
Définition des paramètres : D3
‚
1 = nitrobenzène (chaud)
D2
•
2 = eau (froid) D1 •
λP := 46.4×
W
(paroi) ‚
m× °C
kJ W kg kg
CP1 := 1.38 μ1 := 1× mPa× s λ1 := 0.159 ρ1 := 1200× m1 := 3000×
kg× °C m× °C m
3 h

kJ W kg kg
CP2 := 4.18 μ2 := 0.7× mPa× s λ2 := 0.63 ρ2 := 1000× m2 := 1500×
kg× °C m× °C m
3 h

D1 := 25× mm D2 := 33× mm D3 := 50× mm D4 := 60× mm débits massiques

W W
T1E := 80°C T1S := 30°C T2E := 20°C h1 := 1410× h2 := 1935×
2 2
m × °C m × °C
Température de sortie de l'eau : ( )
m1× CP1× T1E - T1S = m2× CP2× T2S - T2E ( )
m1× CP1 3000 ´ 1.38
T2S := T2E +
m2× CP2
( )
× T1E - T1S = 53.0 °C 20 + ( 80 - 30) ×
1500 ´ 4.18
= 53.014

80°C 80°C Seule l'alimentation à contre-courant est possible !

53°C
53°C
30°C
30°C

Δθ
( )
20°C 20°C
Φ = m1× CP1× T1E - T1S = U× S× Δθ =
1
U× S
Débit de chaleur transféré

( )
Φ := m1× CP1× T1E - T1S = 57.5 kW
3000
3600
× 1380× ( 80 - 30) = 57500

Calcul du coefficient de transfert global la surface d'échange de référence est S2=pD2Le

1 1 1 1 æ D2 ö 1 W
= + + ln ç ÷ U := = 644.5
U × S2 h1× S1 h2× S2 2π Le× λP æ D2 ö 2
è D1 ø 1
+
D2
+ × ln ç
D2
m °C
÷
h2 h1× D1 2× λP
è D1 ø
-1
æ 1 + 33 0.033
× ln æç ö÷ ÷ö
33
ç + = 644.457
è 1935 25 ´ 1410 2 ´ 46.4 è 25 ø ø
Φ = U× S× Δθ = U× π× D2× Le× DTML
DTML =
ΔT 0 - ΔT L
DTML :=
( T1E - T2S) - ( T1S - T2E) = 17.1 °C
æ ΔT 0 ö æ T1E - T2S ö
ln ç ÷ ln ç ÷
è ΔT L ø è T1S - T2E ø
( 80 - 53) - ( 30 - 20)
= 17.1
ln æç
80 - 53 ö
÷
è 30 - 20 ø
57500 Φ
= 50.3 Le := = 50.3 m
644.5× π ´ 0.033× 17.1 U× π× D2× DTML

MH 1/2 15/05/2012
univdocs.com TF06_convection_07.xmcd

On peut également utiliser la méthode du NUT (Nombre dUnités de Transfert)

On dit NTU en anglais.

On commence par définir les débits de capacité 31

thermique (appelés aussi débits calorifiques) : C1 := m1× CP1 C1 = 1.15 ´ 10 ×W
K
débit massique × capacité thermique massique
C = m·cP = ρ·V·cP 31
C2 := m2× CP2 C2 = 1.742 ´ 10 ×W
K

La puissance échangée s'écrit alors : (

Φ := m1× CP1× T1E - T1S ) Φ = 57.5× kW

On détermine le plus petit et le plus grand

de ces paramètres
(
Cmin := min C1 , C2 ) (
Cmax := max C1 , C2 )

Cmin
Ainsi que le rapport des débits de capacité CR := CR = 0.66
thermique Cmax

On évalue le débit énergétique maximal

échangé
Φmax := Cmin× T1E - T2E( ) Φmax = 69.0× kW

Φ
Calcul de l'efficacité ε := ε = 83.3× %
Φmax

Calcul du NUT.
Attention les relations entre ε et NUT dépendent de la disposition
contre-courant

1 æ 1 - ε× CR ö
NUT := × ln ç ÷ = 2.922
1 - CR è 1 - ε ø
U× S
NUT =
Cmin
NUT× Cmin
Lta := = 50.3 m
U × π× D 2

On retrouve bien les valeurs obtenues avec la DTML

MH 2/2 15/05/2012
univdocs.com TF06_convection_08.xmcd

TF06 - Convection - Exercice 8

Définition des paramètres :
kJ kg
CPa := 1.5 Ta1 := 174°C Ta2 := 45°C Ta3 := 88°C ma := 4500×
kg× °C h
(côté gauche) (côté droite)
kJ W W
CPe := 4.18 Te1 := 80°C Te2 := 20°C U1 := 1160× U2 := 731×
kg× °C 2
m × °C
2
m × °C

Calcul du débit d'eau de refroidissement ( )

ma× CPa× Ta1 - Ta2 = me× CPe× Te1 - Te2 ( )
me := ma×
(
CPa× Ta1 - Ta2 ) me = 4500×
1.5× ( 174 - 45)
me = 3471.9×
kg
(
CPe× Te1 - Te2 ) 4.18× ( 80 - 20) h

ma× CPa Te1 - Te2

=
remarque : me× CPe Ta1 - Ta2

Température de l'eau entre les 2 réservoirs ( )

ma× CPa× Ta1 - Ta3 = me× CPe× Te1 - Te3 ( )
ma× CPa Te1 - Te2
(
Te3 := Te1 - Ta1 - Ta3 × )
me× CPe
d'après la remarque
ci-dessus:
(
Te3 := Te1 - Ta1 - Ta3 ×
Ta1 - Ta2
)
= 40 °C

80 - 20
80 - ( 174 - 88) × = 40
174 - 45

Φ = U× S× Δθ = ma× CPa× ΔT a (
Φ1 := ma× CPa× Ta1 - Ta3 = 161.25 kW )
Φ
S=
U× DTML
(
Φ2 := ma× CPa× Ta3 - Ta2 = 80.625 kW )

4500× 1500× ( 174 - 88) 4500× 1500× ( 88 - 45)

S1 = S2 =
( 174 - 80) - ( 88 - 40) ( 88 - 40) - ( 45 - 20)
3600× 1160× 3600× 731×
ln æç
174 - 80 ö
ln æç
88 - 40 ö
÷ ÷
è 88 - 40 ø è 45 - 20 ø

DTML1 :=
( Ta1 - Te1) - ( Ta3 - Te3) DTML1 = 68.443× °C S1 :=
Φ1 2
= 2.03 m
é ( Ta1 - Te1) ù U1× DTML1
ln ê ú
ë ( Ta3 - Te3) û

DTML2 :=
( Ta3 - Te3) - ( Ta2 - Te2) DTML2 = 35.258× °C S2 :=
Φ2 2
= 3.13 m
é ( Ta3 - Te3) ù U2× DTML2
ln ê ú
ë ( Ta2 - Te2) û

MH 1/1 15/05/2012
univdocs.com TF06_convection_09.xmcd

TF06 - Convection - Exercice 9 Échangeur thermique T1E := 85× °C T2E := 5× °C

ep := 2× mm Le := 36× m

D1 := 2.5× cm D4 := 4.2× cm

D2 := D1 + 2× ep D3 := D4 - 2× ep

D2 = 2.9× cm D3 = 3.8× cm
L L
V1 := 60× V2 := 30×
V1 et V2 sont les débits volumiques : min min

Définition des paramètres (eau à 45°C):

kg -3 kJ W W
ρ := 990× μ := 0.65× 10 × Pa× s cP := 4.18× λ := 0.602× λP := 45.2×
m
3 kg× °C m× °C m× °C
conductivité
Différentes définitions du débit d'énergie échangée entre les 2 fuides : thermique du tube

(
Φ = ρ× V1× cP× T1E - T1S ) (
Φ = ρ× V2× cP× T2S - T2E ) Φ = U× S× Δθ
Δθ =
ΔT 0 - ΔT L

æ ΔT 0 ö
æ S2 ö 1 ln ç ÷
1 1 1 1 U =
= + × ln ç ÷ + avec S = S2
1 D2 D2 æ D2 ö è ΔT L ø
U× S h1× S1 2× π Le× λP S1 h2× S2
è ø + + × ln ç ÷
h2 h1× D1 2× λP
S1 := π× D1× Le S2 := π× D2× Le
è D1 ø Prandtl :

4 × ρ× V 1 cP× μ
Calcul du coefficient de transfert h1 : Reynolds : Re1 := = 77570 turbulent Pr := = 4.513
π× μ × D 1 λ
Nusselt ( Re1-1000) × Pr
8× ( 0.79× ln ( Re1 ) - 1.64)
2
Nu1 := = 396.3 Nu1× λ
1+
12.7× Pr( 2¸3
-1 ) h1 :=
D1
= 9.542×
kW
2
( 0.79× ln( Re1) - 1.64) × 8 m × °C

Définition des critères adimensionnels dans l'espace annulaire (entre D2 et D3) ρ× u2× DH
Re =
Calcul du coefficient de transfert h2 : μ
4× Ω æç π× D 2 π× D2 ö÷
2
Définition du diamètre hydraulique DH : DH = 3
Pmouillé 4× ç - ÷
DH := D3 - D2 DH = 0.900× cm DH =
è 4 4 ø = D -D
π× D 3 + π× D 2 3 2

[ debit ] V2
u2 = = 4× V 2 D3 - D2 4 × ρ× V 2
[ section ] 2 2
Re = ρ× × Re2 := = 14472
π× D 3 π× D 2
4
-
4
π× æ D 3 - D 2 ö
è
2 2
ø
μ (
π× μ × D 3 + D 2 )
turbulent
( Re2-1000) × Pr h× De 4× Ω
8× ( 0.79× ln ( Re2 ) - 1.64)
2 Nusselt Nu = De =
Nu2 := = 93.6 λ Ptransfert

1+
12.7× Pr ( 2¸3
-1 ) De : diamètre équivalent
( ( )
0.79× ln Re2 - 1.64 × 8 ) æç π× D 2 π× D 2 ö÷
3 2
ou diamètre thermique
4× ç - ÷ Nu2× λ
De :=
è 4 4 ø
= 2.1× cm h2 := = 2.709×
kW
π× D 2 De 2
m × °C
1 kW
U := = 1.858× 2
D2 D2 æ D2 ö 2 S := S2 = 3.28 m
1 m × °C
+ + × ln ç ÷
h2 h1× D1 2× λP 3 -1
è D1 ø U× S = 6.093 ´ 10 × W× K

MH 1/3 10/06/2013
univdocs.com TF06_convection_09.xmcd

Étude à co-courant : ΔT 0 := T1E - T2E ΔT L = T1S - T2S C1 := ρ× V1× cP C2 := ρ× V2× cP

( )
Φ = C1× T1E - T1S = C2× T2S - T2E( ) 3
C1 = 4.138 ´ 10 × W× K
-1 3
C2 = 2.069 ´ 10 × W× K
-1

Φ=
( T1E - T1S) + ( T2S - T2E) Φ=
ΔT 0 - ΔT L
Φ=
ΔT 0 - ΔT L 1
+
1
=
1 æ ΔT 0 ö
× ln ç ÷
1
+
1 1
+
1 1 æ ΔT 0 ö C1 C2 U× S
è ΔT L ø
C1 C2 C1 C2 × ln ç ÷
U× S ΔT L
è ø
ΔT 0 - ΔT L
ΔT L := ΔT 0× exp éê-æç
1 ö
× U× Sùú
1
+ Φ := = 109.0× kW
C1 C2 ÷ ΔT L = 1.0× °C 1 1
ë è ø û +
C1 C2

Φ Φ
T1S := T1E - T2S := T2E +
C1 C2

T1E = 85× °C T1S = 58.7× °C co-courant :

T2E = 5× °C T2S = 57.7× °C Φ = 109.0× kW

DTML = 17.9× °C

Étude à contre-courant : ΔT 0 = T1E - T2S ΔT L = T1S - T2E

( ) (
Φ = C1× T1E - T1S = C2× T2S - T2E ) (A)

Φ=
( T1E - T1S) - ( T2S - T2E) = ( T1E - T2S) - ( T1S - T2E) Φ = U × S×
( T1E - T2S) - ( T1S - T2E) (B)
1
-
1 1
-
1 æ T1E - T2S ö
C1 C2 C1 C2 ln ç ÷
è T1S - T2E ø
T1E - T2S
= exp éêU× S× æç
1 öù
Ke := exp éêU× S× æç
1 1 1 öù
(B) - ÷ú Posons : - ÷ú Ke = 0.229379
T1S - T2E
ë è C1 C2 øû ë è C1 C2 øû
T2S - T2E C1 C1
(A) = Posons : Kv := Kv = 2
T1E - T1S C2 C2

(B) Ke× T1S + T2S = T1E + Ke× T2E en soustrayant (B) de (A) :
(A) Kv× T1S + T2S = Kv× T1E + T2E ( Kv - Ke) × T1S = ( Kv - 1) × T1E + ( 1 - Ke) × T2E

( Kv - 1) × T1E + ( 1 - Ke) × T2E

T1S := T2S := Kv× T1E + T2E - Kv× T1S
Kv - Ke

( 2 - 1) ´ 85 + ( 1 - 0.229379) ´ 5
= 50.2 2 ´ 85 + 5 - 2 ´ 50.2 = 74.6
2 - 0.229379

( )
Φ := C1× T1E - T1S = 144.1× kW

T1E = 85× °C T1S = 50.2× °C contre-courant :

T2S = 74.6× °C T2E = 5× °C Φ = 144.1× kW

DTML = 23.6× °C

MH 2/3 10/06/2013
univdocs.com TF06_convection_09.xmcd

On peut également utiliser la méthode du NUT (Nombre dUnités de Transfert)

On dit NTU en anglais
Les calculs sont considérablement simplifiés.

3 -1
On commence par définir les débits de capacité C1 := ρ× V1× cP C1 = 4.138 ´ 10 × W× K
thermique (appelés aussi débits calorifiques) :
débit massique × capacité thermique massique 3 -1
C = m·cP = ρ·V·cP C2 := ρ× V2× cP C2 = 2.069 ´ 10 × W× K

La puissance échangée s'écrit alors : (

Φ = C1× T1E - T1S = C2× T2S - T2E ) ( )
On détermine le plus petit et le plus grand
de ces paramètres
Cmin := min C1 , C2 ( ) (
Cmax := max C1 , C2 )
On en déduit le débit énergétique maximal
échangé
Φmax := Cmin× T1E - T2E ( ) Φmax = 165.5× kW

Cmin
Ainsi que le rapport des débits de capacité CR := CR = 0.5
thermique Cmax

U× S
On en déduit le nombre d'unités de transfert (NUT NUT := NUT = 2.945
Cmin

Toutes ces valeurs ci-dessus sont les mêmes en co-courant et à

Φ
contre-courant. ε=
Rappelons par ailleurs la définition de l'efficacité, dont les valeurs sont Φmax
différentes selon la disposition de l'échangeur..

Étude à co-courant : εcoc :=

(
1 - exp é-NUT× 1 + CR ù
ë û ) εcoc = 0.659
1 + CR

εcoc = 65.9× %
Φ := εcoc× Φmax Φ = 109.0× kW

Φ Φ
T1S := T1E - T2S := T2E +
C1 C2

T1S = 58.7× °C T2S = 57.7× °C DTML = 17.9× °C

Étude à contre-courant : εctc :=

(
1 - exp é-NUT× 1 - CR ù
ë û ) εctc = 0.87
(
1 - CR× exp é-NUT× 1 - CR ù
ë û )
εctc = 87.0× %
Φ := εctc× Φmax Φ = 144.1× kW

Φ Φ
T1S := T1E - T2S := T2E +
C1 C2

T1S = 50.2× °C T2S = 74.6× °C DTML = 23.6× °C

On retrouve bien les valeurs calculées avec la méthode DTML !

MH 3/3 10/06/2013
univdocs.com TF06_convection_10.xmcd

D1 := 22× mm
TF06 - Convection - Exercice 10
D2 := 25× mm
Définition des paramètres :
1 = benzène (chaud) W
2 = eau (froid) λP := 45×
m× K

kJ W tonne
CPb := 1.9× hb := 1700× TbE := 80× °C TbS := 30× °C mb := 4.5×
kg× K m ×K
2 h

kJ W
CPe := 4.185× he := 850× TeS := 50× °C TeE := 20× °C débit massique
kg× K 2
m ×K
T1E := TbE T1S := TbS
Installation à contre-courant :
T2S := TeS T2E := TeE

Débit de chaleur transféré ( )

Φ = mb× CPb× TbE - TbS = me× CPe× TeS - TeE = U× S× Δθ( ) Φ = 118.75× kW

débit d'eau CPb æ TbE - TbS ö tonne

me := mb× ×ç ÷ = 3.405
CPe TeS - TeE h
è ø
Calcul du coefficient de transfert global la surface d'échange de référence est S2=pD2Lt

1 1 1 1 æ D2 ö U :=
1
= 531.8
W
= + + ln ç ÷ D2 D2 æ D2 ö 2
U × S2 hb× S1 he× S2 2π Lt× λP 1 m °C
è D1 ø + + × ln ç ÷
he hb× D1 2× λP
è D1 ø
-1
æ 1 + 25 0.025 æ 25 ö ö
Φ = U× S× Δθ ç + × ln ç ÷ ÷ = 531.8
Méthode DTML è 850 1700 ´ 22 2 ´ 45 è 22 ø ø

ΔT 0 - ΔT L Φ 2
Δθ = DTML = ( 80 - 50) - ( 30 - 20) Δθ = 18.2× °C S := = 12.266 m
æ ΔT 0 ö = 18.20 U× Δθ
ln ç ÷ æ
ln ç
80 - 50 ö
÷ Φ
è ΔT L ø è 30 - 20 ø Lt := = 156 m
U× π× D2× Δθ

Méthode NUT Comme on connait toutes les températures, on connait Φ et on en déduit donc l'efficacité ε, on en
déduit alors NUT et donc la surface d'échange.

3 -1 3 -1
C1 := mb× CPb = 2.375 ´ 10 W× K C2 := me× CPe = 3.958 ´ 10 W× K

(
Φ := C1× T1E - T1S = 118.75 kW ) (
Cmin := min C1 , C2 = 1 C1 ) (
Cmax := max C1 , C2 = 1 C2 )
Cmin
CR :=
Cmax
= 0.6 (
Φmax := Cmin× T1E - T2E = 142.5 kW )

Φ C1 T1E - T1S Φ 1 æ 1 - ε× CR ö
ε= = × ε := = 83.3 % NUT := × ln ç ÷ = 2.747
Φmax Cmin T1E - T2E Φmax 1 - CR è 1 - ε ø

U× S NUT× Cmin NUT× Cmin

2
NUT = S := = 12.3 m Lta := = 156 m
Cmin U U × π× D 2

MH 1/2 18/05/2012
univdocs.com TF06_convection_10.xmcd

Système 2 calandres - 4 passes Φ = U× S× CF× Δθ = U× π× D2× Lte × CF× DTML

TbS - TbE TeE - TeS

intérieur = benzène (chaud) P := P = 0.833 R := R = 0.6
extérieur = eau (froid) TeE - TbE TbS - TbE

facteur de correction : CF := 0.75

Φ 2 1 S2
S2 := = 16.4 m Lt4 := × = 52.1 m
U× CF× Δθ 4 π× D 2

La méthode NUT est plus complexe car le calcul de NUT à partie de ε est plutôt compliqué. (5 étapes)
ε = 83.3× %

éê æ 1 - ε× CR ö ùú
2× ç C R - ÷
ê è 1-ε ø ú
ê 1 + CR + ú
1 - ε× CR
ê -1
ú
ê 1-ε ú
ê - 1ú
2
2 ê 1 + CR ú
NUT2_04 := × ln ê ú = 3.596
1 + CR
2 ê æ 1 - ε × C R ö ú
ê 2× ç C R - ÷ ú
è 1-ε ø
ê 1 + CR + ú
ê 1 - ε× CR ú
ê -1 ú
1-ε
ê + 1ú
ê 1 + CR
2 ú
ë û

NUT2_04 = 3.596

U× S
NUT = 1 NUT2× Cmin Lt4
Cmin Lta4 := × = 51.1 m
4 U × π× D 2 = 1.018
Lta4

Soit un écart de moins de 2%

MH 2/2 18/05/2012
univdocs.com TF06_convection_11.xmcd

TF06 - Convection - Exercice 11 Réservoir de vapeur cylindrique TK º 273.15× K

Définition des paramètres :
5
bar º 10 × Pa
Lr := 5× m Dr := 50× cm Mv := 19.3× kg Ta := 15°C

2
π× D r × L r Vr 3
-3 m
Volume spécifique : Vr := Vr = 981.7477 L Vvs := Vvs = 50.868 ´ 10
4 Mv kg

( )
3
m 2.038 3 -1
Relation de Bertin Pvs× Vvs = 2.038× bar× Pvs := × bar× m × kg Pvs = 40.065× bar
kg Vvs

4
4
Pvs
Relation de Duperay Tvs = 100× 0.981× Pvs Tvs := 100× 0.981× × °C Tvs = 250.4× °C
bar

Tvs = 523.5 K
Évaluation de la puissance électrique nécessaire P = ΦCV + ΦR

hc× D Ta = 288.1 K
ΦCV = hc× Slat× Tvs - Ta ( ) Nu =
λ
n
= C× ( Gr× Pr) = C× Ra
n

Ta + Tvs
On prend la moyenne arithmétique de T a et Tvs Tm := = 405.8 K
2
et on interpole :

37.90 - 25.90 æ Tm ö -6 2 -1
25.90 + ×ç - 400÷ = 26.60 νm := 26.60× 10 ×m ×s
500 - 400 è K ø

0.04038 - 0.03365 æ Tm ö W
0.03365 + ×ç - 400÷ = 0.03404 λm := 0.03404×
500 - 400 è K ø m× °C

0.680 - 0.689 æ Tm ö
0.689 + ×ç - 400÷ = 0.688 Prm := 0.688
500 - 400 è K ø
3
g× β× D 1 -2
Calcul du nombre de Grashof Gr = × ΔT β := g := 9.81× m× s
2 Tm
ν
3
g× β× Dr
Gr :=
2
( )
× Tvs - Ta = 1.005 ´ 10
9
Ra := Gr× Prm = 6.916 ´ 10
8
νm
régime laminaire

.25
Nu× λm W
Nu := 0.53× Ra = 85.947 hc := = 5.851
Dr 2
m ×K

-8 -2 -4 (
ΦCV := hc× π× Dr× Lr× Tvs - Ta = 10.817 kW )
σ := 5.67× 10 × W× m ×K

ΦR = σ× Slat× æ Tvs - Ta ö
4 4
è ø

ΦR := σ× π× Dr× Lr× é Tvs

ë ( ) 4 - ( Ta) 4ùû = 30.384 kW
Finalement, la puissance électrique nécessaire :

P := ΦCV + ΦR = 41.202 kW

MH 1/1 18/05/2012