Equilibrio multifásico de mezclas binarias y ternarias
de carbonato de dietilo (DEC), agua y etanol

Isabel Cristina Arango Montoya; Aída Luz Villa*

Grupo Catálisis Ambiental, Departamento de Ingeniería Química, Universidad de Antioquia,
Carrera 53 Nº61-30 Medellín, Colombia.
*alvilla@udea.edu.co

Fecha Recepción: 22 de diciembre de 2011
Fecha Aceptación: 30 de abril de 2012



Resumen

En este trabajo se reporta la solubilidad del DEC en agua, y la curva de solubilidad del sistema ternario carbonato de dietilo (DEC) + agua + etanol entre 296,15 y 338,15K. Se confirma la presencia de azeótropos binarios con las graficas del equilibrio líquido vapor x-y de los constituyentes binarios a 353,15K y 363,15K y la ausencia de azeótropos ternarios mediante la regla azeotrópica. Los datos experimentales se correlacionan con los modelos termodinámicos NRTL y UNIQUAC, encontrándose que ambos modelos representan satisfactoriamente el sistema. Los parámetros de interacción binaria para el sistema ternario y sus sistemas binarios se reportan en función de la temperatura.

Palabras clave: solubilidad, azeótropos, aditivos, DEC.

Multiphase equilibrium for binary and ternary mixtures
containing diethyl carbonate (DEC), water and ethanol

Abstract

The solubility of DEC in water and the solubility curves of the diethyl carbonate (DEC) + water + ethanol ternary system between 296.15 and 338.15K are reported. The results confirm the presence of binary azeotropes using plots x-y of the binary constituents at 353.15K and 363.15K and the absence of ternary azeotropes by the azeotropic rule. The experimental data were correlated with the NRTL and UNIQUAC models, showing that both models provide a good representation of the system. The binary interaction parameter for the ternary system and their binary systems are reported in function of temperature in this work.

Keywords: solubility, azeotropes, additive, DEC.



Introducción

Los aditivos de combustible son alcoholes y éteres que contienen oxígeno, y logran aumentar el octanaje, mejorar la combustión y reducir las emisiones de CO y NOx. Recientemente se ha generado interés por el desarrollo de aditivos para combustibles que mejoren la eficiencia de la reacción de combustión y que no representen un peligro para el ambiente. El carbonato de dietilo (DEC) es un compuesto de gran interés para la industria química y en particular para la industria de los combustibles porque puede ser utilizado como aditivo, presenta un alto contenido de oxígeno (40,7%p), comparado con otros aditivos como el MTBE (18,2%p) y el etanol (34,8%p), mayor contenido de energía, menor presión de vapor y mayor distribución en la gasolina que en el agua, respecto al DMC [1]. Dado el uso potencial del DEC como aditivo, se han realizado estudios sobre su síntesis y purificación. Una ruta promisoria para la producción de DEC es la síntesis directa a partir de CO2 y etanol. En esta reacción se obtiene una mezcla de DEC + agua + etanol [2], siendo entonces necesario la separación de estos compuestos para la purificación del DEC.
Un recurso preliminar muy útil, son las curvas residuales que se realizan en el diagrama ternario y que muestran las líneas por las cuales se sigue una destilación a una composición determinada, indicando la evolución del equilibrio. Éstas sirven para predecir, sin cálculos intensivos, la viabilidad de la separación y se puede realizar un mapa aproximado a partir de puntos singulares como los azeótropos y puntos de ebullición de componentes puros. El sentido de las flechas es el de puntos de ebullición crecientes, es decir, las curvas parten desde un componente puro o un azeótropo de menor punto de ebullición y llegan a un componente puro o azeótropo de mayor punto de ebullición, definiendo de esta forma fronteras de destilación [2].
Los sistemas que contienen DEC cuyo equilibrio líquido-líquido (ELL) se ha reportado son: agua + ácido propiónico + DEC [3], DEC + MTBE [4], dietil y dimetil carbonato con agua, etanol, heptano y metanol a 298,15K y 1atm de presión [5] ajustándose los datos adecuadamente con los modelos UNIQUAC modificado y UNIQUAC extendido. Se ha reportado también estudios del ELL para sistemas similares al de interés como el DMC + agua + metanol, a diferentes temperaturas y su ajuste termodinámico con los modelos UNIQUAC y NRTL [6]. Para el sistema DMC + alcano + metanol, se reportaron estudios termodinámicos con los modelos UNIQUAC y NRTL [7]. Para el Equilibrio líquido vapor (ELV), el par binario etanol + agua ha sido ampliamente estudiado, para el sistema DEC + etanol se han reportado parámetros de interacción binaria para los modelos Wilson y UNIQUAC [8]; sin embargo, no hay reportes del ELV del par binario DEC + agua.
Se han reportado los puntos azeotrópicos, en porcentaje masa, de los sistemas binarios DEC + agua a 70%p de DEC (364,15K, 101,32kPa) y de la mezcla etanol + agua a 96%p de etanol (351,35K, 101,32kPa) [9]. El objetivo de este trabajo es obtener el diagrama x-y del ELV del sistema binario DEC + agua, la solubilidad del DEC en agua a diferentes temperaturas, evaluar la presencia de azeótropos y reportar los parámetros de interacción binaria para el sistema ternario y sus constituyentes binarios.


Metodología

Solubilidad
La solubilidad se determinó entre 296,15 y 341,15K, preparándose mezclas de agua + DEC con 2mL de agua y 2mL de DEC, cuya temperatura se mantuvo constante en un baño de agua. Después de 2h de agitación y 3h en reposo, se analizó una muestra de la fase acuosa y otra de la fase orgánica mediante un cromatógrafo de gases (Agilent technologies 5975c VL MSD), equipado con un detector de conductividad térmica (TCD) que se mantuvo a 523,15K. Para los análisis se utilizó la columna capilar megabore DB-Wax (30m x 0,320mm x 0,50μm), la temperatura del puerto de inyección fue 498,15K, el gas de arrastre fue helio con un flujo de 2mL/min y la partición de 100:1. Después de 1min a 313,15K, la temperatura del horno se aumentó a 343,15K con una velocidad de 293,15K/ min y posteriormente a 303,15K/min hasta alcanzar una temperatura de 393,15K, manteniéndose esta temperatura por 1min. El volumen inyectado fue de 1,0μL y las muestras se analizaron por triplicado, utilizándose curvas de calibración en el rango de concentración requerido, para determinar la composición en peso.

Curva binodal
Para el sistema DEC + agua + etanol, se determinó la curva binodal a 83,09kPa siguiendo el método del punto de nube [6]. Se prepararon mezclas binarias heterogéneas de composición conocida de DEC y agua y se adicionó etanol hasta alcanzar una mezcla homogénea. Las pruebas se realizaron por triplicado para cada temperatura.

Determinación del ELLV y el ELV
Las concentraciones en el ELV isotérmico se determinaron en un HS acoplado a un cromatógrafo de gases, el cual consiste en un cromatógrafo de gases (Agilent 5975c VL MSD) y un headspace automuestreador (Agilent G1888), HSGC. En el sistema HSGC se producen tres etapas sucesivas: regulación de la temperatura con un termostato de precisión (± 0,2°C), la presurización y la inyección, las especificaciones del método utilizado para el headspace se muestran en la Tabla 1. Se evaluaron los equilibrios de los pares binarios DEC + agua, etanol + DEC, agua + etanol y del sistema ternario DEC + agua + etanol, por triplicado a 353,15 y 363,15K. El factor de respuesta de cada componente se determinó por medio del método MHE (del término en inglés Multiple Headspace Extraction), para esto se introdujo una muestra (6μL) en un vial, que fue completamente evaporada. Para determinar el ELV, aproximadamente 4mL de la muestra se adicionaron a un vial de vidrio de 20mL, luego las muestras se llevaron al HS con un tiempo de equilibrio de 1h. El ELV del sistema binario metanol (1) + dimetil carbonato (2) se midió con este procedimiento y se comparó con reportes de la literatura [10], encontrándose una desviación máxima de yi de 1,35%, lo que indica la confiabilidad del procedimiento.

Correlación de los datos experimentales
Los datos experimentales se correlacionaron con los modelos termodinámicos UNIQUAC [11] y NRTL [12]. Los parámetros de interacción binaria se calcularon minimizando las diferencias entre el coeficiente de distribución experimental y calculado para el sistema líquido-líquido, Ecuación 1 [13], mientras que para el sistema líquido-vapor la función de minimización (OF) incluye los coeficientes de actividad experimental y calculado, Ecuación 1. [14]. El algoritmo de optimización utilizado en ASPEN PLUS V7.1 fue el de Britt-Lueck. Los parámetros estructurales para el modelo UNIQUAC, R y Q, se obtuvieron de la literatura [4].

Donde i = componentes; k = 1, 2, 3,… n puntos, Ki,exp es el coeficiente de distribución obtenido de los datos experimentales para el componente i y Ki,calc es el coeficiente de distribución calculado para el componente i.

Donde i = componentes; j = fases, Γijexp es el coeficiente de actividad experimental para el componente i y Γijcal es el coeficiente de actividad calculado para el componente i.
La calidad de las correlaciones para ambos equilibrios fue evaluada por la desviación estándar, RMSD, obtenida por la diferencia entre la fracción mol experimental y calculada, Ecuación 3.


Resultados y Discusión

Solubilidad
La solubilidad de DEC en agua y de agua en DEC, determinada experimentalmente entre 296,15 y 338,15K, se presenta en la Tabla 2. Estas composiciones son consistentes con los datos de la curva binodal a las temperaturas estudiadas mostrada en la Figura 1. La solubilidad de DEC en agua reportada en la literatura a 293,15K es de 1,8g 100mL [15]. La solubilidad de ambos componentes no se afecta apreciablemente por el aumento de la temperatura.

Curva binodal
Las curvas binodales a las temperaturas estudiadas, 296,15K y 338,15K, y a presión atmosférica se presentan en la Figura 1. Todas las concentraciones del sistema terniario son expresadas en fracción mol. La incertidumbre u(x1), se calculó evaluando la desviación estándar de los datos experimentales y se dividió entre 1,96 [16],

asumiendo un nivel de confianza del 95% y una distribución normal de los datos. La incertidumbre u(x1) varió entre 0,0018 y 0,0203.
El sistema presenta un comportamiento de fase tipo I, donde dos componentes (agua y DEC) son parcialmente miscibles y dos pares de componentes son completamente miscibles (agua + etanol y etanol + DEC).

Azeótropos
Los puntos azeotrópicos binarios corresponden a composiciones de la fase vapor y líquida iguales y se identificaron de dos formas: i) haciendo la función de la Ecuación 4 igual a cero; ii) determinando las composiciones que corresponden al cruce de la línea de 45° mostradas en la Figura 2. Estos resultados confirman la presencia de dos puntos azeótropicos binarios para los sistemas (agua + etanol) y (agua + DEC) y son comparados con los reportados en la literatura obtenidos isobáricamente, Tabla 3.

Donde N1 es el número de nodos para los componentes puros en la mezcla, S1 es el número de puntos de silla para los componentes puros, N2 es el número de nodos entre los azeótropos binarios, S2 es la cantidad de puntos de silla dentro de los azeótropos binarios, N3 and S3 son el número de nodos y sillas ternarias respectivamente.
La combinación de los puntos singulares proporciona una restricción con la que se puede estimar las separaciones posibles en una destilación continua [18]. Para determinar la cantidad de sillas y nodos presentes en cada uno de los sistemas, se realiza un mapa aproximado de curvas residuales: primero se etiqueta en el diagrama ternario las temperaturas de ebullición normales de los componentes puros y se ubican los azeótropos binarios, posteriormente se dibujan las flechas en los lados del triángulo entre cada par de especies adyacentes en sentido creciente de temperaturas de ebullición, por último se establecen los puntos nodos y sillas con la ayuda de la Figura 3 [17].
El diagrama mostrado en la Figura 3 muestra que los componentes puros agua y DEC son nodos estables, resultando N1 = 2, mientras que el etanol es un punto de silla, dando como resultado S1 = 1. Además el sistema binario agua + etanol presenta un azeótropo con un nodo inestable, lo que conduce a N2 = 1, mientras que el sistema binario agua + DEC presenta un azeótropo que es un punto de silla, por lo que se obtiene S2 = 1. Bajo esta situación, la regla azeotrópica, Ecuación 5.
Se cumple si y solo si N3 = S3. Sin embargo, aún no hay evidencia experimental de que exista más de un punto azeótropico ternario en un sistema de tres componentes [19]. Esto implica que N3 = 0 y S3 = 0, lo cual es consistente con la observación experimental de este estudio, la ausencia de azeótropos ternarios.

Correlaciones
Los parámetros de interacción binaria τij, Tabla 4 y Tabla 5, se determinaron en función de la temperatura para ambos modelos con las constantes aij y bij, Ecuación 6 para el modelo UNIQUAC y Ecuación 7 para el modelo NRTL.

Los modelos NRTL y UNQUAC proporcionan una buena representación de las composiciones de equilibrio del sistema estudiado tanto cuantitativamente como cualitativamente. Las desviaciones de las composiciones calculadas respecto a las experimentales para el sistema ternario líquido-líquido agua (1) + DEC (2) + etanol (3) y para el equilibrio líquido-vapor para sus constituyentes binarios se muestran en las Tablas 6 y 7, respectivamente.


Conclusiones

Se presentan los resultados experimentales de la solubilidad del sistema parcialmente miscible DEC + agua a diferentes temperaturas, y los resultados experimentales para la curva binodal del sistema ternario agua + etanol + DEC entre 296,15 y 341,15K y presión atmosférica.
Se comprobó la presencia de un azeótropo de mínimo punto de ebullición para los sistemas binarios etanol + agua y agua + DEC. Para el sistema ternario no se hallaron azeótropos ternarios, y se confirmó mediante la regla azeotrópica. Se reportan los parámetros de interacción binaria para el modelo NRTL y el modelo UNIQUAC en función de la temperatura, representando ambos modelos satisfactoriamente el sistema.


Agradecimientos

Los autores agradecen el apoyo financiero del Ministerio de Agricultura y Desarrollo Rural de Colombia, Química Básica y la Universidad de Antioquia a través del proyecto 045-2007D3608- 66-07.


Referencias

[1] Briggs DN, Lawrence KH, Bell AT. An investigation of carbon-supported CuCl2/PdCl2 catalysts for diethyl carbonate synthesis. Appl. Catal. A: Gen. 2009;366(1):71-83.

[2] Fan M, Zhang P, and Ma X. Study on Wackertype catalysts for catalytic synthesis of diethyl carbonate from ethyl nitrite route. Fuel. 2007;86(5-6):902-5.

[3] Bilgin M, Birman I. Separation of propionic acid by diethyl carbonate or diethyl malonate or diethyl fumarate and the synergistic effect of phosphorus compounds and amines. Fluid Phase Equilib. 2010;292(1-2):13-19.

[4] Chen Y, Fu M, Zhang Y, Chen J. Quaternary Liquid–Liquid Equilibria for Fuel Additive Systems Containing Diethyl Carbonate or 1,1-Dimethylethyl Methyl Ether at 298.15K and Atmospheric Pressure. J. Solution Chem. 2009;38(8):1029-41

[5] Chen Y, Cao C. Liquid−Liquid Equilibria for Three Quaternary Systems of Water, 2-Propanol, and Dimethyl Carbonate with Isooctane or Toluene or Heptane at 298.15K. J. Chem. Eng. Data 2009;54(6):1793-8.

[6] Choe J, Song KH. Liquid−Liquid Equilibrium for the Ternary Systems Nonane + Dimethyl Carbonate + Methanol and Decane + Dimethyl Carbonate + Methanol at 298.15K. J. Chem. Eng. Data. 2002;47(5):1253-5.

[7] Lee LS, HC Shen. Azeotropic Behavior of a Water + n-Propanol + Cyclohexane Mixture Using Cyclohexane as an Entrainer for Separating the Water + n-Proponal Mixture at 760 mmHg. Ind. Eng. Chem. 2003;42(23):5905-14.

[8] Rodríguez A, Canosa J, Domínguez A, Tojo J. Isobaric Phase Equilibria of Diethyl Carbonate with Five Alcohols at 101.3kPa. J. Chem. Eng. Data. 2002;48(1):86-91.

[9] Pacheco MA, CL Marshall. Review of Dimethyl Carbonate (DMC) Manufacture and Its Characteristics as a Fuel Additive. Energ Fuel. 1997;11(1):2-29.

[10] Ma XB, Liu XG, Li ZH, Xu GH. Vapor-liquid equilibria for the ternary system methanol + dimethyl carbonate + dimethyl oxalate and constituent binary systems at different temperatures. Fluid Phase Equilib. 2004; 221(1-2):51-6.

[11] Abrams DS, Prausnitz JM. Statistical thermodynamics of liquid mixtures: A new expression for the excess Gibbs energy of partly or completely miscible systems. AIChE J. 1975;21(1):116-28.

[12] Renon H, Prausnitz JM. Local compositions in thermodynamic excess functions for liquid mixtures. AIChE J. 1968;14(1):135-44.

[13] Barrera R, Villa AL, Montes de Correa C. Liquid− Liquid Equilibrium for the Water + Acetonitrile + Limonene System at Different Temperatures. J. Chem. Eng. Data. 2005;50(4):1353-6.

[14] Pereiro AB, Rodríguez A, Canosa J, Tojo J. Vapor-liquid equilibria for systems of diethyl carbonate and ketones and determination of group interaction parameters for the UNIFAC and ASOG methods. Fluid Phase Equilib. 2005;235(1):83-91.

[15] Chemical Safety Information from Intergovernmental Organizations. Disponible en: http://www.inchem.org. Acesso en mayo de 2011.

[16] Chirico RD, Frenkel M, Diky VV, Marsh KN, Wilhoit RC. ThermoMLAn XML-Based Approach for Storage and Exchange of Experimental and Critically Evaluated Thermophysical and Thermochemical Property Data. 2. Uncertainties. J. Chem. Eng. Data. 2003;48(5):1344-59.

[17] Doherty MF, Perkins JD. On the dynamics of distillation processes--III : The topological structure of ternary residue curve maps. Chem. Eng. Sci. 1979;34(12):1401-14.

[18] Modla G, Kopasz A, Lang P. Feasibility of Separation of Ternary Mixtures by Pressure Swing Batch Distillation. In Computer Aided Chemical Engineering, C.A.O.d.N. Rita Maria de Brito Alves and B. Evaristo Chalbaud, Editors. Salvador: Elsevier;2009.p.1017-22.

[19] Wen CC, Tu CH. Vapor-liquid equilibria for binary and ternary mixtures of ethanol, 2-butanone, and 2,2,4-trimethylpentane at 101.3kPa. Fluid Phase Equilib. 2007;258(2):131-9.