Physique : Thermodynamique

Systèmes ouverts et écoulements

Une turbine, une tuyère, un radiateur de voiture : la matière y entre et en ressort sans cesse. Pour les comprendre, on ne suit plus un morceau de fluide : on regarde une boîte fixe, et on compte ce qui franchit sa frontière.

Un système ouvert échange de la matière avec l'extérieur. On le délimite par un volume de contrôle, une frontière imaginaire et fixe. En régime permanent, rien ne change à l'intérieur au cours du temps : ce qui entre est exactement compensé par ce qui sort. Deux lignes suffisent alors.

ṁe = ṁs   et   q + wu = Δ(h + V²/2 + g z)
Tout est compté par kilogramme de fluide qui traverse, et positivement quand c'est reçu par le fluide : q est la chaleur, wu le travail utile (celui d'un arbre ou d'une hélice), h l'enthalpie massique, V la vitesse et z l'altitude.
1

Combien de matière traverse la boîte ?

Le débit massique ṁ est la masse qui franchit une section par seconde : ṁ = ρ·A·V, avec ρ la masse volumique, A l'aire de la section et V la vitesse. La matière ne disparaît pas : la masse contenue ne varie que si les débits d'entrée et de sortie diffèrent, dm/dt = ṁe − ṁs.

Chaque point représente la même masse de fluide. Les deux traits rouges sont les limites du volume de contrôle : le compteur note les points qui les franchissent.
Débit massique ρ₁A₁V₁
Vitesse de sortie
Contrôle : ρ₂A₂V₂ / ρ₁A₁V₁
Points entrés / sortis (compteur)
Le pouce sur le tuyau. Réduire la section accélère l'eau, parce que le même débit doit passer. Pour un gaz qui se dilate, il faut aussi tenir compte de ρ : à section constante, le gaz accélère quand même. La conservation porte sur la masse, pas sur le volume. Sur le compteur, l'écart entre entrées et sorties est le nombre de points en transit : il reste stable, rien ne s'accumule.
Réservoir de 1 m² de section, percé d'un orifice. Le débit de sortie dépend de la hauteur : V = √(2·g·h). Le temps est accéléré.
En haut : hauteur d'eau. En bas : débit entrant (bleu) et sortant (orange). Pointillés : valeurs du régime permanent.
Hauteur d'eau / temps écoulé
Débits entrant / sortant
dm/dt
Régime permanent : hauteur
Contrôle : masse − masse initiale − ∫(ṁₑ − ṁₛ)dt
Le régime permanent n'arrive pas par décret. Il s'établit tout seul : si le réservoir est trop plein, il se vide plus vite qu'il ne se remplit, et inversement. Pour le reste de la page, on suppose ce transitoire terminé, ce qui simplifie tout : dm/dt = 0, et ṁe = ṁs.
2

Pourquoi l'enthalpie remplace-t-elle l'énergie interne ?

Pour entrer dans la boîte, un kilogramme de fluide doit pousser celui qui le précède. Il reçoit du fluide en amont un travail P₁·v₁, et en sortie il cède P₂·v₂ à celui qui suit (v est le volume massique, en m³/kg). C'est le travail de transvasement. Il accompagne toujours le fluide : autant le ranger avec son énergie interne. On obtient h = u + P·v.

Le kilogramme bleu entre en poussant le fluide devant lui, poussé par celui qui le suit. Sa taille suit son volume massique : RT/P. Le kilogramme orange sort de la même façon.
Énergie d'un kilogramme, origine de u à 0 K. Bleu : u. Violet : P·v. Orange : V²/2 (ici, z reste constant). Les deux barres du milieu sont q et wu : vers le haut si le fluide reçoit.
Δh = cp·ΔT / Δu = cv·ΔT
Δ(V²/2)
Travail utile reçu wu
Par le bilan en u : Δu + Δ(V²/2) − q + P₂v₂ − P₁v₁
Verdict
Rien de nouveau, juste du rangement. C'est le premier principe appliqué au kilogramme qui traverse : Δu + Δ(V²/2 + g z) = q + wu + P₁v₁ − P₂v₂. En passant P·v à gauche, u + P·v apparaît tout seul. Attention : dans cette page, le gaz est parfait, donc P·v = R·T ne dépend pas de P. Le tableau montre l'ordre de grandeur de P·v : pour un liquide, il est très petit, et c'est ce qui rend le pompage si peu coûteux.

Le bilan compte trois sortes d'énergie, toutes par kilogramme. Sont-elles du même ordre ? Faites glisser les curseurs jusqu'à égaliser les barres. Les équivalences sont souvent surprenantes.

cp·ΔT
V²/2
g·Δz
Vitesse donnant autant que ΔT
Dénivelé donnant autant que ΔT
Un degré vaut une chute de cent mètres. Pour l'air, 1 K représente autant d'énergie qu'une vitesse d'environ 45 m/s ou qu'un dénivelé d'environ 100 m. Dans un échangeur ou un compresseur, on néglige donc souvent les termes de vitesse et d'altitude. Dans un barrage ou une tuyère, ce sont au contraire eux qui font tout le travail.

Une vanne, ou un simple étranglement, n'échange presque ni chaleur (trop rapide) ni travail (rien ne tourne). Si l'on néglige aussi la variation de vitesse, il reste h₂ = h₁. Pour un gaz parfait, h ne dépend que de T : rien ne change. Les gaz réels, eux, se refroidissent ou se réchauffent : c'est l'effet Joule-Thomson. Toutes les valeurs viennent d'une équation d'état réelle.

Le fluide passe de la pression P₁ à la pression de sortie. Couleur selon la température.
Température de sortie selon P₁, à T₁ fixée. Pointillés verts : le gaz parfait, T₂ = T₁. Point rouge : votre réglage.
Pression de sortie
T₂ / variation de température
Part de vapeur en sortie
Enthalpie h₁ = h₂ (origine CoolProp)
Le froid du détendeur. Dans un gaz réel, les molécules s'attirent : les écarter coûte de l'énergie, qui se prend sur l'agitation. C'est le principe du détendeur d'un réfrigérateur et de la liquéfaction des gaz (essayez le diazote à 130 K). Le dihydrogène fait l'inverse à température ambiante : il se réchauffe, car les répulsions y dominent. Une vanne conserve h, mais produit de l'entropie : c'est une perte nette, jamais compensée.
3

Que produisent une turbine et que coûte un compresseur ?

Dans ces machines, la chaleur échangée et les variations de vitesse et d'altitude sont petites. Le bilan se réduit à wu = h₂ − h₁ : le travail se lit directement sur l'enthalpie. La référence est la machine isentropique (adiabatique et réversible), et le rendement isentropique dit ce qu'on en obtient vraiment.

Le corps de turbine s'élargit parce que la vapeur se dilate énormément. La vapeur cède de l'enthalpie à l'arbre.
Diagramme (h, s), dit de Mollier. Cloche grise : zone liquide et vapeur. Bleu : isobares P₁ et P₂. Vert : détente idéale. Rouge : détente réelle, qui finit à droite car l'entropie augmente.
h₁ / h₂ réelle
Travail fourni / puissance
Isentropique : h₁ − h₂s / ∫v·dP
Sortie : T₂ et titre en vapeur
Entropie produite / volume massique sortie ÷ entrée
Chaque étage comprime et chauffe l'air (de gauche à droite), puis un refroidisseur le ramène à la température d'entrée. L'échelle de couleurs va du froid au chaud.
Plan (v, P), en échelles logarithmiques. Bleu : le trajet réel. Vert : l'isotherme, limite idéale. Orange : l'adiabatique en un seul étage.
Travail total pour 1 à 6 étages (rouge : votre choix), à votre rapport de pression. Pointillés verts : travail isotherme, R·T·ln r / η.
Travail total wu
Température après chaque étage
Chaleur évacuée par les refroidisseurs
Contrôle : wu − q − cp·(Tsortie − T₁)
Deux étages : meilleur rapport du 1ᵉʳ, à l'aveugle / √r
Un étage idéal : cp constant / air réel (CoolProp)
Comprimer froid coûte moins cher. Le travail idéal vaut ∫v·dP : un gaz chaud occupe plus de volume, il coûte plus à comprimer. D'où les refroidisseurs entre les étages, avec le même rapport pour chacun. Les étages rapprochent le travail de la limite isotherme, mais ne l'atteignent qu'à l'infini. Un moteur à gaz, lui, ne refroidit pas : la chaleur retirée serait autant de combustible à refournir.

Pour élever la pression, le travail isentropique vaut ∫v·dP, l'aire à gauche de la courbe. Or un liquide a un volume massique des centaines de fois plus petit qu'un gaz. Toute la différence entre une centrale à vapeur et une turbine à gaz est là.

Volume massique selon la pression, en échelle logarithmique verticale pour voir les deux. L'aire vraie (en échelle linéaire) est le travail affiché.
Comprimer 1 kg d'air à 20 °C
Pomper 1 kg d'eau à 20 °C
Rapport des deux
Part du travail reprise : cycle à vapeur (P₁ = 100 bar)
Part reprise : turbine à gaz (r = 16, 1 400 K, 85 %)
Pourquoi on condense la vapeur. Une centrale à vapeur prend la peine de tout liquéfier avant de remonter la pression : sa pompe ne reprend alors qu'une petite part de ce que fournit sa turbine. Une turbine à gaz doit comprimer un gaz, et son compresseur reprend une grande part du travail. C'est pourquoi les rendements de ses machines comptent tant.
4

Que font une tuyère et un échangeur ?

Ces deux organes n'ont aucune pièce mobile : wu = 0. Une tuyère, sans échange de chaleur, change l'enthalpie en vitesse : V₂²/2 − V₁²/2 = h₁ − h₂. Un échangeur, sans variation sensible de vitesse, fait passer de l'enthalpie d'un fluide à l'autre : ṁc·(he − hs)chaud = ṁf·(hs − he)froid.

Tant que le gaz va moins vite que le son (nombre de Mach M < 1), resserrer le conduit l'accélère. Au-delà de M = 1, c'est l'inverse : le gaz se détend si vite que c'est en élargissant le conduit qu'on l'accélère. D'où la forme en sablier, proposée par Laval pour ses turbines à vapeur.

Chaque point transporte la même masse. Couleur : température. Trait noir : le col, section la plus étroite. Trait rouge : onde de choc, le cas échéant.
Bleu : pression, rapportée à celle du réservoir. Orange : nombre de Mach. Pointillés rouges : pression extérieure.
Régime
Sortie : Mach / pression / température
Vitesse de sortie
Contrôle : h + V²/2 en sortie, écart à h₀
Contrôle : ρAV en sortie / au col
Débit / poussée (col : 100 cm², réservoir : 10 bar)
Entropie produite
Le col bloque le débit. Baissez la pression extérieure : le débit augmente, puis cesse d'augmenter dès que le col atteint M = 1. L'information ne peut plus remonter contre un écoulement supersonique. Plus bas encore, un choc droit se forme dans le divergent : la vitesse y chute et de l'entropie y est créée, sans que l'énergie h + V²/2 change. La tuyère n'est optimale que dite « adaptée », quand sa pression de sortie égale la pression extérieure.
Rouge en haut : fluide chaud. Bleu en bas : fluide froid. Les flèches blanches montrent le sens de chacun ; les flèches entre les deux, la chaleur transmise.
Températures le long de l'échangeur, intégrées pas à pas. À contre-courant, on ajuste la sortie du fluide froid par dichotomie jusqu'à respecter son entrée.
Sorties : chaude / froide
Puissance : côté chaud / côté froid
UA·ΔTlm (écart logarithmique)
Efficacité : intégrée / formule ε-NUT
Entropie produite
Pourquoi le contre-courant gagne. À débits et conductance égaux, l'écart de température y reste bien réparti, alors qu'en co-courant il se consume dès l'entrée. Le fluide froid peut même ressortir plus chaud que le fluide chaud n'en sort, ce qui est impossible en co-courant. Moins l'écart est grand partout, moins on crée d'entropie : on retrouve le « bain après bain » du second principe. Agrandir l'échangeur rapporte pourtant de moins en moins.
5

Que se passe-t-il quand on assemble tout ?

Un turboréacteur met bout à bout entrée d'air, compresseur, chambre de combustion, turbine et tuyère. Chaque organe suit le même bilan, et on peut placer le volume de contrôle où l'on veut, même autour du moteur entier. Et tout cela suppose un écoulement stable : voyons aussi comment il s'établit.

Coupe simplifiée. Le cadre pointillé rouge est le volume de contrôle choisi. Couleur : température totale. Le compresseur et la turbine sont sur le même arbre.
h + V²/2 par kilogramme de gaz, d'une station à l'autre. Le compresseur la fait monter, la chambre encore plus, la turbine redescendre, et la tuyère la garde mais la convertit en vitesse (violet).
Vitesse d'éjection / poussée (100 kg/s d'air)
Combustible : f = ṁc/ṁair / débit
Rendement thermique / cycle de Brayton idéal
Travail compresseur / turbine (même arbre)
Contrôle global : entrée + combustible − sortie
Le même bilan, partout. Cliquez sur chaque organe : ṁ se conserve (la chambre ajoute le combustible), et q + wu = Δ(h + V²/2) s'écrit de la même façon. Autour du moteur entier, le compresseur et la turbine disparaissent du bilan : il ne reste que le combustible qui entre et un gaz plus chaud et plus rapide qui sort. Ici, cp reste constant et sans pertes de pression : c'est un modèle de principe, optimiste par rapport à un vrai moteur. Avec des machines parfaites (100 %), le rendement thermique rejoint exactement celui du cycle de Brayton idéal, qui ne dépend que du rapport de pression total, entrée d'air comprise.

Voici la tuyère de la planche 4 (rapport de sections 1,6), non plus décrite par ses formules de régime permanent, mais calculée pas à pas par les lois de conservation de la masse, de la quantité de mouvement et de l'énergie, sur 150 petits volumes. Le gaz part du repos, et le réservoir s'ouvre : tout ce qui entre et sort est compté.

Bleu : pression calculée. Pointillés verts : solution de régime permanent de la planche 4, obtenue tout autrement. Le dégradé en haut montre le nombre de Mach.
Débit entrant (bleu) et sortant (orange), rapportés au débit bloqué ṁ*. Tant qu'ils diffèrent, la masse contenue varie.
Temps depuis l'ouverture
Débits entrant / sortant (÷ ṁ*)
Contrôle masse : M − M₀ − ∫(ṁₑ − ṁₛ)dt (÷ M₀)
Contrôle énergie (÷ E₀)
Choc : calculé / régime permanent
Sortie : h + V²/2, écart à h₀
Débit de sortie ÷ ṁ* : calculé / formule

Quelques repères

La masse

Dans une boîte fixe, dm/dt = ṁe − ṁs. En régime permanent, ρ·A·V ne change pas d'une section à l'autre : on peut donc en tirer les vitesses.

L'enthalpie

h = u + P·v range avec le fluide le travail qu'il faut pour le faire entrer et sortir. Bilan : q + wu = Δ(h + V²/2 + g z).

Les machines

Turbine et compresseur : wu = Δh. Tuyère : Δh se change en vitesse. Échangeur : Δh passe d'un fluide à l'autre. Vanne : Δh = 0.

Pour aller plus loin

Plusieurs entrées et sorties. Le bilan général s'écrit : dE/dt = Q̇ + Ẇu + Σe ṁ·(h + V²/2 + g z) − Σs ṁ·(h + V²/2 + g z). En régime permanent, dE/dt = 0. L'entropie suit la même logique : Σs ṁ·s − Σe ṁ·s = Σ Q̇/T + Ṡp.

Le travail réversible ∫v·dP. Pour un écoulement réversible, sans variation de vitesse ni d'altitude, wu = ∫v·dP. C'est différent du travail −∫P·dv d'un système fermé : la différence est justement le transvasement. Pour un liquide de densité constante, et sans chaleur ni travail, on retrouve l'équation de Bernoulli : P/ρ + V²/2 + g z = constante.

Les modèles. Tuyère : écoulement quasi unidimensionnel d'un gaz parfait à γ = 1,4, isentropique sauf au travers du choc droit. Hors de la tuyère, les chocs obliques et les faisceaux de détente ne sont pas représentés. Échangeur : U et les capacités constants, pas de pertes vers l'extérieur. Turboréacteur : cp = 1 005 J·kg⁻¹·K⁻¹, γ = 1,4 partout, pouvoir calorifique de 43 MJ/kg, vitesses internes négligées, pas de perte de pression, tuyère adaptée.

La simulation. Volumes finis à deux étages en temps (schéma HLL avec reconstruction linéaire limitée). Le choc y est étalé sur quelques mailles, d'où un léger écart de position. L'entrée suit le réservoir et la sortie impose la pression extérieure tant que le gaz y reste subsonique, par la méthode des caractéristiques.

Les données. Tables de vapeur d'eau, vanne et air réel : CoolProp 8.0.0. Les rendements usuels, les débits proposés, et la masse volumique de l'eau à 1 000 kg/m³ sont donnés à titre d'exemple. Le tableau des centrales de la page voisine donne des rendements réels complets.

Dans la même série : L'enthalpie et la calorimétrie, Le second principe et l'entropie, Moteurs, réfrigérateurs et pompes à chaleur, Changements d'état et diagrammes de phase.

J-C Gaillard