Laouar Imed

‫وزارة اﻟﺘﻌﻠﯿــﻢ اﻟﻌﺎﻟـﻲ واﻟﺒﺤـﺚ اﻟﻌﻠﻤـﻲ‬
BADJI MOKHTAR- ANNABA UNIVERSITY ‫ﻋﻨﺎﺑ��ﺔ‬- ‫ﺟﺎﻣﻌـــــ��ـﺔ ﺑﺎﺟــــ��ـﻲ ﻣﺨﺘــ��ـﺎر‬

UNIVERSITÉ BADJI MOKHTAR- ANNABA
FACULTÉ DES SCIENCES DE L’INGÉNIORAT

DÉPARTEMENT DE GÉNIE MÉCANIQUE
MÉMOIRE
PRÉSENTÉ EN VUE DE L’OBTENTION DU DIPLÔME DE MASTER
INTITULÉ
EEeeeee
OPTIMISATION DES PERFORMANCES
ÉNERGÉTIQUES D’UNE INSTALLATION
CRYOGÉNIQUE DE LIQUÉFACTION DE
L’HÉLIUM
DOMAINE: SCIENCES ET TECHNOLOGIE
FILIЀRE : GÉNIE MÉCANIQUE
SPÉCIALITÉ : ÉNERGÉTIQUE
PRÉSENTÉ PAR :
LAOUAR IMED
DIRECTEUR DU MÉMOIRE :
HAOUAM ABDALLAH (MC-A) U.B.M. ANNABA
DEVANT LE JURY
PRÉSIDENT: MZAD H. (Pr) U.B.M. ANNABA
EXAMINATEUR: MECHIGHEL F. (MC-A) U.B.M. ANNABA
EXAMINATEUR: HAOUAM A. (MC-A) U.B.M. ANNABA
Année : 2018/2019
OPTIMISATION DES
PERFORMANCES ENERGETIQUES
D’UNE INSTALLATION
CRYOGENIQUE DE
LIQUEFACTION DE L’HELIUM
‫ﺑﺳم ﷲ اﻟرﺣﻣن اﻟرﺣﯾم‬
Au nom d’Allah le Tout Miséricordieux le Très Miséricordieux
Dédicaces
Je dédie ce mémoire à mes très

chers parents, à ma sœur et à mon
frère, à ma famille ainsi qu’à tous
les enseignants qui m’ont inculqué
éducation et savoir tout au long de
mon cursus scolaire et universitaire.
Remerciements
Je tiens à exprimer ma profonde gratitude en
remerciant Monsieur HAOUAM Abdallah pour avoir
accepté de m’encadrer. Un grand merci pour ses
précieux conseils, pour son aide ainsi que pour sa
gentillesse et sa grande disponibilité.
Je tiens à remercier également le Pr MZAD Hocine
Président du jury et le Dr MECHIGHEL Farid membre
du jury pour avoir examiné ce travail.
Résumé
Parmi tous les gaz existants, le plus intrigant et le plus intéressant à étudier est le gaz noble
connu sous le nom de « Hélium », lequel possède des propriétés qui le rendent unique. En
effet, l’hélium a la température d’ébullition la plus basse et est capable de se transformer en
un superfluide, contrairement aux autres gaz. Sa liquéfaction est un procédé ambigu au cours
duquel sa température doit être abaissée à 4,2 K, frôlant le zéro absolu. Le travail réalisé dans
ce mémoire comprend le calcul des caractéristiques thermodynamiques des différents cycles
associés à diverses installations cryogéniques destinées à la liquéfaction de l’hélium.
L’optimisation des performances thermiques nous a conduit à dimensionner une installation
de liquéfaction de l’hélium fonctionnant suivant un cycle combinant les processus de détente
isenthalpique et de détente isentropique avec le refroidissement préalable de l’hélium par
l’azote à l’état liquide saturé à 80 K et l’hydrogène liquide à 20 K. Enfin, une simulation des
différents cycles thermodynamiques par le logiciel HYSYS V7.1 a confirmé l’exactitude des
calculs effectués et une interprétation adéquate des résultats obtenus.
‫ﻣﻠﺧص‬
‫ ﻓﺈن أﻛﺛر اﻟﻐﺎزات إﺛﺎرة ﻟﻠﻔﺿول وإﺛﺎرة ﻟﻼھﺗﻣﺎم ھو اﻟﻐﺎز اﻟﻧﺑﯾل‬، ‫ﻣن ﺑﯾن ﺟﻣﯾﻊ اﻟﻐﺎزات اﻟﺣﺎﻟﯾﺔ‬
‫ اﻟﮭﯾﻠﯾوم ﻟدﯾﮫ أدﻧﻰ درﺟﺔ‬، ‫ ﻓﻲ اﻟواﻗﻊ‬.‫ ﻟدﯾﮫ ﺧﺻﺎﺋص ﺗﺟﻌﻠﮫ ﻓرﯾدا ﻣن ﻧوﻋﮫ‬."‫اﻟﻣﻌروف ﺑﺎﺳم "اﻟﮭﯾﻠﯾوم‬
‫ ﻓﺈن ﺗﺳﯾﯾل‬، ‫ وﻣﻊ ذﻟك‬.‫ ﻋﻠﻰ ﻋﻛس اﻟﻐﺎزات اﻷﺧرى‬، ‫ﺣرارة ﻏﻠﯾﺎن وﻗﺎدر ﻋﻠﻰ اﻟﺗﺣول إﻟﻰ ﺳﺎﺋل ﻓﺎﺋﻖ‬
.‫ ﺑﺎﻟﻘرب ﻣن اﻟﺻﻔر اﻟﻣطﻠﻖ‬K 4.2 ‫اﻟﮭﯾﻠﯾوم ﺗﻌﺗﺑر ﻋﻣﻠﯾﺔ ﺻﻌﺑﺔ ﺣﯾث ﯾﺟب ﺧﻔض درﺟﺔ ﺣرارﺗﮫ إﻟﻰ‬
‫ﯾﺷﻣل اﻟﻌﻣل اﻟﻣﻧﺟز ﻓﻲ ھذه اﻟرﺳﺎﻟﺔ ﺣﺳﺎب اﻟﺧﺻﺎﺋص اﻟدﯾﻧﺎﻣﯾﻛﯾﺔ اﻟﺣرارﯾﺔ ﻟﻠدورات اﻟﻣﺧﺗﻠﻔﺔ اﻟﻣرﺗﺑطﺔ‬
‫ ﻗﺎدﻧﺎ ﺗﺣﺳﯾن اﻷداء اﻟﺣراري إﻟﻰ ﺗﺻﻣﯾم ﻣﺣطﺔ ﺗﺳﯾﯾل اﻟﮭﯾﻠﯾوم ﺗﻌﻣل ﻋﻠﻰ‬.‫ﺑﺎﻟﺗرﻛﯾﺑﺎت اﻟﻣﺑردة اﻟﻣﺧﺗﻠﻔﺔ‬
‫ ﻣﻊ اﻟﺗﺑرﯾد اﻟﺳﺎﺑﻖ ﻟﻠﮭﯾﻠﯾوم ﺑﺎاﻟﻧﯾﺗروﺟﯾن‬isentropic ‫ وﻋﻣﻠﯾﺔ اﻟﺗﻣدد‬isenthalpic ‫دورة ﺗﺟﻣﻊ ﺑﯾن اﻟﺗﻣدد‬
‫ أﻛدت ﻣﺣﺎﻛﺎة اﻟدورات‬،‫أﺧﯾرا‬
ً . K 20 ‫ واﻟﮭﯾدروﺟﯾن اﻟﺳﺎﺋل ﻋﻧد‬K 80 ‫ﻓﻲ اﻟﺣﺎﻟﺔ اﻟﺳﺎﺋﻠﺔ اﻟﻣﺷﺑﻌﺔ ﻋﻧد‬
‫( دﻗﺔ اﻟﺣﺳﺎﺑﺎت اﻟﺗﻲ ﺗم إﺟراؤھﺎ‬HYSYS V7.1 ‫اﻟدﯾﻧﺎﻣﯾﻛﯾﺔ اﻟﺣرارﯾﺔ اﻟﻣﺧﺗﻠﻔﺔ )ﺑﺎﺳﺘﻌﻤﺎل ﺑرﻧﺎﻣﺞ‬
.‫وﺗﻔﺳﯾر وافٍ ﻟﻠﻧﺗﺎﺋﺞ اﻟﺗﻲ ﺗم اﻟﺣﺻول ﻋﻠﯾﮭﺎ‬
Abstract
Among all the existing gases, the most intriguing and interesting to study is the noble gas
known as the "Helium". It has properties that make it unique. Indeed, helium has the lowest
boiling temperature and is able to be transformed into super fluid, unlike other gases.
However, the liquefaction of helium is an ambiguous method in which its temperature must
be lowered to 4.2 K brushing absolute zero. The work done in this specification comprises
calculating the thermodynamic characteristics of the different liquefaction helium cycles
associated with various cryogenic installations. Optimization of thermal performance has led
us to design a cryogenic plant operating according to a cycle combining the isenthalpic
expansion process and the isentropic expansion of the pre-cooling helium with nitrogen in the
liquid saturated condition at 80 K and liquid hydrogen at 20 K. Finally, a simulation by
HYSYS V7.1, of different thermodynamic cycles confirmed the accuracy of the calculations
made and an adequate interpretation of the results obtained.
NOMENCLATURE
Abréviations utilisées
ES : Etranglement simple
RP : Refroidissement préalable
FR : Fluide réfrigérant
EC: Echangeur de chaleur
PMB : Point Mort Bas
PMH : Point Mort Haut
BP : Basse Pression
HP : Haute Pression
IRM : Imagerie par Résonance Magnétique
Symboles utilisés
x : Coefficient de liquéfaction
∆HT : Effet isothermique d’étranglement
q in : Imperfection de l’échangeur de chaleur
q a : Apport de chaleur
∆T: Différence de température
Cp : Chaleur spécifique à pression constante
Cv : Chaleur spécifique à volume constant
q 0 totale: Quantité de chaleur enlevée du gaz de travail par le liquide réfrigérant
q x supp : Puissance frigorifique supplémentaire nécessaire au refroidissement préalable
Go : Quantité du liquide réfrigérant
h: Enthalpie massique
T : Température
S : Entropie
Wsupp : Travail supplémentaire
r: Constante spécifique des gaz
ηk : Rendement du compresseur
Wk : Travail du compresseur
β: Coefficient correcteur
Wtotal : Travail total
Wo : Travail spécifique
q x : Puissance frigorifique du cycle
εréel : Coefficient frigorifique réel
εth : Coefficient frigorifique théorique
η: Rendement thermodynamique
LISTE DES FIGURES
Figure I.1 : Diagramme Pression-Température des phases de l’hélium

Figure II.1 : Diagramme Pression-Température de l’eau
Figure II.2 : Cycle idéal de liquéfaction
Figure II.3 : Transformations thermodynamiques inverses de liquéfaction d’un gaz
Figure II.4 : Schéma du procédé et cycle de Brayton inverse associé au processus
Figure II.5 : Schéma du cycle à étranglement simple de Linde
Figure II.6 : Diagramme du cycle à étranglement simple (E.S) de Linde
Figure II.7 : Schémas des cycles à étranglement E.S à un seul étage de R.P
Figure II.8 : Diagramme du cycle à étranglement avec un étage de R.P
Figure II.9 : Schéma du processus combiné avec branchement du détendeur au niveau initial
de température
Figure II.10 : Schéma du processus combiné avec branchement au niveau intermédiaire de
température
Figure II.11 : Schéma du processus combiné avec branchement au niveau inférieur de
température
Figure II.12 : Schéma du cycle combiné associé au processus avec branchement du
détendeur au niveau inférieur de température
Figure II.13 : Schéma du procédé de liquéfaction en cascade
Figure II.14 : Cycle associé au processus de liquéfaction en cascade
Figure III.1 : Schéma du processus de liquéfaction à étranglement simple et cycle associé au
processus à E.S
Figure III.2 : Schéma du processus de liquéfaction et cycle à étranglement simple avec un
étage de R.P
Figure III.3 : Schéma du processus de liquéfaction et cycle à étranglement avec deux étages
de R.P : azote liquide saturé à 80K et azote sous vide à 65K
de R.P : azote liquide saturé à 80K et hydrogène liquide à 20K
Figure III.5 : Cycle de liquéfaction de l’hélium en coordonnées T-S : R.P avec azote liquide
saturé à 80 K et hydrogène liquide à 20 K (P 2 =15atm).
Figure III.8 : Evolution de la consommation d’énergie spécifique en fonction de la pression
de refoulement
Figure III.9 : Evolution du rendement thermodynamique et du coefficient de liquéfaction de
l’hélium en fonction de la pression de refoulement
Figure III.10 : Schéma du processus combiné et du cycle associé au processus de
liquéfaction de l’hélium avec un branchement du détendeur à T=20K et un taux de détente
Ʈ=25
Figure III.11 : Représentation du cycle combiné en coordonnées T-S à 2 étages R.P avec
azote liquide saturé à 80 K et hydrogène liquide à 20 K utilisant la détente adiabatique et de
l’étranglement isenthalpique avec (P 2 = 25atm).
Figure IV.1 : Compresseur mono-vis
Figure IV.2 : Compresseur double-vis
Figure IV.3 : Liaison mécanique dans un compresseur à piston
Figure IV.4 : Echangeur à co-courants
Figure IV.5 : Echangeur à contre-courants
Figure IV.6 : Echangeur à courants croisés
Figure IV.7 : Echangeur à tubes et calandre
Figure IV.8 : Types d’échangeurs de chaleur à faisceaux tubulaires et calandre
Figure IV.9 : Vues diverses d’échangeurs de chaleur à faisceaux tubulaires et calandre
Figure IV.10 : Arrangement des tubes au niveau d’un échangeur de chaleur pour les pas carré
et triangulaire
Figure IV.11 : Variations de température du gaz de travail dans l’échangeur E1 du cycle à
E.S avec R.P
Figure IV.12 : Arrangement des tubes de l’échangeur tubulaire (pas carré)
Figure IV.13 : Echangeur tubulaire E1 du cycle hélium à étranglement simple avec R.P
Figure. V. 1 : Représentation des dix premières étapes de la simulation des cycles
cryogéniques
Figure. V. 2 : Installation du cycle à ES pour la liquéfaction de l’azote avec P 2 =200 bars
Figure. V. 3 : Installation du cycle à ES+RP pour la liquéfaction de l’azote avec P 2 =200bars
Figure. V.4 : Installation du cycle combiné pour la liquéfaction de l’azote avec un
branchement du détendeur au niveau initial de température et un taux de détente Ʈ=200
branchement du détendeur au niveau intermédiaire de température et un taux de détente
Ʈ=200
branchement du détendeur au niveau inférieur de température et un taux de détente Ʈ=200
Figure. V.7 : Installation du cycle à ES+RP pour la liquéfaction de l’hydrogène avec
P 2 =160atm
Figure. V.8 : Installation du cycle combiné pour la liquéfaction de l’hydrogène avec un
Ʈ=120
Figure. V. 9 : Installation du cycle à ES+RP pour la liquéfaction de l’hélium avec P 2 =15atm
Figure. V. 12 : Installation du cycle combiné pour la liquéfaction de l’hélium avec un
branchement du détendeur à T=65K et un taux de détente Ʈ=25
Figure. V. 15 : Allure des courbes de variation du flux de chaleur, de la différence de
température et du coefficient d’échange global en fonction de la température au niveau de
l’échangeur LNG-100 du cycle combiné avec branchement du détendeur à T=20K pour la
liquéfaction de l’hélium.
LISTE DES TABLEAUX

Tableau I.1 : Principaux constituants de l'air sec
Tableau I.2 : Température de saturation des principaux composants de l'air
Tableau I.3 : Propriétés physiques de l’air sec
Tableau I.4 : Propriétés physiques de l’azote
Tableau I.5 : Propriétés physiques de l’argon
Tableau I.6 : Propriétés physiques de l’oxygène
Tableau I.7 : Propriétés physiques de l’hydrogène
Tableau I.8 : Propriétés physiques de l’hélium
Tableau I.9 : Propriétés des isotopes stables de l’hélium
Tableau I.10 : Capacité de production de la centrale à oxygène de l’entreprise SIDER EL
HADJAR-Annaba (2018-2019)
Tableau III.1 : Points figuratifs du cycle à 2 étages de liquéfaction de l’hélium : R.P avec
azote liquide à 80 K et R.P avec hydrogène liquide à 20 K (P 2 =15atm).
azote liquide à 80 K et R.P avec hydrogène liquide à 20 K (P 2 =20atm)
Tableau III.4 : Tableau récapitulatif des résultats des calculs des caractéristiques
thermodynamiques de l’installation cryogénique à 2 étages R.P avec azote liquide à 80 K et
R.P avec hydrogène liquide à 20 K.
azote liquide à Patm à 80 K et R.P avec hydrogène liquide à 20 K utilisant la combinaison de
la détente adiabatique et de l’étranglement isenthalpique (P 2 =25atm)
Tableau III.6 : Caractéristiques thermodynamiques du cycle combiné pour la liquéfaction de
l’hélium
Tableau IV.1 : Caractéristiques du compresseur à pistons choisi
Tableau IV.2 : Propriétés thermiques du fluide chaud
Tableau IV.3 : Propriétés thermiques du fluide froid
Tableau IV.4 : Normes des tubes en fonction du diamètre de la calandre
Tableau IV.5 : Tableau récapitulatif des caractéristiques de l’écoulement des fluides chaud et
froid
Tableau V.1 : Points figuratifs du cycle à ES pour la liquéfaction de l’azote avec P 2 =200bars
Tableau V.2 : Points figuratifs du cycle à ES+RP pour la liquéfaction de l’azote avec
P 2 =200bars
Tableau V.3 : Points figuratifs du cycle combiné pour la liquéfaction de l’azote avec un
Ʈ=200
Tableau V.6 : Points figuratifs du cycle à ES+RP pour la liquéfaction de l’hydrogène avec
P 2 =160atm
Tableau V.7 : Points figuratifs du cycle combiné pour la liquéfaction de l’hydrogène avec un
Ʈ=120
Tableau V.8 : Points figuratifs du cycle à ES+RP pour la liquéfaction de l’hélium avec
P 2 =15atm
P 2 =20atm
P 2 =25atm
Tableau V.11 : Points figuratifs du cycle combiné pour la liquéfaction de l’hélium avec un
Tableau V.13: Liquéfaction de l’azote à P 2 =200 bars
Tableau V.14 : Liquéfaction de l’hydrogène à 160 et à 120 bars
Tableau V.15 : Comparaison des résultats calculés et simulés des cycles à E.S avec deux (02)
étages de refroidissement préalable utilisant de l’azote liquide saturé à T=80K et l’hydrogène
liquide à T=20K pour la liquéfaction de l’hélium
Tableau V.16 : : Comparaison des résultats calculés et simulés des cycles combinés pour la
liquéfaction de l’hélium avec deux étages de refroidissement préalable utilisant de l’azote
liquide saturé à T=80K et l’hydrogène liquide pour un branchement du détendeur à T=20K.
Sommaire
INTRODUCTION GENERALE ..................................................................................................................... 1
Chapitre I. GENERALITES SUR LES GAZ INDUSTRIELS .............................................................................. 2
I.1. Introduction ............................................................................................................................... 2
I.2. Les gaz industriels ...................................................................................................................... 2
I.3. Définitions et propriétés de quelques gaz permanents ............................................................ 2
I.3.1. L’air ..................................................................................................................................... 2
I.3.2. L’Azote ................................................................................................................................. 3
I.3.3. L’Argon ................................................................................................................................ 4
I.3.4. L’Oxygène ............................................................................................................................ 5
I.3.5. L’Hydrogène ........................................................................................................................ 5
I.3.6. L’Hélium .............................................................................................................................. 6
I.4. Localisation de l’hélium .............................................................................................................. 7
I.4.1 Provenance de l’hélium ........................................................................................................... 7
I.4.2. Utilisation et particularité de l’hélium ................................................................................. 7
I.4.3. Particularité de l’hélium ...................................................................................................... 9
I.5. Les principaux producteurs de gaz industriels [21-24] ........................................................... 10
I.5.1. Air Products....................................................................................................................... 10
I.5.2. Linde AG ............................................................................................................................ 10
I.5.3. Air Liquide ......................................................................................................................... 10
I.5.4. Praxair ............................................................................................................................... 10
I.5.5. L’entreprise nationale Sonatrach...................................................................................... 11
I.5.6. La centrale à oxygène de l’entreprise nationale SIDER EL HADJAR-Annaba .................... 11
Chapitre II. PROCEDES DE LIQUEFACTION DES GAZ ET CYCLES THERMODYNAMIQUES ASSOCIES ...... 12
II.1. Généralités sur la liquéfaction des gaz ....................................................................................... 12
II.2. Les cycles de base pour la liquéfaction des gaz ......................................................................... 13
II.2.1. Cycle idéal de liquéfaction................................................................................................... 13
II.2.2. Cycle théorique de Brayton inverse .................................................................................... 14
II.2.3. Cycle à étranglement simple ou cycle de LINDE ................................................................. 14
II.2.4. Le cycle à étranglement avec refroidissement préalable (E.S+R.P) .................................... 17
II.2.5. Le cycle combiné ou cycle de CLAUDE ................................................................................ 20
II.2.6. Cycle en cascade .................................................................................................................. 25
Chapitre III. ETUDE THERMODYNAMIQUE DES INSTALLATIONS ET CYCLES CRYOGENIQUES DE
LIQUEFACTION DE L’HELIUM................................................................................................................. 27
III.1. Cycle à étranglement simple de LINDE ou cycle E.S ............................................................... 27
III.2. Cycle à étranglement simple ES avec un étage de refroidissement préalable (R.P) à l’aide de
l’azote à 80K ...................................................................................................................................... 29
III.3 Cycle à étranglement simple avec deux étages de refroidissement préalable utilisant de l’azote
à 80K et de l’azote sous vide à 65 K .................................................................................................. 31
III.4. Cycle à étranglement simple avec deux étages de R.P fonctionnant avec de l’azote à 80 K et
de l’hydrogène à 20 K ........................................................................................................................ 33
III.4.1. Calcul des caractéristiques thermodynamiques de l’installation cryogénique avec
P 2 =1MPa ........................................................................................................................................ 34
III.4.2. Calcul des caractéristiques thermodynamiques de l’installation cryogénique avec
P 2 =1.5MPa ..................................................................................................................................... 34
III.4.3. Calcul des paramètres de l’installation avec P 2 =2MPa .................................................. 39
III.4.4. Calcul des paramètres de l’installation avec P 2 =2.5MPa ................................................ 43
III.4.5. Représentation graphique de l’évolution des caractéristiques thermodynamiques de
l’installation cryogénique à 2 étages ............................................................................................. 47
III.5. Cycle combiné pour la liquéfaction de l’hélium avec deux étages de refroidissement
préalable R.P utilisant de l’azote à 80K et de l’hydrogène à 20 K avec branchement du détendeur à
T=20K ................................................................................................................................................. 48
III.5.1. Calcul des caractéristiques thermodynamiques de l’installation cryogénique avec P 2 =25
atm ................................................................................................................................................ 49
III.6. Comparaison du cycle combiné avec le cycle à étranglement simple à deux étages de R.P
fonctionnant avec de l’azote liquide à 80 K et de l’hydrogène liquide à 20 K à la pression P 2 =25atm
........................................................................................................................................................... 54
Chapitre IV. DIMENSIONNEMENT DES PRINCIPAUX ORGANES DE L’INSTALLATION CRYOGENIQUE DE
LIQUEFACTION DE L’HELIUM................................................................................................................. 55
IV.1 Introduction ............................................................................................................................... 55
IV.2 Les compresseurs ..................................................................................................................... 55
IV.2.1 Les différents types de compresseurs................................................................................. 55
IV.2.2 Fonctionnement des compresseurs volumétriques à pistons ......................................... 55
IV.2.3 Fonctionnement des compresseurs à vis ......................................................................... 56
IV.2.4. Principaux organes des compresseurs à pistons ............................................................. 58
IV.3 Calcul thermodynamique des caractéristiques du compresseur à pistons de l’installation de
liquéfaction de l’hélium..................................................................................................................... 59
IV.3.1 Choix du taux de compression optimal ............................................................................... 59
IV.3.2. Calcul du débit massique ................................................................................................. 60
IV.3.3. Calcul du volume de gaz aspiré par le compresseur ...................................................... 60
IV.3.4. Calcul du volume de gaz balayé par le compresseur....................................................... 60
IV.3.5. Calcul de la cylindrée ....................................................................................................... 60
IV.3.6. Choix des dimensions des pistons du compresseur ........................................................ 60
IV.3.6 Vérification des dimensions choisies................................................................................ 61
IV.4. Dimensionnement de l’échangeur de chaleur ....................................................................... 61
IV.4.1. Définition ......................................................................................................................... 61
IV.4.2. Classification des échangeurs de chaleur ........................................................................ 61
IV.4.3. Les échangeurs à faisceaux tubulaires et calandre ......................................................... 63
IV.5 Calcul thermique de l’échangeur tubulaire E1 du cycle hélium à étranglement simple avec
R.P...................................................................................................................................................... 66
IV.5.1. Calcul du flux thermique ................................................................................................. 66
IV.5.2. Propriétés des fluides chaud et froid .............................................................................. 66
IV.5.3. Calcul du rapport des taux de capacités thermiques ...................................................... 67
IV.5.4. Calcul du flux thermique maximal ................................................................................... 67
IV.5.5. Calcul de l’efficacité de l’échangeur ................................................................................ 68
IV.5.6. Calcul de la différence de température logarithmique moyenne 𝚫𝚫𝚫𝚫𝚫𝚫𝚫𝚫 ...................... 68 17T
IV.5.7. Correction de la différence de température logarithmique moyenne 𝚫𝚫𝚫𝚫𝚫𝚫𝚫𝚫 ............... 68 17T
IV.5.8. Calcul du coefficient d’échange interne (coté tubes)...................................................... 69

IV.5.9. Calcul du coefficient d’échange externe (coté calandre) ................................................ 70
IV.5.10. Calcul du coefficient d’échange global .......................................................................... 71
IV.5.11. Calcul de la surface d’échange ...................................................................................... 72
IV.5.12. Calcul de la longueur des tubes ..................................................................................... 72
IV.6 Rapport entre l’efficacité de l’échangeur et le coefficient de liquéfaction de l’installation
cryogénique ....................................................................................................................................... 73
IV.6.1 Cycle à étranglement simple avec R.P ................................................................................ 73
IV.6.2 cycle combiné...................................................................................................................... 74
Chapitre V. SIMULATION DES CYCLES CRYOGENIQUES DE LA LIQUEFACTION DE L’HELIUM EN
UTILISANT LE LOGICIEL HYSYS V7.1 ....................................................................................................... 76
V.1. Etapes de simulation ............................................................................................................... 76
V.2. Liquéfaction de l’azote T=77K ................................................................................................. 79
V.2.1 Cycle à étranglement simple pour la liquéfaction de l’azote ............................................ 79
V.2.2. Cycle à étranglement simple avec refroidissement préalable utilisant un fluide
frigorigène (à T=230K) pour la liquéfaction de l’azote .................................................................. 80
V.2.3. Cycle combiné pour la liquéfaction de l’azote avec branchement du détendeur au niveau
initial à T=300K .............................................................................................................................. 82
V.2.4. Cycle combiné pour la liquéfaction de l’azote avec refroidissement préalable pour un
branchement du détendeur au niveau intermédiaire à T=230K ................................................... 84
V.2.5. Cycle combiné pour la liquéfaction de l’azote avec refroidissement préalable à 230K pour
un branchement du détendeur à un niveau intermédiaire à T=160K........................................... 86
V.3. Liquéfaction de l’hydrogène T=20K ......................................................................................... 87
V.3.1. Cycle à étranglement simple avec deux étages de refroidissement préalable (utilisant de
l’azote liquide saturé à la pression atmosphérique à 80K et de l’azote sous vide 65K pour la
liquéfaction de l’hydrogène .......................................................................................................... 88
V.3.2. Cycle combiné pour la liquéfaction de l’hydrogène avec un étage de refroidissement
préalable utilisant de l’azote avec un branchement du détendeur à T=80K ................................ 90
V.4. Cycles cryogéniques pour la liquéfaction de l’hélium ............................................................. 92
V.4.1. Cycle à étranglement simple avec deux étages de refroidissement préalable utilisant de
l’azote liquide saturé à 80K et de l’hydrogène liquide à 20K pour la liquéfaction de l’hélium..... 92
V.4.2. Cycle combiné pour la liquéfaction de l’hélium avec un étage de refroidissement
préalable utilisant de l’azote liquide pour un branchement du détendeur à T=65K .................... 98
V.4.3. Cycle combiné pour la liquéfaction de l’hélium avec deux étages de refroidissement
préalable utilisant de l’azote 80K et de l’azote sous vide 65K pour un branchement du
détendeur à T=45K ........................................................................................................................ 99
V.4.4. Cycle combiné pour la liquéfaction de l’hélium avec deux étages de refroidissement
préalable utilisant de l’azote à 80K et de l’hydrogène pour un branchement du détendeur à
T=20K ........................................................................................................................................... 101
V.5. Interprétation des résultats de la simulation des différents cycles cryogéniques de
liquéfaction des gaz utilisés : azote, hydrogène et hélium ......................................................... 104
CONCLUSION GENERALE ..................................................................................................................... 106
REFERENCES ........................................................................................................................................ 108
ANNEXES.............................................................................................................................................. 111
Introduction générale
INTRODUCTION GENERALE
L’étude thermodynamique des gaz industriels : air, azote, oxygène, hydrogène, argon,
méthane et hélium est un domaine très passionnant aussi bien en sciences physiques qu’en
énergétique, car ils peuvent être utilisés sous différentes formes (gazeuses et liquides).
Certains gaz sont employés comme fluides frigorigènes ou gaz de travail dans les installations
cryogéniques.
Ils sont ainsi utilisés dans différents domaines : industriel, médical et même au quotidien…ils
possèdent des propriétés physiques intéressantes. Parmi ces gaz, l’Hélium, exploité aussi bien
à l’état liquide que gazeux ; il est connu comme étant un gaz noble dont les multiples
caractéristiques sont très remarquables, suscitant ainsi l’intérêt de nombreux chercheurs.
A l’état liquide, l’hélium est considéré comme étant la matière la plus froide sur terre. Il peut
également adopter une seconde forme d’état liquide et dans ce cas est appelé « Superfluide ».
Cependant, sa liquéfaction demeure un véritable défi du fait que sa température doit avoisiner
le « zéro absolu ».
Ce mémoire a pour objectif le calcul thermodynamique des cycles cryogéniques permettant de

choisir les installations les plus performantes, nécessaires à la liquéfaction de l’Hélium. Ce
travail est structuré en cinq (05) chapitres comme suit :
- Chapitre I : Généralités sur les gaz industriels

- Chapitre II : Procédés de liquéfaction des gaz et cycles thermodynamiques associés
- Chapitre III : Etude thermodynamique des installations et cycles cryogéniques de
liquéfaction de l’hélium
- Chapitre IV : Dimensionnement des principaux organes de l’installation cryogénique
de liquéfaction de l’hélium.
- Chapitre V : Simulation des cycles des installations cryogéniques de liquéfaction de
l’hélium.
Enfin, le mémoire s’achève par une conclusion générale dans laquelle le travail réalisé est
synthétisé, faisant ressortir les principaux résultats des étapes franchies pour la détermination
d’une installation cryogénique aux performances optimales pour la liquéfaction de l’hélium.
1
Chapitre I. Généralités sur les gaz industriels
Chapitre I. GENERALITES SUR LES GAZ INDUSTRIELS
I.1. Introduction
Ce chapitre sera consacré essentiellement à l’étude des propriétés physiques des gaz
industriels avec une attention particulière à l’hélium, fluide à liquéfier au sein d’une
installation cryogénique qui sera déterminée ultérieurement.
I.2. Les gaz industriels

Les gaz industriels connus sous le nom de « gaz atmosphériques » sont des gaz qui composent
l’air que nous respirons. Utilisés dans plusieurs domaines, ils ont des propriétés différentes
les unes par rapport aux autres, ils sont très utiles notamment dans le domaine industriel, dans
notre vie quotidienne mais aussi dans le domaine médicale. [1-4]
I.3. Définitions et propriétés de quelques gaz permanents

I.3.1. L’air [5-8]
Tableau I.1 : Principaux constituants de l'air sec
Gaz Pourcentage % en
constituants volume
L’air que nous respirons représente, avec la terre,
l’air sec
le feu et l'eau, l'un des éléments constitutifs du
Azote 78,09
monde. C’est un mélange de gaz, d’humidité et de
Dioxygène 20,95
particules polluantes (Bactéries, poussières…etc.).
Argon 0,93
Dioxyde de 0,035 Constituant l'atmosphère de la Terre, il est incolore
carbone (invisible) et inodore. Cet air peut être dépourvu
Néon 1,8 10-3 d’humidité et d’impuretés, on parle donc d’air sec
Hélium 5,24 10-4 filtré qui pourra être utilisé dans plusieurs
Krypton 1,0 10-4 domaines (agent de pression pour les outils
Hydrogène 5,0 10-4 pneumatiques, dans les moteurs à piston) ou bien
Xénon 8,0 10-6 utiliser ces composant en les séparent. Les
Ozone 1,0 10-6 températures de saturation des principaux
Radon 6,0 10-18 composants de l’air sont représentées sur le tableau
I.2.
Tableau I.2 : Température de saturation des principaux composants de l'air
Gaz N2 O2 Ar Ne He H2
Température -195,79 -183 -185,85 -246,05 -268,93 -252,76
d'ébullition
(°C)
2
Tableau I.3 : Propriétés physiques de l’air sec
Masse molaire (g/mol) 28,9653

Masse volumique (kg/m3) 1,293
Température de fusion (K) 56.95=−216,2°C
Température d’ébullitions (K) 78.85=−194,3 °C
Température critique (°C) -140.73
Pression critique (bar) 37.74
Chaleur latente de vaporisation (air) : L v = 204.1 kJ/kg
Conductivité thermique (air) : 0.0234 W/m. K
Les principaux constituants de l’air (l’azote, l’oxygène et l’argon) peuvent êtres séparés par
liquéfaction en une succession d’étapes. [9]
Etape 1 : Préparation de l’air
• Aspiration et filtration de l’air ;

• Compression de l’air atmosphérique jusqu'à P2= 5 à 6 bars ;
• Refroidissement de l’air après compression dans des tours de refroidissement et de
lavage ;
• Elimination de toute trace d’hydrocarbures et d’humidité dans des tamis moléculaires ;
Etape 2 : refroidissement cryogénique et séparation des gaz de l’air
• Passage de l’air qui a été préparé dans des échangeurs de chaleur et des turbines afin
d’obtenir des températures cryogéniques de l’ordre de -130°C ;
• L’air et les gaz refroidis seront liquéfiés et séparés dans des colonnes dites : colonnes
de séparation ;
• L’oxygène liquéfié est extrait à -183°C, L’argon est extrait à -186°C, L’azote est
extrait à -196°C ;
• Stockages et/ou distribution;
I.3.2. L’Azote [10]
L'azote est l'élément chimique de numéro atomique 7 découvert en 1772 par DANIEL
Rutherford, de symbole N (du latin nitrogenium). C’est le composant principale de l’air il est
incolore, inodore et ininflammable. Il est plus léger que l’eau quand il est à l’état liquide.
Utilisé essentiellement pour créer une atmosphère inerte pour les processus chimiques.
3
Tableau I.4 : Propriétés physiques de l’azote
Numéro atomique 7
Masse molaire (g/mol) 14.0067
Masse volumique (kg/m3) 0.001251
Température de fusion (K) 63.25
Température d’ébullitions (K) 77.35
Structure électronique K2L5
Configuration électronique 1s22s22p3
Rayon atomique 92pm
Charge de l’ion N3-
Température point triple (°C) -210.00
Pression point triple (bar) 0.125
I.3.3. L’Argon [11]
Utilisé essentiellement pour créer une atmosphère inerte (autour des filaments dans les
ampoules, comme isolant dans les vitrages...), agent protecteur pour le soudage, également les
lasers à argon utilisée en chirurgie pour les problèmes oculaires, aussi pour gonfler les
combinaisons de plongée sous-marine, enfin comme agent extincteur ainsi qu’en technologie
alimentaire …
Tableau I.5 : Propriétés physiques de l’argon
Numéro atomique 18
2 8 8
Structure électronique KLM
Configuration électronique 1s 2s 2p63s23p6
2 2
Rayon atomique 71pm

Température critique (°C) -122
4
I.3.4. L’Oxygène [12]
Découvert en 1774 par Carl Wilhelm Scheele, c’est l’élément que l’on trouve en abondance
sur terre. Utilisé dans plusieurs domaines, en médecine (équipements respiratoires, fabrication
de produits pharmaceutiques), en industrie (combustion dans les hauts fourneaux et
élimination d’impuretés dans les aciéries), dans la vie quotidienne (notamment pour le
chauffage et dans d’autres domaines).
Avec son effet oxydant, respectivement comburant, il doit être utilisé avec précautions et doit
être stocké séparément des matières combustibles.
Tableau I.6 : Propriétés physiques de l’oxygène
Numéro atomique 8
Structure électronique K2L6
Configuration électronique 1s22s22p3
Rayon atomique 48pm
Charge de l’ion O2-
I.3.5. L’Hydrogène [13-15]
L’hydrogène est un gaz monoatomique découvert en 1766 par Henry Cavendish, le plus léger
de symbole H, son noyau contient un seul proton et donc il gravite autour de lui qu’un seul
électron. L’hydrogène se présente sous forme de molécules H 2 (dihydrogène), on le trouve en
abondance dans tout l’univers, il est lié à l’oxygène pour former l’H 2 O (vapeur). Dans les
molécules d’eau sous forme d’ion H+(Cation), on peut donc l’extraire par électrolyse. Comme
il se trouve dans les hydrocarbures, Il peut être aussi extrait par reformage du méthane CH 4 .
Les isotopes de l’hydrogène sont le protium 1H (en abondance), le deutérium 2H (2% dans la
nature) et le tritium 3H (de petite quantité).
L’hydrogène est utilisé dans les piles à combustible pour produire de l’électricité (moteur
électrique), pour la combustion (fusées) Il occupe un très grand volume à l’état gazeux, ce qui
n’est pas très pratique. On le comprime pour diminuer l’espace qu’il occupe, il est liquéfié à
T=20.25 K.
5
Tableau I.7 : Propriétés physiques de l’hydrogène
Numéro atomique 1
Structure électronique K1
Configuration électronique 1s1
Rayon atomique 53pm
Charge de l’ion H+, H-
Caractéristique de sa phase gazeuse :
Masse volumique : 0,085 kg/m3
Conductivité thermique : 0.16835 W/m. K
Caractéristique de sa phase liquide :
Masse volumique : 70.97 kg/m3
Chaleur latente de vaporisation : 454,3 kJ/kg
Caractéristique de sa phase solide :
Chaleur latente de fusion: 58,158 kJ/kg
I.3.6. L’Hélium [16-17]
Découvert en 1895 par Pierre Jules César Jansen et Joseph Norman Lockyer,
L’Hélium possède des propriétés très intéressantes et uniques parmi tous les éléments. De
symbole « He » c’est le gaz noble le plus répandu dans l’univers observable. Il est inodore,
incolore et sera le 2ème après l’hydrogène en terme d’abondance, Il est situé dans le coin
supérieur droit du tableau périodique de Mendeleïev. Son numéro atomique est égal à 2.
C’est un gaz inerte utilisé en raison de sa légèreté (par rapport à l’air), du fait que sa masse
molaire est presque 7 fois plus petite que celle de l’air.
M (He) = 4 g/mole
M (Air) = 29 g/mole
6
L’hélium flotte facilement au dessus de l’air, par conséquent la vitesse du son dans l’hélium
est supérieure à celle de l’air.
Par exemple, la vitesse du son dans un air sec (M=28.964 g/mol) à zéro degré Celsius est
de 331.3 m/s. À cette température, pour l'hélium (M=4.003 g/mol) la vitesse du son est de
891.2 m/s.
Tableau I.8 : Propriétés physiques de l’hélium
Numéro atomique 2
Structure électronique K2
Configuration électronique 1s2
Rayon atomique 31pm
I.4. Localisation de l’hélium [16-20]

I.4.1 Provenance de l’hélium
L’hélium n’est pas aussi abondant sur terre : l’atmosphère n’en contient que 0,0005%.
Il peut être retrouvé dans les gisements de gaz naturel en raison de la désintégration
radioactive de l’Uranium et du Thorium dans la croûte terrestre, dans ce cas l’hélium peut
s’accumuler à l’intérieur de la terre avec le gaz naturel où il est extrait par processus
cryogénique, environ 80 à 90% de la production mondiale d'hélium provient de cette
technique.
La Tanzanie pourrait devenir un nouvel acteur du marché de l'hélium, après la récente

découverte d'un énorme gisement naturel de ce gaz noble estimé à 1,5 milliard de m3, soit
l'équivalent de la demande mondiale pendant sept (07) ans.
Autres origines :
• La fusion des isotopes d’hydrogène.
• La transformation spontanée des noyaux d’atomes instables (radioactivité alpha).
• La plupart de l’hélium dans l’univers se trouve dans les étoiles et dans les planètes.
I.4.2. Utilisation et particularité de l’hélium

L’hélium est utilisé sous ses deux formes : liquide et gazeuse. Sous forme de gaz, l’hélium est
principalement utilisé dans l’industrie électronique.
L’hélium est un gaz neutre, c’est cette caractéristique qui intéresse les fabricants.
Ininflammable, il ne réagit avec aucun autre corps. C’est une sorte de gaz protecteur qui
permet de créer une atmosphère neutre dans le processus de fabrication, lors de soudures.
7
À l’état liquide, l’hélium a des propriétés stupéfiantes Il est utilisé pour atteindre de très
basses températures (se liquéfie à -269 °C). Il est ainsi un élément essentiel dans l’imagerie
médicale.
Sans hélium, il n’y aurait pas autant d’IRM, surtout que la demande dans ce secteur (comme
dans celui de l’électronique), ne fait que croître.
Il sert par ailleurs dans le domaine spatial, pour refroidir les lanceurs de satellites, et dans
l'automobile pour la production d'airbags.
Les deux isotopes de l’hélium les plus connus sont :

4 He (2 protons et 2 neutrons) qui est le plus abondant sur terre. Utilisé en :
• En médecine pour scanner les corps humains.

• Gaz de flottabilité pour gonfler les ballons et les dirigeables. (Pour des raisons de
sécurité, 20% d’Oxygène y sont ajoutés : en effet, respirer une grande quantité de ce gaz
peut être dangereux. De plus, l’addition de l’oxygène permet de rendre les ballons plus
sûrs.
• L’hélium liquide peut être utilisé pour refroidir les supraconducteurs, qui à son tour est
refroidi avec l’azote liquide (moins cher). Ainsi, il est utilisé dans des études
scientifiques.
• Sa luminosité dépend de la tension d’électricité et du diamètre de l’ampoule.
Contrairement au Néon et au Xénon, l’Hélium peur briller à distance du fil du générateur.
• Sous forme de plasma, l’hélium peut former des composés extrêmement instables avec
l’H 2 et ils forment même des substances spéciales comme des réseaux piégeant des
molécules.
• Mélange hélium-oxygène comme gaz respiratoire (médecine et plongée).
3 He (2 protons et 1 neutron). Cet isotope est parfait pour mener des réactions de fusion
avec les isotopes de l’Hydrogène mais malheureusement il est très rare sur terre car il
disparaît immédiatement dans l’atmosphère (il est présent dans la lune où il n’y a pas
d’atmosphère). Ce processus de fusion peut aider l’humanité à remplacer le combustible
fossile.
Tableau I.9 : Propriété des isotopes stables de l’hélium

4 3
Propriétés He He
Température d’ébullition (K) 4.2 3.2
Température critique (K) 5.2 3.32
Pression critique (bar) 2.23 1.16
Enthalpie d’évaporation massique à T éb (kJ/kg) 21 8,44
Masse volumique du liquide à Téb et à Patm (kg/m3 ) 125 58
Viscosité du gaz à Téb (pa.s).10-6 1,241 1,0
Conductivité thermique du gaz (300K, 1atm) (W/m. K). 0,152 *
Conductivité thermique du liquide (4,2K, 1atm) (W/m. K). 0,018 *
Chaleur latente de vaporisation (kJ/kg) 20,42 *
8
Enfin l’hélium possède des propriétés très intéressantes et uniques parmi tous les éléments. A
basse température (T= 4,2 K), l’hélium devient le liquide le plus léger dont la densité est près
de 10 fois inférieure à celle de l’eau (d= 0,125 kg/m3) et pour maintenir une température aussi
basse, l’hélium liquide est stocké dans des DEWARS cryogéniques qui sont également
refroidis à l’extérieur avec de l’azote.
A noter qu’il existe une 2ème forme de l’hélium liquide obtenu à T = 2,17 K.
Figure I.1 : Diagramme Pression-Température des phases de l’hélium
Si l’hélium est refroidi davantage (en réduisant la pression dans les vases de stockage) le
liquide finira par atteindre le point lambda λ (T = 2,17 K), d’où la 2ème forme liquide de
l’hélium.
• Le liquide cesse de bouillir et commence à présenter des propriétés complètement

différentes, à de telles températures, la conductivité de l’hélium liquide augmente si
bien qu’elle devient ainsi plus élevée que celle du Cuivre (Cu) et de l’Argent (Ag)
(C’est pourquoi le liquide ne bout pas, car la chaleur se répand immédiatement et
uniformément dans tout le matériau).
• L’Hélium devient un superfluide. Cette superfluidité de l’Hélium liquide à été la 1ère
découverte en 1958.
I.4.3. Particularité de l’hélium

Un élément aussi simple qui à l’état liquide est la matière la plus froide sur terre, possède des
propriétés surprenantes et est irremplaçable. Il peut avoir des propriétés uniques de
superfluidité à des températures extrêmement basses.
9
L’hélium est le seul élément à ne pas posséder un point triple ce qui le rend très spécial.
I.5. Les principaux producteurs de gaz industriels [21-24]

Les acteurs prédominants du marché mondial du gaz industriel sont : Air Products, Linde, Air
liquide et Praxair.
I.5.1. Air Products

Air Products and Chemicals, est un groupe industriel américain qui a été fondé en 1940,
spécialiste des gaz industriels et médicaux. Son siège se trouve aux États-Unis.
Air Products possède des employés issus de divers horizons créant des solutions innovantes
qui profitent à l’environnement, améliorent la durabilité et relèvent les défis auxquels sont
confrontés les clients, les communautés et le monde.
Ils produisent et livrent des gaz industriels et médicaux, leurs effectifs sont de 20000 avec un
chiffre d’affaire de 10,180 milliards de dollars. Il est l'un des principaux fournisseurs
mondiaux de gaz industriels depuis plus de 75 ans. Son cœur de métier, les gaz industriels,
fournit des gaz atmosphériques et de procédé ainsi que des équipements connexes aux
marchés de fabrication dans plusieurs secteurs.
I.5.2. Linde AG
Linde est un groupe international allemand créé en 1878 par Carl Von Linde et ses associés
connu par son nom officiel Linde AG. Spécialisé dans les domaines Gaz et Ingénierie, dont le
siège est à Munich (en Allemagne). Linde est souvent mis en parallèle avec Air liquide avec
un effectif de 63487 et un chiffre d’affaire de 16,655 milliards de dollars, parce qu'ils sont
respectivement le 2ème et le 1er acteur du marché des gaz industriels et sont tous deux basés en
Europe.
Cependant le duo Linde-Praxair de la même famille compte 80 000 salariés, et affiche un chiffre
d’affaires de 24 milliards d’euros, voulant occuper la 1ère place et détrôner Air Liquide avec un chiffre
d’affaires de 20 milliards d’euros.
I.5.3. Air Liquide

Est un groupe industriel français créé en 1902, spécialiste des gaz industriels, la santé,
l'environnement et la recherche. Il est présent dans quatre-vingts (80) pays à travers le monde.
Premier groupe mondial dans le secteur avec un chiffre d’affaire qui dépasse 20 milliards de
dollars avec un effectif de 68000 , son principal concurrent reste Linde, puis viennent
Praxair et Air Products.
I.5.4. Praxair
Créée en 1907 par son fondateur Carl Von Linde qui produit et distribue des gaz purs, utilisés
dans les processus de diverses industries. Praxair est présente dans 30 pays et possède un
chiffre d’affaire similaire à Air Products, depuis près d'une décennie.
10
I.5.5. L’entreprise nationale Sonatrach
Le groupe américain leader mondial dans les gaz industriels (Air Products) et le géant
algérien de l'énergie ont conclu deux accords fin 2018, pour un montant global estimé à 100
millions de dollars. Les deux contrats, concerneront la coentreprise algérienne Helios
(Arzew), détenue à 51 % par Sonatrach et à 49 % par Air Products. Pour la production de
l’hélium en Algérie et la construction de deux unités de séparation des gaz de l'air.
Air Products sera en charge du design et de la construction, Helios sera par la suite
propriétaire et en charge des opérations. Ces projets de production des gaz industriels avec un
partenaire de renommé, permettront de mieux valoriser nos ressources et sécuriser nos besoins
en gaz industriels et ainsi renforcer la position de l’Algérie comme un producteur de référence
à un niveau international, et d’offrir sur le marché algérien des gaz à des prix compétitifs tout
en réduisant les importations.
I.5.6. La centrale à oxygène de l’entreprise nationale SIDER EL HADJAR-Annaba
La centrale à oxygène est considérée comme étant le poumon du complexe sidérurgique

SIDER [EL HADJAR], elle se charge de produire et de stocker l’azote, l’oxygène et l’argon
sous forme gazeuse et liquide. Ces gaz sont par la suite distribués aux unités de production
pour être utilisé majoritairement dans les processus de fabrication de l’acier ou dans le but
d’être commercialisé.
Tableau I.10 : Capacité de production de la centrale à oxygène de l’entreprise SIDER EL

HADJAR-Annaba (2018-2019)
Gaz/Liquide Production stockage

Oxygène gaz (m3/Jours) 250999 180000
Azote gaz (m3/Jours) 191673 108000
Oxygène liquide (l/H) 625 -
Azote liquide (l/H) 417 -
Argon liquide (l/H) 150 -
11
Chapitre II. Procédés de liquéfaction des gaz et cycles thermodynamiques associes
Chapitre II. PROCEDES DE LIQUEFACTION DES GAZ ET CYCLES

THERMODYNAMIQUES ASSOCIES
II.1. Généralités sur la liquéfaction des gaz

Toute substance peut se trouver alternativement dans chacun de ces états (solide-liquide-gaz)
selon les niveaux de température et les pressions auxquelles elle est soumise (Fig. II.1).
Figure II.1: Diagramme Pression-Température de l’eau
La liquéfaction est un changement de la phase gazeuse vers la phase liquide, elle est réalisée
par de divers procédés cryogéniques. La façon de liquéfier un gaz qui a une température
critique supérieure à la température ambiante est d'augmenter sa pression jusqu'à atteindre les
conditions de saturation à l’état liquide. Pour de nombreux gaz dont Tcr < T0 température
ambiante, aucune condensation ne sera possible sans abaisser sa température en dessous de
cette valeur même si l’on augmente la pression du gaz de travail à des pressions infinies. Pour
ces gaz, il faudrait donc faire appel à des procédés industriels d’abaissement de la température
par une utilisation simultanée de la compression et de la détente [25-26].
Les cycles de réfrigération et de liquéfaction cryogéniques mettent en jeu des combinaisons de

compressions isothermes, de refroidissements, et de détentes isenthalpiques et adiabatiques des
fluides. On distingue quatre (04) grandes familles de procédés thermodynamiques cryogéniques :
• Les procédés à détente isenthalpique de Joule-Thomson.

• Les cycles inverses de Brayton à détente isentropique.
• Les procédés mixtes associant une détente isenthalpique et une détente isentropique
(cycle de Claude).
• Les cycles en cascades.
12
II.2. Les cycles de base pour la liquéfaction des gaz
II.2.1. Cycle idéal de liquéfaction
C’est un cycle parfait pour transformer une quantité de gaz entièrement en liquide (Fig. II.2).
Figure II.2 : Cycle idéal pour la liquéfaction d’un gaz
Evolutions thermodynamiques du cycle idéal

1-2 : Compression isotherme
2-C : Détente isentropique
C-3 : Evaporation isotherme et isobare
3-1 : Chauffage isobare
La liquéfaction peut être imaginée en inversant le sens des transformations thermodynamiques
du cycle (Fig. II.3).
Figure II.3 : Transformations thermodynamiques inverses de liquéfaction d’un gaz
Evolutions thermodynamiques inverses :

1-3 : Refroidissement isobare
3-c : Condensation isotherme et isobare
13
II.2.2. Cycle théorique de Brayton inverse
Le cycle de Brayton inverse (Fig. II.4) réalise un effet frigorifique tout en produisant un
travail. Un gaz est comprimé (1-2), refroidi (2-3), détendu (3-4) et réchauffé à son retour
dans l’échangeur (4-1). La température à la fin de la détente étant basse, ce gaz peut être
utilisé pour refroidir une enceinte, soit par contact direct, soit par l'intermédiaire d'un
échangeur [26-29].
1
C
P1 P2
E
4 D
Figure II.4 : Schéma du procédé et cycle de Brayton inverse associé au processus
II.2.3. Cycle à étranglement simple ou cycle de LINDE
Ce cycle réel a été inventé par l’allemand LINDE vers la fin du 17ém siècle. Il repose sur le
P
principe de faire chuter la température d’un gaz grâce à un échangeur de chaleur et à une
valve d’étranglement après sa compression isotherme (Fig. II.5).
14
qk
c
2 1`
q in1
q a1
3
VE
5
4
Figure II.5 : Schéma du cycle à étranglement simple de Linde
L’étranglement du gaz à travers la valve provoquera sa détente isenthalpique, il s’accompagne

d’un refroidissement brusque c’est l’effet de Joule Thomson.
15
Figure II.6: Diagramme du cycle à étranglement simple (E.S) de Linde
Evolutions thermodynamiques du cycle à E.S
1-2 : Compression isotherme du gaz de travail.
2-3 : Refroidissement isobare.
3-4 : Détente isenthalpique.
5-1 : Réchauffement isobare d’une quantité (1-x) du gaz de travail.
𝐪𝐪𝐢𝐢𝐢𝐢 : Pertes de l’échangeur.
𝐪𝐪𝐚𝐚 : Apport de chaleur du milieu extérieur.
16
∆𝐇𝐇𝐓𝐓 : Coefficient d’étranglement isothermique.

𝐓𝐓𝐱𝐱 : Température de condensation du gaz de travail.
La quantité de liquide produite sera caractérisée par un coefficient de liquéfaction :
∆𝐇𝐇𝐓𝐓 − (𝐪𝐪𝐢𝐢𝐢𝐢 + 𝐪𝐪𝐚𝐚 )

𝐱𝐱 =
𝐇𝐇𝟏𝟏 − 𝐇𝐇𝐟𝐟 − 𝐪𝐪𝐢𝐢𝐢𝐢
Rendement du cycle à E.S
• Pour la liquéfaction de l’air (Téb = −194,3°C), le rendement de ce cycle est égal à 5%

• Pour la liquéfaction de l’hélium(Téb = −268,93°C), ce cycle est irréalisable.
II.2.4. Le cycle à étranglement avec refroidissement préalable (E.S+R.P)
Le refroidissement préalable s’effectue par un ou des cycles extérieurs (un ou plusieurs étages
de R.P), pour faire circuler un ou des fluide(s) réfrigérant(s) dans l’installation. Le liquide
réfrigérant rentre en contact avec le gaz de travail sans être mélangé par l’intermédiaire d’un
échangeur de chaleur (E.C).
Le fluide réfrigérant fera chuter la température du gaz de travail en utilisant son énergie
dégagée lors du changement d’état ou en le faisant passer à travers l’échangeur où il
effectuera un échange à contre courant avec le gaz de travail à refroidir (Fig. II.7).
17
Entrée gaz Entrée gaz

qk qk
P1 P1
R R
P2 P2
c c
2 1 8 2 1
qa1 q in E1 qa1 q in E1
1 1
8
2 2
G0 RP G0 RP
7 7
2 6 2 6
qa2 qa2
q in q in
2 E2 2 E2
3 3
VE VE
5 5
4 4
V V
B B
x x
Figure II.7 : Schémas des cycles à étranglement E.S à un seul étage de R.P
18
Figure II.8 : Diagramme du cycle à étranglement avec un étage de R.P
Evolutions thermodynamiques du cycle à E.S + R.P

𝟐𝟐′′ -𝟐𝟐′ : Refroidissement préalable isobare du gaz de travail dû à un échange avec le fluide
réfrigérant.
𝟐𝟐′ -3 : Refroidissement isobare.
5-1 : Réchauffement isobare d’une quantité (1-x) du gaz de travail.
𝐪𝐪𝐢𝐢𝐢𝐢𝐢𝐢 , 𝐪𝐪𝐢𝐢𝐢𝐢𝐢𝐢 : Pertes des échangeurs.
𝐪𝐪𝐚𝐚𝐚𝐚 , 𝐪𝐪𝐚𝐚𝐚𝐚 : Apport de chaleur du milieu extérieur.
∆𝐇𝐇𝐓𝐓𝟐𝟐′ : Coefficient d’étranglement isothermique.
19
La quantité de liquide produite sera :
∆𝐇𝐇𝐓𝐓𝟐𝟐′ − (𝐪𝐪𝐢𝐢𝐢𝐢𝐢𝐢 + 𝐪𝐪𝐚𝐚𝐚𝐚 )

𝐱𝐱 =
𝐇𝐇𝟔𝟔 − 𝐇𝐇𝐟𝐟 − 𝐪𝐪𝐢𝐢𝐢𝐢𝐢𝐢
Rendement du cycle E.S+R.P
• Pour la liquéfaction de l’air(Téb = −194,3°C), le rendement de ce cycle est égal à

15%.
• Pour la liquéfaction de l’hélium(Téb = −268,93°C), le rendement de ce cycle
fonctionnant avec deux étages de R.P est égal à 7,7%.
II.2.5. Le cycle combiné ou cycle de CLAUDE
Pour la liquéfaction des gaz permanents et le cas de l’hélium, la réalisation d’un cycle
combiné permettra d’obtenir des températures qui avoisinent le Zéro (0) absolu.
Le processus associé à ce cycle combine à la fois la détente isenthalpique (valve

d’étranglement) et la détente isentropique (détendeur adiabatique).
Le détendeur peut être branché de différentes manières :
• Au niveau initial de température en fin de compression (Fig. II.9).

• Au niveau intermédiaire de température durant l’écoulement « Aller » du gaz à travers
l’échangeur (Fig. II.10).
• Au niveau inférieur de température (Fig. II.11).
20
Entrée gaz
qk
P1
R
P2
c
qin1 E1
q a1
D
q in2 E2
q a2
VE
B
x
Figure II.9 : Schéma du processus combiné avec branchement du détendeur au niveau

initial de température.
21
Entrée gaz qk
P1
R
c P2
qin1 E1
q a1
D
qin2 E2
qa2
qin3
q a3 VE
V
B
Figure II.10: Schéma du processus combiné avec branchement du détendeur au niveau

intermédiaire de température.
22
Entrée gaz
qk
P1
R
c P2
2
1
qin1 E1
qa 1
3 9
qin2 E2
D
4r
qa 2 5
VE
8
7
Figure II.11: Schéma du processus combiné avec branchement du détendeur au niveau

inférieur de température.
23
Evolutions thermodynamiques du cycle combiné

2-3 : Refroidissement isobare.
3-𝟒𝟒𝐬𝐬 : Détente isentropique dans le détendeur (théorique).
3-𝟒𝟒𝐫𝐫 : Détente réel dans le détendeur (avec création d’entropie).
3-5 : Refroidissement isobare d’une quantité (1-x-G) du gaz de travail.
8-1 : Réchauffement isobare du gaz de travail.
𝐪𝐪𝐢𝐢𝐢𝐢𝐢𝐢 , 𝐪𝐪𝐢𝐢𝐢𝐢𝐢𝐢 : Imperfections des échangeurs.
𝐪𝐪𝐚𝐚𝐚𝐚 , 𝐪𝐪𝐚𝐚𝐚𝐚 : Apport de chaleur du milieu extérieur.
La quantité de liquide produite sera :
(𝟏𝟏 − 𝐆𝐆) (𝚫𝚫𝐇𝐇𝐓𝐓𝐓𝐓 − ∑𝐪𝐪𝟐𝟐 )

𝐱𝐱 =
(𝐡𝐡𝟗𝟗 − 𝐡𝐡𝐟𝐟 − 𝐪𝐪𝐢𝐢𝐢𝐢𝐢𝐢 )
T1 T2 2 1
Enthalpie
P2 P1
/
9
5
4
Tx f
4 r, 8
7 4s
Figure II.12: Schéma du cycle combiné associé au processus avec branchement du détendeur
au niveau inférieur de température.
24
Rendement du cycle combiné
• Pour la liquéfaction de l’air(Téb = −194,3°C), le rendement de ce cycle est égal à

25%.
• Pour la liquéfaction de l’hélium(Téb = −268,93°C), le rendement de ce cycle est égal
à 44%.
II.2.6. Cycle en cascade
Il est possible d'utiliser un processus de réfrigération en cascade, où les fluides subissent des
changements d’états (liquide - gaz) provoqués par des échanges thermiques dans des
échangeurs de chaleurs (Fig. II.13).
La particularité de ce processus est que deux (02) circuits sont couplés thermiquement par
leurs évapo-condenseur.
Condenseur qc
7 6
Cascade B
VE C
Evapo-Condenseur
8 5
3 2
VE C
Cascade A
4 1
Evaporateur
qx
Figure II.13: Schéma du procédé de liquéfaction en cascade
25
Evolutions thermodynamiques du cycle en cascade
Cascade A
1-2 : Compression isentropique dans un compresseur

2-3 : Désurchauffe isobare+condensation isotherme et isobare dans un évapo-condenseur
3-4 : Détente isenthalpique dans une V.E
4-1 : Evaporation isotherme et isobare dans un évaporateur avec une source chaude q x
entrante.
Cascade B
5-6 : Compression isentropique dans un compresseur
6-7 : Désurchauffe isobare+condensation isotherme et isobare dans un condenseur avec q c
sortant
7-8 : Détente isenthalpique dans une V.E
8-5 : Evaporation isotherme et isobare dans un évapo-condenseur
T
qc 6
7 2
B
3
8 5
A 1
4 qx
S
Figure II.14: Cycle associé au processus de liquéfaction en cascade
26
Chapitre III. Etude thermodynamique des installations et cycles cryogéniques de liquéfaction
de l’Hélium
Chapitre III. ETUDE THERMODYNAMIQUE DES INSTALLATIONS ET

CYCLES CRYOGENIQUES DE LIQUEFACTION DE L’HELIUM
III.1. Cycle à étranglement simple de LINDE ou cycle E.S
T[K]
H2 H1
q
k 300 2 1
∆T1
285
Entrée 1
Hélium
c
2 1
qin1 P2
P1
Enth
q a1 alpie
[kJ/kg ]
3
VE
5
4
4,2
f
(I1)
f
x S[kJ/kg K]
Figure III.1: Schéma du processus de liquéfaction à étranglement simple et cycle associé au

processus à E.S
L’équation du bilan du contour (I1) représenté dans le schéma (Fig. III.1) nous permet de
déterminer la quantité de gaz liquéfié (x) : 𝐒𝐒̇ = 𝐄𝐄̇
h 2 + q a1 = x hf + (1 − x)h1 ′
h1 ′ − h2 − q a1
x=
(h1 ′ − hf )
q in1 = h1 − h1 ′ h1 ′ = h1 − q in1
27
de l’Hélium
Notre équation du bilan s’écrit alors :
h1 − h2 − q in1 − q a1
x=
(h1 − hf − q in1 )
∆HT1 − q in1 − q a1
x=
(h1 − hf − q in1 )
∆𝐇𝐇𝐓𝐓𝐓𝐓 − ∑ 𝐪𝐪𝟏𝟏
𝐱𝐱 =
(𝐡𝐡𝟏𝟏 − 𝐡𝐡𝐟𝐟 − 𝐪𝐪𝐢𝐢𝐢𝐢𝐢𝐢 )
Calcul de l’effet isothermique d’étranglement ∆𝐇𝐇𝐓𝐓𝐓𝐓
Hypothèses :
On assume que les apports de chaleur au niveau de chaque étage sont de l’ordre de 3% de la
valeur de l’effet isothermique d’étranglement et que les imperfections des échangeurs de
chaleur sont de l’ordre de 5% de la température du bout chaud de l’échangeur de chaleur
considéré.
On opte au départ pour une pression optimale en fin de compression de 20 atm.
Pour des pressions P1 = 1atm et P2 = 20atm:
∆HT1 = h1 − h2 R
∆HT1 = 1574 − 1580 = −𝟔𝟔 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 < 0
Calcul du coefficient de liquéfaction
q in1 = Cp (T1 − T1 ′) = 5.93 (15) q a1 = −0,18 kJ/kg = 3%∆HT
q in1 = 77,895 kJ/kg ∆T1 = 15K = 5%T1
∆HT1 − ∑ q1 −6 − (−0.18 − 77,895)

x= =
(h1 − hf − q in1 ) (1574 − 9,901 − 77,895)
𝐱𝐱 = −𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎
Le coefficient de liquéfaction (x) obtenu est négatif alors qu’il doit être compris entre 0 et 1,
donc ce cycle à E.S ne pourra pas être réalisé pour liquéfier l’hélium.
28
de l’Hélium
III.2. Cycle à étranglement simple ES avec un étage de refroidissement

préalable (R.P) à l’aide de l’azote à 80K
Dans le but d’améliorer le cycle précédent, on utilisera un étage de R.P fonctionnant avec de
l’azote liquide à 80K (Fig. III.2).
T[K]
Entrée
Hélium H2 H1
qk
2 1
P1 300
∆T1
R 285
P2 1
c
8 2 1/
q a1 q in1 E1
2
2 6
2
80 ∆T2 6
Enthalpie
G0
RP P2 P1
[kJ/kg ]
N2 (80K) 7
2 6
qa2 qin2 E2
3
VE
4 5
V 4,2
f
f
S[kJ/kg K]
x
Figure III.2: Schéma du processus de liquéfaction et cycle à étranglement simple avec un

étage de R.P
29
de l’Hélium
Pour des pressions :P1 = 1atm et P2 = 20 atm :
∆HT2 = h6 − h2 ′
∆HT2 = 431 − 435,7 = − 𝟒𝟒, 𝟕𝟕 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 < 0 R
q in2 = Cp (T6 − T6 ′) = (5,195)(4) = 20,78 kJ/kg
q in2 = 20,78 kJ/kg
q a2 = − 0,141 kJ/kg
∆T2 = 4K
∆HT2 − ∑q 2 −4,7 – (−0,141 + 20,78)

x= =
h6 − hf − q in2 (431 − 9,901 − 20,78)
𝐱𝐱 = −𝟎𝟎. 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎
Le coefficient de liquéfaction (x) bien qu’exprimé au niveau de l’étage inférieur, il demeure

toujours négatif.
Donc le cycle à étranglement avec 1 étage de R.P fonctionnant avec du N 2 liquide à 80K ne
peut pas être réalisé pour liquéfier l’hélium.
30
de l’Hélium
III.3 Cycle à étranglement simple avec deux étages de refroidissement

préalable utilisant de l’azote à 80K et de l’azote sous vide à 65 K
On ajoute un deuxième étage de refroidissement préalable au processus précédent comportant
de l’azote sous vide à 65K (Fig. III.3).
Entrée
Hélium qk T[K]
P1
R H2 H1
P2 2 1
c 300
∆T1
8 2 1/ 285
1
13
q in1 E1
q a1
2
RP
G0
N2(80K) 7
2 6 2//
2/ 6
q in2 E2 80 ∆T2 6/
q a2
9 Enthalpie
65 ∆T3 10 11 [kJ/kg ]
9 11/
P2
G0 12
P1
N2 (65K)
Sous vide 10 11
q in3
E3
q a3
3
VE 5
4 4,2
f
V
B
f S[kJ/kg K]
de R.P : azote liquide saturé à 80K et azote sous vide à 65K
31
de l’Hélium
Pour des pressions P1 = 1atm et P2 = 20atm:
∆HT3 = h11 − h10
∆HT3 = 353 − 356,85 = − 𝟑𝟑, 𝟖𝟖𝟖𝟖 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 < 𝟎𝟎
q in3 = Cp (T11 − T11 ′) = 5,197(3,25) q a3 = −0,1155 kJ/kg
q in3 = 16,89 kJ/kg ∆T3 = 3,25 K
∆HT3 − ∑q 3 −3,85 – (−0.1155 + 16,89)

x= =
h11 − hf − q in3 (353 − 9,901 − 16,89)
𝐱𝐱 = −𝟎𝟎. 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎
Le cycle à étranglement avec 2 étages de R.P fonctionnant avec du N 2 liquide à la pression

atmosphérique et à 80 K et du N 2 sous vide à 65 K ne permet pas de liquéfier l’hélium.
L’hélium a une température de saturation très basse 4,2 K, on utilisera l’hydrogène liquide à
20 K comme gaz de refroidissement préalable en plus de l’azote liquide saturé à 80K à la
pression atmosphérique (Fig. III.4).
32
de l’Hélium
III.4. Cycle à étranglement simple avec deux étages de R.P fonctionnant

avec de l’azote à 80 K et de l’hydrogène à 20 K
T[K]
Entrée
Hélium qk H2 H1
P1 2 1
300
R
P2
∆T1
c 285
1
8 2 1/
13
q in1 E1
q a1
2 (I3)
RP
G0 2//
N2 (80K) 7 2/ 6
2 6 80 ∆T2 6/
q in2 E2 Enthalpie
q a2
[kJ/kg ]
(I4)
P2
9
P1
G0 12
H2 (20K)
(I2)
10 11 9
20 10 11
q in3 E3 ∆T3
q a3
11/
3
3
(I5)
VE 5
4,2
4 f 4 5
(I1)
B S[kJ/kg K]
f
x
de R.P : azote liquide saturé à 80K et hydrogène liquide à 20K
33
de l’Hélium
III.4.1. Calcul des caractéristiques thermodynamiques de l’installation

cryogénique avec P2=1MPa
∆HT3 = h11 − h10
∆HT3 = 118,6 − 114,2 = 𝟒𝟒, 𝟒𝟒 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
Bilan énergétique du contour (I1) :
q in3 = Cp ∆T3 = Cp (T11 − T11 ′) = 5,25 (1) q a3 = 0,132 kJ/kg
q in3 = 5,25 kJ/kg ∆T3 = 1 K
∆HT3 − ∑q 3 4,4 – (0.132 + 5,25)

x= =
h11 − hf − q in3 (118,6 − 9,901 − 5,25)
𝐱𝐱 = −𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 < 0
On obtient un coefficient de liquéfaction x<0

cryogénique avec P2=1.5MPa
∆HT3 = h11 − h10
∆HT3 = 118,6 − 112,2 = 𝟔𝟔, 𝟒𝟒 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
q in3 = Cp (T11 − T11 ′) = 5,25(1) q a3 = 0,192 kJ/kg
q in3 = 5,25 kJ/kg ∆T3 = 1K
∆HT3 − ∑q 3 6,4 – (0.192 + 5,25)

x= = = 0,009261
h11 − hf − q in3 (118,6 − 9,901 − 5,25)
𝐱𝐱 = 𝟎𝟎, 𝟗𝟗𝟗𝟗𝟗𝟗𝟗𝟗%
34
de l’Hélium
q a1 = − 0,12 kJ/kg q in1 = 77,895 kJ/kg
q a2 = − 0,105 kJ/kg q in2 = 20,780 kJ/kg
q a3 = 0,192 kJ/kg q in3 = 5,250 kJ/kg
Détermination de la quantité de chaleur totale 𝐪𝐪𝟎𝟎 𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭 extraite du gaz de travail par
R.P
Bilan énergétique du contour (I2): Ṡ = Ė
h2 + GoN2 h7 + q a1 + q a2 + GoH2 h12 + (1 − x)h11 ′ = GoN2 h8 + GoH2 h13 + h10 +

(1 − x)h1 ′
q 0 totale = q oN2 + q oH2 = (GoN2 h8 − GoN2 h7 ) + (GoH2 h13 − GoH2 h12 )
q 0 totale = (h2 − h10 )(1 − x)(h11 ′ − h1 ′) + q a1 + q a2
q 0 totale = (1578 − 112,2) + (1 − 0,009261)(113,35 − 1496,105) − 0,12 − 0,105
𝐪𝐪𝟎𝟎 𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭 = 𝟗𝟗𝟗𝟗, 𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
Calcul de 𝐡𝐡𝟐𝟐 ′′ à la sortie de l’échangeur E1
Bilan énergétique du contour (I3) : Ṡ = Ė
h2 + q a1 + (1 − x)h6 ′ = h1 ′′ + (1 − x)h1 ′
h2 ′′ = h2 + q a1 + (1 − x)(h6 ′ − h1 ′)
h2 ′′ = 1578 − 0,12 + (1 − 0,009261)(410,22 − 1496,105)
𝐡𝐡𝟐𝟐 ′′ = 𝟓𝟓𝟓𝟓𝟓𝟓, 𝟎𝟎𝟎𝟎 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐓𝐓𝟐𝟐 ′′ = 𝟗𝟗𝟗𝟗 𝐊𝐊
Calcul de 𝐡𝐡𝟗𝟗 à la sortie de l’échangeur E2
Bilan énergétique du contour (I4): Ṡ = Ė
h2 ′ + q a2 + (1 − x)h11 ′ = h9 + (1 − x)h6 ′
h9 = h2 ′ + q a2 + (1 − x)(h11 ′ − h6 ′)
h9 = 434,5 − 0,105 + (1 − 0,009261)(113,35 − 410,22)
𝐡𝐡𝟗𝟗 = 𝟏𝟏𝟏𝟏𝟏𝟏, 𝟐𝟐𝟐𝟐 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐓𝐓𝟗𝟗 = 𝟐𝟐𝟐𝟐𝟐𝟐
35
de l’Hélium
Calcul de 𝐡𝐡𝟑𝟑 (détente isenthalpique)
h3 = h4 = x hf + (1 − x)h5
h3 = h4 = 0,009261(9,901) + (1 − 0,009261)30,691
𝐡𝐡𝟑𝟑 = 𝐡𝐡𝟒𝟒 = 𝟑𝟑𝟑𝟑, 𝟓𝟓𝟓𝟓 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐓𝐓𝟑𝟑 = 𝟕𝟕, 𝟔𝟔 𝐊𝐊
azote liquide à 80 K et R.P avec hydrogène liquide à 20 K (P 2 =15atm).
Points P T H S
[atm] [K] [kJ/kg] [kJ/kg K]
1 1 300 1574,000 31,610
1/ 1 285 1496,105 31,340
2 15 300 1578,000 25,990
2// 15 93 502,066 19,891
2/ 15 80 434,50 19,110
6 1 80 431,00 24,750
6/ 1 76 410,22 24,470
9 15 24 140,274 12,645
10 15 20 112,20 11,570
11 1 20 118,60 17,520
11/ 1 19 113,35 17,245
3 15 7,6 30,50 5,152
4 1 4,2 30,50 8,213
5 1 4,2 30,69 8,510
f 1 4,2 9,901 3,551
36
de l’Hélium
10 11
∆T3
11/
f
4 5
Calcul des caractéristiques thermodynamiques de l’installation cryogénique
Travail du compresseur :
P 15
rT1 ln �P2 � 2,076 (300) ln( 1 )
1
Wk = =
ηk 0,65
𝐖𝐖𝐤𝐤 = 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐, 𝟕𝟕𝟕𝟕𝟕𝟕𝟕𝟕𝟕𝟕𝟕𝟕 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
37
de l’Hélium
Travail supplémentaire :
q 0 totale 95,62569
Wsupp = =
q x supp 1,05
𝐖𝐖𝐬𝐬𝐬𝐬𝐬𝐬𝐬𝐬 = 𝟗𝟗𝟗𝟗, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
Travail total fourni pour le fonctionnement de l’installation :
Wtotal = Wk + β Wsupp = 2594,728716 + 1,2(91,0720857)
𝐖𝐖𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭 = 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐, 𝟎𝟎𝟎𝟎𝟎𝟎 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
Travail spécifique :
Wtotal 2704,015219
Wo = =
x 0,009261
𝐖𝐖𝐨𝐨 = 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐, 𝟕𝟕𝟕𝟕 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤

qx
Coefficient frigorifique réel : εréel =
Wtotal
q x = x (h1 − hf ) = 0,009261(1574 − 9,901)
𝐪𝐪𝐱𝐱 = 𝟏𝟏𝟏𝟏, 𝟒𝟒𝟒𝟒𝟒𝟒𝟒𝟒 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
14,4851
ε𝐫𝐫é𝐞𝐞𝐞𝐞 =
2704,015219
𝛆𝛆𝐫𝐫é𝐞𝐞𝐞𝐞 = 𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎
Coefficient frigorifique théorique :
(h1 − hf )
εth =
T1 (s1 − sf ) − (h1 − hf )
1574 − 9,901
εth =
300(31,61 − 3,551) − (1574 − 9,901)
𝛆𝛆𝐭𝐭𝐭𝐭 = 𝟎𝟎, 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐
Le rendement thermodynamique :
ε𝐫𝐫é𝐞𝐞𝐞𝐞 0,005356885
η= = = 0,0235
εth 0,228215
𝛈𝛈 = 𝟐𝟐, 𝟑𝟑𝟑𝟑%
38
de l’Hélium
III.4.3. Calcul des paramètres de l’installation avec P2 =2MPa
Calcul de l’effet isothermique d’étranglement ∆𝐇𝐇𝐓𝐓𝐓𝐓 R
∆HT3 = h11 − h10
∆HT3 = 118,6 − 110,5 = 𝟖𝟖, 𝟏𝟏 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 > 0
q in3 = Cp (T11 − T11 ′) = 5,25(1) q a3 = 0,243 kJ/kg
q in3 = 5,25 kJ/kg ∆T3 = 1K
∆HT3 − ∑q 3 8,1 – (0.243 + 5,25)

x= = = 0,0252
h11 − hf − q in3 (118,6 − 9,901 − 5,25)
𝐱𝐱 = 𝟐𝟐, 𝟓𝟓𝟓𝟓%
q a1 = − 0,18 kJ/kg q in1 = 77,895 kJ/kg
q a2 = − 0,141 kJ/kg q in2 = 20,780 kJ/kg
q a3 = 0,243 kJ/kg q in3 = 5,250 kJ/kg
R.P
𝐪𝐪𝟎𝟎 𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭 = (h2 − h10 ) + (1 − x)(h11 − h1 ′) + q a1 + q a2
q0 totale = (1580 − 110,5) + (1 − 0,0252)(113,35 − 1496,105) − 0,18 − 0,141
𝐪𝐪𝟎𝟎 𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭 = 𝟏𝟏𝟏𝟏𝟏𝟏, 𝟐𝟐𝟐𝟐 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
Calcul de 𝐡𝐡′′𝟐𝟐 à la sortie de l’échangeur E1
h2 ′′ = h2 + q a1 + (1 − x)(h6 ′ − h1 ′)
h2 ′′ = 1580 − 0, 18 + (1 − 0, 0252) (410, 22 − 1496, 105)
𝐡𝐡𝟐𝟐 ′′ = 𝟓𝟓𝟓𝟓𝟓𝟓, 𝟑𝟑𝟑𝟑 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐓𝐓𝟐𝟐 ′′ = 𝟗𝟗𝟗𝟗, 𝟑𝟑𝟑𝟑𝟑𝟑 𝐊𝐊
h9 = h2 ′ + q a2 + (1 − x)(h11 ′ − h6 ′)
39
de l’Hélium
h9 = 435,7 − 0,141 + (1 − 0,0252)(113,35 − 410,22)
𝐡𝐡𝟗𝟗 = 𝟏𝟏𝟏𝟏𝟏𝟏, 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐓𝐓𝟗𝟗 = 𝟐𝟐𝟐𝟐 𝐊𝐊
h3 = h4 = x hf + (1 − x)h5 = 0, 0252 (9,901) + (1 − 0, 0252) 30,691
𝐡𝐡𝟑𝟑 = 𝐡𝐡𝟒𝟒 = 𝟑𝟑𝟑𝟑, 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐓𝐓𝟑𝟑 = 𝟕𝟕, 𝟑𝟑𝟑𝟑𝟑𝟑
points P T H S
1 1 300 1574,000 31,610
1′ 1 285 1496,105 31,340
2 20 300 1580,000 25,390
2′′ 20 96,37 521,299 19,480
2′ 20 80 435,700 18,510
6 1 80 431,000 24,750
6′ 1 76 410,220 24,470
9 20 26 146,170 12,430
10 20 20 110,500 10,230
11 1 20 118,600 17,520
11′ 1 19 113,350 17,245
3 20 7,34 30,167 4,600
4 1 4,2 30,167 8,138
5 1 4,2 30,691 8,510
f 1 4,2 9,901 3,551
40
de l’Hélium
10 11
∆T3
11/
f
4 5
à 80 K et hydrogène liquide à 20 K (P 2 =20atm).
41
de l’Hélium
P 20
rT1 ln(P2 ) 2,076 (300) ln( 1 )
1
wk = =
ηk 0,65
𝐰𝐰𝐤𝐤 = 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐, 𝟑𝟑𝟑𝟑 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤

q 0 totale 121,269426
Wsupp = =
q x supp 1,05
𝐖𝐖𝐬𝐬𝐬𝐬𝐬𝐬𝐬𝐬 = 𝟏𝟏𝟏𝟏𝟏𝟏, 𝟒𝟒𝟒𝟒 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
Wtotal = Wk + β Wsupp = 2870,3724 + 1,2(115,4947)
𝐖𝐖𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭 = 𝟑𝟑𝟑𝟑𝟑𝟑𝟑𝟑, 𝟗𝟗𝟗𝟗 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
Wtotal 3008,96604
Wo = =
x 0,0252
𝐖𝐖𝐨𝐨 = 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏, 𝟒𝟒𝟒𝟒 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤

qx
Wtotal
q x = x (h1 − hf ) = 0,0252(1574 − 9,901)
𝐪𝐪𝐱𝐱 = 𝟑𝟑𝟑𝟑, 𝟒𝟒𝟒𝟒 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
39,41529
ε𝐫𝐫é𝐞𝐞𝐞𝐞 =
3008,96604
ε𝐫𝐫é𝐞𝐞𝐞𝐞 = 𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎
(h1 − hf )
εth =
T1 (s1 − sf ) − (h1 − hf )
1574 − 9,901
εth =
300(31,61 − 3,551) − (1574 − 9,901)
𝛆𝛆𝐭𝐭𝐭𝐭 = 𝟎𝟎, 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐
42
de l’Hélium
εréel 0,0130992
η= = = 0,05739
εth 0,228215
𝛈𝛈 = 𝟓𝟓, 𝟕𝟕𝟕𝟕 %
III.4.4. Calcul des paramètres de l’installation avec P2=2.5MPa
Calcul de l’effet isothermique d’étranglement ∆𝐇𝐇𝐓𝐓𝐓𝐓 R
∆HT3 = h11 − h10
∆HT3 = 118,6 − 109,3 = 𝟗𝟗, 𝟑𝟑 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
q in3 = cp (T11 − T11 ′ ) = 5,25(1) q a3 = 0,279 kJ/kg
q in3 = 5,25 kJ/kg ∆T3 = 1K
∆HT3 − ∑q3 9,3 – (0.279 + 5,25)

x= = = 0,0364527
h11 − hF − qin3 (118,6 − 9,901 − 5,25)
𝐱𝐱 = 𝟑𝟑, 𝟔𝟔𝟔𝟔%
q a1 = − 0,21 kJ/kg q in1 = 77,895 kJ/kg
q a2 = − 0,177 kJ/kg q in2 = 20,780 kJ/kg
q a3 = 0,279 kJ/kg q in3 = 5,250 kJ/kg
R.P
𝐪𝐪𝟎𝟎 𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭 = (h2 − h10 ) + (1 − x)(h11 − h1 ′) + q a1 + q a2
q 0 totale = (1581 − 109,3) + (1 − 0,0364527)(113,35 − 1496,105) − 0,21 − 0,177
𝐪𝐪𝟎𝟎 𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭 = 𝟏𝟏𝟏𝟏𝟏𝟏, 𝟗𝟗𝟗𝟗 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
Calcul de 𝐡𝐡𝟐𝟐 ′′ à la sortie de l’échangeur E1
h2 ′′ = h2 + q a1 + (1 − x)(h6 ′ − h1 ′)
43
de l’Hélium
h2 ′′ = 1581 − 0,21 + (1 − 0,0364527)(410,22 − 1496,105)
𝐡𝐡𝟐𝟐 ′′ = 𝟓𝟓𝟓𝟓𝟓𝟓, 𝟒𝟒𝟒𝟒𝟒𝟒𝟒𝟒𝟒𝟒 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐓𝐓𝟐𝟐 ′′ = 𝟗𝟗𝟗𝟗, 𝟔𝟔𝟔𝟔𝟔𝟔 𝐊𝐊
h9 = h2 ′ + q a2 + (1 − x)(h11 ′ − h6 ′)
h9 = 436,9 − 0,177 + (1 − 0,0364527)(113,35 − 410,22)
𝐡𝐡𝟗𝟗 = 𝟏𝟏𝟏𝟏𝟏𝟏, 𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐓𝐓𝟗𝟗 = 𝟐𝟐𝟐𝟐, 𝟗𝟗𝟗𝟗𝟗𝟗𝟗𝟗 𝐊𝐊
h3 = h4 = x hf + (1 − x)h5
h3 = h4 = 0,0364527 (9,901) + (1 − 0,0364527)30,691
𝐡𝐡𝟑𝟑 = 𝐡𝐡𝟒𝟒 = 𝟐𝟐𝟐𝟐, 𝟗𝟗𝟗𝟗𝟗𝟗𝟗𝟗 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐓𝐓𝟑𝟑 = 𝟔𝟔, 𝟗𝟗 𝐊𝐊
points P T H S
1 1 300 1574,000 31,610
1/ 1 285 1496,105 31,340
2 25 300 1581,000 24,930
2// 25 98,63 534,488 19,143
2/ 25 80 436,900 18,050
6 1 80 431,000 24,750
6/ 1 76 410,220 24,470
9 25 26,933 150,674 12,120
10 25 20 109,300 10,350
11 1 20 118,600 17,520
11/ 1 19 113,350 17,245
3 25 6,9 29,933 4,1312
4 1 4,2 29,933 8,085
5 1 4,2 30,691 8,510
f 1 4,2 9,901 3,551
44
de l’Hélium
10 11
∆T3
11/
f 5
4
Figure III.7: Cycle de liquéfaction de l’hélium à 2 étages en coordonnées T-S : R.P avec
azote liquide saturé à 80K et hydrogène liquide à 20 K (P 2 =25atm)
P 25
rT1 ln(P2 ) 2,076 (300) ln( 1 )
1
Wk = =
ηk 0,65
𝐖𝐖𝐤𝐤 = 𝟑𝟑𝟑𝟑𝟑𝟑𝟑𝟑, 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
45
de l’Hélium
q 0 totale 138,963153
Wsupp = =
q x supp 1,05
𝐖𝐖𝐬𝐬𝐬𝐬𝐬𝐬𝐬𝐬 = 𝟏𝟏𝟏𝟏𝟏𝟏, 𝟑𝟑𝟑𝟑 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
Wtotal = Wk + β Wsupp = 3084,17825 + 1,2 (132,34586)
𝐖𝐖𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭 = 𝟑𝟑𝟑𝟑𝟑𝟑𝟑𝟑, 𝟗𝟗𝟗𝟗 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
Wtotal 3242,99328
W0 = =
x 0,0364527
𝐖𝐖𝐨𝐨 = 𝟖𝟖𝟖𝟖𝟖𝟖𝟖𝟖𝟖𝟖, 𝟒𝟒𝟒𝟒𝟒𝟒 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤

qx
Wtotal
q x = x (h1 − hf ) = 0,0364527(1574 − 9,901)
𝐪𝐪𝐱𝐱 = 𝟓𝟓𝟓𝟓, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
57,01563
εréel =
3242,99328
εréel = 𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎
(h1 − hf )
εth =
T1 (s1 − sf ) − (h1 − hf )
1574 − 9,901
εth =
300(31,61 − 3,551) − (1574 − 9,901)
𝛆𝛆𝐭𝐭𝐭𝐭 = 𝟎𝟎, 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐
ε𝐫𝐫é𝐞𝐞𝐞𝐞 0,0175811
η= =
εth 0,228215
𝛈𝛈 = 𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎 = 𝟕𝟕, 𝟕𝟕%
46
de l’Hélium
Tableau III.4 : Tableau récapitulatif des résultats des calculs des caractéristiques
thermodynamiques de l’installation cryogénique à 2 étages R.P avec azote liquide à 80 K et
R.P avec hydrogène liquide à 20 K.
P2 Wk X W0 q0 𝛆𝛆 réel
R 𝛈𝛈 th
R
[atm] [kJ/kg] [%] [kJ/kg] [kJ/kg] [%]

15 2594,73 0,93 291966,90 95,53 0,0053 2,35
20 2870,37 2,52 119403,41 121,27 0,0131 5,74
25 3084,18 3,65 88966,49 139,03 0,0176 7,70
III.4.5. Représentation graphique de l’évolution des caractéristiques

thermodynamiques de l’installation cryogénique à 2 étages
350
300 15; 291,9669042
250
W0 [MJ/kg]
200
150
20; 119,4034 W0 [MJ/kg]
100 25; 88,966488
50
0
0 5 10 15 20 25 30
P [atm]
Figure III.8 : Evolution de la consommation d’énergie spécifique en fonction de P 2
10
8 25; 7,7
6 20; 5,74
4 rendement [%]
25; 3,65
15; 2,35 20; 2,52 x [%]
2
15; 0,93
0
0 5 10 15 20 25 30
P [atm]
Figure III.9: Evolution du rendement thermodynamique et du coefficient de liquéfaction de

l’hélium en fonction de la pression P 2
On remarque une augmentation du rendement en abaissant la température et en augmentant la

pression de refoulement du compresseur.
47
de l’Hélium
III.5. Cycle combiné pour la liquéfaction de l’hélium avec deux étages de

refroidissement préalable R.P utilisant de l’azote à 80K et de l’hydrogène à
20 K avec branchement du détendeur à T=20K
Entrée T[K]
Hélium qk
P1 H2 H1
R 300 2 1
P2 ∆T1
c 1
15 2 1
17 q in1 E1
3
qa1 (I4)
3 80 ∆T2
4 13
RP
J0 13
N2 (80K) 14 P2
4 13
q in2 E2 5
qa2 P1
(I6) 6
20 ∆T3 12
5 12
Enthalpie
J0 RP (I5)
[kJ/kg ]
H2(20K) 16
6 12
G
1-G
q in3
1-x
E3
WD D
9 7 11
(I2) ∆T4
7 11 11
q a3
8 6s
1-G-x
1-G
q in4
E4
(I3)
8 10 4,2 f 10
q a4 9
VE
9
V
(I1) S[kJ/kg K]
f
x
Figure III.10 : Schéma du processus combiné et du cycle associé au processus de
liquéfaction de l’hélium avec un branchement du détendeur à T=20K et un taux de détente
Ʈ=25
48
de l’Hélium
azote liquide à Patm à 80 K et R.P avec hydrogène liquide à 20 K utilisant la combinaison de
la détente adiabatique et de l’étranglement isenthalpique (P 2 =25atm)
points P T H S
1 1 300 1574,000 31,610
1′ 1 285 1496,105 31,340
2 25 300 1581,000 24,930
3 25 138,2 740,800 20,880
4 25 80 436,900 18,050
13 1 80 431,000 24,750
13′ 1 76 410,220 24,470
5 25 36,96 207,079 13,864
6 25 20 109,300 10,350
12 1 20 118,600 17,520
12′ 1 19 113,350 17,254
6s 1 5,5 39,930 10,397
7 25 9 39,740 5,335
11 1 9 60,060 13,260
11′ 1 8 54,540 12,610
8 25 5,115 24,080 3,117
9 1 4,2 24,080 6,760
10 1 4,2 30,691 8,510
f 1 4,2 9,901 3,551

cryogénique avec P2=25 atm
Détermination de la quantité de gaz liquéfié

Bilan énergétique du Contour (I1) :
(1 − G)h7 + (1 − G)q a4 = x hf + (1 − G − x)h11 ′
h7 − Gh7 + q a4 − Gq a4 = xhf + h11 ′ − Gh11 ′ − xh11 ′
𝐆𝐆(𝐡𝐡𝟏𝟏𝟏𝟏 ′ − 𝐡𝐡𝟕𝟕 − 𝐪𝐪𝐚𝐚𝐚𝐚 ) + 𝐱𝐱(𝐡𝐡𝟏𝟏𝟏𝟏 ′ − 𝐡𝐡𝐟𝐟 ) = 𝐡𝐡𝟏𝟏𝟏𝟏 ′ − 𝐡𝐡𝟕𝟕 − 𝐪𝐪𝐚𝐚𝐚𝐚
Bilan énergétique du Contour (I2) :
h6 + (1 − G)q a3 + Gh11 ′ + (1 − G − x)h11 ′ = (1 − x)h12 ′ + Gh6 + (1 − G)h7
h6 + q a3 − Gq a3 + Gh11 ′ + h11 ′ − Gh11 ′ − xh11 ′ = h12 ′ − xh12 ′ + Gh6 + h7 − Gh7
𝐆𝐆(𝐡𝐡𝟕𝟕 − 𝐡𝐡𝟔𝟔 − 𝐪𝐪𝐚𝐚𝐚𝐚 ) + 𝐱𝐱(𝐡𝐡𝟏𝟏𝟏𝟏 ′ − 𝐡𝐡𝟏𝟏𝟏𝟏 ′) = 𝐡𝐡𝟏𝟏𝟏𝟏 ′ + 𝐡𝐡𝟕𝟕 − 𝐡𝐡𝟔𝟔 − 𝐡𝐡𝟏𝟏𝟏𝟏 ′ − 𝐪𝐪𝐚𝐚𝐚𝐚
∆HT3 = h12 − h6 = 118,6 − 109,3
49
de l’Hélium
∆𝐇𝐇𝐓𝐓𝐓𝐓 = 𝟗𝟗, 𝟑𝟑 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐪𝐪𝐚𝐚𝐚𝐚 = 𝟎𝟎, 𝟐𝟐𝟐𝟐𝟐𝟐 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
∆𝐇𝐇𝐓𝐓𝐓𝐓 = h11 − h7 = 60,06 − 39,74
∆𝐇𝐇𝐓𝐓𝐓𝐓 = 𝟐𝟐𝟐𝟐, 𝟑𝟑𝟑𝟑 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐪𝐪𝐚𝐚𝐚𝐚 = 𝟎𝟎, 𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
𝐪𝐪𝐢𝐢𝐢𝐢𝐢𝐢 = 𝟓𝟓, 𝟒𝟒𝟒𝟒𝟒𝟒 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
L’équation du contour (I1) s’écrit :
G(h11 ′ − h7 − q a4 ) + x(h11 ′ − hf ) = h11 ′ − h7 − q a4
G(54,54 − 39,74 − 0,6096) + x(54,54 − 9,901) = 54,54 − 39,74 − 0,6096
(14,1904)G + (44,639)x = 14,1904
𝐆𝐆 = 𝟏𝟏 − 𝟑𝟑, 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏 𝐱𝐱
L’équation du contour (I2) s’écrit :
G(h7 − h6 − q a3 ) + x(h12 ′ − h11 ′) = h12 ′ + h7 − h6 − h11 ′ − q a3
= (1 − 3,1457 x)(h7 − h6 − q a3 ) + x(h12 ′ − h11 ′) = h12 ′ + h7 − h6 − q a3 − h11 ′
(1 − 3,1457 x)[(39,74) − 109,3 − 0,279] + x(113,35 − 54,54)
= 113,35 + 39,74 − 109,3 − 0,279 − 54,54
(1 − 3,1457 x)(−69,84) + 58,81 x = −11,029
−69,84 + 219,695 x + 58,81 x = −11,029
−69,84 + 278,505 x = −11,029
(1 − G) (ΔHT4 − ∑q 4 )
x = 0,22645 =
(h11 − hf − q in4 )
𝐱𝐱 = 𝟎𝟎, 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 𝐱𝐱 = 𝟐𝟐𝟐𝟐, 𝟔𝟔𝟔𝟔%
G = 1 − 3,1457 x G = 1 − [3,1457 (0,22645)]
𝐆𝐆 = 𝟎𝟎, 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 kg
Détermination de la quantité de chaleur totale 𝐪𝐪𝟎𝟎 extraite du gaz de travail par R.P
h2 + q a1 + J0N2 h14 + q a2 + J0H2 h16 + (1 − x)h12 ′

= h6 + (1 − x)h1 ′ + J0N2 h15 + J0H2 h17
50
de l’Hélium
h2 + q a1 + q a2 + (1 − x)h12 ′ = h6 + (1 − x)h1 ′ + J0N2 (h15 − h14 ) + J0H2 (h17 − h16 )
q 0 = J0N2 (h15 − h14 ) + J0H2 (h17 − h16 )

= h2 + q a1 + q a2 + (1 − x)h12 ′ − (1 − x)h1 ′ − h6
q 0 = J0N2 (h15 − h14 ) + J0H2 (h17 − h16 ) = (h2 − h6 ) + (1 − x)(h12 ′ − h1 ′) + q a1 + q a2
q 0 = (1581 − 109,3) + (1 − 0,2265)(113,35 − 1496,105) − 0,21 − 0,177
𝐪𝐪𝟎𝟎 = 𝟒𝟒𝟒𝟒𝟒𝟒, 𝟔𝟔𝟔𝟔 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
Calcul de 𝐡𝐡𝟑𝟑 sortie de l’échangeur E1
Bilan énergétique du contour (I4):
h2 + q a1 + (1 − x)h13 ′ = (1 − x)h1 ′ + h3
h3 = h2 + (1 − x)(h13 ′ − h1 ′) + q a1
h3 = 1581 + (1 − 0,22645)(410,22 − 1496,105) − 0,21
𝐡𝐡𝟑𝟑 = 𝟕𝟕𝟕𝟕𝟕𝟕, 𝟖𝟖𝟖𝟖 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐓𝐓𝟑𝟑 = 𝟏𝟏𝟏𝟏𝟏𝟏, 𝟐𝟐 𝐊𝐊
Calcul de 𝐡𝐡𝟓𝟓 à la sortie de l’échangeur E2
h4 + q a2 + (1 − x)h12 ′ = (1 − x)h13 ′ + h5
h5 = h4 + (1 − x)(h12 ′ − h13 ′) + q a2
h5 = 436,9 + (1 − 0,22645)(113,35 − 410,22) − 0,177
𝐡𝐡𝟓𝟓 = 𝟐𝟐𝟐𝟐𝟐𝟐, 𝟎𝟎𝟎𝟎𝟎𝟎 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐓𝐓𝟓𝟓 = 𝟑𝟑𝟑𝟑, 𝟗𝟗𝟗𝟗 𝐊𝐊
Calcul de 𝐡𝐡𝟖𝟖 à la sortie de l’échangeur E4
(1 − G)h7 + (1 − G)q a4 + (1 − G − x)h10 = (1 − G)h8 + (1 − G − x)h11 ′
(1 − G)h8 = (1 − G)(h7 + q a4 + h10 − h11 ′) + x(h11 ′ − h10 )
x(h11 ′ − h10 )
h8 = (h7 + q a4 + h10 − h11 ′) + � �
(1 − G)
51
de l’Hélium
0,22645(54,54 − 30,691)
h8 = (39,74 + 0,6096 + 30,901 − 54,54) + � �
1 − 0,2876
𝐡𝐡𝟖𝟖 = 𝐡𝐡𝟗𝟗 = 𝟐𝟐𝟐𝟐, 𝟎𝟎𝟎𝟎 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐓𝐓𝟖𝟖 = 𝟓𝟓, 𝟏𝟏𝟏𝟏 𝐊𝐊
6 12
∆T3
12
7
∆T4 11
11
8 6s
f 10
9
Figure III.11: Représentation du cycle combiné en coordonnées T-S à 2 étages R.P avec
azote liquide saturé à 80 K et hydrogène liquide à 20 K utilisant la détente adiabatique et de
l’étranglement isenthalpique avec (P 2 = 25atm)
52
de l’Hélium
P 25
rT1 ln(P2 ) 2,076 (300) ln( 1 )
1
Wk = =
ηk 0,65
𝐖𝐖𝐤𝐤 = 𝟑𝟑𝟑𝟑𝟑𝟑𝟑𝟑, 𝟏𝟏𝟏𝟏 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
q 0 totale 401,6828
Wsupp = =
q x supp 1,05
𝐖𝐖𝐬𝐬𝐬𝐬𝐬𝐬𝐬𝐬 = 𝟑𝟑𝟑𝟑𝟑𝟑, 𝟓𝟓𝟓𝟓 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
Travail du détendeur :
WD = ηméc G(h6 − h11 ′) = 0,9 (0,2876) (109,3 − 54,54)
𝐖𝐖𝐃𝐃 = 𝟏𝟏𝟏𝟏, 𝟏𝟏𝟏𝟏 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
Travail total:
Wtotal = Wk + β Wsupp − WD = 3084,17825 + 1,2(382,555) − 14,174
𝐖𝐖𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭 = 𝟑𝟑𝟑𝟑𝟑𝟑𝟑𝟑, 𝟎𝟎𝟎𝟎 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
Wtotal 3529,07
W0 = =
x 0,22645
𝐖𝐖𝐨𝐨 = 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏, 𝟑𝟑𝟑𝟑 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤

qx
Coefficient frigorifique réel : εréel = W
total
q x = x (h1 − hf ) = 0,22645(1574 − 9,901)
𝐪𝐪𝐱𝐱 = 𝟑𝟑𝟑𝟑𝟑𝟑, 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏 𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤
354,1902
εréel =
3529,07
𝛆𝛆𝐫𝐫é𝐞𝐞𝐞𝐞 = 𝟎𝟎, 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏

(h1 − hf )
εth =
T1 (s1 − sf ) − (h1 − hf )
53
de l’Hélium
1574 − 9,901
ε𝐭𝐭𝐭𝐭 =
300(31,61 − 3,551) − (1574 − 9,901)
ε𝐭𝐭𝐭𝐭 = 𝟎𝟎, 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐
εréel 0,10036
η= =
εth 0,228215
𝛈𝛈 = 𝟎𝟎, 𝟒𝟒𝟒𝟒𝟒𝟒𝟒𝟒 = 𝟒𝟒𝟒𝟒, 𝟗𝟗𝟗𝟗%
III.6. Comparaison du cycle combiné avec le cycle à étranglement simple à

deux étages de R.P fonctionnant avec de l’azote liquide à 80 K et de
l’hydrogène liquide à 20 K à la pression P2=25atm
Tableau III.6 : Caractéristiques thermodynamiques du cycle combiné pour la liquéfaction de

l’hélium
Cycles Wk X W0 q0 𝛆𝛆 réel
R 𝛈𝛈 th
R
[kJ/kg] [%] [kJ/kg] [kJ/kg] [%]

ES+R.P 3084,18 3,65 88966,49 139,029 0,0176 7,70
Combiné 3084,18 22,65 15584,325 401,683 0,1004 43,97
Commentaire :
- Le coefficient de liquéfaction est nettement amélioré, il passe de 3.65 à 22,65 %

- La consommation d’énergie spécifique diminue énormément d’environ 89000 à 16000
kJ/kg
- Le rendement du cycle combiné est nettement supérieur à celui du cycle à
étranglement avec R.P (presque 6 fois plus grand).
54
Chapitre IV. Dimensionnement des principaux organes de l’installation cryogénique de
liquéfaction de l’Hélium
Chapitre IV. DIMENSIONNEMENT DES PRINCIPAUX ORGANES DE

L’INSTALLATION CRYOGENIQUE DE LIQUEFACTION DE L’HELIUM
IV.1 Introduction
Ce chapitre sera consacré à la détermination des principales caractéristiques techniques des
principaux équipements de l’installation cryogénique de liquéfaction de l’hélium ayant
enregistré les meilleures performances thermodynamiques à savoir celle du cycle combiné à
deux étages de refroidissement préalable. On se limitera au choix du compresseur ainsi qu’au
dimensionnement d’un échangeur de chaleur.
IV.2 Les compresseurs [30-32]

Les compresseurs sont des appareils indispensables dans une installation cryogénique, ils sont
utilisés afin d’augmenter la pression du gaz de travail, pour cela ils agissent sur le volume ou
la vitesse du gaz à comprimer.
IV.2.1 Les différents types de compresseurs

Les compresseurs sont choisis selon le type de gaz de travail (gaz pur, mélange gazeux…etc.)
ainsi que les pressions à atteindre.
Pour une application cryogénique on distingue deux (02) types de compresseurs souvent
utilisés :
- Les compresseurs volumétriques

La compression est réalisée par une action directe sur le volume d’une masse gazeuse, car la
réduction du volume est accompagnée par un accroissement de pression.
Les compresseurs fonctionnant avec ce principe sont les compresseurs à pistons et les
compresseurs à vis.
- Les turbo-compresseurs
Les turbo-compresseurs sont des machines tournantes, ils agissent sur la vitesse du gaz de
travail en transformant l’énergie cinétique en énergie de pression grâce à des roues à aubes
calées sur un arbre (doté d’une certaine vitesse de rotation).
Les compresseurs fonctionnant avec ce principe peuvent être de types centrifuges, axiaux ou
mixtes.
IV.2.2 Fonctionnement des compresseurs volumétriques à pistons

Les compresseurs à pistons sont très performants, leur conception dépend de la taille de
l’installation (petite, moyenne ou grande puissance) avec des possibilités de raccordement
hermétique, semi hermétique ou ouvert au moteur d’entrainement.
55
Les trois (03) étapes réalisées dans ce type de machines sont (l’admission, la compression et
le refoulement).
L’admission :
Le clapet d’admission s’ouvre. Le gaz de travail BP (Basse Pression) est aspiré dans le
cylindre jusqu'à ce que le piston descende au PMB (Point Mort Bas).
La compression :
Lorsque le piston commence à remonter vers le haut, le clapet d’aspiration se referme, dans ce
cas, la pression du gaz de travail augmente HP (Haute Pression).
Le refoulement :
A un certain niveau lors de la remontée du piston, le clapet de refoulement s’ouvre. Le piston

monte ensuite jusqu’au PMH (Point Mort Haut), le gaz de travail comprimé est quant à lui
évacué vers un collecteur.
IV.2.3 Fonctionnement des compresseurs à vis

Les compresseurs à vis sont à la fois, performants et très silencieux. Ils sont utilisés au niveau
des grandes installations avec des possibilités de raccordement semi hermétique ou ouverts
aux moteurs d’entrainement.
Avec ce type de compresseurs, la réduction du volume du gaz de travail se réalise suite à la

rotation de vis.
Il en existe deux types : les mono rotors-pignons (une seule vis) et les bi-rotors (deux vis).
• Les compresseurs à vis mono rotors-pignons :
Très peu utilisés industriellement, ces compresseurs sont constitués d’une vis centrale (arbre)
avec une forme spéciale afin d’être reliée à deux satellites (Fig. IV.1). Comme pour les
compresseurs à pistons, les 3 étapes réalisées dans ce type de machines sont : l’admission, la
compression et le refoulement.
56
Cannelure
Satellite
Rotor
Figure IV.1: Compresseur mono-vis
L’admission :
Lorsque la vis commence à tourner, le gaz de travail est véhiculé grâce aux cannelures de la
vis jusqu'à son arrivée dans le volume délimité par le carter.
La compression :
Ici vient le rôle des deux satellites, entrainés par le rotor, ils diminuent le volume du gaz de
travail qui les traverse. Il en découle l’augmentation de la pression.
Le refoulement :
En fin de compression, le gaz est refoulé par la rotation de la vis vers un collecteur de
refoulement.
• Les compresseurs à vis double rotors :
Ce compresseur est constitué de deux rotors (male et femelle) ayant des profils très spéciaux
formant des lobes (vis male) et des alvéoles (vis femelle).
Ces derniers tournent à de grandes vitesses et nécessitent une bonne lubrification (afin
d’assurer l’étanchéité et le refroidissement des pièces mécaniques).
La conception de plusieurs étages sera avantageuse pour des taux de compression encore plus
élevés.
57
HP
BP
Figure IV.2: Compresseur double-vis
L’admission :
Dès que les rotors se mettent en mouvement, le gaz de travail est transporté le long de la vis
depuis l’orifice d’aspiration jusqu'à remplissage des espaces inter lobaires.
La compression :
Le volume du gaz se réduit avec la rotation des vis, du fait que les lobes du rotor males
viennent s’aplatir sur les alvéoles des rotors femelles contenant le gaz de travail.
Le refoulement :
Le refoulement commence dés que le gaz de travail atteint l’orifice d’évacuation.
IV.2.4. Principaux organes des compresseurs à pistons
Vilebrequin :
Le vilebrequin est un dispositif mécanique qui permet la transformation du mouvement de

rotation en un mouvement linéaire rectiligne au piston par l’intermédiaire d’une bielle, et
inversement.
La bielle :
Elle assure la transmission du mouvement au piston grace aux articulations de ces extrêmités.
Le piston :
Le piston est la pièce coulissante dans un cylindre. Le mouvement de va et vient du piston

fait balayer le gaz de travail et permet donc de diminuer le volume tout en augmentant la
pression.
58
Les clapets :
Les clapets sont des pièces d’une très grande importance dans un compresseur à piston. Ces
pièces sont exposées au danger de rupture et à des perturbations à cause des différences de
pression et de température.
Quand le piston est au PMH, le clapet de refoulement est fermé (aucune différence de
pression de part et d’autre du clapet). Ici il doit y’avoir un espace nuisible au dessus du
cylindre contenant une petite quantité de gaz HP afin de maintenir le clapet d’aspiration fermé
jusqu'à un certain niveau (Fig. IV.3). En descendant vers PMB, le piston arrive à un certain
niveau ou la BP est supérieur à la pression qui règne dans le cylindre. Ici le clapet
d’admission s’ouvre automatiquement.
Gaz de travail avant Gaz de travail après

compression compression
Clapets
d’aspiration et
de refoulement
Piston
espace
nuisible
Cylindre
Bielle
Vilebrequin
Figure IV.3 : Liaison mécanique dans un compresseur à pistons
IV.3 Calcul thermodynamique des caractéristiques du compresseur à

pistons de l’installation de liquéfaction de l’hélium
Pour un rendement optimal de l’installation, il faudra dimensionner un compresseur à pistons

capable de réaliser un taux de compression optimal
IV.3.1 Choix du taux de compression optimal

P2 25
τ= =
P1 1
𝛕𝛕 = 𝟐𝟐𝟐𝟐
59
IV.3.2. Calcul du débit massique

En choisissant un compresseur à pistons d’une puissance :P = 1,5 MW et connaissant le
travail de compression: W25 atm = 3084178,25 J/kg), on détermine le débit massique du gaz
hélium par :
P = ṁ W
P 1,5 . 106
ṁ = =
W 3084178,25
𝐦𝐦̇ = 𝟎𝟎, 𝟒𝟒𝟒𝟒𝟒𝟒 𝐤𝐤𝐤𝐤/𝐬𝐬 = 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏, 𝟖𝟖𝟖𝟖𝟖𝟖 𝐤𝐤𝐤𝐤/𝐡𝐡
IV.3.3. Calcul du volume de gaz aspiré par le compresseur

ṁ 1750,872
Vaspiré = =
ρHe 0,16
𝐕𝐕𝐚𝐚𝐚𝐚𝐚𝐚𝐚𝐚𝐚𝐚é = 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏, 𝟗𝟗𝟗𝟗𝟗𝟗 𝐦𝐦𝟑𝟑 /𝐡𝐡
IV.3.4. Calcul du volume de gaz balayé par le compresseur

Pour un rendement volumétrique estimé à ηv = 0,81
Vaspiré
ηv = = 0,81
Vbalayé
Vaspiré 10942,947
Vbalayé = =
ηv 0,81
𝐕𝐕𝐛𝐛𝐛𝐛𝐛𝐛𝐛𝐛𝐛𝐛é = 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏, 𝟖𝟖𝟖𝟖𝟖𝟖 𝐦𝐦𝟑𝟑 /𝐡𝐡
IV.3.5. Calcul de la cylindrée

Vbalayé = Cy . N
Pour une vitesse de rotation choisie N= 370 tr/min = 22200 tr/h
Vbalayé 13509,811
Cy = =
N 22200
𝐂𝐂𝐲𝐲 = 𝟎𝟎, 𝟔𝟔𝟔𝟔𝟔𝟔 𝐦𝐦𝟑𝟑
IV.3.6. Choix des dimensions des pistons du compresseur

Pour une telle cylindrée on proposera :
Un alésage : 𝐃𝐃 = 𝟎𝟎, 𝟓𝟓𝟓𝟓𝟓𝟓 𝐦𝐦
Une course du piston : 𝐂𝐂 = 𝟎𝟎, 𝟕𝟕𝟕𝟕𝟕𝟕 𝐦𝐦
Un nombre de cylindres : 𝐧𝐧 = 𝟒𝟒
60
IV.3.6 Vérification des dimensions choisies

π D2 (3,14) (0,510)2
Cy = . Course. n = . 0,745 . 4
4 4
𝐂𝐂𝐲𝐲 = 𝟎𝟎, 𝟔𝟔𝟔𝟔𝟔𝟔 𝐦𝐦𝟑𝟑
Tableau IV.1 : Caractéristiques du compresseur à pistons dimensionné
𝐕𝐕𝐛𝐛𝐛𝐛𝐛𝐛𝐛𝐛𝐛𝐛é 𝐏𝐏𝐟𝐟𝐟𝐟𝐟𝐟𝐟𝐟𝐟𝐟𝐟𝐟 Vitesse Puissance Cylindrée

[ 𝐦𝐦𝟑𝟑 /𝐡𝐡] [𝐚𝐚𝐚𝐚𝐚𝐚] [𝐭𝐭𝐭𝐭/𝐦𝐦𝐦𝐦𝐦𝐦] [MW] [𝐦𝐦𝟑𝟑 ]
13509,811 25 370 1,5 0,609
IV.4. Dimensionnement de l’échangeur de chaleur [33-35]

IV.4.1. Définition
Les échangeurs de chaleur sont utilisés dans un grand nombre d'applications dans plusieurs
domaines (industrie, transport, quotidien…) pour des refroidissements ou des réchauffements
et pour des changements de phases.
Ce sont des dispositifs capables de récupérer de l’énergie thermique d’un fluide à un autre en
présence d’une différence de température.
Les échangeurs de chaleur sont surtout choisis en fonction des propriétés physiques des
fluides, de leur agressivité ainsi qu’en fonction de l’encombrement.
Les plus utilisés sont les échangeurs tubulaires, ils restent les plus adaptés vu les conditions de
pression et de température qu’ils peuvent supporter.
IV.4.2. Classification des échangeurs de chaleur

Les échangeurs de chaleur peuvent être classés en fonction
• De la conception de l’échangeur (tubulaires, à plaques…etc.)
• Du mode de transfert de chaleur (conduction, convection)
• Des critères de performance d’un échangeur (efficacité, coefficient d’échange global,
compacité)
• Des paramètres métallurgiques
• Du contact entre les fluides en circulation (direct ou à mélange)
• Du sens de l’écoulement des fluides.
On distingue :
• Les échangeurs à courants parallèles (co-courants) :
Les fluides chauds et froids circulent parallèlement dans le même sens.
61
T
Tce
Tcs
Tcs
Tfs
Tfe Tfs
Tfe
Tce x
Figure. IV.4 : Echangeur à co-courants
• Les échangeurs à contre-courants :

Les fluides chauds et froids circulent parallèlement dans des sens opposés.
T
Tce Tce
Tfs
Tfe Tfs Tcs
Tcs Tfe
x
Figure. IV.5 : Echangeur de chaleur à contre-courants
• Les échangeurs à courants croisés :

Les fluides chauds et froids circulent perpendiculairement.
62
Fluide
chaud
Fluide
froid
Figure. IV.6 : Echangeur de chaleur à courants croisés
• Les échangeurs avec plusieurs types d’écoulements :

Les fluides chauds et froids circulent dans tout les sens.
Tce
T
Tce
Tfe
Tf inter Tfs Tcs
Tfs Tf inter
Tfe
Tcs x
Figure. IV.7 : Echangeur de chaleur à tubes et calandre
IV.4.3. Les échangeurs à faisceaux tubulaires et calandre

Pour faire face aux problèmes d’encombrement avec de simples échangeurs et pour résister
aux applications de conditions extrêmes tout en aboutissant à d’importantes surfaces
d’échange, les chercheurs ont du développer d’autres géométries d’échange d’où l’apparition
des échangeurs à faisceaux tubulaires et calandre.
63
Figure IV.8: Vues diverses d’échangeurs de chaleur à faisceaux tubulaires et calandre
Ces échangeurs sont actuellement les plus répandus. Un premier fluide traverse le faisceau de
tubes tandis qu’un second circule dans la calandre (enveloppe). C’est ainsi que le transfert de
chaleur s’effectue sans que les fluides ne soient mélangés.
Ces échangeurs conçus par le standard américain TEMA existent en plusieurs modèles en
fonction du type de boites avant, de retour et des types d’enveloppe. Ce qui fera la différence
c’est le nombre de passages côté tubes et côté calandre (simple passe, 4-1, 4-2, avec
chicanes) permettant d’avoir plusieurs natures d’écoulement dans un seul échangeur (faisant
sa particularité).
Tcs Tfe Tce Tfe
Tcs Tfs
Tce Tfs
Echangeur 2-1 Echangeur 4-2
Figure. IV. 9 : Types d’échangeurs de chaleur à faisceaux tubulaires et calandre
64
IV.4.3.1. Principaux organes des échangeurs à tubes et calandre
Boite d’échangeur :
Ce sont des boites placés aux extrémités des tubes afin de distribuer le fluide d’un côté et de
le récupérer de l’autre. Leurs formes dépendent du nombre de passage du fluide côté tubes.
Calandre:
C’est une enveloppe qui englobe le faisceau tubulaire faisant circuler l’un des deux fluides en
un seul passage.
Ici des chicanes (plaques métalliques) peuvent êtres ajoutées pour que le fluide effectue
plusieurs passages côté calandre. Ces plaques peuvent également augmenter la résistance
mécanique de l’échangeur.
Plaques tubulaires :
Ce sont des plaques percées supportant les tubes à leurs extrémités par dudgeonnage ou par
soudure. Elles représentent les parties délicates des échangeurs à faisceau tubulaire à cause
des différences de pression et de température qui règnent entre le fluide dans la calandre et le
fluide dans les boites (tubes).
Le faisceau tubulaire:
Ce sont des tubes normalisés de même longueur avec des épaisseurs théoriquement
négligeables disposés (arrangés) en pas carré ou en pas triangulaire.
Pas
Pas carré
triangulaire
Figure. IV. 10 : Arrangement des tubes au niveau d’un échangeur de chaleur pour les
pas carré et triangulaire
65
IV.5 Calcul thermique de l’échangeur tubulaire E1 du cycle hélium à

étranglement simple avec R.P
T [K]
Δ T2 300
ΔTmax
285
98,63
Δ T1
76
Figure. IV. 11 : Variations de température du gaz de travail dans l’échangeur E1 du cycle à

E.S avec R.P
IV.5.1. Calcul du flux thermique

La puissance échangée est :
Ф = Фc = Фf
Ф = ṁ c Cpc ΔTc = ṁ f Cpf ΔTf
Ф = ṁ c Cpc (Tce − Tcs )
Ф = (0,486 )(5,194)(300 − 98.63)
Ф = 𝟓𝟓𝟓𝟓𝟓𝟓, 𝟕𝟕𝟕𝟕𝟕𝟕 𝐤𝐤𝐤𝐤
IV.5.2. Propriétés des fluides chaud et froid

Propriétés du fluide chaud :
Qv = 10942,947 m3 /h = 3,039 m3 /s
𝐦𝐦̇ 𝐜𝐜 = 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏 𝐤𝐤𝐤𝐤/𝐡𝐡 = 𝟎𝟎, 𝟒𝟒𝟒𝟒𝟒𝟒 𝐤𝐤𝐤𝐤/𝐬𝐬
Tableau IV.2 : Propriétés thermiques du fluide chaud
𝐓𝐓𝐦𝐦𝐦𝐦𝐦𝐦𝐦𝐦𝐦𝐦𝐦𝐦𝐦𝐦 𝐏𝐏𝐟𝐟𝐟𝐟𝐟𝐟𝐟𝐟𝐟𝐟𝐟𝐟 𝛒𝛒𝐜𝐜 𝐂𝐂𝐩𝐩𝐩𝐩 𝐂𝐂𝐯𝐯𝐯𝐯 𝛍𝛍𝐜𝐜 𝛌𝛌𝐜𝐜

[𝐊𝐊] [𝐚𝐚𝐚𝐚𝐚𝐚] [𝐤𝐤𝐤𝐤/𝐦𝐦𝟑𝟑 ] [𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐊𝐊] [𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐊𝐊] [𝐏𝐏𝐏𝐏. 𝐬𝐬] [𝐖𝐖/𝐦𝐦 𝐊𝐊]
199,315 25 5,911 5,194 3,124 15,32 10-6 0,1199
66
Cċ = ṁ c Cpc = (0,486)(5,194. 103 )
𝐂𝐂𝐜𝐜̇ = 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐, 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 𝐖𝐖/𝐊𝐊 = 𝟐𝟐, 𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓 𝐤𝐤𝐤𝐤/𝐊𝐊
Propriétés du fluide froid :
Ф = Фf = ṁ f Cpf ΔTf = ṁ f Cpf (Tfs − Tfe )
Ф 508,7125 . 103
ṁ f = =
Cpf ΔTf 5,193. 103 (285 − 76)
𝐦𝐦̇ 𝐟𝐟 = 𝟎𝟎, 𝟒𝟒𝟒𝟒𝟒𝟒 𝐤𝐤𝐤𝐤/𝐬𝐬
Tableau IV.3 : Propriétés thermiques du fluide froid
𝐓𝐓𝐦𝐦𝐦𝐦𝐦𝐦𝐦𝐦𝐦𝐦𝐦𝐦𝐦𝐦 𝐏𝐏𝐟𝐟𝐟𝐟𝐟𝐟𝐟𝐟𝐟𝐟𝐟𝐟 𝛒𝛒𝐟𝐟 𝐂𝐂𝐩𝐩𝐩𝐩 𝐂𝐂𝐯𝐯𝐯𝐯 𝛍𝛍𝐟𝐟 𝛌𝛌𝐟𝐟

[𝐊𝐊] [𝐚𝐚𝐚𝐚𝐚𝐚] [𝐤𝐤𝐤𝐤/𝐦𝐦𝟑𝟑 ] [𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐊𝐊] [𝐤𝐤𝐤𝐤/𝐤𝐤𝐤𝐤 𝐊𝐊] [𝐏𝐏𝐏𝐏. 𝐬𝐬] [𝐖𝐖/𝐦𝐦 𝐊𝐊]
180,5 1 0,26 5,193 3,116 14,5 10-6 0,1179
Ċf = ṁ f Cpf = (0,469)(5,193. 103 )
𝐂𝐂𝐟𝐟̇ = 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 𝐖𝐖/𝐊𝐊 = 𝟐𝟐, 𝟒𝟒𝟒𝟒𝟒𝟒𝟒𝟒𝟒𝟒𝟒𝟒𝟒𝟒 𝐤𝐤𝐤𝐤/𝐊𝐊
IV.5.3. Calcul du rapport des taux de capacités thermiques

Le rapport des débits de capacité thermique des deux fluides est :
Cmin Ċf
R= =
Cmax Cċ
Nous avons :
Фc = Фf Cċ ΔTc = Cḟ ΔTf
Ċ min Cḟ
R= =
Ċ max Cċ
ΔTc (300 − 98,63)

R= =
ΔTf (285 − 76)
𝐑𝐑 = 𝟎𝟎, 𝟗𝟗𝟗𝟗𝟗𝟗
IV.5.4. Calcul du flux thermique maximal

Фmax = Cmin ΔTmax = Ċf ΔTmax
Фmax = 2434,0311 (300 − 76)
Ф𝐦𝐦𝐦𝐦𝐦𝐦 = 𝟓𝟓𝟓𝟓𝟓𝟓, 𝟐𝟐𝟐𝟐 𝐤𝐤𝐤𝐤
67
IV.5.5. Calcul de l’efficacité de l’échangeur

L’efficacité de l’échangeur est :
Ф 508,7125
ε= =
Фmax 545,222966
𝛆𝛆 = 𝟎𝟎, 𝟗𝟗𝟗𝟗𝟗𝟗
IV.5.6. Calcul de la différence de température logarithmique moyenne 𝚫𝚫𝚫𝚫𝚫𝚫𝚫𝚫

ΔT2 − ΔT1
ΔTLM =
ΔT
ln(ΔT2 )
1
ΔT2 = 300 − 285 = 15 K
ΔT1 = 98,63 − 76 = 22,63 K
𝚫𝚫𝚫𝚫𝚫𝚫𝚫𝚫 = 𝟏𝟏𝟏𝟏, 𝟓𝟓𝟓𝟓 𝐊𝐊
IV.5.7. Correction de la différence de température logarithmique moyenne 𝚫𝚫𝚫𝚫𝚫𝚫𝚫𝚫

Pour un échangeur avec 1 passage côté calandre et 2 passages côté tubes :
𝐅𝐅𝐅𝐅𝐅𝐅𝐅𝐅𝐅𝐅𝐅𝐅𝐅𝐅 = 𝟎𝟎, 𝟗𝟗
ΔTLM𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐é = (ΔTLM)(Facteur de correction)
ΔTLMcorrigé = (18,554)(0,9)
𝚫𝚫𝚫𝚫𝚫𝚫𝚫𝚫𝐜𝐜𝐜𝐜𝐜𝐜𝐜𝐜𝐜𝐜𝐜𝐜é = 𝟏𝟏𝟏𝟏, 𝟔𝟔𝟔𝟔𝟔𝟔 𝐊𝐊
Dans ce cas de figure, il est nécessaire de dimensionner une série d’échangeurs à faisceaux
tubulaires (2-1) dans lesquels le fluide froid circule du côté des tubes disposé en pas carré
tandis que le fluide chaud circule du côté de la calandre.
d
Un Pas
(Pt)
Espace libre
(Pd)
Figure. IV. 12 : Arrangement des tubes au niveau d’un échangeur tubulaire (pas carré).
68
Tableau IV.4 : Normes des tubes en fonction du diamètre de la calandre
Diamètre de la calandre 𝐃𝐃𝐜𝐜 [m] 0,59

Nombre de tubes 𝐍𝐍𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭 42
Diamètre extérieur des tubes 𝐝𝐝𝐞𝐞𝐞𝐞𝐞𝐞 [m] 0,019
Diamètre intérieur des tubes 𝐝𝐝𝐢𝐢𝐢𝐢𝐢𝐢 [m] 0,0122
Pas des tubes 𝐏𝐏𝐭𝐭 [m] 0,0254
Epaisseur du tube e [m] 0,0034
Espace libre entre deux tubes 𝐏𝐏𝐝𝐝 [m] 0,0064
Espacement des chicanes 𝐋𝐋𝐛𝐛 [m] 0,18
Section d’un tube S [m2] 1,168.10-4
Section de 21 tubes 𝐒𝐒𝟐𝟐𝟐𝟐 𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭 [m2] 2,453. 10-3
IV.5.8. Calcul du coefficient d’échange interne (coté tubes)

La vitesse du fluide chaud est :
ṁ c = ρc . Vc . S21tubes
ṁ c 0,486
Vc = =
ρc . S21tubes (5,911)(2,453 . 10−3 )
𝐕𝐕𝐜𝐜 = 𝟑𝟑𝟑𝟑, 𝟓𝟓𝟓𝟓𝟓𝟓 𝐦𝐦/𝐬𝐬
Le nombre de Reynolds est alors égal à :
ρc . Vc . dint
Rec =
µc
(5,911)(33,518)(0,0122)
Rec =
15,32 . 10−6
𝐑𝐑𝐑𝐑𝐜𝐜 = 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏, 𝟕𝟕𝟕𝟕𝟕𝟕 > 10000 Le régime est turbulent.
Le nombre de Prandtl est :
µc . Cpc (15,32 . 10−6 )(5,194. 103 )

Prc = =
λc 0,1199
𝐏𝐏𝐏𝐏𝐜𝐜 = 𝟎𝟎, 𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔
Le nombre de Nusselt est calculé en utilisant la formule de (Dittus et Boelter) soit :
Nuc = 0,0241 Rec 0,8 Prc 0,4
Nuc = 0,0241 (157775,702)0,8 (0,66365)0,4
𝐍𝐍𝐍𝐍𝐜𝐜 = 𝟐𝟐𝟐𝟐𝟐𝟐, 𝟓𝟓𝟓𝟓𝟓𝟓
69
Le coefficient d’échange intérieur est :
Nuc . λc (294,595)(0,1199)
hc = =
dint 0,0122
𝐡𝐡𝐜𝐜 = 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐, 𝟐𝟐𝟐𝟐𝟐𝟐 𝐖𝐖/𝐦𝐦𝟐𝟐 𝐊𝐊
IV.5.9. Calcul du coefficient d’échange externe (coté calandre)

On utilise la méthode de Kern. La vitesse permettant de calculer le nombre de Reynolds est
déterminée à l’aide de la section disponible comprise entre deux chicanes et entre deux tubes
soit :
Pd
Sd = Dc . Lb � �
Pt
0,0064
Sd = (0,59)(0,18) � �
0,0254
𝐒𝐒𝐝𝐝 = 𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 𝐦𝐦𝟐𝟐
D’où la vitesse du fluide froid est :
ṁ f = ρf . Vf . Sd
ṁ f 0,469
Vf = =
ρf . Sd (0,26)(0,026759 )
𝐕𝐕𝐟𝐟 = 𝟔𝟔𝟔𝟔, 𝟒𝟒𝟒𝟒𝟒𝟒 𝐦𝐦/𝐬𝐬
Le nombre de Reynolds est par définition :
ρf . Vf . Déq
Ref =
µf
4 .S
D’où : Déq = pm
πd 2
4 �Pt 2 − � 4ext ��
Déq =
π . dext
3,14 . 0,0192
4 �0,02542 − � ��
4
Déq =
3,14 (0,019)
𝐃𝐃é𝐪𝐪 = 𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 𝐦𝐦
Le nombre de Reynolds est alors égal à :
70
(0,26)(67,411)(0,024255778)
Ref =
14,5 . 10−6
𝐑𝐑𝐑𝐑𝐟𝐟 = 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐, 𝟏𝟏𝟏𝟏𝟏𝟏 > 10000 Le régime est turbulent.
Le nombre de Prandtl est :
µf . Cpf (14,5. 10−6 )(5,193. 103 )

Prf = =
λf 0,11798
𝐏𝐏𝐏𝐏𝐟𝐟 = 𝟎𝟎, 𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔
Le nombre de Nusselt est donné par la corrélation de Kern:
Nuf = 0,36 Ref 0,55 Prf 1/3
Nuf = 0,36 (29319,147)0,55 (0,63823)1/3
𝐍𝐍𝐍𝐍𝐟𝐟 = 𝟖𝟖𝟖𝟖, 𝟕𝟕𝟕𝟕𝟕𝟕𝟕𝟕
Le coefficient d’échange extérieur est :
Nuf . λf (88,7622)(0,11798)
hf = =
De 0,024255778
𝐡𝐡𝐟𝐟 = 𝟒𝟒𝟒𝟒𝟒𝟒, 𝟕𝟕𝟕𝟕𝟕𝟕 𝐖𝐖/𝐦𝐦𝟐𝟐 𝐊𝐊
Tableau IV.5 : Tableau récapitulatif des caractéristiques de l’écoulement des fluides chaud et
froid
T T
entrée sortie 𝐦𝐦̇ V Re Pr Nu h
(K) (K) (kg/s) (m/s) (W/m2K)
Fluide 300 98,63 0,486 33,518 157775 0,66 294,595 2895,241
chaud
Fluide 76 285 0,469 67,411 29319,147 0,64 88,7622 431,739
froid
IV.5.10. Calcul du coefficient d’échange global

Le coefficient d’échange global est donné par la relation :
1 1 dint dext 1 dint dint

= + ln � �+ � � + R si + R se ( )
U hc 2. λp dint hf dext dext
Le coefficient d’échange global est calculé en supposant un échangeur neuf (sans

encrassement) et une conductivité thermique λp = 50 W/m. K :
71
1 1 0,0122 0,019 1 0,0122

= + ln � �+ � �
U 2895,241 2(50) 0,0122 431,739 0,019
𝐔𝐔 = 𝟓𝟓𝟓𝟓𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎 𝐖𝐖/𝐦𝐦𝟐𝟐 𝐊𝐊
IV.5.11. Calcul de la surface d’échange

On en déduit la surface d’échange :
Ф = U S ΔTLMcorrigé
Ф 508,713 . 103
S= =
U . ΔTLMcorrigé (385,835)(16,699)
𝐒𝐒 = 𝟓𝟓𝟓𝟓, 𝟒𝟒𝟒𝟒𝟒𝟒 𝐦𝐦𝟐𝟐
IV.5.12. Calcul de la longueur des tubes

La longueur des tubes dans la calandre est :
S = π dint Ltotale Ntubes
S
Ltotale =
π dint Ntubes
78,955
Ltotale =
(3,14)(0,0122)(21)
𝐋𝐋𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭𝐭 = 𝟕𝟕𝟕𝟕, 𝟒𝟒𝟒𝟒𝟒𝟒 𝐦𝐦
Ltotale 98,146
Lpassage = =
Npasses 16
𝐋𝐋𝐩𝐩𝐩𝐩𝐩𝐩𝐩𝐩𝐩𝐩𝐩𝐩𝐩𝐩 = 𝟒𝟒, 𝟒𝟒𝟒𝟒𝟒𝟒 𝐦𝐦
72
Figure. IV. 13 : Echangeur tubulaire E1 du cycle hélium à étranglement simple avec R.P
IV.6 Rapport entre l’efficacité de l’échangeur et le coefficient de

liquéfaction de l’installation cryogénique
IV.6.1 Cycle à étranglement simple avec R.P

h10 − h11 ′ + q a3
𝐱𝐱 =
hf − h11 ′
h11 ′ − h5
𝛆𝛆 =
h11 − h5
h11 ′ = ε(h11 − h5 ) + h5
73
On introduit h11 ′ dans l’expression du coefficient de liquéfaction (x) :
h10 + q a3 − ε(h11 − h5 ) − h5
x=
hf − ε(h11 − h5 ) − h5
Pour une efficacité parfaite de l’échangeur 𝛆𝛆 = 𝟏𝟏𝟏𝟏𝟏𝟏%
109,3 + 0,279 − 1(118,6 − 30,691) − 30,691

x=
9,901 − 1(118,6 − 30,691) − 30,691
𝐱𝐱 = 𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 = 𝟖𝟖, 𝟐𝟐𝟐𝟐 %
Pour une efficacité de l’échangeur 𝛆𝛆 = 𝟗𝟗𝟗𝟗 %
109,3 + 0,279 − 0,95(118,6 − 30,691) − 30,691

x=
9,901 − 0,95(118,6 − 30,691) − 30,691
𝐱𝐱 = 𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 = 𝟒𝟒, 𝟒𝟒𝟒𝟒 %
109,3 + 0,279 − 0,9(118,6 − 30,691) − 30,691

x=
9,901 − 0,9(118,6 − 30,691) − 30,691
𝐱𝐱 = 𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 = 𝟎𝟎, 𝟐𝟐𝟐𝟐 %
Résultat : la quantité de gaz liquéfiée est tributaire de l’efficacité de l’échangeur, plus cette
efficacité est élevée, plus le coefficient de liquéfaction est important.
IV.6.2 cycle combiné

(1 − G)(h6 − h10 ′ + q a4 )
x=
(hf − h10 ′)
h10 ′ − h9
𝛆𝛆 =
h10 − h9
h10 ′ = ε(h10 − h9 ) + h9
(1 − G)[h6 − h9 − 𝛆𝛆(h10 − h9 )] + q a4
x=
[hf − h9 − 𝛆𝛆(h10 − h9 )]
Pour une efficacité parfaite de l’échangeur 𝛆𝛆 = 𝟏𝟏𝟏𝟏𝟏𝟏%
(1 − 0,29)[45,2 − 30,69 − 1(65,52 − 30,69)] + 0,61

x=
[9,90 − 30,69 − 1(65,52 − 30,69)]
𝐱𝐱 = 𝟎𝟎, 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 = 𝟐𝟐𝟐𝟐, 𝟕𝟕 %
74
(1 − 0,29)[45,2 − 30,69 − 0,95(65,52 − 30,69)] + 0,61

x=
[9,90 − 30,69 − 0,95(65,52 − 30,69)]
𝐱𝐱 = 𝟎𝟎, 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 = 𝟐𝟐𝟐𝟐, 𝟐𝟐𝟐𝟐 %
(1 − 0,29)[45,2 − 30,69 − 0,9(65,52 − 30,69)] + 0,61

x=
[9,90 − 30,69 − 0,9(65,52 − 30,69)]
𝐱𝐱 = 𝟎𝟎, 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 = 𝟐𝟐𝟐𝟐, 𝟔𝟔𝟔𝟔 %
Résultat : La quantité de gaz liquéfiée est tributaire de l’efficacité de l’échangeur, plus cette
efficacité est élevée, plus le coefficient de liquéfaction est important.
NB : Le détendeur permet la compensation de l’inefficacité de l’échangeur.
75
Chapitre V. Simulation des cycles cryogéniques de la liquéfaction de l’Hélium en utilisant le
logiciel HYSYS v7.1
Chapitre V. SIMULATION DES CYCLES CRYOGENIQUES DE LA

LIQUEFACTION DE L’HELIUM EN UTILISANT LE LOGICIEL HYSYS V7.1
V.1. Etapes de simulation

Les étapes permettant d’obtenir une bonne simulation des cycles cryogéniques sont les
suivantes :
76
logiciel HYSYS v7.1
77
logiciel HYSYS v7.1
9
Double clic 10
Figure. V. 1 : Représentation des dix premières étapes à suivre pour la simulation des
cycles cryogéniques de la liquéfaction
78
logiciel HYSYS v7.1
Une fois ces dix étapes réalisées, il sera nécessaire d’entrer :
• Les pressions de l’ensemble des points figuratifs du cycle cryogénique de l’installation

(après avoir relier tous les composants de l’installation).
• Le débit et la composition chimique du gaz au point (1) entrée du cycle (la
composition = 1 pour un gaz pur à 100%).
• Les températures de l’ensemble des points figuratifs du cycle cryogénique de
l’installation (du point 1 jusqu’au dernier point de l’écoulement « aller»).
V.2. Liquéfaction de l’azote T=77K

Pour le refroidissement préalable de l’hélium, l’azote utilisé doit être à l’état liquide et sa
liquéfaction peut être effectuée soit par des « cycles à étranglement simple » (pouvant être
améliorés par des « cycles avec refroidissement préalable ») soit par des « cycles combinés »
afin d’obtenir un rendement optimal.
V.2.1 Cycle à étranglement simple pour la liquéfaction de l’azote
Figure. V. 2 : Installation du cycle à ES pour la liquéfaction de l’azote avec P 2 =200 bars
Le coefficient de liquéfaction x est :
ṁ f = (1 − x)ṁ c
ṁ c − ṁ f ṁ f 3353
x= =1− = 1−
ṁ c ṁ c 3600
𝐱𝐱 = 𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 = 𝟔𝟔, 𝟖𝟖𝟖𝟖%
79
logiciel HYSYS v7.1
Tableau V.1 : Points figuratifs du cycle à ES pour la liquéfaction de l’azote avec P 2 =200 bars
V.2.2. Cycle à étranglement simple avec refroidissement préalable utilisant un

fluide frigorigène (à T=230K) pour la liquéfaction de l’azote
Figure. V. 3 : Installation du cycle à ES+RP pour la liquéfaction de l’azote avec P 2 =200 bars
80
logiciel HYSYS v7.1
Le coefficient de liquéfaction est :
ṁ f = (1 − x)ṁ c
ṁ f 3102
x= 1− =1−
ṁ c 3600
𝐱𝐱 = 𝟎𝟎, 𝟏𝟏𝟏𝟏𝟏𝟏𝟏𝟏 = 𝟏𝟏𝟏𝟏, 𝟖𝟖𝟖𝟖
Tableau V.2 : Points figuratifs du cycle à ES+RP pour la liquéfaction de l’azote avec P 2 =200
bars
81
logiciel HYSYS v7.1
V.2.3. Cycle combiné pour la liquéfaction de l’azote avec branchement du

détendeur au niveau initial à T=300K
Figure. V. 4 : Installation du cycle combiné pour la liquéfaction de l’azote avec un

ṁ f = (1 − G − x)ṁ c
ṁ f 953
x = (1 − G) − � � = (1 − 0,52) − � �
ṁ c 3600
𝐱𝐱 = 𝟎𝟎, 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 = 𝟐𝟐𝟐𝟐, 𝟓𝟓𝟓𝟓%
82
logiciel HYSYS v7.1
83
logiciel HYSYS v7.1
V.2.4. Cycle combiné pour la liquéfaction de l’azote avec refroidissement

préalable pour un branchement du détendeur au niveau intermédiaire à T=230K

Ʈ=200
ṁ f = (1 − G − x)ṁ c
ṁ f 510,6
x = (1 − G) − � � = (1 − 0,596) − � �
ṁ c 3600
𝐱𝐱 = 𝟎𝟎, 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 = 𝟐𝟐𝟐𝟐, 𝟐𝟐𝟐𝟐%
84
logiciel HYSYS v7.1
Ʈ=200
85
logiciel HYSYS v7.1
V.2.5. Cycle combiné pour la liquéfaction de l’azote avec refroidissement

préalable à 230K pour un branchement du détendeur à un niveau intermédiaire à
T=160K

ṁ f = (1 − G − x)ṁ c
ṁ f 293
x = (1 − G) − � � = (1 − 0,654) − � �
ṁ c 3600
𝐱𝐱 = 𝟎𝟎, 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 = 𝟐𝟐𝟐𝟐, 𝟒𝟒𝟒𝟒%
86
logiciel HYSYS v7.1
V.3. Liquéfaction de l’hydrogène T=20K

Pour le refroidissement préalable de l’hélium, l’hydrogène utilisé doit être à l’état liquide et sa
liquéfaction peut être effectuée soit par des « cycles à étranglement simple avec
refroidissement préalable» soit par des « cycles combinés » afin d’obtenir un rendement
optimal.
87
logiciel HYSYS v7.1
V.3.1. Cycle à étranglement simple avec deux étages de refroidissement

préalable (utilisant de l’azote liquide saturé à la pression atmosphérique à 80K et
de l’azote sous vide 65K pour la liquéfaction de l’hydrogène
Figure. V. 7 : Installation du cycle à ES+RP pour la liquéfaction de l’hydrogène avec

P 2 =160atm
ṁ f = (1 − x)ṁ c
ṁ c − ṁ f ṁ f 3780
x= =1− = 1−
ṁ c ṁ c 4932
𝐱𝐱 = 𝟎𝟎, 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 = 𝟐𝟐𝟐𝟐, 𝟑𝟑𝟑𝟑%
88
logiciel HYSYS v7.1
Tableau V.6 : Points figuratifs du cycle à ES+RP pour la liquéfaction de l’hydrogène avec
P 2 =160atm
89
logiciel HYSYS v7.1
V.3.2. Cycle combiné pour la liquéfaction de l’hydrogène avec un étage de

refroidissement préalable utilisant de l’azote avec un branchement du détendeur
à T=80K
Figure. V. 8 : Installation du cycle combiné pour la liquéfaction de l’hydrogène avec un

Ʈ=120
ṁ f = (1 − G − x)ṁ c
ṁ f 845,2
x = (1 − G) − � � = (1 − 0,50) − � �
ṁ c 4932
𝐱𝐱 = 𝟎𝟎, 𝟑𝟑𝟑𝟑𝟑𝟑𝟑𝟑 = 𝟑𝟑𝟑𝟑, 𝟖𝟖𝟖𝟖%
90
logiciel HYSYS v7.1
Tableau V.7 : Points figuratifs du cycle combiné pour la liquéfaction de l’hydrogène avec un
Ʈ=120
91
logiciel HYSYS v7.1
V.4. Cycles cryogéniques pour la liquéfaction de l’hélium
V.4.1. Cycle à étranglement simple avec deux étages de refroidissement

préalable utilisant de l’azote liquide saturé à 80K et de l’hydrogène liquide à 20K
pour la liquéfaction de l’hélium
ṁ f = (1 − x)ṁ c
ṁ c − ṁ f ṁ f 1730,9
x= =1− = 1−
ṁ c ṁ c 1751
𝐱𝐱 = 𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 = 𝟏𝟏, 𝟏𝟏𝟏𝟏%
92
logiciel HYSYS v7.1
P 2 =15atm
93
logiciel HYSYS v7.1
Pour obtenir de meilleurs résultats, on augmente la pression à P 2 =20atm
ṁ f = (1 − x)ṁ c
ṁ c − ṁ f ṁ f 1699,63
x= =1− = 1−
ṁ c ṁ c 1751
𝐱𝐱 = 𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 = 𝟐𝟐, 𝟗𝟗𝟗𝟗%
94
logiciel HYSYS v7.1
P 2 =20atm
95
logiciel HYSYS v7.1
Pour ce cycle, l’installation de base est préservée et la pression augmente jusqu’à atteindre un
maximum de 25 atm.
ṁ f = (1 − x)ṁ c
ṁ c − ṁ f ṁ f 1694,09
x= =1− = 1−
ṁ c ṁ c 1751
𝐱𝐱 = 𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 = 𝟑𝟑, 𝟐𝟐𝟐𝟐%
96
logiciel HYSYS v7.1
P 2 =25atm
97
logiciel HYSYS v7.1
V.4.2. Cycle combiné pour la liquéfaction de l’hélium avec un étage de

refroidissement préalable utilisant de l’azote liquide pour un branchement du
détendeur à T=65K

Le coefficient de liquéfaction x est :
𝐱𝐱 = 𝟎𝟎%
98
logiciel HYSYS v7.1
V.4.3. Cycle combiné pour la liquéfaction de l’hélium avec deux étages de

refroidissement préalable utilisant de l’azote 80K et de l’azote sous vide 65K pour
un branchement du détendeur à T=45K

ṁ f = (1 − G − x)ṁ c
ṁ f 1381
x = (1 − G) − � � = (1 − 0,166) − � �
ṁ c 1751
𝐱𝐱 = 𝟎𝟎, 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 = 𝟒𝟒, 𝟓𝟓𝟓𝟓%
99
logiciel HYSYS v7.1
100
logiciel HYSYS v7.1
V.4.4. Cycle combiné pour la liquéfaction de l’hélium avec deux étages de

refroidissement préalable utilisant de l’azote à 80K et de l’hydrogène pour un
branchement du détendeur à T=20K

ṁ f = (1 − G − x)ṁ c
ṁ f 805,7
x = (1 − G) − � � = (1 − 0,2876) − � �
ṁ c 1751
𝐱𝐱 = 𝟎𝟎, 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 = 𝟐𝟐𝟐𝟐, 𝟐𝟐𝟐𝟐%
101
logiciel HYSYS v7.1
102
logiciel HYSYS v7.1
Figure. V. 15 : Allure des courbes de variation du flux de chaleur, de la différence de

température et du coefficient d’échange global en fonction de la température au niveau de
l’échangeur LNG-100 du cycle combiné avec branchement du détendeur à T=20 K pour la
liquéfaction de l’hélium.
103
logiciel HYSYS v7.1
V.5. Interprétation des résultats de la simulation des différents cycles

cryogéniques de liquéfaction des gaz utilisés : azote, hydrogène et hélium
Tableau V.13 : Liquéfaction de l’azote à P 2 =200 bars
Cycle cryogénique Coefficient

X (%)
Cycle ES 6.86
ES+RP 13.83
Combiné (branchement à 300 K) 21.53
Les résultats obtenus montrent que :
• L’azote peut être liquéfié avec toutes les installations cryogéniques étudiées (cycles à
ES, à ES+RP et les cycles combinés).
• Les cycles à étranglement simple peuvent être améliorés grâce à un refroidissement
préalable ainsi que des cycles combinés.
• Le détendeur peut être branché à n’importe quel niveau de température pour la
liquéfaction de l’azote (initial-intermédiaire-final).
Tableau V.14 : Liquéfaction de l’hydrogène à 160 et à 120 bars
Cycle cryogénique Coefficient

X (%)
Cycle ES 0
ES+RP (P 2 =160 bars) 23.36
Combiné (détendeur à T=80 K) (P 2 =120atm) 32.86
Les résultats obtenus ont montré que :
• L’hydrogène ne peut être liquéfié avec des cycles cryogéniques à étranglement simple.
Cela dit :
• Il a été démontré que celui-ci pouvait être liquéfié avec des cycles à étranglement
simple avec deux étages de refroidissement préalable utilisant de l’azote (à T=80K) et
de l’azote sous vide à (T=65K) ou alors avec des cycles combinés.
• Quant au détendeur, il ne peut être branché à n’importe quel niveau de température
pour la liquéfaction de l’hydrogène.
104
logiciel HYSYS v7.1
Tableau V.15 : Comparaison des résultats calculés et simulés des cycles à étranglement
simple avec deux (02) étages de refroidissement préalable utilisant de l’azote liquide saturé à
T=80K et l’hydrogène liquide à T=20K pour la liquéfaction de l’hélium
P 2 (atm) 15 20 25
Calculé 0.93 2.52 3.65

X(%) Simulé 1.14 2.93 3.25
∆𝑥𝑥 0,21 0,41 0,4
Calculée 93.53 121.27 139.0
𝑄𝑄0 (kJ/kg) Simulée 94.48 122.0 135.5
∆𝑄𝑄0 0,95 0,73 3,5
Tableau V.16 : Comparaison des résultats calculés et simulés des cycles combinés pour la
liquéfaction de l’hélium avec deux étages de refroidissement préalable utilisant de l’azote
liquide saturé à T=80K et l’hydrogène liquide pour un branchement du détendeur à T=20K
Paramètres Résultats des Résultats de la Ecarts

Calculs Simulation
X (%) 22,65 25,226 ∆𝑥𝑥 =2,576
q 0 (kJ/kg) 401,683 404,09236 ∆𝑄𝑄0 =2,409
Interprétation des résultats de la simulation :
• On voit que les résultats de la simulation sont presque similaires à ceux calculés :
Cette légère différence est due au fait que le logiciel ne dispose que de l’isotope
Hélium3 (ayant une température de liquéfaction égale à T=3,2K).
• L’hélium ne peut être liquéfié avec des cycles cryogéniques à étranglement simple, il
peut être liquéfié avec des cycles à étranglement simple à deux étages de
refroidissement préalable utilisant de l’azote liquide à T=80 K et l’hydrogène liquide à
T=20 K ou bien à l’aide des cycles combinés.
• Le détendeur quant à lui doit être branché à un niveau de température inférieur à 65 K
pour avoir un taux de détente Ʈ=25 permettant l’obtention de bons rendements.
105
Conclusion générale
CONCLUSION GENERALE
Ce projet de fin d’étude nous a permis d’étudier de près la liquéfaction des gaz ainsi que le
fonctionnement des différents cycles cryogéniques tout en se familiarisant avec les gaz les
plus connus et les plus utilisés dans le domaine de la cryogénie.
Au niveau de la centrale à oxygène, il a été possible d’observer le déroulement de la

liquéfaction par le procédé de la séparation des gaz de l’air (l’azote, l’oxygène et l’argon) où
nous avons pu voir de près le fonctionnement de ce cycle industriel (Processus détaillés
précédemment dans le chapitre I).
Le but majeur de cette étude fut le calcul thermodynamique et la simulation des cycles
cryogéniques de la liquéfaction de l’hélium à 4.2 K.
Il a été démontré tout au long de ce travail que :
- Le refroidissement préalable utilisant l’azote liquide et l’hydrogène liquide est

indispensable pour la liquéfaction de l’hélium car atteindre une très basse température 4.2 K
nécessite l’intervention d’agents cryogéniques externes fonctionnant suivant des cycles
extérieurs opérant eux-mêmes des températures très basses tels que l’azote liquide à 80 K et
l’hydrogène liquide à 20 K.
- La pression finale de l’hélium dans le cycle à étranglement simple avec refroidissement

préalable utilisant de l’azote liquide saturé à 80 K et de l’hydrogène liquide à 20 K, donnant
de bons résultats est P 2 =25atm avec un rendement thermodynamique de l’installation
cryogénique de η = 7,7% contre à un rendement de η = 5,74% et η = 2,35% pour des
pressions de 20atm et 15atm respectivement.
- Le coefficient de liquéfaction s’améliore considérablement lorsqu’on passe du cycle ES+RP

au cycle combiné, il passe de 3.65% à 22.65%.
- La quantité de chaleur extraite par RP dans le cas du cycle ES+RP est de consommation
d’énergie 139 kJ/kg, elle est plus importante pour le cycle combiné et vaut 401.68 kJ/kg.
- La consommation d’énergie spécifique diminue pour l’ensemble des cycles examinés

lorsque le taux de compression augmente. Elle est de 89 MJ/kg pour le cycle ES+RP à
25atm et vaut 15.6 MJ/kg pour le cycle combiné. Ce qui constitue un gain énergétique et
économique très important
- Le cycle approprié pour la liquéfaction de l’hélium est le cycle combiné avec un

branchement du détendeur à uns niveau de température strictement inférieur à 65K,
température de référence pour le refroidissement préalable par l’azote sous vide
atmosphérique, car d’une part :
- Le cycle combiné avec branchement du détendeur à T=45 K (température d’inversion

supérieur de l’hélium), utilisant un refroidissement préalable avec de l’azote liquide
106
Conclusion générale
avait donné un η = 4,53% (dans ce cas, le détendeur a remplacé le refroidissement

préalable fonctionnant à l’hydrogène liquide).
D’autre-part :
- Le cycle combiné avec branchement du détendeur à T=20 K (utilisant un

refroidissement préalable avec de l’azote et de l’hydrogène liquide) avait quant à lui
donné un très bon rendement thermodynamique η = 43,98% (x6 plus grand).
- Les principaux équipements : compresseur et échangeur de chaleur ont été dimensionnés

pour assurer les performances désirées de l’installation cryogénique.
- Au cours de cette étude, la simulation des cycles cryogéniques par le logiciel HYSYS V.7
fut un excellent moyen permettant de confirmer la fiabilité et la précision dans la
détermination de l’ensemble des caractéristiques thermodynamiques calculées pour les
diverses installations cryogéniques examinées.
107
Références
REFERENCES
[1] PETER WILLIAM ATKINS Elements of Physical Chemistry, Second Edition, Oxford
University Press 1996.
[2] http://lepuitdesciences-com.over-blog.com/2015/04/les-etats-de-la-matiere.html
[3] WILLIAM RAMSAY, Les gaz de l'atmosphère: histoire de leur découverte
[4] https://www.lindegas.fr/fr/products_and_supply/gases_atmospheric/index.html
[5] HUBERT REEVES et JANINE THIBAULT « L’air au quotidien » approche théorique et

expérimentale, Association Graine de chimiste P58
[6] https://www.futura-sciences.com/sciences/definitions/chimie-air-4452/
[7] https://fr.wikipedia.org/wiki/Liqu%C3%A9faction
[8] Techniques de l'ingénieur, opération unitaire traitement des gaz P18
[9] Manuel opératoire de « AIR LIQUIDE », la séparation des gaz de l’air, centre de
perfectionnement TRAINING CENTER
[10] LENNTECH, chemical properties of nitrogen-Healt effects of nitrogen-Envirenmental

effects of nitrogen, Site Web: https://www.lenntech.com/periodic/elements/n.htm
[11] LENNTECH, chemical properties of argon-Health effects of argon-Environmental

effects of argon, Site Web: https://www.lenntech.com/periodic/elements/ar.htm
[12] LENNTECH, chemical properties of oxygen-Health effects of oxygen-Environmental

effects of oxygen, Site Web: https://www.lenntech.com/periodic/elements/o.htm
[13] LENNTECH, chemical properties of hydrogen-Health effects of hydrogen-

Environmental effects of hydrogen, Sites Web:
https://www.lenntech.com/periodic/elements/h.htm https://fr.wikipedia.org/wiki/Hydrogène
[14] A. MIRAOUI et B. BLUNIER, 20 questions sur la pile à combustible- l’hydrogène

vecteur énergétique de demain, © Edition technip, Paris, 2009
[15] Techniques de l'ingénieur, Cryogénie, mise en œuvre, Propriétés physiques usuelles des
fluides cryogéniques B 2 381
[16] LENNTECH, chemical properties of helium-Health effects of helium-Environmental

effects of helium, Site Web: https://www.lenntech.com/periodic/elements/he.htm
[17] B. LECRERCQ et B. DIU, La physique mot à mot Hélium
108
Références
[18] L’Hélium et ses caractéristiques, la découverte de l’Hélium, Site Web : http://lhelium.e-

monsite.com/pages/page/sommaire.html
[19] S.W.VAN SCIVER, Helium Cryogenics, International Cryogenics Monograph series,©

Springer Science Business Media,2012
[20] Techniques de l’ingénieur Froid cryogénie et applications industrielles BE9817 P44 -P46
[21] https://fr.wikipedia.org/wiki/Air_liquide
[22] https://fr.wikipedia.org/wiki/Linde_AG
[23] https://www.info-chimie.fr/air-products-et-sonatrach-voient-encore-plus-grand-en-
algerie,95219
[24] http://www.aps.dz/economie/79930-gaz-industriels-sonatrach-
[25] http://freveille.free.fr/Eau_liquide_vapeur.html
[26] La thermodynamique appliquée aux systèmes énergétiques, Site Web:

https://direns.mines-paristech.fr/Sites/Thopt/fr/co/cycles-combines.html
[27] Techniques de l’ingénieur, Froid cryogénie et applications industrielles « Liquéfaction de

l’hélium et réfrigération à l’hélium de moyennes et fortes puissances » P35
[28] G. VENKATARATHNAM, Cryogenic Mixed Refrigerant Processes, The International

Cryogenics General Editors 2008. 249 Pages
[29] MILIND. D. ATREY, Manuel Cryogenic Engineering Department of Mechanical

Engineering
[30] Manuel d'instruction pour compresseur d'oxygène à piston sec type 4D375-3A
Installation: Société Nationale de Sidérurgie d'El-Hadjar, Algérie
[31] http://www.ef4.be/fr/pompes-a-chaleur/composants/le-compresseur.html
[32] M. JULIEN, Hydraulique industrielle, Théorie et exercices © 2017 les Editions JFD
P.374
[33] https://energie.wallonie.be/fr/differents-types-d-echangeurs-de
chaleurs.html?IDC=8049&IDD=97759
[34] R. LELEU, Transferts de chaleur, Paris, Éditions techniques de l'ingénieur, 1992.
[35] Techniques de l’ingénieur, Echangeurs de chaleur-Dimensionnement thermique BE9517
V1, 2014
Consultation des mémoires :
109
Références
Etude et optimisation de la maintenance de la COX du complexe Sider EL-Hadjar, TAHAR

Ahlam, Université Badji Mokhtar Annaba, 2017.
Calcul thermique d’une installation cryogénique de liquéfaction de l’hydrogène,CHAIB

RASSOU Sabrina, Université Badji Mokhtar Annaba, 2018.
Etude thermodynamique et simulation du cycle de liquéfaction d’azote, AMIRAT Mohamed

Cherif, Université Badji Mokhtar Annaba, 2018.
110
Annexes
ANNEXES
111
Annexes
Annexe I : Représentation des cycles de liquéfaction de l’Azote à E.S et à (E.S+R.P)
2 2 11
Diagramme T-S 1’ 1’
Azote [65K-300K]
Cycle à E.S
Cycle à E.S avec R.P
3
4 8
8’
3
5
ff
57
6 4
112
Annexes
Annexe II : Représentation du cycle combiné de liquéfaction de l’Azote avec

branchement du détendeur au niveau initial de température
2,3' 1
∆T1
1’
Diagramme T-S
Azote [65K-300K]
Cycle combiné avec
détendeur à T=300K
4
∆T2
8
8’
5
f
6 7
3's
113
Annexes
Annexe III: Représentation du cycle combiné de liquéfaction de l’Hydrogène avec

Diagramme T-S
Hydrogène
[0K-150K]
Cycle combiné avec
détendeur à T=80K
4 10
ΔT2
10’
ΔT3 5 9
6 9’
f
9
7 4S
114
Annexes
Annexe IV : Représentation du cycle combiné de liquéfaction de l’Hélium avec

2 1
1’
Diagramme T-S
Hélium
[15K-300K]
Cycle combiné avec
détendeur à T=45K
4 15
15’
5
6 14
14’
7 13
13’
8 12
12’
115
Annexes
Diagramme T-S
Hélium
[3K-25K]
Cycle combiné avec
détendeur à T=45K
8
12
12’
7s
f 11
10
116
Annexes
Annexe V : Représentation des cycles de liquéfaction de l’Hélium les plus rentables

(Cycle combiné avec branchement du détendeur à T=20K et le cycle à ES+RP)
2 2 1 1
1’ 1’
Diagramme T-S
Hélium
[15K-300K]
Cycle combiné avec
détendeur à T=20K
2’’
4 2’ 6 13
6’ 13’
6 12
12’
117
Annexes
Diagramme T-S
Hélium
[3K-25K]
Cycle combiné avec
détendeur à T=20K
6 10 11 12
11’ 12’
7
11
11’
3
8 6S
ff 5 10
9 4
118

Laouar Imed

Transféré par

Informations du documentcliquez pour développer les informations du document

Informations du documentcliquez pour développer les informations du document

Droits d'auteur :

Formats disponibles

Laouar Imed

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Laouar Imed

Transféré par

Droits d'auteur :

Formats disponibles

‫وزارة اﻟﺘﻌﻠﯿــﻢ اﻟﻌﺎﻟـﻲ واﻟﺒﺤـﺚ اﻟﻌﻠﻤـﻲ‬

BADJI MOKHTAR- ANNABA UNIVERSITY ‫ﻋﻨﺎﺑ���ﺔ‬- ‫ﺟﺎﻣﻌـــــ���ـﺔ ﺑﺎﺟــــ���ـﻲ ﻣﺨﺘــ���ـﺎر‬

FACULTÉ DES SCIENCES DE L’INGÉNIORAT

Au nom d’Allah le Tout Miséricordieux le Très Miséricordieux

Je dédie ce mémoire à mes très

Figure I.1 : Diagramme Pression-Température des phases de l’hélium

LISTE DES TABLEAUX

IV.5.7. Correction de la différence de température logarithmique moyenne 𝚫𝚫𝚫𝚫𝚫𝚫𝚫𝚫 ............... 68 17T

IV.5.8. Calcul du coefficient d’échange interne (coté tubes)...................................................... 69

Ce mémoire a pour objectif le calcul thermodynamique des cycles cryogéniques permettant de

- Chapitre I : Généralités sur les gaz industriels

Chapitre I. GENERALITES SUR LES GAZ INDUSTRIELS

I.2. Les gaz industriels

I.3. Définitions et propriétés de quelques gaz permanents

Tableau I.1 : Principaux constituants de l'air sec

Tableau I.2 : Température de saturation des principaux composants de l'air

Tableau I.3 : Propriétés physiques de l’air sec

Masse molaire (g/mol) 28,9653

Chaleur latente de vaporisation (air) : L v = 204.1 kJ/kg

Conductivité thermique (air) : 0.0234 W/m. K

Etape 1 : Préparation de l’air

• Aspiration et filtration de l’air ;

Etape 2 : refroidissement cryogénique et séparation des gaz de l’air

I.3.2. L’Azote [10]

Tableau I.4 : Propriétés physiques de l’azote

I.3.3. L’Argon [11]

Tableau I.5 : Propriétés physiques de l’argon

Rayon atomique 71pm

I.3.4. L’Oxygène [12]

Tableau I.6 : Propriétés physiques de l’oxygène

I.3.5. L’Hydrogène [13-15]

Tableau I.7 : Propriétés physiques de l’hydrogène

Caractéristique de sa phase gazeuse :

Masse volumique : 0,085 kg/m3

Conductivité thermique : 0.16835 W/m. K

Caractéristique de sa phase liquide :

Masse volumique : 70.97 kg/m3

Chaleur latente de vaporisation : 454,3 kJ/kg

Caractéristique de sa phase solide :

Chaleur latente de fusion: 58,158 kJ/kg

I.3.6. L’Hélium [16-17]

Tableau I.8 : Propriétés physiques de l’hélium

I.4. Localisation de l’hélium [16-20]

La Tanzanie pourrait devenir un nouvel acteur du marché de l'hélium, après la récente

I.4.2. Utilisation et particularité de l’hélium

Les deux isotopes de l’hélium les plus connus sont :

• En médecine pour scanner les corps humains.

Tableau I.9 : Propriété des isotopes stables de l’hélium

Figure I.1 : Diagramme Pression-Température des phases de l’hélium

• Le liquide cesse de bouillir et commence à présenter des propriétés complètement

I.4.3. Particularité de l’hélium

I.5. Les principaux producteurs de gaz industriels [21-24]

I.5.1. Air Products

I.5.3. Air Liquide

I.5.5. L’entreprise nationale Sonatrach

I.5.6. La centrale à oxygène de l’entreprise nationale SIDER EL HADJAR-Annaba

La centrale à oxygène est considérée comme étant le poumon du complexe sidérurgique

BADJI MOKHTAR- ANNABA UNIVERSITY ‫ﻋﻨﺎﺑ��ﺔ‬- ‫ﺟﺎﻣﻌـــــ��ـﺔ ﺑﺎﺟــــ��ـﻲ ﻣﺨﺘــ��ـﺎر‬