Cours Convection M2 PDF
Cours Convection M2 PDF
Cours Convection M2 PDF
Master 2 GdP
Ph. Marty
2012-13
1 Introduction 3
4 Convection naturelle 20
4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
4.2 Convection naturelle sur plaque plane verticale chauffée à flux constant . . . . . . . . . . . . . . . . . . . . . 22
4.3 Convection naturelle entre plaques verticales parallèles (cheminée) . . . . . . . . . . . . . . . . . . . . . . . . 23
4.3.1 Condition d’existence de ce régime . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
4.3.2 Equations du régime établi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
4.3.3 Cas de deux plaques de même température . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
4.4 Autres géométries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
4.4.1 Paroi plane inclinée par rapport à la verticale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
4.4.2 Cylindre vertical . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
4.4.3 Cylindre horizontal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
4.4.4 Sphère . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
4.4.5 Plaques horizontales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
1
Nomenclature
T température (K)
R résistance thermique (K.W −1 )
C Chaleur massique (J.Kg −1 .K −1 )
q sources de chaleur volumiques (W.m−3 )
Q chaleur, energie (J)
h coefficient d’échange convectif (W.m−2 )
g gravité (m.s−2 )
P pression (P a)
t temps (s)
~
V vitesse du fluide (m/s)
Nu Nombre de Nusselt
Gr Nombre de Grashof
Ra Nombre de Rayleigh
Pr Nombre de Prandtl
Pe Nombre de Péclet
Ri Nombre de Richardson
ϕ
~ densité de flux de chaleur (W.m−2 )
Φ Flux de chaleur (W )
λ conductivité thermique (W.m−1 .K −1 )
2
∇ = ∆ operateur laplacien
ρ masse volumique (kg.m−3 )
λ
α = ρC diffusivite thermique (m2 .s−1 )
β dilatabilité volumique (K −1 )
2
Chapter 1
Introduction
L’objet de chapitre est de rappeler comment l’écriture sous forme adimensionnelle de l’équation de transport de la chaleur
fait apparatre les nombres sans dimensions caractéristiques de la convection dans les fluides. On y ra pelle aussi ce qui
caractérise les différents régimes de convection : forcée, naturelle et mixte. Ce chapitre se termine par la liste des principaux
nombres sans dimension utiles en convection thermique.
Lorsque le champ de vitesse est imposé, le champ de température est totalement dépendant de celui-ci. Cette situation
est celle de la convection forcée dans laquelle la vitesse est donc insensible aux variations de température dans le fluide. La
température T obéit alors à une équation de transport:
DT
ρCP = λ∇2 T + ΦS
Dt
où:
DT ∂T ~ .∇)T
~
= + (V
Dt ∂t
Le dernier terme représente un éventuel terme source pouvant résulter d’une réaction chimique, de l’effet Joule si le fluide
est métallique ou de la contribution des forces visqueuses si l’écoulement est supersonique (ΦS est en W.m−3 ).
Lorsque le champ de vitesse est crée par le champ de température, on dit que la convection est naturelle, et la vitesse
obéit à:
- En convection forcée:
Si la distribution de vitesse n’est pas connue, on peut la chercher par résolution des équations de Navier-Stokes. En choisis-
sant les variables adimensionnelles suivantes:
T − Tref + U x p t
T+ = ;U = ; x+ = ; p+ = 2 ; t+ =
∆T Uref Lref ρUref Lref /Uref
l’équation de la vitesse devient:
DV~+ ~ + + 1 ∇2 V~+
+
= −∇p
Dt Re
U L
où le nombre de Reynolds est tel que Re = refν ref .
La solution obtenue peut ensuite être ensuite introduite dans l’équation de T qui, en écriture adimensionnelle, devient:
ρCP ∆T DT + λ∆T
. +
= 2 ∇2 T + + Φ+
S
Lref /Uref Dt Lref
soit encore:
DT + 1 2 +
+
= ∇ T + Φ+
S
Dt Pe
3
Uref Lref
où P e désigne le nombre de Péclet: P e = α avec α = λ
ρCP diffusivité thermique du fluide. Le terme Φ+
S est tel que
Φ L
Φ+ S ref
S = Uref ∆T ρCP .
- En convection naturelle:
L’utilisation des mêmes grandeurs adimensionnelles transforme l’équation de Navier-Stokes ainsi:
Si l’écoulement est rapide ( Re, Gr >> 1 ), on peut postuler un équilibre entre les forces d’inertie et la poussée d’Archimède
soit:
ρU 2 p
≈ ρgβ∆T → Uref = gβ∆T Lref
Lref
de sorte que Navier-Stokes devient:
DT + 1
= ∇2 T + + Φ+
S
Dt+ Gr1/2 P r
où le nombre de Prandtl est défini par P r = αν .
Si l’écoulement est lent (faible Grashof et Reynolds), un choix judicieux de l’échelle de vitesse consiste à équilibrer les
forces visqueuses et d’Archimède ce qui donne µ LU2 ≈ ρgβ∆T . L’échelle de vitesse est alors:
gβ∆T L2
Uref =
ν
- Convection mixte:
Lorsque de la convection naturelle se superpose à de la convection forcée, la question se pose de savoir si un des deux champs
de vitesse peut être négligé ou si les deux doivent être pris en considération.
On
p peut par exemple estimer le rapport des deux vitesses attendues pour chacun des modes de convection pris isolément,
soit gβ∆T Lref pour la convection naturelle et U0 pour la convection forcée. On forme ainsi le nombre de Richardson:
gβ∆T Lref
Ri =
U02
Si Ri >> 1, alors la convection naturelle domine alors que si Ri << 1, c’est la convection forcée qui prévaut.
- Le nombre d’Eckert: pour un écoulement de gaz à grande vitesse (nombre de Mach > 0.3), la puissance des forces
U02
visqueuses génère de la chaleur, ce qui se traduit par un terme source dans l’équation de T . Le nombre d’Eckert Ec = CP ∆T
mesure l’importance de cet effet: si Ec << 1, la contribution des forces visqueuses à l’échauffement du fluide peut être
négligée.
Pour l’ensemble des problèmes convectifs, les échanges de chaleur en paroi se mesurent à l’aide du nombre de Nusselt:
ϕreel
Nu =
λ∆T /Lref
où le numérateur désigne le flux de chaleur qui passe effectivement à travers la paroi et le dénominateur le flux qui circulerait
si seule la conduction agissait.
En convection forcée, le nombre de Nusselt est de la forme: N u = f (Re, P r).
En convection naturelle, il s’écrit: N u = f (Gr, P r) ou encore N u = f (Ra, P r)
4
hl
Nombre de Nusselt Nu = λ h : coefficient de convection l : N u traduit la qualité de
longueur caractéristique λ : con- l’échange thermique : une aug-
ductivité thermique du fluide mentation de ce nombre traduit
une contribution importante de
l’écoulement sur l’échange de
chaleur avec la paroi
ν
Nombre de Prandtl Pr = a ν : viscosité cinématique a : dif- P r compare l’aptitude du fluide
fusivité thermique du fluide à diffuser la quantité de mouve-
ment par le biais de sa viscosité à
son aptitude à diffuser la chaleur
par le biais de sa diffusivité ther-
mique
Ud
Nombre de Reynolds Re = ν U vitesse moyenne de
l’écoulement, d dimension
caractéristique et ν viscosité
cinématique du fluide
Ud
Nombre de Péclet Pe = α U vitesse moyenne de une valeur élevée de P e traduit
l’écoulement, d dimension une distorsion importante du
caractéristique et α = ρCλP champ de température due à
diffusivité thermique du fluide l’écoulement par rapport à ce
qu’il serait si seule la diffusion
était présente
gβ∆T l3
Nombre de Grashof Gr = ν2 β coefficient de dilatabilté une augmentation de Gr traduit
du fluide , l dimension car- une augmentation de l’intensité
actéristique, g gravité et ν de la convection naturelle
viscosité cinématique du fluide
gβ∆T l3
Nombre de Rayleigh Ra = νa ν viscosité cinématique , a diffu- Ra = Gr.P r pour de l’air ou
sivité thermique du fluide des fluides de nombre de Prandtl
proche de l’unité, Ra et Gr sont
très proches
gβ∆T l
Nombre de Richardson Ri = U2 en convection mixte, Ri >> 1
traduit l’importance de la con-
vection naturelle par rapport à la
convection forcée
5
Chapter 2
Ce chapitre présente quelques résultats relatifs à la convection interne, c’est à dire en conduite. On traitera d’abord le cas
des écoulements laminaires, puis turbulents.
• le champ de vitesse est un écoulement de Poiseuille supposé non pertubé par la convection naturelle. Il est en régime
établi.
• le champ de température est supposé aussi établi, cela suppose qu’on soit suffisament éloigné de la zone d’entrée où
prend naissance une couche limite thermique
La longueur,L, d’établissement sera estimée plus précisément par la suite. Elle est telle que :
L UD
≈ 0, 05.ReD avecReD =
D ν
Un simple bilan thermique sur une tranche dx donne:
dP = ϕ2πRdx = ṁCP dT
d’où l’expression de l’accroissement de température le long de la conduite:
dT 2πRϕ
= = Cste
dx ρQCP
La température croit donc linéairement le long du tube.
Par ailleurs, on peut très simplement montrer que la solution de l’équation de la chaleur donne un profil de T qui est un
polynome de degré 4. On en déduit le nombre de Nusselt :
48
Nu = = 4, 36
11
Toutefois on peut déduire sans calculs l’ordre de grandeur de la différence de température entre la paroi et l’axe en disant
que cet écart doit être tel qu’il permette le passage du flux ϕ imposé par la paroi:
TP − Taxe ϕR
ϕ≈λ → TP − Taxe ≈
R λ
Ceci est confirmé par la figure 2.3 qui montre les profils de T dans les directions radiales et axiales.
6
r
U
x
dx
ϕP
ϕP
Début du chauffage
Figure 2.2: Etablissement du régime thermique à l’entrée d’une conduite circulaire chauffée sous flux constant
θ = TP - T ϕ/ϕ0 T
TP - Ta
Θ=1
paroi
Région de Région de
développement de régime
la couche limite thermique établi
thermique TP (x)
Tm (x)
η = r/R
0 x0 x
0,617 0,816
Figure 2.3: Conduite circulaire chauffée sous flux constant: profils de température et de flux radial de la chaleur.
7
Convection laminaire en canal rectangulaire chauffé à flux constant La corrélation de Shah et London (1978) donne le
nombre de Nusselt en fonction du rapport d’aspect γ du canal pour un chauffage à flux constant sur les 4 faces du canal :
ou γ est le rapport du plus petit coté au plus grand et est donc inférieur à 1. Pour un canal carré (γ = 1) on obtient ainsi :
N u = 3, 64. Le nombre de Nusselt est défini ici comme : N u = hD λ .
h
Figure 2.4: Nombre de Nusselt et facteur de friction f (f tel que ∆P = f DLh 12 ρV 2 ) pour des écoulements laminaires
pleinement développés en conduites de section quelconque (d’après Incropera and De Witt). La première colonne représente
le rapport d’aspect de la conduite, la seconde le nombre de Nusselt lorsqu’un flux uniforme est imposé à la paroi, la troisième
le nombre de Nusselt lorsque la paroi est main tenue à température uniforme.
8
Ecoulement turbulent en tube circulaire lisse infiniment long :
Pour L/D > 60 ; 0, 7 < P r < 100 ; 104 < ReD < 1, 2.10
5
Formule de COLBURN : N uD = 0, 023.P r1/3 .Re0,8
D N uD = ϕ
λ∆T /D
Les propriétés du fluide sont évaluées à ( Tparoi + Tmélange)/2
On peut aussi utiliser la corrélation tirée d’expériences
de DITTUS et BOELTER :
N uD = 0, 0243.Re0,8.P r0,4 si T p > T m ϕ
N u = λ ∆T
N uD = 0, 0265.Re0,8.P r0,3 si T p < T m D
L
Dans ces deux dernières formules, N uD,∞ est celui d’un tube infiniment long ( D > 60).
La modélisation numérique permet d’avoir des informations plus précises sur le profil transverse de vitesse et de température
dans la conduite mais ce type de renseignements n’est souvent nécessaire que lors du développement d’un nouveau produit
et rarement lors du dimensionnement d’un projet industriel.
Pour des conduites de section quelconque , on peut, à défaut de données fiables, reprendre les formules valables en tubes
circulaires et remplacer D par le diamètre hydraulique DH de la conduite qu’on souhaite étudier.
9
Chapter 3
Contrairement au chapitre précédent qui décrit les échanges thermiques à l’intérieur d’une conduite on s’intéresse ici aux
échanges lorsque le fluide s’écoule à l’extérieur d’un corps. On traitera ainsi le cas des écoulements sur plaque plane, laminaires
ou turbulents, ainsi que quelques exemples d’écoulements autour d’obstacles.
∂u ∂u 1 ∂p ∂2u
u +v =− +ν 2
∂x ∂y ρ ∂x ∂y
et :
∂p
=0
∂y
∂ ∂
Ces équations ont été établies en supposant que ∂x << ∂y et que la pression ne varie pas suivant la direction Oy.
Un raisonnement analogue fournit une équation simplifiée de la température:
∂T ∂T ∂T ∂T ∂2T
>> →u + v. =α 2
∂y ∂x ∂x ∂y ∂y
10
y U0
δ (x)
u(x,y)
AAAAAAAAAAAAAA l
x
Figure 3.1: Développement d’une couche limite hydrodynamique sur une plaque plane
La résolution combinée de ces 3 équations (qdm suivant x, pression et température) permet le traitement d’une couche
limite thermique sur paroi en présence d’un écoulement. Cela reste encore complexe et l’objet de ce cours est plutt de discuter
les ordres de grandeur régissant les échanges de chaleur et le nombre de Nusselt associé. Il est pour cela utile de se rappeler
qu’une couche limite résulte de la compétition entre deux effets qui sont:
- la diffusion (de quantité de mouvement ou de chaleur selon qu’on s’intéresse à une couche limite cinématique ou thermique).
Celle-ci est dominante dans la direction y normale à la paroi qui présente les plus forts gradients.
- la convection , c.a.d. le transport, par l’écoulement lointain qui a lieu principalement suivant Ox. √
Pour la couche cinématique, pendant un temps t, elle aura diffusé transversalement d’une distance δ = νt. Pendant ce
temps, la convection l’aura déplacée de x = U0 t. L’élimination de t dans ces deux relations fournit:
1/2
δ ν
= = Rex−1/2
x U0 x
L’expression précédente fait apparatre le nombre de Reynolds local Rex = Uν0 x basé sur l’abscisse le long de la plaque :
il varie donc e tout point de la plaque depuis la valeur 0 en x = 0 jusqu’à sa valeur maximale en x = l où il vaut Rel = Uν0 l .
Comme il a été vu en cours de Mécanique des Fluides de M 1 , cette dernière valeur doit être inférieure à 105 environ pour
1/2
νx √
que l’écoulement reste laminaire. L’expression de δ peut aussi s’écrire δ = U0
qui montre l’évolution en x de la couche
limite.
Le traitement par ordres de grandeur de problèmes thermiques nécessite de considérer 2 cas extrêmes caractérisés par P r << 1
ou P r >> 1 (voir Figure 3.3).
δ << δt
11
On peut alors estimer le flux de chaleur extrait de la paroi:
TP − T0
ϕp ≈ +λ.
δt
et en le normalisant par : λ TP x−T0 , on obtient:
x
N ux =≈
δt
soit:
N ux ≈ Re1/2
x .P r
1/2
avec:
U0 .x
Rex =
ν
δt << δ
La couche limite thermique est donc incluse dans une région de plus basse vitesse que U0 , variant de 0 à la paroi jusqu’à une
valeur que l’on peut estimer à δδt .U0 si on linéarise le profil de vitesse de U . Pendant un temps t, on écrit donc:
√ δt
δt = αt et x= .U0 t
δ
d’où:
x1/2 ν 1/6 α1/3
δt ≈ 1/2
U0
On peut former le rapport des épaisseurs des 2 couches:
δ p
≈ 3 ν/α = P r1/3
δt
Le nombre de Nusselt s’écrit maintenant:
x
N ux = ≈ Re1/2
x .P r
1/3
δt
Seul l’exposant du nombre de Prandtl a changé par rapport au cas P r << 1. La figure 3.4 tirée de Eckert and Drake, fournit
les résultats de l’analyse exacte (traitement complet des équations de la couche limite que nous avons évoqué dans le premier
pararaphe).
Uinfini
T0
AAAAAAAAAAAAAA
T( y) δ t (x)
TP
AAAAAAAAAAAAAA
AAAAAAAAAAAAAA
x
12
y y
U0
U0
T0
u(y)
δt δ
AAAAAAAAAAA
ut
T0
AAAAAAAAAAAAAAAAAAAAAA δ δt
AAAAAAAAAAAAAAAAAAAAAA
TP TP
Figure 3.3: Profils de température et de vitesse en fonction du nombre de Prandtl. a)- P r << 1 la couche limite se trouve
presque entièrement située dans une zone où la vitese vaut U0 b)- P r >> 1 la couche limite thermique se trouve placée dans
une région où la vitesse évolue de 0 à la paroi jusqu’à une valeur Ut en y = δt .
Figure 3.4: Couche limite laminaire sur plaque plane: Nusselt number vs Prandtl number (Eckert & Drake - Analysis of heat
& mass transfer)
.
13
3.1.4 Expressions exactes - Flux total échangé (en laminaire)
Le calcul complet tenant compte du profil de vitesse réel dans la couche limite est en fait possible pour certains écoulements
et notamment celui sur plaque plane. On trouve :
1/2
• Pour P r < 0.1 : N ux = 0, 565.Rex . P r1/2
(Les propriétés du fluide sont à prendre à la température du film T p+T 2
∞
.)
RL
ϕ.dx R L N ux .λ ∆T
x .dx
RL
Nombre de Nusselt global : N uL = λφ/L ∆T =
0
λ∆T = 0 λ∆T = 0 N ux dx
x
L
Uo .x 1/2 dx
RL 1/2 1/2
= 1, 13.ReL P r1/2
soit : N uL = 0
0, 565.P r ν .x avec
U0 .L
ReL = ν
1/2
• Pour P r > 0.1 : N ux = 0, 332.Rex . P r1/3
1/2
de même : N uL = 0, 664.ReL . P r1/3
Le tableau suivant synthétise les résultats :
1/2
P r << 1 N uL = 1, 13.ReL P r1/2
1/2
P r >> 1 N uL = 0, 664.ReL P r1/3
Nombre de Nusselt pour un écoulement laminaire d’un fluide sur une plaque plane
• u+ = 5.5 + κ1 Lny + pour y + ∈ [100; 1000] où κ = 0.41 est la constante de Karman.
On montre donc que le profil de T s’écrit:
ρCP U ∗ (TP −T )
• dans la zone visqueuse (y + ∈ [0; 10]): T + = P r.y + avec T + = ϕP
où U ∗ est la vitesse de frottement classiquement définie, TP et ϕP la température et le flux à la paroi. T désigne la
température de mélange de l’écoulement.
• dans la zone Log on trouve: T + = 2.195 Lny + + 13.2 P r − 5.66
On voit donc que l’allure générale des profils de T et de U sont très semblables, que ce soit dans la zone visqueuse ou
dans la zone Log.
14
• Si Re > 40 , il existe une allée de tourbillons de Von Karman à l’aval du cylindre
• Si Re ր encore, la couche limite qui se développe à partir du point d’arrêt à θ = 0 devient turbulente et décolle à θd
Re < 3.105 → θd = 80◦
Re > 3.105 → θd = 120◦ et la couche limite est turbulente même entre θ = 0 et θ = θd
Figure 3.5: Ecoulement autour d’un cylindre. Ref.: Favre-Marinet, cours de l’ENSHMG
• Si Re petit, Nu est petit (de quelques unités), surtout à l’aval où les recirculations sont peu actives.
• Pour Re > 426 000, Nu augmente beaucoup (plusieurs centaines) avec un max vers θd = 120◦
15
Figure 3.6: Distribution azimutale du coefficient d’échange autour d’un cylindre dans un écoulement d’air- a) Re petit, b)
Re intermédiaire. Ref.: Eckert and Drake. Attention: l’écoulement va de droite à gauche sur la Fig.a et de bas en haut sur
la Fig.b.
Figure 3.7: Distribution azimutale du coefficient d’échange autour d’un cylindre dans un écoulement d’air- Re >> 1. Ref.:
Eckert and Drake (l’écoulement va de gauche à droite).
16
Figure 3.8: Convection forcée autour d’un cylindre dans l’air. Nombre de Nusselt moyen en fonction de Re d’après Hilpert
(1933).
ϕ̄ φ V∞ .d
avec N u = (T1 −T2 ) et ϕ̄ = S et Re = ν
λ d
Application à la mesure de vitesse par fil chaud : La sonde est un fil de platine de diamètre très petit d ≈ 5µm. Un
asservissement électronique le maintient à température constante par apport calorifique d’un courant engendrant un effet
Joule. La mesure du courant i est fonction du flux φ échangé par le fil.
?
fil i V0
V0 d
Calcul de Re dans l’air à V0 = 100 m/s : la viscosité cinématique ν ≈ 10−5 m2 /s : on en déduit Re = ν ≈ 50 (on est donc
√
en régime laminaire) et le nombre de Nusselt dépend de Re−1/2 d’où Nu ≈ V0
φ/πdL φ
Par ailleurs Nu = λ∆T /d = πL.λ∆T
√
Donc le flux φ est proportionnel à V0 : la réponse de la sonde sera non- linéaire de sorte que la sonde nécessitera un
étalonnage préalable réalisé dans un petite soufflerie par exemple où la vitesse est très bien connue.
φ = hS∆T
17
φ hS∆T hd
N ud = = =
λ ∆T
d .S ∆T.S
λ d λ
On pose N ud = C.Rem d .P r
0,35
avec Red = U∞ .d/ν
La figure 3.9 donne les valeurs de C et m pour diverses géométries.
Figure 3.9: Coefficients d’échanges en fonction de la géométrie en convection forcée externe. (d’après Taine et Petit). La
flèche indique le sens de l’écoulement.
18
3.3.3 Obstacle sphérique
d
'$
~∞
U
- TP
&%
T∞
19
Chapter 4
Convection naturelle
4.1 Introduction
Rappel de la définition : il s’agit de la description du mouvement d’un fluide engendré par les forces d’Archimède dues aux
variations de la masse volumique avec la température. Il y a donc couplage de la dynamique et de la thermique. Le champ
de vitesse transporte la chaleur et, du fait de la dépendance de ρ avec T , agit sur la distribution de masse volumique ; en
retour, les modifications de masse volumique créent du mouvement par le biais de la poussée d’Archimède.
La convection naturelle est un phénomène important dans de nombreux procédés industriels :
• refroidissement d’appareils électriques
• collecteurs solaires
• chauffage des locaux
• centrales nucléaires
Les équations du mouvement en présence de convection naturelle ont été présentées dans le chapitre 1 : leur forme
adimensionnelle a fait apparatre les nombres de Grashof et de Rayleigh :
g β ∆ T . l3 g β ∆ T . l3
Gr = et Ra = Gr.P r =
ν2 αν
où β désigne la dilatabilité thermique du fluide :
1 ∂ρ
β=− ( )P
ρ ∂T
1
Pour un gaz parfait on montre facilement que β = T. Par exemple, pour de l’air à 20C , assimilé à un gaz parfait, on
1
obtient : β = 293 ≈ 3.10−3 .
Avant de présenter les résultats d’échanges thermiques par convection naturelle le long d’une paroi, il faut préciser le
critère qui caractérise la transition laminaire/turbulent d’un écoulement ascendant. La figure 4.1 illustre cette question et
montre que lorsque le nombre de Grashof local atteind 109 environ, l’écoulement devient turbulent.
Régime laminaire Avant de donner les expressions exactes du nombre de Nusselt il est intéressant, comme nous l’avions
fait en convection forcée, de comprendre l’effet du nombre de Prandtl.
• Lorsque P r >> 1 cela correspond donc à un fluide ayant plus de facilité à diffuser la quantité de mouvement que la
chaleur. Il existe donc une couche limite thermique d’épaisseur δt (cf. Figure 4.2) dans laquelle la température est plus
chaude qu’à l’infini où elle vaut T0 . Cette couche limite thermique est donc le siège de forces d’Archimède qui élèvent
le fluide et il y a création d’un jet pariétal d’épaisseur δc supérieur à δt .
• Lorsque P r << 1 l’on s’attendrait à ce que la couche limite thermique soit beaucoup plus épaisse que la couche
cinématique. Cela signifierait alors qu’il existerait une zone de fluide immobile dans laquelle la température est pourtant
supérieure à T0 . Cela est évidement impossible car la poussée d’Archimède contrarierait cet équilibre instable : on
observe en fait que les couches limites cinématiques et thermiques ont la même épaisseur (cf. Figure 4.3).
20
AA
AA
AA
AA
AA
AA
AA Gr L = 10 9
AA
AA
AA
Figure 4.1: Convection naturelle le long d’une plaque verticale: transition à la turbulence
AAA
AAA
z
AAA
AAA T1
u
T0
g
AAA
AAA
Τ
AAA
L
y
AAA
δt δc
AAA
AAA
AAA
Figure 4.2: Convection naturelle le long d’une plaque plane verticale pour P r >> 1.
21
AAA
AAA
z
AAA
AAA
u
T1
T0
g
AAA
AAA
Τ
AAA
L
δc y
AAA
δ1
δt
AAA
AAA
AAA
Figure 4.3: Convection naturelle le long d’une plaque plane verticale pour Pr << 1.
En résolvant le problème de façon analytique, pour un nombre de Prandtl quelconque et pour une paroi de température
imposée, on trouve une formule unique :
1/2 1/4
hz z Pr
N uz = λ = 0, 508. (0.95+P r)1/4
.Grz
1/2 1/4
hL Pr
N uL = λ = 0, 677. (0.95+P r)1/4
.GrL
Régime Turbulent : les formules ci dessus ne sont valables qu’en régime laminaire c’est à dire quand Gr < 109
Pour Gr > 109 , le régime est turbulent et MAC ADAMS propose la corrélation :
4.2 Convection naturelle sur plaque plane verticale chauffée à flux constant
Cette situation se rencontre lorsque la paroi contient par exemple une source de chaleur électrique dont la puissance est
évacuée par la paroi : le flux à la paroi est alors imposé. Les résultats suivants seront alors utilisés:
N uz = 0.6(Grz∗ P r)1/5
Pour 2.1013 < Grz∗ P r < 1016 (régime turbulent):
N uz = 0.568(Grz∗ P r)0.22
Le nombre de Grashof local modifié s’écrit dans ce cas :
gβz 3 ϕP z gβϕP z 4
Grz∗ = =
ν2 λ ν2λ
Dans cette expression, ϕP désigne le flux à la paroi.
22
4.3 Convection naturelle entre plaques verticales parallèles (cheminée)
• Si la longueur L>>2e, on aura un régime hydrodynamique établi.
• On cherche alors un champ de vitesse de la forme :
Ux = u(y)
Uy = 0
AA AAA
x
TP
AA
AA AAA
AAA
TP
AA
AA AAA
AAA y
AA
AA AAA
AAA
T0 AA 2e
AAA T0
∂T
L’équation de la chaleur se simplifie puisqu’ ici on a ∂x =0
2
→ α ∂∂yT2 = 0 → ∂T
∂y = cste = 0
→ T = Tp partout ; le fluide est isotherme. Il subit la poussée d’Archimède par rapport au fluide externe à TO . Cette
poussée sera contrée par la viscosité.
Pour la poussée : p* (x = 0)= p* (x = L) car c’est le fluide externe qui communique la pression ambiante hydrostatique
aux 2 extrémités de la cheminée. On a donc :
∂p∗
= 0 d’où N-S se simplifie :
∂x
23
h i
∂2U gβ(Tp −TO ) gβe2 y 2
∂y 2 =− ν → u= 2ν (Tp − TO ) 1 − e
Ce flux a été nécessaire pour faire passer le fluide de T0 en entrée à Tp en sortie. Il a en fait été donné dans la zone
d’établissement (dans les couches limites), seules régions où le gradient de température à la paroi n’est pas nul. Cela explique
l’absence de L dans l’expression de φ .
Nombre de Nusselt : N uL = φ/2L.l Tp −T0 (facteur 2 car il existe 2 plaques)
λ. L
La présence de L au dénominateur indique que c’est un flux conductif axial qui sert de référence; c’est en effet selon x
qu’il existe un gradient, pas selon y.
On trouve ainsi :
Ra2e hL
N uL = 24 = λ
gβ∆T (2e)3
avec Ra2e = αν (∆T = Tp − T0 )
AAAAAA
AAAAAA
AAAAAA
AAAAAA
g θ
TP AAAAAA
AAAAAA
AAAAAA
T0
Il suffit de remplacer g par g cos θ dans les formules trouvées en paroi verticale (voir Figure 4.5).
−1/4 L
δ = LRaL << d soit : 1/4
<< 1
dRaL
Introduisons le nombre de Rayleigh, Rad , basé sur le diamètre du cylindre :
d3
Rad = RaL . , d’où :
L3
24
δ
d
L
L d3/4
1/4
d Ra1/4 L3/4 << 1 → Rad . Ld >> 1
d
Si cette hypothèse est satisfaite on pourra utiliser les formules obtenues pour plaque plane :
1/4
Pour l’air : N uL = 0, 547.RaL , en laminaire , Raǫ(104 , 109 )
1/3
(P r = 0.7) N uL = 0, 13.RaL , en turbulent , Raǫ(109 , 1013 )
4.4.4 Sphère
1/4
On prend N ud = 2 + 0, 45Grd .P r1/3 , Grd ∈ [1, 106 ].
25
Figure 4.7: Coefficient d’échange en convection naturelle: parois verticales ou horizontales N u = C.Ran . Source: Eyglunent,
Thermique théorique et pratique, Hermès
26
Formulaire
Continuité en cartésien :
∂ρ ∂ρvx ∂ρvy ∂ρvz
+ + + =0
∂t ∂x ∂y ∂z
Continuité en cylindrique:
∂ρ
∂t+ 1r ∂(ρrv
∂r
r)
+ r1 ∂(ρv
∂θ
θ)
+ ∂(ρv
∂z
z)
=0
Equations de Navier Stokes en coordonnées cartésiennes (x, y, z)
2
∂ 2 vx ∂ 2 vx
∂vx ∂vx ∂vx ∂vx ∂p ∂ vx
ρ + vx + vy + vz =− + ρgx + µ + +
∂t ∂x ∂y ∂z ∂x ∂x2 ∂y 2 ∂z 2
2
∂ 2 vy ∂ 2 vy
∂vy ∂vy ∂vy ∂vy ∂p ∂ vy
ρ + vx + vy + vz =− + ρgy + µ + +
∂t ∂x ∂y ∂z ∂y ∂x2 ∂y 2 ∂z 2
2
∂ 2 vz ∂ 2 vz
∂vz ∂vz ∂vz ∂vz ∂p ∂ vz
ρ + vx + vy + vz =− + ρgz + µ + +
∂t ∂x ∂y ∂z ∂z ∂x2 ∂y 2 ∂z 2
Equations de Navier Stokes en coordonnées cylindriques (r, θ, z)
v2 1 ∂ 2 vr ∂ 2 vr
∂vr ∂vr vθ ∂vr ∂vr ∂p 1 ∂ ∂vr vr 2 ∂vθ
ρ + vr + − θ + vz =− + ρgr + µ (r )− 2 + 2 − +
∂t ∂r r ∂θ r ∂z ∂r r ∂r ∂r r r ∂θ2 r2 ∂θ ∂z 2
1 ∂ 2 vθ ∂ 2 vθ
∂vθ ∂vθ vθ ∂vθ vr vθ ∂vθ 1 ∂p 1 ∂ ∂vθ vθ 2 ∂vr
ρ + vr + + + vz =− + ρgθ + µ (r )− 2 + 2 − +
∂t ∂r r ∂θ r ∂z r ∂θ r ∂r ∂r r r ∂θ2 r2 ∂θ ∂z 2
1 ∂ 2 vz ∂ 2 vz
∂vz ∂vz vθ ∂vz ∂vz ∂p 1 ∂ ∂vz
ρ + vr + + vz =− + ρgz + µ (r )+ 2 2
+
∂t ∂r r ∂θ ∂z ∂z r ∂r ∂r r ∂θ ∂z 2
Contraintes visqueuses dans un fluide newtonien en coordonnées cartésiennes:
∂vx
τxx = 2µ
∂x
∂vy
τyy = 2µ
∂y
∂vz
τzz = 2µ
∂z
27
∂vx ∂vy
τxy = τyx = µ +
∂y ∂x
∂vy ∂vz
τyz = τzy = µ +
∂z ∂y
∂vz ∂vx
τzx = τxz =µ +
∂x ∂z
Contraintes visqueuses dans un fluide newtonien en coordonnées cylindriques:
∂vr
τrr = 2µ
∂r
∂vθ
τθθ = 2µ
∂θ
∂vz
τzz = 2µ
∂z
∂ vθ 1 ∂vr
τrθ = τθr = µ r ( ) +
∂r ∂r r ∂θ
∂vθ 1 ∂vz
τθz = τzθ = µ +
∂z r ∂θ
∂vz ∂vr
τzr = τrz = µ +
∂r ∂z
Equations de l’ énergie en coordonnées cartésiennes (avec ρ et λ constants):
2
∂2T ∂ 2T
∂T ∂T ∂T ∂T ∂ T
ρCP + vx + vy + vz =λ + + +q+D
∂t ∂x ∂y ∂z ∂x2 ∂y 2 ∂z 2
1 ∂2T ∂2T
∂T ∂T vθ ∂T ∂T 1 ∂ ∂T
ρCP + vr + + vz =λ (r )+ 2 2 + +q+D
∂t ∂r r ∂θ ∂z r ∂r ∂r r ∂θ ∂z 2
avec q source de chaleur interne (électrique, chimique, nucléaire, etc... en W/m3 ) et D puissance des forces visqueuses
(en W/m3 également mais souvent negligeables) telle que:
en cartésien :
" 2 2 2 # " 2 2 2 #
∂vx ∂vy ∂vz ∂vy ∂vx ∂vz ∂vy ∂vx ∂vz
D = 2µ + + +µ + + + + +
∂x ∂y ∂z ∂x ∂y ∂y ∂z ∂z ∂x
en cylindrique :
" 2 2 2 # 2 2 2
∂vr 1 ∂vθ vr ∂vz ∂ vθ 1 ∂vr 1 ∂vz ∂vθ ∂vr ∂vz
D = 2µ + + + +µ r + +µ + +µ +
∂r r ∂θ r ∂z ∂r r r ∂θ r ∂θ ∂z ∂z ∂r
28
Bibliographie
Plusieurs des figures de ce document sont tirées des ouvrages suivants qui pourront être consultés avec intérêt.
SOLLIEC C., VIAZZO S. : Cours d’Echanges Thermiques de l’Ecole des Mines de Nantes
29