Energía 1
OSTWALD BUNTE
Dado un combustible del cual se conoce el Análisis Ultimo, la expresión de la
fórmula para un quilo de combustible seco es la siguiente:
Ca H b Oc Nd Se
12 1 16 14 32
Donde a, b, c, d y e son los kilogramos de cada componente por kilogramo de
combustible (0 < a, b, c, d, e < 1).
A los efectos del desarrollo que sigue se deberá tener en cuenta:
1) Si existe carbono sin quemar C´, el carbono quemado será:
C quemado = C fuel − C '
2) Cuando el combustible contenga S, lo expresaremos como C equivalente:
C equivalente =
e
32
× 12 con [C equivalente ] = kmol
kg
S
×
kg C
kmol C
fuel
e
C total = C equivalente + C fuel = a + × 12
32
Teniendo en cuenta estas dos consideraciones se define el coeficiente a* de la
siguiente manera:
e
a* = a + × 12 − C'
32
Con este valor rescribimos la fórmula para el combustible de la siguiente manera:
Ca* H b Oc Nd
12 1 16 14
donde no consideramos el S como un elemento separado y lo consideramos junto
con el C. Compruebe que el número de moles de C corregido es igual al número de
moles de C más el de S (menos los moles de C no quemados).
El fundamento para el procedimiento anterior de corrección de la fórmula del
combustible se verá en este repartido más adelante.
Si se escribe la ecuación correspondiente a la combustión completa y con exceso,
se verá que de la realización del balance elemental de masas se obtienen 4
ecuaciones y se tienen 6 incógnitas α, β, γ, ε, ν, µ.
IIMPI – Facultad de Ingeniería – UdelaR 1
Energía 1
C a* H b Oc N d + α (O2 + 3,76 N 2 ) = βCO2 + γH 2O + εCO + υO2 + µN 2
12 1 16 14
a*
C) = β +ε
12
b
H) = γ
2
c γ ε
O) +α = β + +υ +
32 2 2
d
N) + 3,76α = µ
28
De este modo el sistema queda indeterminado y por lo tanto con infinitas soluciones
posibles. Con el fin de obtener más datos se recurre a un análisis de gases de
combustión realizados con un aparato Orsat.
Este análisis nos da los siguientes resultados (porcentaje en volumen o molar) sobre
base seca:
CO2%, O2%, CO% y se calcula el porcentaje de N2 por diferencia.
En realidad el valor de CO2% es de CO2 + SO2 debido a las características del
equipo. Esta es la razón por la cuál se corrige el valor de C incluyendo el S.
Expresando estos resultados como fracciones molares del total de componentes
secos tendremos:
X CO 2 = β
%CO 2
= = (CO2 )
β + ε +ν + µ 100
%O 2
=ν
X O2
β + ε + ν + µ 100 = (O2 )
=
=ε
%CO
X CO
β + ε + ν + µ = 100 = (CO )
Entonces se tienen además de las 4 ecuaciones de balance, 3 valores de Orsat:
(CO2), (O2), (CO), un total de 7 ecuaciones y 6 incógnitas por lo cual el sistema es
incompatible o tiene solución. Ya que sabemos que el fenómeno físico ocurre, en
estas condiciones deberá existir una solución. En ese caso una ecuación es
combinación lineal de las otras y por lo tanto cada coeficiente y el término
independiente de una de ellas es CL de los coeficientes correspondientes de las
ecuaciones restantes.
IIMPI – Facultad de Ingeniería – UdelaR 2
Energía 1
Reordenando las 7 ecuaciones:
a*
= β +ε
12
b
=γ
2
c γ ε
= β + +υ + −α
32 2 2
d
= µ − 3,76α
28
0 = [(CO 2 ) − 1]β + (CO 2 )ε + (CO2 )ν + (CO2 ) µ
0 = (O2 ) β + (O2 )ε + [(O2 ) − 1]ν + (O2 ) µ
0 = (CO ) β + [(CO ) − 1]ε + (CO )ν + (CO ) µ
La siguiente tabla muestra los coeficientes que multiplican a cada variable y el
término independiente (columnas) en cada una de las ecuaciones (filas).
α β γ ε ν µ T.I.
0 1 0 1 0 0 a*/12
0 0 1 0 0 0 b/2
-1 1 ½ ½ 1 0 c/32
- 3.76 0 0 0 0 1 d/28
0 (CO2-1) 0 (CO2) (CO2) (CO2) 0
0 (CO) 0 (CO-1) (CO) (CO) 0
0 (O2) 0 (O2) (O2-1) (O2) 0
a b c d
T.I.) 0 = A + B +C + D + E×0+ F ×0+G×0
12 2 32 28
µ ) 0 = A × 0 + B × 0 + C × 0 + D × 1 + E × (CO2 ) + F × (CO ) + G × (O2 )
α ) 0 = A × 0 + B × 0 + C × (−1) + D × (−3,76) + E × 0 + F × 0 + G × 0
Todas las ecuaciones del sistema tienen que verificarse, en particular, la
correspondiente a la variable µ:
µ ) 0 = A × 0 + B × 0 + C × 0 + D × 1 + E × (CO2 ) + F × (CO ) + G × (O2 )
⇒ E × (CO 2 ) + F × (CO) + G × (O 2 ) + D = 0
Esta es una ecuación general de un plano en las coordenadas (CO2), (O2), (CO).
IIMPI – Facultad de Ingeniería – UdelaR 3
Energía 1
Conclusión: Los 3 valores leídos en Orsat (CO2), (O2) y (CO) no son
independientes entre sí.
Es decir, los valores obtenidos de un análisis Orsat de humos en los que únicamente
se encuentre CO como producto de combustión incompleta1 no serán cualesquiera y
su representación en la terna (O2), (CO2), (CO) determinará un plano.
Buscaremos definir de qué plano se trata encontrando los valores de D, E, F y G.
Determinación del plano
Si se toma D ≠ 0 no quedan incluidos aquellos planos que pasan por el origen, lo
cual es correcto ya que se sabe que el conjunto de valores CO2=0, CO=0 y O2=0
correspondientes al punto (0, 0, 0) carece de significado físico en una combustión.
Dividiendo entonces la ecuación general por –D y renombrando obtendremos:
P × (CO 2 ) + Q × (CO) + R × (O 2 ) = 1
Para determinar este plano bastará conocer 3 puntos particulares.
Tomaremos 3 puntos particulares que dependen únicamente del combustible con el
cual se trabaja. Estos son los correspondientes a las siguientes combustiones:
1. Combustión completa sin exceso (estequiométrica).
2. Combustión donde se utiliza todo el O2 para transformar todo el C* en CO (sin
CO2 ni O2 en humos).
3. Combustión con exceso tendiendo a infinito.
1. Combustión completa sin exceso
(CO ) = (O 2 ) = 0 ⇒ (CO 2 ) = (CO 2 )max corresponde al punto (0, CO2max, 0)
Sustituyendo en la ecuación del plano se obtiene
P= 1
(CO 2 )max
2. Combustión donde todo el C* se transforma en CO y (CO2) = (O2) = 0
(CO 2 ) = (O 2 ) = 0 ⇒ (CO) = (CO)max corresponde al punto (0, 0, COmax)
Q= 1
(CO )max
3. Combustión con exceso tendiendo a infinito
1
También puede haber hollín como producto de combustión incompleta, aunque en la ecuación de combustión
planteada aquí está implícito en la corrección del contenido de carbono (a*).
IIMPI – Facultad de Ingeniería – UdelaR 4
Energía 1
Si el aire entrado tiende a infinito, la fracción molar de humos producidos por la
combustión como el CO2, CO y H2O será despreciable con respecto a la del O2 y el
N2, por lo tanto:
(CO) = (CO 2 ) = 0 ⇒ (O 2 ) = (O 2 )max corresponde al punto (O2max, 0, 0)
R= 1
(O 2 )max
Sustituyendo los valores de P, Q y R hallados en la ecuación del plano se obtiene:
(CO 2 ) (CO) (O 2 )
+ + =1
(CO 2 )max (CO)max (O 2 )max
Ecuación de Ostwald-Bunte
Cálculo de valores máximos
β
1) (0, CO2max, 0) ⇒ ν= ε = 0 ⇒ (CO 2 ) =
β+µ
Los máximos los calculamos con el Análisis Ultimo del combustible, considerando su
combustión completa y sin exceso (Estequiométrica).
a * 12
(CO 2 max ) = ⇒
a * 12 + d 28 + 3,76 (a * 12 + b 4 − c 32 )
1
(CO2 max ) = ⇒
12 d 12 b c
4,76 + + × 3,76 −
28 a * 4 a * 8a *
1
(CO2 max ) =
K
4,76 + 0,43K N + 11,3 K H − O
8
d N b H c O
con K N = = ; KH = = ; KN = =
a* C a* C a* C
ε
2) (0, 0, COmax) ⇒ ν = β = 0 ⇒ (CO) =
µ+ε
IIMPI – Facultad de Ingeniería – UdelaR 5
Energía 1
a * 12
(CO max ) = ⇒
a * 12 + d 28 + 3,76(a * 24 + b 4 − c 32 )
1
(COmax ) = ⇒
3,76 12 d b c
1+ + + 3 × 3,76 −
2 28 a * a * 8a *
1
(CO max ) =
K
2,88 + 0,43K N + 11,3 K H − O
8
3) (O2max, 0, 0)
(O2 max ) = 0,208 ( 20.8% de O2 vol/vol en el aire)
Construcción del Diagrama de Ostwald – Bunte:
A continuación se hallan las expresiones de las rectas de corte del plano de
Ostwald-Bunte (OB) con los planos de CO=0, O2=0 y CO2=0.
a) Corte de OB y CO=0
(CO 2 ) (CO) (O 2 )
+ + =1 (CO 2 ) (O2 )
(CO 2 )max (CO)max (O 2 )max + =1
(CO 2 )max (O2 )max
(CO) = 0
(CO 2 max )
(O 2 )
⇒ (CO 2 ) = 1 -
(O 2 )max
Puntos pertenecientes al segmento CO2max_O2max
(O 2 ), 1 - (CO 2 max ),0
(O 2 )
Los puntos entonces serán de la forma: (O )
2 max
b) Corte de OB y O2 = 0
(CO 2 ) (CO) (O 2 )
+ + =1 (CO 2 ) (CO)
(CO 2 )max (CO)max (O 2 )max + =1
(CO 2 )max (CO)max
(O 2 ) = 0
(CO max )
(CO 2 )
⇒ (CO) = 1 -
(CO 2 )max
IIMPI – Facultad de Ingeniería – UdelaR 6
Energía 1
Los puntos entonces serán de la forma: 0, (CO 2 ) , 1 - (CO max )
(CO 2 )
(CO 2 )max
c) Corte de OB y CO2=0
(CO 2 ) (CO) (O 2 )
+ + =1 (O 2 ) (CO)
(CO 2 )max (CO)max (O 2 )max + =1
(O 2 )max (CO)max
(CO 2 ) = 0
(CO max )
(O 2 )
⇒ (CO) = 1 -
(O 2 )max
Los puntos entonces serán de la forma: (O2 ), 0 , 1 - (CO max )
(O 2 )
(O 2 )max
Figura 1 – Triángulo de Ostwald Bunte en 3D
El corte del plano de Ostwald – Bunte con los tres planos ortogonales forma el
Triángulo de Ostwald Bunte. Solamente se considera el cuadrante con (CO) > 0,
(O2) > 0 y (CO2) > 0 porque contiene los valores físicamente posibles para
estas variables.
Utilización del diagrama:
IIMPI – Facultad de Ingeniería – UdelaR 7
Energía 1
La utilidad del diagrama de Ostwald Bunte reside en que es una representación
gráfica del conjunto de todas las combustiones posibles con un combustible dado2, 3.
Es decir, representando los datos de humos de la combustión en el plano OB,
podemos tener una idea rápida de cómo se está produciendo la misma.
En realidad no se trabaja con el plano visto en 3 dimensiones sino que es más
manejable la proyección del plano de OB sobre el plano CO2 – O2 (o de CO = 0).
CO2m ax
(CO2)
CO = 0
E=0
(O2) O2m ax
Figura 2 – Proyección del Triángulo de Ostwald Bunte en el plano O2 – CO2
La figura muestra la proyección de las rectas de (CO) = 0 y de E = 0.
En el plano, normalmente se representan las rectas de exceso cero y de exceso
constante (distinto de cero), las cuales contienen los puntos de las posibles
combustiones con el exceso dado. Con esto podemos saber simplemente
observando la ubicación del punto en el plano, el exceso de la combustión en
cuestión.
También se representan las rectas de CO constante, es decir, las rectas que
contienen las posibles combustiones con un contenido de CO en humos
especificado. Esto da una idea del grado de completitud de la reacción.
A continuación se desarrollan las ecuaciones de las rectas de exceso constante y de
CO constante.
1) Recta de E = 0
2
En realidad, si en la combustión hay carbono no quemado, entonces hay infinitos planos de Ostwald-Bunte para
un mismo combustible (uno para cada valor de C’), ya que al variar el C* cambian los valores de COmax y
CO2max. Sin embargo esto se puede ver como que al tomar diferentes C* en realidad estamos trabajando con
combustibles diferentes, es decir un combustible con una “fórmula” diferente.
3
Estrictamente hablando no contiene todas las combustiones posibles con dicho combustible, sino que las que
únicamente contengan CO como producto de combustión incompleta en la combustión con aire.
IIMPI – Facultad de Ingeniería – UdelaR 8
Energía 1
Utilizaremos la ecuación de exceso en función de los datos de Orsat:
(CO)
(O2 ) −
E= 2
(N 2 ) (CO)
− (O2 ) −
3,76 2
la cual es correcta para combustibles sin N (o con contenido despreciable de este
elemento) y para combustiones donde se tiene únicamente CO como producto de
combustión incompleta.
Igualando a cero esta ecuación se tiene:
(CO)
(O2 ) − = 0 ⇒ (CO) = 2(O2 )
2
Combinando esta ecuación con la ecuación de OB:
(CO 2 ) (CO) (O 2 )
+ + =1 (CO 2 ) 2(O 2 ) (O 2 )
(CO 2 )max (CO)max (O 2 )max + + =1
(CO 2 )max (CO)max (O 2 )max
(CO) = 2(O 2 )
2 1
⇒ (CO 2 ) = (CO 2 )max 1 − + (O 2 )
(CO)max (O 2 )max
Esta expresión, es la ecuación de una recta en el plano CO2 – O2. Para O2 = 0 ⇒
CO2 = CO2max, es decir que se verifica que el punto de CO2max corresponde al de
combustión completa sin exceso (CO = 0 y O2 = 0).
2) Rectas de E = cte
Combinando las ecuaciones de exceso en función de datos de Orsat y OB y luego
de algunas operaciones se obtiene la siguiente expresión:
E (CO)max 4,76 (CO)max E
− λ (CO 2 ) + 1 + E − λ (O 2 ) − + (CO)max λ = 0
3,76 (CO 2 )max 3,76 (CO 2 )max 3,76
1 1 1
con λ = E - −
3,76 2 2
la cual es la ecuación de una familia de rectas parametrizadas en E. Asignando a E
un valor se obtiene la recta correspondiente a ese valor de exceso.
Se puede verificar que las rectas de exceso constante pasan por el punto: (O2) = 1
(100%), (CO2) = -2 (-200%) independientemente del combustible. Es decir,
cualquiera sea el combustible, las rectas de exceso constante pasan siempre
IIMPI – Facultad de Ingeniería – UdelaR 9
Energía 1
por el punto (1, -2). Evidentemente este punto no corresponde a los datos de una
combustión real.
0,200
0,000
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1
-0,200
-0,400
-0,600
(CO2)
-0,800
-1,000
-1,200
-1,400
-1,600
-1,800
-2,000
(O2)
Figura 3 – Familia de rectas de exceso constante en plano (O2) – (CO2)
Algunas rectas desde exceso = 0 hasta exceso = 2 (200%) para un combustible
de composición: C = 86%, H = 12%, S = 2%.
En realidad, a efectos prácticos, las rectas de exceso constante son prácticamente
paralelas a la de exceso cero, sobre todo las de exceso cercano a cero. Esto se
puede ver en la figura siguiente.
0,20
CO=0
E=0
E = 0,2
(CO2)
E = 0,5
0,10
E = 1,0
E = 2,0
0,00
0 0,05 0,1 0,15 0,2
(O2)
Figura 4 – Familia de rectas de exceso constante en triángulo de OB
Algunas rectas desde exceso constante para un combustible de composición:
C = 86%, H = 12%, S = 2%.
La recta de exceso cero divide el triángulo en dos zonas: oxidante (E > 0) y
reductora (E < 0).
IIMPI – Facultad de Ingeniería – UdelaR 10
Energía 1
La zona oxidante, o de exceso positivo, se encuentra limitada por la recta de exceso
cero y la de (CO) cero. En esta zona caen los puntos de combustiones en que el aire
utilizado es mayor al estequiométrico. Cuanto más lejos de la recta de exceso cero
(hacia los O2 mayores), mayor el exceso.
De la misma manera, la zona reductora o de defecto de aire, es la que se encuentra
entre la recta de (O2) cero y la de exceso cero. Cuanto más lejos (hacia los O2
menores) de la recta de exceso cero, menor el exceso.
3) Rectas de CO = cte
(CO 2 ) (CO) (O 2 )
+ + =1 (CO 2 ) k (O 2 )
(CO 2 )max (CO)max (O 2 )max + + =1
(CO 2 )max (CO)max (O 2 )max
(CO) = k
(O 2 ) k
⇒ (CO 2 ) = (CO 2 )max 1 − +
(O 2 )max (CO)max
Esta es la ecuación de una familia de rectas parametrizada en k. Una característica
es que la pendiente de la curva no depende de k, es decir, todas las curvas de CO
constante son paralelas entre si.
En particular, para (CO) = 0 tenemos:
(O 2 )
(CO 2 ) = (CO 2 )max 1 −
(O 2 )max
Como ya habíamos visto en la sección anterior.
IIMPI – Facultad de Ingeniería – UdelaR 11
Energía 1
0,20
CO=0
CO = 0,02
CO = 0,05
(CO2)
CO = 0,08
0,10
0,00
0 0,05 0,1 0,15 0,2
(O2)
Figura 5 – Familia de rectas de (CO) constante en triángulo de OB
Algunas rectas desde (CO) constante para un combustible de composición:
C = 86%, H = 12%, S = 2%.
Es evidente que a medida que nos alejamos de la recta de (CO) = 0, aumenta la
cantidad de (CO) en humos.
IIMPI – Facultad de Ingeniería – UdelaR 12