Estimación del equilibrio líquido-vapor del sistema binario acetona-cloroformo mediante el modelo termodinámico de Van Laar y Peng Robinson

María Isabel Sandoval Martínez, Samuel Fernando Muñoz Navarro, Humberto Jose Martinez Jiménez

Estimación del equilibrio líquido-vapor del sistema binario acetona-cloroformo mediante el modelo termodinámico de Van Laar y Peng Robinson

Revista ION, vol. 33, núm. 2, 2020

Universidad Industrial de Santander

María Isabel Sandoval Martínez

Universidad Industrial de Santander (UIS), Colombia


Samuel Fernando Muñoz Navarro

Universidad Industrial de Santander (UIS), Colombia


Humberto Jose Martinez Jiménez

Universidad Industrial de Santander (UIS), Colombia


Recibido: 17 Octubre 2019

Aceptado: 25 Mayo 2020

Resumen: La determinación del equilibrio líquido-vapor de diferentes sistemas es información muy importante para la caracterización de las mezclas y los procesos de separación de estas. El equilibrio es estimado de forma experimental; sin embargo, su ejecución tiene altos costos por la precisión que se demanda en dichas pruebas. Una solución viable a la problemática planteada es el uso modelos termodinámicos que permitan obtener el equilibrio de fases. El presente artículo muestra la metodología para la estimación del equilibrio líquido-vapor de un sistema binario no ideal azeotropo, tomando como ejemplo la mezcla de acetona y cloroformo, en la herramienta software Matlab 9.0. Para dar cumplimiento al principal objetivo de la investigación, se elaboró la curva de rocío y burbuja usando la ecuación cúbica de Peng Robinson para calcular los coeficientes de fugacidad, y el modelo termodinámico de Van Laar para los coeficientes de actividad. Las curvas resultantes fueron comparadas con las curvas obtenidas experimentalmente a una temperatura de 50 °C y 35,17 °C, mostrando un buen ajuste, con errores relativos promedios inferiores a 3,9%. Adicionalmente, los resultados fueron comparados con los estimados a través del simulador de procesos químicos Aspen HYSYS y con un programa comercial (VLE Cal) alcanzando porcentajes de error promedio inferiores al 1%.

Palabras clave: Equilibrio Líquido-Vapor, Acetona-Cloroformo, Modelo Van Laar, Coeficiente de Actividad, Coeficiente de Fugacidad.

Abstract: The Determination of the liquid-vapor equilibrium of different systems is very important information for the mixtures characterization and its separation processes. Equilibrium is estimated experimentally; however, its execution has high costs due to the precision demanded in said tests. A good solution to this problem is to employ thermodynamic models that allow obtaining phase equilibrium. This article shows the methodology for estimating the liquid-vapor balance of a non-ideal azeotrope binary system, taking as an example the mixture of acetone and chloroform, in the Matlab 9.0 software tool. To reach the main objective of the research, the dew and bubble curve was developed, using the cubic equation of Peng Robinson to calculate the fugacity coefficients, and the thermodynamic model of Van Laar for the activity coefficients. The resulting curves were compared with the curves obtained experimentally at a temperature of 50 °C and 35.17 °C, showing a well fit, with average relative errors of less than 3.9%. Additionally, the results were compared with those estimated through the Aspen HYSYS chemical process simulator and with a commercial program (VLE Cal), reaching average error percentages of less than 1%.

Keywords: Liquid-Gas Equilibrium, Acetone-Chloroforme, Van Laar Model, Fugacity Coefficient, Activity Coefficient.

Resumo: A determinação do equilíbrio líquido-vapor de diferentes sistemas é uma informação muito importante para a caracterização de misturas e seus processos de separação. O equilíbrio é estimado experimentalmente; no entanto, sua execução apresenta altos custos devido à precisão exigida nos referidos testes. Uma solução viável para o problema proposto é o uso de modelos termodinâmicos que permitem obter o equilíbrio de fases. Este artigo mostra a metodologia para estimar o equilíbrio líquido-vapor de um sistema binário azeotrópico não ideal, tomando como exemplo a mistura de acetona e clorofórmio, na ferramenta de software Matlab 9.0. Para cumprir o objetivo principal da pesquisa, foi desenvolvida a curva de orvalho e bolha, usando a equação cúbica de Peng Robinson para calcular os coeficientes de fugacidade e o modelo termodinâmico de Van Laar para os coeficientes de atividade. As curvas resultantes foram comparadas com as curvas obtidas experimentalmente a uma temperatura de 50 °C e 35,17 °C, mostrando um bom ajuste, com erros relativos médios inferiores a 3,9%. Além disso, os resultados foram comparados com os estimados através do simulador de processo químico Aspen HYSYS e com um programa comercial (VLE cal), atingindo percentuais médios de erro inferiores a 1%.

Palavras-chave: Equilíbrio Líquido-Vapor, Acetona-Clorofórmio, Modelo de Van Laar, Coeficiente de Atividade, Coeficiente de Fugacidade.

La estimación de las condiciones de equilibrio líquido-vapor en mezclas binarias es un tema de interés general para la comprensión del comportamiento de los componentes de una mezcla en equilibrio, que permite la elaboración de diagramas del sistema y que dará las bases para el diseño, análisis y control del proceso de un proceso de destilación, mediante el cual se busca separar las mezclas en componentes puros [1, 2]. El equilibrio líquido-vapor de diferentes mezclas se obtiene comúnmente con pruebas experimentales con el método de destilación, dinámico, de adsorción, entre otros; desarrollados bajo ciertas condiciones de temperatura, presión y concentración [3]. Sin embargo, este tipo de información no puede ser obtenida en todos los casos de forma experimental, ya que los métodos empleados demandan un alto costo de adquisición, relacionados con la necesidad de garantizar puntos con la máxima precisión [2]. Por lo tanto, ha sido necesario elaborar programas computacionales usando ciertos modelos matemáticos para predecir el comportamiento de las mezclas [4,5,6]. El alcance de la presente investigación abarca la programación de modelos termodinámicos y ecuaciones de estado para estimar el equilibrio líquido-vapor en un sistema binario conformado por cloroformo y acetona. Esta mezcla tiene, principalmente, aplicabilidad para la determinación de impurezas mécanicas a muestras de propóleo en laboratorio. Se caracteriza por tener comportamiento no ideal azeótropo a una temperatura de 50 °C [7].

Para ver el equilibrio líquido-vapor de un sistema binario, se construyen curvas P-X-Y, que son gráficas de presión versus la composición de un componente de la mezcla en fase gaseosa y fase líquida. En la elaboración de estas curvas es necesario determinar los coeficientes de fugacidad y actividad. Los coeficientes de fugacidad se pueden obtener mediante ecuaciones de estado [8], tales como: la ecuación Soave-Redlich-Kwong, la cual permite por su simplicidad la determinación del equilibro líquido- vapor (ELV) a altas presiones [9] o la ecuación de Peng Robinson, que fue usada para el desarrollo de este trabajo, ya que predice el equilibrio líquido-vapor con mucha precisión y aproximación a valores experimentales reales, tanto en sistemas binarios como en multicomponentes.

En el cálculo de los coeficientes de actividad se usa la ecuación de Wilson [10], que se basa en el concepto de composición y se caracteriza por ser muy precisa para representar el equilibrio líquido-vapor de sistemas alcohol- hidrocarburo [11] o el modelo de Van Laar, el cual es muy usado por su flexibilidad, simplicidad y capacidad de representar muchos sistemas [12]. Este último fue empleado para el desarrollo de las curvas P-X-Y generadas en este trabajo, puesto que esta ecuación se ha ajustado a los datos de equilibrio líquido - vapor de un sistema binario acetona -cloroformo a una temperatura de 35,17 °C, en trabajos anteriores.

Los resultados de las curvas P-X-Y para un sistema binario cloroformo-acetona, construidas en el presente trabajo a través de la herramienta de software Matlab 9.0, fueron comparadas y mostraron un buen ajuste con errores relativos promedios inferiores a 3,9%, lo que indica la precisión de los modelos termodinámicos y ecuaciones de estado seleccionadas; sin embargo, se recomienda trabajar el mismo sistema con otros modelos para tratar de disminuir el porcentaje de error relativo a valores inferiores al 1%.

Generalidades termodinámicas

Sistemas no-ideales

En la termodinámica existen tres tipos de sistemas: ideales, semi-ideales y no ideales. Los sistemas no ideales se apartan fuertemente en su comportamiento de los modelos que describen los modelos simples y solo se pueden representar mediante métodos de estimación de propiedades termodinámicas que toman en cuenta la interacciones entre partículas. El comportamiento no ideal de un sistema se debe principalmente a su estructura química, la presencia de grupos fuertemente polares en sus moléculas, el gran tamaño de estas, geometría y la posibilidad de ionizarse o asociarse entre sí [12].

En el presente artículo se trabajó con una mezcla no ideal de acetona-cloroformo, la cual tiene grupos polares donde el cloroformo muestra un fuerte desbalance, por el exceso de átomos electronegativos disponibles para formar un puente de hidrógeno; este se debe a que presenta en su composición 3 átomos de cloro frente a un átomo de hidrogeno. Por otro lado, la acetona cuenta con un oxígeno frente a seis hidrógenos [13]. Esta mezcla también es considerada como una disolución azeotrópica, puesto que sus componentes unidos tienen el mismo punto de ebullición, causando que su equilibrio líquido-vapor sea más complejo y por lo tanto, la separación de sus compuestos puros sea imposible de alcanzar mediante método de destilación simple [14,15].

Equilibrio líquido-vapor

Es una condición estática donde no se producen cambios en las propiedades macroscópicas con el tiempo y dinámica a nivel microscópico, ya que las moléculas de la zona líquido están en constante movimiento hacia la zona de vapor y las moléculas de la zona de vapor se mueven hacia la zona de líquido [16].

Para predecir la información respecto al equilibrio líquido-vapor es necesario usar la regla de las fases. Esta permite calcular las fases, grados de libertad o componentes de un sistema para que coexista equilibrio y es expresada mediante la ecuación 1 [17].

(1)

Donde

F=Número de grados de libertad

C=Número de componentes del sistema

P=Número de fases en el sistema

En este trabajo se tuvieron dos grados de libertad, debido a que a que se trabajó con un sistema binario y las fases a estudiar fueron dos (líquido y vapor). Esto indica que se fijaron dos variables para definir completamente el estado de equilibrio. Por lo tanto, en el trabajo se estableció la temperatura y composición para determinar la presión del sistema.

Fugacidad y coeficiente de fugacidad

La fugacidad es una propiedad auxiliar estudiada para describir la presión de un gas real. Cuando se hace referencia a un gas ideal se trabaja con presiones parciales; no obstante, para gases reales se trabaja con el concepto de fugacidad, el cual está relacionado directamente con la tendencia de un sistema a preferir una fase frente a otra [18,19,20]. La fugacidad tiene unidades de presión y depende de esta misma.

El coeficiente de fugacidad es la relación que existe entre la fugacidad y la presión del componente como gas ideal [21]. Este se considera como la corrección del estado ideal y sus valores oscilan entre 0 y 1 [22].

Los modelos más conocidos para la determinación de los coeficientes de fugacidad son: Ecuación Virial, ecuación de Van Der Waals, Readlich Kwong, Readlich Kwong- Soave y Peng Robinson [22,23].

Para este trabajo se empleó la ecuación cúbica de estado de Peng Robinson junto con la regla de mezcla cuadrática, ya que a pesar de haber sido postulada hace muchos años, esta se ajusta a varios sistemas y arroja resultados confiables como para ser tomados en cuenta.

Como un ejemplo se muestra en la figura 1 y 2, los resultados obtenidos de calcular los coeficientes de fugacidad con la ecuación de Peng Robinson para los sistemas binarios metano-agua y acetonitrilo– nitrometano, respectivamente. Las figuras muestran bajos niveles de dispersión, respecto a los obtenidos con pruebas experimentales, lo que indica porcentajes de error promedio inferiores al 1% [24].

Comparación de las fracciones molares calculadas con Peng Robinson y obtenidas experimentalmente para un sistema binario metanol-agua.
Figura 1.
Comparación de las fracciones molares calculadas con Peng Robinson y obtenidas experimentalmente para un sistema binario metanol-agua.


Comparación de las fracciones molares calculadas con Peng Robinson y obtenidas experimentalmente para un sistema binario Acetonitrilo -Nitrometano.
Figura 2.
Comparación de las fracciones molares calculadas con Peng Robinson y obtenidas experimentalmente para un sistema binario Acetonitrilo -Nitrometano.


En la literatura no se encuentra información detallada del modelamiento del equilibrio líquido-vapor de la mezcla acetona-cloroformo con la ecuación de Peng Robinson. Por lo tanto, validar si este modelo representa bien la mezcla estudiada será uno de los principales objetivos y aportes del presente trabajo.

Coeficiente de actividad

El coeficiente de actividad es la medida de la desviación de una disolución respecto al comportamiento ideal. Esta desviación se expresa a través de la función de exceso que depende de la concentración y la temperatura [25].

El coeficiente de actividad se expresa comúnmente a través de ecuaciones tales como: Readlich y Kister (1948) y los modelos de Margules (1985) y Van Laar (1910) [20,26], los cuales se aplican a mezclas binarias, pero no pueden ser extendidos a sistemas multicomponentes [11].

Según la figura 3, la mejor opción para el presente trabajo fue usar el modelo termodinámico de Van Laar, debido a las características de la mezcla cloroformo acetona, ya que esta cuenta con dos componentes polares [27, 22], no electrolitos y la presión de operación se encuentra entre 65 y 80 kPa. Adicionalmente, esta ecuación es muy usada debido a su flexibilidad, simplicidad y capacidad de representar bien muchos sistemas, presentado pequeñas desviaciones con los datos experimentales al modelar los coeficientes de fugacidad de mezclas binarias [28].

Diagrama de flujo para la selección de modelo para la determinación de coeficientes de actividad.
Figura 3.
Diagrama de flujo para la selección de modelo para la determinación de coeficientes de actividad.


Metodología

En esta sección se explica la metodología para la construcción de la curva de burbuja y rocío a partir de los modelos explicados y seleccionados anteriormente, mediante la herramienta de cómputo, Matlab 9.0. El procedimiento está enfocado en un sistema isotérmico, donde se predice el equilibrio líquido-vapor, a través de las presiones para cada una de las relaciones molares establecidas en un sistema binario.

Construcción de la curva de burbuja

En esta primera fase se asumió la composición en fase líquida X1 y X2 (acetona y cloroformo respectivamente) y se supusieron los coeficientes de fugacidad igual a 1(ϕi =1). Posterior a esto, se calculó la presión de vapor (Pi sat) para cada componente, usando la ecuación de Antoine a la temperatura de trabajo a través de la siguiente expresión (ecuación 2).

(2)

Para el uso de la ecuación 2 se necesitaron las constante de Antoine, A,B y C de cada uno de los componentes de la mezcla analizada (cloroformo, acetona), los cuales se presentan en la tabla 1 [12].

Tabla 1.
Constantes de Antoine para cada componente.
Constantes de Antoine para cada componente.


Una vez se definió la presión de saturación para cada componente, se calcularon los coeficientes de actividad (γ1 y γ2) con el modelo de Van Laar a la temperatura de análisis (T) y composición de la fase líquida, a través de las ecuaciones 3 y 4 respectivamente. Las constantes del sistema binario acetona-cloroformo (A 12 y A 21) son obtenidas de la literatura y corresponden a A12 = -0,8643 y A21 = -0,5899.

(3)

(4)

Calculados los coeficientes de actividad de la fase líquida y asumiendo los de fugacidad, se procedió a calcular la presión del sistema (P) usando la ecuación 5.

(5)

Con la presión calculada se determinaron las fracciones de vapor para cada componente a través de la ecuación 6. Posteriormente, se estimaron los coeficientes de fugacidad para cada componente ϕi , usando la ecuación 7.

(6)

(7)

Para usar de manera directa la ecuación 7 fue necesario incluir la constante universal de los gases (R) y calcular el coeficiente de fugacidad de cada componente en la mezcla (φi ), el coeficiente de fugacidad en el punto de saturación (φi sat) y el volumen de cada componente en el punto de saturación (Vi L ).

Como se mencionó anteriormente, el cálculo del coeficiente de fugacidad de los componentes en la mezcla se hizo a través de la ecuación de Peng Robinson (Ecuación 8).

(8)

Para la ecuación 8, por tratarse de una mezcla, los parámetros a y b deben ser calculados a través de la regla de mezcla cuadrática. Para este estudio el cálculo se realizó secuencial mentemediante las siguientes ecuaciones (9,10,11).

(9)

(10)

(11)

(12)

Donde los parámetros ai y bi son propios de cada componente y se calculan con los puntos críticos de la presión y temperatura (Tci y Pci), mediante la ecuación de Peng Robinson para un componente puro con las ecuaciones 13 y 14.

(13)

(14)

Para el caso del cálculo del factor compresibilidad (Z) se resuelve el polinomio de la ecuación cúbica generalizada (ecuación 15). Donde x =-2, u=1 y A’ y B’ conforman los coeficientes principales del polinomio.

(15)

Para el cálculo de los coeficientes de fugacidad en la zona de saturación se usaron las ecuaciones 7,8,9,10,11,12,13,14, asumiendo cada componente como puro y bajo las condiciones de saturación; es decir, los cálculos se hicieron con la presión de saturación de cada componente, calculada previamente con la ecuación 2.

Para el cálculo del volumen molar saturado de cada componente en fase líquida se empleó la ecuación de Raqquet (ecuación 16), la cual incluye el factor de compresibilidad (Zc ), volumen crítico de cada componente (Vc ) y la temperatura reducida (Tr ).

(16)

Con los parámetros calculados previamente, el siguiente paso fue recalcular la presión del sistema con la ecuación 17.

(17)

Adicionalmente se estableció una tolerancia de 0,001 entre la presión inicial (P) y la presión recalculada (Pf), evaluada a través de un proceso iterativo. Una vez se logró cumplir con dicha tolerancia, el proceso terminó y se inició el cálculo de la presión del sistema para una nueva fracción de líquido.

Construcción de la curva de rocío

Inicialmente, se estableció la temperatura del sistema T y la composición en fase vapor Y1 y Y2 , se asumieron los coeficientes de fugacidad y actividad igual a 1 y se calculó la presión del sistema con la ecuación 18.

(18)

Después, se calculó la composición de la fase líquida preliminar con la ecuación 9 y la presión calculada previamente.

(19)

Con las fracciones líquidas determinadas en la ecuación 19, se calcularon los coeficientes de actividad con el modelo de Van Laar en las ecuaciones 3 y 4, esto para calcular nuevamente la presión del sistema con la ecuación 18. Seguido a esto, se determinaron los coeficientes de fugacidad ϕi para cada componente con la ecuación, iniciando la iteración con la presión P.

Con la presión calculada y los coeficientes de fugacidad ϕi , se recalculó la composición de la fase líquida mediante la ecuación 19.

Antes de continuar, fue necesario asegurarse que la suma de las fracciones de líquido fuera igual a 1; de lo contrario, esta debe normalizarse a través de la ecuación 20.

(20)

Con las fracciones normalizadas (X_nor) se calcularon nuevamente los coeficientes de actividad con el modelo de Van Laar y se evaluó la diferencia con los calculados inicialmente. Si la diferencia era mayor a la tolerancia establecida (tolerancia=0,001) era necesario volver a calcular los coeficientes de actividad con las ecuaciones 2 y 3 de lo contrario se podía continuar el proceso de recalculó de la presión del sistema a través con la ecuación 17. Esta presión fue comparada con la calculada inicialmente a través de una tolerancia establecida (tolerancia=0,001). Si esta era alcanzada el proceso iterativo terminaba, obteniendo la presión real del sistema; de lo contrario, el programa debía volver a calcular la presión con la ecuación 17.

Esta rutina fue realizada en Matlab 9.0 a una temperatura de 50 °C y 35,17 °C para obtener las curvas y los resultados que se muestran en la siguiente sección.

Análisis de resultados

En la figura 4 se puede observar el equilibrio líquido-vapor para un sistema binario cloroformoacetona a una temperatura de 50 °C, obtenido del programa elaborado con la metodología anteriormente explicada. En esta figura se muestra un sistema azeótropo a la presión mínima 61,33 kPa donde la fracción molar liquida de la acetona es 0,4. Este comportamiento se debe a que los átomos de cloro del cloroformo forman puentes de hidrógeno con la acetona, produciendo fuertes enlaces de unión y presentando mejor afinidad en mezcla que en estado puro [29]. Por esta razón, se aumentan las interacciones entre las distintas moléculas, particularmente entre el oxígeno del carbonilo en la acetona y el hidrógeno del cloroformo [13]. Esto aleja al sistema de un comportamiento ideal bajo estas condiciones, causando una desviación negativa, la cual se da porque los valores de los coeficientes de actividad son menores a 1.

Diagrama P-X-Y para sistema acetona + cloroformo a 50 °C construido en Matlab 9.0.
Figura 4.
Diagrama P-X-Y para sistema acetona + cloroformo a 50 °C construido en Matlab 9.0.


Con el objetivo de comprobar la exactitud de la metodología y del modelo de Van Laar y Peng Robinson, las curvas fueron comparadas con las obtenidas experimentalmente a una temperatura de 50 °C [12] y 35,17 °C [30], y las elaboradas en el simulador de procesos químicos Aspen HYSYS y el programador comercial VLE-calc. Dichas comparaciones se pueden ver representadas en las figuras 5, 8,11 y 12, respectivamente.

En la figura 5 se muestra la comparación entre en diagrama P-X-Y creada en Matlab 9.0 y el obtenido con pruebas experimentales a 50 °C. Los resultados muestran que, desde la presión de saturación del cloroformo hasta la presión de saturación de la acetona, el porcentaje de error relativo fue inferior al 5%, a pesar del desfase entre estos dos diagramas. Esto indica que la metodología usada, los modelos termodinámicos y las ecuaciones de estados seleccionado para su elaboración, logran representar las tendencias y el comportamiento del equilibrio líquido-vapor en una mezcla de acetona-cloroformo a una temperatura de operación de 50 °C con una desviación baja para la presión de rocío y burbuja, principalmente para las fracciones molares de acetona inferiores a 0,4, pero una desviación alta para valores por encima de dicha fracción.

Comparación entre el diagrama P-X-Y para sistema acetona + cloroformo a 50 °C construido en Matlab 9.0 y obtenidas con pruebas experimentales.
Figura 5.
Comparación entre el diagrama P-X-Y para sistema acetona + cloroformo a 50 °C construido en Matlab 9.0 y obtenidas con pruebas experimentales.


La comparación de los valores de la presión de rocío y burbuja, calculados con los valores experimentales determinados a 50 °C, se muestra en las figuras 6 y 7. La figura 6 presenta altos niveles de desviación entre algunos de los valores de presión de rocío calculados y estimados experimentalmente y un bajo coeficiente de correlación lineal (0,80), lo que se traduce en un porcentaje de error relativo promedio de 3,24%. Por otro lado, la figura 7 muestra un mejor nivel de ajuste de los valores de presión de burbuja con un porcentaje de error relativo promedio de 2,5% y un alto valor del coeficiente lineal (0,95).

Comparación de los valores de la presión de rocío calculados con los valores experimentales, determinados a 50 °C.
Figura 6.
Comparación de los valores de la presión de rocío calculados con los valores experimentales, determinados a 50 °C.


Comparación de los valores de la presión de burbuja calculados con los valores experimentales, determinados a 50 °C.
Figura 7
Comparación de los valores de la presión de burbuja calculados con los valores experimentales, determinados a 50 °C.


A continuación, en las figuras 8 se muestra los resultados de comparar los diagramas PXY con los estimados experimentalmente a 35,17 °C. En esta figura se evidencia que la mezcla binaria de acenotacloroformo, a la temperatura de 35,17 °C, presenta un comportamiento azeotrópico a la presión mínima de 34,10 kPa por las mismas razones explicadas anteriormente. También se puede observar que las curvas construidas tienen menor nivel de desfase con los datos experimentales para esta temperatura, que para el caso de 50 °C, causando menor porcentaje de error.

Comparación entre el diagrama P-X-Y para sistema acetona-cloroformo a 31,7 °C construido en Matlab 9.0 y obtenidas con pruebas experimentales.
Figura 8.
Comparación entre el diagrama P-X-Y para sistema acetona-cloroformo a 31,7 °C construido en Matlab 9.0 y obtenidas con pruebas experimentales.


Lo anterior, también se puede observar en las figuras 9 y 10, donde se presenta la desviación de los valores calculados y obtenidos experimentalmente de la presión de burbuja y rocío, respectivamente. La figura 9 muestra que se tiene un buen ajuste, ya que la mayoría de los datos de la presión de burbuja calculados coinciden con los datos experimentales, causando un porcentaje de error relativo promedio de 1,9% y un alto coeficiente de relación lineal. La figura 10 muestra que los valores de presión de rocío, obtenidos mediante los modelos seleccionados, tienen un menor ajuste con los datos experimentales, causando un porcentaje de error relativo promedio de 2,3%.

Comparación de los valores de la presión de burbuja calculados con los valores experimentales, determinados a 35,17 °C.
Figura 9.
Comparación de los valores de la presión de burbuja calculados con los valores experimentales, determinados a 35,17 °C.


Comparación de los valores de la presión de rocío calculados con los valores experimentales, determinados a 35,17 °C.
Figura 10.
Comparación de los valores de la presión de rocío calculados con los valores experimentales, determinados a 35,17 °C.


El análisis de los resultados anteriormente expuestos permite concluir que a pesar de que los porcentajes de error son bajos y se representan las tendencias, el modelamiento nos muestra los porcentajes de error esperados, según lo observados en la literatura (<1%). También, se puede observar que se tiene mejor ajuste de los valores de presión de burbuja para las temperaturas de 35,17 °C y 50 °C, ya que se alcanzan porcentajes de error del 1,9% y 2,3%, respectivamente. Mientras que los valores de presión de rocío obtenidos alcanzaron porcentajes de error del 2,5% y 3,9%. Adicionalmente, estos porcentajes de error indican que las ecuaciones de Van Laar y Peng Robinson logran representar de mejor manera el sistema de cloroformo-acetona al disminuir la temperatura, debido a que con la temperatura de 35,17 °C se obtienen los menores porcentajes de error relativo. Esto se debe a que, al aumentar la temperatura del sistema, incrementa la probabilidad de que los dos componentes coexistan en estado líquido y bajo esta condición, la ecuación de Van Laar pierde precisión para calcular los coeficientes de actividad [31].

En la figura 11 se muestra la comparación del diagrama P-X-Y para el sistema acetonacloroformo construido en Matlab 9.0 y el arrojado por un programador comercial VLE-CALC. Como se puede ver en esta figura, el comportamiento de las dos gráficas sigue la misma tendencia; sin embargo, el creado en Matlab 9.0 sufre una leve sobrestimación en cuanto a las condiciones baja las cuales se da la formación del azeótropo, pero subestimación en el cálculo de la presión de saturación del componente 1 o acetona. Esto se puede dar porque el programador comercial trabaja con diferentes modelos para la construcción de la curva, no obstante, el porcentaje de error relativo promedio del 1%.

En la figura 12 se compara el diagrama de P-X-Y construido en Matlab 9.0 y el creado en Aspen HYSYS. Aspen HYSYS es uno de los simuladores de procesos químicos más conocidos debido a su alta precisión, cuya función es servir de apoyo para el diseño y modelado en la industria. En esta figura se puede ver que los dos diagramas se ajustan muy bien, principalmente a concentraciones bajas de acetona, donde no se ve desfase entre las curvas. También se puede observar que las condiciones bajo las cuales se da la formación del azeótropo son exactamente las mismas; sin embargo, para fracciones molares de acetona superiores a 0,6, el programa creado en Matlab 9.0 subestima el comportamiento de las curvas del equilibrio líquido vapor y esto lleva a que la presión de saturación de la acetona sea diferente.

Comparación entre el diagrama P-X-Y para sistema acetona + cloroformo a 50 °C construido en MATLAB y el arrojado por el programador ELV-CAL.
Figura 11.
Comparación entre el diagrama P-X-Y para sistema acetona + cloroformo a 50 °C construido en MATLAB y el arrojado por el programador ELV-CAL.


Comparación entre el diagrama P-X-Y para sistema acetona + cloroformo a 50 °C construido en MATLAB y el creado en Aspen HYSYS.
Figura 12.
Comparación entre el diagrama P-X-Y para sistema acetona + cloroformo a 50 °C construido en MATLAB y el creado en Aspen HYSYS.


Bajo los cuatro escenarios de comparación, las curvas tuvieron la misma tendencia y el desfase entre estas fue mínimo; lo que permite asegurar que el programa elaborado para construir las curvas de rocío y curvas de burbuja del equilibrio líquido vapor del sistema binario acetona-cloroformo tiene buen funcionamiento. Sin embargo, se recomienda continuar revisando el mismo sistema con otros modelos y ecuaciones que permitan mejorar el ajuste y alcanzar porcentajes de error inferiores al 1% con respecto a los datos experimentales.

Conclusión

Los resultados obtenidos con la metodología mostrada en este artículo para el modelamiento del equilibrio líquido-vapor, aplicada a un sistema binario (Acetona- cloroformo a 50 y 31,7 °C) no ideal azeótropo, usando el modelo de Van Laar y la ecuación cubica de estado de Peng Robinson, permiten concluir que los porcentajes de error son bajos y logran representar las tendencias; sin embargo, se recomienda continuar revisando el mismo sistema con otros modelos y ecuaciones que permitan mejorar el ajuste con respecto a los datos experimentales.

Con el modelo de Van Laar y la ecuación cúbica de estado de Peng Robinson se tiene mejor ajuste de los valores de presión de burbuja para las temperaturas de 35,17 °C y 50 °C, ya que se alcanzan porcentajes de error promedio del 1,9% y 2,3%, respectivamente. Mientras que los valores de presión de rocío obtenidos alcanzaron porcentajes de error promedio del 2,5% y 3,9%.

Los porcentajes de error indican que las ecuaciones de Van Laar y Peng Robinson logran representar de mejor manera el sistema de cloroformoacetona al disminuir la temperatura, debido a que a 35,17 °C se obtienen los mejores ajustes con los datos experimentales. Esto se debe a que, al aumentar la temperatura del sistema, incrementa la probabilidad de que los dos componentes co-existan en estado líquido y bajo esta condición, la ecuación de Van Laar pierde precisión para calcular los coeficientes de actividad.

Agradecimientos

Un agradecimiento a la Agencia Nacional de Hidrocarburos de Colombia (ANH), Ecopetrol S.A Colombia, Colciencias y al Grupo de investigación Recobro Mejorado de la Universidad Industrial de Santander, por apoyar esta investigación, realizada dentro del proyecto 56937: “Evaluación del uso de flue gas en la inyección cíclica de vapor como método de recobro mejorado, en campos del magdalena medio colombiano”.

Referencias bibliográficas

[1] Chasoy Rojas W. Determinación experimental del equilibrio liquido vapor del sistema etanoagua-glicerina (tesis maestría). Bogotá, Colombia: Universidad Nacional de Colombia; 2012.

[2] Zumalacárregui de Cárdenas L, Pérez Ones O, Hernández Castellanos FA, Cruz Lemus G. Modelación del equilibrio líquido-vapor a presión constante de mezclas etanol-agua utilizando redes neuronales artificiales. Tecnología Química. 2018;38(3):533-53.

[3] Marquez C, Ortega C, Oseguera J, Morales M. Simulación de la separación de una mezcla azeotropica de etanol agua, comparando métodos convencionales y no convencionales (sitio en internet). Researchgate. Disponible en:https://www.researchgate.net/publication/276848759_simulacion_de_la_separacion_de_una_mezcla_azeotropica_de_etanolagua_comparando_metodos_convencionales_y_no_convencionales. Acceso el 2 de diciembre de 2018.

[4] Harding S, Maranas C, Mc Donald, Floudas. Locating All Homogeneous Azeotropes in Multicomponent Mixtures. Ind. Eng. Chem. Res. 1997;36(1):160-78.

[5] Mandagarán B, Campanella E. Cálculo y Predicción de Azeótropos Multicomponentes con Modelos de Coeficientes de Actividad. Inf. Tecnol. 2008;19(5):73-84.

[6] Arango I, Villa A. Equilibrio multifásico de mezclas binarios y ternario s de carbonato de dietilo (DEC), agua, etanol. rev. ion 2012;25(1):43-49.

[7] Wang Y, Chien I. Unique Design Considerations for Maximum-Boiling Azeotropic Systems via Extractive Distillation: Acetone/ Chloroform Separation. Ind. Eng. Chem. Res. 2018;57(38):1-38.

[8] Szarawara J, Gawdzik A. Method of calculation of fugacity coefficient from cubic equations of state. Chem. Eng. Sci. 1989;44(7):1489-94.

[9] Soave G, Coto B. SRK, después de muchos años: ¿Qué se puede mantener y qué se debe cambiar? An. R. Soc. Esp. Quím. 2007:103(4): 14-21.

[10]Ghosh SK, Chopra SJ. Activity Coefficients from the Wilson Equation. Ind Eng Chem Proc Des Dev. 1975;14(3):304-8.

[11] Herrera M, Beltrán JC. Uso de coeficientes de actividad experimentales a dilución infinita para validar simulaciones de proceso. Tecnol. Ciencia Ed. 2012;27(2):69-79.

[12]O’Connell JP, Haile JM. Thermodynamics: Fundamentals for Applications. USA: Reissue edition. Cambridge: Cambridge University Press; 2011.

[13]Bello M, Torres D. Equilibrio liquido-vapor en soluciones no ideales, uso de la ecuación de van Laar para un sistema acetona- cloroformo (Tesis de maestría). México: Universidad Autónomo de Puebla; 2010.

[14]Anguebes-Franseschi F, Bolaños-Reynoso E, Castro-Montoya A, Cantú-Lozano D. Separación de la Mezcla Azeotrópica AcetonaMetanol con Cloruro de Litio. Inf. Tecnol. 2004;15(2):111-6.

[15]Acosta-Esquijarosa J, Rodríguez-Donis I, Lodeiro-Prieto L, Nuevas-Paz L. Separación de los componentes de la mezcla azeotrópica acetona-n -hexano mediante un proceso combinado: extracción líquido-líquido y destilación discontinua. Revista CENIC Ciencias Químicas (sitio en internet). Disponible en: http://www.redalyc.org/resumen.oa?id=181618068004. Acceso el 2 de diciembre de 2018.

[16]Prausnitz J, Lichtenthaler R, Azevedo E. Molecular Thermodynamics of Fluid-Phase Equilibria. New Jersey: Prentice Hall; 1998.

[17]Moran M, Shapiro H. Fundamentos de termodinámica técnica. Barcelona: Reverté; 2012.

[18]Dimian A, Bildea C, Kiss A. Phase equilibria. En: Computer Aided Chemical Engineering. Ámsterdam. Elsevier Science Limited. 2003. p. 201-51.

[19]Rajadell F, Mojilla J. Termodinamica de gases. Fugacidad. En: Termodinámica química. Catellón de la Plana. Universitat Jaume; 2005. p.177-81

[20]Reid R. The Properties of Gases and Liquids. New York: McGraw-Hill; 1987.

[21]Hu J, Wang R, Mao S. Some useful expressions for deriving component fugacity coefficients from mixture fugacity coefficient (Sitio en internet). ResearchGate. Disponible en: https://www.researchgate.net/publication/222208355_Some_useful_expressions_for_deriving_component_fugacity_coefficients_from_mixture_fugacity_coefficient. Acceso el 12 de diciembre de 2018.

[22]Martínez A de L, Camacho JV, Bajo JL. Termotecnia básica para ingenieros químicos: Bases de Termodinámica Aplicada: Ediciones la universidad de Castilla-La Mancha; 2003.

[23]Chakraborty A, Kumar R, Gbechoevi A. Fugacity coefficient estimation of pure Carbon dioxide and water: Approach aimed at comparative study through Van der Waals and Soave-Redlich-Kwong EOS Modelling (Sitio en internet). Researchgate. Disponible en: https://www.researchgate.net/publication/. Acceso el 16 diciembre del 2018.

[24]Corream F. Predicción del equilibrio líquido – vapor por el método UNIFAC – DORTMUND para el cálculo de los coeficientes de actividad mediante el software ‘ELVIQ’ (sitio en internet). Confederación Interamericana en Ingeniería Química. Disponible en: http://www.ciiq.org/varios/peru_2005/Trabajos/V/3/5.3.02.pdf. Acceso el 16 diciembre del 2018.

[25]Marigorta A. Análisis del equilibrio líquido-vapor a 141,3 kPa de mezclas binarias que contienen metanol con alcanos (C5,C6) y con esteres alquílicos (sitio en internet). Universidad de La Laguna. Disponible en: https://dialnet.unirioja.es/servlet/tesis?codigo=838. Acceso el 12 febrero de 2019.

[26]Lerol J, Masson J, Renon H, Fabries J , Sannier H. Accurate Measurement of Activity Coefficient at Infinite Dilution by Inert Gas Stripping and Gas Chromatography. Ind Eng Chem Proc Des Dev. 1977;16(1):139-44.

[27]Rodriguez E, Ortega O del S, Sulroca Y, Delgado J. Introduction of a new mixture of solvents in Thin Layer Chromatography (sitio en internet). Centro Nacional de Toxicología. Disponible en: http://www.medigraphic.com/cgi-bin/new/resumenI.cgi?IDARTICULO=49423. Acceso el 28 de octubre 2018.

[28] Jain T, Singh P. Estimation of Activity Coefficients for binary mixture VLE Data using MATLAB. IARJSET. 2016;3:141-47.

[29]Bedoya F, Vasquez D. Aplicabilidad de la mezcla cloroformo-acetona, equilibrio de fases con ecuaciones de estado. Ingeniería química. Colombia: Universidad de Medellín; 2006. Report No.: 1.

[30]Kamath G, Georgiev G, Potoff J. Molecular Modeling of Phase Behavior and Microstructure of Acetone-Chloroform-Methanol Binary Mixtures. J. Phys. Chem. B. 2005;41(1):19463– 473.

[31]Campbell A, Musbally G. Vapor pressures and vapor- liquid equilibria in the systems: (1) acetone - chloroform, (2) acetone -carbon tetrachloride, (3) benzene - carbon tetrachloride. Can. J. Chem. 1970;48(1):3173- 184.

Notas de autor

*maría.sandoval@correo.uis.edu.co

Información adicional

Cita: Sandoval Martínez MI, Muñoz Navarro SF, Martínez Jiménez HJ. Estimación del equilibrio líquido-vapor del sistema binario acetona-cloroformo mediante el modelo termodinámico de Van Laar y Peng Robinson. rev. ion. 2020;33(2):49-60. doi:10.18273/revion.v33n2-2020004

Secciones
Revista ION
ISSN: 0120-100X
Vol. 33
Num. 2
Año. 2020

Estimación del equilibrio líquido-vapor del sistema binario acetona-cloroformo mediante el modelo termodinámico de Van Laar y Peng Robinson

María IsabelSamuel FernandoHumberto Jose Sandoval MartínezMuñoz NavarroMartinez Jiménez
Universidad Industrial de Santander (UIS)Universidad Industrial de Santander (UIS)Universidad Industrial de Santander (UIS),ColombiaColombiaColombia
Contexto
Descargar
Todas