31/03/2022

Chapitre 5

Initiation aux transferts de chaleur par convection

introduction

L’ étude des transferts de chaleur par convection passe obligatoirement par la résolution
des équations de transport ou encore équations de conservation de :
 masse,
 quantité de mouvement,
 énergie.
associées à leurs conditions aux limites.
Dans ce chapitre on se contente d’introduire le phénomène de convection et de se
focaliser sur la détermination du coefficient d’échange convectif noté h.

1
 31/03/2022

1- Généralités et définitions
1-1 Types de convection

La convection est un phénomène de conductibilité se produisant au sein

d’un fluide en mouvement. (voir chapitre 1)

Au phénomène de transmission par chocs entre molécules, s’ajoute le déplacement

d’énergie calorifique par les molécules se déplaçant suivant le champ des vitesses.

On distingue, deux types de convection selon le mécanisme (cause)

qui génère le mouvement de fluide :

 naturelle : la cause du mouvement est due aux gradients de température.

(Le fluide se met en mouvement sous l’effet des différences de masse
volumique résultant des gradients de température.)

 forcée : le mouvement de fluide est dû, par exemple, à une pompe, un ventilateur,…
le mouvement de fluide est donc indépendant du gradient de température.

NB

L’étude du transfert de chaleur par convection permet de déterminer les échanges

de chaleur se produisant entre un fluide et une paroi.

2
 31/03/2022

Remarque :

La quantité de chaleur échangée par unité de temps dépend de plusieurs

paramètres :

- la différence de température entre la paroi et le fluide (T);

- la vitesse du fluide (V);
- la capacité thermique massique du fluide (c);
- la surface d'échange (S);
- l'état de la surface de solide (lisse, rugueux);
- sa dimension (L, r,…) ;
- etc. . . .

1-2 Régimes d’écoulement

Les transferts de masse et de chaleur sont fortement liés, la considération de la nature

de régime d’écoulement lors d’une étude de convection thermique est inévitable.

A titre d’exemple, considérons l’écoulement d’un fluide dans une conduite.

On distingue deux régimes d’écoulement :

 Ecoulement en régime laminaire

 Ecoulement en régime turbulent

3
 31/03/2022

1-2-1 Ecoulement en régime laminaire

L’écoulement s’effectue par couches de fluide pratiquement indépendantes.

Entre deux filets fluides les échanges de chaleur s’effectuent:

 par conduction uniquement si l’on considère une direction normale aux filets fluides ;

 par convection et conduction (négligeable) si l ’on considère une direction

non normale aux filets fluides.

1.2.2 Ecoulement en régime turbulent

L’écoulement n’est pas unidirectionnel (il est aléatoire).

L’échange thermique dans la zone turbulente s’effectue simultanément par

convection et conduction dans toutes les directions.
On vérifie que la conduction est généralement négligeable par rapport à la convection.

4
 31/03/2022

2- Expression du flux de chaleur

2.1 Loi de Newton
Soit une paroi solide, où on considère un élément de surface ‘’dS’’ porté à une
température Tp, en contact avec un fluide de température T loin de la paroi.

Quelque soit le type de convection (libre ou forcée);

Quelque soit le régime d écoulement de fluide (laminaire ou turbulent),

le flux de chaleur élémentaire qui traverse dS est donné par la relation de Newton:

La quantité de chaleur Q qui travers dS pendant l’intervalle du temps, s’écrit:

h :est le coefficient d’échange par convection, il s’exprime en (W/m2.K)

Q : s’exprime en Joule

Q /dt : en Watt

NB :
Le problème majeur à résoudre avant le calcul du flux de chaleur consiste à déterminer
le coefficient d’échange par convection h, qui dépend de nombreux paramètres :

 caractéristiques de fluide,
 nature de l’écoulement,
 température,
 surface d’échange,...

5
 31/03/2022

2.2 Ordre de grandeur du coefficient convectif pour différentes configurations

Type de convection et configuration géométrique h (W.m-2.K-1)

Convection naturelle:
Plaque verticale de hauteur 0,3 m dans l’air. 4,5
Cylindre horizontal de diamètre 5 cm dans l’air. 6,5
Cylindre horizontal de diamètre 5 cm dans l’eau. 890

Convection forcée:
Courant d’air à 2m/s sur plaque carrée de 2m de coté 12
Courant d’air à 35m/s sur plaque carrée de 0,75m de coté 75
Eau à 0,5 kg/s dans un tube de diamètre 2,5 cm. 3500
Courant d’air à 50m/s perpendiculaire/tube de 5 cm de diamètre 180

Ébullition de l’eau:
Dans un récipient. 2500-35000
En écoulement dans un tube. 5000-10000

Ordres de grandeur du coefficient h (W.m-2.K-1)

Convection libre (air) 5 - 25

Convection libre (eau) 100 - 900
Convection forcée (air) 10 - 500
Convection forcée (eau) 100 - 15 000
Convection forcée (huile) 50 - 2 000
Convection forcée (métaux fondus) 6 000 - 120 000
Convection forcée (eau bouillante) 2 500 - 25 000
Condensation de vapeur d'eau 50 000 - 100 000

6
 31/03/2022

3- Analyse dimensionnelle
L'analyse dimensionnelle est un outil théorique servant à interpréter les problèmes à
partir des dimensions des grandeurs physiques mises en jeu.

Ses objectifs principaux sont :

 Préciser les relations entre les grandeurs mesurables intervenant dans un

phénomène physique.

 Regrouper les grandeurs en produits sans dimension afin d’exprimer les résultats
d’expériences différentes par des relations entre ces produits.

 Etude des similitudes.

3.1 Unités et grandeurs fondamentales

En mécanique des fluides et transferts thermiques par convection, les grandeurs

fondamentales sont :

 Pour une longueur, la dimension est : [Longueur] = L ;

 Pour une masse, la dimension est: [Masse] = M ;
 Pour un temps, la dimension est : [Temps] = T ;
 Pour une température, la dimension est : [Température] = .

Remarques :

 En fonction de L, M, T et , les unités dérivées sont constituées.

 Les grandeurs fondamentales de tout système sont indépendantes les une des autres.
 Le passage d'un système d'unités à un autre n'entraîne que des multiplicateurs
de conversion.

7
 31/03/2022

3.2 Homogénéité dimensionnelle

L’homogénéité dimensionnelle est illustrée par l’égalité des dimensions des termes
physiques des deux côtés d’une égalité. Elle représente la condition de départ de
l’analyse dimensionnelle.

Exemple : Force = Masse x Accélération (ou F = m a)

soit

3.3 Procédure à suivre dans un problème d’analyse dimensionnelle

o Identifier toutes les n variables (dépendants+indépendantes) intervenant dans

le problème étudié.

o Spécifier les dimensions de ces variables en utilisant les dimensions de base (L, M, T et ).

o Choisir les grandeurs fondamentales convenables, soit au nombre m.

o Utiliser une méthode appropriée pour identifier le nombre et la forme des paramètres

sans dimensions (paramètres adimensionnels).

8
 31/03/2022

3.4 Méthodes d’analyse dimensionnelle

Notons qu’il existe deux méthodes d’analyse dimensionnelle :

 la méthode de Rayleigh ;
 le théorème des , ou théorème de Vaschy-Buckingham.

Si le nombre de variables d’un système physique est connu, il est possible d’utiliser l’une
des deux méthodes suivantes pour déduire une relation représentative du phénomène
constituée de ces variables.

3.4.1 Méthode de Rayleigh

C’est une expression utilisée dans le cas où on a trois à quatre variables, dés qu’il y a
plus de variables indépendantes, il sera difficile de trouver une expression pour ce nombre
de variables.

Si la grandeur à exprimer est donnée par:

X=f(X1,X2,X3)

on peut la réécrire sous la forme :

k est une constante et a, b, c sont des puissances réelles à déterminer.

9
 31/03/2022

Exemple 1: Perte de charge dans une conduite

Dans le cas d’une conduite lisse, les variables indépendantes intervenant dans l’expression
de la variation de pression sur une longueur L sont au nombre de quatre (N=4): D (diamètre),
 (viscosité dynamique), V (vitesse moyenne) et  (masse volumique) :

(Eg1)

Les dimensions de base sont au nombre de trois (M, L, T) et ainsi les dimensions des
différentes grandeurs sont :

[P] = M.L-1.T-2, [D] = L, [L] = L,

[] = M.L-1.T-1, [V] = L.T-1, [] = M.L-3

Le membre de gauche de l’égalité (Eg1) et celui de droite doivent avoir la même

dimension. Par conséquent,

M.L-2.T-2 = [L]a [M.L-1.T-1]b [L.T-1]c [M.L-3]d

alors: 1 = b +d ; -2 = a – b + c -3d ; - 2 = -b – c.

On a donc un système à trois équations linéaires et quatre variables inconnues qui

conduit aux solutions :

a = - 1 –b, c = 2 – b ;d =1-b

10
 31/03/2022

La relation Eg1 s’écrit alors :

ou encore

Avec nombre adimensionnel appelé nombre de Reynolds.

Dans le cas de l’écoulement laminaire :

3.4.2 Méthode des π de Vaschy-Buckingham

[Link] Principe de la méthode

La méthode de Rayleigh devient très difficile dés que le nombre de variables dépasse
le nombre de dimensions fondamentales (M, L, T, ). Cette difficulté n’apparait pas dans
la méthode des π dont le principe est présenté dans le paragraphe ci-dessous.

Soit un système composé de ‘’n’’ variables (dépendantes et indépendantes) et

contenant ‘’m’’ dimensions fondamentales (M, L, T, ).

Soit X1, la variable dépendante telle que

X1=f(X2, X3, X4,.., Xn) (Eg2)

ou encore
g(X1,X2, X3, X4,.., Xn)=0

11
 31/03/2022

Si l’équation précédente est dimensionnellement homogène, elle peut être arrangée en

une équation de ‘’ n-m’’ termes adimensionnels.
En effet, si les variables répétées sont X2, X3, X4, on peut écrire :

…….

Chaque équation est résolue par le principe d’homogénéité ensuite on détermine

l’ensemble des π.

Conclusion :

Si 1, 2, 3, …n-m, représentent les nombres adimensionnels issus des quantités X1, X2, …Xn,
alors on peut écrire une équation de la forme :

12
 31/03/2022

[Link] Critères du choix des variables répétées

Le choix des variables répétées doit respecter les conditions suivantes :

 Dans la mesure du possible ; les variables dépendantes ne doivent pas être

utilisées comme variables répétées.

 Chaque variable répétée doit décrire un aspect du système physique :

 propriétés géométriques (longueur, diamètre,...etc.),
 propriétés physiques du fluide (masse volumique, viscosité…etc.) ;
 propriétés cinématique et/ou dynamique de l’écoulement
(vitesse, pression,…etc.).

 Les variables répétées ne doivent pas former un groupe adimensionnel.

 Ensemble les variables répétées doivent avoir le même nombre ‘’m’’ de

dimensions fondamentales.

 Dans la mesure du possible ; deux variables répétées ne doivent pas avoir

la même dimension.

[Link] Les étapes de l’application de la méthode

Pour appliquer la méthode de Vascy-Buckingham, les étapes principales suivantes sont

considérées :

1) Dresser une liste de toutes les quantités physiques Xi (1  i  n). Négliger toute
quantité qui est une fonction directe d’une ou d’autres quantités physiques ( exemple :
Débit et vitesse)

2) Ecrire la fonction : X1=f(X2, X3, X4,.., Xn).

3) Choisir les variables répétées.

4) Ecrire les paramètres i en fonction des exposants inconnus et s’assurer que

toutes les variables Xi sont incluses dans les i.

5) En se basant sur l’homogénéité dimensionnelle, écrire les systèmes d’équations

linéaires dont les inconnus sont les exposants.

6) Résoudre les équations simultanément.

7) Détermination des différents nombres adimensionnels i.

13
 31/03/2022

4- Analyse dimensionnelle appliquée à la convection

4.1 Paramètres influençant sur le coefficient de convection
Les résultats expérimentaux ont montré que le coefficient de convection h dépend,
en général, de 9 grandeurs : géométrique, physiques et cinématique. En effet,

(4-1)

le paramètre  est le coefficient de dilatation à pression constante.

Remarque :

 g (accélération gravitationnelle), T et  n’ont d’influence sur le coefficient d’échange

h que dans le cas de la convection naturelle.

 L’expression (4-1) s’exprime sous la forme de la relation (Eg2) définie ci-dessus.

Les dimensions fondamentales des grandeurs de la relation (4-1) sont regroupées

dans le tableau suivant:

Unités
Grandeurs M L T  Dimensions
(variables)
L 0 +1 0 0 M0.L+1.T0. 0
V 0 +1 -1 0 M0.L+1.T-1. 0
T 0 0 0 +1 M0.L0.T0. +1
 +1 -3 0 0 M+1.L-3.T0. 0
Cp 0 +2 -2 -1 M0.L+2.T-2. -1
 +1 +1 -3 -1 M+1.L+1.T-3. -1
 +1 -1 -1 0 M+1.L-1.T-1. 0
xg 0 +1 -2 -1 M0.L+1.T-2. -1
h +1 0 -3 -1 M+1.L0.T-3. -1

Remarque : Ce tableau forme une matrice du rang m=4.

14
 31/03/2022

4.2 Détermination de h par la méthode de l'analyse dimensionnelle

Pour déterminer le coefficient d’échange convectif h, nous optons pour la méthode

de Vaschy-Buckingham. En effet,

n = 9, le nombre des variables (dépendantes et indépendantes),

m = 4, le nombre des dimension de ces variables (M, L, T et ),
ainsi, nous avons :
n - m = 5 nombres sans dimensions.

Ces cinq (05) nombres i ( 1  i  5) peuvent être déterminés suivants :

D’autre part,

avec ai, bi, ci et di (1  i  5) sont des exposants à déterminer.

15
 31/03/2022

Les variables répétées sont : L, V,  et Cp. Elles ont été choisies selon les critères du choix
mentionnés ci-dessus.

D’un point de vue dimensionnel, les relations précédentes s’écrivent :

A l’aide du principe d’homogénéité des unités, les relations ci-dessus se transforment en

des systèmes d’équations linéaires dont la résolution détermine l’ensemble des i ( i = 1 à 5).

i- Détermination de 1

Soit

Ce qui se traduit par un système de 4 équations linéaires à quatre valeurs inconnues

(a1, b1, c1 et d1). La résolution de ce système donne : (a1 = 1, b1 = 0, c1 = -1 et d1 = 0).

Ainsi, nous avons :

1 est appelé nombre de Nusselt et on le noté :

16
 31/03/2022

ii-Expressions des 2, 3, 4 et 5

Par un raisonnement analogue, les expressions des autres nombres adimensionnels

ont été trouvées. Ainsi

Notons que :

où Pr, Re, Pe=RexPr, Ec, Gr sont respectivement les nombres de :

Prandtl, Reynolds, Peclet, Eckert et Grashof.

Conclusion :

Selon l’analyse dimensionnelle, le coefficient d’échange convectif h représenté par

le nombre de Nusselt dépend, en général, d’autres nombres adimensionnels :

Cette dépendance est fortement liée au type de convection (forcée ou libre).

En pratique :

 Dans le cas de la convection forcée :

 Dans le cas de la convection libre :

Remarque :

Notons que le nombre d’Eckert Ec n’intervient que dans les écoulements dont la vitesse
est proche de celle du son.

17
 31/03/2022

4.3 Signification physique des principaux nombres adimensionnels

4.3.1 Nombre de Nusselt

C’est le rapport de la quantité de chaleur échangée par convection h.S.T à

une quantité de chaleur échangée par conduction .S.T/d

Remarque importante:

Nu est fonction directe du coefficient d’échange convectif h, sa connaissance permet

de déterminer la valeur de h.

4.3.2 Nombre de Grashof

Il compare les forces de gravité par unité de volume et celles de viscosité
par unité de volume :

4.3.3 Nombre d’Eckert

Il Caractérise la dissipation d'énergie par frottement au sein du fluide
(dégradation de l’énergie mécanique en chaleur) :

4.3.4 Nombre de Prandtl

Il caractérise la distribution des vitesses par rapport à la distribution de la température :

4.3.5 Nombre de Reynolds

Il compare les forces d’inertie aux forces de viscosité. Il caractérise également la forme
du profil de vitesse de l’écoulement de fluide :

18
 31/03/2022

4.4 Méthode pratique de calcul de h

4.4.1 Quelques précisions nécessaires

Avant de procéder au calcul de h il faut connaître :

1- la nature de fluide (liquide ou gaz) ;

2- l'intervalle de température de fluide ;

3- le type de convection (naturelle ou forcée) ;

4- le régime d'écoulement (laminaire ou turbulent) ;

5- la nature d’écoulement (interne ou externe) ;

6- la géométrie de la paroi en contact avec fluide ( surface plane, circulaire,…).

4.4.2 Différentes étapes de calcul

i) Dans le cas de la convection forcée :

 Calculer Re et le comparer à Rec (nombre de Reynolds critique approprié) ;

 Si Re < Rec le régime est dit laminaire sinon (Re > Rec) le régime est dit turbulent;
 Utiliser l'une des formules empiriques données pour déterminer h ;
 Calcule de  Q par la formule de Newton et intégrer pour avoir Q

ii) Dans le cas de la convection naturelle :

 Calculer Gr (ou Ra=Gr Pr) et le comparer à Grc (ou Rac) (nombre de

Grashof critique ou nombre de Rayleigh critique approprié) ;
 Si Gr < Grc le régime est dit laminaire sinon (Gr > Grc) le régime est
dit turbulent;
 Utiliser l'une des formules empiriques données pour déterminer h ;
 Calcule de  Q par la formule de Newton et intégrer pour avoir Q

19
 31/03/2022

4.4.3 Expressions empiriques du nombre de Nusselt

Généralement, les nombres de Nusselt moyens utilisés sont donnés par des corrélations
de la forme :

 Dans le cas de la convection forcée :

 Dans le cas de la convection naturelle :

Notons que Ra = Gr Pr est appelé nombre de Rayleigh.

Les constantes C et les exposants n et m dépendent généralement du régime

(laminaire ou turbulent) et du type (interne ou externe) d’écoulement.

4-5 Corrélations (Expressions) du coefficient de convection h

Quelques corrélations pour le calcul des coefficients d’échange par convection sont
données en annexe.

20