Convection F
Convection F
Convection F
EDEDD
PLAN du Cours de Convection forcée
1. INTRODUCTION
EDEDD : 2007 ML 1
Convection forcée
EDEDD
EDEDD : 2007 ML 2
1. Introduction
EDEDD
Exemples de convection
EDEDD : 2007 ML 3
1. Introduction
EDEDD
Objectifs
EDEDD : 2007 ML 4
2. Principes fondamentaux de la convection
EDEDD
2.1. Phénomènes physiques
y
Fluide T∞
convections
sens de mécanique et
l’écoulement diffusion thermique
visqueuse et
conductive
0 Fluide T(y) x
paroi Tp en y = 0
(Tp )y = 0 = (Tfluide )y = 0
∂Tsol ∂T
λ sol = λfluide fluide
∂y y =0 ∂y y =0
EDEDD : 2007 ML 5
2. Principes fondamentaux de la convection
EDEDD
2.1. Phénomènes physiques
EDEDD : 2007 ML 6
2. Principes fondamentaux de la convection
EDEDD
2.1. Phénomènes physiques
● conduction
- flux de chaleur surfacique conductif, dans la paroi, contrôlée par la
conductivité de la paroi
- flux de chaleur surfacique conductif dans le fluide, contrôlée par la
conductivité du fluide : conduction à l’échelle moléculaire
- flux de chaleur : appelé flux convectif à la paroi.
● transport de masse
- l’échange convectif dépend étroitement du transport de masse
- fonction du champ de vitesses
- fonction de la nature de l’écoulement
EDEDD : 2007 ML 7
2. Principes fondamentaux de la convection
EDEDD
2.2. Différentes classes de problèmes en convection
❒ Origine du mouvement
● Convection forcée
mécanisme moteur : mouvement du fluide créé par un organe mécanique
● Convection naturelle
mécanisme moteur : mouvement du fluide créé par une différence de masse
volumique du fluide liée à un gradient de températures
● Convection mixte
mouvement du fluide imposé par un organe extérieur et mouvement du fluide
imposé par un gradient de température du même ordre de grandeur
EDEDD : 2007 ML 8
2. Principes fondamentaux de la convection
EDEDD
2.2. Différentes classes de problèmes en convection
EDEDD : 2007 ML 10
2. Principes fondamentaux de la convection
EDEDD
2.2. Différentes classes de problèmes en convection
u∞ u∞
δ T∞
y ∆ y
EDEDD : 2007 ML 11
2. Principes fondamentaux de la convection
EDEDD
2.2. Différentes classes de problèmes en convection
EDEDD : 2007 ML 12
2. Principes fondamentaux de la convection
EDEDD
2.2. Différentes classes de problèmes en convection
● température
● flux
EDEDD : 2007 ML 13
2. Principes fondamentaux de la convection
EDEDD
2.3. Equations de base 2.3.1. Equations de bilans
∂u
v µdx
∂y
EDEDD : 2007 ML 15
2. Principes fondamentaux de la convection
EDEDD
2.3. Equations de base 2.3.1. Equations de bilans
φ = hS (T p − Tf )
h : coefficient d’échange convectif local ou moyenné sur une surface W/m2K
Tp : température de la paroi
Tf : température caractéristique du fluide
∂Tsol ∂T
ϕ = h( Tp − Tf ) = − λsol = − λfluide fluide
∂y y =0 ∂y y =0
EDEDD : 2007 ML 18
2. Principes fondamentaux de la convection
EDEDD
2.4. Détermination du coefficient d’échange
➨ Méthode analogique
Analogie entre transfert de chaleur et de quantité de mouvement
Calcul du coefficient d’échange à partir du coefficient de frottement
➨ Méthode expérimentale
Etude expérimentale
Généralisation des résultats par introduction de nombres adimensionnels
Lois empiriques pour déterminer h
EDEDD : 2007 ML 19
2. Principes fondamentaux de la convection
EDEDD
2.5. Nombres adimensionnels
❐ Analyse dimensionnelle
ϕ = f ( λ , ρ , c p , µ , V , ∆T , L )
Théorème de Vaschy-Buckingham
7 paramètres indépendants, 4 unités fondamentales 4 NSD
EDEDD : 2007 ML 20
2. Principes fondamentaux de la convection
EDEDD
2.6. Similitude
❐ Similitude en convection
EDEDD : 2007 ML 21
2. Principes fondamentaux de la convection
EDEDD
2.7. Détermination des variables intervenant dans Nu
❐ Géométrie externe
loin de la paroi, gradients de vitesse et de température négligeables
- plaque plane semi-infinie
- cylindre
- sphère
❐ Hypothèses
CF : profils de vitesses indépendants du champ de températures
Propriétés physiques du fluide constantes
Ec<<1 Nux = f(x, Re , Pr) Nu = f (Re, Pr)
Régime laminaire : nombre de Reynolds critique < 4 à 5. 105
∂P ∂ P
= =0 si l’écoulement loin de la paroi a une vitesse constante
∂y ∂x
Fluide T∞
y
sens de
l’écoulement v
V∞ u x
T∞ 0 0
EDEDD : 2007 ML 23
3. Convection forcée externe
EDEDD
3.2. Ecoulement le long d’une plaque plane semi-infinie
∂T ∂T ∂ 2T
u +v = a
∂x ∂y ∂y 2
V x
❐ Epaisseurs des couches limites δ =
4 ,92 x
avec Re x = ∞
Re x υ
δ en x
∆
Fluides usuels : Pr ≥ 0,7 ≈ Pr − 1 3 Pr élevé : ∆ < δ
δ
Pr faible : ∆ > δ
EDEDD : 2007 ML 24
3. Convection forcée externe
EDEDD
3.2. Ecoulement le long d’une pplaque plane semi-infinie
u v u∞ x
Variation adimensionnelle = f ( y u ∞ υx ) Variation adimensionnelle = f ( y u ∞ υx )
u∞ u∞ υ
EDEDD : 2007 ML 25
3. Convection forcée externe
EDEDD
3.2. Ecoulement le long d’une plaque plane semi-infinie
❐ Contrainte pariétale ∂u µu ∞
τ p = µ = 0,332 Re x
∂ y y =0 x
τp
Coefficient de frottement f =
1
ρu ∞
2
2
−1 2 −1 2
f x = 0,664 Re x f L = 1,328 Re L
EDEDD : 2007 ML 26
3. Convection forcée externe
EDEDD
3.2. Ecoulement le long d’une plaque plane semi-infinie
T − Tp
Variation adimensionnelle = f ( y u ∞ υx
T∞ − T p
EDEDD : 2007 ML 27
3. Convection forcée externe
EDEDD
3.2. Ecoulement le long d’une plaque plane semi-infinie
❐ Coefficient d’échange
12
Re x
-Tp imposée Nu x = 0,332 Re1x 2 Pr 1 3 h x = 0,332λ Pr 1 3 = Kx − 1 2
x
L 12
1 ReL
h = ∫ h x dx h = 0,664λ Pr 1 3
L L
0
h L
Nu L = L NuL = 0,664 ReL1 2 Pr 1 3
λ
- flux imposé
12
Re x
12
Nu x = 0,453 Re x Pr 1 3 h x = 0,453λ Pr 1 3 = Kx − 1 2
x
12
NuL = 0,6795 ReL1 2 Pr 1 3 h = 0,906λ
ReL
Pr 1 3
L
Nu ( flux imp )
= 1,36
Nu ( T imp )
EDEDD : 2007 ML 28
3. Convection forcée externe
EDEDD
3.2. Ecoulement perpendiculaire à un cylindre
❐ Géométrie interne
confinement du fluide entre des parois
- tube cylindrique à section circulaire ou non circulaire
- plaques parallèles proches
selon normale à la paroi et loin de la paroi, modification du champ de V et T
❐ Hypothèses
EDEDD : 2007 ML 30
4. Convection forcée laminaire en géométrie interne
EDEDD
4.2. Convection dans un tube à base circulaire
❐ Considérations hydrodynamiques
EDEDD : 2007 ML 31
4. Convection forcée laminaire en géométrie interne
EDEDD
4.2. Convection dans un tube à base circulaire
Zone d’établissement
Le régime est toujours laminaire à l’entrée, pour x croissant : 2 cas
CL se rejoignent avant que le régime turbulent soit atteint
régime laminaire en régime établi
CL turbulentes avant de se rejoindre
régime turbulent en régime établi
Recrit =
(ρV )D = GD = 2200 G : vitesse massique
µ µ
EDEDD : 2007 ML 32
4. Convection forcée laminaire en géométrie interne
EDEDD
4.2. Convection dans un tube à base circulaire
❐ Considérations thermiques
∂P ∂ 2v ∂P ∂ P
2ème composante − +µ =0 − = =0 P = P (x )
∂y ∂y 2
∂y ∂r
∂ ∂ u r ∂P
r =
∂ r ∂r µ ∂ x
r 2 dP
Intégration : u= + a ln r + b sur l’axe en r=0 ln(r) ➜ ∞ a=0
4 µ dx R 2 dP
en r = R u = 0 b=−
4 µ dx
EDEDD : 2007 ML 34
4. Convection forcée laminaire en géométrie interne
EDEDD
4.2. Convection dans un tube à base circulaire
u=−
dP 1
dx 4 µ
(
R2 − r 2 ) u max =
dP R 2
dx 4 µ
r 2 r 2
u = u max 1 −
= 2u ∞ 1 −
2
R R 2
3 y 2
u = u ∞ 1 − u max = 1,5u ∞
2 e
16 64
Coefficient de frottement f = Coefficient de pertes de charge λ =
Re Re
Lh
Longueur d’établissement dynamique = 0,0575 Re D
D
EDEDD : 2007 ML 35
4. Convection forcée laminaire en géométrie interne
EDEDD
4.2. Convection dans un tube à base circulaire
∂u * ∂T *
=0 =0
∂x ∂x
EDEDD : 2007 ML 36
4. Convection forcée laminaire en géométrie interne
EDEDD
4.2. Convection dans un tube à base circulaire
Lt
Longueur d’entrée thermique : = 0,05Pe D = 0,05 ReD Pr
R
EDEDD : 2007 ML 38
4. Convection forcée laminaire en géométrie interne
EDEDD
4.2. Convection dans un tube à base circulaire
Si Pr > 1 Lt >> Lh
Si Pr < 1 Lt << Lh
EDEDD : 2007 ML 39
4. Convection forcée laminaire en géométrie interne
EDEDD
4.2. Convection dans un tube à base circulaire
∂T ∂ 2T 1 ∂T
Equation de conservation de l’énergie : u = a +
∂x ∂r 2 r ∂r
Tm − Tp ∞
G x
Température de mélange : = 8 ∑ n exp − λ2n
T∞ − Tp λ2
n =0 n R Re Pr
EDEDD : 2007 ML 40
4. Convection forcée laminaire en géométrie interne
EDEDD
4.2. Convection dans un tube à base circulaire
ϕx =
R 0
∑ G n e R Re Pr T − T
p ∞
DPe
0,0668
Nu D = 3,66 + L
23 h
DPe
1 + 0,04
L
Lt
Longueur d’entrée thermique : = 0,067Pe D
R
EDEDD : 2007 ML 41
4. Convection forcée laminaire en géométrie interne
EDEDD
4.2. Convection dans un tube à base circulaire
EDEDD : 2007 ML 42
4. Convection forcée laminaire en géométrie interne
EDEDD
4.2. Convection dans un tube à base circulaire
NuD , x = 1,356Gz1x 3 hx
0,104GzL
NuD ,L = 3,66 +
(1 + 0,016GzL )0,8 h
EDEDD : 2007 ML 43
4. Convection forcée laminaire en géométrie interne
EDEDD
4.2. Convection dans un tube à base circulaire
EDEDD : 2007 ML 44
5. Convection forcée en régimeEDEDD
turbulent
5.1. Généralités
❐ Généralités
- Phénomène existant dans la majeure partie des applications en convection
forcée
EDEDD : 2007 ML 45
5. Convection forcée en régimeEDEDD
turbulent
5.1. Généralités
❐ Généralités
- Transition entre les régimes laminaire et turbulent
Convection externe : écoulement sur une plaque plane
EDEDD : 2007 ML 46
5. Convection forcée en régimeEDEDD
turbulent
5.1. Généralités
- Transition entre les régimes laminaire et turbulent
u x
Convection externe Re xcrit = ∞ crit ≈ 5.10 5
υ
Zone de transition : apparition de tourbillons intermittents ou écoulement laminaire
Couche limite turbulente : 3 zones dans la couche limite pour Re > 5. 105
- proche paroi : sous-couche visqueuse (échelle microscopique)
- couche tampon
- zone turbulente (définie statistiquement) :
tourbillons à l’échelle macroscopique (taille ≈ épaisseur de la couche limite, créés par la
contrainte de cisaillement ∂u ∂y )
brassage du fluide tend à homogénéiser V et T sur ∆t
tourbillons ou cellules de turbulence dégénèrent en structures plus petites par
transmission de leur énergie cinétique ➜ dissipation visqueuse (échelle moléculaire)
dissipation de l’énergie cinétique en chaleur : taux de dissipation plus élevé qu’en
laminaire car les vitesses de déformation sont plus élevées
EDEDD : 2007 ML 47
5. Convection forcée en régimeEDEDD
turbulent
5.1. Généralités
-Transition entre les régimes laminaire et turbulent
Convection interne
EDEDD : 2007 ML 48
5. Convection forcée en régimeEDEDD
turbulent
5.1. Généralités
- Transition entre les régimes laminaire et turbulent
Convection interne ρu ∞ D h
ReDh = ≈ 2200
µ
Diamètre hydraulique Dh = 4 sec tion de passage du fluide pour un tube à section constante
périmètre mouillé
- Grandeurs turbulentes
écoulement turbulent = écoulement principal + fluctuations aléatoires
échelle des fluctuations grande/libre parcours moyen des molécules milieu continu
équations de bilans instationnaires difficiles à résoudre
approche statistique
EDEDD : 2007 ML 49
5. Convection forcée en régimeEDEDD
turbulent
5.2. Equations de bilans
grandeur instantanée = grandeur moyennée sur une durée grande devant celle
des fluctuations turbulentes (c’est-à-dire indépendante du temps)
+ grandeur fluctuante (valeur moyenne nulle)
EDEDD : 2007 ML 50
5. Convection forcée en régimeEDEDD
turbulent
5.2. Equations de bilans
∂u ∂u ∂u ∂u 1 ∂P 1 ∂ ∂u 1 ∂ ∂u 1 ∂ ∂u
+u +v +w =− + µ − ρu ' 2 + µ − ρu ' v ' + µ − ρu' w '
∂t ∂x ∂y ∂z ρ ∂x ρ ∂x ∂x ρ ∂y ∂y ρ ∂z ∂z
Equation de l’énergie
∂T ∂T ∂T ∂T ∂ ∂T ∂ ∂T ∂ ∂T
ρc p +u +v +w = − − λ + ρc p u ' T ' − − λ + ρc p v ' T ' − − λ + ρc p w ' T '
∂t ∂x ∂y ∂z ∂x ∂x ∂y ∂y ∂z ∂z
Viscosité turbulente µt ou νt
bilan des forces agissant sur un élément de fluide contrainte tangentielle
∂u
τ =µ − ρu' v' ρu' v' <0 second terme : viscosité supplémentaire
∂y
∂u ρu' v'
Analyse de Prandtl : τ = (µ + µ t ) avec µt = − viscosité turbulente
∂y ∂u ∂y
(Pa/s)
µt
νt = diffusivité thermique turbulente (viscosité cinématique turbulente) (m2/s)
ρ
Conductivité et diffusivité thermiques turbulentes λt ou at
densité de flux thermique : ϕ = −λ ∂T + ρc p v ' T '
∂y
ϕ ∂T v'T '
= −(a + at ) avec at = − diffusivité thermique turbulente
ρc p ∂y ∂T ∂y
(m2/s)
λt = at ρc p conductivité thermique turbulente (W/m.K)
Prt = ν t at
EDEDD : 2007 ML 52
5. Convection forcée en régimeEDEDD
turbulent
5.2. Equations de bilans
❐ Méthodes de fermeture
EDEDD : 2007 ML 54
5. Convection forcée en régimeEDEDD
turbulent
5.3. Ecoulement dans un tube
+ u + = −3,05 + 5 ln y + u + = 5 ,5 + 2,5 ln y +
u u+ = y +
lm = 0,4 yy + / 26
− y + lm = 0,4 y
lm lm = 0,4 y 1 − exp
26
Contrainte tangentielle
νt 1− y R
= −1
ν +
∂u ∂y +
EDEDD : 2007 ML 55
5. Convection forcée en régimeEDEDD
turbulent
5.3. Ecoulement dans un tube
y+
1− y R
Equation de l’énergie T+ = ∫ dy +
at ν + 1 Pr
0
T+
1
T + = y ++
+
ln y + + + 5 ,5
2
T + = −3,05 + ln y + + T =
0,44 0,44
EDEDD : 2007 ML 56
5. Convection forcée en régimeEDEDD
turbulent
5.3. Ecoulement dans un tube
EDEDD : 2007 ML 57
5. Convection forcée en régimeEDEDD
turbulent
5.3. Ecoulement dans un tube
Equation statistique de bilan d’énergie cinétique moyennée sur ∆t : i, j, k sur les 3 axes
∂k ∂k ∂ µ ∂v j
ρ +vj = µ + t
∂k
+ µ t ∂v i ∂v i + − ρε
∂t ∂x j ∂x j σk ∂x j
∂x j ∂x j ∂x i
∂ε ∂ε ∂ µ ∂ε ∂v j
+ C1µ t ε ∂v i ∂v i + − ρC 2 ε
2
ρ + ρv j = µ + t
∂t ∂x j ∂x j σε ∂x j
k ∂x j ∂x j ∂x i k
Grandeurs locales
f x = 0,0576 Re x−0 ,2 δ x = 0,376 x Re x−0,2 Nu x = 0,0288 Pr 1 3 Re 0x,8
f = [0,072 Re − 0 ,2
L ]
L − 0,072 Re x− 0 ,2 + 1,33 Re x crit − 0,5 x crit / L xcrit Recrit = 5.105
crit
(
NuL = 0,036 Pr 1 3 ReL0,8 − 23000 )
EDEDD : 2007 ML 59
5. Convection forcée en régimeEDEDD
turbulent
5.4. Corrélations
Re> 5.105
Coefficient de frottement
ou 1
(
= −0,8 + 0,87 ln Re 4f ) équation de Von Karman
4f
EDEDD : 2007 ML 60
5. Convection forcée en régimeEDEDD
turbulent
5.4. Corrélations
Relation de Dittus-Boelter 0,7 < Pr < 120 104 < Re < 1,2.105 (propriétés à Tm)
0 ,8 n = 0,4 en chauffage
Nu Dh = 0,023 ReD Pr n
h n = 0,3 en refroidissement
EDEDD : 2007 ML 61
5. Convection forcée en régimeEDEDD
turbulent
5.4. Corrélations
Relation de Petukhov 0,5 < Pr < 2000 104 < Re < 5.106 (propriétés à Tm)
( f 8 ) Re Pr
Nu Dh =
(
1,07 + 12 ,7( f / 8 ) 0 ,5 Pr 2 3 − 1 )
f = ( 1,82 log10 Re − 1,64 )−2
Relation de Gnielinski 0,5 < Pr < 2000 104 < Re < 5.106 (propriétés à Tm)
( f 8 )(Re − 1000 ) Pr
Nu Dh =
(
1 + 12 ,7( f / 8 ) 0 ,5 Pr 2 3 − 1 )
f = ( 0,79 ln Re − 1,64 ) −2
EDEDD : 2007 ML 62
5. Convection forcée en régimeEDEDD
turbulent
5.4. Corrélations
D
0 ,7
Nu D ,L = Nu Dh 1 + h
L
EDEDD : 2007 ML 63
5. Convection forcée en régimeEDEDD
turbulent
5.4. Corrélations
EDEDD : 2007 ML 64