INSTITUTO TECNOLOGICO DE ESTUDIOS SUPERIORES DE MONTERREY 2o Parcial de Termodinámica del Equilibrio - B Nombre del alumn@ :Zaira Mendez Ramirez Fecha: 19 de Junio Todos los programas se encuentran al final 1. Para una mezcla binaria de CCl4 y 1-propanol, determinar el diagrama Pxy a 75°C a partir de la ecuación de Roult ideal y Roult modificada con Wilson. Explicar su comportamiento. (20pts) Constantes de Antoine CCL4 (compuesto 1 ) A= 4.10445 4.10445 B=1265.632 C=232.148 1-propanol (compuesto 2 ) A=4.99991 A=4.99991 B=1512.940 C=205.807 Modelo Ideal Equilibrio liquido-vapor del Compuesto 2 (Metodo Ideal) Equilibrio liquido-vapor del Compuesto 1 (Metodo Ideal)
1
1
Composicion de liquido
Composicion de liquido
Composicion de vapor
Composicion de vapor vapor
0.9
0.9
0.8 C
0.8 C
n n
e e
ar ut
ra ut
0.7 ra
0.7 ra
e e
p p m
m e
e T
T
0.6
0.5
0.5
0.4
0.6
0.4
0
0.1
0.2
0.3 0. 4 0. 5 0.6 0.7 Composicion de Compuesto 2
0.8
0. 9
1
0
0. 1
0. 2
0.3 0.4 0. 5 0. 6 0.7 Composicion de Compuesto 1
0. 8
0. 9
1
El comportamiento del modelo ideal para esta mezcla muestra que a temperatura constante la línea de líquido saturado (línea azul) se encuentra por arriba de la del vapor saturado, las gráficas de ambos compuestos solo difieren en la pendiente, para el compuesto 2(1-propanol) presenta pendiente negativa y para el compuesto uno (CCl4) presenta pendiente positiva. Este modelo ideal
no presenta las desviaciones reales que la mezcla puede tener, y de usa fórmula Roult ideal, la cual no toma en cuenta el coeficiente de actividad.
Modelo Wilson Parametros de interaccion binaria l11p =-.6856000000 l21p = 5.019200000 l12p l22p v1p v2p
Equilibrio liquido-vapor del Compuesto 2 (Metodo Wils on) 1 Composicion de liquido Composicion de vapor
= 110.2379000 = -2413.50710 = 0.0959883 = 0.0744555
0.9 Equilibrio liquido-vapor del Compuesto 1 (Metodo Wilson) 1
0.8 C
Composicion de liquido
n
Composicion de vapor e ar tu
0.9
0.7 ar e p m
0.8
e T
C n
0.6 e ra ut
0.7 ar e p
0.5 m e T
0.4
0
0.1
0.2
0.3 0.4 0.5 0.6 0.7 Composicion de Compuesto 2
0.8
0.9
0.6
1 0.5
0.4
0
0.1
0.2
0.3 0.4 0.5 0.6 0.7 Composicion de Compuesto 1
0.8
0.9
1
Los parámetros de interacción binaria se sacaron en el programa aspen,para poder hacer el cálculo de los coeficientes de actividad, que se usan en la fórmula de Roult para sacar las compopsiciones de vapor de los compuestos. Es notable una desviación al modelo ideal anteriormente presentado. Para el compuesto 2(1-propanol) cambia el diagrama en la línea de vapor saturado al acercarse más a la línea de líquido satrurado cuando la composición es menor, en cambio en el diagrama del compuesto 1 (CCl4) cambia la línea del vapor saturado al alejarse de la dell líquido saturado cuando la concentración va aumentando.
2. Los siguientes datos experimentales de equilibrio líquido-vapor corresponden al sistema de acetato de etilo-etanol a una presión de 1atm. El sistema forma un azeótropo con una composición de 54%mol de acetato de etilo, a partir de la ecuación de estado de ideal y NRTL, compruebe el punto azeotrópico realizando los cálculos en Matlab. Grafique los 2 sistemas para analizar sus diferencias. (20pts)
Constantes de Antoine Acetato de etilo (compuesto 1) A=4.13361 B=1195.130 C=212.470 Tvap con Antoine = 77.048°C Etanol A=5.33675 B= 1648.220 C=230.918 Tvap con Antoine = 78.250°C
Ideal
Equilibrio liquido-vapor del Compuesto 2 (Metodo Ideal)
Equilibrio liquido-vapor del Compuesto 1 (Metodo Ideal)
78.4
78.4 Composicion de liquido Composicion de vapor
78.2
78.2
78
78
C n e
C n e 77.8 a r u t a r e p 77.6 m e T
77.8 ar tu ar e p
77.6 m e T
77.4
77.4
77.2
77
Composicion de liquido Composicion de vapor
77.2
0
0.1
0.2
0.3 0.4 0.5 0.6 Composicion de Compuesto 2
0.7
0.8
0.9
77 0.1
0.2
0.3
0.4 0.5 0.6 0.7 Composicion de Compuesto 1
0.8
0.9
1
Se puede observar que los diagramas salen las líneas de vapor y de líquido muy cercanas entre ellas , sin embargo si se agranda la imagen se ve que hay un punto donde se llegan a juntar por la composición 0.5, aunque se tenga un comportamiento extraño si llega a ha verse el punto azeotropico.
NRTL Parámetros de interacción binaria g11 = g21 = g12 = g22 = alfa=
input('Aij:');-0.2431 input('Aji:');-1.1512 input('Bij:');282.9558 input('Bji:');524.4238 input('Cij:');0.3 Equilibrio liquido-vapor del Compuesto 2 (Metodo NRTL)
78.4
Equilibrio liquido-vapor del Compuesto 1 (Metodo NRTL) 78.4
Composicion de liquido
Composicion de liquido Composicion de vapor
Composicion de vapor
78.2
78.2
78
78
e
n
C
C n e 77.8 a r u t a r e p
77.8 ar tu p
e
ar
77.6 m
m e T
e T
77.6
77.4
77.4
77.2
77.2
77
0
0.1
0.2
0.3 0.4 0.5 0.6 Composicion de Compuesto 2
0.7
0.8
0.9
77 0.1
0.2
0.3
0.4 0.5 0.6 0.7 Composicion de Compuesto 1
0.8
0.9
1
Los diagramas con el metodo NRTL usando los parámetros de interacción bianria proporcionados por Aspen son similares a los del método ideal, su forma extraña se debe a que sólo pocos valores se están graficando ya que se tomaron las temperaturas de ebullición del compuesto uno y del compuesto dos para el rango de temperaturas y no los experimentales ya que se esta haciendo la comparación.
3. A partir del ley de Roult ideal y utilizando la ecuación de Harcheler en vez de la de Antoine, realice los siguientes cálculos para el sistema benceno(1)/tolueno(2). (20pts) a) Considere y1 = 0.33 y T = 100°C, encuentre x1 y P b) Considere x1 = 0.33 y P = 120KPa, encuentre y1 y T Compuesto/constante de Harcheler [=] mmHg / K Benceno Tolueno
A
B
C
D
52.100 -5557.61 -5.072 3.610 56.785 -6283.50 -5.681 4.840
En hoja aparte
4. 4. Defina que es equilibrio líquido vapor El equilibrio líquido vapor es un estado de un compuesto o de una mezcla donde se encuentra en transición de líquido a vapor y viceversa y no se puede definir un estado particular, ese es el estado líquido vapor, este estado está definido como el espacio en el diagrama liquidovapor que se encuentra entre las líneas de vapor saturado y líquido saturado. 5. Defina el punto de burbuja Es el punto o la temperatura cuando comienza a cambiar de líquido a vapor y se dice punto de burbuja porque como cuando hierve el agua, es la temperatura a la cual se forma la primera burbuja. 6. Defina el punto de rocío Es un punto a una temperatura donde comienza a cambiar de fase de vapor a fase líquida y se comienza a condensar. 7. Defina presión de vapor Es la presión a la cual un compuesto comienza a entra a la fase de vapor o se comienza a evaporar en cierta condición de temperatura. Se calcula con la ecuación de Antoine 8. Defina que es un azeótropo Un azeótropo se define como una mezcla de dos o más elementos en donde se tiene un comportamiento igual y parece compuesto puro, es muy difícil separarlos ya que tienen la misma temperatura de ebullición, es necesario separarlos por destilación azeotrópica.
EJERCICIO 1 RESULTADOS MODELO IDEAL Elige 1 si el sistema esta a Temperatura constante y 2 si el sistema esta a Presion constante 1 Temperatura del sistema = 75 Introduzca las constantes de Antoine del compuesto 1 A:4.10445 B:1265.632 C:232.148 Introduzca las constantes de Antoine del compusto 2 A:4.99991 B:1512.940 C:205.807 o Presion de Vapor del compuesto 1= 0.964 bar o Presion de Vapor del compuesto 2= 0.409 bar Calculo de Presion en el compuesto 1 en fase liquida Pa= 0.409 Pa= 0.465 Pa= 0.520 Pa= 0.576 Pa= 0.631 Pa= 0.686 Pa= 0.742 Pa= 0.797 Pa= 0.853 Pa= 0.908 Pa= 0.964
Calculo de Presion en el compuesto 2 en fase liquida Pb= 0.409 bar
Pb= 0.465 bar Pb= 0.520 bar Pb= 0.576 bar Pb= 0.631 bar Pb= 0.686 bar Pb= 0.742 bar Pb= 0.797 bar Pb= 0.853 bar Pb= 0.908 bar Pb= 0.964 bar
- Coeficientes de actividad(Gamma)y composiciones en vapor (y): Compuesto 1 0.000000
1.000000
0.207316
0.792684
0.370459
0.629541
0.502187
0.497813
0.610777
0.389223
0.701834
0.298166
0.779286
0.220714
0.845971
0.154029
0.903988
0.096012
0.954923
0.045077
1.000000
0.000000
Compuesto 2
RESULTADOS MODELO WILSON Elige 1 si el sistema esta a Temperatura constante y 2 si el sistema esta a Presion constante1 Temperatura del sistema = 75 Introduzca las constantes de Antoine del compuesto 1 A:4.10445 B:1265.632 C:232.148 Introduzca las constantes de Antoine del compusto 2 A:4.99991
B:1512.940 C:205.807 o Presion de Vapor del compuesto 1= 0.964 bar o Presion de Vapor del compuesto 2= 0.409 bar Calculo de Presion en el compuesto 1 en fase liquida Pa= 0.409 Pa= 0.465 Pa= 0.520 Pa= 0.576 Pa= 0.631 Pa= 0.686 Pa= 0.742 Pa= 0.797 Pa= 0.853 Pa= 0.908 Pa= 0.964
Calculo de Presion en el compuesto 2 en fase liquida Pb= 0.409 bar Pb= 0.465 bar Pb= 0.520 bar Pb= 0.576 bar Pb= 0.631 bar Pb= 0.686 bar Pb= 0.742 bar Pb= 0.797 bar Pb= 0.853 bar Pb= 0.908 bar Pb= 0.964 bar
Introduzca los parametros de interaccion binaria Aij:-.6856000000 Aji:5.019200000
Bij:110.2379000 Bji:-2413.50710 v1:0.0959883 v2:0.0744555 - Parametros de interaccion binaria (Delta):
Delta12 1.241
Delta21 0.336
- Coeficientes de actividad(Gamma)y composiciones en vapor (y): Compuesto 1 Gamma
y
Compuesto 2 Gamma
y
1.565299
0.000000
1.000000
1.000000
1.488225
0.308533
1.002721
0.794841
1.413430
0.523617
1.011969
0.637077
1.341318
0.673592
1.029905
0.512701
1.272403
0.777155
1.059748
0.412478
1.207361
0.847366
1.106505
0.329923
1.147127
0.893940
1.178384
0.260086
1.093085
0.924718
1.289733
0.198657
1.047428
0.946862
1.467772
0.140924
1.013952
0.968246
1.770048
0.079788
1.000000
1.000000
2.337836
0.000000
>>
EJERCICIO 1. Programa en Matlab modelo Ideal X=input('Elige 1 si el sistema esta a Temperatura constante y 2 si el sistema esta a Presion constante'); if (X==1) T1=input('Temperatura del sistema = '); %PRESIONES DE VAPOR A TCTE disp('Introduzca las constantes de Antoine del compuesto 1') A1=input('A:'); B1=input('B:'); C1=input('C:'); disp('Introduzca las constantes de Antoine del compusto 2') A2=input('A:'); B2=input('B:'); C2=input('C:'); %Definir operacion, ecuacion de Antoine
P1= 10^(A1-(B1/(C1+T1))); fprintf(' o Presion de Vapor del compuesto 1= %0.3f bar \n',P1); P2= 10^(A2-(B2/(C2+T1))); fprintf(' o Presion de Vapor del compuesto 2= %0.3f bar \n',P2); %CALCULO DE LAs PRESIONES VARIABLS DEL SISTEMA SEGUN ECUACION P-Pvap2 =(Pvap1-Pvap2)x disp('Calculo de Presion en el compuesto 1 en fase liquida') x =0:0.1:1; for i=1:length(x) %Definir operacion para sacar la Presion del sistema PS1(i)=((P1-P2)*(x(i)))+ P2; fprintf('Pa= %0.3f \n',PS1(i)); end disp(' ') disp('Calculo de Presion en el compuesto 2 en fase liquida') x=0:0.1:1; for i=1:length(x) %Definir operacion para sacar la Presion del sistema PS2(i)= (((P2-P1)*(1-x(i)))+ P1); fprintf('Pb= %0.3f bar \n',PS2(i)); end %Calculo de composiciones de vapor (y) disp(' ') disp('- Coeficientes de actividad(Gamma)y composiciones en vapor (y):') disp(' Compuesto 1 Compuesto 2') disp('') for i=1:length(x) y1(i)= ( x(i)*P1)/ PS1(i) ; y2(i)= ( (1-x(i))*P2)/ PS2(i) ; fprintf(' %f %f xi(i)=1-x(i); end
\n',y1(i),y2(i));
%Graficar Metodo de Ideal Temperatura Constante figure(1) plot(x,PS1,y1,PS1),xlabel('Composicion de Compuesto 1'),ylabel('Temperatura en C'),title('Equilibrio liquido-vapor del Compuesto 1 (Metodo Ideal)'),legend('Composicion de liquido','Composicion de vapor') figure(2) plot(xi,PS2,y2,PS2),xlabel('Composicion de Compuesto 2'),ylabel('Temperatura en C'),title('Equilibrio liquido-vapor del Compuesto 2 (Metodo Ideal)'),legend('Composicion de liquido','Composicion de vapor') end
EJERCICIO 1. Programa en Matlab modelo Wilson X=input('Elige 1 si el sistema esta a Temperatura constante y 2 si el sistema esta a Presion constante'); if (X==1) T1=input('Temperatura del sistema = '); %PRESIONES DE VAPOR A TCTE disp('Introduzca las constantes de Antoine del compuesto 1') A1=input('A:'); B1=input('B:'); C1=input('C:'); disp('Introduzca las constantes de Antoine del compusto 2') A2=input('A:'); B2=input('B:'); C2=input('C:');
%Definir operacion, ecuacion de Antoine P1= 10^(A1-(B1/(C1+T1))); fprintf(' o Presion de Vapor del compuesto 1= %0.3f bar \n',P1); P2= 10^(A2-(B2/(C2+T1))); fprintf(' o Presion de Vapor del compuesto 2= %0.3f bar \n',P2); %CALCULO DE LAs PRESIONES VARIABLS DEL SISTEMA SEGUN ECUACION P-Pvap2 =(Pvap1-Pvap2)x disp('Calculo de Presion en el compuesto 1 en fase liquida') x =0:0.1:1; for i=1:length(x) %Definir operacion para sacar la Presion del sistema PS1(i)=((P1-P2)*(x(i)))+ P2; fprintf('Pa= %0.3f \n',PS1(i)); end disp(' ') disp('Calculo de Presion en el compuesto 2 en fase liquida') x=0:0.1:1; for i=1:length(x) %Definir operacion para sacar la Presion del sistema PS2(i)= (((P2-P1)*(1-x(i)))+ P1); fprintf('Pb= %0.3f bar \n',PS2(i)); end %CALCULO DE LOS COEFICINTES DE ACTIVIDAD disp(' ') disp('Introduzca los parametros de interaccion binaria') l11 = input('Aij:'); l21 = input('Aji:'); l12 = input('Bij:'); l22 = input('Bji:'); v1 = input('v1:'); v2 = input('v2:'); R = 8.314; disp('- Parametros de interaccion binaria (Delta):') disp(' ') disp(' Delta12 Delta21') %Calculo para los parametros de interaccion binaria. A12 = (v1/v2)*exp(- ((l12-l11)/(R*(T1+273)))); A21 = (v2/v1)*exp(- ((l21-l22)/(R*(T1+273)))); fprintf(' %0.3f %0.3f \n',A12,A21); %Calculo para los coeficientes de actividad for i=1:length(x) %Caculo de los coeficientes de actividad (Gamma) gamma1(i) = exp( -log(x(i)+(A12*(1-x(i)))) + (1-x(i))*((A12/(x(i)+(A12* (1-x(i))))) (A21/((1-x(i))+(A21*x(i))))) ); gamma2(i) = exp( -log((1-x(i))+(A21*x(i))) - x(i)*((A12/(x(i)+(A12*(1- x(i))))) (A21/((1-x(i))+(A21*x(i))))) ); end %Calculo de composiciones de vapor (y) disp(' ') disp('- Coeficientes de actividad(Gamma)y composiciones en vapor (y):') disp(' Compuesto 1 Compuesto 2') disp(' Gamma y Gamma y') disp('') for i=1:length(x) y1(i)= ( gamma1(i)*x(i)*P1)/ PS1(i) ; y2(i)= ( gamma2(i)*(1-x(i))*P2)/ PS2(i) ; fprintf(' %f %f %f %f \n',gamma1(i),y1(i),gamma2(i),y2(i)); xi(i)=1-x(i); end
%Graficar Metodo de Wilson Temperatura Constante figure(1) plot(x,PS1,y1,PS1),xlabel('Composicion de Compuesto 1'),ylabel('Temperatura en C'),title('Equilibrio liquido-vapor del Compuesto 1 (Metodo Wilson)'),legend('Composicion de liquido','Composicion de vapor') figure(2) plot(xi,PS2,y2,PS2),xlabel('Composicion de Compuesto 2'),ylabel('Temperatura en C'),title('Equilibrio liquido-vapor del Compuesto 2 (Metodo Wilson)'),legend('Composicion de liquido','Composicion de vapor') end
EJERCICIO 2 RESULTADO MODELO IDEAL Elige 1 si el sistema esta a Temperatura constante y 2 si el sistema esta a Presion constante2 Presion del sistema = 1.013 Introduzca las temperaturas de vapor del compuesto uno y del compuesto dos Tvap1:77.048 Tvap2:78.250 Introduzca las constantes de Antoine del compuesto 1 A:4.13361 B:1195.130 C:212.470 Introduzca las constantes de Antoine del compusto 2 A:5.33675 B:1648.220 C:230.918 1.013 bar
0.966 bar
1.047 bar
1.005 bar
- Composicion en liquido :
Compuesto 1
Compuesto 2
1.000
0.000
Compuesto 1 0.194
Compuesto 2
0.806
- Coeficientes de actividad(Gamma)y composiciones en vapor (y): Compuesto 1 y
Compuesto 2 y
0.999925
0.000075
0.200486
0.799514
RESULTADO MODELO NRTL Elige 1 si el sistema esta a Temperatura constante y 2 si el sistema esta a Presion constante2 Presion del sistema = 1.013 Introduzca las temperaturas de vapor del compuesto uno y del compuesto dos Tvap1:77.048 Tvap2:78.250 Introduzca las constantes de Antoine del compuesto 1 A:4.13361 B:1195.130 C:212.470 Introduzca las constantes de Antoine del compusto 2 A:5.33675 B:1648.220 C:230.918 1.013 bar
0.966 bar
1.047 bar
1.005 bar
- Composicion en liquido :
Compuesto 1
Compuesto 2
1.000
0.000
0.194
0.806
Aij:-0.2431 Aji:-1.1512 Bij:282.9558 Bji:524.5438 Cij:0.3 -0.083011
-0.000312
1.025216
0.999906
-0.082775
-0.000311
1.025143
0.999907
- Coeficientes de actividad(Gamma)y composiciones en vapor (y): Compuesto 1 Gamma
y
Compuesto 2 Gamma
y
1.000000
0.999925
0.920065
0.000069
0.946685
0.189797
0.996750
0.796916
PROGRAMA IDEAL if (X==2) P=input('Presion del sistema = '); disp('Introduzca las temperaturas de vapor del compuesto uno y del compuesto dos') T2=input('Tvap1:'); T3=input('Tvap2:'); T=T2:1:T3; %Defnir rango de Temperatura %PRESIONES DE VAPOR A PCTE disp('Introduzca las constantes de Antoine del compuesto 1') A3=input('A:'); B3=input('B:'); C3=input('C:'); disp('Introduzca las constantes de Antoine del compusto 2') A4=input('A:'); B4=input('B:'); C4=input('C:'); %Definir operacion, ecuacion de Antoine for i=1:length(T) Pvap1(i)= 10^(A3-(B3/(C3+T(i)))); Pvap2(i)= 10^(A4-(B4/(C4+T(i)))); fprintf(' %0.3f bar %0.3f bar \n',Pvap1(i),Pvap2(i));
end %COMPOSICIONES EN LIQUIDO disp(' ') disp('- Composicion en liquido :') for i=1:length(T) %Definir operacion para sacar la composicion del liquido x1(i)= ((P - (Pvap2(i)) )/((Pvap1(i))-(Pvap2(i)))); disp(' ') disp(' Compuesto 1 Compuesto 2') x2(i)= ((P - (Pvap1(i)) )/((Pvap2(i))-(Pvap1(i)))); fprintf(' %0.3f %0.3f \n',x1(i),x2(i)); end %CALCULO DE COMPOSICIONS DE VAPOR(y)(Raoult no ideal) disp(' ') disp('- Coeficientes de actividad(Gamma)y composiciones en vapor (y):') disp(' Compuesto 1 Compuesto 2') disp(' y y') disp('') for i=1:length(T) y3(i)= ( x1(i)*Pvap1(i) )/ P ; y4(i)= ( x2(i)*Pvap2(i) )/ P ; fprintf(' %f %f \n',y3(i),y4(i)); end %Graficar Metodo de Wilson Pcte figure(3) plot(x1,T,y3,T),xlabel('Composicion de Compuesto 1'),ylabel('Temperatura en C'),title('Equilibrio liquido-vapor del Compuesto 1 (Metodo Ideal)'),legend('Composicion de liquido','Composicion de vapor') figure(4) plot(x2,T,y4,T),xlabel('Composicion de Compuesto 2'),ylabel('Temperatura en C'),title('Equilibrio liquido-vapor del Compuesto 2 (Metodo Ideal)'),legend('Composicion de liquido','Composicion de vapor') end
PROGRAMA NRTL if (X==2) P=input('Presion del sistema = '); disp('Introduzca las temperaturas de vapor del compuesto uno y del compuesto dos') T2=input('Tvap1:'); T3=input('Tvap2:'); T=T2:1:T3; %Defnir rango de Temperatura %PRESIONES DE VAPOR A PCTE disp('Introduzca las constantes de Antoine del compuesto 1') A3=input('A:'); B3=input('B:'); C3=input('C:'); disp('Introduzca las constantes de Antoine del compusto 2') A4=input('A:'); B4=input('B:'); C4=input('C:'); %Definir operacion, ecuacion de Antoine for i=1:length(T) Pvap1(i)= 10^(A3-(B3/(C3+T(i)))); Pvap2(i)= 10^(A4-(B4/(C4+T(i)))); fprintf(' %0.3f bar %0.3f bar \n',Pvap1(i),Pvap2(i)); end %COMPOSICIONES EN LIQUIDO disp(' ') disp('- Composicion en liquido :')
for i=1:length(T) %Definir operacion para sacar la composicion del liquido x1(i)= ((P - (Pvap2(i)) )/((Pvap1(i))-(Pvap2(i)))); end disp(' ') disp(' Compuesto 1 Compuesto 2') for i=1:length(T) P= 1.013; %Definir operacion para sacar la composicion del liquido x2(i)= ((P - (Pvap1(i)) )/((Pvap2(i))-(Pvap1(i)))); fprintf(' %0.3f %0.3f \n',x1(i),x2(i)); end %COEFICIENTES DE ACTIVIDAD g11 = input('Aij:'); g21 = input('Aji:'); g12 = input('Bij:'); g22 = input('Bji:'); alfa= input('Cij:'); R = 8.314; for i=1:length(T) %Calculo para los parametros de interaccion binaria t12(i) = (g12-g22)/(R*(T(i)+273)); t21(i) = (g21-g11)/(R*(T(i)+273)); G12(i) = exp(-alfa*t12(i)); G21(i) = exp( alfa*t21(i)); fprintf(' %f %f %f %f \n',t12(i),t21(i),G12(i),G21(i)); end for i=1:length(T) %Caculo de los coeficientes de actividad (Gamma) gamma3(i) = exp((x2(i)^2) * ( ( t21(i)*((G21(i)/(x1(i)+(x2(i)*G21(i))))^2) ) + ( (t12(i)*G12(i))/((x2(i)+(x1(i)*G12(i)))^2)))); gamma4(i) = exp((x1(i)^2) * ( ( t12(i)*((G12(i)/(x2(i)+(x1(i)*G12(i))))^2) ) + ( (t21(i)*G21(i))/((x1(i)+(x2(i)*G21(i)))^2)))); end %CALCULO DE COMPOSICIONS DE VAPOR(y)(Raoult no ideal) disp(' ') disp('- Coeficientes de actividad(Gamma)y composiciones en vapor (y):') disp(' Compuesto 1 Compuesto 2') disp(' Gamma y Gamma y') disp('') for i=1:length(T) y3(i)= ( gamma3(i)*x1(i)*Pvap1(i) )/ P ; y4(i)= ( gamma4(i)*x2(i)*Pvap2(i) )/ P ; fprintf(' %f %f %f %f \n',gamma3(i),y3(i),gamma4(i),y4(i)); end %Graficar Metodo de NRTL Pcte figure(3) plot(x1,T,y3,T),xlabel('Composicion de Compuesto 1'),ylabel('Temperatura en C'),title('Equilibrio liquido-vapor del Compuesto 1 (Metodo NRTL)'),legend('Composicion de liquido','Composicion de vapor') figure(4) plot(x2,T,y4,T),xlabel('Composicion de Compuesto 2'),ylabel('Temperatura en C'),title('Equilibrio liquido-vapor del Compuesto 2 (Metodo NRTL)'),legend('Composicion de liquido','Composicion de vapor') end
t12(i) = a12+(b12/T(i)); t21(i) = a21+(b21/T(i)); g11=-3.447 g21=2.849 g12=1726.751 g22=-423.960
a12=-3.447 a21=2.849 b12=1726.751 b21=-423.960
a12=3.8679 a21=-0.8029 b12=-632.2853 b21=746.2919
t12(i) = (g12-g22)/(R*(T(i))); t21(i) = (g21-g11)/(R*(T(i)));
x11=[0.1399 0.1445] x21=1-x11 x12=[0.9376 0.9236] x22=1-x12 T=[298.57 308.34] alfa=0.3 x11=[0.1399 0.1445] x21=1-x11 x12=[0.9376 0.9236] x22=1-x12 T=[298.57 308.34] alfa=0.3
%%%% a12=-1.42916695702586 a21=-10.662 b12=1227.02048770515 b21=3920.27780496373
