Thermodynamique

Chapitre 2

Phénomènes de diffusion
PC⋆ , Fabert (Metz)

Phénomènes de diffusion
Dans ce chapitre nous allons nous intéresser à la diffusion. Il s’agit là d’un phénomène important
de la thermodynamique car c’est une source importante d’irréversibilité.
La diffusion se caractérise par :
➜ la nécessité d’un support matériel (contrairement à la propagation ou à l’effusion) ;
➜ l’irréversibilité du phénomène ;
➜ la lenteur du phénomène (nous définirons plus loin en quoi ce phénomène est « lent »).
Exemples de diffusion :
➜ taches d’encre sur un buvard (diffusion de particules) ;
➜ conduction thermique (diffusion d’énergie cinétique) ;
➜ conduction électrique (diffusion d’électrons) ;
➜ viscosité des fluides (diffusion de la quantité de mouvement).
Nous séparerons cette étude en deux. Nous commencerons par nous intéresser à la diffusion
thermique avant de regarder le peu de choses qui changent pour la diffusion de particules.

© Matthieu Rigaut 2 / 88 Version du 28 déc. 2013

PC⋆ , Fabert (Metz) TABLE DES MATIÈRES

Table des matières

Biographies succintes 7

I Diffusion thermique 9
I·1 Première approche des transferts thermiques . . . . . . . . . . . . . . . . . . . . . . . 9
I·1·i modes de transfert . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
la convection . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
le rayonnement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
la conduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
I·1·ii point de vue de la thermodynamique . . . . . . . . . . . . . . . . . . . . . . 10
I·1·iii flux thermique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
résultat . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
justification . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
un dernier mot . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
I·1·iv production énergétique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
phénoménologie de base . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
3 grands types de production . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
traduction formelle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
I·2 Bilan thermique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
I·2·i courant thermique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
l’énergie va se déplacer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
généralisation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
avec la convention usuelle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
simplification usuelle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
I·2·ii bilan 1D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
exemple du barreau calorifugé . . . . . . . . . . . . . . . . . . . . . . . . . . 16
approche mésoscopique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
variation dans le temps . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
échange à travers la surface . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
création en volume . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
rassemblement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
résultat connu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
I·2·iii bilan 3D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
situation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
variation dans le temps . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
production en volume . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
échange à travers la surface . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
rassemblement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
I·2·iv pause divergence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
relation fondamentale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
comment ça marche ? . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
interprétation physique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
expression en coordonnées cartésiennes . . . . . . . . . . . . . . . . . . . . . 25
c’est un opérateur vectoriel différentiel . . . . . . . . . . . . . . . . . . . . . 26
c’est une vision 3D d’un théorème connu . . . . . . . . . . . . . . . . . . . . 26
I·2·v petit flash-back . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
passer de la 3D à la 1D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
passer de la 1D à la 3D sur un cas déjà rencontré . . . . . . . . . . . . . . . 27
I·2·vi régime stationnaire sans production . . . . . . . . . . . . . . . . . . . . . . . 27
© Matthieu Rigaut 3 / 88 Version du 28 déc. 2013
PC⋆ , Fabert (Metz) TABLE DES MATIÈRES

en 1D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
continuité du flux en 3D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
I·3 Équation de diffusion thermique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
I·3·i loi de Fourier . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
énoncé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
interprétation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
limites . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
I·3·ii équation en 1D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
équation de diffusion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
cas particulier fréquent : diffusion sans terme de création . . . . . . . . . . . 30
ordres de grandeur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
I·3·iii interprétation de l’équation de diffusion . . . . . . . . . . . . . . . . . . . . . 31
approche technique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
irréversibilité . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
ordre de grandeur . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
régime permanent stationnaire . . . . . . . . . . . . . . . . . . . . . . . . . . 32
I·3·iv équation de diffusion thermique en 3D . . . . . . . . . . . . . . . . . . . . . . 33
I·3·v pause laplacien et nabla . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
définition du laplacien . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
expression du laplacien en coordonnées cartésiennes . . . . . . . . . . . . . . 33
notation nabla . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
I·3·vi nouveau flash-back . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
passer de la 3D à la 1D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
passer de la 1D à la 3D sur un cas déjà rencontré . . . . . . . . . . . . . . . 36
une notation et un nom . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
I·4 Quelques solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
I·4·i méthode . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
I·4·ii choc thermique – solution en régime transitoire en 1D . . . . . . . . . . . . . 37
situation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
graphiquement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
solution et interprétation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
I·4·iii mur de maison – régime stationnaire 1D . . . . . . . . . . . . . . . . . . . . 38
mur simple . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
mur isolé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
I·4·iv barreau radioactif – régime stationnaire en symétrie cylindrique . . . . . . . 44
dispositif et analyse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
mise en équation : méthode du bilan global . . . . . . . . . . . . . . . . . . . 45
mise en équation : méthode sans réfléchir . . . . . . . . . . . . . . . . . . . . 47
comparaison des deux méthodes . . . . . . . . . . . . . . . . . . . . . . . . . 48
I·4·v température dans le sol – régime sinusoïdal forcé en 1D . . . . . . . . . . . . 48
situation et analyse physique . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
analyse technique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
résolution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
interprétation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
I·4·vi résumé des conditions aux limites . . . . . . . . . . . . . . . . . . . . . . . . 53
condition au bord d’un milieu diffusant . . . . . . . . . . . . . . . . . . . . . 53
condition entre deux milieux diffusants . . . . . . . . . . . . . . . . . . . . . 54
I·5 Résistance thermique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
I·5·i définition par l’exemple . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
tout d’abord les circonstances . . . . . . . . . . . . . . . . . . . . . . . . . . 57
© Matthieu Rigaut 4 / 88 Version du 28 déc. 2013
PC⋆ , Fabert (Metz) TABLE DES MATIÈRES

exemple du mur simple . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57

définition retenue . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
interprétation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
I·5·ii analogie électrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
I·5·iii association de résistances thermiques . . . . . . . . . . . . . . . . . . . . . . 59
association parallèle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
association série . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
exemple du mur isolé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
I·5·iv résistance cylindrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
situation, analyse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
méthode : associer des résistances élémentaires . . . . . . . . . . . . . . . . . 63
méthode : partir de la loi globale . . . . . . . . . . . . . . . . . . . . . . . . . 65
comparaison . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66

II Diffusion de particules 67
II·1 Molécules en mouvement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
II·1·i mouvement brownien . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
II·1·ii exemples de phénomènes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
II·2 Bilan moléculaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
II·2·i modélisation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
densité particulaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
courant particulaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
II·2·ii bilan particulaire 1D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
situation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
approche mésoscopique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
variation dans le temps . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
échange à travers la surface . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
production en volume . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
rassemblement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
II·2·iii bilan particulaire 3D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
approche globale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
variation dans le temps . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
échange à travers la surface . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
production en volume . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
rassemblement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
II·3 Équation de diffusion particulaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
II·3·i loi de Fick . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
énoncé . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
interprétation, limites . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
quelques valeurs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
II·3·ii équation de diffusion à 1D . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
partir de l’équation de continuité . . . . . . . . . . . . . . . . . . . . . . . . . 74
approche qualitative . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
interprétation en terme de mouvement . . . . . . . . . . . . . . . . . . . . . 75
II·3·iii équation de diffusion directement en symétrie cylindrique . . . . . . . . . . . 77
situation et analyse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
un système particulier pour un bilan usuel . . . . . . . . . . . . . . . . . . . 78
II·3·iv équation de diffusion directement en symétrie sphérique . . . . . . . . . . . . 81
situation et analyse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
un système particulier pour un bilan usuel . . . . . . . . . . . . . . . . . . . 81
© Matthieu Rigaut 5 / 88 Version du 28 déc. 2013
PC⋆ , Fabert (Metz) TABLE DES MATIÈRES

II·3·v équation de diffusion en 3D . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83

laplacien en coordonnées cylindriques dans un cas particulier . . . . . . . . . 84
laplacien en coordonnées cylindriques dans un cas particulier . . . . . . . . . 84
II·4 Quelques solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
II·4·i tache d’encre – régime transitoire en 1D . . . . . . . . . . . . . . . . . . . . . 84
situation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
solution, représentation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
II·4·ii canal à ion – régime stationnaire 1D . . . . . . . . . . . . . . . . . . . . . . . 85

Fiche de révision 87

© Matthieu Rigaut 6 / 88 Version du 28 déc. 2013

PC⋆ , Fabert (Metz) Biographies succinctes

Biographies succintes

Pierre Simon marquis de Laplace

(1749 Beaumont en Auge – 1827 Paris)

Né dans une famille aisée, Pierre Simon aurait dû devenir éclésiastique mais il
décide de lui-même de monter à Paris avec une lettre de recommandation. Il
fait rapidement ses preuves et Jean le rond d’Alembert lui trouve un poste
rémunéré de professeur de mathématiques. Pierre Simon est élu à l’académie
royale des sciences à 24 ans et à l’académie française en 1816. Entre temps son
prestige international lui permettra de naviguer dans les hautes sphères aussi bien
durant la monarchie que durant la révolution ou sous Napoléon. Pierre Simon
de Laplace est connu pour le « déterminisme laplacien » pour qui tout le futur
et tout le passé est absolument connu à celui qui connaît les lois physiques et tous
les mouvements de toutes les particules.

Jean-Baptiste Joseph Fourier

(1768 Auxerre – 1830 Paris)

Orphelin à 10 ans, Joseph est envoyé dans une école militaire où il brille et
devient professeur à 16 ans. Plongé au cœur de la révolution, il échappe de peu
à la guillotine avec la mort de Robespierre. Fin 1794, il est admis à l’école
Normale supérieure de Paris (qui fermera un an après) puis à la toute nouvelle
école Polytechnique. La renommée de Joseph lui permettra de suivre Napoléon
dans son expédition égyptienne. À son retour, il s’installe à Grenoble comme
préfet de l’Isère et effectue durant son temps libre des travaux en mathématiques
dont son mémoire sur la chaleur. En 1817 il est élu à l’académie des sciences dont
il deviendra le secrétaire perpétuel en 1822.

Robert Brown
(1773 Montrose – 1858 Londres)

Après des études de médecine, Robert, botaniste écossais, a été chirurgien dans
l’armée. C’est durant cette période qu’il apprend l’allemand et s’intéresse à la
botanique. Après, il participe à une expédition de 5 ans près des côtes austra-
liennes, grâce à l’appui du président de la Royal Society, Joseph Banks. À la
mort de ce dernier, en 1820, Robert est nommé conservateur du département de
botanique du British Museum. Robert est un des premiers scientifiques à utiliser
fréquemment le microscope. C’est suite à ses observations qu’il découvre puis dé-
crit le mouvement qui porte son nom et qui ne sera expliqué que bien plus tard
avec l’hypothèse atomiste.
© Matthieu Rigaut 7 / 88 Version du 28 déc. 2013
PC⋆ , Fabert (Metz) Biographies succinctes
George Green
(1793 Nottingham – 1841 Nottingham)

Fils de meunier, George Green ne va à l’école que pendant un an (de 8 à 9

ans) puis aide son père au moulin dont il hérite à la mort de celui-ci. C’est alors
qu’il se met à étudier, seul, les mathématiques et publie un ouvrage en 1828. Un
mathématicien passant par là en achète une copie et encourage George à faire des
études. Ce dernier tergiverse mais accepte finalement. Il rentre à l’âge de 40 ans
à l’université de Cambridge où sa carrière fut brillante. Son diplôme en poche en
1837 il commence à publier quelques ouvrages. Mais il tombe malade en 1840 et
rentre dans son village natal pour y mourir quelques temps plus tard.

Mikahail Ostrogradski
(1801 Pachenna, Ukraine – 1862 Poltava, Ukraine)

Mikhail commença ses études de mathématiques à l’université de Kharkov avant

de les continuer à Paris où rencontra et se lia à de celèbres mathématiciens français
comme Cauchy, Binet, Fourier et Poisson. De retour en Russie il enseigne
à Moscou puis à Saint-Pétersbourg. Ses travaux l’amènent à être élu à l’Académie
américaine des arts et des sciences en 1834, à l’Académie des sciences de Turin
en 1841 et à l’Académie des sciences de Rome en 1853.
Ludwig Boltzmann
(1844 Vienne – 1906 Duino)

Ludwick étudie d’abord à Linz en Autriche avant d’entrer à l’université de Vienne

où il obtient un doctorat en 1867 avec une thèse sur la théorie cinétique des gaz. En
1869 il obtient une chaire de physique théorique puis occupe à partir de 1873 une
chaire de mathématique avant de retourner enseigner la physique expérimentale
en 1876. Les échanges très vifs avec ses collègues qui ne croyaient pas à sa théorie
l’affectent particulièrement. C’est très certainement une des causes de ses deux
tentatives de suicide, la seconde lui étant fatale. Ludwig Boltzmann meure avant
d’avoir vu ses idées pleinement acceptées.

Adolf Eugen Fick

(1829 Kassel – 1901 Blankenberge)

Orphelin à 5 ans, Adolf a été élevé par son oncle, un phisiologiste qui faisait des
recherches sur la vision. Cela va conduire tout naturellement Adolf a suivre des
études de médecine. En 1851 il devient docteur et introduit en 1855 une loi sur
la diffusion moléculaire qui sera démontrée par Albert Einstein en 1905. Ayant
attrapé la tuberculose en 1879 il émigre en Afrique du Sud pour le climat mais
sa femme y contracte la fièvre typhoïde et tous reviennent en Europe. Adolf Fick
est principalement connu pour sa plus célèbre invention : les lentilles de contact.

© Matthieu Rigaut 8 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I – Diffusion thermique

I – Diffusion thermique
I·1 – Première approche des transferts thermiques
I·1·i – modes de transfert
✬ la convection

Le transfert thermique par convection correspond à un transport macroscopique de

matière.

✧ Par exemple s’il faut « chaud » dehors et « froid » dedans, en ouvrant la porte il est possible de
sentir un courant d’air (déplacement de matière) apportant le « chaud » de l’extérieur.
✧ Dans ce type de transport, de la matière en chasse une autre mais de température différente.
✧ A priori dès qu’il y a un fluide, il y a de la convection sauf si le dispositif est tel que le fluide est
confiné comme dans les fenêtres double vitrage 1 .

✧ Ceci étant il est très difficile d’empêcher la convection : les chauffages des maisons et même l’eau utili-
sée pour faire cuire les pâtes « utilisent » la convection car cela permet naturellement d’homogénéiser
le milieu et donc la température.

1. Source : [Link]
[Link].

© Matthieu Rigaut 9 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·1 – Première approche des transferts thermiques

Si rien n’est dit, nous négligerons toujours le phénomène de convection.

K Remarque. Pour étudier la convection, il faut faire de la mécanique des fluides avec un fluide non
homogène.

✬ le rayonnement

Les transferts thermiques par rayonnement correspondent à l’absorption d’ondes

électromagnétiques.

✧ Cela traduit le fait que les ondes électromagnétiques transportent de l’énergie.

✧ Expérimentalement, il suffit de s’allonger au Soleil pour avoir une sensation de « chaud » : c’est
simplement que le corps absorbe de l’énergie reçue par rayonnement.
✧ Un autre chauffage très utile qui utilise le rayonnement est le chauffage par micro-ondes.
✧ Nous n’aborderons a priori pas ce mode de transfert énergétique, sauf à prendre le point de vue de
l’onde électromagnétique.

✬ la conduction

La conduction est un mode de transfert énergétique de proche en proche sans

déplacement de matière.

✧ Une expérience classique consiste à plonger une petite cuillère dans un mug contenant du café (ou
du thé) très chaud.

✧ Au bout de quelques minutes, le bout emmergé de la cuillère est devenu très chaud : l’énergie a
remonté le long du manche.

I·1·ii – point de vue de la thermodynamique

✧ Comme nous le verrons, la conduction (ou la diffusion pour être un peu plus général) c’est un simple
problème de thermodynamique auquel nous rajouterons une loi 2 et c’est tout.
2. La loi de Fourier.

© Matthieu Rigaut 10 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·1 – Première approche des transferts thermiques

✧ Typiquement, nous considérerons un système S idoine (nous verrons comment le choisir) et nous
appliquerons le premier principe en version statique

dU = δQ + δW

✧ Nous dirons alors que la transformation est isochore ce qui permettra de simplifier en

dU = δQ

Dans les phénomènes de diffusion thermique nous ferons toujours l’hypothèse de

transformation isochore.

✧ Si tel n’est pas le cas, c’est clairement plus difficile.

✧ N’oublions pas que cette hypothèse revient à négliger la dilatation des objets suite à des variations
de température, ce que nous avons toujours fait.
✧ Pour peu que les variations de température soient modérées, cette approximation est tout à fait
justifiée : qui a déjà vu son bol augmenter de volume en y versant son café ?

I·1·iii – flux thermique

✧ Du point de vue de la thermodynamique, il ne reste plus qu’à écrire δQ.

✬ résultat

Le transfert thermique δQéch échangé entre deux systèmes s’écrit

δQéch = ϕq × dS × dt où :
➜ dS est l’aire de la surface à travers laquelle se fait l’échange ;
➜ dt est la durée de l’échange ;
➜ ϕq ≷ 0 est le flux surfacique thermique en W.m−2 , c’est un flux surfacique de puissance
algébrique.

✧ Parfois δQ est noté δ 2 Q pour insister sur le fait qu’il provient de deux infiniments petits de nature
différentes (un d’espace et un de temps). Nous ne le ferons pas, sauf si cela permet d’éviter des
collusions de notation.
✧ Nous n’avons certes pas encore tout à fait répondu à la question, mais, pour autant, nous n’avons
pas rien dit.
✧ En effet, même si cela paraît logique et naturel, cette relation impose le fait que le transfert thermique
est proportionnel à la surface d’échange et à la durée d’échange !

✬ justification
✧ Imaginons deux systèmes séparés par une surface fictive.
© Matthieu Rigaut 11 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·1 – Première approche des transferts thermiques
surface d’échange
réelle ou fictive

milieu ① milieu ②

δQ

✧ De l’énergie passe de l’un à l’autre.

✧ Comme l’énergie est une grandeur extensive, ce qui passe à un endroit s’additionne avec ce qui passe
ailleurs ce qui signifie que

δQ ∝ dS
✧ De même l’énergie s’accumule dans le temps puisque rien ne peut la détruire donc ce qui arrive à un
moment s’ajoute à ce qui arrive après d’où

δQ ∝ dt
✧ Par définition, nous poserons ϕq le facteur de proportionnalité, facteur que nous avons interprété
comme le flux surfacique de puissance.

✬ un dernier mot
✧ Cela paraît peut-être évident, mais comme nous aurons plus que l’occasion d’en parler et reparler,
autant insister dès maintenant.

Les échanges thermiques se font à travers la surface délimitant un système.

I·1·iv – production énergétique

✬ phénoménologie de base

La production énergétique se fait au sein du système.

✧ Quelquefois la production énergétique est appelée « production de chaleur » mais c’est ambigu car
cela revient à considérer que la « chaleur » n’est pas une énergie comme une autre et renforce l’idée
commune (et fausse) que la température est proportionnelle à la chaleur.

✬ 3 grands types de production

Un matériau peut produire de l’énergie en son sein suite à trois grands phénomènes :
➜ les réactions chimiques (en incluant les changements de phase) ;
➜ les réactions nucléaires ;
➜ l’effet Joule.

✧ En réalité, nous parlerons, du point de vue de la diffusion thermique, de « production » d’énergie

alors que nous savons bien qu’il est impossible de produire / créer de l’énergie, tout juste pouvons
nous la transformer :
© Matthieu Rigaut 12 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·2 – Bilan thermique

➜ pour les réactions chimiques (ou les changements de phase), l’énergie est initialement sous forme
d’énergie de liaison atomique ;
➜ pour les réactions nucléaires l’énergie est initialement sous forme d’énergie de liaison nucléaire ;
➜ pour l’effet Joule l’énergie est initialement dans le champ électrique responsable du courant.
✧ Dans les trois cas cette énergie est transformée en énergie interne.
✧ Dans ces conditions, i.e. toujours du point de vue de la diffusion, la distinction de l’origine de l’énergie,
à savoir l’écriture sous la forme d’un dUint pour l’énergie de liaison ou d’un δW pour l’effet Joule
ne sert à rien, nous nous conterons de regarder l’énergie qui arrive et non d’où elle vient ce qui nous
permettra de parler dans tous les cas de l’énergie produite.

✬ traduction formelle

L’énergie produite δQprod dans un système de volume δτ s’écrit

δQprod = P × δτ × dt où :
➜ dt est la durée de production ;
➜ P ≷ 0 est la puissance volumique produite en W.m−3 .

✧ Le cas le plus simple de production négative d’énergie est celui d’une réaction chimique endother-
mique.

I·2 – Bilan thermique

I·2·i – courant thermique
✬ l’énergie va se déplacer
✧ Revenons sur la petite cuillère dans le café dont le manche chauffe.

✧ Si la température augmente, c’est que l’énergie cinétique moyenne des atomes qui composent la petite
n
cuillère augmente puisque nous avons 3 hec i = kB T .
2
✧ Cela prouve que de l’énergie s’est déplacée à l’intérieur du manche.

3. n est le nombre de termes quadratiques énergétiques.

© Matthieu Rigaut 13 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·2 – Bilan thermique

Le vecteur ~th représente le déplacement d’énergie dans un problème de diffusion.

milieu diffusant

~th

✧ Considérons quelques cas simples et cherchons ce que peut valoir le transfert thermique δQ1→2 .

milieu ① milieu ② milieu ① milieu ②

~th
~
dS ~th ~
dS

cas 1 cas 2
milieu ① milieu ②

~th ~
dS

cas 3

✧ Dans le premier cas nous avons δQ1→2 = 0 car l’énergie ne passe pas à travers la surface dS, elle la
longe mais ne la traverse pas.
✧ Dans le deuxième cas nous avons δQ1→2 > 0 car l’énergie passe effectivement de 1 vers 2.
✧ Enfin dans le troisième cas nous avons δQ1→2 < 0 car l’énergie va de 2 vers 1.
✧ Imaginons maintenant un quatrième cas.

milieu ① milieu ②

~th ~
dS

cas 4

✧ Cette fois nous avons δQ1→2 > 0 mais quand même inférieur à ce qu’il était quand ~th était normal
à la surface.

✬ généralisation
✧ Finalement, en réunissant toutes ces remarques, nous pouvons généraliser en posant la définition
précise de ~th .

Le vecteur densité surfacique de courant thermique en volume ~th est défini par
~ dt
δQ = ~th · dS où :
~ est le vecteur surface de la surface à travers laquelle se fait l’échange δQ ;
➜ dS
➜ δQ est le transfert thermique dans le sens de dS ~;
➜ dt est la durée du transfert.

© Matthieu Rigaut 14 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·2 – Bilan thermique

✧ De cette définition, nous pouvons tirer immédiatement

En notant ϕq le flux surfacique de 1 vers 2 et ~n12 le vecteur normal dirigé de 1 vers 2,

nous avons
ϕq = ~th · ~n12

✧ Et ainsi

Le vecteur densité surfacique de courant thermique en volume ~th est en W.m−2 .

✧ L’unité de ~th explique son nom :

➜ c’est bien une densité surfacique puisque ~th est par unité de surface ;
➜ ce vecteur correspond bien à de l’énergie qui se déplace dans un volume.
✧ La plupart du temps ~th est appelé « densité de courant thermique volumique » ce qui présente
l’avantage d’être plus léger mais aussi l’inconvénient de laisser sous entendre que l’unité est en m−3
ce qui est faux.
✧ C’est pourquoi, même si la tournure est plus lourde, nous parlerons de « densité surfacique de courant
thermique en volume » ou, pour alléger un peu, de « courant thermique en volume ».
! Les énoncés n’ont pas tous le même vocabulaire donc il faudra systématiquement aller doucement et
réinterpréter physiquement les phénomènes avant d’essayer de les traduire analytiquement.

✬ avec la convention usuelle

Pour un système fermé, comme les vecteurs surface élémentaires sont orientés vers
l’extérieur, nous avons
~ dt
δQreçu = −~th · dS

~
dS

✬ simplification usuelle

Lorsque le courant thermique en volume est uniforme, le transfert thermique algébrique

~ non infinitésimale s’écrit
à travers une surface S
~ dt
δQ = ~th · S

✧ La « démonstration » est immédiate. Comme le transfert thermique est extensif, nous avons
© Matthieu Rigaut 15 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·2 – Bilan thermique

Z Z
δQ = 2
δ QP δQ = ~P dt
~(P ) · dS
P ∈S P ∈S

✧ Comme le courant thermique est uniforme, nous avons ~(P ) = ~ indépendant de P , ce qui nous permet
de le factoriser et ainsi
Z !
δQ = ~ · ~P
dS dt ~ dt
δQ = ~ · S
P ∈S

✧ Ce qui est bien le résultat.

I·2·ii – bilan 1D
✧ Nous n’allons pas tout de suite rechercher l’équation vérifiée par la température mais plutôt celle
vérifiée par l’énergie.
✧ Nous verrons le lien entre énergie et température après.

✬ exemple du barreau calorifugé

✧ Considérons l’exemple canonique suivant.

source source
milieu diffusant
T1 T2

✧ Un barreau en métal est au contact de deux thermostats, un à chaque extrémité.

✧ Nous pouvons supposer que celui-ci est isolé thermiquement sur les côtés de sorte que l’énergie ne
peut se propager que suivant l’axe.
✧ Analyse physique :
➜ le barreau va être le siège de phénomènes de diffusion car il est en contact avec des thermostats
de températures différentes ;
➜ la température ne va donc pas être la même partout et va évoluer dans le temps.
✧ Analyse technique :
➜ comme ici la température varie localement, nous allons faire une approche mésoscopique du
problème ;
➜ pour rester dans les conventions usuelles, nous noterons ~ux l’axe avec 0 à gauche et L à droite.

système étudié
source source
milieu diffusant
T1 T2

~ux x x + δx
O L

© Matthieu Rigaut 16 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·2 – Bilan thermique

✬ approche mésoscopique
✧ Commençons par faire un zoom sur la partie intéressante, ce qui sera notre système S : la tranche
comprise entre x et x + δx.

~(x,t) ~(x + δx,t)

δQ(x) δQ(x + δx)

x x + δx

✧ Faisons un bilan d’énergie sur S entre t et t + dt.

✧ Ce bilan peut se résumer sous la forme

VARIATION dans le temps = ÉCHANGE à travers la surface + CRÉATION en volume

✧ Trois termes donc. . .

✬ variation dans le temps

✧ Ici comme nous avons une transformation isochore, nous allons utiliser l’énergie interne ce qui donne

U (t + dt) − U (t) = CV (T (x,t + dt) − T (x,t))

✧ Et avec une phase condensée (ici un solide) CV = C = ρ c δτ où :

➜ ρ est la masse volumique ;
➜ c est la capacité thermique massique ;
➜ δτ est le volume du système.
✧ Finalement

VARIATION dans le temps = ρ c δτ (T (x,t + dt) − T (x,t))

✬ échange à travers la surface

✧ Regardons le schéma et constatons que :
➜ en tant que cylindre, la surface se décompose naturellement en 3 : la section en x, la section en
x + δx et la surface latérale ;
➜ l’énergie ne passe pas à travers la surface latérale.

section en x section en
x + δx
surface latérale

© Matthieu Rigaut 17 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·2 – Bilan thermique

✧ L’énergie échangée s’écrit donc

δQéch = δQ(x,t) + δQ(x + δx,t)

✧ Or, d’après la définition du courant thermique en volume et parce qu’il est uniforme sur la section,

δQ(x) = +jth (x,t) S dt et δQ(x + δx,t) = −jth (x + δx,t) S dt

✧ Remarquons que le signe n’est pas le même, ce qui peut se voir de deux manières différentes :
➜ soit en constatant que lorsque ~th est dans le sens de ~ux , alors ~th (x,t) « rentre » alors que
~th (x + δx,t) « sort » ;
➜ soit en introduisant les vecteurs surfaces avec S~ (x) = −S ~ux et S
~ (x + δx) = +S ~ux et en utilisant
~
la convention usuelle δQreçue = −~th · S.
✧ Dans tous les cas, il ne faut pas se tromper dans les signes.
✧ Ce terme peut se réécrire

δQéch = +jth (x,t) S dt − jth (x + δx,t) S dt

= S dt jth (x,t) − jth (x + δx,t)

✧ Ce qui donne, avec un petit développement limité à l’ordre 1

∂jth
ÉCHANGE à travers la surface = −S dt (x,t) δx
∂x

✬ création en volume
✧ Par définition même de l’expression de l’énergie produite en volume, nous avons tout de suite

CRÉATION en volume = P (x,t) dt δτ

✬ rassemblement
✧ En reprenant le bilan initial nous avons

∂jth
ρ c δτ (T (x,t + dt) − T (x,t)) = −S dt (x,t) δx + P (x,t) dt δτ
∂x
✧ En remarquant que δτ = S δx et en divisant par δτ dt nous obtenons

(T (x,t + dt) − T (x,t)) ∂jth

ρc =− (x,t) + P (x,t)
dt ∂x
✧ En faisant tendre dt vers 0, nous avons

(T (x,t + dt) − T (x,t)) ∂T

−→ (x,t)
dt ∂t
✧ Nous arrivons alors à

∂T ∂jth ∂T ∂jth
ρc (x,t) = − (x,t) + P (x,t) ρc (x,t) + (x,t) = P (x,t)
∂t ∂x ∂t ∂x
© Matthieu Rigaut 18 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·2 – Bilan thermique

✬ résultat connu

Dans le cas d’une diffusion thermique unidimensionnelle, la conservation de l’énergie se

traduit par
∂T ∂jth
ρc (x,t) + (x,t) = P (x,t)
∂t ∂x

✧ Nous retrouvons là les lois de conservation de l’énergie que nous avions établies dans le cas de la
propagation.
✧ La seule vraie différence est qu’ici il peut y avoir un terme de production.

I·2·iii – bilan 3D
✬ situation
✧ Nous allons cette fois faire un bilan énergétique entre t et t + dt mais pour un système compris dans
un volume V quelconque.

~
dS

dτ

✧ Nous ferons les hypothèses usuelles suivantes :

➜ chaque point du système est immobile donc V ne dépend pas du temps ;
➜ la transformation est isochore (ce qui est un peu lié à l’hypothèse précédente) ;
➜ le corps est une phase condensée.
✧ Le bilan d’énergie sur S entre t et t + dt s’écrit donc, comme tout bilan :

VARIATION dans le temps = ÉCHANGE à travers la surface + CRÉATION en volume

✧ Explicitons les différents termes

✬ variation dans le temps

✧ Tout d’abord nous ne pouvons pas écrire
✭✭✭ ✭
✭✭ ✭ ✭✭✭
✭
U (t + dt) − U (t) = dU VARIATION dans ✭✭
✭le ✭✭✭= C (T (t + dt) − T (t))
temps
✭ ✭
✭✭✭✭
✭✭✭
✧ En effet une telle écriture laisse entendre que la température est définie pour l’ensemble du système
ou, autrement dit, que le système a une température homogène.
✧ C’est bien évidemment faux a priori donc nous ne pouvons pas faire cela.
✧ Commençons alors par exprimer l’énergie interne contenue à t dans le système.
© Matthieu Rigaut 19 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·2 – Bilan thermique

✧ L’extensivité de l’énergie interne pour permet de voir le système comme une réunion de plein de
petits morceaux de masse dmP et donc en notant u(P,t) l’énergie interne massique au point P
ZZZ
U (t) = u(P,t) dmP
V(t)
✧ De même l’énergie interne en t + dt s’écrit
ZZZ
U (t + dt) = u(P,t + dt) dmP
V(t+dt)
✧ La variation d’énergie s’écrit donc
ZZZ ZZZ
U (t + dt) − U (t) = u(P,t + dt) dmP − u(P,t) dmP
V(t+dt) V(t)
✧ Mais comme V (t) = V (t + dt) puisque le système est immobile, nous pouvons réunir les deux sommes
ZZZ ZZZ ZZZ


u(P,t + dt) dmP − u(P,t) dmP = u(P,t + dt) − u(P,t) dmP
V V V
✧ De là nous voyons apparaître une différentielle en faisant un développement limité à l’ordre 1
ZZZ
∂u
U (t + dt) − U (t) = (P,t) dt dmP
V ∂t
✧ Et parce que nous avons une phase condensée nous pouvons écrire, avec dmP = ρ δτP
ZZZ
∂T
VARIATION dans le temps = ρc (P,t) dt dτP
V ∂t
! Même si cela sera souvent le cas, ρ et c ne sont pas forcément uniformes ce qui nous interdit de les
factoriser, i.e. de les sortir du signe somme.

✬ production en volume
✧ Par extensité de l’énergie produite par chaque portion de volume dτP nous pouvons écrire immédia-
tement
ZZZ
CRÉATION en volume = P (P,t) dt dτP
V

✬ échange à travers la surface

✧ Par extensivité de l’énergie reçue, nous pouvons écrire que
ZZ
δQreçu = δ 2 Qreçue (P,t)
S
✧ Remarquons qu’ici la surface est fermée puisque le système est de taille finie.
✧ Nous pouvons réécrire ce terme d’échange avec le vecteur densité surfacique de courant thermique
en volume ~th
ZZ
δQreçu = ~P dt
− ~th (P,t) · dS
S
✧ Ce qui donne
ZZ
ÉCHANGE à travers la surface = − ~P dt
~th (P,t) · dS
S

© Matthieu Rigaut 20 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·2 – Bilan thermique

✬ rassemblement
✧ Le bilan complet donne donc
ZZZ ZZ ZZZ
∂T ~
ρc (P,t) dt dτP = − ~th (P,t) · dSP dt + P (P,t) dt δτP
V ∂t S V

✧ Il y a deux termes dont la sommation est sur le volume et un troisième dont la sommation est sur la
surface. . .
✧ Pour pouvoir aller plus loin, nous allons utiliser le truc qui transforme une sommation sur la surface
en une sommation en volume

Théorème de Green – Ostrogradski

Pour un volume V quelconque délimité par la surface S nous avons
ZZ ZZZ
~
~th (P,t) · dSP = div ~th (P,t) δτP
S V
div ~th est appelé la divergence de ~th

✧ Deux remarques :
➜ ce théorème marche avec n’importe quel champ de vecteur et pas seulement avec ~th qui n’est
là que parce que c’est la première fois que nous rencontrons la notion de divergence ;
➜ le « quelconque » pour le volume est à prendre au sens physique, i.e. « délimité n’importe
quelle surface raisonnable ». Car avec des surfaces pathologiques du type fractale ou bouteille
de Klein ce théorème ne marche plus. Mais fort heureusement cela ne peut pas nous arriver
puisque c’est nous qui choisissons la surface et nous en choisirons systématiquement des surfaces
qui marchent.
✧ Revenons à nos moutons et voyons ce que cela donne
ZZZ ZZZ ZZZ
∂T
ρc (P,t) dt dτP = − div ~th (P,t) δτP dt + P (P,t) dt δτP
V ∂t V V

✧ Comme, encore une fois, tous les volumes de sommation sont les mêmes, nous pouvons tout regrouper
sous le même signe somme ce qui donne

ZZZ !
∂T
ρc (P,t) + div ~
th (P,t) − P (P,t) dt dτP = 0
V ∂t

✧ Ce résultat est vrai quel que soit le volume V ce qui implique que quelle que soit la manière dont
nous sommons l’intégrande, le résultat est nul.
✧ Cela ne peut vouloir dire qu’une chose : que l’intégrande est nulle.

∂T ∂T
ρc + div ~th − P = 0 ρc + div ~th = P
∂t ∂t

La conservation de l’énergie lors d’une diffusion en trois dimension s’écrit

∂T
ρc + div ~th = P
∂t

© Matthieu Rigaut 21 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·2 – Bilan thermique

I·2·iv – pause divergence

✧ Maintenant que nous avons « inventé » un truc appelé divergence, il va falloir regarder d’un peu plus
près ce que c’est.
✧ Il faut tout d’abord savoir que le concept de divergence va être très utilisé tout au long de l’année.
✧ Il est donc indispensable de travailler cette notion tant au nivau physique (comprendre ce qu’il repré-
sente, quand s’en servir) qu’au niveau technique (pour le manipuler correctement dans les équations).

✬ relation fondamentale
✧ C’est la relation qui définit la divergence.

~ à travers la
Pour un volume élémentaire δτ , le flux sortant d’un champ de vecteur A
surface délimitant δτ s’écrit
dΦ , div A ~ × δτ

✧ Autrement dit, nous pouvons transformer une sommation sur une surface en une sommation en
volume

~
dS

δτ

✧ Et cela donne le théorème de Green – Ostrogradski

ZZ ZZZ
Φ= ~ ~
A · dS = ~ × δτ
div A
S V
✧ Cela revient à dire que, pour compter ce qui sort Φ : ZZ
➜ soit il faut regarder ce qui se passe au niveau de la surface ~ · dS
A ~;
ZZZ S

➜ soit il faut regarder ce qui se passe à l’intérieur ~ .

div Aδτ
V

✬ comment ça marche ?
✧ Regardons sur un exemple en 2D comment le théorème fonctionne.
✧ Commençons par prendre un système que nous découperons en 9 cases.

© Matthieu Rigaut 22 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·2 – Bilan thermique

✧ Sur chaque bordure mettons des flèches au hasard

✧ Pour chaque case, comptons les flèches qui sortent

−1 4 0

3 0 −4

1 −2 3

✧ Maintenant comptons le nombre de flèches qui sortent du gros système

8−4= 4
✧ Nous constatons alors que

−1 + 4 + 0 + 3 + 0 − 4 + 1 − 2 + 3 = 4

✧ Il s’agit là ni plus ni moins que du nombre de flèches qui sortent du gros système !
✧ L’explication est simple. Quand nous regardons une flèche séparant deux cases, nous pouvons voir
qu’elle compte obligatoirement positivement pour l’une et négativement pour l’autre

+1 −1

✧ Dans ces conditions, lorsque nous sommons le tout, toutes les flèches à l’intérieur ont une contribution
nulles et il ne reste que les flèches donnant sur l’extérieur.
© Matthieu Rigaut 23 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·2 – Bilan thermique

✬ interprétation physique
✧ Au niveau local, la divergence représente le caractère local d’un champ de vecteur.
✧ Quand le champ « sort », la divergence est positive.

✧ Quand le champ « rentre », la divergence est négative.

✧ Enfin, quand le champ « rentre » autant qu’il sort, la divergence est nulle.

✧ Bref, que la divergence s’appellent divergence et représente la manière dont des lignes de champs
sortent (donc « divergent ») à partir d’un point n’est peut-être pas qu’une heureuse coïncidence.
✧ En électromagnétostatique, nous avons déjà rencontré des champs de vecteurs que nous représentions
sous forme de lignes.
✧ Voici, ci-dessous, deux exemples de cartes de champ, le premier est un champ électrostatique, le
deuxième un champ magnétique.

✧ C’est ainsi que :

➜ pour le champ électrostique, nous savons que les lignes de champs convergent vers les charges
~ non seulement
négatives et divergent à partir des charges positives. Nous aurons donc div E
est relié aux charges, mais en plus est du même signe que les charges ;
➜ les lignes de champs magnétique ne peuvent pas diverger donc div B ~ = 0. C’est là ni plus ni
moins qu’une des 4 équations fondamentales de l’électromagnétisme.

© Matthieu Rigaut 24 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·2 – Bilan thermique

✬ expression en coordonnées cartésiennes

En coordonnées cartésiennes la divergence d’un champ de vecteur ξ~ s’écrit

∂ξx ∂ξy ∂ξz
div ξ~ = + + où :
∂x ∂y ∂z
ξx , ξy et ξz sont les composantes de ξ~ sur respectivement ~ux , ~uy et ~uz .

✧ Pour démontrer ce résultat, partons d’un volume élémentaire δτ et utilisons la relation fondamentale
de la divergence.
z + dz
~uz

x z ~uy
~ux
x + dx
y y + dy

✧ Le flux traversant la surface délimitant δτ peut donc se décomposer en 6 puisqu’il y a 6 surfaces

élémentaires

dΦ = dΦ(x) + dΦ(x + dx) + dΦ(y) + dΦ(y + dy) + dΦ(z) + dΦ(z + dz)

✧ En remplaçant par leurs expressions cela donne

dΦ = ξ~(x,y,z) · (−dy dz ~ux ) + ξ~(x + dx,y,z) · (dy dz ~ux ) · · ·

· · · + ξ~(x,y,z) · (−dx dz ~uy ) + ξ~(x,y + dy,z) · (dx dz ~uy ) · · ·
· · · + ξ~(x,y,z) · (−dx dy ~uz ) + ξ~(x,y,z + dz) · (dy dz ~uz )

✧ Et ainsi en regroupant adroitement

dΦ = ξx (x + dx,y,z) − ξx (x,y,z) dy dz · · ·
 
· · · + ξy (x,y + dy,z) − ξy (x,y,z) dx dz · · ·


· · · + ξz (x,y,z + dz) − ξz (x,y,z) dx dy

✧ En faisant des développement limités, nous voyons alors apparaître 3 dérivées partielles
∂ξx ∂ξy ∂ξz
Φ = dx (x,y,z) × dy dz + dy (x,y,z) × dx dz + dz (x,y,z) × dx dy
∂x ∂y ∂z
✧ Puis, en factorisant
!
∂ξx ∂ξy ∂ξz
Φ= (x,y,z) + (x,y,z) + (x,y,z) × dx dy dz
∂x ∂y ∂z
✧ Comme dx dy dz = δτ , par identification avec la relation fondamentale de la divergence, nous obte-
nons bien le résultat attendu
© Matthieu Rigaut 25 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·2 – Bilan thermique

∂ξx ∂ξy ∂ξz

Φ = div ξ~ × dx dy dz div ξ~ = (x,y,z) + (x,y,z) + (x,y,z)
∂x ∂y ∂z

✬ c’est un opérateur vectoriel différentiel

La divergence est un opérateur vectoriel différentiel et

~ = [ξ]
[div ξ]
longueur

Avec λ1 et λ2 des scalaires quelconques, nous avons

 
div λ1 ξ1 + λ2 ξ2 = λ1 div ξ~1 + λ2 div ξ~2
~ ~

✧ La démonstration basée sur la linéarité de l’opérateur différentielle est laissée au lecteur.

✬ c’est une vision 3D d’un théorème connu

✧ Regardons à nouveau le théorème de Green – Ostrogradski
ZZ ZZZ
~ · dS
A ~= ~ × δτ
div A
S V

✧ Nous voyons que nous pouvons remplacer une somme double par une somme triple à condition de
dériver l’intégrande.
✧ C’est exactement la même philosophie que d’écrire
Z b
f (b) − f (a) = f ′ (x) dx
a

I·2·v – petit flash-back

✬ passer de la 3D à la 1D
✧ Reprenons l’équation de conservation de l’énergie dans le cas de la diffusion en version 3D

∂T
ρc + div ~th = P
∂t
✧ Pour une diffusion unidimensionnelle, nous avons

∂jth
~th = jth (x,t) ~ux div ~th =
∂x
✧ Et ainsi l’équation de conservation de l’énergie se réécrit

∂T ∂jth
ρc + =P
∂t ∂x
✧ Ce qui est bien (et heureusement) ce que nous avions trouvé en raisonnant directement en 1D.
© Matthieu Rigaut 26 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·2 – Bilan thermique

✬ passer de la 1D à la 3D sur un cas déjà rencontré

✧ Reprenons l’équation de conservation de l’énergie que nous avions trouvée dans le cas de la propa-
gation.

∂e ∂P
+ =0
∂t ∂x
✧ Remarquons que cette équation est exactement la même que dans le cas de la diffusion car :
➜ P dans le cas de la propagation représente bien de l’énergie qui bouge, i.e. n’est autre que le
~th ;
➜ le terme de production dans le cas de la propagation est nul.
✧ Nous pouvons donc imager ce que serait l’équation de conservation de l’énergie dans le cas de la
propagation mais en 3D

∂e ~ =0
+ div P
∂t
✧ Les interprétations sont alors :
➜ e l’énergie volumique (et non plus linéique) contenue dans l’onde ;
➜ P ~ un vecteur flux de puissance surfacique.
 En électromagnétisme P ~ n’est autre que le vecteur de Poynting Π.
~

I·2·vi – régime stationnaire sans production

✧ Regardons un cas particulier assez fréquent du cas où la diffusion est envisagée dans un cas station-
naire sans terme de production.

✬ en 1D
✧ Regardons le cas particulier du régime stationnaire.

Un régime est dit stationnaire lorsqu’aucune grandeur ne dépend du temps.

✧ Dans ces conditions la conservation de l’énergie donne

∂jth
0+ =0 jth (x) = Cte
∂x

Dans le cas d’une diffusion unidimensionnelle sans terme de production et en régime

stationnaire, la densité surfacique de courant thermique en volume est uniforme.

✧ Même si elle est rapide à démontrer, c’est une loi que nous pourrons utiliser directement.

✬ continuité du flux en 3D
✧ Sous ces mêmes conditions, la conservation de l’énergie en version 3D donne

La diffusion tridimendionnalle sans terme de production et en régime stationnaire se

traduit par
div ~th = 0

© Matthieu Rigaut 27 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·3 – Équation de diffusion thermique

✧ Autrement dit

Dans le cas diffusion tridimendionnalle sans terme de production et en régime

stationnaire, la densité surfacique de courant thermique en volume est à flux conservatif.

✧ Rappelons que « à flux conservatif » signifie que quelle que soit la surface fermée choisie, le flux est
nul.
✧ Nous pouvons le montrer rapidement grâce au théorème de Green – Ostrogradski
ZZ ZZZ
Φ= ~
~ · dS Φ= div ~ δτ
S V

✧ Et comme div ~ δτ = 0, le résultat tombe.

✧ Là aussi il s’agit d’un résultat rapide à démontrer mais que nous utiliserons parfois immédiatement.

I·3 – Équation de diffusion thermique

✧ Malgré tout ce que nous venons de faire, nous n’avons pas encore tout à fait l’équation de diffusion
thermique.
✧ Physiquement c’est évident parce que nulle part nous avons expliqué de quelle manière l’énergie
bougeait (vite ou pas vite ? dans quelle direction ?)
✧ Techniquement cela se voit car le membre de gauche de l’équation de conservation de l’énergie ne
fait pas intervenir que la température, elle n’est donc pas soluble en l’état.

I·3·i – loi de Fourier

✧ C’est la loi phénoménologique qui relie
➜ ~th qui représente le mouvement de l’énergie ;
➜ T qui est un marqueur de présence d’énergie.

✬ énoncé

loi de Fourier
Pour un matériau usuel
−−→
~th = −λ grad T où :
➜ ~th est le vecteur densité surfacique de courant thermique en volume ;
➜ λ > 0 est la conductivité thermique et dépend du matériau.

L’unité de la conductivité thermique est

[λ] = W.m−1 .K−1

✧ En effet

−−→ W.m−2
[~th ] = W.m−2 et [grad T ] = K.m−1 [λ] = = W.m−1 .K−1
K.m−1
© Matthieu Rigaut 28 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·3 – Équation de diffusion thermique

✬ interprétation

La loi de Fourier traduit le fait que l’énergie se déplace des zones chaudes vers les
zones froides dans le cadre de la conduction thermique.

−−→
✧ En effet, grad T est dirigé vers les hautes températures et la présense du signe − associé au fait que
λ ne peut être que positif impose ce résultat.
✧ En fait la loi de Fourier traduit ce que nous savons du second principe.
✧ Toutefois la loi de Fourier va un peu plus loin en précisant comment l’énergie se déplace.

✬ limites
✧ La loi de Fourier est une loi linéaire faisant apparaître une dérivée première de l’espace (le gradient).
✧ Autrement dit, utiliser la loi de Fourier revient à limiter au premier ordre les effets de la diffusion.
−−→
✧ Nous pouvons alors en déduire la première restriction : il ne faut pas que kgrad T k soit trop grand
sans quoi il faudrait ajouter un terme correctif (non linéaire) du second ordre.
−−→
✧ De plus si kgrad T (t)k varie trop rapidement, il peut y avoir un temps de réponse au niveau moléculaire
−−→
entre ~th et grad T .
✧ Enfin, pour pouvoir utiliser la loi de Fourier il faut que le matériau soit isotrope sinon le λ sur ~ux
ne serait pas le même que le λ sur ~uy .
✧ C’est ainsi que dans le graphique, matériau composé de feuillets de carbone 4 , la conductivité ther-
mique suivant les feuillets est plusieurs centaines de fois plus grande que la conductivité thermique
entre les feuillets.

Sauf précision explicite contraire, nous supposerons a priori que la loi de Fourier est
valide.

✧ Cela ne nous empêchera pas de discuter des résultats obtenus vis à vis de la validité de la loi de
Fourier !

I·3·ii – équation en 1D
✬ équation de diffusion

L’équation de diffusion thermique unidimensionnelle dans un matériau uniforme s’écrit

∂T ∂2T
ρc −λ 2 =P
∂t ∂x
➜ ρ est la masse volumique ;
➜ c est la capacité thermique massique ;
➜ λ est la conductivité thermique ;
➜ P est la production volumique d’énergie.

✧ Pour le montrer, commençons par reprendre l’équation de conservation de l’énergie en 1D.

∂T ∂jth
ρc + =P
∂t ∂x
✧ Introduisons, dans cette équation, la loi de Fourier qui s’écrit, en 1D
4. Voir cours de chimie.

© Matthieu Rigaut 29 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·3 – Équation de diffusion thermique

!
∂T ∂T ∂ ∂T
jth = −λ ρc + −λ =P
∂x ∂t ∂x ∂x
✧ Nous pouvons alors « sortir » la conductivité de la dérivée parce que le matériau est uniforme. Cela
nous conduit à l’équation attendue

∂T ∂2T
ρc −λ 2 =P
∂t ∂x

✬ cas particulier fréquent : diffusion sans terme de création

L’équation de diffusion thermique unidimensionnelle dans un matériau uniforme et sans

terme de source s’écrit
∂T ∂2T λ
=a 2 avec a=
∂t ∂x ρc
➜ ρ est la masse volumique ;
➜ c est la capacité thermique massique ;
➜ λ est la conductivité thermique ;
➜ a est le coefficient de diffusivité thermique en m2 .s−1 .

✧ Cette équation découle de l’équation précédente avec P = 0.

✧ Remarquons que dans cette équation, le coefficient de diffusivité thermique est le seul coefficient
phénoménologique qui apparaît.
✧ Cela a pour première conséquence que techniquement les solutions ne dépendent que de a.
✧ Physiquement a représente le couplage spatio-temporel de la diffusion.
✧ D’une certaine manière, ce n’est pas surprenant que l’équation de diffusion ne dépende que d’un
terme de couplage car c’est ce que nous avions trouvé avec l’équation de propagation qui, elle, ne
dépendait aussi que d’un terme de couplage spatio-temporel, à savoir c.

✬ ordres de grandeur
Y pour la loi de Fourier
✧ Regardons quelques valeurs de conductivité thermique

matériau cuivre eau stagnante verre pyrex polystyrène

1 −1
λ (W.m .K ) 399 0,597 1,4 0,027

✧ Nous pouvons constater que les conductivités thermiques peuvent varier sur 5 ordres de grandeurs
ce qui est considérable.
✧ De plus, conformément à l’intuition, le polystyrène conduit (beaucoup) moins bien la chaleur que le
cuivre ou l’eau.

Y pour l’équation de diffusion

✧ Regardons maintenant quelques valeurs de diffusivité thermique

matériau cuivre eau stagnante verre pyrex polystyrène

a (m .s ) 117.10−6
2 −1
0,143.10−6 0,75.10−6 0,41.10−6

© Matthieu Rigaut 30 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·3 – Équation de diffusion thermique

✧ De manière surprenante nous voyons que la diffusivité du polystyrène est supérieure à celle de l’eau !
✧ Il ne faut surtout pas en conclure que le polystyrène conduit mieux la « chaleur » que l’eau, non,
mais plutôt que le polystyrène est plus rapide à réagir que l’eau.

I·3·iii – interprétation de l’équation de diffusion

✧ Nous allons tout de suite interpréter l’équation de diffusion bien qu’elle ne décrive que le cas unidi-
mensionnel car l’interprétation en 3D est la même.

✬ approche technique
✧ Tout d’abord c’est une équation dont nous n’avons pas, « comme ça » toutes les solutions, contrai-
rement à l’équation de propagation où nous pouvions les écrire a priori.
✧ Toutefois il s’agit d’une équation aux dérivées partielles linéaire ce qui peut parfois faciliter la réso-
lution.
✧ Par exemple dans le cas d’un régime forcé, nous pourrons étudier chaque composante sinusoïdale
séparemment et sommer le tout avec le principe de superposition, une composante s’écrivant

T (x,t) = T0 (x) e j ω t

✬ irréversibilité

L’équation de diffusion traduit le caractère irréversible du phénomène.

✧ En effet, si nous regardons ce qui se passe lorsque nous remontons le temps, i.e. lorsque t −→ t′ = −t,
les dérivées partielles par rapport au temps changent de signe

∂ ∂ ∂
−→ ′ = −
∂t ∂t ∂t
∂
✧ C’est normal : quand quelque chose augmente avec le temps > 0 et que nous remontons le temps,
∂t
∂
cette même chose se met à décroître ′ < 0.
∂t
✧ L’équation de diffusion vue en remontant le temps s’écrit

∂T ∂2T
−ρ c − λ =P
∂t′ ∂x2
✧ Cette équation n’est pas une équation de diffusion.
✧ Cela signifie qu’en remontant le temps nous voyons des phénomènes non diffusif : c’est bien irréver-
sible.
✧ Ce genre d’irréversibilité a déjà été beaucoup vu en thermodynamique de première année et est du
au fait que le « chaud » va vers le froid et non le contraire.

✬ ordre de grandeur
✧ Réécrivons l’équation de diffusion thermique

∂T ∂2T
=a 2
∂t ∂x
© Matthieu Rigaut 31 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·3 – Équation de diffusion thermique

✧ En notant T0 l’ordre de grandeur de la température et τ le temps caractéristique, nous pouvons dire

non seulement d’un point de vue dimensionnel mais aussi parce qu’une dérivée n’est qu’une division
déguisée que

∂T T0
≡
∂t τ
✧ De même pour l’autre terme, comme nous dérivons deux fois
!
∂2T ∂ ∂T ∂2T T0
a 2 =a× a 2 ≡a 2
∂x ∂x ∂x ∂x L

✧ Finalement

T0 T0 √
=a 2 L= aτ
τ L

La distance caractéristique L de diffusion est :

➜ d’autant plus grande
√ que la diffusivité a est grande ;
➜ proportionnelle à τ où τ est le temps caractéristique : c’est un phénomène lent.

✧ Pour passer de L à 10 L, il faut passer de τ à 100 τ , ce qui n’est pas rien !

✧ Pour la propagation, comme nous avons L = c τ , pour passer aller 10 fois plus loin il suffit d’attendre
10 fois plus longtemps.

✬ régime permanent stationnaire

✧ Rappelons ici la différence avec le régime permanent.

Un régime est dit permanent lorsque le régime transitoire est terminé.

✧ Il peut très bien y avoir un régime permanent non stationnaire comme, par exemple, tout ce qui est
régime forcé (sinusoïdal ou non).
✧ Le cas d’un régime permanent sinusoïdal est celui que nous connaissons bien, celui pour lequel nous
nous précipitons sur la notation complexe.
✧ Dans le cas du régime permanent stationnaire (RPS), l’équation de diffusion devient

∂2T
−λ =P
∂x2
✧ Et s’il n’y a pas de source

∂2T
=0
∂x2

Dans le cas d’un régime permanent stationnaire, le champ de température dépend :

➜ de la conductivité (vs. la diffusivité) s’il y a des sources ;
➜ uniquement des conditions aux limites s’il n’y a pas de source.

© Matthieu Rigaut 32 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·3 – Équation de diffusion thermique

I·3·iv – équation de diffusion thermique en 3D

✧ L’équation de conservation de l’énergie en 3D s’écrit
∂T
ρc + div ~th = P
∂t
✧ La loi de Fourier s’écrit, elle
−−→
~th = −λ grad T
✧ En remplaçant dans l’équation de conservation de l’énergie et en « sortant » le λ parce que le matériau
est homogène nous arrivons à
 
∂T −−→ ∂T −−→
ρc + div −λ grad T = P ρc − λ div grad T = P
∂t ∂t
✧ Il s’agit là de l’équation de diffusion thermique mais nous allons la noter autrement.

L’équation de diffusion thermique en 3D s’écrit

∂T
ρc − λ△ T = P
∂t
➜ ρ est la masse volumique ;
➜ c est la capacité thermique massique ;
➜ λ est la conductivité thermique ;
➜ P est la production volumique d’énergie.

I·3·v – pause laplacien et nabla

✬ définition du laplacien

Le laplacien noté △ est défini par

−−→
△ = div grad

✧ Les puristes remarquerons la différence entre ∆ et △ :

➜ ∆ se lit « delta » et est une lettre majuscule grecque représentant en général une variation
(∆U, ∆Ec . . . ) ;
➜ △ se lit « laplacien » et est un triangle représentant un opérateur vectoriel à l’instar de gradient,
divergence. . .
✧ Malgré cela nous éviterons d’utiliser la notation ∆ pour des variations dans des problèmes où inter-
vient le laplacien car, il faut l’avouer, cela prête à confusion.
 −→− → −−→
~ , la définition est associée à
Ceci est la définition du laplacien scalaire, pour le laplacien vectoriel △
la formule bien connue rot rot = grad div − △ ~ .

✬ expression du laplacien en coordonnées cartésiennes

En coordonnées cartésiennes le laplacien de T (x,y,z) s’écrit

∂2T ∂2T ∂2T
△T = + + 2
∂x2 ∂y 2 ∂z

© Matthieu Rigaut 33 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·3 – Équation de diffusion thermique

✧ Nous savons calculer un gradient ? une divergence ? alors nous pouvons calculer un laplacien.
✧ Commençons par le gradient

−−→ ∂T ∂T ∂T
grad T = ~ux + ~uy + ~uz
∂x ∂y ∂z

−−→ ∂T
✧ La composante de grad T sur ~ux est donc (de même pour les autres composantes) ce qui donne
∂x
! ! !
∂ ∂T ∂ ∂T ∂ ∂T
△T = + +
∂x ∂x ∂y ∂y ∂z ∂z

✧ D’où le résultat.

✬ notation nabla
Y kesako ?
✧ Nabla c’est :
➜ une notation (surtout anglo-saxonne) ;
➜ un vecteur ;
−−→
➜ un super moyen mnémotechnique pour retenir les expressions de grad , div , △ en carté-
siennes.
✧ Insistons : le « truc » calculatoire avec nabla ne marche qu’avec les coordonnées cartésiennes, qui-
conque essaiera de l’utiliser avec des coordonnées cylindro-polaire ou sphérique ne pourra qu’échouer.

Par notation le vecteur nabla ∇ ~ est tel que :

−−→ ~ ; ~ · ~th ;
grad T = ∇T div ~th = ∇ △ T = ∇2 T

Y nabla donne le calcul en coordonnées cartésiennes

✧ Nous ne le rappelerons jamais assez : utiliser le vecteur nabla ne marche qu’avec les
coordonnées cartésiennes.
Dans le système de coordonnées cartésiennes, le vecteur nabla peut être vu comme
 

 ∂ 


 ∂x 

 

~ = ∂ 
∇ 

 ∂y 

 

 ∂ 

 ∂z 

−−→ ~ T car
✧ Nous avons alors bien grad T = ∇
© Matthieu Rigaut 34 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·3 – Équation de diffusion thermique

   

 ∂ 

 ∂T 
 

 ∂x 

 ∂x 
 
   

 ∂  × T =  ∂T 
  

 ∂y 

 ∂y 
 
   

 ∂ 

 ∂T 
 
 ∂z   ∂z 

~ · ~th car
✧ De même div ~th = ∇
   

 ∂  
  jth,x



 ∂x  
 


   
 ∂  
 ·  jth,y
 ∂jth,x ∂jth,y ∂jth,z
=
 + +

 ∂y  
 

 ∂x ∂y ∂z
   

 ∂  
  jth,z


 ∂z   

✧ Enfin △ T = ∇2 T car en élevant au carré les composantes de nabla, nous trouvons déjà

∂2 ∂2 ∂2
∇2 = + +
∂x2 ∂y 2 ∂z 2

✧ Et donc
!
2 2 2
2 ∂ ∂ ∂ 2 ∂2T ∂2T ∂2T
∇ T = 2
+ 2+ 2 ×T ∇ T = + + 2
∂x ∂y ∂z ∂x2 ∂y 2 ∂z

I·3·vi – nouveau flash-back

✬ passer de la 3D à la 1D
✧ Nous pouvons (et c’est fort heureux) retrouver l’équation de diffusion en 1D à partir de l’équation
de diffusion en 3D.
✧ En effet, en admettant que la température ne dépende que de x nous avons

∂2T ∂2T ∂2T ∂2T

△T = + + 2 △T =
∂x2 ∂y 2 |{z}
∂z ∂x2
|{z}
=0 =0

✧ Et ainsi l’équation de diffusion donne

∂T ∂T ∂2T
ρc −λ△ T = P ρc −λ 2 =P
∂t ∂t ∂x
© Matthieu Rigaut 35 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

✬ passer de la 1D à la 3D sur un cas déjà rencontré

✧ Nous pouvons maintenant imaginer ce que deviendra l’équation de propagation en 3D.
✧ En effet, le passage 3D −→ 1D est

∂2T
△ T −→
∂x2
✧ Nous pouvons alors imaginer que l’opération 1D −→ 3D est

∂2T
−→ △ T
∂x2
✧ Imaginons une équation de propagation 1D sur, par exemple, une surpression p

∂2p 1 ∂2p
=
∂x2 c2 ∂t2
✧ Cette équation devient, en 3D

1 ∂2p
△p=
c2 ∂t2

✬ une notation et un nom

✧ Nous pouvons réécrire l’équation de propagation sous la forme

1 ∂2p
△ p− =0
c2 ∂t2
✧ Comme il s’agit là d’une équation très fréquente, un vieil usage a introduit une notation spécifique,
peu utilisée de nos jours

p = 0
✧ Comme l’équation est l’équation de d’Alembert, l’opérateur  s’appelle le d’alembertien.
✧ Mais il existe le « laplacien » : △ .
✧ Eh bien l’équation suivante est l’équation de Laplace

△T =0

I·4 – Quelques solutions

I·4·i – méthode
✧ Intéressons-nous maintenant à quelques solutions de l’équation de diffusion.
✧ Pour résoudre l’équation de diffusion thermique, la méthode est simple.
✧ La première question à se poser est : « Le régime est-il stationnaire ou non ? »
✧ Si la réponse est « Oui, le régime est stationnaire » alors l’équation de diffusion devient une équation
uniquement spatiale qui, dans le cas unidimensionnel, est rarement difficile à résoudre.
✧ Si la réponse est « Non, le régime n’est pas stationnaire » alors :
➜ le régime est sinusoïdal forcé, nous pouvons sortir les complexes et résoudre sans difficulté
particulière ;
➜ le régime n’est pas sinusoïdal forcé et il n’y a aucune méthode magique de résolution.

© Matthieu Rigaut 36 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

I·4·ii – choc thermique – solution en régime transitoire en 1D

✬ situation
✧ Modélisons le contact brusque avec un thermostat de la manière suivante :
➜ à t = 0− un milieu occupant un demi-espace est à température uniforme ;
➜ à t = 0 il est mis en contact avec un thermostat de température différente.

milieu semi-infini thermostat milieu semi-infini

T0 T1 T (x,t)
~ux

t<0 t>0

✧ Analyse :
➜ ici nous avons un problème de diffusion ;
➜ il y a deux invariances par translation tant dans la géométrie que dans les contraintes donc la
température ne va dépendre que d’un axe ;
➜ il s’agit ici d’un problème non stationnaire et le régime n’est pas sinusoïdal forcé ;
➜ les grandeurs pertinentes seront donc a (régime non stationnaire) ainsi que T0 comme condition
initiale et T1 comme contrainte.
✧ Analyse technique :
➜ nous allons noter ~ux l’axe normal au plan de séparation ;
➜ nous n’avons pas de méthode particulière de résolution donc il faudra la trouver avec une
méthode particulière.
✧ Les deux conditions se traduisent par :
➜ T (x,t < 0) = T0 ;
➜ T (0,t > 0) = T1 .

✬ graphiquement
✧ Une résolution numérique donne les graphiques suivants.

✧ Les différentes courbes correspondent à des instants répartis linéairement (t = k × t0 avec k entier).
✧ Nous pouvons voir que la température diffuse effectivement vers les x croissants mais de plus en plus
lentement.
© Matthieu Rigaut 37 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

✬ solution et interprétation
✧ Admettons la solution suivante
Z √
x/(2 a t)
T0 − T1 2
T (x,t) = T1 + 2 × √ × e −u du
π 0

✧ Nous pouvons facilement vérifier la 2e condition aux limites

Z 0
T0 − T1 2
T (0,t) = T1 + 2 × √ × e −u du T (0,t) = T1
π 0
| {z }
=0

✧ Rappelons que (cf. cours de maths)

Z ∞ √
−u2 π
e du =
0 2
✧ Nous avons ainsi, pour x > 0 à t = 0 :

Z
√ T0 − T1 ∞
2
x/(2 a t) = ∞ T (x,0) = T1 + 2 × √ × e −u du T (x,0) = T0
π 0

✧ De même à x fixé, si t −→ ∞ cela donne

√ Z 0
T0 − T1 2
x/(2 a t) −→ 0 T (x,0) = T1 + 2 × √ × e −u du T (x,0) = T1
π 0

✧ Autrement dit, en attendant suffisamment longtemps, quel que soit le lieu, la température finira par
atteindre T1 .
x
✧ Enfin remarquons que la température ne dépend que de √ .
2 at
✧ De la sorte, en se fixant une température T ⋆ quelconque (entre T0 et T1 quand même) elle est atteinte
en L à t0 tel que

L √
√ = Cte L∝ t0
2 a t0
✧ Nous retrouvons
√ bien le fait que la diffusion est un phénomène lent puisque la distance est propor-
tionnelle à t.

I·4·iii – mur de maison – régime stationnaire 1D

✧ Il s’agit là d’une situation très classique à tel point que c’est presque l’archétype des phénomènes
de diffusion et c’est pourquoi nous y ferons référence régulièrement.

✬ mur simple
Y situation et analyses
✧ Considérons un mur de maison, suffisamment grand pour pouvoir négliger les effets de bords et donc
être considéré comme infini bien qu’il ne le soit pas.
© Matthieu Rigaut 38 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

intérieur de la mur atmosphère

maison
T0 T1

✧ Ce mur sépare deux milieux (infinis eux aussi) et de températures uniformes et constante.
✧ La question est :
➜ quelle est la répartition de température à l’intérieur du mur en régime permanent ?
➜ quelle est la puissance perdue à travers le mur en régime permanent ?
✧ Analyse physique :
➜ nous sommes face à un problème de diffusion (c’est évident dans ce chapitre mais il faut le
préciser devant un problème « inconnu ») ;
➜ la maison et l’atmosphère vont jouer le rôle de thermostats ;
➜ le mur n’est pas radioactif, n’est pas envahi de réaction chimique et n’est pas parcouru par un
courant électrique donc le terme de production est nul ;
➜ le régime permanent est atteint et les contraites sont constantes (cf. thermostat) donc le régime
permanent est un régime permanent stationnaire ;
➜ les grandeurs pertinentes sont λ (caractéristique matériaux en RPS), e et Smur (géométrie), T0
et T1 (contraintes).
✧ Analyse technique :
➜ nous allons repérer le problème comme sur le schéma ci-dessous ;
➜ comme il y a invariance sur ~uy et ~uz , la température s’écrit T (x).

intérieur de mur atmosphère

la maison
T0 T1
~ux
O e

Y résolution
✧ Dans le cas de la diffusion sans source et en régime permanent stationnaire, l’équation de diffusion
s’écrit

△T =0
✧ Cela donne dans le cas unidimentionnel

d2 T
(x) = 0 T (x) = α x + β
dx2
© Matthieu Rigaut 39 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

✧ α et β sont des constantes d’intégration qui dépendent des contraintes.

À l’échelle mésoscopique, la température est une fonction continue de l’espace.

✧ Nous avons donc

T (0+) = T (0−) T (0+) = T0 et T (0+ ) = α α = T0

✧ De même

T1 − T0
T (e− ) = T (e+ ) T (e− ) = T1 et T (e− ) = T0 + β e β=
e
✧ Finalement

T1 − T0
T (x) = T0 + ×x
e
✧ Nous voyons que le champ de température est affine à l’intérieur du mur.

T0 mur atmosphère

maison

T1
x
e

✧ Nous pouvons constater que λ n’intervient pas, mais c’était prévisible étant donné que la conductivité
n’intervient pas dans l’équation de diffusion.

Y puissance perdue à travers le mur

✧ La puissance thermique perdue est celle qui va de la maison vers l’atmosphère ce qui nous impose le
~ dans le sens de +~ux .
sens de dS
✧ La puissance thermique perdue s’écrit donc
ZZ
Pperdue = Φth,à travers le mur Pperdue = ~P
~th (P ) · dS
P ∈mur

✧ Or la densité surfacique de courant en volume s’écrit, d’après la loi de Fourier

−−→ dT T1 − T0
~th = −λ grad T (x) ~
~th (x) = −λux ~th (x) = −λ ~ux
dx e
✧ De plus pour compter positivement l’énergie qui va de la maison vers l’atmosphère, il faut

~P = +dSP ~ux
dS
✧ Nous avons ainsi
© Matthieu Rigaut 40 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

ZZ ZZ
T1 − T0 T0 − T1
Pperdue = −λ dSP Pperdue =λ × dSP
P ∈mur e e P ∈mur

✧ Et finalement

λ Smur
Pperdue = (T0 − T1 ) ×
e
✧ Ce résultat est bien homogène (pour le voir, utiliser de tête la loi de Fourier) et est cohérent :
➜ si T0 > T1 c’est la maison qui perd de l’énergie ;
➜ plus e est grand plus la perte est faible ;
➜ plus le mur est grand plus la perte est grande.

✬ mur isolé
Y situation et analyses
✧ Considérons maintenant le même mur (sans effet de bord et d’épaisseur e) recouvert d’une couche
isolante d’épaisseur ℓ et de conductivité λi .

intérieur de mur atmosphère

isolant

la maison
T0 T1
~ux
O ℓ e+ℓ

✧ Posons-nous les mêmes questions :

➜ champ de température dans le mur et l’isolant ?
➜ puissance perdue ?
✧ L’analyse physique est analogue à la précédente
➜ c’est un problème de diffusion ;
➜ il n’y a pas de production ;
➜ il y a deux milieux diffusants donc deux équations de diffusion à résoudre et à recoller ;
➜ les grandeurs pertinentes sont les mêmes que précédemment auxquelles il faut ajouter les ca-
ractéristiques de l’isolant λi et ℓ.
✧ Analyse technique :
➜ comme il y a deux milieux diffusants, il y a deux fonctions températures soit Ti (x) et Tm (x) ;
➜ comme les deux milieux sont en régime permanent stationnaire, nous pouvons écrire tout de
suite que les champs de température associés sont linéaires ce qui donne

Ti (x) = A x + B et Tm (x) = C x + D
➜ quatre constantes d’intégration il faudra donc quatre conditions aux limites pour pouvoir
tout résoudre.
© Matthieu Rigaut 41 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

✧ En ce qui concerne la continuité de température, nous en avons 3 : en x = 0, en x = ℓ et en x = ℓ + e.

✧ Il manque une condition, ça sera la continuité du flux thermique en x = ℓ.
✧ N’oublions pas en effet qu’en régime permanent stationnaire la conservation de l’énergie s’écrit

∂jth
0+
=P
∂x
✧ Cette équation montre bien que la fonction jth est dérivable donc qu’elle est continue.

À l’échelle mésoscopique, la densité surfacique de courant thermique en volume est une

fonction continue de l’espace.

✧ Nous avons nos 4 conditions aux limites, nous pouvons les traduire puis résoudre.

Y traduction des conditions aux limites

✧ La continuité de la température en x = 0 donne

T (0−) = Ti (0+) T0 = B
✧ La continuité de la température en x = ℓ + e donne

Tm ((ℓ + e)− ) = T ((ℓ + e)+ ) C (ℓ + e) + D = Text

✧ La continuité de la température en x = ℓ donne

Ti (ℓ− ) = T (ℓ+) Aℓ + B = C ℓ + D
✧ Enfin la continuité du courant thermique donne

jth,i (ℓ− ) = jth,m (ℓ+ )

✧ Soit, avec la loi de Fourier

d2 Ti − d2 Tm +
−λi (ℓ ) = −λm (ℓ ) −λi A = −λm C λi A = λm C
dx2 dx2

Y résolution
✧ Commençons par réécrire les équations obtenues


 C (ℓ + e) + D = Text

B = T0 A ℓ + T0 = C ℓ + D

 λi A = λm C


✧ Nous allons d’abord chercher le flux thermique perdu.

✧ Par continuité du flux thermique, le flux thermique perdu par la maison n’est autre que le flux qui
traverse l’épaisseur ℓ d’isolant.
✧ Or ce flux s’écrit

Φi = −λi × A × Smur
✧ Il nous faut donc trouver A.
✧ Pour ce faire commençons par éliminer C et obtenir un système de deux équations à deux inconnues
© Matthieu Rigaut 42 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·4 – Quelques solutions




 λi
 λm A (ℓ + e) + D = Text


λi
C= A
λm 
 λi


 A ℓ + T0 = Aℓ + D
 λm

✧ Réécrivons ce système sous une forme plus adéquate




 λi


 A (ℓ + e) + D = Text
 λm

!

 λ i


 − 1 ℓ A + D = T0
 λm



✧ En soustrayant les deux équations ainsi obtenues nous arrivons à

 ! 
λi λi✁ Text − T0
A (✁ℓ + e) − ✁ − 1 ℓ = Text − T0 A= λi
λm λ
✁m e+ℓ
λm

✧ Graphiquement, cela donne quelque chose comme représenté ci-dessous.

T0 mur atmosphère
isolant

maison

T1
x
ℓ e+ℓ

Y isolation ?
✧ La puissance perdue par la maison peut s’écrire

Text − T0
Pperdue = −λi A Smur Pperdue = −λi Smur × λi
λm
e+ℓ
✧ Réécrivons cette puissance sous une autre forme pour pouvoir plus facilement la comparer à la
puissance perdue par le mur non isolé

Text − T0
Pperdue = −λm Smur ×
e + λλm ℓ
i

✧ Résultat à comparer avec

Text − T0
Pperdue, non isolé = −λm Smur ×
e
© Matthieu Rigaut 43 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

✧ Donc finalement, isoler un mur de conductivité λm avec un isolant d’épaisseur ℓ et de conductivité

λm
λi revient à augmmenter l’épaisseur du mur de ℓ
λi
✧ Nous retrouvons là quelque chose de très cohérent vis-à-vis de l’isolation thermique :
➜ plus l’épaisseur est grande, plus l’isolation est efficace ;
➜ moins la conductivité de l’isolant est grande, plus l’isolation est efficace ;
➜ plus la conductivité du mur est grande, plus l’isolation est efficace.

I·4·iv – barreau radioactif – régime stationnaire en symétrie cylindrique

✬ dispositif et analyse
✧ Considérons un barreau radioactif cylindrique de rayon R et de hauteur h.

✧ Nous considérerons le barreau suffisamment grand pour pouvoir négliger les effets de bord.
✧ En admettant que la température ambiante soit uniforme et constante égale à T0 , quelle est la
température dans le barreau ?
✧ Analyse physique
➜ ici nous sommes face à un problème de diffusion ;
➜ comme le barreau est radioactif, il y a un terme de source que nous prendrons uniforme et
constant à P0 (en W.m−3 ) ;
➜ il y a une invariance par translation par rapport à l’axe (puisque les effets de bord sont négligés)
et par rotation autour du même axe ;
➜ nous allons considérer le régime permanent stationnaire ;
➜ les grandeurs pertinentes sont R et h (géométrie), λ, P0 (aspect thermique) et T0 (contrainte).
✧ Rappelons que ρ et c ne font pas parti des grandeurs pertinentes car le régime est stationnaire.
✧ Analyse technique
➜ étant donné la géométrie du dispositif, nous allons naturellement utiliser le repérage cylindro-
polaire ce qui permet de dire que, grâce aux invariances, les différentes grandeurs (T , jth ) ne
dépendent que de r ;
➜ dans le cas d’un régime stationnaire avec terme de source (de préférence uniforme), il y a deux
grandes méthodes ;
➙ une mise en équation globale qui permet de simplifier la résolution grâce aux symétries ;
➙ une mise en équation en 3D quelconque suivi de la résolution avec le laplacien en cylindro-
polaire.

© Matthieu Rigaut 44 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

✬ mise en équation : méthode du bilan global

Y un système bien choisi
✧ L’idée est en fait de redémontrer la loi de conservation de l’énergie mais dans le cas particulier du
barreau radioactif.
✧ Le tout étant de bien choisir le système.
✧ Pour choisir le « bon » système, c’est très simple car cela revient à choisir un « bon » volume et une
« bonne » surface qui permettent de faire des calculs de flux facilement.
✧ Et des calculs de flux, cela ressemble beaucoup à l’électrostatique !
✧ En fait, comme pour le théorème de Gauss le bon système sera celui délimité par une surface qui :
➜ respecte les symétries du problème ;
➜ passe par le point M où nous voulons déterminer le champ (de température).
✧ Ici comme la diffusion est radiale nous allons choisir un système contenue dans un cylindre de rayon
r et de hauteur h.
S1

r
r
M M
h h M

Slat

S2 vue de dessus

vue en perspective vue en coupe

Y mise en équation et résolution

✧ Nous allons faire un bilan d’énergie pour ce système entre t et t + dt.
✧ Ce bilan s’écrit

VARIATION dans le temps = ÉCHANGE à travers la surface + CRÉATION en volume

✧ Comme le régime est stationnaire, nous avons tout de suite

VARIATION dans le temps = 0

✧ Commençons par calculer le terme d’échange.
ZZ
Φreçue =− ~P
~th (P ) · dS
P ∈S

✧ En découpant la surface S en trois (1 étant au-dessus et 2 en dessous) et avec l’extensivité du flux

ZZ ZZ ZZ
Φreçue = − ~
~th (P ) · dS1P − ~
~th (P ) · dS2P − ~lat,P
~th (P ) · dS
P ∈S1 P ∈S2 P ∈Sℓ

✧ Or les vecteurs surface élémentaires s’écrivent

© Matthieu Rigaut 45 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

~1P = dSP ~uz ;

dS ~2P = −dSP ~uz
dS et ~lat,P = dSP ~ur
dS
✧ Ce qui donne, avec ~th (P ) = jth (P ) ~ur
ZZ
Φreçue = −0 − 0 − jth (P ) dSP
P ∈Sℓ

✧ Or la surface S est choisie de telle sorte que

rP = Cte = rM = r
not

✧ Et ainsi

ZZ
Φreçue = − jth (rP ) dSP
P ∈Sℓ
ZZ
=− jth (r) dSP
P ∈Sℓ
ZZ
= −jth (r) × dSP
P ∈Sℓ
= −jth (r) × 2 π r h

✧ Le terme de production s’écrit, lui

ZZZ
P= P (P ) δτP
P ∈V

✧ Comme P (P ) = P0 est uniforme, nous avons rapidement

ZZZ
P= P0 δτP
P ∈V
ZZZ
= P0 δτP
P ∈V
= P0 × π r 2 h

✧ En réinjectant ces résultats dans la loi de conservation de l’énergie cela donne

P0 r
0 = Φreçue dt + P dt jth (r) × 2 π r h = P0 × π r 2 h jth (r) =
2
✧ Maintenant nous pouvons utiliser la loi de Fourier qui dit

dT
jth (r) = −λ (r)
dr
✧ Et ainsi T (r) s’écrit

dT P0 r dT P0 r P0 r 2
−λ (r) = (r) = − T (r) = − + Cte
dr 2 dr 2λ 4λ
✧ La constante d’intégration s’obtient grâce à la continuité de la température en r = R

P0 2
T (R) = T0 T (r) = (R − r 2 ) + T0
4λ
© Matthieu Rigaut 46 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

Y vision graphique
✧ Graphiquement la température en fonction de r ressemble à cela

caloporteur
T0
barreau
radioactif
r
R

✧ Nous pouvons constater que :

➜ la température la plus élevée est au centre ;
➜ plus r augmente, plus le courant thermique est grand.

✬ mise en équation : méthode sans réfléchir

Y repartir de l’équation de diffusion thermique
✧ L’idée est de résoudre « mathématiquement » l’équation de diffusion.
✧ Il faut d’abord rétablir dans le cas général l’équation

∂T
ρc − λ△ T = P
∂t
✧ Cette équation se simplifie tout de suite car le régime est permanent stationnaire

P0
−λ △ T = P △T =−
λ
✧ Et maintenant il faut l’expression du laplacien en cylindro polaire pour une fonction qui ne dépend
que de r.
✧ Pour cela en cherchant dans les archives / l’énoncé / la calculatrice, nous trouvons
!
1 d dT
△ T (r) = × r× (r)
r dr dr

! Comme le montre l’expression ci-dessus, contrairement aux coordonnées cartésiennes

d2 T
△ T (r) 6= (r) en cylindro-polaire
dr 2

Y résoudre
✧ Il faut y aller lentement et par étapes.

! !
P0 1 d dT P0 d dT P0
△T =− × r× (r) = − r× (r) = − ×r
λ r dr dr λ dr dr λ

✧ En primitivant cela donne

© Matthieu Rigaut 47 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

dT P0 dT P0 α
r× (r) = − × r2 + α (r) = − ×r+
dr 2λ dr 2λ r
✧ En primitivant une fois de plus nous arrivons à

P0
T (r) = − × r 2 + α ln r + β
4λ
✧ Il y a deux constantes d’intégration, il faut donc deux conditions aux limites pour les trouver.
✧ Ici r = 0 existe car il est dans le domaine de résolution, il faut donc trouver une solution pour r = 0.
✧ Cela implique que α = 0 sans quoi pour r = 0 la température deviendrait infinie, ce qui n’est pas
possible.
✧ Il reste donc

P0
T (r) = − × r2 + β
4λ
✧ À partir de là nous sommes ramenés au cas précédent et la continuité de la température en r = R
permet de conclure.

✬ comparaison des deux méthodes

✧ Vue de loin la 2e méthode paraît plus courte mais c’est oublier que nous n’avons pas retranscrit :
➜ le début qui consiste à retrouver l’équation de diffusition thermique en 3D ;
➜ la fin pour trouver β et tracer T (r).
✧ Malgré cela cette méthode présente quelques caractéristiques qui peuvent être des inconvénients :
➜ cette méthode dépend de la connaissance de △ T (r) qui n’est pas intuitive ;
➜ cette méthode fait apparaître des arguments de cohérence dans la résolution ce qui n’est pas
fréquent ;
➜ des fois il est possible d’avoir un terme en ln r dans une diffusion en symétrie cylindrique.
✧ Considérons un fil de cuivre parcouru par un courant et entouré d’une gaine isolante (électriquement).

✧ La diffusion et radiale et il y a deux milieux diffusif donc il y a deux fonctions T (r) à trouver

Tfil (r) pour 0 6 r 6 R et Tisol (r) pour R 6 r 6 R + a

✧ Dans ces conditions il ne peut pas y avoir de terme en ln r dans Tfil (r) puisque la valeur r = 0
appartient au domaine de résolution mais il peut très bien y avoir ce terme dans Tisol (r) parce que,
justement, r = 0 n’est pas un point de l’isolant.

I·4·v – température dans le sol – régime sinusoïdal forcé en 1D

✬ situation et analyse physique
✧ Considérons une portion de sol suffisamment grande pour pouvoir être considérée comme infinie.
© Matthieu Rigaut 48 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

atmosphère

sol

✧ Nous pouvons mettre de côté le chauffage par rayonnement car la diffusion est un phénomène li-
néaire.
✧ Nous pouvons ainsi nous concentrer sur l’influence thermique de l’air.
✧ La température de l’air varie sinusoïdalement avec des périodes très différentes :
➜ une période correspond à l’évolution saisonnière ;
➜ une autre période correspond à l’évolution quotidienne.
✧ Toujours parce que la diffusion est un phénomène linéaire nous allons nous contenter d’étudier une
seule période T sans préciser laquelle.
✧ Ainsi, pour nous, la température de l’air sera Tair (t) = T0 + θ0 cos (ω t)

Tatmo Tatmo T

θ0

t t

✧ Analyse physique
➜ nous sommes face à un problème de diffusion ;
➜ le sol n’étant pas radioactif 5 , ni le siège de réaction chimique, ni parcouru par un courant, la
diffusion se fera sans terme de production ;
➜ le régime que nous allons chercher est un régime permanent mais non stationnaire ;
➜ la situation est invariante par translation suivant deux axes ;
➜ les grandeurs pertinentes sont a pour l’aspect thermique (et non λ), T0 , θ0 et ω pour les
contraintes. Ici il n’y a pas d’aspect géométrique à moins de désirer compter exactement le flux
qui passe à travers une surface S0 dessinée sur le sol.

✬ analyse technique
✧ Commençons par imposer un repérage de la situation.

atmosphère

sol
~uz

5. Du moins pas à faible profondeur.

© Matthieu Rigaut 49 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

✧ Les invariances imposent donc que la température s’écrit T (z,t) avec z > 0.
✧ Nous pouvons écrire cette température sous la forme

T (z,t) = T0 + Tsol (z,t)

✧ L’avantage de cette écriture et qu’elle fait apparaître le régime permanent stationnaire (i.e. la solution
pour θ0 = 0) de sorte que Tsol (z,t) soit sinusoïdal de valeur moyenne nulle.
✧ Sinusoïdal ? De valeur moyenne nulle ? Utilisons la notation complexe.

Tsol (z,t) = Re(Tsol (z,t))

✧ La notation complexe donne de manière systématique

Tsol (z,t) = f (z) e j ω t

✧ Nous nous irons un peu plus loin en cherchant directement Tsol (z,t) sous une forme un peu plus
particulière

Tsol (z,t) = Tm e j (ω t−k z)

✧ Si cela ne marche pas, nous reviendrons à la forme originelle Tsol (z,t) = f (z) e j ω t mais cela marchera.

✬ résolution
✧ Faisons comme si nous avions retrouvé l’équation de diffusion thermique sans terme de production

∂T 1 ∂T not λ
ρc −λ△ T = 0 △T = avec a=
∂t a ∂t ρc
✧ En une dimension le laplacien donne

∂2T 1 ∂T
(z,t) = (z,t)
∂z 2 a ∂t
✧ Et ainsi en notation complexe

∂2T 1 ∂T
2
(z,t) = (z,t)
∂z a ∂t
✧ En injectant la solution recherchée cela donne

∂2   1 ∂   1
Tm e j (ω t−k z)
= Tm e j (ω t−k z)
−k 2 Tm e j (ω t−k z) = j ω Tm e j (ω t−k z)
∂z 2 a ∂t a

✧ Nous avons finalement

1 ω
−k 2 = jω k2 = − j
a a
✧ Pour trouver k nous allons nous aider de la relation −j = e −j π/2
© Matthieu Rigaut 50 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

1/2
k = k2
 1/2
ω −j π/2
= e
a
r
ω −j π/4
=± e
a
r
ω 1−j
=± × √
a 2

✧ Résultat que nous allons noter

r
ω
avec et
not
k = ε × k0 (1 − j) ε = ±1 k0 =
2a
✧ Nous pouvons maintenant remplacer k dans la solution

Tsol (z,t) = Tm e j (ω t−k z)

= Tm e j (ω t−ε k0 z (1−j))
2
= Tm e j (ω t−ε k0 z) × e +j ε k0 z

= Tm e j (ω t−ε k0 z) × e −ε k0 z

✧ Ici comme z peut tendre vers +∞ et que la température ne peut être infinie, nous avons forcément
ε = +1 et donc

Tsol (z,t) = Tm e j (ω t−k0 z) × e −k0 z

✧ Nous pouvons revenir à la solution réelle (sans oublier l’argument de Tm )

 
Tsol (z,t) = Re Tsol (z,t) Tsol = Tm cos(ω t − k0 z + ϕ) × e −k0 z

✧ Puis à la solution complète

T (z,t) = T0 + Tm cos (ω t − k0 z + ϕ) × e −k0 z

✧ Pour déterminer Tm et ϕ nous allons utiliser la continuité de la température en z = 0.

T (0,t) = T0 + θ0 cos (ω t) T0 + Tm cos (ω t + ϕ) = T (0,t) = T0 + θ0 cos (ω t)

✧ Par identification nous trouvons très vite

Tm = θ0 et ϕ=0

✧ Et la solution finale s’écrit

T (z,t) = T0 + θ0 cos (ω t − k0 z) × e −k0 z

© Matthieu Rigaut 51 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

✬ interprétation
✧ Graphiquement la solution en fonction de x à t fixé est la suivante (pour plusieurs valeurs de t)

T T T
θ0

T0 T0 T0

z z z

✧ En revanche, la solution en fonction de t à x fixé est la suivante (pour plusieurs valeurs de x)

T T T
θ0 θ0 θ0

T0 T0 T0

t t t

✧ Dynamiquement cela ressemble à une onde dont l’amplitude diminue très vite au fur et à mesure de
l’avancement.

Montrer l’animation Maple.

✧ Ici il n’y a aucune dissipation ! L’énergie se conserve et n’est même pas transformée en une autre
forme.
✧ La raison pour laquelle l’amplitude décroît est simple : comme la diffusion se fait dans tous les sens,
un « morceau » d’énergie peut aller vers l’avant ou vers l’arrière.
✧ C’est pourquoi en chaque endroit il y a de l’énergie qui rebrousse chemin, ce qui implique que moins
d’énergie va vers l’avant et donc une amplitude des oscillations de température qui décroît. C’est
inévitable.
✧ Nous voyons aussi que l’effet des variations de température se fait sentir jusqu’à une profondeur
5
d’environ soit
k0
r r
2a a × période
zmax ∼ 5 zmax ∼ 5
ω π
✧ Nous voyons ainsi que la profondeur maximal est directement lié à la période.
✧ En Champagne-Ardenne les caves de vin de champagne sont à 30 m de profondeur et creusées dans un
sol calcaire. En prenant une diffusivité de 0,4.10−6 m2 .s−1 et une période de 365 jours, nous trouvons

zmax = 10 m
© Matthieu Rigaut 52 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

✧ Cette valeur explique pourquoi, quelle que soient les conditions extérieures de température, il fait
toujours 13 ◦ C dans une cave de champagne.
✧ Ci-dessous, une vue de la superbe cave de Moët et Chandon à Epernay 6 .

I·4·vi – résumé des conditions aux limites

✬ condition au bord d’un milieu diffusant

Y paroi adiabatique

✧ Une paroi adiabatique est une paroi qui ne laisse pas passer l’énergie.
✧ Nous pouvons donc traduire techniquement la condition en utilisant la densité de courant en volume.

Au niveau d’une paroi adiabatique, nous avons

~th · ~n = 0 où ~n est le vecteur unitaire normal à la surface

~n
~th
milieu diffusant ~th

Y contact avec un thermostat

✧ Un thermostat un dispositif qui impose sa température donc

6. Source : [Link]
[Link].

© Matthieu Rigaut 53 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

La température de surface d’un milieu en contact avec un thermostat est imposée par le
thermostat.
Tmilieu diffusant (à la surface) = Tthermostat

thermostat à T0

milieu diffusant T0

Y contact avec un fluide en mouvement

✧ Tout le monde a déjà fait l’expérience du ventilateur : quand il y a circulation de fluide près d’un
milieu diffusant, il y a des pertes thermiques.
✧ Ces pertes sont dites conducto-convectives (puisqu’il s’agit d’un échange entre un milieu diffusif et
un milieu convectif) et elles s’écrivent simplement.

Entre un milieu diffusif et un milieu fluide en mouvement, le flux thermique dit

conducto-convectif s’écrit


Φ = ±h S Tfluide − Tsurface où :
➜ h dépend du matériau, du fluide et de la vitesse du fluide ;
➜ S est la surface d’échange ;
➜ Φ est en watt.

fluide en
mouvement
Φ

milieu diffusant

✧ Cette loi n’est pas exigible en tant que telle mais comme elle ressemble énormément à une loi que
nous allons voir dans le paragraphe suivant, autant la connaître de manière à ne pas les confondre.
✧ En ce qui concerne le signe, il faut réfléchir suivant la convention adoptée (qui reçoit Φ ?).

✬ condition entre deux milieux diffusants

Y flux thermique

✧ Généralisons une loi vue lors de la recherche du champ de température dans un mur.
© Matthieu Rigaut 54 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

Entre deux milieux diffusifs, il y a toujours continuité du flux thermique.

Φqui sort de 1 = Φqui rentre dans 2 ou ~th,1 · ~n12 = ~th,2 · ~n12

milieu diffusant 1

~n12
~th,1
milieu diffusant 2
~th,2

✧ Cette loi ne fait que traduire la conservation de l’énergie.

✧ Pour la mettre en défaut, il faudrait qu’il y ait une création (ou une absorption) d’énergie à l’échelle
de l’interface. Ce qui n’est pas possible à l’échelle mésoscopique, mais à l’échelle macroscopique. . .

Y température à l’échelle mésoscopique

✧ Comme nous l’avons déjà signalé sur l’exemple du mur

À l’échelle mésoscopique, la température et la densité surfacique de courant thermique

en volume sont des champs continus de l’espace.

Y température à l’échelle macroscopique

✧ À l’échelle macroscopique, la température n’est pas toujours une fonction continue de l’espace.
✧ Prenons par exemple une maison vue de loin.

atmosphère

Tmaison Tatmo

✧ À cette échelle les murs sont d’épaisseur nulle ce qui implique, à cette échelle, que

T (juste au bord d’un mur) 6= T (juste de l’autre côté du mur)

✧ En fait tout se passe comme si l’épaisseur e avait été réduite à 0.
maison mur atmosphère maison atmosphère

e e

© Matthieu Rigaut 55 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·4 – Quelques solutions

✧ Il peut en être de même pour le terme de production.

✧ Imaginons une paroi source d’énergie et regardons-la de loin.

e e

✧ Nous voyons que nous pouvons alors introduire un terme de puissance surfacique (et non plus volu-
mique).
✧ C’est exactement la même chose qu’en électromagnétisme lorsque nous avons introduit la densité
surfacique de charge.
✧ Ces pourquoi, face à une épaisseur infiniment fine à l’échelle macroscopique :
➜ la température n’est plus continue ;
➜ le courant thermique en volume n’est plus continu si l’épaisseur est source d’énergie.
✧ Il devient alors aussi plus difficile de trouver les flux thermiques, même dans le cas où il n’y a pas de
source.
✧ Reprenons ainsi le mur de maison mais en version « infiniment » fin.

intérieur de
atmosphère
la maison
Tmaison Tatmosphère

✧ Nous savons (car nous l’avons déjà calculé) que le flux thermique qui passe de l’un à l’autre s’écrit,
en régime stationnaire

λ not λ
Φmaison vers atmosphère = × S (Tmaison − Tatmosphère ) = h S (Tmaison − Tatmosphère ) avec h=
e e
✧ Nous voyons alors réapparaître une loi formellement identique à celle qui donne le flux conducto-
convectif mais qui est fondamentalement différente en terme phénoménologique car là ce flux n’est
autre que celui qui traverse un milieu diffusif (ici, le mur).
✧ Cette loi montre aussi que le flux est proportionnel à l’écart en température. Cela nous amène
directement au paragraphe suivant.
© Matthieu Rigaut 56 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·5 – Résistance thermique

I·5 – Résistance thermique

I·5·i – définition par l’exemple
✬ tout d’abord les circonstances
✧ Pour parler de résistance thermique il faut quelques conditions (simples).

Nous pourrons parler de résistance thermique dans un phénomène de diffusion :

➜ en régime permanent stationnaire ;
➜ sans terme de source.

✧ Reste maintenant à savoir ce qu’est cette résistance thermique.

✬ exemple du mur simple

✧ Retrouvons très vite les résultats trouvés sur le mur simple de maison.
✧ Nous savons que la température est une fonction continue de l’espace et que la variation est affine
dans le mur car le régime est permanent et sans terme de source.
✧ Graphiquement cela donne donc ceci

Tmaison mur atmosphère

maison

Tatmosphère
x
e

✧ Analytiquement il n’est pas très difficile de retrouver la fonction

Tatmo − Tmaison
T (x) = Tmaison + ×x
e
✧ Le flux allant de la maison vers l’atmosphère s’écrit donc

Φm→a = jth × S

✧ Avec la loi de Fourier cela donne

dT Tmaison − Tatmo λS
jth = −λ jth = Φm→a = (Tmaison − Tatmo )
dx e e
✧ Autrement dit en généralisant un peu avec des milieux 1 et 2

Φ1→2 = qqch × (T1 − T2 ) T1 − T2 = autre chose × Φ1→2

✧ C’est cette « autre chose » qui est la résistance thermique.

© Matthieu Rigaut 57 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·5 – Résistance thermique

✬ définition retenue

Un milieu diffusant dont les limites sont aux températures T1 et T2 :

T1 − T2
➜ a une résistance thermique Rth telle que Rth = ;
Φ1→2
Φ1→2
➜ a une conductance thermique Gth telle que Gth = .
T1 − T2
Rth et Gth sont indépendants de Φ1→2 et T1 − T2

✧ Pour se souvenir des conventions, c’est facile : les chiffres représentant les milieux sont dans le même
ordre dans Φ1→2 et T1 − T2 .
✧ Sinon, avec un peu de bon sens, le signe se retrouve par cohérence en sachant que le flux thermique
est dirigé des zones chaudes vers les zones froides.

La résistance thermique d’un milieu diffusif unidimensionnel de longeur ℓ, de section S

et de conductivité λ s’écrit
1 ℓ
Rth =
λS

L L

✧ Cette expression ressemble furieusement à celle de la résistance électrique !

✧ En fait ce n’est pas une coïncidence du tout.

✬ interprétation
✧ Reprenons deux milieux quelconques séparés par un milieu diffusant.

thermostat thermostat
milieu diffusant
T1 T2

✧ Nous avons alors

T1 − T2
Φ1→2 =
Rth
✧ Ainsi, si nous voulons favoriser les transferts thermiques, il nous faut diminuer la résistance thermique
alors que si nous voulons au contraire limiter cette diffusion, il faut l’augmenter.

I·5·ii – analogie électrique

✧ Nous pouvons faire une analogie entre les aspects électriques et les aspects thermiques conformément
à ce tableau.
© Matthieu Rigaut 58 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·5 – Résistance thermique

diffusion thermique courant électrique

~th ~
Φ1→2 i1→2
T1 − T2 V1 − V2
1 2 1 2
i1→2

T1 − T2 = Rth Φ1→2 V1 − V2 = R i1→2

✧ De cette manière il est assez facile de se souvenir de la définition de la résistance termique

U = R i1→2 T1 − T2 = Rth Φ1→2

Il est possible de transposer un problème de diffusion en un problème électrique dès lors

que :
➜ le régime est permnanent stationnaire ;
➜ il n’y a pas de terme de source.

✧ Au niveau des lois fondamentales, il est aussi possible de faire des analogie

diffusion thermique courant électrique

−−→ −−→
~th = −λ grad T (loi de Fourier) ~ = −γ grad V (loi d’Ohm)
ZZ
R
Φ1→2 = ~P
~th (P ) · dS ~P
i1→2 = P ∈S ~(P ) · dS
P ∈S

✧ Et comme nous connaissons l’expression de la résistance d’un milieu conducteur de longueur ℓ et de

conductivité γ depuis le cours (de première année) sur le courant électrique, nous avons

1 ℓ 1 ℓ
R= × Rth = ×
γ S λ S

I·5·iii – association de résistances thermiques

✬ association parallèle
Y résultat

Deux milieux diffusifs A et B sont en parallèle lorsque leurs extrémités sont en contact
avec les mêmes thermostats.

Rth,A
T1 T2
Rth,B

© Matthieu Rigaut 59 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·5 – Résistance thermique

Deux milieux A et B diffusifs en parallèle sont équivalents à un seul milieu diffusif de

résistance thermique Rth,éq tel que
1 1 1
= +
Rth,éq Rth,A Rth,B

Rth,A
T1 T2 T1 Rth,éq T2
Rth,B

Y démonstration
✧ Considérons deux milieux en parallèle et demandons-nous s’ils ne seraient pas équivalents (thermi-
quement parlant) à un seul milieu diffusif.

ΦA Rth,A
T1 T2
ΦB Rth,B

✧ L’extensivité du flux nous dit que

Φtot = ΦA + ΦB

✧ Par définition des résistances thermiques Rth,A et Rth,B nous avons ainsi
!
T1 − T2 T1 − T2 1 1
Φtot = + Φtot = + × (T1 − T2 )
Rth,A Rth,B Rth,A Rth,B

✧ Or pour que l’association soit équivalente à une résistance thermique unique, il faut trouver Rth,éq
telle que

1
Φtot = × (T1 − T2 )
Rth,éq

✧ Par identification, nous trouvons bien le résultat

1 1 1
= +
Rth,éq Rth,A Rth,B

✬ association série

© Matthieu Rigaut 60 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·5 – Résistance thermique

Y résultat

Deux milieux diffusifs A et B sont en série lorsqu’ils sont mis bout à bout.

T1 Rth,A Rth,B T2

Deux milieux A et B diffusifs en série sont équivalents à un seul milieu diffusif de

résistance thermique Rth,éq tel que
Rth,éq = Rth,A + Rth,B

T1 Rth,A Rth,B T2 T1 Rth,éq T2

Y démonstration
✧ Considérons deux milieux en série et demandons-nous s’ils ne seraient pas équivalents (thermiquement
parlant) à un seul milieu diffusif.

T′

T1 Rth,A Rth,B T2

✧ En notant T ′ la température au niveau du contact entre A et B, nous pouvons écrire la différence

des températures sous la forme

T1 − T2 = (T1 − T ′ ) + (T ′ − T2 )
✧ Avec les définitions des résistances thermiques Rth,A et Rth,B nous avons ainsi

T1 − T2 = Rth,A × ΦA + Rth,B × ΦB
✧ Et la continuité du flux impose que

ΦA = ΦB = Φtot T1 − T2 = (Rth,A + Rth,B ) × Φtot

✧ Or pour que l’association soit équivalente à une résistance thermique unique, il faut trouver Rth,éq
telle que
1
Φtot = × (T1 − T2 )
Rth,éq
✧ Par identification, nous trouvons bien le résultat

Rth,éq = Rth,A + Rth,B

© Matthieu Rigaut 61 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·5 – Résistance thermique

✬ exemple du mur isolé

✧ Nous pouvons retrouver très vite le flux traversant un mur isolé.
✧ Reprenons le dispositif.

intérieur de mur atmosphère

isolant
la maison
T0 T1

ℓ e

✧ Nous avions trouvé (après « pas mal » de calculs)

Tm − Ta Tm − Ta
Φm→a = λm S × Φm→a =
λm e ℓ
e+ ℓ +
λi λm S λi S

✧ Transposons le problème en problème électrocinétique

Tmaison Tinterface Tatmosphère

Risolant Rmur

✧ Nous avons donc immédiatement

Tm − Ta e ℓ
Φm→a = avec Rth,mur = et Rth,isol =
Rth,mur + Rth,isol λm S λi S

✧ Ce qui n’est autre que le résultat. . . et « immédiat » !

I·5·iv – résistance cylindrique

✬ situation, analyse
✧ Intéressons-nous à une situation de diffusion thermique, mais en symétrie cylindrique, comme par
exemple l’isolation d’un tuyau d’eau chaude par une gaine en mousse.
© Matthieu Rigaut 62 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·5 – Résistance thermique

✧ Ici nous avons bien un milieu diffusif (la gaine isolante) qui ne présente pas de terme de création
d’énergie et qui est au contact de deux thermostats (la fine couche de cuivre elle même thermostatée
par l’eau chaude à la température T1 et l’atmosphère à la température T0 ).
✧ Dans ces conditions, nous pouvons associer une résistance thermique mais comme ce n’est pas une
ℓ
diffusion unidimensionnelle, nous n’avons pas R = et ce d’autant plus qu’il aurait été difficile de
λS
dire ce qu’était la section du milieu.
✧ Quelle est donc l’expression de la résistance thermique de cet isolant ?
✧ Analysons le problème :
➜ le dispositif présente une symétrie de révolution, la diffusion est donc radiale ~th = jth (r) ~ur ;
➜ le terme de production est nul dans la gaine isolante (pas de réaction chimique, pas de réaction
nucléaire, pas de courant électrique) ;
➜ le régime est stationnaire ;
➜ les grandeur pertinentes sont λ (caractéristique thermique), R1 , R2 et h (géométrie) et T1 et
T2 comme contrainte mais si nous contentons de trouver la résistance thermique de la gaine,
nous savons que les températures n’interviendront pas.
✧ Analyse technique :
➜ étant donné la situation (pas de création, RPS), nous pouvons utiliser une approche en terme
de résistance thermique ;
➜ pour déterminer la résistance thermique nous avons deux méthodes ;
➙ nous pouvons procéder par association de résistances thermiques ;
➙ nous pouvons reprendre la définition et calculer le flux en fonction des températures.

✬ méthode : associer des résistances élémentaires

✧ Étant donné que le courant thermique est radial, nous pouvons voir l’ensemble de la gaine comme
l’association de cylindres creux d’épaisseur dr.
© Matthieu Rigaut 63 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·5 – Résistance thermique

dr r

h h

✧ Chaque morceau possède sa propre résistance thermique dRth puisqu’il est en régime permanent et
sans terme de source.
✧ De plus nous voyons que les différentes résistances élémentaires sont mises bout à bout : elles sont
en série.
✧ Dès lors nous pouvons écrire
Z
Rth = dRth

✧ La résistance élémentaire s’écrit

1 longueur traversée
dRth = ×
λ section
✧ Ici nous avons
➜ longueur traversée= dr ;
➜ section= 2 π r h.
✧ Ce qui donne

dr
dRth =
λ2πrh
✧ Maintenant, il suffit de calculer

Z 2
dr
Rth =
1 λ2πrh
Z R2
dr
=
R1 λ 2 π r h
Z R2
1 dr
= ×
2πλh R1 r

✧ La résistance thermique s’écrit donc

1 R2
Rth = × ln
2πλh R1

© Matthieu Rigaut 64 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) I·5 – Résistance thermique

✬ méthode : partir de la loi globale

✧ L’idée est identique à celle qu’il faut avoir quand nous cherchons l’expression d’une résistance équi-
valente alors que nous ne pouvons pas faire d’association

R 2R

3R
A B
2R R

✧ Dans ce cas :
➜ nous imposons V1 − V2 ;
➜ nous calculons i1→2 ;
V1 − V2
➜ nous en déduisons Réq = .
i1→2
✧ Nous pourrions faire la même chose en diffusion mais il est techniquement plus facile d’imposer
un flux et de résoudre en température que le contraire.
✧ Reprenons la situation.

✧ Imposons un flux Φ0 qui va de 1 vers 2.

✧ Parce que le flux est conservé (milieu sans source) nous avons

Φ0 = Φ(r)
✧ Nous trouvons ensuite très vite que
ZZ
Φ(r) = jth (rP ) dSP 2 π r h j (r) = Φ0
P

✧ Si nous ajoutons à cela la loi de Fourier cela donne

dT dT dT Φ0
j (r) = −λ (r) −2 π r h λ (r) = Φ0 (r) = −
dr dr dr 2πrhλ
✧ Nous pouvons alors résoudre
© Matthieu Rigaut 65 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) I·5 – Résistance thermique

Φ0 r
T (r) = − × ln
2πhλ r0
✧ Il ne reste plus qu’à exprimer T1 − T2 sans avoir besoin d’exprimer la constante d’intégration en ln r0
Φ0 R1
T1 − T2 = T (R1) − T (R2) T1 − T2 = − × ln
2πhλ R2
✧ Et nous voyons apparaître l’expression de la résistance thermique.
1 R2
T1 − T2 = Rth × Φ0 Rth = × ln
2πhλ R1
✧ Il s’agit bien là de la même expression que précédemment.

✬ comparaison
✧ La première méthode est certainement plus physique et (donc ?) plus facile mais elle nécessite de
bien voir le découpage et l’association de résistances élémentaires.
✧ La deuxième méthode est plus systématique mais est moins naturelle, surtout parce qu’il faut se
souvenir de l’astuce qui consiste à imposer un flux, et non à partir des températures, ce que tout le
monde fait spontanément.

© Matthieu Rigaut 66 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) II – Diffusion de particules

II – Diffusion de particules
II·1 – Molécules en mouvement
II·1·i – mouvement brownien
✧ Nous savons déjà que les molécules s’agitent dans tous les sens à cause de l’agitation thermique.
✧ Une molécule particulière ou même un objet de taille mésoscopique peut ainsi être influencé par les
chocs incessants de ces particules.
✧ Il en résulte un mouvement saccadé, assez aléatoire, qualifié de « marche au hasard » comme si
chaque « pas » n’avait aucun lien avec le précédent.
✧ Voici ci-dessous une image tirée de Wikipédia 7 illustrant cette marche aléatoire.

Le mouvement brownien est une marche aléatoire d’une particule ou d’un objet
suffisamment petit pour être sensible aux chocs avec les particules du milieu en
mouvement à cause de l’agitation thermique.

✧ Notons que le mouvement brownien existe aussi dans les solides même si, évidemment, c’est beaucoup
plus lent.
✧ C’est cette marche aléatoire qui va être à l’origine de la diffusion particulaire.

II·1·ii – exemples de phénomènes

✧ Dans le quotidien, il y a la tache d’encre sur un tissu ou un buvard : les molécules colorées, entraînées
par le liquide, diffusent à travers les fibres du tissu.
✧ Dans l’industrie, le phénomène de diffusion est utilisé :
➜ dans le dopage du silicium, dopage nécessaire à la fabrication de composants électroniques dont
l’utilité se passe de commentaire ;
➜ dans les réacteurs nucléaires.
✧ Au niveau physiologique l’influx nerveux est régi par des phénomènes de diffusion qui, à l’échelle
microscopique peuvent donc devenir « rapides ».
7. Source : [Link]

© Matthieu Rigaut 67 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) II·2 – Bilan moléculaire

II·2 – Bilan moléculaire

II·2·i – modélisation
✧ Considérons un milieu matériel globalement immobile dans lequel se trouvent des particules intéres-
santes.

✧ Les particules intéressantes sont celles dont nous suivrons la diffusion.

✬ densité particulaire
✧ Pour suivre les particules diffusantes, nous allons utiliser la densité particulaire.

La densité particulaire n⋆ est définie par dN = n⋆ dτ où dN est le nombre de particules

intéressantes dans le volume dτ .
n⋆ est en m−3 .

✧ D’autres fois nous utiliserons la concentration en mol.L−1 ou en mol.m−3 pour suivre le mouvement
des particules intéressantes.
✧ Parfois, mais pas toujours, nous rencontrerons la notation n pour la densité particulaire mais cela
risque d’occasion une collusion de notation entre
➜ la densité particulaire n en m−3 ;
➜ la quantité de matière n en mol ;
➜ la concentration de matière n en mol.m−3 .
✧ Dans la mesure du possible, nous utiliserons la notation n⋆ pour éviter autant que possible les
ambiguités.

✬ courant particulaire
✧ Tout comme il existe une densité surfacique de courant thermique en volume, nous allons introduire
un vecteur ~n caractérisant le mouvement des molécules diffusantes.

Le vecteur densité surfacique de courant de particules en volume ~n est défini par
~ dt
δN = ~n · dS où :
~
➜ dS est le vecteur surface de la surface à travers laquelle passent les particules ;
➜ δN est le nombre de particules qui passent pendant dt dans le sens de dS. ~

II·2·ii – bilan particulaire 1D

✬ situation
✧ Imaginons la situation suivante : dans un tube transparent scellé (pour éviter les phénomènes de
convection) un peu de gaz particulier est introduit à une extrémité.
© Matthieu Rigaut 68 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) II·2 – Bilan moléculaire

✧ Nous savons qu’à terme sa concentration sera homogène mais comment décrire l’évolution de la
concentration ?

✬ approche mésoscopique
✧ Étant donné que le rayon du tuyau est bien plus faible que sa longueur, nous pouvons supposer
l’homogénéité sur une section et donc décrire le nombre de particules uniquement en fonction de
l’abscisse.
✧ Comme le nombres de particules varie en espace et en temps, pour faire un bilan de particules nous
allons utiliser une approche mésoscopique.
✧ Commençons par zoomer.

δN (x) δN (x + δx)

~n (x,t) ~n (x + δx,t)

x x + δx
✧ Nous allons faire un bilan de particules sur S { la tranche entre x et x + δx } entre t et t + dt.
K Remarque. le système S choisi n’est pas fermé. Et pour cause : il y a diffusion de particules. Toutefois
comme nous n’allons pas écrire ni de lois mécaniques ni de lois thermodynamiques dessus, cela ne
sera pas grave.
✧ Le bilan de particules sur S peut se résumer sous la forme

VARIATION dans le temps = ÉCHANGE à travers la surface + CRÉATION en volume

✬ variation dans le temps

✧ Il s’agit là de la variation du nombre de particules intéressantes contenues dans S .
✧ En notant N (t) le nombre de particules contenues dans S à l’instant t nous avons tout de suite, par
définition de la densité particulaire

N (x,t) = n⋆ (x,t) dτ
✧ Ce qui donne

δNvar = N (x,t + dt) − N (x,t) δNvar = (n⋆ (x,t + dt) − n⋆ (x,t)) dτ

✧ Avec un développement limité à l’ordre 1 de n⋆ (x,t + dt) nous avons

⋆ ⋆ ∂n⋆
n (x,t + dt) = n (x,t) + dt × (x,t)
∂t
© Matthieu Rigaut 69 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) II·2 – Bilan moléculaire

✧ Ce qui donne finalement

∂n⋆
VARIATION dans le temps = (x,t) dτ dt
∂t

✬ échange à travers la surface

✧ Dans l’exemple choisi les particules ne peuvent évoluer (diffuser) que le long du tube donc les échanges
ne se font qu’à deux endroits.
✧ Le terme d’échange se réduit donc à la somme de deux termes

δNéch = δNqui rentre en x + δNqui rentre en x + δx

✧ En notant S la section du tube, la notion de densité surface de courant de particules en volume nous
permet d’arriver à (attention aux signes)

δNqui rentre en x = +jn (x,t) S dt et δNqui rentre en x + δx = −jn (x + δx,t) S dt

✧ En rassemblant nous arrivons à

δNéch = jn (x,t) S dt − jn (x + δx,t) S dt

✧ Avec un développement limité à l’ordre 1 de n(x + δx,t) nous avons
∂jn
jn (x + δx,t) = jn (x,t) + δx × (x,t)
∂x
✧ Ce qui conduit à
∂jn
ÉCHANGE à travers la surface = − (x,t) × S δx dt
∂x

✬ production en volume
✧ Comme nous le savons, il n’est pas possible de « créer » des particules, quoique cela puisse encore
être discutable à l’échelle nucléaire.
✧ En revanche il est très facile avec la chimie de transformer une particule inintéressante du substrat
en une particule intéressante qui diffuse.
✧ Comme nous ne comptons que les particules intéressantes, celles-ci sembleront apparaître alors qu’au
fond il n’en est rien.

Les phénomènes à l’origine de la « création » de particules diffusantes sont les réactions

chimiques et nucléaires.

✧ Nous pouvons alors définir un taux de créaction en volume.

Le taux de création volumique de particules σ est tel que

δNcréé = σ dτ dt

✧ Normalement l’unité de σ est m−3 .s−1 si le bilan est particulaire mais si le décompte se fait en mol,
il est possible d’avoir σ en mol.m−3 .s−1 .
✧ Dans ces conditions la partie création s’écrit immédiatement

CRÉATION en volume = σ (x,t) dτ dt

© Matthieu Rigaut 70 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) II·2 – Bilan moléculaire

✬ rassemblement
✧ En rassemblant les trois termes de variation, échange et création nous avons déjà

∂n⋆ ∂j
✚ = − n (x,t) × S δx ✚
(x,t) dτ ✚
dt ✚ + σ (x,t) dτ ✚
dt ✚
dt
∂t ∂x
✧ En remarquand que dτ = S δx nous avons

∂n⋆ ∂j
✚ = − n (x,t) ✚✚ + σ (x,t) ✚✚
dτ
(x,t) ✚ dτ dτ
∂t ∂x
✧ Et donc finalement

∂n⋆ ∂jn
(x,t) + (x,t) = σ (x,t)
∂t ∂x

La loi de conservation des particules, appelée aussi équation de continuité s’écrit, dans le
cas d’une diffusion unidimensionnelle
∂n⋆ ∂jn
(x,t) + (x,t) = σ (x,t)
∂t ∂x

✧ Nous reconnaissons là une équation de conservation avec :

➜ un couplage des variations spatio-temporelles ;
➜ une dérivée simple en temps et en espace.

II·2·iii – bilan particulaire 3D

✬ approche globale
✧ Nous allons cette fois établir l’équation de continuité mais dans un cas plus général.
✧ Cela nous amènera sans surprise au résultat suivant.

La loi de conservation des particules, appelée aussi équation de continuité s’écrit, dans le
cas général
∂n⋆
(x,t) + div ~
n (x,t) = σ (x,t)
∂t

✧ Reprenons la démonstration et commençons par considérer un volume V quelconque mais fixe dans
le temps.

~
dS

dτ

© Matthieu Rigaut 71 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) II·2 – Bilan moléculaire

✧ Faisons un bilan de particules pour ce volume V entre les instants t et t + dt.

VARIATION dans le temps = ÉCHANGE à travers la surface + CRÉATION en volume

✬ variation dans le temps

✧ À un instant t quelconque, la quantité totale N (t) de particules contenues dans le volume V s’écrit,
par extensivité
ZZZ ZZZ
N (t) = dNP N (t) = n⋆ (P,t) dτP
P ∈V P ∈V
✧ La variation du nombre de particules s’écrit donc, en utilisant un développement limité

dN dN
δNvar = N (t + dt) − N (t) et N (t + dt) = N (t) + dt × (t) δNvar = (t) dt
dt dt
✧ Remarquons qu’ici la dérivée temporelle est bien droite.
✧ Avec l’expression de N (t) cela donne
ZZZ !
d
n⋆ (P,t) dτP dt
dt P ∈V

✧ Et comme le domaine d’intégration V ne dépend pas du temps, il est possible de « rentrer » la dérivée
sous le signe somme ce qui donne
ZZZ
∂n⋆
δNvar = (P,t) dτP dt
P ∈V ∂t
✧ Et finalement
ZZZ
∂n⋆
VARIATION dans le temps = (P,t) dτP dt
P ∈V ∂t

✬ échange à travers la surface

✧ Pas de surprise ici : la totalité des particules qui rentrent dans le volume est la somme des particules
qui rentrent par chaque portion de surface
ZZ
δNéch = δ 2 Néch,P
P ∈S
✧ En faisant attention au sens conventionnel du vecteur surface, nous avons
ZZ
2 ~
δ Néch,P = −~n (P,t) · dS dt δNéch = − ~ dt
~n (P,t) · dS
P ∈S
✧ Et en faisant jouer Green – Ostrogradski
ZZZ
δNéch = − div ~n (P,t)dτ dt
P ∈V
✧ Ce qui donne finalement
ZZZ
ÉCHANGE à travers la surface = − div ~n (P,t)dτ dt
P ∈V

© Matthieu Rigaut 72 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) II·3 – Équation de diffusion particulaire

✬ production en volume
✧ En utilisant l’extensivité de la production, nous avons tout de suite
ZZZ ZZZ
2
δNcréé = δ Ncréé,P δNcréé = σ (P,t)dτ dt
P ∈V P ∈V

✧ Ce qui conduit à
ZZZ
CRÉATION en volume = σ (P,t)dτ dt
P ∈V

✬ rassemblement
✧ En rassemblant nous avons d’abord
ZZZ ZZZ ZZZ
∂n⋆
(P,t) dτP dt = − div ~n (P,t)dτ dt + σ (P,t)dτ dt
P ∈V ∂t P ∈V P ∈V

✧ Comme le volume V est le même pour les trois sommes, nous pouvons les regrouper en une seule
ZZZ !
∂n⋆
(P,t) + div ~
n (P,t) − σ (P,t) dτP dt = 0
P ∈V ∂t
✧ Et comme ce résultat est nul quel que soit le volume V c’est que l’intégrande est nul donc

∂n⋆ ∂n⋆
+ div ~n − σ = 0 + div ~n = σ
∂t ∂t
✧ Ce qui est bien le résultat attendu.

II·3 – Équation de diffusion particulaire

II·3·i – loi de Fick
✧ Tout comme pour la diffusion thermique, pour obtenir l’équation de diffusion particulaire il manque
la loi qui explique comment les particules « bougent », i.e. la loi qui explique comment fonctionne
~n .

✬ énoncé
✧ Pas de surprise, elle ressemble à la loi de Fourier.

Loi de Fick
Dans un milieu diffusant, la densité surfacique de courant de particules en volume s’écrit
−−→
~n = −D grad n⋆ où :
➜ n est la densité particulaire en m ;
⋆ −3

➜ D est le coefficient de diffusivité en m2 .s−1 .

✧ Pour l’unité de D nous avons en effet

~ dt 1
δN = ~n · dS [jn ] =
m2 .s
✧ Et ainsi
© Matthieu Rigaut 73 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) II·3 – Équation de diffusion particulaire

−−→ 1 1 m2
[D grad n⋆ ] = [D] × 4 [D] = × m4 [D] =
m m2 .s s

✬ interprétation, limites
✧ L’interprétation est globalement la même que pour la diffusion thermique :
➜ le courant de particule ~n est dirigé vers les zones de densité basse ;
➜ la diffusion tend à uniformiser n⋆ .
✧ Les limites sont elles aussi analogues à celles de la loi de de Fourier :
➜ il ne faut pas un gradient de densité trop important ou d’évolution temporelle trop rapide ;
➜ le milieu qui sert de support à la diffusion doit être isotrope.

✬ quelques valeurs
✧ Dans un gaz au repos (donc sans convection) le coefficient D va de 10−6 à 10−4 m2 .s−1 .
✧ Dans un liquide, là aussi au repos, le coefficient D va de 10−12 à 10−8 m2 .s−1 .
✧ Dans un solide, forcément au repos, le coefficient D va de 10−30 à 10−8 m2 .s−16 .
✧ Nous voyons donc que, suivant le substrat, il peut y avoir une variation de près de 26 ordres de
grandeur en terme de diffusion entre celle qui se fait dans un gaz et celle qui a lieu dans un solide.

II·3·ii – équation de diffusion à 1D

✬ partir de l’équation de continuité
✧ Rappelons l’équation de continuité en 1D

∂n⋆ ∂jn
+ =σ
∂t ∂x
✧ La loi de Fick donne, en 1D

−−→ ∂n⋆
~n = −D grad n⋆ jn = −D
∂x
✧ Et donc en injectant l’un dans l’autre

∂n⋆ ∂ 2 n⋆
−D =σ
∂t ∂x2

L’équation de diffusion unidimensionnelle vérifiée par la densité particulaire n⋆ s’écrit,

en notant D le coefficient de diffusivité
∂n⋆ ∂ 2 n⋆
−D =σ
∂t ∂x2

✧ Avec le cas particulier usuel

L’équation de diffusion unidimensionnelle vérifiée par la densité particulaire n⋆ s’écrit,

en notant D le coefficient de diffusivité et lorsqu’il n’y a pas de terme de source
∂ 2 n⋆ 1 ∂n⋆
=
∂x2 D ∂t

© Matthieu Rigaut 74 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) II·3 – Équation de diffusion particulaire

✬ approche qualitative
✧ Nous voyons tout de suite pourquoi dans le cas de la diffusion thermique nous distinguons la conduc-
tivité λ de la diffusivité a alors que ce n’est pas le cas pour la diffusion particulaire.
✧ En effet dans le cas de la diffusion particulaire il n’y a pas « d’inertie » particulaire alors qu’il existe
une inertie thermique caractérisée par ρ c.
✧ C’est là la « grande » différence entre les deux types de diffusion : dans le cas de la diffusion
particulaire, D sera toujours une grandeur pertinente.
✧ Pour le reste, nous pouvons interpréter cette équation exactement de la même manière que nous
l’avons fait pour la diffusion thermique :
➜ il s’agit d’une loi linéaire ;
➜ il s’agit d’une loi associée à un processus irréversible (à cause de la dérivée première en temps) ;
➜ il s’agit d’une loi associée à un processus lent.

Dans le cas de la diffusion particulaire, la longueur caractéristique de diffusion L est

reliée à la durée caractéristique τ par l’intermédiaire du coefficient de diffusivité D par
√
L = Dτ

✧ Cela vient directement de l’équation de diffusion particulaire sans source

∂ 2 n⋆ 1 ∂n⋆
=
∂x2 D ∂t
✧ En ordre de grandeur nous avons

∂ 2 n⋆ n⋆ 1 ∂n⋆ n⋆
≡ et ≡
∂x2 L2 D ∂t Dτ
✧ Ce qui donne bien

n⋆ n⋆
= L2 = D τ
L2 Dτ

✬ interprétation en terme de mouvement

Y imprévisible au niveau moléculaire – marche au hasard
✧ Comme nous l’avons expliqué au début de cette partie, la diffusion particulaire est liée au mouvement
brownien.
✧ En fait, il est tout à fait possible de retrouver l’équation de diffusion en modélisant le comportement
d’une particule par une marche au hasard 8.
✧ Le fait même que cette marche aléatoire modélise un phénomène de diffusion explique pourquoi ce
dernier est irréversible.
✧ En effet le caractère « aléatoire » de la marche d’une molécule implique qu’il n’est pas possible de
prévoir le pas « juste suivant » de la molécule : elle devient imprévisible.
✧ Sauf que si le mouvement devient imprévisible « vers le futur » c’est qu’il est imprévisible « vers le
passé » et donc qu’il est impossible d’aller spécifiquement et naturellement vers le passé.
✧ Quand un phénomène ne peut pas naturellement revenir sur ses pas, c’est qu’il est intrinsèquement
irréversible !
8. Le modèle le plus simple consiste à prendre un axe infini découpé en cases de même taille a. Une particule saute
à tous les instants n t0 (n entier) d’une case avec une probalité 1/2 pour la droite et 1/2 pour la gauche. Après un
développement limité en a et t0 considérés comme petits, le lecteur pourra montrer que la probabilité de présence
obéit à l’équation de diffusion.

© Matthieu Rigaut 75 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) II·3 – Équation de diffusion particulaire

✧ Le détail « amusant » c’est qu’il est irréversible parce qu’aléatoire (et donc imprévisible) au niveau
microscopique mais il est parfaitement prédictible au niveau méso et macroscopique !
✧ Il s’agit là d’un effet des grands nombres : plus il y a de phénomènes aléatoires, plus leurs compor-
tement d’ensemble est prévisible.

Y prévisible au niveau mésoscopique – vitesse d’ensemble

✧ S’il n’est pas possible de regarder précisément comment chaque particule bouge lors de la diffu-
sion, nous pouvons néanmoins associer à l’ensemble de ces particules une vitesse globale (ou vitesse
d’ensemble) de diffusion.

Lors d’un phénomène de diffusion, les particules diffusantes ont une vitesse moyenne ~vdiff
telle que
~n = n⋆ ~vdiff où :
➜ n⋆ est la densité particulaire ;
➜ ~n est la densité surfacique de courant de particule en volume.

✧ Pour montrer ce résultat, nous allons montrer que le nombre de particules qui traverse une surface
dS durant dt s’écrit

~ dt
δNqui traverse = n⋆ ~vdiff · dS

~
✧ Pour cela considérons une surface dS caractérisé par le vecteur surface dS

~
dS

✧ Comptons le nombre δNqui traverse de particules qui passe à travers dS durant dt.
✧ Pour simplifier considérons que toutes les particules ont la même vitesse ~vdiff .
✧ Dans ces conditions toutes les particules, où qu’elles soient, avanceront de vdiff dt.
© Matthieu Rigaut 76 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) II·3 – Équation de diffusion particulaire

θ ~
dS

dt
v diff

h
✧ Seules celles (en rouge) qui sont dans le cylindre repéré ci-dessus pourront passer. Les autres (en
bleu), ne passeront pas à travers la surface ou ne n’y passeront pas durant dt
✧ Or le volume d’un cylindre ne dépend pas de son inclinaison mais seulement de sa hauteur comme
l’illustre un simple paquet de cartes

✧ Le volume du cylindre dans lequel doivent se trouver les particules vaut donc

dτ = dS × h et h = vdiff dt cos θ dτ = dS × vdiff dt cos θ

✧ Et dans ce volume, le nombre de particules diffusant s’écrit

δNdans dτ = n⋆ dτ δNdans dτ = n⋆ dS vdiff dt cos θ

✧ Donc finalement, le nombre de particules diffusantes qui traversent dS durant dt s’écrit

δNqui traverse = δNdans dτ δNqui traverse = n⋆ dS vdiff dt cos θ

✧ Soit, en notation vectorielle

~ dt
δNqui traverse = n⋆ ~vdiff · dS
✧ Et par identification avec la définition de ~n nous avons

~ dt
δNqui traverse = ~n · dS ~n = n⋆ ~vdiff

II·3·iii – équation de diffusion directement en symétrie cylindrique

✧ Nous allons voir dans cette partie comment trouver l’équation de diffusion dans le cas particulier de
la symétrie cylindrique (diffusion radiale) sans passer par les opérateurs vectoriels.
© Matthieu Rigaut 77 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) II·3 – Équation de diffusion particulaire

✬ situation et analyse

✧ Imaginons une situation de diffusion radiale en symétrie cylindrique.

✧ Des particules ont été introduites dans un barreau de manière à diffuser lentement vers l’extérieur.
✧ Ce genre de procédé peut être utilisé pour la délivrance de produits médicamenteux sur de longues
durées.

✧ Le dispositif est suffisamment grand pour pouvoir négliger les effets de bord.
✧ Analyse physique :
➜ il s’agit d’un phénomène de diffusion ;
➜ la diffusion est radiale donc le courant de particule est porté uniquement par ~ur ;
➜ l’invariance par rotation et l’invariance par translation (ou plutôt la non prise en compte des
effets de bord) permettent d’écrire que les grandeurs ne dépendent pas ni de z ni de θ.
✧ Analyse technique :
➜ ne sachant pas si le terme de création est uniforme ou non, nous allons procéder à une approche
mésoscopique ;
➜ nous allons choisir un volume de contrôle respectant au mieux les symétries du problème.

✬ un système particulier pour un bilan usuel

Y schéma

✧ Choisissons comme système S le volume V { cylindre creux de rayon r, d’épaisseur dr et de hauteur

h }.
© Matthieu Rigaut 78 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) II·3 – Équation de diffusion particulaire

dr
dr

r
h h

vue de dessus

vue en perspective vue en coupe

✧ Sur le schéma ci-dessus nous avons représenté qualitativement le courant de particules.

✧ Nous pouvons faire un bilan de particules pour S entre t et t + dt.

VARIATION dans le temps = ÉCHANGE à travers la surface + CRÉATION en volume

Y variation dans le temps

✧ Comme le volume V est suffisamment bien choisi pour être tel que la densité particulaire y soit
uniforme, nous pouvons dire qu’à l’instant t le nombre de particules qui y sont contenues s’écrit

N (t) = n⋆ (r,t) × V N (t) = n⋆ (r,t) × 2 π r dr h

✧ De même à t + dt

N (t + dt) = n⋆ (r,t + dt) × 2 π r dr h

✧ La variation s’écrit donc


δNvar = N (t + dt) − N (t) δNvar = n⋆ (r,t + dt) − n⋆ (r,t) 2 π r dr h
✧ Et avec un développement limité à l’ordre 1

⋆ ⋆ ∂n⋆ ∂n⋆
n (r,t + dt) = n (r,t) + (r,t) × dt δNvar = (r,t) × dt 2 π r dr h
∂t ∂t
✧ Finalement

∂n⋆
VARIATION dans le temps = (r,t) × dt 2 π r dr h
∂t

Y production en volume
✧ Encore une fois, commme le volume V est uniforme puisqu’à r = Cte , nous avons tout de suite

δNcréé = σ (r,t) × V dt δNcréé = σ (r,t) × 2 π r dr h dt

© Matthieu Rigaut 79 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) II·3 – Équation de diffusion particulaire

✧ Insistons : c’est à l’intérieur de V que la production est uniforme. En prenant en compte la totalité
du barreau le terme de production peut fort bien a priori dépendre de r et t.
✧ Ce qui donne

CRÉATION en volume = σ (r,t) × 2 π r dr h dt

Y échange à travers la surface

✧ Nous avons affaire ici à 4 surfaces :
➜ une couronne au-dessus et en dessous ;
➜ la surface en r ;
➜ la surface en r + dr.
✧ Étant donné que la diffusion est radiale, il n’y a pas de particules qui passent à travers les couronnes
du dessus et du dessous donc il reste

δNreçu = δNreçue en r + δNreçue en r + dr + 0 + 0

✧ Comme il y a uniformité sur chacune des deux surfaces du courant de particules, nous avons

δNreçue en r = +jn (r,t) × S (r) et Nreçue en r + dr = −j (r + dr,t) × S (r + dr)

! Attention ! Ici et contrairement au cas unidimensionnel l’aire de la surface d’échange dépend de la

position !
✧ Cela n’empêche pas de faire un développement limité mais pas directement en jn (r)

∂f
δNreçu = f (r,t) dt − f (r + dr,t) dt avec f (r) = jn (r) S (r) δNreçu = − (r,t) dr dt
∂r
✧ Et ainsi

∂(jn S)
δNreçu = − (r,t) dr dt
∂r
✧ Et comme la surface S (r) s’écrit S (r) = 2 π r h nous obtenons

∂(jn 2 π r h) ∂(jn r)
δNreçu = − (r,t) dr dt δNreçu = − (r,t) 2 πdr dt h
∂r ∂r
✧ Et ainsi

∂(jn r)
= ÉCHANGE à travers la surface = − (r,t) 2 πdr dt h
∂r

Y rassemblement
✧ Le bilan de particules s’écrit donc

∂n⋆ ∂(jn r) ✭ ✘
(r,t) × ✘ 2✘
dt✘ π r✘ ✘
dr✘h=− (r,t) ✭ ✭✭✭
2 πdr 2✟
dt h + σ (r,t) × ✟ dr✘h✘
π r✘ dt
∂t ∂r
✧ Ce qui donne

∂n⋆ 1 ∂(jn r)
(r,t) + (r,t) = σ (r,t)
∂t r ∂r
✧ En faisant intervenir la loi de Fick nous avons ainsi
© Matthieu Rigaut 80 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) II·3 – Équation de diffusion particulaire

!
∂n⋆ ∂n⋆ 1 ∂ ∂n⋆
jn (r,t) = −D (r,t) (r,t) − D × r× (r,t) = σ (r,t)
∂r ∂t r ∂r ∂r

II·3·iv – équation de diffusion directement en symétrie sphérique

✬ situation et analyse
✧ Nous allons considérer une situation similaire à la précédente mais différente dans la géométrie.
✧ Au lieu d’avoir un barreau diffusant, nous allons prendre une sphère diffusante avec une diffusion
radiale.

✧ L’analyse physique est analogue :

➜ il s’agit d’un phénomène de diffusion ;
➜ la diffusion est radiale donc le courant de particule est porté uniquement par ~ur (mais ce n’est
pas le même ~ur qu’en cylindrique) ;
➜ les deux invariances par rotation permettent d’écrire que les grandeurs ne dépendent pas ni de
θ ni de ϕ.
✧ Analyse technique :
➜ ne sachant pas si le terme de création est uniforme ou non, nous allons procéder à une approche
mésoscopique ;
➜ nous allons choisir un volume de contrôle respectant au mieux les symétries du problème.

✬ un système particulier pour un bilan usuel

Y schéma
✧ Choisissons comme système S le volume V { coquille creuse de rayon r et d’épaisseur dr }.

dr

vue en perspective vue en coupe

© Matthieu Rigaut 81 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) II·3 – Équation de diffusion particulaire

✧ Sur le schéma ci-dessus nous avons représenter qualitativement le courant de particules.

✧ Nous pouvons faire un bilan de particules pour S entre t et t + dt comme nous le faisons habituel-
lement.

VARIATION dans le temps = ÉCHANGE à travers la surface + CRÉATION en volume

Y variation dans le temps

✧ Comme le volume V est suffisamment bien choisi pour être tel que la densité particulaire y soit
uniforme, nous pouvons dire qu’à l’instant t le nombre de particules qui y sont contenues s’écrit

N (t) = n⋆ (r,t) × V N (t) = n⋆ (r,t) × 4 π r 2 dr

✧ De même à t + dt

N (t + dt) = n⋆ (r,t + dt) × V4 π r 2 dr

✧ La variation s’écrit donc


δNvar = N (t + dt) − N (t) δNvar = n⋆ (r,t + dt) − n⋆ (r,t) 4 π r 2 dr
✧ Et avec un développement limité à l’ordre 1

⋆ ⋆ ∂n⋆ ∂n⋆ 2
n (r,t + dt) = n (r,t) + (r,t) × dt δNvar = (r,t) × dt 4 π r dr
∂t ∂t
✧ Finalement
∂n⋆ 2
VARIATION dans le temps = (r,t) × dt 4 π r dr
∂t

Y production en volume
✧ Encore une fois, commme le volume V est uniforme puisqu’à r = Cte , nous avons tout de suite

δNcréé = σ (r,t) × V dt δNcréé = σ (r,t) × 4 π r 2 dr dt

✧ Ce qui donne

CRÉATION en volume = σ (r,t) × 4 π r 2 dr dt

Y échange à travers la surface

✧ Ici, la situation est plus simple que pour le barreau puisque le volume considéré n’est constitué que
des deux surfaces en r et en r + dr.
✧ Comme il y a uniformité sur chacune des deux surfaces du courant de particules, nous avons

δNreçue en r = +jn (r,t) × S (r) et Nreçue en r + dr = −j (r + dr,t) × S (r + dr)

! Attention ! Ici aussi, comme pour la diffusion radiale cylindrique, l’aire de la surface d’échange
dépend de la position !
✧ Cela n’empêche toujours pas de faire un développement limité mais pas directement en jn (r)

∂f
δNreçu = f (r,t) dt − f (r + dr,t) dt avec f (r) = jn (r) S (r) δNreçu = − (r,t) dr dt
∂r
© Matthieu Rigaut 82 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) II·3 – Équation de diffusion particulaire

✧ Et ainsi

∂(jn S)
δNreçu = − (r,t) dr dt
∂r
✧ Et comme la surface S (r) s’écrit S (r) = 4 π r 2 nous obtenons

∂(jn 2 π r 2 ) ∂(jn r 2 )
δNreçu = − (r,t) dr dt δNreçu = − (r,t) 4 πdr dt
∂r ∂r
✧ Et ainsi

∂(jn r 2 )
= ÉCHANGE à travers la surface = − (r,t) 4 πdr dt
∂r

Y rassemblement
✧ Le bilan de particules s’écrit donc

∂n⋆ ∂(jn r 2 ) ✘
(r,t) × ✭ 4✭
dt✭
✭
π✭
dr r 2 = − (r,t) ✘ ✘✘✘
4 πdr ✭✭
4 π✭dr
dt + σ (r,t) × ✭
✭
dt r 2
∂t ∂r
✧ Ce qui donne

∂n⋆ 1 ∂(jn r 2 )
(r,t) + 2 (r,t) = σ (r,t)
∂t r ∂r
✧ En faisant intervenir la loi de Fick nous avons ainsi
!
∂n⋆ ∂n⋆ 1 ∂ ∂n⋆
jn (r,t) = −D (r,t) (r,t) − D × r2 × (r,t) = σ (r,t)
∂r ∂t r 2 ∂r ∂r

II·3·v – équation de diffusion en 3D

✧ Aucune surprise.

L’équation de diffusion particulaire s’écrit

∂n⋆
− D △ n⋆ = σ où :
∂t
➜ D est le coefficien de diffusivité ;
➜ n⋆ est la densité particulaire ;
➜ σ est le taux de production volumique.

✧ Pour la redémontrer, il suffit de partir de l’équation de continuité en 3D

∂n⋆
+ div ~n = σ
∂t
✧ Ensuite, introduisons la loi de Fick et utilisons le fait que le matériaux soit homogène

 
−−→ ∂n⋆ −−→ ∂n⋆ −−→
~n = −D grad n⋆ + div −D grad n⋆ =σ + −D div grad n⋆ = σ
∂t ∂t
✧ Et il ne reste plus qu’à reconnaitre la définition du laplacien
−−→ not ∂n⋆
div grad () = △ () − D △ n⋆ = σ
∂t
© Matthieu Rigaut 83 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) II·4 – Quelques solutions

✧ Dans les exemples précédents, nous avons pu mettre en équation la diffusion mais sans utiliser les
opérateurs vectoriels divergence ou rotationnel.
✧ Cela ne signifie pas qu’il existe d’autres équations de diffusion mais bien que nous avons trouvé les
expressions de ces opérateurs dans des géométries particulières.

✬ laplacien en coordonnées cylindriques dans un cas particulier

✧ Dans le cas de la géométrie cylindre avec diffusion radiale, nous avons montré que
!
1 ∂(jn r) ⋆ 1 ∂ ∂n⋆
div ~n = (r,t) et △ n (r,t) = r× (r,t)
r ∂r r ∂r ∂r
✧ Il suffit pour cela de faire l’identification avec
∂n⋆ ∂n⋆ ⋆
(r,t) + div ~
n = σ (r,t) et (r,t) − D △ n (r,t) = σ (r,t)
∂t ∂t

✬ laplacien en coordonnées cylindriques dans un cas particulier

✧ De même dans le cas de la géométrie sphérique avec diffusion radiale, nous avons montré que
!
1 ∂(jn r 2 ) 1 ∂ ∂n⋆
div ~n = 2 (r,t) et △ n⋆ (r,t) = 2 r2 × (r,t)
r ∂r r ∂r ∂r
✧ Il suffit pour cela de faire l’identification avec
∂n⋆ ∂n⋆ ⋆
(r,t) + div ~
n = σ (r,t) et (r,t) − D △ n (r,t) = σ (r,t)
∂t ∂t

II·4 – Quelques solutions

II·4·i – tache d’encre – régime transitoire en 1D
✬ situation
✧ Modélisons une tâche d’encre par une situation 1D.
✧ Imaginons qu’à l’instant t = 0 toutes les particules soient en x = 0 et qu’après elles sont « libres »
de diffuser.
x=0
~ux

✧ Analyse physique :
➜ il s’agit là d’un phénomène de diffusion sans source ;
➜ le régime est non stationnaire, non permanent ;
➜ étant donné que la diffusion se fait en 1D, n⋆ sera une densité particulaire linéique en m−1 ;
➜ les grandeurs pertinente sont N0 le nombre total de particule diffusantes et D la diffusivité du
matériau.
✧ Les contraintes auxquelles devra obéir la solution est :
➜ toutes les particules sont localisées en x = 0 à t =Z 0, i.e. n⋆ (x 6= 0,0) = 0 ;
+∞
➜ le nombre total de particules ne change pas, i.e. n⋆ (x,t) dx = N (t) = N0 = Cte .
−∞

© Matthieu Rigaut 84 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) II·4 – Quelques solutions

✬ solution, représentation
✧ Après des calculs sans grand intérêt physique, nous obtenons
!
N0 x2
n⋆ (x,t) = √ exp −
2 πDt 4Dt
✧ Nous constatons qu’il s’agit là d’une gaussienne (fonction que nous avons déjà rencontrée en optique) :
➜ dont la valeur maximale√ est en x = 0 et décroît avec le temps ;
➜ dont l’écart-type en 2 D t augmente avec le temps.
✧ Voici la représentation de la solution pour différents instants régulièrement espacés.
n⋆

✧ Comme nous le savions, mais c’est toujours

√ mieux de le vérifier, la distance caractérique qui est
ici l’écart-type, est proportionnelle à D t illustrant une fois de plus la lenteur du phénomène de
diffusion.

II·4·ii – canal à ion – régime stationnaire 1D

✧ L’exemple ci-dessus est typique des phénomènes de diffusion au niveau moléculaire.
✧ Nous pouvons modéliser la paroi d’une cellule par un milieu imperméable à certains ions sauf au
niveau de « canaux » dans lequels ils diffusent.

milieu 1 milieu 2

n⋆1 = Cte n⋆2 = Cte

canal

✧ Les milieux à l’intérieur et à l’extérieur de la cellule sont si grands que nous pouvons les considérer
comme des réservoirs de concentration constantes.
✧ Dans ces conditions, si nous supposons le régime permanent stationnaire, comme la diffusion est en
1D, nous pouvons dire que la concentration à l’intérieur du canal est une fonction affine de l’espace
ce qui donne, en respectant la continuité de la concentration en x = 0 et x = e
n⋆2 − n⋆1
n⋆ (x) = n⋆1 + ×x
e
✧ Le flux de particules qui va du milieu ① vers le milieu ② s’écrit
−−→ DS ⋆
Φ1→2 = ~n · (S ~ux ) avec ~n = −D grad n⋆ Φ1→2 = (n1 − n⋆2 )
e
✧ Vérifions quand même la cohérence :
© Matthieu Rigaut 85 / 88 Version du 28 déc. 2013
 PC⋆ , Fabert (Metz) II·4 – Quelques solutions

➜ si n⋆1 > n⋆2 alors les particules vont de ① vers ② ;

➜ plus la diffusivité est grand plus Φ1→2 est grand ;
➜ plus la section S est grande plus Φ1→2 est grand ;
➜ plus la longueur e du canal e est petite plus Φ1→2 est grand.
✧ Pour finir, interprétons microscopiquement ce dernier résultat « plus la longueur e du canal e est
petite plus Φ1→2 est grand ».
✧ N’oublions pas que les particules qui diffusent marchent au hasard.
✧ Ainsi, plus le canal est long, plus il faut qu’une particule fasse de pas dans le même sens pour arriver
à passer de l’autre côté.
✧ C’est comme, sur un série de « pile ou face », faire plus de « pile » que de « face » : moins il y a de
lancers, plus statistiquement l’écart peut être grand.
✧ C’est ainsi que plus la distance à traverser s’allonge moins il est probable que les molécules y diffusent.

© Matthieu Rigaut 86 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) II·4 – Quelques solutions

Phénomènes de diffusion

Au niveau du cours
✬ Programme concerné
✧ Programme de 2e année :
➜ I.E Diffusion de particules – diffusion thermique.

✬ Les définitions
✧ Sont à savoir :
➜ transferts thermiques, par conduction, convection, rayonnement, conducto-convectif ;
➜ coefficient de conduction thermique, de diffusivité thermique ;
➜ flux surfacique, flux thermique, densité volumique de courant thermique ;
➜ création énergétique, courant thermique ;
➜ divergence, laplacien, nabla ;
➜ loi de Fourier, de Fick ;
➜ choc thermique, résistance thermique ;
➜ mouvement brownien, densité particulaire, flux particulaire, création particulaire volumique.

✬ Les grandeurs
✧ Connaître les petites relations suivantes ainsi que leur interprétation :
➜ δ 2 Q = ϕr dt : ϕr en W.m−2 ;
➜ ϕr = ~[Link] ~ : ~ en W ;
−−→
➜ ~th = −λ grad T : λ en W.m.K−1 ;
➜ n en m−3 ;
➜ ~j = n ~v : ~ en m−2 .s−1 ;
−−→
➜ ~ = −D grad n : D en m2 .s−1 .

✬ Les lois
✧ Sont à connaître :
➜ loi de conservation de l’énergie, des particules ;
➜ loi de Fourier, loi de Fick ;
➜ équation de la diffusion thermique, particulaire ;
➜ flux, théorème de Green – Ostrogradski.

✬ la phénoménologie
✧ Savoir :
➜ reconnaître la présence d’un phénomène de diffusion ;
➜ interpréter microscopiquement un phénomène de diffusion particulaire ;
➜ interpréter les lois de Fourier et de Fick ;
➜ le lien entre distance caractéristique et temps caractéristique dans un phénomène de diffusion.

© Matthieu Rigaut 87 / 88 Version du 28 déc. 2013

 PC⋆ , Fabert (Metz) II·4 – Quelques solutions

Au niveau des savoir-faire

✬ petits gestes
✧ Savoir :
➜ traduire les conditions aux limites.

✬ exercices classiques
✧ Savoir :
➜ retrouver l’équation de conservation (énergétique ou particulaire) à une variable (1D linéaire,
à symétrie cylindrique, à symétrie sphérique) ;
➜ retrouver l’équation de diffusion (thermique ou particulaire) à une variable (1D linéaire, à
symétrie cylindrique, à symétrie sphérique) ;
➜ retrouver l’équation de conservation et l’équation de diffusion (thermique ou particulaire) à 3D
(vectoriellement) ;
➜ retrouver les solutions dans le cas du régime permanent à 1D linéaire et avec des symétries
cylindriques ou sphériques.

© Matthieu Rigaut 88 / 88 Version du 28 déc. 2013