TRABAJO FINAL MAQUINAS TÉRMICAS

DISEÑO DE UNA CENTRAL DE GENERACIÓN TÉRMICA CON CICLO RANKINE

SOFÍA MARIANELLA ARCHILA PICO

CAMILA ANDREA MORENO MORA

DIANA DEL PILAR CRUZ JIMENEZ

WILMER ALEJANDRO MORALES GÓMEZ

INGENIERÍA MECÁNICA
FACULTAD TECNOLÓGICA
UNIVERSIDAD DISTRITAL FRANCISCO JOSÉ DE CALDAS
BOGOTÁ
2020
 PLANTEAMIENTO DEL PROBLEMA

Se requiere hacer los cálculos básicos de una central de generación térmica con ciclo Rankine con
recalentamiento intermedio y tres extracciones intermedias hacia tres calentadores dos abiertos y
uno cerrados, completando de esta forma un ciclo regenerativo.
La presión del condensador es de 0,1 bar, y la velocidad de rotación de la turbina es 1.800 rpm. La
turbina cuenta con cuatro tipos de etapas, las que se detallan en la Tabla 1.
Tabla 1. Tipos de etapas en la turbina

Et Tipo Dm (mm) No Alabes Admisión Cant

1 Laval 960 110 Parcial 1
2 Curtis 2R 960 110 Parcial 1
3 Rateau 920 106 Total N1
4 Reacción 940 106 Total N2

La temperatura y presión de salida de la caldera, o de entrada a la turbina, (T1 y P1), corresponde

a los datos suministrados en la Tabla 2; en igual forma se indica la potencia a generar, los tipos de
calentadores, (CA: Calentador Abierto; CC: Calentador Cerrado) para el ciclo regenerativo, y las
presiones a las que se hace cada una de las tres (3) extracciones (Pext1, Pext2 y Pext3),
conformando cuatro (4) secciones (A, B, C y D). También se indica el tipo de etapa correspondiente
a la entrada de las secciones A y B de la turbina.
Tabla 2. Datos complementarios
1 2 3 4 5 6 7 8 9 10 11 12
Eq T1 P1 Pot 1ra 2da 3ra Entr Entr Pext1 Pext2 Pext3
°C (bar) (MW) extr extr extr Sec A Sec B (bar) (bar) (bar)
2 450 65 36 CC1 CA CC2 Curtis Laval 20 5 1,5

Se asume que la central requiere 2% de la potencia nominal para consumo interno. La potencia de
salida debe ser neta y corresponder a la expresada en la cuarta columna de la Tabla 2.
Al momento de hacer la primera extracción, para el ciclo regenerativo, se realiza una extracción
total del vapor para un recalentamiento intermedio que eleva la temperatura hasta la misma
temperatura expresada en la segunda columna de la Tabla 2.
Las cuatro (4) secciones de la turbina estarán conformadas por las etapas indicadas en la Tabla 3.
Tabla 3. Tipos de etapas en cada una de las cuatro (4) secciones de la turbina

Sección A B C D
La indicada en la La indicada en la 50% Etapas Rateau Todas las etapas
Etapas columna 8 de la columna 9 de la + 50% Etapas de son de Reacción
Tabla 2 + Etapas Tabla 2 + Etapas Reacción
Rateau Rateau
Se asume que cada sección de la turbina tiene un rendimiento isoentrópico de 90%.
 Tabla 4.1 Datos necesarios para cálculo de etapas (parte 1)
Etapa Admisión αi Ρopt Βe Km Kf
LAVAL Parcial (15-20)° Cos αi/2 Βi – (0-2)° (0,85-0,9) -
CURTIS (2R) Parcial (15-20)° Cos αi1/4 Βi – (0-2)° (0,85-0,9) (0,85-0,9)
RATEAU Total (18-22)° (0,56-0,60) Βi – (2-5)° (0,85-0,9) -
REACCIÓN Total (20-25)° Cos αi/(2)1/2 αi - -

Tabla 4.2 Datos necesarios para cálculo de etapas (parte 2)

Etapa ηn ηb ηgs Δh´gs
LAVAL 0,93 – 0,97 2Cb*CΔw/Cai 2
ηb* ηn Cb*CΔw
2
CURTIS (2R) 0,93 – 0,97 2Cb*∑CΔw/Cai ηb1y2* ηn Cb*∑CΔw
RATEAU 0,91 – 0,95 2Cb*CΔw/Cai2 ηb* ηn Cb*CΔw
REACCIÓN 0,90 – 0,94 2Cb*CΔw/(Cai2+Cre2-Cri2) ηb* ηn Cb*CΔw
 DESARROLLO ESQUEMÁTICO

Imagen 1. Diagrama de ubicación de equipos de la central

 J
I
G H

E F

P
R
Q

N T

BL

Imagen 2. Diagrama T-S

 J
H I
G

F
E

P
R

O BO

BL

Imagen 2. Diagrama H-S

 MEMORIA DE CÁLCULO

1. Cálculos Central térmica.

1.1 Calculo de los estados termodinámicos.

Usando los datos de la tabla 2 y la aplicación “Steam Tables Free” se calculan los estados
termodinámicos para cada componente, dando como resultado los siguientes datos:

Tabla 5. Estados termodinámicos de la central de generación térmica

h-
ESTADO Presión Temp. s-Entalpía Volumen
Entalpía Calidad
S (bar) (°C) (kJ/kg) específico (m3/kg)
(kJ/kg)
E 65 450,00 3296,30 6,68 0,04758 Sobrecalentado
F 20 291,96 3006,13 6,74 0,12274 Sobrecalentado
G 20 450,00 3357,83 7,29 0,16264 sobrecalentado
H 5 272,23 3007,29 7,36 0,49357 Sobrecalentado
I 1,5 148,40 2769,23 7,41 1,27387 Sobrecalentado
J 0,1 45,83 2388,77 7,54 13,40793 0,918088
K 0,1 45,83 191,83 0,65 0,00101 0
BL 5 45,85 192,33 0,65 0,00101 Subenfriado
M 5 111,37 467,38 1,43 0,00105 Subenfriado
N 5 151,84 640,12 1,86 0,00109 0
BO 65 152,51 646,66 1,86 0,00108 Subenfriado
P 65 280,82 1241,14 3,08 0,00133 0
Q 20 212,37 908,59 2,45 0,00117 0
R 5 151,84 908,59 2,49 0,04845 0,127394
S 1,5 111,37 467,13 1,43 0,00105 0
T 0,1 45,83 467,13 1,51 1,68100 0,115044

1.2 Extracciones de masa

Se asumió un flujo másico total de entrada de 1 kg para calcular las extracciones de masa de la
turbina y se realiza un balance de masas en cada uno de los calentadores:

 Calentador cerrado 1:

Y 1∗hF

1∗hP

1∗hBO

Y 1∗hQ
 Y 1∗hF + 1∗h BO=Y 1∗hQ +1∗h P

Y 1∗hF −Y 1∗hQ=1∗h P−1∗hBO

Y 1∗( hF −h Q )=hP −h BO

h P−hBO
Y 1=
( h F −hQ )
1241,14−646,66
Y 1= =0,2834
( 3006,13−908,59 )
 Calentador abierto:
Y 2∗hH

1∗hN

( 1−Y 1−Y 2 )∗hM

 Y 1∗hR

( 1−Y 1−Y 2 )∗hM +Y 2∗h H +Y 1∗h R=1∗h N

1∗hM −Y 1∗h M −Y 2∗h M + Y 2∗h H +Y 1∗h R =1∗h N

−Y 2∗hM +Y 2∗h H =1∗hN −1∗hM +Y 1∗h M −Y 1∗h R

Y 2∗(−h M +h H )=1∗h N −1∗h M +Y 1∗hM −Y 1∗h R

h N −h M +Y 1∗h M −Y 1∗hR
Y 2=
−h M +h H

640,12−467,38+0,2834∗467,38−0,2834∗908,59
Y 2= =0,0187
−467,38+3007,29

 Calentador cerrado :
Y 3∗hI
 ( 1−Y 1−Y 2 )∗hM

( 1−Y 1−Y 2 )∗hBL

Y 3∗hS

Y 3∗hI + ( 1−Y 1−Y 2 )∗h BL=( 1−Y 1−Y 2 )∗h M +Y 3∗h S

Y 3∗hI −Y 3∗hS =( 1−Y 1−Y 2 )∗h M −( 1−Y 1−Y 2 )∗hBL

Y 3∗(h I −hS )=( 1−Y 1−Y 2 )∗hM −( 1−Y 1 −Y 2 )∗h BL

( 1−Y 1−Y 2 )∗hM −( 1−Y 1−Y 2)∗h BL

Y 3=
h I −h S

( 1−0,2834−0,0187 )∗467,38−( 1−0,2834−0,0187 )∗192,33

Y 3= =0,0833
2769,23−467,13

1.3 Flujo másico

1.4 Calidad de salida del vapor

La calidad del vapor a la salida de la turbina se puede observar en el estado J de la central de

generación térmica, siendo x=0,918088 = 91,80%.

2. Cálculos de la turbina

2.1 Sección A etapas Curtis y Reateau

2.2.1 Primer grupo de alabes móviles etapa Curtis

Según los parámetros definidos en el planteamiento del problema se obtiene:

 Se calcula la velocidad del alabe como:

Cb=W angular∗r curtis

188,49∗0,96 m
Cb= =90,47
2 s
 Se calcula Popt:
 cos (α i 1)
Popt =
4
cos (15)
Popt = =0,241
4
 El valor de la velocidad absoluta de los alabes de entrada es:

Cb
C ai =
P opt
90,47 m
C ai = =374,67
0,241 s

 Calculando la velocidad relativa de entrada a los alabes se obtiene:

Cri=√ cb 2+Cai 2−2∗Cai∗Cb∗cos(αi)

m
Cri=√ 90,47 +374,67 −2∗374,67∗90,47∗cos ( 15 ) =288,23
2 2
s

 Velocidad de cruce axial:

Cfi=Cai∗Sen(αi)
m
Cfi=374,67∗Sen ( 15 )=96,97
s

 Se debe calcular el valor de βi:

Cai∗Sen( αi)
βi=
Cri
374,67∗Sen(15)
βi= =19,66 °
288,23
 Se obtiene el valor de Cre:

Cre=Cri∗K m

m
Cre=288,23∗0,85=245,001
s
 Se calcula el βe:

βe=βi−1
βe=19,66 °−1=18,66 °

 La velocidad absoluta de salida de los alabes es:

 Cae= √Cb 2 +Cre 2−2∗Cre∗Cb∗cos(βe)
m
Cae= √ 90,47 + 245,001 −2∗245,001∗90,47∗cos( 18,66)=161,88
2 2
s
 Se continua con la velocidad de cruce axial de salida:

Cfe=Cae∗Sen( βe)
m
Cfe=161,88∗Sen ( 18,66 ) =51,79
s

 Se calcula la velocidad de volteo de los alabes móviles:

Cw=Cre∗cos ( βe )+Cri∗cos( βi)

m
Cw=245,001∗cos ( 18,66 ) +288,23∗cos ( 19,66 ) =503,55
s
2.2.2 Segundo grupo de alabes móviles etapa Curtis

 Como Cb y Popt son iguales en ambos grupos de alabes móviles se tiene:

m
Cb 2=90,47
s
Popt 2=0,241

 El valor de la velocidad absoluta de los alabes de entrada es:

C ai2 =Cae 1∗K f

m
C ai2 =161,88∗0,85=137,60
s

 Calculando la velocidad relativa de entrada a los alabes se obtiene:

Cri 2=√ Cb 22+Cai 22−2∗Cai 2∗Cb 2∗cos(αi2)

m
Cri 2=√ 90,472 +137,602−2∗137,60∗90,47∗cos ( 22 )=63,51
s

 Velocidad de cruce axial:

Cfi 2=Cai2∗Sen (αi2 )

m
Cfi 2=137,60∗Sen ( 22 )=51,54
s

 Se debe calcular el valor de βi:

 Cai2∗Sen( αi2)
βi 2=
Cri 2
137,60∗Sen(22)
βi 2= =54,25°
63,51
 Se obtiene el valor de Cre:

Cre 2=Cri 2∗K m

m
Cre 2=63,51∗0,85=53,98
s
 Se calcula el βe:

βe 2=βi 2−1
βe 2=54,25° −1=53,25 °

 La velocidad absoluta de salida de los alabes es:

Cae2= √Cb22 +Cre 22−2∗Cre 2∗Cb2∗cos ( βe2 )

m
Cae2= √90,47 2+53,98 2−2∗53,98∗90,47∗cos (53,25)=72,497
s
 Se continua con la velocidad de cruce axial de salida:

Cfe 2=Cae2∗Sen(βe 2)
m
Cfe 2=72,497∗Sen ( 53,25 )=58,08
s

 Se calcula la velocidad de volteo de los alabes móviles:

Cw 2=Cre 2∗cos ( βe2 ) +Cri 2∗cos (βi 2)

m
Cw=53,98∗cos ( 53,25 ) +63,51∗cos ( 54,25 )=69,408
s
 Se asume la eficiencia de las toberas como:

ղ n=0,93

 Calculando la eficiencia de los alabes móviles se obtiene:

2∗Cb∗∑ Cw
ղ b= 2
Cai
 2∗90,47∗(503,55+69,408)
ղ b= 2
=0,739
374,67

 Por lo tanto se obtiene que la eficiencia de la etapa es:

ղ gs =ղ b∗ղ n
ղ gs =0,739∗0,93=0,687