Spécialité : Licence

COURS

TRANSFERTS THERMIQUES

1
TABLE DES MATIÈRES
NOMENCLATURE ............................................................................................................... iv
SIGLES ET ABBREVIATIONS ......................................................................................... vii
LISTE DES TABLEAUX ................................................................................................... viii
LISTE DES FIGURES ......................................................................................................... ix
INTRODUCTION GENERALE ........................................................................................... 1
1. CHAPITRE 1 : INTRODUCTION AUX TRANSFERTS THERMIQUES .............................. 2
INTRODUCTION .................................................................................................................. 3
1.1 Généralités ............................................................................................................................ 3
1.2 Les différents modes de transferts .................................................................................... 4
1.2.1 La conduction ................................................................................................................ 4

1.2.2 Le rayonnement ............................................................................................................ 5

1.2.3 La convection ................................................................................................................ 5

1.3 Combinaison des différents modes de transferts ............................................................ 6

CONCLUSION ...................................................................................................................... 8
2. CHAPITRE 2 : TRANSFERT THERMIQUE PAR CONDUCTION ....................................... 9
INTRODUCTION ................................................................................................................ 10
2.1 Généralités .......................................................................................................................... 10
2.2 La loi de Fourier .................................................................................................................. 11
2.2.1 Définition de quelques grandeurs physiques ......................................................... 11

2.2.2 Enoncé de la loi de Fourier ....................................................................................... 14

2.2.3 Orthogonalité du gradient et de l’isotherme ........................................................... 14

2.3 La conductivité thermique ................................................................................................. 16

2.3.1 Généralités .................................................................................................................. 16

2.3.2 Matériaux anisotropes................................................................................................ 18

2.4 L’équation de la chaleur en conduction .......................................................................... 18

2.4.1 Démonstration ............................................................................................................. 18

2.4.2 Conditions aux limites ................................................................................................ 19

2.4.3 Conductivité et diffusivité thermique ........................................................................ 22

2.4.4 Equation de la chaleur en régime permanent ........................................................ 23

2.5 Exemples de résolution de l’équation de la chaleur en conduction ............................ 24

i
 2.5.1 Cas d’un mur plan ...................................................................................................... 24

2.5.2 Association de murs en série ou en parallèle ........................................................ 25

2.5.3 Mur avec sources internes ........................................................................................ 26

2.5.4 Mur avec conditions de convection aux parois ...................................................... 27

2.5.5 Conduction dans une coque cylindrique ................................................................. 29

2.5.6 Conduction dans une coque sphérique................................................................... 30

2.5.7 Remarques finales sur la conduction ...................................................................... 31

CONCLUSION .................................................................................................................... 32
3. CHAPITRE 3 : TRANSFERT THERMIQUE PAR CONVECTION ..................................... 33
INTRODUCTION ................................................................................................................ 34
3.1 Généralités .......................................................................................................................... 34
3.2 Coefficient d’échange convectif........................................................................................ 36
3.3 Lois de corrélation pour le coefficient d’échange convectif ......................................... 38
3.3.1 Paramètres caractéristiques de la convection ....................................................... 38

3.3.2 Lois de corrélation en convection............................................................................. 42

3.3.3 Méthodologie pour calculer le flux de chaleur en convection .............................. 42

3.3.4 Principales lois de correlation en convection forcée ............................................. 43

CONCLUSION .................................................................................................................... 45
4. CHAPITRE 4 : TRANSFERT THERMIQUE PAR RAYONNEMENT
ELECTROMAGNETIQUE ................................................................................................................. 46
INTRODUCTION ................................................................................................................ 47
4.1 Généralités .......................................................................................................................... 47
4.2 Le corps noir ........................................................................................................................ 50
4.2.1 Définition ...................................................................................................................... 50

4.2.2 La loi de Stefan-Boltzmann ....................................................................................... 51

4.2.3 La loi de Planck ........................................................................................................... 52

4.3 Loi de Wien .......................................................................................................................... 52

4.3.1 Enoncé ......................................................................................................................... 52

4.3.2 Energie émise par le corps noir entre deux longueurs d’onde ............................ 53

4.4 Corps réels .......................................................................................................................... 55

4.5 Loi de Kirchhoff : relation entre absorption et émission ............................................... 56
ii
4.6 Flux radiatif échangé entre deux surfaces noires .......................................................... 57
CONCLUSION .................................................................................................................... 63
CONCLUSION GENERALE ............................................................................................. 64
REFERENCES BIBLIOGRAPHIQUES .......................................................................... 65

iii
NOMENCLATURE

Φ : Flux de chaleur [W]

𝜑 : Densité de flux de chaleur thermique [W.m-2]

T : Température [K]

ℎ : Coefficient d’échange entre la paroi et le fluide ou [W.m-2.K-1]

coefficient de transfert de chaleur

Cp : chaleur massique ou spécifique à pression constante [J.kg-1.K-1]

Cv : chaleur massique ou spécifique à volume constant [J.kg-1.K-1]

𝑡 : temps [s]

λ : conductivité thermique [W.m-1.K-1]

𝑆 : Surface [m2]

𝑒 : Epaisseur du matériau [m]

R th : Résistance thermique [K.W -1]

r Résistance thermique surfacique [K.W -1.m2]

M : masse molaire [kg.mol-1]

μ : viscosité dynamique [kg.m-1.s-1]

𝜈 : viscosité cinématique [m2.s-1]

𝑞 : Densité volumique [W.m-3]

ρ : masse volumique [kg.m-3]

D : diffusivité de masse [m2.s-1]

𝐶 : Concentration du soluté [kg.m-3]

φm : Densité de flux massique [kg.m-2.s-1]

𝛼 : Diffusivité thermique [m2.s-1]

iv
 U : vitesse caractéristique de l’écoulement [m.s-1]

δT : Couche limite thermique

δT : Couche limite dynamique

LV : chaleur latente de vaporisation [J.kg-1]

τν : Temps de diffusion de la quantité de mouvement

τα : Temps de diffusion thermique

Dh : Diamètre hydraulique [m]

V : Volume [m3]

R : constante des gaz parfaits [J.mol-1.K-1]

r : Rayon [m]

L : Longueur [m]

x : Abscisse [m]

y : Ordonnée [m]

z : Ordonnée [m]

H : Hauteur [m]

Re : Nombre de Reynolds

Gr : Nombre de Grashof

Pr : Nombre de Prandtl

Pe Nombre de Péclet

Ra Nombre de Rayleigh

Nu Nombre de Nusselt

m : Masse [kg]

Lref : Longueur caractéristique [m]

v
𝑀0 : Emittance du corps noir [W.m-2]

𝑀 : Emittance du corps réel [W.m-2]

𝜎 Constante de Stefan Boltzmann [W.m-2.K-4]

𝜀 Emissivité monochromatique

ϕi,n : Flux net [W]

Indices

L : Liquide

G : Gaz

V : Vapeur

S : Solide

H : Horizontal

∞ : Ambiant

V : Vertical

vi
SIGLES ET ABBREVIATIONS

% : Pourcentage

CO2 : Dioxyde de carbone

H2O : Eau

U.V : Ultra-violet

I.R : Infra-rouge

vii
LISTE DES TABLEAUX
Tableau 2-1 : Valeurs physiques de la conductivité thermique de quelques corps métalliques
............................................................................................................................................................... 17
Tableau 2-2 : Valeurs physiques de la conductivité thermique de quelques matériaux non-
métalliques........................................................................................................................................... 17
Tableau 2-3 : Valeurs comparées de la conductivité et de la diffusivité thermique de
quelques corps.................................................................................................................................... 23
Tableau 3-1 : Ecoulement dans une conduite de diamètre hydraulique 𝐷ℎ ................................ 43
Tableau 3-2 : Ecoulement autour d’un cylindre de diamètre 𝐷 ...................................................... 44
Tableau 3-3 :Ecoulement autour d’une sphère de diamètre 𝐷 ................................................... 44

viii
LISTE DES FIGURES
Figure 1-1 : Chauffage d’eau .............................................................................................................. 7
Figure 1-2 : Fenêtre à double vitrage ................................................................................................ 7
Figure 2-1 : Interaction entre les molécules ................................................................................... 10
Figure 2-2 : Notations ........................................................................................................................ 12
Figure 2-3 : Exemple .......................................................................................................................... 13
Figure 2-4 : Exemples de courbes isothermes ............................................................................. 13
Figure 2-5 : gradient de température............................................................................................... 15
Figure 2-6: définition d’un tube de flux ............................................................................................ 16
Figure 2-7 : Volume V ........................................................................................................................ 19
Figure 2-8 : Réalisation expérimentale d’une condition aux limites isothermes à T0............... 20
Figure 2-9 : Réalisation expérimentale d’une condition de flux constant : Toutes les faces du
bloc de cuivre autres que la frontière étudiée sont isolées thermiquement .............................. 21
Figure 2-10 : Brique en régime stationnaire (a) et instationnaire (b) .......................................... 23
Figure 2-11: Mur plan......................................................................................................................... 25
Figure 2-12 : Mur plan .......................................................................................................................... 25
Figure 2-13 : Mise en série de résistances thermiques ................................................................ 26
Figure 2-14 : Mur avec sources internes ........................................................................................ 26
Figure 2-15 : Mur avec conditions de convection aux extrémités ............................................... 28
Figure 2-16 : Résistances thermiques équivalente dans le cas d’un mur avec convection ... 28
Figure 2-17 : Conduction dans une coque cylindrique ................................................................. 30
Figure 2-18 : Conduction dans une coque sphérique ................................................................... 30
Figure 3-1 : Nature du régime de l’écoulement ............................................................................. 35
Figure 3-2 : Couche limite ................................................................................................................. 37
Figure 4-1 : Spectre d’un rayonnement électromagnétique......................................................... 48
Figure 4-2 : Divers rayonnements électromagnétiques classées selon leur longueur d’onde 49
Figure 4-3 : Définition des flux réfléchis, absorbés, transmis ...................................................... 49
Figure 4-4: Réalisation expérimentale d’un corps noir ................................................................. 51
Figure 4-5 : Emittance monochromatique du corps noir .............................................................. 54
Figure 4-6 : Courbe donnant la fraction d’énergie émise entre 0 et λT...................................... 55
Figure 4-7 : Courbe donnant la fraction d’énergie émise entre 0 et λT...................................... 56
Figure 4-8 : Echange radiatif entre S1 et S2 ................................................................................... 58
Figure 4-9 : Cavité composée de N surfaces noires ..................................................................... 59
Figure 4-10 : Facteur de forme entre surfaces rectangulaires parallèles .................................. 60
Figure 4-11 : Facteur de forme entre disques parallèles ............................................................. 61
Figure 4-12 : Facteur de forme entre surfaces rectangulaires perpendiculaires ...................... 61
Figure 4-13 : Facteur de forme entre cylindres concentriques .................................................... 62
Figure 4-14 : Facteur de forme entre disques coaxiaux de diamètres différents ..................... 63

ix
INTRODUCTION GENERALE
La thermodynamique permet de prévoir la quantité totale d’énergie qu’un système doit
échanger avec l’extérieur pour passer d’un état d’équilibre à un autre. La thermique
(ou thermocinétique) 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. Les multiples
procédés utilisés dans l'industrie sont très souvent le siège d'échanges de chaleur, soit
parce que c'est le but recherché (fours, coulée, échangeurs, thermoformage, induction,
lits fluidisés, trempe, refroidissement), soit parce que ceux-ci interviennent d'une
manière inévitable (chocs thermiques, pertes de chaleurs, rayonnement).
Cet ouvrage, de parcours Licence, constitue à la fois une introduction aux transferts
thermiques fondée sur la comprehension physique des phénomènes et un document
de base pour l’ingenierie thermique. Les trois modes de transfert (par conduction,
rayonnement et convection) thermiques sont abordés et quelques fois couplés dans
des applications concrètes de façon non approfondie.
La maitrise du premier principe de la thermodynamique et de notions élémentaires de
mathématiques de niveau licence de mécanique ou physique est suffisante pour une
complète assimilation du cours.
Ce cours est focalisé sur la compréhension physique des phénomènes et les transferts
thermiques sont le plus souvent traités à une dimension, de façon à alléger le
formalisme mathématique. Son exploitation permet de dimensionner des systèmes
simples au sein d’applications variées.
Dans le premier chapitre, nous abordons l’introduction aux transferts thermiques en
présentant brièvement les trois modes de transferts thermiques. Le deuxième chapitre
présente la conduction thermique ainsi que quelques exemples de résolution de
problèmes unidimensionnels. Quelques analogies de transferts par rapport à certaines
disciplines ont été introduites. Dans le troisième chapitre est consacré à la convection
thermique et repose essentiellment sur les principaux nombres adimensionnels. Enfin
le quatrième chapitre présente les transferts thermiques par rayonnement
électromagnétique.

1
1. CHAPITRE 1 : INTRODUCTION AUX TRANSFERTS
THERMIQUES

2
INTRODUCTION
Grâce à l'étude de la thermodynamique, vous avons appris que l'énergie peut être
transférée par interactions d'un système avec son environnement. Ces interactions
sont appelées travail et chaleur. Cependant, la thermodynamique traite des états
finaux du processus au cours desquels une interaction se produit et ne fournit aucune
information sur la nature de l’interaction ou sur leur rythme auquel cela se produit.
L'objectif de ce chapitre est d'étendre l'analyse thermodynamique par l'étude des
modes de transfert de chaleur et de poser les bases fondamentales du transfert. Nous
sommes amenés à répondre à ces différentes questions : Qu'est-ce que le transfert
de chaleur ? Comment la chaleur est-elle transférée ? Pourquoi est-ce important ?

1.1 Généralités
De tout temps, les problèmes de transmission d'énergie, et en particulier de la chaleur,
ont eu une importance déterminante pour l'étude et le fonctionnement d'appareils tels
que les générateurs de vapeur, les fours, les échangeurs, les évaporateurs, les
condenseurs, etc, mais aussi pour des opérations de transformations chimiques. En
effet, dans certains systèmes réactionnels, c'est la vitesse des échanges de chaleur
et non la vitesse des réactions chimiques qui détermine le coût de l'opération (cas de
réactions fortement endothermique ou exothermique). D’une façon générale, nous
rencontrons ces situations dans le domaine industriel (moteurs thermiques ou même
électriques, centrales électriques au fuel au gaz, etc.…) ou domestique (chauffage de
l’habitat).
En outre, de nos jours, par suite de l'accroissement relatif du prix de revient de
l'énergie, on recherche dans tous les cas à obtenir le rendement maximal d'une
installation pour une dépense d'énergie minimale. Les problèmes de transfert de
chaleur sont nombreux, et on peut essayer de les différencier par les buts poursuivis
dont les principaux sont :
 L’augmentation de l'énergie transmise ou absorbée par une surface,
 L’obtention du meilleur rendement d'une source de chaleur,
 La réduction ou l'augmentation du passage d'un débit de chaleur d'un milieu à
un autre.
Le potentiel qui provoque le transport et le transfert de l'énergie thermique est la
température. Lorsque deux systèmes sont à des températures différentes, le système
3
le plus chaud cède de la chaleur au plus froid. Il y a échange thermique ou encore
transfert thermique entre ces deux systèmes. Si deux points matériels placés dans un
milieu thermiquement isolé sont à la même température, on peut affirmer qu'il n'existe
aucun échange thermique global entre ces deux points dits en équilibre thermique (il
s'agit bien d'un équilibre thermique car chacun des points matériels émet une énergie
thermique nette de même module, mais de signe opposé).

1.2 Les différents modes de transferts

Un transfert d’´énergie donne lieu à un flux de chaleur qui correspond à un
déplacement de l’énergie du plus chaud vers le plus froid.
Le flux de chaleur dont la densité locale est notée 𝜑
⃗ est une grandeur vectorielle, ce
qui signifie qu’un flux de chaleur est caractérisé non seulement par son intensité mais
aussi par sa direction. Il est défini en chaque point de l’espace et a l’unité d’une densité
surfacique de puissance (W/m2).
Il existe trois modes essentiels de transferts de chaleur : la conduction, le rayonnement
et la convection.

1.2.1 La conduction
La température est une fonction croissante de l’agitation moléculaire dans un corps,
qu’il soit solide, liquide ou gazeux. Considérons pour l’instant un corps solide au sein
duquel la température varie. L’agitation moléculaire (agitation thermique) élevée de la
zone chaude communiquera de l’énergie cinétique aux zones plus froides par un
phénomène appelé conduction de la chaleur. La conduction est un phénomène de
diffusion qui permet donc à la chaleur de se propager à l’intérieur d’un corps solide.
Dans ce cas, la chaleur (transport d'énergie) se propage sans transport de matière. Il
en est de même pour un liquide ou un gaz.
Dans la plupart des cas, nous étudions la conduction dans les milieux solides, puisque
dans les milieux fluides (liquide ou gazeux), il y a souvent couplage avec un
déplacement de matière. Ainsi nous verrons par la suite que pour les fluides, la
convection est un autre mode de transfert de chaleur possible.
La conduction est le seul mécanisme intervenant dans le transfert de chaleur dans un
solide homogène, opaque et compact. Notons enfin que la conduction de la chaleur
n’est pas possible dans le vide puisqu’il n’y a pas de support moléculaire pour cela.

4
 1.2.2 Le rayonnement
La chaleur du soleil frappe pourtant notre planète alors qu’il n’y a aucun support solide,
liquide ou gazeux au-delà de l’atmosphère terrestre. Ceci signifie donc que l’énergie
thermique peut tout de même traverser le vide. Ce mode de transfert s’appelle le
rayonnement. Il correspond à un flux d’ondes électromagnétiques émises par tout
corps, quelle que soit sa température. Comme on l’imagine, le rayonnement
électromagnétique est d’autant plus élevé que sa température est grande. Comme
pour la conduction, ce sont les interactions entre atomes et molécules qui sont à
l’origine de ce rayonnement. Elles peuvent le générer, ce qui diminue leur énergie, ou
encore l’absorber, ce qui l’augmente.
Le rayonnement thermique est le mode de transmission par lequel la chaleur passe
d’un corps à haute température à un autre plus froid sans nécessité de support
matériel. C’est donc le seul mode de transfert de chaleur qui peut se propager dans le
vide. De par sa nature, le rayonnement n’intervient que dans les milieux transparents
(gaz, verre, vide) ou semi-opaque (gaz+fumées de CO2, gaz+vapeur d’eau).
Le rayonnement thermique ne diffère des autres ondes électomagnétiques, comme
les ondes hertziennes par exemple, que par son origine : la température. En effet tout
corps rayonne tant que sa température est différente de 0 K. Il est caractérisé par des
longueurs d'ondes comprises entre 0,1 et 100 µm, il inclut le domaine du visible (ondes
lumineuses ou lumière de 0,4 et 0,8 µm et n'occupe qu'une faible portion du spectre
d'ondes électromagnétiques.

1.2.3 La convection
Un débit ou une circulation de liquide ou de gaz peut transporter avec lui une certaine
quantité d’énergie thermique. Ce transport de chaleur porte le nom de convection
thermique. Ce transport de l’énergie par un écoulement est analogue au transport
d’autres quantités scalaires (non vectorielles): transport d’une concentration de sel par
de l’eau, transport de l’humidité par l’air, etc…. On retiendra donc que dans la
convection, la chaleur se sert du fluide comme véhicule pour se déplacer.
Il existe deux types de transferts convectifs:
 La convection naturelle : les mouvements sont dus aux variations de masse
volumique dans un fluide soumis au champ de pesanteur. Les variations de
masse volumique peuvent être générées par des gradients de température (l’air
5
 chaud est plus léger que l’air froid) et/ou par des gradients de composition.
Cette différence de température créée entre deux points d’un fluide, fait que le
fluide chaud, qui aura une masse volumique plus faible que le fluide froid aura
tendance à monter sous l’effet de la poussée d’Archimède.
 La convection forcée : le mouvement du fluide est forcé par un dispositif
mécanique quelconque (pompe pour un liquide, ventilateur pour de l’air).
 On parlera de convection mixte lorsque les deux types de convection coexistent
dans un système.
Si l’on prend l’exemple d’un chauffage domestique, l’eau chaude qui arrive dans les
radiateurs circule par convection forcée, entretenue par le circulateur (petite pompe
située dans la chaufferie) tandis que l’air des pièces de la maison circule par
convection naturelle depuis le radiateur autour duquel il s’échauffe jusqu’au plafond
vers lequel il s’élève avant de redescendre pour former un circuit fermé.
En convection, on caractérise le flux de chaleur 𝜙 qui est extrait par le fluide de
température 𝑇𝑂 d’une paroi de surface S à la température 𝑇𝑃 par :
𝜙 = ℎ𝑆(𝑇𝑃 − 𝑇𝑂 )
où 𝜙 est en Watt, 𝑆 en m2, 𝑇 en Kelvin et où ℎ désigne le coefficient d’échange entre
la paroi et le fluide (en W.m−2.K−1).

1.3 Combinaison des différents modes de transferts

Dans beaucoup de situations, il y a coexistence de deux (2) ou même trois (3) des
modes de transferts thermiques. Fort heureusement, il est fréquent qu’un mode soit
prépondérant et simplifie l’analyse. Signalons aussi que certains échanges de chaleur
s’accompagnent d’un changement d’état (vaporisation, condensation, fusion, etc…).
Ces phénomènes se comportent alors comme une source (exemple de la
condensation) ou un puits de chaleur (exemple de la vaporisation). Les figures
suivantes illustrent les différents de transferts thermiques.

6
Figure 1-1 : Chauffage d’eau

Figure 1-2 : Fenêtre à double vitrage

7
CONCLUSION
Les trois modes de transfert thermique sont la conduction, la convection et le
rayonnement. La conduction nécessite un milieu matériel et l'énergie est transportée
de proche en proche, généralement dans un solide, sans déplacement de matière
tandis que la convection nécessite un milieu matériel fluide et l'énergie est transportée
par des mouvements de la matière. Le rayonnement ne nécessite pas forcément de
milieu matériel. L'énergie est transportée par des ondes électromagnétiques.

8
2. CHAPITRE 2 : TRANSFERT THERMIQUE PAR
CONDUCTION

9
INTRODUCTION
La conduction est le mode de transfert de chaleur sans déplacement apparent de la
matière sous l'effet d'un gradient de température. C'est le seul mode de transfert
intervenant au sein des milieux solides.
La conduction peut être considérée comme le transfert d'énergie des particules les
plus énergétiques vers les particules les moins énergétiques d'une substance en
raison des interactions entre les particules. Le transfert de chaleur par conduction à
l'intérieur d'un corps s'effectue selon deux mécanismes distincts : une transmission
par la vibration des atomes et une transmission par les électrons libres.

2.1 Généralités
Le transfert par conduction s’effectue par contact. En chauffant l'extrémité d'une tige
métallique, la chaleur se propage dans la tige. On dit qu'il y a conduction lorsque la
chaleur (transport d'énergie) se propage sans transport de matière. On sait que les
molécules/atomes sont en perpétuelle agitation thermique (oscillations autour des
positions d'équilibre (solides, liquides), déplacements désordonnés (gaz)). A cette
agitation thermique est associée une énergie (d'agitation thermique) proportionnelle à
la température T. Les molécules sont en perpétuelle interaction les unes avec les
autres ; au cours de ces chocs, elles échangent de l'énergie ; une molécule "excitée"
peut ainsi perdre un peu de son énergie au profit de ses voisines avec lesquelles elle
interagit (figure suivante).

Figure 2-1 : Interaction entre les molécules

10
Ainsi la molécule 1 va choquer la molécule 2 et globalement lui céder une partie de
son énergie. La molécule 2 va choquer la molécule 3, etc.... Une partie de l'énergie de
la molécule 1 va donc être transférée vers la droite, vers les molécules moins excitées
(donc de température inférieure) et ceci sans déplacement de cette molécule 1. D’où
un transfert de chaleur, dans la matière, sans transfert de matière.
Les molécules effectuent un très grand nombre de chocs, les transferts ci-dessus sont
donc des bilans sur l'ensemble des chocs. Des molécules de même excitation (donc
de même température) échangent de l'énergie lors des chocs, mais le bilan global est
nul (transferts équivalents de chaque côté).
Seule la conduction assure un bon transfert de chaleur à travers les solides. Par
exemple, lorsqu'on chauffe un barreau métallique à l'une de ses extrémités, l'autre
extrémité s'échauffe progressivement. Si l'on chauffe suffisamment longtemps, l'objet
métallique aura la même température en tout point. La chaleur s'est propagée à partir
de l'extrémité chauffée dans tout le reste du matériau. Le barreau métallique a conduit
de la chaleur : cette propriété s'appcelle la conduction thermique. Si l'on arrête
subitement de chauffer l'extrémité du barreau métallique, la température diminuera
progressivement puis le barreau retrouvera sa température initiale en l'occurrence
celle de l'air ambiant. La chaleur transmise à travers les murs ou le plancher d'une
maison se fait par conduction thermique. Les bons conducteurs de chaleur sont
souvent de bons conducteurs électriques. Dans les métaux, la conduction fait
intervenir les électrons libres qui les rendent bons conducteurs de la chaleur. En
revanche dans les isolants, la conduction se fait mal. En résumé, il y a une forte
correspondance entre les propriétés thermiques et électriques des solides. La
conduction s'observe aussi dans des fluides au repos mais elle est beaucoup plus
faible que dans un métal. De plus, elle est souvent dominée par la convection.

2.2 La loi de Fourier

2.2.1 Définition de quelques grandeurs physiques
 Température 𝑇 : elle se définit en chaque point d’un corps liquide, solide ou
gazeux. C’est une fonction scalaire de l’espace et du temps lorsque le problème
en dépend (problème instationnaire). L’unité de température est le Kelvin [K] ou
encore le degré Celsius [oC].

11
  Flux de chaleur 𝜙 : c’est la quantité de chaleur qui traverse une surface S par
unité de temps
𝑑𝑄
𝜙=
𝑑𝑡
𝜙 s’exprime en Watt
 ⃗ : elle représente la puissance qui traverse l’unité de surface.
Densité de flux 𝜑
Pour une surface perpendiculaire au flux de chaleur :
𝑑𝜙
𝜑=
𝑑𝑆
Si le flux est homogène en tout point de la surface alors :
𝜙
𝜑=
𝑆
𝜑 s’exprime en W.m−2.
Pour une surface dont la normale 𝑛⃗ est orientée de manière quelconque par
rapport au flux (confère figure 2) alors :
𝑑𝜙 = 𝜑
⃗ . 𝑛.
⃗⃗⃗ 𝑑𝑆
𝑑𝜙 = |𝜑
⃗ |. 𝑐𝑜𝑠 𝛼. 𝑑𝑆
Le flux à travers une surface quelconque s’écrira donc :

𝜙 = ∫𝜑
⃗ . 𝑛.
⃗⃗⃗ 𝑑𝑆
𝑆

Figure 2-2 : Notations

L’exemple de la figure 2-3 illustre l’importance de l’orientation de la surface par

rapport à la densité de flux : soit un solide carré dont les deux faces (ABCD et
EFGH) sont maintenues à une température constante (respectivement 𝑇 + et 𝑇 − )
et dont les 4 autres sont isolées (par de la laine de verre par exemple). Un flux de
chaleur Φ0 circulera de la face ABCD vers la face EFGH. En tout point du cube,
la densité de flux φ
⃗⃗ est donc parallèle à Ox et vaut :

12
 ϕ0
φ=
S

Le flux à travers les 4 faces isolées est nul. Comme on le verra par la suite, s’il
n’y a pas de sources de chaleur internes au cube, le flux de chaleur se conserve.
Ainsi, si l’on calcule le flux ϕ1 qui traverse la surface diagonale DCFE inclinée à
45o, on trouvera ϕ1 = ϕ0 .

Figure 2-3 : Exemple

 Surface isotherme : si dans un milieu on relève les températures en tout point

à un instant donné et qu’on relie entre eux les points de même température on
obtient des surfaces isothermes en 3D et des lignes isothermes en 2D (voir
figure 2-4).

Figure 2-4 : Exemples de courbes isothermes

13
 2.2.2 Enoncé de la loi de Fourier
Considérons à nouveau le montage de la figure 2-3. Imposons une différence de
température 𝑇 + − 𝑇 − = ∆𝑇0 entre les deux faces non isolées : un flux de chaleur 𝜙0
circulera. Doublons cette différence de température : un flux égal à 2𝜙0 circulera alors.
Pour une valeur donnée de l’écart de température, remplaçons le cube par un matériau
différent. La valeur du flux en sera affectée. Ceci donne l’intuition que le flux qui circule
par conduction est proportionnel à la différence de température et à l’aptitude du
matériau à conduire la chaleur.
En 1811, Fourier propose une formulation locale de cette loi, donc valable en tout
point :
𝜑 ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑇
⃗ = −𝜆𝑔𝑟𝑎𝑑
⃗ en W.m−2 à la conductivité thermique λ du matériau
qui relie la densité de flux 𝜑
(W.m−1.K−1) et au gradient local de température. Le signe − de la loi de Fourier résulte
d’une convention qui rend positif un flux de chaleur s’écoulant du chaud vers le froid,
donc dans le sens d’un gradient négatif.
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑉
Cette loi est l’analogue thermique de la loi d’Ohm : 𝑗 = 𝜎𝐸⃗ = −𝜎𝑔𝑟𝑎𝑑

où 𝑗 est la densité de courant électrique, 𝜎 la conductivité électrique du métal et 𝑉 le

potentiel électrique.
Le flux thermique est l’analogue du courant électrique et la température est l’analogue
1𝑙
du potentiel électrique. Comme on définit la résistance électrique 𝑅𝑒 = 𝜎 𝑆 d’un

conducteur de longueur 𝑙 et de section 𝑆, on définit la résistance thermique par :

1𝑙
𝑅𝑡ℎ =
𝜆𝑆
Telle que :
∆𝑇 = 𝑅𝑡ℎ 𝜙
Avec 𝑅𝑡ℎ en K.W −1

2.2.3 Orthogonalité du gradient et de l’isotherme

Considérons le dessin de la figure 2-5. Entre les points M et M′ distants de ⃗⃗⃗⃗⃗⃗
𝑑𝑀 , la
différence de température vaut :

𝑑𝑇 = ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑔𝑟𝑎𝑑 𝑇. ⃗⃗⃗⃗⃗⃗
𝑑𝑀
Avec :
14
 𝜕𝑇 𝜕𝑇 𝜕𝑇
𝑑𝑇 = 𝑑𝑥 + 𝑑𝑦 + 𝑑𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧
Si l’on place le point M′ sur la même isotherme que M, on obtient :

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑔𝑟𝑎𝑑 𝑇. ⃗⃗⃗⃗⃗⃗
𝑑𝑀 = 0

Ce qui signifie que le gradient de température en chaque point est perpendiculaire à

l’isotherme passant par ce point. On en déduit aussi que les lignes de flux 𝜑
⃗ sont elles
aussi perpendiculaires aux isothermes. Ceci permet de définir la notion de tube de flux.
Sur la figure 2-6 les courbes fermées 𝐶1 , 𝐶2 , 𝐶3 appartiennent aux surfaces isothermes
𝑇1 , 𝑇2 et 𝑇3 et sont liées entre elles par une même ligne de flux. Les lignes de flux étant
par définition des surfaces adiabatiques (aucun flux ne les traverse), les surfaces 𝐶1 ,
𝐶2 , 𝐶3 sont donc traversées par le même flux. Le tube ainsi formé s’appelle un tube de
flux. Il laisse passer un flux 𝜙. On peut définir la résistance thermique 𝑅𝑡ℎ de ce tube
entre les isothermes 𝑇1 et 𝑇2 par exemple, telle que :

𝑇1 − 𝑇2 = 𝑅𝑡ℎ 𝜙

est l’analogue de la loi d’Ohm :

∆𝑉 = 𝑅𝑒 𝐼

Figure 2-5 : gradient de température

15
Figure 2-6: définition d’un tube de flux

2.3 La conductivité thermique

2.3.1 Généralités
Comme on le voit sur les tableaux 2-1 et 2-2, la conductivité thermique 𝜆 des solides
varie de 200 W.m−1K−1 pour le mercure à 0,03 W.m−1K−1 pour la laine de verre. En règle
générale, les métaux qui sont bons conducteurs de l´électricité sont bons conducteurs
thermiques.
Les liquides conduisent en général moins bien la chaleur que les solides, sauf lorsque
ce sont des métaux liquides (mercure, sodium liquide, etc…).
La conductivité thermique des gaz dépend de la pression mais aussi fortement de la
température et de la masse molaire 𝑀 du gaz. La relation suivante donne un bon ordre
de grandeur :

𝑇
𝜆=𝐴√
𝑀

Où 𝐴 est une constante.

16
Tableau 2-1 : Valeurs physiques de la conductivité thermique de quelques corps métalliques

Tableau 2-2 : Valeurs physiques de la conductivité thermique de quelques matériaux non-

métalliques

17
 2.3.2 Matériaux anisotropes
Certains matériaux ont une structure qui rend la conductivité thermique différente selon
la direction de propagation de la chaleur. C’est le cas des matériaux fibreux par
exemple (fibre de verre, fibre de carbone). Dans ces matériaux, qu’on appelle
anisotropes, le flux de chaleur aura donc une direction privilégiée.

2.4 L’équation de la chaleur en conduction

2.4.1 Démonstration
Rappelons d’abord la définition de la chaleur massique. Considérons un solide de
masse m. S’il reçoit une quantité de chaleur 𝑑𝑄 sa température s’élève de 𝑑𝑇 telle
que :
𝑑𝑄 = 𝑚𝐶𝑑𝑇
avec 𝑑𝑄 en Joule, 𝑚 en kg, 𝑑𝑇 en K et où 𝐶 désigne la chaleur massique du solide en
J.kg−1.K−1. Cette grandeur est constante pour un liquide ou un solide tandis que pour
un gaz elle varie en fonction de la pression. On est alors amené à définir 𝐶𝑝 et 𝐶𝑉 selon
que l’apport de chaleur se fait à pression ou à volume constant.
On s’intéresse ici à un solide dont on isole un volume V (cf. figure 2-7), dont la normale
𝑛⃗ est orientée vers l’extérieur. Faisons un bilan de l’énergie échangée par ce volume
V par unité de temps. Puisque nous ne considérons que la conduction comme mode
de transfert, ce bilan traduira le fait que le flux conductif qui est entré à travers la
surface S n’a servi qu’à chauffer ou refroidir la masse de solide contenue dans le
volume V. En prévision de nombreuses situations que l’on verra ultérieurement,
envisageons également la possibilité que ce volume V contient des sources internes
de chaleur, q en W.m−3 comme par exemple une résistance électrique libérant de la
chaleur par effet Joule ou encore une source radioactive pour un problème de
conduction thermique dans un cœur de réacteur. On écrit donc :
𝜕𝑇
− ∫𝜑
⃗ . 𝑛. ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑇. 𝑑𝑆 + ∫ 𝑞𝑑𝑉 = ∫ 𝜌𝐶
⃗⃗⃗ 𝑑𝑆 + ∫ 𝑞𝑑𝑉 = ∫𝜆𝑔𝑟𝑎𝑑 𝑑𝑉
𝑆 𝑉 𝑆 𝑉 𝑉 𝜕𝑡
Le signe − dans le calcul du flux tient au fait que l’on fait un bilan dans le sens opposé
à la normale 𝑛⃗. Par ailleurs, la dérivée partielle de la température est due au fait que
T est une fonction de plusieurs variables :
T (x, y, z, t). L’utilisation du théorème d’Ostrogradski permet d’écrire :

18
 𝜕𝑇
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑇] + 𝑞)𝑑𝑉 = ∫ 𝜌𝐶
∫ (𝑑𝑖𝑣[𝜆𝑔𝑟𝑎𝑑 𝑑𝑉
𝑉 𝑉 𝜕𝑡
Cette relation étant valable quel que soit le volume V, on obtient :
𝜕𝑇
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑇) + 𝑞 = 𝜌𝐶
𝑑𝑖𝑣(𝜆𝑔𝑟𝑎𝑑
𝜕𝑡

Figure 2-7 : Volume V

Plusieurs cas particuliers allègent cette équation :
Régime permanent : la température ne varie pas en fonction du temps. De
𝜕𝑇
ce fait = 0 et l’équation de la chaleur devient :
𝜕𝑡

⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑇) + 𝑞 = 0
𝑑𝑖𝑣(𝜆𝑔𝑟𝑎𝑑
Absence de sources de chaleur : on a ici « équilibre » entre l’apport de
chaleur par conduction dans le solide et son échauffement :
𝜕𝑇
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝑇) = 𝜌𝐶
𝑑𝑖𝑣(𝜆𝑔𝑟𝑎𝑑
𝜕𝑡
Conductivité 𝜆 constante : en dehors des solides anisotropes, 𝜆 est une
constante. On se limitera à ce cas dans ce cours. On peut donc le sortir de
l’opérateur div et la loi de Fourier devient :
𝜕𝑇
𝜆∇2 𝑇 + 𝑞 = 𝜌𝐶
𝜕𝑡

où ∇2 = Δ désigne l’opérateur laplacien.

2.4.2 Conditions aux limites

𝜕𝑇
Même dans un cas très simple sans sources et régime permanent ( 𝜕𝑡 = 𝑞 = 0) la

résolution de l’équation de T, qui devient ∇2 𝑇 = 0 ; nécessite l’adjonction de conditions

aux limites aux frontières du domaine étudié. Ces conditions peuvent être de trois
types :
19
  Conditions aux limites isothermes
La température de la frontière est imposée. Elle peut dépendre du temps et même
varier le long de la frontière. Le cas le plus simple est 𝑇 = 𝑇𝑜 = 𝑐𝑠𝑡𝑒 sur toute la
frontière. Une telle condition est expérimentalement difficile à réaliser sauf si le solide
dans lequel on cherche la distribution de température est en contact avec un milieu
extrêmement conducteur (cuivre par exemple) et que ce milieu est lui-même parcouru
par un fluide qui le maintien à 𝑇𝑜 (cf. figure 2-8).

Figure 2-8 : Réalisation expérimentale d’une condition aux limites isothermes à T0

 Conditions aux limites à flux constant
Une densité de flux de chaleur ⃗⃗⃗⃗
𝜌0 est imposée sur la frontière. Cette condition, de par
la loi de Fourier, impose donc la valeur du gradient de température à la frontière du
𝜕𝑇 𝜌
corps considéré ( = − 𝜆0 ) où 𝑛⃗ est la normale à la frontière. La réalisation
𝜕𝑛

expérimentale d’un flux imposé peut se faire par exemple au moyen de résistances
électriques (cf. figure 2-9). Un cas particulier de flux imposé est 𝜌0 = 0 cela signifie
que la paroi est isolée thermiquement (on dit aussi adiabatique) ce qui se réalise en la
calfeutrant avec un matériau isolant (laine de verre, polystyrène). Toujours de par la
loi de Fourier, on en déduit à la paroi (indice P) :
𝜕𝑇
[ ] =0
𝜕𝑛 𝑃
Ainsi les isothermes seront perpendiculaires à la frontière considérée.

20
Figure 2-9 : Réalisation expérimentale d’une condition de flux constant : Toutes les faces du
bloc de cuivre autres que la frontière étudiée sont isolées thermiquement
 Echanges convectifs sur une paroi
Lorsqu’une paroi échange de la chaleur par convection avec l’extérieur, la densité de
flux ϕ qu’elle échange est proportionnelle à la différence de température entre la paroi
et le milieu fluide extérieur, multipliée par un coefficient d’échange ℎ en W.m−2.K−1 qui
tient compte des différentes propriétés physiques et cinématiques du fluide léchant la
paroi :
𝜑 = ℎ(𝑇𝑝𝑎𝑟𝑜𝑖 − 𝑇𝑓𝑙𝑢𝑖𝑑𝑒 )

En tout point de la frontière, ce flux 𝜑 étant fourni par conduction au travers du solide
étudié, on en déduit le gradient de T à cette frontière :
𝜕𝑇
−𝜆 ( ) = ℎ(𝑇𝑝𝑎𝑟𝑜𝑖 − 𝑇𝑓𝑙𝑢𝑖𝑑𝑒 )
𝜕𝑛 𝑝𝑎𝑟𝑜𝑖
 Echanges radiatifs sur une paroi
Le rayonnement d’une paroi peut être un mode d’échange à prendre en compte,
surtout si sa température est élevée (supérieure à 100 degrés environ). Comme on le
verra, le flux qu’une paroi à la température 𝑇𝑝 échange par rayonnement avec le milieu
externe à la température 𝑇𝑒𝑥𝑡 vaut : 𝜎(𝑇 4 𝑝 − 𝑇 4 𝑒𝑥𝑡 ) où 𝜎 est la constante de Stefan-
Boltzmann. Dans un tel cas, la condition aux limites à la frontière deviendra donc :
𝜕𝑇
−𝜆 ( ) = 𝜎(𝑇 4 𝑝𝑎𝑟𝑜𝑖 − 𝑇 4 𝑒𝑥𝑡 )
𝜕𝑛 𝑝𝑎𝑟𝑜𝑖
21
  Condition de passage entre 2 solides :
Si le problème étudié comporte deux milieux (ou plus) on devra écrire une condition
aux limites à chaque frontière entre deux solides en contact. La conservation de la
chaleur de part et d’autre de la frontière impose :
𝜕𝑇1 𝜕𝑇2
𝜆1 ( )𝑝𝑎𝑟𝑜𝑖 = 𝜆2 ( )
𝜕𝑛 𝜕𝑛 𝑝𝑎𝑟𝑜𝑖
où 𝜆1 et 𝜆2 sont les conductivités thermiques des solides 1 et 2, et 𝑇1 et 𝑇2 sont les
distributions de température dans ces solides.

2.4.3 Conductivité et diffusivité thermique

Considérons la résolution d’un problème de conduction permanent sans sources (𝑞 =
0) dans un milieu homogène. On résout donc :
∇2 𝑇 = 0
Si les conditions aux limites sont isothermes, la conductivité thermique λ du solide
n’intervient pas dans la solution. Elle n’interviendra que si celles-ci prescrivent le flux
𝜕𝑇
qui, on l’a vu, impose la valeur de − 𝜕𝑛 à la frontière.

Considérons maintenant un problème non-permanent de conduction (toujours avec

𝑞 = 0). On résout donc :
𝜕𝑇
𝑎∇2 𝑇 =
𝜕𝑡
𝜆
où 𝑎 = 𝜌𝐶 (en m2.s-1) désigne la diffusivité thermique du solide. Le tableau 2-3 donne
𝑝

les valeurs comparées de 𝜆 et a pour quelques corps. Il est intéressant de comparer

le fer et l’air qui ont des 𝜆 totalement dissemblables mais une diffusivité thermique a
presque égale. On retiendra que la conductivité thermique 𝜆 est l’aptitude à transmettre
un flux sous l’effet d’une différence de température tandis que la diffusivité thermique
a est l’aptitude à égaler plus ou moins vite dans le temps les différences de
température présentes dans le corps étudié.
Sur la figure 2-10, c’est la conductivité thermique λ qui fixera la valeur du flux dans la
brique lorsqu’on la soumet à un écart de température ∆T. C’est par contre la diffusivité
thermique a qui fixera le temps nécessaire à ce que la face de droite ressente l’effet
d’un choc thermique imposé à t=0 sur la face gauche (cf. figure 2-10).

22
Tableau 2-3 : Valeurs comparées de la conductivité et de la diffusivité thermique de quelques
corps

Figure 2-10 : Brique en régime stationnaire (a) et instationnaire (b)

2.4.4 Equation de la chaleur en régime permanent

En régime permanent et dans un matériau isotrope, l’équation de la chaleur est une
équation de Poisson :
𝑞
Δ𝑇 = −
𝜆
L’expression développée varie avec le système de coordonnées :
 Cartésiennes :
𝜕 2𝑇 𝜕 2𝑇 𝜕 2𝑇 𝑞
2
+ 2+ 2 =−
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜆

 Cylindriques :

23
 1𝜕 𝜕𝑇 1 𝜕 2𝑇 𝜕 2𝑇 𝑞
(𝑟 ) + 2 2 + 2 = −
𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧 𝜆
 Sphériques :
1 𝜕 2 𝜕𝑇 1 𝜕 𝜕𝑇 1 𝜕 2𝑇 𝑞
(𝑟 ) + (𝑠𝑖𝑛𝜃 ) + = −
𝑟 2 𝜕𝑟 𝜕𝑟 𝑟 2 𝑠𝑖𝑛𝜃 𝜕𝜃 𝜕𝜃 𝑟 2 𝑠𝑖𝑛2 𝜃 𝜕𝜑 2 𝜆
Les problèmes 1D peuvent souvent se résoudre analytiquement. Les problèmes 3D
nécessitent souvent l’emploi de méthodes numériques traitées par ordinateur. En 1D,
les équations ci-dessus deviennent :
En cartésien :
𝑑2𝑇 𝑞
2
=−
𝑑𝑥 𝜆
En cylindrique :
1𝑑 𝑑𝑇 𝑞
(𝑟 ) = −
𝑟 𝑑𝑟 𝑑𝑟 𝜆
En sphérique :
1 𝑑 2 𝑑𝑇 𝑞
(𝑟 ) = −
𝑟 2 𝑑𝑟 𝑑𝑟 𝜆
2.5 Exemples de résolution de l’équation de la chaleur en
conduction
2.5.1 Cas d’un mur plan
Des températures 𝑇1 et 𝑇2 sont imposées aux bornes d’un mur d’épaisseur e et de
conductivité thermique 𝜆 (cf. figure 2-12). L’équation de la chaleur 1D cartésienne :
𝑑2 𝑇
= 0 admet la solution 𝑇 = 𝐴𝑥 + 𝐵.
𝑑𝑥 2

Les conditions aux limites en 𝑥 = 0 et en 𝑥 = 𝑒 permettent d’éliminer A et B pour

donner :
𝑥
𝑇 = 𝑇1 − (𝑇1 − 𝑇2 )
𝑒
La distribution de température est donc linéaire et les isothermes sont régulièrement
distribuées de 𝑥 = 0 à 𝑥 = 𝑒. La densité de flux en tout point vaut :
𝑑𝑇 𝑇1 − 𝑇2
𝜑 = −𝜆 =𝜆
𝑑𝑥 𝑒

24
et est bien positive si 𝑇1 > 𝑇2 (flux allant du chaud vers le froid). Le flux 𝜙 traversant
𝑇1 −𝑇2
une surface S vaut donc 𝜆 S de sorte que la résistance thermique R définie par
𝑒

𝑇1 − 𝑇2 = 𝑅𝜙 vaut :
1𝑒
𝑅=
𝜆𝑆
Notons enfin que le flux de chaleur 𝜙 est le même quelle que soit l’abscisse à laquelle
on le calcule (de 0 jusqu’à e). Ceci est garanti par la pente constante de la distribution
de température.

Figure 2-12 : Mur plan

2.5.2 Association de murs en série ou en parallèle

Lorsque plusieurs murs sont en série, la résistance totale, 𝑅𝑡 du tube de flux ainsi formé
vaut :
1 𝑒𝑖
𝑅𝑡 = ∑ 𝑅𝑖 = ∑
𝜆𝑖 𝑆𝑖
𝑖 𝑖

Leur mise en parallèle amène à l’expression suivante de 𝑅𝑡 , analogue de celle obtenue

en électricité :
1 1 𝜆𝑖 𝑆𝑖
=∑ =∑
𝑅𝑡 𝑅𝑖 𝑒𝑖
𝑖 𝑖

Cette méthode s’applique à des résistances de sections différentes à condition que

l’on raisonne sur un tube de flux. C’est le cas par exemple sur la figure 2-13 sur laquelle
les frontières latérales des 5 résistances sont isolées.

25
Figure 2-13 : Mise en série de résistances thermiques

2.5.3 Mur avec sources internes

C’est le cas d’une paroi ou d’un plancher chauffé intérieurement par des résistances
électriques (cf. figure 2-14).
 Températures d’extrémités identiques
Appelons T0 la température aux extrémités de la paroi et q (W.m−3) la densité
volumique des sources. L’équation de la chaleur s’écrit :

𝑑2𝑇 𝑞
2
=−
𝑑𝑥 𝜆

Figure 2-14 : Mur avec sources internes

Dont la solution est une parabole :

𝑞 𝑥2
𝑇=− + 𝐴𝑥 + 𝐵
𝜆 2

Les conditions aux limites en x = +e et x = −e donnent finalement :

26
 𝑞 2
𝑇 = 𝑇0 + (𝑒 − 𝑥 2 )
2𝜆
𝑞𝑒 2
L’écart de température au centre du mur vaut donc : Δ𝑇 = La densité de flux vaut :
2𝜆

𝜑 = 𝑞𝑥
Elle est nulle au centre par symétrie et maximum en x = ±e. Les isothermes sont donc
plus serrées à cet endroit.
Le flux 𝜑 = 𝑞𝑥 en x = e est tel qu’il évacue l’ensemble des sources situé entre x = 0 et
x = e.
 Températures d’extrémités différentes
Appelons T1 et T2 les températures en x = −e et en x = e. La solution générale est
inchangée et seules les conditions aux limites changent. On en déduit la solution :

𝑞 2 𝑇1 − 𝑇2 𝑇1 + 𝑇2
𝑇= (𝑒 − 𝑥 2 ) − 𝑥+
2𝜆 2𝑒 2

2.5.4 Mur avec conditions de convection aux parois

La température au loin de l’air qui entoure le mur vaut T1 à gauche et T2 à droite (cf.
figure 2-15).
Appelons ℎ le coefficient d’échange aux parois. La solution générale de l’équation de
la chaleur est inchangée : 𝑇 = 𝐴𝑥 + 𝐵
Les conditions aux limites se traduisent par :
𝑥 = 0 ; 𝜑 = ℎ(𝑇1 − 𝐵) = −𝜆𝐴
𝑥 = 𝑒 ; 𝜑 = ℎ(𝐴𝑒 + 𝐵 − 𝑇2 ) = −𝜆𝐴

27
Figure 2-15 : Mur avec conditions de convection aux extrémités
On peut éliminer A et B de ces deux équations et obtenir la solution finale. On peut
aussi utiliser la notion de résistance thermique (cf. figure 2-16). La résistance
convective 𝑅𝑐𝑣 est telle que 𝑇1 − 𝑇𝑃1 = 𝑅𝑐𝑣 𝜙 mais 𝜙 = ℎ𝑆(𝑇1 − 𝑇𝑃1 ) d’où :
1
𝑅𝑐𝑣 =
ℎ𝑆
1𝑒
tandis que 𝑅𝑐𝑑 = 𝜆 𝑆

Figure 2-16 : Résistance thermique équivalente dans le cas d’un mur avec condition de
convection

La résistance totale entre T1 et T2 vaut donc :

2 𝑒
𝑅𝑡 = 2𝑅𝑐𝑣 + 𝑅𝑐𝑑 = +
ℎ𝑆 𝜆𝑆
𝑇1 −𝑇2
d’où le flux 𝜙 = qui circulera à travers le mur. On en déduit les températures de
𝑅𝑡

paroi :
T1 − T2
Tp1 = T1 −
eh
2+
λ

28
 T1 − T2
Tp2 = T2 +
eh
2+
λ
A l’intérieur du mur la solution linéaire s’écrit :
𝑥
𝑇 = 𝑇𝑃1 − (𝑇𝑃1 − 𝑇𝑃2 )
𝑒
2.5.5 Conduction dans une coque cylindrique
Considérons une coque cylindrique limitée par les rayons R1 et R 2 de températures T1
et T2 , suffisamment longue pour que l’on puisse considérer le transfert de chaleur
uniquement radial (cf. figure 2-17). L’équation de la chaleur s’écrit :
1𝑑 𝑑𝑇
(𝑟 ) = 0
𝑟 𝑑𝑟 𝑑𝑟
𝑑𝑇
d’où : 𝑟 𝑑𝑟 = 𝐴 ce qui donne :

𝑇 = 𝐴𝑙𝑛𝑟 + 𝐵

Les conditions aux limites fournissent les constantes A et B et on en déduit la solution

finale
𝑇 − 𝑇1 ln(𝑟⁄𝑅1 )
=
𝑇2 − 𝑇1 ln(𝑅2 ⁄𝑅1 )
On en déduit le flux à travers une surface cylindrique quelconque de rayon r et de
hauteur L située entre R1 et R 2 :
𝜙 = 𝑆𝜑
𝑑𝑇
𝜙 = 2πrLφ = −2πrLλ
𝑑𝑟
𝑇2 − 𝑇1
𝜙 = −2πLλ
ln(𝑅2 ⁄𝑅1 )
ce qui fournit la valeur de la résistance thermique de la coque :

ln(𝑅2 ⁄𝑅1 )
𝑅𝑐𝑦𝑙 =
2𝜋𝜆𝐿

On pourra vérifier que si R2 et R1 sont très proches, le problème se comporte comme

𝑒
un mur (disparition de l’effet de courbure) et 𝑅𝑐𝑦𝑙 tend vers 𝜆𝑆 .

29
Figure 2-17 : Conduction dans une coque cylindrique

2.5.6 Conduction dans une coque sphérique

L’équation de T en coordonnées sphériques s’écrit :

1 𝑑 2 𝑑𝑇
(𝑟 )=0
𝑟 2 𝑑𝑟 𝑑𝑟

Figure 2-18 : Conduction dans une coque sphérique

𝑑𝑇 𝐴
d’où : 𝑟 2 𝑑𝑟 = 𝐴 et 𝑇 = − 𝑟 + 𝐵. Les conditions aux limites fixent A et B de sorte que

(cf. figure 2-18) :

30
 1 1
𝑇 − 𝑇1 𝑟 − 𝑅1
=
𝑇2 − 𝑇1 1 1
𝑅2 − 𝑅1
Le flux de chaleur vaut alors :
𝜙 = 𝑆𝜑
𝑑𝑇
𝜙 = 4π𝑟 2 φ = −4π𝑟 2 λ
𝑑𝑟
𝑇1 − 𝑇2
𝜙 = 4πλ
1 1
−
𝑅1 𝑅2
et est bien positif si 𝑇1 > 𝑇2 . La résistance thermique d’une coque vaut ainsi :
1 1
𝑅1 − 𝑅2
𝑅𝑠𝑝ℎ =
4𝜋𝜆
2.5.7 Remarques finales sur la conduction
La plupart des exemples traités dans le chapitre se rapportent à un corps solide. Il faut
garder à l’esprit le fait que la conduction existe aussi dans les fluides (liquides ou gaz)
mais que, en général, ceux-ci restent rarement immobiles et, de ce fait, transportent
en se déplaçant un flux de chaleur infiniment plus grand que celui qui circule par
conduction seule. Ce mode de transport est la Convection et sera vu ultérieurement.
 Analogie avec la diffusion de masse
Nous venons de voir qu’un flux de chaleur est proportionnel à la conductivité thermique
et au gradient de température. Considérons un liquide (eau par exemple) contenant
une distribution de sel inhomogène. Du fait des différences de concentration C de sel
dans le liquide, une migration des molécules de sel s’effectuera, elle aussi dirigée dans
le sens des gradients négatifs de C, c’est à dire du plus concentré vers le moins
concentré. Un débit massique 𝑞𝑚 de sel (en kg.s−1) s’en suivra dont la loi de Fick fixe
la densité :
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝐶
𝜌𝑚 = −𝐷𝑔𝑟𝑎𝑑
⃗⃗⃗⃗⃗
𝜌𝑚 représente le vecteur de densité de flux massique (en kg.s−1.m−2), D la diffusivité
où ⃗⃗⃗⃗⃗
de masse (m2.s−1) et C la concentration de soluté (sel) dans le solvant (eau), en kg.m−3.
La divergence de la loi de Fick fournit une équation différentielle très semblable à celle
de la chaleur :
𝜕𝐶
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ 𝐶)
= 𝑑𝑖𝑣(𝐷𝑔𝑟𝑎𝑑
𝜕𝑡
31
où le second membre se résume à D∇2C si D est constant. Cette équation est
𝜕𝑇
l’analogue de = 𝑎∇2 𝑇 pour la conduction de la chaleur sans sources.
𝜕𝑡

 Résolution numérique d’un problème de conduction

La résolution d’un problème industriel de conduction nécessite souvent d’avoir recours
à un outil numérique. Il existe un grand nombre de logiciels de conduction de la chaleur
fonctionnant sur PC ou sur stations de travail.

CONCLUSION
Ce chapitre a permis d’étudier les équations régissant les échanges thermiques par
conduction, essentiellement la loi de Fourier, l’équation de la chaleur, et les diverses
équations différentielles ou intégrales de bilan thermique. Ces équations sont
particulièrement rigoureuses et cadrées.
La conductivité thermique ou la capacité calorifique (évaluée à pression constante)
constituent des paramètres pertinents qu’il faut tenir compte dans l’étude de la
conduction thermique.
Il est à souligner que la conduction thermique est essentiellement présent dans les
milieux solides.

32
3. CHAPITRE 3 : TRANSFERT THERMIQUE PAR
CONVECTION

33
INTRODUCTION
Lorsqu'au sein d'un même fluide (liquide ou gaz) se produit une différence de
température, la densité du fluide n'est plus égale partout. Du fait de la poussée
d'Archimède, tout corps plongé dans un fluide reçoit de la part de ce fluide une poussée
verticale dirigée du bas vers le haut, égale au poids du volume de fluide déplacé, cette
différence de densité produit un mouvement au sein du liquide. Ce mouvement de
brassage, dans lequel les parties les plus chaudes du fluide ont tendance à s'élever et
les parties froides et denses à descendre, s'appelle convection.

3.1 Généralités
La convection est un mode de transfert de chaleur qui met en jeu, en plus de la
conduction, le mouvement macroscopique de la matière. Ce phénomène se produit
au sein des milieux fluides (liquides ou gaz) en écoulement ou entre une paroi solide
et un fluide en mouvement. On distingue deux types de convection :
 Convection naturelle : les mouvements du fluide sont dus aux variations de
masse volumique résultant d’une différence de température sur les frontières
(les variations de masse volumique peuvent être générées par des gradients
de température (l’air chaud est plus léger que l’air froid) et/ou à des gradients
de composition (air d’une pièce chauffé par un radiateur, courants océaniques
ou atmosphériques…)) et d’un champ de forces extérieures (la pesanteur).
 Convection forcée : le mouvement du fluide est provoqué par des actions
mécaniques extérieures (pompe, ventilateur…).
 On parlera de convection mixte lorsque les deux types de convection coexistent
dans un système.
Dans le chapitre précédent, nous avons considéré les échanges par convection
seulement comme une condition aux limites pour traiter des problèmes de conduction
dans les solides (le système étudié était le solide qui échangeait de la chaleur par
convection à sa frontière avec le milieu extérieur). Dans ce chapitre, le système étudié
sera le fluide en mouvement, l’état thermique du solide étant alors pris comme
condition aux limites.
L’étude du transfert de chaleur par convection permet de déterminer les échanges de
chaleur se produisant entre une paroi et le fluide en écoulement. On distingue alors

34
classiquement deux grands types de configurations caractérisant la géométrie du
système :

 Ecoulements externes : typiquement les écoulements autour d’obstacles

(aéronautique, échangeurs…).
 Ecoulements internes : concernent les écoulements dans les tuyaux
(échangeurs) ou dans les locaux (thermique du bâtiment).
L’importance du flux de chaleur échangé par convection va dépendre du régime
d’écoulement sous lequel se produisent les échanges : régime laminaire ou turbulent.
Un écoulement laminaire est un écoulement caractérisé par des lignes de courant bien
identifiables parallèles aux parois. Un écoulement turbulent est caractérisé par des
structures tourbillonnaires qui favorisent le brassage du fluide et donc les échanges
de chaleur. Pour certaines configurations, comme par exemple l’écoulement le long
d’une plaque plane, l’écoulement peut évoluer d’un régime laminaire à un régime
turbulent en passant par une phase de transition.

(a) : régime laminaire-écoulement interne (b) : régime turbulent-écoulement interne

Figure 3-1 : Nature du régime de l’écoulement

35
Le transfert de chaleur se produit :
 Par conduction (diffusion) dans la direction y ; par conduction (généralement
négligeable) et convection dans la direction x (figure a).
 Par convection dans toutes les directions dans la zone turbulente sauf dans la
sous-couche laminaire près des parois (figure b).
La notion de température est directement liée à l’agitation des molécules qui
composent la matière (on parle d’agitation thermique). Plus l’agitation (l’énergie
cinétique) est importante, plus la température est élevée. Ainsi le transfert de chaleur
d’une région chaude vers une autre plus froide correspond à un transfert d’énergie
cinétique lors des chocs entre les molécules. De façon similaire, la viscosité
correspond à une dissipation d’énergie liée au transfert de quantité de mouvement lors
de ces mêmes chocs intermoléculaires. On voit donc que les phénomènes de transfert
de chaleur et de quantité de mouvement sont intimement liés, ce qui a amené à
introduire l’analogie de Reynolds : les profils de vitesse et de température au sein d’un
fluide en mouvement dans un tube et soumis à des échanges de chaleur par
convection sont liés par une relation de similitude.

3.2 Coefficient d’échange convectif

Notion de couche limite : c’est une région de l’espace au sein de laquelle sont
observés les gradients de vitesse (couche limite dynamique) ou les gradients de
température (couche limite thermique). Le développement de la couche limite
dynamique est dû au phénomène de diffusion de quantité de mouvement par
frottement visqueux. Le développement de la couche limite thermique est dû au
phénomène de diffusion d’enthalpie. On note δ𝑑 l’épaisseur de la couche limite
dynamique et δ𝑇 l’épaisseur de la couche limite thermique.

36
Figure 3-2 : Couche limite

Au voisinage des parois, compte tenu des faibles vitesses du fluide, le transfert
d’énergie par diffusion est dominant. Cette couche constitue la principale résistance
au transfert de chaleur entre la paroi et le fluide en mouvement. Ainsi au voisinage
immédiat de la paroi, on pourra définir une résistance thermique locale de conduction,
𝛿𝑇
𝑟, telle que : 𝑟 ≈ où λ est la conductivité thermique, 𝑟 la résistance thermique
𝜆

surfacique du fluide. La densité de flux (flux par unité de surface) échangée entre la
paroi et le fluide s’écrit alors :
∆𝑇 ∆𝑇 𝑇𝑃 − 𝑇∞
𝜑= = =
𝑟 𝛿𝑇 ⁄𝜆 𝛿𝑇 ⁄𝜆
Le transfert de chaleur se produit ensuite par convection dans le fluide et la densité
de flux obéit alors à la loi de Newton :
𝜑 = ℎ∆𝑇
où h désigne le coefficient d’échange convectif (W.m-2.K-1)
On aura alors :
∆𝑇 𝜆
𝜑≈ ≈ ℎ∆𝑇 ⇒ ℎ =
𝛿𝑇 ⁄𝜆 𝛿𝑇
Lorsqu’on parle de densité de flux de chaleur, on s’intéresse au flux de chaleur
échangé localement entre la paroi et le fluide. On y associe un coefficient d’échange
convectif local ℎ. Les conditions d’écoulement pouvant varier d’un point à l’autre de la
paroi, le coefficient d’échange et donc le flux de chaleur échangé, peut aussi varier.
On définit ainsi un coefficient d’échange convectif moyen, ℎ qui correspond au

37
coefficient d’échange local moyenné sur toute la surface de la paroi au contact avec
le fluide :
1
ℎ̅ = ∬ℎ 𝑑𝑆
𝑆 𝑆

Les épaisseurs des couches limites au sein desquelles se produisent les transferts de
chaleur et de quantité de mouvement dépendent d’un grand nombre de paramètres
(nature du fluide, régime d’écoulement, taille et géométrie du système, état de surface
de la paroi…) et sont donc difficiles à caractériser. Le coefficient d’échange convectif,
qui donne accès au calcul du flux de chaleur échangé entre la paroi et le fluide, est
directement lié à ces épaisseurs de couches limites et est ainsi une grandeur
extrêmement difficile à évaluer. Ce coefficient est difficile à calculer précisément mais
on peut toutefois donner des ordres de grandeurs (en W.m-2.K-1) :

 Convection forcée : gaz h~100, liquide h~103 à 105. Application : échangeurs,

refroidissement des circuits électroniques…
 Convection naturelle : gaz h~10, liquide h~102. Application : thermique de
l’habitat, météorologie, mouvements dans le manteau terrestre, courants
océaniques…
Remarque : lorsque la turbulence de l’écoulement augmente, l’épaisseur de la sous-
couche laminaire diminue et donc la résistance thermique décroît. Ainsi, le flux de
chaleur échangé pour un écart de température donné, augmente.

3.3 Lois de corrélation pour le coefficient d’échange convectif

Dans le domaine de l’ingénierie, les coefficients d’échange sont calculés à partir de

lois de corrélations, obtenues soit par l’analyse précise des mécanismes qui
gouvernent les transferts dans les couches limites, ou obtenues à partir
d’expérimentations (lois empiriques).

3.3.1 Paramètres caractéristiques de la convection

Les échanges de chaleur par convection se produisent au sein d’un fluide en
écoulement. On a vu que le coefficient d’échange, qui intervient dans le calcul du flux
de chaleur, était lié à de nombreux paramètres, notamment au régime d’écoulement
(laminaire ou turbulent) et à la nature du fluide. On cherche donc à exprimer ce
38
coefficient d’échange en fonction de grandeurs caractéristiques du régime
d’écoulement et de la nature du fluide. Ainsi on pourra écrire que h est fonction de la
vitesse de l’écoulement, des coefficients de diffusion de quantité de mouvement et de
chaleur (qui contrôlent les épaisseurs de couche limite) et d’une longueur
caractéristique du système.
Dans la pratique, on utilise plutôt des grandeurs sans dimensions.
Le régime d’écoulement en convection forcée est caractérisé à partir d’un
nombre sans dimension : le nombre de Reynolds, qui quantifie l’importance
des forces d’inertie (moteur de l’écoulement) par rapport aux forces visqueuses
(dissipation, frein à l’écoulement). Il s’écrit :
𝑈2
𝜌𝐿
𝑓𝑜𝑟𝑐𝑒𝑠 𝑑′𝑖𝑛𝑒𝑟𝑡𝑖𝑒 𝑟𝑒𝑓 𝑈𝐿𝑟𝑒𝑓
𝑅𝑒 = = =
𝑓𝑜𝑟𝑐𝑒𝑠 𝑣𝑖𝑠𝑞𝑢𝑒𝑢𝑠𝑒𝑠 𝜇 𝑈 𝜈
2
𝐿𝑟𝑒𝑓
Où :
𝑈 est la vitesse caractéristique de l’écoulement (m.s-1)
𝐿𝑟𝑒𝑓 est une longueur caractéristique du système étudié (m)
𝜈 est la viscosité cinématique du fluide (diffusivité de quantité de mouvement) définie
par :
𝜇
𝜈=
𝜌
Où :
𝜇 est la viscosité dynamique du fluide (kg.m-1.s-1)
𝜌 est la masse volumique du fluide (kg.m-3)
𝜈 est ainsi une propriété physique du fluide, qui quantifie la capacité de ce fluide à
diffuser la quantité de mouvement (à atténuer les gradients de vitesse). Elle s’exprime
en m2.s-1.

Le régime d’écoulement en convection naturelle est caractérisé à partir d’un

nombre sans dimension : le nombre de Grashof, qui quantifie l’importance des
forces d’Archimède (moteur de l’écoulement) par rapport aux forces visqueuses
(dissipation, frein à l’écoulement). Il s’écrit :
𝑔𝛽∆𝑇𝐿3𝑟𝑒𝑓
𝐺𝑟𝐿𝑟𝑒𝑓 =
𝜈2
Où :
39
𝑔 est l’accélération de la pesanteur (m.s-2).
𝛽 est le coefficient de dilatation thermique (K-1) :
1 𝜕𝜌
𝛽= )
𝜌 𝜕𝑇 𝑃
∆𝑇 est un écart de température caractéristique du système étudié (K).

Le comportement du fluide vis-à-vis des échanges de chaleur par convection

est caractérisé par le nombre de Prandtl. C’est un paramètre sans dimension
défini par le rapport entre la diffusivité de quantité de mouvement et la diffusivité
thermique :
𝜈
𝑃𝑟 =
𝛼
α est la diffusivité thermique du fluide, définie par :
𝜆
𝛼=
𝜌𝐶𝑝
Où :
𝜆 est la conductivité thermique du fluide (W.m-1.K-1)
𝐶𝑝 est la chaleur spécifique du fluide (J.kg-1.K-1)
𝛼 est ainsi une propriété physique du fluide, qui quantifie la capacité de ce fluide à
diffuser la chaleur (à atténuer les gradients de température). Elle s’exprime en m2.s-1.
Le nombre de Prandtl peut être vu comme un rapport de deux temps caractéristiques :
le temps de diffusion de la quantité de mouvement,
𝐿2𝑟𝑒𝑓
𝜏𝜈 =
𝜈
et le temps de diffusion thermique,
𝐿2𝑟𝑒𝑓
𝜏𝛼 =
𝛼
Ainsi un nombre de Prandtl faible (cas des métaux liquides) signifie que la diffusion de
la chaleur dans le fluide se produit très vite (temps de diffusion très court) de telle sorte
que le champ de vitesse n’a pas le temps d’affecter le champ de température.
Inversement, un nombre de Prandtl élevé signifie que le champ de température dans
le fluide est fortement influencé par le champ de vitesse.
Ordres de grandeur (conditions normales de pression et de température) :
𝑃𝑟𝑎𝑖𝑟 ≈ 0,7 ; 𝑃𝑟𝑒𝑎𝑢 ≈ 7 ; 𝑃𝑟ℎ𝑢𝑖𝑙𝑒 ≈ 1000.

40
  Le nombre de Prandtl peut être combiné au nombre de Reynolds pour former
le nombre de Péclet :
𝑈𝐿𝑟𝑒𝑓
𝑃𝑒𝐿𝑟𝑒𝑓 = 𝑅𝑒𝐿𝑟𝑒𝑓 𝑃𝑟 =
𝛼
 Le nombre de Prandtl peut être combiné au nombre de Grashof pour former le
nombre de Rayleigh :
𝑔𝛽∆𝑇𝐿3𝑟𝑒𝑓
𝑅𝑎𝐿𝑟𝑒𝑓 = 𝐺𝑟𝐿𝑟𝑒𝑓 𝑃𝑟 =
𝜈𝛼
Le flux de chaleur échangé par convection sera caractérisé en le comparant à
un flux de chaleur de référence échangé par conduction. On définit ainsi un
nombre sans dimension, appelé nombre de Nusselt :
𝑓𝑙𝑢𝑥 𝑐𝑜𝑛𝑣𝑒𝑐𝑡𝑖𝑓 𝑙𝑜𝑖 𝑑𝑒 𝑁𝑒𝑤𝑡𝑜𝑛 ℎ𝑆Δ𝑇
𝑁𝑢𝐿𝑟𝑒𝑓 = = =
𝑓𝑙𝑢𝑥 𝑑𝑒 𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑖𝑜𝑛 𝑑𝑒 𝑟é𝑓é𝑟𝑒𝑛𝑐𝑒 𝑙𝑜𝑖 𝑑𝑒 𝐹𝑜𝑢𝑟𝑖𝑒𝑟 𝜆𝑆 Δ𝑇
𝐿𝑟𝑒𝑓
ℎ𝐿𝑟𝑒𝑓
𝑁𝑢𝐿𝑟𝑒𝑓 =
𝜆
Un nombre de Nusselt élevé signifiera donc que les échanges de chaleur par
convection prédominent face aux échanges par conduction.
Le nombre de Nusselt, bien entendu fait apparaître le coefficient d’échange convectif.
On pourra ainsi définir un nombre de Nusselt local à partir du coefficient d’échange
local associé au flux de chaleur échangé localement entre une paroi et le fluide ou
bien un nombre de Nusselt moyen défini à partir du coefficient d’échange moyen
associé au flux de chaleur global sur toute la surface de la paroi.
Exemple : écoulement le long d’une plaque plane de longueur 𝐿
 Nombre de Nusselt local en une position 𝑥 donnée le long de la plaque :
ℎ(𝑥)𝑥
𝑁𝑢𝑥 =
𝜆
 Nombre de Nusselt moyen calculé sur la longueur de la plaque :
ℎ̅𝐿 1 𝐿
𝑁𝑢𝐿 = 𝑜ù ℎ̅ = ∫ ℎ𝑥 𝑑𝑥
𝜆 𝐿 0
Remarque : le nombre de Reynolds caractéristique de l’écoulement en une position 𝑥
donnée le long de la plaque est (𝐿𝑟𝑒𝑓 = 𝑥) :
𝑈𝑥
𝑅𝑒𝑥 =
𝜈

41
 3.3.2 Lois de corrélation en convection
Le transfert de chaleur par convection dépend du régime d’écoulement (laminaire ou
turbulent) et de la nature du fluide.
En convection forcée, on cherchera donc à établir des corrélations qui relient le
nombre de Nusselt aux nombres de Reynolds et de Prandtl :
𝑁𝑢𝐿𝑟𝑒𝑓 = 𝑓(𝑅𝑒𝐿𝑟𝑒𝑓 , 𝑃𝑟)
Les principales corrélations sont présentées en annexe.
En convection naturelle, on cherchera à établir des corrélations qui relient le nombre
de Nusselt aux nombres de Grashof et de Prandtl. Les études montrent que les
corrélations s’écrivent simplement en utilisant le nombre de Rayleigh :

𝑁𝑢𝐿𝑟𝑒𝑓 = 𝑓 (𝑅𝑎𝐿𝑟𝑒𝑓 )

3.3.3 Méthodologie pour calculer le flux de chaleur en

convection

 Calcul du nombre de Reynolds, ReLref (convection forcée) ou du nombre de

Rayleigh, RaLref (convection naturelle), et du nombre de Prandtl Pr.
 Choix de la corrélation.
 Calcul du nombre de Nusselt :

𝑁𝑢𝐿𝑟𝑒𝑓 = 𝑓 (𝑅𝑒𝐿𝑟𝑒𝑓 , 𝑃𝑟) 𝑜𝑢 𝑁𝑢𝐿𝑟𝑒𝑓 = 𝑓(𝑅𝑎𝐿𝑟𝑒𝑓 )

 Calcul du coefficient d’échange (local ou moyen) :

𝜆
ℎ= 𝑁𝑢𝐿𝑟𝑒𝑓
𝐿𝑟𝑒𝑓
 Calcul du flux de chaleur (local ou global) par la loi de Newton.
Remarque : Dans le cas des conduites, la longueur caractéristique est le diamètre
hydraulique, défini par :
𝑆𝑒𝑐𝑡𝑖𝑜𝑛 𝑑𝑒 𝑝𝑎𝑠𝑠𝑎𝑔𝑒 𝑑𝑢 𝑓𝑙𝑢𝑖𝑑𝑒 𝑆
𝐷ℎ = 4 =4
𝑃é𝑟𝑖𝑚è𝑡𝑟𝑒 𝑚𝑜𝑢𝑖𝑙𝑙é 𝑑𝑒 𝑙𝑎 𝑐𝑜𝑛𝑑𝑢𝑖𝑡𝑒 𝑃
Exemples :
 Cylindre de diamètre 𝐷 totalement rempli de fluide :
𝐷2
𝑆=π 𝑒𝑡 𝑃 = 𝜋𝐷 ⇒ 𝐷ℎ = 𝐷
4
42
  Conduite rectangulaire de hauteur ℎ, de largeur 𝐿 totalement remplie de fluide :
𝑆 = ℎ𝐿 𝑒𝑡 𝑃 = 2(ℎ + 𝐿)
ℎ𝐿
𝐷ℎ = 2
ℎ+𝐿
Dans le cas où 𝐿 ≫ ℎ :
ℎ𝐿 ℎ
𝐷ℎ = 2 =2 ≈ 2ℎ
ℎ+𝐿 ℎ
+ 1
𝐿

3.3.4 Principales lois de correlation en convection forcée

𝑁𝑢 = 𝑓(𝑅𝑒, 𝑃𝑟)
Les propriétés thermo-physiques qui interviennent dans les nombres de Reynolds et
de Prandtl sont évaluées à la température moyenne entre la température de
l’écoulement et la température de surface du solide.
 Convection forcée interne
𝑆
Diamètre hydraulique 𝐷ℎ = 4 𝑃

𝑆 : section de passage du fluide

𝑃 : périmètre mouillé de la conduite

Tableau 3-1 : Ecoulement dans une conduite de diamètre hydraulique 𝐷ℎ

Remarque : pour un tube totalement rempli de fluide, une solution analytique peut être
obtenue en régime laminaire :
 Paroi chauffée à température constante : 𝑁𝑢𝐷 = 3,66
43
  Paroi chauffée à flux constant : 𝑁𝑢𝐷 = 4,36
 Convection forcée externe

ReD C m
0,4–4 0,989 0,33
4–40 0,911 0,385
40–4 103 0,683 0,466
4 103 -4 104 0,193 0,618
4 104- 4105 0,027 0,805

Tableau 3-2 : Ecoulement autour d’un cylindre de diamètre 𝐷

Tableau 3-3 :Ecoulement autour d’une sphère de diamètre 𝐷

𝜇∞ est la viscosité dynamique du fluide calculée à la température de l’écoulement à

l’infini (loin de la sphère).
𝜇𝑠 est la viscosité dynamique du fluide calculée à la température de surface de la
sphère.

44
CONCLUSION
La convection thermique est le mode transmission qui implique le déplacement d’un
fluide (liquide, gaz). Dans un fluide, il est pratiquement impossible d’assister à la
conduction pure car le moindre gradient de température entraîne des courants de
convection, c’est-à-dire un transport de masse.
L’étude du transfert de chaleur par convection a permis de déterminer les échanges
de chaleur se produisant entre un fluide et une paroi. La quantité de chaleur échangée
par unité de temps dépend de plusieurs paramètres à savoir la différence de
température entre la paroi et le fluide, la vitesse du fluide, la capacité thermique
massique du fluide, la surface d’échange, l’état de surface du solide ainsi que sa
dimension.

45
4. CHAPITRE 4 : TRANSFERT THERMIQUE PAR
RAYONNEMENT ELECTROMAGNETIQUE

46
INTRODUCTION
Tous les corps, quelque soit leur état (solide, liquide ou gazeux), émettent un
rayonnement de nature électromagnétique. Cette émission d’énergie s’effectue au
détriment de l’énergie interne du corps émetteur. Le rayonnement se propage de
manière rectiligne à la vitesse de la lumière. Il est de radiations de différentes
longueurs d’onde.
Le rayonnement électromagnétique (r.e.m.) obéit aux lois de la physique quantique et
peut être décrit soit par l’aspect corpusculaire soit par l’aspect ondulatoire, selon la
nature des phénomènes qui couplent ce rayonnement à la matière et selon les
énergies mises en jeu. Du point de vue ondulatoire, le rayonnement électromagnétique
(r.e.m.) résulte de la propagation d’une onde électromagnétique (propagation
simultanée d’un champ électrique et d’un champ magnétique) caractérisée par sa
fréquence. L’onde transporte une certaine quantité d’énergie liée à l’intensité du
champ électrique. Du point de vue corpusculaire, le rayonnement électromagnétique
(r.e.m.) est constitué de quantas (photons) d’énergie. Le flux d’énergie est d’autant
plus grand que le nombre de photons est élevé mais aussi que la fréquence est élevée.
C’est en fait l’énergie transportée par chaque photon qui conditionne la forme que
prend l’interaction du rayonnement avec la matière.

4.1 Généralités
Le rayonnement est un mode de transfert particulièrement important dans de
nombreuses situations industrielles à hautes températures (fours par exemple) mais
aussi géophysiques (refroidissement nocturne de la Terre). Le rayonnement est à la
base de nombreuses méthodes de mesures de température (caméra infra-rouge,
pyromètre à filament) qui sont non-intrusives et permettent d’opérer à distance, d’où
l’importance d’étudier soigneusement le rayonnement.
Le rayonnement thermique est de nature électromagnétique en raison de l’agitation de
la matière sous l’effet de la température. Une onde électromagnétique consiste en un
champ électrique ⃗E, un champ magnétique ⃗B et un vecteur d’onde ⃗k, tous trois
perpendiculaires. L’onde se propage dans la direction du vecteur d’onde avec une
C0
vitesse C qui varie selon le milieu qu’elle traverse. Cette vitesse vaut C = où C0 est
n

la vitesse de la lumière dans le vide et n l’indice du milieu. La fréquence d’une onde

47
électromagnétique ne varie pas avec le milieu qu’elle traverse. Elle est liée à C par la
relation :
𝑐
𝜆= = 𝑐𝑇
𝜈

1
où 𝜈 = 𝑇 est la fréquence de l’onde électromagnétique. On voit donc que la longueur

d’onde sera fonction du milieu traversé.

Une onde transporte avec elle de l’énergie. Lorsqu’elle est émise par un corps radiant,
celui-ci perd donc un certain flux d’énergie Φ (Watt). De même, un corps frappé par
une onde électromagnétique reçoit de l’énergie. La distribution énergétique d’un
rayonnement polychromatique (large gamme de λ) varie avec λ (figure 4-1).
On définit ainsi le flux monochromatique Φ (W.m−1) qui caractérise la contribution de
chaque longueur d’onde. On a, bien sûr :

𝜆
Φ = ∫ 𝜙𝜆 𝑑𝜆
0

Figure 4-1 : Spectre d’un rayonnement électromagnétique

Seules les longueurs d’onde comprises entre 0,4µm et 0,8µm sont visibles par l’œil.
La figure 4-2 illustre la variété de rayonnements existants.

48
Les rayons lumineux interagissent avec la matière de sorte que si certains corps la
laissent passer (corps transparents), certains l’absorbent en partie ou en totalité (corps
opaques). Pour un flux incident 𝜙𝑖 , on définit les quantités suivantes (cf. figure 4-3) :
 Flux réfléchi ϕr
 Flux absorbé ϕa
 Flux transmis ϕt
On appelle :
ϕr
 Réflectivité r =
ϕi
ϕa
 Absorptivité a = ϕi
ϕ
 Transitivité t = ϕt
i

Figure 4-2 : Divers rayonnements électromagnétiques classées selon leur longueur d’onde

Figure 4-3 : Définition des flux réfléchis, absorbés, transmis

La conservation de l’énergie impose :

𝑟+𝑎+𝑡 = 1

49
Ces grandeurs varient selon la valeur de λ. Ainsi le verre laisse passer les UV et est
opaque aux Infra-rouges. On doit donc définir rλ , aλ , t λ qui sont les réflectivités,
absorptivités et transitivités monochromatiques. Elles dépendent aussi de la direction
d’émission.

4.2 Le corps noir

4.2.1 Définition
C’est un corps idéal qui absorberait, s’il existait, tout rayonnement qu’il recevrait, quelle
que soit sa longueur d’onde. Pour qu’un corps noir reste en équilibre thermique (sa
température reste constante) il doit émettre également de l’énergie par rayonnement.
Son absorptivité est donc égale à 1.
L’intérêt du corps noir réside dans le fait qu’il sert de référence pour définir les
propriétés radiatives d’un corps réel. Il est important de noter qu’un corps noir n’est
pas forcément de couleur noire. Un corps de couleur noire est donc noir dans le visible
mais peut ne pas l’être pour d’autres longueurs d’onde. On peut réaliser
expérimentalement un corps noir en réalisant le montage de la figure 4-4 où une cavité
est percée d’un petit trou qui piège les rayons entrants. Elle est équipée de résistances
électriques pour réguler sa température. On doit donc retenir qu’un corps noir absorbe
tout rayonnement incident mais émet suivant une loi particulière: la loi de Stefan-
Boltzmann.

50
Figure 4-4: Réalisation expérimentale d’un corps noir

4.2.2 La loi de Stefan-Boltzmann

Définissons d’abord l’émittance du corps noir, M0, qui est la puissance totale émise par
une unité de surface, en intégrant sur toutes les longueurs d’onde. La loi de Stefan-
Boltzmann lie cette grandeur à la température du corps noir :
𝑀0 = 𝜎𝑇 4
où 𝑀0 est en W.m−2, T en Kelvin, et σ est la constante universelle de Stefan-
Boltzmann :

σ=5,670 10−8 W.m−2.K−4

51
 4.2.3 La loi de Planck
Cette loi fixe la contribution respective de chaque longueur d’onde dans l’émission du
corps noir. Elle s’écrit :
𝐶1
𝑀𝜆0 = 𝐶2
𝜆5 (𝑒 𝜆𝑇 − 1)
Avec : C1 = 3,74110−16W.m+2, C2 = 0,0143 m.K et où Mλ0 (en W. m−3 ) désigne
l’émittance monochromatique.

L’intégration de Mλ0 donne évidemment M0 :

∞
∫ 𝑀𝜆0 𝑑𝜆 = 𝜎𝑇 4 = 𝑀0
0

Les figures 25 montrent l’allure de Mλ0 en fonction de T pour diverses températures.

Elles passent toutes par un maximum.

4.3 Loi de Wien

4.3.1 Enoncé
∂M0λ
Cette loi fixe le lieu de ces maximas. En annulant on trouve:
∂λ

𝜆𝑚𝑎𝑥 𝑇 = 2898 µm. K

Où 𝜆𝑚𝑎𝑥 est la longueur d’onde correspondant au maximum d’énergie émise à la

température T. Le déplacement de 𝜆𝑚𝑎𝑥 vers les courtes longueurs d’onde quand T
augmente explique que, lorsqu’un corps chauffe, il commence d’abord par émettre
dans l’infra-rouge et n’est donc pas visible, puis son spectre vient empiéter sur la zone
rouge du spectre visible, puis couvre tout le visible et devient alors blanc. D’où
l’expression « chauffer au rouge » ou « chauffer à blanc ».

Il n’y a quasiment aucune zone commune entre le spectre du rayonnement solaire

(T≈5700 K) et celui d’un corps à température ambiante (300 K). Le soleil émet
principalement dans le visible et dans l’UV tandis qu’un corps à l’ambiante n’émet que
dans l’infra-rouge. Cette propriété est exploitée dans les serres ou les capteurs
solaires.

52
 4.3.2 Energie émise par le corps noir entre deux longueurs
d’onde
Pour un corps à une température donnée, on a souvent besoin d’évaluer la partie de
l’énergie qu’il émet dans une bande spectrale donnée (entre λ 1 et λ2). Cela revient à
calculer la grandeur adimensionnelle :
𝜆
∫𝜆 𝑀𝜆0 𝑑𝜆
2
1
𝐹𝜆1 −𝜆2 = ∞
∫0 𝑀𝜆0 𝑑𝜆 = 𝜎𝑇 4
que l’on peut écrire :
𝜆 𝜆
∫0 𝑀𝜆0 𝑑𝜆 − ∫0 𝑀𝜆0 𝑑𝜆
2 1

𝐹𝜆1 −𝜆2 = = 𝐹0−𝜆2 − 𝐹0−𝜆1

𝜎𝑇 4
Il serait donc utile de disposer d’une table fournissant, pour chaque température, la
valeur de l’intégrale :

𝜆
∫0 𝑀𝜆0 𝑑𝜆 1 𝜆 𝐶1
𝐹0−𝜆 = = ∫ 𝐶2 𝑑𝜆
𝜎𝑇 4 𝜎 0
𝜆5 𝑇 4 (𝑒 𝜆𝑇 − 1)

L’inconvénient est qu’il faudrait une table à 2 entrées (en T et en λ). Pour remédier à
ce problème on utilise 𝜆𝑇 comme variable unique et on pose :

1 𝜆𝑇 𝐶1 𝑑(𝜆𝑇)
𝐹0−𝜆𝑇 = ∫ 𝐶2
𝜎 0
(𝜆𝑇)5 (𝑒 𝜆𝑇 − 1)

53
Figure 4-5 : Emittance monochromatique du corps noir

54
Figure 4-6 : Courbe donnant la fraction d’énergie émise entre 0 et λT
On a ainsi une table avec une entrée unique sur la variable 𝜆𝑇. La courbe de la figure
4-6 montre le résultat obtenu.
On calcule donc l’énergie émise par un corps noir à la température T entre λ 1 et
λ2 par :
Fλ1−λ2 = F0−λ2T − F0−λ1T
Exemple : en supposant que le soleil est un corps noir à 5800 K, quelle est la fraction
de l’énergie qu’il émet dans le domaine visible (0,4–0,8 µm) ?
On calcule λ1T=2320 µm.K et λ2T=4640 µm.K.
Sur la figure 6 on lit :
F0−2320 =0,126 ; F0−4640 = 0,586 d’où :
Fλ2 −λ1 = 0,46
Donc, 46 % de l’énergie solaire émise est rayonnée dans le visible.

4.4 Corps réels

Nous venons d’étudier les lois d’émittance totale M 0 ou monochromatique Mλ0 du corps
noir. Nous ne nous sommes pour l’instant pas encore intéressés à l’importance de la
direction dans laquelle le rayonnement était émis par rapport à la normale à la surface
d’émission. Un corps noir rayonne de façon identique dans toutes les directions de
l’espace : on dit qu’il rayonne de manière diffuse (ou encore qu’il suit la loi de Lambert).
55
Un corps réel ne rayonne jamais vraiment de façon diffuse. Pour simplifier, nous
ignorerons cet aspect directionnel du rayonnement des corps réels étudiés. Par
ailleurs, l’èmittance M d’un corps réel n’est pas égale à celle, M 0 , du corps noir,
considéré comme le radiateur parfait. On pose :
M = εM 0
où ε d´désigne l’émissivité du corps à la température considérée. ε est donc inférieur
à 1 (cf. Figure 4-7). L’émissivité est aussi fonction de la longueur d’onde et il faut définir
l’émissivité monochromatique telle que :
Mλ = ελ Mλ0

Mais, pour simplifier, nous allons faire deux hypothèses :

 Les corps étudiés sont gris : ελ = ε

 Leur émission est diffuse (ou Lambertienne, ou encore isotrope) donc
indépendante de la direction d’émission.

Figure 4-7 : Courbe donnant la fraction d’énergie émise entre 0 et λT

4.5 Loi de Kirchhoff : relation entre absorption et émission

On admettra (sans démonstration) que pour λ donné, T donnée, et une direction
donnée, l’émissivité et l’absorptivité monochromatique directionnelle sont égales :

56
 εOx,λ = aOx,λ
La loi de Kirchhoff énonce donc qu’un corps, à température, longueur d’onde et
direction données, rayonne d’autant plus qu’il absorbe plus :
 Pour un corps gris, on aura : ελ = ε = aλ = a
 Pour un corps noir, on aura : ε = a = 1
L’émissivité des corps naturels dépend de leur nature physico-chimique, de leur état
de surface (poli, rugueux, strié) et varie aussi avec T. La mesure de l’émissivité est
difficile même en laboratoire, et on doit se contenter d’une précision souvent médiocre
sur la valeur de ε.

4.6 Flux radiatif échangé entre deux surfaces noires

Afin d’étudier de façon très générale comment un ensemble de plusieurs surfaces
échangent entre elles de l’énergie par rayonnement, considérons d’abord deux
surfaces comme le montre la Figure 4-8. La surface 𝑆1 envoie un flux total 𝜙1 = 𝑀10 𝑆1
mais seul le flux Φ12 touche la surface S1. Le rapport entre le flux total émis par S1 et
le flux qui frappe effectivement S2 s’appelle facteur de forme entre S1 et S2. Par
définition, ce nombre est compris entre 0 et 1 :
𝜙12
𝐹12 =
𝑀10 𝑆1
Ce facteur est fonction de la forme des surfaces S1 et S2, de leur distance, de l’angle
qu’elles font entre elles, etc. Il se calcule soit par une formule analytique, ce qui est
plus rapide, soit par utilisation d’abaques. De façon équivalente, on peut définir 𝜙12 qui
sera tel que :
𝜙12 = 𝐹12 𝜙1
Réciprocité :
On peut montrer que le flux Φ12 envoyé par S1 vers S2 est égal au flux Φ21 envoyé par
S2 vers S1. On en déduit une relation importante valable quelle que soient S1 et S2 :
F12 . S1 = F21 . S2

57
 S2
F 12
S1

Figure 4-8 : Echange radiatif entre S1 et S2

Relation entre les facteurs de forme dans une cavité fermée :

Soit une cavité fermée telle que montrée sur la Figure 4-9. Le flux Φi émis par Si est en
fait la somme des N flux que Si envoie vers les N−1 autres surfaces plus sur elle-même
si elle est concave, de sorte qu’on peut écrire :
𝜙𝑖 = 𝜙𝑖1 + ⋯ + 𝜙𝑖𝑖 + ⋯ + 𝜙𝑖𝑁

Soit encore :
N N N

ϕi = ∑ ϕij = ∑ Fij ϕi = ϕi ∑ Fij

j=1 j=1 j=1

On en déduit :
N

∑ Fij = 1
j=1

Qui ne traduit rien d’autre que la conservation du flux émis par chaque surface S i.

58
Figure 4-9 : Cavité composée de N surfaces noires
Bilan énergétique dans une cavité fermée : considérons de nouveau la cavité de la
figure 4-9. Le flux Net qu’elle reçoit, ou émet, est la somme algébrique du flux qu’elle
émet et des flux qu’elle reçoit. Afin d´exterminer la valeur de ce flux net, introduisons
la convention suivante :
 Un flux rayonné vers l’extérieur, encore dénommé flux partant, est compté
positivement
 Un flux venant frapper une surface est compté négativement.
Désignons par ϕi,n le flux net (en watt) qui frappe la surface S i. On peut écrire que ce
flux est la somme de ce que la surface Si émet et de ce qu’elle reçoit des N−1 autres
surfaces ainsi que d’elle-même si elle est concave, soit :
N N

ϕi,n = Mi0 . Si − ∑ Fji Mj0 Sj = Mi0 . Si − ∑ Fij Mj0 Si

j=1 j=1

En divisant par Si on obtient la densité de puissance nette φi,n sur Si :

N

φi,n = Mi0 − ∑ Fij Mj0

j=1

On notera les cas particuliers suivants :

 Cas d’une surface isotherme à la température 𝑇𝑖 ;on a alors 𝑀𝑖0 = 𝜎𝑇𝑖4

59
  Cas d’une surface adiabatique (isolée thermiquement) ; alors φi,n = 0
 Cas d’un flux imposé (source de puissance derrière ou dans la paroi par
exemple) alors φi,n = φ0
Abaques donnant les facteurs de forme :
Les Figures suivantes permettent de calculer les facteurs de forme entre surfaces dans
des géométries diverses.

Figure 4-10 : Facteur de forme entre surfaces rectangulaires parallèles

60
Figure 4-11 : Facteur de forme entre disques parallèles

Figure 4-12 : Facteur de forme entre surfaces rectangulaires perpendiculaires

61
Figure 4-13 : Facteur de forme entre cylindres concentriques

62
Figure 4-14 : Facteur de forme entre disques coaxiaux de diamètres différents

CONCLUSION
Le transfert d’énergie par rayonnement résulte des interactions énergétiques entre un
milieu matériel et le champ électromagnétique environnant et se produit par
l’intermédiaire d’ondes électromagnétiques, donc sans support matériel. Lorsque la
propagation se produit dans le vide ou dans un milieu parfaitement transparent, il n’y
a pas de dégradation de l’énergie transportée. La plupart des gaz simples entrent dans
cette catégorie. En revanche, lorsque la propagation des ondes électromagnétiques
s’accompagne d’une diminution de l’énergie, on parle de milieux semi-transparents.
Certains gaz, liquides ou solides entrent dans cette catégorie. Tout corps matériel émet
et absorbe de l’énergie sous forme de rayonnement électromagnétique. Le
phénomène d’émission d’un corps correspond à la conversion d’énergie matérielle en
énergie radiative. Le phénomène d’absorption est la conversion inverse.

63
CONCLUSION GENERALE

64
REFERENCES BIBLIOGRAPHIQUES

65