Taller Final Motores Térmicos. Cantillo - Luquez
Taller Final Motores Térmicos. Cantillo - Luquez
Taller Final Motores Térmicos. Cantillo - Luquez
Presentado por: Brayan David Cantillo Ortega, Kevin Andrés Luquez Vizcaíno.
Cédulas: 1001893257, 1001869827.
Presentado a: MSc. Franklin Consuegra Vargas. 2023-1. Universidad del Atl.
Desarrollo:
1. Considérese en primer lugar un esquema de un ciclo Brayton regenerativo:
𝑇1 = 420 𝐾
𝑃1 = 130 𝑘𝑃𝑎
ℎ1 = 421.26 𝑘𝐽/𝑘𝑔
Estado 2:
𝑃2 = 1110 𝑘𝑃𝑎
𝑇2 = 853 𝐾
Se interpoló entre 840 K y 860 K para determinar el valor de la entalpía en el estado 2
Interpolación h2
T (K) h (kJ/kg)
840 866,08
853 880,5035
860 888,27
ℎ2 = 880.5035 𝑘𝐽/𝑘𝑔
Estado 3:
𝑇3 = 1200 𝐾
ℎ3 = 1277.79 𝑘𝐽/𝑘𝑔
Estado 4: Como se desconocen las propiedades a la salida de la turbina, puede asumirse que en ella
se da una expansión isentrópica y calcular 𝑇4𝑠 como una función de 𝑃3 , 𝑃4 , 𝑇3 y 𝑘 de la siguiente
forma:
𝑇3 𝑃3 (𝑘−1)/𝑘 𝑇3
( )=( ) → 𝑇4𝑠 =
𝑇4𝑠 𝑃4 𝑃 (𝑘−1)/𝑘
( 3)
𝑃4
Pero 𝑃4 = 𝑃1 y 𝑃3 = 𝑃2 . Por lo tanto:
1200 𝐾
𝑇4𝑠 =
1110 𝑘𝑃𝑎 (1.325−1)/1.325
( )
130 𝑘𝑃𝑎
𝑇4𝑠 = 709.136 𝐾
Se interpoló entre 700 K y 710 K para determinar la entalpía isentrópica o ideal en el estado 4
Interpolación h4
T (K) h (kJ/kg)
700 713,27
709,136 723,0914
710 724,02
A partir de la temperatura obtenida, se interpola entre 780 K y 800K para determinar la entalpía real
en 4:
Interpolación h4
T (K) h (kJ/kg)
780 800,03
792,583 813,821
800 821,95
ℎ4 = 813.821 𝑘𝐽/𝑘𝑔
𝑘𝐽 𝑘𝐽
𝑞𝑟𝑒𝑔𝑒𝑛,𝑟𝑒𝑎𝑙 = 830.492 − 880.5035
𝑘𝑔 𝑘𝑔
Tal como puede observarse, el calor en el regenerador es negativo. Esto contradice las condiciones
básicas de funcionamiento de un regenerador, conduciendo a un resultado contraproducente. Esto
puede deberse a las condiciones de entrada y salida del compresor. Obtener un calor negativo, implica
que se requiere un mayor aporte energético en la cámara de combustión, disminuyendo de esta forma
de una manera considerable la eficiencia del ciclo.
𝜂 = 1.1 %
2. En este inciso, se debía realizar un análisis paramétrico variando tanto relaciones de presión
como relaciones de temperatura y manteniendo las presiones de entrada y salida constantes
en un −10% , −20%, +10% , 20%. Este proceso permitió determinar cuáles eran las
condiciones de operación más óptimas en el compresor para un mejor resultado en la
eficiencia del ciclo. Para ello, se realizó una matriz de Excel donde se determinaron cada uno
de los cálculos. Los resultados, se encuentran a continuación.
P2/P1 variable
% Variación qregen,real (kJ/kg) ηtér
-20 29,989 23,5
-10 29,978 23,5
0 -50,012 1,1
10 -6,960 13,5
20 -23,772 8,7
T2/T1 variable
% Variación qregen,real (kJ/kg) ηtér
-20 -4,868 -13,8
-10 -58,394 3,5
0 -50,012 1,1
10 -136,129 22,4
20 -165,253 28,2
Tal como puede observarse variaciones como las del -20% en la relación de temperatura conducen a
valores inconsistentes. Sin embargo, algunas como las de -20% y -10% en la relación de presión o
+10% en la relación de temperatura conducen a resultados más favorables comparados con los del
inciso anterior. Por lo tanto, nos lleva a concluir que es de vital importancia contralar las variables
de entrada y de salida en el compresor, puesto que son de alta relevancia en la eficiencia y en el
rendimiento del regenerador. El Excel puede consultarse en el siguiente hipervínculo: Cálculo Taller
final motores - Cantillo, Luquez
3. Considérese en primer lugar un esquema del proceso:
Variable Valor
P8 (MPa) 12
P8 (kPa) 12000
P10 (MPa) 2
P10 (kPa) 2000
P11 (MPa) 0,7
P11 (kPa) 700
P12 (MPa) 0,4
P12 (kPa) 400
Wneto (MW) 11
Wneto (kW) 11000
ⴖT 0,85
ⴖP 0,95
Determinación de propiedades:
Estado 8: Vapor sobrecalentado
𝑃8 = 12 𝑀𝑃𝑎 = 12000 𝐾𝑝𝑎
𝑇8 = 500°
𝒌𝑱
𝒉𝟖 = 𝟑𝟑𝟒𝟗, 𝟗
𝒌𝒈
𝒌𝑱
𝒔𝟖 = 𝟔, 𝟒𝟗𝟏𝟗𝟖
𝒌𝒈. 𝑲
𝒌𝑱
𝒉𝟗 = 𝟐𝟗𝟒𝟔, 𝟏𝟕𝟒
𝒌𝒈
Ahora calculamos la entropía real
𝒌𝑱
𝑺𝟗 = 𝟔, 𝟔𝟐𝟔
𝒌𝒈. 𝑲
Estado 10: Vapor sobrecalentado
𝑃10 = 2 𝑀𝑃𝑎 = 2000 𝐾𝑝𝑎
𝑇10 = 500°
𝒌𝑱
𝒉𝟏𝟎 = 𝟑𝟒𝟔𝟖, 𝟑
𝒌𝒈
𝒌𝑱
𝒔𝟏𝟎 = 𝟕, 𝟒𝟑𝟑𝟕
𝒌𝒈. 𝑲
𝒌𝑱
𝒉𝟏𝟏 = 𝟑𝟏𝟖𝟕, 𝟎𝟑𝟓𝟔
𝒌𝒈
Ahora calculamos la entropía real, esta se calculó haciendo una doble interpolación ya que no había
tabla de 0,7Mpa.
𝒌𝑱
𝑺𝟏𝟏 = 𝟕, 𝟓𝟏𝟒𝟖
𝒌𝒈. 𝑲
𝒌𝑱
𝒉𝟏𝟐 = 𝟑𝟎𝟔𝟎, 𝟒𝟕𝟔𝟐
𝒌𝒈
Ahora calculamos la entropía real
𝒌𝑱
𝑺𝟗 = 𝟕, 𝟓𝟓𝟓𝟔
𝒌𝒈. 𝑲
𝑘𝐽
𝑆𝑓 = 1,0261
𝑘𝑔. 𝐾
𝑘𝐽
𝑆𝑔 = 7,6691
𝑘𝑔. 𝐾
𝑆13𝑠 −𝑆𝑓
𝑥𝑠 = = 0,9830
(𝑆𝑔 −𝑆𝑓 )
𝒌𝑱
ℎ13 = 𝟐𝟔𝟔𝟔, 𝟐𝟓𝟒𝟒
𝒌𝒈
𝒎𝟑
𝒗𝟏 = 𝟎, 𝟎𝟎𝟏𝟎𝟐𝟔
𝒌𝒈
𝑚3
(0,001026 ) (400 − 40)𝑘𝑃𝑎
𝑘𝑔
𝑊𝑝1 =
0,95
𝑘𝐽
𝑊𝑝1 = 0,3888
𝑘𝑔
Ahora se calcula la entalpia h2 con el trabajo de bomba encontrado anteriormente
ℎ2 = ℎ1 + 𝑊𝑝1
𝑘𝐽 𝑘𝐽
ℎ2 = 317,62 + 0,3888
𝑘𝑔 𝑘𝑔
𝒌𝑱
𝒉𝟐 = 𝟑𝟏𝟖, 𝟎𝟎𝟖𝟖
𝒌𝒈
𝒎𝟑
𝒗𝟑 = 𝟎, 𝟎𝟎𝟏𝟎𝟖𝟒
𝒌𝒈
𝑚3
(0,001084 ) (12000 − 400)𝑘𝑃𝑎
𝑘𝑔
𝑊𝑝2 =
0,95
𝑘𝐽
𝑊𝑝2 = 13,2362
𝑘𝑔
Ahora se calcula la entalpia h4 con el trabajo de bomba encontrado anteriormente
ℎ4 = ℎ3 + 𝑊𝑝2
𝑘𝐽 𝑘𝐽
ℎ4 = 604,66 + 13,2362
𝑘𝑔 𝑘𝑔
𝒌𝑱
𝒉𝟒 = 𝟔𝟏𝟕, 𝟖𝟗𝟔𝟐
𝒌𝒈
c. Determine las fracciones de vapor extraído de la turbina hacia los calentadores de agua de
alimentación.
Balance de energía en el calentador cerrado, se tiene que:
𝑑𝐸 𝑑𝐸
𝐸𝑖𝑛 − 𝐸𝑜𝑢𝑡 = ; 𝑐𝑜𝑚𝑜 𝑠𝑒 𝑒𝑠𝑡𝑎 𝑒𝑛 𝑒𝑠𝑡𝑎𝑑𝑜 𝑒𝑠𝑡𝑎𝑐𝑖𝑜𝑛𝑎𝑟𝑖𝑜 =0
𝑑𝑡 𝑑𝑡
𝐸𝑖𝑛 − 𝐸𝑜𝑢𝑡 = 0
ℎ4 (1) + ℎ11 (𝑦) − ℎ5 (1) − ℎ6 (𝑦) = 0
ℎ5 −ℎ4
𝑦=
ℎ11 − ℎ6
𝐾𝐽
(697 − 617,8962)
𝑘𝑔
𝑦=
𝐾𝐽
(3187,0356 − 697)
𝑘𝑔
𝒚 = 𝟎, 𝟎𝟑𝟏𝟕𝟔𝟖
𝑞𝑛𝑒𝑡𝑜 = 𝑊𝑛𝑒𝑡𝑜
𝐾𝐽
𝑞𝑛𝑒𝑡𝑜 = 1138,691 = 𝑊𝑛𝑒𝑡𝑜
𝑘𝑔
𝑞𝑜𝑢𝑡
𝜂 = 1− = 0,35864
𝑞𝑖𝑛
𝜂 = 35,864%
e. Determine el flujo másico requerido para una entrega de potencia neta de 𝑊𝑊neto.
𝑾̇𝒏𝒆𝒕𝒐 = 𝒎̇ ∗ 𝑾𝒏𝒆𝒕𝒐
𝑊̇𝑛𝑒𝑡𝑜
𝑚̇ =
𝑊𝑛𝑒𝑡𝑜
11000𝐾𝑊
𝑚̇ = 𝐾𝐽
1138,691
𝑘𝑔
𝒌𝒈
𝒎̇ = 𝟗, 𝟔𝟔
𝒔
Cada uno de los cálculos realizados pueden ser consultados en el siguiente Excel: Cálculo Taller
final motores - Cantillo, Luquez
4. En el presente inciso, se solicita realizar un análisis paramétrico de las variables de salida del
ciclo: Eficiencia térmica, tasa de calor y flujo de calor. Para ello, se recurrió a una matriz de
Excel donde se determinaron nuevamente tanto propiedades como cálculos energéticos. Para
ello, se hizo variar las variables de entrada de la tabla de inciso anterior, en un
−20% , −10%, +10%, +20% .
Se mantuvo fijo el valor de los rendimientos isentrópicos tanto de bombas como turbinas para
visualizar cómo cambiaban las variables deseadas en función de las condiciones de entrada. Es
decir, se deseaba estudiar el cambio en las variables de respuesta en condiciones similares. El
archivo de Excel puede ser consultado en el siguiente hipervínculo: Cálculo Taller final motores -
Cantillo, Luquez
Los valores obtenidos, se encuentran en la siguiente tabla:
η vs % variación
R² = 0,9172
41
40 39,894
Eficiencia térmica (%)
39
38
37
36 36,214
35,177 35,864
35,477
35
34
-20 -15 -10 -5 0 5 10 15 20
% Variación
flujo másico vs % variación R² = 0,9864
12 10,599 10,653
9,66
10 8,721
Flujo másico (kg/s)
7,751
8
0
-20 -15 -10 -5 0 5 10 15 20 25
% Variación
1260 1239,06748
1240
Tasa de calor (kJ/kg)
1220
1200
1180
1160 1141,61763
1135,19635 1138,69160
1140
1120 1135,28401
1100
-20 -15 -10 -5 0 5 10 15 20 25
% Variación
Flujo de calor vs % variación R² = 1
13200
14000
12100
12000 11000
9900
Flujo de calor (kW) 10000 8800
8000
6000
4000
2000
0
-20 -15 -10 -5 0 5 10 15 20 25
% Variación
• Conforme se aumentan las presiones en las etapas de turbina, puede observarse un aumento
en la eficiencia térmica del ciclo. Sin embargo, también puede observarse un aumento en el
trabajo consumido en las etapas de bomba.
• Conforme aumenta 𝑊̇𝑛𝑒𝑡𝑜 , es requerido un mayor flujo másico de vapor para garantizar dicho
valor.
• Conforme menor sean las presiones, mayor será la cantidad de calor que debe aportarse
cuando el fluido de trabajo pasa por la caldera y mayor será la cantidad de calor que debe
retirarse en el condensador.
• Conforme se aumenta la potencia neta del ciclo, aumenta también el trabajo neto, puesto que
es un ciclo: 𝑤𝑛𝑒𝑡𝑜 = 𝑞𝑛𝑒𝑡𝑜 por lo tanto
𝑚̇𝑊̇𝑛𝑒𝑡𝑜 = 𝑚̇𝑄̇𝑛𝑒𝑡𝑜
𝑊̇𝑛𝑒𝑡𝑜 = 𝑄̇𝑛𝑒𝑡𝑜