Université Abdelmalek Esaadi
Faculté des Sciences de Tétouan
Master GEE/S2
2021-2022
Transfert de chaleur par
Convection
ZAKIA NGADI
[Link]@[Link]
II. Solution des problèmes de convection
3. Convection libre sur une plaque plane
La convection forcée suppose la présence d’un dispositif pour mettre en mouvement le
fluide (pompe, compresseur…)
En convection libre, le mouvement résulte de l’existence d’un gradient de température
dans le fluide.
En raison de la différence de températures, le champ de densité n'est plus uniforme. Alors,
la flottabilité induira un courant d'écoulement dû au champ gravitationnel et à la variation
de densité.
En général, un transfert de chaleur par convection libre
est généralement beaucoup plus petit qu'un transfert de
chaleur par convection forcée.
Si la température est uniforme, le fluide est au repos.
II. Solution des problèmes de convection
3. Convection libre sur une plaque plane
Considérons la plaque plane verticale présentée ci-dessous. La plaque a une température
uniforme Ts et elle est submergée dans un fluide infini à une température T∞. Le régime
est laminaire permanent et le fluide est considéré newtonien à propriétés constantes.
La dissipation et la génération de l’énergie sont négligeables.
Dans les problèmes de convection naturelle, la masse volumique est une
fonction de la température
La variation de masse volumique d'un fluide avec la température à pression
constante peut être exprimée en termes de coefficient de dilatation volumique
1
( )
T
•Les champs de vitesse et de température sont couplés en convection libre
•Les solutions analytiques nécessitent l'intégration simultanée des équations de
continuité, de quantité de mouvement et d'énergie.
II. Solution des problèmes de convection
3. Convection libre sur une plaque plane
Formulation du problème de convection libre
Considérons le bilan des forces sur un élément de volume de contrôle dans la couche
limite laminaire,
Equation de u v
Conservation 0
de masse x y
Equation de
u u 2u
Conservation u v g (T T ) 2
x y y
de mouvement (selon X)
Equation de
Conservation T T 2T 2
u v k 2 ou u v k 2
de l’énergie x y y x y y
II. Solution des problèmes de convection
3. Convection libre sur une plaque plane
Conditions aux limite de vitesse
Pour la vitesse la condition de non glissement est appliquée à la paroi et la condition de
l’écoulement libre en dehors de la couche limite
y0 u ( x, 0) 0, v( x, 0, ) 0, ( x, 0) 0
y u ( x, ) 0, v( x, ) 0, ( x, ) 0
x0 u (0, y ) 0 (0, y ) 0
II. Solution des problèmes de convection
3. Convection libre sur une plaque plane
Nombre de Grashof et de Rayleigh
L’Adimensionnalisation des équations gouvernantes en remplaçant les variables
dépendantes par des variables indépendantes, l’équation de moment fait apparaître le
nombre de Grashof (Gr) qui représente les effets de la convection naturelle.
g TL3
Gr
2
Un autre nombre adimensionnel important dans les problèmes de convection libre est le
nombre de Rayleigh (Ra)
g TL3
Ra Gr Pr
k
•L’approximation de la couche limite pour la convection libre est valable pour Ra 104
•La transition d’un écoulement laminaire à un écoulement turbulent pour une plaque plane
est considérée par le nombre de Rayleigh de transition donné par Rat 109
II. Solution des problèmes de convection
3. Convection libre sur une plaque plane
Approche de similitude
L’approche repose sur la méthode de similitude pour transformer les trois équations
partielles en deux équations différentielles ordinaires. Dans le cas de convection libre sur
une plaque plane la variable de changement η(x,y) est définie par
1/4
Grx dy
1/4
Grx y d
4 x 4 x
On introduisant la fonction courant Ψ tell que u et v
y x
en utilisant la solution de Blasius de la convection forcée comme guide on peut écrire
v dF ( )
u 2 Grx1/2
y y d
1/4
Gr x
4v x F ( )
4 (Grx )1/4 dF
v
41/4
x d
3 F
x y x
II. Solution des problèmes de convection
3. Convection libre sur une plaque plane
Alors l’équation de moment devient
2
u u 2u
g (T T ) 2
d 3F d 2F dF
u v 3 F ( ) 2 0
x y y d d d
3 2
Équation différentielle ordinaire d’ordre 3
Et l’équation de l’énergie devient
2 2
u v k 2 3Pr F ( ) 0
x y y 2
Équation différentielle ordinaire d’ordre 2
II. Solution des problèmes de convection
3. Convection libre sur une plaque plane
Solution de Pohlhausen
Les conditions aux limites pour la nouvelle équation en fonction de η sont
df (0) df ()
f (0) 0 ; 0
d d
Les conditions aux limites pour la nouvelle équation sont
(0) 1 ; () 0
II. Solution des problèmes de convection
3. Convection libre sur une plaque plane
Solution du problème
•Les équations obtenues sont résolues numériquement. Les solutions sont présentées
graphiquement pour différentes valeurs de Prandtl.
II. Solution des problèmes de convection
3. Convection libre sur une plaque plane
Solution de Pohlhausen
Une fois la distribution de la température est déterminée, on peut calculer
Coefficient de transfert de chaleur h
T ( x, 0) dT d (0)
dy d d y
T ( x, 0)
4 GrL d (0)
1/4
Grx d (0)
1/4
y h
h d
Ts T x 4 d 3 L 4
Nombre de Nusselt
d (0) d (0)
1/4 1/4
hx Gr hL 4 Gr
Nu x x Nu L L
4 d 3 4 d
II. Solution des problèmes de convection
3. Convection libre sur une plaque plane
Solution de Pohlhausen
h et Nu dépendent du gradient de température à la surface qui a son tour dépend du Pr
Nu x 0.6(Pr Rax )1 4 , Pr 0
Nu x 0.503(Pr Grx )1 4 , Pr
Remarques
•L'utilisation de la solution de Pohlhausen nécessite la détermination
des propriétés du fluide.
•Lors des calculs utilisant la solution de Pohlhausen, les propriétés sont
évaluées à la température du film.
Convection libre sur une plaque plane
Application 7
Une plaque carrée de 8cm×8cm est suspendue verticalement dans l'air. La plaque est
maintenue à une température de surface uniforme de 70°C. L'air ambiant est à 10°C.
déterminer les conditions d'écoulement et de transfert de chaleur à l'extrémité x=L.
Plus précisément,
a) Vitesse axiale u à y=0.3cm
b) Température de l’air T à y=0.3cm
c) Épaisseur de la couche limite visqueuse et thermique.
d) Coefficient de transfert de chaleur h et le nombre de Nusselt
e) Taux de transfert de chaleur et le taux total de transfert par la plaque
Propriété de l’air à 40°C
Pr=0.71 ; λ= 0.0271 W/(m.°C) ; =16.96x10-6 m2/s; β=0.0031936K-1
III. Corrélation usuelles en convection
1. introduction
Il existe de nombreuses situations où la détermination analytique du coefficient de
transfert de chaleur h est difficile à obtenir
Dans les cas complexes par exemple géométrie, propriétés variables, écoulement
turbulent, etc. le coefficient de transfert de chaleur est déterminé expérimentalement.
Les résultats expérimentaux sont généralement regroupés et présentés sous forme
d’équations sans dimension.
III. Corrélation usuelles en convection
2. Convection forcé externe
Plaque plane à température uniforme
Régime laminaire x xt
hx
Nux x 0.332(Re x )1/2 (Pr)1/3
hL
Nu L 0.664(Re L )1/2 (Pr)1/3 ; Re L 5 105
Les propriétés du fluide sont
Régime turbulent x xt calculées à Tf
hx x
Nux 0.0296(Re x )4/5 (Pr)1/3
hL
Nu L 0.037(Re L ) 4/5 (Pr)1/3 ;5 105 Re x 107 ; 0.6 Pr 60
Cas particuliers
Plaque plane est suffisamment longue pour que l'écoulement devient turbulent, mais pas assez longue
pour ne pas tenir compte de la région d'écoulement laminaire.
h
1
L 0
xt L
h( x)la min aire dx h( x)turbulent dx
xt
III. Corrélation usuelles en convection
2. Convection forcé externe
Ecoulement autour d’un cylindre ou d’une sphère
1/4
hD 0.4 3.5 ReD 7.6104
Nu D 2 0.4 Re D 0.06 Re D Pr ; 0.7 Pr 380
1/2 2/3
k s
III. Corrélation usuelles en convection
3. Convection forcé interne
Vitesse entièrement développée, température en développement : écoulement laminaire
0.0668( D / L) Re D Pr
Nu D 3.66 ; Re D 2300
1 0.04 ( D / L) Re D Pr
2/3
Vitesse et température entièrement développées : écoulement laminaire
0.14 Les propriétés du fluide sont
Nu D 1.86 ( D / L) Re D Pr
𝑇 +𝑇
1/3
calculées à 𝑇𝑚 = 𝑚𝑒 2 𝑚𝑠
s
III. Corrélation usuelles en convection
3. Convection forcé interne
Écoulement interne turbulent
( f / 8)(Re D 1000) Pr D P
Nu D
1 ( D / L ) 2/3
Avec f
1 12.7( f / 8) (Pr 1)
1/2 2/3
L u 2 / 2
Pour un tube lisse f (0.79ln ReD 1.64)2 ; 104 Re 106
Corrélation de Colburn
Nu D 0.023(Re D ) 4/5 (Pr)1/3 ; 2300 Re D 5 106 ; 0.5 Pr 2000; L / D 60
Les propriétés du fluide sont
𝑇 +𝑇
calculées à 𝑇𝑚 = 𝑚𝑒 2 𝑚𝑠
Pour écoulement dans des tubes non circulaires on utilise le
4𝐴𝑓
diamètre hydraulique au lieu du diamètre 𝐷𝑒 = 𝑃
III. Corrélation usuelles en convection
4. Convection libre externe
Plaque plane verticale à température uniforme
Quel que soit le régime 101 RaL 1012 ; 0 Pr
Les propriétés du fluide sont
1/6
2 calculées à Tf
0.387( Ra )
Nu D 0.825
1 (0.49 / Pr)
9/16 8/27
Plaque plane verticale à flux de chaleur uniforme
1/5
Pr 2
Nu x Grx ;104 Grx* Pr 109 ; 0 Pr
4 9 Pr 10 Pr
1/2
Ts est calculé par la loi de Newton et le
𝛽𝑔𝑞𝑠" 4
nombre de Grashof modifié 𝐺𝑟𝑥∗ = 𝑥
𝜆𝜈2
III. Corrélation usuelles en convection
4. Convection libre externe
Plaque plane inclinée à température uniforme
RaL 109 ; 0 60
2
1/6 g cos (Ts T ) 3
0.387( Ra )
Nu D 0.825 Rax x
1 (0.49 / Pr)
9/16 8/27
k
Les propriétés du fluide sont
calculées à Tf
III. Corrélation usuelles en convection
4. Convection libre externe
Plaque plane horizontale à température uniforme
Surface supérieure de la plaque chaude ou surface inférieure de la plaque froide
Nu L 0.54( RaL )1/4 ; 105 RaL 2 107
Nu L 0.14( RaL )1/3 ; 2 107 RaL 3 1010
Surface inférieur de la plaque chaude ou surface supérieure de la plaque froide
Nu L 0.27( RaL )1/4 ; 3 105 RaL 3 1010
Ces corrélations sont valables pour une variété de formes
As
de plaques ou la longueur caractéristique est Lc
P
III. Corrélation usuelles en convection
4. Convection libre externe
Écoulement autour d’une surface
III. Corrélation usuelles en convection
5. Convection libre interne
Cavité verticale de forme rectangulaire
0.29
Pr Pr
Nu e 0.18 Rae ; 1 e/ L 2, Rae 103
0.2 Pr 0.2 Pr
0.28 1/4
Pr e
Nu e 0.22 Rae ; 2 e/ L 10, Rae 1010
0.2 Pr L
0.3
0.012 e
Nu e 0.42 Rae Pr
1/4
; 10 e/ L 40, 10 4 Rae 107
L
g Tc T f e3
Le flux de chaleur dû à la circulation est
Rae Pr déterminé à partir de la loi de Newton
2
q = ℎ 𝑇𝑐 + 𝑇𝑓
III. Corrélation usuelles en convection
5. Convection libre interne
Cavité fermé horizontale de forme rectangulaire
Nu e 0.06 Rae1/3 Pr 0.074 ; 3 105 Rae 7 109
À des faibles nombres de Rayleigh, le fluide reste stagnant et le
transfert de chaleur à travers la cavité se fait par conduction
jusqu'à Raec 1708
III. Corrélation usuelles en convection
Remarque
Les équations de corrélation présentées sont pour des processus et des configurations limitées. Il
convient de souligner que le traitement ci-dessus est très simplifié. Il existe de nombreuses autres
équations de corrélation pour des sujets tels que la condensation, l'ébullition, l'écoulement à grande
vitesse, la dissipation, les métaux liquides, les géométries à ailettes, les géométries irrégulières, la
microgravité, les fluides non newtoniens, etc.
Notez bien
-Toutes les équations de corrélation ont des limites qui doivent être soigneusement notées avant de les
appliquer.
-Premièrement, la géométrie est un facteur évident. Ainsi, chaque équation est valable pour une
configuration spécifique.
-Deuxièmement, les limites de la plage de paramètres, tels que les nombres de Reynolds, Prandtl, Rayleight
et Grashof, pour lesquels une équation de corrélation est valide, sont déterminées par la disponibilité des
données et/ou la mesure dans laquelle une équation corrèle les données.
-Les équations de corrélation sont basées sur des données déterminées expérimentalement, elles ne
fournissent pas toujours des prédictions très précises de h.