Équations de Stefan-Maxwell

Les équations de diffusion de Stefan-Maxwell décrivent la diffusion dans un milieu multiespèce. Elles ont été établies indépendamment par James Clerk Maxwell (1866) pour les gaz peu denses et Josef Stefan (1871) pour les liquides.

Première version des équations

Il n'est pas possible en général d'exprimer de manière explicite les flux de diffusion en fonction des gradients de concentration. La relation entre ces quantités est donnée par le système linéaire

$\sum _{j\ne i}\frac{x_i{x}_j}{\rho {𝒟}_{ij}}(\frac{{𝐉}_j}{c_j}-\frac{{𝐉}_i}{c_i})=x_i\frac{\mathrm{\nabla }{\mu }_i}{RT}=\mathrm{\nabla }{x}_i,i,j=1,N$

où

• ${𝐉}_i={\rho }_i{𝐕}_{D_i}$ est le flux massique de diffusion pour l'espèce i (Modèle:Unité),
• ${𝐕}_{D_i}$ est la vitesse de diffusion (Modèle:Unité),
• $x_i$ la fraction molaire ou volumique,
• $c_i$ la fraction massique,
• ${𝒟}_{ij}$ le coefficient de diffusion binaire (Modèle:Unité), quantité strictement positive,
• $\rho =\sum _i{\rho }_i$ la masse volumique (Modèle:Unité),
• ${\rho }_i=n_i{m}_i$ où $n_i$ est la densité volumique de particules et $m_i$ leur masse,
• ${\mu }_i$ le potentiel chimique (Modèle:Unité),
• $R$ la constante universelle des gaz,
• $T$ la température,
• $N$ est le nombre d'espèces présentes dans le milieu.

Ce système est d'ordre $N$ mais de rang $N-1$ puisque par définition de la notion de diffusion

$\sum _i{𝐉}_i=0$

Dans le cas d'un système binaire ce système se résout immédiatement et conduit à la loi de Fick

$J_1=-\rho {𝒟}_{12}\mathrm{\nabla }{c}_1$

$J_2=-\rho {𝒟}_{21}\mathrm{\nabla }{c}_2=-J_1$

Ceci implique que

${𝒟}_{12}={𝒟}_{21}$

Cette symétrie est naturelle lorsqu'on regarde l'interaction ij au niveau microscopique.

On peut résoudre formellement le système linéaire :

${𝐉}_i=-\frac{\rho }{{\overline{M}}^2}\sum _{j\ne i}{M}_i{M}_j{D}_{ij}\mathrm{\nabla }{x}_j$

où

• $M$ est la masse molaire,
• $\overline{M}=\sum _i{x}_i{M}_i$ est la masse molaire moyenne.

On obtient une expression multilinéaire explicite où les coefficients de diffusion multicomposants $D_{ij}$ sont à calculer par inversion du système. Dans un calcul de mécanique des fluides cette résolution peut être pénalisante et on utilise diverses approximations permettent d'écrire le système plus simplement :

${𝐉}_i\simeq -{\rho }_i{D}_i\mathrm{\nabla }{c}_i$

où $D_i=\frac{1-w_i}{\sum _{k\ne i}\frac{w_k}{{𝒟}_{ik}}}$

$w_i$ est une pondération que l'on peut prendre égale à $x_i$ ou $c_i$.

Généralisation dans le cas des gaz

Coefficients de diffusion thermique pour l'air à pression normale.

La méthode de Chapman-Enskog permet de généraliser ce système pour un gaz avec la prise en compte des gradients de pression et de température

$\sum _{j\ne i}\frac{x_i{x}_j}{\rho {𝒟}_{ij}}(\frac{{𝐉}_j}{c_j}-\frac{{𝐉}_i}{c_i})=\mathrm{\nabla }{x}_i+(x_i-c_i)\mathrm{\nabla }\ln p+k_i^T\mathrm{\nabla }\ln T$

où

• $p$ est la pression (unité SI :Modèle:Unité),
• $T$ est la température (unité SI :Modèle:Unité),
• $k_i^T$ est coefficient de diffusion thermique multicomposant (sans dimension).

La diffusion par gradient thermique constitue l'effet Soret.

On trouve parfois cette expression écrite sous une forme équivalente

$\sum _{k\ne i}\frac{x_i{x}_k}{\rho {𝒟}_{ik}}(\frac{{𝐉}_k+{𝒟}_k^T\mathrm{\nabla }(\ln T)}{c_k}-\frac{{𝐉}_i+{𝒟}_i^T\mathrm{\nabla }(\ln T)}{c_i})=\mathrm{\nabla }{x}_i+(x_i-c_i)\mathrm{\nabla }(\ln p)$

où ${𝒟}_i^T$ est nommé coefficient de diffusion thermique (unité SI Modèle:Unité : ce n'est pas un coefficient de diffusion).

La relation entre ces deux systèmes d'équations se fait par la relation

$k_i^T=\sum _{j\ne i}\frac{x_i{x}_j}{\rho {𝒟}_{ij}}(\frac{{𝒟}_j^T}{c_j}-\frac{{𝒟}_i^T}{c_i})$

Le rang du système étant N-1 on a :

$\sum _i{𝒟}_i^T=0$

$\sum _i{k}_i^T=0$

Pour un milieu binaire l'équation de Stefan-Maxwell s'écrit :

${𝐉}_2-{𝐉}_1=\frac{\rho {𝒟}_{12}}{x_1{x}_2}(\mathrm{\nabla }{x}_1+k_1^T\mathrm{\nabla }\ln T)$

La solution est :

${𝐉}_1=-{𝐉}_2=-\frac{\rho {𝒟}_{12}}{2x_1{x}_2}(\mathrm{\nabla }{x}_1+k_1^T\mathrm{\nabla }\ln T)$

À l'équilibre (pas de diffusion), il peut exister un gradient de concentration lié au gradient de température :

$\mathrm{\nabla }{x}_2=-\mathrm{\nabla }{x}_1=k_1^T\mathrm{\nabla }\ln T$

Pour les gaz les coefficients de diffusion binaires varient approximativement comme $\frac{T^{\frac{3}{2}}}{p}$. ${𝒟}_i^T$ et $k_i^T$ ont un signe quelconque et une variation en température assez erratique comme indiqué par les courbes jointes.

Extension à un milieu chargé

Le système décrivant la diffusion peut être étendu au cas d'un gaz contenant en faible quantité des particules portant une charge électrique :

$\sum _{j\ne i}\frac{x_i{x}_j}{\rho {𝒟}_{ij}}(\frac{{𝐉}_j}{c_j}-\frac{{𝐉}_i}{c_i})=\mathrm{\nabla }{x}_i+(x_i-c_i)\mathrm{\nabla }\log p+k_i^T\mathrm{\nabla }\log T-\frac{1}{p}(Q_i-c_{i}Q)𝐄$

où

• $Q_i=n_i{Z}_i{q}_e$ est la densité de charge pour l'espèce i,
• $Z_i$ le nombre de charges de la particule i, $Z_i=-1$ pour l'électron,
• $q_e$ la charge de l'électron,
• $Q=\sum _i{x}_i{Q}_i$ la densité de charge totale,
• $𝐄$ le champ électrique.

La simplification de cette expression dans un milieu quasineutre conduit à l'approximation de la diffusion ambipolaire.

Références

Modèle:Portail