Équation du viriel — Wikipédia
L'équation d'état du viriel est une équation d'état utilisée pour décrire le comportement des fluides.
Elle s'écrit le plus souvent comme l'expression de , le facteur de compressibilité, en fonction des puissances de , le volume molaire :
Le coefficient est appelé e coefficient du viriel. Les coefficients sont déterminés expérimentalement pour les fluides réels. En pratique on se limite à et .
Du point de vue théorique, le théorème du viriel (Clausius, 1865) ou bien la physique statistique conduisent à ce type d'équation, les coefficients dépendant du modèle d'interaction entre particules (par exemple le modèle des sphères dures, ou le potentiel de Lennard-Jones).
Les deux formes de l'équation du viriel
[modifier | modifier le code]Il existe deux formes d'équation du viriel[1],[2]. L'une est exprimée en fonction du volume molaire :
l'autre en fonction de la pression :
avec :
- et le e coefficient du viriel, selon la forme retenue,
- la pression,
- la constante universelle des gaz parfaits,
- la température,
- le volume molaire,
- le facteur de compressibilité.
Les coefficients du viriel et ne dépendent que de la température et des fractions molaires des diverses espèces chimiques présentes dans le mélange.
Dans la suite de l'article, il est fait référence à la première forme de l'équation du viriel, en fonction du volume molaire.
Pour passer d'une forme à une autre, on fait apparaitre le facteur de compressibilité dans la fonction de :
Puis on substitue la fonction de :
Après développement et réarrangement des termes par puissance de :
on identifie les divers coefficients :
- etc.
d'où l'on tire :
- etc.
et réciproquement :
- etc.
Notions de mécanique statistique
[modifier | modifier le code]En mécanique statistique, le fluide à décrire est considéré comme un système d'un grand nombre de particules (l’ensemble canonique), chacune des particules étant repérée en coordonnées polaires par son rayon vecteur. Chaque particule possède une certaine quantité de mouvement, , et l'on note la fonction de partition des micro-états des particules.
L’hamiltonien du système de particules s'écrit :
où :
- est le potentiel d'interaction pour les N particules ;
- est la température inverse et par convention .
De par l'indépendance des impulsions et des positions dans le hamiltonien, on peut réécrire la fonction de partition comme :
en posant :
la fonction se calcule aisément :
Par propriété de l'ensemble canonique, on a :
et :
On voit donc, en posant , que :
L'intégrale des configurations Z' décrit donc l'écart à l'idéalité du gaz réel. Pour différents potentiels d'interaction connus (sphères dures, Lennard-Jones, colloides, etc.), le problème est ramené au calcul de Z'.
Lien avec la fonction de corrélation de paire du fluide :
Il est alors possible de montrer que le développement a un lien avec la fonction de corrélation de paire , qui est la probabilité normalisée de trouver une deuxième particule à la distance de la première. Le résultat du développement est :
Cette équation a aussi été trouvée par le théorème du viriel, dans l'article pression cinétique.
Expression des coefficients B2 et B3
[modifier | modifier le code]En poussant les calculs, on peut exprimer les premiers coefficients.
On appelle fonction d'Ursell-Mayer f(r) :
alors :
puis :
Au-delà, les calculs deviennent très complexes et requièrent des techniques de théorie des graphes (cluster-diagrammes). Le calcul de B4 par Boltzmann en 1900 est considéré comme un vrai « tour de force » ; les calculs ne reprendront qu'en 1950.
Développement du viriel dans le cas des sphères dures
[modifier | modifier le code]On appelle gaz de sphères dures (hard spheres en anglais), N sphères de diamètre 2, s'agitant en tous sens, sans autre interaction que l'interdiction stricte de volume exclu. Ce modèle est souvent pris comme référence dans les calculs théoriques ; il est bien adapté aux simulations numériques ; et il représente assez bien l'expérience si on lui adjoint une perturbation représentant l'attraction entre sphères.
Appelons « compacité » la quantité , avec , le volume d'une sphère.
Dans le cas peu dense, le calcul donne :
Boltzmann a calculé théoriquement les premiers coefficients du développement en série :
Depuis, des calculs numériques (méthode de Monte-Carlo) ont donné :
Coefficient | Valeur estimée |
---|---|
28,225 | |
39,82 | |
53,34 | |
68,54 | |
85,8 | |
105,8(+/- 0,5) |
Or Carnahan et Starling[3] (1969) ont remarqué la parenté avec la série entière de coefficient générique . Ceci montre que est une bonne (simple) « fonction-candidat » pour un gaz de sphères dures ; mais elle présente une asymptote pour x = 1, alors que, évidemment, on ne saurait dépasser la compacité maximale : cela limite sa pertinence. Néanmoins, la simplicité de cette équation d'état permet de l'utiliser aisément : l'énergie libre d'excès est . On peut aussi comparer avec la théorie de Percus-Yevick[4]. On peut aussi en déduire une théorie de l'équation d'état pour mélanges binaires.
On sait mieux faire : avec les ordinateurs actuels, on peut trouver les premiers termes de ce développement dit du « viriel » d'un gaz de sphères dures, et on peut approcher l'expression par un approximant de Padé, ce qui donne une meilleure compréhension de l'asymptote, qui n'est évidemment pas x=1. En effet, la compacité maximale est :
Par ailleurs, les simulations numériques sur l'ordinateur Argon révèlent une irrégularité vers 50 %.
Transition de phase liquide-gaz
Si on ajoute une pression interne -a·x2, on retrouve que les sphères peuvent se condenser en une phase liquide. En supprimant le terme en x3 au numérateur de l'équation de Carnahan-Starling, l'équation en V(P,T) est une équation cubique, comme l'équation d'état de van der Waals, et conduit au même type de résultats : existence d'un point critique, et d'une pression de vapeur saturante . La précision est un peu meilleure, sans plus. La théorie du point critique relève nécessairement des méthodes de renormalisation, initiées par D. Parola dans ce cas précis.
Développement du viriel dans le cas du potentiel de Lennard-Jones
[modifier | modifier le code]Le potentiel de Lennard-Jones est :
Dans le cas de faible densité, PV/NkT = 1 + B(T) · N/V, et la fonction B(T) se calcule numériquement en unités réduites (kT/Eo ; r/d) via la formule du viriel :
avec f(r) la fonction d'Ursell-Mayer :
L'allure de B(T) est grossièrement B(T) = b - a/RT, et s'annule donc pour RT = a/b = 3,2 (température de Boyle), mais mieux encore, B(T) possède un léger maximum pour RT ~ 30 (température de Joule). Par ajustement entre cette courbe et la valeur expérimentale de B(T), on en tire les valeurs de Eo et d (Hirschfelder, 1954).
Références
[modifier | modifier le code]- Pierre Infelta et Michael Graetzel, Thermodynamique, Principes et applications, éd. Universal-Publishers, , 488 p. (ISBN 978-1-58112-995-3, lire en ligne), p. 114.
- Jean Vidal, Thermodynamique : application au génie chimique et à l'industrie pétrolière, Paris, Éditions Technip, coll. « Publications de l'Institut français du pétrole. », , 500 p. (ISBN 978-2-7108-0715-5, OCLC 300489419, lire en ligne), p. 96 et 105.
- (en) N.F. Carnahan et K.E. Starling, « Equation of State for Nonattracting Rigid Spheres », Journal of Chemical Physics, no 51, , p. 635-636 (DOI 10.1063/1.1672048, Bibcode 1969JChPh..51..635C, lire en ligne).
- (en) Jerome K. Percus et George J. Yevick, « Analysis of Classical Statistical Mechanics by Means of Collective Coordinates », Phys. Rev., vol. 110, no 1, (DOI 10.1103/PhysRev.110.1, Bibcode 1958PhRv..110....1P, lire en ligne).