Solução Exemplo 4 - Refrigeração em Cascata

Prof. Ricado C. Carpio

Termodinâmica II - IFMG

Exemplo 4. Considere um sistema de refrigeração em cascata de dois estágios (Figura 1) operando entre os limites de pressão de 1,2 MPa e
200 kPa (0,2 MPa) com o refrigerante R134a como fluido de trabalho. A rejeição de calor do cilco inferior para o ciclo superior ocorre em um
trocador de calor contracorrente e adiabático no qual a pressão dos ciclos superior e inferior correspondem a 0,4 MPa e 0,5 MPa
respectivamente. Em ambos os ciclos, o refrigerante é líquido saturado, na saída do condensador, e vapor saturado, na entrada do
compressor, cuja eficiência isentrópica é de 80%. Considerando que a vazão mássica do refrigerante através do ciclo inferior é de 0,15 kg/s,
determine:

a) Escreva as equações de balanço de massa, energia, entropia para cada componente do ciclo em cascata.

b) A vazão mássica do refrigerante através do ciclo superior.

c) A taxa de remoção de calor do espaço refrigerado.

d) O COP desse refrigerador.

 Solução:
In [1]: import numpy as np

import CoolProp.CoolProp as CP

from CoolProp.CoolProp import PropsSI # importa o pacote coolprop para o codigo

CP.set_reference_state('R134a','ASHRAE')

In [60]: # Algoritmo de Interpolação:

def lin_interp(point, vec_x, vec_y):

del_x = vec_x[-1] - vec_x[0]
del_y = vec_y[-1] - vec_y[0]
m = del_y/del_x
b = vec_y[0] - m * vec_x[0]
return point * m + b

Dados:

P1 = 198kP a

P2 = P3 = 1200kP a

T3 = 500 C

P4 = 198 kP a

x4 = 0.20 (Título)

ṁ = 0, 051 kg/s

ηcompressor = 85 %

Tevaporador = 263 K (temperatura do evaporador (fonte fria))

To = Tcondensador = 283 K (temperatura ambiente (fonte quente))

Po = 100kP a (Pressão ambiente)

In [61]: #Dados do Problema:

P5 = 170*1000 # Pascal- Pressão 5

x5 = 0.165

P6 = 170*1000 # Pascal- Pressão 6

x6 = 0.78365

In [62]: #Calculando os valores da entalpia e entropia do Líquido-vapor saturado no ponto 5:

hf = PropsSI('H','P',P5,'Q',0,'R134a')/1000; #entalpia líquido saturado

hg = PropsSI('H','P',P5,'Q',1,'R134a')/1000 #entalpia vapor saturado
hfg = hg - hf
 sf = PropsSI('S','P',P5,'Q',0,'R134a')/1000; #entropia do líquido saturado
sg = PropsSI('S','P',P5,'Q',1,'R134a')/1000; #entropia do vapor saturado
sfg = sg - sf

print('Entalpia do líquido saturado hf = ' "{0:.4f}".format(hf),"KJ/Kg")

print('Entalpia do vapor saturado hg = ' "{0:.4f}".format(hg),"KJ/Kg")
print('Entalpia do líquido-vapor saturado hfg = ' "{0:.4f}".format(hfg),"KJ/Kg")

print('Entropia do líquido saturado sf = ' "{0:.4f}".format(sf),"kJ/Kg K")

print('Entropia do vapor saturado sg = ' "{0:.4f}".format(sg),"kJ/Kg K")
print('Entropia do líquido-vapor saturado sfg = ' "{0:.4f}".format(sfg),"kJ/Kg K")

Entalpia do líquido saturado hf = 78.1550 KJ/Kg

Entalpia do vapor saturado hg = 1415.0533 KJ/Kg
Entalpia do líquido-vapor saturado hfg = 1336.8983 KJ/Kg
Entropia do líquido saturado sf = 0.3226 kJ/Kg K
Entropia do vapor saturado sg = 5.6560 kJ/Kg K
Entropia do líquido-vapor saturado sfg = 5.3334 kJ/Kg K

In [63]: h5 = hf + x5*hfg
s5 = sf + x5*sfg

print("O valor da Entalpia h5 = ","{0:.2f} kJ/kg".format(h5))

print("O valor da Entropia s5 = ","{0:.4f} kJ/kg K".format(s5))

O valor da Entalpia h5 = 298.74 kJ/kg

O valor da Entropia s5 = 1.2026 kJ/kg K

In [67]: #Calculando os valores da entalpia e entropia do Líquido-vapor saturado no ponto 6:

hf = PropsSI('H','P',P6,'Q',0,'134a')/1000; #entalpia líquido saturado

hg = PropsSI('H','P',P6,'Q',1,'134a')/1000 #entalpia vapor saturado
hfg = hg - hf

sf = PropsSI('S','P',P6,'Q',0,'134a')/1000; #entropia do líquido saturado

sg = PropsSI('S','P',P6,'Q',1,'134a')/1000; #entropia do vapor saturado
sfg = sg - sf

print('Entalpia do líquido saturado hf = ' "{0:.4f}".format(hf),"KJ/Kg")

print('Entalpia do vapor saturado hg = ' "{0:.4f}".format(hg),"KJ/Kg")
print('Entalpia do líquido-vapor saturado hfg = ' "{0:.4f}".format(hfg),"KJ/Kg")

print('Entropia do líquido saturado sf = ' "{0:.4f}".format(sf),"kJ/Kg K")

print('Entropia do vapor saturado sg = ' "{0:.4f}".format(sg),"kJ/Kg K")
print('Entropia do líquido-vapor saturado sfg = ' "{0:.4f}".format(sfg),"kJ/Kg K")

CP.PropsSI('D','T',300,'P',101325,'HEOS::R32[0.697615]&R125[0.302385]')

Entalpia do líquido saturado hf = 78.1550 KJ/Kg

Entalpia do vapor saturado hg = 1415.0533 KJ/Kg
Entalpia do líquido-vapor saturado hfg = 1336.8983 KJ/Kg
Entropia do líquido saturado sf = 0.3226 kJ/Kg K
Entropia do vapor saturado sg = 5.6560 kJ/Kg K
Entropia do líquido-vapor saturado sfg = 5.3334 kJ/Kg K
Out[67]: 2.986886779635724

In [68]: h6 = hf + x6*hfg
s6 = sf + x6*sfg

print("O valor da Entalpia h6 = ","{0:.2f} kJ/kg".format(h6))

print("O valor da Entropia s6 = ","{0:.4f} kJ/kg K".format(s6))

O valor da Entalpia h6 = 1125.82 kJ/kg

O valor da Entropia s6 = 4.5021 kJ/kg K

In [69]: h6 - h5

Out[69]: 827.0721625672145

(a) Escrevendo as equações de Balanço de massa, energia, entropia e exergia para cada
processo.

Processo 2 - 3 - Condensador
EBM: ṁ 2 = ṁ3 = ṁ

EBE: ṁ 2 h2
˙
= ṁ3 h3 + Q
cond.

EBEn: ṁ 2 s2
˙ ˙
+ S ger,2−3 = ṁ3 s3 + Q
cond.
/Tcond.

EBEx: ṁ 2 ex2
˙
= ṁ3 ex3 + Q
cond.
˙
(1 − To /Tcond ) + Exd,2−3
 Processo 4 - 1 - Evaporador
EBM: ṁ 4 = ṁ1 = ṁ

EBE: ṁ 4 h4
˙
+ Q
ev.
= ṁ1 h1

EBEn: ṁ 4 s4
˙ ˙
+ S ger,4−1 + Q
ev.
/Tev. = ṁ1 s1

EBEx: ṁ 4 ex4
˙
+ Q
ev.
˙
(To /Tev. − 1) = ṁ1 ex1 + Exd,4−1

Processo 1 - 2 - Compressor
EBM: ṁ 1 = ṁ2 = ṁ

EBE: ṁ 1 h1
˙
+ W ent = ṁ2 h2

EBEn: ṁ 1 s1
˙
+ S g,1−2 = ṁ2 s2

EBEx: ṁ 1 ex1
˙ ˙
+ W ent = ṁ2 ex2 + Exd,1−2

Processo 3 - 4 - Válvula de Expansão

EBM: ṁ 3 = ṁ4 = ṁ

EBE: ṁ 3 h3 = ṁ4 h4 → h3 = h4 (Processo Isoentalpico)

EBEn: ṁ 3 s3
˙
+ S ger,3−4 = ṁ4 s4

EBEx: ṁ 3 ex3 = ṁ4 ex4 + Exd,3−4

Achando as Entalpias e Entropias em cada Ponto do Processo.

Ponto 1: Vapor Saturado

In [4]: ######### Valores de entrada para P1 ######################

# Estado 1: Dados de entrada P1 = 198 kPa = 198000 Pa

# Utilizando o CoolProp ou da tabela A-11 - Moran

h1 = (PropsSI('H', 'P', P1, 'Q', 1, 'R134a')/1000)

s1 = (PropsSI('S', 'P', P1, 'Q', 1, 'R134a')/1000)

print("A Entalpia Vapor Saturado h1 = ","{0:.2f}".format(h1),'kJ/kg')

print("A Entropia do Vapor Saturado s1 = ","{0:.3f}".format(s1),'kJ/kg k')

A Entalpia Vapor Saturado h1 = 244.32 kJ/kg

A Entropia do Vapor Saturado s1 = 0.938 kJ/kg k

Ponto 2: Líquido Comprimido

Achando o valor de h , que se encontra na região de líquido Superaquecido. Será Calculado

2s

com P e s = s :
2 2s 1

Encontrar h : Da tabela A-12 - Moran. Buscamos valores antes e depois de s

2 2s = 0, 938 , para realizar a
interpolação:

s h

0,9527 287,44

0,938 h2s

0,9164 275,52

In [5]: # Interpolando para h2s

s = [0.9164, 0.9527]
h = [275.52, 287.44]
h2s = lin_interp(0.938,s,h)
print("A Entalpia h2s = ","{0:.1f}".format(h2s),'kJ/kg')

A Entalpia h2s = 282.6 kJ/kg

Calculando o valor da entalpia h de: 2

 h2s − h1
ηcomp. =
h2 − h1

(h2s − h1 )
h2 = h1 +
ηcomp

In [6]: # Calculando h2

h2 = h1 + ((h2s - h1)/Eff_c)
print("A Entalpia h2 = ","{0:.1f}".format(h2),'kJ/kg')

A Entalpia h2 = 289.4 kJ/kg

Achando o valor de s , Real. Será Calculado com P e h :

2 2 2

Encontrar s : Da tabela A-12 - Moran. Buscamos valores antes e depois de h

2 2 = 289, 4 , para realizar a interpolação:

h s

298,96 0,9868

289,4 s2

287,44 0,9527

In [7]: # Interpolando para s2

h = [287.44, 298.96]
s = [0.9527, 0.9868]
s2 = lin_interp(289.4,h,s)
print("A Entropia s2 = ","{0:.4f}".format(s2),'kJ/kg K')

A Entropia s2 = 0.9585 kJ/kg K

Ponto 4: Líquido - Vapor

Ponto 3: Líquido Superaquecido para Líquido-Vapor

Achando o valor de s . Será Calculado com h

3 3 = h4 :
Encontrar s ; T : Da tabela A-11 - Moran. Buscamos valores antes e depois de h
3 3 3 = 79, 36 , para realizar a
interpolação:

h s T (K)

71,33 0,2723 288,74

79,36 s3 T3

79,48 0,2999 294,58

In [10]: # Interpolando para s3

h3 = h4
h = [79.48, 71.33]
s = [0.2999, 0.2723]
s3 = lin_interp(h4,h,s)
print("A Entropia s3 = ","{0:.4f}".format(s3),'kJ/kg K')

# Interpolando para T3
h = [79.48, 71.33]
temp = [294.58, 288.74]
T3 = lin_interp(h4,h,temp)
print("Temperatura T3 = ","{0:.2f}".format(T3),'K')
print("Temperatura T3 = ","{0:.2f}".format(T3-273),'C')

A Entropia s3 = 0.2995 kJ/kg K

Temperatura T3 = 294.49 K
Temperatura T3 = 21.49 C

(b) Determinando a taxa de remoção de calor pelo refrigerador e a entrada de energia para o
compressor.

(b1) A Taxa de remoção de Calor:

˙
Qcondensador = ṁ(h2 − h3 )
In [27]: Q_cond = m *(h2 - h3)
print("O Calor de Remoção no Condensador = ","{0:.2f}".format(Q_cond),'KW')

O Calor de Remoção no Condensador = 10.71 KW

(b2) O Trabalho do Compressor:

Ẇ ent = ṁ(h2 − h1 )

In [12]: W_ent = m *(h2 - h1)

print("O Trabalho do Compressor = ","{0:.2f}".format(W_ent),'KW')

O Trabalho do Compressor = 2.30 KW

(c) Encontrar a geração de entropia e as taxas de destruição de exergia para cada componente.

- A Geração de Entropia para o Condensador: Da EBEn:

˙ ˙
ṁ2 s2 + S ger,2−3 = ṁ3 s3 + Q /Tcond.
cond.

˙ ˙
S ger,2−3 = ṁ(s3 − s2 ) + Q /Tcond.
cond.

In [13]: Sger_cond = m*(s3 - s2) + Q_cond/T_cond

print("A Entropia Gerada no Condensador = ","{0:.4f}".format(Sger_cond),'KW/K')

A Entropia Gerada no Condensador = 0.0042 KW/K

- A Destruição da Exergia para o Condensador: Da EBEx:

˙
ṁ2 ex2 = ṁ3 ex3 + Q (1 − To /Tcond ) + Ėxd,2−3
cond

To
˙ ˙
Exd,2−3 = ṁ(ex2 − ex3 ) − Q (1 − )
cond
Tcond

ṁ(ex2 − ex3 ) = (h2 − h3 ) − To (s2 − s3 )

In [14]: # Calculando a variação de Exergia 2-3:

ex23 = (h2 - h3) - To*(s2 - s3)

print("A Variação da Exergia ex23 = ","{0:.2f}".format(ex23),'kJ/Kg')

A Variação da Exergia ex23 = 23.51 kJ/Kg

In [15]: Exd_cond = m*ex23 - Q_cond*(1-(To/T_cond))

print("A Destruição de Exergia no Condensador = ","{0:.2f}".format(Exd_cond),'kW')

A Destruição de Exergia no Condensador = 1.20 kW

- A Geração de Entropia para o Evaporador: Da EBEn:

˙
Q
˙ ev.
ṁ4 s4 + S ger,4−1 + = ṁ1 s1
Tev.

˙
Qev.
˙
S ger,4−1 = ṁ(s1 − s4 ) −
Tev.

Qev. = ṁ(h1 − h4 )

In [16]: # Calculando o Calor do Evaporador

Qev = m*(h1 - h4)

print("O calor do Evaporador = ","{0:.4f}".format(Qev),'kW')

O calor do Evaporador = 8.4131 kW

In [17]: # A entropia gerada no Evaporador

Sger_ev = m*(s1 - s4) - (Qev/T_ev)
print("A Entropia Gerada no Condensador = ","{0:.4f}".format(Sger_ev),'KW/K')

A Entropia Gerada no Condensador = 0.0005 KW/K

- A Destruição da Exergia para o Evaporador: Da EBEx:

 To
˙
ṁ4 ex4 + Q ( − 1) = ṁ1 ex1 + Ėxd,4−1
ev.
Tev.

To
˙ ˙
Exd,4−1 = ṁ(ex4 − ex1 ) + Q ( − 1)
ev.
Tev.

(ex4 − ex1 ) = (h4 − h1 ) − To (s4 − s1 )

In [18]: # Calculando a variação de Exergia 4-1:

ex41 = (h4 - h1) - To*(s4 - s1)

print("A Variação da Exergia ex41 = ","{0:.2f}".format(ex41),'kJ/Kg')

A Variação da Exergia ex41 = 12.67 kJ/Kg

In [19]: Exd_ev = m*ex41 + Qev*((To/T_ev) - 1)

print("A Destruição de Exergia no Evaporador = ","{0:.3f}".format(Exd_ev),'kW')

A Destruição de Exergia no Evaporador = 1.150 kW

- A Geração de Entropia no Compressor: Da EBEn:

˙
ṁ1 s1 + S ger,1−2 = ṁ2 s2

˙
S ger,1−2 = ṁ(s2 − s1 )

In [20]: Sger_comp = m*(s2 - s1)

print("A Entropia Gerada no Compressor = ","{0:.4f}".format(Sger_comp),'KW/K')

A Entropia Gerada no Compressor = 0.0010 KW/K

- A Destruição da Exergia para o Compressor: Da EBEx:

˙ ˙
ṁ1 ex1 + W ent = ṁ2 ex2 + Exd,1−2

Ėxd,1−2 = ṁ(ex1 − ex2 ) + Ẇ ent

(ex1 − ex2 ) = (h1 − h2 ) − To (s1 − s2 )

In [21]: # Calculando a variação de Exergia 1-2:

ex12 = (h1 - h2) - To*(s1 - s2)
print("A Variação da Exergia ex12 = ","{0:.2f}".format(ex12),'kJ/Kg')

A Variação da Exergia ex12 = -39.24 kJ/Kg

In [28]: Exd_comp = m*ex12 + W_ent

print("A Destruição de Exergia no Evaporador = ","{0:.3f}".format(Exd_comp),'kW')

A Destruição de Exergia no Evaporador = 0.296 kW

- A Geração de Entropia na Válvula de Expansão: Da EBEn:

˙
ṁ3 s3 + S ger,3−4 = ṁ4 s4

˙
S ger,3−4 = ṁ(s4 − s3 )

In [23]: Sger_ve = m*(s4 - s3)

print("A Entropia Gerada na Válvula de Expansão = ","{0:.4f}".format(Sger_ve),'KW/K')

A Entropia Gerada na Válvula de Expansão = 0.0006 KW/K

- A Destruição da Exergia para a Válvula de Expansão: Da EBEx:

ṁ3 ex3 = ṁ4 ex4 + Exd,3−4

Ėxd,3−4 = ṁ(ex3 − ex4 )

(ex3 − ex4 ) = (h3 − h4 ) − To (s3 − s4 )

In [24]: # Calculando a variação de Exergia 3-4:

ex34 = (h3 - h4) - To*(s3 - s4)
print("A Variação da Exergia ex34 = ","{0:.2f}".format(ex34),'kJ/Kg')

A Variação da Exergia ex34 = 3.06 kJ/Kg

In [25]: Exd_ve = m*ex34
print("A Destruição de Exergia na Válvula de Expansão = ","{0:.3f}".format(Exd_ve),'kW')

A Destruição de Exergia na Válvula de Expansão = 0.156 kW

(d) Calcular o coeficiente de desempenho energético (COPen) e exergético (COPex) do ciclo.

˙
Q
ev
COPenergetico =
Ẇ ent

˙
Ex Q ev

COPexergetico =
˙
W ent

In [26]: COP_en = Qev/W_ent

COP_ex = (((To/T_ev)-1)*Qev)/W_ent
print("O Coeficiente de Desempenho Energética = ","{0:.2f}".format(COP_en))
print("O Coeficiente de Desempenho Exergético = ","{0:.3f}".format(COP_ex))

O Coeficiente de Desempenho Energética = 3.66

O Coeficiente de Desempenho Exergético = 0.219

Prof. Ricardo C. Carpio