Final
Final
Final
et de la recherche scientifique
************
Université de Carthage
************
Institut National des Sciences
Appliquées et de Technologie
ii
2. Présentation du problème : (exemple 10.3) ................................................................................ 19
3. Résolution du problème : .............................................................................................................. 20
3.1. Résolution manuelle : ............................................................................................................ 20
3.2. Résolution avec HYSYS : ........................................................................................................ 22
3.2.1. Etapes de la résolution : ................................................................................................ 22
3.2.2. Résultats de la simulation HYSYS :................................................................................. 25
3.3. Comparaison des résultats : .................................................................................................. 26
4. Conclusion : ................................................................................................................................... 26
Conclusion générale : ............................................................................................................................ 27
Références bibliographiques : ............................................................................................................... 28
iii
Liste des figures
Figure 1:Principe de fonctionnement d'un compresseur........................................................................ 1
Figure 2 :choix des composés du gaz ...................................................................................................... 3
Figure 3 : Caractéristiques du gaz naturel à l’entrée .............................................................................. 4
Figure 4 : Désignation des caractéristiques du compresseur.................................................................. 4
Figure 5 : Schéma de la simulation .......................................................................................................... 5
Figure 6: Caractéristiques du gaz naturel à la sortie ............................................................................... 5
Figure 7: Puissance consommée par le compresseur d’efficacité 75%................................................... 6
Figure 8: caractéristiques du gaz à la sortie du compresseur ................................................................. 6
Figure 9 : Puissance consommée par le compresseur d’efficacité 10% .................................................. 7
Figure 10:Réprésentation détaillée de l'échangeur de chaleur à étudier ............................................... 9
Figure 11:Évolution des températures du fluide chaud et du fluide froid le long d’un échangeur à co-
courant .................................................................................................................................................. 10
Figure 12 :Évolution des températures du fluide chaud et du fluide froid le long d’un échangeur à
contre-courant ...................................................................................................................................... 11
Figure 13:Exemple d'une abaque de détermination du facteur correctif ............................................ 12
Figure 14: Choix des composés ............................................................................................................. 13
Figure 15: Choix du modèle thermodynamique ................................................................................... 14
Figure 16:Caractérsitques du fluide à l'entrée des tubes...................................................................... 14
Figure 17:Caractéristiques du fluide à l'entrée du calandre ................................................................. 15
Figure 18: Schéma de la simulation....................................................................................................... 15
Figure 19:Résultats affichés par ASPEN HYSYS...................................................................................... 16
Figure 20:Les composants d'une machine frigorifique ......................................................................... 17
Figure 21:Schéma de principe de fonctionnement d'une machine frigorifique ................................... 18
Figure 22:Le cycle frigorifique dans le diagramme de Mollier .............................................................. 19
Figure 23:Diagramme T-S du cycle frigorifique de R134a ..................................................................... 20
Figure 24:Diagramme de Mollier de R134a .......................................................................................... 21
Figure 25:Choix du fluide frigorigène .................................................................................................... 22
Figure 26:Choix du modèle thermodynamique .................................................................................... 22
Figure 27:Ajout du compresseur ........................................................................................................... 23
Figure 28:Caractérisation du fluide ....................................................................................................... 23
Figure 29:Ajout du condenseur ............................................................................................................. 24
Figure 30:Carctérisation du fluide à la sortie du condenseur ............................................................... 24
Figure 31:Ajout du détendeur et de l'évaporateur ............................................................................... 25
Figure 32:Schéma du cycle frigorifique ................................................................................................. 25
Figure 33:Résultats obtenus par HYSYS................................................................................................. 26
Figure 34:Abaque de détermination du facteur de correction ............................................................. 29
iv
Nomenclature
Variable Désignation Unité
Coefficient de compressibilité Ɣ
Efficacité du compresseur ɳ
Travail de l’arbre du Ws KJ/kg
compresseur
Puissance du compresseur Brake kW
Débit massique du gaz qm Kg/h
Fraction molaire x
Débit massique de l’eau qme Kg/s
Facteur de correction F
v
Efficacité du ɳ
compresseur
Production Qf kJ/kg
frigorifique
Travail du W kJ/kg
compresseur
Coefficient de COP
performance
Température de Tk °C
condensation du
fluide frigorigène
vi
Introduction générale :
La résolution de certains problèmes en génie des procédés nécessitent des calculs
complexes et parfois difficiles à réaliser manuellement.
Pour cela on utilise plusieurs méthodes numériques pour résoudre ces problème. Parmi
ces méthodes on cite le logiciel de simulation Aspen HYSYS qui permet de modéliser
mathématiquement des procédés.
On se propose dans ce projet de résoudre trois problèmes différents notamment :
➢ La compression d’un gaz naturel
➢ Etude d’un échangeur de chaleur
➢ Etude d’une machine frigorifique
vii
Problème N°1 : Compression d’un gaz naturel
1. Introduction :
2. Présentation du problème :
Un gaz naturel est comprimé constitué de méthane ,éthane, dioxyde de carbone et diazote
de 30°C et 3bar à 10 bar avec un débit massique est égal à 100kg/h.[2]
Données :[2]
Méthane 80 16
Ethane 10 30
Dioxyde de carbone 5 44
Diazote 5 28
1
Travail demandé :
3. Résolution du problème :
3.1. Résolution manuelle :
Application numérique :
ɣ 𝑹𝑻𝟏 𝑷𝟐
-Ws= ɣ−𝟏* *[(𝑷𝟏)(ɣ/ɣ-1)] [4]
𝑴
Application numérique :
1.31 8.3143∗(30+273.15) 10
-Ws= 1.31−1* *[( 3 )(1.31/(1.31-1))]
19.4
-Ws=181kJ/kg
Efficacité ɳ Brake KW
10% 50.27
75% 6.7
3
Figure 3 : Caractéristiques du gaz naturel à l’entrée
3.2.2 Résultats :
Le schéma de la simulation est présenté comme suit :
4
Figure 5 : Schéma de la simulation
5
Figure 7: Puissance consommée par le compresseur d’efficacité 75%
6
Figure 9 : Puissance consommée par le compresseur d’efficacité 10%
D’après la simulation les puissances consommées par les compresseurs sont données dans le
tableau suivants :
Tableau 3 : Les valeurs de la puissance consommées trouvées par Aspen
3.4 Commentaire :
La résolution avec HYSYS et la résolution manuelle donnent des résultats presque égaux.
7
4 Conclusion :
L’utilisation d’une simulation HYSYS peut donner des résultats différents résultats qu’un
calcul manuel.
Il faut utiliser les données exactes des variables caractéristiques du système étudié pour
avoir des bons résultats.
8
Problème N°2 : Etude d’un échangeur de chaleur
1. Introduction :
Un échangeur de chaleur est un appareil qui permet de réaliser un échange de chaleur
entre deux fluides à températures différentes.
Les échangeurs de chaleur sont couramment utilisés dans le domaine industriel, soit
pour la récupération, soit pour l’évacuation d’énergie. Différentes technologies peuvent
être rencontrées, à savoir, les échangeurs tubulaires simples, les échangeurs à tube et
calandre, les échangeurs à plaque etc.
2. Présentation du problème :
De l’eau chaude à 100°C et un débit de 5kg/s sert à chauffer la glycérine de 15°C à 55°C à
travers un échangeur à tube et calandre ,l’eau sort avec une température est égale à
55°C.[5]
3. Résolution du problème :
3.1. Etude théorique :
On se propose dans cette partie théorique d’expliquer la méthode de différence
logarithmique des températures et la facteurs de correction.
9
➢ celui à co-courants : écoulements parallèles des fluides et dans la même direction
➢ celui à contre-courants : écoulements parallèles des fluides mais dans des directions
inverses
Figure 11:Évolution des températures du fluide chaud et du fluide froid le long d’un échangeur à co-
courant
Dans notre cas, le fluide chaud est l’eau, le fluide froid est la glycérine.
10
Figure 12 :Évolution des températures du fluide chaud et du fluide froid le long d’un échangeur à
contre-courant
Dans ce cas, la différence de température logarithmique moyenne ΔTlm (DTLM) est définie
de la manière suivante :[6]
𝛥𝑇1 − 𝛥𝑇2
𝛥𝑇𝑙𝑚 =
𝛥𝑇
ln 𝛥𝑇1
2
Avec :
Te, Ts: températures entrée/sortie côté calandre
te, ts : températures entrée/sortie côté tube
La valeur du facteur correctif F est déterminé par l’utilisation des abaques tel que l’abaque
présentée ci-dessous
940.5
Application numérique : qmg= 2.48∗(55−15) = 9.45kg/s
qmg= 9.45kg/s
12
Déterminons maintenant la différence logarithmique de température et le facteur de
correction :
T1 − T2 45 − 40
Tlm,CF = = = 42.5 C
ln( T1 / T2 ) ln(45 / 40)
𝑡𝑠 − 𝑡𝑒 55 − 100
𝑃= = = 0.53
𝑇𝑒 − 𝑡𝑒 15 − 100
𝑇𝑒 − 𝑇𝑠 15 − 55
𝑅= = = 0.89
𝑡𝑠 − 𝑡𝑒 55 − 100
Dans ce problème le fluide qui est coté tube est l’eau et le fluide coté calendre est la glycérine
d’où à partir de l’abaque de la figure 34 on a F=0.77
En effet l’aire d’e=échange est donnée par la relation suivante :
A=nπDL avec n est le nombre de tubes
D est le diamètre de tubes
L est la longueur des tubes
A= 10*π*0.015*2= 0.94m2
Calcul de coefficient d’échange global :
𝑄̇ = 𝑈𝐴𝑠 𝐹𝛥𝑇𝑙𝑚,𝐶𝐹
𝑄̇
𝑈=𝐴
𝑠 𝐹𝛥𝑇𝑙𝑚,𝐶𝐹
Application numérique :
940.5
U=0.94∗0.77∗42.5 = 30.6kW/m2.°C
U=30.6kW/m2°C
13
• Ensuite on choisit le modèle thermodynamique qui est dans notre cas « NRTL »
14
Figure 17:Caractéristiques du fluide à l'entrée du calandre
3.3.2. Résultats :
On obtient les résultats suivants :
15
Figure 19:Résultats affichés par ASPEN HYSYS
3.5. Commentaire :
Les résultats obtenus théoriquement sont presque égaux à celle obtenus par HYSYS.
4. Conclusion :
Un problème d’étude d’un échangeur de chaleur est rapide et efficace avec le logiciel HYSYS
puisqu’il donne des valeurs très proches aux valeurs théoriques.
16
Problème N°3 : Etude d’une machine frigorifique
1. Introduction :
La machine frigorifique a pour rôle d'extraire de la chaleur (qu’on note Qf , >0) à la
source froide pour la refroidir d’avantage. Elle le fait tout en recevant un travail W (>0) du
milieu extérieur et en cédant une quantité de chaleur Qc (<0) à la source dite chaude.
17
De plus, il permet de faire circuler le fluide frigorigène dans le circuit fermé.
A la sortie du compresseur, le fluide frigorigène est à l’état vapeur et à haute pression.
➢ Condensation :
A l’entrée du condenseur, le fluide frigorigène est à l’état vapeur et à haute pression.
En passant dans le condenseur, le fluide frigorigène (à haute température) cède son énergie
thermique. De ce fait, le fluide frigorigène se condense et passe à l’état liquide.
A la sortie du condenseur, le fluide frigorigène est à l’état liquide et à haute pression.[9]
➢ Détente :
A l’entrée du détendeur, le fluide frigorigène est à l’état liquide et à haute pression.
Lorsque le fluide frigorigène traverse le détendeur, sa pression ainsi que sa
température diminuent.
Le détendeur permet également de régler le débit de fluide frigorigène parcourant
le circuit fermé.
A la sortie du détendeur, le fluide frigorigène est à l’état liquide et à basse pression.
➢ Evaporation :
A l’entrée de l’évaporateur, le fluide frigorigène est à l’état liquide et à basse pression.
A la sortie de l’évaporateur, le fluide frigorigène est à l’état vapeur et à basse pression.
18
Figure 22:Le cycle frigorifique dans le diagramme de Mollier
La production frigorifique :
Qf= h1-h4
Le coefficient de performance frigorifique COPf qui est le rapport de la puissance
frigorifique 𝑄̇𝑓 sur la puissance mécanique consommé 𝑊̇ :
ℎ −ℎ
COPf=ℎ 1 −ℎ4
2𝑠 1
ℎ1 −ℎ4
COPf= ℎ2 −ℎ1
19
une température égale à 30°C. Le fluide circule avec un débit massique est égal à 0.08kg/s et
le compresseur a une efficacité est égale à 80%.[11]
Travail demandé :
On se propose de :
Calculer la puissance du compresseur
Calculer le coefficient de performance COP de la machine frigorifique
3. Résolution du problème :
3.1. Résolution manuelle :
Afin de déterminer le coefficient de performance on doit faire la lecture graphique des
valeurs des enthalpies à partir le diagramme de Mollier.[12]
20
Figure 24:Diagramme de Mollier de R134a
Application numérique :
272.39−241.35
h2= +241.35= 280.15kJ/kg
0.8
𝑊̇ =3.1kW
Le coefficient de performance est donné par la relation suivante
𝑄 ℎ −ℎ
COP=𝑤𝑓 =ℎ1 −ℎ4
𝐶 2 1
COP=3.86
21
3.2. Résolution avec HYSYS :
3.2.1. Etapes de la résolution :
On choisit le fluide frigorigène à utiliser
On ajoute le compresseur
22
Figure 27:Ajout du compresseur
23
Figure 29:Ajout du condenseur
24
Figure 31:Ajout du détendeur et de l'évaporateur
25
Et les résultats obtenus sont les suivants :
Les résultats obtenus par le calcul théorique sont presque égaux à celle égaux à celle trouvés
par HYSYS.
L’erreur qui peut parvenir à partir la lecture de diagramme peut être remarquable.
4. Conclusion :
L’utilisation de HYSYS pour la détermination de l’efficacité d’une machine frigorifique est
très rapide que le calcul manuel en trouvant des valeurs proches.
26
Conclusion générale :
Ce Projet Personnel Professionnel a été structuré en trois parties. Dans la première partie,
nous avons dimensionner un compresseur du gaz.
Dans la deuxième, nous avons étudié un échangeur de chaleur. Et finalement nous avons fait
l’étude d’une machine frigorifique à compression. Ce travail a montré que la résolution des
problèmes en génie chimique avec des méthodes numériques tel que Aspen HYSYS donne
des résultats précis et en moins du temps que le calcul manuel.
27
Références bibliographiques :
[1] https://fr.wikipedia.org/wiki/Compresseur_m%C3%A9canique
[2] 293535261-Computer-Methods-in-Chemical-Engineering-Nayef-Ghasem-CRC-2012
[3] https://uel.unisciel.fr/chimie/strucmic/strucmic_ch10/co/sevaluer_ch10_02_6.html
[4] https://direns.mines-paristech.fr/Sites/Thopt/fr/res/ComprVolum.pdf
[5] Heat Chap13-054.docx
[6] http://tonepi.com/systemes-energetiques/difference-de-temperature-logarithmique-
moyenne-entre-les-deux-fluides-dtlm.html , consulté le 13/6/2021
[7] http://processs.free.fr/Pages/VersionWeb.php?page=4552 , consulté le 13/6/2021
[8] les propriétés thermodynamiques disponible sur http://www4.ac-nancy-
metz.fr/physique/ancien_site/divers/liq-prop/Liq-prop.htm , consulté le 13/05/2021
[9]https://www.ef4.be/fr/pompes-a-chaleur/technique-generalites/principe-de-
fonctionnement.html ,consulté le 11/6/2021
[10] https://www.ef4.be/fr/pompes-a-chaleur/technique-generalites/cycle.html , consulté le
12/6/2021
[11] Fundamentals-of-Engineering-Thermodynamics.pdf
[12] https://www.abcclim.net/diagramme-mollier.html
28
Annexe:
29