Université de Liège
Faculté des Sciences Appliquées
N. Job, Professeur ordinaire Année académique 2024-2025
Thermodynamique Chimique Appliquée – corps purs : projet
On souhaite installer une pompe à chaleur dans un procédé afin de réchauffer un
fluide. La puissance thermique à fournir au procédé en question, via le condenseur
de la pompe à chaleur, est de 100 kW.
Le schéma de la pompe à chaleur est donné ci-dessous :
Description de la pompe à chaleur
Un fluide Y circule en boucle entre un évaporateur, un compresseur, un condenseur
et une vanne de détente avec un débit F.
- Dans l’évaporateur, le fluide se vaporise à pression constante, ce qui absorbe
une puissance thermique PQ, provenant de la source froide. On suppose que
le fluide n’est pas surchauffé, et sort sous forme de vapeur saturante à T1 = 50
°C ;
- Le compresseur, considéré comme adiabatique, porte la vapeur à la pression
P2. La température de sortie du compresseur, T2, est mesurée égale ou supé-
rieure à 100°C. La compression nécessite l’apport d’un travail W ;
- Dans le condenseur, le fluide est d’abord refroidi, à pression constante,
jusqu’à sa température de saturation. Ensuite, il se condense à température et
pression constante. Il sort sous forme de liquide saturé à la température T3 =
100°C ;
- Dans la vanne de détente, considérée comme adiabatique, le fluide est ra-
mené à la pression P1 et à une température T4. Ce fluide est admis dans
l’évaporateur et le cycle est ainsi bouclé. Comme la transformation est adiaba-
Groupe n° Fluide utilisé :
tique, donc à enthalpie constante, on peut évaluer l’état final de la détente (At-
tention : au terme de la détente, il se peut qu’une partie du liquide se soit va-
porisé).
Le but du projet est d’évaluer le cycle, c’est à dire de calculer le débit F, le travail
minimum du compresseur W (et la puissance correspondante PW), la chaleur extraite
à l’évaporateur Qvap (et la puissance correspondante PQ), et le coefficient de
performance maximum du cycle pour une puissance constante fournie au procédé
de 100 kW.
Le fluide utilisé dans le cycle sera précisé dans la suite de votre travail. Vous ne
disposerez d’aucune donnée sur cette substance. Seul un diagramme enthalpique
du fluide Y vous sera fourni. Les différents fluides sont les suivants :
Code réfrigérant Nom
R11 Trichlorofluorométhane
R12 Dichlorodifluorométhane
R21 Dichlorofluorométhane
R40 Chlorométhane
R113 1,1,2-Trichloro-1,2,2-trifluoroéthane
R114 1,2-Dichloro-1,1,2,2-tétrafluoroéthane
R123 2,2-dichloro-1,1,1-trifluoroéthane
R124 1-Chloro-1,2,2,2-tétrafluoroéthane
R141b 1,1-dichloro-1-fluoroéthane
R142b 1-chloro-1,1-difluoroéthane
R152a 1,1-difluoroéthane
R236ea 1,1,1,2,3,3-hexafluoropropane
R236fa 1,1,1,3,3,3-Hexafluoropropane
R245ca 1,1,2,2,3-pentafluoropropane
R245fa 1,1,1,3,3-pentafluoropropane
R365mfc 1,1,1,3,3-pentafluorobutane
R600 n-Butane
R600a 2-Méthylpropane
R601a 2-Méthylbutane
R717 Ammoniac
RC318 Octafluorocyclobutane
Remarque : on ne peut pas avoir de liquide à l’entrée du compresseur. La vapeur
doit être au minimum saturé.
Pour ce projet, on vous demande :
1) De trouver des données thermodynamiques pour votre fluide : température et
pression critiques, température normale d’ébullition, cp en phase liquide et cp
en phase vapeur, masse moléculaire, enthalpie de vaporisation, facteur
acentrique ainsi que tout autre donnée que vous jugeriez utile. On vous
demande de vérifier ces données à l’aide de deux sources différentes, la
première étant « National Institute of Standards and Technology »
(www.nist.gov). Sur ce site, le lien intitulé « fluid properties » vous permet de
calculer de nombreuses propriétés thermodynamiques dans les conditions de
température et de pression que vous choisissez vous-même.
Groupe n° Fluide utilisé :
2) De présenter le cycle sur un diagramme P=f(T). Pour chaque point du cycle,
indiquer la valeur de la température si celle-ci est connue. Indiquer également
l’état du fluide en chaque point (vapeur saturée, liquide, vapeur surchauffée,
etc.).
3) D’établir, de manière théorique, le bilan d’énergie pour chaque appareillage, et
d’en déduire les relations entre le travail mécanique et l’enthalpie du fluide en
chaque point. Pour vous aider, indiquez clairement, pour chaque
transformation, de quel type de transformation il s’agit et sous quelle condition
particulière (par exemple : détente / compression/… et adiabatique / isobare /
isenthalpique…). Faire apparaître dans vos bilans les 100 kW d’énergie
transmis à la source chaude. Pour faciliter la compréhension, on gardera les
mêmes points que ceux sur le schéma de la page 1.
4) De déterminer, grâce au diagramme enthalpique, quelle seront les pressions
basse et haute de votre cycle. Quelles sont les enthalpies de vaporisation
massique à ces deux températures ainsi que l’enthalpie standard de
vaporisation ?
5) D’utiliser le diagramme enthalpique pour évaluer l’ensemble du cycle (calculer
le débit F de fluide dans le cycle, le travail du compresseur Wmin (et la
puissance correspondante PW), la chaleur échangée Qvap à l’évaporateur (et
la puissance correspondante PQ), et le COP maximum du cycle.
6) Pour chaque point du cycle, on vous demande de déterminer théoriquement
comment calculer l’enthalpie du fluide en supposant que la vapeur peut être
considérée comme un gaz parfait. On prendra h(Tref) = 0 pour Tref = 273,15
K. Quel(s) paramètre(s) est (sont) nécessaire(s) pour calculer ces enthalpies ?
Comment allez-vous déterminer la température T2 à la sortie du
compresseur ainsi que la fraction vapeur à la sortie du détendeur ? Avec
quelle méthode vue au cours, allez-vous le(s) déterminer ?
7) De calculer l’enthalpie du fluide en chaque point du cycle en considérant que,
en phase vapeur, le fluide peut être considéré comme un gaz parfait. On
prendra h(Tref)=0 pour Tref = 273,15 K. Pour Cp = f(T), on utilisera une
régression de type cp(T)=A+B*T utilisant les données fournies pour chaque
fluide. Pour le calcul de l’enthalpie liquide, on prendra l’enthalpie de
vaporisation déterminée graphiquement au point 5. Les résultats obtenus
sont-ils logiques ? Pourquoi ?
8) De calculer F, Qvap, Wmin, PQ et PW toujours en supposant que, en phase
vapeur, le fluide peut être considéré comme un gaz parfait et que le
compresseur à une efficacité isentropique de 75%. Recalculer la température
T2 à la sortie du compresseur;
9) De calculer les paramètres du fluide nécessaires pour vos calculs par des
méthodes de contribution de groupe :
Tc et Pc par Joback-Lydersen ;
Cp par Joback-Lydersen ;
Le facteur acentrique par (i) sa définition et (ii) Constantinou et Gani.
Certaines molécules ne peuvent pas être calculées avec
Constantinou-Gani, dans ce cas, calculer le facteur acentrique
pour une molécule de référence : R21.
10) On vous demande de comparer les valeurs obtenues avec celles disponibles
sur le site du NIST. Pour le Cp, on vous demande de comparer ces valeurs
sur un graphique dans la zone de température désirée.
Groupe n° Fluide utilisé :
11) De calculer les pressions en chaque point du cycle :
a. A partir de la définition du facteur acentrique ;
b. En utilisant la corrélation de Pitzer.
12) D’améliorer les calculs réalisés précédemment en considérant que la vapeur
en entrée et en sortie du compresseur doit être traitée comme un fluide
réel. On calculera les grandeurs résiduelles en utilisant la méthode de Lee-
Kesler. On gardera pour les autres enthalpies les valeurs calculées avec
l’équation des GP, les pressions P1 et P2 seront celles calculées avec
Pitzer ! On prendra la température T2 qui a été réévaluée en tenant
compte de l’efficacité isentropique du compresseur (75%) sauf si T2
calculé <100°C, alors prenez T2=100°C.
13) D’améliorer les calculs réalisés en considérant que la vapeur en entrée et en
sortie du compresseur doit être traitée comme un fluide réel et répond à
l’équation de Soave. On gardera pour les autres enthalpies les valeurs
calculées avec l’équation des GP !
Remarque : pour Soave on arrive à une équation cubique
Z³ - Z² + (A-B-B²)Z-AB= 0
Veuillez préciser les valeurs de A et B
Remarque générale :
1) Pour les points 7,8,9, 10, 11, 12, 13, 14, 15 et 16 votre réponse sera toujours
précédée d’un tableau récapitulatif reprenant la valeur des paramètres
nécessaires à la réalisation de vos calculs
2) En bas de chaque page de votre projet, veuillez préciser votre groupe ainsi
que le fluide utilisé pour vos calculs !
3) La première page de votre rapport sera la page suivante complétée
Groupe n° Fluide utilisé :
Projet de thermodynamique chimique appliquée
Groupe n°….
Nom de l’étudiant 1 :
Nom de l’étudiant 2 :
Fluide Formule
considéré chimique
NIST Autre Joback
source :….
Tb(K) X
Tc(K)
Pc(bar)
NIST ou autre Définition Constantinou
source
Diagramme Définition du Corrélation de Pitzer
enthalpique facteur
acentrique
P2 : Pression haute
du cycle (bar)
P1 : Pression basse
du cycle (bar)
Diagramme Hypothèse Tables de Équation
enthalpique GP Lee de Soave
Kesler
T2 (K)
H1 (kW)
H2 (kW)
H3 (kW)
F (kg/s)
Pw(MW)
PQ(MW)
COP
Groupe n° Fluide utilisé :
