Equilibrio - Quimico Olikara Borman PDF

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 3

UNIVERSIDAD NACIONAL DEL COMAHUE

FACULTAD DE INGENIERIA Maquinas Motrices I

Estrategia de Olikara y Borman para la


determinacion de la composicion de equilibrio de los
productos de la combustion en MCI

Consideremos un combustible general que puede representarse por la formula Cn Hm Oz Nw .


Las especies de mayor importancia que se encuentran en los productos de la combustion en
motores de combustion interna sobre su rango normal de operacion son: CO2 , H2 O, CO, H2 ,
O2 , N2 , NO, OH, O, H y N. Entonces, la reaccion de combustion de este combustible generico
con aire puede escribirse como
nO 2
Cn Hm Oz Nw + ( O2 + 3.773N2 ) =nP (x1 CO2 + x2 H2 O+x3 N2 + x4 O2 + x5 CO
(1)
+ x6 H2 + x7 H+x8 O+x9 OH+x10 NO+x11 N

donde nO2 es el numero de moleculas de O2 requeridas para completar la combustion (n + m/4


z/2), nP es el numero total de moles de los productos, y xi es la fraccion molar de la i-esima
especie en los productos.
Pueden plantearse los siguintes balances atomicos:

Balance de carbono:
n = nP (x1 + x5 ) (2)

Balance de hidrogeno:
m = nP (2x2 + 2x6 + x7 + x9 ) (3)

Balance de oxgeno:
2nO2
z+ = nP (2x1 + x2 + 2x4 + x5 + x8 + x9 + x10 ) (4)

Balance de nitrogeno:
7.546nO2
w+ = nP (2x3 + x10 + x11 ) (5)

Una ecuacion adicional la aporta la restriccion de que la suma de todas las fracciones molares
debe ser igual a la unidad
X11
xi = 1 (6)
i=1

Conocidas la composicion del combustible (n, m, z y w), la relacion de equivalencia, la presion


total de la mezcla (p) y la temperatura (T ), la reaccion de combustion (1) posee 12 incognitas: nP
y las concentraciones de las 11 especies producto consideradas. Tenemos disponibles 4 balances
atomicos y la restriccion (6). Por lo tanto, para resolver completamente el problema, deben
plantearse 7 ecuaciones adicionales apropiadas.

1
UNIVERSIDAD NACIONAL DEL COMAHUE
FACULTAD DE INGENIERIA Maquinas Motrices I

Olikara y Borman1 propusieron un metodo en el cual se consideran en equilibrio las siguientes


7 ecuaciones no redundantes

1 x7 p1/2
H2 = H Kp1 = 1/2
2 x6
1 x8 p1/2
O2 = O Kp2 = 1/2
2 x4
1 x11 p1/2
N2 = N Kp3 = 1/2
2 x3
1 1 x9
H2 + O2 = OH Kp4 = (7)
2 2 1/2 1/2
x4 x6
1 1 x10
O2 + N2 = NO Kp5 = 1/2 1/2
2 2 x4 x3
1 x2
H2 + O2 = H2 O Kp6 = 1/2
2 x4 x6 p1/2
1 x1
CO + O2 = CO2 Kp7 = 1/2
2 x5 x4 p1/2

i Ai Bi Ci Di Ei
1 0.432168 -0.112464105 0.267269101 -0.745744104 0.242484108
2 0.310805 -0.129540105 0.321779101 -0.738336104 0.344645102
3 0.389716 -0.245828102 0.314505101 -0.963730101 0.585643102
4 -0.141784 -0.213308104 0.853461 0.355015104 -0.310227102
5 0.150879101 -0.470959104 0.646096 0.272805101 -0.154444102
6 -0.752364 0.124210105 -0.260286101 0.259556103 -0.162687101
7 -0.415302102 0.148627105 -0.475746101 0.124699103 -0.900227108

Cuadro 1: Coeficientes para la aproximacion de las constantes de equilibrio a presion constante.


Los datos del equilibrio qumico a presion constante para cada una de las once especies indi-
viduales consideradas fueron obtenidos de las Tablas Termoqumicas JANAF.2 Las constantes
de equilibrio a presion constante para cada una de las siete reacciones planteadas en (7) fueron
obtenidas como una funcion de la temperatura utilizando los datos de las Tablas JANAF:
X X
log Kp (T )|reaccion = ni log Kp (T )|JANAF ni log Kp (T )|JANAF (8)
productos reactantes

donde Kp (T ) es la constante de equilibrio a presion constante, p es la presion y ni es el coeficiente


molar de la especie i en la reaccion en el equilibrio. Luego, las constantes de equilibrio para cada
una de las 7 reacciones en (7) obtenidas con la ecuacion (8) fueron fiteadas sobre el rango de
temperaturas de 600 a 4000 K usando la relacion funcional
B
log Kp (T )|reaccion = A ln TA + + C + DTA + ETA2 (9)
TA
siendo TA = T /1000, con T en K. Los coeficientes para las siete reacciones (7) aproximadas por
la ecuacion (9) se listan en el cuadro 1.
El sistema de ecuaciones resultante (2)-(7) se resuelve empleando el metodo iterativo de
Newton-Raphson. Las fracciones molares iniciales para comenzar las iteraciones se estiman a
partir de la composicion de los productos sin disociar. Una vez que se determina la composicion
(es decir, que el metodo iterativo converge), pueden calcularse las propiedades termodinamicas.

2
UNIVERSIDAD NACIONAL DEL COMAHUE
FACULTAD DE INGENIERIA Maquinas Motrices I

Referencias
[1] Olikara C, Borman GL. A computer program for calculating properties of equilibrium com-
bustion products with some applications to I.C. engines. SAE paper 750468, 1975.

[2] JANAF Thermochemical Tables, 2d ed., NSRDS-NB537, U.S. National Bureau of Standards,
June 1971.

También podría gustarte