Equilibrio de Fases2
Equilibrio de Fases2
Equilibrio de Fases2
discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/39379495
Article
Source: OAI
CITATIONS READS
0 208
1 author:
Ana M Blanco-Marigorta
Universidad de Las Palmas de Gran Canaria
28 PUBLICATIONS 333 CITATIONS
SEE PROFILE
All content following this page was uploaded by Ana M Blanco-Marigorta on 25 March 2014.
UNIVERSIDAD DE LA LAGUNA
FACULTAD DE CIENCIAS QUIMICAS
-1-
INDICE
Resumen ................................................................................................................... 6
Abstract .................................................................................................................... 6
1. INTRODUCCION ............................................................................................... 7
1.1. Antecedentes bibliogrficos .......................................................................... 12
1.2. Comparacin de datos bibliogrficos de ELV para las mezclas con metanol . 20
I. PARTE TEORICA
-3-
5.2. Modelos de composicin local ...................................................................... 60
5.2.1. Modelo de Wilson .............................................................................. 60
5.2.2. Modelo NRTL .................................................................................... 61
5.2.3. Modelo UNIQUAC ............................................................................ 63
-4-
11.2. Miscibilidad de la mezcla metanol+n-hexano ......................................... 124
11.3. Resultados de composicin-densidad-volumen de exceso. Correlaciones . 127
11.3.1. Mezclas metanol+n-alcano ............................................................... 128
11.3.2. Mezclas metanol+ster de alquilo .................................................... 131
11.4. Datos de equilibrio lquido-vapor. Test de consistencia. Correlaciones .. 144
11.4.1. Datos de ELV isobrico para las mezclas metanol+alcanos ............. 151
11.4.2. Datos de ELV isobrico para las mezclas metanol+steres de alquilo. 158
11.5. Azetropos ............................................................................................... 192
11.5.1. Azetropos en las mezclas metanol+n-alcanos ................................ 194
11.5.2. Azetropos en las mezclas metanol+steres de alquilo .................... 197
V. CONCLUSIONES Y BIBLIOGRAFIA
VI. APENDICES
-5-
Resumen
En esta Memoria de investigacin se presenta un estudio experimental y terico realizado sobre
los equilibrios lquido-vapor, en condiciones isobricas de 141,3 kPa, de diecisiete mezclas binarias for-
madas por el metanol y dos alcanos (C5, C6) y un conjunto de quince alcanoatos (desde el metanoato al
butanoato) de alquilo (desde el metilo al butilo), con la excepcin del metanoato de metilo. La parte ex-
perimental consiste en la determinacin de los datos de equilibrio, verificando previamente el procedi-
miento experimental con la mezcla metanol+agua. Para todos los sistemas estudiados se obtienen las
curvas patrones de densidad-concentracin, cuya calidad se comprob con las de volmenes de exceso,
los cuales tambin clarificaron el comportamiento de las mezclas.
En el tratamiento de los datos de equilibrio, la fase vapor fue considerada como no-ideal y los
resultados muestran una desviacin positiva de la ley de Raoult; las mezclas metanol+alcano son alta-
mente no-ideales. La consistencia de los valores encontrados se verifica con diversos mtodos y la corre-
lacin de las magnitudes procedentes de los equilibrios lquido-vapor se realiza con los modelos de com-
posicin local Wilson, NRTL y UNIQUAC, as como con la ecuacin de las zetas, que se modifica en este
trabajo utilizando coeficientes dependientes de la temperatura. La modificacin efectuada produce buenos
resultados en la correlacin de datos de equilibrio.
Las mezclas formadas por metanol+(pentano, hexano, etanoato de metilo, propanoato de metilo,
metanoato de etilo, etanoato de etilo, metanoato de propilo) presentan azetropo mnimo a la presin de
trabajo.
La modelizacin terica se realiza mediante la aplicacin de las teoras de contribucin de grupos
ASOG y UNIFAC en dos de sus versiones. La prediccin de los coeficientes de actividad realizada con
dichos modelos es aceptable con los modelos de ASOG y el UNIFAC-modificado, si bien las diferencias
aumentan en las mezclas con steres de alto peso molecular.
Abstract
In this research Report an experimental and theoretical study on vapor-liquid equilibrium is pre-
sented at isobaric conditions of 101.3 kPa, for seventeen binary mixtures formed by methanol and two
alkanes (C5, C6) and a set of fifteen alkyl (from methyl to butyl) alkanoates (from methanoate to buta-
noate), with the exception of methyl methanoate. Experimental part consists of measuring the equilibrium
values, previously checking the experimental procedure for methanol+water mixture. For all systems
studied the standard curves of density-concentration are obtained, and are those checked with the excess
volumes which also clarified the behaviour of the mixtures.
In the treatment of equilibrium data, the vapor phase was considered as non-ideal and the results
show a positive deviation from Raoults law; the methanol+alkane mixtures are highly non-ideal. The
consistency of the obtained values is checked using several methods and the correlations of the different
magnitudes from vapor-liquid equilibria are carried out using the local composition models, Wilson, NRTL
and UNIQUAC and the zetas equation which is modified for this work using temperature depending coe-
fficients. The modification carried out give place to good results in correlation of equilibrium data.
The mixtures formed by methanol+(n-pentane,n-hexane, methyl ethanoate, methyl propanoate,
ethyl methanoate, ethyl ethanoate, propyl methanoate) show a minimum azeotrope at the working pressu-
re.
Theoretical modelisation is carried out by using some group contribution models, ASOG and
UNIFAC using two of its versions. The prediction of the activity coefficients made by those models is
acceptable with the ASOG and modified-UNIFAC models although the differences increase when the high
molecular weight esters are in the mixture.
-6-
INTRODUCCION
el siempre costoso trabajo experimental. Uno de los motivos de este rpido crecimiento
dos alcanzados por dicha materia han permitido un cambio profundo en la dinmica de
Otra de las razones que contribuyen a la importancia del estudio de los fludos es
miento de las fases a partir de las propiedades involucradas. En este sentido, las teoras
-7-
de fluidos han experimentado una evolucin progresiva. En lo que se refiere al equili-
brio lquido-vapor (ELV), objeto de este trabajo de Tesis Doctoral, cabe mencionar, por
Wilson, el NRTL y el UNIQUAC que, debido a sus fundamentos fisicoqumicos, son ca-
ricos ya citados, tambin han alcanzado importancia creciente los modelos de contribu-
cients, de Fredenslund y col. (1975) en sus distintas versiones. En ellos, a partir de cier-
que caracterizan las interacciones entre parejas de grupos. Con dichos parmetros de
interaccin puede predecirse el ELV de mezclas similares que contengan a los grupos
do faltan valores para los parmetros de interaccin que corresponden a algn grupo
lizada que permita su utilizacin en todos los casos ni la adecuada exactitud en la pre-
diccin del comportamiento de los fludos, en particular del equilibrio de fases en siste-
de datos cualificados y por tanto, la modificacin tcnica de los equipos, de forma que
-8-
se puedan determinar con precisin las magnitudes termodinmicas involucradas en el
equilibrio de fases, como son la presin, la temperatura y las composiciones de las fases
coexistentes.
Nuestro grupo de trabajo viene desarrollando desde hace varios aos un amplio
del ELV de sistemas binarios, se han publicado bastantes artculos sobre mezclas de
steres con alcanoles, tanto normales como isomricos, a diferentes presiones. Sin em-
bargo, en todos ellos siempre fue excluido el metanol con objeto de concederle un tra-
a presin constante, de las fases lquida y vapor para un conjunto de sistemas binarios
formados por metanol, como componente comn a todos, con alcanos (C5 y C6) y alca-
noatos (del metanoato al butanoato) de alquilo (de metilo a butilo), todos de bajo peso
ticular con u,v=1,2,3,4, con la excepcin del metanoato de metilo, por la alta volatilidad,
isobricamente para todas las mezclas a 141,3 kPa, ligera sobrepresin de la atmosfrica
normal que facilita el trabajo experimental por su alta volatilidad del metanol, no en-
contrndose datos en la literatura en esas condiciones para mezclas con alcanos que
mencionadas y otras que se expondrn en esta Memoria, es claro que presenta un inters
-9-
de distintos mtodos, el ms ampliamente utilizado emplea fuel-oil residual para la ob-
tencin de gas de sntesis (CO+H2); luego, con este gas y en determinadas condiciones
nol y algunas mezclas de hidrocarburos) e incluso directamente para tal fin. El conoci-
viene el metanol se caracteriza por el carcter fuertemente polar de esta sustancia, con
mente, las mezclas binarias alcanol+alcano, altamente no-ideales, son sistemas estan-
con asociacin por puentes de hidrgeno. Por ello, para este trabajo tambin se ha esco-
-10-
vestigacin es igualmente conseguir medidas de ELV, aportando datos inexistentes a la
cos que van a emplearse en la prediccin y reproduccin de los sistemas en estudio. In-
referencias sobre el ELV de algunas de las mezclas que van a ser objeto de estudio, tanto
ment a la presin de 141,3 kPa. No obstante, todas sern utilizadas para la compara-
establecido varias partes claramente diferenciadas. Una primera seccin recoge la parte
les y los modelos tericos, tanto de correlacin como de estimacin que van a utilizarse.
sis. Con la aplicacin de los diferentes modelos de prediccin a los datos obtenidos se
configura la cuarta seccin de la Memoria. Las ltimas secciones hacen referencia a las
-11-
1.1. Antecedentes bibliogrficos
elegido. Ello permite, por una parte, destacar el inters del trabajo que se pretende, al
observar lo ya realizado por otros autores y por otro, utilizar los resultados aportados por
real de los fludos puros y mezclas. Por ello, esta seccin resulta imprescindible para
situar este trabajo de investigacin y justificar su contenido. Una vez realizada la co-
rrespondientes bsqueda sobre los sistemas binarios elegidos y similares, con datos de
tante comparar los resultados que se presentan para sistemas idnticos y, sobre todo,
determinar nuevos datos para aquellas mezclas que sean necesarias, lo que permite un
formadas por el metanol con alcanos y con los primeros steres de alquilo, que son los
utilizados en este trabajo de Tesis Doctoral. A continuacin se realiza una breve pre-
sentacin de cada uno de los artculos de la tabla, escogidos por orden cronolgico de
aparicin en la literatura y diferenciando segn las mezclas sean de metanol con: a), n-
-12-
Tabla 1-1. Esquema-resumen de la recopilacin bibliogrfica realizada para las mezclas metanol+alcanos,
metanol+steres de alquilo, tanto en condiciones isobricas como isotrmicas.
-13-
1.- Akita y Yoshida (1963). 27.- Nagata y col (1972).
2.- Balashov y col. (1967) (ver Gmehling y Onken, 1977). 28.- Nagata y col (1975) (ver Gmehling y Onken, 1977).
3.- Bekarek (1968) (ver Gmehling y Onken, 1977). 29.- Nagata y col (1976).
4.- Beregovykh y col. (1971) (ver Gmehling y col., 1982a). 30.- Nagata y Ohta (1971).
5.- Bredig y Bayer (1927) (ver Gmehling y Onken, 1977). 31.- Nakanishi y col.(1967)(ver Gmehling y Onken, 1977).
6.- Budantseva y col. (1975) (ver Gmehling y col., 1988). 32.- Ohe y col (1971).
7.- Bushmakin y Kish (1975) (ver Gmehling y col.,1982a). 33.- Oracz y Warycha (1995).
8.- Choi y col. (1981) (ver Gmehling y col., 1988). 34.- Park y col (1973) (ver Gmehling y Onken, 1977).
9.- Crawford y col. (1949) (ver Gmehling y Onken, 1977). 35.- Patlasov y col. (1977) (ver Gmehling y col., 1988).
10.- Dobroserdov y Bagrov (1967) (ver Gmehling y Onken 36.- Polak y Lu (1972b).
1977). 37.- Raal y col. (1972).
11.- Ferguson (1932). 38.- Rajendran y col. (1988).
12.- Figurski y Malanowski (1996). 39.- Rhim y Kwak (1981) (ver Gmehling y col., 1988).
13.- Figurski y von Weber (1976) (ver Gmehling y col., 40.- Scaramucci y Vangell (1984).
1982a). 41.- Serafimov y col. (1966) (ver Gmehling y col., 1988).
14.- Figurski y von Weber (1982). 42.- Severns y col. (1955).
15.- Goral y col. (1983) (ver Gmehling y col., 1988). 43.- Scheller y col (1969) (ver Gmehling y col., 1982a).
16.- Hongo y col. (1994). 44.- Sieg y col (1951) (ver Gmehling y col., 1982a).
17.- Hwang y Robinson (1977)(ver Gmehling y col. 1982a) 45.- Tenn y Missen (1963).
18.- Iguchi (1978) (ver Gmehling y col., 1982a). 46.- Thomas y col. (1991)
19.- Li y col. (1991). 47.- Van Zandijcke y Verhoeye (1974) (ver Gmehling y
20.- Litvinov (1952) (ver Gmehling y Onken, 1977). col.,1982a).
21.- Martn y col. (1994). 48.- Vasileva y col. (1983) (ver Gmehling y col., 1988).
22.- Mato y Cepeda (1984). 49.- Vilim (1961) (ver Gmehling y Onken, 1977).
23.- Murti y van Winkle (1958). 50.- Wilsak y col. (1987).
24.- Nagahama e Hirata (1971) (ver Gmehling y Onken, 51.- Wobst y col (1992).
1977). 52.- Wolff y Hppel (1968) (ver Gmehling y Onken,
25.- Nagata (1962) (ver Gmehling y Onken, 1977). 1977).
26.- Nagata (1969) (ver Gmehling y Onken, 1977). 53.- Zawisza (1985).
Vilim en 1961, ver Gmehling y Onken (1977), obtuvo datos de ELV isobrico de
En 1963 aparecen los datos del ELV isobrico metanol+n-pentano a 99,99 kPa,
Por su parte, Wolff y Hppel (1968), ver Gmehling y Onken (1977), presentan
una amplia relacin de datos experimentales del equilibrio lquido-vapor isotrmico del
Raal y col. (1972) aportaron datos del ELV isobrico, a 101,32 kPa, del sistema
-14-
metanol+n-hexano.
(1988), obtienen datos p-T-x-y del ELV isobrico de los sistemas lquidos binarios meta-
Hwang y Robinson (1977), ver Gmehling y col. (1982a), estudian el ELV del
Al ao siguiente, Iguchi (1978), ver Gmehling y col. (1982a), presenta los datos
Choi y col. (1981), ver Gmehling y col. (1988), determinan los datos correspon-
En ese mismo ao, Rhim y Kwak (1981), ver Gmehling y col. (1988), realizan
Goral y col. (1983), ver Gmehling y col. (1988), obtuvieron para metanol+n-
Zawisza (1985) aporta datos del ELV de este mismo sistema a temperaturas su-
tanol+n-hexano.
-24-
binaria metanol+n-hexano.
Bredig y Bayer en 1927, ver Gmehling y Onken (1977), presentaron los ELV del
Diecisiete aos ms tarde aparecen los datos del ELV isobrico a presin atmos-
frica normal, 101,32 kPa, del sistema formado por metanol+etanotato de metilo, en
En 1951, Sieg y col., ver Gmehling y col. (1982a), determinaron a varias tempe-
Por su parte Litvinov (1952), ver Gmehling y Onken (1977), present el ELV a
etilo, tanto en condiciones isobricas, 101,32 kPa, como isotrmicas, 313,15, 323,15 y
333,15 K.
Nagata en 1962, ver Gmehling y Onken (1977), publica nuevamente datos del
-25-
Al ao siguiente, este mismo sistema y en las mismas condiciones fue tambin
isobricas de baja presin, 26,66 kPa, lo realizan Serafimov y Balashov (1966), ver
ven a presentar el ELV de este sistema binario a la presin indicada antes, a presin at-
De nuevo en 1967, Nakanishi y col., ver Gmehling y Onken (1977), aportan da-
tos de ELV del sistema metanol+etanoato de etilo a una presin cercana a la atmosfri-
Bekarek (1968), ver Gmehling y Onken (1977), vuelve a estudiar los ELV iso-
Nagahama e Hirata, ver Gmehling y Onken (1977), midieron en 1971 las mag-
nitudes caractersticas p-T-x-y del ELV del sistema metanol+etanoato de metilo a una
-26-
presin elevada, 585,63 kPa.
Tambin en 1971, Nagata y Ohta presentan, a 66,66 kPa, el ELV isobrico del
Ohe y col. (1971) aportan datos experimentales del ELV del sistema metanol
En 1972, Nagata y col. obtuvieron datos isotermos del ELV del sistema meta-
el estudio del ELV isobrico del sistema binario metanol+etanoato de etilo en condi-
Bushmakin y Kish (1975), ver Gmehling y col. (1982a) reportan el ELV isobri-
En el mismo ao, 1975, se presentan otros datos de ELV del sistema meta-
Figurski y von Weber (1976), ver Gmehling y col. (1982a), determinan para las
Por su parte, Nagata y col. (1976) estudian el ELV a 318,15 K para la mezcla
-27-
metanol+metanoato de etilo.
Cinco aos ms tarde, Figurski y von Weber (1982) realizan un nuevo estudio
Vasileva y col. (1983), ver Gmehling y col. (1988), presentan datos del ELV,
tanto isobrico, 101,32 kPa, como isotrmico, a 308,15 y 328,15 K, para la mezcla me-
tanol+etanoato de metilo.
En ese mismo ao, Scaramucci y Vangell (1984) estudiaron el ELV del mismo
Martn y col. (1994) publican datos del ELV isotrmico, 298,15 K, del sistema
metanol+etanoato de metilo.
clas en las que interviene el metanol como uno de los componentes, realizados en condi-
ciones isotrmicas. Por el contrario, son escasos los valores para mezclas en condiciones
-28-
isobricas y cuando existen estn realizados principalmente en condiciones atmosfricas
y subatmosfricas. Esto es una seal inequvoca que los cientficos han estado interesa-
metanol en solucin, ms que utilizar datos prcticos para diseo o estudio de procesos
industriales en los que intervenga el metanol. Unicamente se aportan datos del ELV a
metilo realizado por Balashov y col. (1967) a 395,41 kPa y por Nagahama e Hirata
son los que contienen etanoatos de alquilo, con la excepcin del propilo, lo que permite
1.2. Comparacin de datos bibliogrficos de ELV para las mezclas con metanol
rencias con datos experimentales de ELV de las mezclas propuestas para esta investiga-
cin a la misma presin de trabajo. No obstante, puesto que para algunos sistemas bina-
vo con el fin de verificar la calidad de los valores propuestos, eliminndose en esa com-
-29-
Iniciando la comparacin con los sistemas binarios metanol+n-alcanos (C5,C6),
x1metanol+x2n-pentano
Tenn y Missen (1963), 99,99 kPa
0 Budantseva y col. (1975), 101,32 kPa
(ver Gmehling y col.,1988)
y1-x1
x1metanol+x2n-hexano
-0.2 Vilim (1961), 99,32 kPa
(ver Gmehling y Onken, 1977)
Budantseva y col. (1975) 101,32 kPa
pentano (ver Gmehling y col., 1988)
Raal y col. (1972), 101,32 kPa
-0.4
x1metanol+x2n-alcanos
-0.6
0 0.2 0.4 0.6 0.8 1
x1
los de Tenn y Missen (1963) a 99,99 kPa y Budantseva y col. (1975), ver Gmehling y
col. (1988), a 101,32 kPa para la mezcla x1metanol+x2pentano y los datos de la mezcla
x1metanol+x2hexano aportados por Vilim (1961) a 99,32 kPa, Budantseva y col. (1975),
ver Gmehling y col. (1988) a 101,32 kPa y Raal y col. (1972) tambin a 101,32 kPa. La
diferencia media entre los valores de la fraccin molar de la fase vapor, y1, de los dos
-30-
se han representado las correlacines de los datos de Budantseva y col. (1975), ver
Gmeh-ling y col. (1988), para ambas mezclas binarias. No parecen estar justificadas las
bajo, entre los valores presentados por Vilim (1961), ver Gmehling y Onken (1977) a
punto a punto de Fredenslund y col. (1977b), donde se han incluido las modificaciones
diciones isobricas son los que contienen etanoatos de alquilo, con la excepcin del
etanoato de propilo. La Figura (1-2) representa la comparacin entre los valores de (y1-
yen las curvas del ajuste, utilizando la ecuacin polinmica de las zetas (ver aparta-
do 5.1.4.), de los datos de ELV a las presiones de 26,66 kPa de Balashov y col. (1967),
ver Gmeh-ling y Onken (1977), a 101,32 kPa de Mato y Cepeda (1984) y a 585,63 kPa
de Nagahama e Hirata (1971) (ver Gmehling y Onken, 1977). Las discrepancias entre
los datos de las distintas referencias, en idnticas condiciones de presin, 101,32 kPa,
fraccin molar de la fase vapor, y1, no super en ningn caso las 0,009 unidades. La
denslund y col. (19 77b), excepto los de Crawford y col. (1949), ver Gmehling y Onken
-31-
(1977) a 101,32 kPa y de Balashov y col. (1967), ver Gmehling y Onken (1977) a
mostrada en la Figura (1-3), si bien, en este caso, la bibliografa no recoge ningn estu-
presin de 13,33 kPa de Park y col. (1973), ver Gmehling y Onken (1977) y a 101,32
-32-
kPa de Nagata (1962), ver Gmehling y Onken (1977) y de van Zandijcke y Veroheye
(1974), ver Gmehling y col. (1982a). De los conjuntos de datos presentados, nicamente
verifican la consistencia los estudiados por Park y col (1973) a 66,66 kPa, ver Gmehling
y Onken (1977), y por Akita y Yoshida (1963), Murti y van Winkle (1958) y Nagata
(1962), ver Gmehling y Onken (1977) a 101,32 kPa. A esta presin, los datos que pre-
sentan mayores diferencias son los de van Zandijcke y Veroheye (1974), ver Gmehling y
col. (1982a), que frente a los de Nagata (1962), ver Gmehling y Onken (1977) son de
-24-
los trabajos de Beregovykh y col. (1971), ver Gmehling y col. (1982a) y Patlasov y col.
(1977), ver Gmehling y col. (1988), ambos realizados a presin atmosfrica normal. La
Como resumen de este anlisis se deduce que los datos de la bibliografa actual
0.8
Figura (1-4). Comparacin entre los
datos experimentales del equilibrio lqui-
do-vapor isobrico extrados de la biblio-
grafa, para el sistema x1metanol+x2 eta-
0.6
noato de butilo en trminos de la magni-
tud de concentraciones (y1-x1) vs x1. Am-
bos trabajos han sido realizados a la pre-
sin de 101,32 kPa.
y1-x1
x1metanol+x2etanoato de butilo
0
0 0.2 0.4 x1 0.6 0.8 1
mezclas que contienen metanol ya que, adems de existir pocos sistemas estudiados,
-25-
algunos de ellos presentan importantes diferencias. Por ello, esta investigacin puede
con el fin de ampliar la base de datos que pueda ser utilizada para la prediccin e inter-
pretacin correcta del equilibrio lquido-vapor isobrico de sistemas con una misma
naturaleza qumica.
establecerse una vez llevada a cabo la primera fase del mismo, o sea, ya realizada la
dicha seccin se desprende un enorme inters de los investigadores por abordar los estu-
dios de los equilibrios lquido-vapor isobrico de mezclas binarias que contienen meta-
nol, tanto con n-alcanos como con alcanoatos o steres de alquilo, a distintas presiones
sin de trabajo para todas las mezclas de 141,3 kPa, para las que no existen datos en la
literatura consultada y que, al mismo tiempo, supone un lmite cercano para la experi-
aqu y en las secciones anteriores, el conjunto de sistemas binarios elegidos fue el for-
mado por las mezclas del metanol y n-alcanos (CnH2n+2 con n=5,6) o steres de alquilo,
excepcin del metanoato de metilo. Con los valores aportados, no slo se conseguira
-26-
una mejor base de datos que la existente en la actualidad sino tambin llegar a tener ms
con el fin de proponer la modelizacin ms adecuada. Los pasos a seguir que definen el
zados en los estudios del equilibrio lquido-vapor isobrico y que, en general, apoyan su
nes citadas en la literatura, como son aquellas correspondientes a los modelos de com-
ras con el fin de ajustar y reproducir las propiedades mas significativas de las mezclas,
como son el volumen y la entalpa de exceso, vE y hE, as como las magnitudes derivadas
del equilibrio lquido-vapor, la temperatura, T, la energa libre de Gibbs, gE y/o los coe-
los coeficientes de dicha ecuacin. Esta variacin permitir darle una mayor extensin
-27-
en su utilizacin, haciendo posible con ello la estimacin de otras magnitudes, como es
la entalpa, a partir de la correlacin de la funcin de Gibbs obtenida con los datos del
equilibrio isobrico.
del mismo modelo realizada por Weidlich y Gmehling (1987), ver tambin Gmehling y
col. (1993), para la prediccin del equilibrio de fases de estos sistemas. En el modelo
especfico, sin embargo, no sucede lo mismo en el modelo ASOG, donde ha de ser trata-
do como un alcanol ms, a travs del grupo caracterstico (-OH). Aplicando ambos mo-
delos, se compararn los valores predichos para la concentraciones de la fase vapor, y1,
Con la breve descripcin de la labor aqu planteada y los objetivos que se pre-
tenden, parece justificado el contenido de esta Memoria sobre el trabajo realizado, bajo
el ttulo de: Anlisis del Equilibrio Lquido-Vapor a 141,3 kPa de Mezclas Binarias
que contienen Metanol con n-Alcanos (C5, C6) y con Esteres Alqulicos.
0 0 0
-28-
I. PARTE TEORICA
3.1. Introduccin
fases, cada una de las cuales constituye un todo homogneo sin tendencia a cambiar sus
propiedades con el tiempo. Existen dos aspectos muy concretos que se deben destacar,
primero, resulta evidente que la temperatura de las fases debe ser la misma, de lo contra-
rio existira un flujo irreversible de energa calorfica entre ellas y segundo, tambin debe
haber una igualdad de presiones ya que, si no, la inestabilidad ocasionara una compre-
sin de una de las fases y la expansin de la otra con un intercambio irreversible de ener-
ga. Sin embargo, estas dos condiciones no son suficientes para asegurar el equilibrio de
fases.
potenciales mencionadas, llegando a anularse al alcanzar dicho estado. De entre las mag-
nitudes citadas anteriormente, la funcin de Gibbs tiene una atencin especial para el
tratamiento matemtico del equilibrio entre fases. Por ello, en un sistema multifsico la
31
dG=(-SdT+Vdp) 0 (3-1)
del cambio que conduce al equilibrio. As, todos los procesos irreversibles que transcu-
cin de Gibbs.
mediante,
infinitesimales que experimentan las variables de las que depende. As, de acuerdo con
32
G G G
dG = dp + dT + dni (3-3)
p T , ni T p , ni i ni p ,T , n
j i
(3-1), las dos primeras derivadas parciales pueden ser sustituidas, dando lugar a otra ex-
presin anloga para la funcin de Gibbs. Por otro lado, la derivada parcial del ltimo
i, y que, para el caso que nos ocupa, coincide con la propiedad parcial de la correspon-
pero, por otro lado, como: G=xi gi , y d(ng)= gi dni+nid gi , identificando esta ex-
aplicada en este caso a la funcin de Gibbs, pero que puede ser utilizada para cualquier
33
de presin y temperatura constantes en las que resulta muy simplificada:
o sea que, en esas condiciones, las nicas variaciones en el potencial qumico o en la pro-
piedad molar parcial son las debidas a los cambios de composicin. Adems, la ecuacin
(3-6) indica la relacin existente entre las diferentes gi de los compuestos. Otra impor-
tante conclusin de este anlisis, con la funcin del potencial qumico, surge al plantear
la ecuacin (3-3) para un sistema de fases en equilibrio, con T y p constantes, la cual lle-
idni=0 (3-7)
si, como en el caso que nos ocupa, se particulariza esta ecuacin a las fases vapor y lqui-
pero, en la transferencia de masa entre las fases, la ley de conservacin de la materia re-
( iv - il )d niv =0 (3-9)
como los d niv son independientes y toman valores arbitrarios, para que se satisfaga la
34
iv = il (i=1,2,3,...n) (3-10)
equilibrio con respecto a la transferencia de materia entre fases. O sea que, adems de la
nicas variables que se suponen idnticas en todo el sistema, puede considerarse una sola
i ( T , p ) = i ( T , p ) (3-11)
para sentar las bases termodinmicas del tratamiento de sistemas en equilibrio. Sin em-
bargo, el potencial qumico no es una propiedad prctica bajo el punto de vista matemti-
co y por otro lado, en los estudios de fluidos resulta de gran utilidad comparar las propie-
dades de un sistema real con las de otro ideal. Por ello, es conveniente introducir nuevas
35
tes, que sern la conexin entre el mundo abstracto de las Matemticas y de la Termodi-
dGideal=Vdp=RTdlnp (3-12)
de la misma, igualmente vlida, sera sustituir la presin por una nueva funcin denomi-
dGi=RTdlnfi (3-13)
Al identificar las expresiones (3-12) y (3-13) surge una ecuacin diferencial, cuya
do de fluido ideal, o sea, asignar que en p=0 (fi/p)=1 o bien, ln(fi/p)=0. De esta forma
se establece que:
fluido puro. El cociente adimensional (fi/p) que aparece en esta definicin recibe un nom-
36
bre especial, el coeficiente de fugacidad, i, del componente puro i. Estos conceptos, de-
finidos para un componente puro en un fase, pueden ser extendidos a una solucin y a un
compuesto en solucin. As, podrn escribirse los dos conjuntos de relaciones siguientes:
d gi =RTdln f!i (a T constante), lim ( f!i /xi p)=1, ! i = f!i /xi p (3-16)
p0
cuenta las relaciones establecidas hasta ahora, puede utilizarse fcilmente la fugacidad
como criterio de equilibrio de fases, obtenindose de esta forma una relacin similar a la
varias fases debe ser la misma en todas las fases en que est presente en equilibrio. La
37
3.3.2. Actividad. Coeficiente de actividad.
do expresarse por:
a i = ( fi / f i o ) (3-18)
dall (1961), definieron como mezclas ideales aquellas que cumplen la relacin:
fi = f i o xi (3-19)
Segn lo cual, en una mezcla ideal, ai =xi. Por lo tanto, la desviacin de la ideali-
a i f
i = = oi (3-20)
xi f i xi
con un sentido prctico, como tercera condicin, sera mejor utilizar la igualdad de fuga-
38
cidades parciales de los componente individuales del sistema. Es decir, la ecuacin (3-
17) aplicada al caso del equilibrio entre un vapor y un lquido para un componente i, pue-
de escribirse como:
fi v = fi l (3-21)
su vez, pueden combinarse con las correspondientes ecuaciones de estado. Pero cuando
esto no es posible para la fase lquida entonces se introducen los coeficientes de activi-
dad. Si, como es usual en el tratamiento del equilibrio lquido-vapor las concentraciones
ra, a partir de ahora omitiremos los superndices indicadores de las fases, escribiendo
1 RT
p
i = exp vi dp (3-23)
RT 0 p
39
La fugacidad de referencia, fio o fi, en los estudios de ELV a presiones moderadas,
se considera en la mayor parte de los casos, como la del componente puro a la presin y
temperatura de la mezcla. Este valor de fugacidad, fi, est relacionado con la presin de
(v )dp
1 p
f i = pio io exp o
i (3-24)
RT pio
rreccin de Pointing (FP), que tiene en consideracin el hecho de que el lquido se en-
temperatura de la mezcla. Dicho de otra manera, dicho factor establece la correccin rea-
lizada sobre la fugacidad de la fase lquida desde la presin de vapor hasta la presin del
puede escribirse,
1 ( )
vio p pio
p
FP = exp vi
o
dp = exp (3-25)
RT p io RT
Por ahora, la ecuacin establecida en los comienzos, la (3-22), queda como sigue,
cin, io, corrige las desviaciones del comportamiento del vapor saturado con respecto al
40
ro tomando como referencia la presin de saturacin.
1 o RT
pio
i = exp
o
vi dp (3-27)
RT 0 p
sicin de la solucin:
( )exp o RT
yi p vio p pio 1 p
i = exp vi dp (3-28)
xi pio RT RT p io p
pio Bii
io = exp (3-29)
RT
i =
yi p (
Bii vio p pio
exp
)( )
y 2j p ij
exp (3-30)
xi pio RT RT
siendo ij = 2Bij - Bii -Bjj. Ahora, para ultimar el clculo de los i con esta expresin se
41
necesitan valores de los volmenes de compuestos puros en estado lquido saturado as
como la determinacin de los coeficientes de virial, tanto de los compuestos puros como
de las mezclas. En la bibliografa existen varios mtodos para determinar dichas propie-
dades, sin embargo, las secciones siguientes estn dirigidas a presentar nicamente aque-
Por su sencillez y sobre todo por los excelentes resultados que ofrece en nuestros
RT [ + ( )2 /7 ]
vio = c Z RA 1 1 TR (3-31)
pc
precisan datos de las constantes crticas del compuesto y el valor de ZRA, que es una
mentales de densidad a varias temperaturas. Para los compuestos utilizados en este tra-
bajo, steres de alquilo, alcanos y metanol, en la Tabla (A-1) de los apndices se reco-
42
3.4.2. Determinacin de los coeficientes de virial. Correlaciones de Tsonopoulos
Para el clculo de los segundos coeficientes de virial, tanto de las sustancias pu-
ras, Bii, como de las mezclas, Bij, se opt por las estimaciones propuestas por Tsonopou-
los (1974), que constituyen una modificacin del mtodo propuesto por Pitzer y Curl
Bpc
= To (TR ) + T1 (TR ) + T2 (TR ) (3-32)
RTc
y donde, a su vez, cada uno de los sumandos de (3-32) viene dado por las expresiones
empricas:
a b
T2 (T R ) = (3-35)
T R6 T R8
gases polares, segn presenten o no enlaces por puentes de hidrgeno. Cuando esto suce-
de, surgen asociaciones moleculares que originan una mayor complejidad en la influencia
43
De forma simplificada, Tsonopoulos encontr que, para fluidos polares que no posean
10 5 p c2
R = (3-36)
Tc2
Para fluidos con enlaces de hidrgeno, el valor de b ser distinto de cero; por
mricos de a y b para cada uno de los compuestos utilizados en este trabajo. Asimismo,
cientes de virial de mezclas, Bij. En este caso, las magnitudes reducidas se calculan con-
pseudocrticas.
siendo kij una constante caracterstica de cada mezcla binaria, independiente de la tempe-
ratura y cuyo valor es pequeo. En todos los sistemas binarios estudiados en este trabajo
se ha considerado kij=0,01.
RTc ij z c ij
p c ij = z cij = 0.5(z ci + z c j ) y vcij = (1/8)(vc1/3
i
+ vc1/3
j
)3 (3-38)
v c ij
ij = 0.5( i + j ) (3-39)
44
Las constantes aij y bij de la mezcla para el ltimo trmino de la expresin de Tso-
productos puros. Si la mezcla binaria es del tipo polar/no-polar sus valores son nulos.
clculo de las variables involucradas, sino tambin permite el anlisis de los datos con el
fin de valorar su calidad. Las medidas experimentales de equilibrio estn sujetas a erro-
centracin de uno de los componentes. Sin embargo, es necesario efectuar un anlisis ri-
guroso, con las relaciones que plantea la Termodinmica, para demostrar la existencia o
de los equipos o el anlisis de las muestras. Dicho anlisis termodinmico se conoce co-
(T,xi,yi,i) para los equilibrios isobricos y de (p,xi,yi,i) para los equilibrios isotrmicos.
45
cin, aplicando la ecuacin de Gibbs-Duhem a las magnitudes derivadas, i. Dicha ecua-
tiene la forma:
v h
x dln = RT dp RT
i i 2
dT (4-1)
entonces, la expresin final de la ecuacin anterior, la (4-1), aplicada a una mezcla bina-
1 dg E v dp h dT 1
=
RT dx1 RT dx - RT 2 dx + ln (4-3)
1 1 2
que, integrada entre los lmites x1=0 a x1=1, donde la correspondiente funcin de exceso
se anula, produce
T1o p1o
h v
1
ln 1 dx = dT dp (4-4)
0 2 T2o RT 2 p 2o RT
La integral del primer miembro representa la suma algebrica de las reas positiva
cin del componente 1 en la fase lquida, x1, ecuacin de Redllich-Kister (1948), la se-
46
co y el isobrico, la aplicacin de (4-4) da lugar a las siguientes consideraciones. En el
caso isotrmico, la ecuacin (4-4) pierde la primera integral del segundo miembro y en-
reas sea igual al valor de la otra integral. Sin embargo, no es usual encontrar datos que
que nos ocupa en este trabajo de investigacin, que es al mismo tiempo el ms comn en
T1o
h
1
ln 1 dx = dT (4-5)
0 2 T2o RT 2
surge cuando nos lo hay o no son fiables. Para esos casos, Herington (1951) estudi la
forma de las curvas h(x,T) y propuso una aproximacin de dicha integral que requiere
son:
47
donde Tmax es la mayor temperatura de ebullicin presente y Tmin la mnima. De acuerdo
con el criterio de Herington (1951), las condiciones de calidad de los datos de ELV son,
D< J, condicin suficiente para la consistencia; D-J< 10, para la consistencia termo-
tuales en el planteamiento de dicho test. No obstante, dichos mtodos con una valoracin
de las reas, suponen una verificacin global de los datos, no individual. La comproba-
cin individual de cada punto experimental viene dada hoy dia por una serie de clculos
que definen unos procedimientos conocidos como test punto-a-punto. Sin embargo, to-
dos los existentes presentan serias deficiencias pues, en general, no tienen en cuenta los
un test ms coherente que combina ambos, el test punto a punto y el de las reas con uno
nuevo basado en la representacin de (gE/RTx1x2) vs. x1. La ventaja del mtodo es que se
pueden seleccionar los datos de confianza entre varios conjuntos de datos y no se requiere
Wisniak (1993) propone un test punto a punto utilizando relaciones sencillas, co-
Li = T o
k kx sko / s T = g E / s RTw / s = Wi (4-7)
necesidades para las cuales se requieran los datos. El test propuesto por Wisniak (1993)
48
es ms elaborado que el de Herington (1951) pero presenta una mejor base terica y es
ms til que el test propuesto por Kojima (1990) porque no requiere informacin de las
nmica individual de cada punto experimental de ELV, sigue siendo el desarrollado por
Van Ness y col. (1973) en la versin posterior propuesta por Fredenslund y col. (1977b) y
basado en el clculo de las fracciones molares de la fase vapor, yi, a partir de los datos
(T,xi) experimentales, comparndose los mismos con los obtenidos por experimentacin a
travs de las diferencias, yi=yi,exp-yi,cal. Los datos sern consistentes si se verifica que,
yi< 0,01. Este criterio, si bien en su utilizacin rigurosa sera aplicado de forma indivi-
dual a cada uno de los estados de equilibrio, a veces se emplea para valorar globalmente
conjunto de valores.
nacin del ELV, se calcula la cuarta y se compara su valor con el valor medido. As, la
presin total en fase vapor puede ponerse como p=yip. Teniendo en cuenta las expre-
o v o ( p pio )
p= xi pio i exp ln i + i
I
RT
(4-8)
cin (4-2) y luego particularizando las determinaciones y los clculos para una mezcla
49
binaria. De esta forma se plantean las expresiones individuales,
1 1
ln 1 = Q + x2 ln y ln 2 = Q x1 ln (4-9)
2 2
estricta se requiere informacin sobre las entalpas de mezcla, ver ecuacin (4-4), o sea
dQ h dT dQ h dT
ln 1 = Q + x2 + 2
ln 2 = Q x1 + 2
(4-10)
dx1 RT dx1 dx1 RT dx1
nes marcadas por el ELV isobrico a las presiones de trabajo pero, sobre todo, a la tempe-
ratura de ebullicin o cercanas a ella, por tanto, se utiliza un tratamiento alternativo segn
se indica ahora y aplicable al caso que nos ocupa en este trabajo. En condiciones isobri-
presa por:
dQ Q Q dT
= + (4-11)
dx1 x1 T , p T x1 , p dx1
si se identifica esta expresin con la (4-3), se tiene una expresin, que resulta trascen-
50
Q h
= 2 (4-12)
T x1 , p RT
combinando ahora esta relacin con las presentadas anteriormente en las (4-10), se consi-
Q Q
ln 1 = Q + x2 ln 2 = Q x1 (4-13)
x1 p ,T x1 p ,T
sistemas binarios, sin embargo, puede tener el significado matemtico indicado cuando se
Para el asunto que nos ocupa en esta seccin, con la verificacin termodinmica
de los datos de ELV, la representacin de Q=Q(xi) se efecta con los polinomios de Le-
gendre, obteniendo una ptima correlacin de los datos experimentales cuando se mini-
mizan las diferencias de las presiones experimentales y las calculadas por la ecuacin (4-
9). Por ltimo, se determinan las concentraciones de uno de los componentes en la fase
En este trabajo, en la aplicacin del test punto-a-punto con la versin dictada por
Fredenslund y col. (1977b), se introducirn algunas variantes. As, los coeficientes del vi-
51
rial de la correspondiente ecuacin de estado, se calculan por las expresiones empricas
ratura de los volmenes molares de los lquidos puros saturados, con la ecuacin de
Muchas ecuaciones han sido propuestas para correlacionar los coeficientes de ac-
tividad con la composicin y con la temperatura, algunas de ellas con un cierto funda-
mento terico y otras, puramente empricas pero creadas con cierta intuicin. En la ac-
tualidad existen varias clases de correlaciones para los coeficientes de actividad, pero
ninguna de ellas presenta una superioridad absoluta sobre las dems. Una comparacin de
(5-1) se presenta un resumen clarificador de las ms utilizadas por los investigadores, las
cuales sern seguidamente presentadas con cierto detalle. Estrctamente, los coeficientes
dad, a travs de la ecuacin (4-2) que, para una mezcla binaria resulta:
gE
= x1 ln 1 + x 2 ln 2 (5-1)
RT
52
Figura (5-1).Tres de las formas que adoptan las funciones Q=gE/RT y lni vs x1 para sistemas binarios.
cionar los valores de la funcin de exceso de energas de Gibbs tienen normalmente for-
mas empricas que, para una mezcla binaria se expresan genricamente como: gE/RT=
x1x2(x1). Las formas funcionales empleadas pueden estar basadas en algn modelo te-
rico de interaccin entre molculas o puede ser puramente emprico. Las ms antiguas
son la propuesta en 1985 por Margules y en 1910 por Van Laar y la ms popular para co-
tenta que las expresiones para gE o para gE/RT sean lo menos complejas posible, pero con
les, como se ilustran en la Figura (5-1). Tambin conviene que dichas expresiones po-
53
sean el fundamento fsicoqumico apropiado, de tal forma que los parmetros que inter-
caso ser necesario elegir aquellas que mejor se adapten a una situacin concreta. Sin
embargo, despus de una revisin general, las ms importantes pueden ser incluidas en
dos grupos: 1) las basadas en desarrollos en serie de potencias, y 2) las basadas en mo-
La forma ms sencilla de generar una expresin para reproducir las gE/RT vs con-
centracin en una mezcla binaria es con un desarrollo en serie de potencias. Son varias
forma:
54
Q = x1 x2 [(A21 x1 + A12 x2 )(D21 x1 + D12 x2 )x1 x2 +"] (5-2)
donde, en principio, los parmetros A21, A12, D21, D12 ..., son funcin de la temperatura y
la presin. Sin embargo, rara vez se utiliza esta expresin con ms de cuatro parmetros y
[ ( ) ]
ln i = x j2 Aij + 2 Aji Aij xi (5-4)
1 1 A12'
x1 + A21
'
x2 D12
'
x1 + D21
'
x2
= x x + " (5-5)
Q x1 x2 A12 ' ' D' D' 1 2
A21 12 21
donde, tambin aqu, los parmetros A21' , A12' , D21' , D12' , ... son funcin de la presin y de
gE A' A' x x
= ' 12 21 1' 2 (5-6)
RT A12 x1 + A21 x 2
55
y, por lo tanto, los coeficientes de actividad vendrn dados genricamente por:
ln i =
( ) 2
Aij' Aj' i x j2
(5-7)
(A x + A x )
'
ij i
'
ji j
2
[
Q = x1 x2 A + B(x1 x2 )+ C (x1 x2 )2 + D(x1 x2 )3 + " ] (5-8)
gresin de mnimos cuadrados; el nmero de ellos que deba utilizarse depende, a criterio
termodinmica que se estudia y del nmero de puntos tratados. Los coeficientes de acti-
[ ( ) ( ) ] [(
ln i = xi x j A + B(1)i + 1 xi x j + C xi x j 2 + " + x j A x j xi + )
(5-9)
+ B (1) i +1
(6 xi x j 1)+C (xi x j )(8 xi x j 1)+"]
foma adecuada para representar tanto gE/RT como otras magnitudes de exceso caracters-
ticas de una mezcla lquida, como son las entalpas de exceso y volmenes molares de
exceso.
56
5.1.4. Ecuacin de las zetas. Propuesta de modificacin
Surge de una simplificacin, propuesta por Ortega y Alcalde (1992), sobre la ex-
pansin en serie de Wohl (1946) para gE/RT en funcin de las fracciones volumtricas
[
Q = x1 x2 A0 + A1 z + A2 z 2 + A3 z 3 +" ] (5-10)
de los datos (x1,Q), mientras que el valor de k se optimiza para conseguir la mejor corre-
nal deducida y que tiene una forma anloga: Q=z(1-z) Ai zi, que tambin da excelentes
2
ln i = x j2 (A + A z + A z
0 1 2
2
)
+" + xi (1) i +1
(A1 + 2 A2 z +")k z (5-11)
x1
d( g E / T ) hE
= 2 (5-12)
dT T
Si, por simplificar se considera que la entalpa de exceso presenta una variacin
57
lineal con la temperatura, del tipo: hE=+T, entonces, de acuerdo con la (5-12).
gE/T=/T-lnT+ (5-13)
donde es una constante por determinar. Esta consideracin significa que la correspon-
ratura, cosa que, en la literatura actual no se consigue fcilmente y menos an, para las
mezclas que se estudian en este trabajo. Por ello, la modificacin ms simple propuesta
masiado estricta pero permite unos clculos ms sencillos cuando no se dispone de datos
sern los errores provocados por dicha consideracin. O sea que, si ahora hacemos que,
Q gE
= = A0 + A1 z (k ) (5-14)
x1 x2 x1 x2 RT
con la ecuacin (5-12), tambin puede escribirse mediante la siguiente expresin, que se-
58
hE d ( g E / x1 x2 RT )
= T (5-16)
x1 x2 RT dT p, x
hE dA dA z dk
= T 0 T 1 z TA1 (5-17)
x1 x2 RT dT x dT x k dT x
la variacin del parmetro k con la temperatura puede despreciarse, ya que dicha variable
(5-17) se reduce entonces a los dos primeros sumandos; as, dicha expresin pueden ser
hE
= A10 + A11 z (5-18)
x1 x2 RT
dA dA
donde ahora: A01 = T 0 A11 = T 1 (5-19)
dT dT
A01 A11
A01 = y A11 = (5-20)
T T
59
entalpas, los Ai2 se calculan correlacionando los datos de ELV.
Hasta aqu se han mostrado una serie de ecuaciones, basadas todas ellas en desa-
individualmente, es distinto del que resulta si se tiene en cuenta toda la solucin. As,
aparece el concepto de fraccin molar local que se relaciona con la fraccin molar total a
Huggins para soluciones de polmeros. Concretamente, para una mezcla binaria dicha
60
g E / RT = x1 ln( x1 + 12 x 2 ) x 2 ln( x 2 + 21 x1 ) (5-21)
vi ij
donde el parmetro de interaccin ij, es: ij = exp
(5-22)
vj RT
de modo que, ij=ji y ijji. Los vi, vj son los volmenes molares de los componentes
da lugar a la siguiente expresin genrica para los coeficientes de actividad, con el mo-
delo de Wilson:
ij ji
ln i = ln ( xi + ij x j ) + x j (5-23)
x + x x j + ji xi
i ij j
grandes desviaciones positivas respecto de la idealidad, como por ejemplo, los sistemas
inmiscibilidad del ELL. Por vez primera, con la ecuacin de Wilson se logra, describir la
61
NRTL (Non Random Two Liquids), debido a que considera que la composicin interna
del lquido no es aleatoria, sino que las molculas estn semiarregladas. Las energas de
Gibbs de interaccin entre molculas se identifican por gij con un significado fsico simi-
lar al de los ij en Wilson. La ventaja de este modelo NRTL es que puede ser aplicado a
G G g ij
g E / RT = x1 x 2 21 21 + 12 12 siendo: ij = y Gij = exp ( ij ij ) (5-24)
x1 + G21 x 2 x 2 + G12 x1 RT
G ji
2
ij Gij
ln i = x j i
2 + (5-26)
x i + G ji x j +
j 2
( x G x )
j ij i
En el trabajo original, los autores sugieren que ij oscile entre 0,2 y 0,47, sin em-
bargo, en la prctica la cosa es algo distinta, incluso, se han encontrado valores negativos
para mezclas parcial o totalmente inmiscibles. Para el caso que nos ocupa en este trabajo
comprendidos entre 0,3 y 0,5, sin embargo, con el fin de optimizar la representacin de
62
La ecuacin NRTL no ofrece ventajas en sistemas moderadamente no-ideales, so-
bre las otras de Van Laar y Margules, pero realiza una buena representacin de los datos
(1975), propusieron una ecuacin con dos parmetros que, en cierto sentido, es la exten-
(comb) que intenta describir la contribucin entrpica y otra residual (resid) que describe
gE=gE(combinatorial)+gE(residual) (5-27)
las molculas. El trmino residual considera la energa de interaccin a travs de los pa-
g E ( comb ) z
= x1ln 1 + x2ln 2 + q1 x1ln 1 + q2 x2ln 2 (5-28)
RT x1 x2 2 1 2
63
g E ( resid )
= q1 x1ln( 1 + 2 21 ) q2 x2ln( 2 + 1 12 ) (5-29)
RT
El rea de contacto de una molcula i viene dado por el parmetro qi. En esta
y qi a travs de:
xi ri xi qi
i = i =
(
x i ri + x j r j ) y
(
xi q i + x j q j ) (5-30)
uij
ij = exp (5-31)
RT
i z r
ln i = ln + qiln i + j li i l j qiln( i + j j i ) +
xi 2 i rj
(5-32)
ji ij
+ j qi
+ +
i j ji j i ij
z
y donde, li = (ri qi ) (ri 1) . Los parmetros ri y qi se pueden estimar, para cada
2
64
componente, mediante un mtodo de contribucin de grupos, siendo:
ri = n j R j y qi = n j Q j (5-33)
pectivamente, a partir de los correspondientes valores Vwi, Awi de un grupo del tipo de
Vwi Awi
Rj = y Qj = (5-34)
15 ,17 2 ,5 10 9
ca, siendo a veces peor la representacin de los datos con este mtodo que con otras
tribucin de grupos UNIFAC que permite estimar los coeficientes de actividad de ELV.
65
Tabla (5-1). Ecuaciones propuestas por los distintos modelos para la correlacin de Q=gE/RT y ln i.
MODELO Q=gE/RT ln i
Margles
Q = x1 x2 ( A21 x1 + A12 x2 ) [ (
(5-3) ln i = x 2j Aij + 2 A ji Aij x i ) ] (5-4)
van Laar Q=
'
A12 A'21 x1 x 2
(5-6) ln i =
( )x
Aij' A'j i
2 2
j
(A x + A x )
(5-5)
'
A12 x1 + A21
'
x2 ' ' 2
ij i ji j
Redlich- [
Q = x1 x 2 A + B (x1 x 2 )+ C (x1 x 2 )2 + D(x1 x 2 )3 +!] (5-8) ln = x x [A+ B( 1 ) (x x )]+ x [A(x
i i j
i +1
i j j j ) (
x i + B( 1 ) i +1 6 x i x j 1 )] (5-9)
Kister
Q = x x [A + A z + A z ]
Ecuacin 2
+ A3 z 3 +!
( ) z
2
ln i = x 2j A0 + A1 z + A2 z 2 +! + x i ( 1 )i +1 ( A1 + 2 A2 z +!)k
(5-10)
de las 1 2 0 1 2
(5-11)
zetas z = x1 / (x1 + kx2) Ai=Ai1/T+Ai2 x1
Q = x1 ln( x1 + 12 x 2 ) x 2 ln( x 2 + 21 x1 ) (5-21)
ij ji
Wilson ln i = ln ( x i + ij x j ) + x j (5-23)
v ij x i + ij x j x j + ji x i
ij = i exp
vj RT
G G
Q = x1 x 2 21 21 + 12 12 (5-24) 2
1
x + G x x 2 + G12 x 1 G ji ij G ij
ln i = x 2j j i +
21 2
NRTL (5-26)
g ij x + G x ( x +G x )2
ij = y Gij = exp ( ij ij ) i ji j j ij i
RT
i z
Q = Q (combinatorial) + Q (residual) ln i = ln
xi 2 i
r
+ q i ln i + j l j i l j q i ln i + j j i +
r j
( )
1 2 z 1 2
Q( c ) = x1 ln + x 2 ln + q 1 x1 ln + q 2 x 2 ln (5-28) (5-32)
x 2 2 1 2 ji ij
+ j q i
x1
UNIQUAC
i + j ji j + i ij
Q( r ) = q 1 x 1 ln( 1 + 2 21 ) q 2 x 2 ln( 2 + 1 12 ) (5-29)
u ij
l i = (ri q i ) (ri 1)
x i ri xi qi z
i = ; i = ; ij = exp
(
x i ri + x j r j ) (
xi qi + x j q j ) RT
2
-66-
industriales. Desgraciadamente, no siempre es posible contar con buenas bases de datos,
tir de datos existentes, o bien estimar las propiedades apoyndose en una correlacin ade-
cuada. A lo largo de estos ltimos aos se han propuesto varios modelos tericos para
estimar las propiedades de los sistemas multicomponentes, que normalmente, utilizan in-
formacin sobre los compuestos puros y consideran, de manera simple, las interacciones
grupos funcionales (grupos de tomos que no necesariamente coinciden con los grupos
la solucin es la suma de las contribuciones debidas a cada grupo funcional. Estos mto-
otros grupos que forman la molcula, sobre uno dado, es decir, consideran que un deter-
minado grupo funcional contribuye de la misma forma en molculas distintas, sin tener
en cuenta su entorno. Todas estas ideas fueron reflejadas en el trabajo de Wilson y Deal
exceso. Dicha descripcin se realiz en base a los fundamentos que se plantean a conti-
nuacin y que constituyen la parte comn a los modelos, de los que surgen las variantes
67
1. Los coeficientes de actividad se pueden expresar como producto de dos contri-
buciones, una combinatorial, iC, debida a las diferencias de forma y tamao entre las
molculas presentes en la mezcla y otra residual, iR, asociada a las interacciones entre los
lni=lniC+lniR (6-1)
les de cada grupo presente en la solucin, menos la suma de las contribuciones indivi-
duales de los mismos grupos en un estado de referencia formado nicamente por las mo-
k xi
i
Xk = i
[i=1,2,...M(componentes);j=1, 2, ...N (grupos)] (6-3)
j xi
i
i j
68
Basados en estas consideraciones, Wilson y Deal (1962) determinaron los lnk en
ciertas mezclas alcano+alcanol, para los grupos hidroxilo y metileno, procediendo luego
a calcular los coeficientes de actividad para las mezclas consideradas. Wilson y Deal
utilizaron una forma matemtica particular para las contribuciones combinatoriales, la del
la mejor frmula que permita conseguir los resultados ms correctos, se han desarrollado
varios mtodos predictivos. De entre ellos destacamos los dos ms utilizados en los estu-
dios del equilibrio lquido-vapor, el mtodo ASOG (Analytical Solutions of Groups) pro-
puesto por Derr y Deal (1969) y el mtodo UNIFAC (UNIquac Functional-group Activity
Coefficients) de Fredenslund y col. (1975, 1977a); de este ltimo modelo han surgido dos
(1987).
El principal objetivo de estos dos mtodos es utilizar los datos de ELV existentes
en la bibliografa para predecir el equilibrio de fases de otros sistemas para los cuales no
se dispone de datos experimentales. Para ello, con datos que presentan una calidad ade-
cuada, se obtienen los parmetros que caracterizan las interacciones entre parejas de gru-
pos funcionales, utilizndose luego para predecir los coeficientes de actividad de otros
sistemas de los que no se posee informacin experimental, pero que contienen los mis-
Como resumen de las variantes principales que presentan estos dos mtodos, pue-
69
Tabla (6-1). Ecuaciones propuestas por los distintos modelos de contribucin de grupos para la estimacin de lniC y lnk.
-70-
UNIFAC utiliza el potencial de Staverman. Por lo que respecta a los coeficientes de acti-
(1987) tambin modifica los parmetros de volumen, Rj, y de rea, Qj, determinndolos a
res de los parmetros de interaccin entre los grupos, para las mezclas consideradas en
este trabajo, segn los distintos modelos y versiones, aparecen recogidos en los apndi-
ces (A-2) a (A-5). En la Tabla (6-1) se presenta un esquema-resumen con las expresiones
utilizadas por los distintos modelos citados, tanto para calcular la parte combinatorial, iC,
0 0 0
71
II. PARTE EXPERIMENTAL
comerciales con una elevada calidad, evitando as tener que recurrir a tcnicas de purifi-
desgasificacin de todas las sustancias que van a utilizarse con un equipo de ultrasoni-
dos durante varias horas. Para la eliminacin de posibles trazas de humedad se le aadi
a cada sustancia cierta cantidad de tamiz molecular, Fluka, tipo 0,3 nm, y se mantuvo en
la oscuridad durante varios das. Despus de estas operaciones, se midieron sus propie-
o
ebullicin normal, T b,i , o sea, a la presin atmosfrica, cuyos valores se recogen en la
Tabla (7-1) junto a otros de la bibliografa para su comparacin. En general, los valores
lizaron las experiencias del ELV, el cual ser descrito en la seccin 8.1. y las temperatu-
ras se midieron con una precisin de 0,01 K. Las densidades se obtuvieron con un den-
smetro digital de la firma Anton Paar, modelo DMA-55, ver seccin 8.3.1. que tiene
con un bao termosttico Heto Birkerod PT-623. El refractmetro empleado para medir
los n(D,298,15 K), a la longitud de onda de 5893 nm correspondiente a la lnea D del so-
dio, procede de la firma Baush & Lomb y proporciona una lectura de 0,0001 unidades.
75
Tabla (7-1). Caractersticas fsicas de los productos utilizados y comparacin con los datos de la bibliografa.
77
8. SISTEMA EXPERIMENTAL PARA LAS MEDIDAS DE ELV
ELV, conlleva la obtencin de las concentraciones de las fases lquida y vapor, cuando
mtodo empleado para alcanzar los estados de equilibrio, los tipos de recinto ms usa-
dos se clasifican en, los de puntos de burbuja, los de punto de roco, los estticos, los de
vapor a partir de una mezcla lquida en ebullicin que luego se condensa y se recoge
acuerdo con la forma de recirculacin de las fases, los recintos anteriores pueden clasifi-
carse en, circulacin de fase vapor y circulacin de las fases vapor y lquido. Resu-
capacidad de carga de producto y por ello puede considerarse de reducido volumen. Este
equipo fue inicialmente presentado por de Alfonso y col. (1983). Con posterioridad, en
78
nuestro laboratorio se han venido desarrollando algunas modificaciones del mismo, so-
pared (A), por calentamiento elctrico con una resistencia elctrica. El tubo (B), que ac-
T1 que seala la temperatura de ebullicin de la mezcla. Dicho tubo (C) hace de pantalla
para evitar, junto al ensanchamiento del equipo (D) y a la salida lateral (H), el arrastre de
Por un lado, el vapor despus de atravesar (H) y libre de gotas de lquido por lo
que se ha comentado, se condensa en el refrigerante (J) y gotea sobre el colector (L). Por
otro, el lquido cae sobre el embudo (E) que lo canaliza hacia (F). Desde los colectores
lquido rebosan hacia (G), donde se mezclan, regresando el lquido de nuevo haca el
fluida debido a las fuertes fluctuaciones que se producen en la ebullicin. La parte supe-
79
P
T2 H
J
C
D
ML
E Ag
ua
B
K
MV
F
L
Ll1
Ll2
Figura (8-1). Ebullmetro de equilibrio. (A) recipiente invertido de ebullicin; (B) bomba Cotrell; (C)
apantallamiento; (F) colector de fase lquida; (J) refrigerante; (K) cuentagotas; (L) colector de fase vapor;
(P) conexin al sistema de presurizacin.
80
mantiene calorifugado con lana de vidrio y un sostn de yeso. Adems, en la zona (D) se
sita una resistencia elctrica, alimentada por otra fuente de tensin estabilizada, cuyo
efecto se observa por un segundo sensor situado en T2. Una vez establecido el equilibrio,
desde (L) y (F), a travs de las tubuladuras laterales (MV) y (ML) se extraen mediante
agujas hipodrmicas adecuadas, las muestras necesarias para determinar las composi-
ciones de las fases vapor y lquida, respectivamente. Las llaves (Ll1) y (Ll2) permiten el
pe- riencia que produce la observacin continuada del funcionamiento y los distintos
ensayos realizados, cabe puntualizar los siguientes aspectos: a) para un sistema o mezcla
vaporizacin, de tal forma que, a mayor caudal de vaporizacin menor tiempo de opera-
cin, b) cuando las diferencias de las temperaturas de ebullicin de los compuestos pu-
ros son importantes, es necesario un periodo de tiempo ms largo para alcanzar el equi-
T2T1, se produce una condensacin parcial de la fase vapor; sin embargo, cuando la
81
8.2. Equipos auxiliares
cuenta con unos dispositivos auxiliares capaces de controlar y medir las variables ms
rias del estado de equilibrio, como son la temperatura y la presin. En la Figura (8-2) se
presenta el esquema de conexin de todos los equipos auxiliares entre s que constituyen
el sistema experimental montado, as como otra serie de elementos necesarios que com-
ponen el resto de la instalacin, como son manmetros en lnea, llaves de dos y tres v-
as, etc.
vapor isobrico. Tanto para la consecucin de buenos datos de presiones de vapor de las
sustancias puras, o sea, parejas de (T, pio ), como para la realizacin de experiencias a
que alimente al equipo principal. En este caso, gran parte de las experiencias se realiza-
inyectando N2 seco de alta pureza, mientras que, para los estados de bajas presiones se
utiliza una bomba de vaco, (B), segn el esquema detallado en la Figura (8-2).
zaron dos equipos, el primero de ellos (CP) consiste en un controlador de la casa Fisher
82
Figura (8-2). Sistema experimental para la medida de datos de ELV: (A) botella de N2; (B) bomba de vacio; (CP) controlador de presin; (D0, D1, D2)
recipientes amortiguadores; (E, F) llaves; (G) desecante; (H) vlvula de aguja; (I) medidor de presin MKS; (J) sensor de presin; (K) termmetro digital; (L)
fuentes de alimentacin.
83
modelo VK1, cuyo rango de trabajo est comprendido entre 0,3 y 160 kPa, con una pre-
cisin del 0,7%. Conectado en serie con el anterior y con el fin de mantener la presin
dentro de unos intervalos ms pequeos, se situ una vlvula mecnica de aguja (H) de
la casa Edwards modelo VPC1, con una precisin del 0,01%, de la presin de trabajo.
Esto permiti conseguir una estabilidad en la presin, estimada en 0,02 kPa. La medi-
270 B, que presenta una precisin del 0,001% de la lectura mediante un sensor, modelo
390, de la misma firma francesa, con una precisin media de 0,05 % de la lectura. El
depsito regulador D0 amortigua las emboladas de la bomba de vacio, mientras que, los
situados cerca del equipo principal, D1 y D2, sirven de estabilidad a las alteraciones de
do. Por ello, su medida ha de hacerse con la ms absoluta exactitud, poniendo especial
termmetro digital de la firma Comarks Electronics Ltd., modelo 6800, con un rango de
medida de 73-500 K y una lectura de 0,01 K. Los sensores utilizados para la medida,
tanto de T1 como de T2 son resistencias de platino calibradas con las normas ITS-90,
84
con una precisin estimada del 0,01% de la lectura. Peridicamente se contrastaron sus
Las muestras extradas tal y como se indic en la seccin 8.1, con jeringas hipo-
equilibrio se determinan tambin con otras tcnicas analticas, como son la cromatogra-
comparndola en exactitud con la tcnica densimtrica actual, esta ltima fue utilizada
85
das para cada una de las muestras y de esta forma asegurar las correlaciones densidad-
las fases en equilibrio a travs de la densidad es, despus de observar que se ha alcanza-
do el equilibrio, elaborar unas curvas que servirn de patrn para cada sistema en parti-
cular, de la forma =(x). De esta manera, con la densidad de la mezcla se permite cal-
cnica de la firma Anton-Paar (Graz, Austria), modelo DMA-55, que incorpora en una
2 m + v
1/2
T= = 2 (8-2)
V K
86
Haciendo las operaciones pertinentes puede encontrarse una relacin lineal entre
T2=
4 2 m 4 2v
K
+
K
= B + A o sea , =
A
(
1 2
T B ) (8-3)
forma casi directa la densidad de un fludo, conocidas las constantes del aparato.
empleando dos sustancias de densidades bien conocidas a esas temperaturas. Los fluidos
utilizados en este trabajo fueron agua bidestilada y desgasificada con una conductividad
igualmente con ultrasonidos y tratado con tamiz molecular (Fluka, tipo 0,3 nm). La uti-
lizacin de estos productos de calibrado fue justificada por Ortega y col. (1985a). A las
temperaturas utilizadas en el presente trabajo, los datos de densidad para el agua y del
Tabla (8-1). Densidades del agua y nonano (Fluka, ref. 74252), a distintas
temperaturas, empleadas en la calibracin del equipo de densimetra.
Sustancia /(kgm-3)
298,15 K 313,15 K 318,15 K
Agua 997,043 992,214 990,216
87
partir de la ecuacin (8-3), se calculan las constantes A y B, las cuales, una vez introdu-
cidas en el contador digitalizado del equipo, preparan a este para mostrar directamente
za, bien con un termmetro digital Anton Paar DT-100-30 para valores de hasta 303,15
permite medir las temperaturas superiores de 313,15 y 318,15 K, ambos con una preci-
sin de 0,01 K. Los termmetros son continuamente calibrados segn la norma ITS-90.
DT-100-30
25.00
sensor de
temperatura
0.99704
DMA-55
Ultratermostato
Figura (8-3). Esquema de conexin de los distintos elementos que constituyen el equipo de densimetra.
Ultratermostto, Heto Birkerod, termmetro DT-100-30 (Anton Paar) y densimetro digital de tubo vi-
brante DMA-55 (Anton Paar).
88
9. ESTIMACION DE ERRORES EXPERIMENTALES
a una serie de errores que normalmente pueden dividirse en dos grandes grupos, Giam-
bernardino (1971), los errores sistemticos, como aquellos que se repiten constante-
un mismo sentido (pueden ser debidos a la calibracin de los instrumentos, a las condi-
rios, los cuales no dependen de una causa determinada y concreta ni afectan de igual
forma al resultado. Algunos de ellos son los de precisin, debidos a defectos de cons-
fsica medida sea directa o indirecta y, en ese ltimo caso, sea funcin de variables esta-
dsticas o no estadsticas. De forma concreta, en este trabajo nos encontramos con mag-
con una serie de variables, zi , de las cuales unas son estadsticas y otras no. Por lo tanto,
89
para considerar el tratamiento matemtico de los errores experimentales se tomar que
M M M
dM = dz1 + dz2 + ! + dzn (9-2)
z1 z 2 z n
de la funcin y de las variables por incrementos finitos que se identifican con los errores
correspondientes:
con lo cual se obtiene la expresin del error de dicha magnitud indirecta M en funcin
M M M
E(M ) = E ( z1 ) + E ( z2 ) + ! + E ( zn ) (9-4)
z1 z 2 z n
obtenidas en los trabajos de equilibrio lquido-vapor, como son los volmenes de exce-
90
9.1. Errores obtenidos en el clculo de los volmenes de exceso
sin (8-1), por lo tanto, como cantidad indirecta su error puede estimarse mediante:
v E v E v E v E
E( v E ) = E( 1 ) + E( 2 ) + E( ) + E( M 1 ) +
1 2 M 1
(9-5)
v E v E v E
+ E( M 2 ) + E( x1 ) + E( x2 )
M 2 x1 x2
5) en el error global. As, el error en la densidad, tanto de la mezcla como de los com-
puestos puros, ya citado en apartados anteriores fue de 0,02 kgm-3, las masas molecula-
res se estimaron con un error de 0,001 uma y el error medio de las fracciones molares
les, se obtiene para cada volumen de exceso un error experimental que siempre fue
vapor se obtienen a partir de la densidad siendo esta funcin de dos variables, la frac-
91
Por lo tanto, el error cometido en su determinacin podr calcularse por:
E( ) = E( x ) + E( T ) (9-5)
x T
segundo trmino se puede despreciar con lo cual, el error estimado en el clculo de las
1
E( x ) = E( ) (9-6)
x
polinmica sencilla, del tipo = Ai-1 x1i-1. Teniendo en cuenta que la incertidumbre en
cin de la fase lquida de 0,001, y de 0,002 unidades para la fase vapor (incluso te-
_________
Antes de efectuar las medidas de los datos de ELV a la presin de 141,3 kPa
do los resultados experimentales obtenidos para una mezcla que pudiera considerarse
92
patrn, por estar suficientemente representada en la bibliografa y con datos de calidad.
Como el trabajo correspondiente a esta Memoria se refiere a mezclas binarias con meta-
por ello la misma para verificar el funcionamiento del sistema experimental y la meto-
Para el caso que nos ocupa del ELV isobrico, aparecen bastantes referencias
cin de Chemistry Data Series editada por la firma alemana DECHEMA y en la de da-
tos de ELV presentada por Hirata y col. (1975), cuyas referencias concretas se irn men-
cionando en cada ocasin. Son muchos los trabajos publicados a la presin atmosfrica
normal, 101,32 kPa, sin embargo, no se han encontrado referencias a la presin de este
cacin de los datos por el clculo de los vE=(x) para, posteriormente, conseguir las
nes de exceso, los aqu obtenidos sern comparados adecuadamente con los de otros
autores. La validacin de los mismos servir para reafirmar las curvas densidad-
93
10.1. Datos de composicin-densidad-volumen de exceso para la mezcla metanol+
agua. Correlaciones
de estos datos y de las densidades de los compuestos puros, la ecuacin (8-1) permite el
las parejas de (x1,) se correlacionan mediante una ecuacin polinmica sencilla, del
tipo: = Ai-1 x1i-1, que constituye la base del clculo de las composiciones de ambas
fases en el ELV, como ya se indic. Para conseguir mejores ajustes en las fracciones
mnimos cuadrados, se muestran en la Tabla (10-1). Para todas ellas la bondad del
ajuste, r2, fue cercana a la unidad. La verificacin de los datos obtenidos se realiza a
pero donde ahora, kv=v2/v1, es constante, como un cociente entre los volmenes molares
de los compuestos puros a la temperatura de trabajo. Los valores de los coeficientes Vi,
94
Tabla (10-1). Densidad, , y volumen de exceso, vE, para el sistema binario x1metanol+x2agua a la tem-
peratura de 298,15 K.
x1 109 vE x1 109 vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0323 792,22 -118 0,6187 907,10 -951
0,0670 798,22 -241 0,6870 922,56 -873
0,1185 807,20 -403 0,7155 928,91 -824
0,1600 814,58 -519 0,7495 936,40 -752
0,1811 818,53 -578 0,7737 941,61 -692
0,2015 822,17 -624 0,8060 948,86 -610
0,2567 832,82 -757 0,8291 954,25 -550
0,2718 835,47 -777 0,8559 959,91 -460
0,3171 844,24 -855 0,8793 965,44 -390
0,3635 853,32 -914 0,9032 970,78 -305
0,4182 864,57 -973 0,9271 976,48 -224
0,4737 875,77 -992 0,9439 980,63 -166
0,5104 883,76 -1005 0,9663 986,74 -96
0,5512 892,67 -1003 0,9810 991,00 -51
0,5991 902,77 -966 - - -
[0,0-0,2] =997,0-327,35x+748,20x2-1751,79x3
[0,2-0,8] =995,61-245,81x+35,64x2
[0,8-1,0] =996,23-251,72x+42,29x2
mos cuadrados. En la Tabla (10-1) se recogen los valores conseguidos para la correla-
cin, mientras que, en la Figura (10-1a) se presenta la comparacin de los valores expe-
mentales y bibliogrficos, con un error medio inferior en todos los casos al 7%. como se
95
-1100 50
(a) (b)
30
-20
-30
0 0.2 0.4 x 0.6 0.8 1
1
-1050
Figura (10-1). Comparacin de los vE, a 298,15 K, experimentales y
los de la literatura para el sistema x1metanol+x2agua: (a), comparacin
-300 de valores absolutos; grfica superpuesta, comparacin en el intervalo
109 vE/(m3 mol-1 )
96
muestra, validndose de esta forma, los valores de densidad-composicin para la misma
cercano al mximo y el paso de la curva de ajuste para nuestros datos de la Tabla (10-
sentatacin, lo que supone errores inferiores al 2% de de los valores mximos de los vE.
Tabla (10-2). Valores experimentales del ELV del sistema x1metanol+x2agua a diferentes presiones.
97
358,35 0,127 0,480 1,807 1,019 0,091 271,5
356,80 0,154 0,530 1,729 1,008 0,091 270,3
355,65 0,172 0,552 1,681 1,025 0,109 323,3
355,35 0,183 0,565 1,628 1,021 0,106 313,1
353,90 0,210 0,600 1,583 1,029 0,119 349,0
352,74 0,245 0,631 1,485 1,037 0,125 365,2
351,84 0,274 0,653 1,417 1,052 0,132 387,1
351,62 0,278 0,658 1,418 1,052 0,134 390,5
350,75 0,311 0,680 1,348 1,068 0,138 403,5
349,30 0,353 0,716 1,315 1,070 0,140 407,5
348,82 0,394 0,727 1,217 1,119 0,146 422,6
347,50 0,452 0,757 1,158 1,160 0,148 427,5
346,92 0,486 0,774 1,125 1,176 0,141 405,5
346,35 0,513 0,785 1,103 1,208 0,142 410,0
345,35 0,560 0,808 1,078 1,245 0,138 397,0
344,41 0,603 0,829 1,060 1,283 0,134 384,2
343,25 0,671 0,859 1,030 1,334 0,115 327,5
342,20 0,732 0,886 1,013 1,380 0,096 273,6
341,57 0,777 0,905 0,996 1,428 0,077 217,3
340,74 0,826 0,925 0,989 1,481 0,059 168,5
340,15 0,859 0,941 0,989 1,479 0,045 128,4
339,39 0,913 0,961 0,977 1,631 0,022 61,5
338,54 0,963 0,980 0,976 1,985 0,002 5,5
337,83 0,995 0,997 0,986 2,007 -0,010 -29,4
337,42 1,000 1,000 1,000 - 0,000 0,0
141,3 kPa 382,65 0,000 0,000 - 1,000 0,000 0,0
381,30 0,006 0,050 2,714 0,999 0,006 17,5
380,35 0,010 0,081 2,703 1,002 0,012 38,0
379,35 0,016 0,114 2,532 1,004 0,019 59,5
378,58 0,021 0,145 2,484 1,000 0,019 59,1
377,14 0,031 0,178 2,107 1,020 0,042 132,7
375,54 0,043 0,235 2,125 1,015 0,047 145,8
374,41 0,050 0,266 2,140 1,020 0,056 175,6
372,99 0,067 0,311 1,954 1,023 0,066 203,4
369,85 0,100 0,401 1,849 1,029 0,087 268,2
368,70 0,115 0,428 1,763 1,043 0,103 315,6
367,45 0,129 0,460 1,759 1,046 0,112 341,3
366,45 0,140 0,484 1,761 1,049 0,120 366,7
365,47 0,160 0,515 1,689 1,046 0,122 369,9
363,39 0,197 0,565 1,613 1,057 0,139 418,7
362,35 0,221 0,593 1,553 1,062 0,144 433,8
361,56 0,245 0,617 1,494 1,061 0,144 431,6
359,35 0,312 0,674 1,380 1,076 0,151 451,6
358,45 0,352 0,697 1,304 1,097 0,153 456,0
98
356,77 0,431 0,738 1,189 1,154 0,156 462,9
354,55 0,522 0,788 1,129 1,209 0,154 454,3
352,75 0,617 0,835 1,076 1,263 0,134 394,4
351,25 0,698 0,870 1,043 1,338 0,117 342,1
349,95 0,763 0,902 1,035 1,348 0,097 282,4
348,99 0,829 0,927 1,011 1,453 0,073 212,4
348,35 0,867 0,944 1,008 1,468 0,058 166,7
347,35 0,931 0,971 1,000 1,539 0,029 85,1
345,92 1,000 1,000 1,000 - 0,000 0,0
dicha tabla recoge los valores calculados de las composiciones de lquido, x1, y vapor,
y1, determinados con las funciones polinmicas mostradas en la Tabla (10-1), =(x).
lquida, calculados segn se detalla en la seccin 3.4 y los valores de la funcin de ener-
mtodos para los datos de ELV isobricos del sistema metanol+agua a las dos presio-
nes. Asimismo, se aplicaron otros test que evalan de forma global la consistencia de los
datos, como el test de Herington (1951) y la modificacin al mismo propuesta por Wis-
niak (1994). Ambos dieron como resultado consistencia positiva para el sistema reali-
de 141,3 kPa. Los valores de entalpa de mezcla necesarios en la aplicacin del test de
99
375 0.4 500 3
(a) (b)
400
365 0.3 2.5
g /(Jmol )
-1
y1-x1
T/K
E
1.5
345 0.1 100
0 1
335 0.0 -25 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (10-2). Valores experimentales y curvas de ajuste obtenidos para el ELV isobrico del sistema
x1metanol+x2agua a la presin de 101,32 kPa. (a), T vs x1 o y1; (y1-x1) vs x1 (b), gE o i vs x1.
400
375
0.3 2.5
g /(Jmol )
-1
y1-x1
T/K
i
E
1.5
0.1 100
350
0 1
340 0.0 -25 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (10-3). Valores experimentales y curvas de ajuste obtenidos para el ELV isobrico del sistema
x1metanol+x2agua a la presin de 141,3 kPa. (a), T vs x1 o y1; (y1-x1) vs x1 (b), gE o i vs x1.
100
En las Figuras (10-2) y (10-3) se representan los valores experimentales y las
ken (1977), Gmehling y col. (1981), Gmehling y col. (1988)] se recogen valores de la
mezcla metanol+agua, en condiciones isobricas, entre 8 y 506,6 kPa. Por otro lado
Hirata y col. (1975) presentan asimismo una coleccin de datos de ELV para la misma
mezcla, con datos en condiciones ms extremas, entre 810,6 y 1134,8 kPa. Con estas
bases de datos, se considera que existe un nmero suficiente de estudios para validar los
valores obtenidos con nuestro sistema experimental, sobre todo, aquellos correspon-
dientes a la presin de 141,3 kPa, presin de trabajo para la que no existen datos en la
bibliografa consultada.
ciones de presin que aparecen en la literatura, junto a, en rojo, las curvas correspon-
dientes a las dos presiones determinadas para este trabajo. De entre el conjunto de refe-
rencias que existen, se seleccionaron aquellas cuyos datos parecieron tener una mejor
calidad. Al no existir ninguna referencia bibliogrfica con valores a 141,3 kPa, se ob-
(y1-x1) frente a la presin de trabajo, p, ver figura superpuesta en (10-4); una adecuada
correlacin fue (y1-x1)=a+bpc. Esto permiti efectuar una interpolacin de datos, esti-
101
0.5 0.02
(a)
1 : 8.0 kPa
1
2 : 13.3 kPa
2
3 3 : 26.7 kPa 2%
4 4 : 46.7 kPa
(y 1-x1)
5 : 101.3 kPa
0.4 5 6 : 141.3 kPa
0
6 7 : 304.0 kPa
8 : 506.6 kPa
0.3
9
-0.02
0 0.2 0.4 0.6 0.8 1
y1-x1
10
x1
0.02
0.2 (b)
2%
0.50
0.246
(y 1-x1)
(y1-x1)=0.524-0.046p (para x 1=0.2)
0
0.1
y1-x1
-2%
0.25
0 250 p/kPa 1000
0 -0.02
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
x1
Figura (10-4). Curvas de (y1-x1) vs x1 para el ELV isobrico del sistema Figura (10-5). Desviaciones, (y1-x1)=(y1-x1)lit-(y1-x1)exp, para el sistema
x1metanol+x2agua a distintas presiones de: 1-2-3-4, Olevsky y Golubev x1metanol+x2agua entre la curva experimental y los presentados por: (a),
(1956)(Gmehling y Onken, 1977); 5, Ocn y Rebolleda (1958) (Gmehling Ocon y Rebolleda (1958) (ver Gmehling y Onken, 1977) a 101,32 kPa;
y Onken 1977); 6, valores interpolados; 7-8-9-10, Chuem y col. (1967)(Hi- (b), los valores estimados por interpolacin a 141,3 kPa.
rata y col, 1975)( ), curvas de este trabajo a las presiones indicadas.
102
la Tabla (10-2). En la Figura (10-5a) se han representado las diferencias, (y1-x1)=(y1-
como antes se mencion, a la presin de 141,3 kPa, aparece en la Figura (10-5b). Para
ambos casos se apreci un error medio inferior al 2%, lo que permite validar la experi-
0 0 0
103
III. PRESENTACION, TRATAMIENTO Y
ANALISIS DE RESULTADOS
conlleva una cuantiosa e importante experimentacin. Por ello, esta seccin est creada
ra, T, mientras que otras son determinadas indirectamente, como las composiciones de
las fases lquida, xi, y vapor, yi, de acuerdo con la metodologa indicada en la seccin 8
de esta Memoria. A partir de ellas, otras magnitudes termodinmicas son calculadas para
tudio.
Como nota aclaratoria para la presentacin de los datos valga decir que es cos-
tumbre en los trabajos de ELV presentar las composiciones de las sustancias en equili-
brio a travs de la fraccin molar del componente ms voltil, citado siempre en primer
lugar; sin embargo, con la idea de facilitar la comparacin de las magnitudes para los
estos dos aspectos resultan casi coincidentes, ya que el metanol es ms voltil que los
destacado que el primero, o sea la de referirlos a un compuesto comn, todas las compo-
siciones de lquido y vapor estarn referidas al metanol. Concretamente, los datos expe-
-107-
rimentales presentados en esta seccin corresponden a los determinados para las dos
noato de metilo, nico ster no estudiado en esta investigacin por las dificultades que
entraa, en primer lugar su adquisicin comercial y con gran pureza y en segundo, por su
delicada manipulacin.
Siguiendo el mismo orden del trabajo experimental, se presentan los valores con-
nol para las distintas mezclas binarias indicadas y finalmente, los datos correspondientes
xano, cuyos valores aqu determinados, son asimismo incluidos en un apartado especfi-
co de esta seccin.
Las presiones de vapor de los compuestos puros, pio, y sus correspondientes co-
nitudes que definen el ELV, como son los coeficientes de actividad o la funcin de ener-
-108-
ga de Gibbs, repercutiendo por tanto, en el tratamiento y reduccin de datos termodin-
el equilibrio de cada sistema. Por ello, resulta necesario confirmar los datos y/o las co-
nuevos valores experimentales de presin de vapor para el metanol, el agua, los dos n-
butanoato de etilo y para el etanoato y butanoato de propilo. Los datos obtenidos direc-
tamente se recogen en las Tablas (11-1) a (11-10), para el resto de los steres no men-
cionados se utilizaron las correlaciones conseguidas con datos experimentales de este la-
la Tabla (11-11).
Antoine:
log(pio/kPa)=A-B(T/K-C) (11-1)
siendo calculadas las constantes A, B y C con una versin del procedimiento de regresin
no-lineal Simplex, modificado por Nelder y Mead (1965), que minimiza la desviacin
estndar, s(pio), de los datos correlacionados. Los valores obtenidos para estas constan-
tes, junto con la bondad de ajuste calculada para cada caso, as como el intervalo de tem-
-109-
Tabla (11-1). Presiones de vapor experimentales, pio, en funcin de la temperatura para el metanol.
-1
pi /kPa
o
-3 - - - - - TRC (1993) 2%
Ambrose y Sprake (1970)
-4
315 325 345 355
T/K
Figura (11-1). Representacin de pio = pio,lit pio,exp , entre las curvas obtenidas con la eq. (11-1) usando los
parmetros A, B, C de la literatura y los conseguidos con valores de la Tabla (11-1) para el metanol.
-110-
Tabla (11-2). Presiones de vapor experimentales, pio, en funcin de la temperatura para el agua.
-0.5%
-1
pi /kPa
TRC (1993b)
-2
345 360 T/K 380 400
Figura (11-2). Representacin de pio = pio,lit pio,exp , entre las curvas obtenidas con la eq. (11-1) usando los
parmetros A, B, C de la literatura indicada y los conseguidos con valores de la Tabla (11-2) para el agua.
-111-
Tabla (11-3). Presiones de vapor experimentales, pio, en funcin de la temperatura para el n-pentano.
-1
pi /kPa
o
1.5%
-2
- - - - - TRC (1993)
Tenn y Missen (1963)
-3
300 310 320 330
T/K
Figura (11-3). Representacin de, entre las curvas obtenidas con la eq. (11-1) usando los parmetros A, B, C,
de la literatura y los conseguidos con valores de la Tabla (11-3) para el n-pentano.
-112-
Tabla (11-4). Presiones de vapor experimentales, pio, en funcin de la temperatura para el n-hexano.
0
pi /kPa
o
1%
-2 - - - - - TRC (1993)
Raal y col. (1972)
-3
330 340 355 365
T/K
Figura (11-4). Representacin de pio = pio,lit pio,exp , entre las curvas obtenidas con la eq. (11-1) usando los
parmetros A, B, C de la literatura y los conseguidos con valores de la Tabla (11-4) para el n-hexano.
-113-
Tabla (11-5). Presiones de vapor, pio, en funcin de la temperatura para el metanoato de etilo.
2
1%
0
p i /kPa
Figura (11-5). Representacin de pio = pio,lit pio,exp , entre curvas obtenidas con la eq. (11-1) usando los pa-
rmetros A,B,C, de la literatura y los conseguidos con valores de la Tabla (11-5) para el metanoato de etilo.
-114-
Tabla (11-6). Presiones de vapor, pio, en funcin de la temperatura para el etanoato de etilo.
0.5
0
p i /kPa
o
Figura (11-6). Representacin de pio = pio,lit pio,exp , entre curvas obtenidas con la eq. (11-1) usando los par-
metros A, B, C, de la literatura y los conseguidos con valores de la Tabla (11-6) para el etanoato de etilo.
-115-
Tabla (11-7). Presiones de vapor, pio, en funcin de la temperatura para el propanoato de etilo.
0.5
0
p i /kPa
-0.5%
o
Figura (11-7) Representacin de pio = pio,lit pio,exp , entre las curvas obtenidas con la eq. (11-1) usando los pa-
rmetros A, B, C, de la literatura y los conseguidos con valores de la Tabla (11-7) para el propanoato de etilo.
-116-
Tabla (11-8). Presiones de vapor, pio, en funcin de la temperatura para el butanoato de etilo.
2.5
0
p i /kPa
o
Figura (11-8) Representacin de pio = pio,lit pio,exp , entre las curvas obtenidas con la eq. (11-1) usando los pa-
rmetros A,B,C de la literatura y los conseguidos con valores de la Tabla (11-8) para el butanoato de etilo.
-117-
Tabla (11-9). Presiones de vapor, pio, en funcin de la temperatura para el etanoato de propilo.
0
p i /kPa
Figura (11-9). Representacin de pio = pio,lit pio,exp , entre las curvas obtenidas con la eq. (11-1) usando los pa-
rmetros A,B,C de la literatura y los conseguidos con valores de la Tabla (11-9) para el etanoato de propilo.
-118-
Tabla (11-10). Presiones de vapor, pio, en funcin de la temperatura para el butanoato de propilo.
5
2%
0
p i /kPa
o
-2%
TRC (1993b)
-7
Frkov y Wichterle (1993)
Ortega y Galvn (1995)
-10
375 395 T/K 425 440
Figura (11-10). Representacin de pio = pio,lit pio,exp , entre las curvas obtenidas con la eq. (11-1) usando los pa-
rmetros A,B,C, de la literatura y los conseguidos con valores de la Tabla (11-10) para el butanoato de propilo.
-119-
Tabla (11-11). Coeficientes A, B y C de la ecuacin de Antoine, rango de temperaturas, T, de aplicacin y desviacin estndar, s(pio), para los com-
puestos utilizados en este trabajo. Valores obtenidos con resultados experimentales de este trabajo y comparacin con los de la literatura.
120
6,18798 1224,67 57,44 260-380 Ambrose y col. (1981)
6,13360 1195,13 60,68 250-380 TRC (1993b)
Propanoato de etilo 6,30110 1368,91 53,45 335-400 0,08 Exp.
6,14396 1274,70 64,16 260-400 Reid y col. (1966)
6,12169 1260,62 65,75 300-380 Boublick y col. (1973)
6,13486 1268,94 64,85 300-380 Riddick y col. (1986)
6,14285 1274,70 64,15 260-400 TRC (1993b)
Butanoato de etilo 6,36030 1493,89 51,13 355-420 0,06 Exp.
6,07304 1358,30 60,15 275-445 Reid y col. (1966)
6,34123 1515,12 44,72 245-400 Gmehling y col. (1982b)
5,90606 1213,54 83,10 325-400 Frkov y Wichterle (1993)
5,27455 921,56 112,77 275-440 TRC (1993b)
Metanoato de propilo 6,07002 1184,43 61,89 325-375 Galvn y col. (1994)
Etanoato de propilo 6,51158 1524,56 36,19 340-405 0,06 Experimentales
6,14104 1282,28 64,55 200-390 Boublick y col (1973)
6,14361 1284,08 64,36 280-405 Ambrose y col. (1981)
6,39229 1411,45 52,84 260-390 Gmehling y col (1982b)
6,07166 1240,55 69,10 325-385 Frkov y Wichterle (1993)
6,05547 1233,46 70,07 260-410 TRC(1993b)
Propanoato de propilo 6,06538 1301,03 75,11 360-420 Ortega y Galvn (1994)
Butanoato de propilo 6,60941 1723,42 42,02 375-440 0,08 Experimentales
6,32481 1543,16 59,45 345-430 Frkov y Wichterle (1993)
5,40455 1019,49 116,55 300-460 TRC (1993b)
5,97273 1331,40 80,07 375-450 Ortega y Galvn (1995)
Metanoato de butilo 6,70820 1627,85 33,17 350-400 Gonzlez y Ortega (1995)
Etanoato de butilo 6,44352 1584,05 42,20 350-450 Gonzlez y Ortega (1995)
Propanoato de butilo 6,57595 1745,90 36,04 350-450 Gonzlez y Ortega (1995)
Butanoato de butilo 6,13153 1500,71 74,43 400-450 Gonzlez y Ortega (1995)
121
Para el rango de temperaturas en las que se ha llevado a cabo las determinaciones
experimentales de cada una de las sustancias, en las Figuras (11-1) a (11-10) se han re-
presentado las diferencias entre las curvas correspondientes, calculadas con la ecuacin
de Antoine, cuando se emplean las constantes obtenidas con los nuevos valores experi-
mentales y las proporcionadas por la literatura, o sea, con los datos que aparecen en la
Tabla (11-11). El error medio en la mayora de los casos es inferior al 2%; sin embargo,
el resumen para cada una de las sustancias sobre la comparacin con lo existente en la
literatura es:
Para el metanol, las mayores diferencias de nuestros datos aparecen, en temperaturas inferiores
a los 330 K con la curva publicada en el TRC (1993b), pero no superiores al 2%.
Para el agua se han obtenido datos experimentales y correlaciones acordes con los de la biblio-
grafa, con errores medios inferiores al 0,5% en todos los casos.
Para el n-pentano, nuestros datos son regularmente inferiores en 1,4% con relacin a los de
Tenn y Missen (1963), determinados hace ms de treinta aos.
Las curvas del n-hexano estn en concordancia a otras publicadas, con errores inferiores al 1%,
en un rango de 35 K.
Para el metanoato de etilo, las diferencias con relacin a otros de la literatura son, asimismo,
inferiores al 1%. Este hecho se repite igualmente para el etanoato y propanoato de etilo.
Los valores determinados para el butanoato de etilo son quiz los que presentan mayores dife-
rencias con relacin a otros ya publicados, concretamente supera el 3,5% en relacin a los del TRC
(1993b); por el contrario, es inferior al 0,4% en relacin a los ms recientes de Frkov y Wichterle
(1993). Idntico anlisis cabra realizar con el butanoato de propilo, incluyendo una diferencia cercana
al 3% con datos obtenidos anteriormente, Ortega y Galvn (1995), en nuestro mismo laboratorio, debido
posiblemente a errores de medicin en la temperatura del equilibrio.
La diferencia ms significativa del etanoato de propilo es la de Gmehling y col. (1982b), ligera-
mente superior al 1%.
La ecuacin de Antoine, la (11-1), establece para una sustancia pura una relacin
go, es interesante para estas magnitudes utilizar una base comn adimensional basada en
-122-
0 -0.5
(a) (b)
-0.5
steres de metilo
steres de etilo
steres de propilo
-1.3 steres de butilo
n-hexano
o
o
log p r
log p r
metanoato etilo
metanol
n-pentano -2.0
metanol
-2
T r=0.7
T r=0.7
-2.5 -2.8
1.1 1.3 1.6 1.8 1.2 1.3 1.6 1.8
1/Tr 1/Tr
Figura (11-11). Representacin de log pr vs 1/Tr para el metanol y: (a), n-alcanos; (b), steres de alquilo.
Para simplificar, en la grfica (b) slo se muestran las lneas de algunos steres: etanoato y propanoato de
metilo; del metanoto al butanoato de etilo; metanoato, etanoato y butanoato de propilo; y butanoato de bu-
tilo. Los valores de los coeficientes de log pro=a+b/Tr para cada sustancia son:
Compuesto a b Compuesto a b
Metanol 3,829 -3,773 Butanoato de etilo 3,555 -3,477
n-Pentano 3,045 -3,004 Metanoato de propilo 3,334 -3,242
n-Hexano 3,223 -3,164 Etanoato de propilo 3,474 -3,406
Metanoato de metilo 3,344 -3,260 Propanoato de propilo 3,562 -3,457
Propanoato de metilo 3,376 -3,294 Butanoato de propilo 3,703 -3,607
Butanoato de metilo 3,506 -3,411 Metanoato de butilo 3,625 -3,505
Metanoato de etilo 3,302 -3,203 Etanoato de butilo 3,488 -3,435
Etanoato de etilo 3,388 -3,322 Propanoato de butilo 3,624 -3,565
Propanoato de etilo 3,519 -3,437 Butanoato de butilo 3,641 -3,585
con magnitudes reducidas. As, en la Figura (11-11a,b) se representa log pro vs 1/Tr para
cada una de las sustancias de este trabajo y sus respectivas ecuaciones lineales. Tanto pa-
ra los n-alcanos como para los steres se observa una cierta convergencia de las lneas
-123-
correspondientes a los estados de saturacin de dichos compuestos hacia un nico punto.
nito, permite determinar la recta de presiones de vapor con un solo valor para un miem-
bro de la serie. En el caso de las distintas familias de los steres presentados en este tra-
bajo, las coordenadas (1/Tr;log pro) de dicho punto del infinito se estimaron situados en:
(1,024; 0,009); steres de butilo, (0,946; 0,263). De ello, hasta podra concluirse que to-
dos los steres poseen prcticamente el mismo punto del infinito; calculando un punto
comn para este conjunto que se establecera en: log pr=0,169, 1/Tr=0,95. Para los n-
1/Tr=1,112. Sera necesario disponer de datos de otros hidrocarburos similares para de-
ducidas permite obtener de manera sencilla el factor acntrico de las sustancias, tal como
fue definido por Pitzer (1955), = log pro (a Tr = 0,7 ) 1 . Dichos valores han sido de-
terminados y presentados en la Tabla (A-1) del apndice; la comparacin con los refe-
-124-
binaria metanol+n-hexano a la temperatura de 298,15 K, se realizaron nuevas determina-
Upper Critical Solubility Point, y con los mismos productos que van a ser utilizados en
este trabajo. Dichas medidas se realizaron en una clula de vidrio de, aproximadamente,
La temperatura fue obtenida con el equipo descrito en la seccin 8.2.2, siendo la preci-
sin de esta magnitud mejor que 0,02 K. Las fracciones molares, x1, en metanol, se de-
terminaron por pesada con un error inferior a 210-4 unidades. Los valores experimen-
del tipo, T= Bi-1x1i-1. El valor del UCSP, estimado matemticamente, correspondi a las
-125-
310 4
(a) (b)
305
T
0
-2
T/K
(11-12), delimita las fases fluidas existentes. La Figura (11-12b) recoge la comparacin
medio entre nuestra curva experimental y los valores aportados por los distintos autores
del UCSP eligindose, para asegurar la homogeneidad de las fases en las determinacio-
-126-
11.3. Resultados de composicin-densidad-volumen de exceso. Correlaciones
determinarn por densimetra, por ello, previamente deben obtenerse datos para las cur-
vas patrn de las mezclas objeto de este trabajo, x1metanol+x2n-alcano o+x2ster de al-
quilo y sus respectivas correlaciones =(x1). Dichas curvas sern validadas observando
ca sencilla del tipo = Ai-1x1i-1. Con el fin de disminuir los errores en el reclculo de las
centraciones en tres zonas, [0,0-0,2], [0,2-0,8] y [0,8-1,0]. Los valores de los coeficien-
tes respectivos se determinaron por un mtodo de mnimos cuadrados. Por lo que res-
pecta a los volmenes de exceso se emple la ecuacin (10-1) para correlacionar las
parejas (x1,vE), definiendo el parmetro kv como el cociente entre los volmenes molares
de los compuestos puros a la temperatura de trabajo, v2/v1, segn proponen Ortega y Al-
calde (1992). El nmero de coeficientes, Vi, se fija optimizando el grado del polinomio
-127-
11.3.1. Mezclas metanol+n-alcano
presentados como x1--vE para cada una de las mezclas binarias x1metanol+x2n-alcano.
, y los correspondientes volmenes de exceso, vE, de dicha mezcla binaria fueron de-
(298,15 0,01) K. En las Tablas (11-13) a (11-15) se recogen tambin los coeficientes
11.3. En la Figura (11-13) aparecen las curvas experimentales de vE vs x1 para los siste-
mas x1metanol+x2n-alcano(C5 o C6) a las temperaturas indicadas, junto con las curvas
este trabajo. De dicha revisin se tom la curva recomendada para los datos de Tenn y
Missen (1963) para la mezcla metanol+n-pentano, la cual presenta una buena concor-
dancia con los experimentales nuestros a fracciones molares elevadas de metanol, sin
embargo, la discrepancia es algo mayor a bajas concentraciones, con un error medio glo-
bal inferior al 9%. Para el sistema metanol+n-hexano, Costas y col. (1987) aportan datos
experimentales, con un error medio en dicho intervalo inferior al 4%. Puesto que no se
-128-
Tabla (11-13). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2n-pentano a
la temperatura de 298,15 K y valores de los coeficientes de las distintas correlaciones.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0411 623,35 167 0,5159 663,48 377
0,0830 625,40 264 0,6014 675,53 355
0,1181 627,38 317 0,6991 692,67 295
0,2237 634,52 384 0,7629 706,35 243
0,2543 636,91 388 0,8149 719,37 195
0,3628 646,41 396 0,8683 734,85 145
0,4375 654,07 397 0,9476 763,34 55
0,4700 657,78 391
k=0,173* V0=5717 V1=-7616 V2=3175 s(vE)=5
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =621,78+36,15x1+92,22x12
[0,2-0,8] =614,83+99,13x1-86,92x12+149,62x13
[0,8-1,0] =316,71+1251,11x1-1579,99x12+798,96x13
* El parmetro k fue obtenido en este caso por optimizacin en la correlacin de vE vs x1.
Tabla (11-14). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2n-hexano a
la temperatura de 298,15 K y valores de los coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0486 655,79 175 0,8032 723,26 412
0,0633 656,29 202 0,8137 725,46 404
0,1041 657,70 282 0,8338 730,26 359
0,1173 658,23 298 0,8702 739,22 322
0,1643 660,20 355 0,9068 749,98 248
0,1653 660,25 355 0,9436 762,76 154
0,1979 661,76 385 0,9867 780,27 50
0,2125 662,37 415
0,2425 663,46 443
kv=3,232 V0=3598 V1=-19687 V2=68330 V3=-85855 V4=37388 s(vE)=13
-129-
Tabla (11-15). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2n-hexano a
la temperatura de 318,15 K y valores de los coeficientes de las distintas correlaciones.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0425 636,80 197 0,5591 667,73 651
0,1003 638,34 398 0,6068 673,00 631
0,1338 639,46 481 0,6457 677,86 605
0,1914 641,78 576 0,6928 684,54 557
0,2476 644,44 638 0,7281 690,24 513
0,3063 647,63 678 0,7809 699,91 443
0,3875 652,88 695 0,8727 721,37 295
0,4526 657,81 693 0,9455 744,71 138
0,5113 662,95 676
kv=3,246 V0=4888 V1=-20852 V2=69274 V3=-97779 V4=49245 s(vE)=9
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =635,97+16,52x1+71,55x12
[0,2-0,8] =630,29+65,99x1-67,31x12+123,02x13
[0,8-1,0] =252,85+1551,71x1-2027,12x12+990,19x13
0.8
Tenn y Missen (1963)
n-C 6(298,15 K)
n-C 6(318,15 K)
0.6
Costas y col. (1987)
Liu y col. (1991)
n-C 5(298,15 K)
x1
0.4
0.2
0
800
600
200
10 v /(m mol )
-1 3 9 E
-130-
temperatura de 318,15 K, en la representacin se comparan con los determinados por Liu
y col. (1991) a una temperatura ligeramente inferior, 313,15 K. Esta curva de la literatu-
ra tiene una forma casi simtrica, hecho bastante inusual en las mezclas donde intervie-
nen alcanoles, debido a las dificultades inherentes a los efectos asociativos del metanol,
que se destaca principalmente en las zonas de baja concentracin del compuesto asocia-
do, tal como se refleja en las otras mezclas presentadas en la Figura (11-13).
(del metanoato al butanoato) de alquilo (del metilo al butilo), a excepcin del metanoato
de metilo, se presentan en las Tablas (11-16) a (11-30). En ellas figuran adems los
vE, en funcin de la concentracin del metanol para cada sistema. Al mismo tiempo se
muestra la comparacin grfica entre las curvas experimentales de este trabajo y las que
correlaciones con las de otros autores pueden valorarse en la forma siguiente; con la de
Grolier y Viallard (1971) para la mezcla x1metanol+x2etanoato de metilo son del 14% y
el 20% con la curva publicada por Polak y Lu (1972b) para el sistema x1metanol+ x2pro-
-131-
Tabla (11-16). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2etanoato de
metilo a la temperatura de 298,15 K y coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0166 925,89 -7 0,5251 877,66 -85
0,0763 921,75 -34 0,5728 871,32 -86
0,1157 918,76 -44 0,6602 858,32 -74
0,1634 914,98 -55 0,7136 849,53 -65
0,2208 910,13 -66 0,7817 837,21 -54
0,2999 902,89 -79 0,8074 832,17 -47
0,3343 899,50 -81 0,8630 820,60 -36
0,4118 891,33 -90 0,9140 808,87 -21
0,4899 882,15 -89 0,9594 797,62 -14
kv=1,962 V0=-408 V1=154 s(vE)=2
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =926,99-65,10x1-51,10x12
[0,2-0,8] =928,36-79,03x1-2,48x12-58,35x13
[0,8-1,0] =999,76-346,57x1+334,59x12-201,15x13
Tabla (11-17). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2propanoa-to
de metilo a la temperatura de 298,15 K y coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0489 906,87 -12 0,5671 866,81 -77
0,1117 903,49 -29 0,6172 860,78 -72
0,1599 900,64 -37 0,6604 855,07 -63
0,1862 899,03 -43 0,6981 849,78 -62
0,2455 895,19 -56 0,7407 843,26 -56
0,3021 891,17 -64 0,8134 830,71 -45
0,3361 888,57 -66 0,8611 821,22 -32
0,3650 886,31 -72 0,9018 812,34 -24
0,3968 883,67 -76 0,9425 802,51 -17
0,4793 876,08 -77 0,9829 791,61 -5
kv=2,379 V0=-290 V1=-93 V2=153 s(vE)=2
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =909,34-48,37x1-37,62x12
[0,2-0,8] =908,95-44,25x1-51,49x12+25,79x13-50,39x14
[0,8-1,0] =1138,63-883,13x1+1008,94x12-477,79x13
-132-
Tabla (11-18). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2butanoato
de metilo a la temperatura de 298,15 K y coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0538 890,25 -9 0,5797 858,35 -62
0,0622 889,91 -9 0,6203 854,29 -58
0,1253 887,37 -26 0,6414 852,05 -57
0,1716 885,34 -35 0,6870 846,83 -55
0,2107 883,48 -39 0,7490 838,65 -42
0,2766 880,16 -54 0,8163 828,26 -32
0,3510 875,86 -59 0,8525 821,85 -26
0,4072 872,25 -64 0,8821 816,10 -22
0,4125 871,90 -66 0,8965 813,11 -20
0,4865 866,45 -66 0,9575 798,80 -8
0,5090 864,63 -65
kv=2,811 V0=-244 V1=-127 V2=262 s(vE)=2
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =892,24-35,34x1-29,33x12
[0,2-0,8] =890,73-19,52x1-85,96x12+89,99x13-85,93x14
[0,8-1,0] =1070,32-737,28x1+921,88x12 -467,94x13
0.6
4
x1
3
0.4
2
0.2
0
-100
-25
-75
0
10 v /(m mol )
-1 3 9 E
-133-
Tabla (11-19). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2metanoato
de etilo a la temperatura de 298,15 K y coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0568 912,54 -45 0,6020 861,86 -127
0,1413 906,77 -71 0,6449 855,95 -121
0,2169 901,16 -93 0,6807 850,81 -120
0,2356 899,67 -97 0,7066 846,73 -109
0,2519 898,39 -104 0,7453 840,58 -107
0,3222 892,48 -124 0,7852 833,78 -101
0,3843 886,60 -125 0,8203 827,33 -92
0,4268 882,39 -131 0,8631 818,88 -76
0,4701 877,79 -133 0,9175 807,06 -50
0,5192 872,21 -133 0,9677 795,23 -34
0,5763 865,24 -132
kv=1,986 V0=-628 V1=567 V2=-781 s(vE)=5
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =915,85-56,58x1-51,56x12
[0,2-0,8] =916,74-67,33x1-9,28x12-50,33x13
[0,8-1,0] =-1634,80+11596,48x1-19978,97x12+15137,12x13-4333,26x14
Tabla (11-20). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2etanoato de
etilo a la temperatura de 298,15 K y coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0650 891,46 -21 0,5810 856,05 -69
0,1155 889,03 -27 0,6077 853,27 -68
0,2045 884,40 -44 0,6567 847,73 -63
0,2521 881,69 -53 0,6903 843,66 -62
0,2665 880,80 -52 0,7218 839,55 -58
0,3009 878,69 -59 0,7856 830,36 -49
0,3376 876,28 -63 0,8426 820,97 -40
0,3794 873,35 -65 0,8867 812,83 -34
0,4274 869,76 -68 0,9262 804,60 -21
0,5327 860,78 -71 0,9654 795,66 -11
kv=2,419 V0=-269 V1=-42 s(vE)=2
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =894,27-41,38x1-33,90x12
[0,2-0,8] =896,45-60,45x1+25,74x12-71,06x13
[0,8-1,0] =990,24-432,90x1+520,73x12-291,19x13
-134-
Tabla (11-21). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2propanoa-to
de etilo a la temperatura de 298,15 K y coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0398 882,68 -8 0,6463 846,46 -39
0,1887 876,93 -31 0,7261 837,57 -34
0,2179 875,61 -32 0,7671 832,24 -28
0,2557 873,86 -37 0,7797 830,46 -24
0,3153 870,87 -44 0,8079 826,38 -24
0,3958 866,26 -45 0,8627 817,39 -18
0,4609 862,06 -50 0,8898 812,38 -14
0,5284 857,08 -49 0,9419 801,51 -9
kv=2,837 V0=-210 V1=88 s(vE)=2
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =884,01-31,75x1-31,01x12
[0,2-0,8] =887,34-59,95x1+51,26x12-87,10x13
[0,8-1,0] =1116,24-908,72x1+1106,89x12-527,53x13
Tabla (11-22). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2butanoato de
etilo a la temperatura de 298,15 K y coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0400 872,49 -6 0,5707 848,78 -31
0,0876 871,15 -12 0,6417 843,10 -23
0,2049 867,37 -22 0,6776 839,86 -23
0,1671 868,64 -16 0,7001 837,63 -19
0,2460 865,87 -26 0,7733 829,47 -17
0,3199 862,88 -31 0,8310 821,65 -14
0,3754 860,36 -34 0,8656 816,13 -9
0,4227 857,97 -34 0,9110 807,92 -7
0,4731 855,15 -33 0,9594 797,39 -2
0,5122 852,77 -33
kv=3,266 V0=-153 V1=110 s(vE)=2
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =873,60-26,20x1-21,04x12
[0,2-0,8] =872,38-13,24x1-68,15x12+76,93x13-76,90x14
[0,8-1,0] =1130,95-1018,41x1+1285,33x12-611,03x13
-135-
1
0.8
Grolier y Viallard (1971)
Akita y Yoshida (1963)
0.6
4
2
3
x1
2
0.4
0.2
0
-140
-180
-40
0
10 v /(m mol )
-1 3 9 E
panoato de metilo. En la Figura (11-14) se present tambin la curva obtenida por Polak
midad de resultados. Se han localizado varias referencias sobre los vE del sistema x1 me-
nishi y Shirai (1970) y Nikam y col. (1996) con las que nuestros valores difieren, respec-
20% y 89%. Menor desviacin, del 4%, se observa frente a los datos aportados por
nuestro laboratorio hace aos, Ortega y col. (1986), para el sistema x1metanol +x2meta-
-136-
Tabla (11-23). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2metanoato
de propilo a la temperatura de 298,15 K y coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0319 897,49 -14 0,4855 868,85 -121
0,0928 894,70 -40 0,5201 865,62 -119
0,1739 890,55 -67 0,5743 860,24 -119
0,2062 888,75 -75 0,6131 856,03 -117
0,2725 884,80 -93 0,6721 849,00 -113
0,3335 880,79 -106 0,6839 847,47 -110
0,3673 878,39 -111 0,8090 828,87 -86
0,4099 875,17 -116 0,9085 809,59 -47
0,4398 872,76 -118 0,9824 791,74 -10
kv=2,407 V0=-447 V1=-135 s(vE)=2
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =898,87-41,61x1-36,07x2
[0,2-0,8] =900,76-59,62x1+23,12x12-73,17x13
[0,8-1,0] =777,16+293,12x1-283,32x12
Tabla (11-24). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2etanoato de
propilo a la temperatura de 298,15 K y coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0287 881,63 -8 0,5223 856,59 -57
0,0729 880,08 -14 0,5583 853,73 -56
0,0817 879,81 -20 0,6040 849,79 -56
0,1896 875,54 -31 0,6473 845,63 -51
0,2427 873,18 -37 0,6811 842,13 -51
0,2910 870,85 -41 0,7127 838,54 -46
0,3077 870,04 -46 0,7488 834,15 -44
0,3468 867,96 -50 0,8074 826,00 -34
0,3804 866,05 -51 0,8809 813,73 -21
0,4132 864,07 -52 0,9197 806,11 -16
0,4419 862,26 -55 0,9574 797,71 -9
kv=2,842 V0=-213 V1=-27 s(vE)=2
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =882,55-31,40x1-29,44x12
[0,2-0,8] =881,70-21,98x1-66,58x12+71,07x13-75,15x14
[0,8-1,0] =1161,28-1060,98x1+1274,47x12-587,96x13
-137-
Tabla (11-25). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2propanoa-to
de propilo a la temperatura de 298,15 K y coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0463 874,24 -8 0,6037 847,51 -24
0,1164 872,17 -18 0,6407 844,38 -21
0,1569 870,86 -22 0,6860 840,15 -20
0,2197 868,65 -27 0,7205 836,56 -18
0,2775 866,37 -29 0,7786 829,68 -15
0,3469 863,32 -31 0,8069 825,88 -14
0,3778 861,83 -30 0,8854 813,28 -8
0,4479 858,09 -29 0,9342 803,48 -5
0,5329 852,78 -28 0,9608 797,32 -4
0,5547 851,25 -27
kv=3,258 V0=-179 V1=344 V2=-308 s(vE)=1
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =875,48-25,55x1-25,28x12
[0,2-0,8] =878,02-48,35x1+44,46x12-79,59x13
[0,8-1,0] =1161,49-1118,31x1+1400,03x12-656,35x13
Tabla (11-26). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2butanoato
de propilo a la temperatura de 298,15 K y coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0623 866,58 -6 0,6327 842,06 7
0,1344 864,75 -8 0,6647 839,51 8
0,1877 863,26 -11 0,7131 835,14 13
0,2307 861,96 -13 0,7658 829,63 17
0,2749 860,51 -14 0,7987 825,68 18
0,3196 858,91 -14 0,8434 819,55 19
0,4195 854,73 -7 0,8866 812,60 15
0,4737 852,09 -5 0,9214 806,02 14
0,5037 850,50 -4 0,9635 796,65 8
0,5830 845,63 3
kv=3,683 V0=-100 V1=401 s(vE)=1
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =867,98-21,13x1-21,51x12
[0,2-0,8] =865,26+6,98x1-124,22x12+165,76x13-123,84x14
[0,8-1,0] =1321,86-1690,25x1+2066,25x12-910,98x13
-138-
Figura (11-16). Volmenes de exceso a298,15 K de las mezclas binarias
1
0.8
x1
no publicados
0.4
0.2
0
-100
-140
30
10 v /(m mol )
-1 3 9 E
por el metanol y los steres de propilo considerados en este trabajo. Para la mezcla
Por ltimo, en relacin a los steres de butilo, las diferencias medias con respecto a los
resultados obtenidos por Ortega y col. (1985b) fueron superiores al 100%, para las mez-
clas del metanol con el metanoato de butilo, debido posiblemente a las impurezas del
metanoato de butilo empleado en ese momento y del 33% en las mezclas con el etanoato
nuestro laboratorio, y tampoco publicado todava, se obtiene una desviacin media del
-139-
Tabla (11-27). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2metanoato
de butilo a la temperatura de 298,15 K y coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0390 886,90 -16 0,5118 862,00 -106
0,0723 885,71 -27 0,5414 859,61 -108
0,1111 884,23 -37 0,5884 855,38 -99
0,1535 882,54 -51 0,6473 849,54 -95
0,1684 881,91 -55 0,6906 844,61 -82
0,2121 880,02 -70 0,7125 841,92 -75
0,2745 877,03 -84 0,7533 836,60 -70
0,2708 877,21 -83 0,7941 830,65 -59
0,3388 873,63 -99 0,8630 819,02 -43
0,4153 868,94 -104 0,8987 811,97 -32
0,4608 865,84 -107 0,9370 803,49 -21
kv=2,826 V0=-414 V1=-92 V2=279 s(vE)=3
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =888,21-32,61x1-28,44x12
[0,2-0,8] =888,88-41,08x1+8,98x12-61,09x13
[0,8-1,0] =1012,48-548,85x1+709,38x12-386,02x13
Tabla (11-28). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2etanoato de
butilo a la temperatura de 298,15 K y coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0617 874,00 -13 0,5647 850,80 -45
0,1193 872,27 -21 0,6045 847,68 -39
0,1202 872,22 -18 0,6414 844,55 -37
0,2107 869,17 -31 0,7204 836,73 -28
0,2356 868,23 -32 0,7809 829,47 -21
0,2843 866,32 -39 0,8003 826,86 -19
0,3738 862,25 -41 0,8559 818,42 -14
0,4599 857,65 -44 0,8923 811,91 -7
0,4906 855,78 -42 0,9466 800,66 -8
kv=3,258 V0=-202 V1=135 s(vE)=2
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =875,67-25,62x1-24,91x12
[0,2-0,8] =874,56-13,16x1-74,50x12+90,17x13-86,95x14
[0,8-1,0] =-498,84+6419,03x1-11393,16x12+8966,80x13-2707,12x14
-140-
Tabla (11-29). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2propanoa-to
de butilo a la temperatura de 298,15 K y coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0257 870,02 -8 0,6404 843,16 3
0,0843 868,62 -16 0,6822 839,55 8
0,1072 868,02 -17 0,7406 833,71 16
0,1650 866,41 -18 0,7994 826,66 21
0,2532 863,63 -20 0,8479 819,70 23
0,3159 861,38 -23 0,8851 813,47 22
0,4185 856,99 -19 0,9085 809,04 22
0,4751 854,13 -16 0,9208 806,55 20
0,5318 850,84 -8 0,9502 800,02 15
0,5840 847,44 -5 0,9771 793,23 9
kv=3,670 V0=-175 V1=586 s(vE)=2
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =870,57-21,12x1-24,95x12
[0,2-0,8] =867,03+12,83x1-142,47x12+186,17x13-133,98x14
[0,8-1,0] =1340,91-1739,01x1+2112,81x12-927,95x13
Tabla (11-30). Densidades, , volmenes de exceso, vE, para el sistema binario x1metanol+x2butanoato de
butilo a la temperatura de 298,15 K y coeficientes de correlacin.
x1 109vE x1 109vE
kgm-3 m3mol-1 kgm-3 m3mol-1
0,0471 863,44 1 0,5625 845,67 24
0,0766 862,82 2 0,6347 841,00 32
0,1643 860,81 2 0,6555 839,49 33
0,1956 860,01 4 0,6874 836,96 35
0,2376 858,86 6 0,7661 829,46 44
0,3395 855,66 12 0,8171 823,43 43
0,3511 855,26 11 0,8572 817,72 46
0,4204 852,62 14 0,8977 810,95 43
0,4725 850,34 18 0,9524 799,53 34
0,5227 847,84 24
kv=4,095 V0=20 V1=154 V2=753 s(vE)=2
Ecuaciones de ajuste densidad-composicin segn el intervalo de concentraciones
[0,0-0,2] =864,37-18,80x1-17,65x12
[0,2-0,8] =863,23-5,13x1-75,00x12+104,30x13-94,30x14
[0,8-1,0] =1154,11-1157,90x1+1507,03x12-716,30x13
-141-
Figura (11-17).Volmenes de exceso a298,15 K de las mezclas binarias
1
0.8
0.6
4
x1
no publicados
Ortega (1985)
0.4
0.2
0
-250
-350
60
10 v /(m mol )
-1 3 9 E
que se refleja en las Figuras (11-14) a la (11-17) no se puede establecer un anlisis pre-
ciso del comportamiento de esta magnitud para las mezclas de metanol+s-teres de al-
quilo, en funcin de la estructura de los compuestos que intervienen. Por ello, se realiz
una nueva representacin, la de la Figura (11-18), con los valores equimolares obtenidos
para cada una de las mezclas y la cadena hidrocarbonada del ster, R1COOR2, cuando,
vara tanto R1 como R2. Tambin se aprecian, en rojo, las diferencias con los publicados
en la bibliografa internacional.
Con la excepcin de las tres mezclas correspondientes a las Tablas (11-26), (11-
29) y (11-30), que presentan volmenes de exceso cercanos a cero y, por tanto, con com-
-142-
40 Figura (11-18). Representacin de los
volmenes de exceso a fraccin molar
u=4 intermedia, vE a x1=0,5, en funcin del
nmero de CH2 existentes en el ster
alqulico en las mezclas binarias:
{x1CH3OH+x2Cu-1H2u-1COOCvH2v+1} y
0
u=3 comparacin con los encontrados en la
literatura que se menciona:
a
10 v (a x1=0,5)/(m mol )
b
u=2 Grolier y Viallard (1971)
c
Nakanishi y Shirai (1970)
3
d
f Polak y Lu (1972a)
e
e
Polak y Lu (1972b)
f
Ortega y col. (1985)
g
a Ortega y col. (1986)
9 E
h
c Nikam y col. (1996)
d i
No-publicados
i
b Nota.- En esta grfica no se ha represen-
-100 b
tado el dato aportado por Ortega (1985)
i para el sistema x1metanol+ x2metanoato
u=1
i
h de butilo, ya que su valor, vE(a
x1=0,5)=-292, presenta una gran desvia-
cin con el experimental aqu aportado,
g
dificultando la visin de esta representa-
-150 cin.
1 2 v 3 4
son negativos, vE<0, o sea, con contracciones de volumen. Estas contracciones disminu-
yen en valor absoluto con el aumento de la cadena del ster, llegando a alcanzar valores
positivos, con expansiones, a altas concentraciones de metanol en los sistemas en los que
go de composiciones.
las mezclas de steres de alquilo con alcanoles, sin embargo, la presencia del metanol,
-143-
comportamiento. Son relevantes los efectos de asociacin de ambos tipos de molculas,
muy elevados en el metanol, mientras que para los steres disminuyen con el aumento de
la cadena o con el peso molecular; en algunos casos, sobre todo en las mezclas donde
intervienen los primeros de la serie, como los metanoatos, tambin la disociacin de am-
bos tipos de molculas permite la formacin de complejos por fuertes enlaces polares,
bin de acuerdo con el tamao y forma de las molculas, existen efectos fsicos, de impe-
moda entre las de ster con efectos de contraccin en el volumen de mezclado). Pues
mencionaba, las altas contracciones de volumen que son mayores de acuerdo con el or-
Preparacin de las tablas: Corresponde a esta seccin la presentacin de los datos obte-
(desde el metilo al butilo), a la presin constante de (141,30,2) kPa. Con el fin de esta-
-144-
blecer una presentacin ordenada, se ha credo conveniente mostrar la elaboracin de las
cada una de ellas, a los dos conjuntos de mezclas antes citados. Todos los datos, tanto los
posiciones de las fases lquida, x1, y vapor, y1, se tuvieron en cuenta todas las considera-
ciones mencionadas en la seccin 3.4. de esta Memoria y que ahora se describen nueva-
mente por su inters. Los coeficientes de actividad de los componentes en fase lquida,
vez, para conseguir dichos coeficientes de virial se utilizaron las correlaciones empricas
de Tsonopoulos (1974), seccin 3.4.2. Por ltimo, los volmenes molares, vio, y las pre-
por la ecuacin de Rackett modificada por Spencer y Danner (1972), que incluye las
toine, obtenida para cada compuesto participante de la mezcla, ver seccin 11-1. Pues
bien, una vez calculados los coeficientes de actividad, la ecuacin (4-2) permite conse-
guir los valores de la funcin de energa de Gibbs en su forma adimensional, gE/RT, para
-145-
cada fraccin molar y, lgicamente, los correspondientes valores absolutos, gE/(Jmol-1),
olvidarse que cada punto corresponde a un estado de equilibrio, o sea, a unas condicio-
que permiten construir las tablas (11-31a) a (11-47a). Las tres primeras columnas co-
constante p, la cuarta es la diferencia entre las concentraciones en las fases vapor y lqui-
do, (y1-x1), del compuesto de referencia y las cuatro restantes a las cantidades calculadas
como se han descrito, 1, 2, gE/RT y gE. En las figuras que aparecen bajo las tablas, des-
tablecidas para ello, ver seccin 4. En este trabajo se utilizaron algunos de los mtodos
propuestos por varios autores y que han sido desarrollados en dicha seccin, como es el
de Herington (1951), uno de los test de reas, y los test punto-a-punto de Van Ness y col.
cia de errores en el planteamiento del test de Herington, por su valoracin global de los
datos, este ser aqu tambin considerado, ya que el clsico test de reas o de Redlich-
-146-
adecuado para los caosos isobricos. En los apartados 11.4.1 y 11.4.2 se resumen los
los modelos matemticos existentes para correlacionar y reproducir los datos de equili-
as como la modificacin realizada para la ecuacin de las zetas, seccin 5.1.4. De esta
forma, todas ellas presentan en comn un hecho, que permiten correlacionar simult-
puedan conseguirse las entalpas de exceso para la mezcla en estudio. Como no es senci-
conocer hE(T), resulta algo difcil reproducir con la misma correlacin los datos de ELV
y hE. Debido al elevado nmero de coeficientes que presentan las mencionadas ecuacio-
lo que ocasiona que las correlaciones de datos de ELV isobrico, con la funcin gE/RT vs
x1 no reproduzcan bien las curvas de hE para un determinado sistema. Por ello, la corre-
lacin de datos puede llevarse a cabo siguiendo varios procedimientos, sin embargo, solo
b), utilizar datos de gE/RT y de hE/RT efectuando una correlacin simultnea de am-
bas magnitudes.
-147-
El primero de los casos, a), se utiliz para todas las ecuaciones, sin embargo, el
segundo, el b), solo fue considerado en la ecuacin de las zetas por la sencillez y me-
jor exactitud en los clculos con esta ecuacin, como se comprende facilmente. Dicho de
otra forma, los modelos de Wilson, NRTL y UNIQUAC dificultan la buena correlacin
simultnea de los datos de ELV y hE, dando lugar a elevados errores de ajuste en uno o
objetivo definida como proponen Holmes y van Winkle (1970) para el conjunto de los
datos de equilibrio:
( ) ( )
2 2
F .O.= 1,exp 1,cal + 2 ,exp 2 ,cal (11-2)
Sin embargo, para conseguir una ms ptima representacin de los datos de ELV
cin ocasiona prdidas de precisin. Por ello, a la condicin expresada por la ecuacin
(Q Qcal )
2
s( Q ) =
exp
(11-3)
N n
-148-
Las Tablas (11-31b) a la (11-47b) recogen los valores de los coeficientes obteni-
do para cada una de las ecuaciones indicadas en las tablas, excepto para la ecuacin de
las zetas para la cual se ha seguido el segundo de los mtodos indicados. Asmismo,
[ ]
1/ 2
s( i )= s( 1 ) 2 + s( 2 ) 2 (11-4)
se utiliz el procedimiento sealado como b), o sea, utilizando para cada mezcla los da-
tos de entalpa de exceso y de equilibrio lquido-vapor de manera conjunta. As, con los
18), A01 y A11 , los cuales, mediante la (5-20), permiten conseguir los de A01 y A11. Los
el caso de los n-alcanos, con una mayor desviacin de la idealidad, la F.O. requiere
Pues bien, con esta metodologa de clculo se han obtenido los valores que apa-
recen en las Tablas (11-31b) a la (11-47b) con expresin de los parmetros estadsticos
-149-
La presentacin y comparacin cualitativa de las correlaciones obtenidas se reali-
za en las Figuras (11-31c) a (11-47c). En cada caso, los smbolos indican los valores
experimentales y las lneas, representadas en distintos colores, negro, azul, rojo y verde,
las curvas calculadas, respectivamente, con las ecuaciones de las zetas, Wilson, NRTL y
UNIQUAC.
de los diferentes estados de equilibrio, como son las temperaturas y las concentraciones,
nes de las fases lquida y vapor y la temperatura normal de ebullicin del compuesto
Qk = x1 x2 Ai [x1 /( x1 + kx2 )]
i
(11-6)
las magnitudes del equilibrio tal como se ha indicado, se recogen tambin en las Tablas
(11-31b) a (11-47b) y todos ellos son obtenidos por un procedimiento de regresin li-
neal determinando el coeficiente k por optimizacin del ajuste. Insertadas en las Figuras
-150-
(11-31c) a (11-47c) se muestran otras grficas, donde se aprecian los valores experi-
ecuacin (11-6).
Una vez hecha la presentacin de cmo se han elaborado las tablas de datos expe-
tes apartados se reflejan de manera particular los dos conjuntos considerados, los siste-
vapor y otros calculados necesariamente para el anlisis de los sistemas binarios x1meta-
32a). En la misma pgina se muestra grficamente la situacin de los puntos de las mag-
-151-
nos comentarios sobre dichos resultados y su interpretacin.
mismo con el test de Fredenslund y col. (1977b). A pesar de los problemas asociativos
que puede presentar la presencia del metanol, no se observa en estos sistemas ninguna
en regiones extremas aparecen valores ligeramente inferiores a la unidad nica seal que
indica la aparicin de algunos efectos consecuencia del proceso de mezclado con el al-
tico de las magnitudes que definen el equilibrio isobrico lquido-vapor, presentando los
resultados de las relaciones utilizadas, ecuacin de las zetas, Wilson, NRTL y UNI-
de las distintas correlaciones. El ajuste con la ecuacin de las zetas para estas mezclas
de gE/RT y hE, y fue por ello necesario introducir un nuevo trmino para el ajuste de en-
talpas, lo que di lugar a dos nuevos coeficientes en la expresin final de gE/RT vs x1,
-152-
cuyo desarrollo es bastante simple. En la siguiente tabla se presentan los coeficientes
obtenidos para el ajuste de las hE vs x1 con los datos sacados de la literatura. La ltima
columna representa el error medio conseguido al recalcular las entalpas con la ecuacin
cin de los datos experimentales utilizando la ecuacin de las zetas. Para el sistema
QUAC. Cabe sealar que las altas desviaciones estndar promedio de los coeficientes de
actividad en ambos sistemas binarios se debe a los elevados valores que muestran los
En estas mismas tablas aparecen los coeficientes obtenidos para los ajustes de las
unos puntos azeotrpicos muy acusados, se consider conveniente realizar las correla-
ciones utilizando dos intervalos, separados por el punto singular, xaz. En las figuras in-
-153-
Tabla (11-31a) . Magnitudes del ELV a 141,3 kPa para la mezcla binaria x1metanol+x2n-pentano.
30
0 0.5
g /RT
y1-x1
T/K
i
E
0.2
10
-0.6
Figura (11-31). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mezcla
binaria x1metanol+x2n-pentano. (a) T o (y1-x1) vs x1 ; (b) gE/RT o i vs x1.
-154-
Tabla (11-31b). Correlaciones del ELV a 141,3 kPa para la mezcla binaria x1metanol+x2n-pentano.
0
30
1
T/K
y -x
1
1.0
-0.6
g /RT
i
E
310 -0.8
0 0.2 0.4 x1 0.6 0.8 1
0.5
10
0 0.9
0 0.2 0.4 x1 0.6 0.8 1
Figura (11-31c). Correlacin de los datos de equilibrio para el sistema x1metanol+x2n-pentano a 141,3
kPa. Para gE/RT o i vs x1: ecuacin de las zetas, Wilson,, NRTL, UNIQUAC.
Figura insertada: T o (y1-x1) vs x1; curva calculada con la ecuacin de las zetas.
-155-
Tabla (11-32a). Magnitudes del ELV a 141,3 kPa para la mezcla binaria x1metanol+x2n-hexano.
30
0.3
0.5
g /RT
y1-x1
T/K
i
E
0.2
10
-0.2
Figura (11-32). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la
mezcla binaria x1metanol+x2n-hexano. (a) T o (y1-x1) vs x1 ; (b) gE/RT o i vs x1.
-156-
Tabla (11-32b). Correlaciones del ELV a 141,3 kPa para la mezcla binaria x1metanol+x2n-hexano.
0.3
30
1
T/K
y -x
1
1.0
-0.2
g /RT
i
E
330 -0.4
0 0.2 0.4 x1 0.6 0.8 1
0.5
10
0 0.9
0 0.2 0.4 x1 0.6 0.8 1
Figura (11-32c). Correlacin de los datos de equilibrio para el sistema x1metanol+x2n-hexano a 141,3
kPa. Para gE/RT o i vs x1: ecuacin de las zetas, Wilson,, NRTL, UNIQUAC.
Figura insertada: T o (y1-x1) vs x1; curva calculada con la ecuacin de las zetas.
-157-
11.4.2. Datos de ELV para las mezclas metanol+steres de alquilo
(141,3 0,2) kPa, de los quince sistemas binarios x1metanol+x2steres de alquilo aqu
Para estas mezclas aparece, a baja concentracin del alcanol, unos mximos en los i que
efecto qumico que se traduce sobre los coeficientes de actividad. La primera curvatura,
ma casi general, se aprecia que los mximos en los coeficientes de actividad estn en las
mezclas con los steres de cadena larga, indicando que, debido a la diferencia de tamao,
existen interacciones por las caractersticas estructurales distintas de los compuestos que
intervienen, y como consecuencia de ello aparece una inflexin en las energas de Gibbs
de exceso.
De acuerdo con el orden establecido en esta seccin, ahora se presentan los re-
-158-
del metanol con los n-alcanos y los steres de alquilo:
do la ecuacin de las zetas y los modelos Wilson, NRTL y UNIQUAC. Los coeficientes
en las Tablas (11-33b) a (11-47b). Adems del anlisis cuantitativo, en cuanto a la pre-
-159-
cisin de las correlaciones efectuadas, se puede realizar una comparacin visual de la
les de las entalpas de exceso para las mezclas del metanol con los quince steres de al-
quilo, que, en su mayor parte fueron determinados en nuestro laboratorio para completar
este trabajo y por tanto, no-publicados en la literatura. Todos ellos fueron medidos a una
La ltima columna de esta tabla presenta los errores medios obtenidos, en por-
centaje, entre la curva experimental que corresponde a la mostrada aqu y la que resulta
-160-
multnea los datos de la funcin adimensional de energa libre y los coeficientes de acti-
vidad, una vez determinado los valores de A01 y A11 a partir de las entalpas de mezcla.
Los coeficientes obtenidos en este procedimiento son los que se recogen en las Tablas
(11-33b) a (11-47b). En las quince mezclas que constituyen este segundo grupo o fami-
lia, no fue necesario como para las mezclas con n-alcanos, por la forma mas regular de
observa la existencia de una mejor correlacin de los datos con la ecuacin de las zetas,
siendo muy similares entre s las ofrecidas por los otros tres modelos, Wilson, NRTL y
cidad de los modelos de composicin local para representar de forma adecuada el com-
portamiento del ELV. Asimismo, para la correlacin de las magnitudes directas, tempe-
do en los sistemas x1metanol+x2n-alcano. En estos casos, para las mezclas que presenta-
ban con azetropo no fue necesario realizar los ajustes en distintos tramos por la forma
muestran los parmetros y las desviaciones estndar obtenidas en cada caso. En las Fi-
-161-
Tabla (11-33a). Magnitudes del ELV a 141,3 kPa para el sistema x1metanol+x2etanoato de metilo.
2.5
0 0.2
g /RT
y1-x1
T/K
i
E
-0.1 0.1
1.5
1
330 -0.2 0 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (11-33). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mezcla
binaria x1metanol+x2etanoato de metilo. (a) T o (y1-x1) vs x1; (b) gE/RT o i vs x1.
-162-
Tabla (11-33b). Correlaciones del ELV a 141,3 kPa para el sistema x1metanol+x2etanoato de metilo.
0.4
1
T/K
y -x
2.5
1
-0.1
330 -0.2
g /RT
i
E
1.5
0.1
1
0 0.9
0 0.2 0.4 x1 0.6 0.8 1
Figura (11-33c). Correlacin de los datos de equilibrio para el sistema x1metanol+x2etanoato de metilo a
141,3 kPa. Para gE/RT o i vs x1: ecuacin de las zetas, Wilson, NRTL, UNIQUAC.
Figura insertada: T o (y1-x1) vs x1; curva calculada con la ecuacin de las zetas.
-163-
Tabla (11-34a). Magnitudes del ELV a 141,3 kPa para el sistema x1metanol+x2propanoato de metilo.
2.5
0.2
g /RT
y1-x1
T/K
i
E
0.1
1.5
1
342 -0.05 0 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (11-34). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mez-
cla binaria x1metanol+x2popanoato de metilo. (a) T o (y1-x1) vs x1; (b) gE/RT o i vs x1.
-164-
Tabla (11-34b). Correlaciones del ELV a 141,3 kPa para el sistema x1metanol+x2propanoato de metilo.
1
T/K
y -x
2.5
342 -0.05
g /RT
i
E
0.2
1.5
0.1
1
0 0.9
0 0.2 0.4 x1 0.6 0.8 1
Figura (11-34c). Correlacin de los datos de equilibrio para el sistema x1metanol+x2propanoato de metilo
a 141,3 kPa. Para gE/RT o i vs x1: ecuacin de las zetas, Wilson, NRTL, UNIQUAC.
Figura insertada: T o (y1-x1) vs x1; curva calculada con la ecuacin de las zetas.
-165-
Tabla (11-35a). Magnitudes del ELV a 141,3 kPa para el sistema x1metanol+x2butanoato de metilo.
0.2
2.7
g /RT
y1-x1
T/K
i
E
0.05
1.5
1
340 0 -0.05 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (11-35). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mez-
cla binaria x1metanol+x2butanoato de metilo. (a) T o (y1-x1) vs x1; (b) gE/RT o i vs x1.
-166-
Tabla (11-35b). Correlaciones del ELV a 141,3 kPa para el sistema x1metanol+x2butanoato de metilo.
1
T/K
y -x
1
0.4 2.7
g /RT
340 0 i
E
0.1 1.5
1
-0.05 0.9
0 0.2 0.4 x1 0.6 0.8 1
Figura (11-35c). Correlacin de los datos de equilibrio para el sistema x1metanol+x2butanoato de metilo a
141,3 kPa. Para gE/RT o i vs x1: ecuacin de las zetas, Wilson,, NRTL, UNIQUAC.
Figura insertada: T o (y1-x1) vs x1; curva calculada con la ecuacin de las zetas.
-167-
Tabla (11-36a). Magnitudes del ELV a 141,3 kPa para el sistema x1metanol+x2metanoato de etilo.
0 0.2 3
g /RT
1
T/K
y -x
i
E
-0.1 0.1
1.5
1
330 -0.2 0 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (11-36). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mezcla
binaria x1metanol+x2metanoato de etilo. (a) T o (y1-x1) vs x1 ; (b) gE/RT o i vs x1.
-168-
Tabla (11-36b). Correlaciones del ELV a 141,3 kPa para el sistema x1metanol+x2metanoato de etilo.
0.4
1
T/K
y -x
1
-0.1
2.5
330 -0.2
0 0.2 0.4 x1 0.6 0.8 1
g /RT
i
E
1.5
0.1
1
0 0.9
0 0.2 0.4 x1 0.6 0.8 1
Figura (11-36c). Correlacin de los datos de equilibrio para el sistema x1metanol+x2metanoato de etilo a
141,3 kPa. Para gE/RT o i vs x1: ecuacin de las zetas, Wilson, NRTL, UNIQUAC.
Figura insertada: T o (y1-x1) vs x1; curva calculada con la ecuacin de las zetas.
-169-
Tabla (11-37a). Magnitudes del ELV a 141,3 kPa para el sistema x1metanol+x2etanoato de etilo.
0.2 3
g /RT
1
T/K
y -x
i
1
0.1
1.5
0
1
342 -0.05 0 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (11-37). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mez-
cla binaria x1metanol+x2etanoato de etilo. (a) T o (y1-x1) vs x1 ; (b) gE/RT o i vs x1.
-170-
Tabla (11-37b). Correlaciones del ELV a 141,3 kPa para el sistema x1metanol+x2etanoato de etilo.
0.4
1
T/K
y -x
1
0 2.5
g /RT
342 -0.05
0 0.2 0.4 x1 0.6 0.8 1
i
E
1.5
0.1
1
0 0.9
0 0.2 0.4 x1 0.6 0.8 1
Figura (11-37c). Correlacin de los datos de equilibrio para el sistema x1metanol+x2etanoato de etilo a
141,3 kPa. Para gE/RT o i vs x1: ecuacin de las zetas, Wilson, NRTL, UNIQUAC.
Figura insertada: T o (y1-x1) vs x1; curva calculada con la ecuacin de las zetas.
-171-
Tabla (11-38a). Magnitudes del ELV a 141,3 kPa para el sistema x1metanol+x2propanoato de etilo.
0.3 0.2
g /RT
y1-x1
T/K
i
E
0.1
0.1
1
340 0 0 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (11-38). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mezcla
binaria x1metanol+x2propanoato de etilo. (a) T o (y1-x1) vs x1 ; (b) gE/RT o i vs x1.
-172-
Tabla (11-38b). Correlaciones del ELV a 141,3 kPa para el sistema x1metanol+x2propanoato de etilo.
0.3
0.4
1
T/K
y -x
1
0.1
2.5
g /RT
340 0
0 0.2 0.4 x1 0.6 0.8 1 i
E
1.5
0.1
1
0 0.9
0 0.2 0.4 x1 0.6 0.8 1
Figura (11-38c). Correlacin de los datos de equilibrio para el sistema x1metanol+x2propanoato de etilo a
141,3 kPa. Para gE/RT o i vs x1: ecuacin de las zetas, Wilson, NRTL, UNIQUAC.
Figura insertada: T o (y1-x1) vs x1; curva calculada con la ecuacin de las zetas.
-173-
Tabla (11-39a). Magnitudes del ELV a 141,3 kPa para el sistema x1metanol+x2butanoato de etilo.
0.4 0.2
g /RT
1
T/K
y -x
i
E
0.2 0.1
1
340 0 0 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (11-39). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mezcla
binaria x1metanol+x2butanoato de etilo. (a) T o (y1-x1) vs x1 ; (b) gE/RT o i vs x1.
-174-
Tabla (11-39b). Correlaciones del ELV a 141,3 kPa para el sistema x1metanol+x2butanoato de etilo.
0.4
0.4
1
T/K
y -x
1
0.2 3.5
340 0
0 0.2 0.4 x1 0.6 0.8 1
g /RT
i
E
0.1
1.5
1
0 0.9
0 0.2 0.4 x1 0.6 0.8 1
Figura (11-39c). Correlacin de los datos de equilibrio para el sistema x1metanol+x2butanoato de etilo a
141,3 kPa. Para gE/RT o i vs x1: ecuacin de las zetas, Wilson, NRTL, UNIQUAC.
Figura insertada: T o (y1-x1) vs x1; curva calculada con la ecuacin de las zetas.
-175-
Tabla (11-40a). Magnitudes del ELV a 141,3 kPa para el sistema x1metanol+x2metanoato de propilo.
0.2 3.5
0.2
g /RT
y1-x1
T/K
i
E
0.1
1.5
0
1
340 -0.05 0 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (11-40). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mezcla
binaria x1metanol+x2metanoato de propilo. (a) T o (y1-x1) vs x1; (b) gE/RT o i vs x1.
-176-
Tabla (11-40b). Correlaciones del ELV a 141,3 kPa para el sistema x1metanol+x2metanoato de propilo.
0.2
1
T/K
y -x
1
3.5
0.4
0
g /RT
340 -0.05 i
E
0.2
1.5
1
0 0.9
0 0.2 0.4 x1 0.6 0.8 1
-177-
Tabla (11-41a). Magnitudes del ELV a 141,3 kPa para el sistema x1metanol+x2etanoato de propilo.
2.5
0.2
0.3
g /RT
y1-x1
T/K
i
E
0.15
1.5
1
345 0 -0.05 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (11-41). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mez-
cla binaria x1metanol+x2etanoato de propilo. (a) T o (y1-x1) vs x1; (b) gE/RT o i vs x1.
-178-
Tabla (11-41b). Correlaciones del ELV a 141,3 kPa para el sistema x1metanol+x2etanoato de propilo.
0.3
0.4
1
T/K
y -x
1
0.15
345 0
0 0.2 0.4 x1 0.6 0.8 1
g /RT
i
E
0.1
0
1
-0.05 0.9
0 0.2 0.4 x1 0.6 0.8 1
Figura (11-41c). Correlacin de los datos de equilibrio para el sistema x1metanol+x2etanoato de propilo a
141,3 kPa. Para gE/RT o i vs x1: ecuacin de las zetas, Wilson, NRTL,UNIQUAC.
Figura insertada: T o (y1-x1) vs x1; curva calculada con la ecuacin de las zetas.
-179-
Tabla (11-42a). Magnitudes del ELV a 141,3 kPa para el sistema x1metanol+x2propanoato de propilo.
0.2
0.4
g /RT
1
T/K
y -x
i
E
0.2
1.5
0
1
340 0 -0.02 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (11-42). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mezcla
binaria x1metanol+x2propanoato de propilo. (a) T o (y1-x1) vs x1; (b) gE/RT o i vs x1.
-180-
Tabla (11-42b). Correlaciones del ELV a 141,3 kPa para la mexcla x1metanol+x2 propanoato de propilo.
0.4
0.4
1
T/K
y -x
1
0.2
340 0
0 0.2 0.4 x1 0.6 0.8 1
g /RT
i
E
0.1
1.5
0 1
-0.02 0.9
0 0.2 0.4 x1 0.6 0.8 1
-181-
Tabla (11-43a). Magnitudes del ELV a 141,3 kPa para el sistema x1metanol+x2butanoato de propilo.
0.5 0.2
2.5
g /RT
y1-x1
T/K
i
E
0.2
1.5
1
340 0 -0.05 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (11-43). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mez-
cla binaria x1metanol+x2butanoato de propilo. (a) T o (y1-x1) vs x1; (b) gE/RT o i vs x1.
-182-
Tabla (11-43b). Correlaciones del ELV a 141,3 kPa para el sistema x1metanol+x2butanoato de propilo.
0.5
1
T/K
y -x
1
0.4
0.2
2.5
g /RT
340 0
0 0.2 0.4 0.6 0.8 1 i
E
x1
0.2
1.5
0
1
-0.05 0.9
0 0.2 0.4 x1 0.6 0.8 1
Figura (11-43c). Correlacin de los datos de equilibrio para el sistema x1metanol+x2butanoato de propilo
a 141,3 kPa. Para gE/RT o i vs x1: ecuacin de las zetas, Wilson, NRTL, UNIQUAC.
Figura insertada: T o (y1-x1) vs x1; curva calculada con la ecuacin de las zetas.
-183-
Tabla (11-44a). Magnitudes del ELV a 141,3 kPa para el sistema x1metanol+x2metanoato de butilo.
0.4
0.2 3
g /RT
y1-x1
T/K
i
E
0.1 1.5
0
1
340 0 -0.02 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (11-44). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mezcla
binaria x1metanol+x2metanoato de butilo. (a) T o (y1-x1) vs x1; (b) gE/RT o i vs x1.
-184-
Tabla (11-44b). Correlaciones del ELV a 141,3 kPa para el sistema x1metanol+x2metanoato de butilo
0.4
1
T/K
y -x
1
0.4 3
0.1
g /RT
i
E
340 0
0 0.2 0.4 x1 0.6 0.8 1
0.2
1.5
0 1
-0.02 0.9
0 0.2 0.4 x1 0.6 0.8 1
Figura (11-44c). Correlacin de los datos de equilibrio para el sistema x1metanol+x2metanoato de butilo a
141,3 kPa. Para gE/RT o i vs x1: ecuacin de las zetas, Wilson, NRTL, UNIQUAC.
Figura insertada: T o (y1-x1) vs x1; curva calculada con la ecuacin de las zetas.
-185-
Tabla (11-45a). Magnitudes del ELV a 141,3 kPa para el sistema x1metanol+x2etanoato de butilo.
2.5
0.2
0.4
g /RT
y1-x1
T/K
i
E
0.2
1.5
0
1
340 0 -0.02 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (11-45). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mez-
cla binaria x1metanol+x2etanoato de butilo. (a) T o (y1-x1) vs x1; (b) gE/RT o i vs x1.
-186-
Tabla (11-45b). Correlaciones del ELV a 141,3 kPa para el sistema x1metanol+x2etanoato de butilo.
0.4
1
T/K
y -x
1
0.4 0.2
2.5
g /RT
340 0
i
E
0.2
1.5
0 1
-0.02 0.9
0 0.2 0.4 x1 0.6 0.8 1
Figura (11-45c). Correlacin de los datos de equilibrio para el sistema x1metanol+x2etanoato de butilo a
141,3 kPa. Para gE/RT o i vs x1: ecuacin de las zetas, Wilson, NRTL, UNIQUAC.
Figura insertada: T o (y1-x1) vs x1; curva calculada con la ecuacin de las zetas.
-187-
Tabla (11-46a). Magnitudes del ELV a 141,3 kPa para el sistema x1metanol+x2propanoato de butilo.
0.5
0.2 2.5
g /RT
1
T/K
y -x
i
E
0.1
0.2
1.5
0
1
340 0 -0.03 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (11-46). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mez-
cla binaria x1metanol+x2propanoato de butilo. (a) T o (y1-x1) vs x1; (b) gE/RT o i vs x1.
-188-
Tabla (11-46b). Correlaciones del ELV a 141,3 kPa para el sistema x1metanol+x2propanoato de butilo.
0.5
1
T/K
y -x
1
0.4 0.2
2.5
340 0
g /RT
i
E
0.2
1.5
0 1
-0.03 0.9
0 0.2 0.4 x1 0.6 0.8 1
-189-
Tabla (11-47a). Magnitudes del ELV a 141,3 kPa para el sistema x1metanol+x2butanoato de butilo.
0.3
0.6
2.5
g /RT
y1-x1
T/K
i
E
0.1
0.2
1.5
0
1
340 0 -0.03 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (11-47). Representacin de los valores experimentales obtenidos en el ELV a 141,3 kPa de la mez-
cla binaria x1metanol+x2butanoato de butilo. (a) T o (y1-x1) vs x1; (b) gE/RT o i vs x1.
-190-
Tabla (11-47b). Correlaciones del ELV a 141,3 kPa para el sistema x1metanol+x2butanoato de butilo.
0.6 0.6
1
T/K
y -x
1
3
0.2
g /RT
i
340 0
0 0.2 0.4 0.6 0.8 1
E
x1
0.2
1.5
0 1
-0.03 0.9
0 0.2 0.4 x1 0.6 0.8 1
Figura (11-47c). Correlacin de los datos de equilibrio para el sistema x1metanol+x2butanoato de butilo a
141,3 kPa. Para gE/RT o i vs x1: ecuacin de las zetas, Wilson, NRTL, UNIQUAC.
Figura insertada: T o (y1-x1) vs x1; curva calculada con la ecuacin de las zetas.
-191-
11.5. Azetropos
azetropos en las mezclas binarias son bien conocidas, pero nos parece acertada la pre-
sentacin aqu de las mismas, separadas formalmente tal como lo realiza Hiaki y col.
-192-
ya mencionadas, cuyos coeficientes aparecen las Tablas (11-31b) a (11-47b) se obtie-
nen, para cada mezcla, un conjunto de ecuaciones cuya resolucin analtica proporciona
los valores de Taz, y de xaz, correspondientes al punto singular. As por la identidad de las
ecuaciones de correlacin de temperaturas con las composiciones de las fases vapor, Ty,
Ty T - y1Tbo,1 y2Tbo, 2
=1 =1 (11-12)
Tx T - x1Tbo,1 x2Tbo, 2
surgiendo tambin de aqu un valor de x1 que verifique la citada condicin. Por ltimo, la
dTx x2 2
=0 Tb,1
o
Tbo,2 +( 1 2 x1 ) + kz = 0 (11-14)
dx1 x1
-193-
donde: =(Aizi)x y =(d/dz)
y cuarta, las otras dos verifican con mayor precisin el clculo de la composicin y tem-
peratura del azetropo. Los resultados analticos deben ser, al menos en la prctica, coin-
las publicaciones realizadas por Horsley (1973) y ms recientemente por Gmehling y col.
trado valores a la presin de trabajo para los sistemas aqu considerados. No obstante,
-194-
El anlisis de la situacin de los puntos azeotrpicos, comparndolos con otros
relacionan todos los valores encontrados en la literatura para las mezclas de n-alcanos y
la que pertenecen.
Tabla (11-48). Coordenadas azeotrpicas estradas de la literatura para la mezclas binarias formadas por
x1metanol+x2n-alcanos.
-195-
0 0
(a) (b) 30
30
9 9
log p r =-3.153+8.332x1-4.389x1
2
30 30
-0.5 9 -0.5 9
30 30
log p ro=-7.519+15.346x1-7.363x1
2
9 9
n-hexano
o
o
log p ro
log p r
16
16
16 16
experimental
16
3 n-pentano 3
16 1 16
16 1
metanol 16
-2 -2 10
10
10 10
T r=0.7 10 10
-2.5 -2.5
1.1 1.3 1.6 1.8 0.1 0.3 0.6 0.8
1/T r x1
Figura (11-48). Representacin de: (a), log pro vs 1/Tr y (b), log pro vs x1 para los puntos azeotrpicos de las mezclas x1metanol+x2n-pentano () o
x1metanol+x2n-hexano () y las correspondientes correlaciones en coordenadas reducidas (a), para x1metanol+x2n-pentano, log pro=3,276-3,106/Tr; para
x1metanol+x2n-hexano, log pro=3,682-3,424/Tr.
-196-
En la Figura (11-48a,b) se han representado, tanto los valores experimentales
las presiones de vapor de las sustancias puras, ver tambin Figura (11-11a), como de los
puntos azeotrpicos, (x1,log pro) y (1/Tr,log pro). En principio, de acuerdo con las grficas
y las cifras tabuladas, los valores obtenidos en este trabajo parecen concordantes con los
de algunos otros autores, aunque sea en diferentes condiciones de las magnitudes medi-
bles intensivas. En las grficas (b) se muestran tambin las correlaciones que se pueden
mar, por extrapolacin, las coordenadas azeotrpicas cuando no se dispone de datos ex-
perimentales. Como ejemplo, los valores calculados para la presin de trabajo, 141,3
Para el conjunto de las quince mezclas formadas por metanol y los distintos ste-
-197-
res de alquilo, tan slo cinco de ellas presentaron azetropo a la presin de trabajo. Di-
(1994) ofrece numerosas referencias acerca de la localizacin del azetropo en los siste-
la prediccin realizada por Lecat (1926), ver Gmehling y col. (1994), a 101,3 kPa. En
nadas reducidas (1/Tr,log pr) y (x1,log pr) aparece la comparacin de todos los datos re-
cogidos de la literatura junto con los determinados para este trabajo, observndose una
adecuada concordancia entre ellos. Las correlaciones incluidas en la segunda de las figu-
ras, Figura (11-49b), permiten estimar el punto singular en otras condiciones diferentes.
Como ejemplo, valga decir que para la presin considerada aqu se obtiene, en el clculo
-198-
Tabla (11-49). Coordenadas azeotrpicas de las mezclas binarias x1metanol+x2steres de alquilo en literatrura.
p/kPa T/K x1 ref. p/kPa T/K x1 ref. p/kPa T/K x1 ref.
x1metanol+x2etanoato de metilo
13,3 281,7 0,228 1 61,9 313,2 0,346 13 101,3 327,0 0,332 17
24,6 293,2 0,288 2 62,6 314,1 0,304 3 101,3 327,0 0,330 18
26,7 294,7 0,263 1 66,7 316,1 0,302 14 101,3 327,0 0,328 19
27,6 295,2 0,253 3 66,7 316,5 0,290 10 101,3 327,1 0,332 20
31,5 298,2 0,260 4 70,1 316,9 0,312 3 101,3 327,2 0,352 21
32,0 298,2 0,248 5 72,8 318,2 0,323 6 101,3 327,2 0,350 22
34,1 299,8 0,267 3 73,3 318,2 0,310 4 101,3 327,2 0,344 23
39,4 303,2 0,280 2 80,0 320,4 0,322 9 101,3 327,2 0,337 24
39,7 303,3 0,275 3 80,0 321,0 0,306 10 101,3 327,2 0,337 25
40,0 303,6 0,280 1 82,1 321,0 0,323 3 101,3 327,2 0,323 26
48,4 307,9 0,288 3 86,7 323,1 0,310 10 101,3 327,2 0,320 27
48,5 308,2 0,290 6 88,5 322,9 0,340 11 101,3 327,3 0,322 10
48,8 308,2 0,285 4 88,5 323,2 0,343 11 106,7 328,8 0,327 10
49,5 308,2 0,296 7 89,2 324,4 0,332 3 107,2 328,2 0,348 7
52,9 309,9 0,301 8 93,3 325,0 0,316 10 395,4 370,2 0,462 28
53,3 310,2 0,305 9 101,3 326,6 0,340 7 445,8 372,2 0,486 29
53,3 310,4 0,300 1 101,3 326,7 0,336 3 585,7 382,9 0,485 26
53,3 310,9 0,273 10 101,3 326,7 0,338 8 787,5 396,4 0,529 28
56,0 311,4 0,298 3 101,3 326,7 0,352 1 790,3 393,2 0,550 29
59,4 313,2 0,309 2 101,3 326,9 0,350 15 877,5 398,3 0,545 26
59,6 312,9 0,340 11 101,3 327,0 0,347 16 1134,8 408,2 0,611 29
59,6 313,2 0,337 12 101,3 327,0 0,332 9 1172,3 410,2 0,585 26
x1metanol+x2propanoato de metilo
19,7 298,2 0,647 30 101,3 335,2 0,733 31 101,3 335,6 0,713 16
x1metanol+x2metanoato de etilo
53,3 307,3 0,270 32 81,3 318,2 0,300 33 101,3 323,7 0,316 32
66,7 312,8 0,286 32 86,7 319,5 0,301 32 106,7 325,1 0,321 32
80,0 317,4 0,296 32 93,3 321,5 0,308 32 - - - -
x1metanol+x2etanoato de etilo
13,3 291,1 0,608 33 75,8 328,2 0,693 36 101,3 335,4 0,715 38
19,9 298,2 0,634 34 80,0 329,3 0,695 32 101,3 335,4 0,684 16
26,7 304,7 0,643 33 86,7 331,2 0,700 32 101,3 335,5 0,710 39
40,0 312,9 0,675 11 93,3 333,1 0,704 32 101,3 335,5 0,710 22
40,3 313,2 0,630 13 92,1 333,2 0,699 35 101,3 335,4 0,718 33
40,9 313,2 0,650 35 100,9 334,3 0,709 24 101,3 335,5 0,707 35
53,3 319,4 0,670 32 101,3 335,2 0,709 25 101,3 335,6 0,689 40
61,9 323,2 0,686 35 101,3 335,2 0,708 32 101,3 335,7 0,712 41
66,7 324,5 0,694 33 101,3 335,2 0,710 37 101,3 336,0 0,725 42
66,7 324,6 0,682 32 - - - - 106,7 336,5 0,711 32
x1metanol+x2metanoato de propilo
101,3 335,1 0,735 16 - - - - - - - -
1 Serafimov y Balashov (1966) 15 Bushmakin y Kish (1957) 29 Horsley (1973)
2 Bekarek (1968) 16 Lecat (1949) 30 Polak y Lu (1972)
3 Figurski y von Weber (1991) 17 Kudryavtseva y col. (1966) 31 Eastman Data (1973)
4 Figurski y von Weber (1976) 18 Mato y Cepeda (1984) 32 Lee y col. (1986)
5 Martin y col. (1994) 19 Nagata (1969) 33 Park y col. (1973)
6 Nagata (1972) 20 Crawford y col (1949) 34 Staroske y Figursky (1992)
7 Vasileva y col. (1983) 21 Wuyts (1924) 35 Murti y van Winkle (1958)
8 Gmehling (1994) 22 Tanaka (1985) 36 Nagata y col. (1975)
9 Dobroserdov y Bagrov (1967) 23 Hannotte (1926) 37 van Zandijcke y col. (1974)
10 Lee y col. (1986) 24 Ryland (1899) 38 Nakanishi y col. (1967)
11 Bredig y Bayer (1927b) 25 Durrans (1920) 39 Akita y Yoshida (1963)
12 Bredig y Bayer (1927a) 26 Nagahama e Hirata (1971) 40 Rajendran y col. (1988)
13 Litvinov (1952) 27 Wang y Cai (1993) 41 Nagata (1962)
14 Nagata y Ohta (1971) 28 Balashov y col. (1967) 42 Ohe y col. (1971)
Recopilaciones de (1-4,6-28,30,32-39,41,42) en Gmehling y col. (1994); (29,31) en Horsley (1973)
-199-
-0.5 -0.5
(a) (b)
log p or =-4.746+11.409x 1-7.850x21
26
29
26
28
29
metanoato de propilo 26
etanoato de metilo
propanoato de metilo 29
etanoato de etilo 28
-1.3 -1.3
experimental
2
metanoato etilo
log p r o=-4.039+8.356x1-3.690x1
o
log p ro
log p r
32
32
log p r o=-21.505+48.276x1-29.076x12
16
32
log p r =-18.68+42.20x1-25.81x12
metanol 32
32
-2.0 -2 32
35
33
30
o
T r=0.7
33
1
-2.8 -2.8
1.2 1.3 1.6 1.8 0.2 0.3 0.6 0.8
1/T r x1
Figura (11-49). Representacin de: (a), log pro vs 1/Tr y (b), log pro vs x1 para los azeotropos de las mezclas: x1metanol+x2etanoato de metilo (!); x1metanol
+x2propanoato de metilo (); x1metanol+x2metanoato de etilo (); x1metanol+x2etanoato de etilo (); x1metanol+x2metanoato de propilo () y las lneas
de ajuste. Correlaciones en (a): log pr=a+b/Tr, para x1metanol+x2etanoato de metilo, log pr=3,469-3,397/Tr; para x1metanol +x2propanoato de metilo, log
pr=3,956-3,697/Tr; para x1metanol+x2metanoato de etilo, log pr=3,328-3,271/Tr; para x1metanol+x2etanoato de etilo, log pr=3,996 -3,739/Tr. En (a) se incluyen
las lneas de algunos compuestos puros, indicadas tambin en la Figura (11-11).
-200-
error inferior al 0,3 % en todos los casos; sin embargo, la estimacin de las composicio-
nes de los azetropos se obtienen con errores variables, desde un 3% para el sistema
x1metanol+x2etanoato de metilo hasta un valor inferior al 0,4% para los otros sistemas
considerados
141,3 kPa; por ello, la comparacin que se pretende se ve limitada a un anlisis cualita-
res a esta, en aquellas figuras que fueron mostradas en el capitulo primero, concreta-
sultado final es la Figura (12-1), donde se destaca que los valores obtenidos para la
mezcla binaria x1metanol+x2n-pentano son algo inferiores a los presentados por Tenn y
Missen (1963) a 99,99 kPa y por Budantseva y col. (1975), ver Gmehling y col. (1988)
a 101,32 kPa, si bien la presin experimental es superior a la de ambos autores. Sin em-
-201-
bargo, para la mezcla x1metanol+x2n-hexano los datos experimentales son algo superio-
res a los de Budantseva y col. (1975), ver Gmehling y col. (1988), y Raal y col. (1972) a
101,32 kPa; con los de Vilim (1961), ver Gmehling y Onken (1977), a 99,99 kPa, se
metanol.
recogidos en la Figura (1-2) deben ser reducidos con el fin que puedan introducirse los
clarificando as la visualizacin. En este caso, los valores del ELV y las curvas experi-
0.5
Figura (12-1). Comparacin de las curvas expe-
n-hexano rimentales y de la bibliografa para las cantidades
(y1-x1) vs x1 obtenidas para las mezclas x1metanol
+x2n-alcanos.
0.2 x1metanol+x2n-pentano
------
Tenn y Missen (1963), a 99,99 kPa
Budantseva y col. (1975), a 101,32 kPa
0 (ver Gmehling y col. 1988)
experimentales, a 141,3 kPa
y1-x1
x1metanol+x2n-hexano
-202-
0.1 0.25
585.7 kPa 101.3 kPa
395.4 kPa
0.05
0 0.15
y1-x1
y1-x1
101.3 kPa
66.7 kPa
-0.15 0
80 kPa
26.7 kPa
53.3 kPa
Figura (12-2). Comparacin de las cantidades Figura (12-3). Comparacin de las cantidades
(y1-x1) experimentales y bibliogrficas para la (y1-x1) experimentales y bibliogrficas para la
mezcla x1metanol+x2etanoato de metilo. mezcla x1metanol+x2etanoato de etilo.
mentales obtenidas aqu, siguen la variacin que se observa con los datos de otros auto-
presin. El sistema reportado por Balashov y col. (1967), ver Gmehling y Onken
(1977), a 395,4 kPa muestra una fuerte anomala con respecto al orden establecido con
los de otros autores, sobre todo a bajas concentraciones de metanol, presentando tam-
bin una consistencia negativa con el test punto a punto de Fredenslund y col. (1977b).
-203-
Una representacin similar se realiza en la Figura (12-3) para el sistema x1meta-
sin similar de las curvas de (y1-x1) con la de las mezclas anteriores en casi todo el in-
tervalo de concentraciones, exceptuando la zona rica en ster. En ella, los valores expe-
rimentales son algo inferiores a los aportados por Nagata (1962), Gmehling y Onken
0.6
va en la aplicacin del test mencionado
en el parrfo anterior.
-204-
los datos de Beregovykh y col. (1971), ver Gmehling y col. (1982a), a 101,32 kPa y
algo menos con los de Patlashov y col. (1977), ver Gmehling y col. (1988), a 101,32
kPa. Como ya se indic, el test punto a punto de Fredenslund y col. (1977b) no verific
este trabajo de investigacin son de buena calidad para todo el conjunto de los diecisiete
relacin a los datos aportados por la literatura para otras condiciones de trabajo.
del ELV y de acuerdo con el tipo de sistema, para observar as la variacin de las mis-
vs x1 o y1, destaca lo marcado de los puntos azeotrpicos que ofrecen unas curvas ca-
ractersticas para este tipo de mezclas, con distribucin cuasi-lineal en la zona central de
-205-
355 0.5
(a) (b)
y1-x1
T/K
hexano hexano
-0.3
pentano
pentano
310 -0.7
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
0.7 40 800
(c) (d)
hexano
30 600
0.5
pentano
h /(Jmol )
-1
g /RT
i
E
pentano
0.2 hexano
200
10
0 0.9 0
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (12-5a-d). Representacin de las curvas calculadas mediante la ecuacin de las zetas para los
sistemas binarios x1metanol+x2n-alcanos,siendo: (a), T vs x1 o y1; (b), (y1-x1) vs x1; (c), gE/RT o i vs x1;
/d), hE/(Jmol-1) vs x1. En negro valores de la bibliografa, en rojo valores estimados en este trabajo.
-206-
ciones de metanol. Por ltimo, con la modificacin realizada en la ecuacin de las zetas
se destaca nuevamente el hecho por el cual se permite, a partir de una correlacin nica,
reproducir, adems de los datos de ELV la curva de entalpas de exceso para cada mez-
cla, ver Figura (12-5d). La estimacin en estos casos es muy aceptable, con diferencias
fa, ver ltima columna de la pgina 153. Si bien la correlacin no reproduce exacta-
mente las mesetas observadas en la zona central de concentraciones, las curvas obteni-
das conjuntamente con datos de equilibrio isobrico o bien los valores encontrados,
Por ltimo, indicar que la consecuencia del comportamiento de las mezclas me-
destacada auto-asociacin de las molculas del metanol frente a la inercia de los hidro-
carburos saturados y la diferencia del tamao molecular entre las sustancias que consti-
acuerdo con los valores entlpicos y los altos valores de 1 a baja fracciones molares de
metanol, se produce por la rotura de puentes de hidrgeno del alcanol frente a los impe-
Sistemas metanol+steres de alquilo: De la misma manera que para las mezclas meta-
-207-
460 0.8
(a) u=4 (b)
steres butilo
u=4 steres propilo steres butilo
steres etilo u=2 steres propilo
430 steres metilo steres etilo
steres metilo
0.5
u=1
u=2
y1-x1
T/K
u=1
u=1
370
u=1
330 -0.2
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (12-6). Representacin de las curvas calculadas mediante la ecuacin de las zetas para los quince sistemas binarios considerados en este trabajo:
x1CH3OH+x2Cu-1H2u-1COOCvH2v+1, siendo: (a), T vs x1 o y1; y (b), (y1-x1) vs x1
-208-
metanoato etanoato propanoato butanoato
1
metilo
y1
0
1
etilo
y1
0
1
propilo
y1
0
1
butilo
y1
0
0 x1 10 x1 10 x1 10 x1 1
Figura (12-7). Representacin en forma matricial del conjunto de curvas de y1 vs x1, obtenidas para los
quince sistemas formados por: x1metanol+x2alcanoatos (del metanoato al butanoato) de alquilo (desde
el metilo al butilo).
-209-
en dichas figuras se aprecia un orden creciente de las curvas de T y de las diferencias de
(y1-x1) a medida que aumenta la cadena carbonada del alcanoato de alquilo, R1-COO-
R2, o sea, tanto de R1 como de R2. En ambas representaciones se destaca un claro agru-
pamiento de las curvas que corresponden a steres de igual nmero total de tomos de
mao. Otra interesante representacin, en base a las concentraciones de las fases y que
junta para todos los sistemas binarios en estudio. As, la Figura (12-7) constituye una
por las razones apuntadas en la introduccin. Se observa que, en todas las direcciones de
la matriz, existe una variacin uniforme de la magnitud representada con la cadena del
(12-6), para aquellas mezclas donde intervienen alcanoatos con el mismo nmero total
riormente, ver Gonzlez (1996), para otro conjunto de sistemas y se pretende ir am-
res representados, permitir posiblemente mejores logros en futuros estudios, sobre to-
cantidades procedentes casi directamente del equilibrio, como son las concentraciones y
-210-
0.3 4.5 0.3 4.5
(a) (b)
3.5 3.5
0.2 0.2
g /RT
g /RT
i i
E
E
metanoato butilo etanoato butilo
metanoato propilo etanoato propilo
metanoato etilo etanoato etilo
0.1 0.1 etanoato metilo
2 2
1 1
0 0.9 0 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
3.5 3.5
0.2 0.2
g /RT
g /RT
i i
E
1 1
0 0.9 0 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (12-8). Representacin de las curvas calculadas mediante la ecuacin de las zetas para gE/RT o
i vs x1 en los sistemas binarios: (a), x1metanol+x2metanoatos de alquilo; (b), x1metanol+x2etanoatos de
alquilo; (c), x1metanol+x2propanoatos de alquilo; (d), x1metanol+x2butanoatos de alquilo.
-211-
1500 1500
(a) (b)
1100 1100
h /(Jmol )
h /(Jmol )
-1
-1
E
E
400 400
0 0
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
1500 1500
(c) (d)
1100 1100
h /(Jmol )
h /(Jmol )
-1
-1
E
400 400
0 0
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (12-9). Representacin de las curvas calculadas para hE/(Jmol-1) vs x1 en los sistemas binarios: (a)
x1metanol+x2metanoatos de alquilo; (b) x1metanol+x2etanoatos de alquilo; (c) x1metanol+x2propanoatos
de alquilo; (d) x1metanol+x2butanoatos de alquilo. En trazo contnuo se presentan las curvas experimen-
tales y en discontnuo los estimados en este trabajo con la ecuacin de las zetas.
-212-
tamao de la cadena alqulica del ster, R2, sin embargo, no existe una clara variacin
carbonos del ster y, por tanto, la diferencia entre las temperaturas de ebullicin de los
nol, existe la hiptesis de que se forman nuevos complejos por fuertes atracciones di-
polo-dipolo, con efecto exotrmico en la entalpa de mezcla, si bien, la energa neta del
proceso es positiva o endotrmica, ver Figura (12-9ad). Otras tcnicas analticas, como
do que las diferencias de tamao provocan interacciones por las caractersticas estructu-
rales de los compuestos que intervienen, siendo una consecuencia las inflexiones en la
energas de Gibbs de exceso. Posiblemente, para estas mezclas cabra esperar unos re-
sultados que revelaran una mayor discrepancia del comportamiento ideal, sin embargo,
Los resultados de la estimacin de las entalpas de exceso para las quince mez-
valoracin global positiva, si bien ahora, aparecen diferencias medias mayores que para
-213-
el caso de las mezclas con n-alcanos. De cualquier manera, el sistema ms discrepante
con los experimentales es el que contiene metanoato de propilo, con un error medio
inferior al 15%; globalmente, el error medio para el conjunto de quince sistemas es infe-
rior al 6%, resultado que puede considerarse bueno, ver resumen en la ltima columna
0 0 0
-214-
IV. MODELIZACION DE LAS MECLAS
CON METANOL
en el apartado 6 de esta Memoria. El modelo UNIFAC ser aplicado en dos de sus ver-
las parejas OH/CH2 para las mezclas binarias alcanol-alcano y OH/COO para las mez-
zado es el metanol. Los valores de los parmetros se presentan en la Tabla (A-2) del
apndice.
resto de los alcanoles y proponiendo los parmetros CH3OH/CH2 para las mezclas me-
grupo propuestas por Bondi (1968) y normalizados de acuerdo con (5-34); sus valores
ling (1987), considera las mismas interacciones moleculares que la versin original de
-217-
UNIFAC, para las mezclas de metanol con alcanos y con steres de alquilo, incluyendo
la diferenciacin de los metanoatos, con su grupo especfico ya indicado antes. Las Ta-
para los distintos grupos y los valores correspondientes a las interacciones consideradas
fase vapor, yi, que corresponde a cada pareja de (T,xi). En este clculo, el modelo esti-
ma, de acuerdo con los valores de la interaccin, la presencia o no del punto azeotrpi-
co, xaz. La comparacin con los resultados experimentales permitir verificar la aplica-
muestra en la Tabla (13-1) donde se dan los errores medios obtenidos en la prediccin
as como la localizacin del punto azeotrpico de los dos sistemas estudiados de meta-
nol+alcanos.
de la prediccin realizada con los modelos ASOG y UNIFAC-modificado que son los
-218-
Tabla (13-1). Resultados de la estimacin de i, y1 y xaz para las mezclas binarias x1metanol+x2n-alcanos.
modelos que mejor estimacin produce en los coeficientes de actividad y en las con-
bos sistemas binarios, los modelos ASOG y UNIFAC-modificado ofrecen mejores pre-
dicciones que el modelo UNIFAC-original en todas las cantidades, i, gE/RT e yi, consi-
Igual que en el caso de las mezclas con alcanos, en la Tabla (13-2) se muestran
res de alquilo, con respecto a los datos experimentales. En ella se recoge tambin la es-
-219-
Tabla (13-2). Resultados de la estimacin de i,, y1 y xaz para las mezclas x1metanol+x2steres de alquilo.
timacin de los modelos sobre los puntos azeotrpicos de cada uno de los sistemas.
la funcin de Gibbs y los coeficientes de actividad con los modelos que mejor reprodu-
cian los datos experimentales, o se los que ofrecen una mejor prediccin de los coefi-
cientes de actividad y las concentraciones de la fase vapor, que son los de ASOG y
Modelo ASOG. La prediccin realizada por este modelo es aceptable en los sistemas
binarios formados por los steres de metilo y etilo, con un error medio en la estimacin
respectivamente. Para los otros dos conjuntos de steres, de propilo y de butilo, la esti-
-220-
0.7 40 0.7 40
(a) (b)
30 30
0.5 0.5
g /RT
g /RT
i i
E
E
0.2 0.2
10 10
0 0.9 0 0.9
0 0.2 0.4 x1 0.6 0.8 1 0 0.2 0.4 x1 0.6 0.8 1
Figura (13-1). Comparacin de las curvas de gE/RT y i obtenidas mediante la ecuacin de las ze-
tas() y los modelos de contribucin de grupos ASOG () y UNIFAC-modificado () para los
sistemas: (a) x1metanol+x2n-pentano y (b) x1metanol+x2n-hextano.
ster y las diferencias en las curvas de gE/RT son cada vez mayores. La estimacin de la
composicin del azetropo se realiza con una diferencia inferior a 0,017 unidades en
todos los casos, sin embargo, el mtodo estima la existencia de un punto singular que
perimentales al aumentar el nmero de carbonos del ster, con errores medios que van
-221-
metanoato etanoato propanoato butanoato
5
0.3
0.2
metilo
i
g /RT
E
0.1
0 0.9
0.3 5
0.2
etilo
i
g /RT
E
0.1
0 0.9
0.3 5
propilo
0.2
g /RT
i
E
0.1
0 0.9
0.35 5
4
0.25
butilo
g /RT
i
E
0.1 2
0 0.9
0 x1 1 0 x1 1 0 x1 1 0 x1 1
Figura (13-2). Representacin en forma matricial del conjunto de curvas gE/RT y 1 vs x1 obtenidas me-
diante la ecuacin de las zetas () y los modelos de contribucin de grupos ASOG () y UNI-
FAC-modificado () para los quince sistemas formados por: x1metanol+x2alcanoatos (del metanoa-
to al butanoato) de alquilo (desde el metilo al butilo).
-222-
desde el 2,1% para el sistema metanol+etanaoto de metilo al 19,4% para metanol+buta-
presentando un error medio global del 5%. Para los steres de butilo las diferencias me-
dias fueron inferior al 2% y superior al 10% para los butanoatos de metilo y etilo. En
nal de UNIFAC con una cierta disminucin en el error medio tanto de i como de yi para
la mayora de las mezclas binarias estudiadas. El error medio global en el clculo de las
concentraciones de la fase vapor, y1, fue del 4,2% y del 8,5% en los coeficientes de acti-
vidad. Aunque los sistemas que contienen butanoato de metilo y propanoato de etilo no
presentan azetropo experimental, esta versin del modelo lo estima. Para los dems
sistemas, la situacin del azetropo se predice con una desviacin inferior a 0,035 uni-
es de 0,06 unidades.
0 0 0
-223-
V. CONCLUSIONES
Y BIBLIOGRAFIA
14. CONCLUSIONES
15. BIBLIOGRAFIA
14. CONCLUSIONES
Uno los objetivos que se plantearon para el desarrollo de este trabajo, como parte
de un amplio proyecto que viene realizndose en nuestro grupo de investigacin
sobre las mezclas steres de alquilo+alcanoles, fue la de completar los estudios
con el metanol, el cual nunca haba sido empleado en nuestras investigaciones,
dndole con ello un carcter especial al distinguirlo en su comportamiento de los
dems alcanoles. Por ello, el primer paso fue la de recopilar y anlizar lo publi-
cado en la bibliografa hasta el momento, con el fin de convenir o no la redeter-
minacin de nuevos datos o bien crear nuevas condiciones de trabajo. Esta reco-
pilacin, cuya importante labor deseamos destacar, merece la primera de las
conclusiones, ya que permiti elegir adecuadamente la presin de trabajo en
141,3 kPa y las mezclas o sistemas para el desarrollo del mismo.
-227-
Con un ebullmetro de pequea capacidad, chequeado previamente con la mez-
cla metanol+agua a las presiones de 101,32 y 141,3 kPa, se determinaron expe-
rimentalmente los datos de equilibrio lquido-vapor a la presin de (141,30,2)
kPa de dos sistemas binarios formados por el metanol y dos n-alcanos (C5,C6), y
del metanol con quince steres de alquilo, desde el metanoato al butanoato, con
excepcin del metanoato de metilo. La inclusin de los valores obtenidos en los
de la literatura existente, valid el sistema y la tcnica experimental utilizada.
-228-
Se calcularon los valores de la energa de Gibbs en su forma adimensional y los
coeficientes de actividad, que se correlacionaron simultneamente utilizando las
ecuaciones derivadas de los modelos de composicin local, Wilson, NRTL y
UNIQUAC y una modificacin realizada sobre la ecuacin utilizada por nuestro
grupo, las zetas, en las que se incluy una variacin con la temperatura para los
coeficientes de su forma polinmica. El anlisis conjunto realizado para todas
las mezclas permite valorar positivamente la utilizacin de la ecuacin de las
zetas, ya que ofrece mejores o similares correlaciones que los otros tres modelos.
Una de las conclusiones ms positivas es la adecuada reproduccin que se hace
con dicha ecuacin modificada, de las entalpas de exceso de cada mezcla, con
errores que, salvo excepciones, no superaban el 10%. Futuras modificaciones se
preveen para esta ecuacin, cuya extensin permitir conseguir una ptima co-
rrelacin para el tratamiento de datos de ELV y otras magnitudes.
-229-
ms acusada en el mtodo UNIFAC, donde se llegan a establecer diferencias en
los coeficientes de actividad cercanas al 20%. Con este mtodo, a pesar de tener
un grupo especfico para la interaccin de los metanoatos, en la estimacin de
los i las diferencias se elevaron hasta el 12% en la mezcla metanol+metanoato
de butilo. La estimacin de los puntos azeotrpicos con los distintos modelos fue
un tanto irregular y los valores ms cercanos fueron los de ASOG, si bien se
predice la presencia de puntos singulares no existentes en la experimentacin,
para tres de las mezclas estudiadas. Para las mezclas con n-alcanos, la estima-
cin de los coeficientes de actividad y de las concentraciones de la fase vapor
con todos los modelos puede considerarse buena.
Por ltimo, la realizacin del trabajo permitir dirigir mejor las futuras investi-
gaciones de nuestro grupo de Termodinmica. En este sentido, sera conveniente
en primer lugar, la modificacin de algunos aspectos de la tcnica experimental,
bsicamente el propio ebullmetro de equilibrio. Otro objetivo importante es la
finalizacin de una matriz de sistemas binarios sobre las mezclas de alcano-
les+steres de alquilo, con el fin de su correcta y amplia utilizacin de las teo-
ras de contribucin de grupos. En este tema concreto puede decirse que, acep-
tando las estimaciones del modelo ASOG, se recomendara para futuros trabajos
conseguida una mayor base de datos, plantear la redeterminacin de nuevos va-
lores de las dos versiones del modelo UNIFAC con las interacciones correspon-
dientes al grupo carboxilato; aconsejando la simplificacin de dicho grupo y
planteando la interaccin alcanol/ester como OH/COO, ya que, si bien el grupo
funcional es el mismo, la interaccin difiere si se trata de un ester saturado o in-
saturado.
-230-
15. BIBLIOGRAFA
- ALESSI, P.; FERMEGLIA, M. y KIKIC, I. (1989). J. Chem. Eng. Data 34: 236.
- ALFONSO de, C.; EZAMA, R.; LOSADA, P.; CALAMA, M.A.; LLANAS, B.;
PINTADO, M. y SAENZ DE LA TORRE, A.F. (1983). An. Qum. Esp., 79:
243.
- ARCE, A.; BLANCO, A.; PEREZ, J.C. y SOTO, A. (1994). J. Chem. Eng. Data, 39:
95.
-231-
- DERR, E.L. y DEAL. C.H. (1969). I. Chem. E. Symposium Series. No. 32, 3: 40.
- GALVAN, S.; ORTEGA, J.; SUSIAL, P. y PEA, J.A. (1994). J. Chem. Eng. Jpn.,
27: 529.
- GMEHLING, J.; LI, J. y SCHILLER, M. (1993). Ind. Eng. Chem. Res., 32: 178.
-232-
I/5.
- GROLIER, J.-P.E. (1977). Int. DATA Ser., Ser. B., Thermodynamic Properties of
Aqueous Organic Systems, p.4.
- HIAKI, T.; TOCHIGI, K. y KOJIMA, K. (1986). Fluid Phase Equilibria, 26: 83.
- HOLMES, M.J. y VAN WINKLE, P. (1970). Ind. Eng. Chem., 62: 21.
- HONGO, M.; TSUJI, T.; FUKUCHI, K. y ARAI, Y. (1994). J. Chem. Eng. Data,
39: 688.
- HRADETZKY, G. y BITTRICH, H.-J. (1986). Int. DATA Ser., Sel. Data Mixtures,
Ser. A, (3): 216.
-233-
ASOG Method. Kodansha Ltd., Tokyo.
- KOJIMA, K.; MOON, H.M. y OCHI, K. (1990). Fluid Phase Equilibria, 56: 269.
- MARGULES, M. (1895). Sitzber Akad. Wiss. Wien. Math. Naturen. Kl., II, 104:
1243.
- NAGATA, I.; OHTA, T.; OGURA, M. y YASUDA, S. (1976). J. Chem. Eng. Data,
21: 310.
-234-
- ORACZ, P. y WARYCHA, S. (1995). Fluid Phase Equilibria, 108: 199.
- ORTEGA , J. (1995). Int. DATA Ser., Sel Data Mixtures, Ser. A. 23(3): 154.
- PATEL, N.C. y SANDLER, S.I. (1985). J. Chem. Eng. Data, 30: 218.
- PEA, J.A., (1985). Tesis Doctoral, Escuela Superior de Ingenieros Industriales. Univ.
Politcnica de Las Palmas.
- PITZER, K.S. y CURL, R.F. (1957). J. Am. Chem. Soc. 79: 2369.
- RAAL, J. D.; CODE, R. K. y BEST, D. A. (1972). J. Chem. Eng. Data, 17: 211.
-235-
Gases and Liquids, 2nd. ed; Mac Graw-Hill, New York.
- SPENCER, C.F. y DANNER, R.P. (1972) J. Chem. Eng. Data, 17: 236.
-236-
- VAN LAAR, J.J. (1910). Z. Phys. Chem., 72: 723.
- VAN NESS, H.C.; BYERS, S.M. y GIBBS, R.E. (1973). AIChE. Journal, 19: 238.
- WILSON, G.M. y DEAL, C.H. (1962). Ind. Eng. Chem. Fundamentals., 1: 20.
0 0 0
-237-
VI. APENDICES
16. APENDICES
Tabla (A-1). Valores de propiedades y constantes caractersticas de los compuestos puros utilizados en este trabajo: pc, presin crtica; Tc, temperatura
crtica; vc, volumen crtico; ZRA, constante de la ecuacin (3-31); , factor acntrico; a y b , constantes de la ecuacin de Tsonopoulos
.
pc /MPa (1) Tc /K (1) 109 vc (1) exp
(2) (1) (3) (3)
Producto ZRA a b
/m3mol-1
Metanol 8,0959 512,58 117,8 0,2318 0,566 0,561 0,08780 0,0560
n-Pentano 3,3688 469,65 312,3 0,2685 0,249 0,247 0 0
n-Hexano 3,0123 507,43 370,1 0,2635 0,305 0,298 0 0
Etanoato de Metilo 4,6914 506,80 228,0 0,2550 0,325 0,313 -0,01120 0
Propanoato de Metilo 4,0044 530,60 282,0 0,2556 0,353 0,331 -0,00868 0
Butanoato de Metilo 3,4734 554,50 340,0 0,2562 0,381 0,367 -0,00689 0
Metanoato de Etilo 4,7420 508,40 229,0 0,2586 0,285 0,273 -0,01440 0
Etanoato de Etilo 3,8301 523,25 286,0 0,2538 0,361 0,358 -0,00936 0
Propanoato de Etilo 3,3620 546,00 345,0 0,2545 0,394 0,391 -0,00738 0
Butanoato de Etilo 3,0600 571,00 421,0 0,2570 0,419 0,412 -0,00614 0
Metanoato de Propilo 4,0631 538,00 285,0 0,2589 0,318 0,298 -0,01080 0
Etanoato de Propilo 3,3640 549,40 345,0 0,2524 0,394 0,392 -0,00750 0
Propanoato de Propilo 3,1100 578,00 389,0 0,2584 0,376 0,377 -0,00630 0
Butanoato de Propilo 2,8000 594,00 442,0 0,2581 0,448 0,450 -0,00510 0
Metanoato de Butilo 3,5099 559,00 336,0 0,2563 0,384 0,382 -0,00873 0
Etanoato de Butilo 3,1604 579,15 389,0 0,2579 0,417 0,420 -0,00733 0
Propanoato de Butilo 2,7996 594,00 442,0 0,2574 0,475 0,469 -0,00542 0
Butanoato de Butilo 2,5402 616,00 494,0 0,2555 0,485 0,481 -0,00458 0
(1) (2) (3)
Daubert y Danner (1984); Spencer y Danner (1972); Tsonopoulos (1974)
-241-
Tabla (A-2). Parmetros de interaccin mkl y nkl para el modelo ASOG, Tochigi y col., (1990).
Tabla (A-3). Parmetros de volumen y superficie para UNIFAC-original, Hansen y col., (1991)
Tabla (A-4). Parmetros de interaccin amn del modelo UNIFAC-original, Hansen y col., (1991).
-242-
Tabla (A-5). Parmetros de volumen y superficie para UNIFAC-modificado, Gmehling y col., (1993)
Tabla (A-6). Parmetros de interaccin para amn,i para UNIFAC-modificado, Gmehling y col. (1993).
-243-
A-7. Conceptos y definiciones matemticas. Simbologa y unidades
Magnitudes termodinmicas:
H
C p = Capacidad calorfica a presin constante, [JK-1].
T p
M
M i = Propiedad molar parcial, [variable].
ni T , p ,n j i
G
i = Potencial qumico, [Jmol-1].
i T , p ,n j i
n
Parmetros estadsticos:
(Y Ycal )
2 Desviacin estndar. N, nmero de datos experimenta-
s( Y ) = les; n, nmero de coeficientes de la correlacin, [s][Y]
exp
N n
[
s( Yi ) = s( Y1 )2 + s( Y2 )2 ]
1/ 2
Desviacin estndar promedio de Y, [s] [Y].
1 N Yi ,exp Yi ,cal
e ( Yi ) = Y
N i =1
.100 Error relativo medio porcentual, [-].
i ,exp
1 2
e ( Yi ) = e ( Y j )
2 j =1
Valor medio del error de Y, [ e ][Y].
1 N
( Y )= Yi ,exp Yi ,cal
N i =1
Desviacin media, [ ] [Y].
-244-
-245-