Escuela Superior Politécnica del Litoral

Facultad de Ingeniería en Electricidad y Computación

Distribución Eléctrica I – ELEG1032

PARALELO: 1
GRUPO: 2

TÍTULO: MODELAMIENTO DE UNA CENTRAL

TERMOELÉCTRICA DE VAPOR

INTEGRANTES: Duarte Campos, Franklin Fernando

Freire Plaza, Carlos André
Salgado Ochoa, Jocelyne Selena

FECHA DE ENTREGA: Guayaquil, 14 de diciembre de 2020

NOMBRE DEL PROFESOR: Ing. Johnny Rengifo

TÉRMINO ACADÉMICO: 2020-2S

ÍNDICE DE CONTENIDO
1. Objetivos ...........................................................................................................................3
1.1. Objetivo general.......................................................................................................3
1.2. Objetivos específicos .............................................................................................3
2. Problemática .....................................................................................................................3
4. Simulaciones ....................................................................................................................8
5. Análisis de resultados..................................................................................................12
6. Conclusiones..................................................................................................................15
7. Referencias .....................................................................................................................16

INDICE DE FIGURAS
Figura 1. Ciclo Rankine modelado para el problema planteado. .................................... 4
Figura 2. Ciclo Rankin del problema planteado en el software CYCLEPAD ™️. ............ 8
Figura 3. Resultados de la zona 1: Heater1-Turbina1..................................................... 8
Figura 4. Resultados de la zona 2: Turbina1- Heater2.................................................... 9
Figura 5. Resultados de la zona 3: Heater2 -Turbina2.................................................... 9
Figura 6. Resultados de la zona 4: Turbina2 - Cooler ................................................... 10
Figura 7. Resultados de la zona 5: Cooler - Pump ........................................................ 10
Figura 8. Resultados de la zona 6: Pump – Heater1..................................................... 11
Figura 9.Resultados del sistema. ................................................................................... 11
Figura 10. Presión de la zona 4 y 5 donde se encuentra el condensador. ................... 12
Figura 11. Potencia neta de salida del sistema. ............................................................ 13
Figura 12. Eficiencia térmica y BWR del sistema. ......................................................... 13
Figura 13. Reducción del flujo másico en el sistema y el efecto en la potencia de salida
de un 50% menos........................................................................................................... 14
1. Objetivos
1.1. Objetivo general
Simular el comportamiento de una central termoeléctrica de vapor
usando el software CYCLEPAD™️.
1.2. Objetivos específicos
• Simular el comportamiento de la termoeléctrica de vapor
usando un ciclo de Rankine mejorado dado que se toma en
cuenta las eficiencias isentrópicas de la turbina y de la
bomba.
• Aplicar conocimientos adquiridos en las clases teóricas para
el manejo del software CYCLEPAD™️ y así corroborar que
los resultados obtenidos analíticamente sean los correctos.
• Investigar el método utilizado en las centrales
termoeléctricas de vapor para regular la potencia entregada
de la turbina.
2. Problemática
Una planta de energía de vapor opera en el ciclo de recalentamiento
Rankine. El vapor entra en la turbina de alta presión a 12.5 MPa y 550 °C
con un caudal másico de 7.7 kg/s y sale al recalentador a 2 MPa. El vapor
se recalienta hasta 450°C antes de que se expanda en la turbina de baja
presión. Las eficiencias isentrópicas de la turbina y de la bomba son de
85 y 90 por ciento, respectivamente. El vapor sale del condensador como
un líquido saturado. Si el contenido de humedad del vapor a la salida de
la turbina no debe superar el 5 por ciento, utilizando CYCLEPAD
determine:
A) la presión del condensador
B) la potencia neta de salida
C) la eficiencia térmica, el BWR y el heat-rate en BTU/kWh.
D) Investigue el método utilizado para regular la potencia entregada una
turbina de vapor en una central eléctrica simule la operación con un 50%
de potencia respecto al literal (b). Para este punto determine la eficiencia
térmica.
 3. Resolución analítica de la problemática

4.
5.
6.
7.
8.
9.
10.
11.
12.

Figura 1. Ciclo Rankine modelado para el problema planteado.

Datos
𝐾𝑔⁄
𝑃3 = 12.5 𝑀𝑃𝑎, 𝑇3 = 550 ℃, 𝑚̇ = 7.7 𝑠
𝑃4 = 2 𝑀𝑃𝑎
𝑇5 = 450 ℃
𝑥6 = 95%
𝐸𝑙 𝑣𝑎𝑝𝑜𝑟 𝑠𝑎𝑙𝑒 𝑑𝑒𝑙 𝑐𝑜𝑛𝑑𝑒𝑛𝑠𝑎𝑑𝑜𝑟 𝑐𝑜𝑚𝑜 𝑙𝑖𝑞𝑢𝑖𝑑𝑜 𝑠𝑎𝑡𝑢𝑟𝑎𝑑𝑜
Desarrollo
La entalpía y entropía en el punto 3 se la obtiene de la tabla de vapor
sobrecalentado, dado que no hay una tabla con una presión de 12.5 MPa se
procede a utilizar la herramienta para interpolar que ofrece CYCLEPAD.
𝐾𝐽
ℎ3 = 3475.05 ⁄𝑘𝑔

𝐾𝐽
𝑠3 = 6.6289 ⁄𝑘𝑔 ∗ 𝐾

Dado que la turbina es un proceso isentrópico tenemos: 𝑠3 = 𝑠4

𝐾𝐽
𝑠4 = 6.6289 ⁄𝐾𝑔 ∗ 𝐾

Para obtener la entalpia en el punto 4 se utilizará la tabla de vapor

sobrecalentado a una presión de 2 MPa y nos ayudaremos con el dato de
entropía hallado previamente. Dado que en la tabla no está el valor puntual
nuevamente usaremos la herramienta de interpolación de CYCLEPAD, con eso
obtenemos:
𝐾𝐽
ℎ4 = 2946.89 ⁄𝑘𝑔

Dado que la caldera es un proceso isobárico la presión en punto 4 y en el punto

5 son iguales, junto a la temperatura que nos dan como dato ya tenemos las dos
propiedades termodinámicas necesarias para buscar en la tabla de vapor
sobrecalentado la entalpia y la entropía. Al recurrir a la tabla nos percatamos que
hay que interpolar nuevamente ya que a una presión de 2 MPa no están los
valores de entropía y entalpia a 450 °C, usando CYCLEPAD se tiene:
𝐾𝐽
ℎ5 = 3357.5 ⁄𝑘𝑔

𝐾𝐽
𝑠5 = 7.2845 ⁄𝑘𝑔 ∗ 𝐾

La entropía en el punto 6 será la misma que la de entrada de la turbina dado a

que es un proceso isentrópico. Además, conocemos que tenemos una mezcla
de líquido vapor en el punto 6, ya que el dato es la calidad de dicha mezcla.
𝐾𝐽
𝑠6 = 7.2845 ⁄𝑘𝑔 ∗ 𝐾

𝑥6 = 95%
En este punto debemos de trabajar con la tabla de agua saturada. A falta de una
presión con la cual trabajar se procede a identificar los datos de entropía (𝑠𝑔 𝑦 𝑠𝑓 )
con los que se obtiene una calidad del 95%.
𝑠6 − 𝑠𝑓 7.2845 − 1.0753
𝑥6 = = = 0.9498%
𝑠𝑔 − 𝑠𝑓 7.6122 − 1.0753

Con esto nos podemos dar cuenta que los datos con los que se obtiene la calidad
más próxima es a 80 °C y una presión de 0.04739 MPa. La entalpia en dicho
punto se la calcula de la siguiente manera:
ℎ6 = 𝑥6 ∗ (ℎ𝑔 − ℎ𝑓 ) + ℎ𝑓 = 0.95 ∗ (2643.7 − 334.91) + 334.91
𝐾𝐽
ℎ6 = 2528.26 ⁄𝑘𝑔

A la salida del condensador tendremos la misma presión que a la entrada dado

que es un proceso isobárico, por ende 𝑃6 = 𝑃1
𝑃1 = 0.04739 𝑀𝑃𝑎 = 47.39 𝐾𝑃𝑎
Para el punto 1 buscamos en la tabla de agua saturada la presión obtenida y
dado que se conoce la fase (liquido saturado) podemos determinar la entalpia
directamente:
𝐾𝐽
ℎ1 = 334.91 ⁄𝑘𝑔
Para el punto 2, el cual es la salida de la bomba se calcula la entalpia utilizando
el trabajo de compresión. Se debe de tener en cuenta que la presión en dicho
punto será igual a la que se tiene en punto 3 ya que la caldera es un proceso
isobárico.
𝑊̇𝐵
= ℎ2 − ℎ1 ≈ 𝑣1 ∗ (𝑃2 − 𝑃1 )
𝑚̇
ℎ2 ≈ 𝑣1 ∗ (𝑃2 − 𝑃1 ) + ℎ1 ≈ 1.0291 ∗ (12.5 − 0.04739) + 334.91
𝐾𝐽
ℎ2 = 347.72 ⁄𝑘𝑔

Para la obtención de la potencia neta de salida, la eficiencia térmica el BWR y el

Heat Rate debemos de tener en cuenta la eficiencia isentrópica que nos dan
como datos.
𝜂𝑇 = 85%
𝜂𝐵 = 90%
Recalculando la entalpia a la salida de la turbina utilizando la eficiencia
isentrópica se tiene lo siguiente:
ℎ3 − ℎ4
𝜂𝑇 =
ℎ3 − ℎ4𝑠
ℎ4 = ℎ3 − 𝜂𝑇 ∗ (ℎ3 − ℎ4𝑠 ) = 3475.05 − 0.85 ∗ (3475.05 − 2946.89)
𝐾𝐽
ℎ4 = 3026.114 ⁄𝑘𝑔

Recalculando la entalpia a la salida de la bomba utilizando la eficiencia

isentrópica se tiene lo siguiente:
ℎ2𝑠 − ℎ1
𝜂𝐵 =
ℎ2 − ℎ1
ℎ2𝑠 − ℎ1 347.72 − 334.91
ℎ2 = + ℎ1 = + 334.91
𝜂𝐵 0.90
𝐾𝐽
ℎ2 = 349.14 ⁄𝑘𝑔

Para calcular el trabajo neto se debe de tener en cuenta las entalpias

recalculadas:
̇
𝑊𝑛𝑒𝑡𝑜 ̇
𝑊𝑇𝑢𝑟𝑏𝑖𝑛𝑎𝑠 𝑊̇𝐵 ̇
𝑊𝑇1 ̇
𝑊𝑇2 𝑊̇𝐵
= − =( + )−
𝑚̇ 𝑚̇ 𝑚̇ 𝑚̇ 𝑚̇ 𝑚̇
̇
𝑊𝑛𝑒𝑡𝑜
= ((ℎ3 − ℎ4 ) + (ℎ5 − ℎ6 )) − (ℎ2 − ℎ1 )
𝑚̇
 ̇
𝑊𝑛𝑒𝑡𝑜
= ((3475.05 − 3026.114) + (3357.5 − 2528.26)) − (349.14 − 334.91)
𝑚̇
̇
𝑊𝑛𝑒𝑡𝑜 𝐾𝐽
= 1263.946 ⁄𝑘𝑔
𝑚̇
Una vez calculado el trabajo neto procedemos a calcular la potencia generada
̇
𝑊𝑛𝑒𝑡𝑜
𝑃𝑔 = ∗ 𝑚̇ = 1263.946 ∗ 7.7 = 9732.38 𝐾𝑊
𝑚̇
La eficiencia térmica, el heat rate y el BWR del ciclo la obtenemos de la siguiente
manera:
Eficiencia
̇
𝑊𝑛𝑒𝑡𝑜 ̇
𝑊𝑛𝑒𝑡𝑜
1263.946
𝜂 = 𝑚̇ = 𝑚̇ =
𝑄𝑖𝑛 ̇ (ℎ3 − ℎ2 ) + (ℎ5 − ℎ4 ) (3475.05 − 349.14) + (3357.5 − 3026.114)
𝑚̇
𝜂 = 0.3656 ∗ 100 = 36.56%
Heat rate
̇
𝑄𝑖𝑛 𝐾𝐽
1 𝑚 ̇ ∗ 𝑚̇ 3457.296 ∗ 7.7 26621.18 [ ⁄𝑠]
ℎ𝑒𝑎𝑡 𝑟𝑎𝑡𝑒 = = = =
𝜂 𝑊𝑛𝑒𝑡𝑜̇ 1263.946 ∗ 7.7 𝐾𝐽
∗ 𝑚 ̇ 9732.38 [ ⁄𝑠]
𝑚̇
𝐾𝐽 0.9478 𝐵𝑇𝑈 3600 𝑠
26621.18 [ ⁄𝑠 ] ∗ ∗ 1ℎ 90833595.85 [𝐵𝑇𝑈⁄ℎ]
1 𝐾𝐽
ℎ𝑒𝑎𝑡 𝑟𝑎𝑡𝑒 = =
9732.38[𝐾𝑊 ] 9732.38 [𝐾𝑊 ]

ℎ𝑒𝑎𝑡 𝑟𝑎𝑡𝑒 = 9333.13 [𝐵𝑇𝑈⁄𝐾𝑊ℎ]

BWR
𝑊̇𝐵
𝑚̇ (349.14 − 334.91)
𝐵𝑊𝑅 = =
𝑊̇ 𝑊̇ ((3475.05 − 3026.114) + (3357.5 − 2528.26))
( 𝑇1 + 𝑇2 )
𝑚̇ 𝑚̇
𝐵𝑊𝑅 = 1.113 ∗ 10−2 ∗ 100 = 1.113%
4. Simulaciones

Figura 2. Ciclo Rankin del problema planteado en el software CYCLEPAD ™️.

Figura 3. Resultados de la zona 1: Heater1-Turbina1.

Figura 4. Resultados de la zona 2: Turbina1- Heater2.

Figura 5. Resultados de la zona 3: Heater2 -Turbina2

Figura 6. Resultados de la zona 4: Turbina2 - Cooler

Figura 7. Resultados de la zona 5: Cooler - Pump

Figura 8. Resultados de la zona 6: Pump – Heater1

Figura 9.Resultados del sistema.

5. Análisis de resultados
Para el análisis del sistema se decidió llevarlo a un modelo combinado
donde se modelo la turbina de alta presión con su respectiva eficiencia
isentrópica y asimismo la bomba mientras que para la turbina de baja
presión se la considero ideal para aprovechar la propiedad isentrópica de
la turbina y tener la entropía en la entrada del condensador, de lo contrario
no se contaba con las respectivas 2 propiedades termodinámicas
necesarias para realizar el modelamiento real. En el literal A se pide
obtener la presión del condensador dado que este equipo en el ciclo
Rankin tiene como característica un comportamiento isobárico, la presión
que tiene en la zona 4 se mantiene en la zona 4 con un valor de 47.42
[kPa] como se puede observar en la figura 9.

Figura 10. Presión de la zona 4 y 5 donde se encuentra el condensador.

En el literal B se pide obtener la potencia de salida neta del sistema la

cuál tiene un valor de 9747 [kW] como se observa en la figura 10 cabe
recalcar que la potencia de salida neta es la resta de la potencia de salida
con la potencia que consume la bomba que es de 110 [W].
 Figura 11. Potencia neta de salida del sistema.

En el literal C se pide obtener la eficiencia térmica, el BWR y el Heat - rate

los cuales se los puede observar en la figura 11 donde la eficiencia térmica
es de 36.58% que es un valor acercado a la realidad de la eficiencia
térmica que suelen tener este tipo de centrales, el BWR el cual es la
relación del trabajo de entrada de la bomba con respecto al trabajo que
entrega la turbina y tiene un valor de 1.11% y el Heat – rate se lo calcula
por medio del inverso de la eficiencia térmica y tiene un valor de 9333.13
[BTU/kWh].

Figura 12. Eficiencia térmica y BWR del

sistema.
En el literal D para regular la potencia de la turbina de vapor se tomó la decisión
de realizarlo por medio de un bloqueo del flujo de vapor, es decir realizar una
reducción de la capacidad de paso del flujo de la turbina de vapor y el campo en
el perfil de vapor efectivo en las palas de la turbina de vapor. Realizar estos
cambios se obtienen una reducción del flujo de vapor, una potencia reducida y
una degradación de la eficiencia de la turbina. Ejemplos comunes son los
depósitos de cobre y aluminio en turbinas de vapor de alta presión y los depósitos
de sílice en turbinas de presión media y baja.
̇
𝑊𝑛𝑒𝑡𝑜
= ((ℎ3 − ℎ4 ) + (ℎ5 − ℎ6 )) − (ℎ2 − ℎ1 )
𝑚̇
Entonces por unidades se tiene que el trabajo neto másico sobre el flujo másico
𝐾𝐽 𝐾𝑔
es [𝐾𝑔]y si se realiza la multiplicación del flujo másico reducido a la mitad [ ]
𝑠

𝐾𝐽 1 𝐾𝑔 1
Se obtiene [ ]∗ [ ] = ∗ [𝑘𝑊]
𝐾𝑔 2 𝑠 2

Figura 13. Reducción del flujo másico en el sistema y el efecto en la potencia de salida de un 50% menos.
6. Conclusiones
• Se pudo concluir que los datos encontrados analíticamente son
idénticos a los valores arrojados que da el software CICLEPAD, a
pesar de que no se consideró las dos turbinas reales.
• El sistema modelado nos arroja una eficiencia de 36.58% lo cual
es un dato válido dado que se encuentra en el rango de 30-40%
que normalmente se observa en este tipo de centrales
termoeléctricas a pesar de que el modelo utilizado consta de una
turbina considerada ideal.
• El seguimiento del heat rate en este tipo de centrales es primordial
ya que el ajuste inmediato sobre esta variable ayuda a controlar
cualquier anomalía que se presenten al operar la maquina fuera de
sus parámetros de operación, además mejorar el heat rate implica
mejorar la eficiencia de la planta ya que ambas variables se
encuentran ligadas.
• El heat rate es una variable importante ya que es un indicativo de
la cantidad de combustible a utlizarse para generar un KWh, tal y
como se aprecia previamente se obtuvo de valor de 9333.13
[BTU⁄KWh], por lo que se puede concluir que el precio para proveer
de combustible a la central analizada es bastante alto
• Para regular la potencia de salida de una turbina existen varios
métodos, algunos usan válvulas para modificando el delta de
presión y así conseguir reducir la potencia de salida, pero vale
recalcar que el grado de reducción es máximo de 1000 [kW] por
válvula instalada en la entrada de la turbina, se concluyó que el
mejor método para regular la potencia y reducirla al 50% era
realizar una reducción del flujo másico.
7. Referencias

[1] A. Almasi, «Processing Magazine,» 29 Julio 2016. [En línea]. Available:

https://www.processingmagazine.com/home/article/15586711/steam-turbine-flow-
operation. [Último acceso: 14 Diciembre 2020].

[2] A. Fuentes, «Predictiva21,» [En línea]. Available: https://predictiva21.com/gerencia-

combustible-centrales-termicas-vapor/. [Último acceso: 14 Diciembre 2020].

[3] M. Solorzano, M. Lopez y R. Soto, «Slideshare,» 26 Noviembre 2013. [En línea]. Available:
https://es.slideshare.net/MonicaSolorzano/turbinas-de-vapor-presnentacion-final.
[Último acceso: 14 Diciembre 2020].