Equilibrio Quimico Olikara Borman

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

UNIVERSIDAD NACIONAL DEL COMAHUE

FACULTAD DE INGENIERIA

M
aquinas Motrices I

Estrategia de Olikara y Borman para la


determinaci
on de la composici
on de equilibrio de los
productos de la combusti
on 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 combusti
on 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
Cn Hm Oz Nw +

nO 2
( O2 + 3.773N2 ) =nP (
x1 CO2 + x
2 H2 O+
x3 N2 + x
4 O2 + x
5 CO

+x
6 H2 + x
7 H+
x8 O+
x9 OH+
x10 NO+
x11 N

(1)

donde nO2 es el n
umero de moleculas de O2 requeridas para completar la combustion (n + m/4
z/2), nP es el n
umero total de moles de los productos, y x
i es la fraccion molar de la i-esima
especie en los productos.
Pueden plantearse los siguintes balances atomicos:
Balance de carbono:
n = nP (
x1 + x
5 )

(2)

m = nP (2
x2 + 2
x6 + x
7 + x
9 )

(3)

Balance de hidr
ogeno:

Balance de oxgeno:
z+

2nO2
= nP (2
x1 + x
2 + 2
x4 + x
5 + x
8 + x
9 + x
10 )

(4)

7.546nO2
= nP (2
x3 + x
10 + x
11 )

(5)

Balance de nitr
ogeno:
w+

Una ecuaci
on adicional la aporta la restriccion de que la suma de todas las fracciones molares
debe ser igual a la unidad
11
X
x
i = 1
(6)
i=1

Conocidas la composici
on del combustible (n, m, z y w), la relacion de equivalencia, la presi
on
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 restricci
on (6). Por lo tanto, para resolver completamente el problema, deben
plantearse 7 ecuaciones adicionales apropiadas.
1

UNIVERSIDAD NACIONAL DEL COMAHUE


FACULTAD DE INGENIERIA

M
aquinas Motrices I

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


7 ecuaciones no redundantes

i
1
2
3
4
5
6
7

Ai
0.432168
0.310805
0.389716
-0.141784
0.150879101
-0.752364
-0.415302102

1
H2 = H
2

1
O2 = O
2

1
N2 = N
2

1
1
H2 + O2 = OH
2
2

1
1
O2 + N2 = NO
2
2

1
H2 + O2 = H2 O
2

Kp6 =

1
CO + O2 = CO2
2

Kp7 =

Bi
-0.112464105
-0.129540105
-0.245828102
-0.213308104
-0.470959104
0.124210105
0.148627105

Kp1 =
Kp2 =
Kp3 =
Kp4 =
Kp5 =

Ci
0.267269101
0.321779101
0.314505101
0.853461
0.646096
-0.260286101
-0.475746101

x
7 p1/2
1/2

x
6

x
8 p1/2
1/2

x
4

x
11 p1/2
1/2

x
3
x
9

1/2 1/2

(7)

x
4 x
6
x
10

1/2 1/2

x
4 x
3
x
2
1/2

6 p1/2
x
4 x
x
1
1/2

x
5 x
4 p1/2
Di
-0.745744104
-0.738336104
-0.963730101
0.355015104
0.272805101
0.259556103
0.124699103

Ei
0.242484108
0.344645102
0.585643102
-0.310227102
-0.154444102
-0.162687101
-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 individuales consideradas fueron obtenidos de las Tablas Termoqumicas JANAF.2 Las constantes
de equilibrio a presi
on constante para cada una de las siete reacciones planteadas en (7) fueron
obtenidas como una funci
on 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)
reactantes

productos

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 reacci
on 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 composici
on de los productos sin disociar. Una vez que se determina la composici
on
(es decir, que el metodo iterativo converge), pueden calcularse las propiedades termodinamicas.
2

UNIVERSIDAD NACIONAL DEL COMAHUE


FACULTAD DE INGENIERIA

M
aquinas Motrices I

Referencias
[1] Olikara C, Borman GL. A computer program for calculating properties of equilibrium combustion 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