Université de Caen UFR des Sciences 1

Master 1 IMM mention Ingénierie Mécanique (M1)

Transferts de chaleur et de masse
TD2 - Convection Thermique I1

Exercice II.1 : On considère l’écoulement de Couette avec la plaque mobile (notée

plaque 1 avec vitesse = U ) est maintenue à une température uniforme (T1 ) et la plaque
immobile (notée plaque 2) est isolée thermiquement. Déterminer la température de plaque
2 ; exprimer le résultats en fonction des propriétés du fluide, la température et la vitesse
de plaque mobile, T1 et U .

Exercice II.2 : Des mesures expérimentales pour le coefficient de transfert thermique par
convection d’un écoulement croisé sur une
barre de section carrée (comme montré ci-
L = 0, 5 m
contre), ont conduit au résultats suivants :

h1 = 50 W/m2 .K lorsque V1 = 20 m/s

2
h2 = 40 W/m .K lorsque V2 = 15 m/s
Air, V
On admet que la forme fonctionnelle pour le
nombre de Nusselt est N u = CRem P rn , où
C, m et n sont des constantes.
1. Que serait le coefficient de transfert thermique par convection pour une barre similaire
si L = 1 m quand V = 15 m/s ?

2. Que serait le coefficient de transfert thermique par convection pour une barre similaire
si L = 1 m quand V = 30 m/s ?

3. Aurait-on obtenu les mêmes résultats si le coté de carrée est utilisé au lieu du diagonal
en tant que longueur caractéristique ?

Exercice II.3 : L’air atmosphérique s’écoule parallèlement (U∞ = 15 m/s, T∞ =

15◦ C) à un élément chauffant ayant une surface plane (d’aire S = 0, 25 m2 ) destinée à être
maintenue à une température uniforme de 140◦ C. On trouve que la force de traı̂née induite
par la surface de l’élément est 0, 25 N.
1
Ces exercices sont tirés de Chapitre 7 du livre :
Incropera, Dewitt, Bergman, Lavine, Fundamentals of heat and Mass Transfer, Wiley

Adil Ridha Département de Mathématiques et Mécanique

Université de Caen UFR des Sciences 2

Déterminer la puissance électrique nécessaire pour maintenir la surface à la tempéraure

prescrite. Pour les propriétés d’air on peut utiliser le tableau donné à la fin d’exercices.

Exercice II.4 : Les mesures pris dans une journée d’été indiquent que la température
est 27 ◦ C et l’humidité relative est 30%. On trouve que le taux d’évaporation de l’eau à la
surface d’une lac (dont la température est 27 ◦ C) est 0, 10 kg/hr.m2 .
Déterminer le coefficient de transfert de masse par convection.

Exercice II.5 : On observe que le taux de perte de masse d’eau dans une casserole
(diamètre D = 230 mm) est 1, 5 × 10−5 kg/s à 23 ◦ C.

1. Déterminer le coefficient de transfert de masse par convection sous de telles conditions.

2. Estimer le taux de perte de masse par évaporation lorsque l’humidité relative de l’air
vaut 50%.

3. Estimer le taux de perte de masse par évaporation quand les températures de l’air
et de l’eau sont à 47 ◦ C ; on admettra que le coefficient de transfert de masse par
convection reste constant.

Exercice II.6 : Une schématisation assez simple pour dessaler l’eau consiste en main-
tenant un film liquide mince de l’eau saline à la surface de la plaque inférieure de deux
grandes (infinies) plaques planes parallèles, légèrement inclinées et séparées par une dis-
tance L.

condensation
Air, V
TL
y

T0

x Film mince
d’eau saline

Un écoulement incompressible, laminaire et assez lent existe alors dans l’espacement

entre les deux plaques caractérisé par une composante finie de vitesse dans la direction
des x et de composantes nulles dans les directions y et z. L’évaporation a lieu au film
liquide à la surface inférieure, qui est maintenue à une température élevée T0 , tandis que la
condensation a lieu à la surface supérieure, qui est maintenue à une température réduite TL .

Adil Ridha Département de Mathématiques et Mécanique

Université de Caen UFR des Sciences 3

On désigne la concentration molaire de la vapeur d’eau à la surface inférieure par CA,0 et

celles à la surface supérieure par CA,L . On admettra que les distributions de la température
et de concentration dépendent de y uniquement.

1. Trouver une expression pour la distribution de la concentration molaire de vapeur

d’eau CA (y). Déterminer le taux de production massique de l’eau pure par unité d’aire
de surface. Exprimer votre résultats en fonction de CA,0 , CA,L , L, et le coefficient de
diffusion vapeur-eau DAB .

2. Déterminer l’expression donnant le flux thermique à fournir par unité de surface

pour maintenir la surface inférieure à la température T0 . Exprimer votre résultat en
fonction de CA,L , T0 , TL , L, DAB , hf g ( = hf − hg , la chaleur latente de vaporisation
de l’eau, hf ≡ hL ), et la conductivité thermique λ.

Exercice II.7 : On considère un écoulement croisé du gaz X sur un objet ayant une
longueur caractéristique L = 0, 1 m. Pour un nombre de Reynolds de 1 × 104 , on trouve
que le coefficient de transfert thermique moyen est 25 W/m2 .K. Le même objet est en suite
est mis dans le liquide Y et soumis aux mêmes conditions de l’écoulement. Étant données
les propriétés thermophysiques suivantes, déterminer le coefficient de transfert de masse
moyen.

ν (m2 /s) λ (W/m.K) α (m2 /s)

Gaz X 21 × 10−6 0, 030 29 × 10−6
Liquide Y 3, 75 × 10−7 0, 665 1, 65 × 10−7
Vapeur Y 4, 25 × 10−5 0, 023 4, 55 × 10−5
Mélange de gaz X–vapeur Y Sc = 0, 72

Exercice II.8 : Dans un processus industriel, on fait évaporer l’eau d’un film liquide
mince d’une surface à contour particulier. L’air sec s’écoule sur la surface, et les mesures au
laboratoire conduisent à la correlation suivante pour le coefficient de transfert thermique
par convection :
N uL = 0, 43ReL0,58 P r0,4

1. Que serait le taux d’évaporation pour une unité de surface (1 m2 ) et une longueur
caractéristique L = 1 m quand le fluide est l’air s’écoulement à la vitesse 10 m/s et

Adil Ridha Département de Mathématiques et Mécanique

Université de Caen UFR des Sciences 4

à la température 27 ◦ C ? On admettra que la masse volumique de la vapeur saturée

est ρA,sat = 0, 0077 kg/m3 .

2. Que serait la température du film liquide à l’état stationnaire pour cette configura-
tion ?

Propriétés thermophysiques de l’air à la pression atmosphérique

Température Densité Chaleur Conductivité Viscosité Coefficient no. de
spécifique thermique cinématique d’expansion Prandtl
◦
C (kg/m3 ) cp (kJ/kg K) λ (W/m K) ν (m2 /s)×10−6 β (1/K) ×10−3 Pr
-150 2.793 1.026 0.0116 3.08 8.21 0.76
-100 1.980 1.009 0.0160 5.95 5.82 0.74
-50 1.534 1.005 0.0204 9.55 4.51 0.725
0 1.293 1.005 0.0243 13.30 3.67 0.715
20 1.205 1.005 0.0257 15.11 3.43 0.713
40 1.127 1.005 0.0271 16.97 3.20 0.711
60 1.067 1.009 0.0285 18.90 3.00 0.709
80 1.000 1.009 0.0299 20.94 2.83 0.708
100 0.946 1.009 0.0314 23.06 2.68 0.703
120 0.898 1.013 0.0328 25.23 2.55 0.70
140 0.854 1.013 0.0343 27.55 2.43 0.695
160 0.815 1.017 0.0358 29.85 2.32 0.69
180 0.779 1.022 0.0372 32.29 2.21 0.69
200 0.746 1.026 0.0386 34.63 2.11 0.685
250 0.675 1.034 0.0421 41.17 1.91 0.68
300 0.616 1.047 0.0454 47.85 1.75 0.68
350 0.566 1.055 0.0485 55.05 1.61 0.68
400 0.524 1.068 0.0515 62.53 1.49 0.68

Adil Ridha Département de Mathématiques et Mécanique

Université de Caen UFR des Sciences 1

Master 1 IMM mention Ingénierie Mécanique (M1)

Transfert de Chaleur et de Masse
Corrigée de TD2 - Convection Thermique I

Solution de l’exercice - II.1:

Hypothèses :

1. Écoulement incompressible et laminaire en régime permanent.

2. Écoulement de Couette : écoulement entre deux plaques planes engendré par le mou-
vement d’une plaque relativement à l’autre avec ∂p/∂x = 0.

3. La plaque mobile est à tempéraure uniforme et la plaque immobile est isolée ther-
miquement.

4. Propriétés thermophysiques constantes.

Analyse :

Le profil de vitesse pour un écoulement de Couette est donné par

y
u(y) = U
L
Pour l’écoulement de Couette, l’équation de la conservation d’énergie se réduit à :
 2
∂ 2T ∂u
0=λ 2 +µ
∂y ∂y

Afin d’obtenir la distribution de température, on intègre deux fois cette équation

 2  2
∂ 2T µ U ∂T µ U
= − , =− y + C1
∂y 2 λ L  ∂y λ L
2
µ U
T (y) = − y 2 + C1 y + C2
2λ L

Ensuite,on applique les conditions aux limites associées :

∂T
= 0, T (y = L) = Tl , constante
∂y y=0

D’où
µ 2
C1 = 0, C2 = TL + U
2λ

Adil Ridha Département de Mathématiques et Mécanique

Université de Caen UFR des Sciences 2

La température de la plaque immobile (y = 0) est

µ 2
T (0) = TL + U
2λ
Et le flux thermique à la plaque mobile :
∂T µ 2
ϕ = −λ = U
∂y y=L L


Solution de l’exercice - II.2:

Hypothèses :

1. Correlation de la forme N u = CRem P rn , C, m et n étant des constantes.

2. Propriétés thermophysiques constantes.

Analyse :

1. D’après les expériences et les conditions L = 1 m,et V = 15 m/s, il en suit que P r

aussi bien que C, m et n son des constantes. Donc

hL ∝ (V L)m

Ainsi, en utilisant les résultats expérimentaux on peut déterminer m :

m m
h1 L1 V1 L1 50 × 0, 5 20 × 0, 5
 
= =⇒ =
h2 L2 V2 L2 40 × 0, 5 15 × 0, 5
D’où on tire m = 0, 782. Il vient alors que pour L = 1 m et V = 15 m/s, l’on obtient :
m 0,782
L1 V L 0, 5 15 × 1, 0
 
2
h = h1 = 50 W/m .K × = 34, 3 W/m2 .K
L V 2 L2 1, 0 20 × 0, 5

2. Pour la condition L = 1 m et V = 30 m/s, on trouve

m 0,782
L1 V L 0, 5 30 × 1, 0
 
2
h = h1 = 50 W/m .K × = 59, 0 W/m2 .K
L V 2 L2 1, 0 20 × 0, 5

3. Si on utilise le coté de carré comme longueur caractéristique, seule la constant C

changerait. Les coefficients m et n restent inchangés.

Solution de l’exercice - II.3:

Adil Ridha Département de Mathématiques et Mécanique

Université de Caen UFR des Sciences 3

Hypothèses :

1. Régime permanent.

2. L’analogie de Chilton–Colburn s’applique au problème.

3. La surface inférieure est isolée thermiquement.

Propriétés thermophysiques :

Du tableau de propriétés de l’air à Tf = 77◦ C et 1 atm : ρ = 0, 995 kg/m3 , cp =

1009 J/kg.K, P r = 0, 70, ν = 20, 22 × 10−6 m2 /s.

Analyse :

Le nombre de Reynolds est :

U∞ L 40 m/s × 0, 2 m
ReL = = = 4, 0 × 105 < 5, 0 × 105
ν 20, 22 × 10−6 m2 /s
ce qui montre que l’écoulement sur toute la plaque est laminaire. La vitesse U∞ est con-
stante et par conséquent le gradient de pression est nul.
Les valeurs moyennes de la contrainte parietale et le coefficient de frottement sont
respectivement donnés par :
FD 0, 25 N
τs = = = 1 N/m2
S 0, 25 m 2

et
τs 1 N/m2
Cf =1 = 3 = 8, 93 × 10−3 .
2
ρU 2
∞ 0, 5 × 0, 995 kg/m × (15 m/s) 2

On a de l’analogie de Chilton–Colburn
h C f −2/3
St = = Pr
ρU∞ cp 2
et par conséquent on obtient :

h = 0, 995 kg/m3 × 15 m/s × 1009 J/kg.K 8, 93 × 10−3 /2 (0, 70)−2/3 = 85 W/m2 .K

Donc, le flux thermique est :

Φ = hS (Ts − T∞ ) = 85 W/m2 .K × 0, 25 m2 × (140 − 15)◦ C = 2, 66 kW

Solution de l’exercice - II.4:

Adil Ridha Département de Mathématiques et Mécanique

Université de Caen UFR des Sciences 4

Schématique :

L’air T∞ = 27◦ C
ρA,s
Humidité relative, φ = = 0, 3 n′′A = 0, 1 kg/m2 .h
ρA,sat
ρA,∞ ρA,s

Eau, Ts = 27◦ C

Hypothèses :

1. Équilibre thermodynamique à la surface vapeur–liquide.

2. Conditions isothermes.

3. Comportement de gaz parfait pour la vapeur d’eau.

4. L’air est à la pression atmosphérique standarde.

Propriétés :

On obtient du tableau de vapeur d’eau saturée à 300 K :

pA, sat = 0.0357 bar, ρA, sat = 1/vg = 0, 025786 kg/m3 ,

Analyse :

L’équation pour le coefficient de transfert de masse (évaporation) peut s’écrire sous la

forme :
n00A
hm =
ρA,s − ρA,∞
où
ρA,s = ρA, sat ,

qui est la masse volumique de saturation à la température de l’eau. La définition de la

humidité relative (φ) se réduit, dans le cas où on suppose que la vapeur se comporte
comme gaz parfait, à :
ρA,∞
φ= =⇒ ρA,∞ = φρA, sat
ρA, sat

Adil Ridha Département de Mathématiques et Mécanique

Université de Caen UFR des Sciences 5

Il vient alors que :

n00A
hm =
ρA, sat (1 − φ)
Donc, on obtient in insérant les valeurs numériques :

0, 1 kg/m2 .h × 1/3600 s/h

hm = = 1, 54 × 10−3 m/s
0, 025786 kg/m3 (1 − 0, 3)

Solution de l’exercice - II.5:

Schématique :

Air, T∞ = 23◦ C
ou T∞ = 47◦ C ṁévap = 1, 5 × 10−5 kg/s
φ∞ = 0 ou 0, 5
ρA,s

Eau, D = 0, 23 m
Ts = 23◦ C ou 47◦ C

Hypothèses :

1. La vapeur d’eau est saturée à la surface de liquide et peut être traitée comme gaz
parfait.

Propriétés :

Du tableau de la vapeur d’eau saturée, on obtient à T = 296 K :

−1
ρA, sat = vg−1 = 48, 58m3 /kg = 0, 0206kg/m3

et à T = 320 K :

−1
ρA, sat = vg−1 = 13, 868m3 /kg = 0, 072kg/m3

Adil Ridha Département de Mathématiques et Mécanique

Université de Caen UFR des Sciences 6

Analyse :

1. Étant donnée que l’évaporation est un processus de transfert de masse, l’équation

pour la taux d’évaporation s’écrit alors sous la forme :

ṁévap = hm S (ρA, sat − ρA,∞ )

et le coefficient de transfert de masse est

ṁévap 1, 5 × 10−5 kg/s
hm = = = 0, 0175 m/s
(πD2 /4) (ρA, sat − ρA,∞ ) (π/4)(0, 23 m)2 0, 0206kg/m3
avec Ts = T∞ = 23◦ C et φ∞ = 0.

2. Si l’humidité relative de l’air ambiant est augmentée à 50%, le rapport des taux
d’évaporation est
ṁévap (φ∞ = 0, 5) hm S [ρA, s (Ts ) − φ∞ ρA, s (T∞ )] ρA, s (T∞ )
= = 1 − φ∞
ṁévap (φ∞ = 0) hm SρA, s (Ts ) ρA, s (Ts )
Donc
0, 0206kg/m3
" #
ṁévap (φ∞ = 0, 5) = 1, 5 × 10−5 kg/s 1 − 0, 5 = 0, 75 × 10−5 kg/s.
0, 0206kg/m3

3. Si la température de l’air ambiant est augmentée de 23◦ C à 47◦ C, avec φ∞ = 0

maintenu pour les deux cas, il vient que le rapport des taux d’évaporations est :
ṁévap (Ts = T∞ = 47◦ C) hm SρA,s (47◦ C) ρA,s (47◦ C)
= . = .
ṁévap (Ts = T∞ = 23◦ C) hm SρA,s (23◦ C) ρA,s (23◦ C)
Donc,
0, 072kg/m3
ṁévap (Ts = T∞ = 47◦ C) = 1, 5 × 10−5 kg/s 3 = 5, 24 × 10
−5
kg/s.
0, 0206kg/m


Solution de l’exercice - II.6:

Schématique :

y=L
CA,L , TL

Air (B), ṁB

Vapeur d’eau (A) To > TL
y
NA′′

CA,o , To
x

Adil Ridha Département de Mathématiques et Mécanique

Université de Caen UFR des Sciences 7

Hypothèses :

1. Régime permanent.

2. Écoulement incompressible, établi et à propriétés constantes.

3. Forces du volume négligeables.

4. Pas de réaction chimique.

5. Tous formes de travail, y compris de dissipation visqueuse, sont négligeables.

Analyse :

1. L’écoulement est entièrement établi aussi bien pour la concentration de vapeur que
pour les champs de vitesse et températures. Alors,
∂CA
= 0, ou CA (x, y) = CA (y).
∂x
L’écoulement est aussi en régime permanent ∂CA /∂t = 0, ṄA = 0, v = 0 avec
constant coefficient de diffusion DAB . Ains, l’équation pour la conservation d’espèces
se réduit à :
d2 CA
= 0.
dy 2
dont l’intégration (deux fois) conduit à :

CA (y) = C1 y + C2

Il vient en appliquant les conditions aux limites :

CA (y = 0) = CA,o =⇒ C2 = CA,o
CA,o − CA,L
CA (y = L) = CA,L =⇒ CA,L = C1 L + C2 , d’où C1 = −
L
Ainsi, on trouve pour la distribution d’espèces
y
CA (y) = CA,o − (CA,o − CA,L )
L
D’après la loi de Ficks, le flux d’espèces est donné par :
∂CA CA,o − CA,L
NA00 = NA,s
00
= −DAB = DAB
∂y y=0 L
En multipliant par la masse moléculaire de la vapeur d’eau, MA , on obtient le taux
temporel de production d’eau par unité de surface :
CA,o − CA,L
n00A = MA NA00 = MA DAB
L

Adil Ridha Département de Mathématiques et Mécanique

Université de Caen UFR des Sciences 8

2. La quantité de chaleur fournie à la surface inférieure doit être égale à la somme de

chaleur latent et de flux thermique sensible de surface :

ϕ = ϕlat + ϕsen
dT
 
00
ϕ = nA,s hf g + −λ .
dy y=0

La distribution de la température peut être obtenue en utilisant l’équation d’énergie,

qui, pour les conditions prescrites, se réduit à :

d2 T
=0
dy 2
dont la solution générale est donnée par

T (y) = C1 y + C2 .

Les conditions aux limites conduisent au calcul des constantes :

T (y = 0) = 0 =⇒ C2 = To
TL − T0
T (y = L) = TL =⇒ C1 =
L
On obtient alors la distribution de température :

y
T (y) = To − (To − TL ) .
L
D’où :
dT (To − TL )
−λ =λ .
dy y=0 L
Finalement, on tire

CA,o − CA,L (To − TL )

ϕ = MA DAB hfg + λ .
L L

Adil Ridha Département de Mathématiques et Mécanique

Université de Caen UFR des Sciences 9

Solution de l’exercice - II.7:

Schématique :

Sec, hL = 25 W/m2 .K
X(gaz) ou liquide Y , hm,L
Tsat = T∞
T∞
ReL = 1 × 104

L = 0, 1 m

Hypothèses :

1. Analogie de transfert de chaleur et masse s’applique au problème.

2. La vapeur Y se comporte comme gaz parfait.

Propriétés :

(Données) ν (m2 /s) λ (W/m.K) α (m2 /s)

Gaz X 21 × 10−6 0, 030 29 × 10−6
Liquide Y 3, 75 × 10−7 0, 665 1, 65 × 10−7
Vapeur Y 4, 25 × 10−5 0, 023 4, 55 × 10−5
Mélange de gaz X–vapeur Y : Sc = 0, 72

Analyse :

L’analogie entre le trasferts de chaleur et de masse peut être exprimer sous la forme :

hL L hm,L L
N uL = = F (ReL , P r), ShL = = F (ReL , Sc)
λ DAB
Les conditions de l’écoulement sont les mêmes dans les deux cas. Vérifions les valeurs de
P r et Sc. Pour P r, les propriétés pour gaz X(B) :

νB 21 × 10−6 m2 /s
Pr = = = 0, 72
αB 29 × m2 /s

Adil Ridha Département de Mathématiques et Mécanique

Université de Caen UFR des Sciences 10

tandis que Sc = 0, 72 pour le mélange de gaz X(gaz)–vapeur Y (A). Il vient alors que pour
cette situation
hL L hm,L L DAB
N uL = = ShL = ou hm,L = hL .
λ DAB λ
En se rendant compte de

νB 21 × 10−6 m2 /s
DAB = = = 29, 67 × 10−6 m2 /s
Sc 0, 72
et puis en substituant les valeurs numériques, l’on obtient :

29, 67 × 10−6 m2 /s
hm,L = 25 W/m2 .K × = 0, 0243 m/s
0, 03 W/m.K

Solution de l’exercice - II.8:

Schématique :

ϕévap
Air
ϕconv Film d’eau, Ts
T∞ = 300 K
NuL = 0, 43Re0,58
L Pr
0,4
φ∞ = 0
V = 10 m/s S = 1m2 , L = 1 m
p = 1 atm.

Hypothèses :

1. Régime permanent.

2. Propriétés constantes.

3. Transfert thermique par rayonnement négligeable.

4. L’analogie entre les transferts de chaleur et de masse s’applique au problème.

Propriétés :

Tableau des propriétés d’air à 300 K, 1 atm. :

λ = 0, 0263 W/m.K, ν = 15, 89 × 10−6 m2 /s, P r = 0, 707

Adil Ridha Département de Mathématiques et Mécanique

Université de Caen UFR des Sciences 11

Tableau d’eau à Ts ≈= 280 K :

vg = 130, 4 m3 / kg, hfg = 2485 kj/kg.

Tableau d’air-eau à T = 298 K : DAB = 0, 26 × 10−4 m2 /s.

Analyse :

1. Le taux d’évaporation est

ṁévap = nA = hm S [ρA,sat (Ts ) − φ∞ ρA,sat (T∞ )] = hm SρA,sat (Ts )

ù D’après l’analogie entre les transferts de chaleur et de masse

ShL = 0, 43Re0,58 0,4

L Sc .

On a
VL 10m/s1 m ν 15, 89 × 10−6 m2 /s
ReL = = = 6, 29×105 , Sc = = = 0, 61
ν 15, 89 × 10 m /s
−6 2 DAB 0, 26 × 10−4 m2 /s

ShL = 0, 43(6, 29 × 105 )0,58 (0, 61)0,4 = 814

DAB 0, 26 × 10−4 m2 /s
hL = ShL = (814) = 0, 0212 m/s
L 1m
ρA,sat (Ts ) = vg−1 (Ts ) = 0, 0077 kg/m3 .

Donc

ṁévap = 0, 0212 m/s × 1 m2 × 0, 0077 kg/m3 . = 1, 63 × 10−4 kg/s.

2. Le bilan d’énergie à la surface conduit à ϕconv = ϕévap , ou

00
00
ṁévap hfg
hL (T∞ − Ts ) = ṁévap hfg =⇒ Ts = T∞ −
hL
avec
N uL = 0, 43(6, 29 × 105 )0,58 (0, 61)0,707 = 864
λ 0, 0263W/m.K
hL = N uL = 864 = 22, 7 W/m2 . K.
L 1m
Donc,
1, 63 × 10−4 kg/s (2, 465 × 106 J/kg)
Ts = 300 K − = 282, 2 K.
22, 7 W/m2 . K


Adil Ridha Département de Mathématiques et Mécanique