UNIVERSIDAD NACIONAL DEL CENTRO DEL PERÚ

FACULTAD DE INGENIERÍA QUÍMICA

ESCUELA PROFESIONAL DE INGENIERÍA QUÍMICA AMBIENTAL

DIAGRAMAS TERMODINAMICO DE EQUILIBRIO LIQUIDO

VAPOR “ACETATO ETILICO - ETANOL”

CATEDRATICO:
Operaciones De Transferencia De Masa En Equilibrio
I
CATEDRATICO:
Dr. Oré Vidalón Salvador
REALIZADO POR:
Flores Mondragón, Jazmin
SEMESTRE:
VI

HUANCAYO - PERÚ
2019
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

MARGULES ISOTERMICO

%MARGULES ISOTERMICO
%Equilibrio Liq Vap ACETATO ETILICO - ETANOL
%CONSTANTES DE ANTOINE PARA HALLAR PRESIONES DE VAPOR
%Para el componente 1: acetato etílico
A1=7.09808;
B1=1238.710;
C1=217;
%Para el componente 2: etanol
A2=8.0449;
B2=1554.3;
C2=222.65;
T=40; %°C
%valores de las fracciones molares
x1=[0.001 0.101 0.201 0.301 0.401 0.501 0.601 0.701 0.801
0.999];
x2=[0.999 0.899 0.799 0.699 0.599 0.499 0.399 0.299 0.199
0.001];
%calculo de las presiones de vapor
pvp1=10^(A1-B1/(C1+T))
pvp2=10^(A2-B2/(C2+T))
%CONSTANTES PARA MARGULES
A12=0.413;
A21=0.4816;
disp(' x1 x2 pvp1 pvp2 landa1
landa2 P ,mmHg y1 y2 ');
disp('========================================================
================================')
for i=1:10
landa1(i)=10^(((1-x1(i))^2)*(A12+(2*x1(i)*(A21-A12))));
landa2(i)=10^(((x1(i))^2)*(A21+(2*(1-x1(i))*(A12-A21))));
P(i)=x1(i)*pvp1*landa1(i)+x2(i)*pvp2*landa2(i);
y1(i)=x1(i)*pvp1*landa1(i)/P(i);
y2(i)=x2(i)*pvp2*landa2(i)/P(i);
disp([x1(i) x2(i) pvp1 pvp2 landa1(i) landa2(i) P(i) y1(i)
y2(i)])
end
%GRAFICANDO
x=[0 1];
y=[0 1];
subplot(1,2,1),plot(x1,y1,'-
r',x,y);legend('MARGULES');title('Diagrama de acetato
etilico(1)-etanol(2)a 40°C ; x1 vs y1')
subplot(1,2,2),plot(x1,P,y1,P);hold on ;legend('MARGULES-
x1,t','MARGULES-y1,t');title('Diagrama de etil acetico(1)-
etanol(2)a 40°C; P vs y1, x1')

DIAGRAMAS TERMODINAMICOS 2
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

CORRIENDO EL PROGRAMA:
pvp1 = 189.7565
pvp2 = 134.0105

x1 x2 pvp1 pvp2 landa1 landa2 P, mmHg y1 y2

====================================================================
0.0010 0.9990 189.7565 134.0105 2.5841 1.0000 134.3669 0.0036 0.9964

0.1010 0.8990 189.7565 134.0105 2.2130 1.0085 163.9070 0.2588 0.7412

0.2010 0.7990 189.7565 134.0105 1.9110 1.0352 183.7327 0.3967 0.6033

0.3010 0.6990 189.7565 134.0105 1.6671 1.0838 196.7432 0.4840 0.5160

0.4010 0.5990 189.7565 134.0105 1.4721 1.1594 205.0787 0.5462 0.4538

0.5010 0.4990 189.7565 134.0105 1.3181 1.2697 210.2163 0.5961 0.4039

0.6010 0.3990 189.7565 134.0105 1.1992 1.4262 213.0160 0.6420 0.3580

0.7010 0.2990 189.7565 134.0105 1.1105 1.6463 213.6832 0.6913 0.3087

0.8010 0.1990 189.7565 134.0105 1.0488 1.9565 211.5935 0.7534 0.2466

0.9990 0.0010 189.7565 134.0105 1.0000 3.0234 189.9722 0.9979 0.0021

DIAGRAMAS TERMODINAMICOS 3
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

DIAGRAMA ACETATO ETILICO – ETANOL x1 vs y1

DIAGRAMA ACETATO ETILICO – ETANOL P vs x1, y1

DIAGRAMAS TERMODINAMICOS 4
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

MARGULES ISOBARICO
function MARGULES_ISOBARICO
%Equilibrio Liq Vap ACETATO ETILICO - ETANOL
%CONSTANTES DE ANTOINE PARA HALLAR PRESSIONES DE VAPOR
%Para el componente 1: acetato etílico
A1=7.09808;
B1=1238.710;
C1=217;
%Para el componente 2: etanol
A2=8.0449;
B2=1554.3;
C2=222.65;
P=760;%mmHg
R=1.987;%cal/kmol
%declarar a la temperatura una variable
syms T
%valores de las fracciones molares
x1=[0 0.001 0.101 0.201 0.301 0.401 0.501 0.601 0.701 0.801 0.901
0.999];
x2=[1 0.999 0.899 0.799 0.699 0.599 0.499 0.399 0.299 0.199 0.099
0.001];
%calculo de las presiones de vapor
pvp1=10^(A1-B1/(C1+T))
pvp2=10^(A2-B2/(C2+T))
%calculo de las temperaturas iniciales
T10=B1/(A1-log10(P))-C1;
T20=B2/(A2-log10(P))-C2;
%CONSTANTES PARA MARGULES
A12=0.2010;
A21=0.1430;
for i=1:12
landa1(i)=10^(((1-x1(i))^2)*(A12+(2*x1(i)*(A21-A12))));
landa2(i)=10^(((x1(i))^2)*(A21+(2*(1-x1(i))*(A12-A21))));
y1(i)=x1(i)*landa1(i)*pvp1/P;
y2(i)=x2(i)*landa2(i)*pvp2/P;
end
%Determinación de las temperaturas de equilibrio
disp(' t ,ºC x1 Y1 Landa1 Landa2');
disp('====================================================')
f=y1+y2-1;
for i = 1:12
xo=(T10+T20)/2;%valor inicial para la temperatura
e=10^-4;%criterio de convergencia
distancia=1;
while distancia > e
fxi=subs(f(i),sym('T'),xo);
derivada=diff(f(i),sym('T'));
Dxfxi=subs(derivada,sym('T'),xo);
x_1=single(xo-(fxi/Dxfxi));
distancia=single(abs(subs(f(i),sym('T'),x_1)));
xo=x_1;
end
Tm(i)=x_1;
t(i)=Tm(i);
DIAGRAMAS TERMODINAMICOS 5
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

Y1(i)=single(subs(y1(i),sym('T'),Tm(i)));
G1(i)=single(subs(landa1(i),sym('T'),Tm(i)));
G2(i)=single(subs(landa2(i),sym('T'),Tm(i)));
disp([t(i) x1(i) Y1(i) G1(i) G2(i)])
end
%REALIZANDO LOS GRAFICOS
x=[0 1];
y=[0 1];
subplot(1,2,1),plot(x1,Y1,'-
r',x,y);legend('MARGULES');title('Diagrama de ACETATO ETILICO-ETANOL
a 760 mmHg; x1 vs y1')
subplot(1,2,2),plot(x1,t,Y1,t);hold on; legend('MARGULES-
x1,t','MARGULES-y2,t');title('Diagrama ACETATO ETILICO-ETANOL a 760
mmHg;T vs x1,y1')

CORRIENDO EL PROGRAMA:
pvp1 = 10^(44363/6250 - 1361976048442409/(1099511627776*(T + 217)))

pvp2 = 10^(80449/10000 - 15543/(10*(T + 4453/20)))

t ,ºC x1 Y1 Landa1 Landa2

========================================
78.3326 0 0 1.5885 1.0000

78.3155 0.0010 0.0017 1.5867 1.0000

76.9324 0.1010 0.1446 1.4223 1.0058

76.0238 0.2010 0.2550 1.2985 1.0222

75.4312 0.3010 0.3476 1.2054 1.0479

75.0680 0.4010 0.4312 1.1361 1.0819

74.8904 0.5010 0.5116 1.0854 1.1231

74.8815 0.6010 0.5931 1.0493 1.1705

75.0429 0.7010 0.6794 1.0249 1.2227

75.3884 0.8010 0.7739 1.0099 1.2781

75.9415 0.9010 0.8799 1.0022 1.3348

76.7143 0.9990 0.9987 1.0000 1.3894

DIAGRAMAS TERMODINAMICOS 6
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

DIAGRAMA ACETATO ETILICO – ETANOL

x1 vs y1

DIAGRAMA ACETATO ETILICO – ETANOL

T vs x1, y1

DIAGRAMAS TERMODINAMICOS 7
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

WILSON ISOTERMICO
function ISOTERMICO_WILSON
%Equilibrio Liq Vap ACETATO ETILICO - ETANOL
%CONSTANTES DE ANTOINE PARA HALLAR PRESSIONES DE VAPOR
%Para el componente 1: acetato etílico
A1=7.09808;
B1=1238.710;
C1=217;
%Para el componente 2: etanol
A2=8.0449;
B2=1554.3;
C2=222.65;
T=40; %°C
R=1.987; %cal/K*mol
%valores de las fracciones molares
x1=[0 0.001 0.101 0.201 0.301 0.401 0.501 0.601 0.701 0.801 0.901
0.999];
x2=[1 0.999 0.899 0.799 0.699 0.599 0.499 0.399 0.299 0.199 0.099
0.001];
%calculo de las presiones de vapor
pvp1=10^(A1-B1/(C1+T))
pvp2=10^(A2-B2/(C2+T))
%CONSTANTES DE WILSON
landa11=-62.43;
landa22=822.03;
v1=102.1;
v2=60.356;
A12=(v2/v1)*exp(-(landa11)/(R*(T+273.15)));
A21=(v1/v2)*exp(-(landa22)/(R*(T+273.15)));
disp(' P ,mmHg x1 Y1 Landa1 Landa2');
disp('====================================================')
for i=1:12
lamba1(i)=(exp((1-x1(i))*((A12/(x1(i)+A12-x1(i)*A12))-
(A21/(A21*x1(i)+1-x1(i))))))/(x1(i)+A12-A12*x1(i));
lamba2(i)=(exp((-x1(i))*((A12/(x1(i)+A12-x1(i)*A12))-
(A21/(A21*x1(i)+1-x1(i))))))/(1-x1(i)+A21*x1(i));
P(i)=x1(i)*pvp1*lamba1(i)+x2(i)*pvp2*lamba2(i);
y1(i)=x1(i)*lamba1(i)*pvp1/P(i);
y2(i)=x2(i)*lamba2(i)*pvp2/P(i);
disp([P(i) x1(i) y1(i) lamba1(i) lamba2(i)])
end
%REALIZANDO LOS GRAFICOS
x=[0 1];
y=[0 1];
subplot(1,2,1),plot(x1,y1,'-r',x,y);legend('WILSON');title('Diagrama
ACETATO ETILICO - ETANOL 760 mmHg; x1 vs y1')
subplot(1,2,2),plot(x1,P,y1,P); legend('WILSON-x1,t','WILSON-
y1,t');title('Diagrama ACETATO ETILICO - ETANOL a 760 mmHg; P vs
x1,Y1')

DIAGRAMAS TERMODINAMICOS 8
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

CORRIENDO EL PROGRAMA:
pvp1 = 189.7565
pvp2 = 134.0105

P ,mmHg x1 Y1 Landa1 Landa2

====================================================
134.0105 0 0 2.6484 1.0000

134.3782 0.0010 0.0037 2.6435 1.0000

164.1349 0.1010 0.2591 2.2187 1.0094

183.5206 0.2010 0.3945 1.8980 1.0379

196.2210 0.3010 0.4810 1.6523 1.0872

204.5082 0.4010 0.5440 1.4622 1.1616

209.7567 0.5010 0.5958 1.3145 1.2679

212.7018 0.6010 0.6436 1.2004 1.4176

213.5055 0.7010 0.6941 1.1141 1.6299

211.6063 0.8010 0.7558 1.0522 1.9378

205.1669 0.9010 0.8447 1.0136 2.4020

189.9856 0.9990 0.9978 1.0000 3.1233

DIAGRAMAS TERMODINAMICOS 9
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

DIAGRAMA ACETATO ETILICO – ETANOL

x1 vs y1

DIAGRAMA ACETATO ETILICO – ETANOL

P vs x1, y1

DIAGRAMAS TERMODINAMICOS 10
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

WILSON ISOBARICO
function WILSON_ISOBARICO
%Equilibrio Liq Vap ACETATO ETILICO - ETANOL
%CONSTANTES DE ANTOINE PARA HALLAR PRESSIONES DE VAPOR
%Para el componente 1: acetato etílico
A1=7.09808;
B1=1238.710;
C1=217;
%Para el componente 2: etanol
A2=8.0449;
B2=1554.3;
C2=222.65;
P=760;%mmHg
R=1.987; %cal/K*mol
%declarar a la temperatura una variable
syms T
%valores de las fracciones molares
x1=[0 0.001 0.101 0.201 0.301 0.401 0.501 0.601 0.701 0.801 0.901 0.999];
x2=[1 0.999 0.899 0.799 0.699 0.599 0.499 0.399 0.299 0.199 0.099 0.001];
%calculo de las presiones de vapor
pvp1=10^(A1-B1/(C1+T));
pvp2=10^(A2-B2/(C2+T));
%calculo de las temperaturas iniciales
T10=B1/(A1-log10(P))-C1;
T20=B2/(A2-log10(P))-C2;
%CONSTANTES DE WILSON
landa11=-178.81;
landa22=844.69;
v1=102.1;
v2=60.356;
A12=(v2/v1)*exp(-(landa11)/(R*(T+273.15)));
A21=(v1/v2)*exp(-(landa22)/(R*(T+273.15)));
for i=1:12
lamba1(i)=(exp((1-x1(i))*((A12/(x1(i)+A12-x1(i)*A12))-
(A21/(A21*x1(i)+1-x1(i))))))/(x1(i)+A12-A12*x1(i));
lamba2(i)=(exp((-x1(i))*((A12/(x1(i)+A12-x1(i)*A12))-(A21/(A21*x1(i)+1-
x1(i))))))/(1-x1(i)+A21*x1(i));
y1(i)=x1(i)*lamba1(i)*pvp1/P;
y2(i)=x2(i)*lamba2(i)*pvp2/P;
end
%Determinación de las temperaturas de equilibrio
disp(' t ,ºC x1 Y1 Lamba1 Lamba2');
disp('====================================================')
f=y1+y2-1;
for i = 1:12
xo=(T10+T20)/2;%valor inicial para la temperatura
e=10^-4;%criterio de convergencia
distancia=1;
while distancia > e
fxi=subs(f(i),sym('T'),xo);
derivada=diff(f(i),sym('T'));
Dxfxi=subs(derivada,sym('T'),xo);
x_1=single(xo-(fxi/Dxfxi));
distancia=single(abs(subs(f(i),sym('T'),x_1)));
xo=x_1;
end
Tm(i)=x_1;
t(i)=Tm(i);
Y1(i)=single(subs(y1(i),sym('T'),Tm(i)));

DIAGRAMAS TERMODINAMICOS 11
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

G1(i)=single(subs(lamba1(i),sym('T'),Tm(i)));
G2(i)=single(subs(lamba2(i),sym('T'),Tm(i)));
disp([t(i) x1(i) Y1(i) G1(i) G2(i)])
end
%REALIZANDO LOS GRAFICOS
x=[0 1];
y=[0 1];
subplot(1,2,1),plot(x1,Y1,'-r',x,y); legend('WILSON'); title('Diagrama
ACETATO ETILICO - ETANOL 760 mmHg; x1 vs y1')
subplot(1,2,2),plot(x1,t,Y1,t); legend('WILSON-x1,t','WILSON-
y1,t');title('Diagrama ACETATO ETILICO - ETANOL a 760 mmHg; P vs x1,Y1')

CORRIENDO EL PROGRAMA:
t ,ºC x1 Y1 Lamba1 Lamba2
=============================================
78.3326 0 0 1.3113 1.0000

78.3229 0.0010 0.0014 1.3111 1.0000

77.3930 0.1010 0.1330 1.2881 1.0010

76.5470 0.2010 0.2522 1.2623 1.0047

75.8009 0.3010 0.3601 1.2335 1.0126

75.1720 0.4010 0.4576 1.2016 1.0273

74.6802 0.5010 0.5460 1.1666 1.0528

74.3506 0.6010 0.6270 1.1290 1.0966

74.2183 0.7010 0.7027 1.0897 1.1726

74.3475 0.8010 0.7776 1.0507 1.3113

74.9171 0.9010 0.8627 1.0167 1.5902

76.6900 0.9990 0.9979 1.0000 2.2426

DIAGRAMAS TERMODINAMICOS 12
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

DIAGRAMA ACETATO ETILICO – ETANOL

x1 vs y1

DIAGRAMA ACETATO ETILICO – ETANOL

P vs x1, y1

DIAGRAMAS TERMODINAMICOS 13
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

VAN LAAR ISOTERMICO

function VAN_LAAR_ISOTERMICO
%Equilibrio Liq Vap ACETATO ETILICO - ETANOL
%CONSTANTES DE ANTOINE PARA HALLAR PRESSIONES DE VAPOR
%Para el componente 1: acetato etílico
A1=7.09808;
B1=1238.710;
C1=217;
%Para el componente 2: etanol
A2=8.0449;
B2=1554.3;
C2=222.65;
T=40; %°C
%valores de las fracciones molares
x1=[0 0.001 0.101 0.201 0.301 0.401 0.501 0.601 0.701 0.801 0.901 0.999];
x2=[1 0.999 0.899 0.799 0.699 0.599 0.499 0.399 0.299 0.199 0.099 0.001];
%calculo de las presiones de vapor
pvp1=10^(A1-B1/(C1+T))
pvp2=10^(A2-B2/(C2+T))
%CONSTANTES PARA VAN LAAR
A12=0.4151;
A21=0.4833;
disp(' P ,mmHg x1 y1 Lamba1 Lamba2');
disp('====================================================')
for i=1:12
lamba1(i)=10^((A12)/(1+(x1(i)/(1-x1(i)))*(A12/A21))^2);
lamba2(i)=10^((A21)/(1+((1-x1(i))/(x1(i)))*(A12/A21))^2);
P(i)=x1(i)*pvp1*lamba1(i)+x2(i)*pvp2*lamba2(i);
y1(i)=x1(i)*lamba1(i)*pvp1/P(i);
y2(i)=x2(i)*lamba2(i)*pvp2/P(i);
disp([P(i) x1(i) y1(i) lamba1(i) lamba2(i)])
end
%GRAFICANDO
x=[0 1];
y=[0 1];
subplot(1,2,1),plot(x1,y1,'r',x,y);legend('MARGULES');title('Diagrama
ACETATO ETILICO - ETANOL 760 mmHg; x1 vs y1')
subplot(1,2,2),plot(x1,P,y1,P);legend('MARGULES-x1,t','MARGULES-
y1,t');title('Diagrama ACETATO ETILICO - ETANOL a 760 mmHg; P vs x1,Y1')

CORRIENDO EL PROGRAMA:
pvp1 = 189.7565
pvp2 = 134.0105
P ,mmHg x1 y1 Lamba1 Lamba2
====================================================
134.0105 0 0 2.6008 1.0000
134.3694 0.0010 0.0037 2.5965 1.0000
164.7226 0.1010 0.2576 2.2144 1.0150
186.1608 0.2010 0.3910 1.9085 1.0588

DIAGRAMAS TERMODINAMICOS 14
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

201.1050 0.3010 0.4727 1.6642 1.1321

211.2408 0.4010 0.5295 1.4701 1.2380
217.6226 0.5010 0.5755 1.3173 1.3816
220.7353 0.6010 0.6196 1.1992 1.5704
220.5074 0.7010 0.6702 1.1110 1.8150
216.2744 0.8010 0.7374 1.0493 2.1294
206.6838 0.9010 0.8374 1.0124 2.5324
189.9740 0.9990 0.9979 1.0000 3.0372

DIAGRAMA ACETATO ETILICO – ETANOL

x1 vs y1

DIAGRAMAS TERMODINAMICOS 15
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

DIAGRAMA ACETATO ETILICO – ETANOL

T vs x1, y1

DIAGRAMAS TERMODINAMICOS 16
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

VAN LAAR ISOBARICO

function VAN_LAAR_ISOBARICO
%Equilibrio Liq Vap ACETATO ETILICO - ETANOL
%CONSTANTES DE ANTOINE PARA HALLAR PRESSIONES DE VAPOR
%Para el componente 1: acetato etílico
A1=7.09808;
B1=1238.710;
C1=217;
%Para el componente 2: etanol
A2=8.0449;
B2=1554.3;
C2=222.65;
P=760;%mmHg
R=1.987;%cal/kmol
%declarar a la temperatura una variable
syms T
%valores de las fracciones molares
x1=[0 0.001 0.101 0.201 0.301 0.401 0.501 0.601 0.701 0.801 0.901 0.999];
x2=[1 0.999 0.899 0.799 0.699 0.599 0.499 0.399 0.299 0.199 0.099 0.001];
%calculo de las presiones de vapor
pvp1=10^(A1-B1/(C1+T));
pvp2=10^(A2-B2/(C2+T));
%calculo de las temperaturas iniciales
T10=B1/(A1-log10(P))-C1;
T20=B2/(A2-log10(P))-C2;
%CONSTANTES PARA VAAN_LAAR
A12=0.3311;
A21=0.3972;
for i=1:12
lamba1(i)=10^((A12)/(1+(x1(i)/(1-x1(i)))*(A12/A21))^2);
lamba2(i)=10^((A21)/(1+((1-x1(i))/(x1(i)))*(A12/A21))^2);
y1(i)=x1(i)*lamba1(i)*pvp1/P;
y2(i)=x2(i)*lamba2(i)*pvp2/P;
end
%Determinación de las temperaturas de equilibrio
disp(' t ,ºC x1 Y1 Lamba1 Lamba2');
disp('====================================================')
f=y1+y2-1;
for i = 1:12
xo=(T10+T20)/2;%valor inicial para la temperatura
e=10^-4;%criterio de convergencia
distancia=1;
while distancia > e
fxi=subs(f(i),sym('T'),xo);
derivada=diff(f(i),sym('T'));
Dxfxi=subs(derivada,sym('T'),xo);
x_1=single(xo-(fxi/Dxfxi));
distancia=single(abs(subs(f(i),sym('T'),x_1)));
xo=x_1;
end
Tm(i)=x_1;
t(i)=Tm(i);
Y1(i)=single(subs(y1(i),sym('T'),Tm(i)));
G1(i)=single(subs(lamba1(i),sym('T'),Tm(i)));
G2(i)=single(subs(lamba2(i),sym('T'),Tm(i)));
disp([t(i) x1(i) Y1(i) G1(i) G2(i)])
end

DIAGRAMAS TERMODINAMICOS 17
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

%REALIZANDO LOS GRAFICOS

x=[0 1];
y=[0 1];
subplot(1,2,1),plot(x1,Y1,'r',x,y);legend('VAN LAAR'); title('Diagrama
ACETATO ETILICO - ETANOL a 40°C; x1 vs y1')
subplot(1,2,2),plot(x1,t,Y1,t);hold on; legend('VAN LAAR-x1,t','VAN LAAR-
Y1,t');title('Diagrama ACETATO ETILICO - ETANOL a 40°C; P vs x1,Y1')

CORRIENDO EL PROGRAMA:
t ,ºC x1 Y1 Lamba1 Lamba2
====================================================
78.3326 0 0 2.1434 1.0000

78.3007 0.0010 0.0023 2.1407 1.0000

75.5813 0.1010 0.1839 1.8916 1.0130

73.6283 0.2010 0.3052 1.6837 1.0504

72.2601 0.3010 0.3915 1.5111 1.1119

71.3615 0.4010 0.4583 1.3690 1.1989

70.8681 0.5010 0.5154 1.2535 1.3139

70.7671 0.6010 0.5710 1.1617 1.4608

71.1022 0.7010 0.6330 1.0913 1.6450

71.9973 0.8010 0.7114 1.0410 1.8733

73.6998 0.9010 0.8229 1.0104 2.1545

76.6830 0.9990 0.9977 1.0000 2.4919

DIAGRAMAS TERMODINAMICOS 18
 “UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU”
FACULTAD DE INGENIERIA QUIMICA

DIAGRAMA ACETATO ETILICO – ETANOL

x1 vs y1

DIAGRAMA ACETATO ETILICO – ETANOL

T vs x1, y1

DIAGRAMAS TERMODINAMICOS 19