Proyecto final, ciclo Rankine con recalentamiento

Juan Camilo Orozco

Mariana Martínez

Universidad tecnológica de Pereira

Termodinámica II

Pereira-Risaralda

2022
Un ciclo Rankine de vapor de agua opera con un calentador regenerativo bajo las siguientes condiciones:

▪ Entrada del vapor a la turbina: 8+X MPa y 400+(Yx10) oC

▪ Presión del condensador: 7+X kPa

▪ Eficiencia isentrópica de la turbina: (80+Y) %

▪ Presión de extracción: 1 MPa

▪ Potencia neta requerida del ciclo: (100+ 10xX +Y) MW

El proceso de combustión se realiza con gas natural (considerado como metano) o carbón y se establece
un

exceso de aire de 10+X %. La combustión presenta productos ideales (CO2, H2O, N2 y O2).

En el condensador, el agua de enfriamiento ingresa a 25 oC y sale a 40 oC para ser dirigida hacia una
torre de

enfriamiento que debe retornar el agua a la fuente hídrica a 25 oC. El aire atmosférico se encuentra a 25
oC y

humedad relativa del 40%. Se estima que el aire debe abandonar la torre con una temperatura 35 oC y
una

humedad relativa del 100%.

Las condiciones atmosféricas son T0=25 oC y P0=101,325 kPa.

Determinar:

▪ El flujo másico de vapor que debe producir la caldera (kg/h)

▪ El flujo másico de vapor de la extracción en la turbina (kg/h)

▪ La eficiencia térmica del ciclo (%)

▪ La exergía destruida en la caldera (kW)

▪ La exergía destruida en la turbina (kW)

▪ La exergía destruida en el condensador (kW)

▪ El costo del combustible por hora ($/h)

▪ Relación aire-combustible

▪ Determine las emisiones de CO2 por hora (t/h)

▪ Estime el costo del kWh neto producido para satisfacer los costos de combustible ($/kWh)
▪ Determine el caudal volumétrico de aire atmosférico de entrada a la torre de enfriamiento (m3/h)

▪ Determine el flujo másico de agua de reposición de la torre de enfriamiento (kg/h)

▪ Representar el diagrama de flujo de la planta y el ciclo de vapor en los diagramas T-s y h-s.

Datos

 Entrada del vapor a la turbina: 15 [MPa] y 460 [‫ﹾ‬C]

 Presión del condensador 14 [kPa]
 Eficiencia isentrópica de la turbina (86) %
 Presión de extracción 1MPa
 Potencia neta requerida del ciclo (176) [MW]
 El proceso de combustión se realiza con gas natural (metano) o carbón como se establece un
exceso de aire de 17%
 La combustión se presenta productos ideales (CO2, H2O, N2, y O2)
 En el condensador, el agua de enfriamiento ingrese a 25 [C] y sale a 40 [C] para ser dirigida
hacia una torre de enfriamiento que debe retornar el agua a la fuente hídrica a 25 [C]
 El aire atmosférico se encuentra a 25 [C] y humedad relativa del 40%
 Se estima que el aire debe abandonar la torre con una temperatura a 35 [C] y una humedad
relativa del 100%
 Condiciones atmosféricas son T0= 25 [C] Y P0= 101,325 kPa

Ilustración ciclo Rankine

 1) Se inicio el estudio del ciclo a la entrada de la
turbina debido a que es ahí donde tenemos
todos los datos, y por ende podemos hallar
tanto entropía como entalpia real, teniendo
en cuenta que el punto 6 es vapor
sobrecalentado hallamos los datos en la tabla
A-6

Entrada del vapor a la turbina: 15 [MPa] y

460 [‫ﹾ‬C]

h6 = 3188.48 [kJ/kg]

S6=6.18432 [kJ/kg*K]

2) Pasamos al estado 1 es la entrada al condensador, y lo tomamos como un proceso isoentrópico

S6=S1

S1=6.18432 [kJ/kg*K]

P1=14 [KPa]

Estado: Se encuentra
en zona de mezcla

Sf@14[KPa]= 0.73376
[kJ/kg*K]

Sfg@14[KPa]= 7.30168
[kJ/kg*K]

Se debe utilizar la
ecuación de calidad

X1=S1-SF@6/SFG@6

X1=0.74
Procedemos a hallar la entalpia en condiciones ideales (h1i=hf+xhfg) haciendo uso de la tabla A-5

hf@14[KPa]= 219.114 [kJ/kg]

hfg@14[KPa]= 2376.26 [kJ/kg]

h1 = hf + Xhfg

h1i= 1977.5464 [KJ/kg]

ηt= wideal/wreal

86%=h6-h1r/h6-h1i

h1r= 2147.077104 [KJ/kg]

3) Pasamos al estado 7 que es la salida de la turbina, y lo tomamos como un proceso isoentrópico

S6=S7
S7=6.18432 [kJ/kg*K]

P7=1 [MPa]

P7=1000 [kPa]

X7=S7-SF@6/SFG@6

X7=0.909 (Se encuentra en zona de mezcla)

h7 = hf + Xhfg

h7i = 2593.78 [KJ/kg]

Ahora se puede hallar la entalpia en condiciones reales

ηt= wideal/wreal

86%=h6-h7r/h6-h7i

h7r=2677.038 [KJ/kg]
4) Proseguimos con el estado 2, entrada de la bomba P1=P2 en el que tenemos una presión de [14 kPa]
y la calidad es X2=0 debido a que se encuentra en liquido saturado

hf@P2=h2=219.114 [kJ/kg]

vf@P2=v2=0.0010132 [m³/kg]

Sf@P2=S2=0.73376 [kJ/kg*K]

5) Ahora se pasa al estado 3 salida de la bomba 2 P7=P3 en liquido saturado hallo entalpia desde
trabajo de bomba 2

P3=1000 [kPa]

S3=S2

Wbomba1,entrada=v2(P3-P2)

Wbomba1,entrada=0.0010132[m³/kg](1000-14) [kPa]

Wbomba1,entrada=0.9990152 [kJ/kg]

h3=h2+ Wbomba1,entrada

h3=220.104 [kJ/kg]

6) Estado 4

Liquido saturado
Presión de 1000 [kPa]
 hf@1000[kPa]=h4=762.51 [kJ/kg]
Sf@1000[kPa]=S4=2.1381 [kJ/kg*K]
vf@1000[kPa]=v4=0.001127 [m³/kg]

7) Estado 5
P5=P6=15 [MPa]
S5=S4=2.1381 [kJ/kg*K]
Wbomba2, entrada= v4 (P5-P4)
Wbomba2, entrada= 15.778 [kJ/kg]
h5=h4+ Wbomba2, entrada
h5=778.29 [kJ/kg]

▪ Representar el diagrama de flujo de la planta y el ciclo de vapor en los diagramas T-s y h-s.
▪ El flujo másico de vapor que debe producir la caldera (kg/h)

Ẇneto
Ẇneto = ṁ(Wneto) →ṁ = Wneto

Wneto = Wturbina, salida - Wbomba, entrada

Wturbina, salida = ηt * Ws,turbina, salida

Wbomba, entrada = (1-y)*Wbomba I, entrada + Wbomba II, entrada

Primera ley CAA = myh7+(1-y) mh3– mh4 = 0

Y = 0.2283

Wneto = ηt * [ (h5-h4) + (1-y)(h6-h1) ] - [ (1-y)*[V2(P3-P2)] + V4(P5-P4)]

→Wneto = 0,86[ (778.29-762.51)kJ/kg +[ (1 -0.2283)( 0.0010132 m³/kg (1000-14))kJ/kg)] +[(0.001127

m³/kg(15000-1000)kPa)]

→Wneto = 30.1197 kJ/kg

kJ
176000
Ẇneto s kg 3600 s
ṁ= Wneto = 30.1197 kJ = 5843.35 s * 1 h = 21036060 kg/h
kg

▪ El flujo másico de vapor de la extracción en la turbina (kg/h)

Wneto
mv=
Wturbina ( real ) −Wbomba(ideal)

W Turbina de alta

Primera ley turbina de alta

W turbina de alta = m*(h6-h7)

W Turbina de baja

Primera ley turbina de baja

W Turbina de baja = m*(h7-h1)

∑w = 959.31 [KW]

Wbomba1

Primera ley bomba 1

Wbomba1 = m(h2-h3) =0.99 [KW]

Wbomba2

Primera ley bomba 2

Wbomba2= m(h4-h5) =15.78 [KW]

∑W bomba=16.544 [KW]

176000 [ KPa ]
Flujo masico de vapor = = 186.685 [kg/s] = 672066 [kg/h]
( 959.31−16.544 ) [KW ]

▪ La eficiencia térmica del ciclo (%)

Wneta
ηt =
Calor adiconado

Calor adicionado = calor de caldera = m (h6-h5) = 449946.32 [KW]

176000 [kPa]
ηt = = 39.11%
449946.31[ KW ]

▪ La exergía destruida en la caldera (kW)

S6r=6.18432

S5r=2.1381

Th=1500 [K]

q= 1 [kg/s] (h6-h5) [kJ/kg] =2410.19

q
Exergia caldera = To*(s6-s5) - ( )
th
= 726.95 [kJ/kg] * 186.685 [kJ/kg] = 135710.66kJ/s
▪ La exergía destruida en la turbina (kW)

Exergía turbina de alta=wreal-wreversible=9923,66 kw

W real=m*(h6r-h7r)=95155.21 [KW]

W reversible=m(h6r-h7r)-to*(s6r-s7r) m=105078.87

S1r=6.502 [kJ/kgK]

S6r=6.18432 [kJ/kgK]

S7r=6.3627 [kJ/kgK]

m=186.685

Exergía turbina de baja= wreal-wreversible=5980.33 [KW]

Wreal=m*(h7r-h1r)*(1-y)=83935.04KW

Wreversible=m(h7r-h1r)(1-y)-to*(s7r-s1r)(1-y)m=89915.37 [KW]

Primera ley CAA = myh7+(1-y) mh3– mh4 = 0

Y = 0.2283

Exergía turbina=15903.99

▪ La exergía destruida en el condensador (kW)= m(186.685kg/s)*(1-y)*To*(s1-s2)+Qcondensador/TL

=456332,91 kW.

temperatura de baja=456332.91KW

TL=298K

S2=0.7334

Qcondensador = m*(h1-h2)*(1-y)=1448.51 KJ/kg

m=1kg/s

▪ El costo del combustible por hora ($/h)

1 m³ de gas equivale a 0.657 kg

Para la industria el costo del m³ en el 2022 está más o menos en 4000 $ COP

Kg
∗1 m 3
h
∗4000 COP
0.657 Kg
30442.39 =185341796.04 COP / h
1 m3

▪ Relación aire-combustible Se utilizo una combustión completa con un exceso de aire del 17%

CH4+2*(O2+3.76N2) = CO2+2H2O+2*3.76*N2

Con el exceso de aire queda: CH4+2*1.17*(O2+3.76N2) = CO2+2H2O+8.7984 N2+0.34*O2

Formula relación aire combustible = (2.34*4.76*28.97) /16.043=20.1134 kg aire/kg CH4

▪ Determine las emisiones de CO2 por hora (t/h)

Sabemos que 1 kmol de ch4 produce un kmol de co2, y también vemos que 16.043 kg de ch4 producen
44 kg de co2, entonces simplemente sabiendo el flujo másico se hace lo siguiente

(8,98994 kg ch4 /s* 44 kgco2)/ 16.043 kgch4 = 24.656 kgco2/s*3600 s/h = 88761,86 kgco2/h

88761,86 kgco2/h* 1 ton/1000kg = 88, 761 ton co2/h

▪ Estime el costo del kWh neto producido para satisfacer los costos de combustible ($/kWh)

Si con 176000 kPa se produce un costo de 120842.82 $/s, se hace una regla de tres de la que se obtiene
que: 1kW = 0,686 $/s para luego: 0,686 $/ s * 3600 s/ h = 2469.6 $/h

▪ Determine el caudal volumétrico de aire atmosférico de entrada a la torre de enfriamiento (m3/h)

Qcondensador= m(186.685)*q
condensador(1448.51kj/kg)=270415.09kw

Primera ley en la torre de enfriamiento

Qcond=mH2o*(h3-h4)

H3=167.53kj/kg

H4=liquido sat=104.83 kj/kg

Despejando Mh2o=4312.84 kg/s

Ma en torre de enfriamiento= (Mh20*(H3-H4))/(H2-H1)+(W2-W1)*H4

Mreemplazo= Ma en torre(W2-W1)

h1= 45kJ/kg

h2= 127.5 kJ/kg

Estos valores de entalpia se sacaron con sus respectivos procedimientos, el primero se sacó usando la
carta psicrométrica:

Y la entalpia 2 al ser un valor muy alto en la carta psicométrica se utilizó una fórmula para hallarla, que
es:

H 2=CpT 2+W 2⋅ Hg 2
Siendo el Cp del aire: 1,005 kj/kg.K a aprox. 300K

Hg2 como la entalpía de vapor saturado de agua, sacada de tablas

Y W2, calculada de la fórmula:

 Pv
W 2=0,622⋅
PT −Pv
Y Pv siendo la presión de saturación del agua a 35°C de 5.6291 según la tabla anterior y PT la presión
atmosférica que es 101.325 kpa da un valor de:

kg vapor
W 2=0,036
kg aire
Utilizando el mismo proceso para W1, pero siendo la presión de vapor la presión de saturación a 25°C:

kg vapor
W 1=0,007882
kg aire
Reemplazando da que:

Ma de torre de enfriamiento es: 3164,68 kg/s = 11392848 kg/h

V °=Ma ⋅ V
3 3
kg m m
V °=11392848 ⋅0,844 =9615563,71
h kg h

▪ Determine el flujo másico de agua de reposición de la torre de enfriamiento (kg/h)

Mreemplazo= Ma en torre(W2-W1)

Mreemplazo= 3164,68 kg/s*(0,036-0,007882)

Mreemplazo= 88,984 kg/s