Étude Cryogénique : Liquéfaction de l'Air par le Procédé Claude
📝 Présentation du Système de Liquéfaction (Procédé Claude)
L'industrie cryogénique moderne repose sur des cycles thermodynamiques d'une complexité fascinante, conçus pour atteindre des températures proches du zéro absolu.
Historiquement, l'ingénieur Georges Claude a bouleversé ce domaine technique en concevant un procédé novateur de liquéfaction de l'air. Contrairement au simple procédé de Linde, qui se base uniquement sur une détente isenthalpique à travers une vanne, le procédé Claude intègre une détente avec production de travail extérieur au cœur même de l'installation.
D'un point de vue thermodynamique, cette ingénieuse modification permet d'abaisser drastiquement la température du fluide tout en récupérant une puissance mécanique précieuse.
Dans ce système ouvert en écoulement stationnaire, nous allons modéliser une section critique de l'installation fonctionnant en régime permanent. L'air y est rigoureusement assimilé à un gaz parfait diatomique, ce qui simplifiera l'application de la première loi de Joule.
Chronologiquement, le parcours du fluide s'articule en trois phases distinctes. Premièrement, le gaz traverse un compresseur idéal. Afin d'éviter une surchauffe destructrice, cette compression est supposée strictement isotherme (\( T_2 = T_1 \)) grâce à un refroidissement continu.
Ensuite, l'air à très haute pression pénètre dans un échangeur de chaleur isobare parfait, où il subit un refroidissement intense à pression constante (\( P_3 = P_2 \)).
Finalement, ce gaz comprimé et ultra-froid est injecté dans une turbine de détente adiabatique et réversible. Il s'y détend jusqu'à retrouver sa pression atmosphérique initiale (\( P_4 = P_1 \)). C'est ici que l'effondrement thermique se produit, générant un précieux travail de transvasement.
Votre mission d'ingénierie consiste à évaluer rigoureusement les transferts d'énergie de cette séquence cryogénique.
Concrètement, vous devrez déterminer l'énergie massique consommée par le compresseur, la puissance thermique colossale évacuée par l'échangeur, et prédire la température extrême atteinte en sortie de turbine.
En conclusion, il sera impératif de valider la faisabilité théorique du cycle en réalisant un bilan entropique complet pour quantifier l'irréversibilité du modèle.
"Par convention internationale, nous appliquerons la stricte convention thermodynamique égoïste.
Ainsi, toute forme d'énergie (Travail massique \(w\) ou Chaleur massique \(q\)) qui franchit la frontière vers l'intérieur du système est comptabilisée positivement (\(> 0\)). À l'inverse, toute énergie dissipée ou cédée par le fluide vers l'extérieur sera assortie d'un signe négatif (\(< 0\)).
L'énergie cinétique et potentielle macroscopiques sont négligées devant les variations d'enthalpie."
Afin de structurer rigoureusement notre démarche analytique et nos bilans massiques, il est crucial de sanctuariser les variables d'état d'entrée fixées par le cahier des charges industriel.
Dans cette section, nous recensons exhaustivement les pressions de consigne, les températures des sources, ainsi que les propriétés intrinsèques de la matière traitée.
📚 Cadre Théorique et Équations de Comportement
L'air sec, dans cette plage de pressions et de températures (avant la zone de liquéfaction finale), est mathématiquement assimilé à un gaz répondant aux lois idéales.
L'absence de dissociation moléculaire permet de considérer ses capacités thermiques comme parfaitement constantes.
| DÉFINITION DES NŒUDS THERMODYNAMIQUES (États 1, 2, 3 et 4) | |
| Pression atmosphérique d'admission (État 1) | \( P_1 = 1.0 \times 10^5 \text{ Pa} \) (\(1 \text{ bar}\)) |
| Température ambiante d'aspiration (État 1) | \( T_1 = 300 \text{ K} \) |
| Température de sortie du compresseur (Isotherme) | \( T_2 = T_1 = 300 \text{ K} \) |
| Pression de refoulement du compresseur (État 2) | \( P_2 = 40.0 \times 10^5 \text{ Pa} \) (\(40 \text{ bar}\)) |
| Pression dans l'échangeur (Isobare) | \( P_3 = P_2 = 40.0 \times 10^5 \text{ Pa} \) |
| Température cible en sortie d'échangeur (État 3) | \( T_3 = 160 \text{ K} \) |
| Pression d'échappement de la turbine (État 4) | \( P_4 = P_1 = 1.0 \times 10^5 \text{ Pa} \) |
| CONSTANTES INTRINSÈQUES DU FLUIDE (Air Diatomique) | |
| Capacité thermique massique à pression constante | \( c_p = 1005 \text{ J.kg}^{-1}\text{.K}^{-1} \) |
| Constante massique spécifique de l'air | \( r = 287 \text{ J.kg}^{-1}\text{.K}^{-1} \) |
| Coefficient polytropique isentropique (Indice de Laplace) | \( \gamma = \frac{c_p}{c_v} = 1.40 \) |
E. Méthodologie de Résolution (Cadre Thermo Industriel)
En thermodynamique des systèmes fluides ouverts, la rigueur méthodologique est le rempart contre les erreurs de signe.
Nous allons scinder notre étude du cycle de Claude en suivant précisément le parcours massique de l'air.
Étape 1 : Détermination de l'état final
Résolution analytique de l'état 4 en sortie de machine.
Étape 2 : Bilans d'Échange (Travail & Chaleur)
Calcul du travail de transvasement et des transferts thermiques échangés lors de la transformation complète.
Étape 3 : Entropie et Réversibilité
Utilisation du Second Principe pour lier les grandeurs d'état et valider théoriquement le cycle.
Étape 4 : Synthèse et Performances
Calcul final des grandeurs dérivées : rendement d'un cycle et taux de récupération de l'usine.
Liquéfaction de l'Air par le Procédé Claude
🎯 Objectif Scientifique
Avant d'évaluer les transferts d'énergie, nous devons impérativement déterminer les variables d'état manquantes aux nœuds du système.
Dans notre procédé, l'état final 4 (sortie de la turbine) est la clé de la liquéfaction. La température \(T_4\) est inconnue. L'objectif absolu de cette étape est de la calculer rigoureusement pour verrouiller la matrice des états thermodynamiques.
Dans la turbine (3 → 4), le gaz subit une détente extrêmement rapide de \(40 \text{ bar}\) à \(1 \text{ bar}\).
Puisqu'elle est modélisée comme parfaitement isolée et sans frottement interne, cette transformation est adiabatique et réversible, c'est-à-dire isentropique.
Nous pouvons donc lier directement la pression et la température via les lois des gaz parfaits pour trouver l'état final.
Pour un gaz parfait subissant une évolution isentropique, les variables d'état sont intrinsèquement couplées par le coefficient de Laplace \(\gamma\) (rapport des chaleurs spécifiques \(c_p/c_v\)).
La relation \(T^{\gamma} P^{1-\gamma} = \text{constante}\) démontre qu'une chute abrupte de pression engendre inévitablement un effondrement de la température.
📋 Paramètres de l'étape
| Paramètre | Valeur |
|---|---|
| Température d'admission turbine (\(T_3\)) | \( 160 \text{ K} \) |
| Pressions amont et aval | \( P_3 = 40 \text{ bar} \) et \( P_4 = 1 \text{ bar} \) |
| Exposant adiabatique de l'air | \( \frac{\gamma-1}{\gamma} = \frac{1.4-1}{1.4} \approx 0.2857 \) |
Soyez implacable sur la fraction de pression : la turbine détend le gaz. La pression finale (\(1 \text{ bar}\)) est au numérateur, et la pression initiale (\(40 \text{ bar}\)) au dénominateur.
Le ratio \((1/40)\) est inférieur à 1, ce qui garantit mathématiquement que la température s'effondrera !
📝 Calculs Détaillés
1. Évaluation Numérique de la Température finale (\(T_4\))Nous appliquons l'équation de Laplace avec l'exposant préalablement calculé pour éviter les erreurs de parenthèses sur la calculatrice.
Le résultat est spectaculaire. La température plonge à environ \(55.8 \text{ K}\) (soit \(-217 \text{ °C}\)), pulvérisant le seuil de liquéfaction de l'air. L'état final 4 est désormais complètement défini.
La valeur est parfaitement cohérente pour une installation cryogénique.
Une température finale largement inférieure à \(77 \text{ K}\) (point d'ébullition de l'azote à pression atmosphérique) confirme que nous sommes entrés dans la zone diphasique visée par l'invention de Georges Claude.
Oublier de convertir les températures en Kelvin dans la formule de Laplace !
La loi \(T^{\gamma} P^{1-\gamma} = \text{constante}\) n'a de sens physique absolu qu'avec l'échelle absolue des températures (Kelvin).
❓ Question Fréquente
Peut-on utiliser la loi de Charles ou de Gay-Lussac ici ?
Absolument pas, car ni le volume, ni la pression, ni la température ne sont constants lors du passage dans la turbine. Seule l'entropie l'est (transformation isentropique), d'où l'usage exclusif des relations de Laplace.
🎯 Objectif Scientifique
Cette étape constitue le cœur du bilan mécanique. Nous allons balayer chaque organe de la machine (Compresseur, Échangeur, Turbine) pour quantifier rigoureusement toutes les énergies massiques transférées.
L'objectif est d'obtenir le travail de transvasement et la chaleur échangée pour chaque transformation afin d'établir la matrice énergétique de l'usine.
Pour le compresseur (1→2) : Il est isotherme. Donc \(\Delta h = 0\), et \(q_{\text{comp}} = -w_{\text{comp}}\). Le calcul s'opère via l'intégrale \(+\int v dP\).
Pour l'échangeur (2→3) : Il est sans pièce mobile. Donc \(w_{\text{ech}} = 0\), et \(q_{\text{ech}}\) se déduit directement de \(\Delta h = c_p \Delta T\).
Pour la turbine (3→4) : Elle est adiabatique. Donc \(q_{\text{turb}} = 0\), et le travail fourni \(w_{\text{turb}}\) équivaut à la chute d'enthalpie \(\Delta h = c_p \Delta T\).
Dans une machine en écoulement continu (système ouvert), on utilise la fonction Enthalpie massique (\(h\)) au lieu de l'énergie interne (\(u\)).
Le bilan s'écrit \(\Delta h = w_{\text{tr}} + q\), où \(w_{\text{tr}}\) est le travail de transvasement (utile), évalué analytiquement par \(\int v dP\) pour les transformations réversibles.
📋 Données d'Entrée
| Paramètre | Valeur |
|---|---|
| Températures opératoires | \( T_1 = T_2 = 300 \text{ K}, \, T_3 = 160 \text{ K}, \, T_4 = 55.82 \text{ K} \) |
| Capacité et Constante de l'air | \( c_p = 1005 \text{ J.kg}^{-1}\text{.K}^{-1}, \, r = 287 \text{ J.kg}^{-1}\text{.K}^{-1} \) |
Gardez rigoureusement l'ordre des soustractions (\(T_{\text{finale}} - T_{\text{initiale}}\)).
Cela garantit que les signes thermodynamiques ("égoïstes") sortiront tout seuls sans avoir besoin de deviner si le système a perdu ou gagné de l'énergie.
📝 Calculs Détaillés
1. Compresseur IsothermeCalcul de la puissance mécanique exigée et de la chaleur rejetée à l'eau de refroidissement.
Calcul de la chaleur massivement évacuée par le gaz lors du sous-refroidissement.
Calcul de la puissance mécanique générée par l'arbre de transmission.
La convention égoïste est respectée. Le compresseur exige un apport colossal d'énergie (travail positif). La turbine restitue une fraction de cette puissance (travail négatif). Les chaleurs sont négatives, prouvant l'intense refroidissement requis.
Le travail de la turbine (\(104.7 \text{ kJ/kg}\)) représente près d'un tiers du travail requis par le compresseur (\(317.6 \text{ kJ/kg}\)).
C'est la signature de l'ingéniosité de Claude : cette énergie mécanique n'est pas perdue, elle peut être réinjectée sur l'arbre moteur pour économiser l'électricité !
Une erreur tragique très courante consiste à calculer le travail du compresseur en utilisant la formule des systèmes fermés \(-\int P dv\).
Dans une machine à flux continu, les travaux de poussée (refoulement et aspiration) sont colossaux, il faut impérativement utiliser \(+\int v dP\).
❓ Question Fréquente
Pourquoi l'échangeur ne fournit-il aucun travail ?
Parce qu'il ne contient aucune pièce mobile. Le fluide s'y écoule de manière tubulaire sans repousser de pales ni de piston. Mathématiquement, la différentielle du travail de transvasement \(dP\) est nulle (isobare), rendant le travail intégralement nul.
🎯 Objectif Scientifique
Il ne s'agit plus de calculer de l'énergie, mais de juger la cohérence de notre modélisation théorique.
L'objectif magistral est de quantifier la création d'entropie (\(s_{\text{cr}}\)) sur l'ensemble du cycle (de l'état 1 à 4) afin de prouver mathématiquement que les hypothèses de réversibilité absolue de l'énoncé ont été rigoureusement respectées.
Décortiquons le Second Principe : La variation totale d'entropie d'un système (\(\Delta s\)) se compose méticuleusement de deux flux : l'entropie échangée (\(s_{\text{ech}}\)) via les transferts thermiques, et l'entropie sombrement créée (\(s_{\text{cr}}\)) par les frictions internes.
Pour un modèle réversible comme exigé ici, l'Inégalité de Clausius stipule fermement que cette entropie créée doit être nulle (\(s_{\text{cr}} = 0\)).
Pour s'affranchir du chemin suivi par le fluide, on utilise l'identité fondamentale liant les variables d'état : \(dh = Tds + vdP\).
C'est la clé de voûte permettant de dériver la variation exacte \(\Delta s\) d'un gaz parfait en se contentant de connaître les pressions et températures extrêmes.
📋 Données d'Entrée
| Paramètre | Valeur |
|---|---|
| Températures terminales | \( T_1 = 300 \text{ K}, \, T_4 = 55.82 \text{ K} \) |
| Pressions terminales | \( P_1 = P_4 = 1 \text{ bar} \) |
Regardez attentivement les pressions : la turbine détend le gaz exactement à la pression atmosphérique initiale (\(P_4 = P_1\)).
De ce fait, le ratio de pression est de 1, et son logarithme naturel \(\ln(1)\) vaut 0. Le terme de pression s'annihile complètement !
📝 Calculs Détaillés
1. Calcul global de la variation de fonction d'état (\(\Delta s_{14}\))On applique uniquement le logarithme des températures extrêmes de l'installation.
Interprétation : L'entropie du fluide s'est écroulée de \(1690 \text{ J/K}\). L'air extrêmement froid et partiellement liquéfié est dans un état beaucoup plus "ordonné" que le gaz chaud turbulent de départ.
2. Calcul de l'Entropie Échangée globale (\(s_{\text{ech}}\))On somme les flux entropiques du compresseur et de l'échangeur isobare.
Note : La légère différence décimale (\(0.4\)) est uniquement due aux arrondis successifs sur les décimales de \(T_4\).
3. Bilan final : L'Entropie Créée (\(s_{\text{cr}}\))L'inégalité de Clausius est respectée dans son expression la plus pure.
L'obtention mathématique de la borne exacte \(0\) confirme notre hypothèse fondatrice : nous avons évalué des compresseurs, échangeurs et turbines purement réversibles.
Une confusion fatale sépare souvent les étudiants : \(\Delta s\) peut tout à fait être négatif !
Ce n'est pas \(\Delta s\) qui doit être obligatoirement positif, mais bien l'entropie CRÉÉE (\(s_{\text{cr}}\)). Ici, le fluide a rejeté de l'entropie (\(s_{\text{ech}} < 0\)).
❓ Question Fréquente
Que se passerait-il dans une usine réelle ?
Dans la réalité, les frottements mécaniques dans la turbine et les pertes de charges dans l'échangeur créeraient de l'entropie (\(s_{\text{cr}} > 0\)). La température finale \(T_4\) atteinte serait alors fatalement plus élevée (moins froide) que notre valeur idéale de \(55.8 \text{ K}\).
🎯 Objectif Scientifique
La dernière phase consiste à évaluer la pertinence industrielle globale de l'invention de Claude.
Il s'agit de calculer le travail net (\(w_{\text{net}}\)) consommé par l'ensemble de l'installation, ainsi que le taux de récupération d'énergie (\(\chi\)), qui définit l'efficacité du procédé par rapport à une technologie classique sans turbine.
Dans ce couplage complexe, l'effondrement thermique s'accompagne de la production d'un travail moteur par la turbine (\(w_{\text{turb}} < 0\)).
Si nous couplons mécaniquement cette turbine à l'arbre du compresseur principal, ce travail récupéré viendra en déduction du travail monstrueux requis à la compression (\(w_{\text{comp}} > 0\)).
Le bilan global de puissance mécanique s'obtient par l'addition stricte (algébrique) de tous les travaux de transvasement traversant la frontière macroscopique du système.
Le taux de récupération est simplement le ratio entre la valeur absolue de l'énergie produite et l'énergie initialement injectée.
📋 Données d'Entrée
| Paramètre | Valeur |
|---|---|
| Travail unitaire de compression (\(w_{12}\)) | \( +317612 \text{ J.kg}^{-1} \) |
| Travail unitaire de détente (\(w_{34}\)) | \( -104701 \text{ J.kg}^{-1} \) |
Assurez-vous de bien conserver le signe négatif du travail de la turbine lors de l'addition du travail net. C'est l'essence même du gain !
📝 Calculs Détaillés
1. Évaluation du Travail Net (\(w_{\text{net}}\))Opérations de déduction énergétique.
Interprétation : L'usine doit puiser environ \(212.9 \text{ kJ}\) d'électricité par kilogramme d'air traité sur le réseau national, au lieu des \(317.6 \text{ kJ}\) initialement requis par le seul compresseur.
2. Évaluation du Taux de Récupération (\(\chi\))Calcul du pourcentage géant d'économie d'énergie.
Interprétation : L'ajout de la turbine de Claude permet de récupérer et d'autofinancer massivement 33 % de l'effort de compression !
3. Résultats FinauxLe procédé Claude prouve mathématiquement sa suprématie. En plus de précipiter l'air dans la zone cryogénique, l'installation allège sa facture d'électricité d'un tiers absolu par le recyclage mécanique de la détente.
Calculer un "rendement" de machine motrice. Ce système reste une installation frigorifique (récepteur).
On parle donc de coefficient d'efficacité ou, dans le cas de l'arbre, de taux de récupération de l'énergie. Ne le confondez pas avec le théorème de Carnot réversible !
❓ Question Fréquente
Comment cette énergie est-elle physiquement récupérée ?
Dans les immenses installations industrielles (Air Liquide), l'arbre de la turbine est physiquement couplé, via un réducteur de vitesse, à l'arbre du compresseur électrique. La turbine l'aide littéralement à tourner.
📄 La Copie Parfaite (Rapport d'Ingénierie)
Voici le résumé académique de la résolution, présenté selon les standards d'excellence exigés dans la rédaction des rapports industriels ou de concours aux grandes écoles d'ingénieurs.
CRYOGÉNIE
- Système Thermodynamique : Masse unitaire de 1 kg d'air circulant en écoulement stationnaire continu. Le système est radicalement ouvert.
- Équation d'État : Le fluide est modélisé comme un Gaz Parfait idéal. Équation de comportement : \(Pv = rT\).
- Capacités Thermiques : L'air ne présente pas de dissociation moléculaire, les chaleurs massiques \(c_p\) et \(c_v\) sont rigoureusement constantes.
- Conventions de Signes : Approche égoïste. Travail et Chaleur comptabilisés positivement s'ils franchissent la frontière vers le fluide.
| Nœud d'Écoulement | Pression P (bar) | Température T (K) | Observations Techniques |
|---|---|---|---|
| 1. Entrée du procédé | 1.0 | 300 | Ambiant - Aspiration compresseur |
| 2. Fin compression isotherme | 40.0 | 300 | Pression maximale du cycle atteint |
| 3. Sortie échangeur isobare | 40.0 | 160 | Gaz fortement refroidi, densité crête |
| 4. Échappement turbine | 1.0 | 55.8 | Ultra-froid, passage en phase mixte (liquide/vapeur) |
Le cycle, exempt de frottements ou de chocs thermiques brutaux modélisés, se conforme admirablement au canon de la réversibilité thermodynamique absolue.
L'objectif industriel du procédé Claude est validé : nous avons propulsé de l'air sec dans l'abîme cryogénique (\(55.8\text{ K}\)) au moyen d'un effort mécanique global maitrisé.
Posez votre question