Predicción del equilibrio líquido-vapor a partir del equilibrio líquido-líquido
Enviado por Andres David Aguirre Cudris • 6 de Marzo de 2019 • Tarea • 601 Palabras (3 Páginas) • 183 Visitas
PREDICCIÓN DEL EQUILIBRIO LÍQUIDO VAPOR A PARTIR DE EQUILIBRIO LÍQUIDO – LÍQUIDO
Basados en los datos experimentales presentados en el artículo “Vapor-Liquid-Liquid Equilibrium (VLLE) Data for the Systems Ethyl acetate + Water, Toluene + Water and Toluene + Ethyl acetate + Water at 101.3 kPa. Using Modified Equilibrium Still”, para el sistema binario Acetato de etilo – agua, se predijeron las composiciones del equilibrio líquido vapor del sistema.
Se sabe que en el equilibrio líquido – líquido debe cumplirse la igualdad de fugacidades para cada componente en cada fase.
[pic 1]
[pic 2]
En este caso con i=1,2; por tratarse de un sistema binario.
Inicialmente se hace necesaria la elección de un método de coeficientes de actividad. En este caso se utilizó el método UNIQUAC, el cual se planteó de la siguiente manera.
[pic 3]
[pic 4]
[pic 5]
Con:
COMPONENTE | r | q | q' |
ACETATO DE ETILO (1) | 3,48 | 3,12 | 3,12 |
AGUA(2) | 0,92 | 1,40 | 1,00 |
Los datos de equilibrio líquido - líquido se presentan a continuación:
P/atm | T/°C | P1sat/atm | P2sat/atm | X1' | X2' | X1'' | X2'' |
1 | 24,80 | 0,12 | 0,03 | 0,9500 | 0,0500 | 0,0187 | 0,9813 |
1 | 30,77 | 0,16 | 0,04 | 0,9500 | 0,0500 | 0,0176 | 0,9824 |
1 | 34,27 | 0,19 | 0,05 | 0,9357 | 0,0643 | 0,0172 | 0,9828 |
1 | 41,61 | 0,27 | 0,08 | 0,9165 | 0,0835 | 0,0168 | 0,9832 |
1 | 49,68 | 0,37 | 0,12 | 0,8849 | 0,1151 | 0,0170 | 0,9830 |
1 | 57,90 | 0,51 | 0,18 | 0,8291 | 0,1709 | 0,0159 | 0,9841 |
1 | 65,33 | 0,67 | 0,25 | 0,8211 | 0,1789 | 0,0172 | 0,9828 |
Según la ecuación
,[pic 6]
al ser para una temperatura dada, entonces se debe cumplir la igualdad:[pic 7]
[pic 8]
Sacando provecho de esto, se pudieron hallar los parámetros del modelo de actividad escogido, para así poder predecir más tarde los datos del ELV.
Para calcular los parámetros se planteó la anterior igualdad, y el modelo para los coeficientes de actividad desarrollado en términos de sus parámetros y la composición de la fase correspondiente. Se pusieron a variar los parámetros de manera que se cumpliera la igualdad, buscando que:
[pic 9]
Después de realizado este procedimiento se obtuvieron los siguientes parámetros:
u12 [cal/mol] | u21 [cal/mol] |
1736,058 | 94,6590 |
Contando ya con los parámetros para el modelo de coeficientes de actividad, se calcularon los mismos para cada punto y se comprobó la igualdad presentada anteriormente, para cada componente. Estos fueron los resultados.
T/°C | X1' | X2' | X1'' | X2'' | Gamma1' | Gamma1'' | Gamma2' | Gamma2'' |
24,80 | 0,9500 | 0,0500 | 0,0187 | 0,9813 | 1,0120 | 69,0050 | 16,8788 | 1,0036 |
30,77 | 0,9500 | 0,0500 | 0,0176 | 0,9824 | 1,0114 | 69,2321 | 16,2773 | 1,0032 |
34,27 | 0,9357 | 0,0643 | 0,0172 | 0,9828 | 1,0177 | 69,0362 | 14,3144 | 1,0031 |
41,61 | 0,9165 | 0,0835 | 0,0168 | 0,9832 | 1,0272 | 68,0447 | 12,1216 | 1,0029 |
49,68 | 0,8849 | 0,1151 | 0,0170 | 0,9830 | 1,0466 | 66,1481 | 9,7452 | 1,0030 |
57,90 | 0,8291 | 0,1709 | 0,0159 | 0,9841 | 1,0922 | 65,9551 | 7,2309 | 1,0026 |
65,33 | 0,8211 | 0,1789 | 0,0172 | 0,9828 | 1,0969 | 62,9019 | 6,8583 | 1,0030 |
...