Obtención de la curva de equilibrio x-y de la mezcla etanol (1) / agua (2)

Castillo Ludeña, J.J. (20154119B); López Alejos, J.J. (20150238G); Silva De la Vega, R.M. (20150411K).
Escuela de Ingeniería Petroquímica, FIP-UNI.
Laboratorio de Operaciones Unitarias II, 23/04/2019
mcastillol@fip.uni.edu.pe; jlopeza@fip.uni.edu.pe; rsilvad@fip.uni.edu.pe

Resumen: El presente laboratorio se lleva a cabo cuando las presiones de vapor de los componentes no son
muy parecidas; es por ello, que se busca el aprendizaje del buen manejo del equipo y, además, realizar el
adecuado balance de materia para la obtención de los resultados adecuados. Se utiliza la destilación de una
mezcla en el cuál uno de los componentes sea más volátil que el otro, teniendo como resultados valores
experimentales de la fracción de vapor experimental del etanol y las temperaturas de burbuja concluyendo así
que según la Ley de Raoult Modificada se obtendrán datos más cercanos a los reales (experimentales).

Abstract: The present laboratory is carried out when the vapor pressures of the components are not very
similar; it is for this reason, that the learning of the good handling of the equipment is looked for and, in
addition, to realize the suitable balance of matter for the obtaining of the suitable results. The distillation of a
mixture in which one of the components is more volatile than the other is used, resulting in experimental
values of the experimental steam fraction of the ethanol and the bubble temperatures thus concluding that
according to the Modified Raoult Law they will be obtained data closer to the real ones (experimental).

1. INTRODUCCIÓN DF, unen las mezclas en el equilibrio en D y F. Hay

Con el nombre de destilación se entiende la un número infinito de dichas líneas de unión para
separación de los componentes de una mezcla este diagrama. Una mezcla en la curva inferior,
líquida al estado de sustancias puras por como en el punto D, es un líquido saturado; una
vaporización parcial de la misma, de tal manera que mezcla en la curva superior, como en F, es un vapor
la composición del vapor obtenido sea distinta de la saturado. Una mezcla en E es una mezcla de dos
composición del líquido de partida, resultando fases, que consta de una fase líquida de
distinta también la composición del líquido residual composición en D y una fase vapor de composición
[1]. en F, en proporción tal que la composición
promedio de toda la mezcla se representa mediante
Equilibrio Líquido - Vapor: E [2].
Inicialmente, se van a considerar las mezclas
binarias, las cuales serán llamadas “ordinarias”, para
indicar (a) que los componentes líquidos se
disuelven en cualquier proporción para formar
soluciones homogéneas, las cuales no son
necesariamente ideales, y (b) que no hay
complicación alguna de puntos de ebullición
máximos o mínimos. El componente A de la mezcla
binaria A-B se considerará como el más volátil; esto
significa que la presión de vapor de A puro a
cualquier temperatura es mayor que la presión de
vapor B puro [2].

Equilibrio a presión constante: Considérese

primero una sección típica a presión constante
(Fig.1.) la cual la curva superior proporciona la
relación entre la temperatura y la composición del
vapor, mientras que la curva inferior proporciona la
relación entre la temperatura y la composición del
líquido. Las mezclas de líquido y vapor en el
equilibrio están a la misma temperatura y presión, Fig.1. Temperatura vs Fracción mol de A.
de forma que las líneas de unión, como la línea
Las composiciones vapor-liquido en el equilibrio desaparece en este sistema a presiones menores de
pueden mostrarse también sobre un diagrama de 70 mmHg.
composición (x vs y*), como en la Fig.2. Por
ejemplo, el punto P sobre el diagrama representa la
línea de unión DF [2].
Para realizar el balance de masa en el destilador se
tiene en la zona de enriquecimiento:

Fig.2. Diagrama de composición (x vs y*)

Ley de Raoult: Para mezclas ideales [1].

Fig.3. Zona de enriquecimiento

Ec…1
Balance Total: G n+ 1=Ln+ D
Volatilidad relativa: Una medida numérica de lo Balance del componente más volátil:
anterior se conoce como el factor de separación o, y n+ 1 Gn+1 =Ln x n + D x D
particularmente en el caso de la destilación, la Condición de reflujo total:
volatilidad relativa α. Esta es la relación entre la L
R=∞ =1=cte
relación de concentraciones de A y B en una fase y G
en la otra, y es una medida de la posibilidad de Además, x D =x 0
separación [2]. Principio de flujo equimolar y vaporización, siendo
r y n dos platos consecutivos:
Lr−1=Ln y Gn+1=Gr

Ec…2 L
Además, R=  L=RD
Mezclas azeotrópicas: Cuando las desviaciones D
positivas del ideal son lo suficientemente grandes y Por lo tanto, en el balance del componente más
cuando las presiones de vapor de los dos volátil y del balance total se tiene:
componentes no están muy alejadas entre sí, las y n+ 1 G=L x n+ D x D
curvas de presión total a temperatura constante, G=L+ D
pueden aumentar a través de un máximo a cierta Reemplazando
concentración. Una de las más importantes es el R 1
azeótropo etanol-agua, el cual a 1 atm aparece a y n+ 1= xn + x
R+1 R+1 D
89.4% en mol de etanol y 78.2 ºC. El azeotropismo Para condición de reflujo total se tendrá que:
 y n+ 1=x n

1.1 Objetivos:
 Obtención del diagrama de McCabe –
Thiele para una mezcla de etanol (1) / agua
(2).
 Comparación de los valores experimentales
obtenidos con los valores teóricos según la
Ley de Raoult, Ley de Raoult Modificada Fig.4. Equipo de destilación
(Van Laar) y según la volatilidad.
 Comprender el manejo de un equipo de Los elementos necesarios para la realización de esta
destilación. práctica son:
 Equipo UDCC o UDDC.
2. METODOS Y MATERIALES  Solución de etanol y agua destilada al 1:1
La unidad de destilación avanzada (UDCC/UDDC) molar.
está compuesta por un calderín sobre el que pueden Para el desarrollo de esta práctica siga el siguiente
adaptarse dos tipos de columna intercambiables (de procedimiento:
platos y de relleno), un sistema de reflujo, un I. Introduzca la solución problema en el
depósito para la recepción del destilado, una bomba calderín a través de la boca apropiada
de vacío y una bomba para efectuar la alimentación utilizando un embudo para facilitar el
en continuo. El control se realiza con un PC y sus proceso.
correspondientes periféricos, que permiten el II. Encienda el circuito de refrigeración.
registro de los datos. El programa de control es III. Ejecute el programa SACED-UDDC y
sencillo e intuitivo. seleccione un nombre para el archivo de
La columna de 8 platos tipo Brunn presenta tomas datos.
de muestras y de temperatura en cada uno de los IV. Encienda la manta calefactora y programe
platos. como punto de consigna un valor superior
El vapor que llega a la cabeza de columna es a la temperatura de burbuja e inferior a la
enviado a un condensador total. El caudal de agua temperatura de rocío de la mezcla a la
de refrigeración que atraviesa el condensador se presión de operación (1 atm). Para esto,
regula y se indica en un medidor de flujo. La caída medirá el índice de refracción de la
de presión en la columna puede ser medida con un muestra en el calderín y en base a la curva
sensor de presión diferencial. de calibración determinará la fracción
La columna puede trabajar en discontinuo (UDDC) molar de la muestra. Con la curva de
o en continuo (UDCC). Para efectuar la equilibrio (anexo 5.1) determinará las
alimentación en continuo se dispone de una bomba temperaturas de burbuja y rocío para esa
que puede inyectar el alimento directamente en el composición.
calderín o en cualquiera de los platos. La V. Sitúe la válvula de reflujo en la posición en
temperatura de alimentación es regulable mediante que todo el vapor condensado vuelve a la
una resistencia con control PID. columna.
Las temperaturas del sistema se miden mediante VI. Permita que las temperaturas en cada plato
termopares situados en posiciones estratégicas y se se estabilicen.
visualizan en el computador. VII. A un tiempo determinado en el cronómetro
tome una muestra en el calderín y en cada
uno de los platos lo más rápidamente
posible (idealmente, las muestras deben
tomarse simultáneamente). Tape, etiquete y
enfríe rápidamente la muestra en cada uno
de los tubos. Analice las muestras mediante
el refractómetro.
VIII. Registre el tiempo en el cronómetro y el 0.6414 0.7209 0.8034 0.7290 0.8840
tiempo indicado en el programa SACED- 1.0000 1.0000 1.0000 1.0000 1.0000
UDDC (este tiempo servirá para
determinar el registro de temperaturas en
cada plato). Además, registre el número de
refractómetro, la temperatura de medición
y el índice de refracción de cada muestra.
Con el índice de refracción corregido a
25°C obtenga la fracción molar de cada
muestra.
IX. Repita estas mediciones cada 15 minutos
hasta completar tres series de resultados.

3. RESULTADOS Y DISCUSION Tabla 2. Comparación de los valores de Tb

experimental con los teóricos obtenidos
3.1 Resultados según Ley de Raoult, Ley de Raoult
Modificada (Van Laar) y según la volatilidad
Tablas para la 1ra corrida de datos a P constante.
Tb, exp Tb, T (ºC)
x1, Exp: Fracción líquida experimental del etanol x1
(ºC) LR LRM α
[adimensional]. 0.0000 100.00 100.08 100.10 N.P.
y1, exp: Fracción vapor experimental del etanol 0.0389 99.06 98.75 92.70 115.83
[adimensional]. 0.0740 90.92 97.61 89.00 112.79
y1, T: Fracción vapor teórica del etanol 0.0764 89.45 97.53 88.80 112.60
[adimensional]. 0.1879 85.71 94.22 83.70 104.83
Tb, exp: Temperatura de burbuja experimental [ºC]. 0.3193 79.78 90.81 81.40 97.98
Tb, T: Temperatura de burbuja teórica [ºC].
0.4207 82.25 88.46 80.40 93.75
LR: Método del cálculo de y1, T según la Ley de
0.4611 81.44 87.59 80.00 92.26
Raoult.
0.6414 81.29 84.03 78.80 86.55
LRM: Método del cálculo de y1, T según la Ley de
Raoult Modificada. 1.0000 78.30 78.26 78.30 78.26
α: Método del cálculo de y1, T según la volatilidad
Tabla 3. Comparación de los valores de y1
Tabla 1. Comparación de los valores de y1 experimental con los teóricos obtenidos según
experimental con los teóricos obtenidos según Ley de Raoult, Ley de Raoult Modificada (Van
Ley de Raoult, Ley de Raoult Modificada (Van Laar) y según la volatilidad para la 2da corrida de
Laar) y según la volatilidad para la 1ra corrida de datos a P constante.
datos a P constante. y1, T
x1 y1, exp
y1, T LR LRM α
x1 y1, exp 0.0000 0.0000 0.0000 0.0000 0.0000
LR LRM α
0.0000 0.0000 0.0000 0.0000 0.0000 0.0197 0.0188 0.0432 0.1613 0.0789
0.0389 0.0740 0.0834 0.2620 0.1471 0.0188 0.0420 0.0413 0.1557 0.0756
0.0740 0.0764 0.1525 0.3760 0.2539 0.0420 0.0168 0.0898 0.2753 0.1574
0.0764 0.1879 0.1570 0.3820 0.2605 0.0168 0.2484 0.0370 0.1422 0.0679
0.1879 0.3193 0.3433 0.5290 0.4965 0.2484 0.3711 0.4279 0.5666 0.5847
0.3193 0.4207 0.5154 0.5990 0.6665 0.3711 0.3770 0.5726 0.6191 0.7154
0.4207 0.4611 0.6228 0.6370 0.7558 0.3770 0.4240 0.5788 0.6213 0.7206
0.4611 0.6414 0.6606 0.6520 0.7847 0.4240 0.4306 0.6259 0.6386 0.7583
 1.0000 1.0000 1.0000 1.0000 1.0000
TXY a P cte
Experimen-
tal XT
112 Raoult TX
Van Laar
107 TX
Volatilidad
102 TX
Experimen-
97 tal YT

T (ºC)
Raoult TY
Tabla 4. Comparación de los valores de Tb 92
Van Laar
experimental con los teóricos obtenidos TY
87
según Ley de Raoult, Ley de Raoult
Modificada (Van Laar) y según la volatilidad 82
para la 2da corrida de datos a P constante.
Tb, exp Tb, T (ºC) 77
x1
(ºC) LR LRM α 72
0.0000 100.00 100.08 100.10 N.P. 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
0.0197 99.99 99.40 95.73 117.65 x, y
0.0188 98.19 99.43 95.89 117.73 Gráfico 2. Comparación de las gráficas T-X-Y a P
0.0420 97.90 98.65 92.32 115.54 constante obtenidos según Ley de Raoult, Ley de
0.0168 93.56 99.50 96.27 117.94 Raoult Modificada (Van Laar) y según la volatilidad
0.2484 83.99 92.59 82.40 101.44 para la 1ra corrida de datos.
0.3711 79.88 89.58 80.83 95.72
0.3770 77.03 89.44 80.77 95.48
0.4240 77.94 88.39 80.34 93.63 XY a P cte
1.0000 78.30 78.26 78.30 78.26
1
Gráficos 0.8 Experimental
Raoult
0.6
Van Laar
XY a P cte
y

0.4 Volatilidad
1
45°
0.9 0.2
0.8
0.7 Exper- 0
0.6 imental 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
0.5 Raoult x
y

0.4
0.3 Van Laar Gráfico 3. Comparación de las gráficas XY a P
0.2 Volatilidad constante obtenidos según Ley de Raoult, Ley de
0.1 45° Raoult Modificada (Van Laar) y según la volatilidad
0 para la 2da corrida de datos.
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x
Gráfico 1. Comparación de las gráficas XY a P
constante obtenidos según Ley de Raoult, Ley de
Raoult Modificada (Van Laar) y según la volatilidad
para la 1ra corrida de datos.
  Se recomienda dejar el equipo durante un
TXY a P cte tiempo prologado para que este ocurra en
Experimen- un sistema lo más estacionario posible.
tal XT  Se recomienda revisar la preparación de la
112
Raoult XT solución a destilar.
107 Van Laar
XT
 Al momento de tomar las muestras de cada
Volatilidad plato del destilador se recomienda tapar
102
XT rápidamente con papel aluminio para que el
97 Experimen- etanol no volatilice.
tal YT  Dejar enfriar la muestra lo suficiente (25 –
T (ºC)

92 Raoult YT
31 ºC) para poder poner al refractómetro y
Van Laar
87 YT no dañar el equipo.
 Tener en cuenta las temperaturas para la
82 formación del azeótropo. Caso que se
77
llegara a esta situación ya no se podría
destilar y se tendría que reducir la presión
72 para lograr la separación de los
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
componentes de la mezcla binaria.
x, y
Gráfico 4. Comparación de las gráficas T-X-Y a P
constante obtenidos según Ley de Raoult, Ley de
Raoult Modificada (Van Laar) y según la volatilidad
para la 2da corrida de datos.

3.2 Discusión
En las 4 gráficas mostradas se puede notar que un
gran aproximamiento a los datos reales
(experimentales) es brindada por la Ley de Raoult
Modificada, con errores menores en comparación a
los otros métodos (Ley de Raoult y Volatilidad).
Para dicha mezcla a a 1 atm aparece el azeótropo a
89.4% en mol de etanol y 78.2 ºC. No consiguiendo
así la formación del azeótropo con los datos
obtenidos; teniendo así, la separación de la mezcla
etanol (1) / agua (2) sin ningún inconveniente por
esa parte.

4. CONCLUSIONES
 Se concluye que el método de Vaan Laar es
el que mejor se ajusta en cuanto a los datos
experimentales tanto para las gráficas T-X-
Y como para las XY a P constante.
 Se concluye que según las gráficas 1 y 3 se
obtiene una gráfica de McCabe -Thiele
ideal según el método de la volatilidad y la
Ley de Raoult.
 Para este tipo de mezclas se tiene que tener
mucho cuidado con su azeótropo,
concluyendo así que el manejo del equipo y
las lecturas a tomar deben ser muy
cuidadosas.

5. RECOMENDACIONES
6. BIBLIOGRAFIA

[1] J. Ocon garcía y G. Tojo Barreiro, Problemas de ingeniería Química Operaciones Básicas
Tomo II, vol. II, Aguilar.
[2] R. E. Treybal, Operaciones de Transferencia de Masa, McGRAW-HILL, pp. 307-369.

1.3569 37.74
1.3600 47.2
APÉNDICES 1.3626 57.27
1.3643 67.54
Para la composición en fase liquida ( x 1) medida
1.365 76.27
en cada plato, se determinó la composición en
1.3651 85.66
fase vapor ( y 1) en equilibrio y la temperatura de
burbuja. Del cuadro a
Para la Primera base de Datos
a) En base a los resultados obtenidos en el IR ( 25 ° C )=1.339 → %m=9.375
plato superior (curva experimental). Para
esto, se asumirá que las composiciones
1
obtenidas se encuentran en equilibrio. x c=
Tomamos muestras en el calderín y los platos y las Me
∗100−%m
analizamos en el refractómetro M H 2O
1+
%m
Plato IR(T) T
1
Calderín 1.3363 31 xc=
46
1 1.3417 31 ∗100−9.375
18
2 1.3425 30 1+
9.375
3 1.3539 31
4 1.3592 31 x c =0.0389
5 1.3612 31
6 1.3617 31 De la misma forma para los platos
7 1.3619 32
8 1.3618 32 x 1=0.0740 → y c =0.0740

Las temperaturas del calderín y los platos están

IR ( 25 ° C )=IR (T ) + 4.5∗10−4∗(T −25) reportadas en el programa SADEC-UDDC.

Tomando como referencia el calderín T c =99.0617 ° C y T 1=90.9185 ° C

IR ( 25 °C )=1.3363+ 4.5∗10−4∗(31−25) a) Mediante el programa Octave, en base a la

Ley de Raoult hallaremos el y i y T i.
IR ( 25 ° C )=1.339
Para eso se haremos los algoritmos
necesarios.
Tabla anexos
IR(25°C) %m %GRAFICA TXY (P=Cte)
1.3392 9.67 P=1.01325;
1.3456 18.58 %Para una mezcla de Etanol - Agua
1.3522 28.29
 Antoine=[5.33675 1648.220 xlabel('z: COMPOSICION');
230.918;5.11564 1687.537 ylabel('T: ºC')
230.170];
[Tsat]=TempSaturacion(Antoine,P) En donde para el Calderín

x = [0 0.0389 0.0740 0.0764 0.1879 0.3193 x c =0.0389

0.4207 0.4611 0.6414 0.7 0.8 0.9 1]';
n=length(x); Corriendo a ese valor de x c experimental
calculamos
%Modelo Van Laar
B1=1.54; y c =0.0834 y T c =98.750 ° C
B2=0.97;
TVL1=zeros(n,1);yVL1=zeros(n,1);
for i=1:n a) Mediante Excel , en base a la volatilidad
[TVL,yVL]=TVanLaar(x(i),P,B1,B2,Antoi relativa α que lo calcularemos por el
ne); método de mínimos cuadrados utilizando
TVL1(i)=TVL; los valores que se adjuntaron en la guía de
yVL1(i)=yVL; laboratorio como anexo 5.1
end Partiremos de la ecuación brindada para
calcular un y supuesto
%Tabla de Datos
fprintf(' x TVL yVL\n'); α x1
for i=1:length(x)
y¿=
1+ ( α−1 ) x 1
fprintf('%6.3f %6.1f %6.3f\
Entonces para un x 1 tendremos un y supuesto
n',x(i),TVL1(i),yVL1(i));
end De la tabla
x 1=0.0028
plot(x,yVL1,"b",x,x,"r");
xlabel('x'); y 1=0.032
ylabel('y')
Asumiendo un α =3
%INGRESE LAS CONSTANTES DE
ANTOINE α x1
Antoine = [5.33675 1648.220 y¿=
1+ ( α−1 ) x 1
230.918;5.11564 1687.537
230.170];
% ********* GRAFICA P-x-y ********* → y ¿=0.00835
%COMPOSICION 2
x = [0 0.0389 0.0740 0.0764 0.1879 0.3193 Restaremos y lo elevamos al cuadrado ( y ¿ − y 1) y
0.4207 0.4611 0.6414 1]'; así haremos los demás valores que se presentan en
n = length( x ); la tabla.
%INGRESE PRESION DE TRABAJO EN Al final con el método de mínimos cuadrados
BAR ajustaremos los valores a una condición más teórica.
P = 1.01325; Dando así obtendremos un α =4.2610
%TABLA T-x-y
fprintf(' x TRaoult YRaoult\n'); Ahora para un x c =0.0389 hallaremos un y c
[Tburi,Yburi] = TY( P , Antoine , x );
for i=1:n y c =0.14712
fprintf('%8.4f %8.6f %8.4f\
n',x(i),Tburi(i),Yburi(i)); Para el cálculo de la Tbur lo haremos por Ley de
end Raoult
plot(x,Tburi,"b",Yburi,Tburi,"b");
Psat ¿( y ¿¿ c∗P)/ x 1 ¿
 n=length(x);
Con esto calculamos el Psat y aplicando la ley de
Antoine hallamos el T bur =115.828 °C %Modelo Van Laar
B1=1.54;
a) Mediante el programa Octave, en base a la B2=0.97;
TVL1=zeros(n,1);yVL1=zeros(n,1);
Ley de Raoult Modificada hallaremos el y i
for i=1:n
y T i. [TVL,yVL]=TVanLaar(x(i),P,B1,B2,Antoine);
Se harán los algoritmos necesarios TVL1(i)=TVL;
%GRAFICA TXY (P=Cte) yVL1(i)=yVL;
P=1.01325; end
%Para una mezcla de Etanol - Agua
Antoine=[5.33675 1648.220
230.918;5.11564 1687.537 %Tabla de Datos
230.170]; fprintf(' x TVL yVL\n');
[Tsat]=TempSaturacion(Antoine,P) for i=1:length(x)
fprintf('%6.3f %6.1f %6.3f\
x = [0 0.0389 0.0740 0.0764 0.1879 0.3193 n',x(i),TVL1(i),yVL1(i));
0.4207 0.4611 0.6414 0.7 0.8 0.9 0.95 1]'; end
n=length(x);
plot(x,TVL1,"r",yVL1,TVL1,"r");
%Modelo Van Laar xlabel('Composición');
B1=1.54; ylabel('T: ºC');
B2=0.97;
TVL1=zeros(n,1);yVL=zeros(n,1); En donde para el Calderín
for i=1:n
[TVL,yVL]=TVanLaar(x(i),P,B1,B2,Antoine); x c =0.0389
TVL1(i)=TVL;
yVL1(i)=yVL;
Corriendo a ese valor de x c experimental
end
calculamos
%Tabla de Datos
fprintf(' x TVL yVL\n'); y c =0.262 y T c =92.7 °C
for i=1:length(x)
fprintf('%6.3f %6.1f %6.3f\
n',x(i),TVL1(i),yVL1(i));
end

plot(x,yVL1,"r",yVL1,yVL1,"b");
xlabel('Composición x'); Para la segunda base de datos
ylabel('Composición y'); a) En base a los resultados obtenidos en el
plato superior (curva experimental). Para
esto, se asumirá que las composiciones
%GRAFICA TXY (P=Cte) obtenidas se encuentran en equilibrio.
P=1.01325; Tomamos muestras en el calderín y los platos y
%Para una mezcla de Etanol - Agua las analizamos en el refractómetro.
Antoine=[5.33675 1648.220
230.918;5.11564 1687.537
230.170];
[Tsat]=TempSaturacion(Antoine,P)

x = [0 0.0389 0.0740 0.0764 0.1879 0.3193

0.4207 0.4611 0.6414 1]';
 x 1=0.0188 → y c =0.0188

Las temperaturas del calderín y los platos están

reportadas en el programa SADEC-UDDC.

T c =99.99 ° C y T 1=98.19 ° C

b) Plato IR(T) T
Calderín 1.3331 31
1 1.3334 30
2 1.3377 29
3 1.3335 29
IR ( 25 ° C )=IR (T ) + 4.5∗10−4∗(T −25) 4 1.3568 31
5 1.3613 29
Tomando como referencia el calderín 6 1.3614 29
7 1.3617 30
IR ( 25 ° C )=1.3331+ 4.5∗10−4∗(31−25) 8 1.3618 30
Mediante el programa Octave, en base a la Ley de
IR ( 25 ° C )=1.3358 Raoult hallaremos el y i y T i.

IR(25°C) %m Para eso se haremos los algoritmos

1.3358 4.8852598 necesarios.
1.33565 4.6747874
1.3395 10.076913 %GRAFICA TXY (P=Cte)
P=1.01325;
1.3353 4.183685
%Para una mezcla de Etanol - Agua
1.3595 45.782541
Antoine=[5.33675 1648.220
1.3631 60.126154 230.918;5.11564 1687.537
1.3632 60.733846 230.170];
1.36395 65.291538 [Tsat]=TempSaturacion(Antoine,P)
1.36405 65.899231
x = [0.019702, 0.018828, 0.042008,
Del cuadro a 0.016799, 0.248362, 0.37109, 0.37704,
0.423996 ]';
IR ( 25 °C )=1.3358 → %m=4.885 n=length(x);

1 %Modelo Van Laar

x c= B1=1.54;
Me
∗100−%m B2=0.97;
M H 2O TVL1=zeros(n,1);yVL1=zeros(n,1);
1+
%m for i=1:n
[TVL,yVL]=TVanLaar(x(i),P,B1,B2,Antoi
1 ne);
xc= TVL1(i)=TVL;
46
∗100−4.885 yVL1(i)=yVL;
18
1+ end
4.885

x c =0.0197 %Tabla de Datos

fprintf(' x TVL yVL\n');
De la misma forma para los platos for i=1:length(x)
 fprintf('%6.3f %6.1f %6.3f\
n',x(i),TVL1(i),yVL1(i)); Entonces para un x 1 tendremos un y supuesto
end De la tabla
x 1=0.0028
plot(x,yVL1,"b",x,x,"r");
xlabel('x'); y 1=0.032
ylabel('y')
Asumiendo un α =3
%INGRESE LAS CONSTANTES DE
ANTOINE
α x1
y¿=
Antoine = [5.33675 1648.220 1+ ( α−1 ) x 1
230.918;5.11564 1687.537
230.170]; → y ¿=0.00835
% ********* GRAFICA P-x-y *********
%COMPOSICION 2
Restaremos y lo elevamos al cuadrado ( y ¿ − y 1) y
x = [0 0.0389 0.0740 0.0764 0.1879 0.3193
así haremos los demás valores que se presentan en
0.4207 0.4611 0.6414 1]';
la tabla.
n = length( x );
Al final con el método de mínimos cuadrados
%INGRESE PRESION DE TRABAJO EN
ajustaremos los valores a una condición más teórica.
BAR
Dando así obtendremos un α =4.2610
P = 1.01325;
%TABLA T-x-y
fprintf(' x TRaoult YRaoult\n'); Ahora para un x c =0.019702 hallaremos un y c
[Tburi,Yburi] = TY( P , Antoine , x );
for i=1:n y c =0.0788
fprintf('%8.4f %8.6f %8.4f\
n',x(i),Tburi(i),Yburi(i)); Para el cálculo de la Tbur lo haremos por Ley de
end Raoult
plot(x,Tburi,"b",Yburi,Tburi,"b"); Psat ¿( y ¿¿ c∗P)/ x 1 ¿
xlabel('z: COMPOSICION');
ylabel('T: ºC') Con esto calculamos el Psat y aplicando la ley de
Antoine hallamos el T bur =117.6488 °C
En donde para el Calderín
d) Mediante el programa Octave, en base a la
x c =0.019702
Ley de Raoult Modificada hallaremos el y i
x c experimental y T i.
Corriendo a ese valor de
Se harán los algoritmos necesarios
calculamos
%GRAFICA TXY (P=Cte)
y c =0.0432 y T c =99.396° C
P=1.01325;
%Para una mezcla de Etanol - Agua
c) Mediante Excel , en base a la volatilidad Antoine=[5.33675 1648.220
relativa α que lo calcularemos por el 230.918;5.11564 1687.537
método de mínimos cuadrados utilizando 230.170];
los valores que se adjuntaron en la guía de [Tsat]=TempSaturacion(Antoine,P)
laboratorio como anexo 5.1
Partiremos de la ecuación brindada para x = [0.019702, 0.018828, 0.042008, 0.016799,
calcular un y supuesto 0.248362, 0.37109, 0.37704, 0.423996]';
n=length(x);
α x1
y¿= %Modelo Van Laar
1+ ( α−1 ) x 1
B1=1.54; ylabel('T: ºC');
B2=0.97;
TVL1=zeros(n,1);yVL=zeros(n,1); En donde para el Calderín
for i=1:n
[TVL,yVL]=TVanLaar(x(i),P,B1,B2,Antoine); x c =0.019702
TVL1(i)=TVL;
yVL1(i)=yVL; Corriendo a ese valor de x c experimental
end calculamos
%Tabla de Datos y c =0.1613 y T c =95.729° C
fprintf(' x TVL yVL\n');
for i=1:length(x)
fprintf('%6.3f %6.1f %6.3f\
n',x(i),TVL1(i),yVL1(i));
end

plot(x,yVL1,"r",yVL1,yVL1,"b");
xlabel('Composición x');
ylabel('Composición y');

%GRAFICA TXY (P=Cte)

P=1.01325;
%Para una mezcla de Etanol - Agua
Antoine=[5.33675 1648.220
230.918;5.11564 1687.537
230.170];
[Tsat]=TempSaturacion(Antoine,P)

x = [0 0.0389 0.0740 0.0764 0.1879 0.3193

0.4207 0.4611 0.6414 1]';
n=length(x);

%Modelo Van Laar

B1=1.54;
B2=0.97;
TVL1=zeros(n,1);yVL1=zeros(n,1);
for i=1:n
[TVL,yVL]=TVanLaar(x(i),P,B1,B2,Antoine);
TVL1(i)=TVL;
yVL1(i)=yVL;
end

%Tabla de Datos
fprintf(' x TVL yVL\n');
for i=1:length(x)
fprintf('%6.3f %6.1f %6.3f\
n',x(i),TVL1(i),yVL1(i));
end

plot(x,TVL1,"r",yVL1,TVL1,"r");
xlabel('Composición');