Cycle Frigorifique à Compression de Vapeur — Simulateur & Guide Technique

La suite technique interactive ultime pour analyser le Cycle Frigorifique à Compression de Vapeur. Calculez les points d'état thermodynamiques, résolvez les diagrammes de phase $P-h$ et $T-s$, évaluez le travail du compresseur et comparez les performances du R134a, R22, R410A, R717 (Ammoniac), R744 (CO₂) et R1234yf.

❄️ Simulateur VCRC — Cycle Frigorifique à Compression de Vapeur

Physics-based single-stage cycle solver · SI units internally · Educational-grade approximations clearly labeled

⚙️ Paramètres du Cycle

📊 Tableau des Points d'État

ÉtatP (kPa)T (°C)h (kJ/kg)s (kJ/kg·K)xv (m³/kg)

📈 Diagramme P–h log P

📘 Panneau de Formules Docs

Q_e = ṁ(h₁ − h₄)
Refrigeration effect: heat absorbed in evaporator per unit time.
W = ṁ(h₂ − h₁)
Actual compressor work input.
Q_c = ṁ(h₂ − h₃)
Condenser heat rejection. Conservation: Q_c = Q_e + W.
COP_R = Q_e / W , COP_HP = Q_c / W
Refrigeration & heat-pump coefficients of performance.
h₂ = h₁ + (h₂s − h₁)/η_is
Real compressor outlet enthalpy from isentropic outlet h₂s.
h₃ = h₄
Isenthalpic throttling across the expansion valve.
ṁ = Q_e,load / (h₁ − h₄)
Mass flow rate sized to meet the specified cooling load.
COP_Carnot = T_e / (T_c − T_e)
Theoretical maximum COP between the two temperature reservoirs (K).
η_II = COP_R / COP_Carnot
Second-law (exergetic) efficiency of the real cycle.

🔄 Schéma Animé du Cycle

🧊
ÉVAPORATEUR
T₄→T₁
⚙️
COMPRESSEUR
T₁→T₂
🌀
CONDENSEUR
T₂→T₃
🔻
DÉTENDEUR
T₃→T₄
Aspiration (vapeur froide)Refoulement (vapeur chaude)Liquide (tiède)Diphasique (froid)

🧠 Analyse d'Ingénierie

COP (réfrigération)
COP (pompe à chaleur)
COP de Carnot
Rendement Exergétique
%
Débit Massique
kg/s
Travail Compresseur
kW
Rejet Chaleur Cond.
kW
Taux de Compression
Temp. Refoulement
°C
Production Frigorifique
kJ/kg

📈 Diagramme T–s

🧮 Outils Avancés

Export current results for reports or spreadsheets.

Comment Utiliser le Simulateur VCRC

Suivez ces instructions étape par étape pour modéliser, résoudre et analyser les performances du cycle frigorifique.

Étape 1 — Choisir une Configuration ou un Fluide Frigorigène

Sélectionnez un préréglage rapide (Réfrigérateur Domestique, Climatisation A/C R410A, Ammoniac Industriel R717 ou CO₂ Basse Température) ou choisissez un fluide spécifique dans le menu déroulant des Fluides Frigorigènes. Le simulateur met à jour automatiquement les constantes physiques ($T_c, P_c, \omega, M, c_{p,f}, c_{p,g}, h_{fg,ref}$).

Étape 2 — Ajuster les Températures d'Évaporation et de Condensation

Utilisez le curseur de Temp. d'Évaporation ($T_e$) pour régler la température d'ébullition basse pression (plage : $-50^\circ\text{C}$ à $+15^\circ\text{C}$) et le curseur de Temp. de Condensation ($T_c$) pour le côté haute pression (plage : $+10^\circ\text{C}$ à $+70^\circ\text{C}$). Le simulateur calcule instantanément les pressions de saturation $P_1$ et $P_2$ via le modèle Lee–Kesler.

Étape 3 — Régler la Surchauffe, le Sous-refroidissement et le Rendement Isentropique

Ajustez la Surchauffe à l'Aspiration ($\Delta T_{SH}$) à la sortie de l'évaporateur ($T_1 = T_e + \Delta T_{SH}$) et le Sous-refroidissement du Liquide ($\Delta T_{SC}$) à la sortie du condenseur ($T_3 = T_c - \Delta T_{SC}$). Réglez le Rendement Isentropique ($\eta_{is}$) du compresseur.

Étape 4 — Examiner le Tableau d'État et les Diagrammes Thermodynamiques

Consultez le Tableau des Points d'État pour obtenir les valeurs numériques exactes de Pression ($P$), Température ($T$), Enthalpie ($h$), Entropie ($s$), Titre en Vapeur ($x$) et Volume Spécifique ($v$). Observez les tracés sur les diagrammes $P-h$ et $T-s$.

Étape 5 — Comparer les Fluides, Balayage Paramétrique et Exportation

Ouvrez le panneau d'Outils Avancés pour comparer deux fluides côte à côte, générer des courbes de tendance de $\text{COP}_R$ en fonction de $T_e$ ou $T_c$, et exporter les résultats en CSV ou PNG.

Sultan Saudagar — Ingénieur Mécanicien
Écrit et Révisé par
Sultan SaudagarIngénieur Mécanicien & Spécialiste en Thermodynamique

Ingénieur mécanicien spécialisé dans la conception de systèmes thermiques, la mécanique des fluides numérique (CFD), le génie frigorifique et le génie climatique.

1. Qu'est-ce que le Cycle Frigorifique à Compression de Vapeur ?

Le Cycle Frigorifique à Compression de Vapeur (VCRC) est le procédé thermodynamique le plus utilisé et le plus commercialisé de l'histoire pour la production de froid artificiel et le pompage de chaleur. Présent dans des millions de réfrigérateurs ménagers, d'entrepôts frigorifiques commerciaux, de climatiseurs split et de groupes d'eau glacée industriels (chillers), le cycle absorbe de la chaleur dans une source froide et la rejette vers un milieu chaud grâce au travail mécanique ou électrique fourni au compresseur.

Contrairement aux systèmes à absorption ou thermoélectriques, le Cycle Frigorifique à Compression de Vapeur repose sur le changement de phase continu d'un fluide frigorigène circulant entre liquide et vapeur. En vaporisant le fluide à basse pression (et basse température de saturation) et en le condensant à haute pression, le système atteint des Coefficients de Performance ($\text{COP}$) élevés compris généralement entre $2.5$ et plus de $6.0$.

Pourquoi les Ingénieurs et Étudiants Utilisent des Simulateurs VCRC ?

La conception d'équipements frigorifiques exige de trouver un compromis thermodynamique optimal. Le calcul manuel des propriétés d'état dans la zone diphasique est long et sujet aux erreurs. Les ingénieurs utilisent des simulateurs interactifs pour :

  • Cartographier les Points d'État : Calculer instantanément l'enthalpie ($h$), l'entropie ($s$), la température ($T$), la pression ($P$), le titre en vapeur ($x$) et le volume spécifique ($v$).
  • Visualiser les Tracés Thermodynamiques : Tracer la boucle du cycle sur les diagrammes logarithmiques Pression-Enthalpie ($P-h$) et Température-Entropie ($T-s$) avec la courbe de saturation.
  • Quantifier les Irréversibilités Réelles : Évaluer l'impact du rendement isentropique du compresseur ($\eta_{is}$), de la surchauffe d'aspiration ($\Delta T_{SH}$) et du sous-refroidissement ($\Delta T_{SC}$).
  • Comparer les Fluides Frigorigènes : Évaluer les Hydrofluorocarboures (HFCs comme le R134a, R410A) par rapport aux fluides naturels (R717 Ammoniac, R744 CO₂) et aux Hydrofluorooléfines à faible GWP (HFOs comme le R1234yf).

Exemple d'Ingénierie Industrielle : Groupe Frigorifique de Supermarché

Considérons une centrale frigorifique de supermarché nécessitant $3.5\text{ kW}$ de puissance frigorifique nette pour maintenir des produits congelés à $-10^\circ\text{C}$. Le rejet de chaleur s'effectue dans l'air extérieur à $+40^\circ\text{C}$ avec du R134a. La pression d'aspiration s'élève à $201.7\text{ kPa}$ et la pression de refoulement à $1016.6\text{ kPa}$ (taux de compression de $5.04$). Avec $5\text{ K}$ de surchauffe, $3\text{ K}$ de sous-refroidissement et un rendement isentropique de $75\%$, le débit massique est de $\dot{m} = 0.0235\text{ kg/s}$, consommant $1.168\text{ kW}$ de puissance électrique pour un $\text{COP}_R = 2.996$.

2. Comment Fonctionne le Simulateur VCRC ? (Moteur Physique et Mathématique)

Notre simulateur intègre un moteur de calcul de cycle subcritique monotagé fonctionnant en unités SI ($K, \text{kPa}, \text{kJ/kg}, \text{kJ/kg}\cdot\text{K}, \text{kg/s}$). Le moteur combine la thermodynamique des états correspondants avec des modèles d'échelle de propriétés.

Les Quatre Processus Majeurs du Cycle VCRC

  1. Processus 1 → 2 (Compression Réelle Non Idéale) : La vapeur surchauffée basse pression sortant de l'évaporateur à l'état 1 entre dans le compresseur. Le travail mécanique transforme le gaz en vapeur surchauffée haute pression à l'état 2 ($s_2 > s_1$).
  2. Processus 2 → 3 (Rejet de Chaleur Isobare au Condenseur) : Le gaz haute pression entre dans le condenseur, se désurchauffe jusqu'à la température de saturation $T_c$, se condense en liquide et se sous-refroidit de $\Delta T_{SC}$ jusqu'à l'état 3 en rejetant la chaleur $Q_c$.
  3. Processus 3 → 4 (Détente Isenthalpique au Détendeur) : Le liquide sous-refroidi haute pression se détend de manière irréversible à travers un détendeur jusqu'à la pression d'évaporation $P_1$. L'enthalpie reste constante ($h_4 = h_3$), provoquant une chute de température jusqu'à $T_e$ et une vaporisation partielle (titre $x_4$).
  4. Processus 4 → 1 (Absorption de Chaleur Isobare à l'Évaporateur) : Le fluide diphasique froid absorbe la chaleur $Q_e$ de l'enceinte refroidie, se vaporise complètement et gagne la surchauffe $\Delta T_{SH}$ jusqu'à l'état 1.

Modèles et Formules Fondamentales

1. Corrélation de Pression de Saturation de Lee–Kesler

La pression de saturation $P_{sat}(T)$ est calculée à partir de la température réduite $T_r = T / T_c$, de la pression critique $P_c$ et du facteur acentrique $\omega$ :

$$\ln P_r = f^(0)(T_r) + \omega \cdot f^(1)(T_r)$$ $$f^(0)(T_r) = 5.92714 - \frac{6.09648}{T_r} - 1.28862 \ln T_r + 0.169347 T_r^6$$ $$f^(1)(T_r) = 15.2518 - \frac{15.6875}{T_r} - 13.4721 \ln T_r + 0.43577 T_r^6$$

2. Corrélation de Chaleur Latente de Watson

La chaleur latente de vaporisation $h_{fg}(T)$ s'échelonne depuis l'état de référence $T_{ref} = 233.15\text{ K} (-40^\circ\text{C})$ selon :

$$h_{fg}(T) = h_{fg,ref} \cdot \left[ \frac{T_c - T}{T_c - T_{ref}} \right]^{0.38}$$

3. Intégration des Propriétés de Phase

  • Liquide Saturé :
    $$h_f(T) = c_{p,f} (T - T_{ref}), \quad s_f(T) = c_{p,f} \ln\left(\frac{T}{T_{ref}}\right)$$
  • Vapeur Saturée :
    $$h_g(T) = h_f(T) + h_{fg}(T), \quad s_g(T) = s_f(T) + \frac{h_{fg}(T)}{T}$$
  • Vapeur Surchauffée :
    $$h(T,P) = h_g(T_{sat}) + c_{p,g} (T - T_{sat}), \quad s(T,P) = s_g(T_{sat}) + c_{p,g} \ln\left(\frac{T}{T_{sat}}\right)$$

3. Fonctionnalités Interactives du Simulateur VCRC

🎛️ Curseurs Temps Réel

Re-calcul immédiat des tables d'état, diagrammes et COP lors de la modification des paramètres.

🔄 Schéma SVG Animé

Schéma animé en direct représentant le compresseur, le ventilateur et la circulation du fluide.

📈 Diagrammes $P-h$ & $T-s$

Graphiques haute définition avec courbe de saturation, cycle tracé et chemin isentropique.

⚖️ Comparateur de Fluides

Tableau d'évaluation comparative côte à côte de deux fluides frigorigènes.

📊 Balayage Paramétrique

Génération de courbes d'évolution du $\text{COP}_R$ selon la température de condensation ou d'évaporation.

⬇️ Exportation Multi-format

Téléchargez les tableaux au format CSV ou enregistrez les diagrammes sous forme d'images PNG.

4. Paramètres d'Entrée (Référence Technique)

Paramètre d'EntréeSymboleUnitésPlage TypiqueSignification Physique et Impact
Fluide FrigorigèneR134a, R22, R410A, R717, R744, R1234yfDétermine la masse molaire, les coordonnées critiques, la chaleur latente et les impacts écologiques (ODP, GWP).
Température d'Évaporation$T_e$$^\circ\text{C}$ / $\text{K}$$-50^\circ\text{C}$ à $+15^\circ\text{C}$Fixe la pression d'ébullition $P_1$. Une température $T_e$ plus haute augmente la densité à l'aspiration et améliore le $\text{COP}_R$.
Température de Condensation$T_c$$^\circ\text{C}$ / $\text{K}$$+10^\circ\text{C}$ à $+70^\circ\text{C}$Fixe la pression haute $P_2$. Une température $T_c$ plus basse réduit le travail du compresseur et améliore le $\text{COP}_R$.
Surchauffe à l'Aspiration$\Delta T_{SH}$$\text{K}$ / $^\circ\text{C}$$0\text{ K}$ à $25\text{ K}$Garantit la présence de vapeur $100\%$ sèche à l'entrée du compresseur ($T_1 = T_e + \Delta T_{SH}$), évitant les coups de liquide.
Sous-refroidissement du Liquide$\Delta T_{SC}$$\text{K}$ / $^\circ\text{C}$$0\text{ K}$ à $20\text{ K}$Refroidit le liquide en dessous de la saturation ($T_3 = T_c - \Delta T_{SC}$), augmentant la puissance frigorifique sans travail supplémentaire.
Rendement Isentropique$\eta_{is}$$\%$$40\%$ à $95\%$Rapport entre le travail isentropique idéal et le travail réel du compresseur ($\eta_{is} = \frac{h_{2s} - h_1}{h_2 - h_1}$). Un meilleur rendement réduit la puissance $W$.
Charge Frigorifique$Q_{e,load}$$\text{kW}$$0.5\text{ kW}$ à $50\text{ kW}$Puissance thermique à absorber dans l'enceinte. Détermine directement le débit massique $\dot{m}$.

5. Paramètres de Sortie et Performance

Grandeur de SortieSymboleUnitésFormule MathématiqueInterprétation en Ingénierie
Production Frigorifique Spécifique$q_e$$\text{kJ/kg}$$q_e = h_1 - h_4$Chaleur nette absorbée dans l'évaporateur par kilogramme de fluide circulant.
Débit Massique$\dot{m}$$\text{kg/s}$$\dot{m} = \frac{Q_{e,load}}{h_1 - h_4}$Débit de fluide frigorigène nécessaire pour satisfaire la charge thermique.
Puissance du Compresseur$W_{comp}$$\text{kW}$$W = \dot{m}(h_2 - h_1)$Consommation d'énergie mécanique/électrique réelle fournie à l'arbre du compresseur.
Chaleur Rejetée au Condenseur$Q_c$$\text{kW}$$Q_c = \dot{m}(h_2 - h_3)$Puissance thermique totale rejetée au milieu extérieur. Vérifie $Q_c = Q_e + W$.
COP Frigorifique$\text{COP}_R$$\text{COP}_R = \frac{Q_e}{W} = \frac{h_1 - h_4}{h_2 - h_1}$Rapport entre le froid utile produit et l'énergie consommée. Metric clé d'efficacité.
COP Pompe à Chaleur$\text{COP}_{HP}$$\text{COP}_{HP} = \frac{Q_c}{W} = \text{COP}_R + 1$Performance lors du fonctionnement en mode chauffage. Supérieur au $\text{COP}_R$ de $1.0$.
Limite de Carnot$\text{COP}_{Carnot}$$\text{COP}_{Carnot} = \frac{T_e}{T_c - T_e}$Limite théorique maximale d'un cycle réversible fonctionnant entre $T_e$ et $T_c$ (en Kelvin).
Rendement Exergétique$\eta_{II}$$\%$$\eta_{II} = \frac{\text{COP}_R}{\text{COP}_{Carnot}} \times 100$Efficacité selon le second principe mesurant l'écart par rapport à la perfection de Carnot.
Taux de Compression$PR$$PR = \frac{P_2}{P_1}$Rapport entre la pression de condensation et la pression d'évaporation.
Température de Refoulement$T_2$$^\circ\text{C}$$T_2 = T_c + \frac{h_2 - h_g(T_c)}{c_{p,g}}$Température maximale en sortie de compresseur. Une valeur trop élevée dégrade l'huile.

6. Équations d'Ingénierie et Dérivations

6.1 Bilan Massique et Énergétique à l'Évaporateur

D'après le Premier Principe de la thermodynamique en régime permanent :

$$\dot{Q}_e = \dot{m} (h_1 - h_4) \implies \dot{m} = \frac{Q_{e,load}}{h_1 - h_4}$$

6.2 Compression Réelle et Rendement Isentropique

Pour une compression isentropique idéale ($s_{2s} = s_1$), l'enthalpie idéale de sortie $h_{2s}$ s'écrit :

$$T_{2s} = T_c \cdot \exp\left[ \frac{s_1 - s_g(T_c)}{c_{p,g}} \right]$$ $$h_{2s} = h_g(T_c) + c_{p,g} (T_{2s} - T_c)$$

L'enthalpie réelle $h_2$ est déterminée par le rendement isentropique $\eta_{is}$ :

$$\eta_{is} = \frac{h_{2s} - h_1}{h_2 - h_1} \implies h_2 = h_1 + \frac{h_{2s} - h_1}{\eta_{is}}$$ $$W_{comp} = \dot{m} (h_2 - h_1)$$

6.3 Bilan de Chaleur au Condenseur

$$Q_c = \dot{m} (h_2 - h_3)$$ $$Q_c = Q_e + W_{comp} \implies \dot{m}(h_2 - h_3) = \dot{m}(h_1 - h_4) + \dot{m}(h_2 - h_1)$$

6.4 Détente Isenthalpique et Titre en Vapeur

La détente dans le détendeur est isenthalpique ($h_4 = h_3$). Le titre en vapeur $x_4$ est donné par :

$$x_4 = \frac{h_4 - h_{f,e}{h_{fg,e}}$$

6.5 Coefficients de Performance ($\text{COP}_R$ & $\text{COP}_{HP}$)

$$\text{COP}_R = \frac{Q_e}{W_{comp}} = \frac{h_1 - h_4}{h_2 - h_1}$$ $$\text{COP}_{HP} = \frac{Q_c}{W_{comp}} = \text{COP}_R + 1$$

6.6 Limite de Carnot et Rendement Exergétique

$$\text{COP}_{Carnot} = \frac{T_e}{T_c - T_e}, \quad \eta_{II} = \frac{\text{COP}_R}{\text{COP}_{Carnot}}$$

7. Exemple Numérique Résolu Étape par Étape

Énoncé du Problème

Un groupe frigorifique de supermarché fonctionnant au R134a opère avec une température d'évaporation $T_e = -10^\circ\text{C}$ ($263.15\text{ K}$) et de condensation $T_c = 40^\circ\text{C}$ ($313.15\text{ K}$). Surchauffe $\Delta T_{SH} = 5\text{ K}$, sous-refroidissement $\Delta T_{SC} = 3\text{ K}$, rendement isentropique $\eta_{is} = 75\%$, et charge frigorifique $Q_e = 3.5\text{ kW}$.

Résultats Calculés

  • Pressions : $P_1 = 201.7\text{ kPa}$, $P_2 = 1016.6\text{ kPa}$, taux de compression $PR = 5.04$.
  • État 1 : $h_1 = 254.35\text{ kJ/kg}$, entropie $s_1 = 0.9672\text{ kJ/(kg\cdot K)}$.
  • État 2 : $h_2 = 314.04\text{ kJ/kg}$, température de refoulement $T_2 = 80.47^\circ\text{C}$.
  • États 3 & 4 : $h_3 = h_4 = 109.34\text{ kJ/kg}$, titre en vapeur $x_4 = 0.322$.
  • Production frigorifique spécifique : $q_e = 145.01\text{ kJ/kg}$.
  • Débit massique : $\dot{m} = 0.02414\text{ kg/s}$ ($86.89\text{ kg/h}$).
  • Puissance compresseur : $W = 1.441\text{ kW}$.
  • COP Frigorifique : $\text{COP}_R = 2.429$.
  • COP Carnot : $\text{COP}_{Carnot} = 5.263$, Rendement Exergétique $\eta_{II} = 46.15\%$.

8. Tableau Comparatif des Fluides Frigorigènes Industriels

Fluide FrigorigèneType Chimique$T_c$ ($^\circ\text{C}$)$P_c$ ($\text{kPa}$)$M$ ($\text{g/mol}$)$h_{fg,ref}$ ($\text{kJ/kg}$)ODPGWPSécurité ASHRAEDomaine d'Application Typique
R134aHFC$101.06$$4059$$102.03$$225.9$01430A1Climatisation automobile, réfrigérateurs domestiques, chillers
R22HCFC$96.15$$4990$$86.47$$233.2$0.0551810A1Anciennes installations de climatisation (en élimination)
R410AMélange HFC$71.35$$4901$$72.58$$275.0$02088A1Climatisation split moderne, pompes à chaleur
R717 (NH₃)Ammoniac Naturel$132.35$$11333$$17.03$$1390.0$00B2LEntrepôts frigorifiques industriels, agroalimentaire
R744 (CO₂)CO₂ Naturel$30.98$$7377$$44.01$$322.0$01A1Supermarchés commerciaux en cycle transcritique
R1234yfHFO$94.70$$3382$$114.04$$191.0$0<1A2LNouvelle génération de climatisation automobile

9. Foire Aux Questions (FAQ Frigorifique)

Q1 : Quelle est la différence entre un réfrigérateur et une pompe à chaleur ?

Thermodynamiquement, les deux fonctionnent exactement selon le même cycle frigorifique à compression de vapeur. La différence réside dans l'objectif utile : le réfrigérateur extrait de la chaleur $Q_e$ pour refroidir une enceinte, tandis que la pompe à chaleur valorise la chaleur $Q_c$ rejetée au condenseur pour chauffer un espace. Ainsi, $\text{COP}_{HP} = \text{COP}_R + 1$.

Q2 : Pourquoi l'augmentation de la température de condensation réduit-elle le COP ?

Augmenter $T_c$ élève la pression haute $P_2$, ce qui accroît le travail de compression $(h_2 - h_1)$. Parallèlement, l'enthalpie à l'entrée de l'évaporateur $h_4 = h_3$ augmente, diminuant l'effet frigorifique $(h_1 - h_4)$. Plus de travail et moins de froid réduisent fortement le $\text{COP}_R$.

Q3 : Qu'est-ce qu'un coup de liquide en compression ?

L'aspiration de gouttelettes de liquide incompressibles dans le cylindre du compresseur provoque un choc mécanique violent qui détruit les clapets et les embiellages. Une surchauffe de $3\text{ K}$ à $8\text{ K}$ sécurise le fonctionnement.

10. Références Académiques & Ouvrages

  • Duminil, M. (2015). Thermodynamique appliquée aux systèmes frigorifiques. Éditions Pyc.
  • Rapin, P., & Jacquard, P. (2018). Formulaire du froid (15e éd.). Dunod.
  • ASHRAE. (2021). ASHRAE Handbook — Fundamentals. Atlanta, GA.

Outils de Simulation Associés

Découvrez d'autres simulateurs thermodynamiques et mécaniques en ligne gratuits

Tous les calculs du simulateur VCRC s'exécutent localement dans votre navigateur.