Poly 3 TD T Thérmique

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 44

TD de Transfert de Chaleur.

Série N°1 Généralités

I - 1°) Exprimer en Kelvin, degrés Fahrenheit et degrés Rankine (température absolue dans le système
anglo-saxon) les températures de 0°C, 50°C, 100°C, -17.78°C, -273.15°C.
2°) Déterminer la température à laquelle le nombre qui l’exprime est le même en °C et °F. Même question
en K et en °R.
II - 1°) On rencontre dans la littérature anglo-saxonne, la chaleur massique exprimée en Btu/lbF (Btu :
British thermal unit, lbf : pound force). Calculer sa valeur dans le S.I. ainsi qu’en C.G.S. On donne
l1b=453.5g, lBtu=l055 J.
2°) En utilisant les facteurs de conversion entre W et Btu/h, m et ft, K et R, exprimer la constante de Stefan-
Boltzmann σ = 5.67 × 10−8 W / m2 . K4 et le coefficient de h (W/m2.°C) en unité anglo-saxonne Btu / h. ft 2 . R 4 .
Les facteurs de conversion de W, m, et K sont donnés dans le tableau suivant :
1 W = 3.41214 Btu / h
1 m = 3.2808 ft Facteur de conversion
1 K = 1.8 R
3) En utilisant les facteurs de conversion entre °C et °F, le coefficient de conversion entre W et de Btu/h, m
et ft, exprimer le coefficient de h (W/m2.°C) en unité anglo-saxonne (Btu/h.ft2.°F)

III - Une résistance électrique de forme cylindrique (D=0,4cm, L=1,5cm) sur un circuit imprimé dissipe une
puissance de 0,6 W. En supposant que la chaleur est transférée de manière uniforme à travers toutes
les surfaces. Déterminer :
(a) la quantité de chaleur dissipée par cette résistance au cours d'une période de 24 heures,
(b) le flux de chaleur,
(c) la fraction de la chaleur dissipée par les surfaces du haut et du bas.

IV - Un réservoir contient 3m3 d’eau chaude à Ti=80°C. Il est parfaitement calorifugé sauf sur une partie
dont la surface est S=0.3m2. On constate qu’où bout de ∆t=5 heures, la température de l’eau a baissé de
0.6°C quand la température ambiante est de 20°C. En supposant que la capacité calorifique du réservoir est
de 103 kcal/°C.
1) Calculer:
1°/ la quantité de chaleur perdue en 5 heures,
2°/ le flux de chaleur à travers le couvercle,
3°/ la densité de flux thermique à travers le couvercle,
4°/ la résistance thermique du couvercle,
5°/ le coefficient global de transmission thermique.
On donnera les résultats dans les systèmes M.K.H et S.I.
2) Que se passerait-il au bout de lj, 10j ?

V) - On remplit en 5 minutes une baignoire de 500 litres avec de l’eau chaude prélevée dans un réservoir
supposé à température constante de 50°C. La canalisation de diamètre extérieur de 16 mm et intérieur de 14
mm a une longueur de 10 m.
1) Si la chute de température entre le réservoir et le robinet est de 2°C. Calculer le flux de chaleur
perdu par la canalisation pendant le remplissage de la baignoire ainsi que les densités de flux de chaleur
correspondants aux surfaces intérieur et extérieur de la canalisation.
A quel pourcentage de perte de chaleur cela correspond-il si on suppose que l’eau froide était à 15°C avant
d’être chauffée dans le réservoir?
2) La température de l’eau dans la canalisation revient à la température de 20°C au bout de 30
minutes. Calculer le flux de chaleur perdu dans ces conditions.
3) Quelles quantités d’eau chaude faudrait-il prélever en une seule fois pour que les pertes en énergie
entre le réservoir et le robinet ne représentent que 10%, 20%, 50%?
Unités thermiques
S.I. :
Température :
T(K)= T(°C)+273,15 ∆T(K) = ∆T(°C)
Quantité de chaleur :
1 Joule = 0,239 cal = 0,948.10-3 Btu
Flux thermique :
1 W = 0,239 cal.s-1 = 0,86 kcal.h-1 =3,41 Btu.h-1
Densité de flux :
1 Wm-2 =23,6.10-6 cal.cm-2.s-1 = 0,86 kcal .m-2.h-1 =0,317 Btu.ft-2.h-l
Coefficient global de transmission thermique :
1 W.m-2K-l = 23,9.10-6 cal.cm-2.s-1.°C-1 = 0,86 kcal.m-2.h-1.°C-1= 0,176 Btu.ft-2.h-1.°F-1
Résistance thermique :
1 K.W-1 = 4,1855°C.s.ca1-1= 1,163°C.h.kca1-1 = 0,526°F.h.Btu-1

C.G.S. : Centimètre, gramme et seconde.


Quantité de chaleur :
1 cal = 4,1855 J
Flux thermique :
1 cal.s-1= 4,1855 W
Densité de flux :
1 cal.cm-2.s-1 = 4,1855.104 W.m-2
Coefficient global de transmission thermique :
1 cal.cm-2.s-l.°C-1= 4,1855.10-4 W.m-2.K-1
Résistance thermique :
1°C.s.cal-1 = 0,239 K.W-1

ANGLO-SAXON :
Température :
T(°R) = 1,8 T(K) T(°F) = 1,8 T(°C)+32 ∆T(°F) = 1,8 ∆T(°C)
Quantité de chaleur :
1 Btu = 1055 J (Btu = British Thermal Unit)
Flux de chaleur :
1 Btu.h-1 = 0,293 W
Densité de flux :
1 Btu.ft-2.h-1 = 3,15 W.m-2 (1 ft = 0,3048 m) (Feet : pieds en français)
Coefficient global de transmission thermique :
1 Btu.ft-2.h-1.°F-1 = 5,68 W.m-2K-1
Résistance thermique
1°F.h.Btu-1= 1,9 K.W-1
Solution
I) 1°) T (K ) = T (°C ) + 273.15
T (° F ) = 1.8 * T (°C ) + 32
T (° R ) = 1.8 * T (°C ) + 491.67

T(°C) T(K) T(°F) T(°R)


100 373.15 212 631.67
50 323.15 122 581.67
0 273.15 32 491.67
-17.78 255.37 0 459.67
-273.15 0 -459.67 0
2°) T (° F ) = 1.8 * T (°C ) + 32 avec T(°F)=T(°C)=x
32
x=−
= −40
0.8
0 °R=0 K
1Btu 1055 J J 4187 cal cal
II) 1°) = = 4187 = =1
1lb * 1F 1°C kg.K 4185 g .°C g.°C
0.4535kg *
1.8

Btu
3.41214
W h Btu
2°) σ = 5.67 * 10 -8 2 4 = 5.67 × 10 −8 2 4
= 0.171 * 10 -8
m .K (3.2808 ft) (1.8 R) h.ft 2 .R 4
3°) Le facteur de conversion entre °C et °F est donné par : 1°C = 1.8°F
Btu
3.41214
W h Btu
1 2 = = 0.1761
m .°C (3.2808 ft) (1.8 °F)
2
h.ft 2 .° F
III)
Hypothèses : La chaleur est transférée de façon uniforme à travers toutes les surfaces.
Analysis (a) la chaleur dissipée par cette résistance au cours d'une période de 24 heures est :
Q = Q& ∆t = (0.6 W)(24 h) = 14.4 Wh = 51.84 kJ (1 W.h = 3600 W.s = 3.6 kJ)
total total
(b) Le flux de chaleur à la surface de la résistance est
πD 2 π (0.4 cm)2
Stotal = 2 + πDL = 2 + π (0.4 cm)(1.5 cm) = 0.251 + 1.885 = 2.136 cm 2
4 4
Q& total 0.60 W
q& = = = 0.2809 W/cm 2
S total 2.136 cm 2
(c) En supposant que le coefficient de transfert de chaleur est uniforme, le transfert de chaleur est
proportionnelle à la surface. Ensuite, la fraction de la chaleur dissipée par le haut et le bas des surfaces de la
résistance devient
Qhaut −bas S haut −bas 0.251
= = = 0.118 or (11.8%)
Qtotal S total 2.136
Discussion : Le transfert de chaleur par le haut et le bas des surfaces sont faibles par rapport à celle transféré
par la surface cylindrique. C’est le cas des murs au sens thermique du terme.
IV) 1/ Hypothèse : les températures de l’eau et du réservoir sont homogène est égale à une température
moyenne Teau(t) à chaque instant t.

1°) La quantité de chaleur perdue en 5 heures est :


Q (t = 5h) = ((m.c )eau + (m.c )réservoir )(Teau (t = 5h) − Teau (t = 0) )
(( ) ( ))
Q(t = 5h) = 3 *10 3 *1kcal / °C + 10 3 kcal / °C (− 0.6°C ) = −2.4 10 3 kcal (Système MKH)
Q(t = 5h) = −2.4 10 3 kcal = −2.4 10 3 * 4.185kJ = −10044kJ (Système
International)

2°) calcul du flux de chaleur à travers le couvercle


Q(t = 5h) − 2.4 10 3 kcal kcal
φ (t = 5h) = = = −480 (Système MKH)
∆t 5h h
4185 J
φ (t = 5h) = −480 = −480 *1.1625 W = −558 W (SI)
60 * 60s

3°) calcul du densité de flux thermique à travers le couvercle


kcal
− 480
φ (t = 5h) h = −1600 kcal
ϕ (t = 5h) = = 2
(Système MKH)
S 0.3m m 2 .h
kcal 4185 J W W
ϕ (t = 5h) = −1600 2 = −1600 2
= −1600 * 1.1625 2 = −1860 2 (SI)
m .h 3600 m .s m m

4°) résistance thermique du couvercle


Par analogie thermique électrique :
∆T
∆U = Rélectrique * I ∆T = Rthermique * φ ⇒ Rthermique =
φ
∆T Tair − Teau (t = 5h) (20 − (80 − 0.6 ))°C °C.h
Rthermique = = = = 0.12375 (Système MKH)
φ φ (t = 5h) − 480
kcal kcal
h
°C.h °C.3600.s 1 °C °C
Rthermique = 0.12375 = 0.12375 * = 0.12375 * = 0.1064 (SI)
kcal 4185.J 1.1625 W W

5°) le coefficient global de transmission thermique K

φ φ 1
φ = K * S * ∆T ⇒ K= = =
S * ∆T S * Rthérmique * φ Rthérmique * S
1 1 kcal
K= = = 26.93 2 (Système MKH)
Rthérmique * S 0.12375 * 0.3 m .h.°C
kcal W
K = 26.93 = 26.93 * 1.1625 = 31.32 2 (SI)
m .h.°C
2
m .K
2/ Si en suppose que Rthérmique = Cte
Teau (t ) − Tair dT d (Teau (t ) − Tair ) d (Teau (t ) − Tair ) dt
φ (t ) = = − mc = −mc → =−
Rthérmique dt dt Teau (t ) − Tair Rthérmique .mc
t

(T (t ) − T ) = (T (0) − T )e

Rthérmique .mc
eau air eau air avec Rthérmique .mc = 0.12375 * 410 3 = 495h −1
t (h)

Teau (t ) = 20 + 60 * e 495

t= 1j=24h Tair(24h)=77°C
t=10j=240h Tair(240h)=56.9°C
V
1) Soit Qp la chaleur perdue pendant le remplissage la baignoire ( ∆t = 5 min )
QP = 500(kg ) * 1(kcal / kg.°C ) * 2(°C ) = 1000kcal = 4185 kJ
Soit Q la chaleur nécessaire pour chauffer 500 litres d’eau :
Q = 500 *1 * (50 − 15) = 17.5.10 3 kcal = 73.2.10 3 kJ
Le flux perdu est :
Q 4185
φP = P = = 13.95 kW = 12.10 3 kcal.h −1
∆t 5 * 60
Densité de flux rapporté à la surface intérieure ou extérieure de la canalisation :
φP φp φP
ϕ = = = L=10m
S 2.π .R.L π .D.L
Densité de flux intérieure de la canalisation:
φ 13.95 kW
ϕ int = P = = 31.7 kW .m − 2
S int π .14.10 .10 −3

Densité de flux extérieure de la canalisation :


φ 13.95 kW
ϕ ext = P = = 27.8 kW .m − 2
S ext π .16.10 −3.10
Q m.C.∆T perter ∆T perte 2
Pourcentage de perte : P = = = = 5.7%
Q m.C.∆Tchauffage ∆Tchauffage (50 − 15)
2) Chaleur perdue pendant le refroidissement de l’eau dans la canalisation Qr :
Hypothèse : à chaque instant t la température de l’eau dans la canalisation est homogène.
(meau du canalisation .C )eau dans la canalisation .∆θ Di2
φr = avec : meau du canalisation = ρ .π . .L et
∆t 4
kJ
C = 4185
kg

meau du canalisation = 10 .3,14.


3 (14.10 )
−3 2
.10 = 1,54 kg
4
1,54 * 4185 * (50 − 20 )
φr = = 107 W
30 * 60
Les densités de flux :
Densité de flux intérieure de la canalisation:
φ 107 kW
ϕ int, r = r = = 243 W .m − 2
S int π .14.10 .10−3

Densité de flux extérieure de la canalisation :


φ 107 kW
ϕ ext , r = r = = 213 W .m − 2
S ext π .16.10 .10 −3

3) soit y les pertes en énergie :


En note :
Qperdue : quantité d’énergie perdue,
Qconsommée : quantité d’énergie consommée (chauffage par effet joule),
Qutilisée : quantité d’énergie utilisée dans la baignoire.

Qconsommée = Qutilisée − Q perdue Seul le terme perte est négatif.

Q perdue Qutilisée
y= = 1−
Qconsommée Qconsommée
Qutilisée = x * c * (Tbaignoire − Teau froide ) = x * c * (48 − 15) = 33 * x * c
Qconsommée = (x + meau du canalisation )* c * (Tréservoire − Teau froide ) = (x + 1.54 ) * c * (50 − 15)
Qutilisée 33 * x
y = 1− = 1−
Qconsommée 35 * ( x + 1.54 )
53,865(1 − y )
x=
35 * y − 2
Résultats :
y 0.1 0.2 0.5
x(l) 32.3 8.6 1.7
www.almohandiss.com

TD : Transferts de Chaleur
Série N°: 2
Conduction de la chaleur en régime stationnaire
Problèmes 1 D.
I - Un mur de béton de 15 cm d'épaisseur sépare une pièce à la tem pérature Ti = 20 °C de
l'extérieur où la température est Te = 5 °C.
-2
On donne : hi = 9.1 W.m .K-1
-2
he = 16.7 W.m .K-1
-1 -1
λ = 1,74 W.m .K .
Calculer :
- la résistance thermique totale
- la densité de flux
- les températures interne et externe du mur.

II - Le mur d'un local est constitué de trois matériaux différents :

- du béton d'épaisseur e1 = 15 cm à l'extérieur (conductivité thermique λ1 = 0,23 W.m-1.K-1),


- un espace e 2 = 5cm entre les deux cloisons rem pli de polystyrène expansé (conductivité
thermique λ2 = 0,035 W.m-1.K-1),
- des briques d'épaisseur e3 = 5 cm à l'intérieur (conductivité thermique λ3 = 0,47 W.m-1.K-1).
1) On a mesuré en hiver, les températures des parois intérieures θi et extérieure θe qui étaient
θi = 25 °C et θe = - 8 °C.
a) Donner la relation littérale, puis calcule r la résistance therm ique du mur pour un
mètre carré.
b) Donner la relation littérale , puis calculer le flux thermique dans le m ur pour un
mètre carré.
c) Calculer la quantité d e chaleur transmise par jour à travers un m ètre carré de m ur,
pour ces températures. En déduire la quantité de chaleur transm ise, par jour, à travers 10 m 2
de mur.
d) Tracer la courbe de variation de tem pérature θ = f(e) à travers le mur, de paroi
intérieure à paroi extérieure.
2) Les résistances thermiques superficielles interne et externe du mur ont respectivement pour
valeur : 1/hi = 0,11 m2.K.W-1 et 1/he = 0,06 m2.K.W-1
a) A quels types de transfert thermique ces données se rapportent-elles ?
b) Calculer les températures ambiantes extérieure θae et intérieure θai.

III- La paroi d’un four électri que industriel est constituée de plusieurs m atériaux comme
l’indique le schéma ci-dessous.
www.almohandiss.com

Données numériques.
Température ambiante intérieure : θi = 1092 °C.
Température ambiante extérieure : θe = 32 °C.
Surface intérieure du four : S = 8,00 m2.
Résistance superficielle interne pour un mètre carré de paroi : 1/hi = ri = 0,036 m2.K.W-1
Résistance superficielle externe pour un mètre carré de paroi : 1/he = re = 0,175 m2.K.W-1

Caractéristiques des divers matériaux :


Matériau Épaisseur Conductivité thermique
Brique à feu e1 = 230 mm λ1 = 1,04 W.m-1.K-1
Brique réfractaire e2 = 150 mm λ2 = 0,70 W.m-1.K-1
Laine de verre e3 = 50 mm λ3 = 0,07 W.m-1.K-1
Acier e4 = 3 mm λ4 = 45 W.m-1.K-1

1 . Exprimer littéralement puis calculer la résistance thermique globale R de un mètre carré


de paroi.
2 . Exprimer littéralement puis calculer la densité de flux thermique ϕ (puissance thermique
par unité de surface) traversant la paroi.
3 . Déterminer 1es températures au niveau des diverses interfaces de 1'intérieur vers
l'extérieur θsi, θ1, θ2, θ3, θse.
4 . En admettant que la transmission de la chaleur est uniforme sur l'ensemble des parois du
four, calculer la puissance électrique p nécessaire à son fonctionnement à vide.
5 . Calculer le coût de fonctionnement journalier du four sachant que le prix du kW.h est 1,50
Dh.

IV - Soit la sec tion droite d’une résistance électrique constituée d’un cœur en graphite,
entouré d’une enveloppe de verre, elle-m ême enrobée de m icanite (mélange d e mica et d e
résine phénolique agissant comme isolant électrique et therm ique). On de mande de
déterminer l’épaisseur optimum de micanite en vue d’assurer le refroidissement maximum de
l’enveloppe de verre sachant que 40% de l’éner gie électrique dissipée da ns la résistance est
perdue par convection - rayonnem ent avec l’ambiance à 20°c avec un coefficient d’échange
h=17W/m²K. La conductibilité thermique de la micanite est λ=0,1W/mK. Les caractéristiques
électriques de la résistance sont : Pe=1W, Re=106Ω. Sa longueur est de 30mm et son diamètre
est de 1mm.
1 . Quelle est la température à l’interface graphite - verre?
2 . Comparer les résultats avec et sans isolant.
www.almohandiss.com

Solution de la série N°: 1


Transferts de chaleur par conduction

I : ϕ = hi(Ti – TiP) = (λ/e)(TiP – TeP) = he(TeP – Te) = (Ti – Te)/R


Avec R = 1/hi + e/λ = 1/he = 0,11 + 0,15/1,74 + 0,06 = 0,2562 m2.K.W-1
De la dernière égalité :
ϕ = 15/0,2562
ϕ = 58,546 W.m-2
De la première égalité :
TiP = Ti - ϕ/hi = Ti - ϕ.Ri = 20 – 58,546*0,11
TiP = 13,6 °C
De la même façon :
TeP = Te + ϕ/he = Te + ϕ.Re = 5 + 58,546*0,06
TiP = 8,5 °C
II :
1) a) R = Rbe + Rpoly + Rbr = e1/λ1 + e2/λ2 + e3/λ3
= 15.10-2/0,23 + 5.10-2/0,035 + 5.10-2/0,47
R = 2,187 K.W-1.m2
b) ϕ = (θi – θe)/R
ϕ = 15,088 W.m-2
c) Q = ϕ.t = 15,088*24*3600
Q = 1,304.106 J.m-2
Pour 10 m2 de mur :
Q’ = 10Q
Q’ = 13,04.106 J.m2
d) ϕ = (θ1 - θe).λ1/e1
θ1 = θe + e1.ϕ/λ1 = - 8 +9,84*
θ1 = 1,84 °C
De la même façon :
ϕ = (θi - θ2).λ3/e3
θ2 = θi – e3.ϕ/λ3 = 25 – 15,088.5.10-2/0,47
θ2 = 23,39 °C
2) a) Transferts thermiques intérieur et extérieur par convection.
b) ϕ = he(θe - θae)
θae = θe - ϕ/he = - 8 – 15,088*0,06
θae = - 8,9 °C
ϕ = hi(θai - θi)
θai = θ1 + ϕ/hi = 25 + 15,088*0,11
θae = 26,7 °C
III- 1) R = Ra + Rlv + Rbr + Rbf + 1/hi + 1/he
RT = 1,36 m2.K.W-1
2) ϕ = (θI - θe)/R
ϕ = 779 W.m-2
3) ϕ = hi(θi - θsi) θsi = θi - ϕ/hi = 1064 °C
ϕ = λ1(θsi - θ1)/e1 θ1 = θsi - ϕ.e1/λ1 = 892 °C
ϕ = λ2(θ1 - θ2)/e2 θ2 = θ1 - ϕ.e2/λ2 = 725 °C
ϕ = λ3(θ2 - θ3)/e3 θ3 = θ2 - ϕ.e3/λ3 = 169 °C
ϕ = he(θse - θe) θse = θe + ϕ/he = 168 °C
www.almohandiss.com

4) Puissance = Φ = ϕ.S = 6,23 kW


5) E = P.t = 149,6 kWh
coût = 255 Dh
IV

Micanite

r2=1mm
Te T3 T2 Verre

Graphite

Qconv + ray = 0.40 PE = 0.4W


( N .B. PE = I ² RE = I ² 106 ⇒ I = 10 −3 A faible courant )

Avec isolation

⎡ ln r3 r2 1 ⎤
T2 − Te = Qc −v ⎢ + ⎥
⎣ 2π kenv L 2 π rext hext L ⎦
= 0.4(93.984 + 53.073) + 58.8 ⇒ T2 ≅ 79°c

Sans isolant

1
T2 − Te = Qc −v = 0.4 x312 = 124.8 ⇒ T2 ≅ 145°c
2 π rext hext L
www.almohandiss.com

Un tube à vapeur de longueur L= 15 ft (pied), de rayon intérieur r1= 2in (pouces), de rayon extérieur
r2= 2.4in (pouces) et de conductivité thermique λ= 7,2 Btu/h.ft.°F. Le
fluide dans le tuyau est à une température moyenne de 250°F, et la
moyenne du coefficient de convection de chale ur sur la surface
intérieure est h=1.25 Btu/h.ft2.°F. Si la température moyenne à surface
extérieur est T2= 160°F, (a) exprim er l'équation de la c haleur et les
conditions aux limites dans le cas d ’un régime stationaire et pour une
conduction unidimensionnelle, (b) d éduire la v ariation de la
température dans la co nduite par la réso lution de l’équation
différentielle.

Hypothèses : 1- la conduction thermique est stable et


unidimensionnelle car le tuyau est long par rapport à son épaisseur, et il ya une symétrie thermique sur la ligne
centrale. 2- La conductibilité thermique est constante. 3- Il n'ya pas de génération de chaleur dans le tuyau.
Propriétés : La conductivité thermique est donnée à be λ = 7.2 Btu/h⋅ft⋅°F.
Analysis (a) Notant que le transfert de chaleur est une dimension radiale r de la direction, la formulation
mathématique de ce problème peut être exprimé comme

T =160°F
d ⎛ dT ⎞
⎜r ⎟=0
dr ⎝ dr ⎠ Vapeur
dT (r1 ) 250°F
et −k = h[T∞ − T (r1 )] h=1.25
dr
T (r2 ) = T2 = 160° F L = 15 ft
(b) L'intégration de l'équation différentielle une fois en fonction de r donne
dT
r = C1
dr
En multipliant les deux côtés de l'équation ci-dessus par dr/r , l’integration conduit à :
dT C1
=
dr r
T (r ) = C1 ln r + C2
où C1 et C2 sont des constantes arbitraires. Appliquant les conditions aux limites :
C1
r = r1: −k = h[T∞ − (C1 ln r1 + C2 )]
r1
r = r2: T (r2 ) = C1 ln r2 + C2 = T2

D’où les solutions de C1 et C2 sont :


T2 − T∞ T2 − T∞
C1 = and C 2 = T2 − C1 ln r2 = T2 − ln r2
r k r2 k
ln 2 + ln +
r1 hr1 r1 hr1
T 2 − T∞ r
T (r ) = C1 ln r + T2 − C1 ln r2 = C1 (ln r − ln r2 ) + T2 = ln + T2
r2 k r2
ln +
r1 hr1
(160 − 250)°F r r
= ln + 160°F = −24.74 ln + 160°F
2.4 7.2 Btu/h ⋅ ft ⋅ °F 2.4 in 2.4 in
ln +
2 (12.5 Btu/h ⋅ ft 2 ⋅ °F)(2 / 12 ft )
(c) Le flux de chaleur à travers le tuyau est
www.almohandiss.com

dT C T − T∞
φ = −k .S = − k (2πrL) 1 = −2πLk 2
dr r r k
ln 2 +
r1 hr1
(160 − 250)°F
= −2π (15 ft)(7.2 Btu/h ⋅ ft ⋅ °F) = 16,800 Btu/h
2.4 7.2 Btu/h ⋅ ft ⋅ °F
ln +
2 (12.5 Btu/h ⋅ ft 2 ⋅ °F)(2 / 12 ft )
Transfert de chaleur par les murs composés.
Un mur de 4 m de haut et 6 m de long est composé de deux plaques d'acier (λa= 15 W / m.° C) de 2 cm
d'épaisseur chaqune, séparés par 1 cm d'épaisseur et 20 cm de large ur des barres d'acier espacé de 99 cm .
L'espace entre les plaques d'acier est rempli d’isolant de fibre de verre (λi= 0,035 W/m°C). Si la différence de
température entre la surface intérieure et celle de l’extérieure du mur est 22°C :
1) déterminer le flux de chaleur échangé à travers le mur,
2) déterminer le flux de chaleur échangé a travers le mur si on ignore les barres d'acier entre les plaques,
car ils n'occupent que 1 pour cent de surface d’échage.
2 cm 20 cm 2 cm

99 cm

1 cm

Le mur est construit de deux grandes plaques d’aciers séparés par 1 cm d'épaisseur des barres d'acier espacé de
99 cm. L'espace restant entre les plaques d'acier est rempli d'isolant en fibre de verre. Le flux de chaleur à travers
la paroi du mur est à d éterminer, et il est à évaluer si les barres d'acier entre les plaques peuvent être ignorées
dans l'analyse, car ils n'occupent que 1 pour cent de la surface d’échange de chaleur.
Hypothèses :
1 Le transfert de chaleur est stationnaire.
2 Le transfert de chaleur àtravers le mur peut être estimé à une dimension.
3 les conductivités thermiques sont constantes.
4 Les surfaces de la paroi sont maintenues à des températures constantes.
Nous considérons que l'analyse de 1 m de haut et 1 m de large portion de l a paroi qui est
représentative de tout le mur.
Résistance thermique du réseau et les résistances sont

R1 R2 R4
T1 R3 T2
www.almohandiss.com

L 0.02 m
R1 = R4 = Racier = = = 0.00133 °C/W
λacier .S acier (15 W/m.°C)(1 m 2 ) 2 cm 20 cm 2 cm

L 0.2 m
R2 = Rbarre = = = 1.333 °C/W
λacier S barre (15 W/m.°C)(0.01 m 2 )
L 0.2 m
R3 = Risolant = = = 5.772 °C/W
λisolant S isolant (0.035 W/m.°C)(0.99 m 2 )
1 1 1 1 1 99 cm
= + = + ⎯→ Req = 1.083 °C/W

Reqv R2 R3 1.333 5.772
Rtotal = R1 + Reqv + R4 = 0.00133 + 1.083 + 0.00133 = 1.0856 °C/W 1 cm
Le flux de chaleur pour une surface de 1 m2 est :
∆T 22 °C
φ= = = 20.26 W
Rtotal 1.0857 °C/W
Le flux total echangé à traver le mu rest :
φ total = (4 × 6)Q& = 24(20.26 W) = 486.3 W
Si les barres d'acier entre les plaques sont ignorées dans l'analyse, le flux de chaleur pour une surface de 1 m2 est
∆T ∆T 22 °C
φ= = = = 3.81 W
Rtotal R1 + Rinsulation + R4 (0.00133 + 5.772 + 0.00133)°C/W
Le flux de chaleur qui traverse les barres d'acier entre les plaques est de (20.26-3.81)/20,26=81,2% du transfert
de chaleur qui travers les le mur, malgré les faibles espace qu'ils occupent, et bien sûr, leur effet ne peut être
négligé. Le raccordement des bars servent de "ponts thermiques".
www.almohandiss.com

TD de Transfert de Chaleur.
Série N°3.

Milieu thermiquement mince


A Déterminer les relations de la longueur caractéristique d'un mur d'épaisseur 2L, un
très long cylindre de rayon ro et une sphère de rayon ro .
B Pour réchauffer du lait pour bébé, la mère verse le lait dans une mince paroi de verre
dont le diamètre est de 6 cm. La hauteur du lait dans le verre est de 7 cm. Elle place ensuite le
verre dans une grande casserole remplie d'eau chaude à 60°C. Le lait est agité en permanence,
de sorte que sa température est uniforme en tout temps. Si le coefficient de transfert de
chaleur entre l'eau et le verre est 120W/m2.°C, déterminer le temps qu’il faut pour que le lait
se réchauffe de 3°C à 38°C. En admet que les propriétés du lait sont les mêmes que ceux de
l'eau. Dans ces conditions, le lait peut être traité comme un milieu thermiquement mince,
Pourquoi ?
Donnés : La conductivité thermique, densité et la chaleur spécifique de l’eau à 20°C sont λ =
0.607 W/m.°C, ρ = 998 kg/m3 et Cp = 4.182 kJ/kg.°C.

Ailette
Prenons une cuillère en acier inoxydable (λ= 15 W/m.°C),
partiellement immergé dans l'eau bouillante à 93°C dans une
cuisine à 24°C. Le manche de la cuillère est une section de
0,2cm*1,3cm, et s'étend dans l'air de 18cm de la surface libre de
l'eau. Si le coefficient de transfert de chaleur à la surface de la
cuillère exposée à l’air est h=17 W/m2.°C, déterminer la
différence de température à la surface du manche de la cuillère.
Indiquez vos hypothèses.
Facteur de forme
Deux tuyaux parallèles, de 5m de long, dans une épaisse couche de
béton, circule dans l’un de l'eau chaude et dans l’autre de l’eau
froide comme le montre la Figure. Les diamètres des tuyaux sont 5
cm et la distance entre leurs axes est de 30 cm. Les températures des
surfaces chaude et froide des tuyaux sont 70°C et 15°C
respectivement. La conductivité thermique du béton est λ=0,75
W/m.°C. Déterminer le flux de chaleur échangé entre les tuyaux.

Exemple de résolution de problème de conduction 2D par différences finies


Un long barrage en béton (λ=0,6 W/m.°C, coefficient d’absorption de
rayonnement : αs = 0,7 m2/s), de section triangulaire dont la surface
exposée est soumise à des flux de chaleur solaire q&s = 800 W/m2 et de
chaleur convectif avec l’air a une température T0=25°C et un coefficient de
transfert de chaleur h0=30 W/m2.°C. Les deux mètres du haut du barrage,
section verticale, est l'objet de convection de l'eau à Ti=15°C avec un
coefficient de transfert de chaleur hi=150 W/m2.°C. Le flux de chaleur
échangé à travers la surface de la base, de deux mètres de long, est
considérée comme négligeable. Utiliser la méthode de différence finie d'un maillage
∆x=∆y=1m et en supposant que le transfert de chaleur est stationnaire à deux dimensions,
déterminer les températures du début, milieu et du bas de la surface exposée du barrage.

1/4
www.almohandiss.com

Méthode de séparation des variables


Problème :( Session Normale : SMP4- juin 2008 )
On se propose d’étudier le comportement thermique d’une dalle épaisse de béton
soumise à un flux de chaleur. La dalle forme le plancher du rez-de-chaussée du bâtiment à
chauffer ; elle est en béton dense, de propriétés thermo physiques homogènes (conductivité
thermique, masse volumique, chaleur massique).
On s’intéresse ici au comportement thermique de la dalle d’épaisseur e, que l’on peut
caractériser en régime stationnaire par sa réponse à une densité de flux de chaleur ϕ0(t)
imposée sur sa face inférieure(x=0). Les dimensions transversales de la dalle sont très grandes
devant e, de sorte que le problème peut être supposé unidimensionnel. La direction
perpendiculaire à la surface S de la dalle est repérée par la cote verticale x à partir du fond
(Figure 1). La température de l’air ambiant dans le bâtiment est homogène, de valeur
moyenne Ti. La face inférieure de la dalle est parfaitement isolée thermiquement.
x
Air ambiant (h,Ti)
e

Dalle

ϕ(0,t)

0
Isolant Figure 1.

1. Rappeler la loi de Fourier unidimensionnelle liant la densité du flux de chaleur ϕ(x, t)


et la température T(x, t). Quelle est la dimension de ϕ ? En déduire l’unité de la conductivité
thermique λ.
2. En effectuant un bilan thermique de la tranche élémentaire [x, x + dx] de la dalle entre un
instant t quelconque et t + dt, établir les équations différentielles vérifiées par ϕ(x,t) et par
T(x,t).
λ d 2T ( x, t ) dT ( x, t )
3. En déduire que l’équation de diffusion thermique est = . Exprimer la
ρ .C dx 2 dt
diffusivité thermique a de la dalle, à l’aide des grandeurs thermo physiques de la dalle. Quelle
est l’unité de a ?
4. L’échange de chaleur par convection entre la face supérieure de la dalle et l’air ambiant est
bien représenté par la loi de Newton, ϕ(e, t) = h [T(e, t) − Ti] où h est considéré comme
constant et uniforme sur toute la surface. Quelle est l’unité de h ?
5. On suppose qu’un flux constant ϕ0 est maintenu en x = 0. Après un temps t très long, on
peut considérer le régime comme établi, de sorte que ϕ (x) et T(x) sont indépendants du
temps.
a) Montrer que le profil thermique de la dalle est alors donné par
⎛ e 1⎞ ϕ
T ( x) = Ti + ϕ 0 ⎜ + ⎟ − 0 x .
⎝λ h⎠ λ
b) Calculer l’écart de température T(e) − Ti, pour ϕ0 = 20,1 W· m−2 et h = 6, 7 SI.
λ=1,75 SI, e=0,26m, a=9.10-7SI.

2/4
www.almohandiss.com

6. En première approximation le flux de chaleur émis par le circuit solaire au fond de la dalle
peut être approché par une fonction sinusoïdale du temps de sorte que
ϕ (0, t ) = ϕ 0 (δ + sin(ω.t ) où δ est fonction du numéro du jour dans l’année, t le temps écoulé
depuis le lever du soleil, ω la pulsation du cycle jour-nuit de période J = 1 jour solaire (J = 24
heures).
Dans la suite des calculs, on se placera aux cas ou δ=0.
a) Quelle est la nature du régime thermique de la dalle ?
d 2T ( x, t ) 1 dT ( x, t )
b) Pour résoudre l’équation de diffusion thermique : = , on
dx 2 a dt
choisit la méthode de séparation des variables. Dans ce cas, montrer que cette
a d 2 g 1 df
équation peur être réécrite sous la forme suivante : = .
g dx 2 f dt
c) On admettant des solutions de formes exponentielles complexes comme

A. exp(i.ω.t ) et B. exp(i.k .x ) et sachant que i = −


(1 − i )2 , montrer que le système
2
ω
admet comme solution complexe k = ± (1 − i ) .
2a
d) On ne retient que la solution complexe pour k qui possède une partie réelle
ω
négative k = − (1 − i ) , justifier ce choix.
2a
e) Montrer que la solution finale complexe
⎡ ω ⎤ ⎡⎛ ω ⎞⎟⎤
est : T ( x, t ) = T0 . exp ⎢− x ⎥. exp ⎢i⎜⎜ ωt − x ⎥.
⎣ 2a ⎦ ⎣⎢ ⎝ 2a ⎟⎠⎦⎥
( )
f) Déduire la température réelle T ( x, t ) = Re T ( x, t ) et discuter le résultat.

3/4
www.almohandiss.com

Facteurs de forme de conduction de certains Systèmes Bidimensionnels φ = Fc λ (T1 − T2 )

Système Schéma Conditions Facteur de forme (Fc)

Sphère isotherme
enterré dans un milieu
semi-infini.

Cylindre horizontal
isotherme de longueur
L enterré dans un
milieu semi-infini

Cylindre vertical dans


un milieu semi-infini

Conduction entre deux


cylindres de longueur
infinie L en moyenne

Cylindre circulaire
horizontal de longueur
L à mi-chemin entre
les plans parallèles
d'égale longueur et
largeur infinie

Cylindre circulaire de
longueur L centré dans
un carré solide de
longueur égale

Excentrique cylindre
circulaire de longueur
L dans un cylindre de
longueur égale

Conduction par les


bords de murs

4/4
www.almohandiss.com

Solution :

Milieu thermiquement mince

A
V 2 L(a.b ) 2ro
Lc ,mur plan = = =L 2ro
S 2(a.b )
V πr h r
2
a
Lc ,cylindre = = o = o
S 2πro h 2
V 4πro / 3 ro
3

Lc , sphère = = =
4πro
2 b
S 3 2L

B
Hypothèses : 1 Le contenant en verre de forme cylindrique
Eau
d'un rayon de r0 = 3 cm. 2 Les propriétés thermiques du lait
60°C
sont prises comme ceux de l'eau. 3 Les propriétés thermiques
du lait sont constantes à la température ambiante. 4 Le
coefficient de transfert de chaleur est constant et uniforme sur Lait
toute la surface. 5 Le nombre de Biot, dans ce cas est 3°C
importante (beaucoup plus grand que 0,1). Toutefois, le
système global d'analyse est encore applicable puisque le lait
est agité en permanence, de sorte que sa température reste
uniforme en tout temps.
Analyse : La longueur caractéristique et le nombre de Biot du verre de lait sont :
V πro L 2
π (0.03 m) 2 (0.07 m)
Lc = = = = 0.01050 m
S 2πro L + 2πro 2 2π (0.03 m)(0.07 m) + 2π (0.03 m) 2
hLc (120 W/m 2 .°C)(0.0105 m)
Bi = = = 2.076 > 0.1
λ (0.607 W/m.°C)
Pour les raisons citées ci-dessus, nous pouvons utiliser le système global d'analyse afin
de déterminer le temps nécessaire pour que le lait se réchauffe de 3°C à 38°C.
hS h 120 W/m 2 .°C
b= = = = 0.002738 s -1
ρC pV ρC p Lc (998 kg/m 3 )(4182 J/kg.°C)(0.0105 m)
T (t ) − T∞ 38 − 60 -1
= e −bt ⎯
⎯→ = e −( 0.002738 s ) t ⎯
⎯→ t = 348 s = 5.8 min
Ti − T∞ 3 − 60
Par conséquent, il faudra environ 6 minutes pour réchauffer le lait de 3°C à 38°C.

Ailette
Hypothèses : 1 La température de la partie immergée de la cuillère est égale à la température
de l'eau. 2 La température de la cuillère varie le long de la cuillère T (x). 3 Le transfert de
chaleur à partir de l'extrémité de la cuillère est négligeable. 4 Le coefficient de transfert de
chaleur est constant et uniforme sur toute la surface de la cuillère. 5 Les propriétés
thermiques de la cuillère sont constantes. 6 Le transfert de chaleur par rayonnement et
supposé négligeable.

5/4
www.almohandiss.com

Notant que la section transversale de la cuillère est constante et x a pour origine la surface
libre de l'eau. La variation de température le long de la cuillère peut être exprimée comme :
T ( x ) − T∞ cosh a ( L − x )
= h, T∞ 
Tb − T∞ cosh aL
Où 0.2cm
Tb
p = 2(0.002 + 0.013) = 0.030 m 1.3cm
L=18cm
S c = (0.002)(0.013) = 0.000026 m 2

hp 17 * 0.03
a= = = 36 m -1
λ .S c 15 * 0.000026
La température de la pointe de la cuillère est déterminée par :
cosh a( L − L)
T ( L) = T∞ + (Tb − T∞ )
cosh aL
cosh 0 1
= 24 + (93 − 24) = 24° + (93 − 24) = 24.2°C
cosh(36 × 0.18) 326
Par conséquent, la différence de température dans l'ensemble du poignée de la cuillère
exposée est :
∆T = Tb − T (L) = (93 − 24.2) = 68.8°C
Facteur de forme
.
Hypothèse : le transfert de chaleur est en deux dimensions (pas de changement dans la
direction axiale).
La conductivité thermique du béton est λ= 0.75 W/m⋅°C.
Le facteur de forme pour cette configuration est donné dans le tableau :
2πL T1 = 60°C
Fc =
⎛ 2
4 z − D1 − D2
2 2

cosh −1 ⎜⎜ ⎟

T2 = 15°C
⎝ 2 D D
1 2 ⎠
2π (8 m)
= = 9.078 m
−1 ⎛ 4(0.4 m) − (0.05 m) − (0.05 m ) ⎞
2 2 2
cosh ⎜⎜ ⎟⎟
⎝ 2 ( 0 . 05 m )( 0 . 05 m ) ⎠ D = 5 cm
z = 40 cm
L=8m
Le flux de transfert de chaleur échangé entre les tuyaux est :
φ = Fc λ (T1 − T2 ) = (9.078 m)(0.75 W/m.°C)(60 − 15)°C = 306 W

Problèmes Bi-dimension
Nœud 5 sur la frontière isolée peut être considéré comme un nœud de l'intérieur qui donne :
Tdroit + Tgauche + Thaut + Tbas − 4Tnoeud = 0 .
Utilisation de la méthode du bilan d'énergie et en prenant la direction de tous les transferts de
chaleur vers le nœud, les équations aux différences finies pour les nœuds ont été obtenus
comme suit :
l T2 − T1 l/2
l
Nœud 1: hi (Ti − T1 ) + λ + [α s q& s + h0 (T0 − T1 )] = 0
2 2 l sin 45

6/4
www.almohandiss.com

l T1 − T2 l T4 − T2 T − T2
Nœud 2: hi l (Ti − T2 ) + λ +λ + λl 3 =0
2 l 2 l l
T −T T −T
Nœud 3: λl 2 3 + λl 5 3 +
l
[α s q& s + h0 (T0 − T3 )] = 0 1 • h o, T o
l l sin 45
l l T2 − T4 l T5 − T4 Eau 2 3 qs
Nœud 4: hi (Ti − T4 ) + λ +λ =0 • •
2 2 l 2 l hi, Ti
Nœud 5: T4 + 2T3 + T6 − 4T5 = 0
•4 •5 •6
l T5 − T6 l/2
Nœud 6: λ + [α s q& s + h0 (T0 − T6 )] = 0 Isolée
2 l sin 45
Où l = 1 m, λ = 0.6 W/m⋅°C, hi =150 W/m2⋅°C, Ti =15°C, ho = 30 W/m2⋅°C, T0 =25°C, αs =
0.7, et q&s = 800 W/m2 .
Le système de six équations à six inconnues constitue la formulation de la différence finie du
problème. Les températures des six nœuds sont déterminées par la résolution des six
équations ci-dessus.
T1 = Thaut =21.3°C, T2 =15.1°C, T3 = Tmilieu =43.2°C, T4 =15.1°C, T5 =36.3°C, T6 = Tbas
=43.6°C

Discussion Notez que la température la plus élevée se produit à un endroit plus éloigné de
l'eau, comme prévu.

7/4
www.almohandiss.com

Solution de l’épreuve
1. Rappelle de la loi de Fourier unidimensionnelle liant la densité du flux de chaleur ϕ(x,
t) et la température
dT ( x, t )
ϕ ( x, t ) = −λ
dx
Dimension de ϕ(W/m2)
N .m kg.m.s −2
= kg.s −3 → [ϕ ] = M .T −3
J N
= = =
s.m 2 s.m 2 s.m s.m
L’unité de la conductivité thermique λ
W
s
m2 = J
°C m 2 °C
2. Effectuant le bilan thermique de la tranche élémentaire [x, x + dx] de la dalle entre un
instant t quelconque et t + dt :
[ϕ x .s( x) − ϕ x+ dx s( x + dx)]dx = ρ .S ( x).dx.C Tt +dt + Tt
dt
Or s ( x) = s( x + dx) = s

[ϕ x − ϕ x+ dx ].S ( x).dx = ρ .S ( x).dx.C Tt +dt + Tt


dt
dT ( x, t )
dϕ ( x, t ) = − ρ .C
dt
3. L’équation de diffusion thermique :
dT ( x, t )
ϕ ( x, t ) = −λ
dx
d T ( x, t )
2
dT ( x, t ) λ d 2T (x, t ) dT (x, t )
−λ = − ρ .C =
dx 2 dt ρ .C dx 2 dt
La diffusivité thermique a:
λ
a=
ρ .C
L’unité de a est :
W .m 3 .Kg .K
[a] = = m 2 .s −1
m.K .Kg .J

4. L’unité de h :
W
m 2 .K

5. a) Le régime est stationnaire l’équation de diffusion thermique devient


λ d 2 T ( x, t )
=0
ρ .C dx 2
T ( x) = b + c.x
La solution de l’équation est :
dT ϕ
ϕ ( x ) = −λ = −λ.c = ϕ 0 c=− 0
À x=0 φ(x=0)=φ 0
dx d’ou λ

8/4
www.almohandiss.com

ϕ (e) = ϕ 0 = h.(T (e) − Ti )


À x=e
ϕ0
T (e) = + Ti
h
ϕ0
T ( e) = b − .e
λ
ϕ0 ϕ
+ Ti = b − 0 .e
h λ
ϕ0 ϕ0
+ e + Ti = b
h λ
ϕ ϕ ϕ
T ( x) = 0 + 0 e + Ti − 0 .x
h λ λ
⎛1 e ⎞ ϕ
T ( x) = Ti + ϕ 0 ⎜ + ⎟ − 0 .x
⎝h λ⎠ λ
b)
⎛1 e ⎞ ϕ ϕ
T (e) = Ti + ϕ 0 ⎜ + ⎟ − 0 .e = Ti + 0
⎝h λ⎠ λ h
ϕ 20.1
T (e) − Ti = 0 = = 3.1°C
h 6.7

6. a) le régime est périodique


b) pour résoudre l’équation de diffusion thermique :
d 2T ( x, t ) 1 dT ( x, t )
=
dx 2 a dt On choisit la méthode de séparation des variables.
On écrie :
T ( x, t ) = g ( x). f (t )
d T ( x, t )
2
d g ( x) 1 dT ( x, t ) 1
2
df (t )
2
= f (t ). 2
= = g ( x).
dx dx a dt a dt
2
d g ( x) 1 df (t )
f (t ). 2
= g ( x).
dx a dt
2
a d g ( x) 1 df (t )
2
=
g ( x) dx f (t ) dt
c) On admettant des solutions de formes exponentielles complexes comme :
f (t ) = A. exp(i.ω.t ) et g ( x) = B. exp(i.k .x )
a d 2 g ( x) 1 df (t )
2
=
L’équation g ( x) dx f (t ) dt devient :
1
B.i 2 .k 2 . exp(i.ω.x ) = A.i.ω. exp(i.ω.t )
a
B. exp(i.k .x ) A. exp(i.ω.t )

− a.k 2 = i.ω

Sachant que i=−


(1 − i)
2

9/4
www.almohandiss.com

a.k 2
=
(1 − i )2 .ω k 2
=
(1 − i ) 2


2 2.a
D’ou
ω
k=± .(1 − i )
2.a
d) Deux solution possible pour g(x)

⎛ ω ⎞ ⎛ ω ⎞⎟
g ( x) = B. exp⎜⎜ i.(1 − i ) .x ⎟⎟ = B. exp⎜⎜ (i + 1) x⎟
⎝ 2.a ⎠ ⎝ 2.a ⎠
⎛ ω ⎞⎟
Re[g ( x)] = B. exp⎜⎜ + x
⎝ 2.a ⎟⎠
⎛ ω ⎞⎟
lim Re [g ( x ) ] = lim B . exp ⎜
⎜ + x = +∞ pas de convergence
x →∞ x →∞ ⎝ 2.a ⎟⎠
Et
⎛ ω ⎞⎟ ⎛ ω ⎞⎟
g ( x) = B. exp⎜⎜ − i.(1 − i ) .x ⎟ = B. exp⎜⎜ − (i + 1) x⎟
⎝ 2.a ⎠ ⎝ 2.a ⎠
⎛ ω ⎞⎟
Re[g ( x)] = B. exp⎜⎜ − x⎟
⎝ 2 .a ⎠
⎛ ω ⎞⎟
lim Re [g ( x ) ] = lim B . exp ⎜
⎜ − x = 0 Convergence. On retient cette solution
x →∞ x →∞ ⎝ 2.a ⎟⎠
complexe qui possédant une partie réelle négative.

e) la solution finale est :


⎛ ω ⎞⎟ ⎛ ω ⎛ ω ⎞⎟ ⎞⎟
T ( x, t ) = g ( x). f (t ) = B. A. exp⎜⎜ − (i + 1) .x ⎟. exp(i.ω.t ) = B. A. exp⎜ − .x + i⎜⎜ ω.t − .x.⎟
2.a ⎜ 2.a 2.a ⎟
⎝ ⎠ ⎝ ⎝ ⎠⎠
À (x=0,t=0) :
T (0,0) = B. A = T0
⎛ ω ⎛ ω ⎞⎟ ⎞⎟
T ( x, t ) = T0 . exp⎜ − .x + i⎜⎜ ω.t − .x.

⎝ 2.a ⎝ 2.a ⎟⎠ ⎟⎠
f)

[ ] ⎛
T ( x, t ) = Re T ( x, t ) = T0 . exp⎜⎜ −
2
ω ⎞⎟ ⎛⎜
.a
.x ⎟. cos⎜ ω.t −
ω ⎞⎟
2.a
.x ⎟
⎝ ⎠ ⎝ ⎠

⎛ ω ⎞⎟ 2.π 2.π
exp⎜⎜ − .x ⎟ traduit l’amortissement, T = la période temporelle, Λ = la période
⎝ 2.a ⎠ ω ω
2.a
spatiale ou longueur d’onde et υ = 2.a.ω la vitesse de propagation.

10/4
www.almohandiss.com

Solution TD 3:

Milieu thermiquement mince

A
V 2 L(a.b ) 2ro
Lc ,mur plan = = =L 2ro
S 2(a.b )
2
V πr h r a
Lc ,cylindre = = o = o
S 2πro h 2
3
V 4πro / 3 ro
Lc , sphère = = 2
=
S 4πro 3 2L
b

B
Hypothèses : 1 Le contenant en v erre de form e cylindrique
Eau
d'un rayon de r 0 = 3 cm. 2 Les propriétés therm iques du lait
60°C
sont prises comme ceux de l'eau. 3 Les propriétés therm iques
du lait son t constantes à la températu re ambiante. 4 Le
coefficient de transfert de chaleu r est constant et uniforme sur Lait
toute la su rface. 5 Le nom bre de Biot, dans ce cas est 3°C
importante (beaucoup plus gran d que 0,1). Toutefois, le
système global d' analyse est encore applicable puisque le lait
est agité en permanence, de sorte que sa températu re reste
uniforme en tout temps.
Analyse : La longueur caractéristique et le nombre de Biot du verre de lait sont :
2
V πro L π (0.03 m) 2 (0.07 m)
Lc = = = = 0.01050 m
S 2πro L + 2πro 2 2π (0.03 m)(0.07 m) + 2π (0.03 m) 2
hLc (120 W/m 2 .°C)(0.0105 m)
Bi = = = 2.076 > 0.1
λ (0.607 W/m.°C)
Pour les raisons citées ci-dessus, nous pouvons utiliser le système global d'analyse afin
de déterminer le temps nécessaire pour que le lait se réchauffe de 3°C à 38°C.
hS h 120 W/m 2 .°C
b= = = = 0.002738 s -1
ρC pV ρC p Lc (998 kg/m 3 )(4182 J/kg.°C)(0.0105 m)
T (t ) − T∞ 38 − 60 -1
= e −bt ⎯
⎯→ = e −( 0.002738 s ) t ⎯
⎯→ t = 348 s = 5.8 min
Ti − T∞ 3 − 60
Par conséquent, il faudra environ 6 minutes pour réchauffer le lait de 3°C à 38°C.

Ailette
Hypothèses : 1 La température de la partie immergée de la cuillère est égale à la température
de l'eau. 2 La tem pérature de la cuillère varie le long de la cuillère T (x). 3 Le transfert de
chaleur à partir de l' extrémité de la cuillère est néglige able. 4 Le coefficient de transfert de
chaleur est constant et uniform e sur toute la surface de la cui llère. 5 Les propriétés
thermiques de la cuillère sont constantes. 6 L e transfert de chaleur par rayonnem ent et
supposé négligeable.
www.almohandiss.com

Notant que la section transversa le de la cuillère est constan te et x a po ur origine la surface


libre de l'eau. La variation de température le long de la cuillère peut être exprimée comme :
T ( x ) − T∞ cosh a ( L − x )
= h, T∞ 
Tb − T∞ cosh aL
Où 0.2cm
Tb
p = 2(0.002 + 0.013) = 0.030 m 1.3cm
2 L=18cm
S c = (0.002)(0.013) = 0.000026 m
hp 17 * 0.03
a= = = 36 m -1
λ .S c 15 * 0.000026
La température de la pointe de la cuillère est déterminée par :
cosh a( L − L)
T ( L) = T∞ + (Tb − T∞ )
cosh aL
cosh 0 1
= 24 + (93 − 24) = 24° + (93 − 24) = 24.2°C
cosh(36 × 0.18) 326
Par conséquent, la différence de tem pérature dans l'ensemble du poignée de la cuillère
exposée est :
∆T = Tb − T ( L) = (93 − 24.2) = 68.8°C
Facteur de forme
.
Hypothèse : le transfert de chaleur est en de ux dimensions (pas de changem ent dans la
direction axiale).
La conductivité thermique du béton est λ= 0.75 W/m⋅°C.
Le facteur de forme pour cette config uration est donné dans le tableau :
2πL T1 = 60°C
Fc =
⎛ 2 2
4 z − D1 − D2
2

cosh −1 ⎜⎜ ⎟

T2 = 15°C
⎝ 2 D D
1 2 ⎠
2π (8 m)
= = 9.078 m
−1 ⎛ 4(0.4 m) − (0.05 m) − (0.05 m ) ⎞
2 2 2
cosh ⎜⎜ ⎟⎟
⎝ 2 ( 0 . 05 m )( 0 . 05 m ) ⎠ D = 5 cm
z = 40 cm
L=8m
Le flux de transfert de chaleur échangé entre les tuyaux est :
φ = Fc λ (T1 − T2 ) = (9.078 m)(0.75 W/m.°C)(60 − 15)°C = 306 W

Problèmes Bi-dimension
Nœud 5 sur la frontière isolée peut être considéré comme un nœud de l'intérieur qui donne :
Tdroit + Tgauche + Thaut + Tbas − 4Tnoeud = 0 .
Utilisation de la méthode du bilan d'énergie et en prenant la direction de tous les transferts de
chaleur vers le nœud, l es équations aux diffé rences finies pour les nœuds ont été obtenus
comme suit :
l T2 − T1 l/2
l
Nœud 1: hi (Ti − T1 ) + λ + [α s q& s + h0 (T0 − T1 )] = 0
2 2 l sin 45
www.almohandiss.com

l T1 − T2 l T4 − T2 T − T2
Nœud 2: hi l (Ti − T2 ) + λ +λ + λl 3 =0
2 l 2 l l
T −T T −T
Nœud 3: λl 2 3 + λl 5 3 +
l
[α s q& s + h0 (T0 − T3 )] = 0 1 • h o, T o
l l sin 45
l l T2 − T4 l T5 − T4 Eau 2 3 qs
Nœud 4: hi (Ti − T4 ) + λ +λ =0 • •
2 2 l 2 l hi, Ti
Nœud 5: T4 + 2T3 + T6 − 4T5 = 0
•4 •5 •6
l T5 − T6 l/2
Nœud 6: λ + [α s q& s + h0 (T0 − T6 )] = 0 Isolée
2 l sin 45
Où l = 1 m, λ = 0.6 W/m⋅°C, hi =150 W/m2⋅°C, Ti =15°C, ho = 30 W/m2⋅°C, T0 =25°C, αs =
0.7, et q&s = 800 W/m2 .
Le système de six équations à six inconnues constitue la formulation de la différence finie du
problème. Les tem pératures des six nœuds sont déterm inées par la résolution des six
équations ci-dessus.
T1 = Thaut =21.3°C, T2 =15.1°C, T3 = Tmilieu =43.2°C, T4 =15.1°C, T5 =36.3°C, T6 = Tbas
=43.6°C

Discussion Notez que la température la plus élevée se produit à un endroit plus éloigné de
l'eau, comme prévu.
www.almohandiss.com

Solution de l’épreuve
1. Rappelle de la loi de Fourier unidimensionnelle liant la densité du flux de chaleur ϕ(x,
t) et la température
dT ( x, t )
ϕ ( x, t ) = −λ
dx
Dimension de ϕ(W/m2)
N .m kg.m.s −2
= kg.s −3 → [ϕ ] = M .T −3
J N
= = =
s.m 2 s.m 2 s.m s.m
L’unité de la conductivité thermique λ
W
s
m2 = J
°C m 2 °C
1. Effectuant le bilan thermique de la tranche élémentaire [x, x + dx] de la dalle entre un
instant t quelconque et t + dt :
[ϕ x .s( x) − ϕ x+ dx s( x + dx)]dx = ρ .S ( x).dx.C Tt +dt + Tt
dt
Or s ( x) = s( x + dx) = s

[ϕ x − ϕ x+ dx ].S ( x).dx = ρ .S ( x).dx.C Tt +dt + Tt


dt
dT ( x, t )
dϕ ( x, t ) = − ρ .C
dt
2. L’équation de diffusion thermique :
dT ( x, t )
ϕ ( x, t ) = −λ
dx
d T ( x, t )
2
dT ( x, t ) λ d 2T (x, t ) dT (x, t )
−λ = − ρ .C =
dx 2 dt ρ .C dx 2 dt
La diffusivité thermique a:
λ
a=
ρ .C
L’unité de a est :
W .m 3 .Kg .K
[a] = = m 2 .s −1
m.K .Kg .J

3. L’unité de h :
W
m 2 .K

4. a) Le régime est stationnaire l’équation de diffusion thermique devient


λ d 2 T ( x, t )
=0
ρ .C dx 2
T ( x) = b + c.x
La solution de l’équation est :
dT ϕ
ϕ ( x ) = −λ = −λ.c = ϕ 0 c=− 0
À x=0 φ(x=0)=φ 0
dx d’ou λ

1/4
www.almohandiss.com

ϕ (e) = ϕ 0 = h.(T (e) − Ti )


À x=e
ϕ0
T (e) = + Ti
h
ϕ0
T ( e) = b − .e
λ
ϕ0 ϕ
+ Ti = b − 0 .e
h λ
ϕ0 ϕ0
+ e + Ti = b
h λ
ϕ ϕ ϕ
T ( x) = 0 + 0 e + Ti − 0 .x
h λ λ
⎛1 e ⎞ ϕ
T ( x) = Ti + ϕ 0 ⎜ + ⎟ − 0 .x
⎝h λ⎠ λ
b)
⎛1 e ⎞ ϕ ϕ
T (e) = Ti + ϕ 0 ⎜ + ⎟ − 0 .e = Ti + 0
⎝h λ⎠ λ h
ϕ 20.1
T (e) − Ti = 0 = = 3.1°C
h 6.7

5. a) le régime est périodique


b) pour résoudre l’équation de diffusion thermique :
d 2T ( x, t ) 1 dT ( x, t )
=
dx 2 a dt On choisit la méthode de séparation des variables.
On écrie :
T ( x, t ) = g ( x). f (t )
d T ( x, t )
2
d g ( x) 1 dT ( x, t ) 1
2
df (t )
2
= f (t ). 2
= = g ( x).
dx dx a dt a dt
2
d g ( x) 1 df (t )
f (t ). 2
= g ( x).
dx a dt
2
a d g ( x) 1 df (t )
2
=
g ( x) dx f (t ) dt
c) On admettant des solutions de formes exponentielles complexes comme :
f (t ) = A. exp(i.ω.t ) et g ( x) = B. exp(i.k .x )
a d 2 g ( x) 1 df (t )
2
=
L’équation g ( x) dx f (t ) dt devient :
1
B.i 2 .k 2 . exp(i.ω.x ) = A.i.ω. exp(i.ω.t )
a
B. exp(i.k .x ) A. exp(i.ω.t )

− a.k 2 = i.ω

Sachant que i=−


(1 − i)
2

2/4
www.almohandiss.com

a.k 2
=
(1 − i )2 .ω k 2
=
(1 − i ) 2


2 2.a
D’ou
ω
k=± .(1 − i )
2.a
d) Deux solution possible pour g(x)

⎛ ω ⎞ ⎛ ω ⎞⎟
g ( x) = B. exp⎜⎜ i.(1 − i ) .x ⎟⎟ = B. exp⎜⎜ (i + 1) x⎟
⎝ 2.a ⎠ ⎝ 2.a ⎠
⎛ ω ⎞⎟
Re[g ( x)] = B. exp⎜⎜ + x
⎝ 2.a ⎟⎠
⎛ ω ⎞⎟
lim Re [g ( x ) ] = lim B . exp ⎜
⎜ + x = +∞ pas de convergence
x →∞ x →∞ ⎝ 2.a ⎟⎠
Et
⎛ ω ⎞⎟ ⎛ ω ⎞⎟
g ( x) = B. exp⎜⎜ − i.(1 − i ) .x ⎟ = B. exp⎜⎜ − (i + 1) x⎟
⎝ 2.a ⎠ ⎝ 2.a ⎠
⎛ ω ⎞⎟
Re[g ( x)] = B. exp⎜⎜ − x⎟
⎝ 2 .a ⎠
⎛ ω ⎞⎟
lim Re [g ( x ) ] = lim B . exp ⎜
⎜ − x = 0 Convergence. On retient cette solution
x →∞ x →∞ ⎝ 2.a ⎟⎠
complexe qui possédant une partie réelle négative.

e) la solution finale est :


⎛ ω ⎞⎟ ⎛ ω ⎛ ω ⎞⎟ ⎞⎟
T ( x, t ) = g ( x). f (t ) = B. A. exp⎜⎜ − (i + 1) .x ⎟. exp(i.ω.t ) = B. A. exp⎜ − .x + i⎜⎜ ω.t − .x.⎟
2.a ⎜ 2.a 2.a ⎟
⎝ ⎠ ⎝ ⎝ ⎠⎠
À (x=0,t=0) :
T (0,0) = B. A = T0
⎛ ω ⎛ ω ⎞⎟ ⎞⎟
T ( x, t ) = T0 . exp⎜ − .x + i⎜⎜ ω.t − .x.

⎝ 2.a ⎝ 2.a ⎟⎠ ⎟⎠
f)

[ ] ⎛
T ( x, t ) = Re T ( x, t ) = T0 . exp⎜⎜ −
2
ω ⎞⎟ ⎛⎜
.a
.x ⎟. cos⎜ ω.t −
ω ⎞⎟
2.a
.x ⎟
⎝ ⎠ ⎝ ⎠

⎛ ω ⎞⎟ 2.π 2.π
exp⎜⎜ − .x ⎟ traduit l’amortissement, T = la période temporelle, Λ = la période
⎝ 2.a ⎠ ω ω
2.a
spatiale ou longueur d’onde et υ = 2.a.ω la vitesse de propagation.

3/4
www.almohandiss.com

TD de Transfert de Chaleur.
Série N°4 : Rayonnement thermique

I) Utilisation des abaques.

Pour la configuration illustrée sur la Figure, avec S3 la surface plane annulaire entre les deux
cylindres. Déterminer les expressions F13, F31, F32 et F23 en fonction de F11, F12 et les trois surfaces S1,
S2 et S3. F33 est le facteur de forme entre les deux surfaces annulaires opposées. S1 est la surface
intérieure du cylindre extérieur et S2 la surface externe du cylindre intérieur.
Pour L=20cm, r1=20cm et r2=10cm, déterminer F11, F12 et F13 .

II) Règle de Hottel

S1 + S 2 − S3
Démonter que : F12 =
2S1

R  S2 S
III) Méthode des cordes croisées : Facteur de forme entre 2 surfaces non
concaves quelconques (sans obstacle)

Démonter que : F12 =


(S PR + SQS ) − (S RQ + S PS )
2.S1
P
Q  S1
IV) Application de la méthode des cordes croisées : Deux cylindres
parallèles infiniment longs d'un diamètre D sont situés sur une distance d D 
l’un de l'autre. Déterminer le facteur de forme F12 entre ces deux
cylindres. D 
(2) 

(1) 
V) une lampe de 100 W est alimentée sous 220 V. Elle est constituée d 
d’un filament de tungstène placé au centre d’une ampoule sphérique de 8
cm de diamètre à l’intérieur de laquelle on fait le vide.
Pour que la lumière soit assez blanche, il est nécessaire que la température du filament soit de 2600 K.
1) Déterminer le diamètre et la longueur du filament si le facteur d’émission total
hémisphérique du tungstène est de 0.3 (résistivité du tungstène 88 µΩcm).
2) Déterminer la puissance rayonnée dans le visible (entre 0.4 et 0.7 µm) si le facteur
d’émission spectral hémisphérique est de 0.45 dans ce domaine.
3) Quelle est la puissance absorbée par l’ampoule en supposant que le verre est parfaitement
transparent jusqu’à 2.7 µm et se comporte comme un corps noir au-delà, en admettant que
le facteur d’émission spectral hémisphérique du tungstène est de 0.2 dans tout le domaine
au-delà 2.7 µm.
4) Déterminer la température de l’ampoule si le facteur d’émission total hémisphérique du
verre est de 0.93 en négligeant les pertes de chaleur par convection naturelle.

SMP4 TD N°4 : Rayonnement thermique  Page 1 

 
www.almohandiss.com

VI) Flux net d’une surface grise


Soit trois surfaces planes (1,2 et 3) grises de dimensions semi infini fermées (voir figure) :
- Calculer les facteurs de forme suivant (voir figure) : F11, F12, F13, F21, F22, F23, F31, F32, F33.
- Donner la radiosité de chaque surface.
- Donner la définition du flux net d’une surface 1.
- Calculer le flux net de la surface 1.

ε3=0.5,  T3=500°C, L3=0.4m 
ε2=0.5,  T2=500°C, L2=0.3m 

ε1=0.15,  T1=100°C, L1=0.5m 

VII) Deux boucliers fins d’émissivités ε3=0.10 et ε4=0.15 T1 = 600 K


sont, des deux côtés, placés entre deux très grandes plaques parallèles,
qui sont maintenues aux températures uniformes T1=600K et à T2=300K
et ont les émissivités ε1=0.6 et ε2=0.7, respectivement.
ε3 = 0.10 
Déterminez les flux nets échangé entre les deux plaques (1 et 2) avec
et sans les boucliers par unité de surface des plaques. Déterminez les
températures des boucliers T3 et T4 à l’équilibre thermique. ε4 = 0.15
T2 = 300 K

VIII) Un thermocouple protégé par le papier d'aluminium


d’émissivité ε2=0.15 est utilisé pour mesurer la température de gaz  
chauds coulant dans un conduit dont les murs sont maintenus à
Tw=380 K. Le thermomètre montre une lecture de température
Tth=530 K.
En supposant l’émissivité de la jonction de thermocouple ε1=0.7 et le
coefficient de transfert de chaleur de convection h=120 W/m2 °C, Thermocouple 
déterminez la température réelle du gaz. Qu’il est la lecture du Air, Tf   Tth = 530 K 
thermomètre si aucune protection de radiation n'a été utilisée ? ε1 = 0.7 
Tw = 380 K  ε2 = 0.15 

SMP4 TD N°4 : Rayonnement thermique  Page 2 

 
www.almohandiss.com

Facteur de forme - Deux cylindres concentriques de longueur finie. 

SMP4 TD N°4 : Rayonnement thermique  Page 3 

 
www.almohandiss.com

Fraction d’énergie F0-λT rayonnée par un corps noir entre 0 et λ

SMP4 TD N°4 : Rayonnement thermique  Page 4 

 
www.almohandiss.com

Solution

I)

II)  

S1 F12 + S1 F13 = S1 (F11 = 0)


S 2 F21 + S 2 F23 = S 2 (F22 = 0)  Conservation de l’énergie (système géométriquement fermé) 
S3 F31 + S3 F32 = S3 (F33 = 0)

SMP4 TD N°4 : Rayonnement thermique  Page 5 

 
www.almohandiss.com

 
S1 F12 = S 2 F21 , S1 F13 = S 3 F31 , S 2 F23 = S3 F32 Réciprocité 

S1 + S 2 − S3
F12 =  
2S1

III)

Pour PQRS (enceinte fermée S1 convexe) : S1 F12 + S1 F1− PR + S1 F1−QS = S1

S1 + S PR − S RQ
Pour PQR (plan de fermeture) : S1 F1− PR =
2

S1 + SQS − S PS
Pour PQS (plan de fermeture) : S1 F1−QS =
2

S RQ + S PS − (S PR + SQS )
D' où F12 =
2.S1

IV) Le facteur de forme entre les deux cylindres en face l’un de l'autre pour d / D> 3 est 
déterminée par à l’utilisation de la méthode des cordes croisées.

F1− 2 =
∑ diag − ∑ cotés  
2Si

2 d 2 + D 2 − 2d
F1− 2 =  
2(πD / 2)

D’ou  F1− 2 =
(
2 d 2 + D2 − d
 
)
πD

SMP4 TD N°4 : Rayonnement thermique  Page 6 

 
www.almohandiss.com

V) 
1)  Déterminer le diamètre et la longueur du filament si le facteur d’émission total hémisphérique du 
tungstène est de 0.3 (résistivité du tungstène 88 µΩcm). 
M T0 = ε Tu .σ .T 4 = 0.3 * 5.67 10 −8 * 2600 4 = 7.77 10 5 Wm −2  

φ = S * M T0 = 100W  
100 100
S = π .D.L = =
M T 7.7710 5  
0

D.L = 4.1 10 −5 m 2
La résistance électrique du fil est : 
V 2 220 2
Re = = = 484 Ω  
P 100
L
Re = ρ = 484 Ω
π .D 2
4  
L
2
= 43197 10 4 m −1
D
D’où : 
L = 0.898m
 
D = 4.56 10 −5 m
2) la puissance rayonnée dans le visible (entre 0.4 et 0.7 µm) si le facteur d’émission spectral 
hémisphérique est de 0.45 dans ce domaine : 
 
0.7
⎛ 0.7 0.4

PVisible = S * ∫ M λ0.T dλ = S * ε visible ⎜⎜ ∫ M λ0.T dλ − ∫ M λ0.T dλ ⎟⎟ = S * ε visible (F0−0.7.T − F0−0.4.T )σ .T 4  
0.4 ⎝0 0 ⎠
0.05
λ1 = 0.4 µm et T = 2600 K → λ1 .T = 0.4 * 2600 = 1040µm.K → F0−0.4T = = 0.0005
100  
4.18
λ2 = 0.7 µm et T = 2600 K → λ2 .T = 0.7 * 2600 = 1820µm.K → F0−0.7T = = 0.0418
100
 
PVisible = π * D * L * ε visible (0.0418 − 0.0005)5.67 10 −8 * 2600 4
 
PVisible = π * 4.56 10 −5 * 0.898 * 0.45 * (0.0418 − 0.0005)5.67 10 −8 * 2600 4 = 6.2 W
Le calcul peut être fait autrement : 
PTotale = ε * S * σ * T 4 = P = 100 W
La puissance totale rayonnée par le fil est :   P 100  
S *σ *T 4 = =
ε 0.3

SMP4 TD N°4 : Rayonnement thermique  Page 7 

 
www.almohandiss.com

ε Visible
PVisible = S * ε visible (F0−0.7.T − F0−0.4.T )σ .T 4 = P (F0−0.7.T − F0−0.4.T )  
ε
 
3) la puissance absorbée par l’ampoule en supposant que le verre est parfaitement transparent 
jusqu’à 2.7 µm et se comporte comme un corps noir au‐delà, en admettant que le facteur d’émission 
spectral hémisphérique du tungstène est de 0.2 dans tout le domaine au‐delà 2.7 µm : 
ε Inf
PInf = S * ε visible (1 − F0− 2.7.T )σ .T 4 = P (1 − F0−2.7.T )  
ε
80.97
λ3 = 2.7 µm et T = 2600 K → λ3 .T = 2.7 * 2600 = 7020µm.K → F0− 2.7T = = 0.8097  
100
ε Inf
PInf = P (1 − F0−2.7.T ) = 100 * 0.2 (1 − 0.8097) = 12.7 W  
ε 0.3
4) Régime stationnaire : la puissance absorbée par le verre est égale à la puissance 
rayonnée par le verre vers l’extérieure, d’où : 
SVerre * ε Verreσ * TVerre
4
= PInf = 12.3 W   

Avec εVerre facteur d’émission  total hémisphérique du verre 
12.3 W 12.3
4
.TVerre = =
SVerre * ε Verre * σ 0.93 * 5.67 10 −8 * 4 * π * 4 10 − 2 ( )
2
 
.T 4
Verre = 1.16 10 K
10
⇒ Tverre = 328 K
 
VII) 
Hypothèses :   
‐  Les surfaces sont opaques, diffuse, et grises.  
T1 = 600 K
‐ Convection de transfert de chaleur n'est pas pris en considération.   ε1 = 0.6 
 Les émissivités des surfaces sont : ε1 = 0.6, ε2 = 0.7, ε3 = 0.10 et ε4 = 0.15.  
T2 = 300 K 
Le flux net entre les surfaces 1 et 2 sans boucliers est donné par :
ε2 = 0.7 

1 − ε1 1 1− ε2
S1ε 1 S1 F12 S 2ε 2  

σ (T1 4 − T2 4 )
φ12,sans boucliers =  
1 − ε1 1 1− ε2
+ +
S1ε 1 S1 F12 S 2 ε 2
Or F12=1, S1=S2=S 
D’où : 

SMP4 TD N°4 : Rayonnement thermique  Page 8 

 
www.almohandiss.com

σ (T1 4 − T2 4 )
ϕ12,sans boucliers =
1 1
+ −1
ε1 ε2
(5.67 × 10 −8 W/m 2 ⋅ K 4 )[(600 K ) 4 − (300 K ) 4 ]
=  
1 1
+ −1
0.6 0.7
= 3288 W/m 2
Le flux net entre les surfaces 1 et 2 avec boucliers est donné par :

  T1 = 600 K
ε1 = 0.6 
σ (T1 4 − T2 4 )
ϕ12,avec boucliers =
⎛1 1 ⎞ ⎛1 1 ⎞ ⎛ 1 1 ⎞
⎜⎜ + − 1⎟⎟ + ⎜⎜ + − 1⎟⎟ + ⎜⎜ + − 1⎟⎟ ε3 = 0.10 
⎝ ε1 ε 2 ⎠ ⎝ ε3 ε3 ⎠ ⎝ε4 ε4 ⎠
−8
(5.67 × 10 W/m ⋅ K )[(600 K ) − (300 K ) 4 ]
2 4 4
= ε4 = 0.15
⎛ 1 1 ⎞ ⎛ 1 1 ⎞ ⎛ 1 1 ⎞
⎜ + − 1⎟ + ⎜ + − 1⎟ + ⎜ + − 1⎟ T2 = 300 K
⎝ 0.6 0.7 ⎠ ⎝ 0.10 0.10 ⎠ ⎝ 0.15 0.15 ⎠ ε2 = 0.7 
= 206 W/m 2
 
La température d'équilibre thermique de boucliers sont déterminés à partir de : 
σ (T1 4 − T3 4 )
ϕ13 =
⎛1 1 ⎞
⎜⎜ + − 1⎟⎟
⎝ ε1 ε 3 ⎠  
4
(5.67 × 10 −8 W/m 2 ⋅ K 4 )[(600 K ) 4 − T3 ]
206 W/m 2 = ⎯
⎯→ T3 = 549 K
⎛ 1 1 ⎞
⎜ + − 1⎟
⎝ 0.6 0.10 ⎠
σ (T4 4 − T2 4 )
ϕ 42 =
⎛ 1 1 ⎞
⎜⎜ + − 1⎟⎟
  ⎝ε4 ε2 ⎠  
4
(5.67 × 10 −8 W/m 2 ⋅ K 4 )[T4 − (300 K ) 4 ]
206 W/m 2 = ⎯
⎯→ T4 = 429 K
⎛ 1 1 ⎞
⎜ + − 1⎟
⎝ 0.15 0.7 ⎠
 
VIII La température des gaz chauds dans un conduit est mesurée par un thermocouple. La 
température du gaz est à déterminer et  sans rapport avec un écran de protection. 

SMP4 TD N°4 : Rayonnement thermique  Page 9 

 
www.almohandiss.com

Thermocouple 

Tth = 530 K 
Air, Tf  
ε1 = 0.7 
Tw = 380 K 
ε2 = 0.15 
 
Hypothèses : Les surfaces sont opaques, diffuse et  grises. 
Le flux net échangé entre la surface du thermocouple et la surface du conduit est : 
 
1 − ε1 1 1 − ε 2 ,1 1 − ε 2,3 1 1− ε3
S 1ε 1 S 1 F12 S 2 ε 2 ,1 S 2 ε 2,3 S 2 F2 ,3 S 3ε 3  

σ (Tth 4 − TW 4 )
φThermocouple − w, Avec Ecran =
1 − ε1 1 1 − ε 2,1 1 − ε 2,3 1 1− ε3
+ + + + +
S1ε 1 S1 F12 S 2 ε 2,1 S 2 ε 2,3 S 2 F2,3 S 3ε 3
 
S1 surface de la sonde
S1 ≈ S 2 surface de l ' écran
S 3 surface du conduit
Notant que F12=1, F23=1 et S1/S3 tend ver zéro
D’où : 
φThermocouple _ W , Avec Ecran σ (TTh 4 − TW 4 )
ϕ Thermoucouple−écran = =
S1 ⎛1⎞ ⎛ 1 ⎞
⎜⎜ ⎟⎟ + ⎜⎜ 2 − 1⎟⎟
⎝ ε1 ⎠ ⎝ ε 2 ⎠
(5.67 × 10 −8 W/m ⋅ K )[(530 K ) − (380 K ) 4 ]
2 4 4
=  
⎛ 1 ⎞ ⎛ 1 ⎞
⎜ ⎟ + ⎜2 − 1⎟
⎝ 0.7 ⎠ ⎝ 0.15 ⎠
= 239.2 W/m 2
A l’équilibre thermique la densité de flux échangée par le thermocouple par convection est 
égale à celle échangée par rayonnement. 
ϕ convction,Thermocouple = ϕ rayonnement, thermocouple
  h(Tgaz − Tth ) = 239.2 W/m 2  
120 W/m 2 ⋅ °C(Tgaz − 530) = 239.2 W/m 2 ⎯
⎯→ Tgaz = 531.99 K

SMP4 TD N°4 : Rayonnement thermique  Page 10 

 
www.almohandiss.com

Sans écran du thermocouple la température du gaz est : 

Thermocouple 
Air, Tf  

Tw = 500 K  Tth = 850 K 

 
 
ε thσ (Tth 4 − Tw 4 )
Tgaz = Tth +
h  
(0.7)(5.67 × 10 −8 W/m 2 ⋅ K 4 )[(530 K ) 4 − (380 K ) 4 ]
= 530 K + = 549.2 K
120 W/m 2 ⋅ °C

SMP4 TD N°4 : Rayonnement thermique  Page 11 

 
www.almohandiss.com

VI 
• Calculer les facteurs de forme suivant (voir figure) : F11, F12, F13, F21, F22, F23, F31, F32, F33.
F11= F22= F33=0 surfaces planes 
S1 + S 2 − S3
Application de la Règle de Hottel : F12 =
2S1

S1 + S 2 − S 3 L1 + L2 − L3 0.5 + 0.3 − 0.4 0.4


F12 = = = = = 0.4
2 S1 2 L1 2 * 0.5 1

S1 + S 3 − S 2 L1 + L3 − L2 0.5 + 0.4 − 0.3 0.6


F13 = = = = = 0.6
2 S1 2 L1 2 * 0.5 1

S 2 + S 3 − S 21 L2 + L3 − L1 0.3 + 0.4 − 0.5 0.2 1


F23 = = = = =
2S 2 2 L2 2 * 0.3 0.6 3

L2 0.3 1
S 2 F23 = S 3 F3.2 ⇒ F32 = F23 = = 0.25
L3 0.4 3

L1 0.5 2
S 2 F21 = S1 F1.2 ⇒ F21 = F12 = 0.4 =
L2 0.3 3

L1 0.5
S 3 F31 = S1 F1.3 ⇒ F31 = F13 = 0.6 = 0.75
L3 0.4

• Radiosité des différentes surfaces :

J 1 = ϕ1 + (1 − ε 1 )[F11 J 1 + F12 J 2 + F13 J 3 ]  


M 30
J 2 = ϕ 2 + (1 − ε 2 )[F21 J 1 + F22 J 2 + F23 J 3 ]  
1− ε3
 
J 3 = ϕ 3 + (1 − ε 3 )[F31 J 1 + F32 J 2 + F33 J 3 ]   S 3ε 3
J3 
 
1 1
  S1 F13 S 2 F23
  1 − ε1 1 1− ε 2
0 S1ε 1 S 2ε 2
M 1
J1  S1 F12 J2  M 20
 

SMP4 TD N°4 : Rayonnement thermique  Page 12 

 
www.almohandiss.com

• Définition du flux net d’une surface 1 : 

Le flux net de la surface 1 est égal à la différence entre celui le flux émis et celui absorbé.

S1ε 1
φ1,nette =
1 − ε1
( )
M 10 − J 1  

• Calcule du flux net de la surface 1.

⎛ − ϕ1 ⎞ ⎡ −1 (1 − ε 1 )F12 (1 − ε 1 )F13 ⎤⎛ J 1 ⎞
⎜ ⎟ ⎢ ⎜ ⎟
⎜ − ϕ 2 ⎟ = ⎢(1 − ε 2 )F21 −1 (1 − ε 2 )F23 ⎥⎥⎜ J 2 ⎟  
⎜ − ϕ ⎟ ⎢ (1 − ε )F (1 − ε 3 )F32 − 1 ⎥⎦⎜⎝ J 3 ⎟⎠
⎝ 3⎠ ⎣ 3 31

⎛ − ε 1σ * T14 ⎞ ⎡ − 1 0.34 0.51⎤⎛ J 1 ⎞


⎜ ⎟ ⎢ 1 1 ⎥⎜ ⎟
⎜ − ε 2σ * T2 ⎟ = ⎢ −1
4
⎜ J2 ⎟ 
⎜ − ε σ *T 4 ⎟ ⎢ 3 6 ⎥⎥⎜ ⎟
⎝ 3 3 ⎠ ⎣0.375 0.125 − 1 ⎦⎝ J 3 ⎠

⎛ 164 ⎞ ⎡ 1 − 0.34 − 0.51⎤⎛ J 1 ⎞


⎜ ⎟ ⎢ 1 1 ⎥⎜ ⎟
⎜16466 ⎟ = ⎢ − 1 − ⎥⎜ J 2 ⎟  
⎜16466 ⎟ ⎢ 3 6 ⎜ ⎟
⎝ ⎠ ⎣− 0.375 − 0.125 1 ⎥⎦⎝ J 3 ⎠

[K ] = [M ][J ]  
[M ]−1 [K ] = [M ]−1 [M ][J ] = [I ][J ] = [J ]  
⎛ J 1 = 25536 ⎞
⎜ ⎟
Après inversion de la matrice :  ⎜ J 2 = 29934 ⎟  
⎜ J = 29785 ⎟
⎝ 3 ⎠

φ1,nette ε1
S1
=
1 − ε1
(M 10 − J 1 ) = −4477 Wm − 2  

La surface 1 reçoit de la chaleur.

SMP4 TD N°4 : Rayonnement thermique  Page 13 

 
www.almohandiss.com

Solution de L’exercice VI de du TD N°4 
• Calculer les facteurs de forme suivant (voir figure) : F11, F12, F13, F21, F22, F23, F31, F32, F33.
F11= F22= F33=0 surfaces planes 
S1 + S 2 − S3
Application de la Règle de Hottel : F12 =
2S1

S1 + S 2 − S 3 L1 + L2 − L3 0.5 + 0.3 − 0.4 0.4


F12 = = = = = 0.4
2 S1 2 L1 2 * 0.5 1

S1 + S 3 − S 2 L1 + L3 − L2 0.5 + 0.4 − 0.3 0.6


F13 = = = = = 0.6
2 S1 2 L1 2 * 0.5 1

S 2 + S 3 − S 21 L2 + L3 − L1 0.3 + 0.4 − 0.5 0.2 1


F23 = = = = =
2S 2 2 L2 2 * 0.3 0.6 3

L2 0.3 1
S 2 F23 = S 3 F3.2 ⇒ F32 = F23 = = 0.25
L3 0.4 3

L1 0.5 2
S 2 F21 = S1 F1.2 ⇒ F21 = F12 = 0.4 =
L2 0.3 3

L1 0.5
S 3 F31 = S1 F1.3 ⇒ F31 = F13 = 0.6 = 0.75
L3 0.4

• Radiosité des différentes surfaces :

J 1 = ϕ1 + (1 − ε 1 )[F11 J 1 + F12 J 2 + F13 J 3 ]  


M 30
J 2 = ϕ 2 + (1 − ε 2 )[F21 J 1 + F22 J 2 + F23 J 3 ]   1− ε3
 
S 3ε 3
J 3 = ϕ 3 + (1 − ε 3 )[F31 J 1 + F32 J 2 + F33 J 3 ]  
J3 
 
1 1
  S1 F13 S 2 F23

  1 − ε1 1 1− ε 2
0 S 2ε 2
M 1
S1ε 1 J1  S1 F12 J2  M 20
 

SMP4 TD N°4 : Rayonnement thermique  Page 1 

 
www.almohandiss.com

• Définition du flux net d’une surface 1 : 

Le flux net de la surface 1 est égal à la différence entre celui le flux émis et celui absorbé.

S1ε 1
φ1,nette =
1 − ε1
( )
M 10 − J 1  

• Calcule du flux net de la surface 1.

⎛ − ϕ1 ⎞ ⎡ −1 (1 − ε 1 )F12 (1 − ε 1 )F13 ⎤⎛ J 1 ⎞
⎜ ⎟ ⎢ ⎜ ⎟
⎜ − ϕ 2 ⎟ = ⎢(1 − ε 2 )F21 −1 (1 − ε 2 )F23 ⎥⎥⎜ J 2 ⎟  
⎜ − ϕ ⎟ ⎢ (1 − ε )F
⎝ 3⎠ ⎣ 3 31 (1 − ε 3 )F32 − 1 ⎥⎦⎜⎝ J 3 ⎟⎠

⎛ − ε 1σ * T14 ⎞ ⎡ − 1 0.34 0.51⎤⎛ J 1 ⎞


⎜ ⎟ ⎢ 1 1 ⎥⎜ ⎟
⎜ − ε 2σ * T2 ⎟ = ⎢
4
−1 ⎜ J2 ⎟ 
⎜ − ε σ *T 4 ⎟ ⎢ 3 6 ⎥⎥⎜ ⎟
⎝ 3 3 ⎠ ⎣0.375 0.125 − 1 ⎦⎝ J 3 ⎠

⎛ 164 ⎞ ⎡ 1 − 0.34 − 0.51⎤⎛ J 1 ⎞


⎜ ⎟ ⎢ 1 1 ⎥⎜ ⎟
⎜16466 ⎟ = ⎢ − 1 − ⎥⎜ J 2 ⎟  
⎜16466 ⎟ ⎢ 3 6 ⎜ ⎟
⎝ ⎠ ⎣− 0.375 − 0.125 1 ⎥⎦⎝ J 3 ⎠

[K ] = [M ][J ]  
[M ]−1 [K ] = [M ]−1 [M ][J ] = [I ][J ] = [J ]  
⎛ J 1 = 25536 ⎞
⎜ ⎟
Après inversion de la matrice :  ⎜ J 2 = 29934 ⎟  
⎜ J = 29785 ⎟
⎝ 3 ⎠

φ1,nette ε1
S1
=
1 − ε1
(M 10 − J 1 ) = −4477 Wm − 2  

La surface 1 reçoit de la chaleur.

SMP4 TD N°4 : Rayonnement thermique  Page 2 

Vous aimerez peut-être aussi