EQUILIBRIO EN LAS REACCIONES QUÍMICAS

EQUILIBRIO EN LAS REACCIONES QUÍMICAS

ESTUDIANTE: Baltazar Alania, Abigail Cristina                                                                           VI-IQ

CATEDRA: Dr. Salvador Teódulo Oré Vidalón

CATEDRATICO:  Termodinámica de los Procesos Químicos

FECHA DE ENTREGA: 12-01-2021 

% PROGRAMA PARA CALCULAR EL COEFICIENTE DE ACTIVIDAD MEDIANTE UNIFAC

% PARA LA MEZCLA ETANOL - AGUA A UNA TEMPERATURA DE 89°C Y FRACCION MOLAR

% DE X1=0.2

function UNIFAC

clc

clear all

disp('COEFICIENTE DE ACTIVIDAD PARA EL SISTEMA ')

%DATOS PROPUESTOS DEL PROBLEMA

%Fraccion molar

x(1,1)=0.2;

x(2,1)=1-x(1);

%Temperatura

T=83; %°C

T=T+273;%K

z=10; %número de coordinacion

NSG=4; %numero de subgrupo

NC=2; %numero de componentes

%PARAMETROS DE LOS GRUPOS FUNCIONALES SACADAS DE LA REF.1

R=[0.9011 0.6744 1.0000 0.92000];%Parametros de volumen de los subgrupos

Q=[0.848 0.540 1.200 1.400];%Parametros de area de los subgrupos

vk1=[1;1;1;0];%Numeros de grupos funcionales de los grupos en el etanol

vk2=[0;0;0;1];%Numeros de grupos funcionales de los grupos en el agua

a=[0 0 986.5 1318; 0 0 986.5 1318; 156.4 156.4 0 353.5; 300 300 -229.1 0];%Parametros de interaccion binaria

%ECUACIONES EN LA PARTE COMBINATORIAL SACADAS DE LA REF. 2

%CALCULANDO r

for i=1:NSG

    r(1,i)=vk1(i).*(R(i).');

    r(2,i)=vk2(i).*(R(i).');

end

r=[sum(r(1,:));sum(r(2,:))];%Parametro de volumen de los componentes puros

%CALCULANDO q

for i=1:NSG

    q(1,i)=vk1(i).*(Q(i).');

    q(2,i)=vk2(i).*(Q(i).');

end

q=[sum(q(1,:));sum(q(2,:))];%Parametros de area de los componentes puro

%CALCULANDO FI

for j=1:NC

    FI(j,1)=(r(j)*x(j))/sum(r.*x);%Fraccion del segmento de los componentes

end

%CALCULANDO THETA

for j=1:NC

    THETA(j,1)=(q(j)*x(j))/sum(q.*x);%Fraccion de area de los componentes

end

%CALCULANDO L

for i=1:NC

    L(i,1)=(z/2)*(r(i)-q(i))-(r(i)-1); %Definida por la ecuacion

end

%REEMPLAZANDO EN LA ECUACION DE LA PARTE COMBINATORIAL

for i=1:NC

lngammaC(i,1)=log(FI(i)/x(i))+(z/2)*q(i)*log(THETA(i)/FI(i))+L(i)-(FI(i)/x(i))*sum(x.*L);%Coefcientes de actividad combinatorial del componente

end

%ECUACIONES EN LA PARTE RESIDUAL SACADAS DE LA REF. 2

%CALCULO DE lnTAUK

vm=[vk1 vk2];%Numero de grupos funcionales en los componentes

for m=1:NSG

    for j=1:NC

    X(m,j)=x(j)*vm(m,j);

    end

    X1=X(:,1)+X(:,2);

end

Xm=X1/sum(X1);%Fraccion mol de grupo m en la mezcla

for m=1:NSG

    for n=1:NSG

        NI(n)=(Xm(n)*Q(n).');

    end

    THETA(m,1)=(Xm(m)*Q(m).')/sum(NI);%Fraccion de area del grupo m

end

PSI=exp(-a/T);

%REEMPLAZANDO EN LA ECUACION DE lnTAUK

for k=1:NSG

    for m=1:NSG

        A1(m,1)=THETA(m).*PSI(m,k);

        for n=1:NSG

            B12(n,1)=THETA(n).*PSI(n,m);

...