Détermination Expérimentale des Coefficients du Viriel d'un Gaz Réel
📝 Présentation du Système et Contexte Industriel
Dans le cadre stratégique du dimensionnement de cuves de stockage sous très haute pression pour l'industrie aérospatiale et la cryogénie avancée, les ingénieurs thermiciens se heurtent systématiquement aux limites physiques fondamentales. En effet, le modèle du gaz parfait échoue lourdement à prédire le comportement volumétrique des fluides dès que la pression dépasse quelques dizaines de bars.
Pour garantir la sécurité absolue et l'optimisation massique des réservoirs contenant de l'Argon pur (\(\text{Ar}\)), il est impératif de basculer vers une modélisation rigoureuse. Cette approche doit obligatoirement intégrer les interactions intermoléculaires à courte et moyenne distance, en employant le célèbre développement du Viriel.
Définition du Système (\(\Sigma\)) : Le système thermodynamique étudié est un système fermé constitué d'une quantité de matière fixe de gaz Argon. Ce fluide est confiné dans l'enceinte cylindrique de volume variable d'un banc d'essai piézométrique expérimental. La frontière du système est donc strictement délimitée par les parois internes indéformables du cylindre et la face inférieure mobile du piston.
Milieu Extérieur : L'environnement en interaction avec notre système est double. Il comprend, d'une part, un bain thermostatique de très haute précision agissant comme une source de chaleur idéale (capacité thermique infinie). D'autre part, il inclut le dispositif mécanique de la presse hydraulique, agissant comme source de travail des forces de pression.
Nature de la Transformation et États : Le gaz Argon subit une compression isotherme quasi-statique. Nous étudions son évolution depuis un État Initial (\(1\)) (caractérisé par une haute pression modérée) vers un État Final (\(2\)) (caractérisé par une très haute pression). La transformation étant quasi-statique, le système traverse une infinité d'états d'équilibre, garantissant l'homogénéité de la température et de la pression à chaque instant.
Exploiter rigoureusement les données d'état expérimentales (pression et volume molaire) mesurées sur une isotherme unique. Le but est d'isoler mathématiquement et d'évaluer les coefficients du viriel \(B(T)\) (traduisant les interactions à deux corps) et \(C(T)\) (traduisant les interactions à trois corps) pour l'Argon gazeux.
"Convention thermodynamique égoïste respectée : tout transfert reçu par le système est compté positivement. Par ailleurs, la difficulté réside ici dans le refus formel d'appliquer l'équation des gaz parfaits. La modélisation Gaz Réel implique une déviation forte de l'idéalité."
Afin d'alimenter notre modèle mathématique et de procéder aux bilans, l'équipe d'ingénierie a plongé le cylindre dans un thermostat. Le système thermodynamique est ainsi fermement verrouillé à une température d'étude spécifique et invariante.
Dans ces conditions de thermostatisation parfaite, la transformation a été stoppée à deux reprises pour relever les variables d'état d'équilibre. Ces relevés associent, pour chaque palier de la compression, la pression absolue exercée sur le fluide et le volume molaire résultant. Ce sont ces données brutes de laboratoire qui constitueront le socle exclusif de notre résolution.
📚 Modèles Mathématiques Cibles
Équation d'État (Gaz Réel) Développement en série du Viriel tronqué à l'ordre 2| PARAMÈTRE DE CONTRÔLE ISOTHERME | |
| Température du bain thermostatique | \( T = 300.0 \text{ K} \) |
| POINT DE MESURE N°1 (Haute Pression Modérée) | |
| Pression mesurée | \( P_1 = 4.863 \times 10^6 \text{ Pa} \) |
| Volume molaire mesuré | \( V_{\text{m},1} = 5.000 \times 10^{-4} \text{ m}^{3} \cdot \text{mol}^{-1} \) |
| POINT DE MESURE N°2 (Très Haute Pression) | |
| Pression mesurée | \( P_2 = 11.910 \times 10^6 \text{ Pa} \) |
| Volume molaire mesuré | \( V_{\text{m},2} = 2.000 \times 10^{-4} \text{ m}^{3} \cdot \text{mol}^{-1} \) |
| CONSTANTES UNIVERSELLES | |
| Constante des gaz parfaits | \( R = 8.314 \text{ J} \cdot \text{K}^{-1} \cdot \text{mol}^{-1} \) |
E. Méthodologie de Résolution (Cadre Thermo)
En thermodynamique des fluides réels, la rigueur méthodologique est primordiale pour éviter les impasses mathématiques. Nous allons structurer la résolution de ce système d'équations non-linéaires selon une démarche d'isolation séquentielle.
Étape 1 : Évaluation Quantitative de la Déviation Idéale
Calcul rigoureux du Facteur de compressibilité \(Z\) pour chaque état expérimental afin de statuer sur la nature dominante des forces intermoléculaires existantes dans l'Argon.
Étape 2 : Formulation et Linéarisation du Modèle
Écriture du développement du Viriel en volume molaire. Création d'une variable auxiliaire \(Y\) permettant de linéariser le système à deux inconnues \(B(T)\) et \(C(T)\).
Étape 3 : Résolution Analytique du Système
Manipulation algébrique matricielle ou par substitution pour isoler algébriquement puis numériquement le coefficient du troisième ordre \(C(T)\), puis le second \(B(T)\).
Étape 4 : Interprétation Physique et Validation
Analyse du signe des coefficients obtenus. Mise en relation profonde avec la théorie du covolume de Van der Waals et l'attraction dipolaire à moyenne distance de l'Argon.
Détermination Expérimentale des Coefficients du Viriel d'un Gaz Réel
🎯 Objectif Scientifique
Avant même de manipuler le développement asymptotique du viriel, il est rigoureusement indispensable de quantifier la déviation à l'idéalité de notre gaz d'Argon pour chaque point de mesure.
Calculer le facteur de compressibilité permet de juger immédiatement l'erreur monumentale qu'aurait induite l'utilisation naïve du modèle du gaz parfait à ces niveaux de pression.
La thermodynamique macroscopique postule que tout fluide pur à l'équilibre voit son état entièrement défini par une relation liant \(P\), \(V\) et \(T\). Le principe méthodologique fondamental appliqué ici consiste à utiliser le modèle d'état idéal comme étalon de référence absolu, afin de jauger mathématiquement l'amplitude des interactions intermoléculaires agissant dans le fluide réel.
Dans un système gazeux confiné soumis à une très haute pression, le volume libre offert aux déplacements particulaires chute drastiquement. Par conséquent, les molécules sont contraintes de se rapprocher considérablement.
Mon rôle d'analyste est d'observer si ce rapprochement forcé favorise préférentiellement les forces attractives dipolaires de Van der Waals (ce qui tendrait à diminuer la pression relative) ou la répulsion stérique des nuages électroniques (ce qui l'augmenterait brutalement).
Le facteur de compressibilité analytique \(Z\) est une fonction d'état adimensionnelle. Elle est formellement définie comme le rapport exact entre le volume molaire réel du gaz \(V_{\text{m, réel}}\) et le volume molaire qu'il occuperait s'il était rigoureusement idéal (\(V_{\text{m, parfait}} = \frac{R T}{P}\)) aux mêmes conditions de température et de pression.
Par définition, un gaz parfait possède strictement \(Z = 1\), quel que soit l'état d'équilibre du système.
📋 Paramètres de l'étape
| Paramètre | Valeur pour le Point 1 | Valeur pour le Point 2 |
|---|---|---|
| Pression mesurée \(P\) | \( 4.863 \times 10^6 \text{ Pa} \) | \( 11.910 \times 10^6 \text{ Pa} \) |
| Volume molaire \(V_{\text{m}}\) | \( 5.000 \times 10^{-4} \text{ m}^{3} \cdot \text{mol}^{-1} \) | \( 2.000 \times 10^{-4} \text{ m}^{3} \cdot \text{mol}^{-1} \) |
Assurez-vous rigoureusement que le terme d'énergie thermique \(R T\) est bien calculé en Unités du Système International (soit en Joules par mole). Ainsi, le numérateur (des Pascals multipliés par des \(\text{m}^{3} \cdot \text{mol}^{-1}\)) sera parfaitement homogène à des Joules par mole, garantissant que \(Z\) restera une grandeur strictement adimensionnelle.
📝 Exécution du Calcul Numérique Séquentiel
L'application de notre équation de définition à chacun des deux points piézométriques expérimentaux exige au préalable l'évaluation de la constante d'énergie thermique de l'isotherme : \(R T = 8.314 \times 300.0 = 2494.2 \text{ J} \cdot \text{mol}^{-1}\).
Calcul N°1 : Déviation pour l'État 1En injectant le couple de pression modérée et le volume associé, nous isolons la valeur adimensionnelle caractéristique du premier équilibre :
Le résultat analytique du premier point indique que \(Z_1\) est strictement inférieur à \(1\). Cela démontre de manière univoque que le gaz d'Argon est "plus facile" à compresser qu'un gaz idéal à cette pression spécifique de \(48 \text{ bars}\).
Calcul N°2 : Déviation pour l'État 2Procédons de manière parfaitement symétrique pour le second point, situé dans un régime de pression nettement plus oppressant :
Le facteur \(Z_2\) chute encore plus bas sur l'axe des ordonnées ! La déviation à l'idéalité s'amplifie considérablement. L'attraction intermoléculaire mutuelle des atomes d'Argon soulage la pression perçue par les parois par rapport à la prédiction d'un modèle non-interagissant.
Les deux valeurs macroscopiques calculées étant strictement inférieures à l'unité (\(Z < 1\)), nous actons thermodynamiquement que le fluide réel est dominé par des forces cohésives attractives sur cette plage isotherme.
La logique est respectée : l'intensification de la pression (de \(4.8\) à \(11.9 \text{ MPa}\)) force la proximité moléculaire, exacerbant ainsi l'effet réducteur de pression de l'attraction dipolaire.
L'erreur fatale très classique en examen consiste à confondre le volume molaire \(V_{\text{m}}\) (en \(\text{m}^{3} \cdot \text{mol}^{-1}\)) avec le volume massique global \(v\) (en \(\text{m}^{3} \cdot \text{kg}^{-1}\)). L'utilisation aveugle de \(v\) imposerait de diviser par la masse molaire \(M\) de l'Argon et d'utiliser la constante spécifique \(r = R/M\). La vérification de l'homogénéité dimensionnelle doit constituer votre boussole absolue.
🎯 Objectif Scientifique
Il est désormais crucial de relier mathématiquement la déviation expérimentale calculée à l'équation théorique microscopique du Viriel.
Le but exclusif de cette section est de transformer un modèle asymptotique complexe en un système linéaire élémentaire de deux équations à deux inconnues. Cela rendra la solution soluble algébriquement pour extraire les termes \(B\) et \(C\) sans approximation hasardeuse.
Le principe sous-jacent à l'équation du Viriel repose sur la mécanique statistique canonique. Il postule que les écarts à l'idéalité macroscopique peuvent être rigoureusement décomposés en une somme infinie d'interactions microscopiques classées par grappes de particules (paires, triplets, etc.). Chaque coefficient traduit la contribution thermodynamique d'un de ces clusters de collision.
Le développement analytique canonique se présente sous la forme d'une série de puissances de la densité molaire (\(\frac{1}{V_{\text{m}}}\)), ce qui est structurellement non-linéaire.
Étant donné que nous recherchons précisément deux paramètres constants, mon astuce d'ingénieur consiste à regrouper l'intégralité des termes expérimentalement connus du côté gauche de l'égalité. Ainsi, je conserverai une expression polynomiale purifiée des inconnues à droite, évitant un inextricable enchevêtrement algébrique.
Fondée par Heike Kamerlingh Onnes, l'équation d'état du Viriel est l'unique modèle d'état disposant d'une justification théorique exacte. Le second coefficient du viriel \(B(T)\) quantifie thermodynamiquement l'interaction stricte par paires de molécules (attraction/répulsion à deux corps).
Le troisième coefficient du viriel \(C(T)\) évalue l'impact beaucoup plus rare des collisions simultanées impliquant trois molécules. Cette expansion est formellement infinie, mais judicieusement tronquée à haute pression.
📋 Paramètres de l'étape
| Paramètre Analytique | Valeur pour le Point 1 | Valeur pour le Point 2 |
|---|---|---|
| Facteur de Compressibilité \(Z\) | \( 0.97486 \) | \( 0.95501 \) |
| Volume molaire \(V_{\text{m}}\) | \( 5.000 \times 10^{-4} \text{ m}^{3} \cdot \text{mol}^{-1} \) | \( 2.000 \times 10^{-4} \text{ m}^{3} \cdot \text{mol}^{-1} \) |
Puisque la température est fermement verrouillée par la régulation de notre bain thermostatique expérimental (\(T = 300 \text{ K}\)), les fonctions variables \(B(T)\) et \(C(T)\) se figent instantanément en constantes scalaires pures pour la totalité de ce problème.
Afin d'alléger drastiquement notre formalisme algébrique, nous les noterons sobrement \(B\) et \(C\) dans la suite de l'étude.
📝 Exécution de la Linéarisation Numérique
La manipulation algébrique qui suit vise à transformer l'équation asymptotique en une équation de droite de type affine \(Y = a X + b\).
Pour y parvenir, nous devons d'abord soustraire l'unité (qui représente la part du gaz idéal) aux deux membres de l'équation.
Ensuite, pour forcer le coefficient \(B\) à devenir une simple ordonnée à l'origine isolée, nous allons multiplier l'intégralité de l'expression par le volume molaire \(V_{\text{m}}\). Cette manœuvre génère notre nouvelle variable macroscopique \(Y\).
Calcul N°1 : Création de la variable intermédiaire \(Y\)Définissons formellement la grandeur auxiliaire continue \(Y\) en appliquant les opérations algébriques décrites :
Nous devons impérativement quantifier numériquement \(Y_1\) et \(Y_2\) en exploitant la matrice des données thermodynamiques issues de la phase 1.
Le système complexe asymptotique du viriel est désormais méticuleusement converti en un système affine \(Y_i = a X_i + b\) d'une robustesse absolue. La grandeur auxiliaire de calcul \(Y\) possède par construction la stricte même dimensionnalité qu'un volume molaire classique.
Vérifions le postulat de linéarisation : \(Z\) étant un scalaire sans dimension, l'expression \((Z-1)\) l'est également. La multiplication par \(V_{\text{m}}\) garantit formellement que \(Y_1\) et \(Y_2\) s'expriment intrinsèquement en \(\text{m}^{3} \cdot \text{mol}^{-1}\).
Ce verrouillage dimensionnel valide la viabilité de l'équation mathématique linéarisée avant son traitement matriciel.
Une hécatombe classique lors des épreuves de concours réside dans l'arrondi prématuré du facteur \(Z\) à deux malheureux chiffres significatifs. La soustraction opérée par \((Z-1)\) est une opération très sensible : toute troncature décimale excessive ici anéantira instantanément la précision de l'extraction de \(B\) et \(C\). La conservation d'au moins cinq chiffres significatifs de l'étape \(1\) est impérieuse.
🎯 Objectif Scientifique
L'apothéose de notre démonstration thermodynamique se déploie à cette étape précise. L'objectif mathématique implacable est de manipuler le système linéaire \(2 \times 2\) mis en évidence précédemment afin d'extraire les coefficients du viriel sous leur forme pure.
Nous dériverons d'abord formellement le terme de la pseudo-pente \(C(T)\), puis l'ordonnée dynamique à l'origine \(B(T)\).
Le principe mathématique central est fondé sur l'indépendance linéaire stricte des deux équations d'état acquises sous deux niveaux de pression radicalement distincts.
Cette décorrélation expérimentale assure que le déterminant du système de substitution est non nul, garantissant ainsi l'existence et l'unicité des constantes d'interaction recherchées.
Pour éliminer chirurgicalement la constante inconnue \(B\) de notre champ visuel, la soustraction pure et simple des deux équations d'état linéarisées s'impose comme l'opération la plus élégante et la moins risquée.
L'opération différentielle \((Y_2 - Y_1)\) dissipera le coefficient des interactions exclusives à deux corps, mettant spectaculairement en lumière le comportement marginal à trois particules, entièrement absorbé par le terme \(C\).
La méthode par combinaison linéaire est largement privilégiée en ingénierie face aux substitutions itératives. Elle amoindrit la propagation vicieuse des incertitudes expérimentales lors du traitement des données thermodynamiques, car elle centralise les rapports différentiels en une seule équation mère.
📋 Synthèse des Entrants pour le Calcul Final
| Grandeur Théorique | Valeur Préparée |
|---|---|
| Grandeur modifiée \(Y_1\) | \( -12.57 \times 10^{-6} \text{ m}^{3} \cdot \text{mol}^{-1} \) |
| Grandeur modifiée \(Y_2\) | \( -8.998 \times 10^{-6} \text{ m}^{3} \cdot \text{mol}^{-1} \) |
| Densité molaire inverse \(1/V_{\text{m},1}\) | \( 2000 \text{ mol} \cdot \text{m}^{-3} \) |
| Densité molaire inverse \(1/V_{\text{m},2}\) | \( 5000 \text{ mol} \cdot \text{m}^{-3} \) |
Prenez un soin obsessionnel lors de la soustraction des termes \(Y\). Les deux valeurs étant négatives, l'opération \((Y_2 - Y_1)\) fera inévitablement intervenir une collision de signes moins, générant un terme de soustraction positif.
Calculez les dénominateurs inversés séparément pour éviter une erreur de parenthèse fatale sur la machine de calcul.
📝 Exécution Numérique : Isolation du coefficient de choc tertiaire
Nous initions la substitution finale en prêtant une attention clinique à l'ordre de grandeur des puissances de dix négatives régissant les dénominateurs.
Calcul N°1 : Évaluation du 3ème Coefficient \(C(300 \text{ K})\)L'obtention d'une valeur strictement positive pour ce coefficient est capitale. Elle traduit une force thermodynamique de répulsion inévitable se déclenchant lors de rarissimes collisions statistiques à trois particules, imposée par l'impénétrabilité quantique des cœurs électroniques massifs de l'Argon.
📝 Exécution Numérique : Isolation du coefficient binaire
Calcul N°2 : Évaluation du 2nd Coefficient \(B(300 \text{ K})\)En rapatriant la valeur structurelle de \(C\) fraichement validée vers l'équation de la grandeur du premier palier expérimental, nous déverrouillons enfin \(B\) :
Le signe fermement négatif de la constante \(B\) démontre, sans la moindre équivoque théorique, que la composante dipolaire d'attraction globale (force quantique de dispersion de London) surpasse totalement l'effet de volume propre individuel (le covolume stérique) à la température tempérée du laboratoire de \(300 \text{ K}\).
Ces ordres de grandeur nanométriques corroborent magnifiquement les tables d'état officielles du NIST (National Institute of Standards and Technology) pour le fluide Argon. À titre purement prospectif, si notre thermostat était paramétré bien au-dessus de la Température de Boyle (soit \(T \approx 410 \text{ K}\) pour l'Ar), la force cohésive s'effondrerait.
La valeur macroscopique de \(B\) s'annulerait mathématiquement puis deviendrait positive, signalant la domination écrasante du comportement de type "sphères dures répulsives". L'analyse est donc un triomphe.
L'erreur structurelle redoutée concerne l'analyse dimensionnelle du troisième ordre. Assurez-vous par une analyse expresse que le coefficient \(C\) possède impérativement les unités d'un volume molaire élevé au carré : \(\text{m}^{6} \cdot \text{mol}^{-2}\).
Toute dimension résiduelle de pression indiquerait l'oubli de la variable \((Z-1)\) lors de la linéarisation initiale du système d'état.
🎯 Objectif Scientifique
Maintenant que l'équation d'état du fluide réel est rigoureusement calibrée avec ses coefficients du Viriel, nous devons accomplir l'ultime tâche de l'ingénieur thermicien.
L'objectif est de calculer l'énergie mécanique réelle (le travail \(W\)) requise par le compresseur pour faire passer l'Argon de l'État 1 à l'État 2, et mettre en évidence le gain ou la perte par rapport au modèle idéal.
Pour une transformation de compression quasi-statique (où la pression interne du gaz s'équilibre avec la force extérieure à chaque instant infinitésimal), le travail élémentaire des forces de pression s'écrit formellement \(\delta W = -P \text{d}V_{\text{m}}\).
Le signe négatif traduit la convention égoïste : une diminution de volume (\(\text{d}V_{\text{m}} < 0\)) entraîne un travail reçu positif.
La difficulté algébrique réside dans le fait que la Pression \(P\) n'est pas constante. Elle augmente violemment au fil de la compression. Nous devons impérativement substituer la variable \(P\) par son expression mathématique tirée du développement du Viriel, puis intégrer ce polynôme complexe entre les bornes volumiques initiales et finales.
Le fait que la transformation soit isotherme garantit que \(T\), \(B(T)\) et \(C(T)\) peuvent être sortis de l'intégrale comme de simples constantes scalaires, ce qui simplifie grandement l'intégration.
📋 Synthèse des Entrants pour l'Intégration
| Grandeur Thermodynamique Requise | Valeur Préparée |
|---|---|
| Énergie thermique \(R \cdot T\) | \( 2494.2 \text{ J} \cdot \text{mol}^{-1} \) |
| Volume molaire initial (Borne inférieure) \(V_{\text{m},1}\) | \( 5.000 \times 10^{-4} \text{ m}^{3} \cdot \text{mol}^{-1} \) |
| Volume molaire final (Borne supérieure) \(V_{\text{m},2}\) | \( 2.000 \times 10^{-4} \text{ m}^{3} \cdot \text{mol}^{-1} \) |
| Second coefficient du Viriel \(B\) | \( -14.95 \times 10^{-6} \text{ m}^{3} \cdot \text{mol}^{-1} \) |
| Troisième coefficient du Viriel \(C\) | \( 1.191 \times 10^{-9} \text{ m}^{6} \cdot \text{mol}^{-2} \) |
Lors de l'application numérique, l'évaluation du terme d'interaction binaire implique le coefficient \(B\) qui est strictement négatif. Prenez garde à la violente collision de signes lors de la soustraction des bornes de l'intégrale : \(- B \cdot (1/V_{\text{m},2} - 1/V_{\text{m},1})\).
De plus, isolez systématiquement le calcul du travail idéal (\(W_{\text{id}}\)) avant de lui appliquer les corrections croisées du viriel. Cela permet d'isoler l'erreur de calcul si votre résultat final vous semble physiquement aberrant.
📝 Exécution Numérique du Bilan Énergétique
Pour saisir l'impact de nos coefficients, nous allons calculer le travail mécanique qu'aurait exigé un gaz parfait (\(W_{\text{id}}\)), puis lui ajouter les corrections soustractives induites par les forces intermoléculaires réelles.
Calcul N°1 : Évaluation du Travail Idéal de Référence (\(W_{\text{id}}\))Le premier terme du crochet correspond à l'énergie de compression d'un fluide idéal :
Développons la différence des crochets avec \(B = -14.95 \times 10^{-6}\) et \(C = 1.191 \times 10^{-9}\). Notez la double inversion des bornes générée par la soustraction analytique :
Le travail réel calculé est formellement positif, ce qui honore la convention thermodynamique (le système fermé de laboratoire reçoit de l'énergie de la presse hydraulique). Cependant, il est radicalement inférieur au travail idéal prévu.
L'inégalité fondamentale \(W_{\text{réel}} < W_{\text{id}}\) confirme magistralement notre conclusion d'état issue de l'Étape \(3\) : les forces cohésives attractives dominantes (traduites par le signe \(B < 0\)) "aident" physiquement le piston à comprimer le fluide.
Conséquence directe pour le dimensionnement : le compresseur cryogénique industriel dépensera moins d'énergie électrique que ce que prévoyait l'équation naïve des gaz parfaits (soit une économie substantielle de \(\approx 80 \text{ Joules}\) par mole traitée).
📄 La Copie Parfaite (Ce qu'il faut écrire)
Voici le résumé académique de la résolution selon le formalisme exigé en thermodynamique, tel qu'attendu lors d'un audit de conformité d'ingénierie gazeuse.
THERMO MÉTABOLIQUE
- Système \(\Sigma\) : Gaz Argon étudié sous cloche isotherme. Système considéré comme fermé et monophasique (vapeur sèche).
- Modèle Théorique : Rejet absolu du gaz idéal. Utilisation du Développement asymptotique du Viriel, borné au \(3^{\text{ème}}\) ordre. L'équation d'état imposée est \(Z = 1 + \frac{B}{V_{\text{m}}} + \frac{C}{V_{\text{m}}^2}\).
- Objectif : Résolution du système linéaire à l'aide des observations piézométriques pour calibrer empiriquement \(B\) et \(C\).
Posez votre question