Microsoft PowerPoint - TDde Thermique - td2021

ENSAJ TD
Exercices de transferts de chaleur
Table des matières

• Conduction
– conduction en régime permanent
– conduction en régime variable
• Convection
• Rayonnement
• Problèmes combinés
2
I. Conduction
en régime permanent
Prof. EL OTMANI 1
ENSAJ TD
Qu ’est qu’un régime permanent?
Baignoire
Pompe
Conduction: qu’est qu’il faut savoir?

• Concerne principalement les solides
• Mode de transmission:
• vibration moléculaire
• vibration du réseau
• déplacement d’électrons libres
• Lois fondamentales
Générales 1D
Q=qA
Q = Qin − Qout = −  q.next dA T f − Ti
q = −k
q = − k ∇T e
QV + k∆T = ρ c Tɺ QV + k
d ²T
= ρ c Tɺ
dx ²
5
C1: Loi de Fourier
Soit une plaque de conductibilité k =34W/mK, de

surface A = 0,2m² et d ’épaisseur e = 4cm.
La température de la face supérieure de la plaque est
égale à 380K; celle de la face inférieure vaut
360K.
On demande de calculer le débit de chaleur dans

cette plaque.
Prof. EL OTMANI 2
ENSAJ TD
q
Tsup
k
Tinf
∂T T − T sup − 20
• Gradient : = inf = = − 500 K / m
∂n e 0 , 04
∂T
• Flux : q = −k = 34 . 500 = 17 kW / m ²
∂n
• Débit : Q = q A = 17 . 0 , 2 . 10 3 = 3400 W = 3, 4 kW
C2 Loi de Fourier / milieu Composite

Soit une paroi composée d ’une plaque de cuivre de
conductibilité kc=372W/mK et d’épaisseur ec=3mm placé
entre deux plaques identiques en inox, de conductibilité
kI=17W/mK et d’épaisseur eI=2mm. Les contacts entre ces
différentes plaques sont parfaits.
La température de la face gauche de cette paroi est Tg=400°c,

tandis que celle de droite vaut Td=100°c.
On demande de déterminer les températures des deux faces

de la plaque de cuivre, et de vérifier les résultats par
l ’analogie électrique.
8
Inox Cuivre Inox
Tg=400 °c Td=100 °c
kI kC kI
eI eC eI
q I = qC = q I
1 2
Conservation du flux:
∆TI ∆TC
avec qI = −k I et qC = − kC
eI eC
 ∆TI = 145,03°c TC g = 254,97°c qI = 1,232 106W

∆Ttot = Tg − Td = 300°c 
∆TC = 9,94°c TC = 245,03°c qC = 1,232 106W
d
Prof. EL OTMANI 3
ENSAJ TD
C4
– Une paroi composée de 2 milieux homogènes de
conductibilité k1 et k2, et d’épaisseur L1 et L2
respectivement.
– Les conditions aux limites sont les suivantes:
• à gauche: T = T0
= h(T − T f ) (flux de conduction sortant du

dT
• à droite: − k
dx
domaine 2 = flux de convection
dans le fluide)
• contact parfait entre les deux milieux.
– On demande la distribution de T dans chacun des
milieux et le flux thermique, ainsi que le contrôle par
analogie électrique.
10
k1 = 0,8 W / cm °c
k1 k2 h
k 2 = 0,2 W / cm °c
h = 0,5 W / cm 2 °c
T0 = 600 °c
T0 Tf
T f = 200 °c
L1 = 4 cm
L2 = 3 cm
L1 L2 A = 1 cm 2 (aire transversale)
11
QV = 0 pas de source ni de puit

QV + k∆T = ρ c Tɺ  ɺ
 T = 0 régime établi
• a) Solution directe:
milieu 1: d 2T1 T1 = a x + b
=0
dx 2
milieu 2: d 2 T2
=0 T2 = c x + d
dx 2
• Conditions aux limites:

en x = 0 : T1 ( x = 0 ) = T0 b = T0 = 600°c
en x = L1 : T1 ( x = L1 ) = T2 (x = L1 ) a L1 + T0 = c L1 + d
en x = L1 : q1 ( x = L1 ) = q2 ( x = L1 ) − a k1 = −c k 2
dT
en x = L1 + L2 : − k 2 2
dx L1 + L2
(
= h T2 L + L − T f
1 2
) [
− k 2 c = h c (L1 + L2 ) + d − T f ]
12
Prof. EL OTMANI 4
ENSAJ TD
• Solution:
T1 = −22,73 x + 600
T2 = −90,91 x + 872,73 q1 = −k1a = 18,18 W / cm²
q2 = −k 2 c = 18,18 W / cm²
T1 (L1 ) = T2 (L1 ) = 509,1°c
qconv = h(236,36 − 200) = 18,18 W / cm²
T2 (L1 + L2 ) = 236,36°c
• Graphiquement:
T (°c)
700
600
509,1°c
Tf = 200°c
500
400
300
236,4°c
200
100
x (cm)
0 1 2 3 4 5 6 7 8
13
Contact non parfait

x
Gradient de
température
TiI TiII T
qi = qi
I II
qi = hc Ti − Ti ( I II
)
14
Variantes:
a) Au lieu d ’un contact parfait entre les deux
milieux, on considère le cas d ’une séparation
entre les milieux par une conductance
d ’interface cint = 0,1W/cm²°c
Les conditions en x = L1 sont maintenant remplacées par:
q1 ( x = L1 ) = qint = cint (T1 (L1 ) − T2 (L1 ))

15
Prof. EL OTMANI 5
ENSAJ TD
C ’est à dire: − a k1 = cint (a L1 + b − c L1 − d )

− a k 2 = cint (a L1 + b − c L1 − d )
On trouve:
T1 = 600 − 15,625 x
T1 (L1 ) = 537,5°c
q1 = 12,5W / cm²
T2 = 662,5 − 62,5 x
q2 = 12,5W / cm²
T2 (L1 ) = 412,5°c
T2 (L1 + L2 ) = 225°c
q1 = q2 = qint = 0,1(537,5 − 412,5) = 12,5W / cm²

= qconv = 0,5(225 − 200) = 12,5W / cm²
Remarque: Résistance entre les deux milieux -> flux de chaleur plus petit.
16
Variantes:
b) La première paroi est le siège d ’une distribution
de sources volumiques Qv = 50 W/cm³.
Les conditions aux limites sont:
> condition de Dirichlet à gauche T = T0
> condition de Newton (convection) à droite
= h(T − T f )
dT
−k
dx
> contact parfait entre les deux milieux.
On demande les distributions de T et de flux
thermique dans chacun des milieux, et le contrôle
17
par analogie électrique.
Solution directe:
+ QV = 0(= ρ cTɺ )
• Milieu 1 d 2T1
k1
dx 2
Q x²
T1 = T1 + T1 = − V + a x + b
part hom
2k1
dT1
q1 = −k1 = QV x − a k1
dx
• Milieu 2 T2 = c x + d
d ²T2
=0 dT
dx ² q2 = − k 2 2 = − k 2 c
dx
4 constantes à identifier a, b, c, d.
18
Prof. EL OTMANI 6
ENSAJ TD
Conditions aux limites: (4)

1) x = 0 T1 ( x = 0) = T0  b = T0
2) x = L1
− QV L1 ²
T1 ( x = L1 ) = T2 (x = L1 )  + a L1 + T0 = c L1 + d
2k1
3) x = L1
q1 ( x = L1 ) = q2 ( x = L1 )  QV L1 − a k1 = −c k2
4) x = L1 + L2
− k2
dT2
dx L1 + L2
(
= h T2
L1+ L2
− Tf ) [
 − c k 2 = h c (L1 + L2 ) + d − T f ]
19
QV x ²
Solutions: T1 = − + ax + b = −31.25 x ² + 198.86 x + 600 (°c)
2k1
q1 = QV x − ak1 = 50 x − 159.09 (W / cm ²)
T2 = cx + d = −204.55 x + 1713.65 (°c)
q2 = −ck 2 = 40.91 (W / cm²)
T1 ( L1 ) = T2 (L1 ) = 895 .45 °c Temp . de jonction

 159 .09 
T1 max = T1 dT1 = 0 = T1  x =  = 916 .36 °c max de temp .
dx  50 
T paroi = T2 (L1 + L2 ) = T1 ( x = 7 ) = 281 .8 °c
q1 ( x = 0) = −159.09W/cm²
q1 ( x = 0) ≠ q1 ( x = L1 )
q1 ( x = L1 ) = q2 ( x = L1 ) = 40.91W/cm²
( )
q2 ( x = L1 + L2 ) = h TL1 + L2 − T f = 0.5 (281.8 − 200) = 40.9W/cm²
20
Analogie électrique
Électricité Thermique
i Q
V T
Rélectrique Rthermique
i Q
Vi Vf Ti Tf
Rélec Rtherm
Vi − V f = i Rélectrique Ti − T f = Q Rthermique
en un noeud :  i j = 0 en un noeud :  Q j = 0
j j
21
Prof. EL OTMANI 7
ENSAJ TD
Électricité Thermique
Source de courant? Source volumique?
Injecteur de courant i Injecteur de chaleur Q
Remarque: Analogie électrique --->> discrétisation
Calcul des résistances

thermiques
Ti − T f T f − Ti Ti − T f
En conduction: = Q = Aq = − Ak∇T ∇T = =−
R L L
L
R=
kA
(
En convection: q = h Ti − T f ) R=
1
hA
22
C2 Loi de Fourier / milieu Composite

Soit une paroi composée d ’une plaque de cuivre
de conductibilité kc=372W/mK et d ’épaisseur
eI=2mm. Les contacts entre ces différentes
plaques sont parfaits.
La température de la face gauche de cette paroi est
Tg=400°c, tandis que celle de droite vaut
Td=100°c.
On demande de déterminer les températures des
deux faces de la plaque de cuivre, et de vérifier
les résultats par l ’analogie électrique.
23
Inox Cuivre Inox
Tg=400 °c Td=100 °c
kI kC kI
eI eC eI
q I = qC = q I
1 2
Conservation du flux:
∆TI ∆TC
avec qI = −k I et qC = −kC
eI eC
 ∆TI = 145,03°c TC g = 254,97°c qI = 1,232 106W

∆Ttot = Tg − Td = 300°c 
∆TC = 9,94°c TC = 245,03°c qC = 1,232 106W
d
24
Prof. EL OTMANI 8
ENSAJ TD
C2: Confirmation par l’analogie électrique

RI RC RI
Tg Tcg Tcd Td
∆T e
Q = qA = R=
R kA
Tg − Tc Tg − Td
g
=
RI RT
RT = 2 RI + RC OK
Tc − Td Tg − Td
d
=
RI RT
25
C4
– Une paroi composée de 2 milieux homogènes de
conductibilité k1 et k2, et d’épaisseur L1 et L2
respectivement.
– Les conditions aux limites sont les suivantes:
• à gauche: T = T0
= h(T − T f ) (flux de conduction sortant du

dT
• à droite: − k
dx
domaine 2 = flux de convection
dans le fluide)
• contact parfait entre les deux milieux.
– On demande la distribution de T dans chacun des
milieux et le flux thermique, ainsi que le contrôle par
analogie électrique.
26
k1 = 0,8 W / cm °c
k1 k2 h
k 2 = 0,2 W / cm °c
h = 0,5 W / cm 2 °c
T0 = 600 °c
T0 Tf
T f = 200 °c
L1 = 4 cm
L2 = 3 cm
L1 L2 A = 1 cm 2 (aire transversale)
27
Prof. EL OTMANI 9
ENSAJ TD
• C4: Contrôle par l ’analogie électrique
T0 TL1 TL1+L2 Tf
R1=L1/k1A R2=L2/k2A Rext=1/hA
T0 − TL1 TL1 − T f
nœud TL1: = TL1 = 509 ,1°c
L1 L2 1
+
k1 A k 2 A hA
T − TL1 + L2 TL + L − T f
nœud TL1+L2: 0 = 1 2 TL1 + L2 = 236 ,36 °c
L1 L2 1
+
k1 A k 2 A hA
28
Variantes:
a) Au lieu d ’un contact parfait entre les deux
milieux, on considère le cas d ’une séparation
entre les milieux par une conductance
d ’interface cint = 0,1W/cm²°c
Les conditions en x = L1 sont maintenant remplacées par:

29
• Variante par l ’analogie électrique:

Ta Tb1 Tb2 Tc Tf
Rab1 Rb1b2 Rb2c Rcf
L1 1 L2 1
Rab1 = Rb1b2 = Rb2c = Rcf =
k1 A cA k2 A hA
Raf = Rab1 + Rb1b2 + Rb1c + Rcf Rab2 = Rab1 + Rb1b2 ...
Les équations aux nœuds s ’écrivent:
Ta − Tb1 Ta − T f
= Tb1 = 537,5°c
Rab1 Raf
Ta − Tb2 Tb1 − Tb2
= Tb2 = 412,5°c
Rab2 Rb1b2
Ta − Tc Ta − T f
= Tc = 225°c
Rac Raf
30
Prof. EL OTMANI 10
ENSAJ TD
Variantes:
b) La première paroi est le siège d ’une distribution
de sources volumiques Qv = 50 W/cm³.
Les conditions aux limites sont:
> condition de Dirichlet à gauche T = T0
> condition de Newton (convection) à droite
= h(T − T f )
dT
−k
dx
> contact parfait entre les deux milieux.
On demande les distributions de T et de flux
thermique dans chacun des milieux, et le contrôle
31
par analogie électrique.
Variante b: contrôle par l’analogie électrique

R1=L1/k1A R2=L2/k2A Rext=1/hA
T0 TL1 TL1+L2 Tf
Q = Qv L1 A/2
Q = Qv L1 A
TL1 − T0 TL1 − T f
QLA
– en TL1: + − v 1 =0 TL1 = 895,45°c
R1 Rext + R2 2
TL1 + L2 − T0 TL1 + L2 − T f
– en TL1+L2: + = 0 TL1 + L2 = 281,8°c
R1 + R2 Rext
32
T (°c) Milieu Milieu

1 Tmax 2
1000
895.45°c
Qv = 50
800
600
509.1°c
Qv = 0
400 TP = 281.8°c
Tf = 200°c
200
Q (W/m²) 1 2 3 4 5 6 7 x (cm)
40 q2 = 40.91
20 q = 18.18
0
-20
-40
-60
-80 q1
-100
-120
-140 x (cm)
-160
33
Prof. EL OTMANI 11
ENSAJ TD
C6: Les ponts thermiques
Soit une paroi relativement isolante traversée par

des tirants métalliques bons conducteurs disposés
avec un pas horizontal 2b et un pas vertical 2h.
Quelle est l ’influence de ces ponts thermiques sur

le transfert conductif à travers la paroi?
Rem: on considérera l ’écoulement 1D suivant l ’axe des

tirants
34
kp
d
Ti km Te 2h
L 2b
A0 = b h Section d ’échange « témoin »
π d²
Ap = Section occupée par le pont thermique au sein de A0
16
ε = Ap A0 Coefficient d ’occupation du pont par rapport à A
Am = (1 − ε )A0 Section de la paroi non perturbée au sein de A0

L
R0 =
Résistances thermiques: k m A0 mur sans tirant
L
Rm =
k m Am mur avec tirant
L
Rp = tirant
k p Ap
35
Débit de chaleur Rm
Q0
sans tirant T −T T −T Ti Te T Te
Q0 = i e = i e k m A0 i
R0 L R0 Rp
avec tirant Ti − Te Ti − Te Ti − Te
Qeff =
Rm
+
Rp
=
L
[
A0 (1 − ε )k m + ε k p ]
Influence des ponts thermiques
Qeff kp  kp  peut être très différent de 1
= (1 − ε ) + ε = 1 + ε  − 1
Q0 km k
 m  car kp/km>>1
p.ex. k = 43W / mK (acier )  k
p p
= 1075
k m = 0.04W / mK (isolant ) km
Qeff
= 1 + 1074 ε 5 10 −4 1.537 (54% de perte)
Q0
ε = 1 10−3 Qeff / Q0 = 2.074 (107% de perte)
⋮ ⋮ ⋮
N.B. On suppose dans le modèle simplifié que l’écoulement reste 1D malgré
l’inhomogénéité provoquée par les tirants. En réalité, il y a un transfert de chaleur
latéral dont on n ’a pas tenu compte
36
Prof. EL OTMANI 12
ENSAJ TD
Analogie électrique - suite

• Tube cylindrique (écoulement radial)
dT k Ti − Te
r q = −k =
ln  dx r  re  r 
T − Ti r ln  ln  e 
=  i  ri  ∆T r
Te − Ti  re  T −T R= =  i
ln  Q = q 2π r l = 2π k l i e Q 2π k l
 ri  r 
ln e 
 ri 
• Paroi sphérique (écoulement radial)
dT k Ti − Te
q = −k =
1 1 dx r ² 1 − 1 1 1
− −
T − Ti ri r ri re ∆T ri re
= R= =
Te − Ti 1 − 1 T −T
Q = q 4π r ² = 4π k i e Q 4π k
ri re 1 1
−
ri re
37
C7. Influence d ’un isolant sur le refroidissement

d ’une résistance électrique
Soit la section 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 micanite (mélange de mica et de résine phénolique
agissant comme isolant électrique et thermique). On demande 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 ’énergie électrique dissipée dans la résistance est perdue
par convection - rayonnement avec l ’ambiance à 20°c avec un
coefficient d ’échange h=17W/m²K. La conductibilité thermique
de la micanite est k=0,1W/mK. Les caractéristiques électriques de
la résistance sont : PE=1W, RE=106ohms. Sa longueur est de
30mm et son diamètre est de 1mm.
Quelle est la température à l ’interface graphite - verre?

Comparer les résultats avec et sans isolant.
38
Micanite
r2=1mm
Te T3 T2 Verre
Graphite
= 5.88 103 m = r3 (cf . th. p 71)

kenv 0.1
rcr = =
hext 17
épaisseur de micanite assurant le maximum de déperdition
Qconv + ray = 0.40 PE = 0.4W
( N .B. PE = I ² RE = I ² 10 6  I = 10 −3 A faible courant )
 ln r3 r2 1 
T2 − Te = Qc −v  + 
 2π k env L 2 π rext hext L 
= 0.4(93.984 + 53.073) + 58.8  T2 ≅ 79°c
sans isolant: T − T = Q 1
c −v = 0.4 x312 = 124.8  T2 ≅ 145°c
2 π rext hext L
2 e
39
Prof. EL OTMANI 13
ENSAJ TD
I. Conduction
en régime variable
40
Conduction en régime variable

• Théorie
• Analyse trans. où la conduction est négligeable
– cas du chauffage convectif
• Problème du thermocouple
• Cas du changement de phase
– cas du chauffage convectif avec source interne
– cas du chauffage par flux superficiel constant
• Distribution de temp. dans un milieu semi-infini
– concept de profondeur de pénétration thermique
– Étude du régime trans. à travers une paroi d ’ép. Fini
• Discrétisation en régime variable
41
• Analogie électrique
Théorie
• Équation de base:
∂  ∂T 
QV +  k ij  = ρ c ∂T
∂xi  ∂x  ∂t
 j 
Sources/puits Accumulation/
conduction déperdition
• Cas du milieu isotrope homogène:

QV ∂T
+ α ∆T = Δ : Laplacien
ρc ∂t
k
α= diffusivité thermique
ρc
42
Prof. EL OTMANI 14
ENSAJ TD
Analyse transitoire où la conduction est négligeable

• Milieu caractérisé par une seule température
Q = ρ c Tɺ V
Qin − Qout + Qgen = −  q.next dA +  QV dV

A V
– A: surface externe du milieu (surface d’échange effectif)

– V: volume du milieu
43
Cas du chauffage convectif

Q = Qin =  h (T∞ − T ) dA = h (T∞ − T ) A si h, T uniformes
A ∞
spatialement
h T∞
T(t) h LC hV plausibilité de
Bi = = < 0.1
V k kA l’analyse moyenne
A
Bilan de chaleur
ρ c V Tɺ + h A (T − T∞ ) = 0
1
posons θ = T − T∞  θɺ + θ = 0
τ
τ = ρ cV
1
=
[
JK m m −1
=
[J ] = [s ]
−3 3
] temps caract. = capacitance
[
h A Wm −2 K −1m 2 [W ] ] therm. x résistance therm. superficielle
44
Solution:
−t T − T∞ −t
θ = θ0 e τ ou =e τ
T0 − T∞
θ T
T∞
θ0
T0
τ t τ t
V
N.B. = t h A = t h
t
LC = Longueur caractéristique
τ ρ c V ρ c LC A
hL kt θ
= C = B F = e − Bi F0
k ρ c LC 2 θ
i 0
0
45
Prof. EL OTMANI 15
ENSAJ TD
C10.1 Problème du thermocouple

Soit une plaque métallique (e=2.5cm, Ad=5x5 cm²)
de caractéristiques thermiques ρ c=1.2 JK-1cm-3,
k=8Wcm-1K-1.
Ad e
Initialement à T0=500°c et plongé dans une ambiance

T∞=1000°c.
Quelle est la réponse thermique de la plaque utilisée

comme thermocouple si h=0.6Wcm-2K-1 et si on
néglige le transfert à travers l’épaisseur?
46
• Hypothèse: - chauffage convectif

- conduction dans la plaque négl. (Bi < 0.1)
- Température uniforme dans le barreau
• V = 25 x 2.5 cm³ τ=
1.2 x 5
= 2 .5 s
Aéchange = 20 x 2.5 cm² 0 .6 x 4
 T = T∞ + (T0 − T∞ ) e −0.4 t = 1000 − 500 e −0.4 t

t = 0 (s) 1 1,277 2 3 4 5 6 7
T = T0 = 500 (°c) 664,8 700 775,3 849,3 899 932,3 954,6 969,6
• i.e. Plausibilité de l’étude moyenne

h V 0 .6 5
Bi = = = 0.09375 (limite!!!)
kA 8 4
47
C10.2 Cas du changement de phase

Supposons qu’à Tc=700°c, la plaque précédente
entre en fusion en consommant une chaleur latente
de changement d ’état ρ L =70JK-1cm-3.
Comment la réponse thermique sera-t-elle affectée si

on contient le liquide dans le même volume
qu’initialement, et que celui-ci conserve les
caractéristiques thermiques ρ c initiales?
48
Prof. EL OTMANI 16
ENSAJ TD
Lors du changement de phase, on a

ρ ceff = ρ c + ρ L δ (T − Tc )
(aug. apparente de la capa. calorifique en T=Tc)
et Hɺ = Q variation d’enthalpie = chaleur apportée
si le changement de phase se déroule de ti à tf , la

température ayant atteint Tc=700°c, on aura
tf tf
H f − H i =  Q dt soit Vρ L =  h(T∞ − Tc ) Aech dt
ti ti
Dès lors le changement de phase aura la durée,

ρ LV 70 x 5
∆t = t f − ti = = = 0.486 s
h (T∞ − Tc ) Aech 0.6 x 300 x 4
49
Vu les hyp. , le thermocouple donnera la réponse

T = 1000 − 500 e −0.4t pour 0 < t < 1.277s
T = 700 pour 1.277s < t < 1.763s
T = 1000 − 500 e −0.4 (t − 0.486) pour t > 1.763s
Tinf =1000. τ = 2.5s Réponse sans
950. changement de
900. phase
T
(°c) 850.
800.
750.
Réponse avec
Tc = 700.
changement de
∆t = 0.49s
650. phase
600.
550.
t (s)
T0 = 500°c
1.00 2.00 3.00 4.00 5.00 6.00 7.00
50
C10.3 Détermination d’un coefficient d’échange moyen

Soit un barreau d’aluminium de forme cylindrique (diamètre
2.5cm) initialement à une température uniforme de
T0=0°c (le volume du bain est suffisant pour que la
présence du barreau n’affecte pas sensiblement sa
température).
5 minutes après le début de l’opération, on relève pour le

barreau une température de section droite (T supposée
uniforme et ne variant pas avec la longueur) de 5°c.
Évaluer le coefficient d’échange métal-bain en supposant

que l’échange se fait par une loi de Newton?
q = h(T∞ − T )
51
Prof. EL OTMANI 17
ENSAJ TD
En la vertu de l’étude transitoire moyenne

−t soit t θ
θ = θ0 e τ = ln 0
τ θ
ρ c V ρ c V θ
or τ= h= ln 0
hA tA θ
pour l’aluminium k = 230 Wm −1 K −1
V π R 2l R
ρ = 2700 kg m −3 = =
A 2π R l 2
c = 900 J kg −1 K −1
40 − 0
h=
0.9 x 2.7 6
5 x 60
(
10 0.625 10 −2 ln
5−0
)= 105.3Wm − 2 K −1
plausibilité de l’analyse moyenne

h LC 105.3 x R 2 OK
Bi = = = 2.86 10 −3 << 0.1
k 230
52
Cas du chauffage convectif avec source interne

Q = Qin + Qgén =  h (T∞ − T ) dA +  QV dV si h, T∞ et QV uniformes
= h (T∞ − T ) A + QV V
A V
spatialement
= ρ c Tɺ V
De même qu’auparavant, posons

1 Q
θ = T − T∞  θɺ + θ = V
τ ρc
solution homogène: θ h = ε e −t /τ
solution particulière: θ = QV τ = QV V = QV LC
p
ρc Ah h
or θ (0) = θ 0
QV LC QV LC
θ = ε e − t /τ + ε =θ 0 −
h h
53
θ = θ 0 e − t /τ +
QV LC
h
(
1 − e −t / τ )
Terme transitoire relatif à la Terme forcé provenant de l ’existence
condition initiale (influence de la source interne (influence
diminuant avec le temps) grandissant avec le temps)
QV LC QV V
limθ = = ∀θ 0 régime permanent équilibré
t →∞ h hA
Tɺ = 0  h(T∞ − Téq ) A + QV V = 0
aussi Q V
T∞ − Téq = V
hA
54
Prof. EL OTMANI 18
ENSAJ TD
θ0
θ (t )
QV V
équilibre
hA
θ 0 e − t /τ
t
QV V
hA
(
1 − e − t /τ )
Q V
N.B. si θ 0 = V θ = θ 0 toujours car T0 = Téq
hA
QV V
θ0 > θ diminue avec t
hA
QV V
θ0 < θ augmente avec t
hA
55
Cas du chauffage par flux superficiel constant

Q = Qin =  qe dA = qe A si qe uniforme
A
ρ c Tɺ V = qe A qe A
T0
ρ cV
t
qe A
T = T0 + t augmentation non bornée de T
ρ cV
(on n’atteindra jamais un état d’équilibre dans ces

conditions)
56
Distribution de température dans un milieu

semi-infini en régime variable
Concept de profondeur de pénétration thermique
Soit un milieu dont la conduction interne n’est pas négligeable
Si apports ou pertes de chaleur à ses frontières
→ échauffement ou refroidissement commençant par la frontière
→ gradient de température
→ variation de température dans le temps
Zone d’influence du phénomène limitée dans un espace défini (proche
des frontières au départ) pour une période de temps définie
Cette zone d ’influence est la profondeur de pénétration ou couche
limite (interne) de conduction thermique
57
Prof. EL OTMANI 19
ENSAJ TD
C11. Cas d ’un milieu semi-infini siège d’un écoulement mono

dimensionnel de chaleur selon x (cas d’un plan d’extension infinie
limitant le demi-espace: x est la coordonnée normale au plan).
Le milieu est initialement à la température uniforme T0 et on porte
soudainement la température du plan à la valeur Ts (≠T0), qui est
maintenue pour les instants ultérieurs.
L’équation de la chaleur s’écrit:

x
∂ T ∂T
2
α =
∂x 2 ∂t TS T0
k
α= ( diffusivité ) QV = 0 T = T ( x, t )
ρc
58
Les conditions initiales et aux limites sont:
T (x,0 ) = T0 x>0
T (0, t ) = Ts t >0
lim x→∞ T ( x, t ) = T0 t >0 Pas de perturbation
à grande distance
Si on pose:
T − T0
θ=
TS − T0
Ces conditions deviennent:
θ ( x ,0 ) = 0 x > 0
θ (0, t ) = 1 t > 0
lim x→∞ θ ( x, t ) = 0 t > 0
59
L’équation à satisfaire devient:
∂ 2θ ∂θ
T α =
∂x 2 ∂t
TS
T0
x
x ≅ δ (t )
En t > 0 chaleur diffusée dans le milieu jusqu ’à une profondeur
de δ(t) à partir de la frontière
T (0, t ) = Ts
T (δ (t ), t ) = T0
La température à une profondeur x sera fonction du rapport x/δ(t)

60
Prof. EL OTMANI 20
ENSAJ TD
Remplacer la dépendance en les deux variables x et t, par une

dépendance vis-à-vis d’une seule variable η=x/δ(t)
 η = η2
 T ( x1 , t1 ) = T (x2 , t 2 ) si x1 δ (t1 ) = x2 δ (t 2 ) puisque  1
θ = θ (η )
TS
t2
t1
T0
x1 x2 x
61
Solution: Disparition de x (x = η.δ):

∂θ dθ ∂η 1 dθ d 2θ dθ
= = α + η δ δɺ =0
∂x dη ∂x δ dη dη 2 dη
∂ 2θ 1 d 2θ
= Disparition de la dép. temp.:
∂x 2 δ 2 dη 2
d 2θ
∂θ d θ ∂η x δɺ dθ α
= =− 2 dη 2
∂t d η ∂t δ dη = −δ δɺ = constante ∆ − C
dθ
η
L’équation adimentionnelle dη
de la chaleur devient ainsi Premier membre n’a pas de dép.
α d 2θ x δɺ dθ explicite en t
=− 2
δ 2 dη 2 δ dη Second membre dép. explicitement
de t puisque δ= δ(t)
On cherche à faire disparaître de
l’équation toute dépendance Ces deux membres sont donc tous
explicite en x ou en t pour que η deux égaux à une constante,
soit la seule variable indépendante appelée -C
62
L’équation à satisfaire devient:
 δ δɺ = C ⇔
1 d 2
( )
δ =C α
d 2θ
+ Cη
dθ
=0
2 dt dη 2 dη
Équation linéaire du second ordre
 δ 2 = 2 C .t + C0 Posons p = dθ , on a: α dp + Cη p = 0
dη dη
avec δ (t = 0 ) = 0  C0 = 0 dp C
Soit = − η dη
p α
il reste donc : δ = 2 C.t C 2

 ln p = − η + ln C1
2α
63
Prof. EL OTMANI 21
ENSAJ TD
Notons maintenant que le choix de la constante C est arbitraire:

Si C devient C’, alors
C′ C′
δ devient δ ′ = 2 C ′.t = δ η devient η ′ = η
C C
Et l’équation différentielle est:
d 2θ dθ
α + C ′η ′ =0
dη ′2 dη ′
dθ C ′ dθ d 2θ C ′ d 2θ
= et =
dη ′ C dη dη ′2 C dη 2
Ainsi l’équation devient:
C ′ d 2θ C C ′ dθ
α +η C′ =0
C dη 2 C′ C dη
d 2θ dθ Équation différentielle
Ou encore: α +η C =0
dη 2 dη de départ
64
Le choix de la constante C étant donc arbitraire, on choisit:
dθ
p = C1 exp (− η 2 ) =
C
=1 
2α dη
η
Et donc: θ = C2 + C1  exp(− u 2 )du

0
En η = 0, θ = 1  C2 = 1
∞
En η = ∞, θ = 0  0 = 1 + C1  exp − u 2 du ( )
0
π
On retrouve ici une intégrale de Poisson qui vaut , donc C1 = − 2
2 π
Il vient finalement: η
 exp(− u )du = 1 − erf (η ) = erfc(η )

2
θ = 1− 2
π 0
Où erfc(η) est la fonction erreur de Gauss

(table 2.2, pC36)
65
 exp(− u )du = 1 − erf (η ) = erfc(η )

2
θ = 1− 2
π 0
Il est utile de noter que:
puisque δ (t ) = 2 C t = 4 α t
x x
η= =
δ (t ) 4α t
La profondeur de pénétration δ est donc proportionnelle à la racine
carrée du produit de la diffusivité thermique et du temps
erfc(η)
1
0,8
0,6
0,4
0,2
η
0
0 1 2
66
Prof. EL OTMANI 22
ENSAJ TD
The complementary Error Function

η
erfc (η ) = 1 −  exp (− u ) du
2 2
π 0
éta erfc (éta) éta erfc (éta) éta erfc (éta) éta erfc (éta)
0 1 0,55 0,4367 1,1 0,1198 2,1 0,00298
0,05 0,9436 0,6 0,3961 1,2 0,08969 2,2 0,00186
0,1 0,8875 0,65 0,358 1,3 0,06599 2,3 0,001143
0,15 0,832 0,7 0,32222 1,4 0,04772 2,4 0,000689
0,2 0,7773 0,75 0,2889 1,5 0,0339 2,5 0,000407
0,25 0,7237 0,8 0,2579 1,6 0,02365 2,6 0,000236
0,3 0,6714 0,85 0,2293 1,7 0,01621 2,7 0,000134
0,35 0,6206 0,9 0,2031 1,8 0,01091 2,8 0,000075
0,4 0,5716 0,95 0,1791 1,9 0,00721 2,9 0,000041
0,45 0,5245 1 0,1573 2 0,00468 3 0,000022
0,5 0,4795
67
Notons que la profondeur de pénétration effective est matière à

interprétation: on décide d’une tolérance en-deçà de laquelle θ peut être
considéré comme nul (l’information venant de x=0 est considérée
comme négligeable au-delà de cet xlim).
pour xlim = δ → et erfc(1) = 0,157

On obtient environ 16% de la variation en surface.
Abaisser la tolérance conduit aux valeurs suivantes:
Tolérance 9% erfc(ηlim ) = 0,09 ηlim ≅ 1,2 xlim = 1,2 δ = 2,4 α .t

68
Application:
Déterminer la profondeur d’enfouissement d’une canalisation dans le
sol pour que des conditions de gel pendant 48heures (Tair=-10°c) ne
provoquent pas l’obstruction de celle-ci.
Les conditions nominales avant la période de gel sont Tair=Tsol=5°c
air Tair
k sol = 0,52Wm −1 K −1
sol sec
d ρ sol = 2050 kg m −3
csol = 1840 J kg −1 K −1
Tc
TC − T0 0−5
= = erfc(ηlim )  ηlim = 0,68
1
θ= =
Ts − T0 − 10 − 5 3
 k  0,52
α sol =   = = 1,38 10 −7 m 2 s −1
 ρ c  sol 2050 1840
d lim = η lim δ (t = 48h ) = 0,68 4α sol (48 3600) = 0,21m au moins
69
Prof. EL OTMANI 23
ENSAJ TD
C.12 Étude du régime transitoire à travers

une paroi d’épaisseur finie
A. Conditions pour pouvoir appliquer le concept de profondeur de
pénétration à une plaque d’épaisseur finie
T ( x,0) = T0 T (0, t ) = TS t >0
x
Ts 1
1° phase δ <L paroi ≡ milieu infini
2
t > tlim
2° phase δ >L paroi finie
t < tlim
T0
δ face
arrière
L isolée
70
On obtient le temps limite d’observation de la paroi finie à partir de la

solution du milieu 1/2 infini par:
A) choix d’une tolérance ε telle que si θ(x=L) ≤ ε, la solution du
milieu semi-infini est utilisable, par exemple ε=0,01.
θ lim = ε = erfc(ηlim )  η lim = 1,8 =

B) L
δ (tlim )
L L2
C) 2 α t lim = tlim = 0,077
1,8 α
Au-delà de tlim, on doit passer à l’étude de la paroi finie
N.B. tlim est fonction des conditions aux limites; ainsi si la même
paroi présentait 2 faces T=TS pour t>0, le temps limite
deviendrait tlim/4 car l ’épaisseur caractéristique serait divisée
par 2.
71
B. Étude de la paroi finie sollicitée par chauffage convectif

T
T∞ x T∞ x t∞
t>0 k
α= t2 > t1
ρc
t1
h T0
T0 face t0 face
arrière arrière
L isolée k/h L isolée
La paroi est initialement à T0 uniforme. En t=0, on porte

soudainement sa face gauche en contact avec un fluide à température
T∞ (différent de T0) avec lequel elle échange de la chaleur par
convection via un coefficient d ’échange h connu. Sa face droite est
isolée à tout moment.
72
Prof. EL OTMANI 24
ENSAJ TD
Notons ψ = T − T∞ écart de température entre paroi et fluide

∂ 2ψ 1 ∂ψ ∂ψ ∂T
L’équation à résoudre s’écrit = car =
∂x 2 α ∂t ∂t ∂t
Avec les conditions aux limites ∂ψ =0

paroi isolée ou axe de
symétrie
∂x x =0
(pas de flux au travers)
∂ψ h convection
=− ψ
∂x x=L
= h(T − T∞ )
k dT
x=L
−k
dx x= L
Et la condition initiale ψ ( x,0) = ψ 0 = T0 − T∞
73
On cherche une solution par séparation de variable du type:
ψ ( x, t ) = u (x ) v(t ) de sorte que

1
v u′′ = u vɺ
α
′′ ɺ
ou u = 1 v = −λ2 1° membre = f (x ) seul Cste par
f ( x ) = f (t ) =
u α v 1° membre = f (t ) seul identification
u ′′ + λ2u = 0  u = A sin (λx ) + B cos(λx )

vɺ + α λ2 v = 0  v = C e −α λ Vérifiant v(t=∞) = 0 ∀C
2
t
car T≡T∞ pour t = ∞ (équilibre thermique)
D’où la forme générale de la solution
ψ (x, t ) = e −α λ t [C1 sin (λx ) + C2 cos(λx )]

2
74
ψ (x, t ) = e −α λ t [C1 sin (λx ) + C2 cos(λx )]

2
∂ψ
= e −α λ t λ C1 = 0  C1 = 0
2
∂x x =0
∂ψ
= e −α λ t [− λ C2 sin (λ L )] = − ψ L = − e −α λ t C2 cos(λ L )
2 h h 2
∂x x= L k k
λk λL
Soit l’équation transcendante cotg (λ L ) = =
h Bi
régissant les valeurs propres
λn admissibles (sinon C2=0 et
la solution est triviale)
λL
cotg (λ L ) Bi
N .B. si Bi → ∞
π
λL λn L → (2n − 1)
λ 1L λ 2L λ 3L 2
(conditions de Dirichlet)
convection = T à la paroi en
x=L est imposée égale à T∞
0 π 2π 3π
75
Prof. EL OTMANI 25
ENSAJ TD
∞
ψ ( x, t ) =  Cn e −α λ t cos(λn x )
2
En final, on a: n
n =1
avec λn solution de (λn L ) tg (λn L ) = Bi (n=1,2,3,...∞)

C’est donc le mode spatial Cos(λ1x) qui perturbera le plus longtemps
(λ1 étant la plus petite racine de l’équation en valeur absolue)
Pour des conditions de Dirichlet λn = (2n − 1) π

2
Les constantes Cn sont obtenues à partir de la condition initiale:
∞
ψ ( x,0 ) = ψ 0 =  Cn cos (λn x ) soit, de manière équivalente:
n =1
L ∞ L
ψ 0 cos(λm x ) dx =  Cn  cos(λn x ) cos(λm x ) dx

0 n =1 0
sin(λ m L )  = 0 si m ≠ n entiers
ψ0  
λm  = 1 L + 1 sin ( 2 λm L )  si m = n
 2 2 λm

  
76
2ψ 0 sin (λm x )
 Cm =
Lλm + sin (λm x ) cos(λm x )
αt
Et en posant δm=λmL , F0 = = nombre de Fourier
L2
 x
sin  δ m 
ψ ∞
 L
= 2 e −αF δ
2
0 m
ψ0 n =1 δ m + sin (δ m ) cos(δ m )
avec δ m solution de δ m .tg (δ m ) = Bi
77
C13. Discrétisation en régime variable

(écoulement mono dimensionnel de chaleur)
On remplace la connaissance du champ de température instantané T(x,t)
par celle en (n+1) points de discrétisation spatiale à (n+1) stations
temporelles, de manière que:
 xi = i ∆x
 → T ( xi , t k ) = Ti k
 tk = k ∆t
78
Prof. EL OTMANI 26
ENSAJ TD
Ti-1k+1 Tik+1 Ti+1k+1

tk+1
tk ∆t
Ti-1k Tik Ti+1k
...
t2 ∆x
t1
x
x1 x2 ... xi-1 xi xi+1 ... xn=L
(i −1) ∆x i ∆x
(i +1) ∆x
79
Déterminons tout d’abord l’approximation finie de l’équation

transitoire de base
∂ ²T 1 ∂T
=
∂x ² α ∂t
La dérivée temporelle peut s ’évaluer à partir de
∂T
tk
T ( xi , t k + ∆t ) = T ( xi , t k ) + ∆t + ϑ ( ∆t ²) (Taylor )
∂t xi
∂T T (xi , tk + ∆t ) − T (xi , tk )
tk
soit = + ϑ (∆t )
∂t xi ∆t
∂T Ti k +1 − Ti k
k
ou = + ϑ (∆t ) Différence finie avant

∂t i ∆t (dans le temps)
80
La dérivée spatiale s’évalue de manière similaire par développement

en série de Taylor au voisinage du point xi, en tk:
∂T ∂ ²T ∂ ³T
k k k
1 1
Ti +k 1+1 = Ti k + ∆x + ∆x 2 + ∆x 3 + ϑ ( ∆x 4 )
∂x i 2 ∂x ² i 6 ∂x ³ i
k k k
∂T 1 ∂ ²T 1 ∂ ³T
Ti −k1+1 = Ti k − ∆x + ∆x 2 − ∆x 3 + ϑ ( ∆x 4 )
∂x i 2 ∂x ² i 6 ∂x ³ i
∂ ²T
k
2 4
Ti +k1+1 + Ti −k1+1 = 2 Ti k + ∆x + ϑ (∆x )
∂x ² i
∂ ²T Ti +k1 + Ti −k1 − 2 Ti k
k
+ ϑ ( ∆x 2 )
Différence finie centrée
=
∂x ² i ∆x 2 (sur l’espace)
81
Prof. EL OTMANI 27
ENSAJ TD
Par substitution de ces quotients finis dans l’équation de base, on a:
Ti +k1 + Ti −k1 − 2 Ti k 1 Ti k +1 − Ti k
=
∆x 2 α ∆t
Ou de manière explicite, la relation de récurrence:
α ∆t α ∆t
Ti k +1 = (Ti +k1 + Ti −k1 ) + (1 − 2 )Ti k
∆x 2 ∆x 2
a b
Elle montre que la température au point i à l’instant k s’obtient à partir

des températures en ce point et en ses voisins immédiats, à l’instant
précédent (valeurs héréditaires)
82
On montre que ce type de relation de récurrence ne peut s’utiliser

aveuglément, c-à-dire pour des choix arbitraires de ∆x et ∆t. En fait, la
réponse n’est bien conditionnée (au sens d’erreurs numériques qui
reste toujours bornées) que si
α ∆t 1
≤
∆x 2 2
Soit pour un coefficient d’influence de la température au point i en tk
sur la température en ce même point au temps tk+1 restant non négatif:
Ti k +1 = a (Ti +k 1 + Ti −k1 ) + b Ti k b≥0
Tik+1
tk+1
a a
b
tk
Ti-1k Tik Ti+1k+1
83
Méthode graphique de Schmidt

Un cas particulier intéressant consiste précisément à choisir le cas
limite b=0, autrement dit négliger l’influence du nœud sur lui-même à
l’instant ultérieur; on obtient ainsi
1 α ∆t 1
b=0 a= donc =
2 ∆x 2 2
1 k
soit Ti k +1 = (Ti +1 + Ti −k1 )
2
La température en un point i s’obtient par moyenne des températures
de ses voisins immédiats à l’instant précédent.
Tik+1
tk+1
1/2 1/2
tk
Ti-1k Tik Ti+1k+1
84
Prof. EL OTMANI 28
ENSAJ TD
Une telle relation de récurrence est susceptible d’une interprétation

géométrique simple: la température en xi à tk+1 est l’intersection de la
droite verticale x=xi avec la droite joignant les points représentatifs
Tki-1 et Tki+1
Tk(x)
Tik+1 Tki+1
connu
Tki-1
connu ∆x ∆x
xi-1 xi xi+1
85
Notons qu ’une fois le choix de la maille spatiale déterminée (∆x), les

instants auxquels la méthode de Schmidt donne les nouveaux résultats
sont complètement déterminés - ils sont séparé par
∆x 2
∆t =
2α
Invariablement, sauf si l’on change la maille spatiale pour certaines

parties de l’évolution transitoire.
86
Application
Soit une paroi d’épaisseur L initialement à 20°c. En t > 0, on porte

soudainement la température de sa face gauche à Tg=40°c et celle de
sa face droite à Td=30°c. Déterminer l’évolution du profil de T à
partir de la méthode de Schmidt.
L
On choisit ∆x = ,
3
les pas de temps où la réponse apparaîtra sont

∆x 2 L2
donnés par ∆t = = si α est la diffusivité thermique du milieu
2 α 18α
87
Prof. EL OTMANI 29
ENSAJ TD
Solution algébrique
t T1 = Tg T2 T3 T4 = Td
0 20 20 20 30
∆t 40 20 20 30
2∆t 40 30 25 30
3∆t 40 32,5 30 30
4∆t 40 35 31,25 30
5∆t 40 35,625 32,5 30
… … … … …
∞ 40 36,66 33,3 30
La solution finale étant la droite entre Tg=40°c et Td=30°c

Tg − Td
∆T par tranche ∆x = = 3,33°c
88
3
Solution graphique
40°c
30°c
T10 T20 T30 T40

20°c
L/3 2L/3 L
89
T
T11
40°c
30°c T41
T(∆t)
T10 T20= T21 T30= T31 T40

20°c
L/3 2L/3 L
90
Prof. EL OTMANI 30
ENSAJ TD
T
T11 = T12
40°c
30°c T41 = T42

T22
T32 T(∆t)
T10 T20= T21 T30= T31 T40

20°c
L/3 2L/3 L
91
T
T11 = T12
40°c
30°c T41 = T42

T22
T(2∆t)
T32 T(∆t)
T10 T20= T21 T30= T31 T40

20°c
L/3 2L/3 L
92
T
T11 = T12 = T13
40°c
T23
T33
30°c T41 = T42 = T43
T22
T(2∆t)
T32 T(∆t)
T10 T20= T21 T30= T31 T40

20°c
L/3 2L/3 L
93
Prof. EL OTMANI 31
ENSAJ TD
T
T11 = T12 = T13
40°c
T23
T(3∆t)
T33
30°c T41 = T42 = T43
T22
T(2∆t)
T32 T(∆t)
T10 T20= T21 T30= T31 T40

20°c
L/3 2L/3 L
94
T
T11 = T12 =….
40°c = T1∞
T2∞
T(∞∆t)
T3∞
T23
T(3∆t)
T33 T41 = T42 =….
30°c T22 = T4∞
T(2∆t)
T32 T(∆t)
T10 T20= T21 T30= T31 T40

20°c
L/3 2L/3 L
95
Remarque importante:
Les méthodes de récurrence explicite exposées ci-avant, et en
particulier la méthode de Schmidt, donnent naissance à un modèle
discret de la réalité qui a la propriété de présenter une capacité
thermique sous-estimée: le modèle est ainsi plus sensible aux
variations de température (il s’échauffe plus vite localement en cas
d’apport de chaleur et se refroidit plus vite dans le cas inverse) que le
milieu réel. Les températures obtenues par le modèle sont donc des
bornes supérieures aux températures réelles en phase continue
d’échauffement, et des bornes inférieures à celles-ci en phase continue
de refroidissement.
96
Prof. EL OTMANI 32
ENSAJ TD
C.14 Détermination d’un régime thermique

transitoire par la méthode de Schmidt
On désire analyser les distributions de température en régime

transitoire pour une paroi de roche granitique (k=1,6W/mK et
ρc=2,165.106J/m3K) initialement à T0=20°c uniforme.
Son épaisseur est de 9cm et en t=0, elle est soudainement exposée à un
flux qe=100W/m2 sur la surface gauche, tandis qu’un échange
convectif la maintient en contact avec un fluide restant à Tf=20°c, avec
h=10W/m2K.
97
T0 T1 T2 T3
L
qe=100W/m2 ∆x =
3
T4= Tf=20°c
h=10W/m2K
L=9cm
x=0 x=L
On prend ∆x=L/3=3.10-2m.
La méthode de Schmidt donnera en
conséquence des distributions de T à des instants successifs séparés de
∆x 2 9 10 −4
∆t = = = 608s
2α 2 7,4 10 −7
98
Les conditions aux limites donnent les relations à satisfaire à tout

instant à gauche et à droite de la paroi:
dT − qe 100
en x = 0 = =− = −62,5 K / m
dx k 1,6
dT TL − T f TL − 20
en x = L − = = K /m
dx k 0,16
h
Solution algébrique:
Les relations à satisfaire sont donc
T1 − T0
k k k
dT
= = −62,5  T0 = T1 + 62,5 ∆x = T1 + 1,875 ∀k > 0
k k k
a)
dx x =0 ∆x
b) Ti = Ti −1
k
( k −1
+ Ti +1
k −1
)2 ∀k > 0 i = 1,2
T − T3 T − T4
k
0,16 T2 + 0,6
k k k k k
dT
c) − = 2 = 3  T3 = ∀k > 0
k
dx x= L ∆x k 0,19
h
99
Prof. EL OTMANI 33
ENSAJ TD
t T0 T1 T2 T3 T4
0 20 20 20 20 20
∆t 21,875 20 20 20
2 ∆t 22,813 20,938 20 20
3∆t 23,281 21,406 20,469 20,395
4 ∆t 23,75 21,875 20,900 20,758
5∆t 24,2 22,325 21,316 21,109
… … … … … ⋮
∞ 35,625 33,75 31,875 30 20
La solution finale s’obtient par T∞ = ax + b
dT q
−k = qe  a = − e = −62,5
avec
dx 0 k T∞ = 35,625 − 62,5 x
= h[a L + b − 20]  b = 35,625
dT
−k
dx L
100
T
35
a ) T1 =
k (T 0
k −1
+ T2
k −1
)
2
Droite à pente
b) T0 = T1 + 1,875 const. passant
k k
pour k > 0
(T ) par T k
k −1 k −1
+ T3 1
c)T2 =
k
30 1
2
d ) T3 =
k (
0,16 T2 + 0,6 Droite
k
passant par
)
0,19 T k et T k 2 4
25
T01
T4k, k > 0
1,875
20
T00 T10, T11 T20, T21 , T22 T30, T31 , T32 Tf=20°c
L=9cm k/h=16cm
101
T
35
30
25
T4k, k > 0
20
T0 0 T10, T11 T2 0, T21 , T22 T30, T31 , T32
Tf=20°c
L=9cm k/h=16cm
102
Prof. EL OTMANI 34
ENSAJ TD
T
35
30
25
T 02
T 12 T4k, k > 0
20
T0 0 T10, T11 T2 0, T21 , T22 T30, T31 , T32
Tf=20°c
L=9cm k/h=16cm
103
T
35
30
25
T 03
T4k, k > 0
20
T00 T10, T11 T20, T21 , T22 T30, T31 , T32 Tf=20°c
L=9cm k/h=16cm
104
T
35
30
25
T 03
T 13
T23 T33 T4k, k > 0
20
T0 0 T10, T11 T2 0, T21 , T22 T30, T31 , T32
Tf=20°c
L=9cm k/h=16cm
105
Prof. EL OTMANI 35
ENSAJ TD
T
35
30
25
T4k, k > 0
20
T0 0 T1 0, T11 T2 0, T21 , T22 T30, T31 , T32
Tf=20°c
L=9cm k/h=16cm
106
T
35
30
25
T 14
T24 T34 T4k, k > 0
20
T00 T10, T11 T20, T21 , T22 T30, T31 , T32 Tf=20°c
L=9cm k/h=16cm
107
T
T 0∞ a ) T1 =
k (T 0
k −1
+ T2
k −1
)
35 2
T1 ∞
Droite à pente
b) T0 = T1 + 1,875 const. passant
k k
T2∞ pour k > 0

(T ) par T k
k −1 k −1
+ T3 1
c)T2 =
k 1
T3∞ 2
30
d ) T3 =
k (
0,16 T2 + 0,6 Droite
k
passant par
)
0,19 T k et T k 2 4
25
T 04
T 14
T24 T4k, k > 0
T34
20
T00 T10, T11 T20, T21 , T22 T30, T31 , T32 Tf=20°c
L=9cm k/h=16cm
108
Prof. EL OTMANI 36
ENSAJ TD
C.15 Rentrée d’un engin dans l’atmosphère

Donner, par la méthode de Schmidt, une évaluation du temps
pendant lequel la protection thermique d’un engin peut endurer les
conditions thermiques qui accompagnent la rentrée dans l’atmosphère,
sachant que la température maximale admissible pour le matériau
utilisé est 1100°c.
Cette protection thermique est assimilée à une plaque d’épaisseur
L=2,4cm, en acier spécial (k=42W/mK, c=420J/kgK, ρ=8000kg/m3),
initialement à la température uniforme de 40°c. L’écoulement de
chaleur est supposé unidimensionnel à travers l’épaisseur de la paroi.
Il n ’y a pas d’échange thermique coté intérieur. Coté extérieur, la
température au large est prise égale à 2200°c, tandis qu’on prend pour
h la valeur de 2800W/m2K. On utilise un pas de discrétisation
unidimensionnelle ∆x=L/4.
109
Tracer également le graphique des répartitions successives de

température dans la protection.
Comparer avec la réponse obtenue par l’utilisation des abaques en
conduction, en x=0 et x=L, et commenter.
N.B. L’échange thermique par rayonnement est supposé inclus
dans l’évaluation du coefficient d’échange.
110
L k 42
∆x = = 6 10 −3 m ; α = = = 1,25 10 −5 m 2 / s
4 ρ c 420 8000
→ ∆t =
(6 10 −3 )
2
= 1,44s
2 1,25 10 −5
Conditions aux limites:
dT
à droite: " = 0"
dx
T∞
 T4 = T3 n>0
n n
t>0
h à gauche: qconv = qcond
(
h T∞ − T0 =
n n
) k
∆x
(
T0 − T1
n n
)
T0 T1 T2 T3 T4
∆x  T0 = 628,571 + 0,71429 T1
n n
n>0
L
111
Prof. EL OTMANI 37
ENSAJ TD
Par l’utilisation des abaques:
Le problème = cas du mur dont une des faces est isolée et l ’autre
portée soudainement en contact avec un
fluide à T∞=2200°c
→ on utilise l’abaque p C42 (cas d’une paroi plane)
h l 2800 24 10 −3 −1
Bi = = = 1,6  Bi = 0,625
k 42
θ T − T∞ 1100 − 2200
= = = 0,5093
θ 0 T0 − T∞ 40 − 2200
La t° max. critique apparaît d’abord en surface → prendre abaque x = 1

L
→ On lit F0=0,19
t=
F0 L2
=
(
0,19 24 10 −3 )
2
= 8,755s
α 1
80000
112
A cet instant où la face externe atteint 1100°c:

t° face isolée:
x −1
abaque = 0 ; F0 = 0,19 ; Bi = 0,625
L
θ T − 2200
 = = 0,98
θ 0 40 − 2200
 T = 83,2°c
N.B. Méthode de Schmidt: modèle plus sensible aux variations de

température que le milieu réel: il s’échauffe plus vite localement (lors
d’un apport de chaleur)
113
C.16 conduction unidimentionnelle en régime variable

Analogie électrique
Position du problème - Discrétisation de la paroi
Puisqu’on considère un transfert de chaleur en régime variable, il
convient de tenir compte du temps et de la caractéristique inertie
(représenté par une capacité)
Soit la paroi représentée par le schéma équivalent à n cellules R-C-R:
Ts1 T1 Tj Tn Ts2
... ...
Qs1 R1 R2 Rj Rj+1 Rn Rn+1 Qs2
C1 Cj Cn
On convient de numéroter, dans un ordre croissant à partir de la face 1,

les n+1 résistances Rj, les n capacités Cj et les n t° internes Tj.
Les variables externes sont la t° et le débit de chaleur superficiels Tsf et
Qsf (f=1,2)
114
Prof. EL OTMANI 38
ENSAJ TD
 Ts1 − T1
 Qs1 =
Les débits de chaleur Qs1 et Qs2 sont donnés par:  R1
T −T
Qs 2 = n s 2
 Rn +1
En écrivant la conservation de débit aux différents nœuds de ce réseau,

on obtient le système:
 ɺ T1 − Ts1 + T1 − T2 = 0
 C1T1 +
 R1 R2
 ɺ T j − T j −1 T j − T j +1
C jT j + + =0 j = 2,..., n − 1
 Rj R j +1
 ɺ Tn − Ts 2 Tn − Tn −1
 CnTn + R +
Rn
=0
 n +1
115
Application
Reprenons le problème de la paroi de roche granitique (cfr. p C51).
Représentons cette paroi par un schéma à 2 cellules R-C-R:
L L L 1
T0 3k T1 3k T2 3k T3 h T4=20°c
qe
ρcL ρcL On travaille par unité de surface. La
capacité totale est ρcV. Donc, par unité
2 2
de surface ρcL. Et on partage en 2 →
ρcL/2 à chaque nœud.
On en déduit le système d’équations:

 ρ c L dT1 3k
 2 dt + L (T1 − T0 ) + L (T1 − T2 ) = 0 (1)
3k
 T0 − T1
 ρ c L dT2 + 3k (T − T ) + T2 − T4 = 0 et qe =
 2 1
L 1
(2) L
2 dt L + 3k

 3k h
116
L 0,09
d ' où T0 = qe + T1 = 100 + T1  T0 = 1,875 + T1
3k 3 1,6
ρ c L2 dT1
(1)  − 1,875 + T1 = T2
6k dt
ρ c L  ρ c L2 d 2T1
dT  3k  ρ c L2 dT1 
 + 1  +  − 1,875 
2  6k dt 2 dt  L  6k dt 
dans (2 )  1  ρ c L2 dT1 
+ 
 − 1,875 − T4  = 0
L 1  6k dt 
+
3k h
 ρ c L2 
(ρ c )2 L3 d 2T1 + dT1  ρ c L + ρ c L + 6k 

12k dt 2
dt  2 2 L 1
+
→  3k h 

+
T1 3k
= 1,875 +
L 1 L
1
L 1
(1,875 + 20 )
+ +
3k h 3k h
117
Prof. EL OTMANI 39
ENSAJ TD
ou encore :
 
 2  12k
 
ρcL  dT1 + (ρ c ) L T = 284,21 12k
2 3
d 2T1  12k
+ +
dt 2  ρ c L2 1  dt L 1 1
+ (ρ c )2 L3
 L 1  3k h
 + 
 3k h 
d 2T1 dT
+ 1,18 10 −3 1 + 4,73 10 −8 T1 = 1,597 10 − 6
dt 2 dt
−5 −3 1,597 10 − 6
 T1 = C1 e − 4,15 10 t + C2 e −1,1398 10 t +
4,73 10 −8
33, 75
 −5 ρ c L  −3  ρ c L2 
2
− 4 ,15 10−5 t
T = 31,875 + C1 e 1 − 4,15 10  + C2 e −1,1398 10 t 1 − 1,1398 10 −3 
 2  6k   6k 
−5 −3
= 31,875 + 0,9242 C1 e − 4,15 10 t − 1,082 C2 e −1,1398 10 t
118
C1 , C2 ?
en t = 0 : T1 = T2 = 20°c
 C1 + C2 + 33,75 = 20
 
31,875 + 0,9242C1 − 1,082C2 = 20
 C = −13,33
  1 −5 −3
C2 = −0,4148  T1 = −13,33 e −4,15 10 t − 0,4148 e −1,1398 10 t + 33,75

−5 −3
 T2 = 31,875 − 12,3 e − 4,15 10 t + 0,449 e −1,1398 10 t
 T0 = T1 + 1,875
T1 = 20,96°c T1 = 33,75°c
après 20 min (1200 s ) : T2 = 20,29°c pour t → ∞ : T2 = 35,625°c
T0 = 22,84°c T0 = 31,875°c
119
120
Prof. EL OTMANI 40
ENSAJ TD
II. Convection
121
Convection: qu’est qu’il faut savoir?

• Concerne principalement les fluides
• Mode de transmission: transport de particules
• naturel: ≠ de densité généré par ≠ de T
• forcé: action mécanique (≠ de pression d’altitude, vent,…)
• Lois fondamentales
Générales 1D
Q = Qin − Qout = −  q.next dA Q=q A

q = − h ∆T e q = − h (T f − Ti )
122
Convection forcée:
Cas de la plaque
– Nombre de Reynolds
UL
• U = vitesse du fluide au large (m/s)
ReL =
• L = longueur de la plaque (m)
• ν = viscosité cinématique du fluide (m²/s)
ν
Transition entre régimes laminaire et turbulent: Re = 3,5 10 5
cr
L
123
Prof. EL OTMANI 41
ENSAJ TD
Convection forcée:
Cas de la plaque
UL
ReL =
ν
cr
L
– Cas de la paroi isotherme (TP ≠ T∞) 0,6≤ Pr ≤50
• Laminaire: Nu L =
hL
= 0,664 ReL Pr 1/ 3
1/ 2
kf
• Turbulent: Nu L =
hL
= 0,036 ReL Pr 1/ 3
0,8
kf
124
Convection forcée:
Cas de la plaque
UL
ReL =
ν
cr
L
– Cas de la paroi isotherme (TP ≠ T∞) 0,6≤ Pr ≤50
• Laminaire: Nu L =
hL
= 0,664 ReL Pr 1/ 3
1/ 2
kf
• Turbulent: Nu L =
hL
= 0,036 ReL Pr 1/ 3
0,8
kf
– Cas de la paroi à flux imposé qp (W/m²)
• Laminaire: qp L
Nu L = = 0,6795 ReL
1/ 2
Pr1/ 3
(T p − T∞ ) k f
Nu = Nombre de NUSSELT Pr = Nombre de PRANDTL = ν
(coefficient d’échange réduit) α
125
Convection forcée:
Cas de la conduite
• Uav = vit. moy. d’une sect. dro. d’écoulement
U av D
• D = diamètre du tube ou ReD =
Dh diamètre hydraulique équivalent ν
• ν = viscosité cinématique du fluide
N.B.:
Dh = 4 A avec P
A= section de passage
P = périmètre « mouillé »
Transition entre régimes laminaire et turbulent: Re Dcr = 2100 à 2300
126
Prof. EL OTMANI 42
ENSAJ TD
Convection forcée:
Cas de la conduite
U av D
N.B.: avec A= section de passage
Dh = 4 A P = périmètre « mouillé »
P
– Cas de la paroi isotherme (TP ≠ T∞)
• Laminaire: hD
Nu L = = 3,658
kf
= 0,023 ReD Pr 1/ 3 si lisse
0 ,8
• Turbulent: Nu L
f = f(ε)
= f ReD Pr1/ 3 si rugueux ε = rugosité de paroi
8
127
Convection forcée:
Cas de la conduite
U av D
N.B.: avec A= section de passage
Dh = 4 A P = périmètre « mouillé »
P
– Cas de la paroi isotherme (TP ≠ T∞)
• Laminaire: hD
Nu L = = 3,658
kf
= 0,023 ReD Pr 1/ 3 si lisse
0 ,8
• Turbulent: Nu L
f = f(ε)
= f ReD Pr1/ 3 si rugueux ε = rugosité de paroi
8
– Cas de la paroi à flux imposé qp (W/m²)
qp D
• Laminaire: Nu L = = 4,364
(T p − Ts ) k f
128
Convection libre
– Nombre de GRASHOF
• g = acc. gravité
• β = coef. d’expension therm. g β ∆T L3
• ∆T = écart t° nominal Tp-Tf GrL =
• L = longueur caractéristique ν2
• ν = viscosité cinématique
1
Remarque: pour des gaz parfaits, β = ; Tm = temp. moy.
Tm
Sinon, tables dans le bouquin d’exercices!
129
Prof. EL OTMANI 43
ENSAJ TD
Convection libre
– Nombre de RALEIGH Ra L = GrL Pr
130
Convection libre
– Plaque verticale isotherme
0 , 25
 Pr 
Nu L = 0,678 RaL0, 25  
 0,952 + Pr 
131
Convection libre
– Plaque verticale isotherme
0 , 25
 Pr 
Nu L = 0,678 RaL0, 25  
 0,952 + Pr 
– Corrélations empiriques
Plaque Cylindre Sphere
= c1 RaL laminaire 10 4 < RaL < 10 9
0 , 25
Verticale Horizontale
Nu L Ins table Stable
= c 2 RaL RaL > 109
0 , 333
turbulent C1 0,59 0,54 0,27 0,47 0,49
C2 0,13 0,14 - 0,1 -
132
Prof. EL OTMANI 44
ENSAJ TD
H1. Convection forcée

La paroi d’un wagon frigorifique est constituée d’une tôle
d’aluminium tapissée sur sa face interne d’une couche
d’isolant.
La température du contenu du wagon est supposée uniforme
et maintenue à Tint=0°c.
Le wagon se déplace dans un environnement à Text=27°c
uniforme et à pression atmosphérique.
Le gain de chaleur en régime permanent est mesuré égal à
qe=27,5W/m².
On demande d’évaluer le coefficient d’échange convectif
externe, les températures de paroi du wagon et le
coefficient d’échange interne, en supposant l’écoulement
externe uniforme le long du wagon (L=12m).
133
Aluminium
e = 3 mm
k = 238 W/mK
27,5 W/m² = qe
U=40km/h Paroi
Text=27°c
Isolant
L=12m e = 5 cm
k = 0,06 W/mK
Tint=0°c
134
Text Text Text

1
Rext = Rext
hext
Talu-ext Tp-ext
e
Ralu = alu
kalu
Ralu + Riso Rext + Ralu
Talu-in + Riso + Rint
eiso
Riso =
kiso
Tiso Tp-in
1
Rint = Rint
hint
Tin Tin Tin
135
Prof. EL OTMANI 45
ENSAJ TD
hext
A l’extérieur, on est en présence d’un écoulement forcé le
long des parois planes (plaques)
U = 40km / h = 11,11m / s UL
Re L = = 8,51106 > 3,5 105
L = 12m ν air
300 K
ν air
300 K
= 1,566 10 −5 m² / s  Turbulent
Ne connaissant pas la temp. de la plaque,
on prend les prop. de l’air à 300K
L’écoulement étant turbulent, on a:

1  ν
Nu L = 0,036 Pr 3
Re 0,8 = 1,1253 10 4 avec Pr =  = 0,711
 α
Le coefficient d’échange convectif externe se déduit de
27° c
NuL kair 1,1253 0,026
hext = = = 24,5W / m ² K
L 12
136
Le gain de chaleur des parois s’identifie au chauffage convectif, soit
(
qe = qconv = hext Text − Tp
ext
) , dont on tire
q 27,5
= Text − e = 27 − ≅ 25,9°c
ext
Tp
hext 24,5
A travers la paroi composite, ce flux est véhiculé par conduction

Tp − Tp 25,9 − Tp
ext int int
∆T
qe = qcond = = = = 27,5 W / m²  T p = 3°c
int
 rés. ealu + eiso 0,003 + 0,05

k alu kiso 238 0,06
qe 27,5
On en déduit finalement hint = = = 9,2 W / m ² K
Tp − Tint
int
3
Convection interne libre (s’il n’y a pas de ventilation)
La quant. de chal. que la mach. frigo. Int. devra évacuer s’obtient par
Qout = Qin = qe S Où S est la surf. lat. du wagon exp. à la conv.
137
H2. Convection libre

Une plaque d’aluminium à chaud (T0=80°c) est
introduite dans un local rempli d’air calme à 20°c.
La plaque est suspendue verticalement (hauteur
0,5m, largeur 1m, épaisseur 5mm) et on demande
quelle sera sa température après 1 minute?
138
Prof. EL OTMANI 46
ENSAJ TD
On évalue les propriétés de l’air à la température moyenne:

Tp + T∞
0
80°c + 20°c
Tm = = = 50°c = 323K
2 2
On réalise une interpolation dans les tables en considérant
pour cette lecture que 323K = 325K
k300 + k350 29,7 + 26,1 −3
km = = 10 = 27,9 10 −3W / mK
2 2
γ +γ 1,566 + 2,062 −5
γ m = 300 350 = 10 = 1,814 10 −5 m² / s
2 2
α +α 2,203 + 2,920 −5
α m = 300 350 = 10 = 2,562 10 −5 m² / s
2 2
γ m Pr300 + Pr350
Pr = = = 0,708  0,6 < 0,708 < 50  paroi isotherme
αm 2
139
On calcule le nombre de GRASHOF:

g L3 9,81 0,53
GrL = β ∆T =
1
(80 − 20) = 0,692 109
γ2 (1,814 10 )
−5 2 323
1
β= pour les gaz parfaits
Tm
Ra L = GrL Pr = 0,49 109 < 109  laminaire
0, 25
On tire NuL = 0,678 Ra L 0, 25  Pr 

 = 81,52
 0,952 + Pr 
Le coefficient d ’échange en convection naturelle vaut dès lors

Nu L k air 81,52 0,0279
h= = = 4,55 W / m ² K
L 0,5
La perte de chaleur initiale de la plaque s’identifie à
Qconv = h(T p − T∞ ) S = 4,55 60 2 0,5 1 = 273W

140
On est dans les conditions pour réaliser une analyse transitoire

moyenne où la plaque est considérée isotherme à tout instant car
l
h −3
Bi = 2 = 4,55 2,5 10 = 4,7 10 −5 << 0,1 → cond >> conv
k alu 240
Ce qui se comprend intuitivement par la bonne conductibilité de

l’aluminium et le faible échange convectif. Le temps caractéristique
de refroidissement est d’après la page C25:
ρ cV 2707 905 2,5 10 −3
τ= = = 1346 s très long
hS 4,55
Après 1 minute, on aura d ’après la même page
T1min − T0
θ = θ 0 e − t /τ = e − 60 /1346 T1min = 20 + 60 0,956 = 77,4°c
T∞ − T0
141
Prof. EL OTMANI 47
ENSAJ TD
H3. Convection en conduite -

écoulement forcé
Le combustible d’un réacteur nucléaire est
constitué de barreaux d’uranium métallique
de diamètre Dint=3cm. Leur refroidissement
uav=20m/s est assuré par circulation de CO2 gazeux à
Tf=327°c et p=10atm, à la vitesse moyenne
de 20m/s, dans une gaine cylindrique
concentrique de diamètre Dext=6cm.
Dint=3cm
Dext=6cm La température de parois des barreaux en
régime permanent est Tp=560°c.
Tp=560°c
On demande de vérifier que l’évacuation de
chaleur est suffisante pour empêcher le
Tf=327°c dépassement de la température de chgt.
p=10 atm d’état de l’uranium (Tmax=660°c)
142
Le diamètre hydraulique de la conduite de section annulaire s’écrit

4 A A = section de passage du fluide
Dh =
P P = périmètre « mouillé » (surface latérale en contact avec le fluide)
π (D 2
− Dint
2
)
Dh = 4 4
ext
= Dext + Dint = 3 10 − 2 m
π (Dext + Dint )
Pour un gaz à pression non atmosphérique, les caractéristiques de
densité se déduisent de celles à pression atmosphérique par la relation
des gaz parfaits (en supposant que celle-ci est d’application)
p
ρ
= RT R=
ℜ 0 Cste des gaz parfaits
M
=
poids moléculair e
JK −1kg −1 [ ]
Pour deux états (p1, ρ1, T1) et (p2, ρ2, T2), on a donc
ρ CO 300 K / 1atm
≅ 1,8kg / m³
ρ 2 p2 T1 Sachant que 2
=
ρ1 p1 T2 ρ CO 600 K / 10 atm
= 1,8
10 300
= 8,99kg / m ³
On aura 2
1 600
143
On tire de la table A.4. pT. 10:

µ CO 600K = 26,83 10−6 kg m −1s −1
2
µ 26,83 −6
γ CO 600 K
= = 10 = 2,99 10 −6 m² / s
2
ρ 8,99
D’autre part, la table A.4. Donne, si on néglige les effets de pression,
= 0,043W / mK = 1,076 103 J / kgK
600 K 600 K
kCO2 CpCO2
k 0,043
D’où α = = 10 −3 = 4,45 10 −6 m ² / s
ρ Cp 8,99 1,076
Et Pr = ν = 2,99 = 0,673
α 4,45
144
Prof. EL OTMANI 48
ENSAJ TD
L’écoulement est défini par:

uav D 20 0,03 6
Re D = = 10 = 0,2 106
ν 2,99
Et est tout à fait turbulent. En considérant les parois mouillées

comme parfaitement lisses, on a
Nu D = 0,023 Pr1/ 3 Re D = 0,023 (0,876 ) 1,741 10 4 = 350,8

0,8
Le coefficient d’échange en convection forcée est donc

NuD kCO2 350,8 0,043
h= = = 502,8 W / m² K
D 0,03
Le débit de chaleur emporté par convection à la périphérie du
barreau d’uranium, par mètre de hauteur, est donc
qconv = h ∆T S = 502,8 (560 − 327 ) π Dint = 11,041 kW / m
145
On considère que cette chaleur à évacuer correspond à l’irradiation

neuronique uniformément répartie dans le barreau comme une source
de chaleur par unité de volume Qv donnée par:
Où V est le volume du barreau correspondant à
Qgen = QvV = Qconv
1 mètre de hauteur
Qconv 11,041 3
Qv = = 2
10 = 15,62 106 W / m³
V D
π int
4
La distribution de température dans le barreau en régime permanent
correspond donc à la solution de l’équation de la chaleur en
coordonnées cylindriques avec une distribution de sources volumiques
Qv, soit: Q d 2T 1 dT Qv
∆T + v = 0 où 2
+ + =0
kur dr r dr kur
La solution homogène est: Th = a ln r + b et

La solution particulière est: Tpart = − Qv r 2
4kur
146
Les conditions aux limites donnent:

a=0
T (r = 0 ) finie = −a ∞ + b 
b fini
 D  Qv
T  r = int  = Tp  b− rint = T p
2
 2  4kur
Qv   r 2 
T = Tp + rint 1 −   
2
Dès lors:
  int  
4kur r
La température maximale est atteinte au centre du barreau
2
15,62 10 6 0,0152
Tmax = T (r = 0 ) = T p +
Qv rint
= 560 + = 581,4°c
4 kur 4 41
Admissible!
147
Prof. EL OTMANI 49
ENSAJ TD
H4. Échanges convectifs mixtes

Une chambre dont l’ambiance est à Tint=18°c est séparé d’un vide non
ventilé sous-toiture par un plafond de plâtre (kp=0,5W/mK, ep=2cm).
Les échanges convectifs de part et d’autre du plafond (cambre et grenier
non ventilés) sont régis par un coefficient d’échange en convection
naturelle hint=8,2W/m²K, et la surface du plafond est L=10m x l=6m.
La toiture présente une pente de 30° et le plâtre du toit est parallèle au
grand coté. La couverture est constituée de tuiles (kt=0,7W/mK,
et=3cm) et à l’extérieur règne une température uniforme hivernale
Text=0°c avec un vent continu soufflant à 30km/h parallèlement au faîte
du toit.
On demande les pertes calorifiques de la chambre à travers le système
plafond-grenier-toiture en régime permanent, et l’influence sur celles-ci
d’un tapis de laine de verre (kl=0,03W/mK, el=0,06m) déposé sur le sol
du grenier
148
Text=0°c
U=30km/h
Tuiles
kt=0,7W/mK
hint
et=3cm
30°
Plâtre
Chambre
kp=0,5W/mK
hint=8,2W/m²K
ep=2cm
Tint=18°c
L=10m
l=6m
149
hext?
Pour la convection extérieure (convection forcée - paroi plane):
300
10
UL 36
Re L = = = 6,2 106 > 3,5 105  turbulent
ν air
0° c
1,123 + 1,566 −5
10
2
On en déduit
1 0° c  ν
NuL = 0,036 Pr 3
Re 0,8 = 8755 avec Prair  =  = 0,718
 α
Et
Nu L k air 8755 0,0242
hext = = = 21,187 W / m² K
L 10
150
Prof. EL OTMANI 50
ENSAJ TD
Modèle thermique de la toiture 1

Rext =
hext
Text = 0°c Text = 0°c
et
Rtoit =
kt
hext = 21,187W/m²K
Tv inconnu
kt=0,7W/mK 1
1 Rt =
hint = 8,2W/m²K Rint = ht
et=3cm hint
Tv inconnu
Calcul du coefficient de transmission résultant (par unité de surface):

1 1 et 1 1 0,03 1
Rt = = + + = + + = 0,212 m ² K / W  ht = 4,72W / m ² K
ht hext kt hint 21,187 0,7 8,2
l
St = 2 L 2 = 69,282m²  Q = h (T − T ) S = 326,8 T (W )
cos 30°
t t v ext t v
151
Modèle thermique du plafond sans isolant Tv inconnu

1
Rint =
hint
Tv inconnu
kp=0,5W/mK
hint = 8,2W/m²K 1
ep Rp =
ep=0,02m R plafond = hp
kp
Tint = 18°c
hint = 8,2W/m²K 1
Rint =
hint
Tint = 18°c
1 ep 2 0,02 2
Rp = = + = + = 0,284 m ² K / W  h p = 3,52W / m ² K
h p k p hint 0,5 8,2
S p = L l = 60m ²  Q p = h p (Tint − Tv ) S p = 3801,6 − 211,2 Tv (W )
152
Équilibre thermique sans isolant
Text = 0°c
1
Qt Rt =
ht
Tv inconnu
1
Qp Rp =
hp
Tint = 18°c
Q0 = Qt = Q p  326,8 Tv = 3801,6 − 211,2 Tv
Tv = 7,066°c  Q0 = 326,8 7,066 = 2,309kW
153
Prof. EL OTMANI 51
ENSAJ TD
Modèle thermique du plafond avec isolant Tv inconnu

1
Rint =
Tv inconnu hint
hint = 8,2W/m²K kiso=0,5W/mK eiso
Riso = 1
eiso=0,02m kiso Rp =
hp
ep
R plafond =
kp
Tint = 18°c kp=0,5W/mK
hint = 8,2W/m²K ep=0,02m 1
Rint =
hint
Tint = 18°c
1 e p eiso 2 0,02 0,06 2 ′

Rp = = + + = + + = 2,284 m ² K / W  h p = 0,438W / m ² K
h p k p kiso hint 0,5 0,03 8,2
′ ′ ′ ′
S p = L l = 60m²  Q p = h p  Tint − Tv  S p = 472,87 − 26,27 Tv (W )
 
154
Équilibre thermique avec isolant
Text = 0°c
′ 1
Qt Rt =
ht
Tv ’inconnu
′ ′ 1
Qp Rp =
′
hp
Tint = 18°c
′ ′ ′ ′ ′
Q0 = Qt = Q p  326,8 Tv = 472,87 − 26,27 Tv
′ ′
Tv = 1,339°c  Q0 = 326,8 1,339 = 0,438kW
′
Q0 0,438
= = 0,19 (réduction de plus de 5 fois)
Q0 2,309
155
156
Prof. EL OTMANI 52
ENSAJ TD
III. Rayonnement
157
Rayonnement: qu’est qu’il faut savoir?

• Se produit avec ou sans support physique
(plus c’est vide, mieux ça marche!)
• Mode de transmission:
rayonnement électromagnétique (lumière)
• Intensité fonction de
• la température (absolue) en Kelvin
• des propriétés de la surface
• Lois fondamentales (si visibilité totale entre les surf.)
Générales 1D
Q=q A
Q = Qin − Qout = −  q.next dA
*
(
q = −σ ε if T f − Ti
4 4
)
(
q = −σ ε12 T2 − T1 e
* 4 4
) σ = 5,67 10 W / m ² K
−8 4
158
R1. Influence des échanges radiatifs

Faire l’étude d’une paroi de conductibilité k=6W/mK et d’épaisseur
l=0,3m, dans les conditions suivantes:
Te k Tp
σ ε (TR4) σ = 5,67 10 −8W / m² K 4

ε = 0,8
TR
σ ε (Tp4) Te = 600 K
TR = 1200 K
ε
l
159
Prof. EL OTMANI 53
ENSAJ TD
qcond = qray
k 4
(
− (Tp − Te ) = σ ε Tp − TR
4
) Visibilité totale (FpR=1)
( ) (T
l
= h* (Tp + TR )(Tp − TR ) h* = σ ε T p + TR + TR )
2 2
p
Eq. non linéaire pour Tp

k k
σ ε T p 4 + Tp + C = 0 C = −σ ε TR − Te
4
l l
( )
f Tp
Solution par Newton - Raphson (N-R)

Essai T p 0  ( )
f Tp ≠ 0
0
On cherche une correction ∆T p telle que
( ) ( )
f T p + ∆T p ≅ f T p + ∆T p
0 0 ∂f
=0  ∆T p =
− f Tp ( ) 0
∂Tp 0 ∂f
Tp
∂Tp 0
Tp
160
( )= σ εT
k k
+ Te − σ ε T p + T p
4 4
− f Tp
0
R
 ∆T p = l l
∂f k
4σ ε Tp +
3
∂T p l
0
Tp
  k
Si le milieu est mauvais conducteur  σ ε T >> T 
4
 l 
Tp → TR
k 4
Si le milieu est peu isolant  T >> σ ε T 
l 
Tp → Te
4
T  T 
106058,5 − 2000  p  − 4,536 p 
 100   100 
Pour les constantes choisies ∆T p = 3
 Tp 
20 + 0,18144  
 100 
N.B. Attention au « domaine d’attraction » de la méthode N-R

161
Solution itérative
(0 ) Te + TR
i = 0 Tp = = 900 K ∆Tp = 382,86 K!
2
= 1282,86 K (> TR !) ∆Tp = −105,32 K
(1)
i =1 Tp
(2 )
i = 2 T p = 1177,53K ∆Tp = −14,87 K
(3 )
i = 3 T p = 1162,67 K ∆Tp = −0,27 K
(4 )
i = 4 Tp = 1162,4 K
Rem: plus systématiquement, Tp

n +1
= Tp −
n ( )
f Tp
n
(N-R)
f ′(T )
n
p
qcond =
k
(Te − Tp ) = −20 (562,4) = −11,2479 kW / m²
l
 T p  4  T  4 
(
qray = σ ε T p − TR = −4,536 
4 4
)  −  R   = −11,2479 kW / m²
100   100  
 
(le milieu est mauvais cond. Donc Tp → TR) Vérif: qcond = qray
162
Prof. EL OTMANI 54
ENSAJ TD
Variantes:
Faire la même étude
a) avec un matériau cent fois meilleur conducteur, càd k=600W/mK
b) à basse température: TR=600K; Te=0K
k=600W/mK
a b
Te + TR
(0 ) Te + TR
(0 )
Tp = = 900 K ∆T p = −251,236 K Tp = = 300 K ∆T p = −296,518 K
2 2
(1) (1)
Tp = 649 K ∆T p = −5,6138 K Tp = 3,4819 K ∆T p = −0,5425K
(2 ) (2 )
T p = 643K ∆T p = −1,75 10 − 3 K T p = 2,939 K ∆T p = −4,349 10 −10 K
(3 ) (3 )
Tp = 643,14 K Tp = 2,939 K
qcond = q ray = −86,28 kW / m ² qcond = qray = −5,878 kW / m ²
le milieu est bon cond. le milieu est bon cond.

Donc Tp → Te Donc Tp → Te
163
R2. Les écrans radiatifs

Préliminaire
Flux échangé par 2 surfaces grises, opaques,
isothermes, d’étendue infinie (visibilité totale)
T1 T2
( 
)
 T  
4 4
*  T 
q12 = σ ε12 T1 − T2 = 5,67 ε12  1  −  2   W / m²
* 4 4
 100   100  

ε1ε 2 1 ε1 ε2
ε12* = = Coefficient d’émission
ε1 + ε 2 − ε1ε 2 1 + 1 − 1
résultant
ε1 ε2
ε
N.B. si ε1 = ε 2 = ε , → ε12* =
2 −ε
ε 1 >> ε 2 , → ε 12 ≅ ε 2 (le plus faible!)
*
164
Influence d’un écran (gris, opaque, isotherme, d ’étendue infinie)

q1e = σ ε1e T1 − Te
*
( 4 4
) (
qe 2 = σ ε e 2 Te − T2
* 4 4
)
T1 T2
T + T2
4 4
q1e = qe 2  Te = 1 si ε1 = ε 2 = ε ∀ε e
4
2
(système « symétrique » thermiquement)
au niveau des propriétés ε εe ε
T1 − T2
4 4
qeff = q1e = qe 2 = σ ε1e
*
2
1 ε1e ε εe 2−ε
*
1
= q12 = q12 q12=q sans écran
2 ε12* 2 ε + εe − ε εe ε
1 2 −ε 2−ε
N.B. si ε e >> ε → qeff ≅ ε q12 ≅ q12 ≅ q12 ( pas efficace)
2 ε 2
1
εe = ε → qeff = q12 (réduction de moitié )
2
1 2 −ε
ε e << ε → qeff ≅ εe q12 << q12 (très efficace)
2 ε
Ex: plans acier oxydé courant ε=0,8; écran en cuivre poli εe=0,02; qeff=0,015 q12
165
Prof. EL OTMANI 55
ENSAJ TD
R3. Modèle de la bouteille thermos

On idéalise la géométrie d’une bouteille thermos par une double
enveloppe cylindrique de verre rendu réflecteur au sein de laquelle règne
le vide (d entre paroi = 0,5 10-3 m).
On demande le taux de refroidissement du volume d’eau chaude a 60°c
contenu dans le thermos, lui-même placé dans une ambiance à 20°c, les
pertes de chaleur se faisant par rayonnement. Le couvercle et la base du
thermos sont supposés parfaitement adiabatiques.
Tint=60°c
=333K
φ = 10cm = 2 R
h = 25cm
S = 2 π R h = 0,07854 m²
Text=20°c
V = π R 2 h = 1,9634 10 −3 m³ ou litre
h =293K
166
Déperdition calorifique latérale par rayonnement si Tint=333K

homogène et Text=293K homogène:
*
(
Qray = Sε12 σ Tint − Text
4 4
)
On contrôle que le facteur de vue pour 2 cylindres concentriques de
rayons très voisins et de longueur non négligeable est voisin de 1.
L ’abaque du cours théorique donne pour
h 25
γ=
rext 5,00
= = 1,01 et β= = = 5,05
rint 4,95 rint 4,95
 Fext −int = 0,99
Dès lors ε12* = ε = 0,087 Si ε=0,16 pour le verre réflecteur

2 −ε
(
Qray = 0,07854 (0,087 ) 5,67 3,344 − 2,934 = 1,908W )
167
Ce débit de chaleur quittant le système constitué par le thermos, on a
Qray = − ρ CV Tɺ V Ou Tɺ est le taux moyen de refroidissement

(-) car perdu par le fluide du contenu supposé homogénéisé à tout
− Qray instant (pas d’effet de conduction ou de
Tɺ = convection)
ρ C pV
Pour l’eau à 320K: ρ = 989,3 kg / m ³

Table p T8
C p = 4177 J / kgK
− 1,908
Tɺ = = −2,35 10 − 4 K / s = −0,85 K / h
8113
168
Prof. EL OTMANI 56
ENSAJ TD
R4. Dimensionnement
d’un panneau radiatif
On désire déterminer la surface d’un panneau radiatif, considéré
comme un corps noir, présentant une température superficielle en
régime Tp=80°c, et destiné à compenser par son rayonnement
l’inconfort de l’ambiance d’une chambre à Tch=16°c, en étant
suspendu au plafond de la chambre, à une distance h=2m de la tête
d’un sujet.
Le coefficient de convection du sujet avec l’ambiance de la
chambre est pris égal à hch=4W/m²K et la notion de confort s’établit
pour Tt=32°c au minimum.
(Tt désignant la température superficielle de la tête du sujet,
considérée comme isolée, en équilibre thermique avec les seuls flux
extérieurs)
169
Ap Tp = 80°c
Air calme h=2m

Tch=16°c
hch=4W/m²K
Ttmin = 32°c
At
170
a) débit de chaleur apporté à At

(panneau et chambre = émetteurs parfaits)
Fp-t facteur de vue panneau-tête =

( ) (
Qin = σ Tp Ap Fp −t + σ Tch Ach Fch −t
4 4
) fraction énergie émise par le panneau
qui atteint la tête
Fch-t idem pour chambre-tête
b) débit de chaleur perdu par At

(considéré aussi comme émetteur parfait)
( )
Qout = σ Tt At + hch (Tt − Tch ) At
4
Rayonnement + convection
Pas de facteur de vue car c’est

la quantité totale perdue
171
Prof. EL OTMANI 57
ENSAJ TD
En vertu des relations de réciprocité pour deux surfaces échangeant

de la chaleur par rayonnement:
Ai Fij = A j F ji et F j
ij =1 ∀i , on a
Ap Fp −t = At Ft − p , Ach Fch −t = At Ft −ch

Le système tête + chambre + plafond
Ft −t + Ft − p + Ft − ch = 1 étant fermé géométriquement.
La tête du sujet est considérée comme convexe ou plane → Ftt=0
Par équilibre thermique de la tête, on a ainsi

Qin = Qout
et en vertu des relations ci-dessus, on tire
Qin = σ T p At Ft − p + σ Tch At (1 − Ft − p ) = Qout = σ Tp At + h(Tt − Tch ) At

4 4 4
172
σ Tt 4 − σ Tch 4 + h(Tt − Tch )

 Ft − p =
(
σ Tp 4 − Tch 4 )
Tt = 32 + 273 = 305 K 
 490,7 − 395,5 + 4 16
Tp = 80 + 273 = 353K   Ft − p = = 0,328
484,9
Tch = 16 + 273 = 289 K 
On doit donc dimensionner le panneau de manière telle que le

facteur de vue tête-panneau atteigne au moins la valeur
prescrite.
173
Si l’on se sert des tables du cours théorique, on a par exemple:

a) surface élémentaire dAt au centre d’un panneau rectangulaire a, b à
une distance c=2m du plan du panneau, pour lesquels une évaluation
grossière donne: (cf. abaque p. R18)
a
0,328 x = ≅ 1,4
FdA − A = = 0,082 si c On résonne sur la tête
t p
4 b voyant un quart de panneau
y = ≅ 0,4
c
a = 1, 4 2 = 2,8 
soit  A = 4 a b = 8,96 m ²
b = 0, 4 2 = 0,8
p
b) Disque At de rayon r1=20cm face à un disque Ap de rayon r2 (à

définir) à une distance L=2m (centre à centre). L’abaque p. R18
donne (grossièrement):
L 2 r
γ= = = 10 et F12 = 0,328 si β = 2 ≅ 0,75
r1 0,2 L
 r2 = 1,5m et Ap = π r2 = 7,07 m²
2
174
Prof. EL OTMANI 58
ENSAJ TD
On peut aussi réaliser un calcul exact du facteur de vue entre une

facette dAt au centre du disque Ap à une distance h de celui-ci. La
forme théorique donne pour 1/4 de disque:
Rp
r
φ dAp
1 cos θ t cos θ p 1 2 Rp
h2
 r dr dφ
π d
θp= θ F = dAp =
h dAt −
Ap
4
π Ap
d 2
tp 0 0
4
t− p
dt-p
h
θt= θ car cos θ t = cos θ p = = cos θ
dt − p
dAt
175
On obtient, du fait que
h h
r = h tgθ , dr = dθ , d t − p =
cos 2 θ cos θ
Rp
θ max
1 tgθ 1
F =  h4
cos 4 θ dθ
cos 2 θ h 4
Ap
dAt −
4
2 0
θ max h
1 1 θmax
=  tgθ cos θ dθ = sin 2 θ max
2
2 0
4
FdAt − Ap = sin 2 θ max = 0,328  θ max = 34,94°

R p = h tgθ max = 2 0,699 = 1,397 m
Ap = π R p = 6,13 m²
2
176
R5. Détermination des facteurs de vue par la

règle de Hottel
(échanges bi-dimensionnels; surfaces grises - isothermes)
A. Surfaces non concaves
A1
A1 F12 + A1 F13 = A1 (F11 = 0 )  Conservation de
A2 F21 + A2 F23 = A2 (F22 = 0 ) l’énergie (système A2 A3
A3 F31 + A3 F32 = A3 (F33 = 0) géométriquement fermé)
A1F12 = A2 F21 , A1F13 = A3 F31 , A2 F23 = A3 F32 Réciprocité
A1 + A2 − A3
F12 = etc....
2 A1
177
Prof. EL OTMANI 59
ENSAJ TD
B. Surfaces concaves ou localement concaves: plan de fermeture
Surface plane ou convexe traversée par une A3

rainure à faces planes ou convexes A1 et A2 de
mêmes propriétés et mêmes températures. A1 A2
On rend la surface non concave par un plan de ε1 = ε 2 = ε
fermeture A3 de mêmes propriétés ε et T, on
aura pour l’enceinte ainsi formée (A1- A2- A3): T1 = T2 = T
A3 e3 (T ) = A1 F13 e1 (T ) + A2 F23 e2 (T )
 e3 (T ) = ε e 0 (T )
Or e1 (T ) = e2 (T ) = ε e 0 (T ) Même propriété que la rainure
D ' où A3 e3 (T ) = ( A3 F13 + A2 F23 )ε e 0 (T ) Conclusion: pour les calculs

d’échange avec l’extérieur le plan
= A3 (F31 + F32 ) ε e 0 (T ) de fermeture A3 est équivalent à
1 car F33 =0 ( plan) la rainure (A1+A2).
178
C. Entre 2 surfaces non concaves quelconques (sans obstacle)
R A2 S
Pour PQRS
(enceinte fermée
A1F12 + A1 F1− PR + A1 F1−QS = A1
A1 convexe)
Pour PQR P
Q A1
A1 + APR − ARQ
A1F1− PR = (plan de fermeture)
2
Pour PQS
A1 + AQS − APS
A1F1−QS = (plan de fermeture)
2
179
ARQ + APS − (APR + AQS )

D ' où A1 F12 =
2
Règle de Hottel: le facteur de vue entre une surface émettrice (1) et

une surface réceptrice (2) en géométrie bidimensionnelle s ’obtient par
le quotient de (la somme des segments croisés moins la somme des
segments non-croisés tendus entre les extrémités de ces surfaces) par
deux fois la surface émettrice.
En termes de longueurs des arcs, on a aussi
= droite
= arc développé
____

____

____ ____
RQ + PS −  PR + QS 
 =  diag −  cotés
Seulement valable pour
F12 =  des surfaces PLANE
____
2 PQ 2A i ou CONVEXE
180
Prof. EL OTMANI 60
ENSAJ TD
Cas des obstacles

Soit la situation suivante, où 3 est un obstacle semi-infini
2
b
c’’
a g
d e c’ 3
1
1’ 1’’
Décomposons 1 en 1’ et 1’’
181
2 2
b
c’’ c’’
a g g
d e 3 d e c’ 3
1’ 1’’
g + d − (a + e + c′′) F1′′ 2 =
1
 [sin θ (s ) − sin θ (s )]ds
F1′ 2 = Hottel 2 L1′′
1 2
2 L1′ L1′′
g + d − (a + e + c′′) =
1
[(b − g ) + (e + c′′ − c′ − c′′)]
F21′ = Réciprocité
2 L1′′
2L2
=
[b + e − (g + c′)] Théorie p191
2 L1′′
θ1 θ2
ds b + e − ( g + c′ )
F21′′ = Réciprocité
2L2
182
b + d − (a + c′ + c′′)
Au total, F21 = F21′ + F21′′ =
2L2
L2 b + d − (a + c′ + c′′) b + d − (a + c′ + c′′)
F12 = =
L1 2 L2 2 L1
2
b Tout revient donc à utiliser la règle
c’’
a d’Hottel classique en considérant la
corde coupée par l’obstacle comme
d c’ 3 un élastique qui contourne celui-ci
1
1’ 1’’
Remarque sur les surfaces concaves: A2

A2’
F1−2 = F1−2′
mais F21 < F2′−1 car F22 ≠ 0 A1

183
Prof. EL OTMANI 61
ENSAJ TD
Applications
Deux facettes dans des plans normaux
A1 F12 = A2 F21
P
 
h 1  _____ _____ _____ _____ 
1 =  PR + QS − QR − PS 
Q R S 2 
 
___
RS
2
=
1 2
a + h2 + b − (a + b )2 + h 2 
a b 2  
Si a = 0 ( facettes jointives ) A1 F12 = A2 F21 =

1
2
[h + b − b 2 + h2 ]
184
Deux cylindres parallèles de même diam. d distants de a (centre à centre)
a
A
B C
d
_____ _____
A1 F12 = A2 F21 = 2 arc( AB ) + 2 BC − a arc ( AB ) =
d d 1 2
arcsin BC = a −d2
2 a 2
d
A1 F12 = A2 F21 = d arcsin + a2 − d 2 − a A1 = A2 = πd
a
1 a
2
arcsin +   − 1 − 
d a
F12 = F21 =
π a d d
 
185
R6. Détermination des facteurs de vue pour les géométries suivantes

1
F11 = F22 = F33 = F44 = 0 2
Fc . opposés càd F12 , F21 , F34 , F43 :

4 3 a
2a 2 − 2a
Fc . opposés = = 2 − 1 = 0,4142
2a 1
Fc . adjacents càd F13 , F31 , F14 , F41 , F23 , F32 , F24 , F42 : a
2a − a 2 2
Fc . adjacents = = 1− = 0,2929
2a 2
N .B. : on a bien F
j
ij = 1 ∀i
2 2
ex : i = 1 : F11 + F12 + F13 + F14 = 0 + 2 − 1 + 1 − +1− =1
2 2
186
Prof. EL OTMANI 62
ENSAJ TD
2
2
F11 = F22 = F33 = F44 = 0
4 3 a
2a 2 − 2a 2 −1
F12 = = x
2 2
 2x   2x 
2a 1 +   1+   1
a
 a   a  x≤
a 2
2
 2x 
a + a 1 +   − 2a
 a  1 2 −1
F13 = = − = F14
 2x 
2 2  2x 
2
2a 1 +   2 1+  
 a   a 
2 a 2 − 2a
F21 = = 2 −1
2a
2
F23 = F32 = F24 = F42 = 1 −
2
187
2
 2x 
a + a 1+   − a 2
1− 2 1
2
 a   2x 
F31 = F41 = = + 1+  
2a 2 2  a 
2
 2x 
2a 2 − a − a 1 +   2
 a  1 1  2x 
F34 = F43 = = 2− − 1+  
2a 2 2  a 
N .B. :" vérifications":
2 −1 1 2 −1 1 2 −1
*) F11 + F12 + F13 + F14 = 0 + + − + − =1
 2x 
2 2  2x 
2 2  2x 
2
1+   2 1+   2 1+  
 a   a   a 
2 2
F21 + F22 + F23 + F24 = 2 − 1 + 0 + 1 − +1− =1
2 2
188
1− 2 1
2 2
 2x  2 1 1  2x 
F31 + F32 + F33 + F34 = + 1+   +1− +0+ 2 − − 1+   = 1
2 2  a  2 2 2  a 
= F41 + F42 + F43 + F44
* *) lim F12 = 2 − 1
x →0
1 2 −1 2
lim F13 = lim F14 = − = 1−
x →0 x →0 2 2 2
1− 2 1 2
lim F31 = lim F41 = + = 1−
x →0 x →0 2 2 2
1 1
lim F34 = lim F43 = 2 − − = 2 − 1
x →0 x →0 2 2
189
Prof. EL OTMANI 63
ENSAJ TD
3 2
4 3 a
1’
F1′ 2 = F21 = F34 = F43 = 2 − 1 x
2 a
F1′3 = F1′ 4 = F31 = F41 = F23 = F32 = F24 = F42 = 1 − 1 x≤
2 2
a
a 2 −1
F12 = F1′ 2 =
2 2
 2x   2x 
a 1+   1+  
 a   a 
2
1−
a 2
F13 = F14 = F1′ 3 =
2 2
 2x   2x 
a 1+   1+  
 a   a 
190
2
1−
2 −1 2 1
F11 ≠ 0 : F11 = 1 − F12 − F13 − F14 = 1 − −2 = 1−
2 2 2
 2x   2x   2x 
1+   1+   1+  
 a   a   a 
F22 = F33 = F44 = 0
N .B. :
*) on a bien A1 F12 = A2 F21 ; A1 F13 = A3 F31

2
2 −1 = a ( 2 −1) 2  2
 2x  2 1− = a  1− 
a 1+   .  2x 
 a  2 a 1+   . 2
 2 
 2x   a 
1+ 
2
  2x 
 a  1+ 
 a 
2
* *) lim F12 = 2 − 1 ; lim F13 = lim F14 = 1 − ; lim F11 = 0
x→0 x→0 x→0 2 x →0
191
R7. Rayonnement d’un feu de forêt
Lors d ’une incendie de forêt en terrain plat, on observe un front de

flamme rectiligne de hauteur h. On assimile le front et le sol à deux
plans orthogonaux avec ε = α = 1 (surfaces noires). On demande:
a) la fraction de l ’énergie rayonnée par le front qui irradie le sol non
dévasté jusqu ’à une distance l du front de flamme;
b) la densité locale de débit de chaleur transmise au sol devant la
flamme si la température de celle-ci est de 1200°c;
c) l ’influence sur a) d ’une inclinaison de 30° du front de flamme par
suite d ’un vent violent.
flamme
h
sol
l
192
Prof. EL OTMANI 64
ENSAJ TD
a A f F fs = As Fsf =
1
2
[
h + l − h2 + l 2 ]
[ 1 l
] 
2
1 l
F fs = h + l − h 2 + l 2 = 1 + − 1 +   
2h 2 h h 
 
l 0, 25 0,50 0,75 1 2 5 10 ∞
h
F fs 0,110 0,191 0,250 0,293 0,382 0,450 0,475 0,5
N.B. La moitié de l’irradiation totale du sol est confinée dans la zone

0 ≤ x ≤ l = 0,75h
b Densité totale
de flux de chaleur
=
dQ fs d
dl
=
dl
A f F fs σ T f
4
( )
[ σ Tf
]
σ Tf  
4 4
d l
= h + l − h2 + l 2 = 1 − 2 2 
dl 2 2  h +l 
193
σ Tf 4
T f = 1200 + 273 = 1473K  = 133.464kW / m 2
2
l 0 0,25 0,50 0,75 1,00

h
Densité σT f 4
1 0,757 0,553 0,400 0,293

locale 2
c
 1 l 
2
1  _____
l l
F fs = h + l − AB  = 1 + − 1 +   −   
B 2h   2  h  h   h 

flamme
30°
h
l 0,25 0,50 0,75 1,00 ∞
h
F fs 0,174 0,317 0,424 0,5 0,75
sol
A
194
R8. Échange mutuel entre facettes grises en

enceinte fermée en équilibre radiatif
Rappel
1) toute facette i de l’enceinte est le siège d’un double phénomène:
a) Réception de rayonnement
= (1 − α i )Ei + α i Ei si opaque
Ei = ρ i Ei + α i Ei 
= (1 − ε i )Ei + ε i Ei si grise
Ei éclairement de la facette αi absorptivité
ρi réflectivité
εi émissivité
b) Émission de rayonnement
ei = ε i σ Ti
4
195
Prof. EL OTMANI 65
ENSAJ TD
2) On définit la radiosité de la facette comme la totalité du flux

énergétique qui la quitte
Ri = ρ i Ei + ε i σ Ti
4
= (1 − ε i )Ei + ε i σ Ti
4
3) Les radiosités associées aux facettes d’une enceinte fermée en

équilibre sont régies par les équations (Ai désignant l’aire de la facette i):
R j − Ri σ Ti − Ri
4

j ≠i rij
+
rii
=0 Pour la facette i à Ti connue
R j − Ri Pour la facette i qui perd le débit

 j ≠i rij
+ Qi = 0
de chaleur connu Qi
1− εi 1 1
avec rii = et rij = =
ε i Ai Ai Fij Aj Fji
196
En conclusion, on associe à chaque facette i:

a) si Ti connue:
- un potentiel thermique σTi4
-un nœud fictif de radiosité de potentiel thermique Ri, séparé
du premier par la résistance thermique propre de la facette rii
(qui s’annule dans le cas d’une facette noire pour laquelle Ri=
σTi4)
b) si Qi connu
-on injecte directement le débit de chaleur Qi au droit du nœud
de radiosité de potentiel Ri
Dans les deux cas a) et b), tous les nœuds de radiosité sont couplés
mutuellement par des résistances thermiques rij, fonctions des facteurs
de vue des facettes les unes par rapport au autres.
197
Applications
1. Échange de chaleur entre 2 plans indéfinis
T1 T2 r11 r12 r22
σT14 R1 R2 σT24
ε1 ε2
Q12 = Aq12 =
différence potentiel thermique σ T1 − T2
=
( 4 4
) F12 = F21 = 1
 résistances en série r11 + r12 + r22
=
σ T1 − T2 ( 4
=σ ε12
4
* ) A T1 − T2
4 4
( )
1 − ε1 1 1 − ε 2
+ + ε1ε 2 (ε 1 +ε 2 −ε1ε 2 )
Aε1 A Aε 2
198
Prof. EL OTMANI 66
ENSAJ TD
2. Échange entre 2 cylindres concentriques
ε1, T1 r11 r12 r22

r2
σT14 R1 R2 σT24
r1
ε2, T2
Q12 =
(
σ T1 − T2
4 4
)  F12 = 1

;  AF r
 A = 2π r1 L
;  1
r11 + r12 + r22 F = 1 12 = 1  A2 = 2π r2 L
 21 A2 r2
Q12 =
(
σ T1 − T2
4
=
σ ε1 A1
4
) T1 − T2
4 4
( )
1 − ε1 1 1 − ε 2 εA
+ + 1 + 1 1 (1 − ε 2 )
Aε1 A Aε 2 ε 2 A2
199
3. Échange au sein d ’un four industriel
On considère la section droite d’un four

industriel (carré de coté a=3m) dont les 1100°c 3
parois réfractaires présentent toutes une
émissivité ε=0,5 et en régime équilibré 1370°c 1370°c
les températures suivantes:
4 2
plafond : 1100°c
260°c 1
parois : 1370°c a
sol : 260°c
On demande de calculer le débit de chaleur perdu par le sol quand les

parois du four sont en équilibre radiatif, le vide étant réalisé dans le
four ou le fluide contenu étant neutre du point de vue radiatif (ni
émission ni atténuation de rayonnement)
200
Le modèle analogue au four se présente comme suit:

r33
σT34 plafond : T3 = 1100 + 273 = 1373K
R3 Les deux parois latérales sont
r23 regroupées sous le numéro 2.
r22
r13 σT24 parois : T2 = 1370 + 273 = 1643K
R2
r12
R1
r11
σT14 sol : T1 = 260 + 273 = 533K
Évaluons par unité de longueur de four normale au dessin:

1− ε 0,5
r11 = r33 = = = 0,333 m −1
ε A 0,5 (3x1)
1− ε
r22 = = 0,166 m −1 (2 parois )
ε 2A
201
Prof. EL OTMANI 67
ENSAJ TD
Pour les facteurs de vue, on utilise la méthode de Hottel
3 F11 = F33 = 0 ( planes )

a 2 2a 2 − 2 a
F13 = F31 = F22 = = 2 − 1 = 0,414
2a
4 2
1 F12 = F32 = 1 − F13 = 0,586
Ce qui, partant de 1, n’arrive pas en 3, arrive
a forcément sur les parois latérales
1 1
r12 = r23 = = = 0,569 m −1
A1 F12 3 0,586
Dès lors
1 1
r13 = = = 0,805 m −1
A1 F13 3 0,414
202
On peut alors directement écrire le système aux radiosités:
 1
 (R1 − R2 ) + 1 (R1 − R3 ) + 1 R1 = 1 σ T14
 r12 r13 r11 r11
1
 (R2 − R1 ) + (R2 − R3 ) +
1 1 1
R2 = σ T2
4
 r12 r23 r22 r22

 1 (R − R ) + 1 (R − R ) + 1 R = 1 σ T 4
r 3 1
r23
3 2
r33
3
r33
3
 13
C’est un système linéaire de 3 équations à 3 inconnues R1, R2, R3.
Une fois résolu, on trouve par exemple
Q1 =
1
r11
(σ T14 − R1 ) Débit de chaleur à
travers le sol
203
Alternativement pour le système simple à 3 facettes qui nous occupe,

on peut transformer les schémas de résistance comme suit:
a) Passage triangle - étoile (pour le cœur du système)
R3
R3
r3
r23
r2
r13 R2 R2
r12
r1
R1
R1
r12 r13 r1 = r3 = 0,236 m −1

r1 = etc… 
r12 + r13 + r23 r2 = 0,166 m −1
204
Prof. EL OTMANI 68
ENSAJ TD
b) Système en étoile final

Chaque branche de l’étoile est complétée par la résistance propre de
facette, si bien qu’on obtient:
σT34
r3’ = r3+r33 = 0,236+0,333 = 0,569 m-1
σT24
r2’ = r2+r22 = 0,166+0,166 = 0,333 m-1
r1’ = r1+r11
σT14
= 0,236+0,333
= 0,569 m-1
205
Alternativement pour le système simple à 3 facettes qui nous occupe,

on peut transformer les schémas de résistance comme suit:
a) Passage triangle - étoile (pour le cœur du système)
σT34
σT34
r 3’ ℜ 23
r 2’
σT24 ℜ 13 σT24
r 1’ ℜ 12
σT14
σT14
−1 −1 −1 −1
 r ′   r ′   r ′  =  r ′  = 1,757 m
ℜ12
−1
= 1  2 etc …   
1
 3
−1 −1 −1 −1
 r  +  r ′  +  r ′ 
′ ′
 r  = 3,0 m
 1
 2
 
3
 2 
206
−1 −1
ℜ12 = ℜ 23 = 0,809 m
−1
ℜ13 = 0,474 m
−1
(
Qsol = Q13 + Q12 = ℜ13 σ T3 − σ T1 + ℜ12 σ T2 − σ T1
4 4
) −1
( 4 4
)
( )
= 0,474 5,67 13,73 − 5,33 + 0,809 5,67 16,43 − 5,334
4 4
( 4
)
= (93,34 + 330,56) = 423,9 kw / m
207
Prof. EL OTMANI 69
ENSAJ TD
R9. Un mât télescopique pour antenne déployé à partir de la

soute d’une navette spatiale est constitué d’un cylindre en matériau
composite de 45cm de diamètre moyen. Les surfaces interne et externe
du cylindre sont assimilables à des corps gris opaques possédant une
émissivité ε=0,8. On suppose que l’écoulement du chaleur se fait
essentiellement dans une section droite du cylindre et que la
conduction circonférentielle est négligeable. Le vide interplanétaire est
assimilable à un corps noir à T∞=3K.
A) Déterminer la température uniforme d’équilibre des parois
du mât lors de la phase de déployment, qui se fait à l’abri de
l’irradiation solaire et dissipe en régime permanent une puissance
uniformément répartie de 250 w par mètre de hauteur du mât.
208
B) En négligeant la résistance thermique propre des parois et

en l’absence d’obstacle interne, déterminer les températures moyennes
d’équilibre T1 et T2 des deux faces (exposée et non exposée) du mât
après déployment (Qgén=0) lorsqu’il est exposé à l’irradiation solaire
définie par une flux incident qs égal à la constante solaire et dirigé
normalement à l’axe du mât.
T = 3K
ε = 0,8 q s = 1353 w / m 2
ε = 0,8
φ = 45 cm T2 T1
Q gén = 250 w / m
(A) (B)
209
1ère partie
 Téq   T 4 4 
Qgén = ε σ *   −  ∞  π φ par mètre de hauteur
 100   100  
4 4
 Téq  Qgén T 
  = + ∞ 
 100  ε σ π φ  100 
*
4
250  3 
= +  = 38,986
0,8 5,67 π 0,45  100 
Téq = 249,9 K = −23,3 °c

210
Prof. EL OTMANI 70
ENSAJ TD
2ème partie
Pour se ramener à un problème thermique unidimensionnel, il faut calculer q0,
équivalent radial du flux incident qs.
A) Flux irradiant p-u-surface
Rem: On considère que les deux coquilles, exposées et q0
non-exposées, sont chacune à température uniforme.
q s = 1353 w / m 2
R R
θ
ε =α π
2
π R q0 = α R  qs cosθ dθ
π
−
2
α 0,8
q0 = 2 qs = 2 1353 = 699 w / m 2
π π
211
B) Facteur de vue entre 2 1/2 coquilles cylindriques

(Hottel ou cours théorique)
2 2R 2
F12 = = = 0,637
2 πR π
C) Modèle thermique
πφ Surface d’une demi
S= =π R coquille n’intervient pas
2
1− ε 1− ε 1 1− ε 1− ε r1ext = 1
r2ext = 1 r22 = r22 = r12 = r11 = r11 =
ε ε F12 ε ε
σ T∞ 4 R2ext R2int R1int R1ext σ T∞ 4

σ T2 4 σ T14
q0
212
Équilibre thermique par unité de surface
(
q0 − ε σ T1 − T∞ =
4 4
) ε σ F12
ε + 2(1 − ε ) F12
(
T1 − T2 = ε σ T2 − T∞
4 4 4 4
) ( )
Flux net entrant Flux sortant
Flux échangé entre 1 et 2
surface 1 surface 2
T1 T2 T∞
Posons u1 = u2 = u∞ =
100 100 100
( )
 689 − 4,536 u1 − u ∞ = 2,739 u1 − u 2 = 4,536 u 2 − u ∞
4 4
( 4 4
) ( 4 4
)
7,275 u1 − 2,739 u 2 = 689 + 4,536 u∞ 4 4 4
syst de 2 éqs à 2 inc. 
 7,239 u1 − 2,275 u 2 = −4,536 u ∞
4 4 4
T2 = 253,9 K = −19,3°c
u2 = 41,55
4
 T1 = 324,1K = +51,0°c
u1 = 110,35
4
∆T12 = 70,3°c
213
Prof. EL OTMANI 71
ENSAJ TD
Contrôle
(
q1 = q0 − ε σ T1 − T∞ ≅ 188,5 w / m 2
4 4
)
q12 = ε * F12 σ (T 1
4
)
− T2 ≅ 188,5 w / m 2
4
(
q2 = ε σ T2 − T∞ ≅ 188,5 w / m 2
4 4
)
ok
214
Remarque importante:
On aurait pu faire l’exercice en injectant, d’une façon plus
systématique, le flux q0 au nœud de radiosité R1ext. Mais, il convient
alors de ne plus introduire de facteur d’absorption α dans le calcul de
q0: 689
q0 = = 861,25 w / m 2
0,8
1− ε 1− ε 1 1− ε 1− ε r1ext = 1
r2ext = 1 r22 = r22 = r12 = r11 = r11 =
ε ε F12 ε ε
σ T∞ 4 R2ext R2int R1int σ T∞ 4

σ T2 4 σ T14 R1ext
q0
Q
215
( )(=) R 1 −−εσT σT14 − σT∞ 4

ext 4
q0 + σT∞ − R1 = =Q
4 ext 1 1
1 (2 ) 1 2(1 − ε ) 1
+ +
ε ε ε 0,637
(1)  q0 + σT∞ − R1
4 ext
= 4 R1 ( ext
− σT1
4
)
 R1 = 172,25 + 0,8 σT1
ext 4
 R1
ext
=
1
5
(q0 + σT∞ + 4σT1
4 4
) (3)
(2) (
 4 R1 − σT1 = σT1 − σT∞ 0,301218
ext
) (
4 4
) 4
− σT1 = (σT − σT∞ 0,0753 ) (4)

ext 4 4 4
 R1 1
(3) et (4)  172,25 + 0,8σT1 − σT1 = (σT − σT∞ 0,0753 )

4 4 4 4
1
 172,25 = 0,2753σT1
4
 T1 = 324,1 K
216
Prof. EL OTMANI 72
ENSAJ TD
(3)  R1
ext
= 672,733 w / m 2 m 2
R1 − σ T1
( )
ext 4
 Q= = 4 R1 − σ T1 = 188,5 w / m 2
ext 4
1− ε
ε
σ T2 4 − σ T∞ 4
pour T2 , on a  Q =
1−ε
Q
 σ T2 = + σ T∞
4 4
ε
 T2 = 253,9 K
217
On peut démontrer qu’il est normal d’avoir les mêmes résultats par les
deux méthodes.
Soit A’ la grandeur associée à la première méthode.
Soit A’’ la grandeur associée à la deuxième méthode.
Première méthode:
ε  ′ext ′4  ′ext
 R1 − σ T1  = σ T∞ − R1
4
1− ε  
′ext 1 ε
σ T1′
4
 R1 = σ T∞ +
4
1− ε 1− ε
′ext ′4
 R1 = (1 − ε )σ T∞ + ε σ T1
4
ε  ′ext ′4 
donc, Q′ = α q0 +  R1 − σ T1 
1− ε  
ε2 ε
σ T1′ − σ T1′
4 4
Q′ = ε q0 + ε σ T∞ +
4
1− ε 1− ε
218
Seconde méthode:
″ext ε  ″ext ″4 
q0 + σ T∞ − R1 =  R1 − σ T1 
4
1− ε  
″ext 1 ε
σ T1″
4
 R1 = q0 + σ T∞ +
4
1− ε 1− ε
″ext ″4
 R1 = (1 − ε )q0 + (1 − ε )σ T∞ + ε σ T1
4
1  ″ext ″4 
et , Q′′ =  R1 − σ T1 
1− ε  
ε2 ε
σ T1″ − σ T1″
4 4
Q′′ = ε q0 + ε σ T∞ +
4
1− ε 1− ε
Donc, on a bien Q’=Q’’ car, quelque soit l’une des méthodes, T1 est
donné par la même équation:
σ T14 − σ T∞ 4 ′ ″
Q=  T1 = T1 = T1
1 2(1 − ε ) 1
+ +
ε ε F12
219
Prof. EL OTMANI 73
ENSAJ TD
220
IV. Problèmes combinés
221
G1. Un cylindre, de rayon r=2cm, est chauffé par une résistance

électrique dissipant une énergie de 1 kW par mètre. Ce cylindre est
entouré d ’une enveloppe cylindrique concentrique en acier
(k=45W/mK) dont le rayon intérieur vaut 10cm, et dont les surfaces
sont assimilables à des corps gris possédant une émissivité égale à 0,8,
sont séparés par un vide tel que les échanges convectifs sont refroidie
par un fluide à T∞=50°c, dont le coefficient de convection vaut
1kw/m²K.
On demande de calculer les températures de paroi T(1), T(2), T(3), en
régime permanent.
fluide
2cm 10cm
(3) (2) (1)
50°c acier
4mm
vide
222
Prof. EL OTMANI 74
ENSAJ TD
Considérons d’abord l’enveloppe cylindrique extérieure (analogie

électrique) :
 
ln  3  
r
 1 r
T(2 ) − T∞ = Q  +  2 
 2 π r3 l h 2π k l 
ln  3 
r
1  r2   
T∞ hA 2π k l  ln  0,104  
Q  0,1 
+ 
1
=50°c T(3) T(2) =  
2 π l 0,104 1000 45
 
 
Q
= 0,010847
2π l
Or, conservation du débit de chaleur:

Débit de chaleur rayonné = débit de chaleur donné au fluide
T(2 ) − 50 T(3 ) − 50
T(2 ) − T∞ 1000 = 1000 l =
Qr = 1 1
R 0,010487
2π hA
et Qr = 1000 l
T( 2 ) = 51,669 °c T(3) = 51,53 °c
223
Échange radiatif entre deux cylindres concentriques:

cfr p R26:
Q12 = 1000 l =
σ ε1 A1
ε1 A1
(
T(1) − T(2 )
4 4
)
1+ (1 − ε 2 )
ε 2 A2
d ' où 1000 =
5,67 10 −8 0,8 2π 0,02
1 + (0,2)
2
T(1) − T(2 )
4 4
( )
d ' où T(1) = 663,31K
= 390,16°c
Rq.: si ε1=1, alors T(1)=650,9K<663,31K ce qui est logique car le

corps noir émet un maximum
224
G2. Refroidissement d’une plaque
Une plaque d’acier d’épaisseur 1mm, ayant subi un traitement

thermique (Tpl=500°c), est suspendue horizontalement au centre d’un
hall d’usine. Les dimensions d’une section droite du hall, et la position
de la plaque sont données sur le dessin.
On demande de calculer le débit de chaleur perdu par la plaque par
convection (corrélations p.153), et par rayonnement, en faisant
l’approximation d’un problème purement bidimensionnel (extension
très grande normalement au plan du dessin) et en considérant que l’air
est neutre du point de vue radiatif.
La plaque présente une émissivité ε=0,8. Les températures des parois
du hall - considérées comme des surfaces noires - sont les suivantes:
plafond: 5°c, mur de gauche: 10°c, à droite ouverture avec air
extérieur à 0°c, sol: 15°c. L’air dans le hall est à 10°c.
225
Prof. EL OTMANI 75
ENSAJ TD
air 10°c air 0°c

6m
0,5m
15m
226
A. Plaque horizontale soumise à la convection libre

g β ∆T L3
nombre de Grashof : GrL =
ν2
L = long . caractéris tique = 0,5m
avec ∆T = 500 − 10 = 490
1
β = coeff . d ' exp . thermique =
Tm
On évalue les prop. de l’air à la t° moyenne:
T pl + T∞
Tm = = 255°c = 528 K
2
k − 0,03951 528 − 500 k = 0,04124
table A4 → = →
0,0426 − 0,03951 50 (interpolat ion linéaire )
α = 5,96944 10−5
ν = 4,16648 10 −5
ν
Pr = 0,698 =
α
227
D’où 1
9,81 490 0,53
GrL = 528 = 6,555 108
(4,16648 10 −5
2
)
Ra L = GrL Pr = 6,555 108 0,698
= 0,4576 109 < 109 → Laminaire
face sup. : instable : c1 = 0,54
 
 face inf. : stable : c1 = 0,27
a) face stable
Nu L = 0,27 Ra L = 39,489
0 , 25
Nu k fl 39,489 0,04124
h= = = 3,257 W / m² K
L 0,5
Q = hA(Tpl − T∞ ) = 797,98W / m
228
Prof. EL OTMANI 76
ENSAJ TD
a) face instable
Nu L = 0,54 Ra L = 78,978
0 , 25
Nu k fl78,978 0,04124
h= = = 6,514 W / m² K
L 0,5
Q = hA(Tpl − T∞ ) = 1,596 kW / m
D’où: débit de chaleur perdu par la plaque par convection

(par unité de l):
Q = 2,3939 kW / m
229
B. Échange radiatif plaque - parois du hall

Détermination des facteurs de vue: 1
2 a 3
b
4
0,5 + L1 + L1 − L2 − L1
Fa −2 = = 0,5 + L1 − L2
2 0,5
L1 = 32 + 7,252 = 7,84618 L2
avec L1
L2 = 3 + 7,75 = 8,31039
2 2
6
Fa −2 = Fa −3 = Fb −2 = Fb −3 = 0,03579 0,5
2 L − 2 L1 7,25
Fa −1 = 2 = 0,92842 = Fb − 4
2 0,5
Fa −4 = Fb −1 = 0
230
plafond : T1 = 278 K à droite : T3 = 273K

mur à gauche : T2 = 283K sol : T4 = 288 K
(plafond) σT14
ra1
Ra
ra2 ra3
raa
σT2 4 σT34
σTpl4
(gauche) (droite)
rbb
rb2 rb3
Rb
rb4
(sol) σT44
231
Prof. EL OTMANI 77
ENSAJ TD
avec: 1− ε 1 − 0,8
raa = = = 0,5 = rbb
ε Aa 0,8 0,5 1
1
ra1 = rb 4 = = 2,15421
Aa Fa1
1
ra 2 = ra 3 = rb 2 = rb 3 = = 55,87787
Aa Fa − 2
Du graphique des résistances, on tire:

1
raa
(σTpl − Ra =
4
) 1
ra1
Ra − σT1 +
4
( 1
ra 2
)
Ra − σT2 +
4 1
ra 3
(
Ra − σT3
4
) ( )
1 1 1 1  Tpl 4 T14 T2 4 T3 4 
Ra  + + +  =σ + + + 
r r r r
 aa a1 a 2 a 3   raa ra1 ra 2 ra 3 
 
 Ra
232
Avec Ra, on trouve le débit de chaleur perdu par la face a par

rayonnement (par unité de long.):
σ Tpl 4 − Ra
Qa ray = = 7,9622 kW / m
raa
De même avec la face b, d’où:
σ Tpl 4 − Rb
Qb ray = = 7,9431 kW / m
rbb
Donc, le débit de chaleur perdu par la plaque par rayonnement:

15,905 kW / m
Débit total de chaleur perdu par la plaque:

Q = Qconv + Qray = 18,299 kW / m
233
Ou: échange radiatif, 2ème méthode:

Débit de chaleur perdu par la plaque (par unité de longueur):
Qout = σ 0,8 Ta 0,5 + σ 0,8 Tb 0,5 = 16,195 kW / m

4 4
Débit de chaleur reçu par la plaque (par unité de longueur):
[
Qin = ε σ T1 A1 F1−a + σ T4 A4 F4−b + σ T2 ( A2 F2−a + A2 F2−b )
4 4 4
]
+ ε [σ T ( A3 F3− a + A3 F3−b )]
4
3
= 290,06W / m
Q = Qout − Qin = 16,195 − 0,29 = 15,905 kW / m
234
Prof. EL OTMANI 78
ENSAJ TD
G3. Quel est le flux de chaleur q (w/m²) qu’il faut au moins évacuer par
les serpentins de refroidissement logés dans la glace d’une patinoire,
pour éviter que celle-ci ne fonde alors qu’elle est recouverte
provisoirement d’un plancher de bois (kb=0,2w/mK) pour une
représentation théâtrale.
Le plancher (eb=2,5cm) est en contact avec de l’air (Text=20°c,
hext=5w/m²K) et reçoit simultanément, par rayonnement dû à
l’éclairage, un flux de chaleur net qe=150w/m².
Il est d’autre part séparé de la surface de glace par un espace d’air de
2cm présentant une conductance thermique cair=1,3w/m²K et supporté
par des lisses longitudinales de section carrée, faites du même bois,
avec un pas de 60cm.
Les contacts lisses-plancher et lisses-glace sont supposés parfaits.
L’écoulement de chaleur est supposé parfaitement vertical.
235
Les serpentins de refroidissement sont enfuis à une profondeur d=3cm

sous la surface de la glace (kg=2,215w/mK). Ce sont de longs conduits
cylindriques de diamètre φ=2cm, disposés parallèlement aux lisses du
plancher, avec un pas de 3cm et parcourus par du fréon 12 à TF=-13°c.
Leurs parois sont supposées parfaitement lisses et de résistance
thermique négligeable.
Déterminer le débit maximum de fluide nécessaire pour une patinoire
de largeur l=6m transversalement aux serpentins, en supposant que
l’écoulement est parfaitement développé et que l’écoulement de chaleur
reste mono dimensionnel et que les serpentions captent la totalité du
flux de chaleur incident (pas de flux de perte au-delà des serpentins).
236
Text= 20°c
qe=150 w/m²
Air hext= 5w/m²K
Bois eb=2,5cm kb=0,2w/mK

Air + cair=1,3w/m²K
lisse eair=2cm
Glace d=3 cm kg=2,215w/mK
Serpentins
3 cm
... φ=2 cm
TF= -13°c
60 cm
237
Prof. EL OTMANI 79
ENSAJ TD
Schéma électrique équivalent:
Text= 20°c 1
a) Rext = = 0,2 m 2 K / w
Rext hext
e 2,5 10 − 2
qe b) Rb = b = = 0,125 m 2 K / w
Tbsup kb 0,2
Rb  p. u. surface plancher
Tbinf
Réq.
Tgsup= 0°c
238
c) Réq. résulte de deux résistances thermiques en

parallèle pour 0,6m de longueur et 1m de
profondeur:
Rb’ Rair
 ′
 Rb′ =
eb
= 5K / w
1
=
1
+
1
avec 
 kb 2 10 −2 1 ( )
Réq R ′ Rair R = 1
= 1,326 K / w
( )
b
 cair 58 10 − 2 1
air
1 1 1
= + = 0,954 w / k
Réq 5 1,326
2
→ Réq = 1,048K / w = 0,629 m K / w
0,6m de longueur m de largeur
239
d) Détermination de la température du plancher
qe + qconv = qcond
ext
− Tg
sup sup
T −T
sup
Tb
150 + ext b =
Rext Rb + Réq.
250 − 5 Tb = 1,326 Tb
sup sup
 Tb = 39,5°c
sup
e) Flux atteignant la glace
q = qe + qconv = qcond = 52,4 w / m 2

ext
(la convection naturelle a pour effet de réduire l’effet

du rayonnement par un facteur 2/3)
240
Prof. EL OTMANI 80
ENSAJ TD
2ème partie: Le flux, calculé ci-dessus, est conduit jusqu’à une

profondeur d = 3cm
− Ts
sup
Tg
q = kg en w / m 2
d
Où Ts désigne la température moyenne de la glace au droit des
serpentins
qd 52,4 3 10 −2
Ts = Tg − =− = −0,7°c
sup
kg 2,215
Pour une largeur de 60cm, on a 20 serpentins longitudinaux et
l’ensemble du flux conduit dans la glace est repris par le fréon, donc:
Qconv = q 0,6
= h 20 S serp (TS − TF )
avec S serp = π φ par mètre longitudin al (surf . lat .)
241
q 0,6
d ' où : h = = 2,03 w / m 2 K
20 π 2 10 −2 12,3
L’écoulement de fréon est laminaire. Si on en doute, on le suppose
turbulent, dès lors (pour des parois isothermes):
h φ 2,03 0,02
Nu D = = = 0,465
kf 0,0874
1
  0 ,8
 Nu D 
Re D =  1/ 3
 = 25,7 << 2100
 0,023 Pr 
 3, 4 à 13°c 
242
Pour un écoulement laminaire, le coefficient d’échange réel s ’obtient

par (pour des parois isothermes)
hφ 3,658 k f
Nu D = = 3,658  hr = = 15,985 w / m 2 K
kf φ
La vitesse moyenne maximum d’écoulement est
2100 ν f
uav / lim = = 2,37 10 − 2 m / s Pour rester en laminaire
φ
Le débit de fréon pour 6 m de largeur (débit massique) est
π φ2
mɺ F = 10 20 uav / lim ρ = 2,14 kg / s
4 1438kg / m ³
243
Prof. EL OTMANI 81
ENSAJ TD
Contrôle final:
Avec le coefficient d ’échange hr, la température superficielle de glace
Tgsup est-elle inférieure à 0°c?
− Tg − TF
sup sup sup
Text − Tb
sup
Tb Tg
→ 150 + = =
Rext Rb + Req Rg + Rserp
eg 0,03
Rg = = = 0,0135 m 2 K / w
kg 2,215
1 20 hr π φ
Rserp = = 0,0299 m 2 K / w avec h* = = 33,48W / m² K
h* 0,6
= 37,4°c
sup
Tb
 (ok!)
= −10,26°c
sup
Tg
244
V. Exercice résolus
bouquin de théorie
245
Exercice 7:
Un réservoir métallique sphérique (D=0,5m) est testé en laboratoire en
vue de stocker de l’azote liquide à 77K pour une application en
cryogénie spatiale. Le réservoir est pourvu d’une enveloppe isolante de
25mm d’épaisseur, à base de poudre de silice (k=2 10-3w/mK).
L’ambiance du laboratoire est de l’air calme à 300K. La chaleur latente
de vaporisation et la densité de l’azote liquide valent respectivement
200kJ/kg et 804kg/m³.
Quel sont, en régime permanent unidimensionnel, l’apport de chaleur à
l’azote et ses taux de vaporisation massique et volumique?
On négligera la résistance thermique de la paroi métallique du réservoir
ainsi que celles des contacts azote-réservoir et réservoir-enveloppe. On
supposera (en le vérifiant a posteriori) que l’écart de température entre
la paroi externe de l’enveloppe et l’ambiance est au plus de ∆T=1K. On
utilisera la corrélation de l’ingénieur en régime laminaire pour la
convection naturelle.
246
Prof. EL OTMANI 82
ENSAJ TD
air
azote
D = 0,5m
Tazote = 77 K paroi isolant
eiso = 25mm réservoir
Tamb = 300 K
Lazote = 200 KJ / kg
ρ azote = 804kg / m3
247
Exercice 6:
Une conduite d’air chaud d’un échangeur à tube présente une section
semi-circulaire de R=20mm. La face plane est maintenue à Tb=1000K
tandis que le pourtour circulaire est parfaitement isolé (température Tp).
Le débit d’air (à pression atmosphérique) est de 10-2kg/s et
l’écoulement parfaitement développé avec une température moyenne de
section Ts=400K.
Du point de vue radiatif, les parois internes sont grises avec ε=0,8.
Quelle est la quantité de chaleur emportée par le fluide par mètre
courant de canalisation pour maintenir Tb à sa valeur, et que vaut Tp?
Que devient Ts dans le sens de l’écoulement? Justifier.
248
Air Conduite
Ts=400K Tp
mɺ = 10 −2 kg / s ε=0,8
p=1atm
R=20mm
Bas
Tb=1000K
ε=0,8
249
Prof. EL OTMANI 83
ENSAJ TD
Exercice 1:
Une cuillère à café (matériau conducteur métallique k=40w/mK,
barreau de L=10cm, P=5cm, A=0,5cm²) est plongée pour moitié dans
du café chaud, ce qui contribue à lui apporter une quantité de chaleur
superficielle uniformément répartie qs=100w/m². L’autre partie de la
cuillère se refroidit à l’air ambiant (Ta=20°c, h=6w/m²K).
On demande, à partir d’un modèle discret basé sur des résistances
thermiques et injecteurs, la distribution de température longitudinale de
la cuillère, sachant que ses deux extrémités sont adiabatiques. Contrôler
si possible analytiquement.
Conseil: Prendre 3 nœuds le long de la cuillère et concentrer les
échanges convectifs transversaux et/ou l’injection de chaleur en ces
nœuds.
250
Ta=20°c
h=6w/m²K
L/2
qs=100
w/m² P
L/2
A
251
Exercice 3:
On contrôle thermiquement le fini de surface (dépôt d’un revêtement)
d’un collecteur solaire concave (surface 1) par irradiation infrarouge à
partir d’un panneau chauffant plan (surface 2), le tout dans un local
dont l’ensemble des parois est à Tp=20°c.
Les caractéristiques thermiques des surfaces et la géométrie sont
précisées sur la figure.
On demande le flux net de chaleur transmis à la surface du collecteur,
en régime permanent et par unité de longueur, en supposant le problème
2-D, les effets convectifs sont négligeables et pas d’atténuation du
rayonnement par l’air ambiant.
252
Prof. EL OTMANI 84
ENSAJ TD
Collecteur:
T1=300°c
ε=0,5
L=1,5m
Tp Tp=20°c
h=1m
l=1m
panneau:
T2=700°c
ε=0,9
253
Exercice 10:
Un mur d’habitation (longueur 10m, hauteur 2m50) est constitué d’une
couche de brique (épaisseur 12cm), tapissée intérieurement d’une
couche de laine de roche (épaisseur 6cm), elle-même recouverte de
plaques de plâtre (épaisseur 1,5cm).
L’air à l’extérieur est à une température uniforme de 2°c et circule
horizontalement le long de la paroi de brique à une vitesse uniforme de
50km/h. L’ambiance interne de l’habitation est à une température
uniforme de 27°c.
Donner la perte de chaleur a travers le mur en supposant l’écoulement
de chaleur unidimensionnel transversal en régime permanent, avec et
sans effets convectifs. Donner les différentes températures aux parois et
aux jonctions internes
254
h=2,5m
Air
u=50km/h Air
Ta=2°c Brique Isolant Plâtre Ta=27°c
eb=0,12m ei=0,06m ep=0,015m
kb=0,3w/mK ki=0,03w/mK kp=0,5w/mK
L=10m
255
Prof. EL OTMANI 85
ENSAJ TD
VI. Exercice d’examen

exercice n°1 - janvier 2001
256
Une plaque d’aluminium pur (h=2m, e=5mm, L=5m) à T0=20°c

uniforme est intensément et uniformément échauffé par effet Joule
(chauffage inductif) d’une quantité Q=16,2Mw pendant ti=3s,
insuffisant pour permettre un refroidissement substantiel (qu’on
négligera donc), et conduisant à un nouvel état thermique T1 supposé
uniforme, à calculer.
La plaque est ensuite plongée soudainement (trempe) pour moitié dans
un bain d’eau à Te=0°c, et refroidie pour l’autre moitié par convection
forcée (souffle d’air horizontal à 10m/s) et par rayonnement (ε=0,15)
avec l’ambiance du laboratoire à Ta=20°c.
Donner une évaluation du temps nécessaire pour retrouver un état
thermique stable à Tm (à définir) dans la plaque (on néglige la
conduction) en évaluant les coefficients d’échanges convectifs à partir
de T1. Cette évaluation est-elle optimiste où pessimiste? Justifier.
257
Par comparaison, donner la distribution finale de T sur la hauteur de la

plaque par résistances thermiques, en négligeant les pertes radiatives
(ceci est-il justifié?).
Dire comment elle se compare à Tm. Dans tous les cas, on considérera
que les distributions de températures ne varient pas significativement
dans le sens de l’épaisseur et de la longueur, et on négligera les pertes
par les bords de la plaque. Il n’y a pas d’interaction radiative sensible
entre la partie supérieure de la plaque et les bords du récipient d’eau.
258
Prof. EL OTMANI 86
ENSAJ TD
10m/s
Ta=20°c
ε=0,15
e=5mm
Q=16,2Mw
ti=3s
h=2m
L=5m
Te=0°c
259
260
Prof. EL OTMANI 87

Microsoft PowerPoint - TDde Thermique - td2021

Transféré par

Droits d'auteur :

Formats disponibles

Microsoft PowerPoint - TDde Thermique - td2021

Transféré par

Informations du document

Titre original

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

Microsoft PowerPoint - TDde Thermique - td2021

Transféré par

Droits d'auteur :

Formats disponibles

ENSAJ TD

Exercices de transferts de chaleur

Table des matières

Qu ’est qu’un régime permanent?

Conduction: qu’est qu’il faut savoir?

C1: Loi de Fourier

Soit une plaque de conductibilité k =34W/mK, de

On demande de calculer le débit de chaleur dans

C2 Loi de Fourier / milieu Composite

La température de la face gauche de cette paroi est Tg=400°c,

On demande de déterminer les températures des deux faces

Inox Cuivre Inox

 ∆TI = 145,03°c TC g = 254,97°c qI = 1,232 106W

= h(T − T f ) (flux de conduction sortant du

QV = 0 pas de source ni de puit

• Conditions aux limites:

Contact non parfait

Les conditions en x = L1 sont maintenant remplacées par:

q1 ( x = L1 ) = qint = cint (T1 (L1 ) − T2 (L1 ))

C ’est à dire: − a k1 = cint (a L1 + b − c L1 − d )

q1 = q2 = qint = 0,1(537,5 − 412,5) = 12,5W / cm²

Conditions aux limites: (4)

T1 ( L1 ) = T2 (L1 ) = 895 .45 °c Temp . de jonction

Remarque: Analogie électrique --->> discrétisation

Calcul des résistances

C2 Loi de Fourier / milieu Composite

Inox Cuivre Inox

 ∆TI = 145,03°c TC g = 254,97°c qI = 1,232 106W

C2: Confirmation par l’analogie électrique

= h(T − T f ) (flux de conduction sortant du

• C4: Contrôle par l ’analogie électrique

R1=L1/k1A R2=L2/k2A Rext=1/hA

Les conditions en x = L1 sont maintenant remplacées par:

q1 ( x = L1 ) = qint = cint (T1 (L1 ) − T2 (L1 ))

• Variante par l ’analogie électrique:

Rab1 Rb1b2 Rb2c Rcf

Variante b: contrôle par l’analogie électrique

T (°c) Milieu Milieu

C6: Les ponts thermiques

Soit une paroi relativement isolante traversée par

Quelle est l ’influence de ces ponts thermiques sur

Rem: on considérera l ’écoulement 1D suivant l ’axe des

Am = (1 − ε )A0 Section de la paroi non perturbée au sein de A0

Analogie électrique - suite

C7. Influence d ’un isolant sur le refroidissement

Quelle est la température à l ’interface graphite - verre?

= 5.88 103 m = r3 (cf . th. p 71)

Conduction en régime variable

• Cas du milieu isotrope homogène:

Analyse transitoire où la conduction est négligeable

Qin − Qout + Qgen = −  q.next dA +  QV dV

– A: surface externe du milieu (surface d’échange effectif)

Cas du chauffage convectif

C10.1 Problème du thermocouple

Initialement à T0=500°c et plongé dans une ambiance

Quelle est la réponse thermique de la plaque utilisée

• Hypothèse: - chauffage convectif

 T = T∞ + (T0 − T∞ ) e −0.4 t = 1000 − 500 e −0.4 t

• i.e. Plausibilité de l’étude moyenne