John Sebastián Giraldo López johnsgiraldo@gmail.

com

Juan Camilo González Várelas juancamo_3@[Link]

Daniel Alejandro Valencia H. davh82@[Link]

SEGUNDO PARCIAL

Daniel Ricardo Vélez leinad_ri@[Link]

COMBUSTION DEL GAS

PRO
FESOR ANDRÉS AMELL ARRIETA

PRIMER PUNTO

Parte A

Un horno de fusión de cobre opera con gas natural de Guajira, un

diagnostico para optimizar su operación encontró que su potencia térmica
es PT, la temperatura con que salen los gases de combustión es T h, la
temperatura de la pared del horno es TP, la relación de CO/CO2 en humos
seco es RO. La sección transversal del horno es a x b y su longitud L en
metros, la temperatura ambiente es 30ºC, encontrar:

1. La composición de los humos seco.

2. La composición de los humos húmedo.

3. El tipo de combustión con que se opera el horno.

4. El factor de aireación.

5. La fracción de aire utilizado y no utilizado.

6. La eficiencia de combustión con base al PCi.

7. La eficiencia útil con base al PCi.

8. El consumo de gas natural en m3st / h.

Tabla 1. Valores conocidos

D K
atos
Valo 300 20 1 0 30 2 0,99 1,97 7 9 1 0 8,5
r 0 100 ,25 ,52 5 5 ,47 ,44 ,29

Solución:
1. Se estimaron la composición de los humos secos; lo primero fue hallar la
manera como se utilizó el aire (cuanto reacciona con el combustible y
cuanto no reacciona) para caracterizar que tipo de combustión se está
llevando a cabo, para lo anterior se cuento con las siguientes formulas:

(1)

(2)

Dado que (fracción de aire utilizado) la combustión puede ser

Reductora o Mixta.

Ahora se debe hallar :

(3)

(4)

(5)

(6)

A continuación se presentan las ecuaciones para el cálculo de cada uno de

los términos de la anterior ecuación.

(7)

(8)

(9)

(10)

(11)

(12)

Una vez calculados estos valores se procede a hallar los porcentajes de

cada uno de los componentes respecto al volumen de humos secos.

(13)

(14)
 (15)

(16)

(17)

2. Ahora se desea determinar el volumen de humos húmedos y su

composición

(18)

(19)

(20)

Ahora se determinan de la misma forma que el punto anterior los

porcentajes de cada componente de los humos húmedos.

(21)

(22)

(23)

(24)

(25)

(26)

3. Ya se puede determinar el tipo de combustión con la cual opera en

horno, tenemos que:

(27)

Se puede observar de la anterior ecuación que horno opera con una mezcla
rica con un factor de aireación de ; dado que y , el
tipo de combustión se denomina Mixta.

4. Retomando la ecuación 27 el factor de aireación es

5. La fracción de aire utilizado (ecuación 2) es y la de no utilizado

(ecuación 6) es

6. La eficiencia de combustión se calcula con la siguiente fórmula:

 (28)

Siendo

(29)

(30)

Siendo : Calor latente de vaporización del agua

: Porcentaje de agua presente en los humos

: Caudal de gas,

(31)

Siendo

(32)

(33)

(34)

Retomando la ecuación 29, se reemplazan todos los valores y se tiene que:

(35)

Ahora se reemplaza en la ecuación 28. Entonces:

7. La eficiencia útil se define como:

(36)

Falta por determinar las pérdidas de calor por las paredes del horno, las
cuales se definen como:
 (37)

Para las paredes laterales,

53,56 kw (38)

Para la pared superior,

(39)

Para la pared inferior,

(40)

(41)

Retomando la ecuación 36 y remplazando los valores se tiene que:

(42)

8. El consumo de gas natural es:

(43)

Parte B

Si después del diagnostico el horno se optimiza, incorporando un sistema de

calor autoregenerativo, cambiando el aislamiento y ajustando la relación
aire combustible, se encuentra que la temperatura de los humos es Th’, la
temperatura de pared TP’ y la relación CO/CO2 en humos seco es RO’ y el
factor de aireación es 1.2, estimar nuevamente los anteriores numerales.

D PT VC VH VN K CO R Th Tp L*a*b n Y( Va
ato Potencia O2 2O 2 2 0 N2)
s termica
V 300kw 0.99 1.97 7 2 8. 0 10 5 2.5*2* 1. 1 9
alor 5 5 .47 .52 5 0 0 2 2 .29 .44

1.1)Estimación de la composición de humos secos

Partiendo de la primera parte del algoritmo se determina el contenido de

los gases de combustión en los humos secos para el segundo caso.
 RO * CO2 = CO = 0 (1)

VCO2
CO2 = = 0.995 (2)
1 + RO

VH 2 O * K
H 2O = = 1.975 (3)
K + RO

VH 2 O * RO
H2 = =0 (4)
K + RO

RO: Se define como la relación entre CO/ CO2

K: Constante de reacción del gas guajira (2.52)

1.2)Con el fin de calcular la composición de los humos húmedos se utiliza la

ecuación general de eficiencia :

n = n1 + n2 (5)
Para la cual se calcula n1, correspondiente a la fracción de aire
utilizado

 VH 2 O * K   2 + RO 
 + VCO2 − 2 yO2
 K + RO   1 + RO  (6)
n1 = = 1.00
0.42Va

Despejando de (20)
n2 = 0.2 La cual corresponde a la fracción de aire no utilizada
en la combustión

Después de hacer los respectivos cálculos para las fracciones de aire

utilizado y no utilizado se prosigue con la parte del algoritmo para calcular
mediante un balance de masa la cantidad de nitrógeno proveniente del
aire que reacciona (n1. Va), del aire que no reacciona (n2, Va) y del gas
combustible (YN2), y la cantidad de oxigeno básicamente proviene del aire
que no reacciona.

N 2 = 0.79(n1 + n2 ) * Va + yN 2 = 10.23
O2 = 0.21n2Va = 0.396

De los cálculos obtenidos en la primera parte se observa la presencia de

los siguientes gases en los humos CO2, N2.02; debido a esto se concluye
que la combustión es de tipo oxidante ; luego se calcula la cantidad de
humos secos mediante la siguiente ecuación con el fin de poder hallar las
respectivas composiciones:

1.1)
 VHO = CO2 + N 2 + O2 = 11.63 (7)

CO2 N2 O2
%CO2 = = 8.55 %N 2 = = 88 %0 2 = = 3 .4
VHO V HO V HO

1.2) Composición en los humos húmedos (VHH): La cantidad de humos

húmedos corresponde a la sumatoria de los humos secos para una
combustión tipo oxidante y la cantidad de agua calculada en la primera
parte.

V HH = V HO + H 2 O = 13.605 (8)

CO2 CO2 N2
%CO2 = = 7.31 %H 2O = = 14.51 %N2 = = 75.7
VHH VHH V HH

O2
%O2 = = 2.91
V HH

1.6) Determinación de pérdidas de calor por los humos:

Las pérdidas de calor en los humos se debe a los calores sensibles de cada
uno de los gases de combustión , y a la perdida por calor sensible del agua,
ahora las perdidas por compuestos combustibles sin quemar no se tienen
en cuenta ya que las fracciones de CO e H2 son cero para este caso.

LH 2O = H 2 O × h fg × QPCS

kWh
0.55
( )
m H3 2O st
; Calor de evaporización del agua a 15ºC (hfg)

PT
QPCS = = 27.4m 3 / h = 7.611*10 −3 m 3 / s ; Caudal basado en el poder
PCS
calorífico superior

PCS=Poder calorífico superior (39.40MJ/m3)

mH3 2O kWh
3
m gas
LH 2O = 1.975 × 0.55 3 × 27.4 = 29.76kW
m gas (
mH 2O st) h
(24)
 ∆H = VHH × ( Cpvol ) st × ( ∆T ) × QPCS
Pérdidas debidas al calor latente de los gases de combustión

kJ
( Cp vol ) st = 1.4
(m )
3
hh st ׺ C
; Capacidad calorífica volumétrica a condiciones

estándar

3
3
mhh kJ −3
mgas
∆H = 13.605 3 × 1.4 3 × (100 − 15)º C × 7.61*10
mgas mhh st ׺ C ( ) s
∆H = 12.32kW

 PT − ( LH 2O + ∆H ) 
ncombution =   *100% = 85.6%

 P T 

1.7) Eficiencia útil:

La eficiencia útil con base en el poder calorífico inferior se calcula teniendo
en cuenta las perdidas tanto en los productos de combustión como el calor
transferido a los alrededores por convección y radiación

5  T1 + 273  4  T 2 + 273  4 
Q p = h × (T1 − T2 ) 4
+ 4.88 × ε ×   −  
 100   100  
Q p : Calor transferido a los alrededores (kcal / m 2 h)
h : Ceoficiente de convección natural
Pared superior = 2.8; Pared lateral = 2.2; Pared inferior = 1.5
T1 : Temperatura exterior de la pared del horno (º C )
T2 : Temperatura del aire ambiente alrededor del horno (º C )
ε : Emisividad de la superficie exterior de la pared del horno

L = 2.5m, a = 2m, b = 2m

Para las paredes laterales:

5  100 + 273  4  30 + 273  4 

Q p = 2.2 × (100 − 30) 4
+ 4.88 × 0.8 ×   −  
 100   100  
kcal
Q p = 872
m2h
 kcal kcal kJ kJ 1h
Q p = 2 × (2m × 2.5m) * 872 2
= 8720.6 × 4.1868 = 36511.1 × = 10.14kW
m h h kcal h 3600s

kcal kcal kJ kJ 1h
Q p = 2 × (2m × 2.m) * 872 2
= 6976 × 4.1868 = 29207.1 × = 8.11kW
m h h kcal h 3600s

Pared superior:
5  100 + 273  4  30 + 273  4 
Q p = 2.8 × (100 − 30) 4
+ 4.88 × 0.8 ×   −  
 100   100  
kcal
Q p = 993.55
m2h

kcal kcal kJ kJ 1h
Q p = (2m × 2.5m) * 993.55 2
= 4967.7 × 4.1868 = 29799 × = 5.77 kW
m h h kcal h 3600s

Pared inferior:

5  100 + 273  4  30 + 273  4 

Q p = 1.5 × (100 − 30) + 4.88 × 0.8 × 
4
 −  
kcal  100
 kcal   kJ 100   kJ 1h
Q p = (2m × 2.5m) * 730,33 2 = 3651,65 × 4.1868 = 15288.7 × = 4.24kW
kcal m h h kcal h 3600 s
Q p = 730.33 2
m h
Las pérdidas por todas las paredes son:

(Q )
p total = (10.14 + 8.11 + 5.77 + 4.24) kW = 28.26 KW

Perdidas de calor total:

QT = LH 2 O + ∆Hi + (QP ) total = (12.43 + 29.76 + 28.26) KW = 70.45 KW

 P − QT 
nutil =  T  * 100% = 76.5%
 PT 

1.8) Consumo de gas natural: Utilizando la ecuación de potencia térmica

PT = Q * PCI (11)

300 KW 1KJ / s 3600s

Q= * * = 30.48m 3 / h
MJ 1000 KJ 1KW 1h
35.43 3 *
m MJ

Análisis y conclusiones:

La determinación de las composiciones en los humos secos se

determinaron utilizando la siguiente grafica:
%CO = 0 %O2 = 4% %CO2 = 9%
% N 2 = 87%

%H 2 = 0

n1=1.00 n 2=0.2

Los resultados desarrollados gráficamente son bastante similares con los

obtenidos por medio del algoritmo, por tanto se concluye un correcto
cálculo de los porcentajes utilizando esta herramienta.

SEGUNDO PUNTO

A continuación se presenta la reacción del combustible:

Ecuaciones de Balance:

Reacciones de disociación:
 a)

b)

c)

Constantes de equilibrio asociadas:

a)

b)

c)

Debido a que el proceso es adiabático, entonces el balance de energía

queda:

Debemos definir un sistema de ecuaciones usando las ecuaciones de

balance y las ecuaciones asociadas a las constantes de equilibrio:
Para hallar la temperatura de llama adiabática, cuya expresión se mostro
anteriormente se requieren conocer los coeficientes estequiométricos de las
especies en equilibrio y de los valores de las capacidades caloríficas, por lo
tanto se debe de seguir el siguiente procedimiento:

1. Suponer el valor de la temperatura de llama adiabática T.

2. Con dicha temperatura buscar en tablas los valores de las constantes

de equilibrio para cada reacción de disociación y además buscar el
valor de Cp para cada especie.

3. Resolver el sistema de ecuaciones no lineal planteado anteriormente

y hallar los respectivos valores de los coeficientes estequiométricos.

4. Reemplazar los valores de los coeficientes estequiométricos y de los

Cp en la expresión de la temperatura de llama adiabática y hallar su
valor

5. Comparar el valor de la Temperatura supuesto en el numeral 1 con el

hallado en el numeral 4.

6. Si los valores comparados en el numeral 4 no son iguales, entonces

se vuelve a repetir el procedimiento y, la nueva temperatura que se
va a suponer en el numeral 1 es la hallada en el numeral 4.

TERCER PUNTO

Encontrar las expresiones para calcular las tasa de reacción de las especies
presentes en las reacciones.

C2H6 + 3O2 ↔ 3H2O + 2CO (1)

CO+ 1/2O2 ↔ CO2 (2)

Tenemos que la tasa neta de reacción de una especie dada (k) es la suma
de las tasas de reacción de la especie en cada reacción (j).

(1)

Siendo,

(2)

Donde Wk es el peso molecular.

υkj es es la resta de los coeficientes estequiométricos de la especie k en la
reacción j (υkj=υkj’’-υkj’).
Qj es la tasa de progreso de la reacción j.

Donde [Xk] es la concentración molar de la especie k y está dada por ρk/Wk.

Kfj y Krj son las constante de velocidad de reacción directas e inversas de la

reacción j y se calculan de las siguiente manera.

Con Afj siendo el coeficiente pre exponencial.

T la temperatura.
βj exponente de la temperatura.
Ej la energía de activación
R constante de gases ideales.

De donde se despeja Krj.

Kp se puede obtener a partir de tablas o calcular a partir de la función de

Gibbs que también se obtiene en tablas.

Donde ∆G= (υege+ υfgf - υaga - υbgb) con gi siendo la función de Gibbs para
una temperatura dada y siendo una reacción de la forma υA + υB ↔ υE + υF,
siendo la función de Gibbs expresada en términos de entalpias y entropías
de formación, se obtienen estos valores a partir de tablas termodinámicas.

Para una temperatura de 1000K se tiene que la función de Gibbs para las
especies de de la reacción es:
 h0 s0 h298 h1000
C2 H 6 -84680 229.49
O2 0 205.04 8682 31389
H 2O -241820 188.83 9904 35882
CO -110530 197.65 8432 30355
CO2 -393520 213.18 903 34455

gC2H6=-32890kj/kmol gO2=227747kj/kmol gH2O=-27012kj/kmol

gCO=109043kj/kmol gCO2=-393520kj/kmol

Obteniendo por lo tanto para cada una de las reacciones:

∆G1=-513301 ∆G2=-377244

Calculamos entonces Kp; a partir de los ∆G se obtiene entonces para cada

reacción,

Kp1= 1.5x10-27 Kp2=1.97x10-20

A partir de tablas buscamos los coeficientes para dada una de las

reacciones, teniendo entonces:

Para

C2H6

CO

Ahora teniendo los valores de Kf y Kr para las reacciones (1) y (2) a partir de
Kp y de una temperatura de reacción, hallamos su respectiva tasa de
progreso Qj.

ρC2H6= 1.282kg/m3 ρO2=1.354kg/m3 ρH2O=0.082kg/m3

ρCO=1.184kg/m3 ρCO2=1.87kg/m3

WC2H6=30kg/mol WO2=32kg/mol WH2O=18kg/mol

WCO=28kg/mol W CO2=44kg/mol
[X]C2H6=0.043mol/m3 [X]O2=0.042mol/m3 [X]H2O=0.004mol/m3 [X]CO=
0.042mol/m3 [X] CO2=0.042mol/m3

Para la reacción (1)

Para la reacción 2

Teniendo entonces para cada una de las especies de las reacciones (1) y (2)
la tasa de reacción de las mismas: