Université Sultan Moulay SMP_S6_ Parcours Energétique

Slimane M33 : Transfert de Chaleur

Faculté Polydisciplinaire Année Universitaire 2017– 2018
Béni Mellal
M. LAMSAADI
Examen
Session ordinaire (durée 2h30mn)

Exercice 1 : (8 points)
Une ailette cylindrique pleine, de rayon R et de longueur L, est installée sur une pièce
métallique dont la température est égale à Tp. L'air environnant est à la température T∞ et le
coefficient d’échange convectif à la surface de l'ailette est h. On suppose que :
- le transfert conductif est unidirectionnel ( T = T ( x ) ) à travers l’ailette ;
- le régime est permanent sans source interne ;
- h, T∞ et Tp sont les mêmes pour le système sans ailette que pour celui avec ailette ;
- l’ailette est homogène et isotrope ;
- le bout de l’ailette est adiabatique.
1) Faire un bilan thermique sur un élément de volume approprié de l'ailette et déterminer l'équation
différentielle que doit vérifier la température au sein de cette ailette.
2) Déterminer le profil de température au sein de l'ailette, T ( x ) .
3) Exprimer le flux de chaleur, Q0, évacué par l’ailette
4) Calculer le rendement et l’efficacité de l’ailette. On donne :
R = 10mm, L= 0.5 m, Tp = 90°C, T∞ = 30°C, h = 35 W/(m2.K) et λ = 385W/(m.K)

T∞ , h

L adiabatique
λ
Tp 0 x
R

T∞ , h

1/2
Exercice 2 : (8 points)
On considère une enceinte fermée constituée par quatre surfaces S1, S2, S3 et S4 formant un
cylindre de rayon R et de longueur L (voir la figure ci-dessous) (S3 est la surface latérale et (S1+S4)
et S2 sont les deux surfaces de base où S1 est un disque de rayon r). On suppose que les surfaces S1,
S2 et S3 sont des corps noirs maintenus, respectivement, à des températures constantes T1, T2 et T3.
La surfaces S4 est supposée grise, opaque à émission et réflexion diffuses et adiabatique dont
l’émissivité est ε4. On donne : r = 10 cm et R = 120 cm, L = 2m, T1 = 600°C, T2 = 0°C, T3 = 100°C
et ε4 = 0.5
1)- Calculer les différents facteurs de forme mis en jeux.
S2
2)- Calculer le flux total quittant chaque surface Si, Q iq .

3)- Calculer le flux total incident sur chaque surface Si, Q ii .

4)- Déduire la température de la surface S4 et le flux net
perdu ou reçu par chacune des surfaces S1, S2 et S3. L S3
-8 2 4
σ =5.67 10 W/(m .k ) : (constante de Stefan-Boltzmann)
Le facteur de forme entre deux disques parallèles Di
et Dj est exprimé par :

rj S1
S4
(D ) j

ri
Fij =
1
2
[
h − h 2 − 4a 2 b 2 ]
(D i ) a=
rj
L
, b=
L
ri
( )
et h = 1 + 1 + a 2 b 2

Exercice 3 : (convection naturelle le long d’une plaque carrée) (4 points)

On considère une plaque carrée de côté L = 1m, dont l’une de ses deux faces est adiabatique
et l’autre est maintenue à une température constante Tp = 385K. Cette plaque est placée dans l’aire
à pression atmosphérique et à Tf = 315K. Les propriétés physiques de l’air atmosphérique à la
température moyenne sont données par :
ν = 2.076×10-5m2/s, λ = 0.03W/m.°C, Pr = 0.697 et β = 2.86×10-3K-1
Calculer le flux de chaleur échangé par convection naturelle.
1

Le nombre de Nusselt moyen est donné par la corrélation suivante : Nu m = 0.14Ra 3

L

Les valeurs numériques doivent être données avec quatre chiffres significatifs après la virgule
2/2
 Correction d’examen Session Ordinaire 2017-2018

Exercice 1 :
1)
- Régime permanent
- Sans source interne
- λ est constante
- T = T(x) h, T∞

Tp
λ

∆x x
x=0 x=L
h, T∞

Le bilan thermique sur le volume de contrôle ∆v = πR2∆x, donne :

Q cond
e - Q scond - Q sconv = 0

où Q cond
e = πR 2 q x (x) , Q scond = πR 2 q x (x + ∆x) et Q sconv = 2πR∆xh (T ( x ) − T∞ )

πR 2 q ( x ) − πR 2 q ( x + ∆x ) − 2πR∆xh (T ( x ) − T∞ ) = 0

dq ( x )
⇒ − πR 2 ∆x − 2πR∆xh (T ( x ) − T∞ ) = 0
dx
d 2 T(x )
⇒ λπR 2 ∆x − 2πR∆xh (T ( x ) − T∞ ) = 0
dx 2
d 2 T ( x ) 2h
⇒ − (T(x ) − T∞ ) = 0
dx 2 Rλ
2h
Si on pose θ( x ) = T( x ) − T∞ et m 2 = , l’équation précédente s’écrit :
Rλ
d 2 θ( x )
2
− m 2 θ( x ) = 0
dx
Il s’agit d’une équation différentielle linéaire, du second ordre à coefficients constants, qui a
pour solution générale :
θ( x ) = C1e mx + C 2 e − mx
où C1 et C2 sont des constantes d’intégration à déterminer à partir des conditions aux limites
imposées.

3/2
 - En x = 0 : T(x = 0) = Tp ⇒ θ( x= 0) = Tp - T∞ = θp.
∂T ∂θ
- En x = L : =0 ⇒ =0
∂x ∂x
L’application de ces dernières conditions, permet d’obtenir :
C1 + C 2 = θ p

C1e mL − C 2 e − mL = 0

 θp
C1 =
1 + e 2 mL
⇒ 
C = θp
 1 + e − 2 mL
2

θp θp e mx + e m ( x − 2 L ) + e − mx + e − m (x − 2 L )
⇒ θ( x ) = e mx + e − mx ⇒ θ( x ) = θ p
1 + e 2 mL 1 + e − 2 mL 2 + e 2 mL + e − 2 mL

Coh (mx ) + Coh (m( x − 2L))

⇒ ⇒ θ( x ) = θ p
1 + Coh (2mL)

Coh (mx ) + Coh (m( x − 2L))

⇒ T( x ) = (Tp − T∞ )
2h
+ T∞ avec m =
1 + Coh (2mL) Rλ
2)
Le flux de chaleur, Q0, évacué par l’ailette peut être calculé à partir du flux entrant par
conduction ou du celui dû au coefficient d’échange h avec le milieu extérieur :
L
dT( x ) 
Q 0 = −λπR  ou Q 0 = ∫ h 2πR [T( x ) − T∞ ]dx
2

dx  x =0 0

dθ( x ) 
L
⇒ Q 0 = −λπR 2  ou Q 0 = ∫ 2hπRθ( x )dx
dx  x =0 0

Sh (2mL)
Q 0 = θ p λπR 2 m
1 + Coh (2mL)
3)
Si l’ailette est idéale (c'est-à-dire, elle était faite d’un matériau infiniment conducteur (de
conductivité infinie)), alors elle serait isotherme à T = Tp. Le flux de chaleur évacué par l’ailette
serait donc :
L L
Q id = ∫ h 2πR (Tb − T∞ )dx = ∫ h 2πRθ p dx = 2hπRLθ p
0 0

4/2
 Le rendement d’ailette, η, est défini comme le rapport du flux réellement évacué par elle
(Q0) et celui qui serait évacué si elle est idéale (Qid). Donc, le rendement de l’ailette est exprimé
par :
Q0 1 Sh (2mL)
η= = ⇒ η = 0.456
Q id mL 1 + Coh (2mL)

L’efficacité, εailette, d’une ailette est définie comme le rapport entre le flux réellement évacué
par l’ailette (Q0) et celui qui serait évacué sans ailette à travers la surface S (Qsa) :
λm Sh (2mL)
Q sa = hπR 2 (Tb − T∞ ) = hπR 2 θ p ⇒ ε ailette = ⇒ ε ailette = 45,6
h 1 + Coh (2mL)

Exercice 2 :

a) Les quatre surfaces constituent une enceinte fermée. Donc on écrit :

4

∑F
j=1
ij =1

⇒ F11 + F12 + F13 + F14 = 1 , F21 + F22 + F23 + F24 = 1 , F31 + F32 + F33 + F34 = 1 , F41 + F42 + F43 + F44 = 1

Les surfaces S1, S2 et S4 sont planes, donc : F11 = F22 = F44 = 0 .

Les surfaces S1 et S4 sont dans le même plan, alors : F14 = F41 = 0

⇒ F12 + F13 = 1 , F21 + F23 + F24 = 1 , F31 + F32 + F33 + F34 = 1 , F42 + F43 = 1

F12 =
1
2
[ ] ( )
h − h 2 − 4a 2 b 2 , a = 0.6 , b = 20 et h = 1 + 1 + a 2 b 2 = 545 ⇒ F12 = 0.2643

⇒ F13 = 0.7357

S1 r2 S r2
F21 = F12 = 2 F12 ⇒ F21 = 0.0018 , F31 = 1 F13 = F13 ⇒ F31 = 0.0015 ,
S2 R S3 2RL

F2 (1+ 4) =
1
2
[ ]
h − h 2 − 4a 2 b 2 ⇒ F2 (1+ 4 ) = 0.2194 ⇒ F2 (1+ 4 ) = F21 + F24 ⇒ F24 = 0.2176

F23 = 1 − F21 − F24 = 0.7806

S2 R2
F42 = F24 = 2 F24 = 0.2191 , F43 = 1 − F42 = 0.7809
S4 R − r2

S2 S
F32 = F23 = 0.2342 F34 = 4 F43 = 0.2326 , F33 = 1 − F31 − F32 − F34 = 0.5317
S3 S3

5/2
F11 = F22 = F44 = F14 = F41 = 0 , F12 = 0.2643 , F13 = 0.7357 , F21 = 0.0018 , F31 = 0.0015 ,

F24 = 0.2176 , F23 = 0.7806 , F42 = 0.2191 , F43 = 0.7809 , F32 = 0.2342 , F34 = 0.2326 , F33 = 0.5317

2)- flux total quittant chaque surface Si, Q iq .

S1, S2 et S3 sont des surfaces noires : Q iq = Q ei = Si σTi4

⇒ Q1q = 1035.1157W , Q q2 = 1426.8628W , Q 3q = 16568.1938W

S4 n’est pas noire : Q q4 ≠ Q ei : Q q4 = S 4 J 4

Q 4 net = ∑ F4 jS 4 (J 4 − J j ) = 0 ⇒ F42 S 4 (J 4 − J 2 ) + F43S 4 (J 4 − J 3 ) = 0

4

j=1

F42 J 2 + F43 J 3
⇒ J4 = = F42 J 2 + F43 J 3
F42 + F43 = 1

Q q2 Qq
⇒ J 4 = F42 + F43 3
S2 S3

Q q2 Qq
⇒ Q q4 = S 4 F42 + S 4 F43 3 ⇒ Q q4 = F24 Q q2 + F34 Q 3q = 4164.2472W
S2 S3

4 4
3)- Flux total incident sur chaque surface Si, Q ii : Q ii = ∑ Q qji = ∑ Fji Q qj
j=1 j=1

Q1i = F21Q q2 + F31Q 3q ⇒ Q1i = 27.4206 W

Q i2 = F12 Q1q + F32 Q 3q + F42 Q q4 ⇒ Q i2 = 5066.2386W

Q i3 = F13 Q1q + F23 Q q2 + F33 Q 3q + F43 Q q4 ⇒ Q i3 = 13936.5130 W

Q i4 = F24 Q q2 + F34 Q 3q ⇒ Q i4 = 4164.2472 W = Q q4

1 1
σT 4 − J 4  J 4  Q 4
q
4)- - Q 4 net = i = 0 ⇒ T4 =  4  =  4  ⇒ T4 = 357.5747K
1 − ε4 σ  σS 4 
ε 4S 4

- Q1net = Q1q − Q1i = 1007.6951W

- Q 2 net = Q q2 − Q i2 = -3639.3759W

- Q 3net = Q 3q − Q i3 = 2631.6808W

6/2
Exercice 3 : (convection naturelle le long d’une plaque carrée)
h L Nu m λ
Nu m = m → h m =
λ L

Nu m = 0.14Ra 1L/ 3 → Nu m = 0.14(Pr GrL ) = 0.14 Pr 1 / 3 GrL

1/ 3 1/ 3

λ gβ(Tp − Tf )L3
⇒ hm = c Pr 1 / 3 GrL avec GrL = = 5.7 x10 8
1/ 3

L ν 2

0.03
hm = 0.14(0.697)1 / 3 (5.7 x10 8 )1 / 3 = 6.18W / m 2 .°C
0 .5
⇒ Q = h m S∆T = 432.6W

7/2
 Correction d’examen Session Ordinaire 2017-2018

Exercice 1 :
1) Le bilan thermique sur le volume de contrôle ∆v = πR2∆x, donne :
d 2 T ( x ) 2h
πR 2 q ( x ) − πR 2 q ( x + ∆x ) − 2πR∆xh (T( x ) − T∞ ) = 0 ⇒ − (T(x ) − T∞ ) = 0
dx 2 Rλ
d 2 θ( x ) θp θp
− m 2 θ( x ) = 0 ⇒ θ( x ) = C1e mx + C 2 e − mx ⇒ θ( x ) = e mx + e − mx
dx 2
1+ e 2 mL
1 + e −2 mL
Coh (mx ) + Coh (m( x − 2L)) Coh (mx ) + Coh (m( x − 2L))
⇒ θ( x ) = θ p ⇒ T( x ) = (Tp − T∞ ) + T∞
1 + Coh (2mL) 1 + Coh (2mL)

2h
m=
Rλ
L
dT( x )  Sh (2mL)
2) Q 0 = −λπR 2  ou Q 0 = ∫ h 2πR [T( x ) − T∞ ]dx ⇒ Q 0 = θ p λπR m
2

dx  x =0 0
1 + Coh (2mL)
L
Q0 1 Sh (2mL)
3) Q id = ∫ h 2πR (Tb − T∞ )dx =2hπRLθ p ⇒ η = = ⇒ η = 0.456
0
Q id mL 1 + Coh (2mL)

λm Sh (2mL)
Q sa = hπR 2 (Tb − T∞ ) = hπR 2 θ p ⇒ ε ailette = ⇒ ε ailette = 45,6
h 1 + Coh (2mL)
Exercice 2 :

a) S1, S2 et S4 sont planes: F11 = F22 = F44 = 0 , S1 et S4 sont dans le même plan: F14 = F41 = 0

F12 =
1
2
[ ] 1
[ ]
h − h 2 − 4a 2 b 2 ⇒ F12 = 0.2643 , F2 (1+ 4) = h − h 2 − 4a 2 b 2 ⇒ F2 (1+ 4 ) = 0.2194
2
F13 = 0.7357 , F21 = 0.0018 , F23 = 0.7806 , F24 = 0.2176 , , F31 = 0.0015 , F32 = 0.2342 ,

F33 = 0.5317 , F34 = 0.2326 , F42 = 0.2191 , F43 = 0.7809

2)- Q iq = Q ei = Si σTi4 , i = 1-3

⇒ Q1q = 1035.1157W , Q q2 = 1426.8628W , Q 3q = 16568.1938W

Q q2 Qq
Q q4 = S 4 J 4 , Q 4 net = ∑ F4 jS 4 (J 4 − J j ) = 0 ⇒ J 4 = F42 J 2 + F43 J 3 ⇒ J 4 = F42
4
+ F43 3
j=1 S2 S3

Q q2 Qq
⇒ Q q4 = S 4 F42 + S 4 F43 3 ⇒ Q q4 = F24 Q q2 + F34 Q 3q = 4164.2472W
S2 S3

Q 4 net = 0 ⇒ Q q4 = Q i4 ⇒

Q q4 = Q14
q
+ Q q24 + Q 34
q
+ Q q44 = F14 Q1q + F24 Q q2 + F34 Q 3q + F44 Q q4 = F24 Q q2 + F34 Q 3q

8/2
 4 4
3)- Flux total incident sur chaque surface Si, Q ii : Q ii = ∑ Q qji = ∑ Fji Q qj
j=1 j=1

Q1i = F21Q q2 + F31Q 3q ⇒ Q1i = 27.4206 W

Q i2 = F12 Q1q + F32 Q 3q + F42 Q q4 ⇒ Q i2 = 5066.2386W

Q i3 = F13 Q1q + F23 Q q2 + F33 Q 3q + F43 Q q4 ⇒ Q i3 = 13936.5130 W

Q i4 = F24 Q q2 + F34 Q 3q ⇒ Q i4 = 4164.2472 W = Q q4

1 1
σT 4 − J 4  J 4  Q 4
q
4)- - Q 4 net = i = 0 ⇒ T4 =  4  =  4  ⇒ T4 = 357.5747K
1 − ε4 σ  σS 4 
ε 4S 4

- Q1net = Q1q − Q1i = 1007.6951W

- Q 2 net = Q q2 − Q i2 = -3639.3759W

- Q 3net = Q 3q − Q i3 = 2631.6808W

Exercice 3 : (convection naturelle le long d’une plaque carrée)

h L Nu m λ
Nu m = m → h m =
λ L
Nu m = 0.14Ra 1L/ 3 → Nu m = 0.14(Pr GrL ) = 0.14 Pr 1 / 3 GrL
1/ 3 1/ 3

λ gβ(Tp − Tf )L3
⇒ h m = c Pr GrL avec GrL = = 5.7 x10 8
1/ 3 1/ 3

L ν 2

0.03
hm = 0.14(0.697)1 / 3 (5.7 x10 8 )1 / 3 = 6.18W / m 2 .°C
0 .5
⇒ Q = h m S∆T = 432.6W

9/2