Cours de Transfert Thermique I 2020 2021-1

Institut National Supérieur des Sciences et Techniques d’Abéche
INSTA
Tronc Commun : Génie Electrique, Génie Energétique et Génie Mécanique
Transfert Thermique I
Chargé de Cours et TD : M. DJALDI TABDI NGAMSOU

Ing. Energétique et Froid Industriel
Enseignant-Chercheur
Institut National Supérieur des Sciences et Techniques d’Abéché - INSTA
Table des matières

Liste de figures ................................................................................................................................................. 4
Liste de tableaux ............................................................................................................................................... 5
Chapitre1 Introduction au Transfert Thermique ................................................................................................ 6
1. Thermodynamique et transfert thermique .................................................................................................. 6
2. Les différents modes fondamentaux de transfert thermique ........................................................................ 6
2.1. Conduction ...................................................................................................................................... 6
2.2. Convection ....................................................................................................................................... 6
2.3. Rayonnement ................................................................................................................................... 6
3. Définitions et notions fondamentales ......................................................................................................... 7
3.1. Température ..................................................................................................................................... 7
3.2. Champ de Température .................................................................................................................... 7
3.3. Surface Isotherme ............................................................................................................................ 7
3.4. Gradient de Température .................................................................................................................. 7
3.5. Quantité de Chaleur ou Energie-chaleur ............................................................................................ 8
3.6. Flux de Chaleur ................................................................................................................................ 8
3.7. Densité de flux de chaleur ................................................................................................................ 8
4. Bilan thermique ........................................................................................................................................ 8
4.1. Etablissement d’un bilan thermique .................................................................................................. 8
4.2. Expression des flux de chaleur .......................................................................................................... 9
Chapitre 2 Transfert de Chaleur par Conduction............................................................................................. 12
1. Généralités............................................................................................................................................. 12
2. Notions fondamentales............................................................................................................................ 12
2.1. Loi de Fourier ................................................................................................................................ 12
2.2. Equation différentielle de la conduction .......................................................................................... 13
3. Conduction stationnaire unidimensionnelle ............................................................................................. 14
3.1. Corps homogènes ........................................................................................................................... 14
3.2. Corps ayant des structures multicouches ......................................................................................... 21
Chapitre 3 Transfert de Chaleur par Convection ............................................................................................. 26
1. Analyse dimensionnelle .......................................................................................................................... 26
1.1. Dimensions fondamentales ............................................................................................................. 26
1.2. Principe de la méthode ................................................................................................................... 26
1.3. Exemple d’application .................................................................................................................... 27
1.4. Avantages de l’utilisation des grandeurs réduites ............................................................................ 29
2. Généralités et définitions ........................................................................................................................ 30
2.1. Convection naturelle et forcée ........................................................................................................ 30
2.2. Régime d’écoulement ..................................................................................................................... 30
3. Flux de chaleur transféré par convection ................................................................................................ 31
3.1. Analogie de Reynolds .................................................................................................................... 31
3.2. Couches limites dynamique et thermique ........................................................................................ 31
Transfert Thermique I Page 2

3.3. Expression du flux ......................................................................................................................... 32

4. Calcul du flux de chaleur en convection forcée........................................................................................ 32
5. Calcul du flux de chaleur en convection naturelle ................................................................................... 33
5.1. Mécanisme de la convection naturelle............................................................................................. 33
5.2. Signification physique du nombre de Grashof ................................................................................. 34
5.3. Calcul du flux de chaleur en convection naturelle............................................................................ 34
6. Corrélations utilisées en convection forcée ............................................................................................. 35
6.1. Plaques planes................................................................................................................................ 35
6.2. Transfert thermique en conduite en régime laminaire ...................................................................... 35
6.3. Transfert thermique en conduite en régime turbulent ....................................................................... 35
6.4. Convection forcée autour d’un cylindre à température uniforme et constante placé transversalement
par rapport à l’écoulement .......................................................................................................................... 36
6.5. Convection forcée autour d’une sphère solide à température uniforme et constante .......................... 36
6.6. Convection forcée à travers un faisceau de tubes cylindriques de section circulaire alignés par rapport
à l’écoulement ............................................................................................................................................ 36
6.7. Convection forcée à travers un faisceau de tubes cylindriques de section circulaire décalés par rapport
à l’écoulement ............................................................................................................................................ 37
7. Corrélations utilisées en convection naturelle dans les configurations les plus utilisées ........................... 37
7.1. Convection naturelle au-dessus d’une plaque chaude ou en-dessous d’une plaque froide horizontale 37
7.2. Convection naturelle au-dessus d’une plaque froide ou en-dessous d’une plaque chaude horizontale 37
7.3. Convection naturelle au-dessus d’un cylindre horizontal chaud d’axe horizontal.............................. 38
7.4. Convection naturelle au-dessus d’une sphère chaude....................................................................... 38
Chapitre 4 Transfert de Chaleur par Rayonnement .......................................................................................... 39
1. Introduction ........................................................................................................................................... 39
2. Généralités. Définitions .......................................................................................................................... 39
2.1. Nature du rayonnement .................................................................................................................. 39
2.2. Définitions ..................................................................................................................................... 40
3. Lois du rayonnement .............................................................................................................................. 43
3.1. Loi de Lambert............................................................................................................................... 43
3.2. Lois physiques ............................................................................................................................... 44
4. Rayonnement réciproque de plusieurs surfaces ....................................................................................... 46
4.1. Radiosité et flux net perdu .............................................................................................................. 46
4.2. Facteur de forme géométrique ........................................................................................................ 47
4.3. Calcul des flux ............................................................................................................................... 48
4.4. Méthode de résolution .................................................................................................................... 48
4.5. Exemple d’application : Cas de deux plans parallèles infinis ........................................................... 49
4.6. Analogie électrique ........................................................................................................................ 49
Série d’exercices ............................................................................................................................................. 52
Annexes .......................................................................................................................................................... 58
A.1.1 : Fraction d’énergie 𝑭𝟎 − 𝝀𝑻 rayonnée par un corps noir entre 𝟎 et 𝝀 ................................................. 58

A.1.2 : Facteurs de forme géométrique de rayonnement ............................................................................... 59

Bibliographie .................................................................................................................................................. 62
Liste de figures
Figure 1 : Isotherme et gradient de température ................................................................................................. 8
Figure 2 : Système et bilan thermique................................................................................................................ 9
Figure 3 : Schéma du transfert de chaleur par conduction ................................................................................... 9
Figure 4 : Schéma du transfert de chaleur par convection ................................................................................. 10
Figure 5 : Schéma du transfert de chaleur par rayonnement.............................................................................. 10
Figure 6 : Positions relatives du flux surfacique et du gradient de température.................................................. 12
Figure 7 : Signification du flux surfacique ....................................................................................................... 13
Figure 8 : Conduction thermique stationnaire dans un mur plan de conductivité thermique constante................ 14
Figure 9 : Analogie thermoélectrique dans le cas d’un mur plan simple ............................................................ 16
Figure 10 : Conduction thermique stationnaire dans un mur simple en contact avec deux fluides ...................... 16
Figure 11 : Schéma électrique équivalent au transfert thermique entre deux fluides séparés par un mur plan
simple............................................................................................................................................................. 17
Figure 12 : Cylindre creux à surfaces latérales isothermes ................................................................................ 18
Figure 13 : Cylindre creux simple au contact avec deux fluides : distribution de température ............................ 19
Figure 14 : Schéma électrique équivalent d’un Cylindre creux simple au contact avec deux fluides .................. 20
Figure 15 : Mur composite ayant deux surfaces isothermes .............................................................................. 21
Figure 16 : Mur composite en contact avec deux fluides .................................................................................. 23
Figure 17 : Cylindre creux composite de surfaces isothermes ........................................................................... 24
Figure 18 : Schéma de la configuration étudiée ................................................................................................ 27
Figure 19 : Schématisation d’un écoulement laminaire .................................................................................... 30
Figure 20 : Schématisation d’un écoulement turbulent ..................................................................................... 31
Figure 21 : Représentation de l’analogie de Reynolds dans le cas d’un écoulement turbulent dans un tube ........ 31
Figure 22 : Représentation du modèle de Prandtl pour un écoulement turbulent dans une conduite ................... 31
Figure 23 : Représentation du mécanisme de convection naturelle ................................................................... 33
Figure 24 : Principe de l’expérience de William Herschel ................................................................................ 40
Figure 25 : Spectre des ondes électromagnétiques ............................................................................................ 40
Figure 26 : Schéma de l’angle solide ............................................................................................................... 41
Figure 27 : Schéma de définition des angles .................................................................................................... 42
Figure 28 : Schématisation de la répartition d’un flux incident de rayonnement sur un solide ........................... 42
Figure 29 : Représentation simplifiée du pouvoir absorbant monochromatique de la peinture blanche .............. 43
Figure 30 : Schématisation de l’intensité énergétique ....................................................................................... 44
Figure 31 : Schématisation de la luminance et de l’intensité énergétique d’une source isotrope ......................... 44
Figure 32 : Emittance monochromatique d’un corps noir à deux températures différentes ................................. 45
Figure 33 : Schématisation des flux de rayonnement sur une surface ................................................................ 46
Figure 34 : Schéma électrique équivalent du flux radiatif perdu par une surface ............................................... 49
Figure 35 : Schéma électrique équivalent du flux radiatif échangé entre deux surfaces ..................................... 50

Figure 36 : Schéma électrique équivalent du flux radiatif net échangé entre deux surfaces ................................ 50
Liste de tableaux
Tableau 1 : Grandeurs physiques et leurs dimensions (unités) .......................................................................... 27
Tableau 2 : Quelques groupements sans dimensions ........................................................................................ 30
Tableau 3 : Ordre de grandeur du coefficient de transfert de chaleur par convection ......................................... 32

Chapitre1 Introduction au Transfert Thermique

1. Thermodynamique et transfert thermique
La thermodynamique évalue la quantité de chaleur échangée avec l’extérieur pendant un processus subi par un
système passant d’un état d’équilibre à un autre état d’équilibre.
Le transfert thermique (ou de chaleur) appelé aussi thermocinétique, ou thermique, se propose de décrire
quantitativement (dans l’espace et dans le temps) l’évolution des grandeurs caractéristiques du système, en
particulier la température, entre l’état d’équilibre initial et l’état d’équilibre final.
La notion de Chaleur est définie par la thermodynamique comme étant une grandeur de parcours représentant
une forme transitoire d’énergie qui se propage à la suite d’une différence de température soit :
 entre deux parties d’un système quelconque ;
 entre le système et son environnement ;
 entre deux systèmes séparés par un milieu matériel ou pas (vide).
Le transfert thermique traite la description des mécanismes des modes de propagation de la chaleur.
Les méthodes d’étude de la thermodynamique sont valables pour le transfert thermique.
2. Les différents modes fondamentaux de transfert thermique
Le transfert thermique est un processus complexe qui est réalisé par la superposition des trois modes
fondamentaux qui sont : la conduction, la convection et le rayonnement.
2.1. Conduction
La conduction est le processus de propagation de la chaleur au sein d’un milieu opaque, sans déplacement de
matière, par le contact direct entre les particules d’un corps ou entre corps sous l’influence d’une différence de
température, suite au mouvement de ces particules élémentaires.
La propagation de la chaleur par conduction à l’intérieur d’un corps s’effectue selon deux mécanismes distincts :
une transmission par les vibrations des atomes ou molécules et une transmission par les électrons libres.
2.2. Convection
La convection est le processus de transfert de chaleur entre un solide et un fluide déterminé par le mouvement
des particules élémentaires du fluide (liquide ou gaz) entre des zones ayant des températures différentes.
L’énergie étant transmise par le déplacement du fluide.
Compte tenu des forces qui produisent le mouvement du fluide, on a deux types de convection : forcée et
naturelle ou libre.
 Dans la convection forcée, le mouvement de fluide est causé par l’action des forces extérieures du
processus (par exemple pompe, ventilateur, etc.).
 La convection naturelle a comme origine le mouvement produit par les différences de densité entre les
particules chaudes et celles qui sont froides existant dans un fluide situé dans un champ de forces
massiques (par exemple champ gravitationnel).
2.3. Rayonnement
Le rayonnement thermique constitue une forme particulière de transfert de chaleur dans laquelle le porteur
d’énergie n’est plus représenté par des particules de substance, mais par des ondes électromagnétiques ayant à
l’origine l’agitation atomique à la surface d’un corps. Le mécanisme de rayonnement est créé par l’émission et
l’absorption des ondes électromagnétiques porteuses d’énergie rayonnante et par la transformation de celle-ci en
chaleur avec le changement d’état énergétique (énergie interne) des corps qui la retiennent.

Le rayonnement thermique n’exigeant pas un support matériel, il peut se produire même dans le vide.
3. Définitions et notions fondamentales
3.1. Température
La température représente un paramètre d’état thermodynamique intensif d’un système, traduisant l’état
énergétique de la matière.
Pratiquement, elle exprime la qualité de chaud ou de froid qu’on peut attribuer à un corps quelconque solide,
liquide ou gaz.
Il existe plusieurs échelles de températures ; on utilise l’échelle de températures absolues, ayant le symbole 𝑇,
dont l’unité est le Kelvin[𝐾]. Dans la pratique, on utilise le degré Celsius[°𝐶], sur l’échelle correspondante avec
le symbole 𝜃.
On retient que :
𝑇[𝐾] = 273,15 + 𝜃[°𝐶] (1.01)
3.2. Champ de Température
La propagation de la chaleur par tous les modes de transfert de chaleur a lieu seulement s’il y a une différence de
température entre les divers points d’un corps (ou d’un système des corps).
Les transferts de chaleur sont déterminés à partir de l’évolution dans l’espace et dans le temps de la température :
𝑇 = 𝑇(𝑥, 𝑦, 𝑧, 𝑡) (1.02)
La valeur instantanée de la température en tout point de l’espace est un scalaire appelé champ de température ou
distribution de température.
On distingue deux cas :
 Champ de température indépendant du temps : le régime est dit permanent ou stationnaire :
𝜕𝑇
𝑇 = 𝑇1 (𝑥, 𝑦, 𝑧), =0 (1.03)
𝜕𝑡
 Evolution du champ de température avec le temps : le régime est dit variable ou transitoire :
𝑇 = 𝑇(𝑥, 𝑦, 𝑧, 𝑡) (1.04)
3.3. Surface Isotherme
Une surface isotherme est l’ensemble de tous les points de l’espace qui ont la même température. Elle peut être :
 permanente ou stationnaire, si elle ne modifie pas ses coordonnées spatiales par rapport au temps, ou,
autrement dit, si elle est invariable ;
 transitoire ou variable, si elle modifie d’une manière continue sa position par rapport au temps, c’est-à-
dire si elle est mobile et déformable.
Le point singulier représente un point d’un espace matériel ayant une température qu’on ne trouve pas ailleurs.
3.4. Gradient de Température
La variation de température par unité de longueur est maximale le long de la normale à la surface isotherme.
Cette variation est caractérisée par le gradient de température qui est un vecteur normal à la commune de deux
courbes isothermes voisines, à l’intérieur d’un corps quelconque.
L’expression mathématique est donnée par l’application de l’opérateur ∇ (nabla) au champ de température :
𝜕𝑇
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ = ⃗⃗⃗⃗⃗
𝑔𝑟𝑎𝑑𝑇 𝛻𝑇 = 𝑛
⃗⃗⃗⃗ (1.05)
𝜕𝑛 0
où :
 𝑛0 : vecteur unitaire de la normale commune aux surfaces isothermes ;
⃗⃗⃗⃗

𝜕𝑇
 𝜕𝑛
⃗ .
: dérivée de la température le long de la normalen
Si on considère un système monodirectionnel suivant 0𝑥, on peut écrire :

𝜕𝑇
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ =
𝑔𝑟𝑎𝑑𝑇 𝑖 (1.06)
𝜕𝑥
Unité : [𝐾/𝑚].
Figure 1 : Isotherme et gradient de température

3.5. Quantité de Chaleur ou Energie-chaleur
La quantité de chaleur ou chaleur échangée par un système qui évolue entre deux états d’équilibre (1) et (2) et
son milieu extérieur est :
2 2
𝑄 = 𝑄12 = ∫ 𝛿𝑄 = ∫ 𝑚𝑐. 𝑑𝑇 = 𝑚𝑐̅(𝑇2 − 𝑇1 ) (1.07)
1 1
Avec :
 𝑄 : quantité de chaleur en Joules [𝐽]
 𝑚 mase du système en kilogramme [𝑘𝑔],
 (𝑇2 − 𝑇1 ) différence de température entre les deux états en Kelvin [𝐾],
 𝑐̅ chaleur massique spécifique moyenne [𝐽/𝑘𝑔. 𝐾].
3.6. Flux de Chaleur
Le flux de Chaleur noté 𝛷 représente la chaleur transférée à travers une surface donnée dans une unité de
temps :
𝑑𝑄 𝑄
𝛷= = (1.08)
𝑑𝑡 𝑡
Unité : [𝑊]
Le flux de chaleur ou flux thermique est souvent appelé puissance thermique, employée surtout dans la
terminologie des appareils et appareils thermiques.
3.7. Densité de flux de chaleur
Le flux thermique surfacique ou densité de flux thermique ou densité de flux de chaleur noté 𝜑 à travers une
surface 𝑆 représente le flux thermique total rapporté à l’unité de surface ou la chaleur transférée dans une unité
de temps de surface :
𝛷 1 𝑑𝑄
𝜑= = (1.09)
𝑆 𝑆 𝑑𝑆
Unité : [𝑊. 𝑚−2 ].
4. Bilan thermique
4.1. Etablissement d’un bilan thermique
Pour faire l’étude de transfert thermique, il faut définir un système (S) par ses limites dans l’espace et faire
l’inventaire des différentes densités de flux de chaleur qui influent sur l’état du système et qui peuvent être :

 φst : densité de flux de chaleur stockés dans le système (S) ;

 φg : densité de flux de chaleur générée dans le système (S) ;
 φe : densité de flux de chaleur entrant dans le système (S) ;

 φs : densité de flux de chaleur sortant du système (S).
On applique alors le premier principe de la thermodynamique pour établir le bilan d’énergie :
𝜑𝑠𝑡 + 𝜑𝑠 = 𝜑𝑒 + 𝜑𝑔 (1.10)
Figure 2 : Système et bilan thermique

4.2. Expression des flux de chaleur
Il faut ensuite établir les expressions des différents flux de chaleur. En reportant ces expressions dans le bilan
d’énergie, on obtient l’équation différentielle dont la résolution permet de connaître l’évolution de la température
en chaque point du système.
4.2.1. Conduction
La théorie de la conduction repose sur l’hypothèse de Fourier : la densité de flux de chaleur est proportionnelle
au gradient de température :
⃗ = −𝜆 . ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝜑 𝑔𝑟𝑎𝑑(𝑇) (1.11)
Sous forme algébrique et en coordonnées cartésiennes suivant la direction 0𝑥 la densité de flux de chaleur est :
𝜕𝑇
𝜑 = −𝜆. (1.12)
𝜕𝑥
Et le flux de chaleur correspondant s’écrit :
𝜕𝑇
𝛷 = −𝜆. 𝑆. (1.13)
𝜕𝑥
Avec
 𝜆 : conductivité thermique du milieu [𝑊. 𝑚−1 . 𝐾 −1 ] ;
 𝑥 : variable d’espace dans la direction de propagation du flux [𝑚] ;
 𝑆 : aire de la section de passage du flux de chaleur [𝑚2 ].
Figure 3 : Schéma du transfert de chaleur par conduction
4.2.2. Convection
Le mécanisme de transfert thermique par convection est régi par la loi de Newton :
𝜑 = ℎ(𝑇𝑝 − 𝑇∞ ) (1.14)

Figure 4 : Schéma du transfert de chaleur par convection

Avec :
 ℎ : coefficient de transfert de chaleur par convection [𝑊. 𝑚−2 . 𝐾 −1 ] ;
 𝑇𝑝 : température de surface du solide [𝐾] ;
 𝑇∞ : température du fluide loin de la surface [𝐾].

Remarque :
La valeur du coefficient de transfert de chaleur par convection ℎ est fonction de la nature du fluide, de sa
température, de sa vitesse et des caractéristiques géométriques de la surface de contact solide/fluide.
4.2.3. Rayonnement
Dans les problèmes de conduction, on prend en compte le rayonnement entre un solide et le milieu environnant
et dans ce cas nous avons la relation :
𝛷 = 𝜎. 𝜀𝑝 . 𝑆. (𝑇𝑆4 − 𝑇∞4 ) (1.15)
Figure 5 : Schéma du transfert de chaleur par rayonnement

Avec :
 𝜎 : constante de STEPHAN (σ = 5,67.10-8 𝑊. 𝑚−2 . 𝐾 −4 ) ;
 𝜀𝑝 : facteur d’émission de la surface ;
 𝑇𝑆 : température de surface [𝐾] ;

 𝑇∞ : température du milieu environnant la surface [𝐾] ;
4.2.4. Flux de chaleur lié à un débit massique
Lorsqu’un débit massique 𝑚̇ de matière entre dans le système à la température 𝑇1 et en ressort à la
température 𝑇2 , on doit considérer dans le bilan (1.10) un flux de chaleur entrant correspondant :
𝛷𝑒 = 𝑚̇. 𝑐𝑝 (𝑇1 − 𝑇2 ) (1.16)
Avec :
 𝑚̇ : Débit massique [𝑘𝑔. 𝑠 −1 ] ;
 𝑐𝑝 : Chaleur spécifique [𝐽. 𝑘𝑔 −1 𝐾 −1 ],
 𝑇1 et 𝑇2 : Températures d’entrée et de sortie [𝐾].

4.2.5. Stockage d’énergie dans un solide
Le stockage d’énergie thermique dans un corps (solide) correspond à une augmentation de son énergie interne au
cours du temps (à pression constant), d’où :

𝜕𝑇
𝛷𝑠𝑡 = 𝜌. 𝑉. 𝐶. (1.17)
𝜕𝑡
Avec :
 𝜌 : masse volumique [𝑘𝑔. 𝑚−3] ;
 𝑉 : volume [𝑚3] ;
 𝐶 : chaleur massique [𝐽. 𝑘𝑔 −1 . 𝐾 −1] ;
 𝑇 : température [K] ;
 𝑡 : temps [s].
Les paramètres 𝜌, 𝑉, 𝐶 sont supposés constants. Le produit 𝜌. 𝑉. 𝐶 est appelé la capacitance thermique du corps.
4.2.6. Génération d’énergie
Elle intervient lorsqu’une autre forme d’énergie (chimique, électrique, mécanique, nucléaire, …) est convertie en
énergie thermique :
𝛷𝑔 = 𝑞̇ . 𝑉 (1.18)
Avec :
 𝑞̇ = 𝜔 = 𝑠 : densité volumique d’énergie libérée ou fonction de puissance volumique ou source interne
de chaleur [𝑊. 𝑚−3],
 𝑉 : volume [𝑚3].

Chapitre 2 Transfert de Chaleur par Conduction

1. Généralités
La conduction thermique est un phénomène de propagation d’énergie, qui traduit à l’échelle macroscopique les
échanges directs d’énergie à l’échelle microscopique entre les molécules ou les atomes (ou encore leurs
électrons) à l’intérieur d’un matériau. Ces échanges résultent des mouvements des particules constituant le
système. A l’échelle macroscopique, la conduction entre deux points d’un milieu s’établit lorsque ces deux
points ont une température différente.
Par exemple, considérons un barreau chauffé à une extrémité et initialement à la température ambiante.
L’énergie reçue à une extrémité du barreau se propage le long du barreau, ce qui se traduit par une élévation
progressive de la température du barreau.
Le mécanisme de base qui contrôle le transfert thermique par conduction est différent selon la structure physique
du système : solide, liquide ou gaz.
D’une manière générale, la théorie analytique de la conduction repose sur le schéma du milieu continu, ne tenant
pas compte de la structure moléculaire des systèmes concernés. Cette hypothèse est justifiée pour les cas des
systèmes qui ont des dimensions beaucoup plus grandes que celles de leurs molécules et des distances
intermoléculaires.
Le mécanisme de base qui contrôle le transfert thermique par conduction est différent selon la structure physique
du corps : solide, liquide ou gaz.
Comme on sait que la propagation de la chaleur dans un milieu donné a comme origine la distribution non
uniforme de température, le flux thermique transféré par conduction nécessite que le gradient de température soit
différent de zéro entre les différents points du corps.
L’étude du transfert thermique par conduction dans un milieu a comme principaux objectifs la connaissance de
la distribution de température ou champ de température, dans ce milieu, ainsi que le flux thermique qui le
traverse.
2. Notions fondamentales
2.1. Loi de Fourier
La loi fondamentale de la conduction se repose sur la loi de Fourier. Elle s’énonce ainsi : la densité de flux de
chaleur est proportionnelle au gradient de température :
𝜕𝑇
𝜑 ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ = −𝜆𝛻𝑇
⃗ = −𝜆 . 𝑔𝑟𝑎𝑑𝑇 ⃗⃗⃗⃗⃗ = −𝜆
𝑛
⃗⃗⃗⃗ (2.01)
𝜕𝑛 0
On constate donc que la densité de flux de chaleur est un vecteur colinéaire au gradient de température, de
direction perpendiculaire aux surfaces isothermes (fig. 2.1).
Figure 6 : Positions relatives du flux surfacique et du gradient de température

Le facteur de proportionnalité 𝜆 qui apparaît dans la relation (2.01) est obtenu sur la voie expérimentale et
constitue une propriété physique du corps. C’est un scalaire positif, qui définit l’aptitude de celui-ci à transmettre
la chaleur et est appelé conductivité thermique.
La grandeur scalaire qui représente la valeur de la densité de flux de chaleur est donnée par :
𝜕𝑇
𝜑 = −𝜆. [𝑊/𝑚2 ] (2.02)
𝜕𝑛
En utilisant la loi de Fourier, on peut redéfinir les notions de flux thermique et d’énergie-chaleur d’une manière
plus adaptée aux transferts thermiques. Ainsi, le flux thermique ou la puissance thermique représente l’énergie-
chaleur qui passe par unité de temps à travers une surface isotherme :
𝜕𝑇
𝛷 = ∫ 𝜑. 𝑑𝑆 = −𝜆 ∫ 𝜆. . 𝑑𝑆 [𝑊] (2.03)
𝜕𝑛
𝑆 𝑆
Où 𝑑𝑆 est un élément de la surface isotherme (figure 7).
Figure 7 : Signification du flux surfacique

L’énergie-chaleur totale qui traverse une surface isotherme 𝑆 dans le temps 𝑡 est donnée par l’expression :
𝑡 𝑡
𝜕𝑇
𝑄 = ∫ 𝛷. 𝑑𝑡 = − ∫ ∫ 𝜆. . 𝑑𝑆. 𝑑𝑡 [𝐽] (2.04)
𝜕𝑛
0 0 𝑆
Dans le cas concret d’un corps étudié, il est évident que, pour déterminer toutes les grandeurs impliquant le
gradient de température, il faut connaître le champ de température dans ce corps.
2.2. Equation différentielle de la conduction
Pour déterminer le champ de température dans un corps qui subit un transfert thermique par conduction, il existe
une relation entre les propriétés caractéristiques du corps considéré et les paramètres impliqués dans le
phénomène. Cette relation est exprimée par une équation différentielle :
𝜕𝑇 𝑠 𝑠
= 𝑎. ∇2 𝑇 + = 𝑎. ∆𝑇 + (2.05)
𝜕𝑡 𝜌𝑐𝑝 𝜌𝑐𝑝
Avec :
 𝑎 : diffusivité thermique en [𝑚2 /𝑠],
 𝑠 = 𝑞̇ = 𝜔 : densité volumique d’énergie libérée ou fonction de puissance volumique ou sources
internes de chaleur [𝑊. 𝑚−3 ].
 𝑐𝑝 : chaleur spécifique du système en [𝐽/𝑘𝑔°𝐶] ou [𝐽/𝑘𝑔𝐾]
 𝜌 : masse volumique en [𝑘𝑔/𝑚3]

Cette forme exprime l’équation aux dérivées partielles de la conduction thermique (ou de la chaleur). Elle régit
les variations spatiales et temporelles de la température en chaque point du corps où a lieu un transfert de chaleur
par conduction et les caractéristiques thermo-physiques de celui-ci. Compte tenu des hypothèses, elle est valable
pour un corps indéformable, immobile, homogène et isotrope.
Le facteur de proportionnalité « 𝑎 » s’appelle diffusivité thermique et représente une propriété physique du
corps. Elle caractérise l’aptitude des corps solides à laisser passer la chaleur par conduction ou, à diffuser la
température à l’intérieur de leur substance et s’exprime par [𝑚2 /𝑠].
3. Conduction stationnaire unidimensionnelle
L’équation différentielle de la conduction thermique en régime stationnaire correspond à une forme particulière
𝜕𝑇
de l’équation générale (2.05), dans laquelle le champ de température ne varie pas dans le temps ( 𝜕𝑡 = 0). Pour
les corps ayant des sources internes de chaleur uniformément distribués (𝑠 = 𝑐𝑠𝑡𝑒), cette relation est représentée
par l’équation de Poisson :
𝑠
∇2 𝑇 + =0 (2.06)
𝜆
et dans l’absence de sources internes de chaleur (𝑠 = 0) par l’équation de Laplace :
∇2 𝑇 = 0 (2.07)
On se limitera dans ce qui suit à l’étude du transfert de chaleur par conduction unidirectionnelle en régime
stationnaire sans sources internes de chaleur.
3.1. Corps homogènes
3.1.1. Mur simple à faces isothermes
On se placera dans le cas où le transfert de chaleur est unidirectionnel et où il n’y a pas de génération ni de
stockage d’énergie. Ce cas est représenté par un mur de bâtiment.
On considère un mur d’épaisseur 𝛿, constitué d’un matériau homogène et isotrope de conductivité thermique 𝜆
constante et de grandes dimensions transversales dont les faces extrêmes sont à des températures 𝑇1 et 𝑇2 (𝑇1 >
𝑇2 ).
En effectuant un bilan thermique sur le système (𝑆) constitué par le mur dans la direction des abscisses 𝑥 et dans
le cadre de ces hypothèse, l’équation différentielle de Laplace qui décrit la propagation de la chaleur dans ce mur
prend la forme :
𝑑2 𝑇
=0 (2.08)
𝑑𝑥 2
Figure 8 : Conduction thermique stationnaire dans un mur plan de conductivité thermique constante

Après intégrations successives, et en utilisant des conditions aux limites du premier type :
𝑥 = 0 → 𝑇 = 𝑇1
{ (2.09)
𝑥 = 𝛿 → 𝑇 = 𝑇2
On obtient l’équation du champ de température (ou équation de distribution de la température) dans le mur :
𝑥
𝑇(𝑥) = 𝑇1 − (𝑇1 − 𝑇2 ) (2.10)
𝛿
Unité : [𝐾]
Cette expression montre que la distribution de température dans un mur plan, homogène, isotrope et dans
l’hypothèse 𝜆 = 𝑐𝑠𝑡𝑒 est linéaire.
Connaissant le champ de température, la densité de flux thermique dans ce mur s’obtient en utilisant la loi de
Fourier :
𝜆
𝜑= (𝑇 − 𝑇2 )[𝑊/𝑚2] (2.11)
𝛿 1
⃗ est positif suivant l’axe 0𝑥 qui
Cette relation montre que le flux surfacique est constant et que pour 𝑇1 > 𝑇2 , 𝜑
représente aussi l’axe de propagation de la chaleur.
Le flux thermique total ou la puissance thermique transférée entre les deux faces d’aire 𝑆 du mur est donné par la
relation :
𝜆
𝛷 = 𝜑𝑆 = (𝑇 − 𝑇2 )𝑆 [𝑊] (2.12)
𝛿 1
La relation (2.11) permet de donner une signification physique de la conductivité thermique 𝜆 dans le cas du
mur :
𝜑. 𝛿
𝜆= [𝑊/(𝑚𝐾)] (2.13)
𝑇1 − 𝑇2
qui représente donc le flux thermique qui traverse la surface unitaire d’un mur plan d’épaisseur 𝛿 = 1𝑚, si la
différence de température entre ces faces est 𝑇1 − 𝑇2 = 1𝐾.
Si on écrit la relation (2.11) sous la forme :
𝑇1 − 𝑇2 ∆𝑇
𝜑= 𝛿
= (2.14)
𝑅𝑐𝑑
𝜆
On peut faire une analogie avec la loi d’Ohm en électrotechnique :

∆𝑈
𝐼= (2.15)
𝑅é𝑙
Où ∆𝑈 est la différence de potentiel et 𝑅é𝑙 la résistance électrique.
Dans la figure 2.4 est représentée l’analogie thermoélectrique dans le cas d’un mur simple. En faisant la
comparaison entre les relations qui donnent le flux surfacique et l’intensité du courant électrique, on constate
l’analogie suivante entre les deux types de grandeurs : 𝜑 ~ 𝐼 ; ∆𝑇 ~ ∆𝑈 et 𝑅𝑐𝑑 ~ 𝑅é𝑙 .
D’ici dérive le concept de résistance thermique par conduction qui, pour le mur plan simple a pour expression :
𝛿
𝑅𝑐𝑑 = [(𝑚2 𝐾)/𝑊] (2.16)
𝜆

Figure 9 : Analogie thermoélectrique dans le cas d’un mur plan simple

3.1.2. Mur simple en contact avec deux fluides
Considérons maintenant le cas de transfert thermique dans un mur simple, homogène et isotrope de surface 𝑆,
d’épaisseur 𝛿 et de conducteur 𝜆 = 𝑐𝑠𝑡𝑒, séparant deux milieux fluides en mouvement (fig. 10). La chaleur est
transférée entre le fluide chaud ayant la température 𝑇𝑓1 et le fluide froid ayant la température 𝑇𝑓2 , à travers le
mur. On suppose que ces températures et les coefficients de convection ℎ1 et ℎ2 afférents aux fluides sont
constants et uniformes par rapport aux surfaces de contact.
Figure 10 : Conduction thermique stationnaire dans un mur simple en contact avec deux fluides
Dans ces conditions, on doit déterminer la puissance thermique transférée entre ces deux fluides et les
températures sur les faces du mur. On sait que le flux thermique transféré entre le fluide en mouvement et une
surface solide en contact avec lui est donné par la loi de Newton :
𝜑 = ℎ(𝑇𝑓 − 𝑇𝑠 ) (2.17)
Où :
- 𝑇𝑓 : température du fluide [𝐾] ;
- 𝑇𝑠 : température de la surface solide [𝐾] ;

- ℎ : coefficient de transfert par convection ou coefficient de convection [𝑊/𝑚2 𝐾].
Si la relation (2.17) est écrite sous la forme :
𝑇𝑓 − 𝑇𝑠 ∆𝑇
𝜑= 1 = (2.18)
𝑅𝑠
ℎ
En faisant appel à l’analogie thermoélectrique, on a :

1
𝑅𝑠 = [(𝑚2 𝐾)/𝑊] (2.19)
ℎ
𝑅𝑠 est appellée résistance thermique à la surface ou résistance thermique superficielle.

En régime stationnaire et en absence de sources internes de chaleur dans le mur, le flux thermique unitaire se
conserve, c'est-à-dire que le flux cédé par le fluide chaud au mur, celui traversant le mur et celui reçu par le
fluide froid sont égaux :
𝜆
𝜑 = ℎ1 (𝑇𝑓1 − 𝑇1 ) = (𝑇1 − 𝑇2 ) = ℎ2 (𝑇2 − 𝑇𝑓2 ) (2.20)
𝛿
Ou encore :
𝑇𝑓1 − 𝑇1 𝑇1 − 𝑇2 𝑇2 − 𝑇𝑓2 𝑇𝑓1 − 𝑇𝑓2
𝜑= 1 = 𝛿
= 1 = 1 𝛿 1
(2.21)
ℎ1 ℎ2
+ +
𝜆 ℎ1 𝜆 ℎ2
Si on utilise les notions de résistance thermique à la surface (2.19) et résistance thermique conductive (2.16), on
peut écrire que la résistance thermique totale 𝑅𝑡𝑡 de l’ensemble fluides – mur est égale à :
1 𝛿 1
𝑅𝑡𝑡 = 𝑅𝑠1 + 𝑅𝑐𝑑 + 𝑅𝑠2 = + + [(𝑚2 𝐾)/𝑊] (2.22)
ℎ1 𝜆 ℎ2
La densité de flux thermique transférée entre les deux fluides est alors :
𝑇𝑓1 − 𝑇𝑓2 ∆𝑇
𝜑= 1 𝛿 1
= (2.23)
+𝜆+ℎ 𝑅𝑡𝑡
ℎ1 2
Le schéma électrique équivalent (figure 11) est représenté par un circuit ayant trois résistances thermiques en
série :
1
- 𝑅𝑠1 = ℎ résistance thermique superficielle au transfert de la chaleur entre le fluide chaud et la surface
1
du mur ;
𝛿
- 𝑅𝑐𝑑 = 𝜆 résistance thermique conductive au transfert de la chaleur par conduction à travers le mur ;
1
- 𝑅𝑠2 = ℎ résistance thermique superficielle au transfert de la chaleur entre le mur et le fluide froid.
2
Figure 11 : Schéma électrique équivalent au transfert thermique entre deux fluides séparés par un mur
plan simple
Par analogie avec la loi de Newton (2.17), l’expression de la densité de flux thermique peut être mise sous la
forme :
𝜑 = ℎ𝑔 (𝑇𝑓1 − 𝑇𝑓2 ) (2.24)
Où ℎ𝑔 représente le coefficient global de transfert thermique en [𝑊/𝑚2𝐾]. La comparaison des relations (2.23)
et (2.24) nous permet d’écrire :
1 1
ℎ𝑔 = = 1 𝛿 1
[𝑊/𝑚2 𝐾] (2.25)
𝑅𝑡𝑡 +𝜆+ℎ
ℎ1 2
Du point de vue physique, le coefficient global de transfert thermique représente l’inverse de la résistance
thermique totale au passage de la chaleur entre les deux fluides. Le flux total transféré entre ces deux fluides, si
on connaît l’aire de la surface 𝑆 du mur est :
𝛷 = 𝜑𝑆 = ℎ𝑔 (𝑇𝑓1 − 𝑇𝑓2 )𝑆 [𝑊] (2.26)

Les températures sur les faces extrêmes du mur 𝑇1 et 𝑇2 peuvent être déterminées en utilisant la relation (2.21)
qui exprime l’égalité des flux :
1
𝑇1 = 𝑇𝑓1 − 𝜑.
ℎ1
(2.27)
1 𝛿
𝑇2 = 𝑇𝑓1 − 𝜑. ( + )
{ ℎ1 𝜆
3.1.3. Cylindre creux simple à surfaces latérales isothermes
Ce problème est un modèle du cas pratique de transfert thermique par conduction à travers des éléments de
construction de forme cylindrique qui transportent les différents fluides, des tubes de ventilation, etc.
Considérons un cylindre creux fait d’un matériau homogène et isotrope de conductivité thermique 𝜆 constante,
limité par une surface intérieure de rayon 𝑟1 et d’une surface extérieure de rayon 𝑟2 qui ont des températures
constantes et uniformes 𝑇1 et 𝑇2 (𝑇1 > 𝑇2 ) et de longueur 𝑙 beaucoup plus grande que son diamètre, les
isothermes sont des surfaces cylindriques coaxiales de rayon 𝑟1 ≤ 𝑟 ≤ 𝑟2 . La conduction thermique est donc
unidirectionnelle et présente un gradient de température radial.
Figure 12 : Cylindre creux à surfaces latérales isothermes

Comme dans le cas du mur plan, on doit déterminer le champ de température et le flux thermique qui traverse le
matériau du cylindre, par rapport à un système de coordonnées cylindriques dont l’axe est celui du cylindre.
Les hypothèses faites et en coordonnées cylindriques conduisent à l’équation différentielle de la conduction sous
la forme :
𝑑2 𝑇 1 𝑑𝑇
+ =0 (2.28)
𝑑𝑟 2 𝑟 𝑑𝑟
Par son intégration, et en utilisation les conditions aux limites du premier type :
𝑟 = 𝑟1 , 𝑇 = 𝑇1
{ (2.29)
𝑟 = 𝑟2 , 𝑇 = 𝑇2
on obtient finalement l’équation du champ de température dans le cylindre creux :
𝑟
𝑙𝑛
𝑟1
𝑇(𝑟) = 𝑇1 − (𝑇1 − 𝑇2 ) 𝑟2 [𝐾] (2.30)
𝑙𝑛 𝑟
1
La densité de flux thermique est obtenue à partir de la loi de Fourier :

𝜆 𝑇1 − 𝑇2
𝜑= [𝑊/𝑚2] (2.31)
𝑟 𝑙𝑛 𝑟2
𝑟1
Il présente une variation radiale de type hyperbolique.

Le flux thermique total, ou puissance thermique, est donné par la relation :

𝑇1 − 𝑇2
𝛷 = 𝜑𝑆𝑐 = 1 𝑟 [𝑊] (2.32)
2𝜋𝜆𝑙
. 𝑙𝑛 𝑟2
1
Avec 𝑆𝑐 = 2𝜋𝑟𝑙 est la surface latérale du cylindre de rayon 𝑟 et de longueur 𝑙.

Il est constant. En accord avec cette relation, la résistance thermique conductive du cylindre creux est donnée par
l’expression :
1 𝑟2
𝑅𝑐𝑑 = 𝑙𝑛 [𝐾/𝑊] (2.33)
2𝜋𝜆𝑙 𝑟1
On utilise de manière fréquente dans les applications pratiques une puissance thermique linéaire, rapportée à un
mètre linéaire de cylindre :
𝛷 𝑇1 − 𝑇2
𝛷𝑙 = = 1 𝑟 [𝑊/𝑚] (2.34)
𝑙 𝑙𝑛 2
2𝜋𝜆 𝑟1
3.1.4. Cylindre creux simple en contact avec deux fluides

Pratiquement, ce cas est une modélisation d’une conduite qui transporte un fluide et qui est placé dans l’air.
Soit un cylindre creux qui sépare deux milieux fluides ayant des températures 𝑇𝑓1 et 𝑇𝑓2 (𝑇𝑓1 > 𝑇𝑓2) constantes et
uniformes. Il est constitué d’un matériau homogène et isotrope de conductivité 𝜆 constante. Le rayon intérieur 𝑟1
et le rayon extérieur 𝑟2 sont connus et de longueur 𝑙 beaucoup plus grande que ceux-ci. Entre le fluide chaud
ayant la température 𝑇𝑓1 qui s’écoule, par exemple, à l’intérieur du cylindre et le fluide froid de température 𝑇𝑓2
qui s’écoule à l’extérieur de celui-ci, a lieu un transfert thermique. En considérant que les coefficients de
transfert à la surface ℎ1 et ℎ2 afférents aux fluides sont constants et connus, on cherche à déterminer le flux
thermique transféré entre les deux fluides et les températures sur les surfaces limitant le cylindre.
Le problème est unidimensionnel et, en l’absence de sources internes de chaleur, la conservation du flux se
traduit par l’égalité du flux total cédé par le fluide chaud à la surface intérieure du cylindre par convection avec
celui traversant l’épaisseur du cylindre par conduction et avec celui reçu par le fluide froid par convection à la
surface extérieure du cylindre.
Figure 13 : Cylindre creux simple au contact avec deux fluides : distribution de température
Si on exprime ces flux thermiques par rapport à une surface cylindrique de longueur unitaire (𝑙 = 1), on obtient
les expressions :
𝛷𝑙 = ℎ1 . 2𝜋𝑟1 . (𝑇𝑓1 − 𝑇1 )
𝑇1 − 𝑇2
𝛷𝑙 = 2𝜋𝜆 𝑟 (2.35)
𝑙𝑛 ( 2 )
𝑟1
{𝛷𝑙 = ℎ2 . 2𝜋𝑟2 . (𝑇2 − 𝑇𝑓2 )

Où 𝑇1 et 𝑇2 représentent les températures sur les surfaces intérieure et extérieure du cylindre qui ne sont pas
connues a priori. En mettant les équations (2.35) sous la forme :
𝑇𝑓1 − 𝑇1 𝑇1 − 𝑇2 𝑇2 − 𝑇𝑓2
𝛷𝑙 = 1 = 1 𝑟2 = 1 (2.36)
𝑙𝑛
2𝜋ℎ1 𝑟1 2𝜋𝜆 𝑟1 𝜋ℎ2 𝑟2
On obtient par sommation le flux total rapporté à l’unité de longueur du cylindre :

𝑇𝑓1 − 𝑇𝑓2
𝛷𝑙 = 1 1 𝑟 1 [𝑊/𝑚] (2.37)
2𝜋ℎ1 𝑟1
+ 2𝜋𝜆 𝑙𝑛 𝑟2 + 2𝜋ℎ
1 2 𝑟2
La comparaison de cette expression avec la loi d’Ohm nous permet de définir la résistance thermique totale
linéairedu cylindre creux, en contact avec deux fluides :
1 1 𝑟2 1
𝑅𝑡𝑙 = + 𝑙𝑛 + [𝑚𝐾/𝑊] (2. 38)
2𝜋ℎ1 𝑟1 2𝜋𝜆 𝑟1 2𝜋ℎ2 𝑟2
Le schéma électrique équivalent est représenté dans la figure suivante et il est formé de trois résistances
thermiques montées en série :
Figure 14 : Schéma électrique équivalent d’un Cylindre creux simple au contact avec deux fluides
1
 𝑅𝑠1 = 𝜋ℎ 𝑑 résistance thermique surfacique au transfert de chaleur entre le fluide chaud et la surface
1 1
intérieure du cylindre (ou simplement résistance thermique de la surface intérieure) ;

1 𝑟
 𝑅𝑐𝑑 = 2𝜋𝜆 𝑙𝑛 (𝑟2 ) résistance thermique au transfert de chaleur par conduction à travers le cylindre ;
1
1
 𝑅𝑠2 = 𝜋ℎ résistance thermique surfacique au transfert de chaleur entre la surface extérieure du
2 𝑑2
cylindre et le fluide froid (ou résistance thermique de la surface extérieure).

La relation (2.32) peut s’écrire alors sous la forme :
𝑅𝑡𝑙 = 𝑅𝑠1 + 𝑅𝑐𝑑 + 𝑅𝑠2
D’une manière analogue à la relation (2.25), on peut définir pour ce cas aussi un coefficient global de transfert
thermique linéaire sous la forme :
1 1
ℎ𝑔𝑙 = = 1 1 𝑑 1 [𝑊/𝑚𝐾] (2.39)
𝑅𝑡𝑙 + 2𝜋𝜆 𝑙𝑛 𝑑2 + 𝜋ℎ
𝜋ℎ1 𝑑1 1 2 𝑑2
Qui, introduit dans la relation (2.37), conduit à :

𝛷𝑙 = ℎ𝑔𝑙 (𝑇𝑓1 − 𝑇𝑓2 )[𝑊/𝑚] (2.40)
Le coefficient ℎ𝑔𝑙 caractérise l’intensité du transfert thermique entre deux fluides de températures différentes,
séparés par un cylindre creux. Du point de vue numérique, il est égal au flux total linéaire transféré entre les
deux fluides ayant une différence de température de 1𝐾 entre elles.
Les températures sur les surfaces du cylindre résultent des équations (2.37) et (2.38) :
𝛷𝑙
𝑇1 = 𝑇𝑓1 −
𝜋ℎ1 𝑑1
(2.41)
1 1 𝑑2
𝑇2 = 𝑇𝑓1 − 𝛷𝑙 ( + 𝑙𝑛 )
{ 𝜋ℎ1 𝑑1 2𝜋𝜆 𝑑1

3.2. Corps ayant des structures multicouches

Ces corps ont une structure composée de plusieurs couches présentant un contact parfait ou pas.
Chaque couche est homogène et isotrope et ayant une conductivité thermique constante ainsi les dimensions
géométriques connues.
On se propose de déterminer le flux thermique transféré et les températures des surfaces de contact pour
quelques types de structures utilisées d’une manière usuelle.
3.2.1. Mur multicouches à surfaces isothermes
Considérons un mur composé par 𝑛 couches homogènes et isotropes d’épaisseurs 𝛿𝑗 (𝑗 = 1, 2, … , 𝑛) et ayant les
conductivités constantes 𝜆𝑗 en contact parfait. Les surfaces extrêmes ont les températures 𝑇1 et, respectivement,
𝑇𝑛+1 < 𝑇1 .
Figure 15 : Mur composite ayant deux surfaces isothermes
En régime stationnaire et en absence de sources internes de chaleur, le flux thermique surfacique qui traverse
chaque couche du mur reste constant. Suite à la relation (2.11), on obtient le système suivant :
𝜆1
𝜑1 = (𝑇 − 𝑇2 )
𝛿1 1
𝜆2
𝜑2 = (𝑇2 − 𝑇3 )
𝛿2
………… …………
𝜆𝑗 (2.42)
𝜑𝑗 = (𝑇𝑗 − 𝑇𝑗+1 )
𝛿𝑗
……………………
𝜆𝑛
𝜑𝑛 = (𝑇𝑛 − 𝑇𝑛+1 )
{ 𝛿𝑛
Sachant qu’il y a conservation du flux thermique surfacique, on peut écrire :
𝜑1 = 𝜑2 = ⋯ = 𝜑𝑛 = 𝑐𝑠𝑡𝑒 (2.43)
Cette relation permet de déduire l’expression du flux thermique surfacique transféré par le mur :
𝑇1 − 𝑇2 𝑇2 − 𝑇3 𝑇𝑛 − 𝑇𝑛+1
𝜑= 𝛿1
= 𝛿2
=⋯= 𝛿𝑛
(2.44)
𝜆1 𝜆2 𝜆𝑛
D’où :
𝑇1 − 𝑇𝑛+1 𝑇1 − 𝑇𝑛+1
𝜑 = 𝛿1 𝛿2 𝛿𝑛
= 𝛿𝑗
(2.45)
𝜆1
+
𝜆2
+ ⋯+
𝜆𝑛
∑𝑛𝑗=1
𝜆𝑗
𝛿𝑗
La valeur ∑𝑛𝑗=1 est appelée résistance thermique conductive totale du mur composé.
𝜆𝑗

Le flux total est :

𝑇1 − 𝑇𝑛+1
𝛷 = 𝜑𝑆 = 𝛿𝑗
𝑆 [𝑊] (2.46)
∑𝑛𝑗=1
𝜆𝑗
Où 𝑆 est la surface du mur.

Les températures à l’interface entre deux couches voisines sont obtenues à partir du système d’équations (2.42)
et le flux 𝜑 avec l’équation (2.44) :
𝛿1
𝑇2 = 𝑇1 − 𝜑
𝜆1
𝛿2 𝛿1 𝛿2
𝑇3 = 𝑇2 − 𝜑 = 𝑇1 − 𝜑 ( + )
𝜆2 𝜆 𝜆 (2.47)
………………… … … … … … …1 … …2…
𝑛
𝛿𝑗
𝑇𝑛+1 = 𝑇1 + 𝜑 ∑
𝜆𝑗
{ 𝑗=1
La comparaison entre le transfert thermique par conduction à travers un mur multicouches et un mur simple peut
se faire en introduisant la notion de conductivité thermique équivalente𝝀é𝒒 . Elle représente la conductivité d’un
mur simple et homogène, ayant l’épaisseur et la résistance thermique égales respectivement à l’épaisseur et à la
résistance thermique du mur multicouches, c’est – à – dire :
𝑛
∑𝑛𝑗=1 𝛿𝑗 𝛿𝑗
=∑
𝜆é𝑞 𝜆𝑗
𝑗=1
D’où, on obtient :
∑𝑛𝑗=1 𝛿𝑗
𝜆é𝑞 = 𝛿𝑗
(2.48)
∑𝑛𝑗=1
𝜆𝑗
Si entre les couches les contacts ne sont pas parfait, il faut tenir compte des résistances de contact 𝑅𝑐𝑖𝑗 entre les
couches 𝑖 et 𝑗 (𝑖, 𝑗 = 1, 2, … , 𝑛).

3.2.2. Mur multicouches en contact avec deux fluides
Dans ce cas, le mur multicouches a ses faces extrêmes en contact avec deux fluides de températures 𝑇𝑓1 et 𝑇𝑓2
(𝑇𝑓1 > 𝑇𝑓2 ) uniformes et constantes. Les coefficients surfaciques de transfert thermique sont ℎ1 et ℎ2 (figure
2.14). Sous les hypothèses de régime stationnaire de transfert thermique par conduction et en l’absence de
sources internes de chaleur dans le mur, le flux thermique surfacique transféré entre les deux fluides et qui
traverse le mur se conserve. Ainsi l’on peut écrire :
𝜆1 𝜆2
𝜑 = ℎ1 (𝑇𝑓1 − 𝑇1 ) = (𝑇1 − 𝑇2 ) = (𝑇2 − 𝑇3 ) = ⋯ = ℎ2 (𝑇𝑛+1 − 𝑇𝑓2 ) (2.49)
𝛿1 𝛿2
Ou encore :
𝑇𝑓1 − 𝑇1 𝑇1 − 𝑇𝑛+1 𝑇𝑛+1 − 𝑇𝑓2
𝜑= 1 = 𝛿
= 1 (2.50)
ℎ1
∑𝑛𝑗 𝑗 ℎ2
𝜆 𝑗
Cette relation nous permet d’avoir, par sommation, l’expression finale du flux thermique surfacique :
𝜑= 1 𝛿 1
(2.51)
ℎ1
+ ∑𝑛𝑗 𝜆𝑗 + ℎ
𝑗 2
Où l’expression :

𝑛
1 𝛿𝑗 1
𝑅𝑡𝑡 = + ∑ + [𝑚2 𝐾/𝑊] (2.52)
ℎ1 𝜆𝑗 ℎ2
𝑗
𝑅𝑡𝑡 représente la résistance thermique conductive totale du mur multicouches en contact avec deux fluides.
Dans ce cas, le coefficient de transfert thermique global ℎ𝑔 a la forme :
1 1
ℎ𝑔 = = 𝛿𝑗
[𝑊/(𝑚2 𝐾)] (2.53)
𝑅𝑡𝑡 1
+ ∑𝑛𝑗 𝜆 +ℎ
1
ℎ1 𝑗 2
Qui conduit à une expression plus synthétique du flux thermique surfacique :

𝜑 = ℎ𝑔 (𝑇𝑓1 − 𝑇𝑓2 ) [𝑊/𝑚2 ] (2.54)
Le flux total est égal à :
𝛷 = 𝜑𝑆 = ℎ𝑔 𝑆(𝑇𝑓1 − 𝑇𝑓2 ) [𝑊] (2.55)
Où 𝑆 est la surface totale du mur.
Figure 16 : Mur composite en contact avec deux fluides

Après le calcul du flux thermique surfacique par l’équation (2.51), les températures sur les surfaces externes du
mur et dans les plans de contact entre les couches sont déterminées en utilisant l’expression (2.47). On obtient
ainsi :
1
𝑇1 = 𝑇𝑓1 − 𝜑
ℎ1
𝛿1 1 𝛿1
𝑇2 = 𝑇1 − 𝜑 = 𝑇𝑓1 − 𝜑 ( + )
𝜆 ℎ1 𝜆1 [𝐾] (2.56)
… … … … … …1 … … … … … … … ……
𝑛
1 𝛿𝑗
𝑇𝑛+1 = 𝑇𝑓1 − 𝜑 ( + ∑ )
ℎ1 𝜆𝑗
{ 𝑗
3.2.3. Cylindre creux multicouches de surfaces latérales isothermes

D’une manière analogue au cas du mur, le cas d’un cylindre creux composé de 𝑛 couches (figure 2.17) de
diamètres et conductivités thermiques connus, sous les hypothèses précédentes et en utilisant la relation (2.34) du
flux thermique linéique, permet d’obtenir :
𝑇1 − 𝑇2 𝑇2 − 𝑇3 𝑇𝑛 − 𝑇𝑛+1
𝛷𝑙 = 1 𝑟2 = 1 𝑟3 =⋯= 1 𝑟𝑛+1 (2.57)
𝑙𝑛 𝑙𝑛 𝑙𝑛
2𝜋𝜆1 𝑟1 2𝜋𝜆2 𝑟2 2𝜋𝜆𝑛 𝑟𝑛

D’où :
𝑇1 − 𝑇𝑛+1
𝛷𝑙 = 1 1 𝑟𝑗+1 [𝑊/𝑚] (2.58)
∑𝑛𝑗=1 𝑙𝑛
2𝜋 𝜆𝑗 𝑟𝑗
L’expression
𝑛
1 1 𝑟𝑗+1
𝑅𝑡𝑡 = ∑ 𝑙𝑛 [𝑚𝐾/𝑊] (2.59)
2𝜋 𝜆𝑗 𝑟𝑗
𝑗=1
𝑅𝑡𝑡 représente la résistance conductive totale linéaire du cylindre creux multicouches.

La notion de conductivité thermique équivalente𝜆é𝑞 pour un cylindre multicouches garde la même signification
que celle définie pour un mur : la résistance thermique calculée à partir de la conductivité équivalente doit être
égale à la somme des résistances du cylindre multicouches :
𝑛
1 𝑑𝑛+1 1 1 𝑑𝑗+1
𝑙𝑛 = ∑ 𝑙𝑛
2𝜋𝜆é𝑞 𝑑1 2𝜋 𝜆𝑗 𝑑𝑗
𝑗=1
Et donc
𝑑𝑛+1
𝑙𝑛 𝑑1
𝜆é𝑞 = 1 𝑑𝑗+1
[𝑊/𝑚𝐾] (2.60 )
∑𝑛𝑗=1 𝑙𝑛
𝜆𝑗 𝑑𝑗
Les températures à l’interface de deux couches peuvent être calculées en utilisant les relations (2.57) après avoir
déterminé le flux total linéaire 𝛷𝑙 (équation (2.56)) :
1 𝑑2
𝑇2 = 𝑇1 − 𝛷𝑙 𝑙𝑛
2𝜋𝜆1 𝑑1
1 𝑑2 1 𝑑3
𝑇3 = 𝑇1 − 𝛷𝑙 ( 𝑙𝑛 + 𝑙𝑛 )
2𝜋𝜆1 𝑑1 2𝜋𝜆2 𝑑2 (2.61)
………………………………………
𝑛−1
1 𝑑𝑗+1
𝑇𝑛 = 𝑇1 − 𝛷𝑙 ∑ 𝑙𝑛
2𝜋𝜆𝑗 𝑑𝑗
{ 𝑗=1
Figure 17 : Cylindre creux composite de surfaces isothermes

3.2.4. Cylindre creux multicouches en contact avec deux fluides
En effet, on peut déterminer le flux thermique total, rapporté à l’unité de longueur du cylindre, transféré entre les
deux fluides à travers la structure multicouches de celui-ci en tenant compte dans les relations (2.42) des
résistances conductives de chaque couche. Donc, on peut écrire :

𝑇𝑓1 − 𝑇1 𝑇1 − 𝑇2 𝑇𝑛 − 𝑇𝑛+1 𝑇𝑛+1 − 𝑇𝑓2

𝛷𝑙 = 1 = 1 𝑑2 =⋯= 1 𝑑𝑛+1 = 1 (2.62)
𝜋𝑑1 ℎ1 2𝜋𝜆1
𝑙𝑛 𝑑 2𝜋𝜆𝑛
𝑙𝑛 𝑑𝑛 𝜋𝑑𝑛+1ℎ2
1
D’où, par sommation :

𝛷𝑙 = 1 1 1 𝑑𝑗+1 1
[𝑊/𝑚] (2.63)
∑ 𝑛
+ 𝑙𝑛 +
𝜋𝑑1 ℎ1 2𝜋 𝑗=1 𝜆 𝑗𝑑 𝑗 𝜋𝑑𝑛+1 ℎ2
Ce flux a aussi pour expression :

𝛷𝑙 = = ℎ𝑔𝑙 (𝑇𝑓1 − 𝑇𝑓2 ) (2.64)
𝑅𝑡𝑙
L’expression
𝑛
1 1 1 𝑑𝑗+1 1
𝑅𝑡𝑙 = + ∑ 𝑙𝑛 + [𝑚𝐾/𝑊] (2.65)
𝜋𝑑1 ℎ1 2𝜋 𝜆𝑗 𝑑𝑗 𝜋𝑑𝑛+1 ℎ2
𝑗=1
𝑅𝑡𝑙 représente la résistance thermique totale linéaire du cylindre creux multicouches et son inverse ℎ𝑔𝑙 est le
coefficient de transfert thermique total linéaire de celui-ci.
1 1
ℎ𝑔𝑙 = = 1 1 1 𝑑𝑗+1 1
[𝑊/𝑚𝐾] (2.66)
𝑅𝑔𝑙 + 2𝜋 ∑𝑛𝑗=1 𝜆 𝑙𝑛 + 𝜋𝑑
𝜋𝑑1 ℎ1 𝑗 𝑑𝑗 𝑛+1 ℎ2
Le système d’équation (2.62) permet la détermination des températures (𝑇1 et 𝑇𝑛+1 ) sur les surfaces extrêmes du
cylindre et au contact entre les différentes couches, après le calcul de 𝛷𝑙 avec la relation (2.63) :
1
𝑇1 = 𝑇𝑓1 − 𝛷𝑙
𝜋𝑑1 ℎ1
1 1 𝑑2
𝑇2 = 𝑇𝑓1 − 𝛷𝑙 ( + 𝑙𝑛 )
𝜋𝑑1 ℎ1 2𝜋𝜆1 𝑑1 [𝐾] (2.67)
………………………………………
𝑛−1
1 1 𝑑𝑗+1
𝑇𝑛+1 = 𝑇𝑓1 − 𝛷𝑙 ( +∑ 𝑙𝑛 )
𝜋𝑑1 ℎ1 2𝜋𝜆𝑗 𝑑𝑗
𝑗=1 }

Chapitre 3 Transfert de Chaleur par Convection

1. Analyse dimensionnelle
1.1. Dimensions fondamentales
On peut exprimer les grandeurs physiques en fonction d'un nombre limité de dimensions fondamentales.
Exemples : Vitesse : 𝐿. 𝑇 −1 ; viscosité dynamique : 𝑀. 𝐿−1 . 𝑇 −1 ; force : 𝑀. 𝐿. 𝑇 −2.
Sur ces exemples on voit que le nombre de dimensions fondamentales est de 3 : Masse 𝑀, Longueur 𝐿, Temps 𝑇.
Ces trois dimensions fondamentales ne sont pas toujours suffisantes. Pour les problèmes de transfert de chaleur,
il est nécessaire d'ajouter une 4ème dimension : la température 𝛩 et, lorsque l'échange d'énergie entre grandeurs
mécaniques et grandeurs thermiques ne sera pas mesurable, on ajoutera la quantité de chaleur 𝑄 qui sera
considérée comme une 5ème dimension.
Remarque : Q, homogène à un travail qui s'exprime en fonction des dimensions fondamentales 𝑀, 𝐿 et 𝑇 par
𝑄 = 𝑀. 𝐿. 𝑇 −2 n'est pas une vraie dimension fondamentale.
La méthode d'analyse dimensionnelle, qui repose sur le principe de l'homogénéité dimensionnelle des termes
d'une équation, est connue sous le nom de théorème de Vaschy-Buckingam ou théorème des groupements 𝜋.
1.2. Principe de la méthode
Si 1'on peut représenter mathématiquement une loi physique en exprimant une variable physique 𝐺1 en fonction
d'un certain nombre d'autres variable physiques indépendantes 𝐺2 , …,𝐺𝑛 , c’est à dire si 𝐺1 = 𝑓(𝐺2 , 𝐺3 , … , 𝐺𝑛 ) ou
encore 𝑓(𝐺1 , 𝐺2 , … , 𝐺𝑛 ) = 0, le problème peut être simplifié de la manière suivante :
 On écrit pour chaque variable 𝐺𝑖 , l'équation dimension en fonction des dimensions fondamentales. On
dispose alors de 𝑛 équations qui ont nécessité 𝑝 dimensions fondamentales pour caractériser toutes les
grandeurs physiques.
 On prélève 𝑝 de ces 𝑛 équations que l'on considère comme équations de base. Bien que le choix des
équations prélevées soit arbitraire, il faut toutefois que chaque dimension fondamentale apparaisse au
moins une fois sur l'ensemble des 𝑝 équations.
 Les (𝑛 − 𝑝) équations restantes se présentent alors sous forme de (𝑛 − 𝑝) rapports sans dimensions
appelés groupements 𝜋 qui sont des "grandeurs réduites". On obtient alors une équation réduite :
𝑔(𝜋1 , 𝜋2 , … , 𝜋𝑛−𝑝 )
Un groupement 𝜋 est le rapport d'une équation dimension d'une grandeur physique n'appartenant pas à
l'ensemble des équations de base au produit des équations de base, chacune d'elle étant portée à une certaine
puissance :
[𝐺𝑖 ]
𝜋𝑖 = 𝑒𝑖
[𝐺1 ]𝑎𝑖 [𝐺2 ]𝑏𝑖 … [𝐺𝑝 ]
Pour chaque dimension fondamentale 𝑀, 𝐿, 𝑇, 𝜃, 𝑄 figurant au dénominateur, on fait la somme des exposants
que l'on identifie avec l'exposant de la même dimension figurant dans l'équation dimension de la grandeur
physique du numérateur. On obtient ainsi un système linéaire de 𝑝 équations dont la résolution permet de
déterminer les 𝑝 exposants des équations de base du dénominateur.
Il suffit alors d'écrire le rapport 𝜋 en fonction des grandeurs physiques attachées aux équations dimensions de
départ.

1.3. Exemple d’application

Considérons un fluide en circulation forcée dans une canalisation cylindrique pour lequel on se propose de
déterminer le coefficient de convection ℎ relatif au transfert de chaleur fluide-paroi qui correspond à une
convection forcée :
Figure 18 : Schéma de la configuration étudiée

1.3.1. Détermination des grandeurs physiques
Il faut déterminer tous les paramètres dont dépend la densité de flux de chaleur 𝜑 (liée à ℎ par 𝜑 = ℎ. ∆𝑇), ce
sont ici :
 Les caractéristiques du fluide :
 𝜆 coefficient de conductibilité thermique,
 𝑐𝑝 chaleur massique,
 𝜌 masse volumique,
 𝜇 viscosité dynamique.
 Les caractéristiques de l'écoulement :
 𝑢 vitesse moyenne du fluide.
 La géométrie de la surface d'échange :
 𝐷 diamètre de la conduite.
 L'écart de température paroi-fluide ∆𝑇.
d'où on peut réécrire :
𝑓(𝜆, 𝑐𝑝 , 𝜌, 𝜇, 𝑢, 𝐷, ∆𝑇, 𝜑) = 0
1.3.2. Equation dimension de chaque grandeur
Il faut ensuite écrire l'équation aux dimensions fondamentales 𝑀, 𝐿, 𝑇, 𝛩, 𝑄 de chacune des grandeurs, ce qui
s'écrit :
Tableau 1 : Grandeurs physiques et leurs dimensions (unités)
Grandeurs physiques Symboles Unités S.I. Dimensions
2 [𝑀𝑇 3 𝛩 −1 ]
Coefficients de convection ℎ 𝑊/𝑚 𝐾
Longueur caractéristique 𝐷 𝑚 [𝐿]
Conductivité thermique du fluide 𝜆 𝑊/𝑚𝐾 [𝑄𝑇 −1 −1
𝐿 𝛩 −1 ]
Vitesse du fluide 𝑢 𝑚/𝑠 [𝐿𝑇 −1 ]
Viscosité dynamique 𝜇 [𝑀. 𝑇 −1 . 𝐿−1 ]
Chaleur spécifique du fluide 𝑐𝑝 𝐽/𝑘𝑔𝐾 [𝑄. 𝑀−1 . 𝛩 −1 ]
Masse volumique 𝜌 𝑘𝑔/𝑚3 [𝑀. 𝐿−3 ]
Temps 𝑡 𝑠 [𝛩]
Densité de flux de chaleur 𝜑 𝑊/𝑚2 [𝑄. 𝑇 −1 . 𝐿−2 ]

1.3.3. Détermination des groupements 𝝅

Il faut maintenant choisir 5 équations de base (Toutes les dimensions fondamentales ont été utilisées) de façon à
ce que les 5 dimensions fondamentales figurent au moins une fois dans l'ensemble des équations.
Prenons par exemple : 𝜆, 𝜌, 𝑢, 𝐷, ∆𝑇, il reste 𝜑, 𝑐𝑝 et 𝜇.
On écrit alors les 3 rapports sans dimension correspondants à ces variables sous la forme :
𝜑
𝜋1 = 𝑎1 𝑏1 𝑐1 𝑑1 𝑒1
∆𝑇 𝜆 𝜌 𝐷 𝑢
𝑐𝑝
𝜋2 = 𝑎2 𝑏2 𝑐2 𝑑2 𝑒2
∆𝑇 𝜆 𝜌 𝐷 𝑢
𝜇
𝜋3 = 𝑎 𝑏 𝑐 𝑑 𝑒
∆𝑇 𝜆 𝜌 3 𝐷 3 𝑢 3
3 3
Pour chaque rapport 𝜋, on remplace les grandeurs physiques par leurs équations dimensions ce qui donne par
exemple pour 𝜋1 :
[𝑄. 𝑇 −1 . 𝐿−2 ]
[𝜋1 ] =
[𝛩]𝑎1 [𝑄. 𝑇 −1 . 𝐿−1 . 𝜃 −1 ]𝑏1 [𝑀. 𝐿−3 ]𝑐1 [𝐿]𝑑1 [𝐿. 𝑇 −1 ]𝑒1
Pour chaque dimension fondamentale, on identifie les exposants de puissance entre numérateur et dénominateur
relatifs à une même dimension ce qui conduit au système :
 𝑄 ∶ 1 = 𝑏1
 𝑇 ∶ −1 = −𝑏1 − 𝑒1
 𝐿 ∶ −2 = −𝑏1 − 3𝑐1 + 𝑑1 + 𝑒1
 𝛩 ∶ 0 = 𝑎1 − 𝑏1
 𝑀 ∶ 0 = 𝑐1
Après résolution du système, on obtient les valeurs suivantes :
𝑎1 = 𝑏1 = 1 ; 𝑐1 = 𝑒1 = 0 ; 𝑑1 = −1
Le groupement 𝜋1 s'écrit donc :
𝜑 𝜑. 𝐷
𝜋1 = =
∆𝑇. 𝜆. 𝐷−1 ∆𝑇. 𝜆
Ce qui en remplaçant 𝜑 = ℎ. ∆𝑇 dans l’expression de 𝜋1 conduit à :
ℎ. 𝐷
𝜋1 =
𝜆
En procédant de la même manière pour les groupements 𝜋2 et 𝜋3, on obtient :
𝜌. 𝑢. 𝐷. 𝑐𝑝
𝜋2 =
𝜆
𝜇
𝜋3 =
𝜌. 𝑢. 𝐷
Le théorème de Vaschy-Buckingam nous permet d'affirmer que la relation 𝑓(𝜆, 𝑐𝑝 , 𝜌, 𝜇, 𝑢, 𝐷, ∆𝑇, 𝜑) = 0 entre
les 8 variables peut s'exprimer à l'aide des trois nombres sans dimension 𝜋1, 𝜋2 et 𝜋3 sous la forme :
𝑓(𝜋1 , 𝜋2 , 𝜋3 ) = 0
Ou
𝜋1 = 𝑓(𝜋2 , 𝜋3 )
1.3.4. Signification physique des groupements
 𝜋1 est le nombre de Nusselt, il peut aussi s'écrire :

𝐷
𝜆
𝑁𝑢 = 1
ℎ
C'est donc le rapport de la résistance thermique de conduction par la résistance thermique de convection. Il est
d'autant plus élevé que la convection est prédominante sur la conduction. Il caractérise le type de transfert de
chaleur. Il est une mesure de l’amplification locale du transfert de chaleur dans le fluide due à ses mouvements.
Ce nombre est toujours supérieur à l’unité, à tout le moins, de l’ordre de l’unité. Lorsqu’il est de l’ordre de
l’unité, le transfert convectif est faible et le transfert est du type conduction pure.
 𝜋3 est l’inverse du nombre de Reynolds qui caractérise le régime d’écoulement dans la canalisation :
1 𝜇 𝜈
= 𝜋3 = =
𝑅𝑒 𝜌. 𝑢. 𝐷 𝑢. 𝐷
D’où :
𝜌. 𝑢. 𝐷 𝑢. 𝐷
𝑅𝑒 = =
𝜇 𝜈
 𝜋2 est le nombre de Peclet. On peut aussi 1'écrire :
𝜌. 𝑢. 𝐷. 𝑐𝑝 𝜌. 𝑢. 𝐷 𝜇. 𝑐𝑝
𝑃𝑒 = = .
𝜆 𝜇 𝜆
On fait apparaître un nouveau nombre adimensionnel :
𝜇. 𝑐𝑝 𝜌𝑐𝑝 𝜈
𝑃𝑟 = =
𝜆 𝜆
𝑃𝑟 est appelé nombre de Prandtl. Ce nombre est calculable pour un fluide donné indépendamment des conditions
expérimentales (il ne dépend que de la température) et caractérise l'influence de la nature du fluide sur le
transfert de chaleur par convection.
On préfère donc chercher une relation sous la forme :
𝑁𝑢 = 𝑓(𝑅𝑒, 𝑃𝑟)
1.4. Avantages de l’utilisation des grandeurs réduites
Ils concernent essentiellement la représentation, la comparaison et la recherche des résultats expérimentaux :
 La représentation des résultats expérimentaux est simplifiée on pourra avoir une courbe reliant 2
variables ou un abaque reliant 3 variables réduites au lieu d'une relation liant (3 + p) paramètres.
 La comparaison des résultats expérimentaux est aussi très rapide et aisée, quel que soit le chercheur,
même si le système d'unité utilisé est différent puisque les grandeurs réduites sont sans dimension. La
recherche des résultats expérimentaux est facilitée et ordonnée : s’il suffit de tracer une courbe entre
deux variables réduites, c'est qu'il suffit d'effectuer une seule série d'expériences.
Remarque :
Il faut toutefois bien comprendre que la méthode de l’analyse dimensionnelle qui fournit les grandeurs réduites
ne donne pas la forme de la relation qui les lie, la recherche de cette relation fait l'objet du dépouillement des
résultats expérimentaux.

Tableau 2 : Quelques groupements sans dimensions

Nombre Nombre Nombre Nombre Nombre Nombre de Nombre de
de de de de Peclet de Grashof Rayleigh
Reynolds Prandlt Nusselt Margoulis
Groupements 𝑅𝑒 𝑃𝑟 𝑁𝑢 𝑃𝑒 𝑀𝑎 𝐺𝑟 𝑅𝑎
𝜌𝑢𝐷 𝜇𝑐𝑝 ℎ𝐷 𝜌𝑢𝐷𝑐𝑝 ℎ 𝛽𝑔∆𝑇𝜌 𝐷2 3
𝛽𝑐𝑝 𝑔∆𝑇𝜌 2 𝐷3
= = = = = = =
𝜇 𝜆 𝜆 𝜆 𝜌𝑢𝑐𝑝 𝜇2 𝜆µ
= 𝑃𝑟. 𝐺𝑟
2. Généralités et définitions
Les transferts de chaleur qui s’effectuent simultanément avec des transferts de masse sont dits transferts de
chaleur par convection. Ce mode d’échange de chaleur existe au sein des milieux fluides dans lesquels il est
généralement prépondérant.
Dans tout ce qui suit, on n’étudiera que la convection sans changement de phase.
2.1. Convection naturelle et forcée
Selon la nature du mécanisme qui provoque le mouvement du fluide on distingue :
 La convection libre ou naturelle : le fluide est mis en mouvement sous le seul effet des différences de
masse volumique résultant des différences de températures sur les frontières et d’un champ de forces
extérieures (la pesanteur).
 La convection forcée : le mouvement du fluide est induit par une cause indépendante des différences de
température (pompe, ventilateur,...).
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.2. Régime d’écoulement
Compte tenu du lien entre le transfert de masse et le transfert de chaleur, il est nécessaire de prendre en compte
le régime d’écoulement. Considérons à titre d’exemple l’écoulement d’un fluide dans une conduite :
 En régime laminaire, l’écoulement s’effectue par couches pratiquement indépendantes.
Figure 19 : Schématisation d’un écoulement laminaire

Entre deux filets fluides adjacents les échanges de chaleur s’effectuent donc :
 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.
 En régime turbulent, l’écoulement n’est pas unidirectionnel :

Figure 20 : Schématisation d’un écoulement turbulent

L’échange de chaleur dans la zone turbulente s’effectue par convection et conduction dans toutes les directions.
On vérifie que la conduction moléculaire est généralement négligeable par rapport à la convection et à la «
diffusion turbulente » (mélange du fluide dû à l’agitation turbulente) en dehors de la sous-couche laminaire.
3. Flux de chaleur transféré par convection
3.1. Analogie de Reynolds
De même qu’au niveau moléculaire on explique la viscosité des gaz par la transmission des quantités de
mouvement des molécules lors des chocs intermoléculaires, on explique la transmission de la chaleur par la
transmission d’énergie cinétique lors de ces mêmes chocs.
Cette liaison intime des phénomènes de viscosité et de transfert de chaleur conduit à l’analogie de Reynolds :
dans un écoulement fluide avec transfert de chaleur, le profil des vitesses et le profil des températures sont liés
par une relation de similitude schématisée sur la figure 3.4.
Figure 21 : Représentation de l’analogie de Reynolds dans le cas d’un écoulement turbulent dans un tube
3.2. Couches limites dynamique et thermique
Quel que soit le régime d’écoulement, il demeure une couche limite dynamique dans laquelle l’écoulement est
laminaire et dont l’épaisseur est d’autant plus réduite que le nombre de Reynolds est grand. L’épaisseur de cette
couche limite varie en fonction de nombreux paramètres : nature du fluide, température, rugosité de la paroi, ...
L’analogie de Reynolds montre que le gradient thermique est particulièrement important au voisinage de la
paroi, dans une couche limite thermique qui se développe de manière analogue à la couche limite dynamique.
Quel que soit le régime d’écoulement du fluide, on considère que la résistance thermique est entièrement située
dans cette couche limite thermique qui joue le rôle d’isolant.
Figure 22 : Représentation du modèle de Prandtl pour un écoulement turbulent dans une conduite

3.3. Expression du flux

Quelque soit le type de convection (libre ou forcée) et quel que soit le régime d’écoulement du fluide (laminaire
ou turbulent), le flux de chaleur Ф est donné par la relation dite loi de Newton :
Ф = ℎ. 𝑆. ∆𝑇
Le problème majeur à résoudre avant le calcul du flux de chaleur consiste à déterminer le coefficient de transfert
de chaleur par convection ℎ qui dépend d’un nombre important de paramètres : caractéristiques du fluide, de
l’écoulement, de la température, de la forme de la surface d’échange,...
On trouvera dans le tableau 3 l’ordre de grandeur du coefficient de transfert de chaleur par convection pour
différentes configurations.
Tableau 3 : Ordre de grandeur du coefficient de transfert de chaleur par convection
4. Calcul du flux de chaleur en convection forcée

L’application de l’analyse dimensionnelle montre que la relation liant le flux de chaleur transféré par convection
aux variables dont il dépend peut être recherchée sous la forme d’une relation entre trois nombres
adimensionnels :
𝑁𝑢 = 𝑓(𝑅𝑒, 𝑃𝑟)
Définis par :
 Nombre de Nusselt :
ℎ. 𝐷
𝑁𝑢 =
𝜆
 Nombre de Reynolds :
𝜌. 𝑢. 𝐷
𝑅𝑒 =
𝜇
 Nombre de Prandtl :
𝜇. 𝑐𝑝
𝑃𝑟 =
𝜆
Où D est la dimension caractéristique de la géométrie considérée qui sera par exemple le diamètre hydraulique :
4. 𝑠𝑒𝑐𝑡𝑖𝑜𝑛 𝑑𝑒 𝑝𝑎𝑠𝑠𝑎𝑔𝑒
𝐷ℎ =
𝑝é𝑟𝑖𝑚è𝑡𝑟𝑒 𝑚𝑜𝑢𝑖𝑙𝑙é

pour un écoulement dans un conduit (égal au diamètre intérieur pour un conduit cylindrique), le diamètre
extérieur pour un écoulement extérieur perpendiculaire à un tube, la longueur pour un écoulement à surface libre
sur une plaque…
Le calcul d’un flux de chaleur transmis par convection forcée s’effectue donc de la manière suivante :
 Calculer les nombres adimensionnels de Reynolds et de Prandtl.
 Suivant la valeur de 𝑅𝑒 et la configuration, choisir alors la corrélation.
 Calculer le nombre de Nusselt 𝑁𝑢 à partir de la corrélation choisie.
 Puis calculer le coefficient de transfert par convection ℎ
𝜆. 𝑁𝑢
ℎ=
𝐷
 Déduire le flux thermique transféré :
Ф = ℎ. 𝑆. (𝑇𝑝 − 𝑇∞ )
Pour la convection forcée, les principales corrélations sont données pour les propriétés du fluide (𝑐𝑝 , 𝜌, 𝜆, 𝜇)
calculées à une température moyenne dite température de film :
𝑇𝑝 + 𝑇∞
𝑇𝑚 =
2
5. Calcul du flux de chaleur en convection naturelle
5.1. Mécanisme de la convection naturelle
Considérons un fluide au repos en contact avec une paroi plane à température 𝑇0 . Si l’on porte la paroi à une
température 𝑇 = 𝑇0 + ∆𝑇, le fluide au contact de la paroi va s’échauffer par conduction et la masse du volume
unité va passer de 𝜌0 à 𝜌0 − ∆𝜌.
Figure 23 : Représentation du mécanisme de convection naturelle

Il sera donc soumis à une force ascensionnelle 𝑓 = −∆𝜌. 𝑔. Le principe fondamental de la dynamique permet
d’évaluer l’accélération du fluide :
Pour un volume unité : 𝑚 = 𝜌 : ∆𝜌. 𝑔 = 𝜌. 𝛾 et :
∆𝜌
𝛾= .𝑔
𝜌
En introduisant le coefficient de dilatation cubique 𝛽 du fluide défini par :
1 ∆𝜌
𝛽=− ( )
𝜌 ∆𝑇 𝑃
Il vient que :
𝛾 = −𝛽. 𝑔. ∆𝑇
𝛽. 𝑔. ∆𝑇 est donc le module de l’accélération produite par l’expansion thermique due à la variation ∆𝑇 de la
température 𝑇0 . Ce mouvement du fluide induit par les différences de masse volumique résultantes des gradients
de température va donner naissance aux courants de convection.

Dans le cas d’un transfert de chaleur par convection naturelle le long d’une plaque plane, le coefficient de
convection dépend des caractéristiques du fluide : 𝜆, 𝜌, 𝜇, 𝑐𝑝 , 𝛽, 𝑔, de la paroi caractérisée par la longueur 𝐿 et
de l’écart de température ∆𝑇 aux bornes du film, ce que l’on peut traduire par une relation du type :
𝑓(𝜆, 𝜌, 𝜇, 𝑐𝑝 , 𝛽, 𝑔, 𝐿, ∆𝑇)
Dans le système 𝑀, 𝐿, 𝑇, 𝜃, 𝑄, cette relation entre 8 grandeurs se réduit à une relation entre trois nombres
adimensionnels :
𝑁𝑢 = (𝐺𝑟, 𝑃𝑟)
Définis par :
ℎ. 𝐷
𝑁𝑢 =
𝜆
 Nombre de Grashof :
𝛽. 𝑔. ∆𝑇. 𝜌 2 . 𝐿3
𝐺𝑟 =
𝜇2
𝜇. 𝑐𝑝
𝑃𝑟 =
𝜆
5.2. Signification physique du nombre de Grashof
Lorsque la masse unité du fluide, soumise à l’accélération 𝛽. 𝑔. ∆𝑇 subit une variation d’altitude 𝐿, la
conservation de l’énergie permet d’écrire :
𝑢2
= 𝛽. 𝑔. ∆𝑇. 𝐿
2
𝑢2
représente la variation d’énergie cinétique et 𝛽. 𝑔. ∆𝑇. 𝐿 la variation d’énergie potentielle.
2
On voit donc que le nombre de Grashof peut se mettre sous la forme :

1 𝑢𝐿𝜌 2
𝐺𝑟 = ( )
2 𝜇
Il est donc proportionnel au carré d’un nombre de Reynolds caractérisant l’écoulement. En pratique, en
convection naturelle, le courant qui prend naissance reste laminaire jusqu’à ce que le nombre de Grashof atteigne
une valeur d’environ 109 .
5.3. Calcul du flux de chaleur en convection naturelle
L’application de l’analyse dimensionnelle montre que la relation liant le flux de chaleur transféré par convection
aux variables dont il dépend peut être recherchée sous la forme d’une relation entre trois nombres
adimensionnels :
𝑁𝑢 = (𝐺𝑟, 𝑃𝑟)
Définis par :
ℎ. 𝐷
𝑁𝑢 =
𝜆
 Nombre de Grashof :
𝛽. 𝑔. ∆𝑇. 𝜌 2 . 𝐿3
𝐺𝑟 =
𝜇2

𝜇. 𝑐𝑝
𝑃𝑟 =
𝜆
Le calcul d’un flux de chaleur transmis par convection naturelle s’effectue donc de la manière suivante :
 Calculer les nombres adimensionnels de Grashof et de Prandtl.
 Suivant la valeur de 𝐺𝑟 et la configuration, choisir alors la corrélation.
 Calculer le nombre de Nusselt 𝑁𝑢 à partir de la corrélation choisie.
 Puis calculer le coefficient de transfert par convection ℎ
𝜆. 𝑁𝑢
ℎ=
𝐷
 Déduire le flux thermique transféré :
Ф = ℎ. 𝑆. (𝑇𝑝 − 𝑇∞ )
Pour la convection naturelle, les principales corrélations sont données pour les propriétés du fluide (𝑐𝑝 , 𝜌, 𝜆, 𝜇)
calculées à la température moyenne de film comme en convection forcée.
6. Corrélations utilisées en convection forcée
On se propose de donner quelques relations semi-empiriques du calcul du nombre de Nusselt moyen pour des
situations simples qui constituent des cas de référence. Il ne faut pas oublier que ces corrélations n’ont souvent
qu’un domaine limité de validité pas toujours explicite et qu’il faut les utiliser avec beaucoup de précautions.
6.1. Plaques planes
Le nombre de Reynolds critique : 𝑅𝑒𝑐 = 5. 105
 Ecoulement turbulent 𝑅𝑒 > 5. 105
Le nombre de Nusselt est :
𝑁𝑢 = 0,036. 𝑅𝑒 0,8. 𝑃𝑟1/3
 Ecoulement laminaire 𝑅𝑒 < 5. 105
Le nombre de Nusselt est :
𝑁𝑢 = 0,66. 𝑅𝑒 0,5 . 𝑃𝑟1/3
6.2. Transfert thermique en conduite en régime laminaire
Considérons une conduite parcourue un fluide et échangeant de la chaleur avec l’extérieur à travers sa parois
latérale. Le régime d’écoulement dans la conduite est laminaire.
L’expression du coefficient de transfert thermique convectif ℎ en régime laminaire est :
48 𝜆
𝑁𝑢 =.
11 𝐷
En régime de conduction pure, le nombre de Nusselt correspondant est indépendant de l’écoulement :
48
𝑁𝑢 =
11
6.3. Transfert thermique en conduite en régime turbulent
La valeur du nombre de Nusselt est calculée pour différentes valeurs du nombre de Reynolds 𝑅𝑒 et de Prandlt
𝑃𝑟. Nous avons quelques corrélations semi-empiriques, par exemple :
 la corrélation empirique de Dittus et Boelter :
𝑁𝑢 = 0,023. 𝑅𝑒 0,8 . 𝑃𝑟 𝑛
Avec 𝑛 = 0,3 en refroidissement et 𝑛 = 0,4 en chauffage et 0,7 < 𝑃𝑟 < 100.

 la corrélation de Lyon-Martinelli valable pour les métaux liquides à faible nombre de Prandlt et pour
un nombre de Péclet > 50 :
𝑁𝑢 = 7 + 0,025. 𝑅𝑒 0,8 . 𝑃𝑟 0,8
 la corrélation semi-empirique proposée par Burmeieter :
1
𝑁𝑢 𝑅𝑒 −4
𝑆𝑡 = = 0,0384. 1 1
𝑅𝑒. 𝑃𝑟 1 + 1,5𝑃𝑟 −6 (𝑃𝑟 − 1)𝑅𝑒 −8
6.4. Convection forcée autour d’un cylindre à température uniforme et constante placé
transversalement par rapport à l’écoulement
 Ecoulements laminaires et turbulents
𝑈𝐷
𝑁𝑢 = 0,43 + 0,53. 𝑃𝑟 0,31 . 𝑅𝑒 0,5, pour 1 < 𝑅𝑒 = 𝜈
< 4000
𝑁𝑢 = 0,43 + 0,193. 𝑃𝑟 0,31 . 𝑅𝑒 0,62, pour 4. 103 < 𝑅𝑒 < 4. 104

𝑁𝑢 = 0,43 + 0,0265. 𝑃𝑟 0,31 . 𝑅𝑒 0,81, pour 4. 104 < 𝑅𝑒 < 4. 105
𝐷 désignant le diamètre du cylindre
 Ecoulements turbulents
𝑃𝑟 0,25
𝑁𝑢 = 0,025. 𝑅𝑒 0,6 . 𝑃𝑟 0,38 . ( )
𝑃𝑟𝑟
𝑃𝑟𝑟 désignant le nombre de Prandlt de référence pour le calage de la corrélation.
6.5. Convection forcée autour d’une sphère solide à température uniforme et constante
 Conduction pure : 𝑁𝑢 = 2
 Ecoulements à faible nombre de Reynolds :
𝑁𝑢 = 1 + (1 + 𝑃𝑒)1/3 pour 0 < 𝑃𝑒 < 104 avec 𝑃𝑒 = 𝑅𝑒. 𝑃𝑟 = 𝑈𝐷/𝑎
 Ecoulements à grand nombre de Reynolds :
𝑁𝑢 = 1 + 𝑄. 𝑅𝑒 0,41, pour 0 < 𝑃𝑟 < 100 et 0 < 𝑅𝑒 < 400
𝑁𝑢 = 1 + 0,752. 𝑄. 𝑅𝑒 0,472, pour 102 < 𝑅𝑒 < 2. 103
𝑁𝑢 = 1 + 𝑄. (0,44. 𝑅𝑒 0,5 + 0,034. 𝑅𝑒 0,71 ), pour 2. 103 < 𝑅𝑒 < 105
Avec :
1 1/3
𝑄 = (1 + ) . 𝑃𝑟1/3
𝑃𝑒
Il existe aussi une corrélation dite de Ranz et Marshall
𝑁𝑢 = 2 + 0,6. 𝑅𝑒 0,5 . 𝑃𝑟1/3
6.6. Convection forcée à travers un faisceau de tubes cylindriques de section circulaire alignés par
rapport à l’écoulement
Le nombre de Nusselt global pour chaque tube est donnée par :
𝑛
𝑃𝑟 0,25
𝑁𝑢 = 𝐶. 𝑅𝑒𝑚𝑎𝑥 . 𝑃𝑟 0,36 . ( )
𝑃𝑟𝑟
𝑃𝑟𝑟 désignant le nombre de Prandlt de référence pris à la température de
sortie.
Le nombre de Reynolds est évalué à partir de la vitesse maximale 𝑢𝑚𝑎𝑥 à
travers les interstices du faisceau, soit :

𝑢𝑚𝑎𝑥 𝐷 𝐷+𝑙
𝑅𝑒𝑚𝑎𝑥 = 𝜈
, avec 𝑢𝑚𝑎𝑥 = 𝑈. 𝑙
Les valeurs de 𝐶 et 𝑛 sont :

 𝐶 = 0,27, 𝑛 = 0,63 pour 100 < 𝑅𝑒𝑚𝑎𝑥 < 2. 105
 𝐶 = 0,021, 𝑛 = 0,84 pour 2. 102 < 𝑅𝑒𝑚𝑎𝑥 < 2. 106
Le nombre de Nusselt décroît si le nombre 𝑁 de rangés de tube est inférieur à 20. Le coefficient de correction
varie de 0,8 (𝑁 = 2), 0,9 (𝑁 = 4) et 0,97 (𝑁 = 10).
6.7. Convection forcée à travers un faisceau de tubes cylindriques de section circulaire décalés
par rapport à l’écoulement
Le nombre de Nusselt global pour chaque tube est donnée :
𝑛
𝑃𝑟 0,25
𝑁𝑢 = 𝐶. 𝑅𝑒𝑚𝑎𝑥 . 𝑃𝑟 0,36 . ( )
𝑃𝑟𝑟
𝑃𝑟𝑟 désignant le nombre de Prandlt de référence pris à la température de
sortie.
Le nombre de Reynolds est évalué à partir de la vitesse maximale 𝑢𝑚𝑎𝑥 à
travers les interstices du faisceau, soit :
𝑢𝑚𝑎𝑥 𝐷 𝐷+𝑙 𝐷+𝑙
𝑅𝑒𝑚𝑎𝑥 = , avec 𝑢𝑚𝑎𝑥 = 𝑈. si 2. 𝑙𝑑 > 𝑙 et 𝑢𝑚𝑎𝑥 = 𝑈. 2.𝑙 si 2. 𝑙𝑑 < 𝑙
𝜈 𝑙 𝑑
Les valeurs de 𝐶 et 𝑛 sont :

 𝐶 = 0,35𝑚0,2 𝑛 = 0,6 pour 100 < 𝑅𝑒𝑚𝑎𝑥 < 2. 105 si 𝑚 < 0,5
 𝐶 = 0,4𝑛 = 0,6 pour 100 < 𝑅𝑒𝑚𝑎𝑥 < 2. 105 si 𝑚 > 0,5
 𝐶 = 0,022, 𝑛 = 0,84 pour 2. 105 < 𝑅𝑒𝑚𝑎𝑥 < 2. 106
Le nombre de Nusselt décroît si le nombre 𝑁 de rangés de tube est inférieur à 20. Le coefficient de correction
varie de 0,76 (𝑁 = 2), 0,89 (𝑁 = 4) et 0,97 (𝑁 = 10).
7. Corrélations utilisées en convection naturelle dans les configurations les plus utilisées
7.1. Convection naturelle au-dessus d’une plaque chaude ou en-dessous d’une plaque froide
horizontale
Dans cette configuration, les nombres sans dimension sont construits à partir d’une longueur caractéristique 𝐿
égale au rapport aire de la plaque/périmètre. Si la plaque est isotherme, alors le nombre de Nusselt moyen est
estimé par la corrélation suivante :
𝑁𝑢 = 𝐴. 𝑅𝑎𝑚
Avec :
 A = 0,54, 𝑚 = 0,25 si 𝑅𝑎 est tel que : 2. 104 < 𝑅𝑎 < 8. 106 (régime laminaire) ;
 A = 0,15, 𝑚 = 0,33 si 𝑅𝑎 est tel que : 8. 106 < 𝑅𝑎 < 1011 (régime turbulent).
7.2. Convection naturelle au-dessus d’une plaque froide ou en-dessous d’une plaque chaude
horizontale
Si la plaque est isotherme, alors le nombre de Nusselt moyen est estimé par la corrélation suivante :
Avec :
 A = 0,27, 𝑚 = 0,25 pour 𝑅𝑎 est tel que : 105 < 𝑅𝑎 < 1010

7.3. Convection naturelle au-dessus d’un cylindre horizontal chaud d’axe horizontal
Dans cette configuration, les nombres de Rayleigh et de Nusselt moyen sont construits à partir du diamètre du
cylindre. Si le cylindre est isotherme, alors le nombre de Nusselt moyen est estimé par la corrélation suivante :
Avec :
 A = 0,40, 𝑚 = 0 pour 0 < 𝑅𝑎 < 10−2
 A = 1,02, 𝑚 = 0,148 pour 10−2 < 𝑅𝑎 < 102
 A = 0,85, 𝑚 = 0,188 pour 102 < 𝑅𝑎 < 104
 A = 0,53, 𝑚 = 0,25 pour 104 < 𝑅𝑎 < 109
 A = 0,13, 𝑚 = 0,33 pour 109 < 𝑅𝑎 < 1011
Ces corrélations peuvent être condensées en une seule expression analytique valable pour une gamme assez large
de nombre de Rayleigh :
1
0,387𝑅𝑎6
𝑁𝑢1/2 = 0,6 + 8
9 27
0,559 16
(1 + ( ) )
𝑃𝑟
Pour 104 < 𝑅𝑎 < 1013

7.4. Convection naturelle au-dessus d’une sphère chaude
Dans cette configuration, les nombres de Rayleigh et de Nusselt moyen sont construits à partir du diamètre de la
sphère. Si la sphère est isotherme, le nombre de Nusselt moyen est estimé par la corrélation suivante :
𝑁𝑢 = 2 + 𝐴. 𝑅𝑎𝑚
Pour l’air 𝐴 = 0,43, 𝑚 = 0,25 pour 1 < 𝑅𝑎 < 105.
Plus généralement :
0,589
𝐴= 4
9 9
0,469 16
(1 + ( 𝑃𝑟 ) )
Pour 𝑃𝑟 ≥ 0,7 et 𝑅𝑎 > 1011 .

Chapitre 4 Transfert de Chaleur par Rayonnement

1. Introduction
La matière émet et absorbe de la chaleur en permanence. L’émission de chaleur est due à l’agitation moléculaire.
L’absorption de la chaleur par le milieu élève sa température et accroît son agitation moléculaire et donc
l’émission de chaleur. Le transfert de chaleur par rayonnement suppose l’existence d’un émetteur de
rayonnement, d’un espace dans lequel le rayonnement se propage et d’un récepteur sur ou dans lequel le
rayonnement est transformé. Ainsi, contrairement aux transferts thermiques par conduction ou par convection,
qui ont besoin d’un milieu matériel comme support de transfert thermique entre des points à températures
différentes, le transfert thermique par rayonnement peut avoir lieu dans le vide. Ce transfert thermique est très
fortement dépendant du niveau de la température des corps en présence. Il est très intense dès que les
températures dépassent la centaine de degrés Celsius.
2. Généralités. Définitions
2.1. Nature du rayonnement
Tous les corps, quelque soit son état : solide, liquide ou gazeux, ayant une température supérieure au zéro absolu,
émettent un rayonnement de nature électromagnétique. Cette émission d’énergie s’effectue au détriment de
l’énergie interne du corps émetteur.
Dans le vide, les ondes sont caractérisées par leur longueur d’onde𝜆0 et leur fréquence𝑓𝜆 . La relation qui lie la
longueur d’onde et la fréquence de l’onde est :
𝑐0 = 𝜆0 𝑓𝜆 (4.01)
Où 𝑐0 = 2,998. 108 𝑚/𝑠 est la vitesse de la lumière dans le vide. Dans un milieu matériel d’indice de réfraction
𝑛, la célérité 𝑐 de l’onde est :
𝑐0
𝑐=
𝑛
Et la longueur d’onde correspondante est :
𝜆0
𝜆𝑛 =
𝑛
Suivant leur longueur d’onde, les ondes électromagnétiques sont classées selon la nomenclature définie comme
suit :
 Les ondes de longueur d’onde inférieure à 10−6 𝜇𝑚 appartiennent à la catégorie des rayons cosmiques,
 entre 10−5 à 10−4 𝜇𝑚, ce sont des rayons gamma,
 de 10−4 à 10−2 𝜇𝑚 des rayons X, de 10−2 à 0,38 𝜇𝑚 des ultraviolets.
 La lumière visible par l’œil humain est comprise entre 0,38 à 0,76 𝜇𝑚,
 tandis que les infrarouges se situent entre 0,76 à 1000 𝜇𝑚.
C’est Sir William Herschel qui, grâce à l’utilisation d’un prisme dispersif, découvrit l’existence d’un
rayonnement porteur de calories mais non visible à l’œil, et l’a appelé rayonnement infrarouge. Au-delà de
1000 𝜇𝑚, on rencontre les ondes courtes puis les ondes radio. Dans le rayonnement visible, les différentes
couleurs sont associées à des longueurs d’ondes distinctes.

Figure 24 : Principe de l’expérience de William Herschel

Si l’on déplace le long de l’écran un thermomètre, on mesure la température 𝑇𝑒 caractérisant l’énergie reçue par
l’écran dans chaque longueur d’onde. En construisant la courbe 𝑇𝑒 = 𝑓(𝜆), on obtient la répartition spectrale de
l’énergie rayonnée pour la température 𝑇0 de la source. On constate alors que:
 L’énergie émise est maximale pour une certaine longueur d’onde 𝜆𝑚 variable avec 𝑇0 .
 L’énergie n’est émise que sur un intervalle [𝜆1, 𝜆2] de longueur d’onde caractérisant le rayonnement
thermique.
Figure 25 : Spectre des ondes électromagnétiques

2.2. Définitions
2.2.1. Classification
Les grandeurs physiques seront distinguées selon :
 La composition spectrale du rayonnement :
 Si la grandeur est relative à l’ensemble du spectre elle est dite totale.
 Si elle concerne un intervalle spectral étroit 𝑑𝜆 autour d’une longueur d’onde 𝜆 elle est dite
monochromatique : 𝐺𝜆 .
 La distribution spatiale du rayonnement :
 Si la grandeur est relative à l’ensemble des directions de l’espace, elle est dite hémisphérique.
 Si elle caractérise une direction donnée de propagation, elle est dite directionnelle : 𝐺𝑥
2.2.2. Rappel sur les angles solides élémentaires
L’angle solide sous lequel depuis un point 0 on voit une surface 𝑆 est par définition l’aire de la surface
intersection de la sphère de rayon unité et du cône de sommet 0 s’appuyant sur le contour de la surface 𝑆.
L’angle solide élémentaire 𝑑𝛺 sous lequel est vu d’un point 0 le contour d’une petite surface 𝑑𝑆 (assimilée à une
surface plane) peut être calculé par :
𝑑𝑆. cos 𝛼
𝑑𝛺 =
𝑟2

Figure 26 : Schéma de l’angle solide

Propriétés :
 La valeur d’un angle solide 𝛺est comprise entre 0 et 4π.
 Pour un cône de demi-angle au sommet 𝛼 : 𝛺 = 2𝜋(1 − cos 𝛼).
2.2.3. Définitions relatives aux sources
2.2.3.1. Flux
 On appelle flux d’une source S la puissance rayonnée notée 𝛷 par S dans tout l’espace qui l’entoure, sur
toutes les longueurs d’onde. Le flux 𝛷 s’exprime en 𝑊.
 Le flux envoyé par un élément de surface 𝑑𝑆 dans un angle solide élémentaire 𝑑𝛺 est noté 𝑑2 𝛷.
 Le flux envoyé dans tout l’espace par une surface élémentaire 𝑑𝑆 est noté 𝑑𝛷.
 Le flux envoyé par une surface 𝑆 dans l’angle solide 𝑑𝛺 entourant la direction 0𝑥 est noté 𝑑𝛷𝑥.
Nous avons donc les relations suivantes :
𝑑𝛷 = ∫ 𝑑2 𝛷
𝛺
𝛷 = ∫ 𝑑𝛷 = ∫ 𝑑𝛷𝑥
𝑆 𝛺
2.2.3.2. Emittance énergétique monochromatique

Un élément de surface 𝑑𝑆 émet un certain flux d’énergie par rayonnement dans toutes les directions du demi-
espace. Ce flux est réparti sur un intervalle de longueurs d’ondes. Si l’on considère le flux d’énergie 𝑑𝛷𝜆𝜆+𝑑𝜆
émis entre les deux longueurs d’ondes 𝜆 et 𝜆 + 𝑑𝜆, on définit l’émittance monochromatique d’une source à la
température 𝑇 par :
𝑑𝛷𝜆𝜆+𝑑𝜆
𝑀𝜆𝑇 = [𝑊/𝑚3] (4.02)
𝑑𝑆𝑑𝜆
2.2.3.3. Emittance énergétique totale
C’est la densité de flux de chaleur émise par rayonnement par 𝑑𝑆 sur tout le spectre des longueurs d’ondes. Elle
n’est plus fonction que de la température 𝑇 et de la nature de la source :
𝜆=∞
𝛷
𝑀𝜆 = ∫ 𝑀𝜆𝑇 . 𝑑𝜆 = [𝑊/𝑚2] (4.03)
𝑆
𝜆=0
2.2.3.4. Intensité énergétique dans une direction

On appelle intensité énergétique 𝐼𝑥 le flux par unité d’angle solide émis par une surface 𝑑𝑆 dans un angle solide
𝑑𝛺 entourant la direction 0𝑥 :
𝑑2 𝛷𝑥
𝐼𝑥 = (4.04)
𝑑𝛺

2.2.3.5. Luminance énergétique dans une direction

L’intensité énergétique élémentaire 𝑑𝐼𝑥 dans la direction 0𝑥 par unité de surface émettrice apparente 𝑑𝑆𝑥
s’appelle la luminance énergétique 𝐿𝑥 . En partant de la relation (4.04) :
𝐼𝑥 𝐼𝑥 𝑑2 𝛷𝑥
𝐿𝑥 = = = (4.05)
𝑑𝑆𝑥 𝑑𝑆. cos 𝛼 𝑑𝛺. 𝑑𝑆. cos 𝛼
Figure 27 : Schéma de définition des angles

2.2.4. Définitions relatives à un récepteur
2.2.4.1. Eclairement
C’est l’homologue de l’émittance pour une source. L’éclairement est le flux reçu par unité de surface réceptrice,
en provenance de l’ensemble des directions.
2.2.4.2. Réception du rayonnement par un solide
Quand un rayon incident d’énergie 𝛷𝜆 frappe un corps à la température 𝑇, une partie 𝛷𝜆 . 𝜌𝜆𝑇 de l’énergie
incidente est réfléchie par la surface S, une autre partie 𝛷𝜆 . 𝛼𝜆𝑇 est absorbée par le corps qui s’échauffe et le reste
𝛷𝜆 . 𝜏𝜆𝑇 est transmis et continue son chemin d’où la relation :
𝛷𝜆 = 𝛷𝜆 . 𝜌𝜆𝑇 + 𝛷𝜆 . 𝛼𝜆𝑇 + 𝛷𝜆 . 𝜏𝜆𝑇
Finalement :
𝜌𝜆𝑇 + 𝛼𝜆𝑇 + 𝜏𝜆𝑇 = 1 (4.06)
Figure 28 : Schématisation de la répartition d’un flux incident de rayonnement sur un solide

Où :
 𝜌𝜆𝑇 est le pouvoir monochromatique réfléchissant,
 𝛼𝜆𝑇 pouvoir monochromatique absorbant,
 𝜏𝜆𝑇 pouvoir monochromatique filtrant.
Ces pouvoirs monochromatiques sont fonction de la nature du corps, de son épaisseur, de sa température 𝑇, de la
longueur d’onde 𝜆 du rayonnement incident et de l’angle d’incidence.
Si l’on considère l’énergie incidente sur tout le spectre des longueurs d’onde, on obtient les pouvoirs
réfléchissants 𝜌𝑇 , absorbant 𝛼𝑇 et filtrant 𝜏𝑇 totaux.
2.2.5. Corps noir, corps gris
2.2.5.1. Corps noir
Ce type de corps est un modèle. Les valeurs de la luminance et l’émittance qui lui sont attachées servent de
référence au rayonnement des corps réels.

Le corps noir est le corps qui, à une température donnée, émet le plus de puissance. Le corps noir peut être aussi
le corps qui absorbe le plus de rayonnement de toutes les radiations qu’il reçoit indépendamment de son
épaisseur, de sa température, de l’angle d’incidence et de la longueur d’onde du rayonnement incident, il est
défini par : 𝛼𝜆𝑇 = 1. Cette énergie possède une distribution spectrale donnée par la loi de Planck.
Une surface enduite de noir de fumée est approximativement un corps noir.
Propriétés du corps noir :
 Tous les corps noirs rayonnent de la même manière.
 Le corps noir rayonne plus que le corps non noir à la même température.
 Les corps noirs sont des corps Lambertiens.
2.2.5.2. Corps gris
Un corps gris est un corps dont le pouvoir absorbant 𝛼𝜆𝑇 est indépendant de la longueur d’onde 𝜆 du
rayonnement qu’il reçoit. Il est défini par : 𝛼𝜆𝑇 = 𝛼𝑇 .
En général, on considère les corps solides comme des corps gris par intervalle et on utilise un pouvoir absorbant
moyen vis-à-vis du rayonnement émis pour 𝜆 < 3 𝜇𝑚 (rayonnement émis par des corps à haute température
comme le Soleil) et un pouvoir absorbant moyen vis-à-vis du rayonnement émis pour 𝜆 > 3 𝜇𝑚 (rayonnement
émis par les corps à faible température : atmosphère, absorbeur solaire,...). On pourra à titre d’exemple
considérer les valeurs suivantes pour la peinture blanche :
Figure 29 : Représentation simplifiée du pouvoir absorbant monochromatique de la peinture blanche

3. Lois du rayonnement
3.1. Loi de Lambert
Dans le cas où la source est isotrope, la luminance est indépendante de la direction : 𝐿𝑥 = 𝐿
Or
𝐼𝑛
𝐿𝑛 =
𝑆
Et
𝐼𝛼
𝐿𝛼 =
𝑆 cos 𝛼
De l’égalité 𝐿𝑛 = 𝐿𝛼 on déduit la loi de Lambert pour une source isotrope :
𝐼𝛼 = 𝐼𝑛 . cos 𝛼 (4.07)

Figure 30 : Schématisation de l’intensité énergétique

Ainsi l’indicatrice d’émission est une sphère tangente en 0 à la surface émettrice lorsque celle-ci suit la loi de
Lambert :
Figure 31 : Schématisation de la luminance et de l’intensité énergétique d’une source isotrope

Remarque : Comme pour un cône de demi-angle au sommet a : 𝛺 = 2𝜋(1 − cos 𝛼) et 𝑑𝛺 = 2𝜋. sin 𝛼 . 𝑑𝛼,
lorsqu’un corps suit la loi de Lambert :
𝛼=𝜋/2 𝛼=𝜋/2
𝑑𝛷
𝑀= =𝐿 ∫ cos 𝛼 . 𝑑𝛺 = 2𝜋𝐿 ∫ cos 𝛼 . sin 𝛼 . 𝑑𝛼
𝑑𝑆
𝛼=0 𝛼=0
Soit :
𝑀 = 𝜋𝐿 [𝑊/𝑚2 ] (4.08)
3.2. Lois physiques
3.2.1. Loi de Kirchhoff
𝑀𝜆𝑇
A une température 𝑇 donnée et pour une longueur d’onde 𝜆 donnée, le rapport 𝛼𝜆𝑇
est le même pour tous les
corps.
Pour le corps noir : 𝛼𝜆𝑇 = 1, on obtient :
𝑀𝜆𝑇 = 𝛼𝜆𝑇 𝑀0𝜆𝑇 [𝑊/𝑚3 ] (4.09)
Où 𝑀0𝜆𝑇 l’émittance monochromatique du corps noir.
On peut généraliser cette loi de Kirchoff au cas du corps gris. En effet pour un corps gris 𝛼𝜆𝑇 = 𝛼𝑇 on a :
𝜆=∞ 𝜆=∞ 𝜆=∞
𝑀𝑇 = ∫ 𝑀𝜆𝑇 . 𝑑𝜆 = ∫ 𝛼𝜆𝑇 𝑀0𝜆𝑇 . 𝑑𝜆 = 𝛼𝑇 ∫ 𝑀0𝜆𝑇 . 𝑑𝜆

𝜆=0 𝜆=0 𝜆=0
En appelant 𝑀0𝑇 l’émittance totale du corps noir à la température 𝑇, nous obtenons pour un corps gris :
𝑀𝑇 = 𝛼𝑇 𝑀0𝑇 [𝑊/𝑚2 ] (4.10)
3.2.2. Rayonnement du corps noir
3.2.2.1. Emittance monochromatique
Elle est donnée par la loi de Planck :

𝐶1 𝜆−5
𝑀0𝜆𝑇 = 𝐶
(4.11)
𝑒𝑥𝑝 (𝜆𝑇2 ) − 1
Avec 𝐶1 = 3,742. 10−16 𝑊. 𝑚−2 et 𝐶2 = 1,4385. 10−2 𝑚. 𝐾

La loi de Planck permet de tracer les courbes isothermes représentant les variations de 𝑀0𝜆𝑇 en fonction de la
longueur d’onde pour diverses températures :
Figure 32 : Emittance monochromatique d’un corps noir à deux températures différentes

La longueur d’onde 𝜆𝑚𝑎𝑥 pour laquelle l’émission est maximale varie avec la température de la source :
2897,8
𝜆𝑚𝑎𝑥 = [𝜇𝑚. 𝐾] (4.12)
𝑇
Cette relation (4.12) est appelée la première loi de Wien.
Portée dans l’expression de l’émittance, elle donne la deuxième loi de Wien qui donne le maximum de
l’émittance hémisphérique dans le spectre d’émission :
𝑇 5
𝑀0𝜆𝑚𝑎𝑥 = 0,410 ( ) (4.13)
10
Avec 𝑇 : Température [𝐾].
3.2.2.2. Emittance du corps noir 𝑬𝟎
L’intégration de la loi de Planck pour toutes les longueurs d’onde donne la loi Stefan-Boltzmann reliant
l’émittance totale ou densité de puissance 𝑀0 à la température 𝑇 de la surface du corps noir :
∞
𝑀0 = ∫ 𝑀0𝜆 (𝜆, 𝑇). 𝑑𝜆 = 𝜎𝑇 4 [𝑊/𝑚2] (4.14)

0
Où 𝜎 est la constante de Stefan-Boltzmann 𝜎 = 5,67. 10−8 𝑊/𝑚2𝐾 4

Dans les calculs on écrira souvent :
𝑇 4
𝑀0 = 5,67 ( )
100
3.2.2.3. Fraction de l’émittance dans un intervalle donné de longueurs d’onde [λ 1, λ2]
C’est la fraction du flux émis par l’unité de surface du corps noir à la température 𝑇 entre les longueurs d’ondes
𝜆1 et 𝜆2 :

𝜆 𝜆2 𝜆2 𝜆1
∫𝜆 2 𝑀0𝜆𝑇 . 𝑑𝜆
1 1
1
𝐹𝜆1𝑇− 𝜆2𝑇 = ∞ = 4
∫ 𝑀0𝜆𝑇 . 𝑑𝜆 = (∫ 𝑀0𝜆𝑇 . 𝑑𝜆 − ∫ 𝑀0𝜆𝑇 . 𝑑𝜆)
∫0 𝑀0𝜆𝑇 . 𝑑𝜆 𝜎𝑇 𝜎𝑇 4
𝜆1 0 0
𝜆2 𝜆1
1 1
= ∫ 𝑀0𝜆𝑇 . 𝑑𝜆 − 4 ∫ 𝑀0𝜆𝑇 . 𝑑𝜆
𝜎𝑇 4 𝜎𝑇
0 0
Ce qui peut également s’écrire :

𝐹𝜆1𝑇− 𝜆2𝑇 = 𝐹0− 𝜆2𝑇 − 𝐹0− 𝜆1𝑇
3.2.3. Rayonnement des corps non noirs
3.2.3.1. Facteur d’émission ou émissivité
On définit les propriétés émissives des corps réels par rapport aux propriétés émissives du corps noir dans les
mêmes conditions de température et de longueur d’onde et on les caractérise à l’aide de coefficients appelés
facteurs d’émission ou émissivités. Ces coefficients monochromatiques ou totaux sont définis par :
𝑀𝜆𝑇 𝑀𝑇
𝜀𝜆𝑇 = 𝑒𝑡 𝜀𝑇 = (4.15)
𝑀0𝜆𝑇 𝑀0𝑇
D’après la loi de Kirchoff, on montre que :
𝛼𝜆𝑇 = 𝜀𝜆𝑇 (4.16)
3.2.3.2. Cas des corps gris
Ils sont caractérisés par 𝛼𝜆𝑇 = 𝛼𝑇 soit d’après ce qui précède : 𝜀𝜆𝑇 = 𝜀𝑇
Or : 𝑀𝑇 = 𝜀𝑇 𝑀0𝑇, nous en déduisons l’émittance du corps gris à la température 𝑇 :
𝑀𝑇 = 𝜀𝑇 𝜎𝑇 4 [𝑊/𝑚2] (4.17)
4. Rayonnement réciproque de plusieurs surfaces
Hypothèses :
 Les surfaces considérées sont supposées homogènes, opaques, isothermes et grises.
 Les éclairements sont supposés homogènes et les réflexions diffuses.
4.1. Radiosité et flux net perdu
Le rayonnement qui quitte une surface 𝑆𝑖 est la somme de son émission propre et de la réflexion d’une partie du
rayonnement incident sur cette surface. On appelle radiosité, que l’on note 𝐽𝑖 , l’émittance apparente de la surface
𝑆𝑖 donc :
𝐽𝑖 = 𝜀𝑖 𝜎𝑇𝑖4 + (1 − 𝜀𝑖 )𝐸𝑖 [𝑊/𝑚2 ] (4.18)
Avec 𝐸𝑖 : Eclairement de la surface 𝑆𝑖 (𝑊/𝑚2)
Considérons maintenant la surface 𝑆𝑖 choisie parmi 𝑛 surfaces isothermes et homogènes qui délimitent un
volume :
Figure 33 : Schématisation des flux de rayonnement sur une surface

La densité d’énergie nette perdue par rayonnement par 𝑆𝑖 s’écrit :

𝜑𝑖𝑛𝑒𝑡 = 𝜀𝑖 𝜎𝑇𝑖4 − 𝜀𝑖 𝐸𝑖
En introduisant, d’après (4.18), la radiosité 𝐽𝑖 par :
1
𝐸𝑖 = (𝐽 − 𝜀𝑖 𝜎𝑇𝑖4 )
1 − 𝜀𝑖 𝑖
Nous obtenons :
𝜀𝑖
𝜑𝑖𝑛𝑒𝑡 = (𝜎𝑇𝑖4 − 𝐽𝑖 ) = 𝜀𝑖 (𝜎𝑇𝑖4 − 𝐸𝑖 )[𝑊/𝑚2 ] (4.19)
1 − 𝜀𝑖
4.2. Facteur de forme géométrique
On considère une surface 𝑆𝑖 qui sur toute son étendue a une émission apparente 𝛷𝑖 = 𝑆𝑖 𝐽𝑖 .
La surface 𝑆𝑖 est environnée par un nombre 𝑛 de surfaces et 𝛷𝑖 est envoyé sur toutes ces surfaces (la surface 𝑆𝑖
peut également rayonner vers elle-même si elle est concave). Le flux apparent 𝛷𝑖 peut donc se décomposer de la
manière suivante :
𝛷𝑖 = 𝛷𝑖→1 + 𝛷𝑖→2 + ⋯ + 𝛷𝑖→𝑖 + ⋯ + 𝛷𝑖→𝑛
Calculons 𝛷𝑖→𝑘 qui est la part du flux quittant 𝑆𝑖 qui atteint 𝑆𝑘 .
D’après la formule de Bougouer, le flux 𝑑2 𝛷𝑖→𝑘 envoyé par la surface élémentaire 𝑑𝑆𝑖 vers la surface
élémentaire 𝑑𝑆𝑘 s’écrit :
𝑑𝑆𝑖 . cos 𝛼𝑖 . 𝑑𝑆𝑘 . cos 𝛼𝑘
𝑑2 𝛷𝑖→𝑘 = 𝐿𝑖
𝑟2
Avec 𝐿𝑖 = 𝐽𝑖 /𝜋 comme la surface grise 𝑆𝑖 suit la loi de Lambert.
Nous en déduisons :
cos 𝛼𝑖 . cos 𝛼𝑘
𝛷𝑖→𝑘 = 𝐽𝑖 ∫ ∫ 𝑑𝑆𝑖 . 𝑑𝑆𝑘
𝜋𝑟 2
𝑆𝑖 𝑆𝑘
Le facteur de forme géométrique 𝑓𝑖𝑘 de la surface 𝑆𝑖 par rapport à la surface 𝑆𝑘 est alors défini par la relation :
cos 𝛼𝑖 . cos 𝛼𝑘
𝑆𝑖 𝑓𝑖𝑘 = ∫ ∫ 𝑑𝑆𝑖 . 𝑑𝑆𝑘 (4.20)
𝜋𝑟 2
𝑆𝑖 𝑆𝑘
Il ne dépend que de la géométrie et de la disposition relative des surfaces 𝑆𝑖 et 𝑆𝑘 .

Le flux 𝛷𝑖→𝑘 peut alors s’écrire simplement :
𝛷𝑖→𝑘 = 𝐽𝑖 𝑓𝑖𝑘 𝑆𝑖
Le facteur de forme géométrique 𝑓𝑖𝑘 s’interprète simplement comme la fraction du flux total émis en apparence
par 𝑆𝑖 (𝛷𝑖 = 𝐽𝑖 𝑆𝑖 ) qui atteint la surface 𝑆𝑗 .
Remarques :
 Le 2ème membre de la formule (4.20) de définition de 𝑓𝑖𝑘 est symétrique en 𝑖 et 𝑘, on en déduit la
relation de réciprocité des facteurs de forme :
𝑆𝑖 𝑓𝑖𝑘 = 𝑆𝑘 𝑓𝑘𝑖 (4.21)
 La relation 𝛷𝑖 = 𝛷𝑖→1 + 𝛷𝑖→2 + ⋯ + 𝛷𝑖→𝑖 + ⋯ + 𝛷𝑖→𝑛 peut s’écrire :
𝛷𝑖 = 𝐽𝑖 𝑓𝑖1 𝑆𝑖 + 𝐽𝑖 𝑓𝑖2 𝑆𝑖 + ⋯ + 𝐽𝑖 𝑓𝑖𝑖 𝑆𝑖 + ⋯ + 𝐽𝑖 𝑓𝑖𝑛 𝑆𝑖 = 𝐽𝑖 𝑆𝑖 (𝑓𝑖1 + 𝑓𝑖2 + ⋯ + 𝑓𝑖𝑖 + ⋯ + 𝑓𝑖𝑛 )
Or 𝛷𝑖 = 𝐽𝑖 𝑆𝑖
D’où :
𝑓𝑖1 + 𝑓𝑖2 + ⋯ + 𝑓𝑖𝑖 + ⋯ + 𝑓𝑖𝑛 = 1 (4.22)
Ces deux relations sont utiles pour la détermination des facteurs de formes de plusieurs surfaces en présence.

4.3. Calcul des flux

Le flux 𝛷→𝑖 reçu par la surface 𝑆𝑖 s’écrit :
𝑛
𝛷→𝑖 = 𝐸𝑖 𝑆𝑖 = ∑ 𝛷𝑘→𝑖
𝑘=1
Or 𝛷𝑘→𝑖 = 𝐽𝑘 𝑆𝑘 𝑓𝑘𝑖
D’où :
𝑛 𝑛
𝐸𝑖 𝑆𝑖 = ∑ 𝐽𝑘 𝑆𝑘 𝑓𝑘𝑖 = ∑ 𝐽𝑘 𝑆𝑖 𝑓𝑖𝑘
𝑘=1 𝑘=1
D’après (4.21).
En reportant cette expression dans (4.18), nous obtenons :
𝑛
𝐽𝑖 = 𝜀𝑖 𝜎𝑇𝑖4 + (1 − 𝜀𝑖 ) ∑ 𝐽𝑘 𝑓𝑖𝑘
𝑘=1
Soit encore :
𝑛
𝐽𝑖 1
𝜎𝑇𝑖4 = − ∑(1 − 𝜀𝑖 )𝐽𝑘 𝑓𝑖𝑘
𝜀𝑖 𝜀𝑖
𝑘=1
En utilisant le symbole de Kronecker, nous pouvons écrire :

𝑛
𝐽𝑖 = ∑ 𝛿𝑖𝑘 𝐽𝑘
𝑘=1
D’où :
𝑛
𝜀𝑖 𝜎𝑇𝑖4 = ∑[𝛿𝑖𝑘 − (1 − 𝜀𝑖 )𝑓𝑖𝑘 ]𝐽𝑘 (4.23)

𝑘=1
On écrit cette relation pour toutes les surfaces 𝑆𝑖 dont on connaît les températures. Pour celles dont on connaît
plutôt la densité de flux net perdue 𝜑𝑖𝑛𝑒𝑡 on utilise la relation :
𝑛
𝜑𝑖𝑛𝑒𝑡 = 𝐽𝑖 − 𝐸𝑖 = 𝐽𝑖 − ∑ 𝑓𝑖𝑘 𝐽𝑘
𝑘=1
Qui peut encore s’écrire :

𝑛
𝜑𝑖𝑛𝑒𝑡 = ∑(𝛿𝑖𝑘 − 𝑓𝑖𝑘 )𝐽𝑘 (4.24)

𝑘=1
4.4. Méthode de résolution

Si l’on connaît 𝑝 températures et (𝑛 − 𝑝) densités de flux nets 𝜑𝑖𝑛𝑒𝑡 , on écrit 𝑝 fois l’équation (4.23) et (𝑛 − 𝑝)
fois l’équation (4.24), on obtient ainsi un système linéaire de 𝑛 équations à 𝑛 inconnues : 𝐽1, 𝐽2,…, 𝐽𝑝 ,
𝑇𝑝+1 ,…,𝑇𝑛 .
La résolution de ce système permet de calculer les (𝑛 − 𝑝) températures et les 𝑝 radiosités inconnues. Les 𝑝
densités de flux nets inconnues se calculent ensuite par la relation :
𝜀𝑖
𝜑𝑖𝑛𝑒𝑡 = (𝜎𝑇𝑖4 − 𝐽𝑖 )
1 − 𝜀𝑖
Remarque :

Si une surface est noire (𝜀𝑖 = 1), la relation (4.23) ne peut pas être utilisée. Nous avons alors simplement dans ce
cas la relation 𝐽𝑖 = 𝜎𝑇𝑖4 et l’on résout le système des (𝑛 − 1) équations restantes.
4.5. Exemple d’application : Cas de deux plans parallèles infinis
On suppose que les températures 𝑇1 et 𝑇2 ainsi que les émissivités 𝜀1 et 𝜀2 des deux surfaces 𝑆1 et 𝑆2 sont
connues, on cherche à déterminer le flux net perdu par chacune de ces surfaces.
Nous avons 𝑓11 = 𝑓22 = 0 car les surfaces 𝑆1 et 𝑆2 sont planes et ne peuvent pas rayonner vers elles-mêmes.
Nous en déduisons 𝑓12 = 1 et 𝑓21 = 1 en appliquant la relation ∑𝑛𝑘=1 𝑓𝑖𝑘 = 1 pour 𝑖 = 1 et pour 𝑖 = 2.
La relation (4.23) s’écrit alors de la manière suivante pour i = 1 et i = 2 :
𝐽 − (1 − 𝜀1 )𝐽2 = 𝜀1 𝜎𝑇14
{ 1
−(1 − 𝜀2 )𝐽1 + 𝐽2 = 𝜀2 𝜎𝑇24
D’où :
𝜀1 𝑇14 + 𝜀2 (1 − 𝜀1 )𝑇24
𝐽1 = 𝜎
1 − (1 − 𝜀1 )(1 − 𝜀2 )
Et
𝜀1 𝜀1 𝜀1 𝜀1 𝜀1 𝜀2 (1 − 𝜀1 )
𝜑1𝑛𝑒𝑡 = (𝜎𝑇14 − 𝐽1 ) = 𝜎𝑇14 ( − ) − 𝜎𝑇24
1 − 𝜀1 1 − 𝜀1 1 − 𝜀1 𝜀1 + 𝜀2 − 𝜀1 𝜀2 1 − 𝜀1 𝜀1 + 𝜀2 − 𝜀1 𝜀2
𝜀1 𝜀2
𝜑1𝑛𝑒𝑡 = −𝜑2𝑛𝑒𝑡 = 𝜎(𝑇14 − 𝑇24 )
𝜀1 + 𝜀2 − 𝜀1 𝜀2
Soit finalement :
𝜎(𝑇14 − 𝑇24 )
𝜑1𝑛𝑒𝑡 = −𝜑2𝑛𝑒𝑡 = 1 1 (4.25)
+𝜀 −1
𝜀1 2
4.6. Analogie électrique

4.6.1. Flux net perdu par une surface
𝜀
Nous avons montré que 𝜑𝑖𝑛𝑒𝑡 = 1−𝜀𝑖 (𝜎𝑇𝑖4 − 𝐽𝑖 ) ce qui peut encore s’écrire :
𝑖
𝜎𝑇𝑖4 − 𝐽𝑖
𝛷𝑖𝑛𝑒𝑡 = 1−𝜀𝑖
𝜀𝑖𝑆𝑖
Par analogie, cette relation peut être représentée par le schéma électrique équivalent suivant :
Figure 34 : Schéma électrique équivalent du flux radiatif perdu par une surface
On notera que cette résistance thermique de rayonnement ne dépend que des propriétés physiques de la surface
𝑆𝑖 et qu’elle est nulle pour un corps noir.
4.6.2. Flux net échangé entre plusieurs surfaces
Le flux net perdu par la surface 𝑆𝑖 dans ses échanges radiatifs avec l’ensemble des surfaces environnantes s’écrit
d’après la relation (4.19) : 𝛷𝑖𝑛𝑒𝑡 = (𝐽𝑖 − 𝐸𝑖 )𝑆𝑖 .
Le flux 𝛷𝑖 = 𝐽𝑖 𝑆𝑖 quittant la surface 𝑆𝑖 peut se décomposer de la manière suivante :
𝑛
𝛷𝑖 = 𝛷𝑖→1 + 𝛷𝑖→2 + ⋯ + 𝛷𝑖→𝑛 = ∑ 𝐽𝑖 𝑆𝑖 𝑓𝑖𝑗

𝑗=1

L’éclairement 𝐸𝑖 reçu par la surface 𝑆𝑖 peut se décomposer de la manière suivante :

𝑛 𝑛
𝐸𝑖 𝑆𝑖 = ∑ 𝛷𝑗→𝑖 = ∑ 𝐽𝑗 𝑆𝑗 𝑓𝑗𝑖
𝑗=1 𝑗=1
Le flux net perdu par 𝑆𝑖 peut donc s’écrire :

𝑛 𝑛 𝑛 𝑛
𝛷𝑖𝑛𝑒𝑡 = ∑ 𝐽𝑖 𝑆𝑖 𝑓𝑖𝑗 − ∑ 𝐽𝑗 𝑆𝑗 𝑓𝑗𝑖 = ∑ 𝑆𝑖 𝑓𝑖𝑗 (𝐽𝑖 − 𝐽𝑗 ) = ∑ 𝛷𝑛𝑒𝑡𝑖→𝑗

𝑗=1 𝑗=1 𝑗=1 𝑗=1
Le flux net échangé entre les surfaces 𝑆𝑖 et 𝑆𝑗 s’écrit donc :

𝐽𝑖 − 𝐽𝑗
𝛷𝑛𝑒𝑡𝑖→𝑗 = (𝐽𝑖 − 𝐽𝑗 )𝑆𝑖 𝑓𝑖𝑗 = 1
𝑆𝑖 𝑓𝑖𝑗
Cet échange radiatif peut être représenté par le schéma électrique équivalent suivant :
Figure 35 : Schéma électrique équivalent du flux radiatif échangé entre deux surfaces
On notera que cette résistance thermique de rayonnement est purement géométrique et qu’elle ne dépend pas des
propriétés physiques des surfaces 𝑆𝑖 et 𝑆𝑗 .
4.6.3. Application : Echange entre deux surfaces grises
Si les deux surfaces 𝑆1 et 𝑆2 sont seules en présence, le flux net 𝛷1𝑛𝑒𝑡 perdu par 𝑆1 est égal au flux net 𝛷2𝑛𝑒𝑡
gagné par 𝑆2 . Ce flux est encore égal au flux net 𝛷𝑛𝑒𝑡1→2 échangé entre 𝑆1 et 𝑆2 , nous avons donc les égalités :
𝛷1𝑛𝑒𝑡 = 𝛷𝑛𝑒𝑡1→2 = −𝛷2𝑛𝑒𝑡
Soit :
𝜎𝑇14 − 𝐽1 𝐽1 − 𝐽2 𝐽2 − 𝜎𝑇24
𝛷1𝑛𝑒𝑡 = 1−𝜀1 = 1 = 1−𝜀2
𝜀1 𝑆1 𝑆1 𝑓12 𝜀2 𝑆2
Cet échange radiatif peut être représenté par le schéma électrique équivalent suivant :
Figure 36 : Schéma électrique équivalent du flux radiatif net échangé entre deux surfaces
D’où
𝜎(𝑇14 − 𝑇24 )
𝛷1𝑛𝑒𝑡 = −𝛷2𝑛𝑒𝑡 = 1−𝜀1 1 1−𝜀 [𝑊] (4.26)
+ + 2
𝜀1 𝑆1 𝑆1 𝑓12 𝜀2 𝑆2
4.6.4. Utilisation des schémas analogiques

Dans les systèmes simples, il est plus rapide d’utiliser la technique des schémas analogiques que celle du
système linéaire. Lorsqu’on a établi le schéma analogique, on calcule les différentes résistances du circuit puis
on résout par les techniques habituelles utilisées en électricité : loi d’association des résistances en série et en
parallèle, loi des nœuds,...

4.6.5. Exemple d’application

4.6.5.1. Cas d’une surface 𝑺𝟏 convexe complètement entourée par une surface 𝑺𝟐
La surface 𝑆1 étant convexe elle ne peut pas rayonner vers elle-même donc : 𝑓11 = 0.
La relation 𝑓11 + 𝑓12 = 1 nous permet de déduire : 𝑓12 = 1
La relation (4.26) s’écrit alors :
𝜎(𝑇14 − 𝑇24 ) 𝜎(𝑇14 − 𝑇24 )
𝛷1𝑛𝑒𝑡 = −𝛷2𝑛𝑒𝑡 = 1−𝜀1 1 1−𝜀2 = 1 1 1
+ + + −
𝜀1 𝑆1 𝑆1 𝜀2 𝑆2 𝜀1 𝑆1 𝜀2 𝑆2 𝑆2
D’où :
𝜎𝑆1 (𝑇14 − 𝑇24 )
𝛷1𝑛𝑒𝑡 = −𝛷2𝑛𝑒𝑡 = 1 𝑆1 1
[𝑊] (4.27)
+ ( − 1)
𝜀1 𝑆2 𝜀2
4.6.5.2. Cas particulier où la surface 𝑺𝟏 est « petite » devant la surface 𝑺𝟐

𝑆1
Nous avons dans ce cas : ≈ 0 et la relation (4.27) s’écrit alors :
𝑆2
𝛷1𝑛𝑒𝑡 = −𝛷2𝑛𝑒𝑡 = 𝜎𝜀1 𝑆1 (𝑇14 − 𝑇24 )[𝑊] (4.28)

Vue de la surface S1 (« petit corps »), la surface S2 se comporte alors comme un corps noir.

Série d’exercices
Transfert de Chaleur par Conduction
Exercice 1
Le mur d’un bâtiment est fait de briques ayant une épaisseur 𝛿 = 38 𝑐𝑚 et la conductivité thermique 𝜆 =
0,78 𝑊/𝑚𝐾. Les températures sur les faces limitatrices du mur sont égales à 38°𝐶 respectivement 25°𝐶.
Déterminer :
a) le flux surfacique qui traverse le mur ;
b) la profondeur du mur où la température est 30°𝐶.
Exercice 2
Les apports thermiques à travers un mur en béton de 30 m2 de surface est 690 W. Sachant que le mur a une
épaisseur de 10 cm, et que la température de sa face intérieure est 25°C, calculer la température de la face
extérieure supposée supérieure à Ti.
On donne : 𝜆𝑏é𝑡𝑜𝑛 = 1,75 𝑊/𝑚𝐾.
Exercice 3
Le mur d’un bâtiment est fait de briques rouges ayant une épaisseur de 38 𝑐𝑚 et une conductivité thermique 𝜆 =
0,8 𝑊/𝑚𝐾. La température de l’air intérieur est de 40°𝐶 et celle de l’air extérieur 25°𝐶. Les coefficients de
transfert superficiel correspondants sont ℎ1 = 10 𝑊/𝑚2𝐾 et ℎ2 = 20 𝑊/𝑚2 𝐾.
a) Calculer le coefficient de transfert thermique global.
b) Calculer la résistance thermique totale du mur.
c) Calculer le flux thermique surfacique.
d) Calculer les températures sur les faces limitatrices du mur.
Exercice 4
Le circuit d’une surchauffeur à vapeurs est constitué d’un tuyau d’acier thermorésistant ayant le diamètre 42 ∗
5 𝑚𝑚 et une conductivité thermique 𝜆 = 16,3 𝑊/𝑚𝐾. La température de la surface extérieure du circuit est de
580 °𝐶 et celle de la surface intérieure 450 °𝐶.
Déterminer le flux thermique conductif rapporté par mètre linéaire de tuyau.
Exercice 5
Une conduite cylindrique, de rayons interne 6 cm et externe 8 cm, à une conductivité thermique de 0,5 𝑊/𝑚𝐾
supposée indépendante de la température. La température de la surface interne est de 430°C et celle de la face
externe de 30°C. Calculer le flux de chaleur perdu par mètre de longueur de la conduite.
Exercice 6
Un préchauffeur d’air est fait de tuyau d’acier ayant le diamètre 49 ∗ 3 𝑚𝑚 et la conductivité thermique 𝜆 =
50 𝑊/𝑚𝐾. A l’intérieur du tuyau circulent des gaz brûlés ayant la température 𝜃𝑓1 = 250°𝐶 et à l’extérieur un
courant d’air ayant la température 𝜃𝑓2 = 145 𝐶. Les coefficients de transfert thermique surfacique sont à
l’intérieur ℎ1 = 45 𝑊/𝑚2𝐾 et à l’intérieur ℎ2 = 25 𝑊/𝑚2 𝐾.
Déterminer :
a) le flux thermique transféré par mètre linéaire de tuyau ;
b) le coefficient de transfert thermique global linéaire ;
c) les températures sur les surfaces du tuyau.

Exercice 7
Le mur d’un bâtiment est fait de briques ayant une épaisseur 𝛿1 = 38 𝑐𝑚 et la conductivité thermique 𝜆1 =
0,78 𝑊/𝑚𝐾 sur lequel on a appliqué deux couches de crépis sur ses faces limitatrices, d’épaisseur respective
𝛿2 = 2 𝑐𝑚 et de conductivité thermique 𝜆2 = 0,7 𝑊/𝑚𝐾. Les températures sur les faces extérieure et intérieure
sont 20 °𝐶 et respectivement −10 °𝐶.
a) Calculer le flux surfacique qui traverse le mur.
b) Calculer la distribution de température à l’intérieur du mur.
Exercice 8
Le mur d’un four est composé de trois couches : la première couche est réfractaire avec une épaisseur de
150 𝑚𝑚 et une conductivité thermique 𝜆1 = 0,81 𝑊/𝑚𝐾 ; la deuxième couche est faite en diatomite avec une
épaisseur de 95 𝑚𝑚 et 𝜆2 = 0,3 𝑊/𝑚𝐾 et la troisième couche est faite en brique rouge avec une épaisseur de
250 𝑚𝑚 et 𝜆3 = 0,7 𝑊/𝑚𝐾. La température du gaz à l’intérieur du four est de 1200 °𝐶 et la température de
l’air à l’intérieur de la halle est de 30 °𝐶. Les coefficients de transfert thermique superficiel sont ℎ1 =
35 𝑊/𝑚2𝐾 et ℎ2 = 12 𝑊/𝑚2 𝐾.
Déterminer :
a) la résistance thermique totale ;
b) le coefficient de transfert thermique global ;
c) le flux surfacique transféré ;
d) les températures des surfaces limitatrices du mur et dans le plan de contact entre les couches ;
e) la conductivité thermique équivalente du mur ;
f) l’épaisseur que le mur doit avoir, s’il était construit seulement en brique rouge.
Exercice 9
Une conduite de vapeur ayant un diamètre 160 ∗ 5 𝑚𝑚 et une conductivité thermique 𝜆1 = 50 𝑊/𝑚𝐾 est
couverte par une couche d’isolation thermique ayant une épaisseur 𝛿2 = 100 𝑚𝑚 et une conductivité thermique
𝜆2 = 0,08 𝑊/𝑚𝐾. On connaît les températures sur la surface intérieure de la conduite 𝜃1 = 400 °𝐶 et celle sur
la surface de l’isolation 𝜃3 = 50 °𝐶.
Déterminer :
a) le flux thermique linéaire ;
b) la température dans la surface de contact entre la conduite et l’isolation.

Transfert de Chaleur par Convection

Exercice 1
Certains aménagements urbains à revêtements de bitume, tels que les rues, les parcs de stationnement et les toits
d’immeuble, absorbent très bien le rayonnement solaire et leur température de surface augmente pendant le jour
(les espaces verts absorbent aussi le rayonnement solaire mais l’évapotranspiration évite l’accroissement de
température). En conséquence, la température globale au niveau du sol s’élève au cours de la journée. Ce
phénomène entraîne en fin de journée (en été en particulier) le développement de mouvement de convection
naturelle.
On donne : différence de température entre le sol et l’atmosphère de 10°C, une échelle de longueur typique 𝐿 =
500 𝑚, température atmosphérique 𝑇𝑎 = 293 𝐾, 𝜈 = 1,5. 10−5 𝑚2 /𝑠 et 𝑃𝑟 = 0,72.
a. Estimer la vitesse de convection naturelle.
b. Déterminer le flux thermique surfacique.
Corrélation :
Avec :
 A = 0,54, 𝑚 = 0,25 si 𝑅𝑎 est tel que : 2. 104 < 𝑅𝑎 < 8. 106 (régime laminaire) ;
 A = 0,15, 𝑚 = 0,33 si 𝑅𝑎 est tel que : 8. 106 < 𝑅𝑎 (régime turbulent).
Exercice 2
Evaluer le flux de chaleur transféré dans l’air par un radiateur simple constitué par une plaque verticale de
hauteur 𝐿 = 1𝑚 pour les différentes corrélations suivantes. Commenter les résultats.
On donne : 𝑇𝑎 = 293𝐾 et 𝑇𝑝 = 333𝐾 les températures de l’air et de la paroi du radiateur et les propriétés
physiques de l’air : 𝑃𝑟 = 0,71 ; 𝜈 = 1,72. 10−5 𝑚2 /𝑠 ; 𝑎 = 2,4. 10−5 𝑚2 /𝑠 ; 𝜆 = 0,027𝑊. 𝑚−1 . 𝐾 −1 ; 𝛽 =
3,2. 10−3 𝐾 −1.
Corrélation pour 𝑅𝑎 > 109 écoulement turbulent
 Burmeister :
2 1
0,23. 𝑅𝑎5 . 𝑃𝑟 15
𝑁𝑢 = 2
2 5
(1 + 0,494. 𝑃𝑟 3 )
 Churchill et Chu :
2
1
0,387. 𝑅𝑎 6
𝑁𝑢 = 0,825 + 8
9 27
[ (1 + 0,494. (2,033. 𝑃𝑟)−16 ) ]
 Bayley :
𝑁𝑢 = 0,183. 𝑅𝑎0,31
Exercice 3
On réalise le chauffage de l’eau d’un réservoir d’un chauffe eau à la température de 𝑇𝑎 = 303 𝐾 à l’aide de
conduites horizontales de diamètre extérieur 𝑑 = 25 𝑚𝑚. La température externe des conduites est égale à 𝑇𝑐 =
323 𝐾. On suppose que l’espacement des conduites est grand.

Pour la température de référence (𝑇𝑎 + 𝑇𝑐 )/2, les propriétés physiques de l’eau sont : 𝑃𝑟 = 4,31 ; 𝜈 =
0,66. 10−6 𝑚2 /𝑠 ; 𝑎 = 0,15. 10−6 𝑚2 /𝑠 ; 𝜆 = 0,635 𝑊. 𝑚−1 . 𝐾 −1 ; 𝛽 = 3,87. 10−4 𝐾 −1 .
Calculer le flux de chaleur cédé par la conduite chauffée au réservoir par convection naturelle pour une longueur
de conduite 𝑙 = 1 𝑚.
Exercice 4
On veut calculer le flux de chaleur échangé entre de l'air à 5°C
circulant devant une paroi plane à 71°C à une vitesse de 26,8 m/s. On
admet que la température moyenne du film d'air au voisinage de la
paroi est 38°C. Les caractéristiques correspondantes de l'air sont :
𝜌 = 1,136𝑘𝑔/𝑚3, 𝑐𝑝 = 1000𝐽/𝑘𝑔𝐾, 𝜇 = 1,91. 10−5 𝑃𝑎. 𝑠, 𝜆=
0,027𝑊/𝑚𝐾.
Exercice 5
Calculer le flux de chaleur dégagée par une lampe sphérique de diamètre 50 𝑚𝑚 à incandescence de 40𝑊 qui a
une température moyenne de 60°𝐶. La lampe est exposée sous un air de 27°𝐶 circulant à une vitesse de 0,3𝑚/𝑠.
Quel est le pourcentage de chaleur perdue par :
1. Convection forcée (𝑁𝑢 = 0,37(𝑅𝑒)0,6 ).
2. Convection naturelle (𝑁𝑢 = 0,37(𝐺𝑟. 𝑃𝑟)0,25 ).
3. Qu’en déduisez-vous ?
Caractéristiques de l’air à 43,5°𝐶 : 𝐶𝑝 = 1008𝐽. 𝑘𝑔 −1 . 𝐾 −1 ; 𝑎 = 2,98. 10−5 𝑚2 /𝑠 ; 𝜆 = 0,03𝑊. 𝑚−1 . 𝐾 −1 ; 𝜌 =
1𝑘𝑔. 𝑚−3 ; 𝜈 = 2,079. 10−5 𝑚2 /𝑠.

Exercie 6
Dans un canal de section carrée de coté 𝑎 = 25𝑚𝑚, de longueur 𝑙 = 1,6𝑚, circule de l’eau à la vitesse 𝑈 =
1𝑚/𝑠 et à la température moyenne de 313𝐾. La température de la paroi du canal est 𝑇𝑝 = 363𝐾. Les propriétés
physiques de l’eau sont : 𝜆 = 63,5. 10−2 𝑊. 𝑚−1 . 𝐾 −1 ; 𝜈 = 0,66. 10−6 𝑚2/𝑠 ; 𝑃𝑟 = 4,31.
Le nombre de Nusselt peut être estimé par les diverses corrélations à disposition :
 Corrélation de Miheev avec le nombre de Prandtl de référence 𝑃𝑟𝑟 = 1,95 :
𝑃𝑟 0,25
𝑁𝑢 = 0,37. 𝑅𝑒 0,8 . 𝑃𝑟 0,43 . ( )
𝑃𝑟𝑟
 Corrélation de Dittus et Boelter :
𝑁𝑢 = 0,23. 𝑅𝑒 0,8 . 𝑃𝑟 0,4
On cherche à calculer le coefficient de transfert thermique convectif de la paroi ℎ. Commenter les deux résultats.

Transfert de Chaleur par Rayonnement

Exercice 1
Si nous assimilons notre soleil à un corps noir de température 𝑇 = 5800 𝐾, alors nous constatons que notre œil
est capteur centré sur les puissance spectrales d’émission maximales.
a) Déterminer le maximum de la longueur d’onde.
b) Déterminer le maximum de l’émittance hémisphérique monochromatique.
Exercice 2
Un corps noir de forme sphérique, ayant un rayon de 5 cm, est maintenu à une température de 327°C.
a) Quelle est la puissance dans son rayonnement ?
b) Quelle est la longueur d’onde correspondant au maximum d’énergie rayonnée.
Exercice 3
Le Soleil est assimilable à un corps noir à une température de 6000 K.
1. Calculer la puissance émise par le Soleil.
2. Calculer l’éclairement solaire zénithal en limite d’atmosphère terrestre.
On donne : Rayon du Soleil 700 000 km et distance Terre-Soleil 150 106 km.
Exercice 4
Soit une surface peinte dont l’émissivité spectrale, variant peu avec la température, est donnée par la figure ci-
dessous. Déterminer l’émittance totale si la température de surface 27 °C.
Exercice 5
On se propose de calculer les flux infrarouge échangé entre deux surfaces opaques verticales étant en regard
uniquement l’une de l’autre. Leurs températures de surfaces respectives sont T1 et T2, et leurs émissivités valent
ε1 et ε2.
1. Donner l’expression de l’émission apparente J de chacune des surfaces vers l’autre. On écrira
l’émission apparente d’une surface comme l’émission propre Φ de cette surface ajoutée à la réflexion
de l’éclairement E de cette surface.
2. Donner l’expression simple du flux échangé φ1->2 entre les deux surfaces, et identifier une émissivité
effective ε.
Exercice 6
D’après les estimations, la surface du soleil est 6,1 .1018 m² et la puissance rayonnée est de 3,9.1026 W. En
supposant que l’émissivité de la surface du soleil est égale à 1, calculer la température de surface de
rayonnement du soleil.

Exercice 7
Pour chauffer une pièce d’un appartement, on se sert d’un radiateur cylindrique de rayon 1 cm et de 60cm de
longueur. Ce radiateur rayonne comme un corps noir et émet une puissance de 1,2 kW.
1. Calculer sa température.
2. Calculer la longueur d’onde pour laquelle l’émittance est maximale.
3. Quelle devrait être sa température pour que cette longueur d’onde soit 2 mm ?
4. Quelle serait alors sa puissance dégagée ?
Exercice 8
1. Une paroi P1 dont la température est de 1 000 K rayonne avec un facteur d’émission totale hémisphérique de
0,8 vers une autre plane P2 parallèle à P1 dont le facteur d’émission est de 0,7. Cette dernière est soumise sur sa
face opposée à un refroidissement provoqué par une circulation d’eau. Le flux de chaleur capté par ce
refroidissement est de 28 kW/m² de surface. Déterminer la température d’équilibre de la paroi P2.
2) Pour éviter une détérioration prématurée de P2 on polit sa face recevant le rayonnement afin de diminuer son
facteur d’émission. Quelle devra être la valeur de ce dernier pour que la température de P2 ne soit plus que de
100°C ?

Annexes
A.1.1 : Fraction d’énergie 𝑭𝟎−𝝀𝑻 rayonnée par un corps noir entre 𝟎 et 𝝀

A.1.2 : Facteurs de forme géométrique de rayonnement



Bibliographie
1. BIANCHI Ana-Maria, FAUTRELLE Yves et ETAY Jacqueline, Transferts Thermiques, Agence Universitaire
de la Francophonie, Presses Polytechniques et Universitaires Romandes, 2004.
2. JANNOT Yves, Transferts Thermiques, Ecole des Mines Nancy, 2ème Année, 2012.
3. Gibert MINGUY, Salif GAYE, Transfert de Chaleur, Cours et Problèmes, 2004

Cours de Transfert Thermique I 2020 2021-1

Transféré par

Droits d'auteur :

Formats disponibles

Cours de Transfert Thermique I 2020 2021-1

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Cours de Transfert Thermique I 2020 2021-1

Transféré par

Droits d'auteur :

Formats disponibles

Institut National Supérieur des Sciences et Techniques d’Abéche

Chargé de Cours et TD : M. DJALDI TABDI NGAMSOU

Table des matières

Transfert Thermique I Page 2

3.3. Expression du flux ......................................................................................................................... 32

Transfert Thermique I Page 3

A.1.2 : Facteurs de forme géométrique de rayonnement ............................................................................... 59

Transfert Thermique I Page 4

Transfert Thermique I Page 5

Chapitre1 Introduction au Transfert Thermique

Transfert Thermique I Page 6

Transfert Thermique I Page 7

Si on considère un système monodirectionnel suivant 0𝑥, on peut écrire :

Figure 1 : Isotherme et gradient de température

Transfert Thermique I Page 8

 φst : densité de flux de chaleur stockés dans le système (S) ;

 φe : densité de flux de chaleur entrant dans le système (S) ;

Figure 2 : Système et bilan thermique

Figure 3 : Schéma du transfert de chaleur par conduction

Transfert Thermique I Page 9

Figure 4 : Schéma du transfert de chaleur par convection

 𝑇∞ : température du fluide loin de la surface [𝐾].

Figure 5 : Schéma du transfert de chaleur par rayonnement

 𝑇𝑆 : température de surface [𝐾] ;

 𝑇1 et 𝑇2 : Températures d’entrée et de sortie [𝐾].

Transfert Thermique I Page 10

Transfert Thermique I Page 11

Chapitre 2 Transfert de Chaleur par Conduction

Figure 6 : Positions relatives du flux surfacique et du gradient de température

Transfert Thermique I Page 12

Où 𝑑𝑆 est un élément de la surface isotherme (figure 7).

Figure 7 : Signification du flux surfacique

 𝜌 : masse volumique en [𝑘𝑔/𝑚3]

Transfert Thermique I Page 13

Transfert Thermique I Page 14

On peut faire une analogie avec la loi d’Ohm en électrotechnique :

Transfert Thermique I Page 15

Figure 9 : Analogie thermoélectrique dans le cas d’un mur plan simple

- 𝑇𝑠 : température de la surface solide [𝐾] ;

En faisant appel à l’analogie thermoélectrique, on a :

Transfert Thermique I Page 16

Transfert Thermique I Page 17

Figure 12 : Cylindre creux à surfaces latérales isothermes

La densité de flux thermique est obtenue à partir de la loi de Fourier :

Il présente une variation radiale de type hyperbolique.

Transfert Thermique I Page 18

Avec 𝑆𝑐 = 2𝜋𝑟𝑙 est la surface latérale du cylindre de rayon 𝑟 et de longueur 𝑙.

3.1.4. Cylindre creux simple en contact avec deux fluides

{𝛷𝑙 = ℎ2 . 2𝜋𝑟2 . (𝑇2 − 𝑇𝑓2 )

Transfert Thermique I Page 19

On obtient par sommation le flux total rapporté à l’unité de longueur du cylindre :

intérieure du cylindre (ou simplement résistance thermique de la surface intérieure) ;

cylindre et le fluide froid (ou résistance thermique de la surface extérieure).

Qui, introduit dans la relation (2.37), conduit à :

Transfert Thermique I Page 20

3.2. Corps ayant des structures multicouches

Figure 15 : Mur composite ayant deux surfaces isothermes

Transfert Thermique I Page 21

Le flux total est :

Où 𝑆 est la surface du mur.