TD 2 T Thérmique

m
co
-p ro.
TRANSFERT DE CHALEUR
ari
Volume 2 :
k
ab
Travaux dirigés et exercices

al3
w.
ww
ww
w.
al3
ab
kari
-pro.
co
m
TDEX-2.1
Travaux dirigés et exercices du chapitre 2
m
Pour ces deux exercices, l’approche est similaire à celle des notes de cours p 2.3 mais
faites attention aux expressions du gradient (p 2.1) et ne simplifiez trop vite par r !
co
EX-2.1 (S) Dérivation de l’équation thermique en coordonnées cylindriques
-p ro.
ari
En utilisant le volume de contrôle ci-dessus ∆r ∆z r ∆φ et en supposant une génération
de chaleur q′′′ (W/m3), faites le bilan de chaleur et retrouver l’équation de diffusion
thermique en cylindrique.
k
ab
EX-2.2(S) Dérivation de l’équation thermique en coordonnées spériques

al3
w.
En utilisant le volume de contrôle ci-dessus ∆r r ∆ θ r sin θ ∆φ et en supposant une

génération de chaleur q′′′ (W/m3), faites le bilan de chaleur et retrouver l’équation de
ww
diffusion thermique en sphérique

TDEX-2.2
Solutionnaire
m
2.1 Dérivation de l’équation thermique en coordonnées cylindriques
co
1) Hypothèses : T = Tr ,φ , z
2) Volume de contrôle ∆r ∆z r ∆φ
3) Bilan sur le volume de contrôle
ro.
Ce qui rentre – Ce qui sort + Ce qui est généré = Ce qui s’accumule
a) IN/OUT dans la direction r (au travers des faces perpendiculaires à l’axe des r):
∆z r ∆φ qr" − ∆z r ∆φ qr"
-p r r +∆r
b) IN/OUT dans la direction z (au travers des faces perpendiculaires à l’axe des z):
+ ∆r r ∆φ q"z − ∆r r ∆φ q"z
ari
z z +∆z
c) IN/OUT dans la direction φ (au travers des faces perpendiculaires à l’axe des φ ):
+ ∆r ∆z qφ" − ∆r ∆z qφ"
φ φ +∆φ
k
d) Terme de génération (dans le volume):

+ ∆r ∆z r ∆φ q′′′
ab
e) Terme d’accumulation (dans le volume) :

∂T
= ∆r ∆z r ∆φ ρ CP
∂t
al3
On divise par ∆r ∆z ∆φ et on fait tendre les incréments vers 0
Terme a :
w.
∆z r ∆φ qr" − ∆z r ∆φ qr" r qr" − r qr"

r r +∆r
= r r +∆r
∆r ∆z ∆φ ∆r
rqr" − rqr" ∂
lim r r +∆r
=− ( rqr" )
ww
∆r → 0 ∆r ∂r
TDEX-2.3
m
Terme b :
∆r r ∆φ q"z − ∆r r ∆φ q"z r q"z − r q"z q"z − qz"
z z +∆z
= z z +∆z
=r z z +∆z
∆r ∆z ∆φ ∆z ∆z
co
q"z − qz" ∂ "
lim
∆z → 0
z
∆z
z +∆z
=−
∂z
( qz )
ro.
Terme c :
∆r ∆z qφ" − ∆r ∆z qφ" qφ" − qφ"
φ φ +∆φ φ φ +∆φ
+ =
∆r ∆z ∆φ ∆φ
-p
qφ" − qφ" ∂ "
lim
φ
∆ →0
φ
∆φ
φ +∆φ
=−
∂φ
( qφ )
ari
Terme d :
∆r ∆z r ∆φ q′′′
+ = r q′′′
∆r ∆z ∆φ
Terme e :
k
∂T
∆r ∆z r ∆φ ρ C P
∂t = r ρ C ∂T
ab
P
∆ r ∆ z ∆φ ∂t
L’équation s’écrit donc :

al3
∂ ∂ ∂ " ∂T
−
∂r
( rqr" ) − r ( q"z ) −
∂z ∂φ
( qφ ) + r q′′′ = r ρ CP
∂t
La densité de flux s’exprime par la loi de Fourrier en utilisant l’expression du

gradient en coordonnées cylindriques :
w.
∂ ∂ 1 ∂
∇ = δr +δz + δφ
∂r ∂z r ∂φ
∂T ∂T 1 ∂T
qr" = − k q"z = − k qφ" = − k
ww
∂r ∂z r ∂φ
TDEX-2.4
m
∂ ∂T ∂ ∂T ∂ 1 ∂T ∂T
− r −k − r −k − −k + r q′′′ = r ρ CP
∂r ∂r ∂z ∂z ∂φ r ∂φ ∂t
On divise par r:
co
∂ ∂T ∂ ∂T ∂ 1 ∂T ∂T
− r −k − r −k − −k + r q′′′ r ρ CP
∂r ∂r ∂z ∂z ∂φ r ∂φ ∂t
=
r r
ro.
∂ ∂T ∂ ∂T ∂ 1 ∂T
r k r k k
∂r ∂r ∂z ∂z ∂φ r ∂φ ∂T
+ + + q′′′ = ρ CP
r r r ∂t
(commentaire: dans la dérivée -p

∂
∂φ
1
, le terme ne varie pas avec φ )
r
ari
et on obtient finalement:
1 ∂ ∂T ∂ ∂T 1 ∂ ∂T ∂T
kr + k + 2 k + q′′′ = ρ CP
r ∂r ∂r ∂z ∂z r ∂φ ∂φ ∂t
k
ab
al3
w.
ww
TDEX-2.5
2.2 Dérivation de l’équation thermique en coordonnées spériques
m
1) Hypothèses : T = Tr ,φ ,θ
co
2) Volume de contrôle ∆r r ∆ θ r sin θ ∆φ
ro.
a) IN/OUT dans la direction r (au travers des faces perpendiculaires à l’axe des r):
-p
r ∆ θ r sin θ ∆φ qr" − r ∆ θ r sin θ ∆φ qr"
r
b) IN/OUT dans la direction φ (au travers des faces perpendiculaires à l’axe des φ ):
+ ∆r r ∆ θ qφ" − ∆r r ∆ θ qφ"
r +∆r
ari
φ φ +∆φ
c) IN/OUT dans la direction θ (au travers des faces perpendiculaires à l’axe des θ ):
+ ∆r r sin θ ∆φ qθ" − ∆r r sin θ ∆φ qθ"
θ θ +∆θ
d) Terme de génération (dans le volume):
k
+ ∆r r ∆ θ r sin θ ∆φ q′′′
e) Terme d’accumulation (dans le volume) :
ab
∂T
= ∆r r ∆ θ r sin θ ∆φ ρ CP
∂t
On divise par ∆r ∆θ ∆φ et on fait tendre les incréments vers 0

al3
Terme a :
r ∆ θ r sin θ ∆φ qr" − r ∆ θ r sin θ ∆φ qr" r 2 qr" − r 2 qr"

= ( sin θ )
w.
r r +∆r r r +∆r
∆r ∆ θ ∆φ ∆r
r 2 qr" − r 2 qr" ∂ 2 "
lim
∆r → 0
r
∆r
r +∆r
=−
∂r
( r qr )
ww
Terme b :
TDEX-2.6
m
∆r r ∆ θ qφ − ∆r r ∆ θ qφ
" "
r qφ − r qφ
" "
qφ − qφ
" "
φ φ +∆φ φ φ +∆φ φ φ +∆φ
= == r
∆r ∆ θ ∆φ ∆φ ∆φ
qφ" − qφ" ∂ "
( qφ )
co
φ φ +∆φ
lim
φ
∆ →0 ∆φ
=−
∂φ
Terme c :
∆r r sin θ ∆φ qθ" − ∆r r sin θ ∆φ qθ" r sin θ qθ" − r sin θ qθ" sin θ qθ" − sin θ qθ"
ro.
θ θ +∆θ
= θ θ +∆θ
=r θ θ +∆θ
∆r ∆θ ∆φ ∆θ ∆θ
sin θ qθ" − sin θ qθ" ∂
lim
φ
∆ →0
θ
∆θ
θ +∆θ
=−
∂θ
( sin θ qθ" )
Terme d : -p
ari
∆r r ∆ θ r sin θ ∆φ q′′′
+ = r 2 sin θ q′′′
∆r ∆ θ ∆φ
Terme e :
k
∂T
∆r r ∆ θ r sin θ ∆φ ρ CP
∂t = r 2 sin θ ρ C ∂T
ab
P
∆r ∆ θ ∆φ ∂t
L’équation s’écrit donc :
∂ 2 " ∂ " ∂ ∂T
( r qr ) − r ( qφ ) − r ( sin θ qθ" ) + r 2 sin θ q′′′ = r 2 sin θ ρ CP
al3
− sin θ
∂r ∂φ ∂θ ∂t
La densité de flux s’exprime par la loi de Fourrier en utilisant l’expression du

gradient en coordonnées sphériques:
w.
∂ 1 ∂ 1 ∂
∇ = δr + δθ + δφ
∂r r ∂θ r sin θ ∂φ
∂T 1 ∂T 1 ∂T
qr" = − k qθ" = − k qφ" = − k
∂r r ∂θ r sin θ ∂φ
ww
TDEX-2.7
m
∂ 2 ∂T ∂ 1 ∂T ∂ 1 ∂T
− sin θ r −k − r −k − r ( sin θ )( −k )
∂r ∂r ∂φ r sin θ ∂φ ∂θ r ∂θ
∂T
+ r 2 sin θ q′′′ = r 2 sin θ ρ CP
co
∂t
∂ 2 ∂T 1 ∂ ∂T 1 ∂ ∂T
sin θ r k + r k + r k sin θ
∂r ∂r r sin θ ∂φ ∂φ r ∂θ ∂θ
ro.
∂T
+ r 2 sin θ q′′′ = r 2 sin θ ρ CP
∂t
On divise par r 2 sin θ :
1 ∂
r ∂r
2
k r2
∂T
∂r
+
1 ∂
r sin θ ∂φ
2 2
k
∂T
∂φ -p +
1 ∂
r sin θ ∂θ
2
k sin θ
∂T
∂θ
+ q′′′ = ρ CP
∂T
∂t
k ari
ab
al3
w.
ww
ww
w.
al3
ab
kari
-pro.
co
m
TDEX-3.1
om
Travaux dirigés et exercices du chapitre 3:
Conduction dans les solides en régime permanent
Exercices
Ex-3.1 Profils de température (à l=examen de 1997)
Une paroi de conductivité k sépare un fluide chaud à la températures Ta d'un fluide froid à la température Tb.
o.c
Sur chacune des faces, les coefficients de transfert de chaleur avec les fluides sont respectivement ha et hb.
Pour un régime stationnaire, on étudie la variation de la température le long d'un axe perpendiculaire à la
paroi. Pour chacune des trois situations suivantes A, B ou C, décrites ci-dessous, associer la figure qui lui
correspond. Cocher la case □
A) la conductivité de la paroi est très élevée 1□ 2□ 3□
et dans ce cas, quel est le coefficient le plus élevé (ha ou hb) ?
-pr
B) les deux coefficients ha et hb sont très élevés 1□ 2□ 3□
C) les coefficients sont moyens mais ha > hb 1□ 2□ 3□
k ari
ab
Ex-3.2(S) Mur avec différentes conditions frontières (voir figure notes de cours p 3.1)
Trouver l'expression du profil de température dans un mur d'épaisseur L, de conductivité k et dont la face à
x=0, est soumise à un flux constant q"O alors que l'autre à x=L échange sa chaleur par convection (coefficient
hF) avec un fluide à une température TF. (Refaire toute la démonstration: schéma, hypothèses, bilan,
al3
conditions frontières)
,,
q 1 L
Réponse: T(x)= - 0 x + q 0,,[ + ] + T F
k hF k
Ex-3.3 Epaisseur d'un isolant ( à l'examen de sept. 1989)
Vous souhaitez isoler la paroi intérieure d'un mur d'une maison. Ce mur en brique, d'épaisseur
w.
o
LB=0.10 m, a une conductivité thermique kB=0.72 W/m/ C. En utilisant le concept de résistance thermique (on
présentera une schéma du problème réel ainsi que le schéma correspondant des résistances thermiques),
déterminez l'épaisseur de la couche d'isolant (laine de verre, kI =0.043 W/m/C) à installer pour limiter les
2
pertes thermiques à une valeur de 50 W/m dans les conditions suivantes d'utilisation:
0
- température intérieure de la pièce T1 = 20 C
0
- température extérieure de l'air T2 = -20 C
ww
2
- coefficient de transfert de chaleur dans la pièce h1=10 W/m /C
2
- coefficient de transfert de chaleur à l'extérieur h2=100 W/m /C
Quelles seraient les pertes thermiques sans isolation ?
2
Réponses: épaisseur d'isolant= 2.37 cm, q"=161 W/m
TDEX-3.2
m
Travaux dirigés No 1
Ex-3.4(S) Conduction avec génération de chaleur (à l'examen final 1991)
co
Une sphère métallique de rayon R et de conductivité k est le siège d'une génération de chaleur uniforme q'''
3
(W/m ). Cette sphère échange avec l'environnement à la température TE et le coefficient de transfert de
chaleur à la surface vaut h.
a) Déterminer l'expression du profil de température dans la bille métallique.
b) Sans utiliser l'expression du profil obtenue en a) et par un raisonnement indépendant,
ro.
retrouver la température à la surface.
q′′′ 2 2 Rq′′′
Réponse a): T(r)= - ( r - R )+ +T E
6k 3h
Ex-3.5 Pertes thermiques d'un tube isolé (à l'examen de sept. 1990)

0
ri-p
Un fluide, à une température TF=150 C, s'écoule dans une tube d'acier (rayon intérieur r1=0.023 m, rayon
exterieur r2=0.026 m) de conductivité thermique kA= 12 W/m/K. Le coefficient de transfert de chaleur, hF entre
2
le fluide et la paroi interne du tube est égal à 450 W/m /K. Ce tube est recouvert d'un isolant (rayon externe
0
r3=.050 m) de conductivité thermique kISO=0.05 W/m/K. L'air environnant est à la température TAIR=20 C, et le
2
coefficient d'échange à la surface de l'isolant vaut hAIR=5 W/m /K.
a) Faites le schéma de ce problème et, en utilisant le concept de résistance thermique, faites le schéma
ka
de l'analogie électrique.
b) Trouvez l'expression des pertes thermiques, par unité de longueur de tube, pour les conditions
indiquées ci-dessus.
c) L'augmentation de l'épaisseur d'isolation réduirait-elle les pertes thermiques ?.
(PS: on négligera les résistances de contact) Réponses: b) pertes= 47.5 W/m c) oui, car (dq/dr3) < 0
b
Ex-3.6 Profil dans une résistance électrique (à l=examen de

1997)
3a
Une résistance électrique (épaisseur b, largeur W, hauteur L)

est isolée parfaitement sur une de ses faces. Sur l'autre face,
il y a échange par convection (avec un coefficient h) avec l'air
environnant (à la température T4). En régime permanent, on
suppose que la conductivité k de la résistance est constante
et que la génération de chaleur par effet Joule est uniforme et
.al
,,, 3
vaut q (W/m ).
Résistance
a) Déterminer l'expression du profil de température
dans la résistance électrique (on négligera les effets
de bouts). électrique
w
b) Sans utiliser l'expression du profil obtenue en a),

trouver la température sur la surface en contact avec
l'air.
c) A quelle profondeur dans la plaque électrique
ww
observe-t-on la température maximale ? Tracer sur la figure ci-dessus la forme du profil de

température dans la plaque.
q′′′ 2 1 b
Réponse a: T(x)= - x +q′′′b( + )+ T ∞
2k h 2k
TDEX-3.3
om
Conduction dans les ailettes
Ex-3.7(S) Température de l'anse d'un chaudron
L'anse d'un chaudron est une tige métallique de 10 mm de rayon

et de 0.5 m de longueur. Les deux parties inférieures de l'anse
0
c
sont soudées au chaudron qui est à une température de 90 C.
Sachant que le coefficient de transfert de chaleur avec l'air
2 0
ambiant est de 35 W/(m .K) et que l'air est à 30 C, déterminer la
ro.
température au milieu de l'anse si la tige est faite:
a) en cuivre, k=385 W/(m.K)

b) en acier, k=44 W/(m.K)
(on pourra négliger la courbure de l'anse et le gradient radial de

température)
ri-p
réponses: 66.9 et 35.1 C
Travail complémentaire:
0
c) Avec Excel par exemple, étudier l'effet de h et de k sur la

témperature au milieu de l=anse du chaudron (TC ). Tracer les
2
ka
graphes TC vs h (h variant entre 1 et 100 W/m .K) et TC vs k (k variant de 10 à 400 W/m.K). Commentez.
Ex-3.8(S) Ailette circulaire (à l'examen de 1989)

ab
Une ailette circulaire d'épaisseur W et de rayon exterieur R2 est installée sur un tube de rayon externe R1 dont
la surface est à la température T1. L'air environnant est à la
température T4 et le coefficient de convection à la surface de
l'ailette est h.
a) Faites un bilan de chaleur sur un élement de volume

al3
approprié de l'ailette et obtenez l'équation différentielle

que doit satisfaire la température de cette ailette.
(remarque: on ne demande pas d'intégrer cette équation)
b) Quelles sont les différentes conditions frontières

possibles ?
w.
c) Si on appelle qC, les pertes de chaleur (en watts) au

travers de cette ailette, donnez alors les expressions du
rendement et de l'efficacité de l'ailette.
Réponse a:
2
∂ T 1 ∂T 2h
ww
+ - (T -T∞ )= 0
∂ r 2 r ∂r kW
TDEX-3.4
om
Exercice
Ex-3.9(S) Profil de température dans une ailette fixée entre deux parois
Une ailette de largeur W, d'épaisseur t et de longueur L

est fixée entre deux pièces métalliques dont les Tp
températures sont identiques et égales à TP. Les deux
c
faces de l'ailette ne sont pas soumises aux mêmes
conditions. La face supérieure de l'ailette est soumise à
un courant d'air à une température TA et le coefficient
ro.
hA TA
d'échange vaut hA. La face inférieure quant à elle, est
soumise à un autre courant d'air à la température TB et
le coefficient d'échange vaut hB. Afin de déterminer le
profil axial de température,
hB Tb
(attention: on supposera que les extrémités latérales
ri-p
sont isolées et les régions A et B sont séparées; la
convection ne se fait que sur les faces horizontales ).
ka
a) Faites un bilan de chaleur sur un volume approprié. Posez clairement vos hypothèses et obtenez
l'équation différentielle que doit vérifier la température.
b) Posez les conditions frontières.
c) Obtenez l'expression générale du profil de température ainsi que les relations que doivent vérifier les
ab
constantes d'intégration figurant dans cette expression
Réponse a:
al3
2
d T - [T h A + h B - h ATA + h BT B ] = 0
2
dx kt kt
2
dT 1 d 2U
- [ α T - β ] = 0 ⇔ -U = 0
dx
2
α dx 2
h + hB +h
w.
h
avec α = A , β = AT A BT B et U = α T - β
kt kt
ww
TDEX-3.5
m
Exemple d’examen No 1
GCH-18243 Examen partiel No 1 - 25 pts Aucun document sauf 1 feuille (recto-verso)

19 Septembre 2003, Durée 110 min
co
I) Profil dans un mur de béton ( 6 pts)
Un mur de béton (épaisseur b, largeur W, hauteur

L, conductivité k ) sépare une pièce à la température TINT ,
ro.
de l’extérieur dont la température est T EXT . Les coefficients
de transfert de chaleur sur les faces intérieure et extérieure
de ce mur sont respectivement hINT et hEXT
a) Faites un bilan de chaleur sur un volume de

contrôle pertinent et obtenez l’équation différentielle
que doit vérifier la température dans le mur (on
négligera les effets de bouts) (2 pts)
ri-p
b) Intégrer cette équation, poser les conditions
frontières et déterminer l'expression du profil de température dans le mur (4. pts)
Réponse
T(x)= C1 x + C2 ; C2 = TINT +
k C1
; C1 =
(TEXT − TINT )
hINT  k k 
ka
 +b+ 
 hINT hEXT 
II) Stockage d’air liquide dans un réservoir sphérique (6 pts)

ab
Un réservoir sphérique de rayon r1=1.5 m, contient

de l’air liquide. Le réservoir est isolé par une épaisseur de
0.05 m d'un matériau isolant de conductivité k=0.05
W/(m.K). On peut supposer que la paroi externe du
réservoir reste à température constante T1=80 K. A la
surface de l'isolant, il y a échange par convection avec l'air
al3
ambiant à la température T∞=283 K et le coefficient de

2
convection h est égal à 18 W/(m .K).
a) En utilisant, le concept des résistances thermiques,

calculer la température T2, à la surface de l'isolant en
contact avec l'air. (4.5 pts) (On négligera la résistance
de contact paroi-isolant)
w.
(réponse 272.6K)
b) Est-ce que l’épaisseur de l’isolant vous parait suffisante ? Commenter. (1.5 pt)
On rappelle les expressions des résistances pour la sphère creuse:

ww
1 1 1 1
Rconduction =  −  Rconvection =
sphèrecreuse 4π k  r1 r2  sur une sphère 4π h r22
TDEX-3.6
om
III) Température d’un fil électrique ( 3 pts)
Un fil électrique de 1 mm de diamètre est dénudé sur une longueur de 1 m. Calculer la

température de ce fil (température qu’on supposera uniforme), sachant que la puissance générée par
o
effet joule est de 1.5 watt/m, que la température de l’air est de 20 C et que le coefficient de transfert
2
de chaleur à la surface du fil avec l’air environnant vaut 12. W/m .C.
Réponse : 59.8
o.c
IV) Quiz (2 pts) Cocher la bonne réponse:
• Pour la construction d’une ailette, il est préférable d’utiliser un matériau qui à un faible
conductivité. VRAI FAUX
-pr
• Dans une ailette, l’augmentation de la conductivité favorise la diminution du gradient de la
température. VRAI FAUX
• Dans un bilan de chaleur associé à un problème de conduction, si il y a un terme de génération de

chaleur alors le terme d’accumulation est non nul. VRAI FAUX
• À l’interface entre deux solides, plus la résistance de contact est faible, plus la différence entre les
ari
températures de ces deux surfaces est grande. VRAI FAUX
Réponses : FVFF
V) Ailette en forme de cylindre creux (8 pts)

k
Une ailette, de longueur L, a la forme

d'un cylindre creux. Sa base est fixée sur une
paroi à la température Tp. La face intérieure du
ab
cylindre est refroidie par de l’air à la température

T1 avec un coefficient de convection h1 alors que
la face extérieure est refroidie avec un air à la
température
T2 avec un coefficient de convection h2.
al3
On supposera que le flux de chaleur axial à x=L

sur l’extrémité de l’ailette est négligeable.
a) Faites un bilan de chaleur sur un élément de

volume approprié de l'ailette et obtenez l'équation
différentielle que doit satisfaire la température de
cette ailette. (4.5 pts)
w.
b) Quelles sont les conditions frontières ? (1 pt)
c)Obtenez l’expression du profil de température. (2.5 pts)

ww
TDEX-3.7
m
Solutionnaire Ex-3.2 Mur avec différentes conditions frontières
Trouver l'expression du profil de température dans un mur d'épaisseur L ( de surface A), de conductivité k et
co
dont la face à x=0, est soumise à un flux constant q"O alors que l'autre à x=L échange sa chaleur par
convection (coefficient hF) avec un fluide à une température TF. (Refaire toute la démonstration: schéma,
hypothèses, bilan, conditions frontières)
1) Hypothèses :
ro.
∂
T=T(x) k=Constante régime permanent =0
∂t
2) Volume de contrôle A ∆x
ri-p
Ce qui rentre – Ce qui sort =0
A q"x − Aqx" =0
x x +∆x
On divise par A∆x et on fait tendre ∆x vers 0
ka
A q"x − Aqx" qx" − qx"
x x +∆x x x +∆x
= =0
A∆x ∆x
q"x − q"x ∂q"x
ab
x x +∆x
lim =− =0
∆x → 0 ∆x ∂x
On intègre une première fois :
∂q"x
− = 0 ⇒ qx" = C1
∂x
al3
On utilise la loi de Fourier

∂T
qx" = − k = C1
∂x
Et on intègre une deuxième fois :
∂T C  C 
=− 1 ⇒ Tx =  − 1  x + C2 (1)
w.
∂x k  k 
Les constantes d’intégration C1 et C2 sont obtenues en écrivant les conditions frontières
 ∂T 
à x = 0 ⇒ q"x ⇒ q"x
ww
CF1 : = q0" =  −k "

 = C1 = q0
x =0 x=0
 ∂x  x =0
CF2 : à x = L ⇒ q"x = hF (Tx = L − TF )

x= L
TDEX-3.8
m
 ∂T    C   
q"x =  −k  = C1 = q0" = hF   − 1  x + C2  − TF 
x=L
 ∂x  x = L   k   x= L 
o
 q  
"
q0" = hF  − 0  L + C2 − TF 
 k  
o.c
 1 L
⇒ C2 = q0"  +  + TF
 hF k 
Et en remplaçant C1 et C2 dans (1), on obtient l’expression de la variation de la température dans le mur :
-pr
q0"  1 L
Tx = − x + q0"  +  + TF
k  hF k 
ri
ka
ab
al3
w.
ww
TDEX-3.9
m
Solutionnaire Ex-3.4 Conduction avec génération de chaleur
Une sphère métallique de rayon R et de conductivité k est le siège d'une génération de chaleur uniforme q'''
3
(W/m ). Cette sphère échange avec l'environnement à la température TE et le coefficient de transfert de
o
chaleur à la surface vaut h.
a) Déterminer l'expression du profil de température dans la bille métallique.
o.c
1) Hypothèses :
∂
T=T(r) k=Constante q′′′ = Constante régime permanent =0
∂t
2) Volume de contrôle 4π r 2 ∆r
-pr
4π r 2 qr" − 4π r 2 qr" + 4π r 2 ∆r q′′′ = 0

r r +∆r
ari
On divise par 4π ∆r et on fait tendre ∆r vers 0
4π r 2 qr" − 4π r 2 qr" 4π r 2 ∆r q′′′

r r +∆r
+ =0
4π ∆r 4π ∆r
k
r 2 qr" − r 2 qr" ∂ 2 "

Par définition : lim
∆r → 0
r
∆r
r +∆r
=−
∂r
( r qr ) on obtient donc :
ab
∂ ∂ 2 "
− ( r 2 qr" ) + r 2 q′′′ = 0 ⇒ ( r qr ) = r 2 q′′′
∂r ∂r
r3
On intègre une première fois: ( r 2 qr" ) = q′′′ + C1
al3
3
∂T
On exprime alors la densité de flux par la loi de Fourrier qr" = −k :
∂r
r3
3 q′′′ + C1
 ∂T  r ∂T ∂T r3 q′′′ C q′′′ C
r 2  −k  = q′′′ + C1 ⇒ =− 3 ⇒ = − − 12 = − r − 12
 ∂r  3 ∂r k r2 ∂r r2 3k kr 3k kr
w.
et on intègre une deuxième fois pour obtenir la variation T(r):
2
 q′′′  r C  1
Tr =  −  − 1  −  + C2 (1)
ww
 3k  2 k  r
On utilise les conditions frontières (ou de symétrie) pour déterminer les constantes d’intégration C1 et C2
TDEX-3.10
om
∂T
CF1: à r=0 la température est finie, ou encore =0 ⇒ C1 = 0
∂r
2
 q′′′  r
T = −  + C2 (2)
 3k  2
CF2: à r=R, la chaleur qui arrive à la surface par conduction sort par convection :
∂T
qr" = h (Tr = R − TE )
o.c
= −k
r=R ∂r r=R
et on remplace alors T par l’expression (2):
 q′′′  r 2 
∂  −  + C2 
∂T  3k  2  q′′′ 2r  q′′′   q′′′  r 2  
= −k  
-pr
qr" = −k = − k −  = R = h   −  + C2  − TE 

∂r ∂r  3 k 2  r=R 3   3k  2
r =R
r=R r=R 
r=R
q′′′   q′′′  R 2  q′′′  q′′′  R

2
R = h−  + C2 − TE  ⇒ R+  + TE = C2
ari
3   3k  2  3h  3k  2
q′′′R q′′′R 2
⇒ C2 = + + TE
3h 6k
En remplaçant C1 et C2 dans (1), l’expression de la variation de la température s’écrit alors :
bk
 q′′′  2 q′′′R
 (r − R ) +
2
Tr =  − + TE (3)
 6k  3h
b) Sans utiliser l'expression du profil obtenue en a) et par un raisonnement indépendant retrouver la
température à la surface.
3a
Le bilan de chaleur sur la sphère en entier s’écrit :
Ein − Eout + E généré = Eaccumulé

Ein = 0 Eaccumulé = 0
al
4
Eout = 4π R 2 h (Tr = R − TE ) et Egénéré = π R 3 q′′′
3
w.
4
⇒ −4π R 2 h (Tr = R − TE ) + π R 3 q′′′ = 0
3
4 3 
 π R q′′′ 
Tr = R = TE + 
3  = T + q′′′R
ww
( 4π R h ) E 3h
2
Cette expression est identique à celle obtenue avec l’équation (3) lorsque r=R.
TDEX-3.11
om
Solutionnaire Ex-3.7 Température de l'anse d'un chaudron
L'anse d'un chaudron est une tige métallique de 10 mm de rayon et de 0.5 m de longueur. Les deux
0
parties inférieures de l'anse sont soudées au chaudron qui est à une température de 90 C. Sachant que le
2 0
coefficient de transfert de chaleur avec l'air ambiant est de 35 W/(m .K) et que l'air est à 30 C, déterminer la
température au milieu de l'anse si la tige est faite:
a) en cuivre, k=385 W/(m.K)
c
b) en acier, k=44 W/(m.K)
(on pourra négliger la courbure de l'anse et le gradient radial de température)
ro.
Hypothèses :
On néglige la courbure de l’anse et le problème à résoudre est celui de la variation dans une tige
cylindrique de rayon R et de longueur L=0.5 m
ri-p Convection, h
T∞
avec air à
∆x
ka
x
o
x=0; T=Tp =90 C x=L; T=Tp =90oC
ab
∂
T=T(x) k=Constante Eg=0 régime permanent =0
∂t
al3
Volume de contrôle π R 2 ∆x
Bilan sur le volume de contrôle
Ce qui rentre(conduction) – Ce qui sort (conduction+convection à la surface de la tige)=0

w.
π R 2 qx" x − π R 2 qx" x +∆x

− 2π R∆x h(Tx − T∞ ) = 0
On divise par π R 2 ∆x et on fait tendre ∆x vers

ww
TDEX-3.12
om
2 " 2 "
πR q x x − πR q x x +∆x − 2π R∆x h(Tx − T∞ )
2
=0
π R ∆x
q"x − qx" 2h
x x +∆x
− (Tx − T∞ ) = 0
∆x R
q"x − qx" ∂qx"
o.c
x x +∆x
lim =−
∆x → 0 ∆x ∂x
"
∂q 2 hx
− − (Tx − T∞ ) = 0
∂x R
 ∂T 
∂  −k x  2h
∂T
qx" = − k x ⇒ −  −
∂x
-pr
(Tx − T∞ ) = 0
∂x ∂x R
2
∂ Tx 2 h
− (Tx − T∞ ) = 0 (1)
∂x 2 kR
Faisons le changement de variable suivant
ari
∂θ ∂Tx ∂ 2θ ∂ 2T
θ = Tx − T∞ ⇒ = ⇒ `2 = 2x
∂x ∂x ∂x ∂x
2h
et posons m 2 =
kR
k
2
∂ Tx 2 h ∂ 2θ
− (Tx − T∞ ) = 0 ⇔ − m2θ = 0 (2)
∂x 2 kR ∂x 2
ab
L’équation différentielle (2) est une équation différentielle linéaire du second ordre à coefficients
constants) dont le polynôme caractéristique r 2 − m 2 = 0 admet deux racines +m et –m La solution de
(2) est donc :
θ = (Tx − T∞ ) = C1e m x + C2e − m x ( S1)
al3
Les constantes C1 et C2 sont obtenues en écrivant les conditions frontières :

à x = 0 Tx x=0
= Tp ⇒ θ x =0 = (Tp − T∞ ) = θ p = C1e m 0 + C2 e − m 0
⇒ C1 + C2 = θ p (a)
w.
à x = L Tx x=L
= Tp ⇒ θ x = L = (Tp − T∞ ) = θ p = C1e mL + C2 e − mL
ww
⇒ C1emL + C2 e − mL = θ p (b)
TDEX-3.13
La résolution du système des deux équations a) et b) donne les expressions des constantes C1 et C2 :
m
C1 + C2 = θ p ⇒ C1 = θ p − C2 ( S 2)
on remplace C1 dans b)
co
(θ p − C2 ) emL + C2 e− mL = θ p ⇒ C2 (e − mL − e mL ) = θ p (1 − e mL )
θ p (1 − emL )
⇒ C2 = ( S 3)
(e − mL − emL )
ro.
Application numérique :
2 o o
R=0.01 m, L=.05 m, h= 35 W/m .K T∞= 30 C; Tp= 90 C
k=385
k=44
ri-p C2= 53.638
C2= 58.890
C1= 6.361
C1= 0.109
au centre à z=L/2 T=66.9 C
au centre à z=L/2 T=35.1 C

o
o
ka
ab
al3
w.
ww
TDEX-3.14
m
Solutionnaire Ex-3.8 Ailette circulaire
Une ailette circulaire d'épaisseur W et de rayon exterieur R2 est installée sur un tube de rayon externe
R1 dont la surface est à la température T1. L'air environnant est à la température T4 et le coefficient de
convection à la surface de l'ailette est h.
co
a) Faites un bilan de chaleur sur un élement de volume approprié de l'ailette et obtenez l'équation
différentielle que doit satisfaire la température de cette ailette.
(remarque: on ne demande pas d'intégrer cette équation)
b) Quelles sont les différentes conditions frontières possibles ?
c) Si on appelle qC, les pertes de chaleur (en watts) au travers de cette ailette, donnez alors les
expressions du rendement et de l'efficacité de l'ailette.
ro.
a)Hypothèses : On néglige la variation de la température dans l’épaisseur (W) de l’ailette
∂
T=T(r) k=Constante Eg=0 régime permanent =0
Volume de contrôle
ri-p
2π r W ∆r
∂t
Ce qui rentre(conduction) – Ce qui sort (conduction+convection dessus et dessous)=0

ka
2π r Wqr" − 2π rWqr" − 2 [ 2π r ∆r h(Tr − T∞ ) ] = 0
r r +∆r
On divise par 2π W ∆r et on fait tendre ∆r vers

b
2π rWqr" − 2π rWqr" − 4π r ∆r h(Tr − T∞ )

r r +∆r
=0
3a
2π W ∆r
r qr" − rqr" 2 hr
r r +∆r
− (Tr − T∞ ) = 0
∆r W
r qr" − r qr" ∂rqr"
.al
r r +∆r
lim =−
∆r → 0 ∆r ∂r
∂rqr" 2 hr
− − (Tr − T∞ ) = 0
∂r W
w
ww
TDEX-3.15
m
  ∂T 
∂ ( r )  − k r   ∂ 2T   ∂Tr   ∂ ( r )   ∂ 2Tr 
"
∂T ∂rq  ∂r   ∂Tr 
qr" = − k r ⇒ − =−  r
= k ( r )  2r  + k    = kr  2 +k 
∂r ∂r ∂r  ∂r   ∂r   ∂r   ∂r   ∂r 
co
 ∂ 2Tr   ∂T  2 hr
kr  2  + k  r − (Tr − T∞ ) = 0
 ∂r   ∂r  W
 ∂ 2Tr  1  ∂Tr  2h
 2 +  − (Tr − T∞ ) = 0
 ∂r  r  ∂r  kW
ro.
b) Conditions frontières :
à r=R1 T=T1
à r=R2 une condition parmi les suivantes
1) densité de flux nulle qr"

r = R2
ri-p
= −k
∂Tr
∂r r = R2
=0
∂Tr
2) échange par convection : qr" = −k = h(Tr = R2 − T∞ )
r = R2 ∂r r = R2
3) si R2 très grand Tr=R2 = T∞

ka
c)
1) rendement = (perte réelle)/(perte si matériau infiniment conducteur)
rendement = (perte réelle)/(perte si ailette à T uniforme = T1)
ab
qc
η=
2 π ( R − R12 ) h (T1 − T∞ )
2
2
2) efficacité= (perte réelle)/(perte sans ailette)

al3
qc
ε=
2π R1Wh (T1 − T∞ )
w.
ww
TDEX-3.16
m
Solutionnaire Ex-3.9 Profil de température dans une ailette fixée entre deux parois
Hypothèses :
∂
T=T(x) k=Constante régime permanent =0
co
∂t
Volume de contrôle W t ∆x
Ce qui rentre(conduction) – Ce qui sort (conduction+convection dessus et dessous)=0
ro.
Wt q"x − Wtq"x − W ∆x hA (Tx − TA ) − W ∆x hB (Tx − TB ) = 0
x x +∆x
On divise par Wt ∆x et on fait tendre ∆x vers
Wt q"x − Wtq"x
x x +∆x
ri-p
− W ∆x hA (Tx − TA ) − W ∆x hB (Tx − TB )
=0
Wt ∆x
q"x − q"x hA h
x x +∆x
− (Tx − TA ) − B (Tx − TB ) = 0
ka
∆x t t
q"x − qx" ∂q"x
x x +∆x
lim =−
∆x → 0 ∆x ∂x
"
∂q hA + hB (h T + h T )
ab
x
− − Tx + A A B B = 0
∂x t t
 ∂T 
∂  −k x 
∂T ∂x  hA + hB (h T + h T )
qx" = −k x ⇒ −  − Tx + A A B B = 0
∂x ∂x t t
al3
2
∂ Tx  hA + hB (h T + h T ) 
2
− Tx − A A B B  = 0 (1)
∂x  kt kt 
Faisons le changement de variable suivant

w.
h + h   hATA + hBTB 
U = α Tx − β avec α =  A B  et β =  ( S1)
 kt   kt
∂ 2U ∂ 2Tx ∂ 2Tx 1 ∂ 2U
ww
∂U ∂T
=α x ⇒ = α ⇒ =
∂x ∂x ∂x`2 ∂x 2 ∂x 2 α ∂x`2
TDEX-3.17
L’équation (1) s’écrit donc :
om
1 ∂ 2U ∂ 2U
−U = 0 (2) ⇔ − αU = 0
α ∂x 2 ∂x 2
constants) dont le polynôme caractéristique r 2 − α = 0 admet deux racines α et − α . La solution de (2)
c
est donc :
ro.
αx αx
U = C1e + C2 e − ( S 2)
Les constantes C1 et C2 sont obtenues en écrivant les conditions frontières :
ri-p
à x = 0 Tx
⇒ C1 + C2 = U p
x=0
= Tp ⇒ U x = 0 = (α Tp − β ) = C1e
(a) avec
α0
+ C2 e − α0
U p = (α Tp − β ) ( S 3)
à x = L Tx x=L
= Tp ⇒ U x = L = (α Tp − β ) = U p = C1e αL
+ C2 e − αL
ka
⇒ C1e αL
+ C2 e − αL
=Up (b)
La résolution du système des deux équations a) et b) donne les expressions des constantes C1 et C2 :
C1 + C2 = U p ⇒ C1 = U p − C2 ( S 4)
ab
(U p − C2 ) e αL
+ C2 e − αL
=Up ⇒ C2 ( e − αL
−e αL
) = U p (1 − e αL
)
αL
U p (1 − e )
al3
⇒ C2 = − αL αL
( S 5)
(e −e )
w.
ww
TDEX-3.18
m
Solutionnaire : Travaux dirigés No 3
I) Profil dans un mur de béton ( 6 pts)
Un mur de béton (épaisseur b, largeur W, hauteur
co
L, conductivité k ) sépare une pièce à la température TINT ,
de l’extérieur dont la température est T EXT . Les coefficients
de transfert de chaleur sur les faces intérieure et extérieure
de ce mur sont respectivement hINT et hEXT
I.1 Faites un bilan de chaleur sur un volume de
ro.
contrôle pertinent et obtenez l’équation différentielle que doit
vérifier la température dans le mur (on négligera les effets
de bouts) (2 pts)
1) Hypothèses :
∂
T=T(x)
ri-p
k=Constante
2) Volume de contrôle
régime permanent
A ∆x
∂t
=0

ka
A q"x − Aqx" =0
x x +∆x
On divise par A∆x et on fait tendre ∆x vers 0
A q"x − Aqx" qx" − qx"
x x +∆x x x +∆x
= =0
A∆x ∆x
ab
q"x − qx" ∂q"x

x x +∆x
lim =− =0
∆x → 0 ∆x ∂x
On intègre une première fois :
∂q"x
al3
− = 0 ⇒ q"x = C1′
∂x
On utilise la loi de Fourier
∂T
q"x = − k = C1′
∂x
Et on intègre une deuxième fois :
w.
∂T C′  C′ 
=− 1 ⇒ Tx =  − 1  x + C2 ⇒ Tx = C1 x + C2
∂x k  k 
 
Les constantes d’intégration C’1 et C2 sont obtenues en écrivant les conditions frontières
ww
 ∂T 
CF1 : à x = 0 ⇒ q"x = hINT ( TINT − Tx =0 ) =  − k 
x=0
 ∂x  x = 0
TDEX-3.19
om
 ∂T 
CF2 : à x = b ⇒ q"x =  −k  = hEXT (Tx =b − TEXT )
x =b
 ∂x  x =b
 ∂T 
à x = 0 = hINT (TINT − Tx = 0 ) = hINT (TINT − [ C1 0 + C2 ]) = hINT (TINT − C2 ) =  − k  = − kC1
 ∂x  x =0
k
hINT (TINT − C2 ) = − kC1 ⇒ C2 = TINT + C1
o.c
hINT
à x = b ⇒ hEXT (Tx =b − TEXT ) = hEXT ([C1b + C2 ] − TEXT ) = − kC1
 k   k k 
hEXT  C1b + TINT + C1 − TEXT  = − kC1 ⇒ C1  b + +  = (TINT − TEXT )
 hINT   hINT hEXT 
-pr
C1 =
(TEXT − TINT )
 k k 
b + + 
 hINT hEXT 
II) Stockage d’air liquide dans un réservoir sphérique (6 pts)

ri
Un réservoir sphérique de rayon r1=1.5 m, contient de l’air liquide. Le réservoir est isolé par une épaisseur
de 0.05 m d'un matériau isolant de conductivité k=0.05 W/(m.K). On peut supposer que la paroi externe
ka
du réservoir reste à température constante T1=80 K. A la surface de l'isolant, il y a échange par
convection avec l'air ambiant à la température T∞=283 K et le coefficient de convection h est égal à 18
2
W/(m .K).
1 1 1 1
Rconduction =  −  Rconvection =
4π k  r1 r2  4π h r22
ab
sphèrecreuse sur une sphère
−2
Rconduction = 3.42 ×10 Rconvection = 1.84 ×10−3 RTOTALE = 3.606 ×10−2
sphèrecreuse sur une sphère
∆T∞−1 = RTOTALE q ⇒ q=
∆T
=
(T∞ − T1 ) = ( 283 − 80 ) = 5629.5 W
RTOTALE RTOTALE 3.606 ×10−2
al3
∆T∞− 2 = Rconvection q ⇒ (T∞ − T2 ) = 1.84 ×10−3 × 5629.5 = 10.35 = ( 283 − T2 )

T2 = 283 − 10.35 = 272.6 K
À 272.6 K l’eau gèle, il y aura donc condensation de la vapeur d’eau sur la surface du réservoir puis
congélation, le réservoir va devenir un gros glaçon : l’épaisseur est insuffisante !
w.
III) Température d’un fil électrique ( 3 pts)
Un fil électrique de 1 mm de diamètre est dénudé sur une longueur de 1 m. Calculer la

ww
température de ce fil (température qu’on supposera uniforme), sachant que la puissance générée par
o
effet joule est de 1.5 watt/m, que la température de l’air est de 20 C et que le coefficient de transfert
2
de chaleur à la surface du fil avec l’air environnant vaut 12. W/m .C.
TDEX-3.20
Bilan macroscopique sur le fil :
m
Ein − Eout + Egénéré = Eaccumulé mais Ein = 0 Eaccumulé = 0
Eout = π DLh (Tr = R − Tair ) et Egénéré = 1.5 W
co
1.5
π DLh (Tr = R − Tair ) = 1.5 ⇒ Tr = R = Tair + = 20 + 39.79 = 59.79o C
π × 0.001× 1× 12
IV) Quiz (2 pts) Cocher la bonne réponse:
ro.
• Pour la construction d’une ailette, il est préférable d’utiliser un matériau qui à un faible
conductivité. VRAI FAUX
• Dans une ailette, l’augmentation de la conductivité favorise la diminution du gradient de la

température. VRAI FAUX
-p
Dans un bilan de chaleur associé à un problème de conduction, si il y a un terme de génération de
chaleur alors le terme d’accumulation est non nul. VRAI FAUX
À l’interface entre deux solides, plus la résistance de contact est faible, plus la différence entre les
ari
•
températures de ces deux surfaces est grande. VRAI FAUX
V) Ailette en forme de cylindre creux (8

pts)
Une ailette, de longueur L, a la forme d'un
bk
cylindre creux. Sa base est fixée sur une paroi à la

température Tp. La face intérieure du cylindre est
refroidie par de l’air à la température T1 avec un
coefficient de convection h1 alors que la face
extérieure est refroidie avec un air à la
3a
température
T2 avec un coefficient de convection h2.
On supposera que le flux de chaleur axial à x=L

sur l’extrémité de l’ailette est négligeable.
al
a) Faites un bilan de chaleur sur un élément de

volume approprié de l'ailette et obtenez l'équation différentielle que doit satisfaire la température de cette
ailette. (4.5 pts)
b) Quelles sont les conditions frontières ? (1 pt)

w.
c)Obtenez l’expression du profil de température. (2.5 pts)

∂
T=T(x) k=Constante Eg=0 régime permanent =0
∂t
ww
Volume de contrôle π ( r2 2 − r12 ) ∆x = s ∆x

TDEX-3.21
m
Ce qui rentre(conduction) – Ce qui sort (conduction+convection à la surface de la tige)=0
co
s qx" − sqx" − 2π r1∆x h1 (Tx − T1 ) − 2π r2 ∆x h2 (Tx − T2 ) = 0
x x +∆x
On divise par s∆x et on fait tendre ∆x vers
ro.
s qx" − sqx" − 2π r1∆x h1 (Tx − T1 ) − 2π r2 ∆x h2 (Tx − T2 )
x x +∆x
=0
s ∆x
qx" − qx" 2π r1 h1 (Tx − T1 ) + 2π r2 h2 (Tx − T2 ) qx" − qx" ∂q"
x x +∆x x +∆x
− =0 lim x =− x
∆x s δ x →0 ∆x ∂x
"
qx = −k
∂ 2Tx
∂Tx
∂x
⇒−
∂qx
∂x
"
=−
 ∂T 
∂  −k x 

 2π r h + 2π r2 h2 
∂x
∂x  ∂ 2Tx
=k 2 ⇒ 2 −
∂x ∂x
 2π r1 h1T1 + 2π r2 h2T2 
ri-p
∂ 2Tx 2π r1 h1 (Tx − T1 ) + 2π r2 h2 (Tx − T2 )
ks
=0
2
− Tx  1 1  +  = 0 (1)
∂x  ks   ks 
ka
 2π r h + 2π r2 h2   2π r1 h1T1 + 2π r2 h2T2  2
posons U = Tx  1 1 −  = α Tx − β


ks

 ks 
α2 β
2 2 2 2
∂U ∂T ∂T 1 ∂U
⇒ = α 2 2x ⇒ 2x = 2 2
b
2
∂x ∂x ∂x α ∂x
2
1 ∂U ∂ 2U
(1) ⇔ − U = 0 ⇔ − α 2U = 0 (2)
3a
α 2 ∂x 2 ∂x 2
constants) dont le polynôme caractéristique r 2 − α 2 = 0 admet deux racines +α et −α La solution de
(2) est donc :
U = α 2Tx − β = C1eα x + C2 e −α x ( S1)
.al
Les constantes C1 et C2 sont obtenues en écrivant les conditions frontières:

à x = 0 Tx x=0
= Tp ⇒ U x = 0 = α 2Tx =0 − β = α 2Tp − β = U p = C1eα 0 + C2 e −α 0
⇒ C1 + C2 = U p (a)
w
∂Tx ∂T ∂U
à x = L le flux est négligeable qx" = − k =0⇒ x =0⇒ =0
∂x ∂x ∂x
ww
∂U
⇒ = α C1eα L − α C2 e −α L = 0 (b)
∂x x=L
TDEX-3.22
On résout le système de 2 équations a) et b) à 2 inconnues C1 et C2
m
C1 + C2 = U p ⇒ C1 = U p − C2
o
α (U p − C2 ) eα L − α C2 e−α L = 0 ⇒ C2 ( eα L + e−α L ) = U p eα L
o.c
U p eα L
⇒ C2 = avec U p = α 2Tp − β
(eαL
+e −α L
)
ri -pr
ka
ab
al3
w.
ww
TDEX-4.1
m
Conduction stationnaire en deux dimensions
Exercices
co
Ex-4.1 Convention des flux entrants
En utilisant la convention des flux entrants (page 4.6), retrouver par bilan les équations pour les 3 géométies
ci-dessous (on appellera L la profondeur des pièces et on a )x=)y)
ro.
a) Convection dans un coin intérieur:
h∆x h∆x
ri-p
2( T m-1,n + T m,n+1 )+( T m+1,n + T m,n-1 )+ 2
k
T ∞ - 2(3+
k
)T m,n = 0
ka
b) Convection sur la surface:
h∆x h∆x
(2T m-1,n + T m,n+1 + T m,n-1 )+ 2 T ∞ - 2( + 2)T m,n = 0
k k
ab
al3
c) Convection sur un coin extérieur:
h∆x h∆x
( T m-1,n + T m,n-1 )+ 2 T ∞ - 2( +1)T m,n = 0
k k
w.
ww
TDEX-4.2
m
Solutionnaire 4.1a :
Autour du nœud Tm,n, le volume de contrôle
.co
pertinent est :
3
∆x ∆x L
4
On identifie les 6 faces au travers des quelles
il y aura un flux entrant dans ce volume et
dont les expressions sont les suivantes :
o
 ∆x  (Tm,n −1 − Tm, n )
q1 =  L k -pr q2 = ( ∆x L ) k
(T m −1, n − Tm ,n )
ari
 2  ∆x ∆x
q3 = ( ∆x L ) k
(Tm,n+1 − Tm,n )  ∆x  (Tm +1,n − Tm ,n )
q4 =  L k
∆x  2  ∆x
 ∆x   ∆x 
q5 =  L  h (T∞ − Tm ,n ) q6 =  L  h (T∞ − Tm ,n )
bk
 2   2 
6
Le bilan s’écrit : ∑q
i =1
i =0
3a
L L
k (Tm ,n −1 − Tm ,n ) + L k (Tm −1,n − Tm ,n ) + L k (Tm ,n +1 − Tm ,n ) + k (Tm +1, n − Tm ,n )
2 2
 ∆x   ∆x 
+ L  h (T∞ − Tm ,n ) +  L  h (T∞ − Tm ,n ) = 0
.al
 2   2 
k (Tm ,n −1 − Tm ,n ) + 2 k (Tm −1, n − Tm ,n ) + 2 k (Tm ,n +1 − Tm ,n ) + k (Tm +1, n − Tm ,n ) + 2∆x h (T∞ − Tm, n ) = 0

w
2 ∆x h  ∆x h 
⇒ 2 (Tm −1,n + Tm ,n +1 ) + (Tm ,n −1 + Tm +1, n ) + T∞ − 2  + 3  Tm ,n = 0
k  k 
ww
TDEX-4.3
m
Solutionnaire 4.1b :
Autour du nœud Tm,n, le volume de contrôle pertinent est :
.co
∆x ∆x L
2
On identifie les 4 faces au travers des quelles il y aura un
flux entrant dans ce volume et dont les expressions sont
les suivantes :
q1 = ( ∆x L ) h (T∞ − Tm, n )
o
 ∆x  (Tm, n −1 − Tm ,n )
q2 = 
-pr
L k
 2  ∆x
q3 = ( ∆x L ) k
(Tm−1,n − Tm,n )
∆x
 ∆x  (Tm ,n +1 − Tm ,n )
q4 =  L k
ari
 2  ∆x
4
i =1
i =0
bk
∆x  (Tm ,n −1 − Tm ,n ) (Tm−1,n − Tm,n ) +  ∆x L  k (Tm,n+1 − Tm,n ) = 0

( ∆x L ) h (T∞ − Tm,n ) +  L k + ( ∆x L ) k  
 2  ∆x ∆x  2  ∆x
1 1
∆x h (T∞ − Tm ,n ) + k (Tm ,n −1 − Tm ,n ) + k (Tm −1,n − Tm ,n ) + k (Tm ,n +1 − Tm ,n ) = 0
3a
2 2
2∆x h
k
(T∞ − Tm,n ) + (Tm,n−1 − Tm,n ) + 2 (Tm−1,n − Tm,n ) + (Tm,n+1 − Tm,n ) = 0
( 2Tm−1,n + Tm,n−1 + Tm,n+1 ) + 2∆kx h T∞ − 2  ∆xk h + 2  Tm,n = 0
.al
 
w
ww
TDEX-4.4
m
Solutionnaire 4.1c :
Autour du nœud Tm,n, le volume de contrôle pertinent est :
co
∆x ∆x L
4
On identifie les 4 faces au travers des quelles il y
aura un flux entrant dans ce volume et dont les
expressions sont les suivantes :
ro.
 ∆x L 
q1 =   h (T∞ − Tm ,n )
 2 
 ∆x  (Tm, n −1 − Tm ,n )
q2 =  L k
 2 
 2 
 ∆x L 
k
ri-p ∆x
 ∆x L  (Tm −1,n − Tm ,n )
q3 = 
∆x
q4 =   h (T∞ − Tm ,n )
 2 
ka
4
i =1
i =0
 ∆x  (Tm ,n −1 − Tm ,n )  ∆x L  (Tm −1, n − Tm ,n )  ∆x L 

ab
 ∆x L 
  h (T∞ − Tm ,n ) +  L k + k +  h (T∞ − Tm ,n )
 2   2  ∆x  2  ∆x  2 
2∆xh
k
(T∞ − Tm,n ) + (Tm,n−1 − Tm,n ) + (Tm−1,n − Tm,n ) = 0
al3
(Tm,n−1 + Tm−1,n ) + 2∆kxh T∞ − 2  ∆xh 

+ 1 Tm ,n = 0
 k 
w.
ww
TDEX-4.5
m
Ex-4.2 : Exemple de résolution d==un problème de conduction 2D par différences finies
(exemple de résolution par la méthode de Gauss-Seidel avec Excel)
La chambre d=un four industriel comporte une colonne de brique de conductivité k=1 w/m.K et de section carré
co
de 1 m de côté. Trois faces de cette colonne sont soumises à une température de 500 K alors que la dernière
est soumise à la convection avec l=air environnant à la température de 300 K et le coefficient h vaut 10
2
W/m .K (on appellera H la hauteur de cette colonne)
En utilisant un maillage de 0.25 m, déterminer le profil de température dans la colonne.
ro.
a) On remarque que le problème présente une symétrie et il suffit de le résoudre sur la moitié de la
colonne. Ci-dessous, on définit donc 8 nœuds et le plan qui passe par les points 2-8 correspond au
plan de symétrie.
b) Les nœuds sur les 3 faces de la colonne sont à la température TS=500 K (et ne seront pas identifiés
par un numéro, on utilisera TS).
ri-p
h=10W/m .K
K
2
Tair = 300
7 8 7
ka
5 6 5
TS =500 K TS =500 K
b
3 4 3
3a
1 2 1
TS =500 K
.al
Rappel :dans un bilan de chaleur, si une face délimitant le volume de contrôle correspond à un plan
de symétrie alors le flux au travers de ce plan est nul.
w
c) Bilan sur les nœuds 1 à 8. Pour ce faire, on utilise un autre maillage décalé d’une demi-maille
permettant de définir les volumes de contrôle pertinent autour de chaque nœud (voir pointillé ci-
ww
dessous)
TDEX-4.6
m
- Pour les nœuds 1-3-5, le volume est :
∆x ∆x H
- Pour les nœuds 2-4-6-7, le volume est :
co
∆x
∆x H
2
- Pour le nœud 8, le volume est :
∆x ∆x
ro.
H
2 2
Pour chaque nœud, on fait le bilan de chaleur en

utilisant la convention des flux entrants.
ri-p
Nous allons détailler le bilan uniquement pour un
nœud, le no 4. Au travers des faces a,b,c il y a apport
de chaleur par conduction alors qu’au travers de la
face d le flux est nul puisque il s’agit d’un plan de
symétrie.
ka
 ∆x  T2 − T4
qa =  Hk
 2  ∆x
T −T
qb = ( ∆x H ) k 3 4
∆x
b
 ∆x  T6 − T4
qc =  Hk
 2  ∆x
3a
qd = 0
qa + qb + qc = 0
.al
 ∆x  T −T T − T  ∆x  T −T
⇒ H  k 2 4 + ( ∆x H ) k 3 4 + 
w
H  k 6 4 =0
 2  ∆x ∆x  2  ∆x
T2 − T4 T −T
+ T3 − T4 + 6 4 = 0 ⇒ T2 − T4 + 2T3 − 2T4 + T6 − T4 = 0
ww
2 2
pour le noeud 4 : T2 + 2 T3 + T6 − 4 T4 = 0
(faites par vous-même les bilans sur tous les autres nœuds)
TDEX-4.7
Les équations pour les 8 nœuds sont donc les suivantes:
m
pour le noeud 1: T2 + T3 + 1000 − 4 T1 = 0
pour le noeud 3 : T1 + T4 + T5 + 500 − 4 T3 = 0
pour le noeud 5 : T3 + T6 + T7 + 500 − 4 T5 = 0
o
pour le noeud 2 : 2T1 + T4 + 500 − 4 T2 = 0
pour le noeud 4 : T2 + 2T3 + T6 − 4 T4 = 0
o.c
pour le noeud 6 : T4 + 2T5 + T8 − 4 T6 = 0
pour le noeud 7 : 2T5 + T8 + 2000 − 9 T7 = 0
pour le noeud 8 : 2T6 + 2T7 + 1500 − 9 T8 = 0
Ces équations nous permettent d=exprimer les températures en chaque nœud :
-pr
T2 + T3 + 1000 2T1 + T4 + 500 T1 + T4 + T5 + 500 T2 + 2T3 + T6
T1 = T2 = T3 = T4 =
4 4 4 4
T + T + T + 500 T + 2T5 + T8 2T + T + 2000 2T + 2T7 + 1500
T5 = 3 6 7 T6 = 4 T7 = 5 8 T8 = 6
4 4 9 9
ri
d) Pour résoudre ce système avec la méthode de Gauss-Seidel, nous allons utiliser Excel et son outil
Solveur (voir ci-dessous et récupérer la feuille Excel sur le site du cours).
ka
La procédure est la suivante :
- dans la colonne B, écrire le nom de

chacun des noeuds (TC1, TC2 …TC8)
- dans la colonne C, nous rentrons un

ab
estimé des températures des 8 nœuds (par

exemple 450 K)
- pour chaque cellule de C, nous définissons

un nom de variable (TC1, TC2 …TC8)
commandes Excel : se placer sur une
al3
cellule puis Insertion-Nom-Definir.
ps par défaut le nom apparaissant dans la

cellule à gauche est retenu.
-dans la colonne C, pour chaque nœud

écrire (sous forme litérale) l’équation
w.
correspondante.
dans la colonne E nous calculons le carré

de la différence entre la valeur estimée
(colonne C) et la valeur calculée (colonne
ww
D). La somme de ces carrés est faite dans

la cellule E17.
TDEX-4.8
Se placer sur la cellule E17 et choisir le Solveur dans la rubrique Outils.
m
Si cette fonctionalité n’apparait pas
choisir alors l’onglet Macros
o
complémentaires et installer le Solveur
o.c
-pr
Se placer sur la cellule E17, choisir le
Solveur et dans la fenêtre Paramètres du
solveur indiquer au Solveur que les
valeurs estimées des températures (C8-
C15) sont les cellules à faire varier afin
de minimiser la somme de la cellule E17.
ri
Le Solveur va faire des itérations sur les
valeurs des températures estimées
jusqu’à avoir T estimée= Tcalculée.
ka
Cliquer Résoudre
ab
al3
w.
ww
TDEX-4.9
Après optimisation, on trouve:
m
co
ro.
ri-p
ka
d) Sur le volume en pointillé faites un bilan macroscopique (pour 1 m de

hauteur de colonne) et vérifier que la chaleur entrant par conduction est
égale à celle sortant par convection:
ab
- échange par convection avec l=air:

∆x
out
qconvection = 1 ∆x h (T7 − Tair ) + 1 h (T8 − Tair )
2
(réponse 191.3 W pour 1m)
al3
-apport d’énergie par conduction

w.
in
=k
∆x  TS − T2  (TS − T1 ) + k ∆x (TS − T1 ) + k ∆x (TS − T3 ) + k ∆x (TS − T5 ) + k ∆x  TS − T7 
qconduction   + k ∆x  
ww
2  ∆x  ∆x ∆x ∆x ∆x 2  ∆x 
= k [5TS − 2T1 − T3 − T5 − 0.5T2 − 0.5T7 ] (réponse :191.31 W / m)
ww
w.
al3
ab
kari
-pro.
co
m
TDEX- 5.1
m
Ex-5.1 (S) Conduction en transitoire
co
Une plaque métallique rectangulaire (largeur W, longueur L, épaisseur
b) de conductivité k, est initialement à une température uniforme To.
Soudain ses deux faces sont mises en contact avec deux fluides ayant
des températures T1 et T2 et les coefficients de transfert de chaleur sur
les faces sont respectivement h1 et h2.
ro.
Quelle est l'expression de la température au centre de la plaque en
fonction du temps?
On supposera: (h1 b/k) <<1; (h2 b/k) <<1
Réponse :
( h 1T 1 + h 2T 2 ) - ( h 1 + h 2 )T ( h1+h 2 )t
= e - bρ C P
ri-p
( h 1T 1 + h 2T 2 ) - ( h 1 + h 2 )T 0
Ex-5.2 (S) Régime transitoire en cylindrique (abaque de Heissler)

ka
O
Un cylindre de 2 m de long et de 0.20 m de diamètre, initialement à une température uniforme de 400 C, est
O -5 2
plongé dans un bain d=eau à 50 C. Sachant que la diffusivité thermique du cylindre, α, vaut 10 m /s la
conductivité thermique k vaut 40 W/m.K et que le coefficient de transfert de chaleur à la surface du cylindre h
2
vaut 200 W/m .K , on demande de calculer après 20 minutes d=immersion:
- la température au centre du cylindre
ab
- la température à la surface du cyclindre

- la fraction de la chaleur potentiellement échangeable (Q/Q0) aisni que la chaleur
échangée avec le bain (exprimée en kW.hr)
O O
On négligera les effets de bouts (Réponses: 184.8 C, 156.6 C, 0.656,16.03 kWhr)
al3
w.
ww
TDEX- 5.2
m
Chauffage par induction, régime transitoire avec résistance interne négligeable
co
Ex-5.3 (S)
Une bille métallique de rayon R initialement à une température TO est indroduite dans un dispositif de
'''
chauffage par induction (comme un four micro-onde). La bille est le siège d'une génération de chaleur q
3
(W/m ) et on suppose qu'à la surface de la bille, il y a échange par convection (coefficient h) avec l'air
environnant qui reste à une température constante, TA. On souhaite étudier la variation dans le temps de la
ro.
température de la bille.
a) A quelle condition peut-on négliger le gradient radial de température dans la bille ?
b) On suppose que le gradient radial de température est négligeable, faites un bilan de chaleur et
obtenez l'équation différentielle que doit vérifier la température de la bille ?
c) Obtenez l'expression de la variation dans le temps de la température ?
''' '''
Réponse: q -(3h/R)(T-TA)= C exp(-3ht/RρcP) avec C=q -(3h/R)(TO-TA)
d) Si on attend suffisamment longtemps, quelle sera la température d'équilibre de la bille ? On peut
Réponse:
'''
TEQUILIBRE-TA=q R/(3h) ri-p
répondre à la question sans utiliser la réponse c)
Ex-5.4 (S) Régime transitoire géométrie semi-infinie

-5 2
Une pièce de bronze très épaisse (diffusivité thermique α=0.86 10 m /s et conductivité k=26 W/m.K) est
ka
O O
initialement à une température de 250 C. Soudainement sa surface est exposée à un fluide à 25 C et le
2
coefficient d'échange est de 150 W/m .K. Déterminer la température dans le bronze à 5 cm sous la surface
O
après 10 min d'exposition.(utiliser une résolution graphique) (Réponse: 214 C)
b
Ex-5.5 (S) Profondeur pour enterrer les conduites d’eau (géométrie semi-infinie)
3a
Les conduites d’eau potable enterrées dans le sol doivent être à l’abri du gel. Déterminer la profondeur
d’enfouissement sécuritaire correspondant aux conditions extrêmes suivantes:
o
- le 31 octobre la température du sol est uniforme et égale à 10 C
er o
- le 1 novembre, il y a refroidissement à -15 C, et cette température se maintient pendant 120
2
jours avec un fort vent qui impose un coefficient de transfert de chaleur à la surface du sol de 35 W/m .K
.al
er
- et le 1 mars, il fait beau !
3
Propriétés du sol: k=0.52 W/m/K; D =2050 kg/m ; CP=1840 J/kg/K (réponse : 1.36m)
w
ww
TDEX- 5.3
m
Exercices pour pratiquer les abaques de Heisler. Faites les résolutions graphiques et
comparez les solutions obtenues en utilisant les logiciels disponibles sur le site du cours
co
Ex-5.6 (S) Géométrie plane:
Dans un four, une des parois ayant une épaisseur de 50 cm, est initialement à une température uniforme de
0
60 C. La face extérieure de cette paroi est parfaitement isolée alors que la face intérieure est soudainement
0 2
exposée à une température de 900 C et le coefficient de transfert de chaleur à la surface vaut h=25 W/m /K.
3
Les propriétés de la paroi étant les suivantes: cP=837 J/kg/K; ρ=500 kg/m ; k=1.25 W/m/K, déterminer le
ro.
0
temps requis pour que la température de la face externe de la paroi atteigne 600 C. A cet instant, quelles sont
: a) la température 10 cm en dessous de la face interne b) l'énergie totale reçue par la paroi depuis le début
2 8
du chauffage (pour 1 m de surface de mur) [réponses : 51750-775 -Q/Qo=0.75 et Q=1.32 10 ]
Ex-5.7 (S) Géométrie sphérique:
ri-p
bille est soumise à un traitement thermique en deux étapes.
a)
0
0
Une bille métallique de 5 mm de rayon est initialement à une température uniforme de 400 C. Cette
Dans la première, on la refroidit dans l'air à 20 C pendant un temps tair nécessaire pour avoir une
0
température au centre de la sphère égale à 335 C. Pendant cette étape, le coefficient à la surface vaut
2
hair=10 W/m .K L'hypothèse d'une résistance thermique interne négligeable est-elle applicable à cette étape ?
ka
Calculer le temps tair ?[94 s]
0
b) Dans une deuxième étape, la bille est alors introduite dans un bain d'eau à la température de 20 C.
2
Ce bain est fortement mélangé et le coefficient heau vaut alors 6000 W/m .K Quel est le temps requis pour
0
avoir une température au centre de 50 C ? [3 s]
ab
3
Données: k=20 W/m/K, : cP=1000 J/kg/K; ρ=3000 kg/m
Ex-5.8 (S) Géométrie mixte:

al3
Un cylindre de 10 cm de diamètre et 8 cm d=épaisseur est initialement à une température de 700 K. Pour le

refroidir, on le plonge dans un bain d=huile à 300 K et le coefficient de transfert de chaleur à la surface vaut
2
400 W/m K
Après 10 minutes d=immersion, quelles sont les températures minimale et maximale dans le cylindre ?
3
Données: ρ=8050 kg/m , k=18.6 W/.mK, CP=536 J/kg.K, [Réponses: 316 et 346 K]
w.
ww
TDEX- 5.4
m
Ex-5.9 (S) Résolution d’un régime transitoire avec l’ hypothèse du régime quasi-permanent :
détermination du temps de solidification d'une goutte d'eau
co
On souhaite prédire le temps de solidification,tf, d'une gouttelette sphérique de rayon R de
conductivité k, à la surface de laquelle la chaleur est échangée avec de l'air froid à une température TA, le
coefficient d'échange à la surface de la goutte est constant et vaut h. On supposera que:
- la goutte est initialement à sa température de fusion Tf
- la chaleur sensible de la partie solidifiée est négligeable devant la chaleur de fusion
- les densités des phases liquide et solide sont égales
ro.
Initialement, tout est liquide et à la fin du processus tout est solide. La solidification commence à se faire à
la périphérie puis le front de solidification pénètre à l'intérieur de la goutte. Considérez le schéma de la
goutte à un instant quelconque du processus de refroidissement. Dans ce problème, on va considérer que
la variation du rayon de solidification, Rf, est lente. Pendant un bref intervalle de temps, on va donc
résoudre en stationnaire, le profil de température dans la couche déjà solidifiée.
R
ri-p
glace
ka
Rf Échange avec l’air à TA
et coefficient h
Eau à Tf
ab
al3
a) Par bilan d'énergie en régime stationnaire sur la partie solidifiée, obtenez le profil de température:
C1
T (r ) = + C2
r
w.
b) En utilisant les conditions frontières pertinentes, trouvez les constantes C1 et C2
T f − TA C1
C1 = et C 2 = T f −
1 1 k Rf
− +
ww
R f R hR 2
TDEX- 5.5
m
b) Démontrez que les pertes de chaleur avec l'air qui arrivent par conduction à r=R s’exprime par:
4 π k (T f − T A )
Qf =
1 1 k
− +
co
R f R hR 2
c) En remarquant que durant un instant dt, ces pertes sont associées à l’énergie dégagée lors de la
ro.
solidification d’une petite couche d’eau d’épaisseur dRf , trouver l’équation différentielle reliant les pertes de
chaleur à la variation dans le temps de l'interface liquide-solide (Rf). On appellera ∆Hf la chaleur latente
de solidification. (attention dRf est négatif)
d) Combinez b et c et après intégration, obtenez le temps de solidification de la goutte:
ri-p
Application numérique
2
tf =
ρ ∆H f
k (T f − T A )
3
 R 2 kR 
 +
 6 3h 

0
R=0.01 m, h=300 W/m C, k=2.2 W/m/C, ∆Hf=333 KJ/kg, ρ=1000 kg/m , (Tf-TA)=10 C
(622s)
ka
ab
al3
w.
ww
TDEX- 5.6
om
Solutionnaire Ex-5.1 Conduction en transitoire
Une plaque métallique rectangulaire (largeur W, longueur L, épaisseur b) de conductivité k, est initialement à
une température uniforme To. Soudain ses deux faces sont mises en contact avec deux fluides ayant des
températures T1 et T2 et les coefficients de transfert de chaleur sur les faces sont respectivement h1 et h2.
Quelle est l'expression de la température au centre de la plaque en fonction du temps?

On supposera: (h1 b/k) <<1; (h2 b/k) <<1
o.c
1) Hypothèses :
puisque (h1 b/k) <<1 et (h2 b/k) <<1 alors on peut supposer que T=T(t)
k=Constante
2) Volume de contrôle WLb

ri -pr
– Ce qui sort = Ce qui s’accumule
−WL h1 (T − T1 ) − WL h2 (T − T2 ) = WLb ρ CP
∂T
∂t
ka
On divise par WL
∂T
− ( h1 + h2 ) T + ( h1T1 + h2T2 ) = b ρ CP [1]
∂t
posons : u = − ( h1 + h2 ) T + ( h1T1 + h2T2 )
ab
(h1T1 + h2T2 ) − u ∂T 1 −∂u

⇒T = ⇒ =
( h1 + h2 ) ∂t ( h1 + h2 ) ∂t
1 −∂u
[1] ⇒ u = bρ CP
al3
( h1 + h2 ) ∂t
⇒
∂u (h + h )
=− 1 2 u ⇒
∂u (h + h )
= − 1 2 ∂t
∂t bρ CP u b ρ cP
u
∂u
t
( h1 + h2 ) u ( h1 + h2 )
∫u u t =∫0 − bρCP ∂t ⇒ Ln u0 = − bρCP ( t − 0 )
=
w.
( h1 + h2 )
u −(h1 + h2 ) T + (h1T1 + h2T2 ) −
b ρ cP
t
= =e
u0 −(h1 + h2 ) T0 + (h1T1 + h2T2 )
ww
TDEX- 5.7
m
Solutionnaire Ex-5.2 Régime transitoire en cylindrique (abaque de Heissler)
αt 10−5 × 1200
t = 20 × 60 = 1200 s Fo = = = 1.2
ro2 0.102
co
h ro 200 × 0.10
Biot = = = 0.5 ⇒ Biot −1 = 2.
k 40
Sur la figure page 5.13
To − T∞ T − 50
ro.
pour Fo = 1.2 et Biot −1 = 2. ⇒ θ o* = ≈ 0.38 = o ⇒ To = 0.38 × 350 + 50 = 183o C
Ti − T∞ 400 − 50
Sur la figure 5.12 page 5.14 :
r θ Tr = ro − T∞
pour = 1 et Biot −1 = 2 ⇒ = ≈ 0.78 ⇒ Tr = ro ≈ 0.78 × (183 − 50 ) + 50 = 153.7 o C
ro θo To − T∞
ri-p
(avec le logiciel, les valeurs exactes sont :To= 184.8 et TS=156.6)
Sur la figure 5.13 page 5.14 :
Q
pour Bi 2 Fo = 0.52 × 1.2 = 0.30 et Bi = 0.5 on lit ≈ 0.6
Qo
ka
k k 40
α= ⇒ ρ cP = = −5 = 4.106
ρ cP α 10
⇒ Q = 0.60 × Qo = 0.60 × ρ cP × (π ro2 L ) × (Ti − T∞ )
ab
Q = 0.60 × 4.106 × (π × 0.12 × 2 ) × ( 400 − 50 ) = 52778756.6 joules

1 kWh = 1000 × 3600 = 3.6 106 joules
52778756.6
Q= = 14.7 kWh
3.6 106
al3
O O Q
p.s. : Les réponses exactes avec le logiciel sont : :TO=184.8 C, TS=156.6 C, 0.656, Q =16.03 kWhr
Qo
w.
ww
TDEX- 5.8
om
Solutionnaire Ex-5.3 Chauffage par induction, régime transitoire
a) le gradient radial de température dans la bille est négligeable si Biot< 0.1

b) Hypothèse : Biot< 0.1 ⇒ T = T (t ) uniquement
4
Volume contrôle: tout l’objet = π R3
3
Bilan :
c
Ein − Eout + E généré = Eaccumulé
ro.
4π R 3 4π R 3 ∂T (t )
0 − 4π R 2 h ( T (t ) − T∞ ) + q′′′ = ρ CP
3 3 ∂t
3h ∂T (t )
0 − (T (t ) − T∞ ) + q′′′ = ρ CP
R ∂t
∂T (t )
∂t
=−
3h
ρ cP R
c)
ri-p
(T (t ) − T∞ ) +
q′′′
ρ CP
3h
=−
3h
ρ CP R
 q′′′ 3h
 q′′′
T (t ) +  +
3h
 ρ CP ρ CP R 
 ∂u

T∞ 
3h ∂T (t ) ∂T (t ) ρ CP R ∂u
posons : u = − T (t ) +  + T∞  ⇒ =− ⇒ =−
ρ CP R  ρ CP ρ CP R  ∂t ρ CP R ∂t ∂t 3h ∂t
ka
∂T (t ) ρ CP R ∂u ∂u 3h
⇒ = − =u ⇒ =− ∂t
∂t 3h ∂t u ρ CP R
 u 
ut t
∂u 3h 3h
∫u u ρCP R ∫0 ∂t ⇒ Ln  ut =0  = − ρCP R t
= −
ab
t =0
3h  q′′′ 3h  3h  3h 
− T (t ) +  + T∞  3h − T (t ) +  q′′′ + T∞ 
u ρ CP R  ρ CP ρ CP R  − t
R  R 
⇒ = =e ρ CP R
=
ut = 0 3h  q′′′ 3h  3h  3h 
− T0 +  + T∞  − T0 +  q′′′ + T∞ 
ρ CP R  ρ C P ρ CP R  R  R 
al3
d)
Méthode 1 : la température d’équilibre s’obtient à t=∞
q′′′ 3h
+ T
3h  q′′′ 3h  ρ CP ρ CP R ∞ q′′′R
ut =∞ = 0 ⇒ − T + + T = 0 ⇒ Tt =∞ = = + T∞
ρ CP R t =∞  ρ CP ρ CP R ∞  3h 3h
w.
ρ CP R
Méthode 2 : à l’équilibre le bilan de chaleur donne : − Eout + Egénéré = 0 ⇒ E généré = Eout
4π R 3 Rq′′′
q′′′ = 4π R 2 h (Tt =∞ − T∞ ) ⇒ Tt =∞ = T∞ +
3 3h
ww
TDEX- 5.9
m
Solutionnaire Ex-5.4 Régime transitoire géométrie semi-infinie
La pièce de bronze est très épaisse et on suppose que loin de la surface c’est l’infini et on utilise le
graphique p 5.22.
x 0.05
co
x = 0.05 m t = 10 × 60 = 600 s η = = = 0.348
2 αt 2 0.86 10 −5 × 600
−5
h α t 150 0.86 10 × 600
= = 0.414
k 26
ro.
(T − Ti ) ≈ 0.18 ⇒ T ≈ 0.18 T − T + T = 0.18 × 25 − 250 + 250 = 209.5o C
sur le graphique, on trouve : ( ∞ i) i ( )
(T∞ − Ti )
Le calcul exact avec l’équation [14] p. 5.20 donne:
(T − Ti ) = 0.1774 ⇒ T = 0.1774 × 25 − 250 + 250 = 210.1o C
ri-p
(T∞ − Ti )
Solutionnaire
( )
Ex-5.5 Profondeur pour enterrer les conduites d’eau

Les conduites d’eau potable enterrées dans le sol doivent être à l’abri du gel. Déterminer la profondeur
d’enfouissement sécuritaire correspondant aux conditions extrêmes suivantes:
o
- le 31 octobre la température du sol est uniforme et égale à 10 C
ka
er o
- le 1 novembre, il y a refroidissement à -15 C, et cette température se maintient pendant 120
2
jours avec un fort vent qui impose un coefficient de transfert de chaleur à la surface du sol de 35 W/m .K
er
- et le 1 mars, il fait beau !
3
Propriétés du sol: k=0.52 W/m/K; D=2050 kg/m ; CP=1840 J/kg/K
ab
Réponse :
Il s’agit là d’un cas d’échange à la surface d’un solide semi infini (figure p 5.22). On supposera que la
O er
al3
température initiale du sol est uniforme et égale à 10 C.Au temps t=0 (le 1 novembre), il y a échange à
2 O
la surface avec un coefficient h=35 W/m .K avec de l’air à -15 C et ces conditions durent pendant 120
jours. Pour éviter que l’eau de la conduite gèle, il faut qu’après 120 jours la température à la surface de la
O
conduite (enterrée à une profondeur X que l’on cherche) soit supérieure ou égale à 0 C
w.
ww
TDEX- 5.10
m
Tair=- 15 C Échange avec la surface du sol, h=35 W/m2. oC
o
.co
Air
X=0 Sol
o
-pr
X, profondeur à déterminer
ari
Pour ce problème, on a donc:
Ti = 10 0C T∞ = − 15 0C T ( x, t = 120 j ) = 00 C
T ( x, t ) − Ti (0) − 10 10
= = = 0.40
T∞ − Ti (−15) − 10 25
k
k 0.52
α= = = 1.38 10−7 t = 120 24 3600 = 1.037 107
ρ CP 2050 1840
ab
h α t 35 (1.38 10 )(1.037 10 )
−7 7
= = 80.47
k 0.52
h αt
al3
Sur le graphique, la courbe correspondant à une valeur de = 80.47 est la même que celle
k
h αt
pour = ∞ . Pour une valeur de l’ordonnée de 0.4, l’intersection avec la courbe ∞ se fait à une
k
valeur de l’abscisse
x
≈ 0.58
w.
2 αt
On trouve donc la valeur de

ww
x = 0.58 2 α t = 1.14 1.38 10 −7 1.037 107 = 1.38 m

TDEX- 5.11
m
Solutionnaire Ex-5.6 Géométrie plane
co
Pour résoudre ce problème, il faut remarquer
que la paroi externe du four est parfaitement
islolée alors que la paroi interne est soumise à
de la convection La résolution en transitoire
ro.
dans ce mur d’épaisseur L est équivalente à la
résolution faite au paragraphe 5.2.1 pour un
mur d’épaisseur 2L dont les deux faces sont
soumises à de la convection; au centre il y a
une ligne de symétrie au travers de laquelle le
flux est nul et nous retrouvons ainsi la
géométrie de notre problème du four.
Biot =
h L 25 × 0.50
= = 10 ⇒ Bi −1 = 0.1
ri-p
k 1.25
To − T∞ 600 − 900 300

ka
= = = 0.357
Ti − T∞ 60 − 900 840
Sur la graphique 5.8 page 5.9, on trouve la valeur du nombre de Fourier :

αt Fo L2 k 1.25
Fo = ≈ 0.62 ⇒ t = α= = = 2.98 10−6
b
2
L α ρ CP 500 × 837
2
0.62 × 0.50
3a
t= = 51894 (14h 24 min.54 s )

2.98 10 −6
a)10 cm sous la surface correspond à x=0.40 et on utilise le graphique 5.19 page 5.10
x 0.40 T −T
= 0.80 ⇒ pour Bi −1 = 0.1 on trouve x =0.40 ∞ ≈ 0.415
.al
=
L 0.50 To − T∞
⇒ Tx =0.40 = 0.415 × ( 600 − 900 ) + 900 = 775o C
b) Pour trouver l’énergie reçue par la paroi, on utilise le graphique 5.10 page 5.10
w
Q
Bi 2 Fo = 10 2 × 0.62 = 62 ⇒ ≈ 0.75
Qo
Qo = Qmax = L S ρ CP (T∞ − Ti ) = 0.50 × 1× 500 × 837 × ( 900 − 60 ) = 1.758 108
ww
Q = 1.318 108 joules

TDEX- 5.12
m
Solutionnaire Ex-5.7 Géométrie sphérique:
0
Une bille métallique de 5 mm de rayon est initialement à une température uniforme de 400 C. Cette bille est
soumise à un traitement thermique en deux étapes.
co
0
a) Dans la première, on la refroidit dans l'air à 20 C pendant un temps tair nécessaire pour avoir une
0
température au centre de la sphère égale à 335 C. Pendant cette étape, le coefficient à la surface vaut
2
hair=10 W/m .K L'hypothèse d'une résistance thermique interne négligeable est-elle applicable à cette étape ?
Calculer le temps tair ?]
0
b) Dans une deuxième étape, la bille est alors introduite dans un bain d'eau à la température de 20 C.
ro.
2
Ce bain est fortement mélangé et le coefficient heau vaut alors 6000 W/m .K Quel est le temps requis pour
0
avoir une température au centre de 50 C ?
3
Données: k=20 W/m/K, : cP=1000 J/kg/K; ρ=3000 kg/m
-p
Réponses :
h R 10 0.005
Bi = = = 0.0025
a) Le nombre de Biot vaut : k 20
ari
Biot étant inférieur à 0.1 la résistance thermique interne est donc négligeable
⇒ T = T (t )
Bilan :
EIN − EOUT + Eg = Eaccumulé

k
4π R 3 ∂T
−4π R 2 h (T − Tair ) = ρ CP ( )
ab
3 ∂t
T t
 ρ CP R  ∂T  ρ CP R  ∂T
t
−
 3h  (T − T ) = ∂t ⇒  3h  ∫ (T − T ) = ∫ ∂t
−
air Tt =O air 0
 ρ CP R  Tt − Tair
t = − Ln
al3
 3h  Tt = 0 − Tair
application numérique :
 3000 1000 0.005  335 − 20
t = −  Ln = 500 0.1876 = 93.80 s
 3 10  400 − 20
w.
b) Le nombre de Biot vaut
h R 6000 0.005
Bi = = = 1.5
k 20
ww
TDEX- 5.13
om
Biot étant supérieur à 0.1 la résistance thermique interne n’est plus négligeable, ⇒ T = T (r , t ) et il faut
donc utiliser les abaques de Heisler.
Tcentre à t − Tair 50 − 20 1
θ0 = = = 0.095 Bi = 1.5 ⇒ Bi −1 = = 0.67
Tinitiale − Tair 335 − 20 1.5
o.c
*
Sur la figure p 5.15 pour les valeurs ci-dessus on lit la valeur de t :
αt t*R2
0.84 0.0052 0.84 0.0052
t = 2 ≈ 0.84 ⇒ t =
*
= =
R α  k   20 
   
 ρ CP   3000 1000 
-pr
0.84 0.0052
t= = 3.15 s
( 6.66 10−6 )
ari
Solutionnaire Ex-5.8 Géométrie mixte
La géométrie du problème est celle d’un cylindre court correspondant au cas (i) de la page 5.24. La solution
sera le produit de la solution d’un cylindre infini par un plan.
- La température maximale est celle au centre de l’objet : r=0 et x=0
- La température minimale est celle à la périphérie : r=r0=5 cm et x=L=4 cm
k
On cherche les températures au centre d’un cylindre infini et d’un plan infini après 10 minutes. On utilisera la
feuille Excel plutôt que les graphiques
ab
Au centre du plan :
pour le plan Tx =0 = 451.84 K Tx =0.04 = 404.27 K
pour le cylindre Tr =0 = 386.59 K Tr =0.05 = 354.02 K
à r = 0 et x = 0
al3
Tr =0, x = 0 − T∞ Tx =0 − T∞ Tr =0 − T∞ 451.84 − 300 386.59 − 300

= × = × = 0.082
Ti − T∞ Ti − T∞ Ti − T∞ 700 − 300 700 − 300
Tr =0, x =0 = 400 × 0.082 + 300 = 332.9 K
w.
à r = 0.05 et x = 0.04
Tr =0.05, x =0.04 − T∞ Tx =0.04 − T∞ Tr = 0.05 − T∞ 404.27 − 300 354.02 − 300
= × = × = 0.035
Ti − T∞ Ti − T∞ Ti − T∞ 700 − 300 700 − 300
ww
Tr =0, x =0 = 400 × 0.035 + 300 = 314 K

TDEX- 5.14
om
Solutionnaire Ex-5.9
Résolution d’un régime transitoire avec l’ hypothèse du régime quasi-permanent : détermination
du temps de solidification d'une goutte d'eau
o.c
glace
Rf Échange avec l’air à TA
et coefficient h
Eau à Tf
a) Hypothèses
ri -pr
∂T
ka
pas de génération de chaleur :EG=0, régime quasi permanent; = 0 k=constante
∂t
conduction unidirectionnelle : T = T (r )
Volume de contrôle
4π r 2 ∆r
EIN − EOUT + EG = E ACC ⇒ EIN − EOUT = 0
ab
Bilan de chaleur
4 π r 2 q′′r|r - 4 π r 2 q′′ r|r+∆r = 0 on divise par 4 π ∆r (attention ne pas simplifier par r 2 )

r 2 q′′r|r - r 2 q′′r|r+∆r r 2 q′′r|r - r 2 q′′r|r+∆r ∂ ( r 2 q′′r )
=0 lim| ∆r →0 =− =0
al3
∆r ∆r ∂r
 dT  dT C′ 1
⇒ r 2 qr" = C ′1 = r 2  −k  ⇒ =- 1 2
 dr  dr k r
−C ′1 1 1
⇒T = + C2 = C1 + C2
k r r
w.
b)
1
CF1 : r = R f , T = T f = C1 +C2
Rf
 ∂T 
ww
CF2 : r = R,  −k  = h(Tr = R − TA )
 ∂r  r = R
TDEX- 5.15
m
C1
CF1 ⇒ C2 = T f -
Rf
 ∂T   −C   C   C C 
CF2 ⇒  − k  = − k  2 1  = h(Tr = R − TA ) = h  1 + C 2 − TA  = h  1 + T f - 1 − TA 
co
 ∂r  r = R  R   R   R Rf 
 
 k h h  h ( T f − TA ) ( T f − TA )
C1  2 − +
R  = h ( T f − TA ) ⇒ C1 = =
 R R f   k h h   k 1 1 
 2 − +   2 − + 
R R f 
ro.
R R Rf   hR
c)
 ∂T   kC   k  ( T f − TA )
Q f = 4π R 2 q′′r|r = R = 4π R 2  − k  = 4π R 2  21  = 4π R 2  2 
 ∂r  r = R  R   R  k 1 1 
 2 − + 
Qf =
4 π k (T f − TA )
1 1
− + 2
R f R hR
k
ri-p  hR R Rf 
d)
 
ka
 4 π k (T − T ) 
⇒ − 4π R 2f ∂R f ρ ∆H f =   ∂t
f A
−4π R 2f ∂R f ρ ∆H f = Q f ∂t
 1 1 k 
 R − R + hR 2 
 f 
ab
 1 1 k  k
R 2f  − + 2  ∂R f = −
 R f R hR  ρ ∆ H
(T f − TA ) ∂t
  f
 R f kR f 
2 2
 1 k  k
 R f − + 2  ∂R f =  R f − R 2f  + 2   ∂R f = −
  R hR   ρ ∆H f
(T f − TA ) ∂t
 R hR  
al3
R f =0 t =t f
 2  1 k  k
∫  R f + R f  − R + hR 2   ∂R f = − ρ ∆H (T f − TA ) ∫ ∂t
Rf =R    f t =0
R f =0
 R 2f R 3f  k 1  k
 + 
3  hR 2
−  
R  R =R
=− ( T f − TA ) [t ]0f
t
w.
 2 ρ ∆H f
f
 R2 R3  k 1   R 2 kR R 2  k
− + 
3  hR 2
− 
R
 = −  + − =− (T − TA ) t f
ρ ∆H f f
 2  2 3h 3 
ww
ρ ∆H f  R 2 kR 
tf =  + 
k (T f − TA )  6 3h 
ww
w.
al3
ab
kari
-pro.
co
m
TDEX- 6.1
m
Objectifs: Bien comprendre les contributions d==advection (ou convection) dans les bilans de
chaleur quand il y a mouvement de matière au travers des surfaces délimitant le
co
volume de contrôle.
Ex-6.1 (S) Refroidissement d'un fil de polymère (examen de 1996)
Un fil de polymère (P.), de rayon R, produit par extrusion,

sort de la filière avec une vitesse V. Ce fil doit alors être
ro.
refroidi, de T1 à T2, dans un bain d'eau à la température
Ta. Le coefficient de transfert de chaleur à la surface du fil
est constant est vaut ha.
Négligeant la conduction axiale, le gradient radial

de température à l'intérieur du fil, et l'échange avec l'air
avant l'entrée dans l'eau:
a)
b)
Réponse :
ri-p
déterminer l=expression du profil de température
T(z), avec z distance parcourue dans l=eau.
déterminer la distance, L, que le fil doit parcourir
pour atteindre la température T2.
ka
T(z) - T a ρ C P VR -
Ln T 1 T a
2hZ
= e- ρ C pVR L=
T1 -T a 2h T2 -T a
Ex-6.2 (S) Sphère poreuse (examen 1996)

b
Une sphère creuse (rayon intérieur r1 ,rayon

extérieur r2) est faite d'un matériau poreux. Les
3a
températures des faces intérieure et extérieure sont

3
T1 et T2. Un débit d'air QA (m /s) s'écoule au travers
du matériau poreux de l'intérieur vers l'extérieur. On Débit d’air
appelle k, la conductivité apparente de la paroi
(air+matériau poreux). On souhaite déterminer le
profil de température dans la paroi.
.al
a) Faites un bilan de chaleur sur un élément

de volume approprié et obtenez l'équation
différentielle que doit vérifier la température.
b) Trouvez l'expression du profil de
température
w
Réponse :
Matériau poreux
β β
ww
T(r) - T 1 e- r - e- rsub 1 Q ρ CP
= β β
avec β = A
T 2 -T1 - -
e r 2 - e r1 4πk
TDEX- 6.2
m
Travaux dirigés No 6 (suite)
Ex-6.3 (S) Écoulement et chauffage sur un plan incliné (examen 91)
Un liquide (conductivité k, masse volumique ρ,
co
capacité calorifique CP) s'écoule sur un plan
incliné (longueur L et largeur W) dont la
surface est maintenue à une température
constante TW (voir figure ci-dessous). A
l'entrée du plan le fluide, d=épaisseur δ, est à
une température uniforme TO (<TW ) et on
ro.
suppose que le transfert de chaleur avec l'air
environnant à la surface du liquide est
négligeable. En régime laminaire, le profil de
vitesse sur le plan est donné par:
  x 2 
Épaisseur *
ri-p v Z = V max 1-   
  δ  
La température dans le liquide est donc
fonction de la profondeur et de la distance
parcourue depuis l'entrée sur le plan.
ka
a) Faites un bilan de chaleur (le plus général possible) sur un élement de volume approprié dans le fluide et
obtenez l'équation différentielle que doit vérifier la température du liquide (ne pas résoudre cette équation)
2 2
∂T ∂T ∂T
k + k - ρ C P vz =0
b
2 2
réponse : ∂z ∂x ∂z
b) Simplifier l'équation précédente dans le cas où, suivant la direction z, la convection ou advection est
3a
beaucoup plus importante que la conduction.

2
∂T ∂T
+k 2
- ρ CP v z =0
réponse : ∂x ∂z
c) Quelles sont les conditions frontières nécessaires à la résolution des équations obtenues en a et b
.al
réponses :
pour b)CF1 : z = 0 T = T 0 ; CF2 : x = δ T = Tw ;CF3 : x = 0 - k ∂T/ ∂x = 0 ;
pour a)CF4 z : bilan d ′énergie sur le film entre z = 0 et z
w
ww
TDEX- 6.3
Ex-6.4 (S) Couche limite sur une plaque
m
Une solution approximative au problème de la couche limite
thermique sur une plaque est donnée par:
co
(T - T S ) 3 y 1 y 3
= ( )- ( )
(T ∞ -T S ) 2 δ 2 δ
a v e c δ = 5 x R e -1 /2 P r -1 /3
ro.
a)
et R e = ( ρ V ∞ x )/ µ
Donner l'expression du nombre de Nusselt local à une position x donnée ( NuX=hX x /k) en fonction des
nombres de Reynolds et Prandtl.
3k h xx
hx = Nux = = 0.3 Re 1/2 Pr 1/3
2δ k
Réponse :
ri-p
b) Donner l'expression du coefficient moyen de transfert de chaleur, h, entre l'entrée de la plaque et cette
position x. Comparez les nombres de Nusselt local et moyen (Nu=hx/k)
Réponse :
Nu = 2N u x
Ex-6.5 (S) Paroi poreuse

ka
Une paroi poreuse d'épaisseur b, de largeur W et de hauteur L,
sépare une chambre froide d'une pièce tempérée. Cette paroi poreuse
3
est traversée par un débit d'air DA(m /s) qui s'écoule de la chambre b
froide vers l'extérieur. La face froide de la paroi est à une température
T1 alors que la face extérieure est à la température T2. La conductivité
apparente dans la paroi poreuse (air+ matériau poreux) est k . DA
ab
On souhaite déterminer le profil de température dans le matériau T1 T2

poreux ainsi que la quantité de chaleur qui rentre dans la chambre
froide.
al3
a) Faites un bilan de chaleur sur un élément de volume

approprié dans la paroi
b) Obtenez l'expression du profil de température
T ( z ) − T1 e A z − 1 D ρC
w.
= Ab avec A = A P
Réponse :
T2 − T1 e −1 WLk
c) Déterminer l'expression de la quantité de chaleur, Q, qui rentre dans la chambre froide par la face
intérieure.
ww
A
Q = WLk (T2 − T1 ) Ab
Réponse : e −1
d) Sachant que la quantité, QO, qui rentrerait si il n'y avait pas d'écoulement d'air vaut, en valeur absolue:
TDEX- 6.4
QO=WLkAPP(T2-T1)/b, commentez l'effet du débit d'air, sur la variation du rapport Q/QO
m
Q Ab
= Ab
Réponse :
Q0 e − 1 (Q/Q0) vaut respectivement 0.77, 0.58, 0.31, 0.03 quand Ab vaut
respectivement 0.5, 1, 2, 5 : donc la chaleur qui entre dans la chambre froide décroit quand DA
co
augmente.
Ex-6.6 (S) Refroidissement d'une feuille métallique (examen 2003)
ro.
Une feuille métallique de largeur W et d'épaisseur b, à la
température uniforme T1 sort d'un four avec une vitesse
V. La feuille échange avec l'air ambiant à la température
TA et les coefficients de transfert de chaleur sur les faces
supérieure et inférieure de la feuille valent respectivement
hS et hI. On suppose que le gradient de température dans
ri-p
l'épaisseur de la feuille est négligeable.
On supposera que hSb/k <<1 et hIb/k <<1
a) Faites le bilan de chaleur sur un élément de

volume approprié et obtenez l’équation
différentielle que doit vérifier la température dans
la feuille.
ka
Dans le cas où la vitesse V est élevée, on peut supposer que la conduction axiale est négligeable devant
l'advection axiale. Simplifier l’équation obtenue en a) et obtenez l'expression de la variation de la
température de la plaque avec la position z?
ab
al3
w.
ww
TDEX- 6.5
om
Solutionnaire: 6.1 Refroidissement d'un fil de polymère
a) On peut représenter le problème comme ci-dessous :
o.c
-pr
∂
Hypothèses : Gradient radial négligeable YT=T(z) ; Eg=0 régime permanent =0
∂t
Volume de contrôle π R 2 ∆z
Ce qui rentre (conduction+ advection) –Ce qui sort (conduction+advection+convection)=0
ri
π R 2 q"z z − π R 2 q"z + π R 2V ρ CP (Tz − Tref ) − π R 2V ρ CP (Tz − Tref ) − 2π R∆z h (Tz − Ta ) = 0
ka
z +∆z z z +∆z

terme de conduction axiale négligeable
Comme on suppose que le terme de conduction axiale est négligeable devant les autres termes, le bilan
se réduit à :
ab
π R 2V ρ CP (Tz − Tref ) z − π R 2V ρ CP (Tz − Tref ) z +∆z − 2π R∆z h (Tz − Ta ) = 0

RV ρ CP (Tz − Tref ) − RV ρ CP (Tz − Tref ) − 2∆zh(Tz − Ta ) = 0
z z +∆z
Tz − Tz +∆z Tz − Tz +∆z ∂T
RV ρ CP − 2h (Tz − Ta ) = 0 lim =− z
∆z ∆z ∂z
al3
∆z → 0
∂Tz 2h ∂T  2h 
− − (Tz − Ta ) = 0 ⇒ z = −   (Tz − Ta ) = 0
∂z RV ρ CP ∂z  RV ρ CP 
∂θ ∂Tz ∂θ  2h  ∂θ  2h 
posons θ = (Tz − Ta ) ⇒ = ⇒ = − θ ⇒ = −  ∂z
∂z ∂z ∂z  RV ρ CP  θ  RV ρ CP 
w.
 2h 
T −T  2h  Tz − Ta  2h  Tz − Ta − z
 RV ρ CP 
Ln z a = −   z ⇒ Ln = −   z ⇒ = e
Tz = 0 − Ta  RV ρ CP  T1 − Ta  RV ρ CP  T1 − Ta
ww
b) à z=L, T=T2 :
T2 − Ta  2h   RV ρ CP  T1 − Ta
Ln = − L ⇒ L=   Ln T − T
T1 − Ta  RV ρ C P   2 h 2 a
TDEX- 6.6
m
Solutionnaire 6.2: Sphère poreuse (examen 1996)
Une sphère creuse (rayon intérieur r1 ,rayon extérieur r2) est
faite d'un matériau poreux. Les températures des faces
co
3
intérieure et extérieure sont T1 et T2. Un débit d'air QA (m /s)
s'écoule au travers du matériau poreux de l'intérieur vers
l'extérieur. On appelle k, la conductivité apparente de la Débit d’air
paroi (air+matériau poreux). On souhaite déterminer le profil
de température dans la paroi.
Faites un bilan de chaleur sur un élément de volume
ro.
approprié et obtenez l'équation différentielle que doit vérifier
la température.
a) Faites un bilan de chaleur sur un élément de volume
approprié et obtenez l'équation différentielle que doit
vérifier la température.
b) Trouvez l'expression du profil de température
a)
ri-p
Hypothèses :
T=T(r) k=conductivité apparente = Constante

ka
∂
=0
Eg=0 régime permanent
∂t
Volume de contrôle 4π r 2 ∆r
ab
Ce qui rentre(conduction+ advection) – Ce qui sort (conduction+advection)=0

al3
4π r 2 qr" − 4π r 2 qr" + QA ρ CP (T − Tref ) − QA ρ CP (T − Tref ) =0

r r +∆r r r +∆r
on divise par 4π∆r et on fait tendre ∆r vers 0

w.
4π r 2 qr" − 4π r 2 qr" QA ρ CP (T − Tref ) r − QA ρ CP (T − Tref ) r +∆r

r r +∆r
+ =0
4π∆r 4π∆r
∂ ( r 2 qr" ) QA ρ CP ∂T ∂T
− − =0 mais qr" = − k
ww
∂r 4π ∂r ∂r
TDEX- 6.7
m
 ∂T   ∂T 
∂  r2  ∂  r2 
k 
∂r  Q
= A P
ρ C ∂T
⇔  ∂r  =  QA ρ CP  ∂T = β ∂T
 
∂r 4π ∂r ∂r  4π k  ∂r ∂r
.co
on intègre une première fois, puis une deuxième fois :
 2 ∂T  ∂T ∂r
r  = β T + C1 ⇒ = 2
 ∂r  β T + C1 r
1 1 β
Ln ( β T + C1 ) = − + C2 ⇒ Ln ( β T + C1 ) = − + C2
β r r
o
β
− + C2  1   − β +C2 
β T + C1 = e r
⇒ T =    e r − C1 
 β  
à r=r1 T=T1 et à r=r2 T=T2

-pr
Pour trouver C1 et C2, on utilise les deux conditions frontières :
Pour éviter de calculer C1 puis C2,on calcule directement l’expression (T-T1)/(T2- T1) :
ari
  1   − r1 +C2   −β − 
β β
1   − βr +C2
 −    e − C1  eC2 e − e 1 
r
β  e − C1 r
T − T1    β    
= =  
T2 − T1  1   − rβ +C2   1   − r + C2
β
  −r
β
− 
β
bk
 β   e  −    e 1 − C1  e 2 e 2 − e 1 
C r
2
− C1
   β    
    
 −β − 
β
 
r1
e r
− e
T − T1   

3a
=
T2 − T1  − r β
− 
β
e 2 − e 1 
r
 
 
w .al
ww
TDEX- 6.8
m
Solutionnaire Ex-6.3 Écoulement sur un plan incliné
co
ro.
ri-p
ka
∂
a) Hypothèses : T=T(x,z) ; Eg=0 régime permanent =0
ab
∂t
Volume de contrôle ∆x∆zW
Ce qui rentre (conduction+ advection) –Ce qui sort (conduction+advection)=0
al3
W ∆x q"z − W ∆x qz" + W ∆z qx" − W ∆z qx" + W ∆xVz ρ CP (Tx , z − Tref ) − W ∆xVz ρ CP (Tx , z − Tref ) =0
z z +∆z x x +∆x z z +∆z
W ∆x q"z − W ∆x qz"
z z +∆z
W ∆z q"x − W ∆z qx"
x x +∆x
W ∆xVz ρ CP (Tx , z − Tx , z +∆z )
+ + =0
W ∆x ∆z W ∆x ∆z W ∆x ∆z
q"z − q"z q"x − q"x (T − Tx , z +∆z )
w.
z +∆z x +∆x x,z

z
+ x
+ Vz ρ CP =0 on fait tendre ∆x ∆z vers 0
∆z ∆x ∆z
∂q"z ∂q"x ∂T ∂Tx , z ∂Tx , z
− − − Vz ρ C P x , z = 0 qz" = − k qx" = − k
∂z ∂x ∂z ∂z ∂x
ww
2
∂ Tx , z 2
∂ Tx , z ∂Tx , z 2
∂ Tx , z 2
∂ Tx , z   x 2  ∂T
k +k − Vz ρ C P =0 ⇔ k +k − Vmax 1 −    ρ CP x , z = 0
∂z 2 ∂x 2 ∂z ∂z 2 ∂x 2   δ   ∂z
b)
TDEX- 6.9
m
∂ 2Tx , z ∂ 2Tx , z  x  2
∂T
k 2
+k 2
− Vmax 1 −    ρ CP x , z = 0

∂z
∂x   δ   ∂z
conduction suivant z conduction
co
négligeable suivant x advection suivant z
devant l ' advection
∂ 2Tx , z   x 2  ∂T
k − Vmax 1 −    ρ CP x , z = 0
∂x 2   δ   ∂z
ro.
e er
c) L’équation différentielle obtenue en b) est du 2 ordre par rapport à x et du 1 ordre par rapport à
z, il faut donc 3 conditions frontières :
CF1 : à z = 0 Tx ,0 = T 0 pour 0 ≤ x ≤ δ
CF2 : à x = δ
-p
Tδ , z = Tw pour 0 ≤ z ≤ L
CF3 : échange de chaleur négligeable à la surface libre du liquide

ari
∂Tx , z
à x=0 -k = 0 pour 0 ≤ z ≤ L
∂x x =0
k
e e
L’équation différentielle obtenue en a) est du 2 ordre par rapport à x et du 2 ordre par rapport à z, il faut
donc 4 conditions frontières. On utilisera les 3 conditions ci-dessus.
e
Pour obtenir une 4 condition, plusieurs possibilités:
ab
a) on fait un bilan de chaleur macroscopique sur toute l’épaisseur du film entre l’entrée z=0 et z=L :
x =δ
  x 2  L
 ∂Tx , z 
∫x=0 max   δ   P ( x,0 x, L )
V 1 − ρ C T − T Wdx = ∫z =0 −k  ∂x  Wdz
 
al3
x =δ
b) Le plan est suffisamment long pour que le profil de température suivant z soit établi :
∂Tx , z
à x= L -k = 0 pour 0 ≤ x ≤ δ
∂z z=L
c) Le plan est suffisamment long pour que la température du film soit égale à TW
w.
à x= L Tx , z = L = TW pour 0 ≤ x ≤ δ
ww
TDEX- 6.10
om
Solutionnaire : Ex-6.4
Une solution approximative au problème de la couche limite
thermique sur une plaque est donnée par:
(T - T S ) 3 y 1 y 3
= ( )- ( )
(T ∞ -T S ) 2 δ 2 δ
o.c
a v e c δ = 5 x R e -1 /2 P r -1 /3
a)
et R e = ( ρ V ∞ x )/ µ
Donner l'expression du nombre de Nusselt local à une position x donnée ( NuX=hX x /k) en fonction des
nombres de Reynolds et Prandtl.
-pr
b) Donner l'expression du coefficient moyen de transfert de chaleur, h, entre l'entrée de la plaque et cette
position x. Comparez les nombres de Nusselt local et moyen (Nu=hx/k)
Réponse:
a) à la paroi
ari
∂T
q′′ = −k = hx (Ty =0 − T∞ )
∂y y =0
(T - T S ) 3 y 1 y 3 ∂T
connaissant = ( ) - ( ) on peut calculer et Ty = 0
(T ∞ -T S ) 2 δ 2 δ ∂y
bk
∂T  3 1 1 3y2  3
= (T ∞ -T S ) - 3 
= (T ∞ -T S )
∂y  2 δ 2 δ  y =0 2δ
(Ty = 0 - T S )
= 0 ⇒ Ty =0 = T S
(T ∞ -T S )
3a
∂T 3
−k = hx (Ty =0 − T∞ ) ⇔ −k( T ∞ - T S ) = hx (TS − T∞ )
∂y y =0
2δ
3k
⇒ hx =
al
2δ
3k
x
hx x 2δ 3x
Nu x = = = sachant que δ = 5 x Re -1/2 Pr -1/3
w.
k k 2δ
3x 3x 1/2
Nu x = = = 0.3 Re Pr 1/3
2 ( 5 x Re -1/2
Pr
-1/3
) 10 x Re -1/2
Pr
-1/3
ww
b)
TDEX- 6.11
m
h x k
Nu x = x = 0.3 Re Pr 1/3 ⇒ hx = 0.3 Re Pr 1/3
1/2 1/2
k x
1/2
o
k  ρ V∞ x 
x x x
1 1 k 1
hx = ∫ hx dx = ∫ 0.3 Re Pr 1/3dx = ∫ 0.3 
1/2 1/3
 Pr dx
x0 x0 x x0 x  µ 
o.c
1/2
0.3 k  ρ V∞ 
x
x1/2
 Pr ∫
1/3
hx =  dx
x  µ  0
x
x x
x1/2
-pr
x
∫
0
x
dx = ∫ x −1/2 dx =  2 x1/2  = 2 x1/2
0
0
0.3 k  ρ V∞ 
1/2
 k  ρ V∞ x 
1/2
  k
1/2

⇒ hx =   Pr
1/3
2 x 1/2
= 2  0.3   Pr
1/3
 = 2  0.3 Re 1/3
Pr  = 2 hx
x  µ  x  µ 
ri    x 
ka
hx x 2hx x hx
Nu x = = = 2 x = 2 Nu x
k k k
Nu x == 2 Nu x
ab
al3
w.
ww
TDEX- 6.12
m
Solutionnaire 6.5 : Paroi poreuse
co
ro.
ri-p
ka
ab
∂
al3
Hypothèses : T=T(z) ; Eg=0 régime permanent =0

∂t
Volume de contrôle W L ∆z

Ce qui rentre (conduction+ advection) –Ce qui sort (conduction+advection)=0
w.
WL qz − WL q"z + DA ρ CP (Tz − Tref ) − DA ρ CP (Tz − Tref ) =0

z z +∆z z z +∆z
WL qz" − WL qz" DA ρ CP (Tz − Tz +∆z ) qz" − q"z DA ρ CP (Tz − Tz +∆z )

z z +∆z z z +∆z
+ = + =0
WL∆z WL∆z ∆z WL ∆z
ww
TDEX- 6.13
om
on fait tendre ∆z vers 0
∂q"z  DA ρ CP  ∂Tz ∂qz"  D ρ C  ∂T  D ρC 
− −  = 0 ⇒ = −  A P  z ⇒ qz" = −  A P  Tz + C1
∂z  WL  ∂z ∂z  WL  ∂z  WL 
∂T  D ρC  ∂Tz  DA ρ CP  C1  D ρC 
qz" = − k z = −  A P  Tz + C1 ⇒ =  Tz − =  A P  Tz + C1′ = α Tz + C1′
∂z  WL  ∂z  kWL  k  kWL 
o.c
∂Tz ∂Tz 1
= ∂z ⇒ ∫ = ∫ ∂z ⇒ ln α Tz + C1′ = z + C2
α Tz + C1′ α Tz + C1′ α
C2′ eα z − C1′
⇒ α Tz + C1′ = e ( 2 ) = eα z +α C2 = eα C2 eα z = C2′ eα z
α z +C
⇒ Tz =
α
Pour trouver les constantes d’intégration, on utilise les conditions frontières à z=0 T=T1 et à z=b T=T2 .
-pr
Tz − Tz =0
Plus simplement on va calculer directement l’expression :
Tz =b − Tz =0
C2′ eα z − C1′ C2′eα 0 − C1′

−
Tz − Tz =0 T −T eα z − 1
ari
= z 1= α α =
Tz =b − Tz =0 T2 − T1 C2′eα b − C1′ C2′eα 0 − C1′ eα b − 1
−
α α
αz
Tz − T1 e − 1 DA ρ C P
= avec α =
T2 − T1 eα b − 1 WL
bk
La quantité de chaleur, Q, entre dans la chambre froide à z=0. En ce point, c’est la conduction qui apporte
la chaleur de l’extérieur :
 ∂T  eα z − 1 ∂Tz (T2 − T1 ) α eα z
3a
Q = WLq"z = WL  − k z  Tz = (T2 − T1 ) + T1 −k = −k
 ∂z  z = 0 eα b − 1 ∂z z =0 eα b − 1 z =0
α
Q = −WLk (T2 − T1 ) αb
signe négatif car le flux est dans le sens contraire de l ' axe
e −1
al
Si on compare avec le flux Q0 entrant sans écoulement: (voir Notes de cours paragraphe 3.2)
α
WLk (T2 − T1 )
Q
= e −1 = α b
αb
w.
Q0 WLk (T2 − T1 ) eα b − 1
b
(Q/Q0) vaut respectivement 0.77, 0.58, 0.31, 0.03 quand α b vaut respectivement 0.5, 1, 2, 5 : donc la
quantité de chaleur qui entre dans la chambre froide décroit quand DAaugmente: le flux d’air sortant de la
paroi s’oppose à l’entrée de la chaleur (cet effet est appelé en anglais « transpiration cooling» (voir le livre
ww
Transport Phenomena).
TDEX- 6.14
m
Solutionnaire : Ex-6.6 Refroidissement d'une feuille métallique (examen 2003)
co
température uniforme T1 sort d'un four avec une vitesse
V. La feuille échange avec l'air ambiant à la température
TA et les coefficients de transfert de chaleur sur les faces
supérieure et inférieure de la feuille valent respectivement
hS et hI. On suppose que le gradient de température dans
l'épaisseur de la feuille est négligeable.
ro.
a) Faites le bilan de chaleur sur un élément de
volume approprié et obtenez l’équation
différentielle que doit vérifier la température dans
la feuille.
b) Dans le cas où la vitesse V est élevée, on peut supposer que la conduction axiale est négligeable
ri-p
devant l'advection axiale. Simplifier l’équation obtenue en a) et obtenez l'expression de la variation de
la température de la plaque avec la position z?
Réponses :
a) Hypothèses:
Comme Biot est très inférieur à 1, on supposera que la température ne varie pas dans l’épaisseur de la plaque
ka
et varie uniquement avec la position axiale z : T=T(z)
∂T
=0
k=constante, régime permanent ∂t
Volume de contrôle:
W b ∆z
ab
Bilan d’énergie thermique:
IN OUT
conduction Wbq′′ z − Wbq′′ z +∆z
al3
advection Wb ρ CPV (T − Tref ) − Wb ρ CPV (T − Tref )

z z +∆z
convection dessous − W ∆z hI (T − Ta )
convection dessus − W ∆z hS (T − Ta ) = 0
w.
On divise par
W b ∆z et on fait tendre ∆z vers 0
ww
TDEX- 6.15
m
q′′ z − q′′ z +∆z  Tz − Tz +∆z  ( h + hS ) T − T = 0
+ ρ CPV   − I ( a)
∆z  ∆z  b
∂q′′ ∂T ( h + hS ) T − T = 0
− − ρ CPV − I ( a)
co
∂z ∂z b
∂ 2T ∂T ( h + hS ) T − T = 0
k 2 − ρ CPV − I ( a)
∂z ∂z b
ro.
b) Si la conduction axiale est négligeable devant l’advection alors :
≈0
∂ 2T ∂T ( hI + hS )
− ρ CPV − (T − Ta ) = 0
-p
k 2
∂z ∂z b
∂T ( hI + hS ) ∂T ( h + hS ) T − T
− ρ CPV − (T − Ta ) = 0 ⇔ ρ CPV =− I ( a)
∂z b ∂z b
ari
∂T ( h + hS ) ∂z
=− I ⇔
∂ (T − Ta ) ( h + hS ) ∂z
=− I
(T − Ta ) b ρ CPV (T − Ta ) b ρ CPV
T
∂ (T − Ta ) ( hI + hS ) z
∫ (T − Ta )
=−
b ρ CPV z∫=0
∂z
bk
T1
Ln
(T − Ta ) ( h + hS ) z
=− I
(T1 − Ta ) b ρ CPV
al 3a
w.
ww
TDEX- 6.16
m
Exemple d’examen No 2
GCH-18243, Examen No 2 Notation sur 25 pts Notes de cours autorisées:
co
Date: 24 octobre 2003 Durée 110 min.
si vous utilisez des solutions graphiques, mentionnez les numéros et pages de ces graphiques.
ro.
1) Refroidissement d’une plaque de cuivre (4. pts)
Une plaque de cuivre (longueur L= 8 cm, largeur W= 5 cm, épaisseur b=1 mm) sort d’un
four à une température uniforme de 150 oC. On la plonge immédiatement dans un fluide de
refroidissement qui est à une température de 15 oC et le coefficient de transfert de chaleur à la
surface de la plaque vaut 250 W/m2.K Quel est le temps nécessaire pour refroidir la plaque à 25
o
Données:
ri-p
C ? Justifier toutes vos hypothèses et faites la démonstration du bilan thermique.
Pour le cuivre: k=398.0 W/m/K; =8933 kg/m3; CP=389 J/kg/K
(réponse : 17.5 s)
ka
2) Refroidissement d’une bille de plastique (4. pts)
Une bille de plastique de 35 mm de rayon est initialement à une température uniforme de

85 C. On la place dans un courant d’air frais à 10oC. Sachant que la température au centre de la
o
bille est de 41oC après 184 minutes de refroidissement, trouver la valeur du coefficient de
ab
transfert de chaleur à la surface de la bille durant ce refroidissement.
Données: Pour le plastique: k=0.5 W/m/K; =1170 kg/m3; CP=1930 J/kg/K

(réponse : 2.38)
al3
3) Méthode explicite de différences finies en transitoire (4.5 pts)
Pour le maillage ci-contre d’une pièce métallique d’ épaisseur L et

dans le cas d’un régime transitoire écrivez le bilan thermique
w.
autour du noeud 5 et en utilisant la méthode explicite, obtenez la

relation permettant de relier T5 (t+)t) en fonction des températures
des noeuds au temps t.
(On a )x=)y et on pourra utiliser l’hypothèse des flux entrants)

ww
TDEX- 6.17
m
4) Profondeur pour enterrer les conduites d’eau (4.5 pts)
Les conduites d’eau potable enterrées dans le sol doivent être à l’abri du gel. Déterminer la
profondeur d’enfouissement sécuritaire correspondant aux conditions extrêmes suivantes:
co
- le 31 octobre la température du sol est uniforme et égale à 10 oC
- le 1ernovembre, il y a refroidissement à -15 oC, et cette température se maintient pendant
120 jours avec un fort vent qui impose un coefficient de transfert de chaleur à la surface du sol de
35 W/m2.K
- et le 1er mars, il fait beau !
ro.
Propriétés du sol: k=0.52 W/m/K; D=2050 kg/m3; CP=1840 J/kg/K
(réponse 1.38 m , voir Ex-5.5)
5) Quiz (2 pts): Cochez la case appropriée
ri-p
PROPOSITIONS
le coefficient h sur un objet dépend de la conductivité thermique de l'objet

OUI
VRAI
NON
FAUX
lors d’un processus transitoire de refroidissement dans un objet, plus la résistance thermique
interne est élevée, plus le gradient de température dans l'objet est faible
ka
pour calculer le nombre de Biot associé à un objet, on utilise uniquement les propriétés du
matériau de l'objet
pour une géométrie donnée, le nombre de Nusselt est uniquement fonction du nombre de
Reynolds
ab
6) Refroidissement d'une feuille métallique (6 pts)

al3

température uniforme T1 sort d'un four avec une vitesse V. La
feuille échange avec l'air ambiant à la température TA et les
coefficients de transfert de chaleur sur les faces supérieure et
inférieure de la feuille valent respectivement hS et hI. On
suppose que le gradient de température dans l'épaisseur de la
w.
feuille est négligeable.

a) Faites le bilan de chaleur sur un élément de volume approprié et obtenez l’équation
différentielle que doit vérifier la température dans la feuille.
b) Dans le cas où la vitesse V est élevée, on peut supposer que la conduction axiale est
ww
négligeable devant l'advection axiale. Simplifier l’équation obtenue en a) et obtenez

l'expression de la variation de la température de la plaque avec la position z?
(réponses voir Ex-6.6)
ww
w.
al3
ab
kari
-pro.
co
m
TDEX- 7.1
m
Travaux dirigés et exercices du chapitre 7 : Convection externe
co
Ex-7.1 (S) Calcul du coefficient h sur une sphère à l=aide d=une corrélation
On laisse tomber dans un réservoir profond rempli d'une huile à une température constante TH, une
bille sphérique de rayon R. On néglige la phase d'accélération de la bille et on suppose que la vitesse
terminale de chute VT est atteinte immédiatement.
ro.
Supposant que les propriétés de l'huile sont constantes, quelle est la valeur du coefficient h de transfert de
chaleur à la surface de la sphère ? [réponse : 651.1]
3
Données: la bille: R =0.003 m, VT=0.23 m/s, CPS = 1000 J/kg/K , ρS= 2600 kg/m
l'huile : µL= 0.055 Pa.s, CPL = 1885 J/kg/K
3
ρL= 900 kg/m , kL = 0.145 W/m/K
Nu = (h D/kL) = 2+ 0.6 Re ri-p

0.5
Pr
0.33
, avec Re=ρLVD/µL, Pr=cPL µL/kL
Ex-7.2 (S) Calcul d=un coefficient h à partir d=une solution analytique
a) Faire l’exercice 6.4

ka
b) Cas limite pour l’écoulement autour d=une sphère
Démontrer que pour la conduction dans un fluide autour d'une sphère, on peut écrire Nu=2.
On suppose qu'une sphère de rayon R est plongée dans un fluide stagnant dont la température loin de la
sphère est T. La température à la surface de la sphère est TS. Déterminer l'expression du profil de
ab
température pour r>=R. Calculer le flux à la surface et trouver la valeur du nombre de Nusselt en utilisant
la définition du coefficient h.
Ex-7.3(S) Calcul du coefficient h autour d=un cylindre à l=aide d=une corrélation

al3
Une huile, dont les propriétés sont données à l=exercice 1, s=écoule avec une vitesse de 0.55 m/s
perpendiculairement à un cylindre de 0.05 m de rayon.
Quelle est la valeur du coefficient h à la surface du cylindre en utilisant la relation 7.45 p 7.4.1 ?
[réponse : 210.8]
Comparer avec la valeur obtenue en utilisant l=équation 7.47 [réponse : 246.5]
w.
ww
TDEX- 7.2
m
Solutionnaire 7.1 Calcul d’un coefficient h
hD ρ V D 900 × 0.23 × 0.006
Nu = = 2 + 0.6 Re0.5 Pr 0.33 Re = L T = = 22.58
kL µL 0.055
CPL µ L 1885 × 0.055
.co
Pr = = = 715 Nu = 2 + 0.6 × 22.580.5 × 7150.33 = 26.94
kL 0.145
26.94 × 0.145
h= = 651.1
0.006
o
Solutionnaire Ex-7.2b :
-pr
Démontrer que pour la conduction autour d'une sphère, on peut écrire Nu=2.
On suppose qu'une sphère de rayon R est plongée dans un fluide stagnant dont la température loin de la
sphère est T. La température à la surface de la sphère est TS. Déterminer l'expression du profil de
température pour r>=R. Calculer le flux à la surface et trouver la valeur du nombre de Nusselt en utilisant
la définition du coefficient h.
ari
1) Hypothèses :
∂
T=T(r) k=Constante EG=0 régime permanent =0
∂t
2) Volume de contrôle 4π r 2 ∆r
k

ab

(terme uniquement de conduction si fluide stagnant)
4π r 2 qr" − 4π r 2 qr" =0
r r +∆r
al3
On divise par 4π ∆r et on fait tendre ∆r vers 0
4π r 2 qr" − 4π r 2 qr"
r r +∆r
=0
4π ∆r
w.
r 2 qr" − r 2 qr" ∂ 2 "

Par définition : lim
∆r → 0
r
∆r
r +∆r
=−
∂r
( r qr ) on obtient donc :
∂ 2 " ∂ 2 "
− ( r qr ) = 0 ⇒ ( r qr ) = 0
ww
∂r ∂r
On intègre une première fois: (r q ) = C

2 "
r 1
TDEX- 7.3
m
" ∂T
On exprime alors la densité de flux par la loi de Fourrier: qr = −k
∂r
 ∂T  ∂T C
r 2  −k  = C1 ⇒ = − 12
co
 ∂r  ∂r kr
et on intègre une deuxième fois pour obtenir la variation T(r):
ro.
C1  1
Tr = −  −  + C2 [1]
k  r
On utilise les conditions frontières pour déterminer les constantes d’intégration C1 et C2
CF1: à r=R , TR = TS
CF2: à r=∞,
ri-pTR =∞ = T∞
Calculons directement :
 C1  1    C1  1  
−  −  + C2  −  − −  + C2 
Tr − T∞  k  r   k  ∞  =R
ka
=
TS − T∞  C1  1    C1  1   r
 −  −  + C2  −  − −  + C2 
 k  R   k  ∞ 
2) Calculons le flux à la surface :

ab
∂T 1 k
qr" = − k = (− k ) (TS − T∞ ) R (− 2
) = (TS − T∞ )
∂r r =R r r=R R
Ce flux sortant peut s’exprimer par la loi de refroidissement de Newton:

al3
k k
qr" = h (TS − T∞ ) = (TS − T∞ ) ⇒ h =
R R
k
2R
w.
hD R
Nu = = =2
k k
Commentaire : Dans la section 7.5 des notes de cours, on a mentionné que la plupart des corrélations
ww
pour l’écoulement autour d’une sphère sont de la forme:

Nu = 2 + f ( Re , Pr )
On vient ainsi d’expliquer la valeur limite de 2 dans ces corrélations.

TDEX- 7.4
m
Solutionnaire 7.3
On utilise les expressions de la page 7.6 :
hD ρ LV D 900 × 0.55 × 0.10
o
1
Equation 7.45 : Nu = = C Rem Pr 3 Re = = = 900 Pr = 715
kL µL 0.055
hD
o.c
1
⇒ C = 0.683 m = 0.466 ⇒ Nu = = 0.683 × 9000.466 × 715 3 = 145.4
kL
145.4 × 0.145
h= = 210.8
0.1
-pr
   0.4  2 / 3     Re 5/ 8 
−1/ 4 4/5
hD  0.5
1
Equation 7.47 : Nu = = 0.3 + 0.62 Re Pr × 1 + 
3
   × 1 +   
kL  
  Pr      282000  

   0.4  2 / 3     900 5 / 8 
−1/ 4 4/5
 0.5
1
Nu = 0.3 + 0.62 × 900 × 715 3 × 1 +     × 1 +    = 169.98
   715      282000  
ri
169.98 × 0.145
h= = 246.48
ka
0.1
ab
al3
w.
ww
TDEX- 8.1
m
Travaux dirigés et exercices du chapitre 8 : Convection interne
Ex-8.1 (S) Un tube, de 0.05 m de diamètre interne et de 1 m de longueur, est enveloppé d'une résistance
co
électrique chauffante dont la surface externe est parfaitement isolée. Cette résistance délivre une certaine
énergie à la paroi interne du tube dont la température TS est supposée uniforme. Un débit d'air de 0.001 kg/s
0 0
entre dans le tube à une température de 20 C et en ressort à 50 C.
a) quelle est la densité moyenne de flux de chaleur à la surface interne du tube ?
b) quel est le coefficient moyen de transfert de chaleur dans le tube (attention choisir la bonne
corrélation, table p 8.5.1)?
ro.
c) en utilisant la différence moyenne logarithmique de température, déterminer la température de la
surface.
-7
Propriétés de l'air (unités SI) cP=1007, µ=188. 10 , k=0.0269, Pr=.71
2 2 o
[réponses : 192.3 W/m -3.64W/m .K-89.3 C]
a)
b)
c)
ri-p
Ex-8.2 (S) Une huile s'écoule avec un débit de 0.5 kg/s dans un tube de 25 mm de diamètre et de 10 m de
0
longueur. La température de l'huile à l'entrée est de 25 C et on suppose que la température de la paroi interne
0
du tube reste constante à 100 C.
Quelle corrélation de la page 8.5.1 s'applique au calcul du coefficient h dans ce problème?
Calculer le coefficient de transfert de chaleur par convection à la surface du tube.
Calculer la température de l'huile en sortie ainsi que la chaleur totale échangée
-2
Propriétés de l'huile (unités SI) cP=2076, µ=5.31 10 , k=0.139, ρ=860
ka
2 o
[réponses : 8.56-92.5W/m .K-30.1 C-5293.8W]
Ex-8.3 Refaire le problème 2 en supposant qu'il s'agit d'écoulement complètement établi

[réponses : 8.55-20.3-26.1-1141.8]
-4
Ex-8.4 (S) Un débit d'air de 3.10 kg/s s'écoule dans un canal rectangulaire de dimensions 4 mm x 16
ab
2
mm x 1 m de longueur. Une densité de flux uniforme de 600 W/m est imposée sur la paroi interne du canal.
a) Quelle est la quantité totale de chaleur reçue par l'air?
0
b) Quelle est la température de l'air en sortie sachant qu'il entre à 27 C?
c) Déterminer le coefficient moyen de transfert de chaleur dans le canal
d) Déterminer la température de la paroi interne du canal à la sortie.
al3
-7
Propriétés de l'air (unités SI) cP=1007, µ=184.6 10 , k=0.0263, Pr=.707
Pour cette géométrie, en utilisant le diamètre hydraulique DH, on a Nu=hDH/k=5.33
o 2 o
[ réponses : 24W-106 C-22W/m .K-133 C]
Ex-8.5(S) Un fluide froid à une température Tm,i entre avec un débit massique m, dans un tube de
longueur L, dont la surface interne est à une température constante TS (chaude). Le diamètre du tube décroît
linéairement de D=a à l'entrée à D=b à la sortie [ D est une fonction de x, distance depuis l'entrée:
w.
D(x)=a+(b-a)(x/L)].
Pour cette géométrie, le coefficient de transfert de chaleur local à une distance x est donné par:
1.8
h(x)/h0=[b/D(x)]
a) déterminer l'équation différentielle reliant la variation de la température avec la distance x
b) montrer qu'en sortie, la température du fluide Tm,o vérifie la relation:
1.8 0.2 0.2
(TS-Tm,O)/(TS-Tm,I)=exp{-πhob L [b -a ] / [0.2 m cP(b-a)]}
ww
Remarque: La résolution des problèmes 1-4 est grandement facilitée par le fait que les valeurs des
propriétés physiques sont données dans l'énoncé. En pratique, il faudrait utiliser une méthode itérative pour
déterminer la température pertinente à laquelle on calcule ces propriétés.
TDEX- 8.2
m
Un tube, de 0.05 m de diamètre interne et de 1 m de longueur, est enveloppé d'une résistance électrique
chauffante dont la surface externe est parfaitement isolée. Cette résistance délivre une certaine énergie à la
paroi interne du tube dont la température TS est supposée uniforme. Un débit d'air de 0.001 kg/s entre dans le
co
0 0
tube à une température de 20 C et en ressort à 50 C.
a) quelle est la densité moyenne de flux de chaleur à la surface interne du tube ?
b) quel est le coefficient moyen de transfert de chaleur dans le tube (attention choisir la bonne
corrélation, table p 8.5.1)?
c) en utilisant la différence moyenne logarithmique de température, déterminer la température de la
ro.
surface.
-7
Propriétés de l'air (unités SI) cP=1007, µ=188. 10 , k=0.0269, Pr=.71
Solution :
a)
ri-p
qéchangée = m CP ∆T = 0.001 1007 ( 50 − 20 ) = 30.21
Surface d ' échange = S = π D L = π 0.05 1 = 0.1571
[W ]
 m 2 
qéchangée 30.21
q′′ = = = 192.3 W / m 2 
S 0.1571
ka
b)
ρVd πd2 4m 4 × 0.001
Re = m= ρV ⇒ρV = = = 0.5093
µ 4 π d 2 π × 0.052
ab
0.5093 × 0.05
Re = = 1354.1
188.10 −7
Pour Re<2300, on est donc en régime laminaire. Calculons les longueurs d’entrée hydrodynamique et
thermique (pages 8.1 et 8.2):
al3
entrée
X hydrodynamique
= 0.05 Re = 0.05 × 1354.5 = 67.72
D
⇒ X hydrodynamique
entrée
= 67.72 0.05 = 3.38 m
entrée
X thermique CP µ 1007 ×188 10−7
w.
= 0.05 Re Pr Pr = = = 0.704
D k 0.0269
entrée
X thermique
= 0.05 ×1354.5 × 0.704 = 47.66
D
ww
⇒ X thermique
entrée
= 47.66 × 0.05 = 2.38 m
Ces deux longueurs d’entrée sont plus grandes que la longueur de 1 m du tube; on est donc dans une
région combinée d’entrée (hydrodynamique+ thermique), et c’est l’équation 8.57 qui s’applique.
TDEX- 8.3
om
1 1
hD  Re Pr  3  µ 0.14  1354 × 0.704  3

Nu = = 1.86     ≈ 1.86   1 = 6.74
k  L   µS   1 
 D   0.05 
6.74 × 0.0269
h= = 3.63 W / ( m 2 . o C )  ou W / ( m 2 . K ) 
0.05
o.c
qéchangée q′′ 192.3
qéchangée = AS h ∆TLM ⇒ ∆TLM = = = = 52.97
AS h h 3.63
∆TLM =
(T S ) (
− Tin − TS − Tout )=
(Tout − Tin ) = ( 50 − 20 )
-pr
(T − T )
Ln S in
(T − T ) Ln (TS − Tin )
Ln S in
(TS − Tout ) (TS − Tout ) (TS − Tout )
(T − T ) 30 = 0.566 ⇒ (TS − Tin ) = e0.566 = 1.762
Ln S in =
(TS − Tout ) 52.97 (TS − Tout )
ri
(TS − Tin ) = (TS − Tout ) 1.762 ⇒ Ts (1.762 − 1) = 1.762 Tout − Tin
ka
1.762 Tout − Tin 1.762 × 50 − 20
Ts = = = 89.4 O C
(1.762 − 1) 0.762
ab
al3
w.
ww
TDEX- 8.4
m
Solutionnaire 8.2
a)
co
ρVd πd2 4m 4 0.5
Re = m= ρV ⇒ρV = 2
= = 1018
µ 4 πd π 0.0252
1018 × 0.025
Re = = 479.56 ⇒ régime laminaire
5.31 10−2
ro.
x
longueur d ' entrée hydrodynamique : E ,hydro = 0.05 Re = 23.97 ⇒ xE = 23.97 × 0.025 = 0.599 m
D
−2
C µ 2076 × 5.3110
Prandtl = Pr = P = = 793.
k 0.139
x
-p
longueur d ' entrée thermique : E ,ther = 0.05 Re Pr ⇒ xE ,ther = 0.05 × 479.56 × 793 × 0.025 = 475 m
D
Avec un tube de 10m de longueur, l’hydrodynamique est dans une zone développée (0.599 << 10m) alors
ari
que seule la thermique est en zone d’entrée (10<<475). On utilisera la corrélation 8.56 page 8.7 (laminar
thermal entry length).
b)
D
0.0668   Re Pr
k
Nu =
hD
= 3.66 + L = 16.71 ⇒ h = 92.9 W / m 2 .K
2/3
k  D  
ab
1 + 0.04   Re Pr 
 L  
c) Page 8.5, on a la relation (13)
As h
al3
md CP = −
 ∆T 
Ln  out 
 ∆Tin 
 ∆T  As h
⇒ Ln  out  = −
 ∆Tin  md CP
w.
As h
∆Tout −
⇒ =e md CP
∆Tin
ww
TDEX- 8.5
om
∆Tout 100 − Tout As h 0.785 × 92.9
= As = π DL = π × 0.025 ×10 = 0.785 m 2 = = 0.0702
∆Tin 100 − 25 md CP 0.5 × 2076
100 − Tout
= e −0.0702 = 0.932 ⇒ Tout = 100 − 75 × 0.932 = 30.10 C
100 − 25
o.c
chaleur échangée : Q = md CP (Thuile ,out − Thuile ,in ) = 0.5 × 2076 × ( 30.1 − 25 ) = 5293.8 watts
Solutionnaire 8.4
-pr
a) La surface d’échange vaut :
As = Périmètre × Longueur = 2 ( 0.004 + 0.016 ) × 1 = 0.04 m 2

Q = As q′′ = 0.04 × 600 = 24 W
ari
b)
chaleur échangée : Q = mair CP ,air (Tair ,out − Tair ,in )
Q 24
⇒ Tair ,out = + Tair ,in = −4
+ 27 = 106.40 C
mair C P ,air 3. 10 × 1007
k
c) Page 8.8, on trouve la corrélation pour un tube de section rectangulaire :

ab
b 16 hDH section d ' écoulement

= = 4 ⇒ Nu = = 5.33 ⇒ DH = 4 RH = 4
a 4 k périmètre mouillé
16 × 4 × 10−6 5.33 × 0.0263
DH = 4 = 6.4 10−3 m h= = 21.9 W / m 2 .K
2 × (16 + 4 ) × 10 −3
6.4 10 −3
al3
d) On applique la loi de refroidissement de Newton à la paroi du tube en sortie :
q′′ 600
q′′ = h (Ts − Tair ) ⇒ Ts = + Tair ⇒ Ts = + 106.4 = 133.8 o C
h 21.9
w.
ww
TDEX- 8.6
m
Un fluide froid à une température Tm,i entre avec un débit massique m, dans un tube de longueur L, dont la
surface interne est à une température constante TS (chaude). Le diamètre du tube décroît linéairement de D=a
à l'entrée à D=b à la sortie [ D est une fonction de x, distance depuis l'entrée: D(x)=a+(b-a)(x/L)]. Pour cette
1.8
géométrie, le coefficient de transfert de chaleur local à une distance x est donné par: h(x)/h0=[b/D(x)]
o
a) déterminer l'équation différentielle reliant la variation de la température avec la distance x
b) montrer qu'en sortie, la température du fluide Tm,o vérifie la relation:
1.8 0.2 0.2
(TS-Tm,O)/(TS-Tm,I)=exp{-•hob L [b -a ] / [0.2 m cP(b-a)]}
o.c
ri -pr
ka
ab
Solution :
a) En négligeant la conduction axiale, le bilan d’énergie s’écrit :
m CP (Tm , x − Tref ) − m CP (Tm , x − Tref ) + π Dx ∆x hx (TS − Tm , x ) = 0

x x +∆x
al3
T − T 
m CP  m, x m , x +∆x  + π Dx hx (TS − Tm , x ) = 0
 ∆x 
∂T
− m CP m , x + π Dx hx (TS − Tm , x ) = 0
∂x
∂Tm , x
w.
m CP = π Dx hx ∂x
(TS − Tm , x )
b)
ww
TDEX- 8.7
om
1.8
x  b 
Dx = a + ( b − a ) hx = h0  
L  Dx 
1.8
dTm, x d (TS − Tm , x )  b 
m CP = − m CP = π Dx h0   dx = π h0 b1.8 Dx−0.8 dx
(TS − Tm , x ) (TS − Tm, x )  Dx 
 − m CP  d (TS − Tm, x )
o.c
 1.8 
= Dx−0.8 dx
 π h0 b  (TS − Tm , x )
Dx = a + ( b − a )
x
⇒ dDx =
( b − a ) dx ⇒ dx = L dD
L L (b − a ) x
 m CP  d (TS − Tm , x ) L
-pr
− 1.8 
= Dx−0.8 dDx
 π h0 b  (TS − Tm , x ) ( b − a )
 ( b − a ) m CP  d (TS − Tm, x )
− 1.8 
= Dx−0.8 dDx
 L π h0 b  (TS − Tm , x )
 ( b − a ) m CP  x = L d (TS − Tm , x ) b
1 1
ari
0.2 b
− 1.8  ∫ = ∫ D −0.8
x dDx = 
 Dx
 = ( b 0.2 − a 0.2 )
 L π h0 b  x = o (TS − Tm, x ) 0.2 0.2
a
a
(TS − Tm , x = L )  L π h0 b1.8  0.2

 ( b − a )
0.2
Ln =  −
(TS − Tm , x =o )  0.2 ( b − a ) m CP 
bk
 L π h0 b1.8 
(TS − Tm,out ) −
 0.2( b − a ) m C (
 b0.2 − a 0.2 )
 
=e P
(TS − Tm ,in )
al 3a
w.
ww
ww
w.
al3
ab
kari
-pro.
co
m
TDEX- 9.1
m
Travaux dirigés et exercices du chapitre 9 : Échangeurs de chaleur
Ex-9.1 (S) Introduction au dimensionnement d’un échangeur

0 0 0
On souhaite refroidir 0.1 kg/s d'une huile de 100 C à 60 C. Pour cela, on utilise 0.2 kg/s d'une eau à 30 C. On
o
envisage d'utiliser un échangeur composé de deux tubes concentriques: le diamètre du tube interne, Di, est de
25 mm alors que celui du tube externe, DO est de 45 mm (on pourra négliger l'épaisseur du tube interne ainsi
que l'encrassement). Quelle est la longueur de tubes nécessaire pour refroidir l'huile sachant qu'elle s'écoule
o.c
dans l'espace annulaire et que l'on fonctionne à contre-courant ? Quelle serait cette longueur pour un
échangeur à co-courant
Données:
-6
EAU CPE = 4178 J/kg/K, µE= 725 10 Pa.s, kE=0.625 W/m.K, Pr=4.85
HUILE : CPH = 2131 J/kg/K, µH= 0.0325 Pa.s, kH=0.138 W/m.K
-pr
corrélation dans le tube intérieur:
0.8 0.4
Nu = (h D/k) = 0.023 Re Pr , avec Re=ρVD/µ, Pr=cP µ/k
corrélation dans l'espace annulaire:

Nu=(h D/k)= 5.56 avec D=diamètre hydraulique
Méthodologie:
ri
1) faire un schéma avec les températures
2) calculer la température de l'eau en sortie [réponse : 40.2]
3) calculer la LMTD [réponse : 43.2]
ka
4) calculer les coefficients de transfert de chaleur sur les surfaces interne et externe du tube interne.
[réponse : 2249.5 - 38.4]
5) calculer le coefficient global d'échange U, la surface d'échange nécessaire et la longueur requise.
[réponse :37.75-5.23-66.5
6) calculer la longueur pour un co-courant [réponse 72.3]
ab
Ex-9.2 (S) Méthode des NTU

2 2
1) Un échangeur à contre courant a une surface de 12.5 m et un coefficient global U=400 W/m /C. On
o
l'utilise pour refroidir un débit de 2 kg/s d'une huile à 100 C (CP,H=2000 J/kg/K) avec un débit de 0.48 kg/s
o
al3
d'une eau à TE=20 C (CP,E=4170 J/kg/K). Calculer les températures de sortie des deux fluides.
[réponse 86.6-66.7]
2) Résolution de 1) en utilisant une méthode itérative avec la LMTD et le Solveur d’Excel.

w.
ww
TDEX- 9.2
Ex-9.3 (S): Coefficient moyen d'échange dans un échangeur. [réf: A.P. Colburn Ind. Eng. Chem. 25 873
m
(1933)]
Considérons l’échangeur co-courant ci-dessous
.co
TC
1.
2.
Tc1
o
Tc2
dq
Th1
-pr Th2
ari
Th
Si on suppose que le coefficient global d'échange, U, varie linéairement avec la différence de température :
U − U1 ∆T - ∆T 1
= avec ∆T = Th − Tc [1]
bk
U 2 − U1 ∆T 2 - ∆T 1
Démontrer alors que :

q = U MOYEN A ∆T LM
3a
1 U 1 ∆T 2 - U 2 ∆T 1 ∆ -∆
avec U MOYEN = et ∆T LM = T 1 T 2
∆T LM ∆ ∆
Log U 1 T 2 Log T 1
U 2 ∆T 1 ∆T 2
Début de la démonstration :
.al
dd dq = -m h c Ph dT h = UdA( T h - T c ) dq = m cC Pc dT c = UdA( T h - T c )
U U
-dT h = ( T h - T c )dA dT c = ( T h - T c )dA
m hC Ph m cC Pc
1 1 1 1
w
d( T h - T c ) = -( + )U( T h - T c )dA d( ∆T) = -( + )U( ∆T)dA

m c C Pc m h C Ph m c C Pc m h C Ph
On remplace U par l'expression linéaire (1), on regroupe les ∆T dans le membre de droite et on intègre par
ww
rapport à ∆T. On rappelle que :

1 1 1
= −
x (a + x) a x a ( x + a)
TDEX- 9.3
m
Solutionnaire 9.1 Introduction au dimensionnement d’un échangeur
1)
co
ro.
in out
− Thuile
ri-p
2) la quantité de chaleur, Q, perdue par l’huile est reçue par l’eau:
Q = mhuileCP ,huile (Thuile ) = 0.1× 2131× (100 − 60 ) = 8524 W

Q 8524
Q = meau CP ,eau (Teau
out in
− Teau ) ⇒ Teauout = in
+ Teau = + 30 = 40.2o C
meau CP ,eau 0.2 × 4178
ka
∆T1 − ∆T2 (100 − 40.2 ) − ( 60 − 30 ) 29.8

3) ∆TLN = = = = 43.20 C
Ln
∆T1 (100 − 40.2 ) Ln
59.8
ab
Ln
∆T2 ( 60 − 30 ) 30
4) pour l’eau à l’intérieur du tube :

al3
ρ eauVeau D π D2 4meau
Nu = 0.023 Re0.8 Pr 0.4 Re = meau = ρ eauVeau ⇒ ρ eauVeau =
µeau 4 π D2
4meau D 4 meau 4 × 0.2
Re = = = = 14049.5
π D µeau π Dµeau π × 0.025 × 725 10−6
2
heau D 0.8 0.4 Nu keau 89.98 × 0.625

= 0.023 (14049.5 ) ( 4.85 ) = 89.98 ⇒ heau = = 2249.5W / m 2 .K
w.
Nu = =
keau D 0.025
pour l’huile dans l’espace annulaire (conduite non circulaire), on va utiliser le rayon hydraulique RH
ww
TDEX- 9.4
om
sec tion d ' écoulement
RH =
périmètre mouillé
π
(D 2
2 − D12 )
RH = 4 ( D − D1 )( D2 + D1 )
= 2
π ( D2 + D1 ) 4 ( D2 + D1 )
o.c
RH =
( D2 − D1 )
4
DH = 4 RH = D2 − D1 = 45 − 25 = 20 mm
-pr
hhuile DH Nu khuile 5.56 × 0.138
Nu = 5.56 = ⇒ hhuile = = = 38.4 W / m 2 .K
khuile DH 0.020
5) Comme on néglige la conduction dans l’épaisseur du tube interne et l’encrassement sur les deux faces
du tube interne, la résistance totale au transfert la somme des résistances de convection :
ari
1 1 1 1 1
= + = + = 0.02648 ⇒ U = 37.75
U heau hhuile 2249.5 38.4
Q 8524
Q = UA∆TLN ⇒ A= = = 5.226 m 2
U ∆TLN 37.75 × 43.2
bk
A 5.226
A = π D1 L ⇒ L = = = 66.54 m
π D1 π × 0.025
3a
6) dans le cas d’un co-courant, ∆TLN est alors donné par :

∆T1 − ∆T2 (100 − 30 ) − ( 60 − 40.2 ) 50.2
∆TLN = = = = 39.750 C
Ln
∆T1
Ln
(100 − 30 ) Ln
70
∆T2 ( 60 − 40.2 ) 19.8
al
Q 8524
A= = = 5.68 m 2
U ∆TLN 37.75 × 39.75
w.
A 5.68
A = π D1 L ⇒ L = = = 72.3 m
π D1 π × 0.025
ww
Rappel : le contre-courant est plus efficace que le co-courant

TDEX- 9.5
m
Solutionnaire 9.2 Méthode des NTU
Ceau = meau CP ,eau = 0.48 × 4170 = 2001.6 Chuile = mhuileCP ,huile = 2 × 2000 = 4000
Cmin 2001.6
co
Cmin = min ( Ceau , Chuile ) = Ceau = 2001.6 = = 0.5
Cmax 4000
UA 400 × 12.5
NTU = = = 2.5
Cmin 2001.6
Sur la figure 11.15 page 9.30, on trouve, ε ≈ 0.83
ro.
o i
Cfr ( Tfro - Tfri ) Ceau ( Teau - Teau ) ε Cmin
ε = i
Cmin ( Tch - Tfr )i
= i i
Cmin ( Thuile - Teau )
⇒ Teau
o
=
Ceau
( i
Thuile i
- Teau ) i
+Teau
0.83 × 2001.6
o
Teau = ( 100 - 20 ) +20=86.4 o C
2001.6
ε = ch chi
ri-p
C ( T i - Tcho )
i
Cmin ( Tch - Tfr )
⇒ T o
huile
0.83 × 2001.6
= T i
huile −
i
ε C min ( Thuile
C huile
i
- Teau )
o
Thuile = 100 − ( 100 - 20 ) = 66.7 o C
4000
ka
Solutionnaire 9.3 Coefficient U moyen
dq = -m h c Ph dT h = UdA( T h - T c ) dq = m cC Pc dT c = UdA( T h - T c )
ab
U U
-dT h = ( T h - T c )dA dT c = ( T h - T c )dA
m hC Ph m cC Pc
1 1 1 1
d( T h - T c )= -( + )U( T h - T c )dA d( ∆T)= -( + )U( ∆T)dA
m c C Pc m h C Ph m c C Pc m h C Ph
al3
U − U1 ∆T - ∆T 1
= avec ∆T = Th − Tc
U 2 − U1 ∆T 2 - ∆T 1
∆T - ∆T 1
U = U1 + (U 2 − U1 )
∆T 2 - ∆T 1
1 1 1 1  ∆T - ∆T 1 
) U1 + (U 2 − U1 )
w.
d( ∆T)= -( + )U( ∆T)dA = -( + ∆TdA

m c C Pc m h C Ph m c C Pc m h C Ph  ∆T 2 - ∆T 1 
posons α = -(
1
+
1
) β=
(U 2 − U1 )
ww
m c C Pc m h C Ph ∆T 2 - ∆T 1
 ∆T - ∆T 1   (U − U1 ) ∆T - ∆  ∆TdA
d( ∆T) = α U1 + (U 2 − U1 )  ∆TdA = α U1 + 2 ( T 1) 
 ∆T 2 - ∆T 1   ∆T 2 - ∆T 1 
TDEX- 9.6
m
d ( ∆T ) = α U1 + β ( ∆T - ∆T 1 )  ∆TdA = α [U1 − β ∆T 1 + β ∆T ] ∆TdA
d ( ∆T ) U 
= α dA posons : a =  1 − ∆T 1 
[U1 − β∆T 1 + β∆T ] ∆T
co
β 
1 1 1 1 1
= =
[U1 − β∆T 1 + β∆T ] ∆T β  U1 − ∆  + ∆T  ∆T β [ a + ∆T ] ∆T
 T 1 
 β  
ro.
1 1 1
mais = −
[ a + ∆T ] ∆T a∆T a [ a + ∆T ]
d ( ∆T ) 1  d ( ∆T ) d ( ∆T ) 
=  −  = α dA
[U1 − β∆T 1 + β∆T ] ∆T β  a∆T a [ a + ∆T ] 
∆T2
∫
∆T1
d ( ∆T )
a∆T
∆T
− ∫
2
∆T1
a
d ( ∆T )
[ a + ∆T ] AS
ri-p
1
a
∆T 1
= ∫ αβ dA ⇒ Ln 2 − Ln
∆T1 a
a + ∆T2
a + ∆T1
= αβ AS
∆T2 a + ∆T2
Ln − Ln = aαβ AS
∆T1 a + ∆T1
ka
U  U1 U ( ∆T 2 - ∆T 1) − (U 2 − U1 ) ∆T1 + (U 2 − U1 ) ∆T2
a + ∆T2 =  1 − ∆T 1  + ∆T2 = − ∆T 1 + ∆T2 = 1
β  (U 2 − U1 ) (U 2 − U1 )
ab
∆T 2 - ∆T 1
U ( ∆T2 − ∆T1 )
a + ∆T2 = 2
(U 2 − U1 )
al3
U  U1 U ( ∆T 2 - ∆T 1) − (U 2 − U1 ) ∆T1 + (U 2 − U1 ) ∆T1
a + ∆T1 =  1 − ∆T 1  + ∆T1 = − ∆T 1 + ∆T1 = 1
β  (U 2 − U1 ) (U 2 − U1 )
∆T 2 - ∆T 1
U ( ∆T2 − ∆T1 )
a + ∆T1 = 1
w.
(U 2 − U1 )
U 2 ( ∆T2 − ∆T1 )
a + ∆T2 (U 2 − U1 ) = Ln U 2
ww
Ln = Ln
a + ∆T1 U1 ( ∆T2 − ∆T1 ) U1
(U 2 − U1 )
TDEX- 9.7
m
∆T2 a + ∆T2 ∆T U U ∆T
Ln − Ln = Ln 2 − Ln 2 = Ln 1 2 = aαβ AS
∆T1 a + ∆T1 ∆T1 U1 U 2 ∆T1
co
 1 1 
α = − + 
 m c C Pc m h C Ph 
ro.
q q
q = m c C Pc (Tc 2 − Tc1 ) ⇒ (Tc 2 − Tc1 ) = q = m h C Ph (Th1 − Th 2 ) ⇒ (Th1 − Th 2 ) =
m c C Pc h C Ph
q q
(Tc 2 − Tc1 ) + (Th1 − Th 2 ) = + = −α q (Tc 2 − Tc1 ) + (Th1 − Th 2 ) = (Th1 − Tc1 ) − (Th 2 − Tc 2 ) = ∆T1 − ∆T2
m c C Pc m h C Ph
∆T1 − ∆T2 = −α q ⇒ α =
Ln
U1∆T2
ri-p
∆T − ∆T1
= aαβ AS = 2
∆T2 − ∆T1
q
a β AS ⇒ q =
∆T2 − ∆T1  U1

U1∆T2  β
 (U − U1 )
− ∆T 1  2 AS
U 2 ∆T1 q Ln  ∆T 2 - ∆T 1
U 2 ∆T1
ka
 
1  U1  1  U1 ( ∆T 2 - ∆T 1) ∆T 1 (U 2 − U1 ) 
q=  − ∆T 1  (U 2 − U1 ) AS =  −  (U 2 − U1 ) AS
Ln 1 2  ( 2
U ∆T U − U1 )  U1∆T2  U 2 − U1 U 2 − U1 
 Ln
U 2 ∆T1  ∆T 2 - ∆T 1  U 2 ∆T1
ab
1 U1 ( ∆T 2 - ∆T 1 ) − ∆T 1 (U 2 − U1 )  1
q=
U1∆T2    (U 2 − U1 ) AS = (U ∆T − ∆T 1U 2 ) AS
U1∆T2 1 2
Ln U 2 − U 1  Ln
U 2 ∆T1 U 2 ∆T1
 
al3
 (U ∆T 2 − U ∆T 1 ) 
q= 1 2
 AS = U moyen AS ∆TLN
 Ln 1 2 U ∆T 
 U 2 ∆T1 
(U1∆T 2 − U 2 ∆T 1) ( ∆T 2 − ∆T 1)
w.
q = U moyen AS ∆TLN avec U moyen = et ∆TLN =

U1∆T2 ∆T2
∆TLN Ln Ln
U 2 ∆T1 ∆T1
ww
ww
w.
al3
ab
kari
-pro.
co
m
TDEX- 10.1
m
Travaux dirigés et exercices du chapitre 10 : Convection naturelle
co
Ex-10.1 (S) Convection naturelle à la surface d’un mur
Déterminer les coefficients de convection moyen à la surface d'un mur (vertical !) de 2.5 m de haut à l'intérieur d'une
maison pour les conditions suivantes:
a) Température de la surface=100C température de l'air ambiant=200C
ro.
données pour l'air en SI: ß=1/Tf(K), <=14.82 10-6, k=0.0253, α=20.9 10-6, Pr=0.71
b) Température de la surface=370C température de l'air ambiant=270C

données pour l'air en SI: ß=1/Tf(K), <=16.39 10-6, k=0.0267, α=23.2 10-6, Pr=0.706
(réponses: a: h=3.03 W/m2/K b: h=2.94 W/m2/K)
ri-p
Ex-10.2 (S) Convection naturelle à la surface d’un cylindre
Un cylindre horizontal de 5 mm de diamètre est immergé dans de l'eau à 180C. Si la température de la surface du cylindre
est de 560C, déterminer les pertes de chaleur due à la convection naturelle par mètre de longueur de cylindre
données pour l'eau en SI: ρ= 993., ß=361.2 10-6, <=6.999 10-7, k=0.628, α=1.514 10-7, Pr=4.62,
ka
( réponses : h=1312 W/m2/K et q=783 W/m)
Ex-10.3 (S) Convection naturelle à la surface d’une sphère

ab
Une sphère de 25 mm de diamètre contient une résistance électrique. Déterminer la puissance nécessaire pour maintenir
la température à la surface de la sphère à 94 0C si cette sphère est immergée dans les fluides suivants:
a) air à 20 0C: <=18.91 10-6, k=0.0285, α=26.9 10-6, Pr=0.711, ß=3.03 10-3
b) eau à 20 0C: <=0.497 10-6, k=0.650 , α=0.158 10-6, Pr=3.15, ß=0.504 10-3
al3
c) glycol à 200C: <=5.15 10-6, k=0.260 , α=0.936 10-6, Pr=55., ß=0.65 10-3
(réponses: a: h=10.6 q=1.55 W b: h=1297 q=188 W c: h=230 q=33.4 w)

w.
Ex-10.4 (S) Perte de chaleur par convection naturelle à la surface d’un cylindre
De la vapeur circule dans un tuyau horizontal de 0.3048 m de diamètre. La température de la surface du tuyau est de
2500C alors que la température du milieu environnant est de 150C. Calculer les pertes de chaleur (par mètre de tuyau) due
à la convection naturelle.
ww
Propriétés de l'air
( à Tf=[250+15]/2=132.5 0C) k=0.034 W/ (m.K) β=2.47 10-3 K-1 <=26.26 10-6 m2/s Pr=0.687
réponses : h=7.27 W/m2.K q=1.6 kW /m

TDEX- 10.2
om
Solutionnaire 10.1 :
a)
1 1 1
β= = = = 3.47 × 10−3
Tf 273 + (10 + 20) 288
2
o.c
3
g β (T S - T ∞) L 3 9.81× 3.47 × 10 × ( 293 − 283) × 2.5
−3
Ra = Gr Pr = = -6 −6
=1.717 × 1010
να 14.82 × 10 × 20.9 × 10
2
 
 
 
-pr
1
 0.387 Ra 6  hL 299.4 × 0.0253
Nu L = 0.825 + 8  = 299.4 = ⇒h= = 3.03 W / m 2 .K
 k 2.5
  0.492  16  
9 27
 1 +    
   Pr   

ari
b)
1 1 1
k
β= = = = 3.28 × 10−3
Tf 273 + (27 + 37) 305
2
ab
3
g β (T S - T ∞) L 3 9.81× 3.28 × 10−3 × 10 × 2.5
Ra = Gr Pr = = -6 −6
=1.322 × 1010
να 16.39 × 10 × 23.2 × 10
2
 
al3
 
 1 
 0.387 Ra 6  hL 275.4 × 0.0267
Nu L = 0.825 + 8  = 275.4 = ⇒h= = 2.94 W / m2 .K
 k 2.5
  0.492  16  
9 27
 1 +    
   Pr   

w.
ww
TDEX- 10.3
om
3
g β (T S - T ∞) D 3 9.81× 361.2 × 10−6 × ( 56 − 18 ) × .005
Ra = Gr Pr = = −7 −7
=1.588 × 105
να 6.999 × 10 × 1.514 × 10
2
 
 
o.c
 1 
 0.387 Ra 6  hD
Nu D =  0.60 + 8  = 10.45 =
   k
 0.559  16 
9 27
 1 +    
   Pr   
10.45 × 0.628
-pr
⇒h= = 1312 W .m2 .K q = π DLh (TS − T f ) = π × 0.005 × 13121× 1× (56 − 18) = 783 W / m
0.005
ri
a)
3
g β (T S - T ∞) D 3 9.81× 3.03 × 10 × ( 94 − 20 ) × .025
−3
ka
Ra = Gr Pr = = =67564
να 18.91× 10−6 × 26.9 × 10−6
1
0.589 Ra 4
9.33 × 0.0285
Nu D = 2 + 4
= 9.33 ⇒ h= =10.64 W/m 2 .K
ab
 0.025
 0.469  16 
9 9
1 +   
  Pr  
2
 0.025 
q = 4π R h (TS − T∞ ) = 4π × 
2
 × 10.64 × ( 94 − 20 ) = 1.546 W
 2 
al3
b)
3
g β (T S - T ∞) D 3 9.81× 0.504 × 10−3 × ( 94 − 20 ) × .025
Ra = Gr Pr = = −6 −6
=7.373 ×107
να 0.497 ×10 × 0.156 × 10
1
49.9 × 0.650
w.
4
0.589 Ra
Nu D = 2 + 4
= 49.9 ⇒ h= =1297.4 W/m 2 .K
 0.025
 0.469  16 
9 9
1 +   
  Pr  
ww
2
 0.025 
q = 4π R 2 h (TS − T∞ ) = 4π ×   × 1297.4 × ( 94 − 20 ) = 188.5 W
 2 
TDEX- 10.4
m
c)
3
g β (T S - T ∞) D 3 9.81× 0.650 × 10−3 × ( 94 − 20 ) × .025
= 1.529 × 106
co
Ra = Gr Pr = = −6 −6
να 5.15 × 10 × 0.936 ×10
1
0.589 Ra 4
22.11× 0.260
Nu D = 2 + 4
= 22.11 ⇒ h= =229.98 W/m 2 .K
 0.025
 0.469  16 
ro.
9 9
1 +   
  Pr  
2
 0.025 
q = 4π R 2 h (TS − T∞ ) = 4π ×   ×188.5 × ( 94 − 20 ) = 33.4 W
 2 
Solutionnaire 10.4
g β (T S - T ∞) D 3
-p
9.81× 2.47 × 10−3 × ( 250 − 15 ) × 0.3048
3
ari
Gr = = = 2.338 × 108
−6 2
ν 2
( 26.26 ×10 )
Ra = Gr Pr = 2.338 × 108 × 0.687 = 1.606 × 108
k
 
 
 1 
ab
 0.387 Ra 6  65.15 × 0.034

Nu D =  0.60 + 8  = 65.15 ⇒ h = = 7.27 W / m 2 K
   0.3048 =
 0.559  16 
9 27
 1 +    
   Pr   
al3
q = π DLh (TS − T∞ ) = π × 0.3048 × 1× 7.27 × ( 250 − 15 ) = 1635 = 1.635 kW / m

w.
ww
TDEX- 11.1
om
Travaux dirigés et exercices du chapitre 11: Ébullition et condensation
Ébullition nucléée en piscine
Ex-11.1 (S) Casserole
o.c
Une casserole en cuivre, de 150 mm de diamètre, est utilisée pour faire bouillir de l'eau. Le fond de la
O
casserole est maintenu à une température de TW =115 C.
a) Déterminer la densité de flux à la paroi (comparer votre résultat avec la figure p 10.2.1) et la
5 2
puissance disponible pour faire bouillir l'eau. [réponses : 4.6 10 W/m , 8.16 kW]
b) Quelle est la quantité d'eau évaporée chaque minute ? [réponse : 216g/min.]
c) Quelle est la valeur du rapport du flux réel au flux critique? [réponse : 0.367]
3 3 -4 2
-pr
Données: ρL=957.9 kg/m ; ρV=0.5955 kg/m ; µL=2.79 10 Pa.s, g=9.81 m/s
6
∆hLV=2.257 10 J/kg= chaleur latente de vaporisation.
-3
σ=58.9 10 N/m; PrL=1.76; CPL=4217 J/kg
CSf=0.013 et n=1 (équations p 11.4)
Ex-11.2 (S): Tubes

Des tubes de cuivre de 25 mm de diamètre et de 0.75 m de longueur sont utilisés pour l=ébullition en piscine
ari
2
sous 1 atm. Sachant que les tubes opèrent à 75% du flux critique (utilisée la valeur de 1.26 MW/m trouvée à
l=exercice 1), combien de tubes sont nécessaires pour permettre une vitesse d=évaporation de 750 kg/h ?
O
Quelle est la température à la surface des tubes ? [réponses : 8.45, 117.6 C]
3 3 -4 2
6
bk
-3
σ=58.9 10 N/m; PrL=1.76; CPL=4217 J/kg
Condensation:
Ex-11.3 (S): Condensation en film laminaire sur un cylindre vertical

3a
0
La paroi extérieure d'un cylindre vertical de 80 mm de diamètre et de 1 m de longueur est maintenue à 50 C
grâce un fluide circulant à l'intérieur du cylindre. Ce cylindre est environné d'une vapeur d'eau saturée sous
une pression de 1 atm. Déterminer le coefficient moyen de transfert de chaleur sur la surface, la quantité
totale de chaleur transférée ainsi que la masse de vapeur se condensant sur le cylindre.
3 3 -6 2
Données: ρV=0.596 kg/m ; ρL=975 kg/m ; µL=375 10 Pa.s, g=9.81 m/s
al
0 6
kL=0.668 w/m; TSAT=100 C; ∆hLV=2.257 10 J/kg= chaleur latente
CPL=4193 J/kg (utiliser l'équation p 11.7 pour une plaque verticale)
2
[réponses : 4094W/m /K; 51446 W; 21.4 g/s]
w.
Ex-11.4 (S) Condensation en film laminaire sur un cylindre horizontal

O
Un cylindre horizontal de 50 mm de diamètre dont la température de la surface est de 34 C, est exposé à une
vapeur sous 0.2 bar. Calculer la masse de vapeur, MV,qui se condense sur le tube et la chaleur échangée ,Q,
par m de longueur de tube.
ww
2
[réponses : h=6927 W/m .K, 11.6 g/s.m, 28300 W/m]
Données:
3 3 -6 2
ρV=0.129 kg/m ; ρL=989 kg/m ; µL=577 10 Pa.s, g=9.81 m/s
kL=0.640 w/m; TSAT=333K, ∆hLV=2358 KJ/kg= chaleur latente
CPL=4180 J/kg (utiliser l'équation p 11.7 pour le cylindre)
TDEX- 11.2
om
Ex-11.5: Démonstration pour l==obtention d==un coefficient de condensation en film sur une plaque
En 1916, Nusselt a proposé une analyse du phénomène de condensation en film sur une paroi froide verticale.
En utilisant les suggestions ci-dessous retrouvez l'analyse de Nusselt et obtenez une expression du coefficient
de transfert de chaleur pour la condensation film.
o.c
- Faites un schéma d'une plaque verticale, de largeur W, à la température TW , sur laquelle se
condense une vapeur qui est à sa température de saturation TSAT. Au fur et à mesure que la vapeur
se condense un film de liquide se developpe et s'écoule le long de la plaque. (on appellera x l'axe
vertical et y l'axe horizontal)
- On suppose que l'épaisseur du film, δ(x), varie lentement le long de la plaque et que le profil de
vitesse Vx(y) dans le film est donné par la relation:
-pr
( ρ L - ρ V )gδy y
V x (y) = [1 - ]
µL 2δ
- Déterminez l'expression du débit massique de condensat, Mc, à la cote x

3
ari
[réponse: Mc= W ρ L(ρL-ρV) g δ /(3 µL)]
- Quelle est la variation de ce débit entre les cotes x et x+dx en fonction de la variation de l=épaisseur
du film ?
- On suppose que le transfert de chaleur dans le film liquide est purement conductif et accepte
l'hypothèse du régime quasi-permanent. Trouver l'expression du profil de température et du flux de
chaleur en paroi (à x donné).
k
[réponse (T-Tw) /(TSAT-TW )=y/ δ]

- Trouvez l'expression du flux à la paroi [réponse : q@=(TSAT-TW ) (kL/δ)]
ab
- En considérant que le flux de chaleur à la paroi correspond à la chaleur nécessaire pour condenser la
vapeur (on négligera la chaleur sensible), trouver l'équation différentielle reliant la variation de
l'épaisseur du film avec la cote x et obtenez l=expression de la variation δ(x)
4x k L( T sat - T W )µ L
al3
Réponse: δ 4 =
∆ h LV g ρ L( ρ L - ρ V )
- Trouver les relations pour le coefficient de transfert de chaleur local h(x) et pour le coefficient h moyen
(pour 0< x<L)
Réponse:
w.
∆ h LV g ρ L( ρ L - ρ V )k 3L 1/4 ∆ h LV g ρ L( ρ L - ρ V )k 3L 1/4
hx = [ ] et h Moyen = 0.943[ ]
4x( T sat - T W )µ L L( T sat - T W )µ L
ww
TDEX- 11.3
m
Solutionnaire Ex-11.1 Casserole
Une casserole en cuivre, de 150 mm de diamètre, est utilisée pour faire bouillir de l'eau. Le fond de la
O
casserole est maintenu à une température de TW =115 C.
a) Déterminer la densité de flux à la paroi (comparer votre résultat avec la figure p 10.2.1) et la
co
puissance disponible pour faire bouillir l'eau.
b) Quelle est la quantité d'eau évaporée chaque minute ?
c) Quelle est la valeur du rapport du flux réel au flux critique?
3 3 -4 2
6
ro.
-3
σ=58.9 10 N/m; PrL=1.76; CPL=4217 J/kg
Solution :
a) La corrélation de Rohsenow donne la valeur du flux transféré :
"
q S = µ L∆ h LV [
"
q S = 2.79 10−4
g( ρ L - ρ V )
σ
2.257 106 [
ri-p
c P,L( T S - T SAT ) 3
1/2
] [
C sf ∆ h LV Pr L
n
]
9.81 ( 957.9 - 0.5955 ) 1/2

] [
4217( 115 − 100 )
]
3
0.0589 0.013 2.257 106 1.761

"
q S = 462254 W/m 2 ≈ 4.622 105 W/m 2
ka
"  0.152 
La puissance disponible est donc q = q S As = 4.622 105  π  = 8168.7 W = 8.169 kW
 4 
b
b) à chaque minute, l’énergie fournie permet l’évaporation d’une masse d’eau m telle que :
q ∆t
3a
q ∆t = m ∆ h LV ⇒ m =
∆ h LV
8.169 60
m= = 0.217 kg / min
2.257 106
.al
c) Le flux critique est donné par :

1
"  σ g (ρ L - ρV )  4
q C = 0.149 ∆ h LV ρ V  
 ρ V2 
1
w
" 6  0.0589 9.81 (957.9 − 0.5955)  4

q = 0.149 2.257 10
C
0.5955  
 0.59552
"
q C =1258540.8
ww
"
qS 462254
"
= = 0.3673
q C 1258540
TDEX- 11.4
m
a)
"
q S = 0.75 × qC′′ = 0.75 × 1.26 × 106 = 9.45 ×105 W / m 2
mévaporée × ∆hLV
co
⇒ nt =
"
qtotale = nt × π DL × q S = mévaporée × ∆hLV
π DL × q S"
750
× 2.257 × 106
nt = 3600 = 8.45
π × .025 × 0.75 × 9.45 × 105
ro.
b)
1 3
"  g (ρ L - ρV )  2  C P,L ( T S - T SAT ) 
q = µ L ∆ h LV 
S   n 
 σ   C sf ∆ h LV Pr L 
⇒ T S = T SAT +
C sf ∆ h LV
C P,L

n 
Pr L 

ri-p "
qS
1 
 µ L ∆ h LV  g ( ρ L - ρ V )  
2



1/ 3
  σ  
  
1/ 3
ka
 
 
0.013 × 2.257 × 106 × 1.761  9.45 × 105 
T S = 100 + × 1 
4217   9.81 × ( 957.9 - 0.5955 )  2

−4 6
 2.79 × 10 × 2.257 × 10  −3  
58.9 × 10
ab
   
1/ 3
 9.45 ×105 
T S = 100 + 11.304 ×   = 117.57 o C
 251442 
al3
w.
Solutionnaire Ex-11.3Condensation en film laminaire sur un cylindre vertical

0
La paroi extérieure d'un cylindre vertical de 80 mm de diamètre et de 1 m de longueur est maintenue à 50 C
grâce un fluide circulant à l'intérieur du cylindre. Ce cylindre est environné d'une vapeur d'eau saturée sous
une pression de 1 atm. Déterminer le coefficient moyen de transfert de chaleur sur la surface, la quantité
ww
totale de chaleur transférée ainsi que la masse de vapeur se condensant sur le cylindre.
3 3 -6 2
Données: ρV=0.596 kg/m ; ρL=975 kg/m ; µL=375 10 Pa.s, g=9.81 m/s
0 6
kL=0.668 w/m; TSAT=100 C; ∆hLV=2.257 10 J/kg= chaleur latente
CPL=4193 J/kg (utiliser l'équation p 11.7 pour une plaque verticale)
TDEX- 11.5
m
Solution :
a) Le cylindre étant vertical, on néglige la courbure et on utilise la corrélation pour un plan vertical :
co
1
h L  ρ g ( ρ L - ρ V ) ∆ h LV
′
L 
3 4
Nu = L = 0.943  L 
kL  µ L k L (T SAT - T S ) 
∆ h′ LV = ∆ h LV + 0.68 CP (T SAT - T S ) = 2.257 106 + 0.68 4193 (100 − 50) = 2399562
ro.
∆ h′ LV = 2399562 = 2.399 106
1
 975 9.81 (975-0.596) 2.399 106 13  4
Nu = 0.943   = 6128
 375 10−6 0.668(100-50) 
Nu k L 6128 0.668
-p
hL = = = 4094 W / m2C
L 1
b) la quantité totale de chaleur échangée s’exprime par :
ari
q = h L A (T SAT - T S ) = 4094 (π 0.080 1) (100 − 50) = 51446 W
c)
q = m ∆ h′ LV
k
q 51446
m= = = 0.0214 kg / s = 21.4 g / s
∆ h′ LV 2.399 106
ab
1
 ρ L g ( ρ L - ρ V ) ∆ h LV
′
kL 
3 4
h D = 0.729  
al3
 µ L D (T SAT - T S ) 
∆ h LV = ∆ h LV + 0.68 C PL ( T SAT - T S ) = 2358 × 103 + 0.68 × 4180 ×  333 − ( 273 + 34 )  = 2.4319 × 106
′
1
 989 × 9.81× ( 989 - 0.129 ) × 2.4319 × 103 × 0.643  4
h D = 0.729   = 6927.4 W / m 2o C
 577 × 10−6 × 0.05 ×  333 − ( 273 + 34 )  
w.
q π DL h D (T SAT - T S ) π × 0.05 × 1× 6927.4 × 333 − ( 273 + 34 ) 

MV = = =
∆ h′ LV ∆ h′ LV 2.4319 × 103
M V = 0.0116 kg / s.m = 11.6 g / s.m
ww
q = π DL h D (T SAT - T S )= M V × ∆ h′ LV = 0.0116 × 2.4319 × 103 = 28291 ≈ 28.3 kW / m detube

TDEX- 11.6
m
Demonstration pour l==obtention d==un coefficient de condensation en film sur une plaque
En 1916, Nusselt a proposé une analyse du phénomène de condensation en film sur une paroi froide verticale.
co
En utilisant les suggestions ci-dessous retrouvez l'analyse de Nusselt et obtenez une expression du coefficient
de transfert de chaleur pour la condensation film.
ro.
a) Faites un schéma d'une plaque verticale, de largeur W, à la
température TW , sur laquelle se condense une vapeur qui est à sa
température de saturation TSAT. Au fur et à mesure que la vapeur se
condense un film de liquide se développe et s'écoule le long de la plaque.
(on appellera x l'axe vertical et y l'axe horizontal
ri-p
Réponse ci-contre
ka
b) On suppose que l'épaisseur du film, •(x), varie lentement le long de la
plaque et que le profil de vitesse Vx(y) dans le film est donné par la
relation:
ab
( ρ L - ρ V )gδy y
V x (y) = [1 - ]
µL 2δ
al3
Déterminez l'expression du débit massique de condensat, Mc, à la cote x
Réponse :
Le débit s’obtient par intégration du profil de vitesse à une position x donnée :

w.
ww
TDEX- 11.7
m
δ δ
( ρ L - ρV ) g δ y  y 
M C = ∫ ρ L Vx ( y ) W dy = ∫ ρ L 1-  W dy =
0 0 µL  2δ 
δ δ
ρ (ρ -ρ )gδ W  y  ρL ( ρ L - ρ V ) g δW  y 2 1 y3 
∫0  2δ 
co
MC = L L V y 1- dy =  2 − 2δ 3 
µL µL  0
ρ L ( ρ L - ρ V ) g δ W  δ 2 δ 3  ρ L ( ρ L - ρ V ) g δ W  3δ 2 δ 2 
MC =  2 − 6δ  =  6 − 6
µL   µL  
ro.
ρ ( ρ - ρ ) g δ W  2δ  ρ L ( ρ L - ρ V ) g W δ
2 3
MC = L L V  6 =
µL   3µ L
c) Quelle est la variation de ce débit entre les cotes x et x+dx en fonction de la variation de l=épaisseur du film
?
Réponse :
ri-p
 ρL ( ρ L - ρ V ) g W δ 3 
d 
ρ L ( ρ L - ρ V ) g W d (δ )
3
d MC  3µ L 
= =
dx dx 3µ L dx
ka
d M C ρL ( ρ L - ρ V ) g W 2 dδ ρL ( ρ L - ρ V ) g W 2 dδ
= 3δ = δ
dx 3µ L dx µL dx
ρL ( ρ L - ρ V ) g W 2
d MC = δ dδ
µL
b
3a
d) On suppose que le transfert de chaleur dans le film liquide est purement conductif et accepte l'hypothèse du
régime quasi-permanent. Trouver l'expression du profil de température et du flux de chaleur en paroi (à x
donné).
Réponse :
.al
On résout le problème de la conduction suivant la direction y dans un film liquide stagnant

d’épaisseur δ avec comme conditions frontières
à y=0 T=TW et à y= δ T=Tsat
Le profil est donc linéaire (voir chapitre 3) et s’exprime par :

w
Ty − TW y
=
Tsat − TW δ
ww
∂T 1
q′′y = −k = − k (Tsat − TW )
y =0 ∂y y =0
δ
TDEX- 11.8
m
Commentaire : ce flux est négatif, la chaleur de condensation entre dans la paroi.
e) En considérant que le flux de chaleur à la paroi correspond à la chaleur nécessaire pour condenser la
vapeur (on négligera la chaleur sensible), trouver l'équation différentielle reliant la variation de l'épaisseur du
film avec la cote x et obtenez l=expression de la variation δ (x)
co
Réponse:
Entre les cotes x et x+dx , l’apport d’une quantité de condensat dMC a nécessité l’échange
à la surface Wdx de la paroi de la chaleur latente correspondante soit:
ro.
∆H LV d M C = q′′y y =0 Wdx
ρL( ρ L - ρ V ) g W 2 1
∆H LV δ dδ = k L (Tsat − TW ) Wdx
µL δ
ri-p
∆H LV ρ L ( ρ L - ρ V ) g 3
k L (Tsat − TW ) µ L
∆H LV ρ L ( ρ L - ρ V ) g
δ dδ = dx
δ x
∫ δ dδ = ∫ dx
3
k L (Tsat − TW ) µ L 0 0
ka
4
∆H LV ρ L ( ρ L - ρ V ) g δ
=x
k L (Tsat − TW ) µ L 4
4 x k L ( T sat - T W ) µ L
δ4=
∆ H LV ρ L ( ρ L - ρ V ) g
b
1
 4 x k L ( T sat - T W ) µ L  4
δ = 
 ∆ H LV ρ L ( ρ L - ρ V ) g 
3a
f) Trouver les relations pour le coefficient de transfert de chaleur local h(x) et pour le coefficient h moyen (pour
0< x<L)
.al
Réponse :Le flux à la paroi obtenu en d) s’exprime aussi par la loi de refroidissement de
Newton :
w
ww
TDEX- 11.9
om
kL
q′′y = hx (Tsat − TW ) = (Tsat − TW )
x δ
1
hx =
kL
=
kL
=
(k )
L
4 4
1 1
δ  4 x k L ( T sat - T W ) µ L  4
 4 x k L ( T sat - T W ) µ L  4
   
 ∆ H LV ρ L ( ρ L - ρ V ) g   ∆ H LV ρ L ( ρ L - ρ V ) g 
o.c
1 1
 ∆ H LV ρ L ( ρ L - ρ V ) g k L 4  4
 ∆ H LV ρ L ( ρ L - ρ V ) g k L3  4
hx =   = 
 4 x k L ( T sat - T W ) µ L   4 x( T sat - T W ) µ L 
-pr
Pour la valeur moyenne du coefficient h, on intègre en 0 et L :
1
1  ∆ H LV ρ L ( ρ L - ρ V ) g k L 3  4
L L
1
hL = ∫ hx dx = ∫   dx
L0 L 0 4 x( T sat - T W ) µ L 
1
1  ∆ H LV ρ L ( ρ L - ρ V ) g k L 3  4 L 1
ari
−
hL = 
L 4( T sat - T W ) µ L


∫ x 4 dx
0
L
1
 
1  ∆ H LV ρ L ( ρ L - ρ V ) g k L 3  4  1
 1 
 − +1 
hL =    x 4  
k
L 4( T sat - T W ) µ L    − 1 + 1 
  4  
0
ab
1
1  ∆ H LV ρ L ( ρ L - ρ V ) g k L 3  4
 4 34 
hL =    L 
L 4( T sat - T W ) µ L  3 
1
 4 1   ∆ H LV ρ L ( ρ L - ρ V ) g k L 3  4
1  34 
hL =  1   L 
al3
3 4 4   ( T sat - T W ) µ L  L 
1
 ∆ H LV ρ L ( ρ L - ρ V ) g k L 3  4
hL = 0.943  
 L ( T sat - T W ) µ L 
w.
ww
ww
w.
al3
ab
kari
-pro.
co
m
TDEX- 12.1
m
Travaux dirigés et exercices du chapitre 12: Rayonnement
Ex-12.1 (S) Soleil
co
2
Le flux solaire qui atteint l’atmosphère terrestre a une valeur moyenne de q’’S=1353 W/m . Les diamètres
9 7
du soleil et de la terre sont respectivement de DS=1.39 10 m et de DT =1.27 10 m et la distance des
11
centres de la terre et du soleil est LTS=1.5 10 m.
a) quelle est l’émittance totale du soleil ?
ro.
b) en supposant la surface du soleil comme un corps noir, quelle est sa température ?
c) à quelle longueur d’onde l’émittance spectrale du soleil est-elle maximale ?
d) en supposant que la surface de la terre est un corps noir dont la seule source d’énergie est celle
du soleil, estimer la température moyenne de la surface de la terre.
e) En fait, approximativement 30% du rayonnement solaire est réfléchi par la terre (albedo);
recalculer dans ces conditions la température moyenne de la surface de la terre et commenter.
Solutionnaire
-p
a) L’atmosphère terrestre se trouve à une distance du soleil
ari
LTS = 1.5 1011
Considérons une calotte sphérique de rayon LTS centrée sur le soleil. La surface de cette calotte
sphérique est S = 4 π l 2 l’énergie totale QSOLEIL émise par le soleil passe au travers de cette surface
et elle vaut donc:
2
QSOLEIL = 4 π L2TS qS′′ = 4 π (1.5 1011 ) (1353) = 3.8255 1026 W
k
Toute cette énergie est émise par la surface du soleil dont l’émittance totale Eb est donnée par :
ab
QSoleil 3.8255 10 26
Eb , soleil = 2
= = 6.3025 107 (W / m 2 )
D  6.06991018
4π  S 
 2 
al3
b) En utilisant la loi de Stefan-Boltzman :
1/ 4 1/ 4
E   6.3025107 
EB=σT
4
⇒ TSOLEIL =  B , SOLEIL  = −8 
= 5774 K
 σ   5.67 10 
w.
c) Avec la loi de Wien λ maxT = 2898 ( µ m.K) , on trouve donc :
2898
λ max = = 0.5 µ m
5774
ww
TDEX- 12.2
m
2
d) Le flux solaire q’’S=1353 W/m est celui qui traverse la calotte sphérique de rayon LTS définit en a.
La projection de la terre sur cette calotte sphérique est égale à la section de la terre
π DT 2
co
Sp =
4
alors que la surface de la terre vaut :
2
D 
Sterre = 4π  T  = π DT 2 .
 2 
ro.
Le bilan d’énergie s’écrit donc :
qT′′ Sterre = qS′′ S p
π DT 2
qT′′ π DT 2 = qS′′
⇒ qT′′ =
qS′′
4
4
= 338.25 W / m 2 = σ Tterre
 338.25 
1/ 4
4
ri-p
⇒ Tterre = −8 
= 277.9 K (≈ 5O C )
 5.67 10 
ka
e) uniquement 70% de l’énergie solaire est disponible
1/ 4
 0.7 338.25 
Tterre = −8 
= 254 K (≈ −19O C )
 5.67 10 
ab
o
Commentaire : La température moyenne de la terre est en fait de 15 C à cause des gaz à effets de serre
(H2O, CO2, N2O,CH4, O3) qui absorbent en partie le rayonnement émis par la terre.
Ex-12.2 (S) Efficacité d’une ampoule

al3
Calculer l’efficacité d’une ampoule de 100 W en supposant que son filament assimilé à une fine plaque de
5 mm de long par 2 mm de large est à une température de 2900 K et qu’il émet comme un corps noir. On
définit l’efficacité comme le rapport la puissance émise dans le visible sur la puissance totale de
l’ampoule.
Solutionnaire :
w.
a) On calcule la fraction de l’énergie émise dans le visible : F 0.4- 0.7 = F 0.7 - F 0.4 .
pour λ = 0.4 µ m : λT = 0.4 2900 = 1160
pour λ = 0.7 µ m : λT = 0.7 2900 = 2030
Par interpolation linéaire des valeurs du tableau p. 12.4, on trouve :
ww
TDEX- 12.3
om
F 0.4 − 0.00101 = 1160-1111 ⇒
F 0.4 = 0.00168
0.00252 − 0.00101 1222-1111
F 0.7 − 0.06672 2030-2000
= ⇒ F 0.7 = 0.0716
0.08496 − 0.06672 2111-2000
F 0.4-0.7 = 0.0716 − 0.00168 = 0.07
o.c
La surface du filament qui émet est égale à :
A S = 2. largeur .longueur = 2 5.10−3 2.10−3 = 2 10−5 (m 2 )
On suppose la surface comme un corps noir et la puissance totale pour l’émission dans le visible
vaut :
-pr
4
qvisible = A S σ TS4 F 0.4-0.7 = 2 10 −5 5.68 10−8 ( 2900 ) 0.07 = 5.62 watts
5.62
⇒ η= = 5.6%
100
Ex-12.3 (S)
ri
a) Déterminez les longueurs d’onde correspondant au maximum d’émission spectrale pour les surfaces
suivantes : soleil à 5800 K ; filament de tungsten à 2500 K ; surface métallique à 1500 K ; la peau à 305
surface métallique à 60 K refroidie par de l’azote liquide.
ka
b) Calculer les fractions du rayonnement solaire émises dans les trois domaines : ultraviolet, , visible et
infrarouge
Solutionnaire :
ab
a) Avec la loi de Wien λ maxT = 2898 ( µ m.K) , on obtient donc :
surface soleil tungsten metal peau métal

Température K 5800 2500 1500 305 60
λ max ( µ m ) 0.50 1.16 1.93 9.50 48.30
al3
b) Pour T=5800K, on s’intéresse aux fractions émises respectivement dans les gammes suivantes :
(voir figure p. 12.2)
UV Visible IR
a-b 0.01-0.4 0.4-0.7 0.7-100
5
λaT - λbT 58 - 2320 2320 - 4060 4060 - 5.8 10
w.
F 0 −λ a F 0 −λb - 0 - 0.125 0.125 -0.491 0.491 - 1

(tableau p 12.4)
F λb −λ a 0.125 0.366 0.509
ww
On retient donc que 12.5% du rayonnement solaire est dans l’UV, 36.6% dans le visible et 51% dans l’IR.
Commentaire : Les UV sont absorbés par l'ozone de la stratosphère, les IR sont absorbés par la vapeur
d'eau et le gaz carbonique de la troposphère puis réémis vers le sol et vers l'espace. Le Visible et le
proche Infra-Rouge (0.78-1.4 µm) atteignent la surface.
TDEX- 12.4
Ex-12.4 (S) Thermocouple
m
0
Un gaz s’écoule dans une conduite dont les parois sont à la température de 250 C. Un thermocouple
0
placé dans l’écoulement au centre de la conduite indique une valeur de 150 C . Sachant que le
2 0
coefficient de transfert de chaleur à la surface du thermocouple avec le gaz vaut h=100 W/m . C et que
co
l’émissivité de sa surface est ,=0.45, calculer la température du gaz.
(on supposera le thermocouple parfaitement isolé et le régime permanent).
Solutionnaire
ro.
Le bilan donne sur le thermocouple s’écrit : Ein-Eout=0
Il y a apport d’énergie par radiation et perte par convection :

Ein = εσ (Tpar
4
− Ttc4 ) Eout = h (Ttc − Tgaz )
Ein − Eout = 0
⇒ εσ (Tpar
4
− Ttc4 ) = h (Ttc − Tgaz )
Tgaz = Ttc −
εσ (Tpar
4
− Ttc4 )
ri-p
h
0
Il faut convertir les C en K :
ka
Tpar = 523 K Ttc = 423 K
εσ (Tpar
4
− Ttc4 ) 0.45 5.68 10−8 ( 5234 − 4234 )
Tgaz = Ttc − = 423 − = 412.08 K
h 100
Tgaz = 139.09 O C
ab
Commentaire : Le rayonnement perturbe significativement la mesure de la température du gaz.

L’installation d’un écran (pour éviter les échanges avec la paroi) serait appropriée.
al3
Ex-12.5 (S) Capteur solaire thermique
La surface AS d’un panneau solaire, non encore protégée par une vitre, a une émissivité de 0.1 et une
2
absorptivité du rayonnement solaire de α =0.95. Avec un ensoleillement Gsoleil =750 W/m , la
o o
température TS à la surface du panneau atteint 120 C. La température de l’air ambiant est de 30 C et on
o
peut supposer que la température du ciel dégagé est de -10 C. La surface AS du panneau échange aussi
w.
par convection avec l’air environnant avec un coefficient h alors que la face arrière échange sa chaleur QC
avec le fluide caloporteur.
2
Calculer l’efficacité du panneau pour les deux valeurs du coefficient h=1 et h=5 W/m C
ww
TDEX- 12.5
m
Solutionnaire :
En régime permanent : : Ein-Eout=0
co
- Le flux entrant est celui apporté par le soleil: Ein = AS α Gsoleil
- Il y a trois contributions sortantes :
- les pertes par convection avec l’air Econvection = hAS ( TS − Tair )
- les pertes par rayonnement avec le ciel froid (
Erayonnement = AS ε σ TS4 − Tciel
4
)
ro.
- le transfert au fluide caloporteur QC
(
ASα Gsoleil = AS ε σ TS4 − Tciel
4
)
+ hAS ( TS − Tair ) + QC
QC
α Gsoleil = ε σ (T 4 − Tciel
4
) + h (TS − Tair ) +
 QC

 AS

S
( )
ri-p
AS
 = qC = α Gsoleil − ε σ TS − Tciel − h (TS − Tair )


2
Par m , on récupère donc :
4 4
2
Si h=1 W/m C qC = 514 W / m 2
ka
2
Si h=5 W/m C qC = 154 W / m 2
On définit l’efficacité comme le rapport de l’énergie récupéré sur l’énergie solaire disponible :
514 154
ηh =1 = = 69% η h =5 = = 20%
750 750
ab
Commentaire : Afin d’éviter les pertes, on installe une vitre qui évite le refroidissement par convection et
retient le rayonnement IR émis par la surface du panneau.
al3
w.
ww
TDEX- 12.6
Ex-12.6 (S) Facteur de forme A
m
a) Pour la géométrie ci-contre, calculer le facteur
de forme F1-3 sachant que les rayons des disques
o
1 et 2 sont de 10 cm et 120 cm respectivement et
que L=2m
o.c
b) Si les températures des surfaces A1 et A3
O O
(supposées corps noir) sont de 600 C et 100 C,
quel est l’échange net de chaleur de A1 vers A3 ?
c) Refaire les calculs L=1 m
-pr
Solutionnaire :
a) Pour cette enceinte fermée, tous les rayons
quittant la surface A1 tombent sur les surfaces de
l’enceinte :
4
ri
∑F
k =1
1− k =1
⇒ F1−1 + F1− 2 + F1−3 + F1− 4

ka
La surface A1 est plane donc F11 =0. Les surfaces A1 et A4 sont dans le même plan donc F1-4=0.
⇒ F1− 2 + F1−3 = 1 ⇒ F1−3 = 1 − F1− 2

ab
Tous les rayons quittant A1 vont tomber sur A2 et A3.
On utilise le graphique de la page 12.8 pour déterminer F1-2.
L/r1=2/0.1=20 r2/L=1.2/2=0.6 F1-2≈0.26 Y F1-3=0.74

al3
q1−3 = A1 F1−3 σ (T14 − T34 ) = π 0.12 0.74 5.6810 −8 ( 273 + 600 ) − ( 273 + 100 )  = 736 W
4 4
b)  
q1−3 = 400 W
c) Pour L=1, F1-2≈0.6 Y F1-3=0.4
w.
ww
TDEX- 12.7
Ex-12.7 (S) Facteur de forme B
m
Un disque de diamètre D1 est placé à une distance L1 d’un four tubulaire de diamètre D2 et de longueur L2.
La paroi interne de ce four est à une température de 1000 K. On souhaite calculer la puissance radiative
émise par la paroi interne du four qui sera récupérée par le disque.
co
ro.
ri-p
ka
Solutionnaire :
Il faut estimer le facteur de forme F2-1. On a la relation de réciprocité

ab
A1
A1 F1− 2 = A2 F2 −1 ⇒ F2−1 = F1− 2
A2
On définit les surfaces fictives A3 et A4 aux deux extrémités du four. Les rayons qui quittent A1 et qui
entrent dans le four passent par A4 puis, ils tombent sur A2 ou ils sortent par A3 : on en déduit la relation
al3
des facteurs de formes :
F1− 4 = F1− 2 + F1−3 ⇒ F1− 2 = F1− 4 − F1−3
Les facteurs F1− 4 et F1−3 peuvent s’obtenir avec la figure p 12.8

w.
F1− 4 : L/r1=0.1/0.025=4 r2/L=0.05/0.1=0.5 F1-4≈0.20

F1−3 : L/r1=0.2/0.025=8 r2/L=0.05/0.2=0.25 F1-3≈0.06 F1-2≈-0.14
A   π 0.052 
 ( 5.67 10 1000 )
−8 4
ww
q2→1 = F2 −1 A2 Eb ,2 =  1 F1− 2  A2 Eb ,2 = F1− 2 A1 Eb ,2 = 0.14 

 A2   4 
q2→1 = 15.6 W
TDEX- 12.8
m
Ex-12.8 (S) Perte thermique d’un réservoir de stockage d’azote liquide
Un réservoir d’azote liquide est constitué de deux sphères

concentriques de diamètres de 1m et 1.5 m. Dans
co
l’espace entre les deux sphères, on fait le vide afin
d’isoler la sphère intérieure qui contient l’azote liquide.
Les surfaces des sphères en aluminium sont polies et
l’émissivité, ε , vaut 0.03 Sachant que les températures
des surfaces des deux sphères sont respectivement de
77 K et 273 K, calculer la quantité nette de chaleur reçue
ro.
par la sphère intérieure. Calculer la quantité d’azote
vaporisé par heure. Commenter
Données : pour deux sphères concentriques de surfaces

A1 et A2 (avec rayons r1 < r2), les facteurs de formes sont :
F1− 2 = 1 F1−1 = 0 F2 −1 =
A1
A2
-pF2 − 2 = 1 −
A1
A2
ari
3
Pour l’azote : )HLV = 198 kJ/kg D=808 kg/m
Solutionnaire : L’espace entre les deux sphères est fermée et la relation (28) p. 12.15 s’applique :
E b,i - E b, j 1 − εi 1 1− ε j
qi = avec R i + R i- j + R j = + +
( R i + R i- j + R j ) Ai ε i Ai Fi − j Aj ε j
k
On choisit i=2 et j=1

1− ε2 1 1 − ε1 A1
Rtotale = R 2 + R 2-1 + R1 = + + mais F2 −1 =
ab
A2 ε 2 A2 F2 −1 A1 ε1 A2
1− ε2 1 1 − ε1 1 − ε 2 1 1 − ε1 1 − ε 2 ε1 + 1 − ε 1 1− ε2 1
Rtotale = + + = + + = + = +
A2 ε 2 A A1 A1 ε1 A2 ε 2 A1 A1 ε1 A2 ε 2 A1 ε1 A2 ε 2 A1 ε1
2
A2
al3
E b,2 - E b,1 σ (T24 - T14 )

q i → j = q2 = = = 20.6 W = ∆H LV mazotevaporisé
1− ε2 1 1− ε2 1
+ +
A2 ε 2 A1 ε1 A2 ε 2 A1 ε1
20.6
⇒m= = 1.04 10 −4 kg / s = 0.104 g / s = 347 g / heure ≈ 9kg / jour
198 103
w.
4
volume du réservoir = π R13 = 0.523 m3 = 523 litres ⇒ capacité de 423 kg d ' azote
3
ww
TDEX- 12.9
Ex-12.9 (S) Cuisson d’une pièce peinte
m
La surface peinte d’une pièce métallique
doit subir un traitement thermique à 450K.
co
On utilise le four ci-contre composé de 3
parois de surface identique (largeur W,
hauteur H):
- la surface 1 est la résistance électrique à
une température de 1100 K
ro.
-la surface 2 est la pièce métallique à
traiter à une température de 450 K
- la surface trois est une paroi isolée.
Calculer la puissance nécessaire du four

(par m de hauteur) ainsi que la
température de la surface isolée. ri-p
Solutionnaire :
ka
Pour cette enceinte fermée, l’équation (31) p. 12.16 s’applique
N
Ji- J j E -J
qi = ∑ -1
= i,b i
j =1 ( Ai F i- j )
1- ε i
Aiε i
et on va écrire le système de trois équations à 3 inconnues :les radiosités des 3 surfaces J1, J2 et J3. On a
b
donc :
3a
3
J1- J j E -J J1- J 2 J1- J 3
pour la surface 1 : q1 = ∑ -1
= 1,b 1 = -1
+ -1
[1]
j =1 ( A1 F 1- j )
1-ε1 ( A1F 1- 2 ) ( A1F 1-3 )
A1ε 1
3
J2- J j E -J J 2 - J1 J2- J3
pour la surface 2 : q2 = ∑ -1
= 2 ,b 2 = -1
+ -1
[2]
1- ε 2
.al
j =1 ( A 2 F 2 - j ) ( A2 F 2-1 ) ( A2 F 2-3 )
A2ε 2
3
J3- J j J 3- J1 J3- J2
pour la surface 3 qui est isloée: q 3 = ∑ -1
=0= -1
+ -1
= 0 [3]
j =1 ( A3 F 3- j ) ( A2 F 3-1 ) ( A2 F 3- 2 )
w
Considérant la géométrie et la symétrie du four, le rayonnement émis par une surface est reçu pour moitié
par chacune des deux autres surfaces:
F 1-2 + F 1-3 = 1 F 2-1 + F 2-3 = 1 F 3-1 + F 3- 2 = 1
ww
F 1-2 = F 1-3 = F 2-1 = F 2-3 = F 3-1 = F 3- 2 = 0.5

1- ε 1 1- ε 2
Si on définit β1 = et β 2 = le système d’équations linéaires (1-3) à résoudre s’écrit
ε1 ε2
TDEX- 12.10
donc :
m
 1  E 1
J1 1 +  = 1,b + ( J 2 + J 3 ) avec E1,b = σ T14
 β1  β1 2
 1  E
co
1
J 2 1 +  = 2,b + ( J1 + J 3 ) avec E2,b = σ T24
 β2  β2 2
1
J 3 = ( J1 + J 2 )
2
ro.
Ce système peut être résolu numériquement par itération (méthode de Gauss-Siedel) avec Excel et son
Solveur (méthode déjà utilisée pour le projet no 1) :
E1,b1
(J + J )
+
β1 2 2 3
J1 =
 1 
-p 1 + 
 β1 
E2,b 1
+ (J + J )
β2 2 1 3
ari
J 2 ==
 1 
1 + 
 β2 
1
J 3 = ( J1 + J 2 )
k
2
E1,b + E2,b
On pourra initialiser les valeurs J1, J2 J3 respectivement à E1,b , E2,b ,
ab
2
On obtient :
J1 = 76 472 J 2 = 41 579 J 3 =59 025
E 1,b - J 1
al3
q1 = = 26.17 kW ( A1 = 1 m 2 )
1- ε 1
A1ε 1
1
E  4
J 3 = E 3,b = σ T 4
⇒ T3 =  3,b  = 1010 K
w.
3
 σ 
ww

TD 2 T Thérmique

Transféré par

Droits d'auteur :

Formats disponibles

TD 2 T Thérmique

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

TD 2 T Thérmique

Transféré par

Droits d'auteur :

Formats disponibles

m

Travaux dirigés et exercices

EX-2.2(S) Dérivation de l’équation thermique en coordonnées spériques

En utilisant le volume de contrôle ci-dessus ∆r r ∆ θ r sin θ ∆φ et en supposant une

diffusion thermique en sphérique

3) Bilan sur le volume de contrôle

d) Terme de génération (dans le volume):

e) Terme d’accumulation (dans le volume) :

On divise par ∆r ∆z ∆φ et on fait tendre les incréments vers 0

∆z r ∆φ qr" − ∆z r ∆φ qr" r qr" − r qr"

L’équation s’écrit donc :

La densité de flux s’exprime par la loi de Fourrier en utilisant l’expression du

(commentaire: dans la dérivée -p

3) Bilan sur le volume de contrôle

On divise par ∆r ∆θ ∆φ et on fait tendre les incréments vers 0

r ∆ θ r sin θ ∆φ qr" − r ∆ θ r sin θ ∆φ qr" r 2 qr" − r 2 qr"

L’équation s’écrit donc :

La densité de flux s’exprime par la loi de Fourrier en utilisant l’expression du

On divise par r 2 sin θ :

Ex-3.5 Pertes thermiques d'un tube isolé (à l'examen de sept. 1990)

Ex-3.6 Profil dans une résistance électrique (à l=examen de

Une résistance électrique (épaisseur b, largeur W, hauteur L)

b) Sans utiliser l'expression du profil obtenue en a),

observe-t-on la température maximale ? Tracer sur la figure ci-dessus la forme du profil de

Ex-3.7(S) Température de l'anse d'un chaudron

L'anse d'un chaudron est une tige métallique de 10 mm de rayon

a) en cuivre, k=385 W/(m.K)

(on pourra négliger la courbure de l'anse et le gradient radial de

c) Avec Excel par exemple, étudier l'effet de h et de k sur la

Ex-3.8(S) Ailette circulaire (à l'examen de 1989)

a) Faites un bilan de chaleur sur un élement de volume

approprié de l'ailette et obtenez l'équation différentielle

b) Quelles sont les différentes conditions frontières

c) Si on appelle qC, les pertes de chaleur (en watts) au

Une ailette de largeur W, d'épaisseur t et de longueur L

b) Posez les conditions frontières.

constantes d'intégration figurant dans cette expression

GCH-18243 Examen partiel No 1 - 25 pts Aucun document sauf 1 feuille (recto-verso)

Un mur de béton (épaisseur b, largeur W, hauteur

a) Faites un bilan de chaleur sur un volume de

II) Stockage d’air liquide dans un réservoir sphérique (6 pts)

Un réservoir sphérique de rayon r1=1.5 m, contient

ambiant à la température T∞=283 K et le coefficient de

a) En utilisant, le concept des résistances thermiques,

On rappelle les expressions des résistances pour la sphère creuse:

Un fil électrique de 1 mm de diamètre est dénudé sur une longueur de 1 m. Calculer la

• Dans un bilan de chaleur associé à un problème de conduction, si il y a un terme de génération de

V) Ailette en forme de cylindre creux (8 pts)

Une ailette, de longueur L, a la forme

cylindre est refroidie par de l’air à la température

On supposera que le flux de chaleur axial à x=L

a) Faites un bilan de chaleur sur un élément de

b) Quelles sont les conditions frontières ? (1 pt)

c)Obtenez l’expression du profil de température. (2.5 pts)

3) Bilan sur le volume de contrôle

Ce qui rentre – Ce qui sort =0

On utilise la loi de Fourier

CF1 : = q0" =  −k "

CF2 : à x = L ⇒ q"x = hF (Tx = L − TF )

Et en remplaçant C1 et C2 dans (1), on obtient l’expression de la variation de la température dans le mur :