Cours de Transferts Thermiques

4ière Année ENSAM Meknès

Pr Mhamed MOUQALLID

Version 1

1
Table des matières
Introduction ...................................................................................................................................................... 4
I. Généralités sur les transferts de chaleur ............................................................................................ 4
I.1. Modes de transmission de la chaleur ............................................................................................... 4
I .1.2 – La conduction ............................................................................................................................. 5
I.1.3 – La convection ............................................................................................................................. 5
I.1.3 – Rayonnement .............................................................................................................................. 6
ière
1 partie : Transfert de chaleur par conduction ....................................................................................... 8
I. Généralités ............................................................................................................................................... 8
I. 1. Champ de température .................................................................................................................. 8
I.2. Surfaces isothermes ........................................................................................................................ 8
I.3. Flux thermiques ................................................................................................................................ 8
I.4. Loi de Fourrier .................................................................................................................................. 9
I.5. Conductivité thermique des matériaux........................................................................................ 11
I.6. La loi de Fourrier dans les corps anisotropes ............................................................................ 13
I. Equations du champ thermique ............................................................................................................... 14
2 - Equation thermique de la chaleur ..................................................................................................... 14
2 – 1. Equation thermique de la chaleur ................................................................................................ 16
3 - Conditions aux limites ....................................................................................................................... 16
4 - Solutions unidimensionnelles de l’équation de la chaleur ................................................................ 19
5 - Le problème de la tige mince de section constante .......................................................................... 26
6 - Ailettes de refroidissement ............................................................................................................... 31
7 - Conduction en régime instationnaire ................................................................................................ 36
2ième Partie : Convection ......................................................................................................................... 65
I. Introduction ............................................................................................................................................. 65
II. Analyse dimensionnelle appliquée à la convection................................................................................. 66
II.1 Exemples................................................................................................................................................ 67
A) Pertes de charge d’un tube lisse...................................................................................................... 67
B) Pertes de charge d’un tube rugueux ................................................................................................... 68
II.2 Expression adimensionnelle de la loi de convection ............................................................................. 69
II.3 Signification physique des nombres sans dimension ............................................................................ 71
III. Corrélations pour le calcul du coefficient de convection pour différentes configurations................. 74
3ième Partie : Rayonnement Thermique ...................................................................................................... 77
Introduction ..................................................................................................................................................... 77
I. Grandeurs énergétiques. .............................................................................................................................. 81

2
 I. 1. Flux total .............................................................................................................................................. 81
I. 2. Emittance ou Radiance M .................................................................................................................... 81
I. 3. Intensité I ............................................................................................................................................. 81
I. 4. Luminance............................................................................................................................................ 81
I. 5. Grandeurs monochromatiques ........................................................................................................... 82
I. 6. Eclairement .......................................................................................................................................... 82
I. 7. Relation entre l’éclairement du récepteur et la luminance de l’émetteur. ........................................ 83
II. Lois du rayonnement ................................................................................................................................... 83
II. 1. Loi de Lambert ..................................................................................................................................... 83
II. 2. Réception du rayonnement par un corps ........................................................................................ 84
II. 3. Corps noir ............................................................................................................................................ 84
II. 4. Loi de Planck ........................................................................................................................................ 85
II. 5. Loi de Stephan-Boltzman ..................................................................................................................... 87
II. 6. Facteur d’émission d’un corps non noir .............................................................................................. 88
II. 7. Corps en équilibre : Loi de Kirchhoff ................................................................................................... 88
II. 8. Corps gris ............................................................................................................................................. 90
II. 9. Modèles des « corps gris » et des corps « gris par morceaux » .......................................................... 90
II.8.1. Corps gris de rayonnement uniforme et isotrope ......................................................................... 90
II.8.1. Corps gris par morceau, de rayonnement uniforme et isotrope .................................................. 90
III. Echanges radiatifs entre surfaces séparées par un milieu parfaitement transparents ............................. 91
III. 1. Echanges radiatifs entre surfaces noires. ........................................................................................... 91
III. 1. 1. Facteurs de forme. ...................................................................................................................... 92
III. 1. 2. Puissance nette échangée par rayonnement entre S1 et S2 ....................................................... 93
III. 1. 3. Relations algébriques entre facteurs de forme .......................................................................... 94
III. 1. 4. Evaluation des facteurs de forme .............................................................................................. 95
III. 2. Echanges radiatifs entre surfaces opaques grises séparées par un milieu parfaitement transparent.
................................................................................................................................................................... 100
III. 1. 1. Surfaces opaques grises et diffusantes en émission et en réflexion. ....................................... 100
III. 1. 2. Echanges dans une enceinte vide aux parois grises et diffusantes en émission et réflexion. . 101
Bibliographie.................................................................................................................................................. 119

3
Introduction

La thermique est une discipline qui permet à l’ingénieur de comprendre le fonctionnement d’un
grand nombre d’équipements industriels tels que les générateurs de vapeurs, les fours, les
turbines à gaz et vapeur, les réacteurs nucléaires, l’industrie du froid, le chauffage, la climatisation,
etc….

Nous avons vu qu’en thermodynamique que l’énergie peut être transférée par interaction avec
l’environnement. Cette interaction a été appelée travail et chaleur.
Si la thermodynamique quantifie les quantités d’énergie mises en jeu dans l’évolution d’un
système, elle ne donne aucune indication sur les mécanismes et la dynamique de ces échanges.
C'est la tâche de la science du transfert thermiqu de clarifier les lois de ce processus.
Trois modes de transfert de chaleur peuvent être distingués : conduction, convection, et
rayonnement.

Le transfert de chaleur peut lui aussi se produire suivant deux modes semblables :

- Soit par contact : c’est la conduction thermique.

- Soit à distance : c’est le rayonnement thermique.

On considère un troisième mode de transfert d’énergie calorifique qui est la convection. Dans ce
cas, le phénomène thermique est compliqué par des déplacements de matière et au transfert de
chaleur se superpose le transfert de masse.

Nous allons étudier séparément, la conduction, la convection et le rayonnement, mais en réalité

les trois modes se présentent simultanément.

I. Généralités sur les transferts de chaleur

I.1. Modes de transmission de la chaleur

Deux corps à températures différentes échangent de l’énergie sous forme de chaleur jusqu’à
disparition complète de leur différence de température. Cet échange peut se faire suivant trois
modes de transmission fondamentaux : la conduction, la convection, le rayonnement.

Chacun de ces modes est lié à un phénomène physique bien déterminé ; En effet comme l’énergie
thermique d’un milieu matériel correspond à l’énergie cinétique de ses constituants fondamentaux
ayant une certaine liberté de mouvement (molécules, atomes, électrons libres, etc….). Ceux-ci
pourront échanger tout ou une partie de leur énergie thermique ; c'est-à-dire gagner ou perdre de
l’énergie cinétique :

- Soit par interaction directe avec les particules voisines (choc de molécules par exemple) ce
qui correspond à la conduction.
- Soit par absorption ou émission de radiations électromagnétiques, ce qui correspond
au rayonnement.

4
I .1.2 – La conduction

Le phénomène de conduction pure se rencontre essentiellement dans les solides. La propagation

de la chaleur par conduction à l’intérieur d’un corps se fait suivant deux mécanismes bien
distincts :

- Une transmission par les vibrations des atomes ou molécules (solides non métalliques).
- Une transmission par les électrons libres (métaux).
L’importance relative de chacun des deux mécanismes dépend évidemment de la nature du corps.

Exemple : La propagation de la chaleur le long d’une barre de fer dont une extrémité est soumise à
l’action d’une flamme.

La théorie de la conduction repose sur la loi expérimentale de Fourrier : la densité de flux est
proportionnelle au gradient de température :

Le flux de chaleur qui traverse une plaque d’épaisseur 𝐿 dont les

parois sont portées aux températures 𝑇1 𝑒𝑡 𝑇2 avec 𝑇1 > 𝑇2 est
donné par la loi de Fourrier :

𝑇1 − 𝑇2
= 𝑆
𝐿
 : Conductivité thermique
𝐿 : Epaisseur de la plaque

I.1.3 – La convection

La convection intervient en particulier dans les échanges thermiques entre une paroi et un fluide
en mouvement.

Lorsque le mouvement des fluides est provoqué par une action mécanique extérieure (action d’une
pompe pour la circulation d’un fluide dans une canalisation, agitation mécanique pour un fluide
chauffé dans un réacteur), la convection est dite forcée.
On parlera de la convection naturelle lorsque le mouvement des fluides est dû aux variations de
masse volumique au sein des fluides qui apparaissent sous l’effet des échanges thermiques entre
le fluide et la paroi par suite de différence de température.

Le mouvement résulte de la différence de masse volumique entre les parties chaudes et froides du
fluide.

5
Le flux de chaleur échangé par convection entre une paroi de
température 𝑇𝑃 , de surface 𝑆 et un fluide environnant de
température 𝑇𝑓 est régi par la loi de Newton :

 = 𝑕 𝑆 𝑇𝑃 − 𝑇𝑓

𝑕 : Coefficient de convection
𝑆 : Surface d’échange
𝑇𝑃 : Température de la paroi
𝑇𝑃 : Température du fluide environnant

Le tableau ci-dessous regroupe quelques valeurs typiques du coefficient de convection 𝑕

𝑕 (𝑊 𝑚−2 𝐾 −1 )

Convection libre 5 - 25
Convection forcée gaz 25 - 250
Convection forcée liquide 50 - 20000
Convection avec changement de 2500 - 100000
phase ( ébullition, condensation )
Tableau 1. Valeurs typiques du coefficient de convection

Exemple : Au contact d’un radiateur, la température de l’air croit et sa masse volumique décroît.
Cet air subit alors une poussée d’Archimède du fait de l’air environnant plus froid. Cet air est
remplacé par l’air plus froid qui s’échauffe à son tour et le processus continu. Cet exemple ne fait
intervenir que les forces de pesanteur : convection naturelle.

I.1.3 – Rayonnement

Tous les corps émettent par leur surface de l’énergie sous forme de radiations
électromagnétiques. Si cette énergie rencontre un corps absorbant, elle se transforme en chaleur.

La transmission de chaleur par rayonnement se produit même aux basses températures (il suffit
que la température soit supérieure à 0𝐾) dès que deux corps à des températures différentes sont
séparés par un milieu perméable. Ce transfert de chaleur qui ne nécessite pas de support matériel
se produit même dans le vide.

Exemple : Le rayonnement solaire nous parvient après avoir parcouru une distance considérable
dans le vide spatial.

La loi fondamentale du rayonnement thermique est la loi de Stéphan-Boltzmann :

𝜑 = 𝜍 𝑇4

𝜑 : Flux de chaleur émis par unité de surface du corps noir.

6
Un corps noir est un corps idéal qui émet le maximum d’énergie par rayonnement thermiqueà une
température donnée. C’est aussi un corps qui absorbe la totalité des radiations thermiques qu’il
reçoit.

𝜍 : Constante de Stephan-Boltzmann 𝜍 = 5.67 10 −8 𝑊𝑚−2 𝐾 −4

𝑇 : Température absolue

Quand le corps n’est pas noir il émet un flux par unité de surface :

𝜑 = 𝜀𝜍 𝑇 4

𝜀 : Propriété radiative de la surface (émissivité) avec 0 ≤ 𝜀 ≤ 1

Le flux de chaleur échangé entre une surface 𝑆 à la température 𝑇 de coefficient d’émission 𝜀 et un

milieu ambiant de température 𝑇𝑎 est donné par :

 = 𝜀𝜍𝑆 𝑇 4 − 𝑇𝑎4

7
 1ière partie : Transfert de chaleur par conduction

I. Généralités

I. 1. Champ de température

On peut définir en chaque point 𝑀 𝑥, 𝑦, 𝑧 d’un corps solide, liquide ou gazeux, une température,
fonction scalaire des coordonnées du point M et du temps : 𝑇 𝑥, 𝑦, 𝑧, 𝑡 .

Lorsque la température dépend du temps, on sait que le régime thermique est variable. Dans le
cas contraire, le régime est permanent.

I.2. Surfaces isothermes

Le lieu des points ayant à chaque instant la même température est appelé surface isotherme ou
isotherme 𝑇.

Le vecteur 𝑔𝑟𝑎𝑑 𝑇 défini en tout point M est perpendiculaire à la surface isotherme et dirigée vers
une surface à température supérieure à 𝑇.

I.3. Flux thermiques

Considérons une plaque chauffée de manière uniforme sur toute sa la surface 𝑆.

Soit 𝑑𝑄 la quantité de chaleur échangée entre cette plaque chauffante et l’air ambiant pendant le
temps 𝑑𝑡, on appelle :

- Flux thermique  : la puissance échangée par la surface 𝑆 de la plaque

𝑑𝑄
 =
𝑑𝑡

- Densité de flux thermique : la puissance échangée par une surface unité de la plaque
 𝑑𝑄
𝜑= =
𝑆 𝑆𝑑𝑡

D’une façon générale, pour un élément de surface 𝑑𝑆 orienté par sa

normale 𝑛, le flux élémentaire 𝑑 peut être considéré comme le flux
d’un vecteur 𝜑 à travers 𝑑𝑆 d’où :

𝑑 = 𝜑 ∙ 𝑛 𝑑𝑆

8
Le vecteur 𝜑 représente la densité locale du flux thermique au point M. Il caractérise en chaque
point du milieu, la direction, le sens et l’intensité du flux de chaleur.

L’ensemble des vecteurs 𝜑 constitue un champ de vecteurs analogue à ceux rencontrés dans
d’autres domaines. Par analogie (champs électriques, champs de forces, vitesses,……..).

𝑄 représente une énergie 𝑄 = 𝐽𝑜𝑢𝑙𝑒 (𝐽)

 représente une énergie par unité de temps  = Watts (W)

𝜑 représente une puissance par unité de surface 𝜑 = 𝑊𝑎𝑡𝑡𝑠 𝑚−2

I.4. Loi de Fourrier

Enoncé : Il existe une relation linéaire entre la densité de flux thermique et le gradient de
température. En tout point d’un milieu isotrope, la densité de flux thermique instantanée est
proportionnelle à la conductivité thermique  du milieu et au gradient de température :

𝜑 = −  𝑔𝑟𝑎𝑑 𝑇

C’est une loi qui est basée sur les études expérimentales.

Par convention 𝜑 est compté positivement dans le sens de l’écoulement de la chaleur, c'est-à-dire
vers les températures décroissantes.

𝑔𝑟𝑎𝑑 𝑇 est de sens contraire de 𝜑, mais porté par le même axe.

En un point quelconque du milieu, on a à tout instant :

𝑑𝑇 = 𝑔𝑟𝑎𝑑 𝑇 ∙ 𝑑𝑀

Si on considère un déplacement élémentaire sur l’isotherme passant par 𝑀, on a :

𝑔𝑟𝑎𝑑 𝑇 ∙ 𝑑𝑀 = 0

𝑔𝑟𝑎𝑑 𝑇 est en chaque point normal à la surface isotherme passant par ce point.

De point de vue physique, le gradient de température représente le taux de variation de la

température suivant la direction normale à l’isotherme ou dérivée normale à l’isotherme :

𝜕𝑇
𝜕𝑛 𝑆𝑢𝑟𝑓𝑎𝑐𝑒 𝑖𝑠𝑜𝑡 𝑕𝑒𝑟𝑚𝑒

La conductivité thermique peut varier d’un point à un autre d’un corps, soit directement dans les
solides non homogènes,  𝑥, 𝑦, 𝑧 , soit indirectement à travers la température, dans les solides où
la conductivité thermique varie fortement avec la température  𝑇

9
On pour les projections du vecteur densité de flux 𝜑 sur le système de coordonnées :

𝜕𝑇
𝜑𝑥 = −
𝜕𝑥
𝜕𝑇
𝜑𝑦 = −
𝜕𝑦
𝜕𝑇
𝜑𝑧 = −
𝜕𝑧

Pour un milieu isotrope, la conductivité thermique  est une grandeur scalaire positive fonction du
point 𝑀 et de la température 𝑇,  𝑀, 𝑇 . Le tableau 1 montre des valeurs de la conductivité
thermique de quelques substances.

En général, la conductivité thermique des métaux dits purs (cuivre, fer, zinc, magnésium, etc..)
décroit avec quand la température augmente. Certains métaux tels que le béryllium, le cobalt, le
platine voient au contraire leurs conductivité thermiques croitre avec la température. Beaucoup de
métaux comme le plomb, le bismuth, le manganèse voient leurs conductivité thermique passer
par un minimum.
Pour de nombreux métaux solides (métalliques ou non) et liquides, il est acceptable de supposer
une conductivité thermique constante avec la température, soit une conductivité thermique qui
varie suivant la loi linéaire :
 = 0 1 + 𝑏 𝑇 − 𝑇0

0 est la conductivité thermique à la température 𝑇0 .

𝑏 est une constante à déterminer expérimentalement.

On verra que prendre  = 𝑇 rend l’équation de la conduction non linéaire.

Matériaux  (𝑊𝑚−1 °𝐶 −1 ) Matériaux  (𝑊𝑚−1 °𝐶 −1 )

Argent 419 Plâtre 0.48
Cuivre 386 Amiante 0.16
Aluminium 204 Bois 0.12 – 0.23
Acier doux 45 Liège 0.044 – 0.049
Acier inox 15 Laine de roche 0.038 – 0.041
Glace 1.88 Laine de verre 0.035 – 0.051
Béton 1.4 Polystyrène expansé 0.036 – 0.047
Brique terre cuite 1.1 Polyuréthane 0.030 – 0.045
Verre 1.0 Polystyrène extrudé 0.028
Eau 0.6 Air 0.026
Tableau 2. Conductivité thermique de quelques substances

La figure 1 montre l’évolution de la conductivité thermique de quelques matériaux en fonction de la

température.

10
 Fig. 1. Variation de la conductivité thermique des matériaux en fonction de la température

I.5. Conductivité thermique des matériaux

La conductivité thermique des matériaux homogènes dépend de nombreux facteurs et notamment

de leur composition chimique, de leur état (solide, liquide ou gazeux) de leur structure cristalline
s’il s’agit d’un solide, de leur température et pour les gaz et les vapeurs leur pression.

D’après la relation :

𝜑 = −  𝑔𝑟𝑎𝑑 𝑇

𝜑 = 𝑊𝑎𝑡𝑡𝑠 𝑚−2

 = 𝑊𝑎𝑡𝑡𝑠 𝑚−1 °𝐶 −1

a) Conduction thermique dans les solides

Nous avons signalé deux mécanismes pour la conduction dans les solides :

- une transmission de chaleur par les vibrations des atomes ou molécules que l’on
caractérise par un coefficient de conductibilité thermique 𝑣 .
- Une transmission de chaleur par les électrons libres caractérisé par un coefficient de
conductibilité thermique 𝑒 .
La conductibilité thermique  d’un corps est donnée par :

 = 𝑣 + 𝑒

a-1. Conduction thermique dans les métaux purs

11
Dans les métaux purs, bons conducteurs de l’électricité, le terme 𝑒 est le facteur primordial. La
conductivité thermique et la conductivité électrique sont extrêmement liées par la loi expérimentale
de Wiedman et Franz (1853).

Le rapport de la conductibilité thermique  à la conductibilité électrique 𝜍 est le même pour tous

les métaux à une température donnée.
D’après Lorenz, ce rapport est proportionnel à la température absolue.

= 𝐿0 𝑇
𝜍

𝐿0 = 2.45 10−8 𝑊  𝐾 −2

En raison des impuretés et défauts du réseau cristallin, les valeurs de  sont très variables.

En conséquence, pour les métaux industriels, il ne sera pas possible d’utiliser les valeurs des
conductivités données pour les métaux purs.

Métaux à 20°C  𝐾𝑐𝑎𝑙 𝑕−1 𝑚−1 °𝐶 −1

𝐴𝑔 365
𝐶𝑢 330
𝐴𝑙 196

a-2. Conductivité thermique dans les alliages

Dans les alliages ou la conductivité thermique dépend de la nature des constituants et de la

structure de l’alliage, la contribution relative au terme v est plus importante que les métaux purs.

Les structures des alliages étant moins ordonnée que celle des métaux, la conductivité thermique
sera plus faible.

Acier à 20°C  𝐾𝑐𝑎𝑙 𝑕−1 𝑚−1 °𝐶 −1

Acier au chrome 18 <  < 28
Acier au chrome, nickel 15 <  < 48
Acier doux 40

a-3. Conductivité thermique des liquides

Au point de fusion d’un corps, on observe généralement une chute brutale de la conductivité
thermique par suite essentiellement de la destruction de l’organisation de la structure. Par contre
les métaux fondus représentent une classe de liquides à part : leur conductivité thermique est due
essentiellement aux électrons libres et est 10 à 1000 fois plus grande que celle des liquides non
métalliques.

 𝐾𝑐𝑎𝑙 𝑕−1 𝑚−1 °𝐶 −1

A 30°C eau 0.54
A 30° Benzène 0.14
A 98°C Sodium fondu 74

12
a-4. Conductivité thermique des gaz

La conductivité thermique des gaz s’explique par la théorie cinétique des gaz. On démontre que la
conductivité thermique est donnée par :

1
 = 𝜌 𝑢 𝐶𝑣 𝐿
3

𝑢 ∶ 𝑣𝑖𝑡𝑒𝑠𝑠𝑒 𝑚𝑜𝑦𝑒𝑛𝑛𝑒 𝑑’𝑎𝑔𝑖𝑡𝑎𝑡𝑖𝑜𝑛 𝑑𝑒𝑠 𝑚𝑜𝑙é𝑐𝑢𝑙𝑒𝑠

𝐿: 𝑙𝑖𝑏𝑟𝑒 𝑝𝑎𝑟𝑐𝑜𝑢𝑟𝑠 𝑚𝑜𝑦𝑒𝑛
𝐶𝑣 ∶ 𝑐𝑕𝑎𝑙𝑒𝑢𝑟 𝑚𝑎𝑠𝑠𝑖𝑞𝑢𝑒 à 𝑣𝑜𝑙𝑢𝑚𝑒 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡
𝜌 ∶ 𝑚𝑎𝑠𝑠𝑒 𝑣𝑜𝑙𝑢𝑚𝑖𝑞𝑢𝑒
On peut aussi relier  au coefficient de viscosité dynamique par :
 = 𝑎 𝐶𝑣 𝜇

𝑎 = 2.5 𝑝𝑜𝑢𝑟 𝑙𝑒𝑠 𝑔𝑎𝑧 𝑚𝑜𝑛𝑜𝑎𝑡𝑜𝑚𝑖𝑞𝑢𝑒𝑠

1 < 𝑎 < 2 𝑝𝑜𝑢𝑟 𝑙𝑒𝑠 𝑔𝑎𝑧 𝑝𝑜𝑙𝑦𝑎𝑡𝑜𝑚𝑖𝑞𝑢𝑒𝑠

 𝐾𝑐𝑎𝑙 𝑕−1 𝑚−1 °𝐶 −1

A 0°C Air 0.0192
Ammoniac 0.0165
Oxyde de carbone 0.0121

I.6. La loi de Fourrier dans les corps anisotropes

Le développement précédent de la loi de Fourrier est valable pour les corps anisotropes. Pour les
corps anisotropes, on suppose que les projections du vecteur densité de flux, selon chacune des
directions, sont fonctions des gradients de température suivant les trous directions. On donc :

𝜕𝑇 𝜕𝑇 𝜕𝑇
𝜑𝑥 = − 𝑥𝑥 + 𝑥𝑦 + 𝑥𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑇 𝜕𝑇 𝜕𝑇
𝜑𝑦 = − 𝑦𝑥 + 𝑦𝑦 + 𝑦𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑇 𝜕𝑇 𝜕𝑇
𝜑𝑧 = − 𝑧𝑥 + 𝑧𝑦 + 𝑧𝑧
𝜕𝑥 𝜕𝑦 𝜕𝑧

Il existe donc trois coefficients 𝑖𝑗 suivant chaque direction du système de cordonnées. La
conductivité est alors un tenseur symétrique du second ordre.

𝑥𝑥 𝑥𝑦 𝑥𝑧

𝑖𝑗 = 𝑦𝑥 𝑦𝑦 𝑦𝑧
𝑧𝑥 𝑧𝑦 𝑧𝑧

13
I. Equations du champ thermique

2 - Equation thermique de la chaleur

On considère un repère 𝑂𝑥𝑦𝑧. Choisissons un parallélépipède élémentaire de volume 𝑑𝑥𝑑𝑦𝑑𝑧 et

faisons un bilan thermique. On suppose constante la conductivité thermique du matériau.

Le flux de chaleur à travers la surface 𝑑𝑥𝑑𝑧 est :

𝜕𝑇
−  𝑑𝑥𝑑𝑧
𝜕𝑦 𝑦
En 𝑦 + 𝑑𝑦 on a :
𝜕𝑇 𝜕𝑇 𝜕 2 𝑇
−  𝑑𝑥𝑑𝑧 = −  𝑑𝑥𝑑𝑧 + 𝑑𝑦
𝜕𝑦 𝑦+𝑑𝑦 𝜕𝑦 𝜕𝑦 2
De même pour les flux qui traversent 𝑑𝑧𝑑𝑦 et 𝑑𝑦𝑑𝑧.

Le solde des actions calorifiques sur le volume élémentaire par unité de temps à travers
l’ensemble des parois a pour expression :

𝜕2 𝑇 𝜕2 𝑇 𝜕2 𝑇
 𝑑𝑥𝑑𝑦𝑑𝑧 + +
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2

On suppose que cette action calorifique a pour effet d’accroitre l’énergie interne de ce volume à
raison de :

𝜕𝑇
𝜌𝐶 𝑑𝑥𝑑𝑦𝑑𝑧
𝜕𝑡

𝐶 : Chaleur massique à volume constant du matériau.

𝜌 : masse volumique du matériau.
D’où l’équation de la chaleur :
𝜕𝑇 𝜕2 𝑇 𝜕2 𝑇 𝜕2 𝑇
𝜌𝐶 =  + + + 𝑝
𝜕𝑡 𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2

𝜌𝐶 𝜕𝑇 𝜕2 𝑇 𝜕2 𝑇 𝜕2 𝑇 𝑝
= + + 2 +
 𝜕𝑡 𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 

14
 1 𝜕𝑇 𝜕2 𝑇 𝜕2 𝑇 𝜕2 𝑇 𝑝
= + + 2 +
𝑎 𝜕𝑡 𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 

1 𝜕𝑇 𝑝
= ∆𝑇 +
𝑎 𝜕𝑡 

𝑝 : est la puissance générée par les sources internes.

𝑝 ∶ 𝑊 𝑚3

𝑎=
𝜌𝐶
𝑎 : est la diffusivité thermique du matériau (caractérise le comportement transitoire du matériau).
𝑎 ∶ 𝑚2 𝑠

Air Laine de Bois Verre Béton Acier Cuivre

verre plein doux
𝑎 𝑚2 𝑠 2.2 10−5 2. 10−6 1.5 10−7 7.5 10−7 7 10−7 1.8 10−5 1.2 10−4
Tableau 3. Diffusivité thermique de quelques corps à la température ambiante

Equation de la chaleur en coordonnées cylindriques :

1 𝜕𝑇 1 𝜕 𝜕𝑇 1 𝜕2 𝑇 𝜕2 𝑇 𝑝
= 𝑟 + 2 2+ 2+
𝑎 𝜕𝑡 𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧 

Equation de la chaleur en coordonnées sphériques :

1 𝜕𝑇 1 𝜕 2 𝜕𝑇 1 𝜕 𝜕𝑇 1 𝜕2 𝑇 𝑝
= 2 𝑟 + 2 𝑠𝑖𝑛𝜃 + 2 2 +
𝑎 𝜕𝑡 𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝑠𝑖𝑛𝜃 𝜕𝜃 𝜕𝜃 𝑟 𝑠𝑖𝑛 𝜃 𝜕𝜑2 

Autre démonstration :

Production = Echange + Accumulation

𝜕𝑇
𝑝𝑑𝑣 = 𝜑. 𝑛 𝑑𝑆 + 𝜌𝐶 𝑑𝑣
𝜕𝑡

On utilise la formule de Green Ostogradski

𝜑. 𝑛 𝑑𝑆 = 𝑑𝑖𝑣 𝜑 𝑑𝑣

L’équation devient :

𝜕𝑇
𝑝𝑑𝑣 = 𝑑𝑖𝑣 𝜑 𝑑𝑣 + 𝜌𝐶 𝑑𝑣
𝜕𝑡

𝜕𝑇
𝑝 = 𝑑𝑖𝑣 𝜑 + 𝜌 𝐶
𝜕𝑡

Avec :

𝜑 = − 𝑔𝑟𝑎𝑑 𝑇
15
 𝜕𝑇
𝑝 = 𝑑𝑖𝑣 − 𝑔𝑟𝑎𝑑 𝑇 + 𝜌 𝐶
𝜕𝑡
𝑝 1 𝜕𝑇
∆𝑇 + =
 𝑎 𝜕𝑡

Cas particuliers

 Champ instationnaire
𝜕𝑇
=0
𝜕𝑡

On obtient l’équation de Poisson

𝑝
∆𝑇 + =0

 Champ instationnaire, sans production (𝑝 = 0).

On obtient l’équation de Laplace :

∆𝑇 = 0

Dans les cas où la conductivité thermique n’est pas constante on refait le bilan. Le flux de chaleur
pour la direction 𝑦 :

𝜕 𝜕𝑇
 𝑦 + 𝑑𝑦 =  𝑦 + 𝑑𝑥𝑑𝑧 − 𝑑𝑦
𝜕𝑦 𝜕𝑦

De même pour les autres directions et on obtient l’équation suivante :

𝜕𝑇 𝜕 𝜕𝑇 𝜕 𝜕𝑇 𝜕 𝜕𝑇
𝜌𝐶𝑣 =  +  +  + 𝑝
𝜕𝑡 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧

Pour les coordonnées cylindriques l’équation s’écrit :

𝜕𝑇 1 𝜕 𝜕𝑇 1 𝜕 𝜕𝑇 𝜕 𝜕𝑇
𝜌𝐶𝑣 = 𝑟 + 2  +  + 𝑝
𝜕𝑡 𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝜃 𝜕𝑧 𝜕𝑧

Pour les coordonnées sphériques l’équation s’écrit :

𝜕𝑇 1 𝜕 𝜕𝑇 1 𝜕 𝜕𝑇 1 𝜕 𝜕𝑇
𝜌𝐶𝑣 = 2  𝑟2 + 2  𝑠𝑖𝑛𝜃 + 2 2  + 𝑝
𝜕𝑡 𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝑠𝑖𝑛𝜃 𝜕𝜃 𝜕𝜃 𝑟 𝑠𝑖𝑛 𝜃 𝜕𝜑 𝜕𝜑

2 – 1. Equation thermique de la chaleur

3 - Conditions aux limites

Comme on le sait, la solution d’une équation aux dérivées partielles comme l’équation de la
chaleur dépend dans une large mesure des conditions initiales et des conditions aux limites.
16
Lorsque le champ thermique n’est pas stationnaire, il y a bien lieu de fixer de façon détaillée la
distribution initiale des températures.

à 𝑡 = 0, 𝑇 𝑥, 𝑦, 𝑧, 𝑡 = 𝑇0 𝑥, 𝑦, 𝑧

Dans tous les cas, il faut imposer aux différentes frontières du système étudié une des conditions
aux limites que nous allons passe en revue :

3.a. Température imposée

La température est une fonction donnée du point 𝑀𝑝 sur la paroi et/ou du temps :

𝑇𝑝 = 𝑓 𝑀𝑝 , 𝑡

C’est une condition pratique mais peu réaliste (c’est la condition de Dirichlet 1ière espèce).

Si 𝑇𝑝 est une constante, c’est la condition d’isothermie (réalisée par exemple dans les
thermostats).

3.b. densité de flux imposée

La densité de flux thermique est une fonction donnée du point sur la paroi et/ou du temps, soit :

𝜕𝑇
𝜑 = − = 𝑓 𝑀𝑝 , 𝑡
𝜕𝑛 𝑝

𝜕𝑇
∶ 𝑑é𝑠𝑖𝑔𝑛𝑒 𝑙𝑎 𝑑é𝑟𝑖𝑣é𝑒 𝑛𝑜𝑟𝑚𝑎𝑙𝑒 𝑑𝑒 𝑇, 𝑐𝑎𝑙𝑐𝑢𝑙é𝑒 à 𝑙𝑎 𝑝𝑎𝑟𝑜𝑖
𝜕𝑛 𝑝

(dérivée de la température le long de la normale).

C’est une condition très réaliste (c’est la condition de Newmann 2ière espèce).

Si le corps est thermiquement isolé, le flux est nul en tout point de sa surface.

3.c. Paroi en contact d’un fluide en mouvement

La densité de flux traversant la surface frontière est proportionnelle à la différence de température

entre la paroi et le milieu environnant (généralement fluide) :

𝜑 = 𝑕 𝑇𝑝 − 𝑇∞

𝑕 est appelé coefficient de convection.

𝑇∞ 𝑒𝑠𝑡 𝑙𝑎 𝑡𝑒𝑚𝑝é𝑟𝑎𝑡𝑢𝑟𝑒 𝑑𝑢 𝑓𝑙𝑢𝑖𝑑𝑒 𝑙𝑜𝑖𝑛 𝑑𝑒 𝑙𝑎 𝑝𝑎𝑟𝑜𝑖.

Cette condition suppose que la variation de température de 𝑇𝑝 à 𝑇∞ est confinée dans l’épaisseur
d’une couche fluide très mince, en contact avec la paroi solide. Au delà le fluide n’intervient plus.
C’est la condition courante de transfert de chaleur par convection dite Loi de Newton.
En utilisant la loi de Fourrier, cette condition s’écrit :

𝜕𝑇
− = 𝑕 𝑇𝑝 − 𝑇∝
𝜕𝑛 𝑝

17
C’est la condition de Fourrier de 3ième espèce.

La condition de newton n’implique pas que la température de la paroi soit uniforme, même si celle
du fluide l’est.

3.d. Paroi échangeant de la chaleur avec son voisinage par rayonnement

Un corps de température absolue 𝑇1 placé dans une enceinte à la température absolue 𝑇2

échange avec celle-ci la densité de flux thermique :

𝜑 = 𝐴𝜍 𝑇14 − 𝑇24

A est un paramètre qui tient compte des propriétés radiatives des corps en présence et de la
géométrie du problème.

𝜍 𝑒𝑠𝑡 𝑙𝑎 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 𝑑𝑒 𝑆𝑡é𝑝𝑕𝑎𝑛 − 𝐵𝑜𝑙𝑡𝑧𝑚𝑎𝑛

𝜍 = 5.67 10−8 𝑊/𝑚2 𝐾 4

La condition de transfert par rayonnement s’écrit :

𝜕𝑇
− = 𝐴𝜍 𝑇𝑝4 − 𝑇𝑒4
𝜕𝑛 𝑝

Cette condition n’est pas linéaire puisque la température intervient en puissance quatre.

Pour un écart 𝑇𝑝 − 𝑇𝑒 faible la formule peut être linéarisée :

𝑇𝑝4 − 𝑇𝑒4 = 𝑇𝑝2 + 𝑇𝑒2 𝑇𝑝 + 𝑇𝑒 𝑇𝑝 − 𝑇𝑒

Pour 𝑇𝑝 − 𝑇𝑒 ≤ 100°𝐶

𝜕𝑇
𝜑𝑟 = −  = 𝑕𝑟 𝑇𝑝 − 𝑇𝑒
𝜕𝑛 𝑝

Avec 𝑕𝑟 = 4𝐴𝜍𝑇𝑒∗3

𝑇 ∗ représente la température moyenne entre 𝑇𝑝 𝑒𝑡 𝑇𝑒 .

Par analogie avec la convection, le paramètre 𝑕𝑟 est appelé coefficient d’échange thermique par
rayonnement.

3.e. Transfert à l’interface de deux solides de nature différentes

Lorsque deux solides de conductivités respectives 1 𝑒𝑡 2 possédant une frontière commune S,

la conservation du flux s’écrit :

1 𝑔𝑟𝑎𝑑 𝑇1 = 2 𝑔𝑟𝑎𝑑 𝑇2

C’est la condition de Fourrier de 3ième espèce

18
La seconde condition est donnée dans le cas du contact parfait par l’égalité des températures des
deux corps à l’interface : 𝑇1 = 𝑇2 𝑠𝑢𝑟 𝑆.

4 - Solutions unidimensionnelles de l’équation de la chaleur

Dans de nombreux problèmes, la géométrie du solide est très simple. Si les conditions aux limites
le sont également, l’équation fondamentale se simplifie beaucoup.

4.1. La plaque plane simple

Considérons une plaque plane à faces parallèles, d’étendue

infinie, d’épaisseur 𝐿, constituée d’un matériau homogène
isotrope.

a) Les deux faces de la plaque sont maintenues à des

températures uniformes 𝑇1 𝑒𝑡 𝑇2 .

𝜕𝑇
Si nous limitons à l’étude d’un champ stationnaire 𝜕𝑡
=0
et sans source interne de la chaleur 𝑝 = 0 . Du fait que les axes 𝑂𝑦 𝑒𝑡 𝑂𝑧 sont parallèles aux
plans isothermes, les variables 𝑦 𝑒𝑡 𝑧 ne jouent aucun rôle et l’équation du champ thermique
s’écrit :

𝑑2 𝑇
=0
𝑑𝑥 2

Les conditions aux limites sont :

𝑇 𝑥 = 0 = 𝑇1

𝑇 𝑥 = 𝐿 = 𝑇2

La solution de l’équation est :

𝑇2 − 𝑇1
𝑇 𝑥 = 𝑇1 + 𝑥 (∗∗∗)
𝐿

La loi de répartition de la température est linéaire.

Le vecteur densité de flux est :

𝑇1 − 𝑇2
𝜑 = − 𝑔𝑟𝑎𝑑 𝑇 =  𝑥
𝐿

Si 𝑇1 > 𝑇2 , 𝜑 est suivant l’axe 𝑥 qui l’axe de propagation de la chaleur.

L’équation (***) montre que les plans 𝑥 = 𝑐𝑠𝑡𝑒 sont des plans isothermes.

Si on considère une surface S de cette paroi, le flux de chaleur qui la traverse est :

19
 𝑇1 − 𝑇2
∅= 𝜑𝑆=  𝑆
𝐿
𝐿
𝑇1 − 𝑇2 = ∅
𝑆
Dans l’analogie électrique 𝑇1 − 𝑇2 est analogue à la tension 𝑉1 − 𝑉2 et ∅ à un courant électrique 𝐼
𝑉1 − 𝑉2 = 𝑅𝐼 (Loi d’Ohm).
𝐿
Donc 𝑆
est analogue à une résistance électrique, on définit donc la résistance thermique
conductive de la plaque par :
𝐿
𝑅𝑐𝑑 =
𝑆

b) La plaque sépare deux fluides à des températures

extérieures uniformes.

Entre les parois et les fluides s’établit un échange

convectif. On admet que les coefficients de convection
𝑕1 𝑒𝑡 𝑕2 des fluides 1 et 2 au contact de la plaque
aient une même valeur en tout point de cette paroi.

Le problème étant unidimensionnel, la conservation du flux se traduit par l’égalité des flux cédé par
le fluide chaud au mur, traversant le mur, et reçu par le fluide froid :
𝑇1 − 𝑇2
∅ = 𝜑 𝑆 = 𝑕1 𝑆 𝑇∞1 − 𝑇1 =  𝐿
𝑆 = 𝑕2 𝑆 𝑇2 − 𝑇∞2

On définit ainsi les résistances thermiques convectives :

1 1
𝑅𝑐𝑣1 = 𝑒𝑡 𝑅𝑐𝑣2 =
𝑕1 𝑆 𝑕2 𝑆
D’ou :
𝑇∞1 − 𝑇∞2
∅=
1 𝐿 1
+ 𝑆 +
𝑕1 𝑆 𝑕2 𝑆
Cette définition traduit la loi d’Ohm pour les trois résistances 𝑅𝑐𝑑 , 𝑅𝑐𝑣1 𝑒𝑡 𝑅𝑐𝑣2 en série.

4.2. La paroi plane composée

Si la paroi comprend 𝑛 plaques en contact, d’épaisseurs successives 𝑙1 , 𝑙2 , 𝑙3 , … … . . 𝑙𝑛 et de

conductivités thermiques respectives 1 , 2 , 3 , … … . . 𝑛 et si comme précédemment, nous
admettons que les coefficients de convection des fluides 1 et 2 au contact de la paroi ont même
valeur en tout point de cette paroi.

𝑇1 , 𝑇2 , 𝑇, … … . . 𝑇𝑛 , 𝑇𝑛+1 désignent les températures des surfaces de contact intermédiaires.

La conservation du flux thermique le long d’un tube de courant de section S implique les égalités :

1 𝑛
∅ = 𝑕1 𝑆 𝑇∞1 − 𝑇1 = 𝑇1 − 𝑇2 = ⋯ 𝑇𝑛 − 𝑇𝑛+1 = 𝑕2 𝑆 𝑇𝑛+1 − 𝑇∞2
𝑙1 𝑙𝑛
20
On obtient :
𝑇∞1 − 𝑇∞2
∅=
1 1
+ 𝑛1 𝑅𝑖 +
𝑕1 𝑆 𝑕2 𝑆
Avec
𝑙𝑖
𝑅𝑖 =
𝑖 𝑆

4.3. La paroi cylindrique simple

Considérons une paroi cylindrique de longueur axiale infinie, de

rayon intérieur 𝑅1 et de rayon extérieur 𝑅2 , constitué d’un matériau
homogène et isotrope dont les surfaces cylindriques sont à des
températures uniformes 𝑇1 𝑒𝑡 𝑇2 . Les isothermes sont des surfaces
cylindriques coaxiales, le gradient de température est radial.

Si l’on rapporte cette paroi à un système d’axes de cordonnées

cylindriques

𝑒𝑟 , 𝑒𝜃 , 𝑒𝑧 , l’équation du champ thermique s’écrit (en régime

𝜕𝑇
stationnaire 𝜕𝑡
= 0 et sans source interne de la chaleur 𝑝 = 0 )
s’écrit :

𝑑2 𝑇 1 𝑑𝑇
+ =0
𝑑𝑟 2 𝑟 𝑑𝑟

Les conditions aux limites sont :

𝑇 𝑟 = 𝑅1 = 𝑇1

𝑇 𝑟 = 𝑅2 = 𝑇2

La résolution de cette équation donne :

𝑇2 − 𝑇1 𝑟
𝑇 𝑟 = 𝑇1 + 𝐿𝑛
𝑅 𝑅1
𝐿𝑛 𝑅1
2

21
On déduit de la loi de Fourrier, la densité de flux :

𝑇1 − 𝑇2 1
𝜑 = − 𝑔𝑟𝑎𝑑 𝑇 = − 𝑒𝑟
𝑅
𝐿𝑛 𝑅1 𝑟
2

Le flux thermique ∅ traversant chaque surface cylindrique de longueur L est :

𝑇1 − 𝑇2
∅ = 2𝜋 𝐿
𝑅
𝐿𝑛 𝑅2
1

𝑅
𝐿𝑛 𝑅2
1
𝑇1 − 𝑇2 = ∅
2𝜋 𝐿

D’où la résistance thermique de conduction du cylindre creux :

𝑅2
𝐿𝑛
𝑅1
𝑅=
2𝜋 𝐿

Remarque :

Dans le cas d’un cylindre creux mince d’épaisseur e.

𝑅2 = 𝑅1 + 𝑒

𝑅2 𝑒 𝑒
𝐿𝑛 = 𝐿𝑛 1 + ~
𝑅1 𝑅1 𝑅1

𝑇1 − 𝑇2 1  𝑅1
𝜑𝑟 =  𝑒 = 𝑇1 − 𝑇2
𝑟 𝑒 𝑟
𝑅1

2𝜋𝐿𝑅1
∅ = 𝑇1 − 𝑇2
𝑒

4.4. Conduction de la chaleur dans un tube circulaire composite

Soit un cylindre composé de n matériaux superposés par les

cylindres de rayons 𝑅1 , 𝑅2 , 𝑅3 , … … . . , 𝑅𝑛 et de conductivités
thermiques respectives 1 , 2 , 3 , … … . . 𝑛 . Les résistances
thermiques de ces cylindres sont :

𝑅 𝑅 𝑅
𝐿𝑛 𝑅1 𝐿𝑛 𝑅2 𝐿𝑛 𝑅 𝑛
0 2 𝑛−1
𝑅 1 = , 𝑅 2 = , … … … … 𝑅 𝑛 =
2𝜋 1 𝐿 2𝜋 2 𝐿 2𝜋 1𝑛 𝐿

La résistance thermique du tube est donc :

𝑅
𝑛 𝐿𝑛 𝑅 𝑖
𝑖−1
𝑅=
2𝜋 𝑖 𝐿
𝑖=1

22
Si 𝑇1 𝑒𝑡 𝑇2 sont imposées, on peut calculer comme précédemment le flux thermique par la relation :

𝑇2 − 𝑇1
∅ =
𝑅

4.5. La sphère creuse

C’est un problème au problème du cylindre creux (Cf 4.3). On considère une sphère creuse de
rayon intérieur 𝑅1 et de rayon extérieur 𝑅2 , constitué d’un matériau homogène et isotrope dont les
surfaces cylindriques sont à des températures uniformes 𝑇1 𝑒𝑡 𝑇2 . Les isothermes sont des
surfaces sphériques concentriques, le gradient de température est radial. L’équation du champ
thermique s’écrit :

𝑑2 𝑇 2 𝑑𝑇
+ =0
𝑑𝑟 2 𝑟 𝑑𝑟

Les conditions aux limites sont :

𝑇 𝑟 = 𝑅1 = 𝑇1

𝑇 𝑟 = 𝑅2 = 𝑇2

La résolution de cette équation donne :

1 1
−
𝑟 𝑅1
𝑇 𝑟 = 𝑇1 + 𝑇1 − 𝑇2
1 1
𝑅1 − 𝑅2

On déduit de la loi de Fourrier, la densité de flux :

𝑇1 − 𝑇2 1
𝜑 = − 𝑔𝑟𝑎𝑑 𝑇 =  𝑒
1 1 𝑟2 𝑟
−
𝑅1 𝑅2

Le flux thermique ∅ traversant chaque surface cylindrique de longueur L est :

𝑇1 − 𝑇2 1
∅= × 4 𝜋𝑟 2
1 1 𝑟2
𝑅1 − 𝑅2

𝑅1 𝑅2
∅ = 4 𝜋 𝑇1 − 𝑇2
𝑅2 − 𝑅1

D’où la résistance thermique de conduction de la sphère creuse :

1 𝑅1 𝑅2
𝑅=
4𝜋 𝑅2 − 𝑅1

4.6. Epaisseur critique d’un revêtement isolant

Supposons un tube cylindrique composite ou non de longueur L et de rayon intérieur 𝑅1 et de

rayon extérieur 𝑅2 dont la résistance thermique est 𝑅 0 . Soit 𝑇1 la température de la face intérieure
du tube. Supposons qu’autour de ce tube, soit placé un isolant de rayon extérieur 𝑟 et de
conductivité thermique . Soit 𝑕 le coefficient de convection avec l’air ambiant de température 𝑇𝑎.
23
La résistance thermique du système est alors :
𝑟
𝐿𝑛 𝑅 1
2
𝑅 = 𝑅0 + +
2𝜋 𝐿 2𝜋𝑕𝑟 𝐿

∆𝑇
∅=
𝑅

∆𝑇
∅= 𝑟
𝐿𝑛 𝑅 1
2
𝑅0 + +
2𝜋 𝐿 2𝜋𝑕𝑟 𝐿

On voit que si 𝑟 augmente, les deux termes du dénominateur varient en sens contraire. Le
1
premier à une augmentation de ∅ due à un accroissement de la surface d’échange par
2𝜋𝑕𝑟 𝐿
𝑟
𝐿𝑛
𝑅2
convection. Tandis que le second correspond à une augmentation de la résistance
2𝜋  𝐿

thermique de la couche de calorifuge avec r. Il est donc possible qu’étant donné des effets
opposés qui doivent s’ajouter, il existe une valeur de 𝑟 pour laquelle ∅ soit un extrémum.

𝑑∅ 𝑑𝑅
= 0 𝑜𝑢 =0
𝑑𝑟 𝑑𝑟

On vérifie que ces deux relations sont vraies pour


𝑟 = 𝑟𝑐 =
𝑕

𝑟𝑐 est le rayon critique qui ne dépend que de  et 𝑕.

𝑑2∅
𝑑𝑟 2 𝑟= 𝑟𝑐
est négative, donc l’extrémum est un maximum

A partir de cette courbe, il faut distinguer deux cas :

∎ 𝑟 < 𝑟𝑐

Au fur et à mesure qu’on ajoute un calorifuge sur le fil de rayon 𝑟, la déperdition ∅ commence à
augmenter jusqu’à ce que le rayon extérieur 𝑟 atteint le rayon critique 𝑟𝑐 puis si l’on continue ∅
décroit.

On profite de ce fait pour augmenter au maximum la déperdition des petits fils électriques, en
tentant de se rapprocher de 𝑟 = 𝑟𝑐 .

∎ 𝑟 > 𝑟𝑐 (On augmente l’épaisseur de l’isolant)

24
Au fur et mesure que l’on applique le calorifuge, la déperdition ∅ diminue. Dans la plupart des cas
d’isolation des tuyaux, on trouve 𝑟 > 𝑟𝑐 , car 𝑟 est relativement grand et  est relativement petit
(isolants).

4.7. Conduction de la chaleur dans un tube creux avec production de la chaleur

La loi de propagation de la chaleur devient alors :

𝑑2 𝑇 1 𝑑𝑇 𝑝
+ = −
𝑑𝑟 2 𝑟 𝑑𝑟 
On suppose que la production 𝑝 est constante dans tout le volume.

𝑑2 𝑇 𝑑𝑇 𝑝𝑟
𝑟 + = −
𝑑𝑟 2 𝑑𝑟 
𝑑 𝑑𝑇 𝑝𝑟
𝑟 = −
𝑑𝑟 𝑑𝑟 
𝑑𝑇 𝑝 𝑟2
𝑟 = − + 𝐴
𝑑𝑟 2
𝑑𝑇 𝑝𝑟 𝐴
= − +
𝑑𝑟 2 𝑟

𝑝 𝑟2
𝑇= − + 𝐴 𝑙𝑛 𝑟 + 𝐵
4

𝐴 𝑒𝑡 𝐵 sont des constantes à déterminer par les conditions aux limites.

a) Températures imposées sur les surfaces.

𝑇 𝑟 = 𝑅1 = 𝑇1

𝑇 𝑟 = 𝑅2 = 𝑇2

𝑇2 − 𝑇1 𝑝 𝑅22 − 𝑅12
𝐴= +
𝑅
𝑙𝑛 𝑅2 4  𝑙𝑛 𝑅2
1 𝑅1

𝑇1 𝑙𝑛 𝑅2 − 𝑇2 𝑙𝑛 𝑅1 𝑝 𝑅12 𝑙𝑛 𝑅2 − 𝑅22 𝑙𝑛 𝑅1
𝐵= +
𝑅
𝑙𝑛 𝑅2 4 𝑅
𝑙𝑛 𝑅2
1 1

La densité de flux est :

𝑑𝑇 𝑝𝑟 𝐴
𝜑=− = −
𝑑𝑟 2 𝑟
𝑝𝑟  𝑇2 − 𝑇1 𝑝 𝑅22 − 𝑅12
𝜑= − +
2 𝑟 𝑅
𝑙𝑛 𝑅2 4  𝑙𝑛 𝑅2
1 𝑅1
b) Température imposée sur une face et un flux sur l’autre.

25
 𝑇 𝑟 = 𝑅1 = 𝑇1
𝜑 𝑟 = 𝑅2 = 𝜑2
On obtient :
𝑝 𝑅2  𝐴
𝜑2 = −
2 𝑅2

D’où :
𝑝 𝑅22 𝜑2 𝑅2
𝐴= −
2 
𝑝 𝑅12
𝑇1 = − + 𝐴 𝑙𝑛 𝑅1 + 𝐵
4
𝑝 𝑅12 𝑝 𝑅22 𝜑2 𝑅2
𝐵 = 𝑇1 + − − 𝑙𝑛 𝑅1
4 2 
𝑝 𝑅22 𝜑2 𝑅2 𝑟
𝑇 = 𝑇1 + − 𝑙𝑛
2  𝑅1

4.8. Conduction de la chaleur dans une barre avec production de chaleur.

C’est le cas précédent où 𝑅1 = 0, on a alors 𝐴 = 0

𝑝 𝑟2
𝑇= − +𝐵
4
𝑑𝑇 𝑝𝑟
𝜑 = − =
𝑑𝑟 2

La loi de distribution de la température est parabolique, le maximum étant au centre de la barre.

Les conditions aux limites ne peuvent être que la température sur la face extérieure, c'est-à-dire :

𝑇 𝑟 = 𝑅2 = 𝑇2
𝑝
𝑇 = 𝑇2 + 𝑅2 − 𝑟2
4 2
La température au centre est :
𝑝 2
𝑇0 = 𝑇2 + 𝑅
4 2
D’où l’écart de température :
𝑝 2
𝛿𝑇 = 𝑇0 − 𝑇2 = 𝑅
4 2
Le flux de chaleur sortant du cylindre est :

𝑝 𝑅2
 = 2𝜋 𝑅2 𝐿 = 𝜋 𝑅22 𝐿 𝑝 = 𝑝 𝑉
2

𝑉 est le volume de la barre. Le résultat est logique.

5 - Le problème de la tige mince de section constante

Ce problème présente deux aspects intéressants :

26
 - La tige constitue une modélisation des ailettes ou d’autres dispositifs utilisés pour accroitre
la surface effective d’échange d’un appareil de transmission calorifique sans augmenter
son encombrement.
- Il s’agit d’un exemple de problème où une mise en équation directe tenant compte
d’hypothèses simplificatrices.

A – Mise en équation

Considérons une tige cylindrique ayant une section droite implantée sur une surface plane.

On suppose que la conductibilité thermique ( = 𝑐𝑡𝑒) est constante et que la production p est
nulle (𝑝 = 0).

Contrairement au cas précédent, où on a résolu le système général constitué de l’équation locale

et les conditions aux limites ; ici on va utiliser une formulation globale incluant directement les
conditions de surface dans l’équation.

On va écrire un bilan sur un élément de volume de longueur 𝑑𝑥 couvrant la section totale S de la

barre. La section est suffisamment faible de telle sorte qu’on peut considérer que 𝑇 = 𝑇 𝑥 .

Bilan :
- Flux entrant par le plan situé à l’abscisse 𝑥 est :
𝑑𝑇(𝑥)
− 𝑆
𝑑𝑥
- Flux sortant par le plan situé à l’abscisse 𝑥 + 𝑑𝑥 est :

𝑑𝑇(𝑥 + 𝑑𝑥) 𝑑 𝑑𝑇(𝑥)

− 𝑆 = − 𝑇+ 𝑑𝑥 𝑆
𝑑𝑥 𝑑𝑥 𝑑𝑥

- Flux sortant la tranche de surface 𝑝 𝑑𝑥 est : 𝑕 𝑝𝑑𝑥 𝑇 𝑥 − 𝑇∝

Où 𝑇∝ est la température supposée uniforme du fluide.

𝑝 est le périmètre de la tige.

A l’équilibre (en régime permanent) le flux entrant est égal aux flux sortant.

𝑑𝑇(𝑥) 𝑑 𝑑𝑇(𝑥)
− 𝑆 = − 𝑇+ 𝑑𝑥 𝑆 + 𝑕 𝑝𝑑𝑥 𝑇 𝑥 − 𝑇∝
𝑑𝑥 𝑑𝑥 𝑑𝑥

𝑑2 𝑇(𝑥)
 𝑆 = 𝑕 𝑝 𝑇 𝑥 − 𝑇∝
𝑑𝑥 2

27
C’est une équation différentielle du second ordre à coefficients constants qu’on peut écrire sous la
forme :

𝑑2 𝑇(𝑥)
− 𝛼 2 𝑇 𝑥 − 𝑇∝ = 0
𝑑𝑥 2

Avec

𝑕𝑝
𝛼2 =
𝑆
La solution générale de l’équation est de la forme :

𝑇 𝑥 − 𝑇∞ = 𝐴 𝑒 −𝛼𝑥 + 𝐵 𝑒 𝛼𝑥 𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝐴

Ou
𝑇 𝑥 − 𝑇∞ = 𝐴1 𝐶𝑕 𝛼𝑥 + 𝐵1 𝑆𝑕 𝛼𝑥 𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝐵

Les constantes A et B sont déterminées par les conditions aux limites (les conditions aux
extrémités).

L’équation A est bien adaptée aux problèmes relatifs aux barres très longues. La deuxième B aux
géométries limitées.

B – Cas d’une barre très longue.

A l’extrémité non chauffée d’une barre de très grande longueur 𝑥 → ∞ , la température tend vers
celle du milieu ambiant 𝑇∞ .

𝑥 𝑔𝑟𝑎𝑛𝑑 , > 0 , 𝑒 𝛼𝑥 → ∞ or T doit rester finie  B=0.

A l’extrémité chauffée 𝑥 = 0 , 𝑇 = 𝑇0  𝑇0 − 𝑇∞ = 𝐴
D’où :
𝑇 𝑥 − 𝑇∞
= 𝑒 −𝛼𝑥
𝑇0 − 𝑇∞
𝛼 = 𝐿−1 Caractérise la rapidité avec laquelle l’échauffement s’estompe le long de la barre pour
des abscisses croissantes.

Le flux traversant une section à l’abscisse 𝑥 est donné par la loi de Fourrier :

𝑑𝑇(𝑥)
∅= − 𝑆 = 𝛼  𝑆 𝑇 𝑥 − 𝑇∞
𝑑𝑥

Le flux de dissipation sur toute la surface latérale de la barre peut être calculé par intégration du
flux de convection local :
∝ 1
∅𝑝𝑒𝑟𝑡𝑒𝑠 = 0
𝑕 𝑝 𝑇 𝑥 − 𝑇∝ 𝑑𝑥 = 𝑕 𝑝  𝑆 2 𝑇0 − 𝑇∝

Ce flux (en régime permanent) est égal au flux transféré par conduction dans la section de la barre
à l’abscisse 𝑥 = 0.
1
∅ 𝑥 = 0 = 𝛼  𝑆 𝑇 0 − 𝑇∝ = 𝑕 𝑝  𝑆 2 𝑇0 − 𝑇∝ .

28
C – Cas d’une barre dont l’extrémité non chauffée est adiabatique.

Pour une barre de longueur L, si la puissance traversant l’extrémité non chauffée est nulle ou
négligeable (extrémité isolée) les constantes 𝐴1 et 𝐵1 de la solution générale (équation B) sont
données par :
𝑇 0 = 𝑇0  𝑇0 − 𝑇∞ = 𝐴1
Et

𝑑𝑇(𝑥)
∅ 𝐿 = − 𝑆 = 0  𝐴1 𝛼𝑆𝑕 𝛼𝐿 + 𝐵1 𝛼𝐶𝑕 𝛼𝐿 = 0
𝑑𝑥 𝑥=𝐿
𝐵1
= 𝑡𝑕 𝛼𝐿
𝐴1
D’où :
𝑇 𝑥 − 𝑇∞ 𝐶𝑕𝛼 𝑥 − 𝐿
=
𝑇0 − 𝑇∞ 𝐶𝑕 𝛼𝐿

Le flux dissipé s’évalue comme précédemment à partir du flux conductif à la base de la barre :

𝑑𝑇(𝑥)
∅ 𝑥 =0 = − 𝑆
𝑑𝑥 𝑥=0

∅ 𝑥 = 0 = 𝛼  𝑆 𝑇0 − 𝑇∞ 𝑡𝑕 𝛼𝐿
1
∅ 𝑥 =0 = 𝑕𝑝𝑆 2 𝑇0 − 𝑇∞ 𝑡𝑕 𝛼𝐿

D – Cas d’une barre dont l’extrémité est à une température constante

Pour une barre de longueur 𝐿 dont l’extrémité 𝑥 = 𝐿 est à la température 𝑇1 , les constantes 𝐴1 et
𝐵1 de la solution générale (équation B) sont données par :
𝑇 0 = 𝑇0  𝑇0 − 𝑇∞ = 𝐴1
𝑇 𝐿 = 𝑇1
𝑇 𝐿 − 𝑇∞ = 𝑇1 − 𝑇∞ = 𝐴1 𝑐𝑕 𝛼𝐿 + 𝐵1 𝑠𝑕 𝛼𝐿

𝑇1 − 𝑇∞ − 𝑇0 − 𝑇∞ 𝑐𝑕 𝛼𝐿
𝐵1 =
𝑠𝑕 𝛼𝐿

D’où :

𝑇1 − 𝑇∞ − 𝑇0 − 𝑇∞ 𝑐𝑕 𝛼𝐿
𝑇 𝑥 − 𝑇∞ = 𝑇0 − 𝑇∞ 𝑐𝑕 𝛼𝑥 + 𝑠𝑕 𝛼𝑥
𝑠𝑕 𝛼𝐿

𝑇0 − 𝑇∞ 𝑐𝑕 𝛼𝑥 𝑠𝑕 𝛼𝐿 + 𝑇1 − 𝑇∞ − 𝑇0 − 𝑇∞ 𝑐𝑕 𝛼𝐿
𝑇 𝑥 − 𝑇∞ = 𝑠𝑕 𝛼𝑥
𝑠𝑕 𝛼𝐿

𝑇0 − 𝑇∞ 𝑐𝑕 𝛼𝑥 𝑠𝑕 𝛼𝐿 − 𝑐𝑕 𝛼𝐿 𝑠𝑕 𝛼𝑥 + 𝑇1 − 𝑇∞ 𝑠𝑕 𝛼𝑥
𝑇 𝑥 − 𝑇∞ =
𝑠𝑕 𝛼𝐿

𝑇0 − 𝑇∞ 𝑠𝑕𝛼 𝐿 − 𝑥 + 𝑇1 − 𝑇∞ 𝑠𝑕 𝛼𝑥
𝑇 𝑥 − 𝑇∞ =
𝑠𝑕 𝛼𝐿
29
Le flux dissipé s’évalue comme précédemment à partir du flux conductif à la base de la barre :

𝑑𝑇(𝑥)
=∅ 𝑥 =0 = − 𝑆
𝑑𝑥 𝑥=0

𝑇1 − 𝑇∞
𝑐𝑕 𝛼𝐿 −
1 𝑇0 − 𝑇∞
= 𝑕𝑝𝑆 2 𝑇0 − 𝑇∞
𝑠𝑕 𝛼𝐿

E – Cas d’une barre avec échange convectif à l’extrémité non chauffée.

Lorsque la barre est courte, tout le flux n’est pas dissipé par sa surface latérale et sa section droite
extrême est encore traversée par un flux convectif. On a donc :

𝑇 0 = 𝑇0  𝑇0 − 𝑇∞ = 𝐴1

𝑑𝑇(𝑥)
− 𝑆 = 𝑕 𝑆 𝑇 𝐿 − 𝑇∝
𝑑𝑥 𝑥=𝐿
𝑕
𝑠𝑕 𝛼𝐿 + 𝐶𝑕 𝛼𝐿
𝐵1 = − 𝑇0 − 𝑇∞ 𝛼
𝑕
𝑐𝑕 𝛼𝐿 + 𝑠𝑕 𝛼𝐿
𝛼
D’où la solution :
𝑕
𝐶𝑕𝛼 𝐿 − 𝑥 +
𝑇 𝑥 − 𝑇∞
= 𝛼 𝑠𝑕𝛼 𝐿 − 𝑥
𝑇0 − 𝑇∞ 𝑕
𝐶𝑕 𝛼𝐿 + 𝑠𝑕 𝛼𝐿
𝛼
Le flux total de dissipation est :
𝑑𝑇(𝑥)
∅𝑝𝑒𝑟𝑡𝑒𝑠 = ∅ 𝑥 = 0 = −  𝑆
𝑑𝑥 𝑥=0

𝑕
𝑡𝑕𝛼𝐿 +
∅𝑝𝑒𝑟𝑡𝑒𝑠 = 𝛼  𝑆 𝑇0 − 𝑇∞ 𝛼
𝑕
1+ 𝑡𝑕𝛼𝐿
𝛼
La solution donnant le flux total de dissipation pour une
tige avec de la convection à l’extrémité (cas 2) est
compliquée comparée à la solution donnant le flux total
de dissipation pour une tige avec un flux nul à l’extrémité
(cas 1). Une solution de compromis à été proposé par
Harper-Brown qui consiste à utiliser à utiliser la solution
du 𝑐𝑎𝑠 1 pour le 𝑐𝑎𝑠 2 à condition de prendre une
longueur 𝐿𝐶 telle que :

𝑆
𝐿𝐶 = 𝐿 +
𝑝
𝑒
𝐿𝐶 = 𝐿 + 𝑝𝑜𝑢𝑟 𝑢𝑛𝑒 𝑡𝑖𝑔𝑒 𝑟𝑒𝑐𝑡𝑎𝑛𝑔𝑢𝑙𝑎𝑖𝑟𝑒
2
𝐷
𝐿𝐶 = 𝐿 + 𝑝𝑜𝑢𝑟 𝑢𝑛𝑒 𝑡𝑖𝑔𝑒 𝑐𝑦𝑙𝑖𝑛𝑑𝑟𝑖𝑞𝑢𝑒
4

30
𝑒 𝑒𝑠𝑡 𝑙 ′ é𝑝𝑎𝑖𝑠𝑠𝑒𝑢𝑟 𝑑𝑒 𝑙𝑎 𝑡𝑖𝑔𝑒 𝑟𝑒𝑐𝑡𝑎𝑛𝑔𝑢𝑙𝑎𝑖𝑟𝑒
𝐷 𝑒𝑠𝑡 𝑙𝑒 𝑑𝑖𝑎𝑚è𝑡𝑟𝑒 𝑑𝑒 𝑙𝑎 𝑡𝑖𝑔𝑒 𝑐𝑦𝑙𝑖𝑛𝑑𝑟𝑖𝑞𝑢𝑒
La correction de Harper-Brown ne s’applique pas au calcul de la température mais uniquement à
celui du flux.

6 - Ailettes de refroidissement

En vue d’accroître l’échange de chaleur entre une surface et un fluide, il est d’usage courant de
garnir cette surface d’extensions. Si la surface d’échange est ainsi majorée, l’écart moyen de
température entre la paroi et le fluide se trouve réduit. Une tige mince de section rectangulaire
dont une dimension est prépondérante est le modèle d’une ailette de refroidissement.

Hypothèses simplificatrices :
- Régime permanent
- La longueur de l’ailette est grande par rapport à son épaisseur.
- L’ailette est constituée d’un matériau homogène et isotrope. La conductivité thermique est
constante.
- La température du fluide ambiant est uniforme.
- Il n’y a pas de résistance de contact entre l’ailette et la surface mère.
- Il n’y a pas de source interne de chaleur dans l’ailette.

6.a. Cas des ailettes d’épaisseur constante

Considérons une ailette de longueur L, soumise à la convection sur ses surfaces latérales et sur
son extrémité. L’évolution de la température dans l’ailette s’écrit :

𝑕
𝐶𝑕𝛼 𝐿 − 𝑥 +
𝑇 𝑥 − 𝑇∞
= 𝛼 𝑠𝑕𝛼 𝐿 − 𝑥
𝑇0 − 𝑇∞ 𝑕
𝐶𝑕 𝛼𝐿 + 𝑠𝑕 𝛼𝐿
𝛼

Avec :
𝑕𝑝
𝛼2 =
𝑆
Ou 𝑝 et 𝑆 représentent le périmètre et la surface d’une coupe transversale de l’ailette.

Le flux calorifique au pied de l’ailette à pour expression (c’est le flux de chaleur dans la section
d’encastrement) :
𝑕
𝑡𝑕𝛼𝐿 +
∅0 = 𝛼  𝑆 𝑇0 − 𝑇∞ 𝛼
𝑕
1+ 𝑡𝑕𝛼𝐿
𝛼
En l’absence de l’ailette, le flux calorifique traversant la même surface S sous la différence de
température 𝑇0 − 𝑇∝ vaudrait :

𝑕𝑆 𝑇0 − 𝑇∞

Le recours à des ailettes sera donc avantageux si l’on a :

∅0 > 𝑕𝑆 𝑇0 − 𝑇∞

31
 𝑕
𝑡𝑕𝛼𝐿 +
𝛼  𝑆 𝑇0 − 𝑇∞ 𝛼 > 𝑕𝑆 𝑇0 − 𝑇∞
𝑕
1+ 𝑡𝑕𝛼𝐿
𝛼
𝑕𝑆
𝑝 <1

Pour une ailette de la forme suivante :

𝑝 = 4𝑏 + 2𝑥 ≅ 2𝑥

𝑆 = 2𝑏𝑥

𝑕𝑏
 <1
𝐵𝑖 < 1

𝐵𝑖 est le nombre de Biot.

On constate de la hauteur L est absente de cette condition. Par conséquent, selon que le nombre
de Biot est inférieur ou supérieur à 1, on a intérêt à mettre ou à ne pas mettre d’ailettes, quelque
soit la hauteur de celles-ci.

6.b. Dimensions optimales d’une ailette rectangulaire

On souhaite refroidir une plaque, avec des ailettes dont la masse,

onc le volume est donné.

Doit-on choisir la forme 1 ou la forme 2 ?

La dimension transversale étant donnée, on pose :

𝐵 = 2 𝑒 𝐿 = 𝑐𝑠𝑡𝑒

On suppose que la surface 𝑥 = 𝐿 est adiabatique. L’évolution de la température dans l’ailette

s’écrit :

𝑇 𝑥 − 𝑇∞ 𝐶𝑕𝛼 𝑥 − 𝐿
=
𝑇0 − 𝑇∞ 𝐶𝑕 𝛼𝐿

La puissance par unité de largeur 𝑙 = 1 évacuée par l’ailette est :

1 𝑕𝐿2
∅ = 2𝑒 𝑕/ 𝑒 2 𝑇0 − 𝑇∞ 𝑡𝑕
𝑒
Ou

1 𝑕𝐵2
∅= 𝐾𝑒 2 𝑡𝑕
4 𝑒 3

Il faut à volume donné, évacuer une puissance maximale :

Si on pose :

32
 𝑕𝐵2
𝑢= = 𝛼𝐿
4 𝑒 3

𝜕∅
=0
𝜕𝑢

𝜕 −1
 𝑢 3 𝑡𝑕 𝑢 = 0
𝜕𝑢
Il faut résoudre graphiquement l’équation :

𝑡𝑕 𝑢 = 3𝑢 1 − 𝑡𝑕2 𝑢

 𝑢 = 1.419

Cette valeur de 𝑢 correspond à

𝑡𝑕 𝛼𝐿
= 0.626
𝛼𝐿

Le couple 𝑒, 𝐿 se déduit des relations :

𝑕𝐿2
𝐵 =2𝑒𝐿 𝑒𝑡 =2
𝑒

6.c. Efficacité d’une ailette

Par définition, l’efficacité 𝜀 d‘une ailette est la rapport entre le flux calorifique ∅𝑟é𝑒𝑙 transféré par
l’ailette au fluide et le flux ∅𝑚𝑎𝑥 qui serait transféré par la même ailette si elle était isotherme (ce
qui revient au même, si elle était constitué d’un matériau de conductibilité thermique infinie).

∅𝑚𝑎𝑥 = 𝑕 𝑝 𝐿 𝑇0 − 𝑇∞

 Pour l’ailette dont le flux à l’extrémité est nul

1
𝑕 𝑝  𝑆 2 𝑇0 − 𝑇∝ 𝑡𝑕 𝛼𝐿 𝑡𝑕 𝛼𝐿
𝜖= =
𝑕 𝑝 𝐿 𝑇0 − 𝑇∝ 𝛼𝐿

Cette relation montre que si la longueur de l’ailette croit, l’efficacité diminue. 𝛼𝐿 fait intervenir la
longueur et les propriétés de l’ailette. A propriétés égales, plus l’ailette est longue, moins elle est
efficace. Dans la pratique, la on choisit des ailettes efficaces (𝜖 de 0.8 à 0.9).

 Lorsque le flux n’est pas nul à l’extrémité

𝑕
𝑡𝑕𝛼𝐿 +
𝛼  𝑆 𝑇0 − 𝑇∝ 𝛼
𝑕 𝑕
1+ 𝑡𝑕𝛼𝐿 𝑡𝑕𝛼𝐿 +
𝜖= 𝛼  = 𝛼
𝑕 𝑝 𝐿 𝑇0 − 𝑇∝ 𝐿𝑕
𝛼𝐿 + 𝑡𝑕𝛼𝐿

L’expression s’emploi souvent sous la forme ∅𝑟é𝑒𝑙 = 𝜖 ∅𝑚𝑎𝑥 où 𝜖 est souvent déterminé par des
graphiques (figure 2 et figure 3). La relation montre bien que le flux réel est le produit du flux
maximal par une efficacité.

33
 Fig. 2. Efficacité des ailettes 2D [1]

Fig. 3. Efficacité des ailettes 2D [1]

6.d. Efficacité des surfaces ailettées

La notion d’efficacité dune ailette peut être étendue à une surface ailettée. Considérons une
surface de base 𝑆𝑏 (aire initiale) garnie d’ailettes qui procurent une surface d’échange
additionnelle. On adjoint à la surface 𝑆𝑏 un certain nombres d’ailettes de surface 𝑆𝑎 ).

On définit l’efficacité  de la surface d’échange ailettée par le rapport du flux effectivement

échangé par la surface totale 𝑆𝑇 (𝑆𝑇 = 𝑛𝑆𝑎 + 𝑆𝑏 ) du corps garni de ses ailettes au flux maximum
que pourrait échanger cette même surface si tout le corps était à une température uniforme égale
à la température de base de l’ailette.

L’efficacité d’une surface ailettée  est définie comme le rapport du flux échangé sur toute la
surface, sur le flux qui aurait été échangé si la totalité de la surface était à la température 𝑇0 .

34
 ∅ 𝑇𝑜𝑡𝑎𝑙
 =
𝑕 𝑆𝑇 𝑇0 − 𝑇∞

𝑆𝑇 : Surface d’échange totale (totalité des ailettes+base)

∅ 𝑇𝑜𝑡𝑎𝑙 : Flux thermique au travers de la surface totale 𝑆𝑇 .
Le flux total s’écrit :
∅ 𝑇𝑜𝑡𝑎𝑙 = 𝑕 𝑆𝑏 𝑇0 − 𝑇∞ + 𝑁 𝜖 𝑕 𝑆𝑎 𝑇0 − 𝑇∞

𝑁 : Nombre d’ailettes sur la surface

𝑆𝑎 : Surface d’échange d’une ailette
𝑆𝑏 : Surface de base à la température 𝑇0

𝑆𝑇 = 𝑛 𝑆𝑎 + 𝑆𝑏

𝑕 𝑆𝑏 𝑇0 − 𝑇∞ + 𝑁 𝜖 𝑕 𝑆𝑎 𝑇0 − 𝑇∞ 𝑆𝑏 + 𝑁 𝜖 𝑆𝑎 𝑆𝑇 − 𝑛𝑆𝑎 + 𝑁 𝜖 𝑆𝑎
 = = =
𝑕 𝑆𝑇 𝑇0 − 𝑇∞ 𝑆𝑇 𝑆𝑇

On obtient l’expression de l’efficacité d’une surface ailettée :

𝑁 𝑆𝑎
 = 1− 1−𝜖 𝑆𝑇

6.e. Ailettes triangulaire

Hypothèses :

- Régime permanent
- 𝐿, 𝑙 ≫ 𝑏 𝐿
- On inverse l’axe 𝑂𝑥
- 𝑙 : largeur de l’ailette

Bilan :
𝑏 𝐿
𝑏 𝑥 =𝑥
𝐿

𝑕𝑝 𝑕
𝛼2 = =
𝑆 𝑏
𝑑 𝑑𝑇
𝑆 𝑥 =𝑕𝑝 𝑥 𝑇 𝑥 − 𝑇𝑒
𝑑𝑥 𝑑𝑥

𝑆 𝑥 =2𝑏 𝑥 𝑙

𝑝 𝑥 = 2 2 𝑏 𝑥 + 𝑙 ≅ 2𝑙

L’équation est la suivante :

𝑑2 𝜃 𝑑𝜃
𝑥 + − 𝐿𝛼 2 𝜃 = 0
𝑑𝑥 2 𝑑𝑥
35
La résolution de cette équation se fait par les fonctions de Bessel.

7 - Conduction en régime instationnaire

Dans cette partie, nous étudierons les échanges thermiques par conduction dans le cas ou la
température est fonction non seulement des conditions spatiales ( en particulier x) mais également
du temps et ce, en chaque point du solide.

On distingue parmi les problèmes instationnaires, les problèmes :

- Périodiques : caractérisés par la répétition d’un cycle sinusoïdal ou quelconques (variations

diurnes, moteurs à combustion interne).
- Apériodiques : Mise à température de régime des moteurs.
Démarrage ou arrêt d’installation thermique, de réacteurs nucléaires,
…..etc….

L’étude des phénomènes instationnaires est extrêmement importante dans la pratique et il est
difficile de trouver une activité humaine ou il n’intervient pas.

Il convient d’ajouter que un grand nombre de problèmes thermiques, le régime permanant n’est
jamais atteint ; d’ailleurs théoriquement la durée d’établissement d’un régime permanent est infinie.
En pratique, on considère que le régime est établi quand la température en chaque point du
système ne diffère de la température du régime que d’une certaine valeur préconvenue choisie par
exemple d’après le degré de précision des appareils de mesures utilisées.

Une autre raison qui rend importante l’étude du régime instationnaire est que les gradients de
températures sont beaucoup plus importants, en régime transitoire qu’en régime permanents, il
s’ensuit des dilatations, des tensions internes d’origines thermiques « contraintes thermiques » qui
peuvent être beaucoup plus importantes que celles qui existent en régime permanent.

Le but de notre étude est de déterminer la distribution de température en fonction du temps.

Nous commencerons par l’étude des problèmes où la résistance interne par conduction est
𝐿
négligeable (grande conductivité thermique et petite dimension caractéristique L) pour
𝑆
lesquels la température peut être considérée comme uniforme.

Ensuite, nous étudierons les problèmes où au contraire, la résistance de surface par convection
1
𝑕𝑆
est négligeable ( h =∞) et pour lesquels on impose instantanément en surface, une
température donnée.

L’expérience montre qu’on peut classer les corps du point de vue leur comportement thermique en
corps thermiquement minces et corps thermiquement épais.

Pour les corps thermiquement minces, la température peut y être considérée comme uniforme.
Elle ne dépend que du temps. Tandis que pour les systèmes épais, la température est fonction du
point considéré et du temps.

Nous avons vu que :

𝜕2 𝑇 𝜕2 𝑇 𝜕2 𝑇 1 𝜕𝑇
2
+ 2
+ 2=
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝑎 𝜕𝑡

36
Pour des raisons de simplification, on utilise les méthodes de résolution de cette équation dans le
cas suivant :

𝜕2 𝑇 1 𝜕𝑇
2
=
𝜕𝑥 𝑎 𝜕𝑡

7. 1. Refroidissement ou chauffage d’un solide très conducteur.

On considère un corps de forme quelconque brusquement plongé

dans un fluide de température 𝑇∞ . La température initiale du corps
est 𝑇0 .

Pour que la chaleur se propage dans le milieu, il faut qu’il n’y ait
pas de gradient de température.

En effet,

𝜕𝑇 𝑇𝑆 − 𝑇𝐶
≪
𝜕𝑙 𝐿

TC est la température au centre du solide.

𝑇𝑆 − 𝑇𝐶 ≪ 𝑇𝑆 − 𝑇𝑒

Avec la condition aux limites :

𝜕𝑇
− ≪ 𝑕 𝑇𝑆 − 𝑇𝐶
𝜕𝑙

on a :
𝑇𝑆 − 𝑇𝐶
− = 𝑕 𝑇𝑆 − 𝑇𝑒
𝐿

Donc :
𝑇𝑆 − 𝑇𝐶 𝑕𝐿
= = 𝐵𝑖 ≪ 1
𝑇𝑆 − 𝑇𝑒 

On admet que la conductivité thermique est très grande de telle sorte que l’on puisse négliger les
variations de température du corps.

Le problème est entièrement contrôlé par la résistance de convection entre le fluide et la paroi.

On mesure l’importance relative de ces deux résistances par le nombre de Biot Bi qui est défini
quelque soit la forme géométrique du solide étudié par la relation :

𝑕𝐿
𝐵𝑖 =

𝐿 : dimension caractéristique, obtenue en divisant le volume du solide par la surface d’échange
offerte à l’ambiance.
𝑙
𝐵𝑖 = 
1
𝑕
Un nombre de Biot petit signifie que la résistance interne est petite par rapport à la résistance
externe.
37
Longueur caractéristique :

𝑉𝑜𝑙𝑢𝑚𝑒 𝑑𝑢 𝑠𝑜𝑙𝑖𝑑𝑒
𝐿=
𝑆𝑢𝑟𝑓𝑎𝑐𝑒 𝑑′é𝑐𝑕𝑎𝑛𝑔𝑒

 Pour un mur d’épaisseur 2 échangeant de la chaleur sur les deux faces, chacune d’aire A on a :

2 𝐴
𝐿= = 
2𝐴
La dimension caractéristique est égale à la moitié de son épaisseur.

 Pour un mur d’épaisseur  échangeant de la chaleur sur une seule face d’aire A et isolé sur
l’autre on a :

𝐴
𝐿= = 
𝐴
La dimension caractéristique est égale à son épaisseur.

 Pour un cylindre infini de rayon R on a :

 𝑅2 𝐿 𝑅
𝐿= =
2 𝑅 𝐿 2

 Pour une sphère de rayon R on a :

4
3
 𝑅3 𝑅
𝐿= =
4 𝑅 2 3

Plusieurs auteurs définissent le nombre de Biot du cylindre et de la sphère comme (Attention il y

a confusion ?!) :
𝑕𝑅
𝐵𝑖 =

 Pour un parallélépipède on calculera 𝐵𝑖𝑥 , 𝐵𝑖𝑦 , 𝐵𝑖𝑧 .

Il faut donc calculer le nombre de Biot avant de se lancer dans les calculs compliqués.

L’équation caractéristique du problème s’obtient en écrivant que la quantité de chaleur perdue par
le corps est égale à celle évacuée dans le fluide.

− 𝑕 𝑆 𝑇 − 𝑇∞ 𝑑𝑡 = 𝑐 𝑉 𝑑𝑇

La quantité de chaleur transmise par convection pendant dt est égale à la dimunition de l’énergie
interne du corps. On a donc :
𝑑𝑇
𝑕 𝑆 𝑇 − 𝑇∞ = − 𝑐 𝑉
𝑑𝑡
Avec 𝑇 = 𝑇0 à t=0.

𝑕𝑆 𝑑𝑇
− 𝑑𝑡 =
𝑐 𝑉 𝑇 − 𝑇∞

38
𝑉 est le volume du corps, 𝑆 est la surface extérieure du corps.

Soit
 = 𝑇 − 𝑇∞

𝑕𝑆 𝑑
− 𝑑𝑡 = 𝑐 𝑉
𝑐 𝑉 

Avec  𝑡 = 0 = 𝑇0 − 𝑇∞ = 0

La solution de l’équation est :

 −
𝑕𝑆
𝑐 𝑉 𝑡
= 𝑒
0

En introduisant les paramètres adimensionnels :

𝑕𝐿 𝑕𝑉
𝐵𝑖 = =
 𝑆
𝑎𝑡 𝑡
𝐹𝑂 = =
𝐿2 𝑐𝐿2


= 𝑒 − 𝐵𝑖 𝐹𝑂
0

Notons que cette relation est indépendante de  puisque 𝐵𝑖 𝐹𝑂 est indépendant de  . Ceci montre
que  est quelconque.

La connaissance de 𝐵𝑖 𝑒𝑡 𝐹𝑂 permet de chartériser l’évolution de la température.

Signification des nombres adimensionnels

𝑉
𝑉 𝑆
𝑕𝐿 𝑕
𝐵𝑖 = = 𝑆 =  𝑆 = 𝑅é𝑠𝑖𝑠𝑡𝑎𝑛𝑐𝑒 𝑖𝑛𝑡𝑒𝑟𝑛𝑒
  1 𝑅é𝑠𝑖𝑠𝑡𝑎𝑛𝑐𝑒 𝑒𝑥𝑡𝑒𝑟𝑛𝑒
𝑕𝑆
∆𝑇 ∆𝑇
𝑎𝑡 𝑎𝑡  𝑙 𝑆  𝑆 𝑙 𝐹𝑙𝑢𝑥 𝑡𝑕𝑒𝑟𝑚𝑖𝑞𝑢𝑒 à 𝑡𝑟𝑎𝑣𝑒𝑟𝑠 𝑙𝑎 𝑠𝑢𝑟𝑓𝑎𝑐𝑒 𝑆
𝐹𝑂 = 2 = = = =
𝐿 𝑉 2 𝜌 𝐶 𝑉 ∆𝑇 𝑑 𝑚𝑐𝑇 𝐹𝑙𝑢𝑥 𝑑𝑒 𝑠𝑡𝑜𝑐𝑘𝑎𝑔𝑒 𝑑𝑎𝑛𝑠 𝑙𝑒 𝑣𝑜𝑙𝑢𝑚𝑒 𝑉
𝑆 𝑡 𝑑𝑡

Remarque

L’analogie électrique de ce problème est la décharge d’un condensateur sur un circuit purement
résistant

𝑉 𝑡
= 𝑒 − 𝑅𝐶
𝑉0
𝑐𝑉
La quotient 𝑕𝑆
a la dimension d’un temps et est appelé C « constante du temps du système ». Il
exprime la vitesse de réponse d’un système contenant une capacité simple, vis à vis d’un

39
changement brusque de la température ambiante. Plus C est grand , plus le système réagit
lentement.

Pour 𝑡 = 𝐶
𝑇 − 𝑇∞
= 𝑒 − 1 = 0,368
𝑇 − 𝑇0
Après une période de temps écoulée égale à la constante de temps C la différence de
température 𝑇 − 𝑇∞ tombe à 36,8 % de sa valeur initiale 𝑇 − 𝑇0 .

Exemple : Application à la trempe des métaux

Plaçons dans les conditions suivantes :

𝑕 = 4 104 𝑊 𝑚−2 𝐾 −1 (Convection avec ébullition)

 = 7 103 𝑘𝑔 𝑚−3 (métal)

𝐶 = 4 102 𝐽 𝑘𝑔−1 𝐾 −1

Prenons
𝑉
𝐿= = 0,1m comme dimension typique du corps immergé, on a :
𝑆

𝑐𝑉
𝐶 = = 10 𝑠
𝑕𝑠

Le temps caractéristique de changement de température est donc de 10 secondes.

On peut juger sur cet exemple l’importance du coefficient d’échange.

Si le métal était placé dans l’air ou 𝑕 = 40 𝑊 𝑚−2 𝐾 −1 (convection naturelle lente) le temps
caractéristique est :

𝑐𝑉
𝐶 = = 104 𝑠
𝑕𝑠
Conséquence :

Une conséquence très importante est relative au temps de réponse d’un thermocouple, en régime
instationnaire. Si la soudure du thermocouple est localement associée à un cylindre de rayon R.

𝑉 𝑅
𝐿= =
𝑆 2
𝑐𝑉 𝑐𝑅
𝐶 = =
𝑕𝑠 2𝑕

Il parait que le temps de réponse d’un thermocouple est d’autant plus petit que R est plus petit.
D’où l’intérêt de construire des thermocouples dont le diamètre externe soit aussi petit que
possible (quelques centièmes de mm) de façon à obtenir des temps de réponse très courts.

40
7. 2. Refroidissement ou chauffage d’une plaque plane.

On considère une plaque d’épaisseur e de conductivité 

isolée d’un côté et perdant de la chaleur par l’autre côté.

L’équation de la chaleur est :

𝜕2 𝑇 1 𝜕𝑇
2
=
𝜕𝑥 𝑎 𝜕𝑡
On pose  = 𝑇 − 𝑇∞

Les conditions aux limites sont :

𝜕𝑇
= 0
𝜕𝑥 𝑥=0
𝜕𝑇
− = 𝑕 𝑇 𝑒 − 𝑇∞
𝜕𝑥 𝑥=𝑒

 𝑡 = 0 = 0 = 𝑇0 − 𝑇∞

Ce problème a été résolu par Fourrier (1892) : Théorie analytique de la chaleur. C’est à ce propos
que Fourrier a développé l’étude des séries qui portent son nom.

La méthode de résolution de ces problèmes linéaires consistent à chercher une solution à

variables séparées.

 𝑥, 𝑡 = 𝐹 𝑡 𝐺 𝑥

L’équation de chaleur devient :

𝐹 ′ 𝑡 𝐺 𝑥 = 𝑎 𝐹 𝑡 𝐺" 𝑥
𝐹′ 𝑡 𝐺" 𝑥
= = ∓ 𝜇2
𝑎𝐹 𝑡 𝐺 𝑥

𝜇 est une constante appelée valeur propre ou constante de séparation.

𝐹 ′ 𝑡 = ∓ 𝜇2 𝑎 𝐹 𝑡

𝐺" 𝑥 = ∓ 𝜇2 𝐺 𝑥

La première équation a pour solution :

𝐹 𝑡 = 𝐶1 𝑒𝑥𝑝 ∓ 𝜇2 𝑎𝑡

Ce qui montre qu’il faut choisir le signe moins si nous voulons que  𝑥, 𝑡 décroisse dans le temps
(après un temps infini, on aura 𝑇 = 𝑇∞ ) (plaque plus chaude que le fluide). Donc :
𝐹 𝑡 = 𝐶1 𝑒𝑥𝑝 − 𝜇2 𝑎𝑡
𝐺 𝑥 = 𝐶2 𝑒𝑥𝑝 𝑖 𝜇 𝑥 + 𝐶3 𝑒𝑥𝑝 −𝑖 𝜇 𝑥

La solution générale s’écrit :

 𝑥, 𝑡 = 𝑒𝑥𝑝 − 𝜇2 𝑎𝑡 𝐴 𝑐𝑜𝑠 𝜇 𝑥 + 𝐵𝑠𝑖𝑛 𝜇 𝑥

Les constantes 𝐴 , 𝐵 et 𝜇 doivent être déterminées à l’aide des conditions aux limites.

𝜕
= 0  𝐵=0
𝜕𝑥 𝑥=0
41
La condition sur l’autre face :
𝜕 𝑕 
=− 𝑙  𝜇 = 𝑐𝑜𝑡𝑔 𝜇𝑙
𝜕𝑥 𝑥=𝑙  𝑕

Cette équation n’a pas une seule racine, elle en a toute une suite que nous appellerons 𝜇𝑛

(résolution graphique). La figure ci-dessous montre les intersections de la droite 𝑕
𝜇 avec la
courbe 𝑐𝑜𝑡𝑔 𝜇𝑙 .

La solution générale sera la somme de chacune des solutions, on a donc :

∞

 𝑥, 𝑡 = 𝐴𝑛 𝑒𝑥𝑝 − 𝜇𝑛2 𝑎𝑡 𝑐𝑜𝑠 𝜇𝑛 𝑥

𝑛=1

On ignore encore quelle est la constante 𝐴𝑛 qui correspond à chaque valeur propre 𝜇𝑛 , c’est la
condition initiale qui permet de le savoir.
∞

 𝑡 = 0 = 0 = 𝐴𝑛 𝑐𝑜𝑠 𝜇𝑛 𝑥
𝑛=1
𝑙 ∞ 𝑙
0 𝑐𝑜𝑠 𝜇𝑚 𝑥 𝑑𝑥 = 𝐴𝑛 𝑐𝑜𝑠 𝜇𝑛 𝑥 𝑐𝑜𝑠 𝜇𝑚 𝑥 𝑑𝑥
0 0 0

On utilise ici la propriété d’orthogonalité de la fonction 𝑐𝑜𝑠 𝜇𝑛 𝑥 .

𝑙
𝑐𝑜𝑠 𝜇𝑛 𝑥 𝑐𝑜𝑠 𝜇𝑚 𝑥 𝑑𝑥 = 0 𝑠𝑖 𝑚 ≠ 𝑛
0
𝑙
𝑙 𝑠𝑖𝑛 2𝜇𝑚 𝑙
𝑐𝑜𝑠 2 𝜇𝑛 𝑥 𝑑𝑥 = + 𝑠𝑖 𝑚 = 𝑛
0 2 4𝜇𝑚
Donc
𝑙 𝑙
𝑠𝑖𝑛 𝜇𝑚 𝑙 𝑙 𝑙 𝑠𝑖𝑛 2𝜇𝑚 𝑙
0 = 0 𝑐𝑜𝑠 𝜇𝑚 𝑥 𝑑𝑥 = 𝑐𝑜𝑠 2 𝜇𝑚 𝑥 𝑑𝑥 = = 𝐴𝑚 +
𝜇𝑚 0 0 2 2 4𝜇𝑚
D’où
2𝑠𝑖𝑛 𝜇𝑚 𝑙
𝐶 = 0
𝜇𝑚 𝑙 + 𝑠𝑖𝑛 𝜇𝑚 𝑙 𝑐𝑜𝑠 𝜇𝑚 𝑙

42
Et la solution générale est :
∞
 2𝑠𝑖𝑛 𝜇𝑛 𝑙 𝑐𝑜𝑠 𝜇𝑛 𝑥
= 𝑒𝑥𝑝 − 𝜇𝑛2 𝑎𝑡
0 𝜇𝑛 𝑙 + 𝑠𝑖𝑛 𝜇𝑛 𝑙 𝑐𝑜𝑠 𝜇𝑛 𝑙
𝑛=1

∞
𝑇 𝑥, 𝑡 − 𝑇∞ 2𝑠𝑖𝑛 𝜇𝑛 𝑙 𝑐𝑜𝑠 𝜇𝑛 𝑥
= 𝑒𝑥𝑝 − 𝜇𝑛2 𝑎𝑡
𝑇0 − 𝑇∞ 𝜇𝑛 𝑙 + 𝑠𝑖𝑛 𝜇𝑛 𝑙 𝑐𝑜𝑠 𝜇𝑛 𝑙
𝑛=1

On peut aussi l’écrire sous la forme :

∞
𝑇 𝑥, 𝑡 − 𝑇∞ 2𝑠𝑖𝑛 𝜔𝑖 𝑥 𝑎𝑡
= 𝑐𝑜𝑠 𝜔𝑖 𝑒𝑥𝑝 − 𝜔𝑖 2 2
𝑇0 − 𝑇∞ 𝜔𝑖 + 𝑠𝑖𝑛 𝜔𝑖 𝑐𝑜𝑠 𝜔𝑖 𝑙 𝑙
𝑖=1

𝜔𝑖 = 𝜇𝑖 𝑙 𝑒𝑠𝑡 𝑓𝑜𝑛𝑐𝑡𝑖𝑜𝑛 𝑑𝑢 𝑛𝑜𝑚𝑏𝑟𝑒 𝑑𝑒 𝐵𝑖𝑜𝑡, 𝜔𝑖 = 𝜔𝑖 𝐵𝑖 .

On définit le nombre de Fourrier , 𝐹𝑂 par :

𝑎𝑡
𝐹𝑂 =
𝑙2
L’expression (****) est tabulée sous la forme :

∞
𝑇 𝑥, 𝑡 − 𝑇∞ 𝑎𝑡 𝑥
= 𝐴𝑖 𝑒𝑥𝑝 − 𝜔𝑖 2 2
𝑐𝑜𝑠 𝜔𝑖
𝑇0 − 𝑇∞ 𝑙 𝑙
𝑖=1
C'est-à-dire :

𝑇 𝑥, 𝑡 − 𝑇∞ 𝑎𝑡 𝑥 𝑎𝑡 𝑥
= 𝐴1 𝑒𝑥𝑝 − 𝜔1 2 2 𝑐𝑜𝑠 𝜔1 +𝐴2 𝑒𝑥𝑝 − 𝜔2 2 2 𝑐𝑜𝑠 𝜔2 + ⋯
𝑇0 − 𝑇∞ 𝑙 𝑙 𝑙 𝑙

La contribution du terme de rang j est d’autant plus faible que j est élevé.

Pour 0.1 ≤ 𝐵𝑖 ≤ ∞ , le terme exponentiel s’amortit très vite et on peut généralement se contenter
du 1ier terme de la série ( à condition que le nombre de Fourrier 𝐹𝑂 ne soit pas trop petit) et on
vérifie ensuite que le deuxième terme est négligeable.

On admet généralement que :

𝑇 𝑥, 𝑡 − 𝑇∞ 𝑎𝑡 𝑥
≅ 𝐴1 𝑒𝑥𝑝 − 𝜔1 2 2 𝑐𝑜𝑠 𝜔1 𝑝𝑜𝑢𝑟 𝐹𝑂 ≥ 0.2
𝑇0 − 𝑇∞ 𝑙 𝑙

En revanche proche du choc thermique 𝑡 ~ 0 , la convergence de cette série est mal assurée, on
utilise la méthode de transformation de Laplace.

Le tableau 3 présente les valeurs des coefficients 𝐴1 et de 𝜔1 pour différents valeurs du nombre
de Biot ainsi qu'une représentation graphique de cette solution approchées. Ces graphiques sont
connus sous le nom de graphes (ou abaques) de Heisler, le premier auteur a les avoir publiés.

43
Temps de relaxation :

Calculons le temps t au bout duquel la température de la plaque est partout quasiment égale à
celle du fluide. Plaçons nous au point le plus éloigné du choc thermique 𝑥 = 0 .

𝑇 0, 𝑡 − 𝑇∞ 𝑎𝑡
≅ 𝐴1 𝑒𝑥𝑝 − 𝜔1 2 2
𝑇0 − 𝑇∞ 𝑙

𝐴1 dépend du nombre de Biot 𝐵𝑖

𝐴1 = 1.016 𝑝𝑜𝑢𝑟 𝐵𝑖 = 0.1 𝑐𝑜𝑛𝑣𝑒𝑡𝑖𝑜𝑛 𝑓𝑎𝑖𝑏𝑙𝑒 .
𝐴1 = 1.273 𝑝𝑜𝑢𝑟 𝐵𝑖 = ∞ 𝑐𝑜𝑛𝑣𝑒𝑡𝑖𝑜𝑛 𝑓𝑜𝑟𝑡𝑒 .
𝐴1 varie peu avec le nombre de Biot, on prend 𝐴1 ~ 1
Si 99% de l’échelon initial est relaxé en tout point de la plaque on a :

𝑇 0, 𝑡 − 𝑇∞ 𝑎𝑡 𝑙2
= 0.01  𝜔1 2 = 4.6  𝑡 = 4.6
𝑇0 − 𝑇∞ 𝑙2 𝑎 𝜔1 2

𝑙 = 15𝑐𝑚 𝐶𝑢𝑖𝑣𝑟𝑒 𝑎 = 10 −4 𝑚2 𝑠 −1 𝐵𝑖 = 0.1 𝑡 = 3.23 𝑕𝑒𝑢𝑟𝑒𝑠

𝐵𝑖 = ∞ 𝑡 = 3.23 𝑕𝑒𝑢𝑟𝑒𝑠
𝑙 = 15𝑐𝑚 𝐴𝑐𝑖𝑒𝑟 𝑎 = 5 10 −6 𝑚2 𝑠 −1 𝐵𝑖 = 0.1 𝑡 = 64.6 𝑕𝑒𝑢𝑟𝑒𝑠
𝐵𝑖 = ∞ 𝑡 = 2.5 𝑕𝑒𝑢𝑟𝑒𝑠
𝑙 = 15𝑐𝑚 𝐵é𝑡𝑜𝑛 𝑎 = 2 10 −7 𝑚2 𝑠 −1 𝐵𝑖 = 0.1 𝑡 = 1615 𝑕𝑒𝑢𝑟𝑒𝑠
𝐵𝑖 = ∞ 𝑡 = 63.3𝑕𝑒𝑢𝑟𝑒𝑠

44
 Tableau 3. Valeurs des coefficients 𝐴1 et de 𝜔1
pour différents valeurs du nombre de Biot pour une plaque plane.

45
 Figure 4. Abaques de Heisler pour une plaque plane

7. 3. Echelon de température appliqué à un cylindre de grande longueur

L’équation de chaleur dans le cylindre s’écrit :

1 𝜕𝜃 𝜕 2 𝜃 1 𝜕𝜃
= +
𝑎 𝜕𝑡 𝜕𝑟 2 𝑟 𝜕𝑟

On utilise la méthode de séparation des variables, on pose

𝜃 𝑟, 𝑡 = 𝜑 𝑡  𝑟
𝜕𝜃
= 𝜑′ 𝑡  𝑟
𝜕𝑡
𝜕𝜃
= 𝜑 𝑡 ′ 𝑟
𝜕𝑟
𝜕2 𝜃
= 𝜑 𝑡  ′′ 𝑟
𝜕𝑟 2
D’où
1
𝜑′ 𝑡  ′′ 𝑟 + 𝑟  ′ 𝑟
=𝑎 = ∓ 𝑎 𝑛2
𝜑 𝑡  𝑟
Et l’équation de Bessel :
1
𝑛2  𝑟 +  ′′ 𝑟 + 𝑟
′ 𝑟 =0

Les solutions de cette équation sont les fonctions de Bessel.

La solution de l’équation de chaleur est donnée par :

∞
𝑇 𝑟, 𝑡 − 𝑇∞ 𝑟
= 𝐵𝑖 𝐽0 𝜔𝑖 𝑒𝑥𝑝 − 𝜔𝑖 2 𝐹𝑂
𝑇0 − 𝑇∞ 𝑅
𝑖=1
𝑎𝑡 𝑕𝑅
𝐹𝑂 = 𝑒𝑡 𝐵𝑖 =
𝑅 2 

46
𝐽0 est la fonction de Bessel de première espèce non modifiée d’ordre zéro. Le tableau 5 présente
les valeurs de des fonctions de Bessel.

2 𝐽1 𝜔𝑖
𝐵𝑖 =
𝜔𝑖 𝐽02 𝜔𝑖 + 𝐽12 𝜔𝑖

On admet généralement que :

𝑇 𝑥, 𝑡 − 𝑇∞ 𝑟
≅ 𝐵1 𝐽0 𝜔1 𝑒𝑥𝑝 − 𝜔1 2 𝐹𝑂 𝑝𝑜𝑢𝑟 𝐹𝑂 ≥ 0.2
𝑇0 − 𝑇∞ 𝑅

Le tableau 4 présente les valeurs des coefficients 𝐵1 et de 𝜔1 pour différents valeurs du nombre de
Biot ainsi qu'une représentation graphique de cette solution approchées. Ces graphiques sont
connus sous le nom de graphes (ou abaques) de Heisler (figure 5).

Tableau 4. Valeurs des coefficients 𝐴1 et de 𝜔1

pour différents valeurs du nombre de Biot pour une cylindre.

47
Tableau 5. Valeurs des fonctions de Bessel.

48
 Figure 5. Abaques de Heisler pour une cylindre

7. 4. Echelon de température appliqué à une sphère pleine

L’équation de chaleur dans la sphère s’écrit :

1 𝜕𝜃 𝜕 2 𝜃 2 𝜕𝜃
= +
𝑎 𝜕𝑡 𝜕𝑟 2 𝑟 𝜕𝑟

On utilise la méthode de séparation des variables, on pose

𝜃 𝑟, 𝑡 = 𝜑 𝑡  𝑟
𝜕𝜃
= 𝜑′ 𝑡  𝑟
𝜕𝑡
𝜕𝜃
= 𝜑 𝑡 ′ 𝑟
𝜕𝑟
𝜕2 𝜃
= 𝜑 𝑡  ′′ 𝑟
𝜕𝑟 2
D’où
2
𝜑′ 𝑡  ′′ 𝑟 + 𝑟  ′ 𝑟
=𝑎 = ∓ 𝑎 𝑛2
𝜑 𝑡  𝑟
Et l’équation de Bessel :
2
𝑛2  𝑟 +  ′′ 𝑟 + ′ 𝑟 =0
𝑟

Les solutions de cette équation sont sous la forme

𝑠𝑖𝑛 𝑛𝑟
 𝑟 =
𝑛𝑟
49
La solution de l’équation de chaleur est donnée par :
∞
𝑇 𝑟, 𝑡 − 𝑇∞ 𝑠𝑖𝑛 𝜔𝑖 𝑟∗
= 𝐶𝑖 𝑒𝑥𝑝 − 𝜔𝑖 2 𝐹𝑂
𝑇0 − 𝑇∞ 𝜔𝑖 𝑟∗
𝑖=1

𝑎𝑡 𝑟
𝐹𝑂 = 𝑒𝑡 𝑟 ∗ =
𝑅2 𝑅

Avec :
𝑠𝑖𝑛 𝜔𝑖 − 𝜔𝑖 𝑐𝑜𝑠 𝜔𝑖
𝐶𝑖 = 2
𝜔𝑖 − 𝑠𝑖𝑛 𝜔𝑖 𝑐𝑜𝑠 𝜔𝑖

Les 𝜔𝑖 sont les valeurs propres de l’équation

𝜔𝑖
𝑡𝑔𝜔𝑖 =
1 − 𝐵𝑖

On admet généralement que :

𝑇 𝑥, 𝑡 − 𝑇∞ 𝑠𝑖𝑛 𝜔1 𝑟∗
≅ 𝐶1 𝑒𝑥𝑝 − 𝜔1 2 𝐹𝑂 𝑝𝑜𝑢𝑟 𝐹𝑂 ≥ 0.2
𝑇0 − 𝑇∞ 𝜔1 𝑟∗

Exemple

Une billette de rayon 𝑅 = 1 𝑐𝑚, de conductivité thermique  = 1 𝑊 𝑚−1 𝐾 −1 de température

𝑇0 = 300°𝐶 est plongée brusquement dans un bain de température 𝑇𝑒 = 15°𝐶 avec un coefficient
de convection de 𝑕 = 50 𝑊 𝑚−2 𝐾 −1 . La diffusivité thermique de la bille est 𝑎 = 2.4 10−7 𝑚2 𝑠 −1

1) Calculer le temps tel que 𝑇 0, 𝑡 = 50°𝐶.

2) Quelle est alors la température pariétale de la sphère.

𝑕𝑅
𝐵𝑖 = = 0.5

D’où
𝜔1 = 1.1656, 𝐶1 = 1.144
𝜔2 = 1.1656, 𝐶2 = −0.2216

1)
𝑠𝑖𝑛 𝑥
→ 1 𝑞𝑢𝑎𝑛𝑑 𝑥 → 0
𝑥

50 − 15 2
= 1.144 𝑒𝑥𝑝 − 1.1656 𝐹𝑂
300 − 15
 𝐹𝑂 = 1.6426
 𝑡 = 684𝑠

2)
𝑇 𝑅, 𝑡 − 𝑇∞ 𝑠𝑖𝑛 1.1656 2
= 1.144 𝑒𝑥𝑝 − 1.1656 × 1.6426
𝑇0 − 𝑇∞ 1.1656

 𝑇 𝑅, 𝑡 = 43°𝐶

50
Le tableau 6 présente les valeurs des coefficients 𝐶1 et de 𝜔1 pour différents valeurs du nombre de
Biot ainsi qu'une représentation graphique de cette solution approchées. Ces graphiques sont
connus sous le nom de graphes (ou abaques) de Heisler (figure 6).

Tableau 6. Valeurs des coefficients 𝐴1 et de 𝜔1

pour différents valeurs du nombre de Biot pour une sphère.

51
 Figure 6. Abaques de Heisler pour une sphère.

7. 5. Problèmes tridimensionnels

On considère un parallélépipède initialement de

température uniforme 𝑇0 , plonge dans un milieu de
température 𝑇∞ , le coefficient d’échange par convection de
la face perpendiculaire à la direction 𝑥𝑖 sera noté 𝑕𝑖 .

Le problème monodimensionnel avait la solution suivante :

∞
𝑇 𝑥, 𝑡 − 𝑇∞ 𝑎𝑡 𝑥
= 𝐴𝑖 𝑒𝑥𝑝 − 𝜔𝑖 2 2
𝑐𝑜𝑠 𝜔𝑖
𝑇0 − 𝑇∞ 𝑙 𝑙
𝑖=1

On peut définir des nombres de Biot suivant les trois axes

𝑥1 , 𝑥2 , 𝑥3 par :
𝑕1 𝑙1 𝑕2 𝑙2 𝑕3 𝑙3
𝐵𝑖1 = , 𝐵𝑖1 = , 𝐵𝑖1 =
  
Considérons la fonction 𝐹1 𝑥1 , 𝑡 solution de l’équation de propagation dans la direction 𝑥1 ∶

𝜕 2 𝐹1 1 𝜕𝐹1
− =0
𝜕𝑥1 2 𝑎 𝜕𝑡

Avec les conditions aux limites :

𝐹1 𝑥1 , 𝑡 = 1

𝐹1 𝑥1 , ∞ = 0

52
 𝜕𝐹1
− = 𝑕1 𝐹1 𝑙1 , 𝑡
𝜕𝑥1 𝑙1

On peut considérer de la même façon, les fonctions similaires 𝐹2 𝑥2 , 𝑡 , 𝐹3 𝑥3 , 𝑡 décrivant les

processus monodimensionnels instationnaires.

Nous allons démontrer que le produit 𝐹1 𝑥1 , 𝑡 𝐹2 𝑥2 , 𝑡 𝐹3 𝑥3 , 𝑡 est solution de l’équation

tridimensionnelle.

Partons de :

𝜕 2 𝐹1 1 𝜕𝐹1
− =0
𝜕𝑥1 2 𝑎 𝜕𝑡
Multiplions par 𝐹2 𝐹3 :
𝜕 2 𝐹1 𝜕 2 𝐹1 𝑥1 𝐹2 𝑥2 𝐹3 𝑥3 1 𝜕𝐹1
𝐹2 𝐹3 = = 𝐹2 𝐹3
𝜕𝑥1 2 𝜕𝑥1 2 𝑎 𝜕𝑡
Permutant les indices
𝜕 2 𝐹2 𝐹1 𝐹3 1 𝜕𝐹2
2
= 𝐹1 𝐹3
𝜕𝑥2 𝑎 𝜕𝑡

𝜕 2 𝐹3 𝐹2 𝐹1 1 𝜕𝐹3
2
= 𝐹2 𝐹1
𝜕𝑥3 𝑎 𝜕𝑡
Si 𝐹 = 𝐹1 𝐹2 𝐹3
𝜕2 𝐹 𝜕2 𝐹 𝜕2 𝐹 1 𝜕𝐹1 1 𝜕𝐹2 1 𝜕𝐹3
2
+ 2
+ 2
= 𝐹2 𝐹3 + 𝐹1 𝐹3 + 𝐹2 𝐹1
𝜕𝑥1 𝜕𝑥2 𝜕𝑥3 𝑎 𝜕𝑡 𝑎 𝜕𝑡 𝑎 𝜕𝑡

1 𝜕𝐹
∆𝐹 =
𝑎 𝜕𝑡
La solution s’écrit donc :
∞
𝑇 𝑥1 , 𝑥2 , 𝑥3 , 𝑡 − 𝑇∞ 2 𝑎𝑡 𝑥1
= 𝐴 1𝑖 𝑒𝑥𝑝 − 𝜔1𝑖 2 𝑐𝑜𝑠 𝜔1𝑖
𝑇0 − 𝑇∞ 𝑙1 𝑙1
𝑖=1

∞
2 𝑎𝑡 𝑥2
× 𝐴 2𝑖 𝑒𝑥𝑝 − 𝜔𝑖2 2 𝑐𝑜𝑠 𝜔𝑖2
𝑙2 𝑙2
𝑖=1

∞
2 𝑎𝑡 𝑥3
× 𝐴 3𝑖 𝑒𝑥𝑝 − 𝜔𝑖3 2 𝑐𝑜𝑠 𝜔𝑖3
𝑙3 𝑙3
𝑖=1

Les 𝐴𝑖 𝑒𝑡 𝜔𝑖 dépendent du nombre de Biot.

Application :

Soit un cube en acier de côté 𝑙 = 30𝑐𝑚, de diffusivité thermique 𝑎 = 1.2 10−5 𝑚2 𝑠 −1 , de

température 𝑇0 = 300°𝐶 est plongée brusquement dans un bain de température 𝑇∞ = 15°𝐶 avec
un coefficient de convection de 𝑕 = 50 𝑊 𝑚−2 𝐾 −1 .

1) Calculer le temps tel que 𝑇 0, 𝑡 = 50°𝐶.

53
 𝑙
𝑕
𝐵𝑖 = 2 = 0.15

𝜔1 = 1.1656, 𝐴1 = 1.0236 (C’est des valeurs moyennes entre 𝐵𝑖 = 0.1 et 𝐵𝑖 = 0.2.
3

𝑇 0, 𝑡 − 𝑇∞ 𝑎𝑡 35
= 𝐴1 𝑒𝑥𝑝 − 𝜔1 2 = = 0.1228
𝑇0 − 𝑇∞ 𝑙 2 285
2

𝑎𝑡
𝐴1 𝑒𝑥𝑝 − 𝜔1 2 2 = 0.497
𝑙
2
 𝐹𝑂 = 5.22
 𝑡 = 9790𝑠 = 2.7𝑕𝑒𝑢𝑟𝑒

7. 6. Mur semi infini mis en contact avec un solide ou un fluide de température uniforme et constante.

Notion de mur semi-infini : Il s’agit d’explorer les régions d’un milieu matériel touchées par une
contrainte thermique alors que d’autres régions ne l’ont pas été encore ; ce modèle s’applique en
particulier aux instants qui suivent le « choc thermique ».

7. 6. 1. Mur semi infini mis en contact avec un solide de température différente et constante.
Considérons une plaque d’épaisseur L, de température uniforme 𝑇0 , dont la face 𝑥 = 0 est
brutalement portée puis maintenue à la température 𝑇1 . Comment le champ de température dans
cette plaque évolue-t-il ?

L’équation de la chaleur dans la plaque s’écrit :

𝜕2 𝑇 1 𝜕𝑇
2
=
𝜕𝑥 𝑎 𝜕𝑡

L’une des méthodes permettant d’accéder à la distribution de la température 𝑇 𝑥, 𝑡 est la méthode

de transformation de Laplace.

Par définition :

54
 ∞
𝑇 𝑥, 𝑝 = 𝑇 𝑥, 𝑡 𝑒 −𝑝𝑡 𝑑𝑡
0
𝑝 ∶ 𝑛𝑜𝑚𝑏𝑟𝑒 𝑟é𝑒𝑙 𝑜𝑢 𝑐𝑚𝑝𝑙𝑒𝑥𝑒.

𝑇 𝑥, 𝑝 𝑖𝑚𝑎𝑔𝑒 𝑜𝑢 𝑡𝑟𝑎𝑛𝑠𝑓𝑜𝑟𝑚é𝑒 𝑑𝑒 𝐿𝑎𝑝𝑙𝑎𝑐𝑒.

Appliquons cette définition à l’équation de la chaleur monodimensionnelle :

𝜕2 𝑇 1 𝜕𝑇
=
𝜕𝑥 2 𝑎 𝜕𝑡
∞
𝜕 2 𝑇 −𝑝𝑡 ∞
1 𝜕𝑇 −𝑝𝑡
2
𝑒 𝑑𝑡 − 𝑒 𝑑𝑡 = 0
0 𝜕𝑥 0 𝑎 𝜕𝑡
𝑑2 𝑇 1 −𝑝𝑡 ∞ ∞
−𝑝
− 𝑇𝑒 + 𝑇 𝑥, 𝑡 𝑒 −𝑝𝑡 𝑑𝑡 = 0
𝑑𝑥 2 𝑎 0 0 𝑎

𝑑2 𝑇 1 𝑝
2
+ 𝑇 𝑥, 0 − 𝑇 𝑥, 𝑝 = 0
𝑑𝑥 𝑎 𝑎

Si 𝑇 𝑥, 0 = 𝑇0 𝑡𝑒𝑚𝑝é𝑟𝑎𝑡𝑢𝑟𝑒 𝑖𝑛𝑖𝑡𝑖𝑎𝑙𝑒 .

𝑑2 𝑇 𝑝 𝑝
2
− 𝑇 𝑥, 𝑝 − 𝑇0 = 0 𝑝𝑢𝑖𝑠𝑞𝑢𝑒 𝑇0 = 𝑇0
𝑑𝑥 𝑎 𝑎

Si on pose 𝑇 ∗ = 𝑇 − 𝑇0
𝑑2 𝑇 ∗
− 𝑘2 𝑇∗ = 0
𝑑𝑥 2
Avec
𝑝
𝑘2 =
𝑎
D’où :
𝑇 ∗ 𝑥, 𝑝 = 𝐴 𝑒 −𝑘𝑥 + 𝐵 𝑒 𝑘𝑥
A et B sont déterminés par les conditions aux limites.

Par les transformées inverse de 𝑇 ∗, on obtient donc 𝑇 ∗ et par conséquent

𝑇 − 𝑇0 = 𝑇 𝑥, 𝑡 − 𝑇0 .

Cas d’un saut de température en 𝒙 = 𝟎

𝜕2 𝑇 ∗ 1 𝜕𝑇 ∗ 𝑑2 𝑇 ∗
= →→ →→ 𝑇 𝑑𝑒 𝐿𝑎𝑝𝑙𝑎𝑐𝑒 →→→→ − 𝑘2 𝑇∗ = 0
𝜕𝑥 2 𝑎 𝜕𝑡 𝑑𝑥 2

𝑇 ∗ 𝑥, 0 = 0 →→→→ 𝑇 𝑑𝑒 𝐿𝑎𝑝𝑙𝑎𝑐𝑒 →→→→ 𝑇 ∗ 𝑥, 0 = 0

𝑇1 − 𝑇0
𝑇 ∗ 0, 𝑡 = 𝑇1 − 𝑇0 →→→→ 𝑇 𝑑𝑒 𝐿𝑎𝑝𝑙𝑎𝑐𝑒 →→→→ 𝑇 ∗ 0, 𝑝 =
𝑝

La solution est donc :

𝑇 ∗ 𝑥, 𝑝 = 𝐴 𝑒𝑥𝑝 −𝑘𝑥 + 𝐵 𝑒𝑥𝑝 −𝑘𝑥

𝑥 → ∞ 𝑙𝑎 𝑡𝑒𝑚𝑝é𝑟𝑎𝑡𝑢𝑟𝑒 𝑑𝑜𝑖𝑡 𝑟𝑒𝑠𝑡𝑒𝑟 𝑓𝑖𝑛𝑖𝑒, 𝑑𝑜𝑛𝑐 𝐵 = 0

55
 𝑇1 − 𝑇0
𝑇 ∗ 0, 𝑝 = 𝐴 + 𝐵 =
𝑝

𝑇1 − 𝑇0
𝐴=
𝑝

𝑇1 − 𝑇0 −𝑘𝑥
𝑇 ∗ 𝑥, 𝑝 = 𝑒
𝑝
𝑥 𝑒𝑥𝑝 − 𝑘𝑥
𝑒𝑟𝑓𝑐 →→→ 𝑇 𝑑𝑒 𝐿𝑎𝑝𝑙𝑎𝑐𝑒 →→→
2 𝑎𝑡 𝑝

D’où :
𝑥
𝑇 ∗ 𝑥, 𝑡 = 𝑇 𝑥, 𝑡 − 𝑇0 = 𝑇1 − 𝑇0 𝑒𝑟𝑓𝑐
2 𝑎𝑡

𝑇 𝑥, 𝑡 − 𝑇1 𝑥
= 𝑒𝑟𝑓𝑐
𝑇0 − 𝑇1 2 𝑎𝑡

Rappel : Fonction erreur

La fonction erreur 𝑛𝑜𝑡é𝑒 𝑒𝑟𝑓 𝑢 est définie par :

𝑢
2
𝑒𝑟𝑓 𝑢 = 𝑒𝑥𝑝 −2 𝑑
𝜋 0

La fonction erreur complémentaire 𝑛𝑜𝑡é𝑒 𝑒𝑟𝑓𝑐 𝑢 est définie par :

𝑒𝑟𝑓𝑐 𝑢 = 1 − 𝑒𝑟𝑓 𝑢

La figure 7 montre l’évolution de la fonction erreur.

Figure 7. Evolution de la fonction erreur.

Dans quelle mesure le modèle semi-infini est il justifié ?

Si
𝑥0
𝑒𝑟𝑓 = 0.99
2 𝑎𝑡

56
On peut considérer que le front thermique est en 𝑥0 à l’instant 𝑡0 , cala donne
𝑥0
~1.8
2 𝑎𝑡

Le temps d’apparition du front thermique vaut :

𝑥0 2
𝑡0 ~ 0.08
𝑎

Le modèle semi infini n’est valide que si 𝑥0 < 𝛿

Exemple

 Aluminium 𝑎 = 8.6 10−5 𝑚2 𝑠 −1 , 𝛿 = 5𝑐𝑚  𝑡0 = 2.32 𝑠

 Béton 𝑎 = 4.2 10−7 𝑚2 𝑠 −1 , 𝛿 = 20𝑐𝑚  𝑡0 = 2.1 𝑕

Vitesse de propagation

Nous définissons la vitesse de propagation par :

𝑑𝑥0 1 𝑎
=
𝑑𝑡0 2 0.08𝑡0
La relation montre que la vitesse se propage comme 𝑎, donc 14 fois plus vite dans l’aluminium
que dans le béton (on découvre le sens physique de la diffusivité thermique 𝑎.

7. 6. 2. Mise en contact brusque de deux milieux.

D’après ce qui précède, la densité de flux qui traverse le plan 𝑥 est donnée par :

𝑑𝑇  𝑥2
𝜑 𝑥, 𝑡 = − = 𝑇1 − 𝑇0 𝑒𝑥𝑝 −
𝑑𝑥 2 𝑎𝑡 4𝑎𝑡

On fait apparaitre un nouveau paramètre appelé effusivité 𝑏 du matériau :


𝑏= = 𝜌𝑐
𝑎
Nous verrons plus loin son sens physique.

Le flux qui traverse la face 𝑥 = 0 est donné par :

𝑏
𝜑 𝑥 = 0, 𝑡 = 𝑇1 − 𝑇0
𝜋𝑡

Lors de la mise en contact de deux murs semi-infinis, le plan frontière prend la température 𝑇𝐶
constante. Le corps 1 cède le flux :
𝑏1
𝜑1 = 𝑇𝐶 – 𝑇1
𝜋𝑡
Le corps 2 reçoit 𝑇1 > 𝑇2 le flux
𝑏2
𝜑2 = 𝑇𝐶 – 𝑇2 .
𝜋𝑡
La conservation du flux 𝜑1 + 𝜑2 = 0 donne :
𝑏1 𝑏2
𝑇𝐶 – 𝑇1 + 𝑇𝐶 – 𝑇2 = 0
𝜋𝑡 𝜋𝑡

57
 𝑏1 𝑇1 + 𝑏2 𝑇2
𝑇𝐶 =
𝑏1 + 𝑏2
𝑇1 +𝑇2
 Si 𝑏1 ~ 𝑏2  𝑇𝐶 = 2
 Si 𝑏1 ≫ 𝑏2  𝑇𝐶 = 𝑇1

Ceci montre que le corps de plus grande effusivité tend à imposer sa température à l’autre corps.

Air Laine de Bois Verre Béton Acier Cuivre

verre plein doux
𝑏 𝐽 𝑚−2 𝐾 −1 𝑠 −
1
2 5.5 30 6 102 1.3 103 2 103 1.2 104 3.5 104
Tableau 7. Effusivité thermique de quelques corps à la température ambiante

Exemple

 Contact entre la main à 37°C et le sable à 55°C

𝑇1 + 𝑇2
𝑇𝐶 = = 46°𝐶
2
 Contact entre la main à 37°C et du métal à 55°C
𝑇𝐶 = 𝑇𝑚 é𝑡𝑎𝑙 = 55°𝐶
Il y a risque de brûlure pour ce dernier cas.

L’effusivité représente la capacité du milieu à imposer sa température.

7. 6. 3. Mur semi infini avec flux pariétal imposé et constant.

𝑇 ∗ 𝑥, 𝑡 = 𝑇 𝑥, 𝑡 − 𝑇0 𝑎𝑣𝑒𝑐 𝑇0 = 𝑇 0, 𝑡

𝑑2 𝑇 ∗
− 𝑘2 𝑇∗ = 0
𝑑𝑥 2
𝑝
𝑘2 =
𝑎

Comme précédemment, la solution de l’équation est :

𝑇 ∗ 𝑥, 𝑝 = 𝐴 𝑒𝑥𝑝 −𝑘𝑥 + 𝐵 𝑒𝑥𝑝 −𝑘𝑥

𝑥 → ∞ 𝑙𝑎 𝑡𝑒𝑚𝑝é𝑟𝑎𝑡𝑢𝑟𝑒 𝑑𝑜𝑖𝑡 𝑟𝑒𝑠𝑡𝑒𝑟 𝑓𝑖𝑛𝑖𝑒, 𝑑𝑜𝑛𝑐 𝐵 = 0

Donc :
𝑇 ∗ 𝑥, 𝑝 = 𝐴 𝑒𝑥𝑝 −𝑘𝑥
La condition en 𝑥 = 0 s’écrit :
𝑑𝑇
− = 𝜑𝑤
𝑑𝑥 𝑥=0

∞ ∞
𝑑𝑇
 𝑒 −𝑝𝑡 𝑑𝑡 = 𝑒 −𝑝𝑡 𝜑𝑤 𝑑𝑡
0 𝑑𝑥 𝑥=0,𝑡 0

𝜑𝑤 𝑑𝑇 ∗ 𝑥, 𝑡
= 
𝑝 𝑑𝑥 𝑥=0

58
Or 𝑇 ∗ 𝑥, 𝑝 = 𝑇 𝑥, 𝑝 − 𝑇0 = 𝐴 𝑒𝑥𝑝 −𝑘𝑥

𝑑𝑇 ∗ 𝑥, 𝑡 𝜑𝑤
= −𝑘 𝐴  𝐴 =
𝑑𝑥 𝑥=0
𝑘𝑝

𝜑𝑤
𝑇 ∗ 𝑥, 𝑝 = 𝑒𝑥𝑝 −𝑘𝑥
𝑘𝑝

La transformée de Laplace inverse (voir tableau) montre que :

𝜑𝑤 𝑎𝑡 𝑥2 𝑥
𝑇 𝑥, 𝑡 − 𝑇0 = 2 𝑒𝑥𝑝 − − 𝑥 𝑒𝑟𝑓𝑐
 𝜋 4𝑎𝑡 4𝑎𝑡

On peut noter le rôle de l’effusivité qui contrôle la montée en température de la paroi.

La température de la paroi est donnée par :

2𝜑𝑤 𝑡
𝑇 0, 𝑡 = 𝑇0 +
𝑏 𝜋

Quant le flux est grand, on peut choisir un matériau d’une effusivité grande pour éviter que la
température de la paroi tend vers une température élevée.

Calculons le flux traversant l’unité de surface en 𝑥, à l’instant 𝑡 :

𝑑𝑇 𝜑𝑤 𝑑 𝑎𝑡 𝑥2 𝑥
𝜑 𝑥, 𝑡 = − = − 2 𝑒𝑥𝑝 − − 𝑥 𝑒𝑟𝑓𝑐
𝑑𝑥  𝑑𝑥 𝜋 4𝑎𝑡 4𝑎𝑡

𝑎𝑡 2𝑥 𝑥2 𝑥 2 𝑥2 1
𝜑 𝑥, 𝑡 = −𝜑𝑤 2 − 𝑒𝑥𝑝 − − 𝑒𝑟𝑓𝑐 − 𝑥 − 𝑒𝑥𝑝 ×
𝜋 4𝑎 4𝑎𝑡 4𝑎𝑡 𝜋 4𝑎𝑡 4𝑎𝑡

𝑥
𝜑 𝑥, 𝑡 = 𝜑𝑤 𝑒𝑟𝑓𝑐
4𝑎𝑡

Remarque

La formule est similaire à celle d’un échelon de température fixé :

𝑇 𝑥, 𝑡 − 𝑇1 𝑥
= 𝑒𝑟𝑓𝑐
𝑇0 − 𝑇1 4𝑎𝑡

7. 6. 4. Mur semi-infini mis en contact avec un fluide de température uniforme et constante.

La face 𝑥 = 0 d’un milieu semi-infini de température 𝑇0 uniforme et constante est brutalement mise
en contact avec un fluide de température 𝑇∞ .

Comme précédemment on a :

59
 𝑇 ∗ 𝑥, 𝑡 = 𝑇 𝑥, 𝑡 − 𝑇0 𝑎𝑣𝑒𝑐 𝑇0 = 𝑇 0, 𝑡

𝑑2 𝑇 ∗
− 𝑘2 𝑇∗ = 0
𝑑𝑥 2
𝑝
𝑘2 =
𝑎

Comme précédemment, la solution de l’équation est :

𝑇 ∗ 𝑥, 𝑝 = 𝐴 𝑒𝑥𝑝 −𝑘𝑥 + 𝐵 𝑒𝑥𝑝 −𝑘𝑥

𝑥 → ∞ 𝑙𝑎 𝑡𝑒𝑚𝑝é𝑟𝑎𝑡𝑢𝑟𝑒 𝑑𝑜𝑖𝑡 𝑟𝑒𝑠𝑡𝑒𝑟 𝑓𝑖𝑛𝑖𝑒, 𝑑𝑜𝑛𝑐 𝐵 = 0

Donc :
𝑇 ∗ 𝑥, 𝑝 = 𝐴 𝑒𝑥𝑝 −𝑘𝑥
Il reste à écrire la condition en 𝑥 = 0.
𝑑𝑇
− = 𝑕 𝑇𝑒 − 𝑇 0, 𝑡
𝑑𝑥 𝑥=0
∞ ∞
𝑑𝑇
𝑒 −𝑝𝑡
− 𝑑𝑡 = 𝑒 −𝑝𝑡 𝑕 𝑇𝑒 − 𝑇 0, 𝑡 𝑑𝑡
0 𝑑𝑥 𝑥=0 0

∞ ∞
𝑑 𝑕𝑇𝑒
− 𝑒 −𝑝𝑡 𝑇 𝑥, 𝑡 𝑑𝑡 = − 𝑕 𝑒 −𝑝𝑡 𝑇 0, 𝑡 𝑑𝑡
𝑑𝑥 0 𝑥=0
𝑝 0

𝑑 𝑇 𝑥, 𝑝
− = 𝑕 𝑇𝑒 − 𝑇 𝑥, 𝑝
𝑑𝑥 𝑥=0
𝑥=0

Or 𝑇 ∗ 𝑥, 𝑝 = 𝑇 𝑥, 𝑝 − 𝑇0 = 𝐴 𝑒𝑥𝑝 −𝑘𝑥

Avec
𝑑𝑇 ∗ 𝑥, 𝑡
= −𝑘 𝐴
𝑑𝑥 𝑥=0
et
𝑇 𝑥, 𝑝 𝑥=0 = 𝐴 + 𝑇0
On a :

𝑘 𝐴 = 𝑕 𝑇𝑒 – 𝐴 − 𝑇0

𝑇𝑒 – 𝑇0 𝑕
𝐴=
𝑘 + 𝑕 𝑝

𝑕 𝑒𝑥𝑝 −𝑘𝑥
𝑇 ∗ 𝑥, 𝑝 =
 𝑝 𝑘+𝑕

La transformée de Laplace inverse donne (voir tableau) :

𝑇 𝑥, 𝑡 − 𝑇0 𝑕  𝑥  𝑕 𝑎𝑕2 𝑥 𝑕
= 𝑒𝑟𝑓𝑐 − 𝑒𝑥𝑝 𝑥+ 𝑡 𝑒𝑟𝑓𝑐 + 𝑎𝑡
𝑇𝑒 – 𝑇0  𝑕 4𝑎𝑡 𝑕   2
4𝑎𝑡 

Si on pose :

60
 𝑥 𝑕
𝑢= 𝑒𝑡 𝑋 = 𝑥
4𝑎𝑡 

𝑢 est le nombre de Fourrier local.

𝑋 est le nombre de Biot local.

𝑇 𝑥, 𝑡 − 𝑇0 𝑋2 
= 𝑒𝑟𝑓 𝑢 + 𝑒𝑥𝑝 𝑋 + 1 − 𝑒𝑟𝑓 𝑋 +
𝑇𝑒 – 𝑇0 4𝑢2 2𝑢

Remarque :

𝑆𝑖 𝐵𝑖 → ∞ , Le fluide impose sa température à la paroi et on retrouve :

𝑇 𝑥, 𝑡 − 𝑇0 𝑥
= 𝑒𝑟𝑓𝑐
𝑇𝑒 − 𝑇0 4𝑎𝑡

7. 6. 4. Mur semi-infini en régime sinusoïdal

Ce type de problème concerne :

- Le chauffage des locaux par suite des apports calorifiques diurnes du soleil.
- Le chauffage des parois du cylindre d’un moteur à piston.

Ici on examine le cas d’un mur semi infini soumis à une variation sinusoidale de la température de
sa surface d’abscisse 𝑥 = 0.

Au bout d’un temps suffisamment long, la température évolue autour d’une valeur moyenne
constante en chaque point.

On s’intéresse ici à la fluctuation 𝜃 𝑥, 𝑡 de la température autour de la valeur moyenne. On ne

s’occupe ici que du régime quasi stationnaire : on suppose que le régime périodique pur est établi,
c'est-à-dire que le terme transitoire relatif aux premiers instants s’est éteint. Le solide a
complètement oublié la distribution initiale de température, de laquelle il ne reste plus de trace.

𝜃 𝑥, 𝑡 = 𝑇 𝑥, 𝑡 − 𝑇𝑚

𝑇𝑚 ∶ 𝑇𝑒𝑚𝑝é𝑟𝑎𝑡𝑢𝑟𝑒 𝑚𝑜𝑦𝑒𝑛𝑛𝑒 𝑒𝑛 𝑥 = 0.

Le problème à résoudre est le suivant :

𝜕2 𝜃 1 𝜕𝜃
2
=
𝜕𝑥 𝑎 𝜕𝑡

La température pariétale fluctue à une pulsation 𝜔 donnée autour de la valeur moyenne 𝑇𝑚 .

𝜃 0, 𝑡 = 𝑇 0, 𝑡 − 𝑇𝑚 = 𝜃0 𝑠𝑖𝑛 𝜔𝑡 = 𝜃 𝑡

𝜃 ∞, 𝑡 → 0

Le problème est de nature périodique, on cherche une solution de nature périodique, dans ce cas
de forme sinusoïdale de même pulsation 𝜔 que la variation de la température de la face
extérieure.
61
On cherche une solution de la forme :
𝑥 𝑥
𝜃 𝑥, 𝑡 = 𝜃0 𝑒𝑥𝑝 − 𝑠𝑖𝑛 𝜔 𝑡 −
 

𝜕𝜃 𝑥 𝑥
= 𝜃0 𝜔 𝑒𝑥𝑝 − 𝑐𝑜𝑠 𝜔 𝑡 −
𝜕𝑡  

𝜕𝜃 𝜃0 𝑥 𝑥 1 𝑥 𝑥
= − 𝑒𝑥𝑝 − 𝑠𝑖𝑛 𝜔 𝑡 − + 𝜃0 − 𝑒𝑥𝑝 − 𝑐𝑜𝑠 𝜔 𝑡 −
𝜕𝑥      

𝜕2 𝜃 𝜃0 1 𝑥 𝑥 1 𝑥 𝑥
= − − 𝑒𝑥𝑝 − 𝑠𝑖𝑛 𝜔 𝑡 − + − 𝑒𝑥𝑝 − 𝑐𝑜𝑠 𝜔 𝑡 −
𝜕𝑥 2       
1 𝑥 𝑥 1 𝑥 𝑥
+ − 𝑒𝑥𝑝 − 𝑐𝑜𝑠 𝜔 𝑡 − + 𝑒𝑥𝑝 − 𝑠𝑖𝑛 𝜔 𝑡 −
     

𝜕2 𝜃 𝜃0 2 𝑥 𝑥
= − − 𝑒𝑥𝑝 − 𝑐𝑜𝑠 𝜔 𝑡 −
𝜕𝑥 2    

𝜕2 𝜃 2𝜃0 𝑥 𝑥
= 2 𝑒𝑥𝑝 − 𝑐𝑜𝑠 𝜔 𝑡 −
𝜕𝑥 2
  

D’où

2𝜃0 𝑥 𝑥 𝑥 𝑥
𝑒𝑥𝑝 − 𝑐𝑜𝑠 𝜔 𝑡 − = 𝜃0 𝜔 𝑒𝑥𝑝 − 𝑐𝑜𝑠 𝜔 𝑡 −
2    

2θ0
= θ0 ω
2

2a
 =
ω
La solution est donc :
𝜔 𝜔
𝜃 𝑥, 𝑡 = 𝜃0 𝑒𝑥𝑝 − 𝑥 𝑠𝑖𝑛 𝜔 𝑡 − 𝑥
2𝑎 2𝑎
𝑥 𝑥
𝜃 𝑥, 𝑡 = 𝜃0 𝑒𝑥𝑝 − 𝑠𝑖𝑛 𝜔 𝑡 −
𝛿 𝛿

𝛿 ∶ é𝑝𝑎𝑖𝑠𝑠𝑢𝑒𝑟 𝑑𝑒 𝑝é𝑛é𝑡𝑟𝑎𝑡𝑖𝑜𝑛

Remarque :

 L’amplitude des oscillations décroit rapidement :

- Lorsque l’on s’éloigne de la face du mur.

- Pour les fréquences croissantes.

 Il existe un déphasage entre les oscillations observées en deux points ;ce déphasage croit avec
la fréquence.

62
Application n°1 :

On suppose que les variations journalières de température sinusoïdale. A quelle profondeur L doit-
on pénétrer dans le sol pour que l’amplitude de la température à cette profondeur, soit inférieure
au 1 100 de l’amplitude en surface sur une période de 20 heurs.

Cette profondeur est caractérisée par l’épaisseur de pénétration 𝛿 :

2𝑎
𝛿=
𝜔

 L’épaisseur de pénétration des fluctuations journalières de la température dans le sol est :

Avec la diffusivité thermique dans le sol 𝑎 = 2.4 10_7 𝑚2 /𝑠.

2 × 2.4 10−7 × 24 × 60 × 60
𝛿= = 0.08𝑚
2𝜋

 L’épaisseur de pénétration des fluctuations annuelles de la température dans le sol est :

2 × 2.4 10−7 × 365 × 24 × 60 × 60

𝛿= = 1.55𝑚
2𝜋

𝜃 𝑥=𝐿 𝜔 1
= 𝑒𝑥𝑝 − 𝐿 =
𝜃 𝑥=0 2𝑎 100

𝜔
𝐿 = 4.605
2𝑎

2𝑎
𝐿 = 4.605
𝜔

2𝜋
𝑇=
𝜔

2𝑎𝑇
𝐿 = 4.605
2𝜋

𝐿 = 2.60 𝑎 𝑇

𝐿 = 2.60 𝑎 𝑇

Avec la diffusivité thermique dans le sol 𝑎 = 2.8 10_7 𝑚2 /𝑠.

𝐿 = 2.60 2.8 10_7 × 20 × 60 × 60

𝐿 = 0.37𝑚

63
Application n°2 :

Dans un moteur à combustion interne, la température des cylindres subit une variation sinusoïdale.
Le modèle du mur semi-infini est-il d’application ?

Prenons un moteur à quatre temps tournant à 4000 tr/mn.

Comme, il y a 2 tours de manivelle par cycle complet, il y a donc 2000 explosions par minute (le
cylindre est rempli une fois quand le moteur fait 2 tours).

La fréquence des oscillations est 2000x60 = 120000 cycles/heure.

Avec la diffusivité thermique de l’acier 𝑎 = 0.04 𝑚2 /𝑠. La profondeur pour laquelle on aura une
atténuation de 0.01 sera :

𝜔 1
𝑒𝑥𝑝 − 𝑥 =
2𝑎 100

𝐿𝑛 0.01 2𝑎
𝑥 = −
𝜔

𝐿𝑛 0.01 2 0.04
𝑥 =−
2 𝜋 × 120000

𝑥 = 1.06 𝑚𝑚

Cette valeur est inférieure à l’épaisseur des parois. Par conséquent la théorie précédente est
d’application. Les oscillations ne se feront pas sentir sur la paroi extérieure. Elles seront réduites à
1/100 de la valeur sur la paroi pour une profondeur de 1.06m.

64
2ième Partie : Convection

I. Introduction
On désigne par convection l’échange de chaleur entre un solide et un fluide en contact avec lui.
Cet échange est régi par la loi de Newton :

𝜑 = 𝑕∆𝑇 = 𝑕 𝑇𝑃 − 𝑇∞ 𝑊/𝑚2

∆𝑇 est l’écart de température entre la paroi et le fluide à une distance suffisante de celle-ci.
L’objectif de ce chapitre est de déterminer le coefficient 𝑕.

Les mouvements du fluide ont une influence considérable sur la convection. Si les mouvements
sont engendrés par une cause extérieure (ventilateur, compresseur, pompe,…..) il s’agit de la
convection forcée ; si au contraire ces mouvements proviennent d’écarts de température au sein
du fluide, on parle de la convection naturelle.

Lorsque 𝑕 est connu, la relation précédente permet de déterminer la quantité de chaleur traversant
une surface donnée pendant un temps donné.

Cette écriture ne fournit aucune indication quant au processus responsable du transfert de

chaleur : il s’agit d’une formulation globale.

Le coefficient 𝑕 dépend de :

- Caractéristiques géométriques de la paroi (forme, rugosité)

- Caractéristiques du fluide : 𝜌, 𝜇, 𝐶𝑃 , , … ….
- De l’écoulement du fluide : U (vitesse), régime laminaire, turbulent.
- De la définition de la température extérieure : 𝑇𝑚𝑜𝑦 , 𝑇𝑚𝑖𝑛 , 𝑇𝑚𝑎𝑥

Valeurs typiques du coefficient de convection 𝒉

𝑕 (𝑊 𝑚−2 𝐾 −1 )

Convection libre 5 - 25
Convection forcée gaz 25 - 250
Convection forcée liquide 50 - 20000
Convection avec changement de 2500 - 100'000
phase ( ébullition, condensation )

𝑕 = 𝑕 𝜌, 𝜇, 𝐶𝑃 , , … … ….

Pour préciser cette relation, on doit d’abord recourir à l’expérience. L’analyse dimensionnelle
permet de réduire le nombre d’essais nécessaires et de les organiser de façon optimale.

65
 II. Analyse dimensionnelle appliquée à la convection

A) Méthode Rayleigh (Analyse dimensionnelle aveugle)

L’objet de l’analyse dimensionnelle est d’exprimer une relation fonctionnelle entre m grandeurs
physiques indépendantes 𝑄1 , 𝑄2 … … … … . 𝑄𝑚 sous forme d’une relation mathématique entre 𝑛 < 𝑚
produits adimensionnels 𝜋1 , 𝜋2 … … … … . 𝜋𝑛 . Chacun d’eux étant une combinaison de puissances
du type :
𝑦
𝑄1𝑥 𝑄2 𝑄3𝑧 … …

Cette méthode est appelée méthode de Rayleigh.

La mécanique des fluides est l’un des champs d’application le plus important de l’analyse
dimensionnelle.

Principales difficultés rencontrées :

- La liste des variables 𝑄1 , 𝑄2 … … … … . 𝑄𝑚 intervenant dans un phénomène donné doit être

en principe établie sur base de l’expérience. Cependant, l’influence de certaines variables
sera plus facilement mis en évidence par un modèle mathématique que par l’expérience :
par exemple l’accélération de la pesanteur 𝑔 qu’il est impossible de faire varier
expérimentalement.

- La liste des unités fondamentales n’est pas toujours évidente. On considère en général
comme fondamentales les unités de longueur, de masse, de temps et de température.

- Les produits adimensionnels 𝜋 peuvent être constitués à partir des grandeurs 𝑄 de

nombreuses façons. Seules certains d’entre elles ont été sanctionnées par la pratique.

B) Méthode de Vaschy-Buckingham
C) Dans la méthode précédente, les dimensions étaient définies à partir des grandeurs
fondamentales. En généralisant cette méthode, on peut dire que l’analyse d’un phénomène
physique comprenant Méthode de Vachy-Buckingham
D) Dans la méthode précédente, les dimensions étaient définies à partir des grandeurs
fondamentales Méthode de Vachy-Buckingham

Dans la méthode précédente, les dimensions étaient définies à partir des grandeurs
fondamentales. En généralisant cette méthode, nous pouvons dire que l’analyse d’un phénomène
physique comprenant 𝑛 variables dont les dimensions sont définies à partir de 𝑝 grandeurs
fondamentales, permet de définir 𝑛 − 𝑝 produits sans dimension indépendants. Cette
généralisation constitue le théorème de Vaschy-Buckingham (ou théorème 𝜋).

Enoncé : Soit un phénomène physique comprenant 𝑛 variables, dans les dimensions desquelles
interviennent 𝑝 grandeurs fondamentales. L’équation :

𝑓 𝑥1 , 𝑥2 ,∙∙∙∙∙∙∙∙∙, 𝑥𝑛 =0

qui régit le phénomène peut se mettre sous la forme :

66
 𝐹 𝜋1 , 𝜋2 ,∙∙∙∙∙∙∙∙∙, 𝜋𝑛 =0

𝜋1 , 𝜋2 ,∙∙∙∙∙∙∙∙∙, 𝜋𝑛 sont des produits sans dimension.

𝑥𝑖
𝜋𝑖 = 𝑎 𝑖 𝑏𝑖 𝑐𝑖
𝑥1 𝑥2 𝑥3 ∙∙∙∙∙

D’après le théorème 𝜋, les produits sans dimension peuvent être comme de nouvelles
variables. ???????????????,,

II.1 Exemples

A) Pertes de charge d’un tube lisse

La variation de la pression sur une longueur 𝐿 d’un tuyau lisse s’écrit :

∆𝑃
= 𝑘 𝐷 𝑎 𝜇𝑏 𝑈 𝑐 𝜌 𝑑
𝐿

Il faut que le membre de gauche ait la même dimension que celui de droite.

𝑃 = 𝑘𝑔 𝑚−1 𝑠 −2

𝜇 = 𝑘𝑔 𝑚−1 𝑠 −1

𝑈 = 𝑚 𝑠 −1

𝜌 = 𝑘𝑔 𝑚−3

𝐷 =𝑚

Nous avons donc :

𝑘𝑔 𝑚−2 𝑠 −2 = 𝑚𝑎 𝑘𝑔𝑏 𝑚−𝑏 𝑠 −𝑏 𝑚𝑐 𝑠 −𝑐 𝑘𝑔𝑑𝑏 𝑚−3𝑑

On déduit le système suivant :

1=𝑎+𝑑

−2 = 𝑎 − 𝑏 + 𝑐 − 3𝑑

−2 = −𝑏 − 𝑐

Nous avons trois équations avec quatre inconnus. Pour résoudre le système, on doit fixer un des
exposants. On déduit le système suivant :

𝑑 =𝑎−𝑏

𝑐 =2−𝑏

𝑎 = −1 − 𝑏

67
L’équation (1) s’écrit :

∆𝑃
= 𝑘 𝐷 −1−𝑏 𝜇𝑏 𝑈 2−𝑏 𝜌−1−𝑏
𝐿

∆𝑃 𝜌 𝑈2 𝜇 𝑏
=𝑘
𝐿 𝐷 𝐷𝑈𝜌

∆𝑃 𝜌 𝑈2
=𝑘 𝑅𝑒 −𝑏
𝐿 𝐷

∆𝑃
𝐿 =  = 𝑓 𝑅𝑒
𝜌 𝑈2
𝐷

Rappel : (cours de mécanique des fluides 3ième année)

Pour un écoulement laminaire :

64
 =
𝑅𝑒

B) Pertes de charge d’un tube rugueux

La variation de la pression sur une longueur L d’un tuyau rugueux s’écrit :

∆𝑃
= 𝑘 𝐷 𝑎 𝜇𝑏 𝑈 𝑐 𝜌 𝑑 𝜀 𝑒
𝐿

1=𝑏+𝑑

−2 = 𝑎 − 𝑏 + 𝑐 − 3𝑑 + 𝑒

−2 = −𝑏 − 𝑐

∆𝑃
= 𝑘 𝐷 −1−𝑏−𝑒 𝜇𝑏 𝑈 2−𝑏 𝜌−1−𝑏 𝜀 𝑒
𝐿

∆𝑃 𝜌 𝑈2 𝜇 𝑏 𝜀 𝑒
=𝑘
𝐿 𝐷 𝐷𝑈𝜌 𝐷
𝜀
𝐷
: est la rugosité relative

∆𝑃
𝐿 =  = 𝑓 𝑅𝑒, 𝜀
𝜌 𝑈2 𝐷
𝐷

68
II.2 Expression adimensionnelle de la loi de convection

L’expérience et l’analyse dimensionnelle de la convection conduisent à admettre que le coefficient

de convection 𝑕 dépend en général de 9 paramètres :

𝑕 = 𝑓 𝑈, 𝐿, 𝜇, , 𝜌, 𝐶𝑃 , ∆𝑇, 𝛽, 𝑔

𝛽 est le coefficient de dilatation à pression constante définit par :

1 𝜕𝑉
𝛽=
𝑉 𝜕𝑇 𝑃

Les trois paramètres ∆𝑇, 𝛽, 𝑔 n’ont d’influence sur le coefficient de convection 𝑕 que dans le cas de
la convection naturelle.

Le coefficient de convection 𝑕 s’écrit :

𝑕 = 𝑘 𝑈 𝑥1 𝐿𝑥2 𝜇 𝑥3 𝑥4 𝜌 𝑥5 𝐶𝑃 𝑥6 ∆𝑇 𝑥7 𝛽𝑔 8

𝑘 et les exposants 𝑥𝑖 sont des nombres purs.

Les différentes grandeurs intervenant dans cette relation ont des dimensions M, L,T,  c'est-à-dire
(Masse, Longueur, Temps, Température).

Nous avons donc 8 grandeurs physiques et 4 grandeurs fondamentales qui définissent le système
d’unité. Le nombre de produits sans dimension est 8-4 = 4.

Dimension L : 𝑥1 + 𝑥2 − 𝑥3 + 𝑥4 − 3𝑥5 + 2𝑥6 + 𝑥8 = 0

Dimension M : 𝑥3 + 𝑥4 + 𝑥5 = 1

Dimension T : 𝑥1 + 𝑥3 + 3𝑥4 + 2𝑥6 + 2𝑥8 = 3

Dimension  : 𝑥4 + 𝑥6 − 𝑥7 + 𝑥8 = 1

Ce système de quatre équations à huit inconnues laisse quatre exposants libres, par exemple
𝑥1 , 𝑥6 , 𝑥7 , 𝑥8 , nous trouvons donc le système suivant :

𝑥4 = 1 − 𝑥6 + 𝑥7 − 𝑥8 = 0

𝑥3 = −𝑥1 + 𝑥6 − 3𝑥7 + 𝑥8 = 0

𝑥5 = 𝑥1 + 2𝑥7

𝑥2 = −1 + 𝑥1 + 2𝑥7 + 𝑥8 = 0

La formule (2) s’écrit :

69
 𝑕 = 𝑘 𝑈 𝑥1 𝐿−1+𝑥 1 +2𝑥 7 +𝑥 8 𝜇 −𝑥 1 +𝑥 6 −3𝑥 7 +𝑥 8 1−𝑥 6 +𝑥 7 −𝑥 8 𝜌 𝑥 1 +2𝑥 7 𝐶𝑃 𝑥6 ∆𝑇 𝑥7 𝛽𝑔 8

𝑥7
𝐿𝑈𝜌 𝑥1
𝜇 𝐶𝑃 𝑥6
𝐿2  𝜌2 ∆𝑇 𝐿𝜇𝛽𝑔 𝑥8
 1
𝑕=𝑘
𝜇  𝜇3  𝐿

 𝐿𝑈𝜌 𝜇 𝐶𝑃 𝐿2  𝜌2 ∆𝑇 𝐿𝜇𝛽𝑔
𝑕= 𝐹 , , ,
𝐿 𝜇  𝜇3 

On reconnait dans cette relation les nombres suivants :

Le nombre de Nusselt :

𝑕𝐿
𝑁𝑢 =

Le nombre de Reynolds :

𝐿𝑈𝜌
𝑅𝑒 =
𝜇

Le nombre de Prandtl :

𝜇 𝐶𝑃
Pr =

Les autres nombres ne sont pas d’utilisation courante. Notons cependant que :

𝐿2  𝜌2 ∆𝑇 𝐿𝜇𝛽𝑔 𝛽𝑔𝐿3 𝜌2 ∆𝑇
× =
𝜇3  𝜇2

Est le nombre de Grashof :

𝛽𝑔𝐿3 𝜌2 ∆𝑇
𝐺𝑟 =
𝜇2

D’autre part le produit :

𝐿2  𝜌2 ∆𝑇 𝜇 𝐶𝑃 𝜇2 𝐶𝑃 ∆𝑇
× × =
𝜇 3  2 2
𝐿 𝑈 𝜌 2 𝑈2

Ce nombre est l’inverse d’un nombre sans dimension connu comme le nombre d’Eckert (Ec) :
−1
𝐶𝑃 ∆𝑇
𝐸𝑐 =
𝑈2

Enfin le produit :

𝐿2  𝜌2 ∆𝑇 𝜇2  ∆𝑇
3
× 2 2 2
=
𝜇 𝐿 𝑈 𝜌 𝜇𝑈 2

Ce nombre est l’inverse d’un nombre sans dimension connu sous le nom de Brinkman (Br) :

70
  ∆𝑇 −1
𝐵𝑟 =
𝜇𝑈 2

Nous avons donc le choix entre les deux relations :

𝑁𝑢 = 𝐹1 𝑅𝑒, Pr, 𝐺𝑟, 𝐸𝑐

Ou

𝑁𝑢 = 𝐹2 𝑅𝑒, Pr, 𝐺𝑟, 𝐵𝑟

Notons que le nombre de Eckert et le nombre de Brinkman sont peu utilisés.

La relation 𝑁𝑢 = 𝐹1 𝑅𝑒, Pr, 𝐺𝑟, 𝐸𝑐 représente l’expression la plus générale de la loi de convection.

En pratique, il arrive fréquemment que la convection forcée soit suffisamment intense pour effacer
l’effet des écarts de température sur le mouvement du fluide ; dans ce cas l’influence du nombre
de Grashof devient négligeable et l’on a :

𝑁𝑢 = 𝑓 𝑅𝑒, Pr

II.3 Signification physique des nombres sans dimension

- Nombre de Nusselt :

𝑕𝐿 𝑕∆𝑇
𝑁𝑢 = =
 ∆𝑇
 𝐿

C’est le rapport d’un flux transmis par conduction et d’un flux transmis par conduction.

- Nombre de Reynolds :

𝜌𝑈 2 𝐹𝑜𝑟𝑐𝑒 𝑑′𝑖𝑛𝑒𝑟𝑡𝑖𝑒
𝑅𝑒 = =
𝑈 𝐹𝑜𝑟𝑐𝑒 𝑑𝑒 𝑣𝑖𝑠𝑐𝑜𝑠𝑖𝑡é
𝜇𝐿

Ce nombre joue un rôle important fondamental dans la caractérisation de l’écoulement.

Si 𝑅𝑒 < 2400, on est en régime laminaire.

Pour des vitesses plus élevées, Si 𝑅𝑒 > 2400, le régime est turbulent.

71
- Nombre de Prandtl :

𝜇
𝜌  𝑉𝑖𝑠𝑐𝑜𝑠𝑖𝑡é 𝑐𝑖𝑛é𝑚𝑎𝑡𝑖𝑞𝑢𝑒
Pr = = =
 𝑎 𝑑𝑖𝑓𝑓𝑢𝑠𝑖𝑣𝑖𝑡é 𝑡𝑕𝑒𝑟𝑚𝑖𝑞𝑢𝑒
𝜌𝐶𝑃

Gaz Pr

H2 0.69

Air 0.69

Ar 0.66

CO2 0.75

CO 0.72

He 0.71

N2 0.70

O2 0.70

H2O vapeur 1.06

Tableau 8. Quelques valeurs du nombre de Pr

72
 Liquide T(°C) Pr
Eau 0 13.6
20 7.03
100 1.75
Alcool éthylique 0 21.8
30 13.9
60 12.1
Glycol 20 203
100 25
Glycérine 0 100000
30 5200
Tableau 9. Variations du nombre de 𝑃𝑟 en fonction de la température

Pour les métaux liquides, au contraire Pr est très petit de l’ordre de 0.01.

- Nombre de Peclet :

𝑈𝜌𝐿𝐶𝑃 𝑈𝜌𝐶𝑃 ∆𝑇
𝑃𝑒 = = 𝑅𝑒 × Pr =
 ∆𝑇
 𝐿

𝑓𝑙𝑢𝑥 𝑐𝑎𝑙𝑜𝑟𝑖𝑓𝑖𝑞𝑢𝑒 𝑙𝑖é à 𝑢𝑛 𝑡𝑟𝑎𝑛𝑠𝑝𝑜𝑟𝑡 𝑑𝑒 𝑚𝑎𝑡𝑖è𝑟𝑒

𝑃𝑒 =
𝑓𝑙𝑢𝑥 𝑐𝑎𝑙𝑜𝑟𝑖𝑓𝑖𝑞𝑢𝑒 𝑡𝑟𝑎𝑛𝑠𝑚𝑖𝑠 𝑝𝑎𝑟 𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑖𝑜𝑛

- Nombre de Stanton (ou nombre de Margoulis Ma) :

𝑕 𝑁𝑢 𝑕∆𝑇
𝑀𝑎 = = =
𝜌𝑈𝐶𝑃 𝑅𝑒 𝑃𝑟 𝜌𝑈𝐶𝑃 ∆𝑇

𝑓𝑙𝑢𝑥 𝑐𝑎𝑙𝑜𝑟𝑖𝑓𝑖𝑞𝑢𝑒 𝑙𝑖é à 𝑙𝑎 𝑐𝑜𝑛𝑣𝑒𝑐𝑡𝑖𝑜𝑛

𝑀𝑎 =
𝑓𝑙𝑢𝑥 𝑐𝑎𝑙𝑜𝑟𝑖𝑓𝑖𝑞𝑢𝑒 𝑙𝑖é à 𝑢𝑛 𝑡𝑟𝑎𝑛𝑠𝑝𝑜𝑟𝑡 𝑑𝑒 𝑚𝑎𝑡𝑖è𝑟𝑒

- Nombre d’Eckert (Ec) :

𝑈2 ∆𝑇𝑎𝑑 𝐴𝑐𝑐𝑟𝑜𝑖𝑠𝑠𝑒𝑚𝑒𝑛𝑡 𝑑𝑒 𝑙𝑎 𝑐𝑕𝑎𝑙𝑒𝑢𝑟 𝑒𝑛 𝑎𝑑𝑖𝑎𝑏𝑎𝑡𝑖𝑞𝑢𝑒

𝐸𝑐 = = 2 =
𝐶𝑃 ∆𝑇 ∆𝑇 𝐸𝑐𝑎𝑟𝑡 𝑑𝑒 𝑡𝑒𝑚𝑝é𝑟𝑎𝑡𝑢𝑟𝑒 𝑑𝑒 𝑟é𝑓é𝑟𝑒𝑛𝑐𝑒

- Nombre de Grashof (Gr) :

𝛽𝑔𝐿3 𝜌2 ∆𝑇 𝑓𝑜𝑟𝑐𝑒 𝑎𝑠𝑐𝑒𝑛𝑡𝑖𝑜𝑛𝑛𝑒𝑙𝑙𝑒

𝐺𝑟 = 2
= 𝑅𝑒 ×
𝜇 𝑓𝑜𝑟𝑐𝑒 𝑣𝑖𝑠𝑞𝑢𝑒𝑢𝑠𝑒

73
 III. Corrélations pour le calcul du coefficient de convection pour différentes
configurations

Ce paragraphe rassemble quelques formules fréquemment utilisées dans les calculs de

transmission de chaleur.

Bien que leurs origines soient généralement expérimentales, leurs formes sont souvent inspirées
des travaux théoriques présentées dans le paragraphe IV.

A- Convection le long d’une surface plane

- Ecoulement laminaire
1 1
𝑁𝑢𝑥 = 0.332 𝑃𝑟 3 𝑅𝑒𝑥 2

1 1
𝑁𝑢𝑥 = 0.664 𝑃𝑟 3 𝑅𝑒𝑥 2

𝑥 est la longueur mesurée depuis le bord d’attaque.

𝑁𝑢𝑥 est le nombre de Nusselt local.
𝑁𝑢𝑥 est le nombre de Nusselt moyen entre 0 𝑒𝑡 𝑥.

- Ecoulement turbulent
1
𝑁𝑢𝑥 = 0.0292 𝑃𝑟 3 𝑅𝑒𝑥 0.8
1
𝑁𝑢𝑥 = 0.036 𝑃𝑟 3 𝑅𝑒𝑥 0.8

B- Convection à l’intérieur de tuyaux cylindriques

- Ecoulement laminaire

 Tubes de grandes longueurs :

𝑁𝑢𝐷 = 3.66

 Tubes de courtes longueurs :

1
𝐷 3 𝜇 0.14
1 1
𝑁𝑢𝐷 = 1.86 𝑅𝑒𝐷 3 𝑃𝑟 3
𝐿 𝜇𝑃

𝜇 viscosité calculée à la température qu’aurait le fluide bien mélangé.

𝜇𝑃 représente la viscosité à la température de la paroi.

- Ecoulement turbulent

𝑁𝑢𝐷 = 0.023 𝑅𝑒𝐷 0.8 𝑃𝑟 𝑛

74
𝑛 = 0.3 𝑜𝑢 0.4 selon que le fluide est refroidi ou réchauffé.

Si l’écart de température entre le fluide et la paroi est important, la formule peut être améliorée
comme suit :

𝜇 0.14
1
𝑁𝑢𝐷 = 0.027 𝑅𝑒𝐷 0.8 𝑃𝑟 3
𝜇𝑃

C - Convection à l’extérieur de tubes cylindriques

2 1 1
𝑃𝑜𝑢𝑟 𝑅𝑒𝐷 < 100 𝑁𝑢𝐷 = + 0.56 𝑅𝑒𝐷 2 𝑃𝑟 3
2𝑙
𝑙𝑛
𝐷
𝑃𝑜𝑢𝑟 𝑅𝑒𝐷 > 100 𝑁𝑢𝐷 = 𝐶 𝑅𝑒𝐷 𝑛 𝑃𝑟 0.35

Les valeurs des constantes 𝐶 et 𝑛 sont données en fonction du nombre de 𝑅𝑒𝐷 dans le tableau ci-
dessus :

𝑅𝑒𝐷 𝐶 𝑛
100 à 1000 0.665 0.47
5000 à 50000 0.22 0.6
4 0.026 0.8
> 10
D - Convection à l’extérieur de tubes cylindriques.

En ligne En quinconce

75
 𝑙 𝑙2
La constante 𝐶 est tabulée en fonction de 𝑎 = 𝑑1 et 𝑏 = 𝑑
Pour les cylindres en ligne : 𝐶 𝑣𝑎𝑟𝑖𝑒 𝑑𝑒 0.66 à 1.07.
Pour les cylindres en quinconce : 𝐶 𝑣𝑎𝑟𝑖𝑒 𝑑𝑒 0.91 à 1.26.

La comparaison avec les résultats du cas d’un tube isolé, montre, comme on pouvait
l’attendre, une augmentation du coefficient d’échange moyen. Cela s’explique par le fait que
seuls les tubes de la première rangée se comportent comme des tubes isolés ; pour les autres
tubes le transfert de chaleur par convection est amélioré par la turbulence due aux sillages
respectifs des divers tubes.

D - Convection naturelle

Dans les formules qui sont données ci-dessus, les caractéristiques du fluide sont rapportées à
la température du film :
𝑇𝑃 + 𝑇𝑒
𝑇𝑓 =
2

 Plaque plane verticale, de hauteur 𝐿 et à la température uniforme 𝑃𝑟 > 0.6

1
𝑁𝑢𝐿 = 0.53 𝐺𝑟𝐿 𝑃𝑟 4 𝑝𝑜𝑢𝑟 𝐺𝑟𝐿 𝑃𝑟 < 109

2
𝑁𝑢𝐿 = 0.021 𝐺𝑟𝐿 𝑃𝑟 5 𝑝𝑜𝑢𝑟 109 < 𝐺𝑟𝐿 𝑃𝑟 < 1012

 Long cylindre horizontal de diamètre 𝐷 et à la température uniforme 𝑃𝑟 > 0.6

1
𝑁𝑢𝐷 = 0.52 𝐺𝑟𝐷 𝑃𝑟 4 𝑝𝑜𝑢𝑟 103 < 𝐺𝑟𝐷 𝑃𝑟 < 109

1
𝑁𝑢𝐷 = 0.1 𝐺𝑟𝐷 𝑃𝑟 3 𝑝𝑜𝑢𝑟 𝐺𝑟𝐷 𝑃𝑟 > 109

 Sphère de rayon 𝑅 et à la température uniforme 𝑃𝑟 > 0.6

1
𝑁𝑢𝑅 = 0.53 𝐺𝑟𝑅 𝑃𝑟 4 𝑝𝑜𝑢𝑟 𝐺𝑟𝐷 > 103

 Plaque horizontale carrée de coté 𝐿 et à la température uniforme 𝑃𝑟 > 0.6

- Si la face chaude est tournée vers le haut ou si la face froide (par rapport au fluide) est
tournée vers le bas.
1
𝑁𝑢𝐿 = 0.54 𝐺𝑟𝐿 𝑃𝑟 4 𝑝𝑜𝑢𝑟 105 < 𝐺𝑟𝐿 < 2 107

1
𝑁𝑢𝐿 = 0.14 𝐺𝑟𝐿 𝑃𝑟 3 𝑝𝑜𝑢𝑟 𝐺𝑟𝐷 𝑃𝑟 > 2 107

- Si la face chaude est tournée vers le bas ou si la face froide (par rapport au fluide) est
tournée vers le haut, la convection est moins active.
1
𝑁𝑢𝐿 = 0.68 𝐺𝑟𝐿 𝑃𝑟 5 𝑝𝑜𝑢𝑟 3 105 < 𝐺𝑟𝐿 < 108
1
𝑁𝑢𝐿 = 0.27 𝐺𝑟𝐿 𝑃𝑟 4 𝑝𝑜𝑢𝑟 108 < 𝐺𝑟𝐿 < 3 1010

76
3ième Partie : Rayonnement Thermique

Introduction
Un système S isolé dans l’espace, c’est à dire sans contact matériel avec d’autres corps émet un
rayonnement qui peut se perdre à l’infini où être capté par un autre système S’. Le système S perd
donc continuellement de l’énergie dont une partie est reçu par S’. Tant que l’énergie interne du
corps S n’est pas nulle ; tant que la température absolue n’est pas nulle S continue à émettre un
rayonnement.

Par ailleurs, les objets qui entourent S lui communiquent une partie de leurs énergies ce qui évite à
S de se refroidir complètement. Il arrive un moment où S atteint un état d’équilibre lorsque
l’énergie qu’il capte est égale à celle qu’il emet si par ailleurs, il n’a aucun contact avec d’autres
objets.

Dans la pratique, on constate que le rayonnement peut contenir des radiations visibles, infrarouges
ou ultraviolettes. On l’appelle rayonnement thermique du système.

Tout corps matériels rayonnent ; c'est-à-dire libèrent en permanence sous forme non matérielle
qu’on appelle « rayonnement électromagnétique ».

- Ce phénomène peut être considéré de point de vue ondulatoire (ondes électromagnétiques

soumises aux équations de Maxwell).
- Ce phénomène peut être considéré de point de vue corpusculaire (particules sans masses
ni charges électriques appelées « photons »).

Le rayonnement est constitué de photons caractérisés par leurs directions de propagation  et

par leur fréquence . Chaque photon se propage dans un milieu d’indice de réfraction n transporte
une énergie 𝐸 = 𝑕 et une quantité de mouvement

𝑕
𝑝=𝑛 
𝐶0

𝑕 = 6.625 10−34 𝐽 𝑠 𝑒𝑠𝑡 𝑙𝑎 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 𝑑𝑒 𝑃𝑙𝑎𝑛𝑐𝑘.

𝐶0 = 3 108 𝑚 𝑠 𝑒𝑠𝑡 𝑙𝑎 𝑣𝑖𝑡𝑒𝑠𝑠𝑒𝑑𝑒 𝑝𝑟𝑜𝑝𝑎𝑔𝑎𝑡𝑖𝑜𝑛 𝑑𝑒 𝑙𝑎 𝑙𝑢𝑚𝑖è𝑟𝑒.

Dans la littérature, le photon est généralement caractérisé par sa longueur d’onde dans le vide
𝐶0 𝐶
0 = 
et dans le milieu considéré par  = 𝑛0 plutôt que sa fréquence.

L’échange énergétique entre un système matériel et un champ de rayonnement s’effectue selon

deux processus :

 Emission : conversion de l’énergie matérielle (translation, rotation, vibration) en énergie radiative

(production de photons).

 Absorption : le processus inverse au cours duquel des photons disparaissent en cédant

l’intégralité de leur énergie et de leur quantité de mouvement au milieu matériel.

Les phénomènes de réflexion et de diffusion de la lumière correspondent à un simple échange de

quantité de mouvement (déviation de photons) sans échange énergétique.

77
Lorsque l’émission se fait au détriment de l’énergie thermique du système, on parle de
« rayonnement thermique ». Les photons émis se répartissent de façon continue sur un domaine
spectral qui dépend de la température.

Lorsque l’émission est emprunté à l’excitation du milieu, on parle de » rayonnement de

fluorescence ». Son spectre est caractérisé par des fortes discontinuités (émission de raies
atomiques et de bandes moléculaires.

Exemple : Une flamme riche en combustible est l’illustration de ces deux cas de figures
(émission). Le spectre se compose de bades de rotation-vibration émises par les radicaux libres
excités (OH, CH, CN, NO,………) et du fond continu émis les particules de suies présentes dans la
flamme.

Fluorescence

On considère un tube d’éclairage domestique : les électrons libres accélérés par le champ
électrique produisent des collisions inélastiques qui excitent les atomes, c’est à dire modifient la
structure de leur cortège électronique. L’état excité ainsi produit est généralement instable.
L’atome se désexcite rapidement soit au cours d’une nouvelle collision, sous forme radiative : le
rayonnement de fluorescence ainsi produit atteste du caractère quantique des niveaux d’énergie
atomique ou moléculaire.

Notons que certains niveaux excités appelés « métastables » ont une grande durée de vie : le flux
radiatif produit par leur lente désexcitation est le rayonnement de « phosphorescence ».

Ils existent trois processus radiatifs qui participent au peuplement des niveaux atomiques à savoir :

- L’émission spontanée.
- L’absorption.
- L’émission induite (stimulée).

Emission spontanée :

Un atome M dans l’état d’énergie 𝐸𝑗 se désexcite spontanément vers un état d’énergie 𝐸𝑗 en

émettant un photon d’’énergie 𝑕  𝑖 𝑗 = 𝐸𝑗 − 𝐸𝑖

Absorption :

Un atome M dans l’état d’énergie 𝐸𝑖 absorbe un photon dont l’énergie est celle de la transition
𝑖 → 𝑗.

𝑀𝑖 + 𝑕  𝑖 𝑗 → 𝑀𝑗

Emission induite (ou stimulée) :

Pour rendre les processus de peuplement compatible avec la loi de Planck, Einstein propose un
troisième type d’interaction entre atome et champ électromagnétique.

L’émission se fait dans des conditions très originales : le photon émis est de même fréquence,
même direction et même phase que le rayonnement inducteur. Ce phénomène de résonance est
émis à profit dans la conception des sources LASER (Light Amplification by Stimulated Emission
of Radiation).

78
 Rappel : Angle solide

1- Vecteur élément de surface

𝑑𝑆 = 𝑛 ∙ 𝑑𝑆

2 - Angle solide

C’est la manière de caractériser la portion de l’espace intérieure à un cône.

C’est une généralisation dans l’espace à trois dimensions de la notion d’angle plan.
C’est la rapport entre l’aire découpée sur une sphère 𝛴 et le carré de son rayon.

3 - Angle solide élémentaire

𝑑𝑆 ∙ 𝑢 𝑑𝛴 𝑑𝑆 𝑐𝑜𝑠𝛼
𝑑 = 2
= 2 =
𝑟 𝑟 𝑟2

L’angle solide est sans dimension ; l’unité de l’angle solide est le stéradian.

Exemple

L’angle solide 𝑑 sous lequel on voit une couronne élémentaire d’un disque depuis un point O de
son axe.

79
 𝑑𝑆 𝑐𝑜𝑠𝛼
𝑑 =
𝑙2

𝑑𝑆 = 2 𝜋 𝑟 𝑑𝑟

𝑟 = 𝑎 𝑡𝑔𝛼

𝑎 𝑑𝛼
𝑑𝑟 =
𝑐𝑜𝑠 2 𝛼
𝑎
𝑙=
𝑐𝑜𝑠𝛼
2 𝜋 𝑟 𝑑𝑟 𝑐𝑜𝑠𝛼
𝑑 =
𝑙2

𝑎 𝑑𝛼
2𝜋𝑟 𝑐𝑜𝑠𝛼 2 𝜋 𝑟 𝑎 𝑑𝛼 2 𝜋 𝑎 𝑡𝑔𝛼 𝑎 𝑑𝛼
𝑑 = 𝑐𝑜𝑠 2 𝛼 = =
𝑙 2 2
𝑙 𝑐𝑜𝑠𝛼 𝑙 2 𝑐𝑜𝑠𝛼

𝑑 = 2 𝜋 𝑠𝑖𝑛𝛼 𝑑𝛼

Application

L’angle solide sous lequel on voit le disque (angle 𝜃) :

𝑑 = 2 𝜋 𝑠𝑖𝑛𝛼 𝑑𝛼
𝜃
 = 2 𝜋 𝑠𝑖𝑛𝛼 𝑑𝛼 = 2 𝜋 1 − 𝑐𝑜𝑠𝜃
0

Remarque

- Le stéradian est l’aire de 1𝑚2 , intercepté sur une sphère de rayon 1𝑚, par un cône dont le
sommet est placé au centre de cette sphère.
- Un angle solide sphérique vaut 4𝜋 𝑆𝑡é𝑟𝑎𝑑𝑖𝑎𝑛.
80
I. Grandeurs énergétiques.

I. 1. Flux total
Le flux énergétique étant la quantité d’énergie transportée sous forme de radiations pendant l’unité
de temps ; il caractérise le rayonnement total émis dans toute les directions par une surface
émettrice. On l’appelle ∅. ∅ = 𝑊𝑎𝑡𝑡𝑠 𝑊 .

I. 2. Emittance ou Radiance M
Pour rendre ce flux indépendant des dimensions de la surface émettrice, on introduit l’émittance
énergétique ou la radiance M, rapport du flux émis dans toutes les directions par unité de surface
émettrice.

∅
𝑀=
𝑆

𝑀 = 𝑊 𝑚−2 .

I. 3. Intensité I
L’intensité I caractérise le flux énergétique émis par une surface S dans une direction donnée 𝑂𝑥
(figure). Soit ∅𝑂𝑥 le flux élémentaire contenu dans un angle solide 𝑑 autour de la direction 𝑂𝑥.

Dans la direction 𝑂𝑥, le flux énergétique ∅𝑂𝑥 est tel

que :

𝑑∅𝑂𝑥
∅𝑂𝑥 = = 𝐼𝑂𝑥
𝑑

𝐼 = 𝑊 𝑠𝑡 −1 .

𝑑∅𝑂𝑥 est la portion de flux rayonné dans un angle

solide élémentaire 𝑑 entourant la direction 𝑂𝑥.

𝐼𝑂𝑥 est le flux par unité d’angle solide.

I. 4. Luminance
La surface émettrice S fait un angle avec la normale 𝛼 avec la normale à la direction dans laquelle
on a évalué l’intensité énergétique. La projection de S sur cette normale est 𝑆𝑐𝑜𝑠𝛼, appelée
surface apparente  et l’intensité énergétique par unité de surface apparente définit la
luminance.

𝐼 𝑑∅
𝐿= =
 𝑑 𝑆𝑐𝑜𝑠𝛼
81
 𝐿 = 𝑊 𝑠𝑡 −1 𝑚−2 .

Le flux émis par un élément de surface 𝑑𝑆 dans un angle solide 𝑑 entourant une direction 𝑂𝑥,
incliné d’un angle 𝜃 sur la normale à cette surface, a pour expression :

𝑑2 ∅𝑂𝑥 = 𝐿𝑂𝑥 𝑑 𝑑𝑆𝑐𝑜𝑠𝜃

I. 5. Grandeurs monochromatiques

Toutes les grandeurs introduites précédemment à propose du rayonnement total (c'est-à-dire

l’ensemble des longueurs d’onde rayonnées par une source) peuvent être également rapportées
aux échanges confinés dans un domaine spectral étroit, de longueur 𝑑 (ou 𝑑) centré sur une
longueur d’onde  (ou une fréquence ).

On définit le flux énergétique :

𝑑∅
∅ =
𝑑

De la même manière, on définit :

𝑑𝑀
𝑀 =
𝑑
𝑑𝐿
𝐿 =
𝑑

On définit une relation entre flux énergétique et flux monochromatique ou spectral :

∞
∅= ∅ 𝑑
0

De même :
∞
𝑀= 𝑀 𝑑
0

∞
𝐿= 𝐿 𝑑
0

I. 6. Eclairement
Les notions de flux, d’intensité et de luminance s’appliquent aussi bien au rayonnement incident
sur une surface qu’au rayonnement émis par celle-ci. Par contre la notion d’émittance est
remplacée par l’éclairement.

L’éclairement est le flux reçu par l’unité de surface réceptrice, en provenance de l’ensemble des
directions d’où elle peut recevoir du rayonnement.

𝑑∅
𝐸=
𝑑𝑆

82
 𝐸 = 𝑊 𝑚−2

𝑑𝑆 est l’aire de la surface recevant un flux 𝑑∅.

I. 7. Relation entre l’éclairement du récepteur et la luminance de l’émetteur.

Le flux émis par une surface 𝑑𝑆2 en direction d’une surface réceptrice 𝑑𝑆1 s’écrit :

𝑑2 ∅2 = 𝐿2 𝑑𝑆2 𝑐𝑜𝑠𝜃2 𝑑2

𝑑𝑆1 𝑐𝑜𝑠𝜃1
𝐴𝑣𝑒𝑐 𝑑2 = , 𝑜𝑛 𝑎 ∶
𝐷2
𝑑𝑆2 𝑐𝑜𝑠𝜃2 𝑑𝑆1 𝑐𝑜𝑠𝜃1
𝑑 2 ∅ 2 = 𝐿2
𝐷2

L’éclairement de 𝑑𝑆1 sera donc :

𝑑 2 ∅2 𝑑𝑆2 𝑐𝑜𝑠𝜃2 𝑐𝑜𝑠𝜃1

𝐸= = 𝐿2
𝑑𝑆1 𝐷2

Comme précédemment, on définit un éclairement monochromatique 𝐸 (ou 𝐸 ).

II. Lois du rayonnement

II. 1. Loi de Lambert

Soit 𝐼𝑛 l’intensité énergétique selon la normale à la surface émettrice 𝑆 et 𝐼𝛼 l’intensité énergétique
dans une direction faisant l’angle 𝛼 avec la normale.

On a :

𝐼𝑛
𝐿𝑛 =
𝑆

83
 𝐼𝛼 𝐼𝛼
𝐿𝛼 = =
 𝑆𝑐𝑜𝑠𝛼

On dit que le rayonnement suit la loi de Lambert si la relation : 𝐼𝛼 = 𝐼𝑛 𝑐𝑜𝑠𝛼 est vérifiée.

Par conséquent :

𝐼𝛼 𝐼𝛼 𝐼𝑛 𝑐𝑜𝑠𝛼
𝐿𝛼 = = = = 𝐿𝑛
 𝑆𝑐𝑜𝑠𝛼 𝑆𝑐𝑜𝑠𝛼

La luminance énergétique est la même dans toutes les directions lorsqu’un rayonnement suit la loi
de Lambert (la réciproque est vraie).

D’autre part, on montre qu’un rayonnement suit la loi de Lambert est tel que l’émittance et la
luminance sont proportionnelles :

𝑀 = 𝜋𝐿

II. 2. Réception du rayonnement par un corps

Lorsqu’un rayonnement atteint la surface d’un corps, une fraction 𝑟 de la puissance totale est
réfléchie, c'est-à-dire envoyée sans avoir pénétré dans le matériau, une autre fraction 𝑎 est
absorbée dans la masse du récepteur constituant un apport énergétique pour ce dernier. Une
fraction 𝑡 peut être transmise (traverse le corps et ressort de l’autre côté).

On définit les facteurs de réflexion, transmission et absorption par :

∅𝑟
𝑟=
∅𝑖

∅𝑡
𝑡=
∅𝑖

∅𝑎
𝑎=
∅𝑖

Où ∅𝑟 , ∅𝑡 et ∅𝑎 sont respectivement le flux réfléchi, le flux transmis et le flux absorbé par le

corps. ∅𝑖 est flux incident.

D’après la loi de conservation on trouve :

𝑟+𝑡+𝑎 = 1

Pour une longueur d’onde donnée on a :

𝑟 + 𝑡 + 𝑎 = 1

C’est naturellement le coefficient d’absorption 𝑎 qui est à la base de l’étude du rayonnement.

II. 3. Corps noir

84
On appelle corps noir toute surface dont le facteur d’absorption est égal à l’unité, quelle que soit la
longueur d’onde de la radiation incidente .

𝑎 = 𝑎 = 1 ∀ 𝑙𝑎 𝑙𝑜𝑛𝑔𝑢𝑒𝑢𝑟 𝑑 ′ 𝑜𝑛𝑑𝑒  (𝑟 = 𝑡 = 0).

C’est un concept théorique vers lequel tendent certaines substances comme le noir de platine, le
noir de Bismuth ou le noir du camphre 𝑎 > 0.92.

Les grandeurs relatives au corps noir seront affectées d’un indice supérieur "0" à droite, par
exemple :

𝑀0 désignera l’émittance monochromatique du corps noir à la longueur d’onde .

𝐿0 désignera la luminance monochromatique du corps noir à la longueur d’onde .

Un tel corps apparait à température modéré effectivement noir car il

absorbe toutes les radiations du spectre visible et ne réfléchit donc
rien vers l’œil.

On peut réaliser un corps noir en prenant une cavité recevant du

rayonnement par une faible ouverture percée sur sa paroi. La face
intérieure de cette paroi sera aussi absorbante que possible (figure).

L’énergie radiante qui pénètre à l’intérieur est réfléchie un grand

nombre de fois avant de pouvoir sortir. A chaque réflexion, une
grande fraction de l’énergie est absorbée de manière à ce qui peut
en sortir est extrêmement faible. Donc tout ce qui tombe dans le trou
O est donc absorbé. Pratiquement
∀ 𝑙𝑎 𝑡𝑒𝑚𝑝é𝑟𝑎𝑡𝑢𝑟𝑒 𝑇, ∀ 𝑙𝑎 𝑙𝑜𝑛𝑔𝑢𝑒𝑢𝑟 𝑑 𝑜𝑛𝑑𝑒  , le trou O se comporte
′

comme un corps noir de point de vue absorption.

II. 4. Loi de Planck

Cette loi relie l’émittance monochromatique du corps noir, 𝑀0 à la longueur d’onde  et à la
température 𝑇 𝐾𝑒𝑙𝑣𝑖𝑛 . Elle s’exprime sous la forme :

2 𝜋 𝑕 𝑐 2 −5
𝑀0 = 𝑕𝑐
𝑒 𝑘𝑇 − 1

𝑐 est la vitesse des ondes électromagnétiques dans le milieu où se propage le rayonnement.

𝑐0
𝑐=
𝑛

𝑐0 𝑒𝑠𝑡 𝑙𝑎 𝑣𝑖𝑡𝑒𝑠𝑠𝑒 𝑑𝑒 𝑙𝑎 𝑙𝑢𝑚𝑖è𝑟𝑒 𝑐0 = 3 108 𝑚 𝑠 −1

𝑛 𝑒𝑠𝑡 𝑙 ′ 𝑖𝑛𝑑𝑖𝑐𝑒𝑑𝑒 𝑟é𝑓𝑟𝑎𝑐𝑡𝑖𝑜𝑛 𝑑𝑢 𝑚𝑖𝑙𝑖𝑒𝑢.

𝑕 𝑒𝑠𝑡 𝑙𝑎 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 𝑑𝑒 𝑃𝑙𝑎𝑛𝑐𝑘 𝑕 = 6.625 10−34 𝐽 𝑠

𝑘 𝑒𝑠𝑡 𝑙𝑎 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 𝑑𝑒 𝐵𝑜𝑙𝑡𝑧𝑚𝑎𝑛 𝑘 = 1.3805 10−23 𝐽 𝐾 −1

85
Lorsque le rayonnement se propage dans un milieu dont l’indice de réfraction 𝑛 est égal à l’unité
(c'est le cas du vide ou de l’air), en première approximation, la loi de Planck s’exprime sous la
forme suivante, qui celle que l’on utilise dans la pratique :

𝐶1 −5
𝑀0 = 𝐶2
𝑒𝑇 − 1

Avec 𝐶1 et 𝐶2 sont des constantes.

𝐶1 = 2 𝜋 𝑕 𝑐02

𝑕 𝑐0
𝐶2 =
𝑘

𝐶1 = 3.741 10−16 𝑊 𝑚2

𝐶2 = 0.014388 𝑚 𝐾

Les courbes spectrales du rayonnement noir 𝑀0 = 𝑓  sont du type Gauss avec dissymétrie
(figure 7).

Figure 7. Evolution de l’émittance d’un corps noir en fonction

de la longueur d’onde et de la température.

La chute d’énergie est beaucoup plus lente après le maximum. La croissance avec  (courtes
longueurs d’onde) est beaucoup plus rapide que leur décroissance.

86
Approximation de la loi de Planck

Pour les courtes longueurs d’onde  < 5𝜇𝑚 ∶ 𝑣𝑖𝑠𝑖𝑏𝑙𝑒 𝑒𝑡 𝑝𝑟𝑜𝑐𝑕𝑒 𝑖𝑛𝑓𝑟𝑎𝑟𝑜𝑢𝑔𝑒 , la formule de Planck
peut être approchée par l’expression dite formule de Wien:
𝐶2
𝑀0 = 𝐶1 −5 𝑒 −  𝑇

On néglige 1 du dénominateur de la formule de Planck devant l’exponentielle. Le domaine de

validité de cette expression est :  𝑇 < 14000 𝜇𝑚 𝐾.

Loi de déplacement de Wien

L’abscisse 𝑚 du maximum de 𝑀0 se déplace vers les courtes longueurs d’onde lorsque la
température croit. C'est ce exprime la « loi de déplacement de Wien :

𝑚 𝑇 = 2898 𝜇𝑚 𝐾

Remarque :

Le déplacement de 𝑚 vers les courtes longueurs d’onde lorsque la température croît, explique le
fait que, pendant son échauffement un corps n’émet d’abord que dans l’infrarouge (radiation non
visible) puis commence à rougir lorsque son spectre vient empiéter sur le domaine du visible (en
commençant par le rouge). Pour une température suffisamment élevée, le spectre d’émission du
corps recouvre le domaine du visible ce qui correspond à une émission de lumière blanche.

L’émittance énergétique totale du corps noir à une température donnée :

∞
𝑀0 = 𝑀0 𝑑
0

Est représentée par l’aire comprise sous la courbe spectrale à cette température.

II. 5. Loi de Stephan-Boltzman

L’émittance énergétique totale du rayonnement noir est proportionnelle à la puissance 4 de la

température absolue de l’émetteur noir.

∞ ∞
𝐶1 −5
𝑀0 = 𝑀0 𝑑 = 𝐶2 𝑑 = 𝜍𝑇 4
0 0 𝑒  𝑇 − 1

Cette loi a été établie par Stephan en 1879 et retrouvée à partir de la thermodynamique statistique
par Boltzman en 1884, elle introduit la constante de Stephan 𝜍.

𝜍 = 5.67 10−8 𝑊𝑚−2 𝐾 −4

87
II. 6. Facteur d’émission d’un corps non noir

L’étude du rayonnement d’un corps non noir se fait en comparant le rayonnement obtenu à une
certaine température, avec le rayonnement du corps noir émettant à la même température.
Par définition :
𝑀
𝜖 = 0
𝑀

𝜖 facteur d’émission monochromatique.

𝑀
𝜖=
𝑀0

𝜖 facteur total d’émission.

𝜖 et 𝜖 varient avec la température et l’état de surface.

Substance Etat Emissivité

Aluminium Poli 0.02 – 0.06

Cuivre Poli 0.03 – 0.04

Cuivre Oxydé 0.5 – 0.8

Asphalte - 0.85 – 0.93

Bois - 0.82 – 0.92

Neige - 0.82 – 0.9

Végétation - 0.92 – 0.96

Peau - 0.95

Emissivité moyenne : Exemples

II. 7. Corps en équilibre : Loi de Kirchhoff

Soient deux corps de faibles dimensions et de surface

respectives 𝑆1 et 𝑆2 situés loin l’un de l’autre à l’intérieur
d’une grande enceinte 𝑆 parfaitement isolée et maintenue
à une température uniforme.

Au bout d’un certain temps, il s’établit à la suite des

échanges par rayonnement entre les corps et les parois
rayonnantes de l’enceinte, un équilibre thermique.

Pour la longueur d’onde  il « tombe » par unité de

surface sur chaque corps l’énergie  .
88
Le corps (1) de coefficient d’absorption 𝑎1  absorbe le flux énergétique 1 telle que :


1 = 𝑎1   𝑆1


De même le corps (2) absorbe le flux 2 telle que :



2 = 𝑎2   𝑆2


En tant qu’émetteur, le corps (1) émet à la longueur d’onde  le flux :

′1  = 𝑀1  𝑆1

De même le corps (2) émet le flux ′2  telle que :

′ 2 = 𝑀2  𝑆2


A l’équilibre :

1 = ′1  et 2 = ′2 
 

Soit :

𝑎1   𝑆1 = 𝑀1  𝑆1

𝑎2   𝑆2 = 𝑀2  𝑆2

𝑎1  𝑀1 
=
𝑎2  𝑀2 

𝑎1  𝑎2 
= = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
𝑀1  𝑀2 

Le rapport est constant pur tous les corps. Il est indépendant de la forme et de la nature du corps.
En particulier si le corps (2) est un corps noir, donc :

𝑎2  = 1 et 𝑀2  = 𝑀0 , d’où :
𝑀 𝑀0
=
𝑎 1
𝑀 𝑀0
=
𝑎 1
𝑀 = 𝑎 𝑀0 Loi de Kirchoff

Par comparaison de la définition du facteur d’émission monochromatique et de la loi de Kirchoff,

on déduit que :

𝑎 = 𝜀

Par intégration sur toute la courbe spectrale, on trouve :

89
 𝑎= 𝜀

L’étude du comportement radiatif des matériaux a cependant montré que la relation

𝑎 = 𝜀

Reste valable même en dehors des conditions l’équilibre thermodynamique aussi bien pour les
grandeurs hémisphériques que directionnelles pour une longueur d’onde donné.

II. 8. Corps gris

On appelle corps gris tout corps dont le coefficient d’absorption est indépendant de la longueur
d’onde .

𝑎 = 𝑎 ∀ 

Dans le cas des surfaces grises on a :

𝑀
= 𝑎 ∀
𝑀0

Nous avons donc : 𝜀 = 𝑎 = 𝜀 𝑐𝑎𝑟 𝑙𝑎 𝑙𝑜𝑖 𝑑𝑒 𝐾𝑖𝑟𝑐𝑕𝑜𝑓𝑓 𝑎 = 𝜀

Dans la pratique un grand nombre de matériaux peuvent être considérés comme des corps gris,
on pet donc prendre 𝑎 = 𝜀.

II. 9. Modèles des « corps gris » et des corps « gris par morceaux »

II.8.1. Corps gris de rayonnement uniforme et isotrope

L’émissivité spectrale est indépendante de la longueur d’onde :

𝜀 = 𝜀 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒

La loi de Kirchhoff s’étend à tout le spectre 𝑎 = 𝜀 . Le corps noir est un cas particulier du corps
gris 𝑎 = 𝜀 = 1 , c’est donc à la fois l’émetteur et le récepteur idéal.

II.8.1. Corps gris par morceau, de rayonnement uniforme et isotrope

Le spectre peut être divisé en différents domaines spectraux à l’intérieur desquels l’émissivité
spectrale est constante :

𝜀 = 𝜀0 𝑝𝑜𝑢𝑟 0 <  < 1

𝜀 = 𝜀1 𝑝𝑜𝑢𝑟 1 <  < 2

……………………………

𝜀 = 𝜀𝑗 𝑝𝑜𝑢𝑟 𝑗 <  < ∞

L’émissivité totale hémisphérique s’exprime en fonction de la fraction d’existence du corps noir

𝐹(𝑥).
90
 ∞
𝑀 𝜀 𝑀0 𝑑
0  
𝜖= 0= ∞ 0
𝑀
0
𝑀 𝑑

1 2 0 ∞
𝑀 𝑀0 𝑑 𝑀 𝑑 𝑗
𝑀0 𝑑
0 1
𝜖 = 0 = 𝜀0 + 𝜀0 +∙∙∙∙∙∙∙∙∙∙∙∙∙∙ +𝜀𝑗
𝑀 𝜍𝑇 4 𝜍𝑇 4 𝜍𝑇 4

𝜖 = 𝜀0 𝐹 𝑥1 + 𝜀1 𝐹 𝑥2 − 𝐹 𝑥1 +∙∙∙∙∙∙∙∙∙∙∙∙∙∙ +𝜀𝑗 1 − 𝐹 𝑥𝑗

Avec

0
𝑀0 𝑑
𝐹 𝑥 =
𝜍𝑇 4

Figure 8. Evolution de l’émittance en fonction de la longueur d’onde et de la température

pour corps noir, gris et réel

III. Echanges radiatifs entre surfaces séparées par un milieu parfaitement transparents

III. 1. Echanges radiatifs entre surfaces noires.

On considère deux surfaces 𝑆1 et 𝑆2 maintenues à des températures 𝑇1 et 𝑇2 (uniformes) et
pouvant rayonner mutuellement l’une vers l’autre.

Le flux totala hémisphérique émis par 𝑆1 est :

1 = M10 S1 dont seule une fraction 12 atteint 𝑆2 :

12 = F12 1 = F12 M10 S1

F12 est un nombre sans dimension représentant la fraction du flux total hémisphérique de la
surface 𝑆1 qui atteint 𝑆2 .

F12 est appelé facteur de forme ou facteur d’angle de 𝑆1 vers 𝑆2 .

91
III. 1. 1. Facteurs de forme.
On considère deux éléments 𝑑𝑆1 et 𝑑𝑆2 prises respectivement sur 𝑆1 et 𝑆2 et écrivons l’expression
du flux émis par 𝑑𝑆1 en direction de 𝑑𝑆2 .
Ce flux est contenu dans l’angle solide élémentaire 𝑑𝑆1 sous lequel 𝑑𝑆2 est vu à partir de 𝑑𝑆1 .

𝑑2 12 = 𝐿01 𝑑𝑆1 𝑐𝑜𝑠𝜃1 𝑑12

Or

𝑑𝑆2 𝑐𝑜𝑠𝜃2
𝑑12 =
𝑙2

Et

𝑀10
𝐿01 = (𝐿𝑜𝑖 𝑑𝑒 𝐿𝑎𝑚𝑏𝑒𝑟𝑡)
𝜋

𝑀10 𝑑𝑆1 𝑐𝑜𝑠𝜃1 𝑑𝑆2 𝑐𝑜𝑠𝜃2

𝑑2 12 =
𝜋 𝑙2

Par intégration de 𝑑2 12 sur 𝑆2 on obtient le flux 𝑑12 émis par 𝑑𝑆1 sous lequel 𝑆2 :

𝑑12 = 𝑑 2 12
𝑆2

En intégrant 𝑑12 sur 𝑆1 , on obtient le flux 12 émis par 𝑆1 sur 𝑆2 donc :

𝑑𝑆1 𝑐𝑜𝑠𝜃1 𝑑𝑆2 𝑐𝑜𝑠𝜃2

12 = 𝑑2 12 = 𝑀10
𝑆1 𝑆2 𝑆1 𝑆2 𝜋 𝑙2

1 𝑑𝑆1 𝑐𝑜𝑠𝜃1 𝑑𝑆2 𝑐𝑜𝑠𝜃2

12 = 𝑀10 𝑆1
𝑆1 𝑆1 𝑆2 𝜋 𝑙2

12 = 𝑀10 𝑆1 𝐹12

Avec :

1 𝑑𝑆1 𝑐𝑜𝑠𝜃1 𝑑𝑆2 𝑐𝑜𝑠𝜃2

𝐹12 =
𝑆1 𝑆1 𝑆2 𝜋 𝑙2

92
 12 = 1 𝐹12

On remarque que le facteur 𝐹12 est une quantité purement géométrique.

Le flux émis par 𝑆2 est atteignant 𝑆1 s’écrit de manière similaire :

21 = 𝑀20 𝑆2 𝐹21

1 𝑑𝑆1 𝑐𝑜𝑠𝜃1 𝑑𝑆2 𝑐𝑜𝑠𝜃2

𝐹21 =
𝑆2 𝑆1 𝑆2 𝜋 𝑙2

𝐹21 est la fraction du flux hémisphérique de 𝑆2 qui atteint 𝑆1 .

𝑑𝑆1 𝑐𝑜𝑠𝜃1 𝑑𝑆2 𝑐𝑜𝑠𝜃2

𝑆1 𝐹12 = 𝑆2 𝐹21 =
𝑆1 𝑆2 𝜋 𝑙2

Donc :

12 = 𝑀10 𝑆1 𝐹12 = 𝑀10 𝑆2 𝐹21

Et

21 = 𝑀20 𝑆2 𝐹21 = 𝑀20 𝑆1 𝐹12

Nous avons le choix de calculer soit 𝐹12 ou 𝐹21 . On choisit généralement le plus simple à calculer.

III. 1. 2. Puissance nette échangée par rayonnement entre S1 et S2

Les deux surfaces étant noires, tout flux reçu par l’une ou l’autre est entièrement absorbé (il n’y a
pas de réflexion, ni de transmission).

On obtient le flux net échangée entre 𝑆1 et 𝑆2 noté 12𝑛𝑒𝑡 .

En faisant sur 𝑆1 un bilan entre la puissance émise par 𝑆1 vers 𝑆2 (il s’agit d’une perte d’énergie
pour 𝑆1 ) et la puissance absorbée par 𝑆1 en provenance de 𝑆2 (il s’agit d’un gain d’énergie pour
𝑆1 ) :

12𝑛𝑒𝑡 = 12 − 21 = 𝑀10 𝑆1 𝐹12 − 𝑀20 𝑆2 𝐹21

12𝑛𝑒𝑡 = 𝑆1 𝐹12 𝑀10 − 𝑀20 = 𝑆2 𝐹21 𝑀10 − 𝑀20

Ou encore

12𝑛𝑒𝑡 = 𝑆1 𝐹12 𝜍 𝑇14 − 𝑇24 = 𝑆2 𝐹21 𝜍 𝑇14 − 𝑇24

Si 𝑇1 > 𝑇2 le flux échangé de 𝑆1 à 𝑆2 est positif, cela veut dire que 𝑆1 émet davantage de
rayonnement vers 𝑆2 qu’elle en absorbe en provenance de cette dernière.

93
Avec cette convention : on compte positivement les pertes d’énergie d’une surface par
rayonnement.

III. 1. 3. Relations algébriques entre facteurs de forme

Nous avons vu que pour deux surfaces 𝑆𝑖 à 𝑆𝑗 , on a la relation

appelée règle de réciprocité

𝑆𝑖 𝐹𝑖𝑗 = 𝑆𝑗 𝐹𝑗𝑖

Considérons une enceinte fermée constituée de 𝑛 surfaces

noires individuellement isothermes.

Pour la 𝑖 𝑖è𝑚𝑒 surface S, on peut définir 𝑛 facteurs de formes :

𝐹𝑖1 , 𝐹𝑖2 ∙∙∙∙∙∙∙∙∙∙∙∙∙ 𝐹𝑖𝑖 , 𝐹𝑖𝑗 ∙∙∙∙∙∙∙∙∙∙ 𝐹𝑖𝑛

Soit 𝐹𝑖𝑗 𝑎𝑣𝑒𝑐 𝑗 = 1, … … … … 𝑛

Avec

𝑖𝑗
𝐹𝑖𝑗 =
𝑖

𝐹𝑖𝑗 𝑓𝑟𝑎𝑐𝑡𝑖𝑜𝑛 𝑑𝑢 𝑓𝑙𝑢𝑥 𝑑𝑒 𝑙𝑎 𝑠𝑢𝑟𝑓𝑎𝑐𝑒 𝑆𝑖 qui atteint 𝑆𝑗 .

𝐹𝑖𝑖 est le facteur de forme de l’échange de 𝑆𝑖 avec elle-même. Il n’existe que si 𝑆𝑖 est concave.

Le flux total 𝑖 émis par 𝑆𝑖 est absorbé par toutes les surfaces constituant l’enceinte y compris 𝑆𝑖
si elle est concave.
n n n

1 = ij = i Fij = i Fij

j=1 j=1 j=1

D’où :
n

Fij = 1
j=1

Intérêt pratique

Exemple :

Soient deux surfaces 𝑆𝑖 à 𝑆𝑗 telles que :

𝑆𝑗 = 𝑆𝑗 1 + 𝑆𝑗 2

94
 ij = ij1 + ij2

ij = Fij i = Fij1 i + Fij2 i

Fij = Fi j1+j2 = Fij1 + Fij2

𝑆𝑖 Fij = 𝑆𝑖 Fij1 + 𝑆𝑖 Fij2 (3)

Compte tenu des relations de réciprocité :

𝑆𝑖 𝐹𝑖𝑗 = 𝑆𝑗 𝐹𝑗𝑖

𝑆𝑖 𝐹𝑖𝑗 1 = 𝑆𝑗 1 𝐹𝑗 1 𝑖

𝑆𝑖 𝐹𝑖𝑗 2 = 𝑆𝑗 2 𝐹𝑗 2 𝑖

On remplace dans l’équation (3)

𝑆𝑗 𝐹𝑗𝑖 = 𝑆𝑗 1 𝐹𝑗 1 𝑖 + 𝑆𝑗 2 𝐹𝑗 2 𝑖

Dans certains cas on a : 𝑆𝑗 1 𝐹𝑗 1 𝑖 = 𝑆𝑗 2 𝐹𝑗 2 𝑖 , nous avons donc :

𝑆𝑗 𝐹𝑗𝑖
𝑆𝑗 1 𝐹𝑗 1 𝑖 = 𝑆𝑗 2 𝐹𝑗 2 𝑖 =
2

Ceci constitue un moyen pour calculer les facteurs de forme 𝐹𝑗 1 𝑖 et 𝐹𝑗 2 𝑖 connaissant le facteur de
forme 𝐹𝑗𝑖 (ou 𝐹𝑖𝑗 ).

III. 1. 4. Evaluation des facteurs de forme

On peut évaluer les facteurs de forme de trois manières :

1) Calcul d’une intégrale de surface (calcul long).

2) Détermination à partir de considérations simples de géométrie.
3) Il existe des abaques ou de tables de facteurs de forme pour différentes configurations
courantes (les autres sont déduits par les relations algébriques)

Exemples

a) 𝐹12 = 𝐹21

95
Pour toutes les surfaces parallèles de grandes dimensions par rapport à leur distance :

- Plans parallèles
- Sphères concentriques
- Cylindres coaxiaux très longs

Ces surfaces sont dîtes en influence totale (tout le flux de l’une atteint l’autre).

b) Demi sphère ou demi cylindre très long par rapport à leur plan de base 𝑆1 .

Tout le flux émis par 𝑆1 atteint 𝑆2

𝐹12 = 1

𝑆1 𝐹12 = 𝑆2 𝐹21

𝑆1
𝐹21 =
𝑆2

- Pour la demi-sphère
𝜋𝑅 2 1
𝐹21 = 2
=
2𝜋𝑅 2
- Pour le demi cylindre

2𝑅𝑙 2
𝐹21 = =
𝜋𝑅𝑙 𝜋

c) Pour toute surface convexe 𝑆1 située à l’intérieur d’une surface concave 𝑆2

𝐹12 = 1

𝑆1 𝐹12 = 𝑆2 𝐹21

𝑆1
𝐹21 =
𝑆2

- Cas des deux sphères concentriques 𝑆1 étant la plus petite

𝐹12 = 1
2
𝑅1
𝐹21 =
𝑅2

Puisque 𝐹21 + 𝐹22 = 1, nous avons :

96
 2
𝑅1
𝐹22 = 1 − 𝐹21 = 1 −
𝑅2

𝐹22 est la fraction du flux émis par la surface concave, qui revient directement sur celle-ci.

Si 𝑅1 𝑡𝑒𝑛𝑑 𝑣𝑒𝑟𝑠 𝑅2 on retrouve le premier cas aussi bien pour les cylindres que pour les sphères.

Application

Dans un four parallélépipédique (voir figure), on assimile la voute ( 𝑆1 ), la sole ( 𝑆2 ) et les parois
latérales (𝑆3 ) à des surfaces noires de températures uniformes respectives 𝑇1 , 𝑇2 , 𝑇3 . La surface
𝑆3 regroupe les surfaces 𝑆4 , 𝑆5 , 𝑆6 𝑒𝑡 𝑆7 . La porte a pour surface 𝑆0 et pour température 𝑇0 . ( 𝑇0
est la température ambiante de la salle ou se trouve le four).

Calculer les pertes radiatives par la porte d’enfournement lorsque celle-ci est ouverte pendant le
fonctionnement ? La porte à pour surface 𝑆0 à la température 𝑇0

On donne 𝑇1 = 800𝐾, 𝑇2 = 1200𝐾, 𝑇3 = 1000𝐾.

Le flux net échangé entre 𝑆0 et l’une des trois surfaces 𝑆1 , 𝑆2 , 𝑜𝑢 𝑆3 a pour expression :

0𝑖𝑛𝑒𝑡 = 𝑆0 𝐹0𝑖 𝑀00 − 𝑀𝑖0 𝑎𝑣𝑒𝑐 𝑖 = 1, 2, 3

Le flux échangé entre 𝑆0 et l’ensemble de parois intérieures du four est :

3

0𝑖𝑛𝑒𝑡 = 𝑆0 𝐹0𝑖 𝑀00 − 𝑀𝑖0

𝑖=1

0𝑖𝑛𝑒𝑡 = 𝑆0 𝐹01 𝜍 𝑇04 − 𝑇14 + 𝑆0 𝐹02 𝜍 𝑇04 − 𝑇24 + 𝑆0 𝐹03 𝜍 𝑇04 − 𝑇34

97
Evaluation des facteurs de forme :

On utilise les abaques

 Calcul de 𝐹02 ?

𝑎 = 4𝑚, 𝑏 = 1𝑚, 𝑐 = 3𝑚

D’après l’abaque figure 𝐹15 = 0.135

𝐹02 = 0.35

𝑆0 𝐹02 = 1.05 𝑚2

 Calcul de 𝑆0 𝐹01 ?

La surface 𝑆3 est composée des surfaces 𝑆4 , 𝑆5 , 𝑆6 𝑒𝑡 𝑆7 , on a :

𝐹10 = 𝐹1 7+1 − 𝐹17 = 𝐹15 − 𝐹17

𝐹15 ?

𝑎 = 2𝑚, 𝑏 = 4𝑚, 𝑐 = 3𝑚

D’après l’abaque figure 𝐹15 = 0.135

𝐹17 ?

𝑎 = 1𝑚, 𝑏 = 4𝑚, 𝑐 = 3𝑚

D’après l’abaque figure 𝐹17 = 0.0875

Donc 𝐹10 = 0.0475

𝑆0 𝐹01 = 0.57𝑚2

 Calcul de 𝑆0 𝐹03 ?

𝑆0 𝐹03 = 𝑆0 𝐹04 + 𝑆0 𝐹05 + 𝑆0 𝐹06 + 𝑆0 𝐹07

𝑆0 𝐹04 = 𝑆4 𝐹40 = 𝑆4 𝐹47 𝑐𝑎𝑟 𝐹40 = 𝐹47 (𝑝𝑎𝑟 𝑠𝑦𝑚é𝑡𝑟𝑖𝑒)

𝑆4 𝐹4 7+0 = 𝑆4 𝐹47 + 𝑆4 𝐹40 = 2 𝑆4 𝐹40

Donc

𝑆4 𝐹4 7+0 𝑆4 𝐹45
𝑆0 𝐹04 = =
2 2

  Calcul de 𝐹45 ?

98
𝑎 = 3𝑚, 𝑏 = 4𝑚, 𝑐 = 2𝑚

En lisant sur l’abaque de la figure 2, on trouve

𝐹45 = 0.135

𝑆0 𝐹04 = 0.54 𝑚2

 Calcul de 𝑆0 𝐹05 ?

𝑆0 𝐹05 = 𝑆5 𝐹50 = 𝑆5 𝐹57 𝑐𝑎𝑟 𝐹50 = 𝐹57 (𝑝𝑎𝑟 𝑠𝑦𝑚é𝑡𝑟𝑖𝑒)

𝑆5 𝐹5 7+0 = 𝑆5 𝐹57 + 𝑆5 𝐹50 = 2 𝑆5 𝐹50

𝑆5 𝐹5 7+0
𝑆0 𝐹05 = 𝑆5 𝐹50 =
2

𝑆5 𝐹5 7+0 𝑒𝑠𝑡 𝑓𝑜𝑢𝑟𝑛𝑖 𝑝𝑎𝑟 𝑙 ′ 𝑎𝑏𝑎𝑞𝑢𝑒𝑑𝑒 𝑙𝑎 𝑓𝑖𝑔𝑢𝑟𝑒 ? ? ? ?

𝑎 = 2𝑚, 𝑏 = 3𝑚, 𝑐 = 4𝑚

𝐹5 7+0 = 0.097

D’où 𝑆0 𝐹05 = 0.29 𝑚2

Finalement

𝑆0 𝐹03 = 𝑆0 𝐹04 + 𝑆0 𝐹05 + 𝑆0 𝐹06 + 𝑆0 𝐹07

𝑆0 𝐹03 = 1.38 𝑚2

On aurait pu calculer 𝐹03 directement en utilisant la relation pour une surface fermée :

𝐹01 + 𝐹02 + 𝐹03 + 𝐹07 = 1

𝐹03 = 1 − 𝐹01 + 𝐹02 + 𝐹07

Le flux net est

0𝑖𝑛𝑒𝑡 = 𝑆0 𝐹01 𝜍 𝑇04 − 𝑇14 + 𝑆0 𝐹02 𝜍 𝑇04 − 𝑇24 + 𝑆0 𝐹03 𝜍 𝑇04 − 𝑇34

0𝑖𝑛𝑒𝑡 = 5.67 10−8 0.57 3004 − 12004 + 1.05 3004 − 8004 + 1.38 3004 − 10004

0𝑖𝑛𝑒𝑡 = −168𝑘𝑊

Ce flux est négatif, il s’agit en effet d’un flux gagné par 𝑆0 , donc ceci est conforme avec les
conventions de signe adoptées dans ce cours.

99
III. 2. Echanges radiatifs entre surfaces opaques grises séparées par un milieu parfaitement
transparent.

Nous avons étudié les surfaces dans le paragraphe précédent. Ceci n’est réaliste que si le
coefficient d’émission 𝜀 > 0.9. Nous nous intéressons dans ce paragraphe aux surfaces grises et
diffusantes en émission. Leur réflexions sera soit diffuse soit spéculaire.

III. 1. 1. Surfaces opaques grises et diffusantes en émission et en réflexion.

Pour faire cette étude, nous adoptons les hypothèses

suivantes ;

Chaque surface bénéficie des propretés suivantes :

- L’isothermie
- Propriétés radiatives indépendantes de la longueur
d’onde et de la direction.
- Pour chaque surface on a : 𝑎 = 𝜀 = 1 − 𝑟.
- E flux incident sur toute surface y est uniformément répart.

𝑖 flux incident

𝑎 flux absorbé

𝑟 flux réfléchi

L’éclairement de la surface est :

𝑖
𝐸=
𝑆

L’émittance de la surface est :

𝑀 = 𝜀 𝑀0

On introduit une nouvelle grandeur la radiosité notée 𝑱, qui est la somme de l’émittance de la
surface et du flux réfléchi par l’unité de surface qui est 𝑟 𝐸. La radiosité s’écrit :

𝐽 = 𝜀 𝑀0 + 𝑟 𝐸

Pour une surface opaque, le coefficient de transmission est nul

𝑡 = 0 et le coefficient de réflexion s’écrit :

𝑟 =1−𝑎 = 1− 𝜀

La relation de la radiosité devient :

𝐽 = 𝜀 𝑀0 + 1 − 𝜀 𝐸

Le flux perdu par la surface S sera égal à la différence entre le flux émis et absorbé par cette
surface, soit :
100
 𝑛𝑒𝑡 = 𝜀 𝑀0 − 𝑎𝐸 𝑆

Comme 𝑎 = 𝜀

𝑛𝑒𝑡 = 𝜀 𝑆 𝑀0 − 𝐸

En élimant E de cette relation on obtient :

𝜀𝑆
𝑛𝑒𝑡 = 𝑀0 − 𝐽
1−𝜀

Ou

𝑛𝑒𝑡
= 𝑀0 − 𝐽
𝑆

Ainsi la densité de flux nette perdue par une surface est égale à la différence entre sa radiosité
(densité de flux quittant la surface par émission et réflexion) et son éclairement (densité de flux
arrivant sur la surface).

III. 1. 2. Echanges dans une enceinte vide aux parois grises et diffusantes en émission et
réflexion.

Considérons une enceinte fermée constituée de 𝑛 surfaces grises

individuellement isothermes 𝑆𝑖 𝑎𝑣𝑒𝑐 𝑖 = 1, … … … … 𝑛

La radiosité de 𝑆𝑖 est égal à la somme de l’émittance propre et du

flux en provenance de toutes les surfaces de l’enceinte et réfléchi
par l’unité de surface de 𝑆𝑖 .

Radiosité de 𝑆𝑖 = 𝐽𝑖

Emittance propre de 𝑆𝑖 = 𝜀𝑖 𝑀𝑖0

Flux incident sur 𝑆𝑖 :

𝑛 𝑛

𝑆𝑗 𝐹𝑗𝑖 𝐽𝑗 = 𝑆𝑖 𝐹𝑖𝑗 𝐽𝑗
𝑗 =1 𝑗 =1

Flux incident sur l »unité de surface de 𝑆𝑖 :

𝑛

𝐹𝑖𝑗 𝐽𝑗
𝑗 =1

Flux réfléchi par l’unité de surface de 𝑆𝑖 :

𝑛 𝑛

𝑟𝑖 𝐹𝑖𝑗 𝐽𝑗 = 1 − 𝜀𝑖 𝐹𝑖𝑗 𝐽𝑗
𝑗 =1 𝑗 =1

101
D’où
𝑛

𝐽𝑖 = 𝜀𝑖 𝑀𝑖0 + 1 − 𝜀𝑖 𝐹𝑖𝑗 𝐽𝑗
𝑗 =1

𝐽𝑖 − 1 − 𝜀𝑖 𝐹𝑖𝑗 𝐽𝑗 = 𝜀𝑖 𝑀𝑖0 = 𝜀𝑖 𝜍 𝑇𝑖4

𝑗 =1

Cette équation est à utiliser pour toutes les surfaces 𝑆𝑖 de l’enceinte ayant une température 𝑇𝑖
connue ou imposée.

Pour les surfaces ayant un flux imposé, on écrit à nouveau l’équation de bilan :

Le flux net perdu par l’unité de surface est égal à la différence entre le flux quittant l’unité de
surface 𝑆𝑖 (radiosité) et le flux arrivant sur l’unité de surface 𝑆𝑖 .
𝑛
𝑖𝑛𝑒𝑡
𝜑𝑖 = = 𝐽𝑖 − 𝐹𝑖𝑗 𝐽𝑗
𝑆
𝑗 =1

Application

Dans un four électrique à chambre parallélépipédique (voir

figure). La surface ( 𝑆1 ), est à la température 𝑇1 . La surface
( 𝑆2 ) est à la température 𝑇2 . La surface (𝑆3 ) comprend les
quatre parois latérales adiabatiques à la température 𝑇3 .

Calculer les flux radiatifs 1𝑛𝑒𝑡 , 2𝑛𝑒𝑡 , 𝑇3 .

On donne 𝑇1 = 800𝐾, 𝑇2 = 1200𝐾, 𝜀1 = 0.8, 𝜀2 = 0.9

Solution

Pour les surfaces 𝑆1 𝑒𝑡 𝑆2 qui ont des températures imposées on a les équations :

𝐽1 − 1 − 𝜀1 𝐹11 𝐽1 + 𝐹12 𝐽2 + 𝐹13 𝐽3 = 𝜀1 𝜍 𝑇14 (𝑎)

𝐽2 − 1 − 𝜀2 𝐹21 𝐽1 + 𝐹22 𝐽2 + 𝐹23 𝐽3 = 𝜀2 𝜍 𝑇44 (𝑏)

Pour la surface 𝑆3 , nous avons un flux imposé (comme la surface adiabatique le flux net imposé
est nul) :

3𝑛𝑒𝑡
= 0 = 𝐽3 − 𝐹31 𝐽1 + 𝐹32 𝐽2 + 𝐹33 𝐽3 𝑐
𝑆

Nous avons trois équations à trois inconnues 𝐽1 , 𝐽2 𝑒𝑡 𝐽3 .

102
Evaluation des facteurs de forme :

𝐹32 = 𝐹31 (𝑝𝑎𝑟 𝑠𝑦𝑚é𝑡𝑟𝑖𝑒)

La surface étudiée est une surface fermée :

n

Fij = 1
j=1

𝐹31 + 𝐹32 + 𝐹33 = 1

𝐹33 = 1 − 𝐹31 + 𝐹32 = 2 𝐹31

L’équation (c) devient :

𝐽3 − 𝐹31 𝐽1 − 𝐹32 𝐽2 – 1 − 2 𝐹31 𝐽3 = 0

D’où

𝐽1 + 𝐽2
𝐽3 =
2
En plus nous avons :

𝐹13 = 𝐹23

𝐹12 = 𝐹21

En utilisant l’abaque (figure ???)

𝑎 = 4𝑚, 𝑏 = 1𝑚, 𝑐 = 1𝑚

𝐹12 = 0.34

Les équations (a) et (b) s’écrivent :

2 1 − 𝜀1 1 − 2 𝐹12 𝐽1 − 1 − 𝜀1 1 + 𝐹12 𝐽2 = 2 𝜀1 𝜍 𝑇14

− 1 − 𝜀2 1 + 𝐹12 𝐽1 + 2 − 1 − 𝜀2 1 − 𝐹12 𝐽2 = 2 𝜀2 𝜍 𝑇24

On obtient le système :

1.868 𝐽1 − 0.268 𝐽2 = 3.7159 104

−0.134 𝐽1 + 1.934 𝐽2 = 21.1631 104

𝐽1 = 11.19 104 𝑊/𝑚2

𝐽2 = 3.59 104 𝑊/𝑚2

𝐽3 = 7.39 104 𝑊/𝑚2

Calcul des flux :

𝜀1 𝑆1
1𝑛𝑒𝑡 = 𝜍 𝑇14 − 𝐽1 = −204.6 𝑘𝑊
1 − 𝜀1

103
 𝜀2 𝑆2
2𝑛𝑒𝑡 = 𝜍 𝑇24 − 𝐽2 = +204.6 𝑘𝑊
1 − 𝜀2

Ce résultat étais prévisible puisque la surface 𝑆3 est adiabatique, par conséquent tout le flux cédé
par 𝑆2 est gagné par 𝑆1 .

Calcul des flux :

𝜀3 𝑆3
3𝑛𝑒𝑡 = 𝜍 𝑇34 − 𝐽3 = 0
1 − 𝜀3

𝜍 𝑇34 = 𝐽3
1
𝐽3 4
𝑇3 = = 1068 𝐾
𝜍

La température 𝑇3 a une valeur intermédiaire entre 𝑇1 𝑒𝑡 𝑇2 , ce qui est normal.

Fig. Exemples de courbes des facteurs de formes

104
Application

Deux plans indéfinis et parallèles P1 et P2 dont la distance les séparant est petite devant leurs
dimensions sont maintenus respectivement aux températures T1 et T2. Supposons que les deux
plans sont des corps gris de facteurs d’émissions 1 et  2 ; ils sont séparés par de l’air dont on
néglige les mouvements de convection.
1°) Déterminer par les deux méthodes ci-dessous le flux de chaleur transmis d’une surface à
l’autre par rayonnement.

1ière méthode : Choisir une densité de flux 1 qui part de P1 et une densité de flux 2 qui part de
P2 et qui vont subir une suite de réflexions sur les plans P1 et P2 (Figure 2) 2ième méthode : En
faisant un bilan entre les radiosités J1 et J2 (respectivement) de P1 et P2 et le flux de chaleur
échangé entre les deux plans.

Solution :

Le plan 𝑃1 émet
1 = 𝜀1 𝜍 𝑇14

La totalité du flux atteint le plan 𝑃2

La plan 𝑃2 considéré comme récepteur reçoit une partie

du rayonnement incident égale à 𝜀2 1 (le facteur
d’émission est égale au facteur d’absorption) et réfléchit
le reste du rayonnement soit 1 − 𝜀2 1 .

L’expression du flux absorbé à partir du flux qu’il a lui-

même émis est :

𝜀1 1 − 𝜀2 1 + 𝜀1 1 − 𝜀1 1 − 𝜀2 2 1 +∙∙∙∙∙∙∙
∙ +𝜀1 1 − 𝜀1 𝑛
1 − 𝜀2 𝑛+1
1

= 𝜀1 1 − 𝜀2 1 1 + 1 − 𝜀1 1 − 𝜀2 +∙∙∙∙∙∙∙
𝑛 𝑛
∙ + 1 − 𝜀1 1 − 𝜀2

1 − 1 − 𝜀1 𝑛 1 − 𝜀2 𝑛
= 𝜀1 1 − 𝜀2 1
1 − 1 − 𝜀1 1 − 𝜀2

𝜀1 1 − 𝜀2 1
=
1 − 1 − 𝜀1 1 − 𝜀2
𝑛 𝑛
𝐶𝑎𝑟, 𝑞𝑢𝑎𝑛𝑑 𝑛 → ∞, 𝑙𝑒 𝑡𝑒𝑟𝑚𝑒 1 − 𝜀1 1 − 𝜀2 → 0

Le flux absorbé par 𝑃1 provenant du flux émis par 𝑃2 est :

𝜀1 2
′1 =
1 − 1 − 𝜀1 1 − 𝜀2

Le flux absorbé par 𝑃2 à partir du flux émis par 𝑃1 est :

105
 𝜀2 1
2 =
1 − 1 − 𝜀1 1 − 𝜀2

Le flux absorbé par 𝑃2 provenant du flux émis par 𝑃2 est :

𝜀2 1 − 𝜀1 2
′2 =
1 − 1 − 𝜀1 1 − 𝜀2

Bilan sur le plan 𝑷𝟏

Au total 𝑃1 a donc émis 1 = ε1 σ T14 et absorbé la somme :

𝜀1 1 − 𝜀2 1 𝜀1 2
+
1 − 1 − 𝜀1 1 − 𝜀2 1 − 1 − 𝜀1 1 − 𝜀2

𝜀1 1 − 𝜀2 1 + 𝜀1 2
=
1 − 1 − 𝜀1 1 − 𝜀2

Le flux de chaleur net perdu par 𝑃1 , c'est-à-dire échangé avec 𝑃2 est la différence du flux émis et
du flux absorbé.

𝜀1 1 − 𝜀2 1 + 𝜀1 2
1𝑛𝑒𝑡 = 1 −
1 − 1 − 𝜀1 1 − 𝜀2

𝜀2 1 − 𝜀1 2
1𝑛𝑒𝑡 =
1 − 1 − 𝜀1 1 − 𝜀2
𝜀 1 𝜀2
1𝑛𝑒𝑡 = 1− 1−𝜀 1 1−𝜀 2
𝜍𝑆 𝑇14 − 𝑇24

Méthode basée sur la radiosité

𝐽1 = 𝜑1 + 1 − 𝜀1 𝐽2

𝐽2 = 𝜑2 + 1 − 𝜀2 𝐽1

𝜑2 + 1 − 𝜀2 𝜑1
𝐽2 =
1 − 1 − 𝜀1 1 − 𝜀2

𝜑1 + 1 − 𝜀1 𝜑2
𝐽1 =
1 − 1 − 𝜀1 1 − 𝜀2

Le flux échangé entre les deux surfaces est :

𝜑1 + 1 − 𝜀1 𝜑2 𝜑2 + 1 − 𝜀2 𝜑1
𝜑12 = 𝐽1 − 𝐽2 = −
1 − 1 − 𝜀1 1 − 𝜀2 1 − 1 − 𝜀1 1 − 𝜀2
𝜀2 𝜑1 − 𝜀1 𝜑2
𝜑12 =
1 − 1 − 𝜀1 1 − 𝜀2
𝜀1 𝜀2
𝜑12 = 𝜍 𝑇14 − 𝜍𝑇24
1 − 1 − 𝜀1 1 − 𝜀2
106
Energie solaire reçue par la terre

Supposons que le soleil rayonne comme un corps noir à la température 𝑇 = 5800°𝐶, on se

propose de calculer l’éclairement solaire 𝐸 reçu par mètre carré de surface de terre. On donne :

Le rayon du soleil : 𝑅𝑆 = 696000 𝑘𝑚

La distance Terre-Soleil : 𝐷 = 149637000 𝑘𝑚

On suppose que le coefficient de transmission de l’atmosphère est égal à 𝜏 = 0.75

Solution

Soit 𝑑𝑆 un élément de surface de la surface du soleil.

Evaluons sa contribution à l’éclairement de l’élément de surface 𝑑𝑆′ considéré sur la terre.

𝑑 : angle solide sous lequel on voit 𝑑𝑆′ de la terre, depuis 𝑑𝑆 sur le soleil.

Le flux élémentaire émis par 𝑑𝑆 dans un angle solide 𝑑 entourant la direction de 𝑑𝑆 vers 𝑑𝑆′ est :

𝑑2  𝑑𝑆→𝑑𝑆′ = 𝐿0𝑇 𝑑𝑆 𝑐𝑜𝑠𝜃 𝑑𝑑𝑆→𝑑𝑆′

𝑑𝑆′𝑐𝑜𝑠𝜃′
𝑑𝑑𝑆→𝑑𝑆′ =
𝑑2

𝑀𝑇0 𝑑𝑆 𝑐𝑜𝑠𝜃 𝑑𝑆′𝑐𝑜𝑠𝜃′

𝑑2  𝑑𝑆→𝑑𝑆′ =
𝜋 𝑑2

On intègre sur toute la surface du soleil :

𝑀𝑇0 𝑑𝑆 𝑐𝑜𝑠𝜃 𝑑𝑆′𝑐𝑜𝑠𝜃′

𝑑 𝑆→𝑑𝑆′ =
𝜋 𝑑2

𝐷 ≅ 150 𝑅  𝑑 ≅ 𝐷 et 𝜃′est négligeable quelque soit la position de 𝑑𝑆 ′ .

𝑀𝑇0 𝑑𝑆′ 𝑀𝑇0 𝑑𝑆′

𝑑  𝑆→𝑑𝑆′ = 𝑑𝑆 𝑐𝑜𝑠𝜃 = 𝑑
𝜋 𝐷2 𝜋 𝐷2

Avec :

𝑑 =  = 𝜋 𝑅 2

𝑀𝑇0 𝑑𝑆′
𝑑 𝑆→𝑑𝑆′ = 𝜋 𝑅2
𝜋 𝐷2
2
𝑅
𝑑 𝑆→𝑑𝑆′ = 𝑀𝑇0 𝑑𝑆′
𝐷

𝑑 𝑆→𝑑𝑆′ 𝑅 2
𝐸= = 𝑀𝑇0
𝑑𝑆′ 𝐷
2
𝑅
𝐸 = 𝜍 𝑇4
𝐷

107
Avec :

𝑊
𝜍 = 5.67 10−8 .
𝑚2 𝐾 4

𝑇 = 5800°𝐶.

𝑅𝑆 = 696000 𝑘𝑚

𝐷 = 149637000 𝑘𝑚

L’éclairement est donc : 𝐸 = 1402 𝑊/𝑚2

Si on tient compte de l’atmosphère : 𝜏 𝐸 = 1052 𝑊/𝑚2 . En réalité la valeur qu’on mesure est plus
faible que cette valeur.

108
109
110
111
112
113
114
115
  Corrélations pour le calcul des coefficients de transfert en convection
forcée

116
117
 Corrélations pour le calcul des coefficients de transfert en convection
naturelle

118
Bibliographie

[1] D.P. De witt et F.P. Incropera. Fundamentals of heat and mass transfer. John Wiley and
Sons.
[2] Yunus A. . Cengel and Afshin J. Ghajar. heat and mass transfer. Mc Graw Hill
[3] J. F. Sacadura. Initiation aux transferts thermiques. Edition Tec & Doc .
[4] B. Cheron. Transferts thermiques - Résumé de cours, problèmes corrigés. Edition Ellipses
[5] J. Crabol. Transfert de chaleur Tome 1: Les Principes. Edition Masson
[6] J. Crabol. Transfert de chaleur Tome 2 : Applications industrielles. Edition Masson
[7] J. Crabol. Transfert de chaleur Tome 3 - Corrigés de problèmes Edition Masson
[8] J. Taine, E. Iacona, J.P. Petit. Transferts thermiques Introduction aux transferts d'énergie
Collection. Dunod
[9] Y. Jannot. Transferts thermiques. Ecole des Mines de Nancy.
[Link]/principal/thermique.p

119
120