CORRECTION de la composition : autour de la diffusion

I. Préliminaires

1) Le moteur microscopique commun à tous les phénomènes de diffusion est l’agitation

thermique : à une température fixée, chaque molécule possède une énergie cinétique et peut
subir des chocs avec les molécules voisines. Plus l’agitation thermique est importante, plus les
chocs sont nombreux et les transports diffusifs efficaces (diffusion de particules, de la chaleur,
de la quantité de mouvement).

2) D, le coefficient de diffusion, est homogène à une longueur2/temps.

3) Cette équation, étant non invariante par renversement du temps, est irréversible. Par ailleurs,
cette équation est linéaire (une combinaison linéaire de solutions de cette équation est elle-
même une solution de cette dernière).

Du fait de la première propriété, les systèmes subissant des transformations adiabatiques et des
phénomènes de diffusion ont une entropie qui augmente.

Du fait de la deuxième propriété, on peut utiliser le théorème de superposition pour résoudre

l’équation de diffusion.

A A
4) Il découle directement de l’équation de diffusion : D , où A désigne l’ordre de
 L*
2

grandeur caractéristique du champ scalaire a(M,t). On en déduit : L*  D , relation

caractéristique des phénomènes diffusifs. Ces derniers sont peu efficaces pour des temps longs.

5) On parlera de régime stationnaire lorsque le champ scalaire a(M,t) est indépendant du temps.
Dans ce régime, l’équation de diffusion prend la forme suivante : a  0 . On reconnait
l’équation de Laplace qui admet une solution unique pour des conditions aux limites fixées.

6) En notation complexe, l’équation de diffusion s’écrit j a   D k 2 a ou encore k   j / D ,

2

ce qui correspond à la relation de dispersion cherchée.


j  1 j  1 j 2D
On en déduit : k  e 4
  avec  . Ainsi
D 2 D  
1 j 1 j z
j (t   z) j (t  z) j (t  )
a( M , t )  Ae 
 Ae 
 Ae 
e z /  pour une onde se propageant dans le sens
des z croissants.

1
La partie réelle de k dépendant de la pulsation  de manière non linéaire, la propagation est
dispersive.

Par ailleurs, le vecteur d’onde est complexe et possède une partie imaginaire non nulle négative
lorsque l’onde se propage dans le sens des (Oz) croissant ; la propagation de l’onde
s’accompagne donc d’une atténuation.

L’épaisseur de peau  est la distance caractéristique sur laquelle s’effectue la propagation et

2D
elle dépend de la fréquence :   .


II Etude de la diffusion de chaleur

II.1 La loi de Fourier

7) Une loi phénoménologique est une loi qui est déduite d’observations expérimentales et qui
n’a pas de fondement théorique.

Conditions de validité de la loi de Fourier : gradient thermique pas trop élevés, pas trop rapides
et matériaux isotropes.

8) J. Fourier est un mathématicien et physicien français du début du 19ième siècle (1768-1830).

Il est à l’origine de la décomposition des fonctions périodiques en séries trigonométriques et a
conduit des études sur la propagation de la chaleur.

9) Le coefficient est appelé conductivité thermique. Son unité est : W.K-1m-1 .

Ce coefficient dépend naturellement de la nature du matériau considéré (bois, métal..etc..). Il
peut également dépendre de la température, de la pression, de l’humidité…

10) Ordres de grandeurs :

métal eau air

W.K m )
-1 -1 1× 10 2 1 1× 10 -2

11) On observe que le glaçon sur la poêle fond plus vite que le glaçon sur la plaque de bois.
Ces deux glaçons n’ont donc pas reçu la même quantité de chaleur lors de cette transformation
à pression extérieure constante. Pendant un même intervalle de temps, le glaçon sur la poêle a
reçu une quantité de chaleur plus grande que celui sur le bois. Etant donné que les deux glaçons
sont entourés d’un air à la même température ambiante, on peut supposer, que la poêle est
meilleure conductrice de chaleur que la plaque en bois si on émet l’hypothèse que c’est que la
diffusion thermique qui régit les échanges thermiques. Si la poêle et la plaque ont la même
épaisseur, cette expérience montre que la conductivité thermique du bois est inférieure à la
conductivité thermique du métal.

2
12) On considère un système immobile, constitué du matériau à étudier, délimité par une surface
fixe  et indéformable et de volume V constant. On suppose donc que le système subit une
évolution isochore et le travail des forces de pression, qui correspond au seul travail des forces
extérieures non conservatives, est nul. Cette hypothèse est raisonnable si on considère des
phases condensées (solide ou liquide) ou des systèmes gazeux confinés dans un volume fixé
(air dans un double vitrage). On suppose également que le système a une énergie mécanique
macroscopique constante (énergie cinétique macroscopique nulle, énergie potentielle des forces
extérieures constante).

Appliquons le premier principe de la thermodynamique à ce système pendant une durée dt :

𝑑𝑈 + 𝑑𝐸𝑚𝑚𝑎𝑐𝑟𝑜 = 𝛿𝑄 + 𝛿𝑊𝑒𝑥𝑡 𝑛𝑜𝑛 𝑐𝑜𝑛𝑠𝑒𝑟𝑣𝑎𝑡𝑖𝑣𝑒𝑠 = 𝛿𝑄 + 𝛿𝑊𝑓𝑜𝑟𝑐𝑒𝑠 𝑒𝑥𝑡 𝑝𝑟𝑒𝑠𝑠𝑖𝑜𝑛

T
D’où 𝑑𝑈 = 𝛿𝑄 et ainsi dU   cv dt    jth .d 2 Sdt    div jth d 3 dt .
V
t  V

T
L’égalité devant être vraie quelle que soit la surface  choisie, il vient : cv  div jth  T ,
t
T  
ou encore  T  Dth T avec Dth  .
t cv cv

II.2 Mesure de conductivité thermique d’un solide : méthode de la plaque chaude gardée

13) On peut supposer que le modèle unidimensionnel est valable si la distance caractéristique
e est très inférieure devant la dimension transversale  . Ici  est de l’ordre de 600 mm et
l’épaisseur e varie entre 20 mm et 100 mm donc le rapport e /  varie entre 0,03 et 0,2. Cette
hypothèse semble donc raisonnable.

Par ailleurs cette hypothèse semble justifiée dans la mesure où le dispositif expérimental
comporte un anneau de garde.

On peut imaginer que l’anneau de garde fonctionne de la manière suivante : cet anneau doit
délivrer une densité de flux thermique supérieure à la densité de flux thermique  S délivrée par
la plaque chauffante pour compenser les perte latérales convectives. En pratique, l’anneau de
garde et la plaque chauffante sont séparés par une mince couche d’air.

 2T
14) En régime stationnaire T   0 . On en déduit que la fonction T (x) est une fonction
x 2
affine de x de type T ( x)  a x  b pour l’échantillon du bas et T ( x)  a x  b pour celui du
haut , où a , b , a et b sont des constantes.

3
15) La résistance thermique d’un échantillon de surface s et d’épaisseur e en régime stationnaire
T
est défini par l’expression suivante : Rth  , où  th représente le flux total de chaleur
 th
traversant l’échantillon et T  Te  Ts , l’écart de température auquel il est soumis, dans la
convention récepteur de l’électricité.

T
D’après la question précédente, est une grandeur constante dans un échantillon d’épaisseur
x
e et de surface s, le flux de chaleur dans les échantillons est donc constant et est donné par la
T e
loi de Fourier :  th  s . On en déduit la résistance thermique d’un échantillon : Rth  .
e s
La résistance thermique est d’autant plus élevée que l’épaisseur e est grande, le matériau peu
conducteur thermique et la surface s petite, ce qui est intuitif.

16) Le fait d’utiliser deux échantillons de part et d’autre de la zone de chauffe permet d’une
part d’éviter de faire un calorifugeage sur une seule zone latérale de la zone de chauffe et d’avoir
une plus grande incertitude dans l’évaluation du flux thermique traversant un unique échantillon.

Cela permet d’autre part de réduire l’incertitude sur la mesure de la conductivité thermique
puisque l’on fait une moyenne sur deux échantillons.

17) Si on suppose que la puissance thermique générée par la plaque chauffante est intégralement
transmise aux échantillons et que le régime stationnaire est bien atteint, il vient :

T S T S
  2SS   .
e e

2S
On en déduit alors :   .
T T

e e

Pour obtenir la formule précédente, nous avons fait l’hypothèse que toute la puissance
thermique générée par la plaque chauffante est intégralement transmise aux deux échantillons.

Afin de respecter expérimentalement cette hypothèse, il faut prendre quelques précautions

expérimentales pour négliger les résistances thermiques au niveau des contacts (plaque de
chauffe-échantillons) devant la résistance thermique de l’échantillon : on peut utiliser de la
graisse conductrice pour assurer un meilleur contact thermique et exercer une pression entre les
deux surfaces en contact.

Nous avons négligé les résistances thermiques de contact (plaque de chauffe-échantillons)

devant la résistance thermique des échantillons. Ce dispositif fonctionne donc essentiellement
pour des matériaux isolants.

4
18) On trouve :   3,79  102 W.K-1m-1 .

On trouve une conductivité thermique du même ordre de grandeur que celle de l’air, ce qui est
plutôt attendu pour un matériau dans lequel est emprisonné de l’air (c’est le cas aussi pour
d’autres matériaux isolants utilisés dans le bâtiment tel que la laine de verre, le polystyrène…).

19) En utilisant la question 4), la durée caractéristique du régime transitoire est

e2 e 2 cv
T  , soit de l’ordre de 4.103s ou encore de l’ordre de 1h. Il faut donc attendre
Dth 
plusieurs heures avant d’être surs expérimentalement d’avoir atteint le régime stationnaire. On
peut remarquer ici que nous n’avons évalué que la durée temporelle nécessaire à la
thermalisation des échantillons. En toute rigueur, Il faudrait également tenir compte de la
thermalisation de tout le dispositif, ce qui doit nécessiter beaucoup plus de temps.

En pratique, les expérimentateurs conditionnent les échantillons plusieurs jours.

II.3 Dispositif des cylindres coaxiaux

20) Il existe trois modes de transferts thermiques : le rayonnement, la conduction (ou diffusion
thermique) et la convection. Les transferts par convection sont accompagnés d’un mouvement
macroscopique de la matière et n’existent donc pas dans un solide. Les transferts par conduction
sont associés à un transport d’énergie sans déplacement macroscopique de la matière et ont
pour moteur l’agitation thermique des particules. Le rayonnement est un transfert d’énergie
dans un milieu transparent par l’intermédiaire d’un rayonnement électromagnétique (surtout les
infrarouges).

Bilan :

rayonnement conduction convection

solide Possible si transparent possible impossible
fluide Possible si transparent possible possible

Dans un fluide, contrairement au solide, les transferts thermiques peuvent être assurés par
convection (mouvement macroscopique de la matière possible).

21) On note T* l’échelle caractéristiques de la température et la grandeur l est plutôt la distance

caractéristique de variation de la température (et non de l’écoulement).

Les ordres de grandeurs caractéristiques cherchée sont :

T *
- pour le vecteur densité de courant thermique diffusif :
l
- pour le vecteur densité de courant thermique convectif : cvT *.V .

cvT *. Vl
On en déduit le nombre de Péclet thermique : Peth  Vl  .
T *
Dth

5
22) Ici l = r2-r1.
On utilise un faible espace annulaire pour minimiser le transport de chaleur par convection.

23) Le problème admet une invariance par rotation autour de l’axe (Oz). Comme on néglige les
effets de bord, il y a également invariance par translation parallèlement au même axe. On en
déduit : T (M )  T (r ) .

24) Appliquons, pendant une durée dt, le premier principe de la thermodynamique au volume
compris entre deux cylindres quelconques de rayons a et b (r1<a<b<r2) et de hauteur htot la
hauteur totale des deux cylindres coaxiaux : dU  0  ( (a)   (b))dt , puisque le régime est
stationnaire et puisqu’il n’y pas de production de chaleur.
On en conclut que le flux de chaleur radial traversant une surface cylindrique quelconque de
rayon r de hauteur fixée est constante, c’est-à-dire indépendant du rayon r de la surface.
T T A
Or  (r )   2rhtot , on en déduit  et par conséquent T (r )  A ln(r )  B , où A et B
r r r
sont des constantes réelles.

25) Par continuité des températures, on a T1  A ln(r1 )  B et T2  A ln(r2 )  B , d’où

T T T T
A  1 2 . Il vient  (r )  A2htot   1 2  2htot .
ln(r1 / r2 ) ln(r1 / r2 )

T1  T2
La résistance thermique du fluide contenu dans ce dispositif étant définie par Rth '  , on
 (r )
ln( r2 / r1 )
obtient pour la résistance thermique : Rth '  .
 2htot

 (r ) ln(r2 / r1 )
26) D’après la question précédente :   r1 .
T1  T2 Seff

On détermine le flux de chaleur radial  (r ) , en déterminant la puissance dissipée par effet Joule
dans le cylindre interne, à partir de la connaissance de la résistance chauffante Rchauf et de la
tension à laquelle elle est soumise Uchauf :

 (r ) = Uchauf 2/ Rchauf .

2
r1 ln(r2 / r1 ) U chauf
On en déduit   .
T1  T2 Seff Rchauf

Application numérique :  (r ) =29,91 W et   0,142 W.K-1m-1 .

6
  (r ) U chauf
27) On a 2 et (T1  T2 )  2T .
 (r ) U chauf

2 2 2 2
  (r )   (T1  T2 )   Uchauf   T 
D’où        4   2  .
   (r )   T1  T2   Uchauf   T1  T2 


On trouve  6,5  10 2 d’où   0,142  0,009 W.K-1m-1.


On remarque que l’incertitude sur la température est prépondérante.

II.3 La méthode flash

T T
28) Equation de diffusion de la chaleur dans le matériau :  Dth 2 .
t x

D’où, en réinjectant la solution cherchée dans l’équation de diffusion, il vient :

g ' (t ) f ' ' ( x)
f ( x) g ' (t )  Dth f ' ' ( x) g (t ) ou encore  Dth  cte , où la grandeur cte est une
g (t ) f ( x)
constante indépendante de x et t puisque les termes de droite et de gauche dépendent de deux
variables différentes.

Suivant le signe de la constante cte, on a des solutions différentes (trigonométriques ou

exponentielles) pour la fonction f(x).

Quel que soit le signe de cte, la fonction g(t) a la forme suivante : 𝑔(𝑡) = 𝐷𝑒 𝑐𝑡𝑒.𝑡 . Si cte>0, g(t)
diverge exponentiellement et si cte<0, g(t) décroît vers 0 exponentiellement. Etant donné que
l’on donne à l’échantillon une quantité de chaleur finie, seule la deuxième possibilité a un sens
physique. La constante cte a nécessairement un signe négatif et peut s’écrire : cte   K 2 Dth .

On en déduit la forme générale suivante pour f(x) et g(t) :

D K t 2
f ( x)  B cos(Kx)  C sin( Kx) et g (t )  De th .

29) Les faces avant et arrière étant isolées thermiquement, on peut écrire que le flux thermique
est nul en x=0 et x=L.

Or j th ( M , t )   grad (T )  f ' ( x) g (t )u x et f ' ( x)  KB sin( Kx)  CK cos(Kx) .

Les conditions aux limites s’écrivent donc : f ' (0)  CK  0 et f ' ( L)  KB sin( KL)  0 .

n
On en déduit C=0 et les valeurs quantifiées possibles pour la grandeur K ; K n  .
L

7
On en déduit pour la solution générale de la forme cherchée :
 n 
2

n  Dth   t
T ( x, t )  Ti  Bn cos( x )e  L 
.
L

30) On utilise le théorème de superposition (linéarité de l’équation de diffusion) et on cherche

 n 
2

n   Dth   t
maintenant à trouver une solution de la forme T ( x, t )  Ti   Bn cos( x)e  L  vérifiant
n 0 L
la condition initiale. : T ( x, t  0)  T0  Ti pour 0  x  d et T ( x, t  0)  Ti pour d  x  L .


n
Or T ( x,0)  Ti   Bn cos( x) .
n 0 L

n 
Dans l’expression de T (x,0) , le terme
n 0 L
B
x) correspond à un développement en
n cos(

série de Fourier d’une fonction, paire, 2L-périodique dont les harmoniques sont données par les
coefficients Bn .

Pour déterminer les coefficients Bn, il suffit d’identifier cette dernière fonction T (x,0) à la
fonction x  h(x) , définie dans l’énoncé, en prenant C  T0 .

n
d sin( d)
d L
On en déduit B0  T0 et Bn  2 T0
L L n
d
L

Il vient finalement pour le champ des températures cherché :

 n 2 

d   sin( d)
n
 n  
 Dth   t
T ( x, t )  Ti  T0 1   2 L cos( x )e  L   .
L  n 1 n  L 
 d 
 L 

Ce champ de température vérifie bien l’équation de la chaleur, les conditions aux limites ainsi
que la condition initiale.

d
31) Pour des temps infiniment longs, on obtient : T ( x, t  )  Ti  T0 et la température
L
d
d’équilibre est donc : Teq  Ti  T0 . Cette température est d’autant plus importante que
L
l’épaisseur d de la couche chauffée, l’élévation de température T0 due au flash sont élevées,
ce qui est assez intuitif.

8
32) Si on suppose que l’épaisseur d est suffisamment petite pour faire des calculs à l’ordre le
  n  
2
  Dth 
plus bas non nul en d/L on obtient : T ( L, t )  Ti  T0 1   2(1) n e  L   puisqu’alors
d  t

L  n1 
 
n
sin( d)    n  2 
 Dtht  
L  1 . D’où l’expression cherchée : T ( L, t )  Ti  (Teq  Ti ) 1  2 (1) e  L   .
k
n  n 1 
d  
L

33) Pour des temps longs, on peut ne considérer que le premier terme de la série infinie :
  Dtht   
  2

T ( L, t )  Ti  (Teq  Ti ) 1  2e  L   .
 
 

  Dth t1/ 2   
  2   Dtht1 / 2   
  2

Ainsi T ( L, t1/ 2 )  Ti  (Teq  Ti ) 1  2e   

L
 Ti  (Teq  Ti ) / 2  1  2e L 
 1/ 2
   
   

L2 ln(4) L2
 t1 / 2  . Ainsi t  0,1406 .
 2 Dth 1/ 2
Dth

34) Sur la figure 6, on constate que la courbe expérimentale a l’allure attendue. L’équilibre
thermique est atteint au bout de 0,5s environ. Toutefois, pour vérifier quantitativement que le
modèle théorique explique les observations expérimentales, il faudrait faire un ajustement de la
courbe expérimentale à l’aide d’une fonction de type f (t )  a  bect à l’aide d’un logiciel
adapté et pour des temps t suffisamment longs.

On peut également observer, que proche du flash, correspondant à l’instant t=0, la courbe
expérimentale ne peut pas être modélisée par une fonction f (t )  a  bect puisque la pente de
le courbe expérimentale à l’origine est nulle (la chaleur met un certain temps pour se propager
dans l’épaisseur du matériau). Ceci n’est pas en contradiction avec le modèle théorique
développée dans la question 33) puisque ce dernier n’est valable que pour des temps longs.

On peut observer que pour des temps supérieurs à 0,5s, il y a une petite décroissance de la
température, preuve d’une perte thermique dans le dispositif. L’hypothèse d’adiabaticité du
dispositif peut être remise en question.

35) Sur la figure 5, on mesure t1/ 2  (0,116  0,005) s, soit Dth  (4,7  0,1)  10 6 m2.s-1 . C’est
l’incertitude sur l’épaisseur L de l’échantillon qui limite la précision sur Dth.

36) Quatre sources d’erreurs majeures peuvent être responsables des incertitudes de mesure de
la méthode flash en utilisant le modèle décrit dans ce sujet. Tout d’abord, les pertes thermiques
radiatives et convectives de l’échantillon ont été négligées. L’inhomogénéité du flux thermique
incident doit être considéré. Il faudrait également tenir compte de la durée et la forme exacte de
l’impulsion flash et enfin considérer des effets de non linéarité.

9
De nombreux autres modèles ont été développés ces dernières années, tenant compte de la
dissipation de chaleur dans le dispositif expérimental.

Il s’agit d’une méthode de mesure rapide ne nécessitant pas l’attente de l’installation du régime
stationnaire. Cette méthode est facile à mettre en oeuvre et peut être utilisée pour un grand
nombre de matériaux (isolant, conducteur, solide, liquide, granulaire) sur des grandes gammes
de températures.

37) En appliquant le premier principe de la thermodynamique sur l’ensemble de l’échantillon

entre l’instant t=0- et un instant où l’équilibre thermique final est atteint, il vient :
U  CT  C (Teq  Ti )  Q , où Q est la chaleur reçue par l’échantillon au cours du flash. On
peut effectuer la même mesure avec le même flash, délivrant la même chaleur Q, sur un autre
échantillon de référence, dont on connait la capacité thermique, et en déduire la chaleur Q
fournie par le flash au cours de la mesure. On peut alors déterminer expérimentalement une
mesure de la capacité thermique de l’échantillon étudié.

Connaissant, la diffusivité thermique de l’échantillon, sa capacité thermique, sa masse

volumique et son volume, il est alors possible d’en déduire une mesure de sa conductivité
thermique.

III. Diffusion de la quantité de mouvement

III.1 Force de viscosité et approche microscopique de la diffusion de la quantité de

mouvement

38) Exemples de fluides non newtoniens : ketchup, mayonnaise, dentifrice, peinture, mélange
d’eau et de maïzena, mousse à raser, yaourt, ciment. ….

A l’échelle microscopique, le caractère non newtonien d’un fluide peut s’expliquer par la
présence dans le fluide d’objets de grande taille par rapport à l’échelle atomique mais de
dimension petite par rapport à la dimension de l’écoulement : exemples de particules en
suspension, des macromolécules dans les solutions de polymères, des gouttelettes dans les
émulsions. Ces objets peuvent eux-mêmes former des structures encore plus grandes qui
influencent de manière importante les propriétés de l’écoulement.

39) L’expression de la force de viscosité exercée sur un élément de surface dS  dS u z , situé à

la côte z, par le fluide situé au-dessus de z sur le fluide situé en dessous de z en fonction du
v ( z )
champ des vitesses v x (z ) est : d F   x dS u x .
z

Le coefficient de viscosité dynamique a pour dimension :    MT 1L1 .

On utilise souvent le Poiseuille (Pl) comme unité.

10
La viscosité dépend bien sûr du fluide considéré, de la température et de la pression.

Dans le cas du glycérol pur, l’ordre de grandeur de la viscosité dynamique est 1Pl.

40) Initialement, lorsque la plaque du dessus est fixe, tout le fluide est au repos. Lorsque la
plaque est mise en mouvement brutalement, seule la couche de fluide au-dessous de la plaque
accroche le solide et se met en mouvement. Il y a diffusion de la quantité de mouvement dans
la mesure où c’est l’hétérogénéité du champ des vitesses qui est à l’origine de la mise en
mouvement de toutes les couches de fluide et cette diffusion a pour rôle d’homogénéiser le plus
possible le champ des vitesses, compte tenu des conditions aux limites. Cette diffusion
s’effectue dans la direction (Oz), perpendiculairement à l’écoulement.

En régime stationnaire, on peut comparer ce type d’écoulement au cas de la conduction

unidimensionnelle de la chaleur entre deux plaques parallèles portées à des températures
différentes, le champ de vitesse vx(z) jouant alors le rôle du champ des températures T(z). A la
vx
relation de Fourier se substitue une relation de proportionnalité entre et la force de
z
viscosité surfacique exercée sur un élément de surface dS  dS u z , situé à la côte z, par le fluide
situé au-dessus de z sur le fluide situé en dessous de z.

41) Pendant la durée , le nombre de molécules venant du haut et traversant la surface dS à la

1
cote z est donné par nudS .
6

On en déduite la quantité de mouvement moyenne dans la direction (Ox) traversant pendant 

l’unité de surface dS à la cote z associé aux molécules venant du haut:
1
dP haut x  nudS (mv( z  l ))
6

De même, la quantité de mouvement moyenne dans la direction (Ox) traversant pendant 

l’unité de surface dS à la cote z associée aux molécules venant du bas:

1
dP bas x  nudS (mv( z  l )) .
6

1
On en déduit : dPx  dP x  dP x  nudS (mv( z  l )  mv( z  l ))
haut bas

42) Comme le libre parcours moyen l est très inférieur à la distance caractéristique de variation
du champ de vitesse, on peut écrire :

dFx 1 dPx 1 1 v 1 v
  nu (mv( z  l )  mv( z  l ))  num .2l  numl .
dS  dS 6 6 z z 3 z z

1
On en déduit :   numl .
3

11
43) L’énergie cinétique de la molécule étant associé à trois degrés de liberté, on a
3 1
k BT  mu 2 .
2 2

1 1 3k BT 1 k BT m
D’où   numl  nm . Il vient :   .
3 3 m 2n C 6 C

On constate que la viscosité dynamique, dans ce modèle, est indépendant de la concentration

moléculaire. Cette viscosité augmente lorsque la température augmente. Cette situation est
opposée à celle des liquides pour lesquels, la viscosité diminue lorsque la température croit,
suivant une loi d’Arhénius qui décrit le processus d’activation de saut à travers la barrière de
potentiel d’interaction séparant deux molécules

44) On trouve Pl dans cas du gaz d’hélium à température ambiante.

III.2 Equation de diffusion

45) Le nombre de Reynolds est un nombre sans dimension caractérisant un écoulement et

représentant le rapport du flux convectif de quantité de mouvement sur le flux diffusif de
quantité de mouvement (analogie avec le nombre de Péclet thermique).

U * L*
Re 


46) On part de l’équation de Navier Stokes, en négligeant toutes les autres forces que les forces
 v 
de viscosité :    (v.grad )v   v . Si Re<<1 alors on peut négliger  (v.grad )v devant la
 t 
v v 
force volumique de viscosité et il vient   v ou encore  v  v  Dvisc v .
t t 


Ainsi Dvisc 


U * L*
Il n’est pas tout le temps nécessaire d’avoir la condition Re  <<1 pour obtenir

l’équation de diffusion précédente. Il existe des écoulements (comme les écoulements de
cisaillement) pour lesquels l’accélération convective est nulle.

12
Ordres de grandeurs :

glycérol pur eau air

Pl) 1 1. 10-3 2.10-5
kg.m-3 1.103 1.103 1
Dvisc m2.s-1 10-3 10-6 2.10-5

La diffusion de la quantité de mouvement est un peu plus efficace dans l’air que dans l’eau,
même si l’eau est un fluide plus visqueux que l’air.

  1
47) En reprenant le modèle microscopique : Dvisc    ul .
 nm 3

On trouve le même type de formule que pour le coefficient de diffusion thermique. Ce résultat
illustre l’unicité du mécanisme microscopique (agitation thermique) associé aux phénomènes
de transport diffusif dans les gaz dilués.

III.3 Frottement fluide sur une plaque en mouvement rectiligne uniforme

48) La plaques étant supposée infinie, il y a invariance du système par translation parallèlement
aux axes (Ox) et (Oy). On en déduit que le champ de pression ne dépend que de la coordonnée
z : P(M,t)=P(z,t).

Il en est de même pour le champ des vitesses : v( M , t )  v( z , t ) . L’écoulement étant par ailleurs
de cisaillement, on en déduit : v( M , t )  vx ( z , t )u x .

z
49) On cherche un champ des vitesses de la forme v( M , t )  f (u )u x avec u  un
Dvisct
paramètre sans dimension.

 v 
On écrit l’équation de Navier Stockes :    (v.grad )v    grad ( P)   g  v . En
 t 
 
projection sur (Ox), il vient :

 f (u )   2 f (u ) 1 z 1 1
   Dvisc   f ' (u )  Dvisc f ' ' (u )
 t  z 2
2 Dvisct t Dvisct

1
 f ' (u )u  f ' ' (u )  0
2

13
50) En posant F(u)=f’(u), l’équation différentielle vérifiée par la fonction F est
1
F (u )u  F ' (u )  0 , dont la solution générale est F (u )  F0e  u / 4 .
2

2
u
51)  F ( X )dX  f (u)  f (0) d’où
0

u u/2
v( z, t )     F0e u 2 / 4
du    2  F0e  z dz     F0erf (u / 2) .
2

0 0

Pour déterminer les constantes  et F0, il suffit de faire intervenir les conditions aux limites en z=0..
et z   .

Pour t>0, v(0, t )     erf (0 / 2)    U 0 et

lim z  v( z , t )     F0 lim u  erf (u / 2)  0     F0 .

D’où   U 0 et F0   /  .

z
On en conclut : v( z , t )  U 0 (1  erf ( ))
2 Dvisct

La fonction v(z,t) ne dépend que de la variable u. Tous les profils de vitesses se déduisent donc
les uns des autres par dilatation de l’échelle des longueurs suivant (Oz) d’un facteur 2 Dvisct .
A un instant t donné, l’épaisseur caractéristique sur lequel le champ des vitesses varie, ou encore
la distance sur laquelle le déplacement de la plaque se fait sentir est de l’ordre de d= 2 Dvisct .
Comme tout phénomène diffusif, la diffusion de la quantité de mouvement est peu efficace pour
les temps longs.

52)

t=10 s t=1 h
eau d=6 mm d=12 cm
air d=3 cm d=53 cm

La mise en mouvement d’un fluide par couplage visqueux est d’autant plus rapide que Dvisc est
grande c’est-à-dire que sa viscosité dynamique  est grande et sa masse volumique  petite. De
ce point de vue, comme on l’avait remarqué plus haut, la diffusion de la quantité de mouvement
s’effectue un peu mieux dans l’air que dans l’eau.

14
53) La force df subie par un élément de surface dS de la plaque mobile a pour expression :
v 1 1 1  U0
 u x dS  dSf ' (u ) u x  dSF (u ) u x  dS ux .
z z  0 Dvisct Dvisct Dvisct 

1  U0
Ainsi : d f  dS ux .
Dvisct 

Il s’agit bien sûr d’une force qui s’oppose au mouvement de la plaque. Cette force diminue
lorsque le temps augmente (puisque le gradient du champ des vitesses diminue grâce à la
diffusion de la quantité de mouvement).

54) Ce dispositif est rigoureusement identique à celui de la diffusion de la chaleur dans un demi-
espace à partir d’un plan dont la température est fixée et subit une variation brutale de
température à t=0. Le champ de température est alors l’analogue de v(z,t) et Dvisc est remplacé
par Dther.

Ce dispositif est également l’analogue de la diffusion de particules en régime non stationnaire

entre deux réservoirs contenant deux concentrations de particules différentes.

III.4 Viscosimètre de Couette

55) Le cylindre en rotation met en mouvement la couche de fluide juste au-dessous de lui. Puis
il y a diffusion de la quantité de mouvement dans la direction radiale car les couches de fluide
successives se mettent en mouvement par l’intermédiaire de la force de viscosité. La couche de
fluide en contact avec le cylindre interne met en mouvement ce dernier. Comme ce dernier est
retenu par un couple de force de moment <0 par rapport à l’axe (Oz), en régime permanent,
ce cylindre est à l’équilibre, le couple  compensant la force de viscosité s’exerçant de la part
du fluide sur le cylindre. En mesurant l’angle dont a tourné le cylindre interne, on peut
déterminer la viscosité du fluide étudié.

Le transport diffusif de la quantité de la quantité de mouvement s’effectue radialement. Le

transport convectif de la quantité de mouvement se fait quant à lui, dans la direction de
l’écoulement, c’est-à-dire dans la direction orthoradiale.

56) Il y invariance par rotation autour de l’axe (Oz). Le champ de pression ne dépend donc que
des coordonnées cylindriques r et z : P(M)=P(r,z). De même v( M )  v(r , z ) .

57) On considère le référentiel du laboratoire galiléen. Les forces extérieures s’appliquant sur
le système sont : les forces de viscosité, les forces de pression et la force de pesanteur. Seules
les forces de viscosité ont un moment non nul par rapport à l’axe (Oz). En régime permanent le
moment cinétique du système par rapport à (Oz) est constant.

15
Il vient :
dL  v   v     v 
 0  r 2    2rH   (r  dr ) 2    2 (r  dr ) H  2H  r 3     .
dt r  r  r r  r  r dr r  r  r 

dL
0
dt

  3   v 
D’où r     0 et en intégrant par rapport à la variable r, on obtient :
r  r  r 
A
v (r )    Br .
2r

Pour déterminer les constantes A et B, il suffit ensuite de réinjecter les conditions aux limites :
v ( R1 )  0 et v ( R2 )  R2 .

R2 (r 2  R1 )
2 2

Après calculs, on obtient le champ des vitesses suivant : v (r )  .

r ( R2  R1 )
2 2

58) Appliquons maintenant le théorème du moment cinétique au cylindre interne par rapport à
  v 
l’axe (Oz) lorsque celui-ci est en équilibre :   R1   2R1H  0 .
2

r  r  R1

R R
2 2
A
D’où,   H 2R  4H 2 1 2 2 .
3

R2  R1
3 1
R1

R2  R1
2 2

On en déduit    2 .
4HR1 R2
2

59) Application numérique = 6,12 102 Pl.

60) On veut retrouver l’écoulement laminaire entre deux plaques planes dont celle du dessus
est en mouvement et celle du dessous est immobile, c’est-à-dire l’écoulement de Couette plan
décrit dans la partie III.1, pour lequel le profil des vitesses est linéaire en fonction de z.

On pose donc r= R1+r avec r<<R1.

R2 (r 2  R1 ) R2  R2  R2 
2 2 2 2 2 2 2
R R1
v (r )   (r  1 )  ( R1   (r) - ) 2 (r )
r ( R2  R1 )
2 2
( R2  R1 )
2 2
r ( R2  R1 )
2 2
R1   (r) ( R2 2  R12 )
R2  R2  R
2 2
v (r )  2 (r )  2 (r )  2  (r ) .
( R2  R1 )
2 2
eC ( R1  R2 ) eC

On retrouve bien le même profil des vitesses que celui dans la partie III.1 avec les
correspondances suivantes : R2   V0 , eC  a et  (r )  z .

16
61) On trouve les deux types de dispositifs dans les rhéomètres industriels. A grande vitesse de
rotation, le rhéomètre avec cylindre interne en rotation peut être moins stable à cause des
instabilités de Taylor-Couette.

IV Diffusion de la matière

IV. 1 Loi de Fick et diffusion de particules

62) La loi de Fick permet de relier le gradient de concentration au flux de particules d’une
espèce. Elle s’écrit :

𝑗⃗ = −𝐷mol ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗
grad𝑐

où, comme précisé dans l’énoncé, j s’exprime en en mol.m².s-1, D en m²/s et c en mol/m3.

63) On peut penser aux expériences suivantes :

- Diffusion d’une vapeur d’ammoniac dans un tube placé verticalement. Ceci conduit à un ordre
de grandeur.

- La diffusion d’une solution de glycérol dans de l’eau (expérience avec la nappe laser).

IV. 2 Evaporation d’une goutte d’eau sphérique

Cas a) Evaporation d’une goutte d’eau sphérique en atmosphère humide.

64) En électrostatique dans une région vide de charges, le potentiel électrostatique V vérifie
∆𝑉 = 0 et pour une distribution de charges admettant la symétrie sphérique, en dehors de cette
𝐴
distribution, le potentiel prend la forme : 𝑉(𝑟) = 𝑟 + 𝐵.

Par analogie, la solution de l’équation ∆𝑐 = 0 ayant une symétrie sphérique est de la forme en
𝐴(𝑡)
régime quasi-stationnaire :𝑐(𝑟, 𝑡) = + 𝐵(𝑡),
𝑟

On insère les conditions aux limites 𝑐(𝑟𝑔 , 𝑡) = 𝑐𝑠𝑎𝑡 et 𝑐(∞, 𝑡) = 𝐻𝑐𝑠𝑎𝑡 . On obtient alors :

𝑟𝑔 (𝑡)
𝑐(𝑟, 𝑡) = 𝑐𝑠𝑎𝑡 (1 − 𝐻) + 𝐻𝑐𝑠𝑎𝑡 .
𝑟

65) La variation de masse de la goutte par unité de temps peut se calculer à partir du nombre de
particules qui s’échappent de la surface :

𝐼(𝑡) = 4𝜋𝑟𝑔 ²𝑀𝑗(𝑟𝑔 ) = 4𝜋𝑟𝑔 𝑀𝐷mol 𝑐𝑠𝑎𝑡 (1 − 𝐻)

On observe de façon non totalement intuitive que le taux instantané d’évaporation n’est pas
proportionnel à la surface de la goutte à un instant donné mais plutôt à la racine carrée de la
surface de la goutte.

17
 4
𝑑𝑚 𝑑( 𝜇𝜋)
66) On écrit naturellement = 3
𝑟𝑔3 .
𝑑𝑡 𝑑𝑡

On en déduit :

2 2𝑀 2𝑀
𝑟𝑔2 (𝑡) = 𝑟𝑔0 − 𝐷mol 𝑐𝑠𝑎𝑡 (1 − 𝐻)𝑡 ou 𝑟𝑔2 = 𝐷mol 𝑐𝑠𝑎𝑡 (1 − 𝐻)(𝑡𝑓 − 𝑡).
µ µ

En représentation log-log, on s’attend à ce que la pente de la droite soit égale à 0,5. En effet :

1 2𝑀 1
log (𝑟𝑔 (𝑡)) = 2 log ( 𝐷mol 𝑐𝑠𝑎𝑡 (1 − 𝐻)) + 2 log(𝑡𝑓 − 𝑡).
µ

1000
8
6
rg [µm]

100
2 3 4 5 6 7 8 9
1 10
tf-t [min]

Sur le graphe, la pente de la droite en trait continu est 0.5 comme attendu.

On cherche alors à déterminer la valeur du coefficient de diffusion. On lit que pour 𝑡𝑓 − 𝑡 =

1 min, on a 𝑟𝑔 = 170 µm. On calcule 𝑐sat ~ 0,94 m-3. On trouve 𝐷mol =3,0 m²/s.

67) Pour réaliser une telle expérience, on peut attacher une gouttelette d’eau à un capillaire très
fin. La forme d’une goutte sur un capillaire très fin est essentiellement sphérique car, même en
situation de mouillage total, il y a peu de gain énergétique à s’étaler. Une autre solution est de
placer la goutte en lévitation acoustique.

68) Pour utiliser l’hypothèse d’un profil stationnaire établi, il faut comparer la vitesse de
l’interface 𝑑𝑟𝑔 /𝑑𝑡 à l’échelle de vitesse 𝐷mol /𝑟𝑔 associée à la fabrication du profil de
concentration autour de la goutte. La première vitesse est typiquement de 100.10-6/(10*60)~10-
7
m/s. La seconde vitesse est 10-5/100.10-6=0,1 m/s qui est bien plus rapide. L’approximation
du régime quasi stationnaire était donc bien justifiée.

18
Cas b) Evaporation d’une goutte d’eau sphérique en atmosphère saturée

69) La tension de surface est due au fait que l’environnement moléculaire d’une molécule à la
surface n’est pas la même qu’une molécule dans le volume du liquide. En conséquence, il faut
apporter de l’énergie pour placer une molécule à la surface.

Typiquement, on s’attend donc à ce que la tension de surface soit de l’ordre de l’énergie d’une
liaison intermoléculaire par taille de molécule :
𝐸 𝑘𝑇 4.10−21
𝛾~ 𝑎² ~ 𝑎² ~ (3.10−10 )² ~0.04 J.m-2 .

Pour l’eau, les liaisons sont un peu plus énergétiques (liaisons H plutôt que liaisons de van der
Waals, d’où un nombre un peu plus élevé).

Quelques ordres de grandeur : eau 72 mN/m, éthanol : 22 mN/m, eau savonneuse : 30 mN/m.

70) La différence de pression des deux côtés d’une interface courbée est liée à la courbure
moyenne de cette interface On a la relation :
2𝛾
∆𝑝 = .
𝑅

2𝛾
On en déduit ici que : 𝑝𝑙 − 𝑝𝑣 (𝑟𝑔+ ) = .
𝑟𝑔

𝑑µ
71) On a la relation : (𝑑𝑝) = 𝑣. Si on fait l’hypothèse que le liquide est incompressible :
𝑇

µ𝑙 (𝑝𝑙 , 𝑇) = µ𝑙 (𝑝sat , 𝑇) + 𝑉m (𝑝𝑙 − 𝑝sat )

72) En faisant l’hypothèse que le gaz est un gaz parfait, on obtient l’équation (9). On égalise
les deux potentiels chimiques et on obtient :

𝑝v 2𝛾
𝑅𝑇ln ( ) = 𝑉m (𝑝𝑙 − 𝑝sat ) = 𝑉m (𝑝𝑣 − 𝑝sat + )
𝑝sat 𝑟𝑔
2𝛾
pour une goutte d’eau de 1 µm vaut typiquement 105 Pa, ce qui est beaucoup plus grand que
𝑟𝑔
la pression de vapeur saturante de l’eau (quelques kPa).

𝑝 𝑝v −𝑝sat 2𝛾 2𝛾
D’où : 𝑅𝑇ln (𝑝 v ) ~𝑅𝑇 = 𝑉m ( 𝑟 ). Ainsi : 𝑝v − 𝑝sat = 𝑝sat 𝑉m (𝑅𝑇𝑟 ).
sat 𝑝sat 𝑔 𝑔

73) On voit que 𝑝𝑣 > 𝑝sat donc une goutte d’eau peut s’évaporer, même en présence d’air
totalement humide.

74) En refaisant le même calcul que précédemment, on peut calculer le taux d’évaporation 𝐼(𝑡) :

8𝜋𝛾𝑀𝐷mol 𝑝sat 𝑉m
𝐼(𝑡) =
(𝑅𝑇)²

19
Dans ce cas-là, on trouve 𝐼 = cte. On déduit donc que le temps de vie est :

4
µ 3 𝜇𝜋𝑟𝑔3
Δ𝑡 =
𝐼
On trouve pour 𝐷mol = 1,0. 105 m²/s, Vm = 18.10-6 m3, I=2,1.10-15 kg/m3, d’où une durée de vie
pour une goutte de masse 5,2.10-19 kg : Δ𝑡 = 0,2 ms.

75) Les effets de courbure ne sont des effets qui n’apparaissent que pour des courbures
nanométriques ou très très proche de la saturation. On a donc eu raison de négliger ces effets
de correction de courbure.

IV. 3 Mouvement brownien

76) Une particule colloïdale est une particule dans un solvant de taille comprise entre quelques
nanomètres et quelques micromètres afin de ne pas sédimenter et d’être soumises au
mouvement brownien.

77) Les équations (12) et (13) signifient qu’un premier effet systématique du solvant a été pris
en compte sous la forme d’un terme systématique qui est le terme de friction. Les équations (12)
et (13) signifient alors d’une manière mathématique le caractère aléatoire des collisions des
molécules du solvant sur la particule colloïdale.

Une hypothèse importante du calcul est le fait que les temps caractéristique des chocs sur la
particule est bien plus faible que le temps d’évolution des vitesses.

Temps caractéristique des chocs :

La vitesse moyenne d’une molécule d’eau à 20 °C peut être estimée à partir du théorème
d’équipartition :

1 3
𝑚𝑣 2 = 𝑘𝑇
2 2
Pour une molécule d’eau (𝑀 = 18 g/mol), cela donne 𝑣~640 m/s. Donc, pour une molécule
de rayon 𝑏, le temps typique de collision des molécules entre elles est 𝑏/𝑣. La taille typique
d’une molécule d’eau peut être obtenue en écrivant µ = 1 g/cm3, ce qui conduit à un volume
moléculaire de 1,5.10-29 m3 et donc un rayon de molécule de 0,2 nm.

On trouve donc :
𝑏
𝜏~ 𝑣 ~30 ps

Le temps de décorrélation des vitesses est typiquement le temps de diffusion d’une particule
colloïdale sur sa taille. Pour une particule de 1 µm ce temps vaut typiquement 1 ms. Pour une
particule de 10 nm, ce temps est de 0.001 ms.

20
Sur les temps typiques d’intérêt, les forces exercées par le solvant sur la particule sont des
variables indépendantes décorrelées.

78) On peut supposer que la friction sur la particule est une friction de type Stokes, i.e. :

6𝜋𝜂𝑎 = 𝑚𝜉, où a est le rayon de la particule.

79) On peut résoudre l’équation de Langevin par la méthode de variation de la constante :

1 𝑡
𝑣⃗(𝑡) = 𝑣⃗(0) exp(−𝜁𝑡) + ∫ 𝑓⃗(𝑠)exp (−𝜁(𝑡 − 𝑠)) 𝑑𝑠
𝑚 0

On déduit simplement que

1 𝑡 𝑡
〈𝑣⃗(𝑡). 𝑣⃗(𝑡)〉 = 〈𝑣02 〉 exp(−2𝜁𝑡) + ∫ 𝑑𝑠′ ∫ 𝑑𝑠exp (−𝜁(2𝑡 − 𝑠 − 𝑠′)) 〈𝑓⃗(𝑠). 𝑓⃗(𝑠′)〉
𝑚² 0 0

Car 〈𝑓⃗(𝑠). 𝑣⃗(0)〉 s’annule.

80) On en déduit ainsi en utilisant les propriétés du bruit que la limite à grand t vaut :

3𝑔0
lim 〈𝑣²(𝑡)〉 =
𝑡→∞ 𝑚²𝜁

En utilisant le théorème d’équipartition de la mécanique statistique d’équilibre, on sait que le

terme de droite de cette équation vaut : 3𝑘𝑇/𝑚 soit :

𝑚𝜁𝑘𝑇 = 𝑔0

Ceci signifie qu’il existe une relation simple entre l’amplitude de la force fluctuante 𝑔0 et la
dissipation 𝜁 .Cette relation thermodynamique montre qu’il existe une relation entre la
dissipation caractérisée par la viscosité à l’échelle macroscopique et la force fluctuante
microscopique.

81) On cherche à calculer 〈𝑣⃗(𝑡). 𝑣⃗(𝑡′)〉. Le système est supposé à l’équilibre thermodynamique.
Il est donc invariant par symétrie temporelle. On a donc :

〈𝑣⃗(𝑡). 𝑣⃗(𝑡′)〉 = 〈𝑣⃗(𝑂). 𝑣⃗(𝑡 − 𝑡′)〉.

On en déduit alors pour 𝑡 − 𝑡’ > 0 :

〈𝑣⃗(𝑡). 𝑣⃗(𝑡′)〉 = 〈𝑣0 ²〉 exp(−𝜁(𝑡 − 𝑡 ′ ))

En utilisant le théorème d’équipartition, on tire :

3𝑘𝑇
〈𝑣⃗(𝑡). 𝑣⃗(𝑡′)〉 = exp(−𝜁|𝑡 − 𝑡′|)
𝑚
On en déduit que 𝜁 −1 est le temps typique au-delà duquel la particule perd la mémoire de sa
vitesse. C’est le « temps de relaxation de la vitesse ».

21
82) On développe :

〈(𝑟⃗(𝑡) − 𝑟⃗(𝑡′))²〉 =< 𝑟 2 (𝑡) > +< 𝑟 2 (𝑡′) > −2 < 𝑟⃗(𝑡). 𝑟⃗(𝑡′) >

En notant que :
𝑡 𝑡′
〈𝑟⃗(𝑡). 𝑟⃗(𝑡′)〉 = ∫ 𝑑𝑠 ∫ 𝑑𝑠′〈𝑣⃗(𝑠). 𝑣⃗(𝑠′)〉
0 0

On tire :

6𝑘𝑇
〈(𝑟⃗(𝑡) − 𝑟⃗(𝑡′))²〉 = (|𝑡 − 𝑡′| + 𝜁 −1 (exp(−𝜁|𝑡 − 𝑡′|) − 1))
𝑚𝜁

On en déduit l’allure de graphe suivante :

〈(𝑟⃗(𝑡) − 𝑟⃗(0))²〉

𝑡
𝜁 −1
On observe à temps court (𝑡 ≪ 𝜁 −1 ) un régime balistique. A temps, long (𝑡 ≫ 𝜁 −1 ), un régime
diffusif. On identifie

6𝑘𝑇
6𝐷mol =
𝑚𝜁
𝑘𝑇
D’où 𝐷mol = 𝑚𝜁 .

𝑘𝑇
83) 𝐷mol = 6𝜋𝜂𝑎 = 4,4. 10−11 m.s-2. C’est le bon ordre de grandeur.

22