INSTITUTO POLITÉCNICO NACIONAL.
ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS
EXTRACTIVAS.
DEPARTAMENTO DE FORMACIÓN BÁSICA.
ACADEMIA DE FISICOQUÍMICA BÁSICA.
LABORATORIO DE TERMODINÁMICA DEL EQUILIBRIO DE FASES.
Práctica 4. EQUILIBRIO LIQUIDO - VAPOR DE UNA
SOLUCION NO IDEAL A BAJA PRESION.
INTEGRANTES DEL EQUIPO:
• Aguilar Ortiz Paulina Isabel
• Angeles Soto Jonathan
• García García Samuel
• Ramírez Cordero Isaac
2IM32
Equipo 1.
Fecha de entrega: martes 15 de mayo 2022
Objetivos
1.- Preparar una serie de soluciones de cloroformo (1) y metanol (2).
Medir a cada una su temperatura de burbuja y la composición de los vapores, a
presión atmosférica.
2.- A partir de los datos experimentales, determinar el coeficiente de actividad
(yi) de cada Componente en cada solución.
3.- Investigar el modelo de solución (Margules o Van Laar) al que se ajusta el
sistema.
4.- Corroborar el punto 3, elaborando los cálculos correspondientes y
comparando los resultados con los datos experimentales.
Introducción Teórica
En termodinámica e ingeniería química, el equilibrio vapor-líquido (VLE) describe la
distribución de una especie química entre la fase de vapor y una fase líquida.
La concentración de un vapor en contacto con su líquido, especialmente en
equilibrio, se expresa a menudo en términos de presión de vapor, que será una
presión parcial (una parte de la presión total del gas) si hay otros gases presentes
con el vapor. La presión de vapor de equilibrio de un líquido depende en general en
gran medida de la temperatura. En el equilibrio vapor-líquido, un líquido con
componentes individuales en ciertas concentraciones tendrá un equilibrio de vapor
en el cual las concentraciones o presiones parciales de los componentes de vapor
tienen ciertos valores dependiendo de todas las concentraciones del componente
líquido y de la temperatura. Lo contrario también es cierto: si un vapor con
componentes a ciertas concentraciones o presiones parciales se encuentra en
equilibrio vapor-líquido con su líquido, entonces las concentraciones de
componentes en el líquido se determinarán dependiendo de las concentraciones de
vapor y de la temperatura. La concentración de equilibrio de cada componente en
la fase líquida es a menudo diferente de su concentración (o presión de vapor) en
la fase de vapor, pero existe una relación. Los datos de concentración de VLE
pueden ser determinados experimentalmente, o computados o aproximados con la
ayuda de teorías tales como ley de Raoult, ley de Dalton, y ley de Henry.
Tal información de equilibrio vapor-líquido es útil en el diseño de columnas para
destilación, especialmente destilación fraccionada, que es una especialidad
particular de ingeniero químicos. La destilación es un proceso utilizado para separar
o separar parcialmente los componentes de una mezcla por ebullición
(vaporización) seguida de condensación. La destilación aprovecha las diferencias
de concentración de los componentes en las fases líquida y vapor.
En las mezclas que contienen dos o más
componentes, las concentraciones de cada
componente se expresan a menudo como
fracciones molares. La fracción molar de un
componente dado de una mezcla en una fase
particular (ya sea la fase de vapor o la fase líquida)
es el número de moles de ese componente en esa
fase dividido por el número total de molares de
todos los componentes en esa fase.
Las mezclas binarias son las que tienen dos componentes. Las mezclas de tres
componentes se denominan mezclas ternarias. Puede haber datos VLE para
mezclas con aún más componentes, pero tales datos son a menudo difíciles de
mostrar gráficamente. Los datos VLE son una función de la presión total, como 1
atm o cualquier presión a la que se lleve a cabo el proceso.
Cuando se alcanza una temperatura tal que la suma de las presiones de vapor de
equilibrio de los componentes líquidos se vuelve igual a la presión total del sistema
(de lo contrario es menor), entonces las burbujas de vapor generadas por el líquido
comienzan a desplazar al gas que mantenía la presión total, y se dice que la mezcla
Hierve. Esta temperatura se denomina "punto de ebullición" de la mezcla líquida a
la presión dada. (Se supone que la presión total se mantiene constante ajustando
el volumen total del sistema para adaptarse a los cambios de volumen específicos
que acompañan a la ebullición.) El punto de ebullición a una presión global de 1 atm
se denomina "punto de ebullición normal" de la mezcla líquida.
Desarrollo Experimental
Preparar una serie de Medir a cada solución el Colocar la primera
soluciones de metanol (1) índice de refracción ἠ y solución en el recipiente
e isopropanol (2) En la elaborar una gráfica de ἠ del equipo. Poner el
que la concentración del vs X1 la Cual se llama recipiente sobre la parrilla
metanol (1) varíe de 0,0.1 Curva de Calibración y comenzar agitación
hasta 1 fracción lenta.
Depositar la muestra en Del depósito de Calentar el líquido hasta
un tubo de ensaye condensado que tiene el ebullición y esperar un
previamente etiquetado equipo en la base del momento para asegurare
tapándolo refrigerante tomar una que la temperatura sea
inmediatamente muestra por medio de constante y anotarla. (Es
una pipeta delgada y una la temperatura de burbuja)
perilla
Esperar a que las Con ayuda de la curva de
Repetir los pasos para muestras que contienen calibración determina la
cada solución el condensado estén a la fracción mol del metanol
temperatura ambiente en los vapores en el
para medir el índice de momento en que se dio el
refracción equilibrio líquido-vapor.
Llenar la tabla de datos
experimentales.
Cálculos
Compuesto Sustancia Peso Molecular Densidad
(PM) g/mol g/ml
1 CHCl3 119.38 1.483
Cloroformo
2 CH3OH 32 0.791
Metanol
Tabla de datos experimentales
ml ml X1 X2 ἠ (liq) T(C) ἠ(cond) Y1 Y2
clorof. Metanol
0 1 1.3298 57 1.3299 1.45227563 -0.45227563
0.1 0.9 1.3485 54 1.3670 1.4522677 -0.4522677
0.2 0.8 1.3666 50 1.3927 1.45221472 -0.45221472
0.3 0.7 1.3790 49 1.4049 1.45215271 -0.45215271
0.4 0.6 1.3924 47 1.4186 1.45206219 -0.45206219
0.5 0.5 1.4072 46.5 1.4191 1.45193383 -0.45193383
0.6 0.4 1.4113 47 1.4212 1.45189301 -0.45189301
0.7 0.3 1.4257 47 1.4257 1.4517315 -0.4517315
0.8 0.2 1.4314 48 1.4257 1.45165977 -0.45165977
0.9 0.1 1.4419 50 1.4434 1.45151609 -0.45151609
1 0 1.4439 54 1.440 1.45148702 -0.45148702
Grafica x y n liq
nliq vs N1
1.46
1.44
1.42
y = -0.068x2 + 0.1817x + 1.3309
1.4
R² = 0.9974
nliq
1.38
1.36
1.34
1.32
0 0.2 0.4 0.6 0.8 1 1.2
X1
Y= -0.068X2+0.1817X+1.3309 R2=0.9974
Para calcular y1:
𝑷𝟏 𝒔𝒂𝒕
𝒚𝟏 = 𝒙𝟏
𝑷
𝑌1 𝑒𝑥𝑝 = −0.068(1.3298)2 + 0.1817(1.3298) + 1.3309 = 1.4522756
𝑌2 𝑒𝑥𝑝 = −0.068(1.3485)2 + 0.1817(1.3385) + 1.3309 = 1.4522677
𝑌3 𝑒𝑥𝑝 = −0.068(1.3666)2 + 0.1817(1.3666) + 1.3309 = 1.4522147
𝑌4 𝑒𝑥𝑝 = −0.068(1.3790)2 + 0.1817(1.3790) + 1.3309 = 1.45215271
𝑌5 𝑒𝑥𝑝 = −0.068(1.3924)2 + 0.1817(1.3924) + 1.3309 = 1.45206219
𝑌6 𝑒𝑥𝑝 = −0.068(1.4072)2 + 0.1817(1.4072) + 1.3309 = 1.45193383
𝑌7 𝑒𝑥𝑝 = −0.068(1.4113)2 + 0.1817(1.4113) + 1.3309 = 1.45189301
𝑌8 𝑒𝑥𝑝 = −0.068(1.4257)2 + 0.1817(1.4257) + 1.3309 = 1.4517315
𝑌9 𝑒𝑥𝑝 = −0.068(1.4314)2 + 0.1817(1.4314) + 1.3309 = 1.45165977
𝑌10 𝑒𝑥𝑝 = −0.068(1.4419)2 + 0.1817(1.4419) + 1.3309 = 1.45151609
𝑌11 𝑒𝑥𝑝 = −0.068(1.4439)2 + 0.1817(1.4439) + 1.3309 = 1.45148702
Para y2
𝒚𝟐 = 𝟏 − 𝒚𝟏
y1 = 1 − 1.4522756 = −0.4522756
y2 = 1 − 1.4522677 = −0.4522677
y3 = 1 − 1.4522147 = − 0.4522147
y4 = 1 − 1.45215271 = −0.45215271
y5 = 1 − 1.45206219 = −0.45206219
y6 = 1 − 1.45193383 = −0.45193383
y7 = 1 − 1.45189301 = −0.45189301
y8 = 1 − 1.4517315 = −0.4517315
y9 = 1 − 1.45165977 = −0.45165977
y10 = 1 − 1.45151609 = −0.45151609
y11 = 1 − 1.45148702 = −0.45148702
Calcular las presiones de saturación para cada componente, aplicando la ecuación
de Antoine a la temperatura seleccionada (T experimental)
A B C
Clorofomo 15.9731 26.9679 -46.16
Metanol 18.5875 3626.55 -34.29
𝐵
𝐿𝑛𝑃𝑠𝑎𝑡 = 𝐴 − 𝑇+𝑐
𝑩
[𝑨−𝑻+𝑪]
𝑷𝟏 𝒔𝒂𝒕 = 𝒆
26.9679
P1 sat = e[15.9731−330.15−46.16]=7866621.674
26.9679
P2 sat = e[15.9731−327.15−46.16]=7858650.13
26.9679
P3 sat = e[15.9731−323.15−46.16]=7847765.86
26.9679
P4 sat = e[15.9731−322.15−46.16]=7844997.904
26.9679
P5 sat = e[15.9731−320.15−46.16]=7839404.36
26.9679
P6 sat = e[15.9731−319.65−46.16]=7837993.82
26.9679
P7 sat = e[15.9731−320.15−46.16]=7839404.36
26.9679
P8 sat = e[15.9731−320.15−46.16]=7839404.361
26.9679
P9 sat = e[15.9731−321.15−46.16]=7842210.8
26.9679
P10 sat = e[15.9731−323.15−46.16]=7847765.86
26.9679
P11 sat = e[15.9731−327.15−46.16]=7858650.23
𝑩
𝒔𝒂𝒕 [𝑨−𝑻+𝑪]
𝑷𝟐 =𝒆
3626.55
P1 sat = e[18.5875−330.15−34.29]=561.069258
3626.55
P2 sat = e[18.5875−327.15−34.29]=494.8622042
3626.55
P3 sat = e[18.5875−323.15−34.29]=416.881589
3626.55
P4 sat = e[18.5875−322.15−34.29]=399.0905449
3626.55
P5 sat = e[18.5875−320.15−34.29]=365.419165
3626.55
P6 sat = e[18.5875−319.65−34.29]=357.385929
3626.55
P7 sat = e[18.5875−320.15−34.29]=365.419165
3626.55
P8 sat = e[18.5875−320.15−34.29]=365.419165
3626.55
P9 sat = e[18.5875−321.15−34.29]=381.942602
3626.55
P10 sat = e[18.5875−323.15−34.29]=416.881589
3626.55
P11 sat = e[18.5875−327.15−34.29]=494.862204
X1 T(K) Psat 1 Psat 2
0 330.15 7866621.7 561.06926
0.1 327.15 7858650.1 494.8622
0.2 323.15 7847765.9 416.88159
0.3 322.15 7844997.9 399.09054
0.4 320.15 7839404.4 365.41917
0.5 319.65 7837993.8 357.3859
0.6 320.15 7839404.4 365.4191
0.7 320.15 7839404.4 365.4191
0.8 321.15 7842210.8 381.9426
0.9 323.15 7847765.9 416.8815
1 327.15 7858650.1 494.8622
Para calcular los coeficientes de actividad, se emplea la ley de Raoult modificada
P=585mmHg
𝒆𝒙𝒑 𝒆𝒙𝒑
𝒆𝒙𝒑 𝒚𝟏 𝑷 𝒆𝒙𝒑 𝒚𝟐 𝑷
ɣ𝟏 = 𝒔𝒂𝒕−𝒆𝒙𝒑 ɣ𝟐 = 𝒔𝒂𝒕−𝒆𝒙𝒑
𝒙𝟏 𝑷 𝟏 𝒙𝟐 𝑷 𝟐
Para ɣ1
exp (1.45227563)(585)
ɣ1 = =0
(0)(7866621.674)
exp (1.4522677)(585)
ɣ2 = = 1.08x10−3
(0.1)(7858650.128)
exp (1.45221472)(585)
ɣ3 = = 5.41x10−4
(0.2)(7847765.86)
exp (1.45215271)(585)
ɣ4 = = 3.60x10−4
(0.3)(7844997.905)
exp (1.45206219)(585)
ɣ5 = = 2.70x10−4
(0.4)(7839404.361)
exp (1.45193383)(585)
ɣ6 = = 2.16x10−4
(0.5)(7837993.823)
exp (1.45189301)(585)
ɣ7 = = 1.80x10−4
(0.6)(7839404.361)
exp (1.4517315)(585)
ɣ8 = = 1.54x10−4
(0.7)(7839404.361)
exp (1.45165977)(585)
ɣ9 = = 1.35x10−4
(0.8)(7842210.805)
exp (1.45151609)(585)
ɣ10 = = 1.20x10−4
(0.9)(7847765.86)
exp (1.45148702)(585)
ɣ11 = = 1.08x10−4
(1)(7858650.128)
Para ɣ2
(0.4522756)(585)
ɣ1𝑒𝑥𝑝 = = 0.4715661
(1)(561.06926)
(0. 4522677)(585)
ɣ𝑒𝑥𝑝
2 = = 0.5940522
(0.9)(494.8622)
(0. 4522147)(585)
ɣ𝑒𝑥𝑝
3 = = 0.7932277
(0.8)(416.88159)
(0. 4521527)(585)
ɣ𝑒𝑥𝑝
4 = = 0.9468289
(0.7)(399.09054)
(0. 4520622)(585)
ɣ𝑒𝑥𝑝
5 = = 1.2061782
(0.6)(365.41917)
(0. 4519338)(585)
ɣ𝑒𝑥𝑝
6 = = 1.4795283
(0.5)(357.38593)
(0. 451893)(585)
ɣ𝑒𝑥𝑝
7 = = 1.8085902
(0.4)(365.41917)
(0. 4517315)(585)
ɣ𝑒𝑥𝑝
8 = = 2.4105917
(0.3)(365.41917)
(0. 4516598)(585)
ɣ𝑒𝑥𝑝
9 = = 3.4589094
(0.2)(381.9426)
𝑒𝑥𝑝 (0. 4515561)(585)
ɣ10 = = 6.3360177
(0.1)(416.88159)
𝑒𝑥𝑝 (0. 451487)(585)
ɣ11 = =0
(0)(494.8622)
Se procede a graficar ln 1 x1 t ln 2 y x2 y determinar las constantes de Margules y
Van Laar
𝒆𝒙𝒑 𝒆𝒙𝒑
Cálculo de 𝐥𝐧 ɣ𝟏 y 𝐥𝐧 ɣ𝟐
𝒆𝒙𝒑
Para 𝐥𝐧 ɣ𝟏
𝒆𝒙𝒑
lnɣ𝟏 (0) = 0
𝒆𝒙𝒑
lnɣ𝟐 (1.08x10−3 ) = −6.8298022
𝒆𝒙𝒑
lnɣ𝟑 (5.41x10−4 ) = −7.521599
𝒆𝒙𝒑
lnɣ𝟒 (3.60x10−4 ) = − 7.9267549
𝒆𝒙𝒑
lnɣ𝟓 (2.70x10−4 ) = −8.2137861
𝒆𝒙𝒑
lnɣ𝟔 (2.16x10−4 ) = −8.4368381
𝒆𝒙𝒑
lnɣ𝟕 (1.80x10−4 ) = −8.6193677
𝒆𝒙𝒑
lnɣ𝟖 (1.54x10−4 ) = −8.7736296
𝒆𝒙𝒑
lnɣ𝟗 (1.35x10−4 ) = −8.9075684
𝒆𝒙𝒑
lnɣ𝟏𝟎 (1.20x10−4 ) = −9.0261585
𝒆𝒙𝒑
lnɣ𝟏𝟏 (1.08x10−4 ) = − 9.132925
𝒆𝒙𝒑
Para 𝐥𝐧 ɣ𝟐
𝒆𝒙𝒑
lnɣ𝟏 (0.4715661) = −0.75169598
𝒆𝒙𝒑
lnɣ𝟐 (0.5940522) = −0.52078801
𝒆𝒙𝒑
lnɣ𝟑 (0.7932277) = −0.23164499
𝒆𝒙𝒑
lnɣ𝟒 (0.9468289) = −0.05463683
𝒆𝒙𝒑
lnɣ𝟓 (1.2061782) = 0.18745686
𝒆𝒙𝒑
lnɣ𝟔 (1.4795283) = 0.3917233
𝒆𝒙𝒑
lnɣ𝟕 (1.8085902) = 0.59254765
𝒆𝒙𝒑
lnɣ𝟖 (2.4105917) = 0.87987225
𝒆𝒙𝒑
lnɣ𝟗 (3.4589094) = 1.24095333
𝒆𝒙𝒑
lnɣ𝟏𝟎 (6.3360177) = 1.84625044
𝒆𝒙𝒑
lnɣ𝟗 (0) = 1
Modelo de Van Laar (𝑨𝟏𝟐 == −6.8298022 𝑨𝟐𝟏 = 1.84625044)
𝐴12
[ ]
𝐴12 ∗𝑥1 2
[1+ ]
Sacando 𝜸𝟏𝒄𝒂𝒍 𝛾1𝑐𝑎𝑙 = 𝑒 𝐴21 ∗𝑥2
−6.8298022
[ 2]
−6.8298022∗0
[1+ ]
𝛾1𝑐𝑎𝑙1 = 𝑒 1.84625044∗1 = 0.0010811
−6.8298022
[ 2]
−6.8298022∗0.1
[1+ ]
𝛾1𝑐𝑎𝑙2 = 𝑒 1.84625044∗0.9 = 0.002901
−6.8298022
[ ]
−6.8298022∗0.2 2
[1+ ]
𝛾1𝑐𝑎𝑙3 = 𝑒 1.84625044∗0.8 = 0.0251914
−6.8298022
[ ]
−6.8298022∗0.3 2
[1+ ]
𝛾1𝑐𝑎𝑙4 = 𝑒 1.84625044∗0.7 = 0.143149
−6.8298022
[ ]
−6.8298022∗0.4 2
[1+ ]
𝛾1𝑐𝑎𝑙5 = 𝑒 1.84625044∗0.6 = 0.38122
−6.8298022
[ 2]
−6.8298022∗0.5
[1+ ]
𝛾1𝑐𝑎𝑙6 = 𝑒 1.84625044∗0.5 = 0.628074
−6.8298022
[ ]
−6.8298022∗0.6 2
[1+ ]
𝛾1𝑐𝑎𝑙7 = 𝑒 1.84625044∗0.4 = 0.806672
−6.8298022
[ ]
−6.8298022∗0.7 2
[1+ ]
𝛾1𝑐𝑎𝑙8 = 𝑒 1.84625044∗0.3 = 0.913516
−6.8298022
[ 2]
−6.8298022∗0.8
[1+ ]
𝛾1𝑐𝑎𝑙9 = 𝑒 1.846250444∗0.2 = 0.964426
−6.8298022
[ 2]
−6.8298022∗0.9
[1+ ]
𝛾1𝑐𝑎𝑙10 = 𝑒 1.846250444∗0.1 = 0.993862
−6.8298022
[ 2]
−6.8298022∗1
[1+ ]
𝛾1𝑐𝑎𝑙11 = 𝑒 1.846250444∗0 =0
𝐴21
[ ]
𝐴 ∗𝑥 2
[1+ 21 2 ]
Sacando 𝜸𝟐𝒄𝒂𝒍 𝛾2𝑐𝑎𝑙 = 𝑒 𝐴12 ∗𝑥1
1.846250444
[ ]
1.846250444∗1 2
[1+ ]
𝛾2𝑐𝑎𝑙1 = 𝑒 −6.8298022∗0 = 0
1.846250444
[ ]
1.846250444∗0.9 2
[1+ ]
𝛾2𝑐𝑎𝑙2 = 𝑒 −6.8298022∗0.1 = 2.457636
1.846250444
[ ]
1.846250444∗0.8 2
[1+ ]
𝛾2𝑐𝑎𝑙3 = 𝑒 −6.8298022∗0.2 = 2.73𝑥1038
1.846250444
[ ]
1.846250444∗0.7 2
[1+ ]
𝛾2𝑐𝑎𝑙4 = 𝑒 −6.8298022∗0.3 = 760080.07
1.846250444
[ ]
1.846250444∗0.6 2
[1+ ]
𝛾2𝑐𝑎𝑙5 = 𝑒 −6.8298022∗0.4 = 185.5869
1.846250444
[ ]
1.846250444∗0.5 2
[1+ ]
𝛾2𝑐𝑎𝑙6 = 𝑒 −6.8298022∗0.5 = 32.0596
1.846250444
[ ]
1.846250444∗0.4 2
[1+ ]
𝛾2𝑐𝑎𝑙7 = 𝑒 −6.8298022∗0.6 = 15.5989
1.846250444
[ ]
1.846250444∗0.3 2
[1+ ]
𝛾2𝑐𝑎𝑙8 = 𝑒 −6.8298022∗0.7 = 10.6099
1.846250444
[ ]
1.846250444∗0.2 2
[1+ ]
𝛾2𝑐𝑎𝑙9 = 𝑒 −6.8298022∗0.8 = 8.3609
1.846250444
[ ]
1.846250444∗0.1 2
[1+ ]
𝛾2𝑐𝑎𝑙10 = 𝑒 −6.8298022∗0.9 = 7.1161
1.846250444
[ ]
1.846250444∗0 2
[1+ ]
𝛾2𝑐𝑎𝑙11 = 𝑒 −6.8298022∗1 = 6.336
Calcular las temperaturas de ebullición de los componentes puros a la presión del
sistema 𝑙𝑛 𝑃 𝑠𝑎𝑡 = 𝐴 − 𝐵 𝑇 + 𝐶 𝑑𝑜𝑛𝑑𝑒 𝑇[=]𝐾, 𝑃[=]𝑚𝑚𝐻𝑔 𝑇1 = ( 𝐵1 𝐴1 − 𝑙𝑛( 𝑃) − 𝐶1)
; 𝑇2 = ( 𝐵2 𝐴2 − 𝑙𝑛( 𝑃) − 𝐶2)
Temperatura supuesta en base a los componentes puros, para cada composición
𝑇𝑠𝑢𝑝 1 = (𝑥1 ∗ 𝑇1) + (𝑥2 ∗ 𝑇2)
P=mmHg
585 = 𝛾1 𝑐𝑎𝑙 ∗ 𝑥1 ∗ 𝑃1𝑠𝑎𝑡 + 𝛾2 𝑐𝑎𝑙 ∗ 𝑥2 ∗ 𝑃2𝑠𝑎𝑡 ∗ 𝑇
585 − 𝛾1 𝑐𝑎𝑙 ∗ 𝑥1 ∗ 𝑃1𝑠𝑎𝑡
𝑇=
𝛾2 𝑐𝑎𝑙 ∗ 𝑥2 ∗ 𝑃2𝑠𝑎𝑡
585 − (0.00108107)(0)( 7866621.7)
𝑇1 = =0
(0)(1)( 561.06926)
585 − (0.002901)(0.1)( 7858650.1)
𝑇2 = = −1.5483
(2.457636 )(0.9)( 494.8622)
585 − (0.0251914)(0.2)( 78747765.9)
𝑇3 = = −4.35𝑥10−36
(2.73𝑥1038 )(0.8)( 416.88159)
585 − (0.143149)(0.3)( 7844997.9)
𝑇4 = = −0.001583
(760080.072)(0.7)( 399.09054)
585 − (0.38122)(0.4)( 7839404.4)
𝑇5 = = −29.36407
(185.5869)(0.6)( 365.41917)
585 − (0.628074)(0.5)( 7837993.8)
𝑇6 = = −429.5532
(32.0596)(0.5)( 357.38593)
585 − (0.806672)(0.6)( 7839404.4)
𝑇7 = = −1663.869
(15.5989)(0.4)( 365.41917)
585 − (0.913516)(0.7)( 7839404.4)
𝑇8 = = −4309.459
(10.6099)(0.3)( 365.41917)
585 − (0.964426)(0.8)( 7842210.8)
𝑇9 = = −9472.711
(8.3609)(0.2)( 381.9426)
585 − (0.993862)(0.9)( 784765.9)
𝑇10 = = −23660.49
(7.1161)(0.1)( 416.88159)
585 − (0)(1)( 7858650.1)
𝑇11 = =0
(6.336)(0)( 494.8622)
Temperatura de burbuja
𝛾1 𝑐𝑎𝑙 𝑥1 𝑃10 𝛾2 𝑐𝑎𝑙 𝑥2 𝑃20
𝑦1 𝑐𝑎𝑙 = 𝑦2 =
𝑃 𝑃
Para 𝒚𝟏 𝒄𝒂𝒍
0.00108107 ∗ 0 ∗ 7866621.7
𝑦1 𝑐𝑎𝑙1 = =0
585
0.002901 ∗ 0.1 ∗ 7858650.1
𝑦1 𝑐𝑎𝑙2 = = 3.8970844
585
0.0251914 ∗ 0.2 ∗ 7847765.9
𝑦1 𝑐𝑎𝑙3 = = 67.588447
585
0.143149 ∗ 0.3 ∗ 7844997.9
𝑦1 𝑐𝑎𝑙4 = = 575.89928
585
0.38122 ∗ 0.4 ∗ 7839404.36
𝑦1 𝑐𝑎𝑙5 = = 2043.4446
585
0.628074 ∗ 0.5 ∗ 7837993.823
𝑦1 𝑐𝑎𝑙6 = = 4207.5557
585
0.806672 ∗ 0.6 ∗ 7839404.361
𝑦1 𝑐𝑎𝑙7 = = 6485.9774
585
0.91351 ∗ 0.7 ∗ 7839404.361
𝑦1 𝑐𝑎𝑙8 = = 8569.2221
585
0.964426 ∗ 0.8 ∗ 7842210.805
𝑦1 𝑐𝑎𝑙9 = = 10342.881
585
0.993862 ∗ 0.9 ∗ 7847765.86
𝑦1 𝑐𝑎𝑙10 = = 11999.379
585
0 ∗ 1 ∗ 7857650.128
𝑦1 𝑐𝑎𝑙11 = =0
585
Para 𝑦2 𝑐𝑎𝑙
0 ∗ 1 ∗ 561.06926
𝑦2 𝑐𝑎𝑙1 = =0
585
2.457636 ∗ 0.9 ∗ 494.8622042
𝑦2 𝑐𝑎𝑙2 = = 1.87106334
585
𝑐𝑎𝑙3
2.73𝑥1038 ∗ 0.8 ∗ 416.8815891
𝑦2 = = 1.5544𝑥1038
585
760080.072 ∗ 0.7 ∗ 399.0905449
𝑦2 𝑐𝑎𝑙4 = = 362971.862
585
185.5869 ∗ 0.6 ∗ 365.419165
𝑦2 𝑐𝑎𝑙5 = = 69.5559077
585
32.0596 ∗ 0.5 ∗ 357.3859292
𝑦2 𝑐𝑎𝑙6 = = 9.79286319
585
15.5989 ∗ 0.4 ∗ 365.419165
𝑦2 𝑐𝑎𝑙7 = = 3.89752958
585
10.6099 ∗ 0.3 ∗ 365.419165
𝑦2 𝑐𝑎𝑙8 = = 1.98823631
585
8.3609 ∗ 0.2 ∗ 381.942602
𝑦2 𝑐𝑎𝑙9 = = 1.09175518
585
7.116 ∗ 0.1 ∗ 416.881589
𝑦2 𝑐𝑎𝑙10 = = 0.50710617
585
6.336 ∗ 0 ∗ 494.862204
𝑦2 𝑐𝑎𝑙11 = =0
585
Gráficas de datos Experimentales.
Título del gráfico
1.4525
1.4524
1.4523
1.4522
1.4521
1.452
1.4519
1.4518
1.4517
1.4516
1.4515
1.4514
0 0.2 0.4 0.6 0.8 1 1.2
nliq vs X1
1.2
1 y = 9.1751x - 12.339
0.8 R² = 0.9741
x1
0.6
0.4
0.2
0
1.34 1.36 1.38 1.4 1.42 1.44 1.46
nliq
0 Lineal (0) Polinómica (0)
Observaciones:
La experimentación en laboratorio comenzó con una serie de soluciones de
cloroformo y metanol, se determinó la temperatura de burbuja, es decir, se registró
el punto de temperatura en el cual la mezcla comenzó a sacar burbujas (ebullir).
Dado que se trata de un sistema de equilibrio - vapor en el que las sustancias
involucradas (cloroformo-metanol) son inmiscibles, se aplican las ecuaciones de
Van Laar o Margules, para llevar acabo un ajuste en los datos obtenidos y poder
graficar el diagrama de la fase liquido-vapor.
A través de los diagramas de equilibrio liquido-vapor se comprueba el cambio de
fase que existe en la solución de cloroformo-metanol, así como el punto en el que
ambas soluciones coexisten en equilibrio. Por lo que se puede decir que es una
binaria, en las cuales se presenta el punto azeotrópico.
Con los diagramas también es posible observar que cuando la presión es mayor
que la presión de vapor corresponde al líquido como fase estable y los puntos por
debajo de la línea cuando la presión es menor que la presión de vapor, de forma
que la muestra se evapora entonces corresponden al vapor como fase estable.
Cuando el vapor y el líquido están en equilibrio, sus composiciones no son las
mismas.
Finalmente decir que la experimentación no pudo ser llevada a cabo debido a
problemas con los materiales, pero se trató de dar una descripción lo más exacta
posible y la maestra procuró darnos a entender la manera de interpretar los datos.
Conclusiones:
Durante la experimentación se observa que al iniciar la experimentación se puede
considerar que dos o más fases se presentan cuando hay una superficie
macroscópica de separación entre las fases homogéneas, de manera tal que el
número de fases que pueden identificar depende de cómo se observe la muestra.
Se comprobó que en una mezcla homogénea de dos líquidos volátiles A y B se
pueden presentar dos casos, esto mediante la ecuación de Raoult, en el primero si
las fuerzas intermoleculares entre las moléculas A y B son más débiles que las que
existen entre las moléculas de cada componente se dice que existe una mayor
tendencia para considerarse real.
En el segundo caso, si las moléculas A y B se atraen con mayor fuerza que entre
las del mismo tipo la presión de vapor es menor que la suma de las presiones
parciales y se conoce como desviación negativa.
De igual forma se realizaron los cálculos de γ(1^cal ) y γ(2^cal ) por medio de la
ecuación de Van Laar donde se observa que en ambos casos aumenta
gradualmente y con estos datos realizar el cálculo de las temperaturas.
Por último, se concluye que la práctica se llevó a cabo de manera correcta, al
verificar los cálculos y las gráficas se observa de manera más clara el
comportamiento de las dos soluciones.
