Instituto Politécnico Nacional.

Escuela Superior de Ingeniería e Ciencias Extractivas.

Departamento de Ingeniería Química Petrolera

Laboratorio De Termodinámica Del Equilibrio De Fases.

Profesor: Dr. Octavio Elizalde Solís

Auxiliar: Ing. María Mateos Hernández

Grupo: 2PM31
Equipo: 6

Nombre De Práctica: Equilibrio liquido-vapor de una mezcla binaria. Ley de Raoult

modificada
Numero de práctica: 5
Miembros Del Equipo:

.- González Barragán Jared…………………

.- Melo Lugo Denisse Guadalupe…………….

.- Soto Orihuela Karla Michelle……………………

.- Toral Flores Raúl……………………………….

.- Zavala Díaz Fernando……………………….

Fecha experimental: 27/05/2021

1
 INDICE

 Objetivos……………………………………………………….2
 Introducción……………………………………………………2
 Sección experimental………………………………………...8
 Sección de resultados………………………………………..12
 Observaciones y recomendaciones…………………………34
 Conclusiones…………………………………………………..34
 Bibliografías……………………………………………………37
OBJETIVOS

 Determinar experimentalmente el equilibrio liquido-vapor para una mezcla

binaria a presión atmosférica.
 Preparar soluciones de la mezcla binaria 2-propanol + agua a diferentes
composiciones.
 Obtener la temperatura de burbuja experimental para cada composición.
 Determinar el coeficiente de actividad de cada componente en la solución
utilizando el modelo de Margules o van Laar
 Determinar la temperatura y composición a la cual se forma el azeótropo en la
mezcla.
 Comparar los resultados experimentales con los teóricos a través de un cálculo
de punto de burbuja.
INTRODUCCIÓN
POTENCIAL QUÍMICO
La afinidad química, también conocida como energía molar parcial, es un tipo de
energía que puede ser absorbida o liberada durante
un efecto químico. También puede cambiar con el tiempo. El potencial
de las especies químicas en una mezcla se puede definir como una disminución de
la energía sistémica en respuesta a cambios en el número de moles de cada una
de las especies que tenemos. [1]
FUGACIDAD
Es una propiedad que comúnmente se usa como una unidad de soporte
termodinámico que actúa como un amplificador para convertir la corriente de
calibración al modo de operación real. En el caso de los líquidos, debe entenderse
en función de la temperatura del agua, teniendo en cuenta la presión de vapor. Cada
propiedad esencial debe medirse por el tamaño de las partículas en el vapor de
cada material. Después de esta carga, debe utilizar la propiedad de estado para
evaluar la correlación media. [2]

2
 COEFICIENTE DE FUGACIDAD
El coeficiente de fugacidad es la relación entre fugacidad y presión. Dependiendo
del gas, la presión del gas será igual a la presión del sistema. Para los líquidos, es
necesario considerar la presión de vapor, que se entiende en función de la
temperatura del agua. Ésta es la relación que existe entre el voltaje y la presión del
componente como gas positivo. Esta sustancia se considera positiva y su valor se
encuentra entre 0 y 1. [2]
EQUILIBRIO MULTIFÁSICO
El equilibrio multifásico como su nombre lo dice es un equilibrio termodinámico
donde existen un gran número de fases. En termodinámica se conoce el cálculo
más importante en equilibrio de fases, Esto incluye el número exacto y el tipo de
pasos y los elementos relevantes al equilibrio, de esa manera ahora podemos decir
que la energía de Gibbs en un sistema sea mínimo. [3]
EQUILIBRIO MULTICOMPONENTE
Una balanza de componentes múltiples es un conjunto de balanzas de varios
componentes. En las redes de procesos químicos, los componentes son ciertas
especies químicas presentes en las corrientes como componentes de una mezcla.
Generalmente, las especies se pueden transformar unas en otras mediante
reacciones químicas; por lo que sus cantidades no pueden conservarse
individualmente. [4]
EQUILIBRIO TERMODINÁMICO
Cuando la fuerza resultante que actúa sobre la estabilización se vuelve cero, el
proceso está en estabilización mecánica. Si todas las partes u objetos que lo
producen están a la misma temperatura, están en el rango de temperatura, se intuye
una estabilización térmica, y si no ocurre ninguna reacción química, están en
equilibrio químico. Por lo tanto “Si un sistema que está en equilibrio mecánico,
térmico y químico se dice que está en equilibrio termodinámico”. [5]
EQUILIBRIO TÉRMICO
Se ha observado que cuando un objeto caliente se encuentra con un objeto frío, el
calor se convierte en algo a bajas temperaturas. Los dos elementos alcanzan la
misma temperatura y ambos mantienen la temperatura ambiente sin perder el otro.
Entonces se dice que están en el termómetro. El control de la temperatura es objeto
de la ley cero de la termodinámica. [5]
EQUILIBRIO MECÁNICO
El equilibrio es una condición de desarrollo que satisface una de dos condiciones:
Cuando la suma de potencia y tiempo para cada parte del sistema es cero, el
sistema se convierte en la suma de los equipos de ingeniería estándar. Si la posición

3
del sistema en la cámara fría es donde la gravedad superficial es cero, el sistema
está mecánicamente en equilibrio. [5]
EQUILIBRIO QUÍMICO
Asegúrese de que no haya reacción química ni contaminante. El contenido siempre
responde, pero el estado del sistema de respuesta no cambia con el tiempo. [6]
EQUILIBRIO LIQUIDO-VAPOR
Esto se logra cuando el efecto químico es el mismo en ambos niveles en un
ambiente y presión particulares. El ajuste ocurre cuando se alcanza la presión de
vapor (la presión aplicada cuando el estado de vapor está sincronizado con el nivel
del agua). Además de depender de la temperatura, depende de las propiedades de
los componentes y del diseño de cada componente de las fases que se tenga. [7]
PUNTO DE SATURACIÓN
Desde un punto de vista químico, el significado de los nervios es que el sistema
digestivo no puede disolver un problema a ninguna temperatura. Hay soluciones
que requieren una gran cantidad de soluciones para lograr este problema y otras
que lo hacen en pequeñas cantidades. Todas las sustancias pueden transportar una
variedad de líquidos, gases, proteínas y energía. Una vez que se alcanza el punto
de saturación, no se puede absorber de manera adicional más cosas. [8]
PUNTO DE ROCIÓ
El punto de rocío se conoce como la temperatura y la presión a las que la mezcla
de vapor comienza a condensarse. También significa la presión más baja durante
el turno de tiempo. [7]
PUNTO DE BURBUJA
Esto significa la temperatura y presión a las que se forman por primera vez las
partículas de vapor en la mezcla de agua, es decir, comienzan a hervir. Además, la
clave para que pueda existir un cambio es el tiempo ejerciendo una mayor presión.
[7]

4
 DIAGRAMAS DE FASES PARA MEZCLAS BINARIAS

[9]
LEY DE RAOULT PARA SOLUCIONES NO IDEALES
Para soluciones no lineales, la ley de Raoult se puede ajustar al incluir dos
elementos que describen las interacciones de diferentes moléculas.
El primero es compensar la escasez de líneas de gas o desviarse de la legislación
de gas correcta. Esto se llama transitorio y se explica de la siguiente manera:
𝜑𝑝,𝑖

El segundo factor, comúnmente conocido como coeficiente de actividad, es la

corrección de la interacción de diferentes moléculas en el medio líquido. Esta
situación se describe a continuación.
𝛾𝑖
Para que finalmente se coloque de la siguiente manera: 𝛾𝑖 𝜑𝑝,𝑖 𝑃𝑖 (𝑇) = 𝑥𝑖𝐿 𝛾𝑖 𝑃𝑖0 (𝑇) [10]

MODELOS DE SOLUCION
A que se refiere esta terminología, su definición como tal son por así decirlo bases de
cálculo de la cuales nos guiaremos para resolver por cierto procedimiento especifico,
ahora bien, debemos resaltar médelos más específicos como son los modelos
termodinámicos son desenlaces del tipo matemáticos que nos ayudan a analizar el
comportamiento que en este caso es termodinámico de una sustancia pura o también
de una mezcla donde estén presentes propiedades como la presión, la polaridad, la

5
temperatura, etc., a través de una metodología más precisas hechas por personas o
maquinas como se hace hoy en día. [11]
AZEOTROPO
Esto se conoce como mezcla azeotrópica o mezcla azeotrópica con una mezcla de
compuestos (dos o más sustancias) que están en estado líquido y actúan como si
solo se reunieran. La naturaleza de esta mezcla azeotrópica (mezcla azeotrópica)
le permite utilizar la mezcla como un ingrediente o componente puro (fácil de
procesar y almacenar).
Existen dos tipos de azeótropos los que son positivos y los que son negativos, en
donde los presentan desviaciones positivas y negativas respectivamente, todo con
la ley de Raoult. [12]
IMPORTANCIA Y APLICACIÓN EN LA INDUSTRIA DEL EQUILIBRIO LÍQUIDO-
VAPOR
 Procesamiento de productos naturales o petróleo.
 Recuperación de bifenilo
 Extracción de líquidos (NGLs) del gas natural
 Producción supercrítica de biodiesel
Al aplicar sistemáticamente este sistema para resolver problemas de diseño y
mejorar los métodos de tratamiento químico, puede identificar posibles áreas
funcionales que conducirán a proyectos exitosos. [13]
PRINCIPIO DE OPERACIÓN DE UN REFRACTÓMETRO
Un refractómetro, como su nombre indica, es un instrumento visual que se centra
en el estudio de la distorsión de la luz. Se utiliza un refractómetro para medir el
índice de refracción de neumáticos hechos de material translúcido congelado. El
refractómetro funciona según el principio de respuesta, desviando la luz a medida
que se mueve de un objeto a otro. Un ejemplo de esto ocurre cuando miras un objeto
que se hunde y tratas de alcanzarlo. Esto se debe a la distorsión parcial o parcial
de las ondas de luz que se mueven del agua al aire y viceversa. [14]
EQUILIBRIO LIQUIDO-VAPOR EN SISTEMAS NO IDEALES
La solución es que una cierta presión de vapor no es proporcional al peso molecular
de la porción al nivel del agua de todos los microorganismos. Tenga en cuenta que
no hay diferencia entre sólido y soluble. Las mezclas no ideales presentan
desviaciones en ambas fases, liquidas y gaseosa, dentro de las cuales las más
comunes son los azeótropos [15]

6
 MODELOS SOLUCION
Flory-Huggins
La teoría de Flory-Huggins es un modelo matemático de las propiedades térmicas
de las soluciones de polímeros, teniendo en cuenta los valores de gran tamaño
molecular y ajustando la representación típica para condiciones mixtas. El resultado
es la conversión de energía libre de Gibbs al combinar el polímero con un disolvente.
También incluye la simplificación de supuestos, pero con resultados que ayudan a
interpretar ciertos experimentos. [16]
UNIQUAC
El método UNIQUAC se utiliza para calcular las calorías libres y el exceso de gas y
agua. Se utiliza en sistemas no objetivos con muy buenos resultados. Este modelo
permite el uso de varios parámetros avanzados (objetivo, virial, etc.), pero el original
se usa para un funcionamiento estándar basado en el cálculo de temperatura de
Redlich-Kwong. El modelo de producto es ineficaz. Cambio de temperatura en
comparación con otros tipos de modelaciones por mencionar a Wilson. [16]
UNIFAC
Este es un sistema intermedio para detectar la actividad no conductora de mezclas
no deseadas. UNIFAC utiliza un grupo de trabajo de moléculas que combinan
mezclas líquidas para cuantificar el valor de un proyecto. La interacción de cada
grupo funcional en la celda, así como otras interacciones binarias importantes, se
pueden utilizar para calcular la función de cada pieza de información. Esta
información se puede utilizar para obtener información sobre la estabilidad del agua
que es útil para muchos cálculos de temperatura, como la composición química y
las propiedades de fusión. [16]
Van Laar
La ecuación de Van Laar es un modelo de actividad termodinámica, para explicar
los equilibrios de etapa de mezclas líquidas. Gracias a esto, el modelo perdió la
conexión con las características moleculares y, por consiguiente, debería considerar
como un modelo experimental para corregir los resultados experimentales. [16]
NRTL
El término de NRTL se fundamenta en la conjetura de Wilson de que la
concentración local cerca de una molécula es distinta de la concentración de masa.
Esta diferencia se debería a una diferencia entre la energía de relación de la
molécula central con las moléculas de su propio tipo y las moléculas del otro tipo.
La diferencia de energía además introduce una no aleatoriedad a grado molecular
local. [16]

7
 Wilson
Este modelo se creó por Wilson, posibilita calcular los coeficientes de actividad en
etapa líquida para sistemas enormemente no ideales, en especial sistemas de
alcoholes y agua, continuamente en una sola etapa líquida. La aplicación industrial
más común se encuentra en la estandarización de estándares para sistemas
eficientes no invasivos basados en vapor de agua. Esto es especialmente útil para
operaciones destructivas y, en general, interacciones vapor-agua. [16]
CALCULO DE COEFICIENTES DE FUGACIDAD Y ACTIVIDAD
Fugacidad
Primeramente, para realizar los cálculos de los coeficientes de fugacidad, se puede
realizar por medio de propiedades específicas y su función de estado, haciéndolo
por separado, como en la fase liquida y también para la fase vapor, primeramente,
para la fase vapor en necesario conocer su entalpia, su entropía y su presión, de
esa manera hacemos la relación de fugacidad y la presión, pero del sistema; para
la fase liquida, conociendo la presión o el volumen de saturación. Se hace la relación
con su fugacidad y la presión con respecto del vapor, finalmente se hace una
evaluación de fugacidades por fase y la presión del sistema, que en este caso que
deberá ser diferente de uno.
Actividad
Si tenemos en consideración que la actividad para un sistema no ideal podemos
decir que tiene como proveniente la fracción molar, y debido a que existe una
relación entre la mezcla de gases ideales con su respectiva solución verdadera [17]
Sección Experimental
Materiales:
 Termómetro.
 Embudo.
 Buretas de 50 ml.
 Matraz de bola.
 Vasos de precipitados.
 Frascos.
 Etiquetas adhesivas.
 Servilletas absorbentes de papel.
 Refractómetro.
Reactivos:
 etanol 99.9% pureza
 n-heptano 99.9% pureza

8
Diagrama del equipo experimental

9
 Diagrama de bloques

Preparar 14 frascos con la Agregar los volúmenes de

mezlca binaria de etanol + Agitar cada solución
cada reactivo para cada
n-heptano, para tener solución por medio de buretas dentro de los frascos.
diferentes composiciones en frascos previamente
(fase líquida). etiquetados para su
identificación.

Para cada solución

Para cada solución,
Colocar el termómetro agregar unas gotas en el
verter en el matraz del
en el equipo y calentar refractómetro y obtener
equipo de equilibrio
hasta ebullición. el índice de refracción de
liquido-vapor con la
la fase liquida 𝑛𝑙𝑖𝑞.
ayuda de un embudo.

Registrar la temperatura Colectar muestras de la Analizar cada muestra del

de burbuja 𝑇𝐵 de cada fase vapor condensada destilado con el
solución a presión en frascos y taparlos refractómetro (fase vapor)
atmosférica. inmediatamente. 𝑛𝑣𝑎𝑝.

Repetir los pasos del 4 al 9

para cada solución (frasco).

10
 Tabla De Datos Experimentales

Frasco V1(ml)Et V2(ml)hept nliq T°C nvap

1 0 40 1.3875 81 1.3875
2 0.4 39.6 1.3875 74 1.3813
3 2 38 1.387 65 1.3774
4 4 36 1.3835 63 1.377
5 7 33 1.382 61.5 1.376
6 10 30 1.379 61 1.3755
7 15 25 1.3745 61 1.3755
8 20 20 1.3715 61 1.3755
9 25 15 1.3675 61 1.375
10 30 10 1.3645 61 1.375
11 35 5 1.3615 64 1.374
12 38 2 1.36 68 1.369
13 39 1 1.3595 70 1.3662
14 40 0 1.358 71 1.358

11
 Sección De Resultados

Tabla De Datos Bibliográficos [19]

Datos Bibliográficos
T°C x1 y1
90.5 0.013 0.205
85 0.023 0.33
76.8 0.051 0.49
75.6 0.083 0.535
73.3 0.181 0.58
72.4 0.241 0.59
71.6 0.309 0.605
71.4 0.406 0.625
71.1 0.546 0.635
71.4 0.61 0.648
71.2 0.66 0.653
71.3 0.679 0.66
71.3 0.69 0.666
71.2 0.742 0.678
71.5 0.821 0.7
71.9 0.88 0.732
72.1 0.886 0.732
72.5 0.91 0.76
74.5 0.961 0.847
74.3 0.961 0.837
75.4 0.969 0.881

1.- Calcular la composición de la fase liquida (𝑥𝑖 𝑒𝑥𝑝 ) a partir de los volúmenes
agregados en las buretas.
𝑚 𝑛
𝑚𝑖 = 𝑝𝑖 ∗ 𝑣𝑖 𝑛𝑖 = 𝑃𝑀𝑖 𝑥𝑖 𝑒𝑥𝑝 = 𝑛 𝑖
𝑖 𝑇

Para el componente 1
 Evento 1
𝑚𝑖 = 0.789 𝑔 ∗ 0𝑐𝑚3 = 0𝑔
𝑐𝑚3
 Evento 2
𝑚𝑖 = 0.789 𝑔 ∗ 0.4𝑐𝑚3 = 0.3156𝑔
𝑐𝑚3

12
  Evento 3
𝑚𝑖 = 0.789 𝑔 ∗ 2𝑐𝑚3 = 1.586𝑔
𝑐𝑚3

 Evento 1
0𝑔
𝑛𝑖 = = 𝑜𝑚𝑜𝑙
46.07 𝑔
𝑚𝑜𝑙

 Evento 2
0.3156 𝑔
𝑛𝑖 = = 0.006850445𝑚𝑜𝑙
46.07 𝑔
𝑚𝑜𝑙
 Evento 3
1.586𝑔
𝑛𝑖 = 46.97 = 0.034252225 𝑚𝑜𝑙
𝑔
𝑚𝑜𝑙
 Evento 1
0𝑚𝑜𝑙
𝑥𝑖 𝑒𝑥𝑝 = =0
0.273026644 𝑚𝑜𝑙
 Evento 2
0.006850445𝑚𝑜𝑙
𝑥𝑖 𝑒𝑥𝑝 = = 0.024717747
0.277146823 𝑚𝑜𝑙

 Evento 1

0.034252225 𝑚𝑜𝑙
𝑥𝑖 𝑒𝑥𝑝 = = 0.11665195
0.293627537 𝑚𝑜𝑙

Para el componente 2
 Evento 1
𝑚𝑖 = 0.684 𝑔 ∗ 40𝑐𝑚3 = 27.36𝑔
𝑐𝑚3
 Evento 2
𝑚𝑖 = 0.684 𝑔 ∗ 39.6𝑐𝑚3 = 27.0864𝑔
𝑐𝑚3
 Evento 3
𝑚𝑖 = 0.684 𝑔 ∗ 38𝑐𝑚3 = 25.992𝑔
𝑐𝑚3

 Evento 1

13
 27.36 𝑔
𝑛𝑖 = = 0.273026644𝑚𝑜𝑙
100.21 𝑔
𝑚𝑜𝑙

 Evento 2
27.0864 𝑔
𝑛𝑖 = = 0.270296378𝑚𝑜𝑙
100.21 𝑔
𝑚𝑜𝑙
 Evento 3
25.992𝑔
𝑛𝑖 = 100.21 = 0.259375312 𝑚𝑜𝑙
𝑔
𝑚𝑜𝑙
 Evento 1
0.27302664𝑚𝑜𝑙
𝑥𝑖 𝑒𝑥𝑝 = =1
0.273026644 𝑚𝑜𝑙
 Evento 2
0.270296378𝑚𝑜𝑙
𝑥𝑖 𝑒𝑥𝑝 = = 0.975282253
0.277146823 𝑚𝑜𝑙

 Evento 3

0.259375312 𝑚𝑜𝑙
𝑥𝑖 𝑒𝑥𝑝 = = 0.88334805
0.293627537 𝑚𝑜𝑙

2. - Graficar la relación composición de la fase liquida (𝑥𝑖 𝑒𝑥𝑝 ), eje y, - índice

de refracción de la fase liquida 𝑛𝑙𝑖𝑞 , eje x.

1.2
xi exp VS nliq y = -631.54x2 + 1701.8x - 1145.4

0.8
Xi exp

0.6

0.4

0.2

0
1.355 1.36 1.365 1.37 1.375 1.38 1.385 1.39
Indice De Refraccion

14
3.- Obtener una ecuación para calcular la composición del componente uno en
función del índice de refracción, no considere el frasco inicial y final.

𝑦 = −631.54 𝑥 2 + 1701.8 𝑥 − 1145.4

4.- A partir de la ecuación anterior, calcular la composición del componente uno
𝑦1 exp utilizando el índice de refracción de la fase vapor 𝑛𝑣𝑎𝑝 . La composición del
componente uno en los frascos inicial y final para la fase vapor, es la misma que
los correspondientes en fase liquida.
Componente 1
 Evento 1
𝑦 = (−631.54 (0)2 + 1701.8 (0)) − 1145.4 = 0

 Evento 2
𝑦 = (−631.54 (1.3813 )2 + 1701.8 (1.3813)) − 1145.4 = 0.324531177

 Evento 3
𝑦 = (−631.54 (1.3774)2 + 1701.8 (1.3774)) − 1145.4 = 0.48220583

Componente 2
 Evento 1
𝑦 =1−0= 1
 Evento 2
𝑦 = 1 − 0.324531177 = 0.675468823

 Evento 3
𝑦 = 1 − 0.48220583 = 0.51779417

15
5.- Graficar la relación, eje y- eje x, (𝑇𝐵𝑒𝑥𝑝 − 𝑥1 𝑒𝑥𝑝 , 𝑇𝑏 𝑒𝑥𝑝 − 𝑦1 exp )

Tb exp VS X1, y1
90
80
70
60
Tb exp

50
40 x1
30 y1
20
10
0
0 0.2 0.4 0.6 0.8 1 1.2
X1, Y1

6.- Graficar la relación, eje y- eje x, (𝑇𝐵 𝑏𝑖𝑏𝑙𝑖 − 𝑥1 𝑏𝑖𝑏𝑙𝑖 , 𝑇𝑏 𝑏𝑖𝑏𝑙𝑖 − 𝑦1 bibli )

T°C Bbl-x1, T°C-y1

95

90

85

80
x1
75 y1

70

65

60
0 0.2 0.4 0.6 0.8 1 1.2

16
7.- Graficar la relación, eje y- eje x, (𝑦1 𝑒𝑥𝑝 − 𝑥1 𝑒𝑥𝑝 ) trazando una línea en cuarenta
y cinco grados con respecto al origen.

Y1 exp Vs x1 exp vs 45
1.2

0.8
y1

0.6

0.4

0.2

0
0 0.2 0.4 0.6 0.8 1 1.2
X1

8.- Graficar la presión “eje y” – composición 𝑥1𝑒𝑥𝑝 “eje x” las siguientes series de
datos: presión total, presión parcial del componente uno, presión parcial del
componente 2, presión de vapor del componente 1, presión de vapor del
componente 2.

PT VS Pi1 VS Pi2 VS Pv1 VS Pv2

0.8
0.7
0.6
0.5
x1

0.4
0.3
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1 1.2
PT-Pi1-Pi2-Pv1-Pv2

17
9.- Calcular la presión de vapor de cada componente a partir de la ecuación de
Antoine.
Etanol n HEPTANO
A 8.2133 6.89385
B 1652.05 1264.37
C 231.48 216.636

Componente 1
𝐵
𝑃𝑠𝑎𝑡 = 10𝐴−(𝑇˚𝐶−𝐶 )

 Evento 1
1652.05
𝑃𝑠𝑎𝑡 = 108.2133−(81˚𝐶−231.4) = 844.1150403 𝑃𝑎

 Evento 2
1652.05
8.2133−( )
𝑃𝑠𝑎𝑡 = 10 74˚𝐶−231.4 = 638.6363685𝑃𝑎

 Evento 3

8.2133−1652.05
𝑃𝑠𝑎𝑡 = 10 65˚𝐶−231.48 = 437.6087757 𝑃𝑎

𝑃𝑝𝑎𝑟𝑐𝑖𝑎𝑙 = 𝑥1 ∗ 𝑃𝑠𝑎𝑡

 Evento 1
𝑃𝑝𝑎𝑟𝑐𝑖𝑎𝑙 = 0 ∗ 844.1150403𝑃𝑎 = 0 𝑃𝑎

 Evento 2
𝑃𝑝𝑎𝑟𝑐𝑖𝑎𝑙 = 0.024717747 ∗ 638.6363685𝑃𝑎 = 15.78565203 𝑃𝑎

 Evento 3
𝑃𝑝𝑎𝑟𝑐𝑖𝑎𝑙 = 0.11665195 ∗ 437.6087757𝑃𝑎 = 51.04791724 𝑃𝑎

Componente 2
 Evento 1

18
 𝐵
𝑃𝑠𝑎𝑡 = 10𝐴−(𝑇˚𝐶−𝐶 )

 Evento 1
1264.37
𝑃𝑠𝑎𝑡 = 106.89385−(81˚𝐶−216.636) = 442.3935317 𝑃𝑎

 Evento 2
1652.05
𝑃𝑠𝑎𝑡 = 108.2133−(74˚𝐶−231.4) = 349.5378644 𝑃𝑎

 Evento 3

8.2133−1652.05
𝑃𝑠𝑎𝑡 = 10 65˚𝐶−231.48 = 253.7894952 𝑃𝑎

𝑃𝑝𝑎𝑟𝑐𝑖𝑎𝑙 = 𝑥1 ∗ 𝑃𝑠𝑎𝑡

 Evento 1
𝑃𝑝𝑎𝑟𝑐𝑖𝑎𝑙 = 1 ∗ 442.3935317𝑃𝑎 = 442.3935317𝑃𝑎

 Evento 2
𝑃𝑝𝑎𝑟𝑐𝑖𝑎𝑙 = 0.975282253 ∗ 340.898076𝑃𝑎 = 340.898076 𝑃𝑎

 Evento 3
𝑃𝑝𝑎𝑟𝑐𝑖𝑎𝑙 = 0.88334805 ∗ 253.7894952𝑃𝑎 = 224.1844555 𝑃𝑎

10. calcular el coeficiente de actividad para cada componente a partir

de la ley de Raoult modificada

𝑦𝑖 𝑃 = 𝛾𝑖 𝑥𝑖 𝑃𝑠𝑎𝑡
𝑦𝑖 ∗ 𝑃
𝛾𝑖 =
𝑥𝑖 ∗ 𝑝𝑖𝑠𝑎𝑡
𝑃𝑠𝑖𝑠 = 585 𝑚𝑚𝐻𝑔
FRASCO 2
0.324531177 ∗ 585
𝛾1𝑒𝑥𝑝 = = 12.02679107
0.024717747 ∗ 638.6363685
0.675468823 ∗ 585
𝛾2𝑒𝑥𝑝 = = 1.159141952
0.975282253 ∗ 349.5378644

19
FRASCO 3
0.48220583 ∗ 585
𝛾1𝑒𝑥𝑝 = = 5.525992549
437.6087757 ∗ 0.11665195
0.51779417 ∗ 585
𝛾2𝑒𝑥𝑝 = = 1.351162323
253.7894952 ∗ 0.88334805
FRASCO 3
0.49729134 ∗ 585
𝛾1𝑒𝑥𝑝 = = 3.326992814
401.0897762 ∗ 0.21800844
0.50270866 ∗ 585
𝛾2𝑒𝑥𝑝 = = 1.595530978
235.7029003 ∗ 0.78199156
𝐺𝐸
A partir de la definición de en función de 𝛾 y etapas algebraicas,
𝑅𝑇
obtener los parámetros de iteración (𝐴12𝑒𝑥𝑝 , 𝐴21𝑒𝑥𝑝 ) del modelo de
Margules.
𝐺𝐸
= 𝑥𝑖 𝑙𝑛(𝛾𝑖 ) + 𝑥2 ln(𝛾2 )
𝑅𝑇
FRASCO 2
𝑙𝑛(𝛾𝑖 ) = 2.48713675
𝑙𝑛(𝛾2 ) = 0.147680035
𝐺𝐸
= 2.48713675 ∗ 0.024717747 + 0.975282253 ∗ 0.147680035
𝑅𝑇
= 0.205506133
FRASCO 3
𝑙𝑛(𝛾𝑖 ) = 1.709462878
𝑙𝑛(𝛾2 ) = 0.300965202
𝐺𝐸
= 1.709462878 ∗ 0.11665195 + 0.88334804 ∗ 0.300965202
𝑅𝑇
= 0.465269203
FRASCO 4
𝑙𝑛(𝛾𝑖 ) = 1.202068837
𝑙𝑛(𝛾2 ) = 0.467206582
𝐺𝐸
= 1.202068837 ∗ 0.21800844 + 0.78199156 ∗ 0.467206582
𝑅𝑇
= 0.627412756
 Evento 1
𝑥2 = 0.9752822532 = 0.951175473
2

𝑥 2 ∗ 𝑙 𝑛(𝛾2) 𝑒𝑥𝑝 = (0.951175473) ∗ (0.147680035) =0.14402972

20
  Evento 2
𝑥2 = 0.7803037772 = 0.265857024
2

𝑥 2 ∗ 𝑙 𝑛(𝛾2) 𝑒𝑥𝑝 = (0.300965202) ∗ (0.88334805) = 0.265857024

 Evento 3
𝑥2 = 0.781991562 = 0.6115108
2

𝑥 2 ∗ 𝑙 𝑛(𝛾2) 𝑒𝑥𝑝 = (0.78199156) ∗ (0.467206582) = 0.365351604

Tabla De Resultados

Fra 𝑥1 𝑒𝑥𝑝 𝑥2 𝑒𝑥𝑝 𝑇𝑏𝑒𝑥𝑝 ˚𝐶 𝑦1 𝑒𝑥𝑝 𝑦2 𝑒𝑥𝑝 𝑝1 𝑠𝑎𝑡 𝑝2 𝑠𝑎𝑡 ƴ1 𝑒𝑥𝑝 ƴ2 𝑒𝑥𝑝 𝐿𝑛ƴ1 𝑒𝑥𝑝 𝐿𝑛ƴ2 𝑒𝑥𝑝 𝐺𝐸
sco 𝑅𝑇
1 0 1 0 1 1.12539 0.58980 1 1.3603 0
4164 9534
2 0.02471 0.97 84.954 0.32453 0.675468 0.85144 0.46601 7.33292 0.706747 1.3582145 3.65558E- -0.289256241
7747 528 64737 1177 823 5132 2159 423 133 61 05
225
3
3 0.11665 0.88 73.939 0.48220 0.517794 0.58343 0.33835 2.59988 0.635698 1.3105069 0.0042372 -0.288726911
195 334 94595 583 17 0384 8166 3831 482 8
805
4 0.21800 0.78 64.698 0.49729 0.502708 0.53474 0.31424 1.54553 0.741195 1.1966133 0.02798024 -0.139284789
844 199 53696 134 66 2389 4689 6797 42 18 4
156
5 0.34735 0.65 64.698 0.53412 0.465879 0.50052 0.29708 1.12978 0.883622 0.9842449 0.11373449 -0.038360951
5685 264 53696 096 04 0899 6812 7118 698 4 7
431
5
6 0.45544 0.54 63.975 0.55206 0.447937 0.48953 0.29153 0.94516 1.077003 0.7729304 0.25687749 0.014713141
4034 455 72493 2115 885 3724 898 7925 978 77 4
596
6
7 0.60086 0.39 64.099 0.55206 0.447937 0.48953 0.29153 0.77045 1.580256 0.4786165 0.59077584 0.025945423
9366 913 60034 2115 885 3724 898 4232 853 34 2
063
4
8 0.71502 0.28 64.165 0.55206 0.447937 0.48953 0.29153 0.68309 2.335136 0.2693564 0.99626758 -0.03082957
4665 497 74139 2115 885 3724 898 7385 165 29
533
5
9 0.80701 0.19 64.038 0.56968 0.430312 0.48953 0.29153 0.65081 3.451881 0.1328986 1.43304040 -0.107542196
6595 298 35608 75 5 3724 898 8757 029 65 9
340
5
10 0.88272 0.11 64.340 0.56968 0.430312 0.48953 0.29153 0.61476 5.869135 0.0519245 1.87596456 -0.221931346
858 727 18387 75 5 3724 898 042 852 65 3
142
11 0.94613 0.05 66.199 0.60399 0.396009 0.55863 0.32612 0.62406 12.31030 0.0114598 2.31042071 -0.3108643
0933 386 76909 096 04 8504 156 6503 122 1 8
906
7

21
12 0.97945 0.02 68.894 0.75656 0.243437 0.66344 0.37730 0.76558 20.65000 0.0017054 2.56350861 -0.199418009
4491 054 73786 206 94 9598 598 7387 091 76 8
550
9
13 0.98988 0.01 70.229 0.82820 0.171790 0.72177 0.40521 0.83296 30.11469 0.0004163 2.64635652 -0.146475035
4066 011 15166 9042 958 7542 7428 0713 431 71 8
593
4
14 1 0 1 0 0.75252 0.41978 1 1 0 2.7284
2697 1887

11.- Obtener los parámetros de interacción 𝐴12 𝑒𝑥𝑝 , 𝐴21 𝑒𝑥𝑝 del modelo de
Margules a partir de la siguiente relación y de la tabla anterior.
𝐺𝐸
𝐶=
𝑅𝑇𝑋 2 1 𝑋2
 Evento 1
𝑋1 0.975282253
= = 39.45676209
𝑋2 0.024717747

0.205506133
𝐶= = 344.8869086
0.0247177472 0.975282253
 Evento 2

𝑋1 0.88334805
= = 7.572509897
𝑋2 0.11665195

0.465269203
𝐶 = 0.116651952 =38.7069047
0.88334805

 Evento 3

𝑋1 0.78199156
= = 3.586978372
𝑋2 0.21800844

0.627412756
𝐶= = 16.88125241
0.218008442 0.78199156

22
 constantes de Margules
400 y = 8.7339x - 4.8425
350
300
250
GE/RTx1^2x2

200
150
100
50
0
0 5 10 15 20 25 30 35 40 45
-50
x2/x1

𝐴12 𝑐𝑎𝑙 = 8.7339 𝐴21 𝑐𝑎𝑙 = −4.8425

12.- Comparar los valores de 𝐴12 𝑒𝑥𝑝 , 𝐴21 𝑒𝑥𝑝 con respecto a los bibliográficos:

Para el sistema:
B) Etanol (1)+n-heptano (2)
𝐴12 𝐵𝐼𝐵𝐿𝐼 = 1.3603 𝐴21 𝐵𝐼𝐵𝐿𝐼 = 2.7284
𝐴12 𝑐𝑎𝑙 = 8.7339 𝐴21 𝑐𝑎𝑙 = −4.8425
Se puede observar una gran diferencia entre los datos calculados y los datos
experimentales, que pudieron verse afectados debido a las condiciones
ambientales de la experimentación, así como las diferencia al obtener los datos,
ya que los datos experimentales se obtuvieron a una presión constante y os
bibliográficos a una temperatura constante, por lo que es normal una variación tan
notable entre estos.

23
13.- A partir de la ley de Raoult modificada para soluciones no ideales,
realizar cálculos de temperatura de burbuja para cada composición
utilizando el modelo de Margules y los parámetros de referencia.
Para evento 1
Cálculo de 𝐽1 𝑐𝑎𝑙

𝐵
𝑇𝑠𝑎𝑡 1 = (− − 𝐶)
𝐴 + ln(585)
1652.05
𝑇𝑠𝑎𝑡 1 = (−8.2133 − 231.48) = 71.86305508 ˚𝐶
+ ln(585)

1652.05
𝑇𝑠𝑎𝑡 1 = (−8.2133 − 231.48) = 85.52362037 ˚𝐶
+ ln(1.0189)

1652.05
𝑇𝑠𝑎𝑡 1 = (−8.2133 − 231.48) = 75.0083067 ˚𝐶
+ ln(665.342731)

𝐵
𝑇𝑠𝑎𝑡 2 = − ( − 𝐶)
𝐴 + ln(585)

1652.05
𝑇𝑠𝑎𝑡 2 = − (8.2133 − 231.48) = 89.75210606˚𝐶
+ ln(585)

1652.05
𝑇𝑠𝑎𝑡 2 = − (8.2133 − 231.48) = 89.75210606˚𝐶
+ ln(585)

1652.05
𝑇𝑠𝑎𝑡 2 = − (8.2133 − 231.48) = 89.75210606˚𝐶
+ ln(585)

𝑇 ∗ = (𝑋1 ∗ 𝑇𝑠𝑎𝑡 1 ) + (𝑋 2 ∗ 𝑇𝑆𝐴𝑇 2 )

𝑇 ∗ = (0 ∗ 71.86305508) + ( 1 ∗ 89.75210606) = 89.75210606 ˚𝐶

𝑇 ∗ = (𝑋1 ∗ 𝑇𝑠𝑎𝑡 1 ) + (𝑋 2 ∗ 𝑇𝑆𝐴𝑇 2 )

𝑇 ∗ = (0.11665195 ∗ 71.86305508) + ( 0.88334805 ∗ 89.75210606)

= 87.66531338 ˚𝐶
𝑇 ∗ = (𝑋1 ∗ 𝑇𝑠𝑎𝑡 1 ) + (𝑋 2 ∗ 𝑇𝑆𝐴𝑇 2 )

𝑇 ∗ = (0.21800844 ∗ 71.86305508) + (0.78199156 ∗ 89.75210606)

= 85.85214196 ˚𝐶

24
 𝐿𝑛 ƴ = (𝑥2 2 ∗ 𝐴12 ) + 2(𝐴21 − 𝐴12 ) ∗ (𝑥1 )

𝐿𝑛 ƴ1 = (12 ∗ 1.3603 ) + 2(2.7284 − 1.3603 ) ∗ (1 ) = 1.3603

𝐿𝑛 ƴ2 = (02 ∗ 1.3603 ) + 2(2.7284 − 1.3603 ) ∗ (1 ) = 0

𝐿𝑛 ƴ1 = (0.9752822532 ∗ 1.3603 ) + 2(2.7284 − 1.3603 )

∗ (0.024717747 ) = 1.358214561

𝐿𝑛 ƴ2 = (0.0247177472 ∗ 1.3603 ) + 2(2.7284 − 1.3603 )

∗ (0.975282253 ) = 3.65558𝐸 − 05

𝐿𝑛 ƴ1 = (0.883348052 ∗ 1.3603 ) + 2(2.7284 − 1.3603 ) ∗ (0.11665195 )

= 1.31050698

𝐿𝑛 ƴ2 = (0.116651952 ∗ 1.3603 ) + 2(2.7284 − 1.3603 ) ∗ (0.88334805 )

= 0.0042372

ƴ1 𝑐𝑎𝑙 = 𝐸𝑋𝑃(𝑙𝑛ƴ1 𝑐𝑎𝑙 )

ƴ1 𝑐𝑎𝑙 = 𝐸𝑋𝑃 (𝑙𝑛1.358214561 ) = 3.889243092

ƴ1 𝑐𝑎𝑙 = 𝐸𝑋𝑃 (𝑙𝑛1.31050698 ) = 3.708053144

ƴ1 𝑐𝑎𝑙 = 𝐸𝑋𝑃 (𝑙𝑛1.196613318 ) = 3.308891762

𝐵
𝑇𝑛1 = (− )−𝐶
𝐴 + 𝐿𝑜𝑔(ƴ1 𝑐𝑎𝑙 )

 Evento 1

1652.05
𝑇𝑛1 = (− ) − 231.48 = 85.52362037 ˚𝐶
8.2133 + 𝐿𝑜𝑔(3.889243092 )

 Evento 2

1652.05
𝑇𝑛1 = (− ) − 231.48 = 75.0083067 ˚𝐶
8.2133 + 𝐿𝑜𝑔(3.708053144 )

25
  Evento 3

1652.05
𝑇𝑛1 = (− ) − 231.48 = 68.92076946 ˚𝐶
8.2133 + 𝐿𝑜𝑔(3.308891762 )

Para evento 1
Cálculo de ƴ2 𝑐𝑎𝑙

ƴ2 𝑐𝑎𝑙 = 𝐸𝑋𝑃(𝑙𝑛ƴ2 𝑐𝑎𝑙 )

ƴ2 𝑐𝑎𝑙 = 𝐸𝑋𝑃 (ln 0) = 1

ƴ2 𝑐𝑎𝑙 = 𝐸𝑋𝑃(𝑙𝑛3.65558𝐸 − 05 ) = 1.000036556

ƴ2 𝑐𝑎𝑙 = 𝐸𝑋𝑃 (𝑙𝑛0.0042372 ) = 1.00424619

𝐵
𝑇𝑛1 = (− )−𝐶
𝐴 + 𝐿𝑜𝑔(ƴ1 𝑐𝑎𝑙 )

 Evento 1

1652.05
𝑇𝑛1 = (− ) − 231.48 = 85.52362037 ˚𝐶
8.2133 + 𝐿𝑜𝑔(3.889243092 )
 Evento 2

1652.05
𝑇𝑛1 = (− ) − 231.48 = 75.0083067 ˚𝐶
8.2133 + 𝐿𝑜𝑔(3.708053144 )

 Evento 3

1652.05
𝑇𝑛1 = (− ) − 231.48 = 68.92076946 ˚𝐶
8.2133 + 𝐿𝑜𝑔(3.308891762 )

26
 Tabla De Resultados Calculados
Fra 𝑥1 𝑒𝑥𝑝 𝑥2 𝑒𝑥𝑝 𝑇𝐵 𝑐𝑎𝑙 ƴ1 𝑐𝑎𝑙 ƴ2 𝑐𝑎𝑙 ln(𝑦1 𝑐𝑎𝑙 )ln(𝑦2 𝑐𝑎𝑙 ) 𝑦1 𝑐𝑎𝑙 𝑦2 𝑐𝑎𝑙
sco
1 0 1 89.7521 0 1 1.3603 0 0 1
0606
2 0.02471 0.97528 84.9546 0.16151 0.83852 1.35821 3.65558 0.02471 0.97528
7747 2253 4737 4522 8615 4561 E-05 7747 2253
3 0.11665 0.88334 73.9399 0.47106 0.52895 1.31050 0.00423 0.11665 0.88334
195 805 4595 6047 5689 698 72 195 805
4 0.21800 0.78199 64.6985 0.61190 0.38810 1.19661 0.02798 0.21800 0.78199
844 156 3696 098 7195 3318 0244 844 156
5 0.34735 0.65264 64.6985 0.68625 0.31374 0.98424 0.11373 0.34735 0.65264
5685 4315 3696 6819 6851 494 4497 5685 4315
6 0.45544 0.54455 63.9757 0.70587 0.29412 0.77293 0.25687 0.45544 0.54455
4034 5966 2493 5693 6976 0477 7494 4034 5966
7 0.60086 0.39913 64.0996 0.69758 0.30241 0.47861 0.59077 0.60086 0.39913
9366 0634 0034 5868 6684 6534 5842 9366 0634
8 0.71502 0.28497 64.1657 0.67532 0.32468 0.26935 0.99626 0.71502 0.28497
4665 5335 4139 2427 6226 6429 758 4665 5335
9 0.80701 0.19298 64.0383 0.66130 0.33870 0.13289 1.43304 0.80701 0.19298
6595 3405 5608 0144 2931 8665 0409 6595 3405
10 0.88272 0.11727 64.3401 0.67592 0.32411 0.05192 1.87596 0.88272 0.11727
858 142 8387 5605 9577 4565 4563 858 142
11 0.94613 0.05386 66.1997 0.75391 0.24609 0.01145 2.31042 0.94613 0.05386
0933 9067 6909 5455 8485 981 0718 0933 9067
12 0.97945 0.02054 68.8947 0.86677 0.13322 0.00170 2.56350 0.97945 0.02054
4491 5509 3786 5022 6126 5476 8618 4491 5509
13 0.98988 0.01011 70.2291 0.92528 0.07472 0.00041 2.64635 0.98988 0.01011
4066 5934 5166 57 4057 6371 6528 4066 5934
14 1 0 71.8630 1 0 0 2.7284 1 0
5508

27
 14.- Graficar la relación, eje y- eje x (𝑇𝑏 𝑐𝑎𝑙 − 𝑥1 𝑒𝑥𝑝 , 𝑇𝑏 𝑐𝑎𝑙 − 𝑦1 𝑐𝑎𝑙 )

Tbcal-x1exp, Tbcal-y1cal
95

90

85

80
Tbcal-x1exp
75 Tbcal-y1cal

70

65

60
0 0.2 0.4 0.6 0.8 1 1.2

15.- Graficar la relación, eje y- eje x (𝑦1 𝑐𝑎𝑙 − 𝑥1 𝑒𝑥𝑝 )

y1cal vs x1exp
1.2

0.8

0.6

0.4

0.2

0
0 0.2 0.4 0.6 0.8 1 1.2

28
16. Calcula la desviación entre la composición de la fase vapor experimental 𝑦1 𝑒𝑥𝑝
y la calculada 𝑦1𝑐𝑎𝑙 .
𝜎 = 𝑦1 𝑒𝑥𝑝 −𝑦1𝑐𝑎𝑙

FRASCO 2 𝜎 = 0.33245312 − 0.1615145 = 0.1630167

FRASCO 3 𝜎 = 0.4822058 − 0.471066 = 0.0111398
FRASCO 4 𝜎 = 0.4972913 − 0.611901 = −0.11461

17. Calcular la desviación entre la temperatura de burbuja experimental 𝑇𝐵 𝑒𝑥𝑝 y

calculada 𝑇𝐵 𝑐𝑎𝑙 .
𝜎 = 𝑇𝐵 𝑒𝑥𝑝 − 𝑇𝐵 𝑐𝑎𝑙
FRASCO 1 𝜎 = 81 − 89.752106 = −8.752106
FRASCO 2 𝜎 = 74 − 84.954647 = −10.95465
FRASCO 3 𝜎 = 65 − 73.939946 = −8.939946

Discusión De Resultados
Análisis de la gráfica 1) mediante esta grafica pudimos obtener una ecuación de
ajuste que relaciona los datos del índice de refracción y la fracción mol liquida y que
a su vez mediante la sustitución del índice de refracción vapor nos ayuda a obtener
la fracción mol de la sustancia en fase vapor.
Análisis grafica 2) Por medio de esta grafica nos podemos percatar de cómo es que
a una misma temperatura se pueden presentar distintas fracciones mol de un
componente, dependiendo de la fase en la que se busque, además en ella se
observa como en cierta parte de la gráfica a una cierta temperatura la fracción mol
liquida es igual a la fracción mol vapor, este punto es conocido como azeótropo y
su condición de temperatura seria de 59°C basándonos en la grafica
Análisis 3) en esta grafica se observa el comportamiento que el fenómeno presenta
bibliográficamente y observamos que es muy semejante al generado mediante los
cálculos de la practica por lo cual esta grafica nos sirve como un indicador de que
los cálculos se están realizando correctamente.
Análisis 4) mediante la línea trazada desde el origen hasta 1 nos podemos percatar
de que el comportamiento de la mezcla que estamos trabajando difiere mucho de la
idealidad, pero que a su vez hay algunos puntos donde la curva genera por su
comportamiento atraviesa dicha línea y por ello podemos inferir que hay ciertos
valores donde se pude llegar a tener un comportamiento ideal y a su vez nos ayuda
con el cálculo del azeótropo por el método gráfico.
Análisis 5) mediante esta grafica se puede observar cómo es que la composición de
cada componente afecta a la presión parcial del otro y también se observa cómo
entre más grande sea la fracción del componente uno mayor será la presión total
generada en el sistema.

29
Análisis 6) Ahora bien al realizar esta grafica comprobamos la eficiencia del método
de iteraciones para la obtención de y1 y de tburb puesto que el comportamiento que
sigue es muy parecido al comportamiento experimental y al bibliográfico lo cual
refleja una buena adaptación al fenómeno.
Análisis 7) En esta última grafica de nuevo podemos observar cómo es que el
comportamiento de nuestra mezcla se aleja de la idealidad, pero esta vez una
manera mejor definida en cuanto a las curvas generadas y así como en la ocasión
anterior se puede llegar observa que este comportamiento chocara en determinado
punto con la idealidad.

CUESTIONARIO
¿Para este sistema, es posible la formación de un azeótropo? Demostrar las
condiciones a las cuales ocurre por método gráfico.
A través de la gráfica de T B vs X1, TB vs Y1 se puede observar la formación de un
azeótropo al juntarse en un punto el etanol y el n-heptano, y haciendo un análisis
más detallado de éste, se puede observar gráficamente el comportamiento de un
azeótropo de temperatura de ebullición mínima, esto se visualiza al notar que el
mínimo de la curva T-y coincide con el mínimo de la curva T-x.
En caso de que exista la formación de un azeótropo, ¿por qué es positivo o
negativo?
Es positivo, pues como ya se mencionó en la pregunta anterior, la gráfica tiene el
comportamiento de un azeótropo de temperatura de ebullición mínima.
Determinar por un método analítico las condiciones de formación del
azeótropo a presión constante, realizar los cálculos detalladamente y reportar
los resultados.
¿El sistema se comporta idealmente en ambas fases? explicar.
No, pues este sistema binario en equilibrio liquido-vapor se desarrolla con la ley de
Raoult modificada precisamente porque la fase vapor es no ideal y la fase líquida
es no ideal, esto debido a que nuestro azeótropo presenta desviaciones positivas
de la idealidad.
¿Qué valor tiene la constante de equilibrio “K” en un punto azeotrópico y en
el punto crítico de una mezcla binaria?
El valor de la constante de equilibrio K es igual a uno en ambos puntos.
Con base a la pregunta anterior y con respecto al valor de la constante de
equilibrio “K”, ¿se puede considerar el punto crítico de una mezcla binaria
como un azeótropo? Explicar.

30
Si nos basamos en el valor de la constante de equilibrio podría decirse que sí,
además de que punto crítico cumple con algunas de las características de lo que es
el azeótropo, sin embargo, en un punto azeótropo, existen dos fases en equilibrio
termodinámico formalmente: liquido y vapor, y las puedes ver en el frasco las
composiciones dadas. No existe un solo punto azeótropo, ya que el diagrama T-x a
P constante ves uno, pero si cambias la P encuentras otro, y así sucesivamente.
Cada fase tiene una composición fija en el azeótropo.
Un punto crítico no es un estado de coexistencia de fases. De hecho, la
termodinámica no lo describe, ya que, en ese punto, las densidades de las dos fases
fluctúan enormemente, para, aun valor diferencialmente mayor, de T o P, volverse
una sola fase. Físicamente son distintos.
Determinar por un método analítico las condiciones de formación del
azeótropo a presión constante, realizar los cálculos detalladamente y reportar
los resultados.
Para las condiciones de formación del azeótropo a presión constante se ocupa la
𝑌 𝑌
𝐾𝑖 = 𝑖 ; donde nos dice que 𝐾𝑖 = 𝑖 = 1 es el punto del azeótropo ya que al
𝑋𝑖 𝑋𝑖
𝑃𝑠𝑎𝑡 𝑌𝑖
despejar nos dice que 𝐾𝑖 = 𝑃𝑇
, ya que al final se verá el azeótropo el cual hará
que cada uno de los puntos poseen por uno solo, y eso determinara la obtención de
fórmulas como la de 𝑃𝑇 = 𝑃1 𝑠𝑎𝑡 𝑌1 como ecuación 5 y 𝑃𝑇2 = 𝑃2 𝑠𝑎𝑡 𝑌2 como ecuación
6, ya que al ver cómo actúa la presión como la temperatura para poder obtener
𝑋1 𝑦 𝑋2 y obtener las interacciones a ocupar y la obtención de 𝑌1 , ya que nos ayuda
también a despejar y obtener el azeótropo.

Anexo cálculos de azeótropo

𝑌𝑖 𝑌1 𝑃1 𝑠𝑎𝑡 𝑌1 𝑌2 𝑃2 𝑠𝑎𝑡 𝑌2
𝐾𝑖 = Por lo tanto 𝐾1 = =1= 𝐾2 = =1=
𝑋𝑖 𝑋1 𝑃𝑇 𝑋2 𝑃𝑇

𝑌1
𝐾1 = =1
𝑋1

𝑌2
𝐾2 = =1
𝑋2

𝑃𝑇 = 𝑃1 𝑠𝑎𝑡 𝑌1 ……. 5 𝑥2 + 𝑥1 = 1

𝑃𝑇2 = 𝑃2 𝑠𝑎𝑡 𝑌2 ……6 𝑥2 = 1 − 𝑥1

31
 𝐵
𝑃𝑠𝑎𝑡 = 10𝐴− 𝑇+𝐶
2 [ 𝐴 +2(𝐴 − 𝐴 ) 𝑥 ]
Y= 𝑒 𝑥2 12 21 12 1 = 0.75

𝐴12 = 1.3603
𝐴21 = 2.7284
𝐵1 2
𝐴1 − ∗ 𝑒 𝑥2 [ 𝐴12+2(𝐴21− 𝐴12) 𝑥1]
0.75=10 𝑇+𝐶1

𝐵1 2
𝐴2 − ∗ 𝑒 𝑥2 [ 𝐴12+2(𝐴21− 𝐴12) 𝑥2]
0.75=10 𝑇+𝐶1

Etanol n-heptano
A=8.2133 A=6.89385
B=1652.05 B=1264.37
C=231.48 C=216.636

2 [ 1.3603+2(2.7284− 1.3603) 𝑥]
0.75= 𝑒 (1−𝑥)
0.75= 𝑒 (1−𝑥)(1−𝑥)[ 1.3603+2(2.7284− 1.3603) 𝑥]
2 )[ 1.3603+2(2.7284− 1.3603) 𝑥]
0.75= 𝑒 (1−2𝑥+𝑥
2 )(1.3603+2.7362𝑥)
0.75= 𝑒 (1−2𝑥+𝑥
2 +1.3603𝑥 2 +2.7362𝑥 3 )
0.75= 𝑒 (1.3603+2.7362𝑥)−2.7206𝑥−5.4724𝑥
3)
0.75= 𝑒 (1.3603+0.0156−4.1121+2.7362+𝑥
X= 0.9085602964

1652.05 2
𝑒 ((1−(0.9085602964) ) (1.3603+2(2.7284−1.3603)(0.9085602964))
0.75= 10(8.2133−𝑇+231.48) ∗
1652.05 ((1−(0.9085602964)2) (1.3603+2(2.7284−1.3603)(0.9085602964))
Log 10(8.2133−𝑇+231.48) ∗ 𝑒 = log 0.75
1652.05
(8.2133 − 𝑇+231.48) 𝑒 −2.175052909 = log 0.75
187.6764202
0.9330484806 - = log 0.75
𝑇+231.48
187.6764202
0.9330484806 - = -0.1249387366
𝑇+231.48

32
 187.6764202
- = -1.057987217
𝑇+231.48

-187.6764202= -1.057987217 * 𝑇 + 231.48

419.1564202=-1.057987217T
T= 396.1828803
2
Sustituyendo x en 𝑒 (1−𝑋1 ) [ 𝐴12 + 2(𝐴21 − 𝐴12 ) 𝑥1 ]
2
𝑒 (1−(0.9085602964) ) (1.3603 + 2(2.7284 − 1.3603)(0.9085602964) = 0.75
𝐵1 2
𝐴1 − ∗ 𝑒 𝑥2 [ 𝐴12+2(𝐴21− 𝐴12) 𝑥1]
Sustituyendo en 0.75=10 𝑇+𝐶1

1652.05
10(8.2133−396.1828803+231.48) ∗ 0.75

0.75=0.7459

𝐵1 2
𝐴2 − ∗ 𝑒 𝑥2 [ 𝐴12+2(𝐴21− 𝐴12) 𝑥2]
0.75=10 𝑇+𝐶1

2 [ 1.3603+2(2.7284− 1.3603) 𝑥]
0.75= = 𝑒 (1−𝑥)
2 [ 1.3603+2(2.7284− 1.3603) 𝑥]
In (0.75) = In (𝑒 (1−𝑥) )
In (0.75) = (1 − 𝑥)2 [ 1.3603 + 2(2.7284 − 1.3603) 𝑥 ]𝐼𝑛 (𝑒)
In (0.75) = (1 − 𝑥)2 [ 1.3603 + 2(2.7284 − 1.3603) 𝑥 ]
X2= -0.07022630842
1264.37 2
𝑒 ((1−(−0.07022630842) ) (1.3603+2(2.7284−1.3603)(−0.07022630842))
0.75= 10(6.89385−𝑇+216.636) ∗
1264.37 ((1−(−0.07022630842)2 ) (1.3603+2(2.7284−1.3603)(−0.07022630842))
Log 10(6.89385−𝑇+216.636) ∗ 𝑒 = log 0.75
1264.37
(6.89385 − ) 𝑒 1.162385785 = log 0.75
𝑇+216.363
4042.889906
1.754299233 - = log 0.75
𝑇+216.363
4042.889906
1.754299233 - = -0.1249387366
𝑇+216.636
4042.889906
- = -1.87923797
𝑇+216.636

-4042.889906= -1.87923797* 𝑇 + 216.636

-4259.525906=-1.87923797T
T= 226.624011

33
 2
Sustituyendo x en 𝑒 (1−𝑋2 ) [ 𝐴12 + 2(𝐴21 − 𝐴12 ) 𝑥2 ]
2 ) (1.3603+2(2.7284−1.3603)(−0.07022630842))
𝑒 ((1−(−0.07022630842) = 0.75
𝐵1 2
𝐴2 − ∗ 𝑒 𝑥2 [ 𝐴12+2(𝐴21− 𝐴12) 𝑥2]
Sustituyendo en 0.75=10 𝑇+𝐶1

1264.37
10(6.89385−226.624011+216.636) ∗ 0.75

0.75=0.75064

Observaciones Y Recomendaciones

Se pudo observar durante la práctica el manejo del refractómetro para poder medir
la refracción de la mezcla binaria de cada frasco, se observó cómo se obtenían las
gotas de la condensación y así obtener la fase vapor, después de haber obtenido
esto se realizó el calculó de su índice de refracción de la fase vapor y después
calcular el índice de refracción de la fase líquida.

Conclusiones
Mediante cálculos matemáticos, se determinó el equilibrio líquido-vapor para una
mezcla binaria de etanol + n-heptano, para poder obtener los cálculos se tuvo que
utilizar la Ley de Raoult modificada y la Ley de Margules o Van Laar para obtener
los parámetros de interacción. Se determinaron las corrosiones y los coeficientes
de actividad, donde al momento de graficar cada uno de los puntos de temperatura
de rocío y burbuja se podía observar el comportamiento que cada uno de ellos tenía.
Se pudo comprobar así mismo el punto azeotrópico de la solución de la composición
“x” y “y” donde dio como resultado que son iguales.
Al término de la práctica pudimos darnos cuenta de que se cumple con los objetivos
de la práctica y al final de los cálculos hechos por la Ley de Raoult modificada nos
dio el resultado de 585 mmHg.

Por medio de la práctica anterior se pudo determinar el equilibrio líquido-vapor de la

mezcla etanol y n-heptano, todo esto mediante cálculos matemáticos a través de la

34
ley de Raoult modificada. También pudimos obtener la temperatura y composición
de nuestro azeótropo formado, por medio del método de iteración para calcular las
temperaturas de burbuja y así poder hacer la comparación con la temperatura de
burbuja experimental. A partir de las gráficas obtenidas me pude percatar de que el
azeótropo obtenido es de temperatura de ebullición mínima, por lo tanto, es un
azeótropo positivo, esto debido a que nuestro sistema presenta desviaciones
positivas de la idealidad. Asimismo, se calculó el coeficiente de actividad del etanol
y del n-heptano a partir del modelo de Margules o van Laar para así hacer un
análisis con respecto a A12, A21 bibliográficos. Esta práctica fue muy similar a la
práctica 4, pues en la ley de Raoult modificada al tener el coeficiente de actividad
se puede decir que tiene una mejor precisión y un mejor ajuste en gráficas que la
ley de Raoult, sin embargo, se sabe ya que una es para fases no ideales y la otra
ley para ideales respectivamente. Es necesario recalcar la importancia de la ley de
Raoult modificada para sistemas con n-componentes donde las fases líquido y
vapor se alejan de la idealidad, ya sean desviaciones positivas o negativas y de
igual manera reconocer que los azeótropos en la industria nos permiten realizar
mezclas más estables para transportar y almacenar.

-Melo Lugo Denisse Guadalupe

Esta práctica partió por así decirlo de la anterior ya que se tomaron en consideración
algunos conceptos anteriores incluyendo unos más para poder desarrollar la parte
experimental, como son los métodos de resolución, que si bien existen muchos para
cada tipo de área, en este caso se usó los modelos que existen para termodinámica,
donde de la misma manera se dividían en más como son los métodos especiales
usados en la industria del petróleo, para poder saber las ecuaciones de estado y
finalmente los que nos compete en esta práctica que son los que no ayudan en la
determinación de los coeficientes de actividad, siendo el caso de usar el de
Margules y esto debido a que nuestro sistema se encontraba en no ideal ambas
fases que son el vapor y la liquida, y todo cambia con respecto a eso, aquí se vuelve
hacer el uso del refractómetro, el saber el comportamiento de las sustancias para
su comprobación, ya que de esa manera sabremos si vamos por el camino correcto
y además de una pieza clave que se menciona en la introducción y que viene
repercutiendo por toda la práctica que es el azeótropos, que para en este caso si
son de gran utilidad, aunque se tuvieron que hacer bastantes cálculos, bastantes

35
iteraciones, y bastantes gráficas, aunque principalmente aquí fue el hecho de haber
entendido el tema más que nada.

Al término de esta práctica se logró observar el comportamiento de un sistema en

equilibrio liquido-vapor obtener la temperatura y la composición de nuestro
azeótropo a través del método de iteración para obtener sus distintas temperaturas
de burbuja y compararla con las temperaturas experimentales y obtener sus
porcentajes de error, se logró observar por medio de las gráficas el comportamiento
de un azeótropo ya que en nuestro sistema representa diversas variaciones a un
sistema ideal, también se logró comprobar matemáticamente por medio de la
ecuación de Raoult en nivel de actividad de un sistema así como sus fracciones de
vapor para un sistema con presión contante .

Al realizar esta práctica logre concluir que los diagramas obtenidos a partir de los
datos calculados para el equilibrio de la mezcla etanol con heptano son muy útiles
ya que gracias a ellos se puede dar una estimación de cuáles serán las condiciones
en la cuales se forme un azeótropo, así mismo logre aprender que para la mayoría
de condiciones de temperatura los componentes de una mezcla no ideal se
comportan como 2 compuestos separados pero que ha ciertas condiciones se
pueden comportar como un solo compuesto. además de eso logré aprender que la
ley de Raoult modificada es muy importante ya que nos ayuda a conocer el
coeficiente de actividad y además de eso combinada con el modelo de solución de
Margules nos ayuda a calcular los valores de temperatura de burbuja y fracciones
vapor de una mezcla no ideal a presión constante, también sobre esta ley y el
modelo de margulles aprendí que son muy exactos ya que al comparar los
resultados obtenidos experimentalmente con los que se obtienen mediante los
cálculos de la práctica se observa que la diferencia es mínima lo cual refleja que
hay un buena representación matemática del fenómeno. Así como también
analizando la gráfica de presiones logre inferir que entre más cantidad de
componente 1 haya en la mezcla mayor terminara siendo la presión total dentro del
sistema. Por último logre comprender que el cálculo analítico del azeótropo es una
forma más viable de calcular el azeótropo puesto que muestra mayor exactitud que
el método gráfico.

36
 BIBLIOGRAFIA

[1] Kaplan, T. A. (2006, 16 marzo). The Chemical Potential. Journal of Statistical

Physics. https://link.springer.com/article/10.1007/s10955-005-8067-
x?error=cookies_not_supported&code=46725d85-f752-472d-b468-
60ed7b1a2539

[2] Roo, R. (2019, 11 julio). Termodinámica de las soluciones. SlideShare.

https://www.slideshare.net/RobertoRoo1/termodinmica-de-las-soluciones

[3]García-Sánchez, F., & N. Justo-García, D. (2011, 8 julio). Representación de los

Equilibrios Líquidos Multifásicos de Sistemas Cuaternarios con los Modelos
NRTL y UNIQUAC. Daimler N. Justo-García.
https://scielo.conicyt.cl/pdf/infotec/v23n1/art14.pdf

[4] Alfaro, G. M. B. (2011, 12 diciembre). Equilibrio de fases en sistemas de

multicomponentes. Slideshare.
https://es.slideshare.net/AngBarraza/equilibrio-de-fases-en-sistemas-de-
multicomponentes#:%7E:text=%EF%82%A7Los%20dos%20tipos%20de,qu
e%20se%20presenta%20cada%20especie.

[5] Hyperphysics. (2000). Retrieved 30 March 2021, from http://hyperphysics.phy-

astr.gsu.edu/hbasees/thermo/heatra.html

[6] Hernández Segura, G. INTRODUCCIÓN A LA TERMODINÁMICA. Consultado

el 30 de marzo de 2021 en
http://depa.fquim.unam.mx/~fermor/blog/programas/claset1.pdf

[7]Vargas, F., Vera, F., López, A., Orozco, G., & Raya, M. (1997, 15 septiembre).
Equilibrio líquido vapor - Biosep liq-liq. Biosep liq-liq.
https://sites.google.com/site/biosepliqliq/about/teaching-
stories#:%7E:text=En%20una%20mezcla%20binaria%20de,mostrados%20
en%20la%20figura%201.

[8] E. (2021, 1 marzo). ¿Qué es el punto de saturación? QueSignificado.org.

https://quesignificado.org/que-es-el-punto-de-saturacion/

[9] Vargas, F., Vera, F., López, A., Orozco, G., & Raya, M. (1997, 15 septiembre). Ley de
Henry y Ley de Raoult. Biosep liq-liq.
https://sites.google.com/site/biosepliqliq/about/teaching-stories

[10] Noguera, I. B. (2020, 14 octubre). Ley de Raoult. Ingeniería Química Reviews.

https://www.ingenieriaquimicareviews.com/2020/09/ley-de-raoult.html

37
[11] Ríos A., L. G. (2004). MODELOS TERMODINÁMICOS PARA EL EQUILIBRIO
VAPOR - LIQUIDO A BAJAS PRESIONES: FASE LIQUIDA, MODELO DE
WILSON. Redalyc.org. https://www.redalyc.org/articulo.oa?id=84912053034

[12] Noguera, I. B. (2020b, diciembre 15). ¿Qué es un azeótropo o mezcla azeotrópica?

IngenierÃa QuÃmica Reviews.
https://www.ingenieriaquimicareviews.com/2020/12/que-es-un-azeotropo-o-
mezcla-azeotropica.html

[13] Pereda, S. (2014). Ingenieria del Equilibrio entre Fases. Course Hero.
https://www.coursehero.com/file/p5rapr5h/proceso-puede-dividirse-en-dos-
%C3%A1reas-una-operando-por-encima-de-la-temperatura/

[14]Químico, L. T. P. (2017, 24 agosto). Refractómetro. TP - Laboratorio Químico.

https://www.tplaboratorioquimico.com/laboratorio-quimico/materiales-e-
instrumentos-de-un-laboratorio-quimico/refractometro.html

[15] Van Ness, Hendrick C .; Soczek, CA ; Peloquin, GL ; Machado, RL ,

Propiedades de exceso termodinámico de tres sistemas alcohol-hidrocarburo , J.
Chem. Ing. Datos , 1967, 12, 2, 217-224, https://doi.org/10.1021/je60033a017

[16] Ríos A., L. G. (2004). MODELOS TERMODINÁMICOS PARA EL EQUILIBRIO

VAPOR - LIQUIDO A BAJAS PRESIONES: FASE LIQUIDA, MODELO DE
WILSON. Redalyc.org. https://www.redalyc.org/articulo.oa?id=84912053034

[17]Salazar, F. (2015, 27 octubre). FUGACIDAD Y ACTIVIDAD. Issuu.

https://issuu.com/fsalazar2011/docs/equilibrio_21

[18] B.E. Poling J. M. Prausnitz, J.P.O´Connell, The Properties of Gases and

Liquids. 5a. Edición, Mc Graw-Hill, ISBN: 978-0-07-011682-5.

[19] Journal, C. (1993). Chemical Engineering Journal. Retrieved 26 May 2021,

from https://www.journals.elsevier.com/chemical-engineering-journal

38