Solução Problema 3- Refrigeração por Absorção

Prof. Ricado C. Carpio

Termodinâmica II - IFMG

Exemplo 3 Um sistema de refrigeração de absorção de amônia-água, conforme mostrado na Figura 3, recebe calor de 100 kW de uma fonte de
calor a 450 K para manter um espaço resfriado a 293 K a uma pressão de 200 kPa. Tanto o absorvedor como o condensador liberam calor com
uma taxa de 120 kW. O fluido de trabalho entra no evaporador com um título de 0,165 e uma razão de massa de amônia de 0,909 e sai com uma
temperatura de 280 K e um título de 0,783. Considere a temperatura e a pressão ambiente como sendo 298 K e 100 kPa, respectivamente.

a) Escreva as equações de balanço de massa, energia, entropia e exergia para cada componente do resfriador de absorção e de todo o ciclo.

b) Determine a taxa de transferência de calor que ocorre no evaporador.

c) Determine o COP en e o COP ex deste ciclo de refrigeração.

d) Determine a vazão de massa que passa pelo evaporador.

e) Determine a taxa de geração de entropia deste evaporador.

f) Determine a taxa de destruição de exergia deste evaporador.

Solução:
Dados:

˙
Q = 100kW
ge

˙ ˙
Q + Q = 120kW
ab cn

P5 = P6 = 200kP a

x4 = 0.165 (Título)
x5 = 0.783 (Título)

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

To = Tcondensador = 298 K (temperatura ambiente (fonte quente))
Po = 100kP a (Pressão ambiente)
In [1]: import numpy as np

import CoolProp.CoolProp as CP # importa o pacote coolprop para o codigo

from CoolProp.CoolProp import PropsSI # utiliza o Sistema internacional

CP.set_reference_state('NH3','ASHRAE')
#CP.set_reference_state('NH3','DEF')

In [2]: #Dados do Problema:

P1 = 198*1000 # Pascal - pressão 1

P5 = 200*1000 # Pascal - Pressão 5 antes do evaporador

x5 = 0.165

P6 = 200*1000 # Pascal - Pressão 6 depois do evaporador

x6 = 0.783
T6 = 280 # Kelvin

Qge = 100 # KW - Calor gerador

Qcn = 120 # KW - Calor Absorvedor + Condensador

Wp = 0 # Watt - Trabalho da bomba

To = 298 # kelvin - temperatura no estado morto

Tev = 293 # kelvin - temperatura no evaporador
Tge = 450 # kelvin - temperatura no gerador

a) Escreva as equações de balanço de massa, energia, entropia e exergia para cada componente
do resfriador de absorção e de todo o ciclo.:

Processo 1 - 2 Gerador
MBE : ṁ 1 = ṁ2 + ṁ3

EBE : ṁ 1 h1
˙
+ Q
ge
= ṁ2 h2 + ṁ3 h3

˙
Qge
EnBE : ṁ 1 s1 +
˙
+ S gen,ge = ṁ2 s2 + ṁ3 s3
Ts,ge

T0
ExBE : ṁ 1 ex1
˙
+ Qge (1 −
˙
) = ṁ2 ex2 + ṁ3 ex3 + Exd,ge
Ts,ge
Processo 2 - 4 Condensador
MBE : ṁ 2 = ṁ4

EBE : ṁ 2 h2 = ṁ4 h4 + Q̇
cn

Q̇
EnBE : ṁ
cn
2 s2 + Ṡ gen,cn = ṁ4 s4 +
Tb,cn

T0
ExBE : ṁ 2 ex2 = ṁ4 ex4 + Q̇
cn
(1 − ) + Ėxd,cn
Tb,cn

Processo 4 - 5 Valvula 1
MBE : ṁ 4 = ṁ5

EBE : ṁ 4 h4 = ṁ5 h5

EnBE : ṁ 4 s4
˙
+ S gen,tv1 = ṁ5 s5

ExBE : ṁ 4 ex4
˙
= ṁ5 ex5 + Exd,tv1

Processo 5 - 6 evaporador
MBE : ṁ 5 = ṁ6

EBE : ṁ 5 h5
˙
+ Q
ev
= ṁ6 h6

˙
Q
EnBE : ṁ
ev ˙
5 s6 + + S gen,ev = ṁ6 s6
Ts,gev

T0
ExBE : ṁ 5 ex5
˙
+ Q
ev
(
˙
− 1) = ṁ6 ex6 + Exd,ev
Ts,ev

Processo 6 - 10 Absorvedor
MBE : ṁ 6 + ṁ10 = ṁ7

EBE : ṁ 6 h6
˙
+ ṁ10 h10 = ṁ7 h7 + Q
ab

˙
Q
EnBE : ṁ
ab
˙
6 s6 + ṁ10 s10 = +S gen,ab = ṁ7 s7 +
Tb,ab

T0
ExBE : ṁ 6 ex6
˙ ˙
+ ṁ10 ex10 = ṁ7 ex7 + Exd,ab + Q
ab
(1 − )
Tb,ab

Processo 7 - 8 Bomba
MBE : ṁ 7 = ṁ8

EBE : ṁ 7 h7
˙
+ W p = ṁ8 h8

EnBE : ṁ 7 s7
˙
+ S gen,p = ṁ8 s8

ExBE : ṁ 7 ex7 + Ẇ p = ṁ8 ex8 + Ėxd,p

Processo 9 - 10 Valvula 2
MBE : ṁ 9 = ṁ10

EBE : ṁ 9 h9 = ṁ10 h10

EnBE : ṁ 9 s9 + Ṡ gen,tv2 = ṁ10 s10

ExBE : ṁ 9 ex9
˙
= ṁ10 ex10 + Exd,tv2
 Processo 1 - 3/ 8 - 9 Trocador de Calor
MBE : ṁ 3 + ṁ8 = ṁ1 + ṁ9

EBE : ṁ 3 h3 + ṁ8 h8 = ṁ1 h1 + ṁ9 h9

EnBE : ṁ 3 s3
˙
+ ṁ8 s8 + S gen,he = ṁ1 s1 + ṁ9 s9

ExBE : ṁ 3 ex3
˙
+ ṁ8 ex8 = ṁ1 ex1 + ṁ9 ex9 + Exd,he

Ciclo completo

EBE : Q
˙
ge
˙
+ Q
ev
˙ ˙
+ Wp = Q
ab
˙
+ Q
cn

˙ ˙ ˙ ˙
Q Q Q Q
ge
EnBE :
ev ab cn
˙
+ + S gen,ov = +
Ts,ge Ts,ev Tb,ab Tb,cn

T0 T0 T0 T0
ExBE : Q
˙
ge
(1 −
˙
) + Q
ev
(
˙ ˙
− 1) + W p = Q
ab
(1 −
˙
) + Q
cn
(1 −
˙
) + Exd,ov
Ts,ge Ts,ev Tb,ab Tb,cn

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

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

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

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

sg = PropsSI('S','P',P5,'Q',1,'NH3')/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 = 94.5913 KJ/Kg

Entalpia do vapor saturado hg = 1420.0789 KJ/Kg
Entalpia do líquido-vapor saturado hfg = 1325.4876 KJ/Kg
Entropia do líquido saturado sf = 0.3875 kJ/Kg K
Entropia do vapor saturado sg = 5.5998 kJ/Kg K
Entropia do líquido-vapor saturado sfg = 5.2123 kJ/Kg K

In [4]: 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 = 313.30 kJ/kg

O valor da Entropia s5 = 1.2475 kJ/kg K

In [5]: #Calculando as entalpias e as entropias nos pontos 5 e 6

h5 = PropsSI('H','P',P5,'Q',x5,'NH3')/1000;
s5 = PropsSI('S','P',P5,'Q',x5,'NH3')/1000;

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))

h6 = PropsSI('H','P',P5,'Q',x6,'NH3')/1000;
s6 = PropsSI('S','P',P5,'Q',x6,'NH3')/1000;

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 h5 = 313.30 kJ/kg

O valor da Entropia s5 = 1.2475 kJ/kg K
O valor da Entalpia h6 = 1132.45 kJ/kg
O valor da Entropia s6 = 4.4687 kJ/kg K

b) Determine a taxa de transferência de calor que ocorre no evaporador.

In [6]: #Calculo do calor no evaporador usando a EBE do ciclo completo

Qev = (Qcn - Qge - Wp)

print("O valor do calor no evaporador Qev = ","{0:.2f} kW".format(Qev))

O valor do calor no evaporador Qev = 20.00 kW

c) Determine o COP en e o COP ex deste ciclo de refrigeração.

In [7]: #Calculo do COPen

COPen = Qev / Qge

COPex = (Qev*((To/Tev) - 1))/(Qge*(1 - (To/Tge)))

print("O valor do COPen = ","{0:.3f}".format(COPen))

print("O valor do COPex = ","{0:.4f}".format(COPex))

O valor do COPen = 0.200

O valor do COPex = 0.0101

d) Determine a vazão de massa que passa pelo evaporador.

In [8]: # Sabendo que m5 = m6 podemos usar a equação EBE do evaporador
# e depois encotrar as entalpias para os dois pontos

m5 = np.round (Qev/(h6 - h5),4)

print("O valor da vazão em massa no evaporador = ","{0:.4f} kg/s".format(m5))

O valor da vazão em massa no evaporador = 0.0244 kg/s

e) Determine a taxa de geração de entropia deste evaporador.

In [9]: #Aplicando a EnBE no evaporador

Sev = np.round (m5*(s6 - s5) - (Qev/Tev),4)

print("A geração de entropia no evaporador = ","{0:.3f} kW/K".format(Sev))

A geração de entropia no evaporador = 0.010 kW/K

f) Determine a taxa de destruição de exergia deste evaporador.

In [10]: #Rearanjando a equação ExBE e coloando em função das entalpias e entropias

Exev = m5*((h5 - h6) - To*(s5 - s6)) + Qev*((To/Tev) - 1)

print("A taxa de destruição de exergia = ","{0:.3f} kW/K".format(Exev))

A taxa de destruição de exergia = 3.776 kW/K

Prof. Ricardo C. Carpio