Leçon

 La diffusion moléculaire d'un point de vue microscopique et macroscopique. Les molécules solubles sur le côté gauche de la barrière (ligne violette) diffusent pour remplir le volume complet.

En haut : une seule molécule se déplace aléatoirement.

Au milieu : Le soluté remplit le volume disponible par marche aléatoire.

En bas : au niveau macroscopique, le côté aléatoire devient indétectable. Le soluté se déplace des zones où les concentrations sont élevées vers les zones à concentrations plus faibles. Ce déplacement est décrit par la loi de Fick.


 

La loi de Fick décrit la diffusion de la matière dans un milieu binaire. Elle a été établie par Adolf Fick en 1855.

Reliant le flux de matière au gradient de concentration, elle est analogue à l'équation de la chaleur introduite par Joseph Fourier  en 1822. Ce type de loi nommée loi de diffusion en mathématiques apparaît dans les systèmes décrivant un transport (masse, énergie, etc.) chaque fois que l'on peut séparer les échelles microscopiques d'un phénomène décrit par une équation cinétique comme l'équation de Boltzmann et les échelles du milieu continu macroscopique.

Cette loi, au départ empirique, a été justifiée et généralisée dans le cas d'un milieu multicomposant sous le nom de équations de Stefan-Maxwell d'après les travaux de James Clerk Maxwell pour les gaz en 1866 et Josef Stefan pour les liquides  en 1871.

La loi exprime une relation linéaire entre le flux de matière et le gradient de concentration de celle-ci :

J j = ρ D i j c j {\displaystyle \mathbf {J} _{j}=-\rho {\mathcal {D}}_{ij}\nabla c_{j}}

avec

J j {\displaystyle \mathbf {J} _{j}} flux massique (kg.m-2.s-1),
ρ {\displaystyle \rho } masse volumique (kg.m-3),
D i j {\displaystyle {\mathcal {D}}_{ij}} coefficient de diffusion binaire (m2.s-1),
c j {\displaystyle c_{j}} fraction massique (sans unité).

Propriété

Les quantités contenues dans l'équation sont telles que D j i = D i j {\displaystyle {\mathcal {D}}_{ji}={\mathcal {D}}_{ij}} (symétrie de l'interaction entre les particules i et j) et c i + c j = 1 {\displaystyle c_{i}+c_{j}=1} (par définition de la fraction massique). On en déduit que la diffusion ne transporte pas globalement de masse, elle ne fait que répartir différemment celle-ci :

J i + J j = 0 {\displaystyle \mathbf {J} _{i}+\mathbf {J} _{j}=0}

Cette propriété résulte en fait de la définition de la vitesse d'un fluide comme la vitesse d'ensemble (vitesse barycentrique, appelée généralement "vitesse", sans qualificatif) transportant globalement la masse et de la vitesse de diffusion transportant une composante de celle-ci par rapport au barycentre.

Pour un soluté

On peut exprimer cette loi sous une autre forme pour un milieu incompressible où ρ = 0 {\displaystyle \nabla \rho =0} en divisant par la masse molaire M j {\displaystyle M_{j}} (kg.mol-1) du soluté :

j j = D i j C j {\displaystyle \mathbf {j} _{j}=-{\mathcal {D}}_{ij}\nabla C_{j}}

avec

j j = J j M j {\displaystyle \mathbf {j} _{j}={\frac {\mathbf {J} _{j}}{M_{j}}}}  : flux molaire (mol.m-2.s-1),
C j = ρ c j M j {\displaystyle C_{j}={\frac {\rho c_{j}}{M_{j}}}}  : concentration molaire (mol.m-3).

On note que j i + j j 0 {\displaystyle \mathbf {j} _{i}+\mathbf {j} _{j}\neq 0}

Équation de conservation

Article détaillé : Équation de conservation.

On peut définir une loi de conservation pour une variable extensive ϕ {\displaystyle \phi } entraînée à la vitesse V {\displaystyle \mathbf {V} } et comportant un terme de production volumique S {\displaystyle S} par :

ϕ t + ( ϕ V ) = S {\displaystyle {\frac {\partial \phi }{\partial t}}+\nabla \cdot (\phi \mathbf {V} )=S}

Dans notre cas on prend ϕ = ρ c i {\displaystyle \phi =\rho c_{i}} , V = J i ρ c i {\displaystyle \mathbf {V} ={\frac {\mathbf {J} _{i}}{\rho c_{i}}}} et S = 0 {\displaystyle S=0} , ce qui donne dans le cas général :

t ( ρ c i ) + J i = 0 {\displaystyle {\frac {\partial }{\partial t}}(\rho c_{i})+\nabla \cdot \mathbf {J} _{i}=0}

Dans le cas d'un liquide ou plus généralement d'un fluide incompressible, en divisant par M :

C i t + j i = 0 {\displaystyle {\frac {\partial C_{i}}{\partial t}}+\nabla \cdot \mathbf {j} _{i}=0}

Cette loi de conservation est parfois appelée "seconde loi de Fick". Elle est en tout point analogue à l'équation de la chaleur. On dispose donc pour l'analyse de tout l'arsenal théorique et numérique liée à celle-ci.