Chapitre 3
Chapitre 3
Chapitre 3
III.1. Définitions
Nous appelons convection thermique la transmission d'énergie-chaleur entre une
surface et un fluide se déplaçant le long de la surface. Dans les chapitres précédents nous
avons appliqué la convection comme condition aux limites pour les problèmes de
conduction avec le coefficient de convection donné.
Dans les chapitres suivants nous discutons les mécanismes de transmission par
convection et nous développons des méthodes pratiques pour la détermination du
coefficient de convection h.
La transmission de chaleur par convection est étroitement liée à la mécanique des
fluides car le transport de chaleur se fait ici essentiellement par le mouvement des
particules du fluide dans le voisinage de la surface de contact.
On peut considérer deux types de convection, en fonction des causes qui produisent
le mouvement du fluide: la convection forcée et la convection libre ou naturelle.
La convection forcée: se manifeste quand le mouvement du fluide est une
conséquence des actions extérieures imposées (pompe, ventilateur, vent, etc…). Dans ce
cas , le champ de température est convecté par un écoulement extérieur imposé.
La convection libre ou naturelle: est produite par un mouvement du fluide, engendré
par une différence de masse volumique entre les parties froides et chaudes du fluide. La
distribution de température engendre son propre mouvement en créant des forces
d’Archimède rotationnelles.
Quand les deux types de convection existent simultanément sans que l’une soit
négligeable par rapport à l’autre la convection est dite mixte.
Page 1
Chapitre V : Principes Fondamentaux de la convection thermique
φ = h. S. ∆T (III. 1)
Page 2
Chapitre V : Principes Fondamentaux de la convection thermique
- La température T.
Ce qui permet ensuite d’évaluer les transferts de chaleur. Un nombre suffisant de
conditions aux limites et initiales doit être connu.
Pour calculer les six fonctions indiquées ci-dessous, il est nécessaire de disposer de
six équations. Pour cela nous allons établir des relations traduisent les principes de
conservation de la mécanique et de la thermodynamique. Nous obtiendrons ainsi les
équations de conservation de :
La masse,
La quantité de mouvement (équation vectorielle équivalente à 3 équations
scalaires),
L’énergie.
Enfin comme derrière équation nous utiliserons une relation décrivant l’état du fluide.
III.5. Equations de conservation
Equation de conservation de la masse (Equation de continuité):
Considérons l’écoulement d’un fluide et effectuons un bilan matière sur le système
constitué par l’élément parallélépipédique de fluide de côtés ∆x, ∆y et ∆z et de masse dm :
La conservation de la masse dans ce volume entre les instants t et t+dt peut s’écrire :
Masse entrante + Masse initiale + Masse générée = Masse finale + Masse sortante
∂ρ ∂ρ
+ ∇. [ρU] = ⃗)=0
+ div(ρ. U (III. 2)
∂t ∂t
où : ρ Masse volumique
⃗ Vitesse
𝑈
𝜕𝜌
⃗ )=0
- Pour un écoulement permanent ( 𝜕𝑡 = 0): div(ρ. U
- Pour un écoulement permanent, incompressible (ρ=Cte) : div(U⃗ )=0
Equations de quantité de mouvement :
Il s’agit de la conservation de la quantité de mouvement :
DU ∂
ρ. = ρ. Fi + (σ ) (III. 3)
Dt ∂xj ij
Fi représente la densité massique des forces extérieures (l’exemple le plus courant : la
pesanteur, pour laquelle ⃗F = ⃗g )
Page 3
Chapitre V : Principes Fondamentaux de la convection thermique
DU ∂P ∂τij
ρ. = ρ. Fi + + (III. 5)
Dt ∂xi ∂xj
La signification physique de cette équation apparaît clairement:
Variation de la Forces Forces associés à Contraintes
quantité de = volumiques + la pression par + visqueuses par unité
mouvement par unité par unité de unité de volume de volume
de volume volume
On adoptera les notations suivantes pour les forces de cisaillement (frottement visqueux)
s’appliquant sur les surfaces d’un volume élémentaire dx dy dz :
- xx, yy, zz : composante normale de la force de surface suivant les directions Ox, Oy et
Oz.
- xy, yx, xz, zx, yz , zy : composante tangentielle de la force de surface : le 1er indice
indique la direction normale à la surface considérée, le 2ème indice indique la direction dans
laquelle la composante agit.
La loi de Newton permet alors d’écrire selon chacune des 3 directions :
Du ∂τxx ∂τyx ∂τzx ∂P
ρ = ρ. Fx + + + −
Dt ∂x ∂y ∂z ∂x
Dv ∂τxy ∂τyy ∂τzy ∂P
ρ = ρ. Fy + + + − (III. 6)
Dt ∂x ∂y ∂z ∂y
Dw ∂τxz ∂τyz ∂τzz ∂P
ρ = ρ. Fz + + + −
{ Dt ∂x ∂y ∂z ∂z
-Cas d’un fluide Newtonien
Un fluide est dit newtonien si les contraintes de cisaillement sont proportionnelles aux
gradients de vitesse,
on a dans ce cas :
∂u ∂v ∂u ∂w ∂v ∂w
τxy = τyx = μ ( + ) ; τxz = τzx = μ ( + ) ; τyz = τzy = μ ( + );
∂y ∂x ∂z ∂x ∂z ∂y
∂u 2 ∂v 2 ∂w 2
τxx = μ (2 − ∇. U) ; τyy = μ (2 − ∇. U) ; τzz = μ (2 − ∇. U)
∂x 3 ∂y 3 ∂z 3
En remplaçant le tenseur des contraintes par cette expression dans l’équation générale du
mouvement déduite de la loi de Newton on obtient finalement l’équation de Navier-Stokes:
DU 1 1
= F𝑖 − ∇P + ν. ∇2 U + ν. ∇(∇. U) (III. 7)
Dt ρ 3
Pour un fluide incompressible : ∇. U = 0 et on obtient une forme simplifiée de l’équation
de Navier-Stokes :
DU 1
= F𝑖 − ∇P + ν. ∇2 U (III. 8)
Dt ρ
On obtient alors, pour un écoulement incompressible, les équations de Navier-Stokes :
⃗
DU 1
=F ⃗⃗⃗⃗⃗⃗⃗⃗⃗ p + ν. ΔU
⃗ 𝑖 − grad ⃗ (III. 9)
Dt 𝜌
Page 4
Chapitre V : Principes Fondamentaux de la convection thermique
Avec :
D ∂ ∂ ∂ ∂
= + u + v + w : Dérivée particulaire
Dt ∂t ∂x ∂y ∂z
ρ: Densité
Fx, Fy, Fz : Forces extérieures
P : Pression,
μ: Viscosité dynamique,
u, v, w : Composantes du vecteur vitesse.
Equation de conservation de l’énergie
En appliquant le principe de conservation de l’énergie thermique déduit du premier
principe de la thermodynamique, l’équation de conservation de l’énergie s’exprime comme
suit :
DT dp
ρ. Cp . ⃗⃗⃗⃗⃗⃗⃗⃗⃗ T) + Tβ. + q̇ + ψ
= div(λ. grad (III. 11)
Dt dt
q̇ : sources de chaleur internes : rayonnement thermique, réactions chimiques…
ψ : la fonction de dissipation thermique
1 ∂ρ
β = − ρ (∂T) le taux de dilatation volumique à pression constante.
P
- Dans la plupart des cas usuels λ et μ peuvent considérés comme constants.
- Pour les liquides, en général
dp dT
≪ Cp
dt dt
1
- Pour les gaz parfaits : β = T
DT dp
ρ. Cp . ⃗⃗⃗⃗⃗⃗⃗⃗⃗ T) +
= div(λ. grad +q+ψ (III. 12)
Dt dt
DT
= a. ΔT (III. 13)
Dt
λ
Avec : a = ρ.C la diffusivité thermique.
P
DT ∂T ∂T ∂T ∂T
= +u +v + w : Dérivée particulaire de la temperature
Dt ∂t ∂x ∂y ∂z
Page 5
Chapitre V : Principes Fondamentaux de la convection thermique
Alors toutes les grandeurs physiques sont exprimées dans les dimensions
fondamentales suivantes:
• La masse (kg) M
• La longueur (m) L
• Le temps (s) T
• La température (K) θ
Page 6
Chapitre V : Principes Fondamentaux de la convection thermique
Chaque paramètre (πi ) sera formulé par des expressions potentielles des variables de base
Pour que π reste adimensionnel, il faut éliminer les dimensions, les exposants entre
parenthèses doivent avoir la valeur zéro, ce qui nous définit le système d'équations
linéaires suivant
b+d+e+g= 0
a + b + c − 3d − e + 2f = 0
{ (III. 16)
−3b − c − e − 2f − 3g = 0
−b − f − g + i = 0
Afin de déterminer les paramètres adimensionnels πi nous devons résoudre le système
d'équations linéaires (III.1). Pour la définition de chaque paramètre πi nous pouvons choisir
quatre coefficients de manière arbitraire (généralement on donne la valeur zéro) et calculer
les autres quatre coefficients par l'équation (V.1).
Nous obtenons par la méthode de l'analyse de dimension pour la convection thermique les
paramètres adimensionnels suivants:
Le nombre de Nusselt
g = 1 Pour obtenir une loi de la forme h = f ( . . .)
c = d = 0 Le groupement π trouvé ne dépendra pas de l'énergie cinétique du fluide ρ U2
i = 0 Le groupement π trouvé ne dépendra pas de l'écart de température Tp- T∞
Page 7
Chapitre V : Principes Fondamentaux de la convection thermique
b+d+e+g= 0 b + e = −1
{
a + b + c − 3d − e + 2f = 0
−3b − c − e − 2f − 3g = 0
{
a + b − e + 2f = 0
−3b − e − 2f − 3g = 3
−b − f − g + i = 0 −b − f = 1
⇒ a = 1, b = −1, e = 0, f = 0
π𝑖 = Da . λb . U c . ρd . μe . Cf . hg . (TP − T∞ )i
Avec : a = 1, b = −1, c = 0, d = 0, e = 0, f = 0, g = 1, i = 0
h. D
⇒ 𝜋1 = Nu = (III. 17)
λ
Le nombre de Reynolds
π𝑖 = Da . λb . U c . ρd . μe . C f . hg . (TP − T∞ )i
4 des 8 paramètres peuvent être choisis de manière arbitraire
b+d+e+g= 0 d+e=0
a + b + c − 3d − e + 2f = 0 a + c − 3d − e = 0
{ {
−3b − c − e − 2f − 3g = 0 −c − e = 0
−b − f − g + i = 0 i=0
ρ. U. D
⇒ 𝜋2 = Re = (III. 18)
μ
Nombre de Prandtl Pr
π𝑖 = Da . λb . U c . ρd . μe . Cf . hg . (TP − T∞ )i
Avec a=c=g=i=0
b+d+e+g= 0 b+d+e =0
{
a + b + c − 3d − e + 2f = 0
−3b − c − e − 2f − 3g = 0
b − 3d − e + 2f = 0
−3b − e − 2f = 0
−b − f − g + i = 0 −b − f = 0
Page 8
Chapitre V : Principes Fondamentaux de la convection thermique
μ. C
⇒ 𝜋3 = Pr = (III. 19)
λ
4.[p.(D22 −D21 ]
Espace annulaire : Dh = (p.D1 +p.D2 )
= D2 − D1
Page 9
Chapitre V : Principes Fondamentaux de la convection thermique
β. g. ΔT. ρ2 . D3
Gr =
μ2
β : dilatabilité du fluide en [K-1 ] ,
∆T : différence de température entre fluide et paroi : ∆T = Tparoi– Tfluide.
g. β. ΔT. D3
Ra = Pr. Gr =
a. ν
a = λ / (ρ .Cp) : diffusivité thermique [m²/s] ,
ν = µ / ρ : viscosité cinématique du fluide [m²/s].
⃗
DU 1
= − ⃗⃗⃗⃗⃗⃗⃗⃗⃗ ⃗⃗⃗
gradp + ν. ΔU (III. 21)
Dt ρ
T − Tref U x p t
T∗ = ; U∗ = ; x∗ = ; p∗ = 2 ; t∗ =
∆T Uref Lref ρ. Uref Lref /Uref
⃗∗
DU 1
⃗ p∗ + ΔU
= −∇ ⃗∗ (III. 22)
Dt ∗ Re
Uref . Lref
Re =
ν
La solution obtenue peut ensuite être ensuite introduite dans l´équation de T qui, en
écriture adimensionnelle, devient:
Page 10
Chapitre V : Principes Fondamentaux de la convection thermique
ρ. Cp . ΔT DT ∗ λ. ΔT ∗
. = 2 ΔT + q∗̇ (III. 23)
Lref /Uref Dt ∗ Lref
soit encore:
DT ∗ 1
= Δ𝑇 ∗ + q∗̇
Dt ∗ 𝑃𝑒
Uref .Lref
ou Pe désigne le nombre de Péclet: Pe = a
λ q̇ .𝐿𝑟𝑒𝑓
avec a = ρ.C diffusivité thermique du fluide. Le terme 𝑞 ∗ est tel que q̇ ∗ = 𝑈
P 𝑟𝑒𝑓 .Δ𝑇.𝜌.𝐶𝑝
L’écoulement du fluide est produit par des forces d’Archimède dues aux variations
locales de la masse volumique. Il n’y a pas de dispositif mécanique mettant le fluide en
mouvement. L’échelle de vitesse uo n’est plus imposée, mais dépend des autres échelles.
Les problèmes dynamique et thermique sont couplés et doivent être résolus simultanément.
L’utilisation des mêmes grandeurs adimensionnelles transforme l’équation de
Navier-Stokes ainsi:
⃗⃗ ∗
DU 1 g. β. ΔT. Lref ∗
⃗ p∗ +
= −∇ ⃗⃗ ∗ +
ΔU T (III. 24)
Dt ∗ Re 2
Uref
Il reste toutefois à déterminer Uref qui n’est pas imposée par un mécanisme externe mais
par la convection naturelle elle-même.
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:
2
ρ. Uref
≈ ρgβΔT ⇒ Uref = √gβΔTLref
Lref
de sorte que Navier-Stokes devient:
⃗∗
DU 1
⃗ p∗ + 1 ΔU
= −∇ ⃗ ∗ − ⃗⃗⃗
g∗. T ∗ (III. 25)
Dt ∗
𝐺𝑟 2
g. β. ΔT. L3ref
Gr =
ν2
L´équation de la température devient alors:
DT ∗ 1
∗
= 1/2
Δ𝑇 ∗ + q̇ ∗ (III. 26)
Dt 𝐺𝑟 . 𝑃𝑟
Page 11
Chapitre V : Principes Fondamentaux de la convection thermique
ν
Pr =
a
U
μ ≈ ρ. g. ΔT
L2
L’échelle de vitesse est alors:
g.β.ΔT.L2
Uref = ν
Page 12