LABORATORIO DE INGENIERÍA QUÍMICA I

29 de noviembre del 2022

INFORME CALDERA

Grupo 2
Maria Luisa Castaño Mazuera; Johan David Delgado Sandoval; Alejandro Pabón Huérfano;
Daniel Antonio Serna Gómez.

1. Objetivos1
1.1. General: Estudiar la caldera como intercambiador de calor, su
funcionamiento térmico y su rendimiento calórico.
1.2. Específicos
I. Calcular la cantidad de combustible gastado por ciclo y por unidad de
vapor producido.
II. Calcular el calor cedido por el combustible, en la combustión y en el
paso por el banco de tubos.
III. Determinar los coeficientes peliculares para los gases de combustión
y en el vapor de agua.
IV. Estimar los coeficientes totales de transferencia de calor, Uc, Ud y Rd.
V. Calcular la eficiencia de la caldera.
VI. Graficar temperatura de salida de gases de chimenea y presión de
salida del vapor contra tiempo. Analizar los periodos de trabajo de
una caldera y sus niveles de estabilización.
VII. Comparar el funcionamiento térmico con el de otros equipos.

2. Datos y cálculos.
Para el desarrollo de la práctica se registraron los datos (Tabla 1) que relacionaban las
condiciones de salida de vapor de la caldera, así como las condiciones de operación de
esta. Es de notar que el combustible que utilizó el equipo fue gas natural proporcionado
por Efigas2 y el agua que se usó fue la proveniente del acueducto, la cual se trató en una
columna de intercambio catiónico para reducir la cantidad de minerales que esta
pudiera tener. Las casillas que se encuentran resaltadas indican el momento en que la
caldera se encontraba operando.
 Medición P caldera
Medición T
Tiempo gas Q T vap P vap
Ciclo gas final Caldera
[min] inicial [kg/h] [°C] [Barg]
[m3] [°C] [psig] [Barg]
[m3]
5 170 95 6,55 63 165,1 6
10 170 110 7,58 72 168,8 6,76
15 150 102 7,03 52 168 6,56
1 8387,0 8388,0
20 140 100 6,89 52 166,3 6,24
25 135 90 6,21 48 164,7 5,92
30 135 82 5,65 46 162,4 5,59
5 170 90 6,21 67 165,2 6,03
10 170 110 7,58 54 169,2 6,74
15 145 102 7,03 42 167,3 6,44
2 8388,0 8389,0
20 140 90 6,21 47 165,3 6,06
25 140 83 5,72 50 163,1 5,66
30 130 79 5,45 52 161.3 5,41
5 170 85 5,86 122 163,4 5,8
10 170 100 6,89 104 166,3 6,31
15 175 110 7,58 89 168,4 6,64
3 8389,0 8390,5
20 150 110 7,58 56 168,1 6,58
25 140 100 6,89 56 166,4 6,24
30 140 90 6,21 54 164,7 5,95
Tabla 1. Registro de datos a lo largo de la práctica
Como se puede notar en la tabla, el gasto volumétrico del gas en los dos primeros dos ciclos
fue de 1m3 mientras que en el último ciclo fue de 0.5 m3.
Para determinar el valor de la capacidad calorífica además del peso molecular de la mezcla
del gas es necesario saber la composición de este, pero es de notar que esta varía de
yacimiento a yacimiento 3, por lo cual se tomaron los valores típicos de composición en
volumen con que se encuentra el gas natural domiciliario que es4:

• Metano (CH4): 81.86%.

• Etano (C2H6): 11.61%.
• Propano (C3H8): 1.92%.
• I-Butano (C4H10): 0.23%.
• N-Butano (C4H10): 0.22%.
• Nitrógeno (N2): 0.90%.
• Dióxido de carbono (CO2): 3.18%.
Asumiendo que la fracción de volumen es igual a la fracción molar, se puede determinar la
densidad del gas:
 kg
 =  ( vi  PM ) = 0.775 [1]
m3

2.1. Combustible gastado en cada ciclo y por unidad de vapor producido

Vgastado
Fcombustible = [2]
ttotal

Donde,
Vgastado = V final − Vinicial [3]

Vgastado = 8390.5m3 − 8388m3 = 2.5m3 [4]

De esta manera, reemplazando los valores del Vgastado en [2] se

obtiene que

2.5m3 m3
Fcombustible = = 0.0278 [5]
90 min min
Ahora,
kg
mcombustible = Fcombustible  0.775 [6]
m3
m3 kg kg
mcombustible = 0.0278  0.775 3 = 0.0216 [7]
min m min
De esta manera, el volumen de combustible gastado es
Vc / c = Fcombustible  ti [8]
m3
Vc / c = 0.0278  90 min = 2.502m3 [9]
min
El calor ganado por el agua es igual a:

 1
( ) 
 
(a + b T + c T + d1T 3 ) dT +  (a + b2T + c2T 2 + d 2T 3 + e2T 4 ) dT + v 
Tsat Toc
Q agua = magua   1 1 1
2
2 [10]
 Pmagua Tent Tsat


Donde, el λv a la presión promedio de operación de la caldera de 100psig, es igual a

880 BTU/lb.5
 Y, las constantes de las integrales se pueden observar en la tabla 2.

a1 1.82964*101
b1 4.72118*10-1
c1 -1.33876*10-3
d1 1.31424*10-6
a2 3.40471*101
b2 -9.65064*10-3
c2 3.29983*10-5
d2 -2.04467*10-8
e2 4.30228*10-12

Tabla 2. Constantes de las integrales para la ecuación 101.

Luego, el calor cedido por el combustible es
 
Q c = m combustible  Pc ,total [11]
kJ
Teniendo en cuenta que Pc ,total = 39900 2
kg

kJ
Qc = 0.0216  39900 = 861.84 [12]
min
La presión atmosférica en Manizales es 585mmHg=11.312psi1, de esta manera

Pmax = 168.8 + 11.312 = 180.112 psi [13]

Pmin = 161.3 + 11.312 = 172.612 psi [14]

Como el vapor se encuentra en estado saturado, con [13] y [14], se puede establecer
las temperaturas máxima y mínima de saturación
189.5185C + 187.5913C
Tsat = [15]
2
Tsat = 188.555C [16]

A partir del balance de energía [10] se calcula la cantidad de agua evaporada


 Q agua
m agua = [17]
  ( a1 + b1T + c1T 2 + d1T 3 ) dT + v
1 Tsat

Pmagua Tent
Donde, el calor ganado por el agua se calcula a partir del calor cedido por la combustión y
las pérdidas ocasionadas por convección y por radiación.

El calor por convección se calcula de la siguiente manera

 Qconvección = Ac  ho  (Tcomb − T ) [18]
Y Ac está determinada por

Ac =   Dc  Lc [19]
Ac =  12inch 160inch = 6031.86inch = 3.8915m 2
[20]
Ahora, reemplazando en [18]

BTU BTU
Qconvección = 41.8878ft 2  0.11  (120 − 157.5 )  F = −172.787 [21]
h  ft  F
2
h
Luego, el calor por radiación se calcula así

Qradiación = Ac      (Tcomb 4 − T4 ) [22]

 0.3  (1204 − 157.54 ) K = −27.0067W = −92.0928

W BTU
Qradiación = 3.8915m2  5.67 x10−8 [23]
m K
2 4
h
Las pérdidas totales son la suma del calor perdido por convección y radiación

Q pérdidas = Qconvección + Qradiación [24]

BTU [25]
Qpérdidas = ( −172.787 ) + ( −92.0928 ) = −264.88
h
Así, el calor ganado por el agua es entonces

Qagua = Q pérdidas − Qc
[26]

BTU
Qagua = 264.88 − 861.84 = 596.96 [27]
h
Ahora, reemplazando los valores en [17] se obtiene que

kg [28]
magua = 0.4336
min
Por último, con el flujo másico del agua [28] encontrado y el flujo másico de combustible [7]
se establece la relación buscada

Cantidad combustible gastado mcombustible

= [29]
Cantidad vapor producido magua

2.502  0.4336 [30]

Cantidad vapor producido = = 50.2253
0.0216
 2.2. Calor cedido por el combustible en la combustión, y en el paso por el banco de
tubos

El calor cedido por el combustible es hallado a través de la siguiente expresión:

[31]

Y para encontrar los valores de cada uno de los anteriores términos se usan las siguientes
ecuaciones, de la [32] a la [34]:

Para hallar la entalpía de los productos, es necesario calcular el flujo de aire que entra a la
caldera, lo cual se hace a través del cálculo del oxígeno teórico entrante y de la cantidad en
exceso de aire en la caldera. El aire en exceso en la caldera es cuantificado a través del análisis
de Orsat que se realiza en el segundo ciclo de la caldera, el cual, para este caso es de 25%. Y
para calcular el aire teórico y el aire que entra se usan las siguientes ecuaciones [36] y [37]:

Ahora bien, para calcular la entalpía de los productos es necesario realizar los balances
descritos a continuación en [38], [39] y [40], para los cuales se requieren las fracciones
másicas de los productos de combustión los cuales se obtienen de los resultados obtenidos
en el análisis de Orsat que se muestran en la siguiente figura:
 Figura 1. Resultados obtenidos en el análisis de Orsat.

Para encontrar el flujo de gases de combustión a la salida, es necesario realizar los balances con
base en las siguientes reacciones llevadas a cabo en la combustión del gas natural:

𝐶𝐻4 + 2𝑂2 → 𝐶𝑂2 + 2𝐻2 𝑂

3
𝐶𝐻4 + 𝑂2 → 𝐶𝑂 + 2𝐻2 𝑂
2
𝐶3 𝐻8 + 5𝑂2 → 3𝐶𝑂2 + 4𝐻2 𝑂
 7
𝐶3 𝐻8 + 𝑂2 → 3𝐶𝑂 + 4𝐻2 𝑂
2
5
𝐶2 𝐻6 + 𝑂2 → 2𝐶𝑂 + 3𝐻2 𝑂
2
7
𝐶2 𝐻6 + 𝑂2 → 2𝐶𝑂2 + 3𝐻2 𝑂
2
𝐶4 𝐻10 + 13𝑂2 → 8𝐶𝑂2 + 10𝐻2 𝑂
9
2𝐶4 𝐻10 + 𝑂2 → 4𝐶𝑂 + 5𝐻2 𝑂
2

Tabla 3. Flujo de gases de combustión

Flujo 1 Flujo 2
Flujo 3 Flujo 1 N2 Flujo 3 Flujo 3 Flujo 3 Flujo 3
O2 O2
Ciclo O2 teórico CO2 N2 CO H2O
teórico teórico
(Kmol/s) (Kmol/s) (Kmol/s) (Kmol/s) (Kmol/s) (Kmol/s)
(Kmol/s) (Kmol/s)

1 3,02E-04 6,36E-05 2,36E-04 2,36E-04 6,13E-07 2,36E-04 9,28E-06 2,03E-06

2 3,27E-04 6,45E-05 2,54E-04 2,54E-04 6,63E-07 2,54E-04 9,87E-06 2,18E-06

3 3,96E-04 6,87E-05 2,71E-04 2,71E-04 6,98E-07 2,71E-04 1,06E-06 2,32E-06

Ahora, en cuanto a la entalpía de los reactivos es necesario calcular la capacidad calorífica

(Cp), para las diferentes sustancias que componen el gas natural, como se muestra en la
tabla 4.

Tabla 4. Valores de las constantes a,b,c,d y e para el cálculo del Cp

Cp
Sustancia a b c d e
(J/mol)

-
CH4 34,942 -0,039957 0,00019184 0,00000015 3,932E-11 36,206
3

C2H6 28,146 0,043447 0,00018946 -1,908E-07 5,335E-11 52,918

C3H8 28,277 0,116 0,00019597 -2,327E-07 6,867E-11 74,035

C4H10 6,772 0,34147 -0,00010271 -3,685E-08 2,043E-11 97,680

-4,3116E-
N2 29,342 -0,0035395 1,0076E-05 2,594E-13 29,066
09

CO2 27,437 0,042315 -1,9555E-05 3,9968E-09 -2,99E-13 38,307

 Ahora, para calcular el flujo de cada compuesto presente en el gas natural, se usa la
ecuación [41] y se muestran los resultados en la tabla 5.

[41]

Tabla 5. Flujo de gas natural

Flujo CH4 Flujo C2H6 Flujo C3H8 Flujo C4H10 Flujo N2

Ciclo
(Kmol/s) (Kmol/s) (Kmol/s) (Kmol/s) (Kmol/s)

1 1,98E-05 2,74E-06 7,54E-08 1,14E-08 1,61E-07

2 1,98E-05 2,74E-06 7,54E-08 1,14E-08 1,61E-07

3 2,98E-05 4,11E-06 1,13E-07 1,71E-08 2,41E-07

Finalmente, en la tabla 6 se reporta el calor cedido por el combustible para cada uno de los
3 ciclos realizados durante la práctica.

Tabla 6. Calor cedido en cada ciclo

Flujo de aire Calor cedido por

Ciclo que entra el combustible
(m3/s) (KJ/min)
1 0,013623 1798,432
2 0,014722 1793,128
3 0,018180 1769,324

2.3. Coeficientes peliculares para los gases de combustión y en el vapor de agua

para calcular el coeficiente convectivo de los gases de combustión se hace de la siguiente
manera:
 Di     [42]
hio = cci   
mezcla

 De    
Donde
fi [43]
cci = 2.44  mgas 0.8 
Di1.8
  Cp 
0.4 [44]
fi =  mezcla   kmezcla 0.6
 mezcla 

Donde el Di es de 2 pulgadas y el De el diámetro externo de los tubos con 2,375 pulgadas

Cp K 𝜇
1,453 0,026 0,0136
Remplazando los valores en las ecuaciones tenemos los siguientes resultados.
 BTU BTU  0,7253
fi   
 cp h  ft  F 
 BTU BTU  25,36534212
 kg  
cp h  ft  F
cci   
 min inch 
 
 BTU  21,3603
hio  
 h  ft  F 
2

Y el coeficiente para el vapor del agua por fuera de los tubos se hace de la
siguiente forma:
 Qc    Po  
6.7 0.17 12 10 [45]
 Po   Po 
ho = 0.00658  Pc 0.69
   1.8    + 4  + 10    
 A    Pc   Pc   Pc  
Donde el Qc es el calculado en el objetivo 1 la Po es la presión de vapor del
sistema.
La presión critica se obtuvo de las tablas de vapor a la Temperatura promedio de
la caldera 157,5 °C la cual fue de 5,722559294 bares

Para el área se calculó mediante

[46]

Con la información dispuesta sobre la caldera el resultado fue

A = 3  38 1.80    2
[47]
A = 107, 4149635 ft 2

El valor de ho fue de 1,42665E+22

2.4. Coeficientes totales de transferencia de calor, Uc, Ud y Rd

1 1 1 [48]
= +
Uc hio h0
 Qc [49]
Ud =
A  MLDT
1 1 [50]
Rd = −
Uc Ud

Utilizando los coeficientes de película calculados en el inciso anterior y despejando

Uc

Uc = 21,3602881 [Btu/(hr-ft2°F)]. []

Para calcular el Ud se necesita MLDT que es la temperatura media logarítmica y se

encuentra de la siguiente manera

(Tsat − Tent ) − (Tp max − Tp min )

MLDT =
 (Tsat − Tent ) 
ln   [51]
 (T
 p max − Tp min ) 
MLDT = 118.12 F
Y reemplazando valores tenemos que

Ud = 16,17183 [Btu/(hr-ft2°F)]. [52]

Ya con los dos valores de los coeficientes de transferencia se puede encontrar el

Rd

Rd = −0,01502008207 [Btu/(hr-ft2°F)]. [53]

2.5. Eficiencia de la caldera

La eficiencia de la caldera se define como
Qagua [54
= 100
Qc ]
[
596.96
= 100 = 69.27% 55]
861.84
 Figura 2. Eficiencia teórica de una caldera vs porcentaje de aire. Imagen tomada de:
[Link]

2.6. Gráficas de variación de temperatura y presión a través del tiempo de operación

Figura 3. Gráfica temperatura vs tiempo de toda la operación.

 Figura 4. Gráfica presión vs tiempo de toda la operación.

2.7. Comparar el funcionamiento térmico con el de otros equipos

Comparando la eficiencia de una caldera pirotubular, se tiene que las calderas
pirotubulares operan habitualmente en un rango de eficiencia entre 88% y 92%6 si
se compara con la eficiencia obtenida en la práctica que fue del 69.27%, se
encuentra que el equipo no se encuentra operando en óptimas condiciones y que
se le debería realizar un mantenimiento, para obtener mejores resultados.

3. Análisis de Resultados
• Observando el porcentaje obtenido de la eficiencia de la caldera, se
encuentra que el procedimiento fue realizado de manera correcta y que la
caldera se encuentra operando en buenas condiciones, pero se podría
realizar un mantenimiento para que su funcionamiento mejore.
• Al analizar el volumen de combustible que se gastó en la práctica, se llega a
una conclusión de que es un proceso de costo bajo para la importancia que
tiene ya que gracias a todo el vapor que se produce, se pueden alimentar
otros equipos que se encuentran en el laboratorio.
• En la figura 3 podemos observar que el rango de la temperatura de vapor va
entre 161-169 °C y los mayores valores de temperatura se obtienen mientras
la caldera se encuentra encendida.
• En la figura 4 podemos observar que el rango de la presión de vapor va entre
5.4 - 6.7 bares, los mayores valores de presión se obtienen mientras la
caldera se encuentra encendida.
 • Viendo los anteriores comportamientos de la temperatura y la presión de
vapor, es de suma importancia poder tener los valores de operación para
poder llevar un riguroso análisis de cuánto tiempo se debe mantener
encendida la caldera para conseguir los valores necesarios del vapor que
alimentara el proceso que estemos realizando, debido a que la eficiencia de
los procesos industriales puede llegar a depender en gran medida de este.
• El calor cedido por el combustible es similar en cada ciclo realizado durante
la práctica, en parte esto se explica porque los tiempos de cada ciclo fueron
totalmente iguales, sin embargo, el ciclo 3 es un poco diferente de los dos
anteriores debido a que este presento mayor variación en los datos.
• Debido a que, en el proceso el aire en exceso fue del 25% es posible afirmar
que el valor es medianamente adecuado debido al análisis de los gases de
combustión Orsat, ya que en dicho análisis se encuentran cantidades de CO2
y O2, productos que se obtienen cuando se da una combustión completa,
por otro lado, no hay presencia de monóxido de carbono en los gases, lo cual
es bueno porque el monóxido de carbono aparece cuando la reacción no se
completa con éxito.
• Comparando los resultados obtenidos de Ud y Uc vemos que no son valores
tan lejanos por lo que la transferencia de calor no se esta viendo afectada
por situaciones como falta de mantenimiento del equipo sin embargo es
normal que a través de los años el equipo empiece a transferir menos, esto
debido al desgaste que este va sufriendo.
• Tener en constante monitoreo estos coeficientes transferencia debe ser una
buena practica ingenieril aplicada ya que con estos podemos tener una visión
general de cómo funciona nuestro equipo y tener un análisis mas profundo
al momento en que este presente menos eficiencia o diferentes fallas, se
podría descartar que sea por motivos de ensuciamiento.
• Viendo la figura 2 y sabiendo que el exceso de aire reportado por nosotros
es de 25% el rango que podemos ver de eficiencia es del 60% al 70% sobre
el cual se encuentra la eficiencia calculada en el objetivo cinco, la cual fue de
69%,esto nos lleva a decir que se tiene una buena consistencia de la practica
con lo visto en la realidad ingenieril.
4. Conclusiones.
4.1. Debido a la aplicabilidad y la importancia que representa el uso del vapor
saturado en el desarrollo de la industria, el conocimiento adquirido durante la
realización de esta práctica fue de gran aporte para nuestra formación como
ingenieros.
4.2. Es importante verificar que las propiedades del agua de alimentación a la
caldera estén dentro de los requerimientos establecidos, garantizando de esta
 manera un buen desempeño del equipo y así evitar daños como la corrosión o las
incrustaciones; es por esto que la caldera que se utiliza en el laboratorio tiene una
columna de intercambio catiónico.
4.3. Se notó que la duración de los ciclos de la caldera se veía restringidos por el uso
del vapor en los diferentes equipos del laboratorio, ya que habitualmente cada ciclo
tiene una duración de 45 minutos y los tres ciclos observados tuvieron una duración
de 30 minutos porque en ese momento se encontraban en funcionamiento 2
equipos que requerían del uso del vapor generado en la caldera.
4.4 con la practica se pudieron explorar diversas funciones que debe tener el
ingeniero a la hora de enfrentarse con una caldera y que esta no sea un equipo
extraño para él, sino que calculando los calores, el consumo de combustible, los
coeficientes de transferencia, se pueda sacar provecho del equipo y que no sea algo
raro para el ingeniero a la hora de salir a su mundo laboral.
5. Recomendaciones.
Sería interesante poder realizar la prueba con distintas fuentes de energía diferentes
al gas debido a que esto puede presentar un análisis, más enriquecedor sobre temas
de eficiencia de combustibles y contaminantes generados, realizando diversas
pruebas Orsat7.

Una recomendación que nos parece interesante es que podría ser muy enriquecedor
para los estudiantes poder conocer un poco más de la mecánica del mantenimiento
de los equipos, ya que, si bien no es algo que nos competa como ingenieros, tener
el conocimiento y poder actuar en un caso de emergencia nos haría tener mas
aptitudes para ser un ingeniero versátil, que es algo muy criticado en nuestros días.

6. Referencias.
1. Mesa Rueda FA, Peña Hernandez GA. Manual de Laboratorio de Operaciones Unitarias II.
Manizales: 2006.
2. EFIGAS S.A. E.S.P. [Internet]. [cited 2022 Nov 29];Available from:
[Link]
3. ¿Qué es el Gas Natural? 2014;
4. ¿Qué es? - Grupo Vanti.
5. Tablas de vapor | Armstrong International [Internet]. [cited 2022 Nov 28];Available from:
[Link]
6. CALDERA (máquina). Qué es, tipos, partes, funcionamiento [Internet]. [cited 2022 Nov
18];Available from: [Link]
7. Análisis de Orsat para gases de combustión [Internet]. [cited 2022 Nov 18];Available from:
[Link]