THERMO-ENERGETIQUE

Enoncés de TD
1A GI

Contact : [Link]@[Link]
 Thermo-Energétique

Machines cycliques pour la production d’énergie mécanique

et thermique

1. Production d’énergie mécanique par un cycle à gaz

Une installation de production d’énergie électrique mettant en œuvre un cycle à air est
représentée sur la figure 1. L’air extérieur (point 1) subit une compression en trois étages avec
refroidissement intermédiaire, puis se préchauffe dans un échangeur thermique appelé
‘récupérateur’, avant d’alimenter (point 7) une 1ère chambre de combustion CC-1. Le gaz sortant
de celle-ci (point 8) est détendu dans la turbine ‘haute pression’ T-HP, entraînant les 3
compresseurs.
Après une deuxième chambre de combustion CC-2, le gaz (point 10) est détendu à travers
la turbine basse pression T-BP, entraînant un alternateur. Le gaz d’échappement de cette turbine
(point 11) sert à préchauffer l’air comprimé dans le récupérateur, avant d’être rejeté (point 12)
à l’atmosphère.

Hypothèses :
- On supposera que les compressions et les détentes sont isentropiques et que toutes les
autres transformations (combustion, échange de chaleur) sont isobares.
- Les gaz seront assimilés à de l’air : on négligera les modifications de débit massique et
de composition provoqués par les combustions, celles-ci s’effectuant avec un important excès
d’air.
- On considèrera l’air comme un gaz parfait de constante massique R = 0,287 kJ/(kg.K)
et dont le rapport des chaleurs massiques est :  =1,35.

Données :
- Températures à l’entrée des compresseurs : T1 = T3 = T5 = 25°C
- Températures à l’entrée des turbines : T8 =T10 = 680°C
- Température à l’entrée et à la sortie du récupérateur : T6 = 106°C T7 = 300°C
- Température du gaz en sortie de turbine HP : T9 = 444°C
- Température du gaz rejeté : T12 = 160°C
- Pressions : P1 = P11 = P12 = 1.105 Pa P2 = P3 = 2,5.105 Pa
P4 = P5 = 6.105 Pa P6 = P7 = P8 = 15.105 Pa
P9 = P10 = 5.105 Pa
- Quantités de chaleur apportées dans les chambres de combustions :
QCC1 = 421 kJ/kg QCC2 = 261,5 kJ/kg

2
 Thermo-Energétique

Travail à réaliser :
1) Sur un diagramme entropique (T,S) qualitatif, tracer le cycle parcouru par l’air. Quel est
ce cycle mis en œuvre ?
2) Calculer le travail massique de compression pour chacun des 3 compresseurs.
3) Déterminer la température T11 à la sortie de la turbine basse pression T-BP. En déduire
le travail massique fourni par le gaz au rotor de cette turbine.
4) Déterminer le travail fourni par la turbine T-HP.
5) Effectuer le bilan des travaux massiques sur chacune des deux lignes d’arbres des
turbines T-BP et T-HP (travaux produits par les turbines et absorbés par les
compresseurs et l’alternateur). On négligera les pertes mécaniques de transmission.
Déterminer le rendement thermique de l’installation.
6) Calculer le débit massique d’air (tonnes/h) nécessaires à une installation produisant 5
MW.

Figure 1 : Schéma de l’installation motrice

3
 Thermo-Energétique

15 bar 6 bar 5 bar 2,5 bar 1 bar

680°C

300°C

25°C

2. Usine avec ou sans centrale de cogénération

Une usine consomme en énergie :
- Puissance électrique : 28 MW.
- Puissance thermique : 37,2 MW à T>60°C.

I- Usine élémentaire (sans centrale de production d’énergie)

- L’énergie thermique est fournie par combustion d’un fuel dont le pouvoir calorifique
est de : pc=4,08.104 kJ/kg
- L’énergie électrique est fournie intégralement par le réseau national à partir d’une
centrale électrique fonctionnant au fuel, dont le rendement thermique vaut r = 0,38. Le fuel

4
 Thermo-Energétique

utilisé par cette centrale possède le même pouvoir calorifique que ci-dessus et les pertes en
ligne sont évaluées à 20%.

Dans ces conditions, la consommation de fuel est :

- Pour l’énergie thermique : Pth/Pc = 37,2.106/4,08.107 = 0,911 kg/s = 3282 kg/h de fuel

- Pour l’énergie électrique :

o Pertes en lignes de 20% donc en sortie de la centrale : Pelec = 28/0,8 = 35 MW
o Rendement (r =Pelec/Pth) de la centrale de 0,38 donc Pth à fournir par la
centrale : 35/0.38 = 92,1 MW
o Débit de fuel : 92.1.106/4,08.107 = 2,257 kg/s = 8125 kg/h de fuel

Consommation totale de fuel = 8125 + 3282 =11407 kg/h

II- Usine avec centrale de production d’énergie

En réalité l’usine comporte une centrale à vapeur qui produit une partie du besoin en
énergie électrique et la totalité de l’énergie thermique nécessaire sous forme de vapeur.
L’appoint en énergie électrique est fait par le réseau national.
Dans le circuit de la centrale de l’usine, l’eau liquide à 15°C et 1 bar (point A) est
pompée pour alimenter sous 30 bar (point B) une chaudière chauffée par combustion de
fuel. L’élévation de température lors de ce pompage est faible puisque la température de
l’eau au point B est de 16°C. La vapeur sortant de la chaudière, surchauffée à 480°C (point
D), est détendue jusqu’à 2 bar à travers une turbine (avec un rendement isentropique  =
0,9) entraînant un alternateur.
La vapeur issue de la turbine (point E) est envoyée dans les circuits d’utilisation comme
moyen de chauffage. Dans ces circuits, le condensat demeure utilisable comme fluide
chauffant jusqu’à une température de 60°C (point F).
Dans ce problème et pour toute l’installation, on ne tiendra pas compte des pertes de
charge, pertes thermiques, ainsi que des pertes mécaniques par transmission.

Données :
- Diagramme entropique de l’eau
- Enthalpies : HA = 63 et HF = 251 (en kJ/kg)

Travail à réaliser :
1) Sur le diagramme entropique (T,S), tracer le cycle parcouru par l’eau dans la centrale et
dans les circuits d’utilisation. Compléter le tableau des points de fonctionnement. Quel
est le cycle mis en œuvre ?

5
 Thermo-Energétique

2) Calculer la quantité de chaleur cédée par chaque kg de vapeur de chauffage dans les
circuits d’utilisation.
3) Déterminer le débit massique de l’eau circulant dans la centrale de l’usine.
4) Calculer la puissance électrique produite par cette centrale.
5) Comparer la puissance mécanique nécessaire au pompage de l’eau (trajet A-B) à celle
produite par la détente de la vapeur à travers la turbine (trajet D-E)
6) Calculer le débit de fuel consommé par la chaudière de la centrale de l’usine.
7) Quelle est l’incidence de cette usine sur la consommation nationale de fuel ?

Point Etat T (°C) P (105 Pa) H (kJ/kg)

A Liq. 15 1 63
B Liq. 16 30 66
C Vap. Sat. 234 30 2800
D Vap. Surch. 480 30 3415
E Vap. 170 2 2817
F Liq. 60 2 251

AB: pompage de l’eau.

BD: transformation à travers la chaudière (chauffage de l’eau, vaporisation, surchauffe).
DE : détente de la vapeur à travers la turbine.
EF: transformation dans le circuit d’utilisation comme moyen de chauffage
(refroidissement de la vapeur, condensation, refroidissement du condensat).

152 30 bar 480 C

30 bar
I bar

2 bar
60 %
2 bar

Figure 1 : représentation schématique de l’installation

6
 Thermo-Energétique

&
·
S

·*

&

&

⑧ ⑧ &

-
:

7
 Thermo-Energétique

3. Production de froid par cycle à compression de vapeur – Partie 1

Une installation frigorifique à compression de vapeur, de puissance QE = 100 kW à
l’évaporateur, doit permettre de maintenir la température d’une chambre froide à 4°C. La
température d’évaporation du fluide est fixée à TE = -5°C. Le condenseur est refroidi par de
l’eau à 20°C. La température de condensation est choisie à Tc = 38°C. On ne prendra pas en
compte les pertes de charges dans l’installation.
La machine frigorifique envisagée suit un cycle de Rankine classique : compression,
condensation, détente, évaporation (figure 1).
On suppose que :
- La compression s’effectue en une étape de façon isentropique.
- La détente est également en une étape et se fait de manière isenthalpique
- Le fluide frigorifique liquide sortant du condenseur est saturé. Il en est de même pour
la vapeur sortant de l’évaporateur.

I- Utilisation du diagramme enthalpique du fréon 12 (R-12 ou

dichlorodifluorométhane)
Le Fréon 12 est utilisé comme fluide frigorifique. Nous disposons de son diagramme
enthalpique.

1) Quelles sont les pressions à l’évaporateur et au condenseur ?

2) Tracer le cycle parcouru par le fréon sur le diagramme (P,H) et en déduire la fraction
vaporisée au point 4 et la température T2.
3) Calculer le débit massique de fréon 12 (kg/h) circulant dans l’installation.
4) Déterminer la puissance de compression.
5) En déduire le COP de l’installation.

II- Calcul du cycle du fréon 12 par l’équation d’état gaz parfait

Nous utilisons toujours le fréon 12 comme fluide frigorifique mais nous supposons ne
plus disposer de son diagramme enthalpique. Les grandeurs lues précédemment devront donc
être maintenant calculées. Les données disponibles sur ce corps pur sont :

Loi d’Antoine : log10(P) = A − (B / (T + C))

P = Tension de vapeur (bar) T = température (K)
A = 3,81109 B = 782,072 C = -37,773 Tmin = 154,7 K Tmax = 243,4 K

Les Cpgaz et Cpliq seront considérés constants : Cpgaz = 71,06 J/mol/K ; Cpliq = 112,6 J/mol/K
hvap(-5°C) = 19051,6 J/mol hvap (38°C) = 16250 J/mol
Tc = 384,95 K Pc = 40,71 atm M = 120,9 g/mol

8
 Thermo-Energétique

1) Quelles sont les pressions à l’évaporateur et au condenseur (à calculer avec la loi

d’Antoine)?
2) Calculer la température du point 2 avec l’hypothèse comportement gaz parfait.
3) Calculer les enthalpies spécifiques des points 1, 2 et 3.
4) Calculer le débit massique de fréon 12 (kg/h) circulant dans l’installation.
5) Déterminer la puissance de compression.
6) En déduire le COP de l’installation.
↑

7) Commenter les résultat obtenus.

Eau de refroidissement

T = 4°C

Figure 1 : Cycle de base

9
 /
Thermo-Energétique

E &

-
·
I i

-
--

10
 Thermo-Energétique

EQUATIONS D’ETAT - PROPRIÉTÉS THERMODYNAMIQUES DES

CORPS PURS

4- Calcul d'une enthalpie de vaporisation

On se propose de calculer l'enthalpie de vaporisation du 1-3 butadiène à 140°C, c'est-à-
dire dans des conditions telles que la phase vapeur n'est plus un gaz parfait et que le volume
molaire du liquide n'est plus négligeable devant celui du gaz.
Pour cela, on supposera que l'équation d'état de Soave Redlich et Kwong rend
convenablement compte du comportement du 1-3 butadiène.
1) Tracer, à 140°C, les variations de la pression du 1-3 butadiène, données par l'équation
d'état de Soave Redlich et Kwong, en fonction de son volume molaire (0,1 l/mol < v <
0 , 02 pas
1 l/mol) 0, 1 pas
2) En déduire, pour la tension de vapeur à 140°C, les valeurs des volumes molaires du 1-3
butadiène gaz et liquide.
3) Calculer, par la relation différentielle de CLAPEYRON, l'enthalpie molaire de
vaporisation du 1-3 butadiène. On rappelle que sur un intervalle de température peu
étendu, le logarithme de la tension de vapeur varie linéairement avec 1/T.

5- Fonctions d'état du propylène

1) Pour un fluide obéissant à l'équation d'état de Redlich et Kwong, trouver l'expression de

l'enthalpie molaire résiduelle h-h*.

2) Pour ce même fluide trouver l'expression de l'entropie molaire résiduelle s-s*.
3) En supposant que le propylène suit l'équation d'état de Redlich et Kwong, calculer
l'enthalpie molaire du propylène à 105 °C sous 30 bar, en prenant comme origine des
enthalpies le propylène gaz parfait à 25°C. Nous utiliserons pour cet exercice un Cp gaz
constant pour l’éthylène et égal à 70,5 J/mol/K.
4) On détend en régime permanent de manière isenthalpe (détente libre de JOULE-
THOMSON dans un accident de parcours) du propylène de 105°C, 30 bar à 1 bar.
Quelle est la température du propylène après détente.
5) On effectue maintenant cette détente en régime permanent dans une turbine adiabatique
fonctionnant réversiblement, quel travail utile est récupéré par mole de propylène
détendue.

11
 Thermo-Energétique

6- Calculs de fugacités
1) Calculer la fugacité de l'eau gaz à 251,1°C sous 40 bar par l'équation de PITZER.
2) Calculer la fugacité de l'eau gaz à 251,1°C sous 40 bar par l'équation de PENG et
ROBINSON.
3) Calculer la fugacité de l'eau liquide à 251,1°C sous 50 bar par l'équation de PENG et
ROBINSON.
4) Calculer la fugacité de l'eau liquide à 251,1°C sous 50 bar par la méthode de la tension
de vapeur.

Données complémentaires:
volume molaire expérimental de l'eau liquide à 251.1°C = 0,0225 [Link]-1

 
1 1 1 1
on rappelle: 2   

v  2.b.v  b 2  
2 2 .b  v  1  2 .b  

v  1  2 . b 

7- Calcul d'équilibre liquide-vapeur de corps pur par équation d'état

Le benzène suit l'équation d'état de Soave Redlich et Kwong.
1) Calculer par cette équation d'état les tensions de vapeur du benzène pour les
températures:
80,1°C 103,8°C 142,5°C 178,8°C 221,5°C 249,5°C 272,3°C

et comparer les avec les valeurs expérimentales.

2) Calculer, à l'aide des enthalpies résiduelles, l'enthalpie de vaporisation du benzène à
249,5°C.
3) Calculer cette même enthalpie de vaporisation à 249,5°C en utilisant la relation de
Clapeyron, les valeurs expérimentales des tensions de vapeur et l'équation d'état de
Soave Redlich et Kwong pour déterminer les volumes molaires.
4) Calculer cette même enthalpie de vaporisation à 249,5°C en utilisant l'enthalpie normale
de vaporisation (sous 1 bar) et les chaleurs molaires du benzène gaz et liquide.
On utilisera pour cela le tableur Excel et les fonctions personnalisées de résolution
d'équations du troisième degré mises à disposition.

12
 Thermo-Energétique

Rappel sur l'équation SRK:

 a    b 

m  0,480  1,574.  0,176. 2 ,   1  m. 1  Tr  
2

 .Pr
A  a
T 2
, B  b
Pr
Tr
 
, z 3  z 2  A  B  B 2 .z  A.B  0
r

h  h*  A T   B 
  1  z  .1  m. r . ln 1  
R.T  B    z 

s  s* A Tr  B
 ln  z  B   .m. . ln 1  
R B   z
f A B
ln( )  z  1  ln( z  B )  ln(1  )
P B z

8- Production de froid par cycle à compression de vapeur – Partie 2

III- Calcul du cycle du fréon 12 par l’équation d’état RK
Nous reprenons la machine frigorifique de l’exercice 3. Les hypothèses de fonctionnement sont
toujours les mêmes. Au point d’entrée du compresseur le fréon-12 sera considéré comme un
gaz parfait, la pression étant de 2,5 bar. Par contre, au point de sortie, la pression étant de 9 bar,
l’équation d’état utilisée sera Redlich-Kwong. Nous rappelons que les expressions analytiques
des enthalpies et entropies résiduelles avec l’équation de Redlich-Kwong ont été établies dans
l’exercice 5.

1) Ecrire et résoudre les 2 équations permettant de trouver le volume molaire gazeux

et la température au point 2.
2) Calculer les enthalpies spécifiques des points 1, 2 et 3.
3) Calculer le débit massique de fréon 12 (kg/h) circulant dans l’installation
4) Déterminer la puissance de compression et en déduire le COP de l’installation.
5) Comparer les résultats obtenus avec ceux de la partie I.

13
 Thermo-Energétique

IV- Utilisation de Prophy+ pour le cycle du fréon 12

Nous utilisons toujours le fréon 12 comme fluide frigorifique mais seul le logiciel
Prophy+ est disponible pour réaliser les calculs. Nous utiliserons le modèle thermodynamique
idéal dans un premier temps.
1) Calculer les pressions d’équilibre à 38°C (HP) et -5°C (LP).
2) Calculer les enthalpies spécifiques des points 1 et 3 (Figure 2).
3) Connaissant la pression LP et l’enthalpie h3, calculer la fraction vaporisée au point 4.
Calculer l’entropie au point 1.
4) Calculer la température au point 2, connaissant HP et s1. Analyser les résultats obtenus.
Que proposez-vous ?
5) Déduire le débit de fréon 12 pour obtenir les 100 kW souhaités à l’évaporateur. Calculer
la chaleur échangée au condenseur ainsi que le travail minimum à fournir au
compresseur.
6) Comparer les résultats obtenus avec ceux des parties I et III
7) Reproduire la Figure ci-dessous sur excel après transfert des résultats fournis par
Prophy+:
- La courbe de saturation peut être obtenue avec des calculs pour P variant de 0.1 atm
jusqu’à une valeur proche de Pc (500 points).
- Les isothermes 38°C et - 5°C avec des calculs pour P variant de 0.1 atm jusqu’à une
valeur proche de Pc (500 points).

Diagramme Log(P)=f(H)
LogP

-5°C 38°C

3
HP
2

LP
1
4

H(cal/mole)

14