Matemática Superior aplicada
Sistemas de Ecuaciones No Lineales
Profesor: Dr. Alejandro S. M. Santa Cruz
JTP: Dr. Juan Ignacio Manassaldi
Aux. 1ra: Ing. Amalia Rueda
Ejemplo: Equilibrio líquido vapor no ideal
• El equilibrio de una mezcla binaria ideal se
representa mediante la Ley de Raoult:
P = Pv A P = p A + pB
T
A A+ B
x A Pv A = p A x A Pv A = y A P ¡No se cumple para
xB PvB = pB xB PvB = yB P
mezclas no ideales!
Vapor gas ideal
Ejemplo: Equilibrio líquido vapor no ideal
• Como primera aproximación a la no idealidad,
se utiliza un modelo de actividad.
• A bajas presiones podemos asumir que el vapor
se comporta de manera ideal.
• La condición de equilibrio líquido-vapor a partir
del coeficiente de actividad corresponde a:
A x A PvA = y A P
B xB PvB = yB P
Ejemplo: Equilibrio líquido vapor no ideal
• Comenzamos a armar el sistema de ecuaciones
A x A PvA = y A P
B xB PvB = yB P
• ¿Qué otra ecuación conocen?
x A + xB = 1 → xB = 1 − x A
yB = 1 − y A
• Reemplazamos
A x A PvA = y A P
B (1 − x A ) PvB = (1 − y A ) P
Ejemplo: Equilibrio líquido vapor no ideal
• Planteamos el problema del calculo de la
presión de rocío (dew point pressure).
• Datos: y A T
A x A PvA = y A P 2 ecuaciones
4 incógnitas
B (1 − x A ) PvB = (1 − y A ) P x P
A B A
• ¿Se puede resolver?
• No se puede resolver, se deben agregar
ecuaciones para representar al coeficiente de
actividad.
Ejemplo: Equilibrio líquido vapor no ideal
• Uno de los modelos mas sencillos corresponde
al de Margules con dos sufijos:
ln ( A ) = AxB 2
ln ( B ) = Ax A 2
• Formamos el sistema final:
A x A Pv A = y A P
B (1 − x A ) PvB = (1 − y A ) P
ln ( A ) = A (1 − x A )
2
ln = Ax 2
( B) A
Ejemplo: Equilibrio líquido vapor no ideal
• Sistema final: A x A PvA − y A P = 0
B (1 − x A ) PvB − (1 − y A ) P = 0
ln ( A ) − A (1 − x A ) = 0
2
ln ( B ) − Ax A = 0
2
A x A Pv A − y A P xA
B (1 − x A ) PvB − (1 − y A ) P
P
f ( x) = x=
ln ( A ) − A (1 − x A )
2
A
ln ( ) − Ax 2
B A
B
Ejemplo: Equilibrio líquido vapor no ideal
Caso de estudio 1:
• Mezcla de Agua (A) y Metanol (B)
• Datos: y A = 0.4
T = 365 K
Pv 365
A = 0.75261bar
PvB365 = 2.705 bar
A = 0.6833
Ejemplo: Equilibrio líquido vapor no ideal
• Reemplazamos en el sistema: y A = 0.4
T = 365 K
Pv365
A = 0.75261bar
PvB365 = 2.705 bar
A = 0.6833
A x A 0.75261 − 0.4 P xA
B (1 − x A ) 2.705 − 0.6 P
P
f ( x) = 2 x=
ln ( A ) − 0.6833 (1 − x A)
A
ln ( ) − 0.6833 x 2
B A
B
Ejemplo: Equilibrio líquido vapor no ideal
• Resolución:
A x A 0.75261 − 0.4 P xA
B (1 − x A ) 2.705 − 0.6 P
P
f ( x) = 2 x=
ln ( A ) − 0.6833 (1 − x A )
A
ln ( ) − 0.6833 x 2
B A
B
A 0.75261 −0.4 x A 0.75261 0
− B 2.705 −0.6 0 (1 − xA ) 2.705
1
J ( x ) = 2 0.6833 (1 − x A ) 0 0
A
1
−2 0.6833 x A 0 0
B
Ejemplo: Equilibrio líquido vapor no ideal
• ¿Valor de arranque?
• ¡Ley de Raoult!
x A Pv A = y A P x A 0.75261 = 0.4 P
xB PvB = yB P (1 − x A ) 2.705 = 0.6 P
x A 0.75261 + (1 − x A ) 2.705 = P
x A 0.75261
x A 0.75261 + (1 − x A ) 2.705 =
0.4
x A 0.75261 + (1 − x A ) 2.705 = x A1.881525
x A = 0.705545 → P = 1.3275 bar
Ejemplo: Equilibrio líquido vapor no ideal
x A = 0.705545 → P = 1.3275 bar
0.705545
1.3275
=
(0)
x
1
1
0.7870393 x −x
(1) (0)
= 0.4670967
1.515948
x
(1)
=x
(0)
( ) ( )
−J x
(0) −1
F x
(0)
=
1.0264512 ( ) = 0.1181059
f x
(1)
0.7779367 1.4187191 0.7779346
1.5138665 1.5138641
= =
(3) (4)
x x
1.0342683 1.0342697
1.5121278 1.5121326
f x ( ) = 0.0000143
(3)
f x ( ) = 2.946 10
(3) -11
Ejemplo: Equilibrio líquido vapor no ideal
Caso de estudio 2:
• Mezcla de Agua (A) y Metanol (B)
• Datos: y A = 0.5259
T = 322.91K
Pv322.91
A = 0.122232 bar
PvB322.91 = 0.548176 bar
A = 0.6833
• Valor experimental:
x A* = 0.8782
P = 0.20932 bar
*
Ejemplo: Equilibrio líquido vapor no ideal
• Ley de Raoult: x A = 0.832628 → P = 0.193523 bar
0.832628 0.910879 0.895141
0.193523 0.211951 0.209622
x = →x = =
(0) (4)
(1)
... → x
1 1.001243 1.0075414
1 1.562750 1.728948
0.895141 ( ) = 3.114 10
f x
(4) -10
0.209622
x = Solución
*
1.0075414
1.728948
x A* = 0.8782 Valor experimental
ErrorxA = 1.93%
P* = 0.20932 bar ErrorP = 0.14 %
Ejemplo: Equilibrio líquido vapor no ideal
• Pantear el sistema utilizando el modelo de Van
Laar:
−2
A xA
ln ( A ) = A 1 +
B xB
−2
B xB
ln ( B ) = B 1 +
A xA
Ejemplo: Equilibrio líquido vapor no ideal
Caso de estudio 3:
• Mezcla de Agua (A) y Metanol (B)
• Datos: y A = 0.4
T = 365 K
Pv365
A = 0.75261bar
Pv 365
B = 2.705 bar
A = 0.2439
B = 0.3861