TRANSFERTS THERMIQUES CONVECTIFS

Master 2 GdP

Ph. Marty
2012-13

Isothermes autour dun cylindre chauffe

en presence dun ecoulement dair `a Re = 1260
Ref.: Eckert and Drake.

ES

GENIE
DES PROCED
Master 2
Universite Joseph Fourier, Grenoble
version modifiee le 9 Juillet 2012
Philippe.Marty@hmg.inpg.fr

Contents
1 Introduction

2 Convection forc
ee interne
2.1 Convection forcee laminaire en conduite circulaire chauffee `a flux constant . . . . . . . . . . . . . . . . . . . .
2.2 Convection forcee laminaire dans des conduites de section quelconque . . . . . . . . . . . . . . . . . . . . . .
2.3 Convection forcee turbulente dans un tube quelconque . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

6
6
8
8

3 Convection forc
ee externe
3.1 Convection forcee laminaire sur plaque plane . . . . . . . . . . . . .
3.1.1 Rappel sur les couches limites hydrodynamique et thermique
3.1.2 Echanges thermiques sur plaque plane pour P r << 1 . . . . .
3.1.3 Echanges thermiques sur plaque plane pour P r >> 1 . . . . .
3.1.4 Expressions exactes - Flux total echange (en laminaire) . . .
3.2 Convection forcee turbulente sur plaque plane . . . . . . . . . . . . .
3.3 Ecoulements forces autour dobstacles . . . . . . . . . . . . . . . . .
3.3.1 Obstacles cylindriques . . . . . . . . . . . . . . . . . . . . . .
3.3.2 Obstacles non circulaires . . . . . . . . . . . . . . . . . . . . .
3.3.3 Obstacle spherique . . . . . . . . . . . . . . . . . . . . . . . .

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

10
10
10
11
12
14
14
14
14
17
19

4 Convection naturelle
4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.2 Convection naturelle sur plaque plane verticale chauffee `a flux constant
4.3 Convection naturelle entre plaques verticales parall`eles (cheminee) . . .
4.3.1 Condition dexistence de ce regime . . . . . . . . . . . . . . . . .
4.3.2 Equations du regime etabli . . . . . . . . . . . . . . . . . . . . .
4.3.3 Cas de deux plaques de meme temperature . . . . . . . . . . . .
4.4 Autres geometries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.4.1 Paroi plane inclinee par rapport `
a la verticale . . . . . . . . . . .
4.4.2 Cylindre vertical . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.4.3 Cylindre horizontal . . . . . . . . . . . . . . . . . . . . . . . . . .
4.4.4 Sph`ere . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.4.5 Plaques horizontales . . . . . . . . . . . . . . . . . . . . . . . . .

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.
.
.

20
20
22
23
23
23
23
24
24
24
25
25
25

.
.
.
.
.
.
.
.
.
.

Nomenclature
T
temperature (K)
R
resistance thermique (K.W 1 )
C
Chaleur massique (J.Kg 1 .K 1 )
q
sources de chaleur volumiques (W.m3 )
Q
chaleur, energie (J)
h
coefficient dechange convectif (W.m2 )
g
gravite (m.s2 )
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 Peclet
Ri
Nombre de Richardson

~
densite de flux de chaleur (W.m2 )

Flux de chaleur (W )

conductivite thermique (W.m1 .K 1 )

2
= operateur laplacien

masse volumique (kg.m3 )

diffusivite thermique (m2 .s1 )

= C

dilatabilite volumique (K 1 )

Chapter 1

Introduction
Lobjet de chapitre est de rappeler comment lecriture sous forme adimensionnelle de lequation de transport de la chaleur
fait apparatre les nombres sans dimensions caracteristiques de la convection dans les fluides. On y ra pelle aussi ce qui
caracterise les differents regimes de convection : forcee, 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 impose, le champ de temperature est totalement dependant de celui-ci. Cette situation
est celle de la convection forcee dans laquelle la vitesse est donc insensible aux variations de temperature dans le fluide. La
temperature T obeit alors `
a une equation de transport:
CP

DT
= 2 T + S
Dt

o`
u:

T
DT
~ .)T
~
=
+ (V
Dt
t
Le dernier terme represente un eventuel terme source pouvant resulter dune reaction chimique, de leffet Joule si le fluide
est metallique ou de la contribution des forces visqueuses si lecoulement est supersonique (S est en W.m3 ).
Lorsque le champ de vitesse est cree par le champ de temperature, on dit que la convection est naturelle, et la vitesse
obeit `
a:
DV~
~ + 2 V
~ g(T Tref )~g
= p
Dt
o`
u le dernier terme represente la poussee dArchm`ede par unite de volume fluide.

Ecriture adimensionnelle des equations

- En convection forcee:
Si la distribution de vitesse nest pas connue, on peut la chercher par resolution des equations de Navier-Stokes. En choisissant les variables adimensionnelles suivantes:
T+ =

T Tref +
x
p
U
t
; x+ =
; p+ =
;U =
; t+ =
2
T
Uref
Lref
Uref
Lref /Uref

lequation de la vitesse devient:

DV~+
~ + + 1 2 V~+
= p
+
Dt
Re
U

o`
u le nombre de Reynolds est tel que Re = ref ref .
La solution obtenue peut ensuite etre ensuite introduite dans lequation de T qui, en ecriture adimensionnelle, devient:
CP T DT +
T
.
= 2 2 T + + +
S
+
Lref /Uref Dt
Lref
soit encore:

1 2 +
DT +
=
T + +
S
+
Dt
Pe

o`
u P e designe le nombre de Peclet: P e =

Uref Lref

avec =

CP

diffusivite thermique du fluide. Le terme +

S est tel que

S ref
+
S = Uref T CP .
- En convection naturelle:
Lutilisation des memes grandeurs adimensionnelles transforme lequation de Navier-Stokes ainsi:

DV~+
~ + + 1 2 V~+ + gT L T +
= p
2
+
Dt
Re
Uref
Il reste toutefois `a determiner Uref qui nest pas imposee par un mecanisme externe mais par la convection naturelle
elle-meme.
Si lecoulement est rapide ( Re, Gr >> 1 ), on peut postuler un equilibre entre les forces dinertie et la poussee dArchim`ede
soit:
p
U 2
gT Uref = gT Lref
Lref

de sorte que Navier-Stokes devient:

DU~+
~ + + 1 2 U~+ T + g~+
= p
Dt+
Gr1/2
o`
u le nombre de Grashof est defini par: Gr =
Lequation de la temperature devient alors:

gT L3
.
2

DT +
1
2 T + + +
=
S
Dt+
Gr1/2 P r
o`
u le nombre de Prandtl est defini par P r = .
Si lecoulement est lent (faible Grashof et Reynolds), un choix judicieux de lechelle de vitesse consiste `a equilibrer les
forces visqueuses et dArchim`ede ce qui donne LU2 gT . Lechelle de vitesse est alors:
Uref =

gT L2

- Convection mixte:
Lorsque de la convection naturelle se superpose `
a de la convection forcee, la question se pose de savoir si un des deux champs
de vitesse peut etre neglige ou si les deux doivent etre pris en consideration.
On
p peut par exemple estimer le rapport des deux vitesses attendues pour chacun des modes de convection pris isolement,
soit gT Lref pour la convection naturelle et U0 pour la convection forcee. On forme ainsi le nombre de Richardson:
Ri =

gT Lref
U02

Si Ri >> 1, alors la convection naturelle domine alors que si Ri << 1, cest la convection forcee qui prevaut.
- Le nombre dEckert: pour un ecoulement de gaz `a grande vitesse (nombre de Mach > 0.3), la puissance des forces
U02
visqueuses gen`ere de la chaleur, ce qui se traduit par un terme source dans lequation de T . Le nombre dEckert Ec = CP T
mesure limportance de cet effet: si Ec << 1, la contribution des forces visqueuses `a lechauffement du fluide peut etre
negligee.

Pour lensemble des probl`emes convectifs, les echanges de chaleur en paroi se mesurent `a laide du nombre de Nusselt:
Nu =

reel
T /Lref

o`
u le numerateur designe le flux de chaleur qui passe effectivement `a travers la paroi et le denominateur le flux qui circulerait
si seule la conduction agissait.
En convection forcee, le nombre de Nusselt est de la forme: N u = f (Re, P r).
En convection naturelle, il secrit: N u = f (Gr, P r) ou encore N u = f (Ra, P r)

hl

Nombre de Nusselt

Nu =

h : coefficient de convection l :
longueur caracteristique : conductivite thermique du fluide

Nombre de Prandtl

Pr =

: viscosite cinematique a : diffusivite thermique du fluide

Nombre de Reynolds

Re =

Ud

Nombre de Peclet

Pe =

Ud

U
vitesse
moyenne
de
lecoulement,
d
dimension
caracteristique et viscosite
cinematique du fluide
U
vitesse
moyenne
de
lecoulement,
d
dimension
caracteristique et = CP
diffusivite thermique du fluide

Nombre de Grashof

Gr =

gT l3
2

Nombre de Rayleigh

Ra =

gT l3
a

Nombre de Richardson

Ri =

gT l
U2

coefficient de dilatabilte
du fluide , l dimension caracteristique, g gravite et
viscosite cinematique du fluide
viscosite cinematique , a diffusivite thermique du fluide

N u traduit la qualite de
lechange thermique : une augmentation de ce nombre traduit
une contribution importante de
lecoulement sur lechange de
chaleur avec la paroi
P r compare laptitude du fluide
`a diffuser la quantite de mouvement par le biais de sa viscosite `a
son aptitude `a diffuser la chaleur
par le biais de sa diffusivite thermique

une valeur elevee de P e traduit

une distorsion importante du
champ de temperature due `a
lecoulement par rapport `a ce
quil serait si seule la diffusion
etait presente
une augmentation de Gr traduit
une augmentation de lintensite
de la convection naturelle
Ra = Gr.P r pour de lair ou
des fluides de nombre de Prandtl
proche de lunite, Ra et Gr sont
tr`es proches
en convection mixte, Ri >> 1
traduit limportance de la convection naturelle par rapport `a la
convection forcee

Figure 1.1: Nombres sans dimension utiles en convection

Chapter 2

Convection forc
ee interne
Ce chapitre presente quelques resultats relatifs `
a la convection interne, cest `a dire en conduite. On traitera dabord le cas
des ecoulements laminaires, puis turbulents.

2.1

Convection forc
ee laminaire en conduite circulaire chauff
ee `
a flux constant

situation tr`es courante (echangeurs)

la conduite fournit un flux constant : p (Tparoi > Tfluide )
le champ de vitesse est un ecoulement de Poiseuille suppose non pertube par la convection naturelle. Il est en regime
etabli.
le champ de temperature est suppose aussi etabli, cela suppose quon soit suffisament eloigne de la zone dentree o`
u
prend naissance une couche limite thermique
La longueur,L, detablissement sera estimee plus precisement par la suite. Elle est telle que :
L
0, 05.ReD
D
Un simple bilan thermique sur une tranche dx donne:

avecReD =

UD

dP = 2Rdx = mC
P dT
do`
u lexpression de laccroissement de temperature le long de la conduite:
dT
2R
= Cste
=
dx
QCP
La temperature croit donc lineairement le long du tube.
Par ailleurs, on peut tr`es simplement montrer que la solution de lequation de la chaleur donne un profil de T qui est un
polynome de degre 4. On en deduit le nombre de Nusselt :
Nu =

48
= 4, 36
11

Toutefois on peut deduire sans calculs lordre de grandeur de la difference de temperature entre la paroi et laxe en disant
que cet ecart doit etre tel quil permette le passage du flux impose par la paroi:
TP Taxe
R
TP Taxe
R

Ceci est confirme par la figure 2.3 qui montre les profils de T dans les directions radiales et axiales.

U
x
dx

Figure 2.1: Conduite circulaire chauffee sous flux constant

Couche limite thermique

D

Ta (axe)
t (x)

x
rgime
thermique tabli

Tp

P
Dbut du chauffage

Figure 2.2: Etablissement du regime thermique `a lentree dune conduite circulaire chauffee sous flux constant

/0

= TP - T
TP - Ta
=1

paroi
Rgion de
dveloppement de
la couche limite
thermique

Rgion de
rgime
thermique tabli
TP (x)

Tm (x)

= r/R
0,617

0,816

x0

Figure 2.3: Conduite circulaire chauffee sous flux constant: profils de temperature et de flux radial de la chaleur.

Convection laminaire en canal rectangulaire chauffe `a flux constant La correlation de Shah et London (1978) donne le
nombre de Nusselt en fonction du rapport daspect du canal pour un chauffage `a flux constant sur les 4 faces du canal :
N u = 8, 235(1 2, 0421 + 3, 0853 2 2, 4765 3 + 1, 0578 4 0, 1861 5)
ou est le rapport du plus petit cote au plus grand et est donc inferieur `a 1. Pour un canal carre ( = 1) on obtient ainsi :
h
N u = 3, 64. Le nombre de Nusselt est defini ici comme : N u = hD
.

2.2

Convection forc
ee laminaire dans des conduites de section quelconque

La Figure 2.4 , extraite du livre de Incropera and de Witt, donne le nombre de Nusselt et le coefficient de frottement pour
des ecoulements laminaires dans des conduites de section quelconque.

Figure 2.4: Nombre de Nusselt et facteur de friction f (f tel que P = f DLh 12 V 2 ) pour des ecoulements laminaires
pleinement developpes en conduites de section quelconque (dapr`es Incropera and De Witt). La premi`ere colonne represente
le rapport daspect de la conduite, la seconde le nombre de Nusselt lorsquun flux uniforme est impose `a la paroi, la troisi`eme
le nombre de Nusselt lorsque la paroi est main tenue `a temperature uniforme.

2.3

Convection forc
ee turbulente dans un tube quelconque

Du fait de la presence de tourbillons dans lecoulement, la diffusivite apparente du fluide augmente et les transferts de
chaleur sen trouvent intensifies. On donne ici dabord quelques resultats relatifs aux ecoulements en tubes cylindriques puis
on donnera la marche `a suivre pour des conduites turbulentes de section quelconque.

Ecoulement turbulent en tube circulaire lisse infiniment long :

5
Pour L/D > 60 ; 0, 7 < P r < 100 ; 104 < ReD < 1, 2.10

Formule de COLBURN : N uD = 0, 023.P r1/3 .Re0,8
N uD =
D

T /D

Les proprietes du fluide sont evaluees `

a ( Tparoi + Tmelange)/2
On peut aussi utiliser la correlation tiree dexperiences
 de DITTUS et BOELTER :
N uD = 0, 0243.Re0,8.P r0,4
si T p > T m

N u = T
D
N uD = 0, 0265.Re0,8.P r0,3
si T p < T m
[Proprietes `a prendre `
a (Tp + Tm)/2].
Pour des nombres de Reynolds plus grands, (104 < Re < 5.106) et un nombre de Prandtl proche de celui de lair
(0.5 < P r < 1.5), la formule suivante peut etre employee:
N uD =
Tubes courts:


 0.4

= 0.0214 Re0.8
D 100 P r
T /D

Selon la longueur on utlisera:



L
< 60 N uD,L = N uD,h 1 + 6 D
- pour 20 < D
L

0.7 i
L
- pour 2 < D
< 20 N uD,L = N uD, 1 + D
L

L
> 60).
Dans ces deux derni`eres formules, N uD, est celui dun tube infiniment long ( D

La modelisation numerique permet davoir des informations plus precises sur le profil transverse de vitesse et de temperature
dans la conduite mais ce type de renseignements nest souvent necessaire que lors du developpement dun nouveau produit
et rarement lors du dimensionnement dun projet industriel.
Pour des conduites de section quelconque , on peut, `a defaut de donnees fiables, reprendre les formules valables en tubes
circulaires et remplacer D par le diam`etre hydraulique DH de la conduite quon souhaite etudier.

Chapter 3

Convection forc
ee externe
Contrairement au chapitre precedent qui decrit les echanges thermiques `a linterieur dune conduite on sinteresse ici aux
echanges lorsque le fluide secoule `
a lexterieur dun corps. On traitera ainsi le cas des ecoulements sur plaque plane, laminaires
ou turbulents, ainsi que quelques exemples decoulements autour dobstacles.

3.1

Convection forc
ee laminaire sur plaque plane

On assimile la paroi `a une plaque plane lorsquelle est effectivement plane ou bien lorsque sa concavite est suffisamment faible
pour etre negligee. Dans une situation confinee, le diam`etre hydraulique peut etre suffisamment faible pour que lecoulement
soit laminaire. Lecoulement pr`es dune paroi consiste en une zone lointaine (eloignee de la paroi), de vitesse U0 , et une zone
tr`es mince, o`
u les gradients transverses de vitesse sont eleves, permettant de respecter la condition U = 0 `a la paroi. Cette
zone est appelee couche limite hydrodynamique (voir Figure 3.1).
Cette couche limite est dune importance essentielle dans les transferts thermiques entre le fluide et la paroi : il existe
egalement une zone mince pr`es de la paroi o`
u les variations de temperature sont rapides : cest la couche limite thermique
(voir Figure 3.2). Nous allons montrer, dans le prochain paragraphe, comment les equations du mouvement et de lenergie
pourraient etre simplifiees en vue dune eventuelle resolution analytique. Puis, dans les paragraphes suivants, nous simplifierons encore le probl`eme pour faire ressortir non pas la valeur exacte du nombre de Nusselt mais son ordre de grandeur en
fonction du nombre de Reynolds de lecoulement et du nombre de Prandtl du fluide.

3.1.1

Rappel sur les couches limites hydrodynamique et thermique

Pour simplifier, on supposera :

ecoulement 2D (lecoulement na donc que 2 composantes de vitesse, une le long de la paroi et une dans la direction
normale `a celle-ci.
regime permanent /t = 0
S = 0 : pas de sources internes de chaleur : la puissance thermique associee `a la dissipation visqueuse est negligee.
proprietes physiques (, , , ) constantes
Un traitement rigoureux de la couche limite necessiterait la resolution compl`ete des equations de Navier-Stokes. Leur
complexite a donne `a PRANDTL lidee de les simplifier pour ne retenir que les termes les plus importants . Lidee principale
consiste `a negliger les gradients axiaux /x devant les gradients transverses /y. On obtient ainsi les equations de
PRANDTL de la couche limite:
u
et :

u
1 p
2u
u
+v
=
+ 2
x
y
x
y
p
=0
y

<< y
et que la pression ne varie pas suivant la direction Oy.
Ces equations ont ete etablies en supposant que x
Un raisonnement analogue fournit une equation simplifiee de la temperature:

T
T
T
T
2T
>>
u
+ v.
= 2
y
x
x
y
y
10

U0

(x)

u(x,y)

AAAAAAAAAAAAAA

Figure 3.1: Developpement dune couche limite hydrodynamique sur une plaque plane

La resolution combinee de ces 3 equations (qdm suivant x, pression et temperature) permet le traitement dune couche
limite thermique sur paroi en presence dun ecoulement. Cela reste encore complexe et lobjet de ce cours est plutt de discuter
les ordres de grandeur regissant les echanges de chaleur et le nombre de Nusselt associe. Il est pour cela utile de se rappeler
quune couche limite resulte de la competition entre deux effets qui sont:
- la diffusion (de quantite de mouvement ou de chaleur selon quon sinteresse `a une couche limite cinematique ou thermique).
Celle-ci est dominante dans la direction y normale `
a la paroi qui presente les plus forts gradients.
- la convection , c.a.d. le transport, par lecoulement lointain qui a lieu principalement suivant Ox.

Pour la couche cinematique, pendant un temps t, elle aura diffuse transversalement dune distance = t. Pendant ce
temps, la convection laura deplacee de x = U0 t. Lelimination de t dans ces deux relations fournit:

=
x

U0 x

1/2

= Rex1/2

Lexpression precedente fait apparatre le nombre de Reynolds local Rex = U0 x base sur labscisse le long de la plaque :
il varie donc e tout point de la plaque depuis la valeur 0 en x = 0 jusqu`
a sa valeur maximale en x = l o`
u il vaut Rel = U0 l .
Comme il a ete vu en cours de Mecanique des Fluides de M 1 , cette derni`ere valeur doit etre inferieure `a 105 environ pour
 1/2

x
que lecoulement reste laminaire. Lexpression de peut aussi secrire = U
qui montre levolution en x de la couche
0
limite.
Le traitement par ordres de grandeur de probl`emes thermiques necessite de considerer 2 cas extremes caracterises par P r << 1
ou P r >> 1 (voir Figure 3.3).

3.1.2

Echanges thermiques sur plaque plane pour P r << 1

Dans ce cas la chaleur diffuse plus rapidement que la quantite de mouvement et lepaisseur de la couche limite thermique
sera plus grande que celle de la couche limite cinematique:
<< t
Le meme raisonnement que pour la couche limite cin
ematique peut etre applique `a la Couche limite thermique. Pendant
un temps t elle aura diffuse transversalement de t = t. Pendant ce meme temps elle aura ete convectee `a la vitesse U0
(puisque presque toute la couche thermique se trouve plongee dans une zone se deplacant `a U0 , voir Figure 3.3-a). Ainsi,
on a :
r
x
t
U0
do`
u:
p

/ = P r1/2 << 1
t
11

On peut alors estimer le flux de chaleur extrait de la paroi:

p +.

TP T0
t

et en le normalisant par : TP xT0 , on obtient:

N ux =

x
t

soit:
1/2
N ux Re1/2
x .P r

avec:

U0 .x

Rex =

3.1.3

Echanges thermiques sur plaque plane pour P r >> 1

Le meme raisonnement sapplique au calcul de qui garde la meme expression. Toutefois, celui relatif `a t doit etre revu.
En effet, le nombre de Prandtl est maintenant tel que:
t <<
La couche limite thermique est donc incluse dans une region de plus basse vitesse que U0 , variant de 0 `a la paroi jusqu`
a une
valeur que lon peut estimer `
a t .U0 si on linearise le profil de vitesse de U . Pendant un temps t, on ecrit donc:
t =

do`
u:
t

et

x=

t
.U0 t

x1/2 1/6 1/3

1/2

U0

On peut former le rapport des epaisseurs des 2 couches:

Le nombre de Nusselt secrit maintenant:

3 / = P r1/3
t
N ux =

x
1/3
Re1/2
x .P r
t

Seul lexposant du nombre de Prandtl a change par rapport au cas P r << 1. La figure 3.4 tiree de Eckert and Drake, fournit
les resultats de lanalyse exacte (traitement complet des equations de la couche limite que nous avons evoque dans le premier
pararaphe).

Uinfini
T0

AAAAAAAAAAAAAA
AAAAAAAAAAAAAA
AAAAAAAAAAAAAA
t (x)

T( y)

TP

Figure 3.2: Couche limite thermique sur plaque plane

12

U0

U0
T0
u(y)

AAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAA
ut

T0

TP

TP

<t (Pr<1)
Mtaux liquides

>t (Pr>1)
Liquides ususels

Figure 3.3: Profils de temperature et de vitesse en fonction du nombre de Prandtl. a)- P r << 1 la couche limite se trouve
presque enti`erement situee dans une zone o`
u la vitese vaut U0 b)- P r >> 1 la couche limite thermique se trouve placee dans
une region o`
u la vitesse evolue de 0 `
a la paroi jusqu`
a 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

echang
e (en laminaire)

Le calcul complet tenant compte du profil de vitesse reel dans la couche limite est en fait possible pour certains ecoulements
et notamment celui sur plaque plane. On trouve :
1/2

Pour P r < 0.1 : N ux = 0, 565.Rex

. P r1/2

.)
(Les proprietes du fluide sont `
a prendre `
a la temperature du film T p+T
2
RL
R L N ux . T
RL
.dx
x .dx
0
Nombre de Nusselt global : N uL = /L
= 0
= 0 N ux dx
T =
T
T
x
L

soit : N uL =
ReL =

U0 .L

RL
0

0, 565.P r

1/2

Uo .x 1/2 dx
.x

1/2

Pour P r > 0.1 : N ux = 0, 332.Rex

1/2

de meme : N uL = 0, 664.ReL

1/2

= 1, 13.ReL P r1/2

avec

. P r1/3

. P r1/3

Le tableau suivant synthetise les resultats :

P r << 1
P r >> 1

1/2

N uL = 1, 13.ReL P r1/2
1/2
N uL = 0, 664.ReL P r1/3

Nombre de Nusselt pour un ecoulement laminaire dun fluide sur une plaque plane

3.2

Convection forc
ee turbulente sur plaque plane

Correlations usuelles : Lorsque le nombre de Reynolds local Rex exc`ede la valeur Rex > 3.105 on adopte les resultats suivants
:
)
4/5
0, 5 < P r < 50
N ux = 0, 029.P r1/3 .Rex
4/5
1/3

ecoulement turbulent sur toute la plaque

N uL = 0, 036.P r .ReL
Une etude detaillee du profil de temperature dans la direction transverse `a la plaque rev`ele que celui-ci est tr`es similaire
`a ce qui est observe pour le champ de vitesse (voir cours de Mecanique des Fluides de M2) qui sexprime ainsi:
u+ = y + pour y + [0; 10]
u = 0.41 est la constante de Karman.
u+ = 5.5 + 1 Lny + pour y + [100; 1000] o`
On montre donc que le profil de T secrit:
dans la zone visqueuse (y + [0; 10]): T + = P r.y + avec T + =

CP U (TP T )
P

o`
u U est la vitesse de frottement classiquement definie, TP et P la temperature et le flux `a la paroi. T designe la
temperature de melange de lecoulement.
dans la zone Log on trouve: T + = 2.195 Lny + + 13.2 P r 5.66
On voit donc que lallure generale des profils de T et de U sont tr`es semblables, que ce soit dans la zone visqueuse ou
dans la zone Log.

3.3
3.3.1

Ecoulements forc
es autour dobstacles
Obstacles cylindriques

Cette situation correspond par exemple `

a celle des echangeurs `a tubes et calandre dans lesquels le fluide de la calandre croise
les tubes. Lecoulement autour dun cylindre depend fortement de Re (voir Figure 3.5).
Si Re << 1 ( rare) : ecoulement purement visqueux symetrique
Si Re mais reste inferieur `
a 40, il existe 2 tourbillons attaches derri`ere le cylindre
14

 Si Re > 40 , il existe une allee de tourbillons de Von Karman `a laval du cylindre

Si Re encore, la couche limite qui se developpe `a partir du point darret `a = 0 devient turbulente et decolle `a d
Re < 3.105
Re > 3.105

d = 80
d = 120

et la couche limite est turbulente meme entre = 0 et = d

Figure 3.5: Ecoulement autour dun cylindre. Ref.: Favre-Marinet, cours de lENSHMG
Resultats experimentaux concernant les echanges thermiques (Fig. 3.6 et 3.7):
Si Re petit, Nu est petit (de quelques unites), surtout `a laval o`
u les recirculations sont peu actives.
Pour Re > 426 000, Nu augmente beaucoup (plusieurs centaines) avec un max vers d = 120
Synth`ese des resultats experimentaux : (Hilpert, 1993)
Re [1, 4000]
Re [4000, 40000]
Re [40000, 400000]

N u = 0, 43 + 0, 53.P r0,31 .Re0,5

N u = 0, 43 + 0, 193.P r0,31.Re0,618
N u = 0, 43 + 0, 0265.P r0,31.Re0,805

Ecoulement autour dun cylindre : nombre de Nusselt vs Re.

15

Figure 3.6: Distribution azimutale du coefficient dechange autour dun cylindre dans un ecoulement dair- a) Re petit, b)
Re intermediaire. Ref.: Eckert and Drake. Attention: lecoulement va de droite `a gauche sur la Fig.a et de bas en haut sur
la Fig.b.

Figure 3.7: Distribution azimutale du coefficient dechange autour dun cylindre dans un ecoulement dair- Re >> 1. Ref.:
Eckert and Drake (lecoulement va de gauche `
a droite).

16

Figure 3.8: Convection forcee autour dun cylindre dans lair. Nombre de Nusselt moyen en fonction de Re dapr`es Hilpert
(1933).

avec N u =

(T1 T2 )
d

et =

et Re =

V .d

Application `a la mesure de vitesse par fil chaud : La sonde est un fil de platine de diam`etre tr`es petit d 5m. Un
asservissement electronique le maintient `
a temperature constante par apport calorifique dun courant engendrant un effet
Joule. La mesure du courant i est fonction du flux echange par le fil.

fil


?
i 

V0


Calcul de Re dans lair `
a V0 = 100 m/s : la viscosite cinematique 105 m2 /s : on en deduit Re =

en regime laminaire) et le nombre de Nusselt depend de Re1/2 do`

u Nu V0

V0 d

50 (on est donc

/dL

Par ailleurs Nu = T
/d = L.T

Donc le flux est proportionnel `

a V0 : la reponse de la sonde sera non- lineaire de sorte que la sonde necessitera un
u la vitesse est tr`es bien connue.
etalonnage prealable realise dans un petite soufflerie par exemple o`

3.3.2

Obstacles non circulaires

L`
a encore on introduit un coefficient dechange h tel que:
= hST
de sorte que le nombre de Nusselt devient:
17

N ud =

T
d

.S

hST
hd
=
T.S

0,35
On pose N ud = C.Rem
avec Red = U .d/
d .P r
La figure 3.9 donne les valeurs de C et m pour diverses geometries.

Figure 3.9: Coefficients dechanges en fonction de la geometrie en convection forcee externe. (dapr`es Taine et Petit). La
fl`eche indique le sens de lecoulement.

18

3.3.3

Obstacle sph
erique

~
U

d
'$

TP

&%


Dans les conditions suivantes:

0, 6 < P r < 400

et

Re < 7.104

on a :
N ud =

hd

1/2

= 2 + 0.6 P r1/3 .Red

Ecoulement autour dune sph`ere: nombre de Nusselt en fonction de Re et P r.

19

Chapter 4

Convection naturelle
4.1

Introduction

Rappel de la definition : il sagit de la description du mouvement dun fluide engendre par les forces dArchim`ede dues aux
variations de la masse volumique avec la temperature. Il y a donc couplage de la dynamique et de la thermique. Le champ
de vitesse transporte la chaleur et, du fait de la dependance de avec T , agit sur la distribution de masse volumique ; en
retour, les modifications de masse volumique creent du mouvement par le biais de la poussee dArchim`ede.
La convection naturelle est un phenom`ene important dans de nombreux procedes industriels :
refroidissement dappareils electriques
collecteurs solaires
chauffage des locaux
centrales nucleaires
Les equations du mouvement en presence de convection naturelle ont ete presentees dans le chapitre 1 : leur forme
adimensionnelle a fait apparatre les nombres de Grashof et de Rayleigh :
Gr =

g T . l3
2

o`
u designe la dilatabilite thermique du fluide :

et

Ra = Gr.P r =

g T . l3

1
)P
= (
T

Pour un gaz parfait on montre facilement que =

1
3.103 .
obtient : = 293

1
T.

Par exemple, pour de lair `a 20C , assimile `a un gaz parfait, on

Avant de presenter les resultats dechanges thermiques par convection naturelle le long dune paroi, il faut preciser le
crit`ere qui caracterise la transition laminaire/turbulent dun ecoulement ascendant. La figure 4.1 illustre cette question et
montre que lorsque le nombre de Grashof local atteind 109 environ, lecoulement devient turbulent.
R
egime laminaire Avant de donner les expressions exactes du nombre de Nusselt il est interessant, comme nous lavions
fait en convection forcee, de comprendre leffet du nombre de Prandtl.
Lorsque P r >> 1 cela correspond donc `
a un fluide ayant plus de facilite `a diffuser la quantite de mouvement que la
chaleur. Il existe donc une couche limite thermique depaisseur t (cf. Figure 4.2) dans laquelle la temperature est plus
chaude qu`a linfini o`
u elle vaut T0 . Cette couche limite thermique est donc le si`ege de forces dArchim`ede qui el`event
le fluide et il y a creation dun jet parietal depaisseur c superieur `a t .
Lorsque P r << 1 lon sattendrait `
a ce que la couche limite thermique soit beaucoup plus epaisse que la couche
cinematique. Cela signifierait alors quil existerait une zone de fluide immobile dans laquelle la temperature est pourtant
superieure `a T0 . Cela est evidement impossible car la poussee dArchim`ede contrarierait cet equilibre instable : on
observe en fait que les couches limites cinematiques et thermiques ont la meme epaisseur (cf. Figure 4.3).

20

AA
AA
AA
AA
AA
AA
AA
AA
AA
AA

Gr L = 10 9

Figure 4.1: Convection naturelle le long dune plaque verticale: transition `a la turbulence

AAA
AAA
AAA
AAA
AAA
AAA
AAA
AAA
AAA
AAA
AAA
z

T1

g
L

T0

Figure 4.2: Convection naturelle le long dune plaque plane verticale pour P r >> 1.

21

AAA
AAA
AAA
AAA
AAA
AAA
AAA
AAA
AAA
AAA
AAA
z

T1

g
L

T0

c
t

Figure 4.3: Convection naturelle le long dune plaque plane verticale pour Pr << 1.

En resolvant le probl`eme de facon analytique, pour un nombre de Prandtl quelconque et pour une paroi de temperature
imposee, on trouve une formule unique :

N uz =
N uL =

hz z

hL

1/2

1/4

1/2

1/4

Pr
.Grz
= 0, 508. (0.95+P
r)1/4

Pr
.GrL
= 0, 677. (0.95+P
r)1/4

R
egime Turbulent : les formules ci dessus ne sont valables quen regime laminaire cest `a dire quand Gr < 109
Pour Gr > 109 , le regime est turbulent et MAC ADAMS propose la correlation :


N uL = 0.13(GrL .P r)1/3 valable pour RaL 109 , 1013

4.2

Convection naturelle sur plaque plane verticale chauff

ee `
a flux constant

Cette situation se rencontre lorsque la paroi contient par exemple une source de chaleur electrique dont la puissance est
evacuee par la paroi : le flux `
a la paroi est alors impose. Les resultats suivants seront alors utilises:
Pour 105 < Grz P r < 1011 (regime laminaire):
N uz = 0.6(Grz P r)1/5
Pour 2.1013 < Grz P r < 1016 (regime turbulent):
N uz = 0.568(Grz P r)0.22
Le nombre de Grashof local modifie secrit dans ce cas :
gz 3 P z
gP z 4
=
2
2
designe le flux a
` la paroi.
Grz =

Dans cette expression, P

22

4.3

Convection naturelle entre plaques verticales parall`

eles (chemin
ee)

Si la longueur L>>2e, on aura un regime hydrodynamique etabli.

On cherche alors un champ de vitesse de la forme :
Ux
Uy

4.3.1

= u(y)
= 0

Condition dexistence de ce r
egime

Pour une plaque, nous avons vu que lepaisseur de la couche limite etait donnee par :

1/4
= RaL
L

lorsque P r 1(air)

1/4

Il y aura regime etabli si >> e soit Le .RaL << 1

Sous forme dimensionnelle cela secrit :
e
L1/4

gT

1/4

<< 1

Pour obtenir ce regime, on peut donc e ou T ou L.

4.3.2

Equations du r
egime
etabli

Les hypoth`eses sur le champ de vitesse imposent :

2
N Sx : p
0 g(T T0 ) + yu2 = 0
x = 

2 T
2 T
T : u T
x = x2 + y 2
p* designe la pression generalisee : p = p + 0 gx. En l absence de mouvement et de temperature, la solution sera bien
hydrostatique : p* = constante ; T = T0 .

4.3.3

Cas de deux plaques de m

eme temp
erature

TP

T0

AA
AA
AA
AA
AA
AA
AA
AA

AAA
AAA
AAA
AAA
AAA
AAA
AAA
AAA

TP

2e

T0

Figure 4.4: Convection naturelle entre deux plaques de meme temperature

Lequation de la chaleur se simplifie puisqu ici on a
2
yT2 = 0 T
y = cste = 0

T
x

=0

T = Tp partout ; le fluide est isotherme. Il subit la poussee dArchim`ede par rapport au fluide externe `a TO . Cette
poussee sera contree par la viscosite.
Pour la poussee : p* (x = 0)= p* (x = L) car cest le fluide externe qui communique la pression ambiante hydrostatique
aux 2 extremites de la cheminee. On a donc :
p
= 0 do`
u N-S se simplifie :
x
23

2U
y 2

g(Tp TO )

u=

ge2
2 (Tp

do`
u le debit par unite de largeur :

h
TO ) 1
q/l =

y 2
e
Z

+e

udy =
e

2 ge3
(Tp TO )
3

Bilan thermique entree/sortie :

3

2
= 0 CP ql (TP T0 ) = 0 .Cp 2ge
3 (Tp T0 )
Ce flux a ete necessaire pour faire passer le fluide de T0 en entree `a Tp en sortie. Il a en fait ete donne dans la zone
detablissement (dans les couches limites), seules regions o`
u le gradient de temperature `a la paroi nest pas nul. Cela explique
labsence de L dans lexpression de .
Nombre de Nusselt : N uL = /2L.l
Tp T0 (facteur 2 car il existe 2 plaques)

La presence de L au denominateur indique que cest un flux conductif axial qui sert de reference; cest en effet selon x
quil existe un gradient, pas selon y.
On trouve ainsi :
N uL =
avec Ra2e =

4.4
4.4.1

gT (2e)3

Ra2e
24

hL

(T = Tp T0 )

Autres g
eom
etries
Paroi plane inclin
ee par rapport `
a la verticale

AAAAAA
AAAAAA
AAAAAA
AAAAAA
AAAAAA
AAAAAA
AAAAAA

TP

T0

Figure 4.5: Convection naturelle le long dune paroi inclinee

Il suffit de remplacer g par g cos dans les formules trouvees en paroi verticale (voir Figure 4.5).

4.4.2

Cylindre vertical

Leffet de courbure sera negligeable si lepaisseur de la couche limite est petite devant le diam`etre d du cylindre :
1/4
<< d Pour un gaz (P r 1) par exemple letude sur plaque plane donne : L = RaL
La condition de courbure negligeable devient donc :
1/4

= LRaL

<< d

soit :

L
1/4

dRaL

Introduisons le nombre de Rayleigh, Rad , base sur le diam`etre du cylindre :

Rad = RaL .

24

d3
, do`
u:
L3

<< 1

d
L

Figure 4.6: Convection naturelle le long dun cylindre vertical

d3/4
L
d Ra1/4 L3/4
d

<< 1

Rad . Ld

1/4

>> 1

Si cette hypoth`ese est satisfaite on pourra utiliser les formules obtenues pour plaque plane :
1/4
Pour lair :
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.3

Cylindre horizontal

On utilise encore le diam`etre d comme longueur de reference.

Avec une erreur < 15 % on prendra :
1/4
en laminaire N ud = 0, 53.Rad
Ra[103 , 109 ]
1/3
en turbulent N ud = 0, 10.Rad
Ra[109, 1013 ]

4.4.4

Sph`
ere
1/4

On prend N ud = 2 + 0, 45Grd .P r1/3 , Grd [1, 106 ].

4.4.5

Plaques horizontales

Voir tableau 4.7. La dimension caracteristique `

a introduire dans le calcul du nombre de Rayleigh est la largeur de la plaque.
Plaque chauffant vers le haut
Plaque chaufant vers le bas

N u = 0, 54.Ra1/4
N u = 0, 14.Ra1/3
N u = 0, 27.Ra1/4
N u = 0, 07.Ra1/3

Ra[105, 107 ]
Ra[107 , 1010 ]
Ra[3.105, 3.1010 ]
Ra[3.1010, 1013 ]

25

Figure 4.7: Coefficient dechange en convection naturelle: parois verticales ou horizontales N u = C.Ran . Source: Eyglunent,
Thermique theorique et pratique, Herm`es

26

Formulaire
Gradient dun scalaire en coordonnees cylindriques (r, , z):
~ f
grad

f
r
f
r
f
z

~ =
Divergence dun vecteur en coordonnees cartesiennes: div V
~ =
Divergence dun vecteur en coordonnees cylindriques: div V
Continuite en cartesien :

vy
vz
vx
x + y + z
vz
1 (rvr )

+ 1r v
r r
+ z

vx
vy
vz
+
+
+
=0
t
x
y
z
Continuite en cylindrique:
r)
z)
)
+ 1r (rv
+ r1 (v
+ (v
=0
r

z
Equations de Navier Stokes en coordonnees cartesiennes (x, y, z)



 2
p
vx
vx
vx
2 vx
2 vx
vx
vx
=

+
+
+ vx
+ vy
+ vz
+ gx +
t
x
y
z
x
x2
y 2
z 2



 2
p
vy
vy
vy
2 vy
2 vy
vy
vy
=

+
+
+ vx
+ vy
+ vz
+ gy +
t
x
y
z
y
x2
y 2
z 2



 2
vz
p
vz
vz
vz
2 vz
2 vz
vz

=
+
+
+ vx
+ vy
+ vz
+ gz +
t
x
y
z
z
x2
y 2
z 2

Equations de Navier Stokes en coordonnees cylindriques (r, , z)

vr
vr
v vr
v2
vr
+ vr
+
+ vz
t
r
r
r
z

p
+ gr +
=
r

v
v
v
v v
vr v
+ vr
+
+
+ vz
t
r
r
r
z

1 p
=
+ g +
r

vz
vz
vz
v vz
+ vr
+
+ vz
t
r
r
z

1 vr
1 2 vr
2 v
vr
2 vr

(r
) 2 + 2
+
r r r
r
r 2
r2
z 2


1 2 v
2 vr
v
2 v
1 v

(r
) 2 + 2
+
r r r
r
r 2
r2
z 2

p
+ gz +
z

2 vz
1 2 vz
1 vz
+
(r
)+ 2
2
r r r
r
z 2

Contraintes visqueuses dans un fluide newtonien en coordonnees cartesiennes:

xx = 2

vx
x

yy = 2

vy
y

zz = 2

vz
z

27

xy = yx =

vx
vy
+
y
x

yz = zy =

vz
vy
+
z
y

vz
vx
=
+
x
z

zx = xz

Contraintes visqueuses dans un fluide newtonien en coordonnees cylindriques:

rr = 2

vr
r

vz
zz = 2
z


1 vr
v
r = r = r ( ) +
r r
r


v
1 vz
z = z =
+
z
r


vr
vz
+
zr = rz =
r
z
= 2

Equations de l energie en coordonnees cartesiennes (avec et constants):



 2

T
T
T
T
T
2T
2T
=
CP
+q+D
+
+
+ vx
+ vy
+ vz
t
x
y
z
x2
y 2
z 2
Equations de l energie en coordonnees cylindriques:




T
1 T
T
T
2T
v T
1 2T
CP
=
+q+D
+ vr
+
+ vz
(r
)+ 2 2 +
t
r
r
z
r r r
r
z 2
avec q source de chaleur interne (electrique, chimique, nucleaire, etc... en W/m3 ) et D puissance des forces visqueuses
(en W/m3 egalement mais souvent negligeables) telle que:
en cartesien :
"
"
2 
2 
2 #
2 
2 
2 #
vy
vz
vz
vx
vy
vx
vy
vz
vx
+
+
+
+
+
+
+
+
D = 2
x
y
z
x
y
y
z
z
x
en cylindrique :
"
2 
2 
2 #

2
2
2


vr
 v  1 vr
1 v
vr
vz
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 tirees des ouvrages suivants qui pourront etre consultes avec interet.

BEJAN A: Convection Heat Transfer, Edition Wiley.

INCROPERA F. DEWITT D.: Fundamentals of Heat and Mass Transfer.
PETIT J.P., Transfert de chaleur et de masse, Cours de lEcole Centrale de Paris.
FAVRE MARINET M.: Cours de Convection Thermique, ENSHMG.
EYGLUNENT B., Manuel de thermique, Hermes.
SOLLIEC C., VIAZZO S. : Cours dEchanges Thermiques de lEcole des Mines de Nantes

29