Intégrales de collision

Les intégrales de collision sont des quantités qui interviennent dans le calcul des propriétés de transport des gaz (coefficients de diffusion, viscosités, conductivité) et qui font le lien avec le niveau microscopique décrit par le potentiel d'interaction entre deux particules.

Collision élastique de deux particules

Schéma d'une interaction moléculaire élastique dans le système lié au barycentre.

Les vitesses avant interaction sont Modèle:Math et Modèle:Math dans un référentiel galiléen. Les indices représentent indifféremment une même espèce ou deux espèces différentes. Ces vitesses valent Modèle:Math et Modèle:Math après interaction. On se place dans un système centré sur le barycentre qui a une vitesse constante du fait de la conservation de la quantité de mouvement. Dans ce système qui est donc galiléen la vitesse initiale de la particule Modèle:Mvar est la vitesse relative Modèle:Math. Par symétrie on peut affirmer que la trajectoire sera contenue dans le plan contenant l'origine et Modèle:Math. On choisit un repère tel que $\mathrm{\Omega }=\frac{{𝐠}_{ij}}{||{𝐠}_{ij}||}=[1,0,0]$ (voir figure). Dans ce repère la déviation est Modèle:Mvar, fonction du paramètre d'impact Modèle:Mvar, de la vitesse relative Modèle:Math et du potentiel d'interaction que l'on suppose ne dépendant que de la distance en les deux particules en interaction. Si cette hypothèse est rigoureuse pour l'interaction entre deux atomes, on peut la considérer utilisable pour deux molécules : le potentiel est alors un potentiel moyen statistique.

La direction de sortie d'interaction est définie par ${\mathrm{\Omega }}^{\prime }=\frac{{{𝐠}^{\prime }}_{ij}}{||{{𝐠}^{\prime }}_{ij}||}$. On peut calculer les vitesses finales à partir des considérations suivantes :

• la conservation de la quantité de mouvement dans l'interaction implique

${{𝐯}_i}^{\prime }+{{𝐯}_j}^{\prime }={𝐯}_i+{𝐯}_j$

• la vitesse relative Modèle:Math a un module constant du fait de la conservation de l'énergie, donc Modèle:Math ou

$|{{𝐯}_i}^{\prime }-{{𝐯}_j}^{\prime }|=|{𝐯}_i-{𝐯}_j|$

Les vitesses après l'interaction sont donc

${{𝐯}_i}^{\prime }(b,g)={𝐯}_i-({𝐠}_{𝐢𝐣}\cdot {\mathrm{\Omega }}^{\prime }){\mathrm{\Omega }}^{\prime }$

${{𝐯}_j}^{\prime }(b,g)={𝐯}_j+({𝐠}_{𝐢𝐣}\cdot {\mathrm{\Omega }}^{\prime }){\mathrm{\Omega }}^{\prime }$

De plus la conservation du moment cinétique au cours de l'interaction conduit à Modèle:Mvar.

Le système décrivant la collision est réversible. Le théorème de Liouville permet donc d'écrire

$\mathrm{d}{{𝐯}_i}^{\prime }\mathrm{d}{{𝐯}_j}^{\prime }=\mathrm{d}{𝐯}_i\mathrm{d}{𝐯}_j$

Définition de l'intégrale de collision

La méthode de Chapman-Enskog permet, dans le cas d'un gaz, de retrouver les équations de Navier-Stokes comme solution à l'ordre un d'un développement asymptotique de l'équation de Boltzmann. Les quantités qui apparaissent font intervenir les intégrales suivantes, liées aux collisions des couples de particules Modèle:MathModèle:,Modèle:,Modèle:,.

${\mathrm{\Lambda }}_{ij}^{l,s}=\sqrt{\frac{2\pi kT}{m_{ij}}}\underset{0}{\overset{\mathrm{\infty }}{\int }}\underset{0}{\overset{\mathrm{\infty }}{\int }}{\mathrm{e}}^{-{\gamma }_{ij}^2}{\gamma }_{ij}^{2s+3}(1-{\cos }^l\theta )b\mathrm{d}b\mathrm{d}{\gamma }_{ij},\theta =\theta (b,g_{ij})$

où

$m_{ij}={(m_i^{-1}+m_j^{-1})}^{-1}$	est la masse réduite,
${\gamma }_{ij}=\sqrt{\frac{m_{ij}}{kT}}\frac{g_{ij}}{2}$	est la vitesse relative réduite.

La distance d'approche minimale pour le potentiel choisi Modèle:Math supposé de révolution (au moins en moyenne sur toutes les collisions) est :

$r_m=b{(1-\frac{V(r)}{\frac{1}{2}{m}_{ij}{g}_{ij}^2})}^{-\frac{1}{2}}$

L'angle de déviation est donné par :

$\theta =\pi -2\underset{r_m}{\overset{\mathrm{\infty }}{\int }}{(\frac{1}{r_m^2}-\frac{1}{r^2})}^{-\frac{1}{2}}\frac{\mathrm{d}r}{r^2}$

Le calcul est analytique pour un potentiel sphères dures de diamètres Modèle:Mvar :

${\left[{\mathrm{\Lambda }}_{ij}^{l,s}\right]}_{SD}=\sqrt{\frac{kT}{2\pi m_{ij}}}\frac{(s+1)!}{2}[1-\frac{1+(-1{)}^l}{2(1+l)}]\pi {\sigma }_{ij}^2$

Dans le cas général, par exemple avec un potentiel de Lennard-Jones très utilisé par le passé à cause de sa simplicité et de son réalisme, on définit des intégrales de collision réduites :

$\stackrel{i,j}{\bar{\mathrm{\Lambda }}}=\frac{{\mathrm{\Lambda }}^{i,j}}{{\left[{\mathrm{\Lambda }}_{ij}^{l,s}\right]}_{SD}}$

Par construction la valeur de ces intégrales est proche de l'unité. Elles sont tabulées pour divers potentielsModèle:,.

Propriétés de transport

Les quantités l et s sont relatives aux expressions des coefficients de transport dans lesquelles elles interviennent.

Propriété	collision (i, i)	collision (i, j≠i)
viscosité	${\mathrm{\Lambda }}^{2,2}$	${\mathrm{\Lambda }}^{2,2}$, ${\mathrm{\Lambda }}^{1,1}$
conductivité	${\mathrm{\Lambda }}^{2,2}$	${\mathrm{\Lambda }}^{2,2},{\mathrm{\Lambda }}^{1,s},s=1,2,3$
diffusion	${\mathrm{\Lambda }}^{1,1}$	${\mathrm{\Lambda }}^{1,1}$
diffusion thermique	---	${\mathrm{\Lambda }}^{2,2},{\mathrm{\Lambda }}^{1,s},s=1,2,3$

Ainsi on obtient pour un gaz homogène :

- le coefficient d'auto-diffusion	$𝒟=\frac{1}{\rho }\frac{3}{8}\frac{\sqrt{\pi mkT}}{\stackrel{1,1}{\bar{\mathrm{\Lambda }}}}$
- la viscosité dynamique	$\mu =\frac{5}{16}\frac{\sqrt{\pi mkT}}{\stackrel{2,2}{\bar{\mathrm{\Lambda }}}}$
- la conductivité thermique	$\lambda =\frac{5}{2}{C}_V\mu =\frac{25}{32}\frac{C_V\sqrt{\pi mkT}}{\stackrel{2,2}{\bar{\mathrm{\Lambda }}}}$

Voir aussi

Articles connexes

• Potentiel sphères dures
• Potentiel de Lennard-Jones
• Potentiel de Morse
• Potentiel de Buckingham