INGENIERIA MECANICA

EN MANTENIMIENTO
INDUSTRIAL

Termodinámica II

“Ciclo de refrigeración por compresión mecánica”

NOMBRE: Hans Muena.

José Luis Yañez.
CARRERA: Ingeniería Mecánica en Mantenimiento Industrial, Sección 823.
ASIGNATURA: Termodinámica II.
PROFESOR: Diohanic Bustos.
FECHA: 12/07/2018.
Tabla de contenido
Termodinámica II .................................................................................................................... 1
“Ciclo de refrigeración por compresión mecánica” ................................................................. 1
1 Introducción .................................................................................................................... 3
2 Objetivos .......................................................................................................................... 4
2.1 Objetivo general: .............................................................................................................. 4
2.2 Objetivos Específicos: ....................................................................................................... 4
3 Marco teórico. .................................................................................................................. 5
3.1 Ciclo de refrigeración: ....................................................................................................... 7
3.2 Características técnicas de los equipos: ................................................................................ 8
4 Desarrollo del tema. ......................................................................................................... 9
4.1 Método de Caso. ............................................................................................................... 9
4.1.1 Ciclo de refrigeracion estándar de compresión mecánica. ...................................................... 9
4.1.2 Ciclo de refrigeración de compresión doble en cascada........................................................ 15
5 Discusión........................................................................................................................ 26
6 Conclusión. .................................................................................................................... 27
7 Bibliografía. ................................................................................................................... 28
8 Webgrafía. ..................................................................................................................... 28

Termodinámica II Página 2
1 Introducción

 La elaboración del presente informe trata el análisis de un ciclo de refrigeración por

compresión mecánica de una cámara frigorífica, donde se describen los dispositivos
empleados para producir la energía frigorífica necesaria para el ciclo, en el cual se
utilizara refrigerante R-22.
 Se analizará un ciclo de compresión doble en cascada, con temperaturas de evaporación y
condensación específicas.

Termodinámica II Página 3
2 Objetivos

2.1 Objetivo general:

 Elaborar informe de análisis de ciclo estándar de compresión mecánica y ciclo de

compresión doble en cascada.

2.2 Objetivos Específicos:

 Especificar características técnicas de los equipos.

 Dibujar diagramas P – h y T – s.
 Definir los estados de cada fase del ciclo de compresión mecánica
 Calcular el trabajo de compresión.
 Calcular el calor de condensación.
 Calcular la relación de compresión.
 Calcular caudal volumétrico de refrigerante manejado por el compresor.
 El coeficiente de eficiencia energética de refrigeración.

Termodinámica II Página 4
3 Marco teórico.

Libro: “Fundamentos de la termodinámica técnica”.

Edición: 2ª edición / correspondiente a la 4ª edición original.
Autores: Michael J. Moran.
Howard N. Shapiro.
Constantes:
 Aceleración constante de la gravedad:
9,80665 𝑚/𝑠 2
g ={
32,174 𝑓𝑡/𝑠 2

 Presión atmosférica estándar:

1,01325 𝑏𝑎𝑟
1 atm = {
14,696 𝑙𝑏𝑓/𝑖𝑛.2

Formulas presentes en el informe:

 Trabajo de compresión:
𝐖𝐜
= 𝐡𝟐 − 𝐡𝟏
𝐦

 Calor transferido por unidad de masa de refrigerante:

𝑸𝒄
= 𝒉𝟏 − 𝒉𝟒
𝒎

 Calor de condensación:
𝐐𝐬
= 𝐡𝟐 − 𝐡𝟑
𝐦

Termodinámica II Página 5
  Coeficiente de eficiencia energética de refrigeración:

𝑄𝑒/𝑚 ℎ1 − ℎ4
β= =
𝑊𝑐/𝑚 ℎ2 − ℎ1

 Relación de compresión:

𝑃 𝑐𝑜𝑛𝑑
Rc =
𝑃 𝑒𝑣𝑎𝑝

 Caudal volumétrico de refrigerante:

Vr = m x V1

Termodinámica II Página 6
3.1 Ciclo de refrigeración:

Para comenzar el objetivo de un sistema de refrigeración es mantener una región fría a una
temperatura por debajo a la del entorno. Esto se consigue habitualmente empleando un sistema de
refrigeración que utiliza vapor.

Vamos a seguir al refrigerante a través de cada uno de los equipos del ciclo, comenzando por la
entrada del evaporador. El refrigerante entra en el evaporador como mezcla de líquido y vapor en
el estado 4. En el evaporador parte del refrigerante cambia de fase de líquido a vapor como
consecuencia del calor transferido del foco a la temperatura del refrigerante. La temperatura y la
presión del refrigerante permanecen constantes durante el proceso que va del estado 4 al estado 1.
El refrigerante se comprime adiabáticamente desde el estado 1, donde es una mezcla bifásica
liquido-vapor, hasta el estado 2, donde es vapor saturado. Durante este proceso la temperatura del
refrigerante se incrementa desde la Tf (temperatura de frio) hasta Tc (temperatura de calor), y la
presión también aumenta. El refrigerante pasa del compresor al condensador, donde cambia de
fase de vapor saturado a liquido saturado debido al calor transferido desde el foco de temperatura
Tc. La presión y la temperatura permanecen constantes desde el proceso que va desde el estado 2
al estado 3. El refrigerante vuelve a su estado de entrada en el evaporador después de su expansión
adiabática en una turbina. En este proceso del estado 3 al estado 4, la temperatura decrece desde
Tc a Tf, y hay un descenso de la presión.

Termodinámica II Página 7
3.2 Características técnicas de los equipos:

Proceso Maquina Características

1-2 Compresor Compresión isoéntropica del

refrigerante desde el estado 1
hasta la presión del condensador
en el estado 2.
2-3 Condensador Transferencia de calor desde el
refrigerante que fluye a presión
constante en el condensador. El
refrigerante sale como liquido
en el estado 3.
3-4 Válvula de expansión Proceso de estrangulación desde
el estado 3 hasta la mezcla
bifásica liquido – vapor en el
estado 4.
4-1 Evaporador Transferencia de calor hacia el
refrigerante que fluye a presión
constante en el evaporador hasta
completar el ciclo

Termodinámica II Página 8
4 Desarrollo del tema.

4.1 Método de Caso.

Una cámara frigorífica, demanda una potencia frigorífica de 100.000 (Frig/h) su temperatura de
evaporación debe ser de -30°C y su temperatura de condensación de 40°C. Si se pretende usar en
ambos casos refrigerante R-22. Dibujar diagramas de P-h y T-s. Calcular trabajo de compresión, el
calor de condensación, la relación de compresión, el caudal volumétrico de refrigerante manejado
por el compresor y el coeficiente de eficiencia energética de refrigeración.
Nota:
Para el caso de refrigeración en cascada la temperatura de evaporación del ciclo superior debe ser
inferior a la temperatura de condensación del ciclo inferior, a esta diferencia de temperaturas se le
llama solape. Si se utiliza un solape de 5°C y una temperatura intermedia.
t3 = 2,7°C y t5 = t8 = -2,3°C
La potencia frigorífica suministrada por todos los ciclos debe ser:
Qe = 100.000 frig/h = 116,3 kW.

4.1.1 Ciclo de refrigeracion estándar de compresión mecánica.

Componentes de un sistema de refrigeración por compresión de vapor:

Termodinámica II Página 9
Diagrama de t-s (Temperatura – Entropía) del ciclo de compresión mecánica.

Diagrama de P-h (Presion – entalpia) del ciclo de compresion mecanica:

Termodinámica II Página 10
Estado 1:

-Fase: Vapor saturado.

-Presión 1 (p1): 163,9 kPa.
-Temperatura 1 (t1): -30°C.
-Titulo (x): 1.
-Entalpia (h1): 392,5 kJ/kg.
-Entropía (s1) = sg =1,801 kJ/ (kg x K).
-Volumen específico (℣1): 0,1354 𝑚3 /kg.

Estado 2:

-Fase: Vapor sobre calentado.

-Presión 2 (p2): 1534,1 kPa.
-Temperatura 2 (t2): 77,92°C
-Titulo (x): no existe título porque está fuera de la curva.
-Entalpia (h2): 450,33 kJ/kg.
-Entropía (s2): s1 = sg = 1,801 kJ/ (kg x K).

 Interpolación para temperatura:

s (Entropía) t (°C)
1,806 80
1,801 t2
1,782 70
𝑡2 − 80 70 − 80
=
1,801 − 1,806 1,782 − 1,806

t2 = 77,92°C

 Interpolación para entalpia (h):

t (°C) h (entalpia)
80 452,1
77,92 h2
70 443,6
ℎ2 − 452,1 443,6 − 452,1
=
77,92 − 80 70 − 80

h2 = 450,33 kJ/kg.

Termodinámica II Página 11
Estado 3:

-Fase: Liquido comprimido..

-Presión 3 (p3): 1534,1 kPa.
-Temperatura 3 (t3): 40°C.
-Titulo (x): 0
-Entalpia (h3): 249,8 kJ/kg.
-Entropía (s3): 1,167 kJ/ (kg x K).

Estado 4:

-Fase: Mezcla bifásica.

-Presión 4 (p4): 163,9 kPa.
-Temperatura 4 (t4): -30°C
-Titulo (x): 37,10
-Entalpia (h4) = h3 = 249,8 kJ/kg.
-Entropía (s4): 1,2141 kJ/ (kg x K).
-Entropía del líquido (sl): 0,868 kJ/ (kg x K).
-Entropía del gas (sg): 1,801 kJ/ (kg x K).
-Entalpia del líquido (hl): 165,6 (kJ/kg).
-Entalpia del gas (hg): 392,5 (kJ/kg).

 Calculo de título de mezcla bifásica (x = %):

h4 = hl + x (hg – hl)
249,8 = 165,6 + x (392,5 – 165,6)

 Despejando X:

249,8 − 165,6
x=
226,9

x = 0,371088 x 100
x = 37,10%

 Calculo de entropía (s4):

s4 = sl + x (sg – sl)
s4 = 0,868 + 0,3710 (1,801 – 0,868)
s4 = 1,2141 kJ/ (kg x K)

Termodinámica II Página 12
  Calculo del trabajo de compresión (Wc):
𝑊𝑐
= ℎ2 − ℎ1
𝑚

 Donde m es el flujo másico del refrigerante.

 Tenemos la fórmula del Calor transferido por unidad de masa de refrigerante
(Qc/m):
𝑄𝑐
= ℎ1 − ℎ4
𝑚

 Despejamos m:
𝑄𝑐
=𝑚
(ℎ1−ℎ4)

 Donde Qc es el calor es la potencia frigorífica en kW.

116,3𝑘𝑊
=𝑚
(392,5 − 249,8)

0,8149 𝑘𝑔/𝑠 = 𝑚

 Teniendo el flujo másico se calcula Wc:

𝑊𝑐
= ℎ2 − ℎ1
𝑚

𝑊𝑐 = (ℎ2 − ℎ1) 𝑥 𝑚
𝑊𝑐 = (450,33 − 392,5) 𝑥 0,8149
𝑊𝑐 = 47,12𝑘𝑊

 Calculo del calor de condensación (Qs):

𝑄𝑠
𝑚
= ℎ2 − ℎ3

𝑄𝑠 = (ℎ2 − ℎ3)𝑥 𝑚
𝑄𝑠 = (450,33 − 249,8)𝑥 0,8149
𝑄𝑠 = 163,4𝑘𝑊
Termodinámica II Página 13
  Relación de compresión:

𝑃 𝑐𝑜𝑛𝑑 (𝑃2)
Rc =
𝑃 𝑒𝑣𝑎𝑝 (𝑃4)

1534,1𝑘𝑃𝑎
Rc =
163,9𝑘𝑃𝑎

Rc = 9,359

 Caudal volumétrico de refrigerante:

Vr = m x V1

kg
Vr = 0,8149 x 0,1354𝑚3/kg
s

Vr = 0,1103 𝑚3 /𝑠 = 397,08 𝑚3 /ℎ

 Calculo del coeficiente de eficiencia energética de refrigeración:

𝑄𝑒/𝑚 ℎ1 − ℎ4
β= =
𝑊𝑐/𝑚 ℎ2 − ℎ1

142,71 392,5 − 249,8

β= =
57,82 450,33 − 392,5

β = 2,4681

Termodinámica II Página 14
 4.1.2 Ciclo de refrigeración de compresión doble en cascada.

Ciclo de refrigeración en cascada por compresión de vapor.

Termodinámica II Página 15
 Diagrama P –h del ciclo de compresión doble en cascada:

Diagrama T –s del ciclo de compresión doble en cascada:

Termodinámica II Página 16
Estado 1:

-Fase: Vapor saturado.

-Presión 1 (p1): 163,9 kPa.
-Temperatura 1 (t1): -30°C.
-Titulo (x): 1.
-Entalpia (h1): 392,5 kJ/kg.
-Entropía (s1) = sg =1,801 kJ/ (kg x K).
-Volumen específico (℣1): 0,1354 𝑚3 /kg.

Estado 2:

-Fase: Vapor sobre calentado.

-Presión 2 (p2): 1534,1 kPa.
-Temperatura 2 (t2): 77,92°C
-Titulo (x): no existe título porque está fuera de la curva.
-Entalpia (h2): 450,33 kJ/kg.
-Entropía (s2): s1 = sg = 1,801 kJ/ (kg x K).

Estado 3:

-Fase: Liquido comprimido..

-Presión 3 (p3): 1534,1 kPa.
-Temperatura 3 (t3): -2,7°C.
-Titulo (x): 0
-Entalpia (h3): 197,03 kJ/kg.
-Entropía (s3): 0.986 kJ/ (kg x K).

 Interpolación para entropía:

T (°C) s (Entropia)
0 0,998
-2,7 s3
-10 0,955
ℎ3 − 0,998 0,955 − 0,998
=
−2,7 − 0 −10 − 0

s3 = 0.986 kJ/ (kg x K).

Termodinámica II Página 17
  Interpolación para entalpia (h3):

t (°C) h (entalpia)
0 200,2
-2,7 h3
-10 188,6
ℎ3 − 200,2 188,6 − 200,2
=
−2,7 − 0 −10 − 0

h3 = 197,03 kJ/kg .

Estado 4:

-Fase: Mezcla bifásica.

-Presión 4 (p4): 163,9 kPa.
-Temperatura 4 (t4): -30°C
-Titulo (x): 13,32%
-Entalpia (h4) = h3 = 197,03 kJ/kg.
-Entropía (s4): 0,992 kJ/ (kg x K).
-Entropía del líquido (sl): 0,868 kJ/ (kg x K).
-Entropía del gas (sg): 1,801 kJ/ (kg x K).
-Entalpia del líquido (hl): 165,6 (kJ/kg).
-Entalpia del gas (hg): 392,5 (kJ/kg).

 Calculo de título de mezcla bifásica (x = %):

h4 = hl + x (hg – hl)
197,03 = 165,6 + x (392,5 – 165,6)

 Despejando X:

197,03 − 165,6
x=
226,9

x = 0,1332 x 100
x = 13,32%

 Calculo de entropía (s4):

s4 = sl + x (sg – sl)
s4 = 0,868 + 0,1332 (1,801 – 0,868)
s4 = 0,992 kJ/ (kg x K)

Termodinámica II Página 18
Estado 5:

-Fase: Vapor saturado.

-Presión 5 (p5): 461,98 kPa.
-Temperatura 5 (t5): -2,3°C.
-Titulo (x): 1.
-Entalpia (h5): 404,09 kJ/kg.
-Entropía (s5) = sg =1,75345 kJ/ (kg x K).
-Volumen específico (℣5): 1289,3 𝑚3 /kg.

 Interpolación para presión:

T (°C) P (Presión)
-2 466,5
-2,3 P5
-4 0436,4
𝑝5 − 466,5 436,4 − 466,5
=
−2,3 + 2 −4 + 2

P5 = 461,98 kPa.

 Interpolación para entalpia:

P (presión) h (entalpia)
466,5 404,2
461,98 h5
436,4 403,5
ℎ5 − 404,2 403,5 − 404,2
=
461,98 − 466,5 436,4 − 466,5

P5 = 404,09 kJ/kg..

 Interpolación para entalpia:

P (presión) s (Entropia)
466,5 1,753
461,98 s5
436,4 1,756
𝑠5 − 1,753 1,756 − 1,753
=
461,98 − 466,5 436,4 − 466,5

s5 = 1,7534 kJ/(kg x K).

Termodinámica II Página 19
Estado 6:

-Fase: Vapor sobre calentado.

-Presión 6 (p6): 1534,1 kPa.
-Temperatura 6 (t6): 59,03°C
-Titulo (x): no existe título porque está fuera de la curva.
-Entalpia (h6): 434,027 kJ/kg.
-Entropía (s6): s5 = 1,75345 kJ/ (kg x K).

 Interpolación para temperatura:

s (Entropía) t (°C)
1,756 60
1,7534 t6
1,729 50
𝑡6 − 60 50 − 60
=
1,7534 − 1,756 1,729 − 1,756

t6 = 59,03°C

 Interpolación para entalpia (h6):

T (°C) h (entalpia)
60 434,9
59,03 h6
50 425,9
ℎ6 − 434,9 425,9 − 434,9
=
59,03 − 60 50 − 60

h6 = 434,027 kJ/kg.

Termodinámica II Página 20
Estado 7:

-Fase: Liquido comprimido..

-Presión 7 (p7): 1534,1 kPa.
-Temperatura 7 (t7): 40°C.
-Titulo (x): 0
-Entalpia (h7): 249,8 kJ/kg.
-Entropía (s7): 1,167 kJ/ (kg x K).

Estado 8:

-Fase: Mezcla bifásica.

-Presión 8 (p8): 461,98 kPa.
-Temperatura 8 (t8): -2,3°C
-Titulo (x): 25,40%
-Entalpia (h8) = h7 = 249,8 kJ/kg.
-Entropía (s8): 1,179 kJ/ (kg x K).
-Entropía del líquido (sl): 0,9842 kJ/ (kg x K).
-Entropía del gas (sg): 1,75345 kJ/ (kg x K).
-Entalpia del líquido (hl): 197,25 (kJ/kg).
-Entalpia del gas (hg): 404,095 (kJ/kg).

 Interpolación para entalpia del líquido (hl):

T (°C) hl (entalpia
liquido)
-2 197,6
-2,3 hl
-4 195,3

ℎ𝑙 − 197,6 195,3 − 197,6

=
−2,3 + 2 −4 + 2

hl = 19725 kJ/kg.

Termodinámica II Página 21
  Interpolación para entalpia del gas (hg):

T (°C) hg (entalpia
gas)
-2 404,2
-2,3 hg
-4 403,5

ℎ𝑔 − 404,2 403,5 − 404,2

=
−2,3 + 2 −4 + 2

hg = 404,095 kJ/kg.

 Calculo de título de mezcla bifásica (x = %):

h8 = hl + x (hg – hl)
249,8 = 197,25 + x (404,095 – 197,25)

 Despejando X:

249,8 − 197,25
x=
(404,095 − 197,25

x = 0,2540 x 100
x = 25,40%

 Interpolación para entropía del gas (sg):

T (°C) sg (entropía
gas)
-2 1,753
-2,3 sl
-4 1,756

𝑠𝑔 − 1,753 1,756 − 1,753

=
−2,3 + 2 −4 + 2

sg = 1,75345 kJ/(kg x K).

Termodinámica II Página 22
  Interpolación para entropía del líquido (sl):

T (°C) sl (entropía
liquido)
-2 0,983
-2,3 sl
-4 0,991

𝑠𝑙 − 0,983 0,991 − 0,983

=
−2,3 + 2 −4 + 2

sl = 0,9842 kJ/(kg x K).

 Calculo de entropía (s4):

s4 = sl + x (sg – sl)
s4 = 0,9842 + 0,2540 (1,75345 – 0,9842)
s4 = 1,179 kJ/ (kg x K)

 Calculo del trabajo de compresión (Wc)

𝑊𝑐
= ℎ2 − ℎ1
𝑚

 Donde m es el flujo másico del refrigerante.

 Tenemos la fórmula del Calor transferido por unidad de masa de refrigerante
(Qc/m):
𝑄𝑐
= ℎ1 − ℎ4
𝑚

 Despejamos m: (Balance de energía sobre el evaporador)

𝑄𝑐
= 𝑚a
(ℎ1−ℎ4)

Termodinámica II Página 23
  Donde Qc es el calor es la potencia frigorífica en kW.

116,3𝑘𝑊
= 𝑚𝑎
(392,5 − 197,03 )

0,5949 𝑘𝑔/𝑠 = 𝑚𝑎

 Balance de energía sobre el intercambiador intermedio:

mb (h8-h5) = ma (h3-h2)
ℎ3−ℎ2
mb = ma ( )
ℎ8−ℎ5
197,03−450,33
mb = 0,5949 ( )
249,8−404,09
mb = 0,9766 kg/s

 Trabajo de compresión Wc:

Wc = ma (h2-h1) + mb (h6-h5)
Wc = 0,614 (450,33-392,5) + 0,9766 (425,9-404,09)
Wc = 56,02 kW.

 Calculo del calor de condensación (Qc):

𝑄𝑠 = (ℎ6 − ℎ7)𝑥 𝑚b
𝑄𝑠 = (425,9 − 249,8)𝑥 0,9766
𝑄𝑠 = 179,97𝑘𝑊

Termodinámica II Página 24
  Relación de compresión:
𝑃 𝑐𝑜𝑛𝑑 (𝑃2)
Rc =
𝑃 𝑒𝑣𝑎𝑝 (𝑃4)

1534,1𝑘𝑃𝑎
Rc =
163,9𝑘𝑃𝑎

Rc = 9,359

 Caudal volumétrico de refrigerante:

Vr = (ma x V1) + (𝑚𝑏 𝑥 𝑉5)

Vr = (0,614 x 0,1354) + (0.9766 𝑥 1289,3)

Vr = 1259,2 𝑚3 /𝑠 = 4533,1 𝑚3 /ℎ

 Calculo del coeficiente de eficiencia energética de refrigeración:

𝑄𝑐
β=
𝑊𝑐
116,3
β=
56,02

β = 2,076

Termodinámica II Página 25
5 Discusión

En el informe propuesto se busca si una cámara frigorífica es capaz de satisfacer la demanda de

una potencia de 100.000 [frig/horas] que equivalen a 116,3 kW, para ello lo que se propone es
realizar los cálculos pertinentes para cada caso.
En el ciclo estándar los resultados obtenidos arrojan que la cámara frigorífica propuesta no cumple
con la demanda, ya que la potencia de este es de 47,11 kW (40514,6 [frig/h]) estableciendo que
esta cámara no es acorde a la potencia designada.
En el caso del ciclo doble en cascada al igual que el ciclo anterior tampoco cumpliría con lo
propuesto ya que según los cálculos recabados se obtiene una potencia de 56,02 kW (48177,2
[frig/h]); Además se visualiza que el rendimiento del ciclo doble en cascada es mayor al ciclo
estándar con un rendimiento de 2,076 en el ciclo doble en cascada y 2,4681 en el ciclo estándar de
compresión mecánica.
Sin embargo cabe recalcar que tanto como el primer ciclo como el segundo revelan que esta
cámara frigorífica no es apta para lo que se le designa.

Termodinámica II Página 26
6 Conclusión.

En el informe técnico entregado se realizó la confección de diagramas correspondientes tanto de

P-H y T-S. Además se idéntica cada uno de los estados en ambos ciclos (ciclo estándar de
refrigeración y ciclo doble en cascada)
Finalmente se llega a la conclusión que ambos ciclos no son aptos para la demanda propuesta, ya
que ningún de los dos alcanzan siquiera la mitad de la potencia requerida, al ser así se debe
investigar más a fondo o elegir otra cámara de refrigeración realizando los cálculos y recabando
datos para corroborar en sí, si la cámara nuevamente propuesta es apta para el trabajo y la potencia
que se le dignan a concretar.

Termodinámica II Página 27
7 Bibliografía.

Libro: “Fundamentos de la termodinámica técnica”.

Edición: 2ª edición / correspondiente a la 4ª edición original.
Autores: Michael J. Moran.
Howard N. Shapiro.

8 Webgrafía.

https://lms.inacap.cl/pluginfile.php/1265155/mod_resource/content/1/ES3%20-
%20Método%20de%20Caso%20Termodinámica%20II.pdf

https://es.scribd.com/doc/142456174/FULL-Termodinamica-Tecnica-Moran-Shapiro

https://www.pdfdrive.net/termodinamica-cengel-7th-d33969081.html

https://termoaplicadaunefm.files.wordpress.com/2009/02/tema-2-ciclos-de-refrigeracion1.pdf

http://jfc.us.es/DESCARGAS/TF_GITI/PROBLEMAS_RESUELTOS_TF_2013.pdf

Termodinámica II Página 28