Metodo Margules, Van Lar y Wilson

UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU
FACULTAD DE INGENIERIA QUIMICA

CARRERA PROFESIONAL DEL INGENIERIA QUIMICA
TERMODINAMICA DE LOS PROCESOS QUIMICOS 2
APELLIDOS Y NOMBRES: HURTADO TOCAS ANYELO YONEL

CODIGO: 2016101455B
1. PROGRAMACION CON EL METODO DE VAN LAAR:
van_laar_isotermico
clc, clear all

% ideal isotermico
% Equilibrio Liquido Vapor - Acetona(1) - Metanol (2)
% Datos de la acetona
x1e=[0.0281 0.0670 0.0644 0.0858 0.1046 0.1357 0.1454 0.1663 0.2173
0.2390 0.2787 0.3579 0.4050 0.4480 0.5052 0.5432 0.6332 0.6538 0.6605
0.6945 0.7327 0.7525 0.7752 0.7922 0.9080 0.9088 0.9197 0.9448];
% Parte gaseosa de la acetona
x2e=1-x1e;
% Datos del metanol
y1e=[0.0647 0.1295 0.1407 0.1848 0.2190 0.2637 0.2694 0.3055 0.3633
0.3863 0.4184 0.4779 0.5135 0.5512 0.5844 0.6174 0.6772 0.6849 0.6926
0.7124 0.7383 0.7618 0.7729 0.7876 0.8959 0.8963 0.8941 0.9336];
% Parte gaseosa del metnol
y2e=1-y1e;
% Presiones variables ya que es un oroceso isotermico
Pe=[542.17 564.61 569.56 581.45 592.15 610.13 618.98 628.16 650.74 657.70
675.68 699.07 712.65 722.76 732.37 738.49 748.61 752.18 749.65 752.11
753.53 753.85 757.52 757.97 749.10 750.31 750.47 748.52];
% Temperatura constante en kelvin
T=55+273.15;
%Datos encontrados en tablas termodinamicas
VPA1=-7.45514;
VPB1=1.202;
VPC1=-2.43926;
VPD1=-3.36690;
% Presion y temperatura critica de la acetona
Pc1=47;
Tc1=508.1;
VPA2=-8.54796;
VPB2=0.76982;
VPC2=-3.1085;
VPD2=1.54481;
% Presion y temperatura critica del metanol
Pc2=80.9;
Tc2=512.6;
% Calculo de las presiones de vapor
TR1=1-T/Tc1;
LnPvp1=(VPA1*TR1+VPB1*TR1^1.5+VPC1*TR1^3+VPD1*TR1^6)/(1-TR1);
Pvp1=exp(LnPvp1)*Pc1*750,062;
TR2=1-T/Tc2;
% Aplicacion de la ecuacion de van laar
x1=[0:0.01:1];
% Constantes de van laar conseguido en tablas
A12=0.2635;
A21=0.2801;
lamda1=10.^(A12./(1+(x1*A12)./((1-x1)*A21)).^2)
lamda2=10.^(A21./(1+((1-x1)*A21)./(x1*A12)).^2)
P=Pvp1*lamda1.*x1+Pvp2*lamda2.*(1-x1);
y1=Pvp1*x1.*lamda1./P
x=[0:0.01:1];
y=[0:0.01:1];
% Grafica
subplot(2,1,1), plot(x1e,y1e,'+')
hold on
plot(x1,y1,'g')
plot(x,y)
title('ISOTERMICO: Equilibrio Liquido Vapor - Acetona - Metanol'),
xlabel('x1'),ylabel('y1')
legend('datos experimentales','VAN LAAR','DIAGONAL')
grid
subplot(2,1,2),plot(x1,P,y1,P)
hold on
plot(x1e,Pe,'+',y1e,Pe,'*')
xlabel('x1 y1'),ylabel('P(mmHg)')
legend('DATOS EXPERIMENTALES','DATOS EXPERIMENTALES','P Vs x1','P Vs y1')
grid
ISOTERMICO: Equilibrio Liquido Vapor - Acetona - Metanol

1
0.5
y1
datos experimentales
VAN LAAR
DIAGONAL
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1
800
700
P(mmHg)
DATOS EXPERIMENTALES
600 DATOS EXPERIMENTALES
P Vs x1
P Vs y1
500
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1 y1
2. PROGRAMACION CON LA ECUACION DE MARGULES:
margules_isotermico
clc, clear all

% Datos experimentales encontrados
x1e=[0.0281 0.0670 0.0644 0.0858 0.1046 0.1357 0.1454 0.1663 0.2173
0.2390 0.2787 0.3579 0.4050 0.4480 0.5052 0.5432 0.6332 0.6538 0.6605
0.6945 0.7327 0.7525 0.7752 0.7922 0.9080 0.9088 0.9197 0.9448];
x2e=1-x1e;
y1e=[0.0647 0.1295 0.1407 0.1848 0.2190 0.2637 0.2694 0.3055 0.3633
0.3863 0.4184 0.4779 0.5135 0.5512 0.5844 0.6174 0.6772 0.6849 0.6926
0.7124 0.7383 0.7618 0.7729 0.7876 0.8959 0.8963 0.8941 0.9336];
y2e=1-y1e;
Pe=[542.17 564.61 569.56 581.45 592.15 610.13 618.98 628.16 650.74 657.70
675.68 699.07 712.65 722.76 732.37 738.49 748.61 752.18 749.65 752.11
753.53 753.85 757.52 757.97 749.10 750.31 750.47 748.52];
T=55+273.15; %Kelvin
R=1.987; %cal/mol*K
%Datos de las tablas
VPA1=-7.45514;
VPB1=1.202;
VPC1=-2.43926;
VPD1=-3.36690;
%Presion y temperatura critica 1
Pc1=47;
Tc1=508.1;
VPA2=-8.54796;
VPB2=0.76982;
VPC2=-3.1085;
VPD2=1.54481;
Pc2=80.9;
Tc2=512.6;
% Calculamos las presiones de vapor
TR1=1-T/Tc1;
TR2=1-T/Tc2;
% Programacion con la ecuacion de margules
x1=[0:0.01:1];
% Constantes de margules
A12=0.2762;
A21=0.2877;
lamda1=10.^((1-x1).^2.*(A12+2*x1*(A21-A12)))
lamda2=10.^(x1.^2.*(A21+2*(1-x1)*(A12-A21)))
y1=Pvp1*x1.*lamda1./P;
x=[0:0.01:1];
y=[0:0.01:1];
% Graficando la ecuacion
hold on
plot(x1,y1,'g')
plot(x,y,'r')
legend('datos experimentales','MARGULES','DIAGONAL')
grid
hold on
grid
1
MARGULES
DIAGONAL
0.5
y1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1
800
700
P Vs x1
P(mmHg)
P Vs y1
600
500
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1 y1
3. PROGRAMACION CON LA ECUACION DE WILSON:
wilson_isotermico
clc, clear all

x1e=[0.0281 0.0670 0.0644 0.0858 0.1046 0.1357 0.1454 0.1663 0.2173
0.2390 0.2787 0.3579 0.4050 0.4480 0.5052 0.5432 0.6332 0.6538 0.6605
0.6945 0.7327 0.7525 0.7752 0.7922 0.9080 0.9088 0.9197 0.9448];
x2e=1-x1e;
y1e=[0.0647 0.1295 0.1407 0.1848 0.2190 0.2637 0.2694 0.3055 0.3633
0.3863 0.4184 0.4779 0.5135 0.5512 0.5844 0.6174 0.6772 0.6849 0.6926
0.7124 0.7383 0.7618 0.7729 0.7876 0.8959 0.8963 0.8941 0.9336];
y2e=1-y1e;
Pe=[542.17 564.61 569.56 581.45 592.15 610.13 618.98 628.16 650.74 657.70
675.68 699.07 712.65 722.76 732.37 738.49 748.61 752.18 749.65 752.11
753.53 753.85 757.52 757.97 749.10 750.31 750.47 748.52];
T=55+273.15; %K
R=1.987; %cal/mol*K
%Datos encontrados en tablas
VPA1=-7.45514;
VPB1=1.202;
VPC1=-2.43926;
VPD1=-3.36690;
Pc1=47;
Tc1=508.1;
VPA2=-8.54796;
VPB2=0.76982;
VPC2=-3.1085;
VPD2=1.54481;
Pc2=80.9;
Tc2=512.6;
%Calculamos las presiones de vapor
TR1=1-T/Tc1;
TR2=1-T/Tc2;
%Programacion con el metodo de wilson
x1=[0:0.01:1];
% Constantes de Wilson
lamba12=-214.95; %cal/mol-g
lamba21=664.08; %cal/mol-g
v2=42.4809;
v1=77.3603;
A12=v2/v1*exp(-(lamba12)/(R*T));
A21=v1/v2*exp(-(lamba21)/(R*T));
lamda1=exp(-log(x1+A12*(1-x1))+(1-x1).*((A12./(x1+A12*(1-x1)))-(A21./
(A21*x1+(1-x1)))));
lamda2=exp(-log((1-x1)+A21*x1)-x1.*((A12./(x1+A12*(1-x1)))-(A21./
(A21*x1+(1-x1)))));
x=[0:0.01:1];
y=[0:0.01:1];
%Graficando las ecuaciones
hold on
plot(x1,y1,'g')
plot(x,y,'r')
legend('datos experimentales','WLSON','DIAGONAL')
grid
hold on
grid
1
0.5
y1
WLSON
DIAGONAL
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1
800
700
P(mmHg)
P Vs x1
P Vs y1
500
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1 y1
4. PROGRAMACION CON EL METODO VAAN LAR EN UN SISTEMA ISOBARICO:
disp('***********PROGRAMA SISTEMA ISOBARICO (ACETONA-

METANOL)*****************')
clc, clear all
x1e=[1 0.93 0.867 0.819 0.783 0.773 0.735 0.66 0.594 0.554 0.519 0.407
0.394 0.369 0.324 0.281 0.263 0.229 0.195 0.151 0.1 0.074 0.053 0.024 0];
x2e=1-x1e;
y1e=[1 0.918 0.856 0.812 0.782 0.774 0.745 0.689 0.644 0.618 0.594 0.514
0.504 0.485 0.45 0.41 0.392 0.357 0.319 0.265 0.191 0.148 0.11 0.053 0];
y2e=1-y1e;
%TEMPERATURAS VARIABLES YA QUE ES UN PROCESO ISOBARICO
Te=[56.53 55.67 55.43 55.31 55.24 55.25 55.3 55.39 55.74 55.96 56.15 56.9
57.05 57.29 57.69 58.32 58.49 58.97 59.57 60.41 61.53 62.21 62.79 63.69
64.7];
Te1=273.15+Te
%PRESION CONSTANTE
P=800;
%DATOS DE LA ACETONA
VPA1=-7.45514;
VPB1=1.202;
VPC1=-2.43926;
VPD1=-3.36690;
%PRESION Y TEMPERATURA CRITICA DE LA ACETONA
Pc1=47;
Tc1=508.1;
%DATOS DEL METANOL
VPA2=-8.54796;
VPB2=0.76982;
VPC2=-3.1085;
VPD2=1.54481;
%PRESION Y TEMPERATURA CRITICA DEL METANOL
Pc2=80.9;
Tc2=512.6;
%CALCULAMOS LAS PRESIONES DE VAPOR
TR1=1-Te1./Tc1
LnPvp1=(VPA1*TR1+VPB1*TR1.^1.5+VPC1*TR1.^3+VPD1*TR1.^6)/(1-TR1);
TR2=1-Te1./Tc2;
%PROGRAMACION CON EL METODO DE VAN LAAR
x1=[1 0.93 0.867 0.819 0.783 0.773 0.735 0.66 0.594 0.554 0.519 0.407
0.394 0.369 0.324 0.281 0.263 0.229 0.195 0.151 0.1 0.074 0.053 0.024 0];
%CONSTANTES DE VAN LAAR
A12=0.2635;
A21=0.2801;
lamda1=10.^(A12./(1+(x1*A12)./((1-x1)*A21)).^2)
lamda2=10.^(A21./(1+((1-x1)*A21)./(x1*A12)).^2)
y1=Pvp1.*x1.*lamda1./P
x=[0:0.01:1];
y=[0:0.01:1];
%GRAFICANDO
hold on
plot(x1,y1,'r')
plot(x,y,'g')
title('ISOBARICO: Equilibrio Liquido Vapor - Acetona - Metanol'),
grid
subplot(2,1,2),plot(x1,Te1,y1,Te1)
hold on
plot(x1e,Te1,'+',y1e,Te1,'*')
xlabel('x1 y1'),ylabel('T °K')
legend('DATOS EXPERIMENTALES','DATOS EXPERIMENTALES','T Vs x1','T Vs y1')
grid
ISOBARICO: Equilibrio Liquido Vapor - Acetona - Metanol
1.5
1
y1
0.5 datos experimentales

VAN LAAR
DIAGONAL
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1
340
T Vs x1
T °K
T Vs y1
330
325
0 0.2 0.4 0.6 0.8 1 1.2 1.4
x1 y1
5. PROGRAMACION CON EL METODO MARGULES EN UN SISTEMA ISOBARICO:
disp('*********PROCESO ISOBARICO EQUILIBRIO ACETONA-METANOL***********')

clc, clear all
x1e=[1 0.93 0.867 0.819 0.783 0.773 0.735 0.66 0.594 0.554 0.519 0.407
0.394 0.369 0.324 0.281 0.263 0.229 0.195 0.151 0.1 0.074 0.053 0.024 0];
x2e=1-x1e;
y1e=[1 0.918 0.856 0.812 0.782 0.774 0.745 0.689 0.644 0.618 0.594 0.514
0.504 0.485 0.45 0.41 0.392 0.357 0.319 0.265 0.191 0.148 0.11 0.053 0];
y2e=1-y1e;
%ESTAS SON LAS TEMPERATURAS VARIADAS YA QUE ES UN PROCESO ISOBARICO
Te=[56.53 55.67 55.43 55.31 55.24 55.25 55.3 55.39 55.74 55.96 56.15 56.9
57.05 57.29 57.69 58.32 58.49 58.97 59.57 60.41 61.53 62.21 62.79 63.69
64.7];
Te1=273.15+Te
%Presion constante a 800 mmHg
P=800;
R=1.987; %cal/mol*K
VPA1=-7.45514;
VPB1=1.202;
VPC1=-2.43926;
VPD1=-3.36690;
Pc1=47;
Tc1=508.1;
%DATOS DEL METANOL
VPA2=-8.54796;
VPB2=0.76982;
VPC2=-3.1085;
VPD2=1.54481;
Pc2=80.9;
Tc2=512.6;
TR1=1-Te1./Tc1;
TR2=1-Te1./Tc2;
%PROGRAMACION CON EL METODO DE MARGULES
x1=[1 0.93 0.867 0.819 0.783 0.773 0.735 0.66 0.594 0.554 0.519 0.407
0.394 0.369 0.324 0.281 0.263 0.229 0.195 0.151 0.1 0.074 0.053 0.024 0];
%CONSTANTES DE MARGULES
A12=0.2762;
A21=0.2877;
lamda1=10.^((1-x1).^2.*(A12+2*x1*(A21-A12)))
lamda2=10.^(x1.^2.*(A21+2*(1-x1)*(A12-A21)))
x=[0:0.01:1];
y=[0:0.01:1];
%GRAFICANDO
hold on
plot(x1,y1,'r')
plot(x,y,'g')
legend('datos experimentales','MARGULES','DIAGONAL')
grid
hold on
grid
1.5
MARGULES
1
DIAGONAL
y1
0.5
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1
340
335
T Vs x1
T °K
T Vs y1
330
325
0 0.2 0.4 0.6 0.8 1 1.2 1.4
x1 y1
6. PROGRAMACION CON EL METODO WILSON EN UN SISTEMA ISOBARICO:
disp('*********PROCESO ISOBARICO EQUILIBRIO ACETONA-METANOL***********')

clc, clear all
x1e=[1 0.93 0.867 0.819 0.783 0.773 0.735 0.66 0.594 0.554 0.519 0.407
0.394 0.369 0.324 0.281 0.263 0.229 0.195 0.151 0.1 0.074 0.053 0.024 0];
x2e=1-x1e;
y1e=[1 0.918 0.856 0.812 0.782 0.774 0.745 0.689 0.644 0.618 0.594 0.514
0.504 0.485 0.45 0.41 0.392 0.357 0.319 0.265 0.191 0.148 0.11 0.053 0];
y2e=1-y1e;
Te=[56.53 55.67 55.43 55.31 55.24 55.25 55.3 55.39 55.74 55.96 56.15 56.9
57.05 57.29 57.69 58.32 58.49 58.97 59.57 60.41 61.53 62.21 62.79 63.69
64.7];
Te1=273.15+Te
%Presion constante a 800 mmHg
P=800;
R=1.987; %cal/mol*K
VPA1=-7.45514;
VPB1=1.202;
VPC1=-2.43926;
VPD1=-3.36690;
Pc1=47;
Tc1=508.1;
%DATOS DEL METANOL
VPA2=-8.54796;
VPB2=0.76982;
VPC2=-3.1085;
VPD2=1.54481;
Pc2=80.9;
Tc2=512.6;
TR1=1-Te1./Tc1;
TR2=1-Te1./Tc2;
%PROGRAMACION CON EL METODO DE WILSON
x1=[1 0.93 0.867 0.819 0.783 0.773 0.735 0.66 0.594 0.554 0.519 0.407
0.394 0.369 0.324 0.281 0.263 0.229 0.195 0.151 0.1 0.074 0.053 0.024 0];
%CONSTANTES DE WILSON
lamba12=-214.95; %cal/mol-g
lamba21=664.08; %cal/mol-g
v2=42.4809;
v1=77.3603;
A12=v2/v1*exp(-(lamba12)./(R*Te1));
A21=v1/v2*exp(-(lamba21)./(R*Te1));
lamda1=exp(-log(x1+A12.*(1-x1))+(1-x1).*((A12./(x1+A12.*(1-x1)))-(A21./
(A21.*x1+(1-x1)))));
lamda2=exp(-log((1-x1)+A21.*x1)-x1.*((A12./(x1+A12.*(1-x1)))-(A21./
(A21.*x1+(1-x1)))));
x=[0:0.01:1];
y=[0:0.01:1];
%GRAFICANDO
hold on
plot(x1,y1,'r')
plot(x,y,'g')
legend('datos experimentales','WILSON','DIAGONAL')
grid
hold on
grid
1.5
WILSON
1 DIAGONAL
y1
0.5
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1
340
T Vs x1
T °K
T Vs y1
330
325
0 0.2 0.4 0.6 0.8 1 1.2 1.4
x1 y1
GRACIAS!!!!
clc, clear all
syms T
% ideal isotermico
% Equilibrio Liquido Vapor - BENCENCENO (1) - CICLOPENTANO (2)
% Datos de fraccion molar
x1e=[0;0.2;0.4;0.6;0.8;1];
% Parte gaseosa de la acetona
y2e=1-x1e;
% Presiones variables ya que es un oroceso isotermico
Pc1=48.9;
Pc2=45.1;
% Temperatura constante en kelvin
Tc1=572.3;
Tc2=511.3;
%Datos encontrados en tablas termodinamicas
VPA1=-6.38273;
VPB1=1.31233;
VPC1=3.62863;
VPD1=-1.15399;
% Presion y temperatura critica de la acetona
VPA2=-6.5180;
VPB2=0.33442;
VPC2=-1.11706;
VPD2=-4.50035;
% Calculo de las presiones de vapor
TR1=1-T/Tc1;
TR2=1-T/Tc2;
% Aplicacion de la ecuacion de van laar
x1=[0:0.01:1];
% Constantes de van laar conseguido en tablas
A12=0.1655;
A21=0.1302;
lamda1=10.^(A12./(1+(x1*A12)./((1-x1)*A21)).^2)
lamda2=10.^(A21./(1+((1-x1)*A21)./(x1*A12)).^2)
x=[0:0.01:1];
y=[0:0.01:1];
% Grafica
subplot(2,1,1), plot(x1e,y1,'+')
hold on
plot(x1,y1,'g')
plot(x,y)
grid
hold on
grid

Metodo Margules, Van Lar y Wilson

Cargado por

Copyright:

Formatos disponibles

Metodo Margules, Van Lar y Wilson

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Metodo Margules, Van Lar y Wilson

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD NACIONAL DEL CENTRO DEL PERU

FACULTAD DE INGENIERIA QUIMICA

APELLIDOS Y NOMBRES: HURTADO TOCAS ANYELO YONEL

1. PROGRAMACION CON EL METODO DE VAN LAAR:

clc, clear all

ISOTERMICO: Equilibrio Liquido Vapor - Acetona - Metanol

clc, clear all

3. PROGRAMACION CON LA ECUACION DE WILSON:

clc, clear all

4. PROGRAMACION CON EL METODO VAAN LAR EN UN SISTEMA ISOBARICO:

disp('***********PROGRAMA SISTEMA ISOBARICO (ACETONA-

0.5 datos experimentales

5. PROGRAMACION CON EL METODO MARGULES EN UN SISTEMA ISOBARICO:

disp('*PROCESO ISOBARICO EQUILIBRIO ACETONA-METANOL***')

6. PROGRAMACION CON EL METODO WILSON EN UN SISTEMA ISOBARICO:

disp('*PROCESO ISOBARICO EQUILIBRIO ACETONA-METANOL***')

También podría gustarte