Seguridad Industiral

Difusividad másica del metanol en H2O
Para hallar la difusividad del metanol en agua variando la temperatura se desarrollo el siguiente código en Matlab (R2008b), el cual arroja los seis valores dedifusividad pedidos para una temperatura fija utilizando el método de Hsu-Chen.
1. Método de Hsu-Chen
clear all

xa=[0+eps,0.25:0.25:0.95,0.90,1+eps];
Ra=1.4311;
Rb=0.9200;
Qa=1.432;Qb=1.400;

%Paso 1 H-C

landaa= (Ra^(1/3));
landab= (Rb^(1/3));

for i=1:length(xa)

phia(i)=((xa(i)*landaa)/((xa(i)*landaa)+((1-xa(i))*landab)));phib(i)=(((1-xa(i))*landab)/((xa(i)*landaa)+((1-xa(i))*landab)));

end

for j=1:length(xa)

tethaa(j)=((xa(j)*Qa)/((xa(j)*Qa)+((1-xa(j))*Qb)));
tethab(j)=(((1-xa(j))*Qb)/((xa(j)*Qa)+((1-xa(j))*Qb)));
end

%Paso2 H-C

Aab=-10.7575;
Aba=194.5302;
T=313.13;

tauab= exp((-Aab)/T);
tauba= exp((-Aba)/T);
%Paso 3 H-S

tethaaa=((tethaa)./(tethaa+(tethab.*tauba)));tethaba=((tethab.*tauba)./(tethaa+(tethab.*tauba)));
tethabb=((tethab)./(tethab+(tethaa.*tauab)));
tethaab=((tethaa.*tauab)./(tethab+(tethaa.*tauab)));


%Paso 4 Difusividad diluida
Estas son difusividades diluidas medidasexperimentalmente.

Dabo=2.10e-5;
Dbao=2.67e-5;


%Paso 5 (Continuación del método H-C)

DDab= ((1-xa).*log(Dabo)) + (xa.*log(Dbao)) +(2.*[(xa.*log(xa./phia))+((1-xa).*log((1-xa)./phib))])+...
(2.*xa.*(1-xa).*[((phia./xa).*(1-(landaa/landab)))+((phib./(1-xa)).*(1-(landab/landaa)))])+...
((1-xa).*Qa.*[((1-((tethaba).^2)).*log(tauba))+((1-((tethabb).^2)).*tauab.*log(tauab))])+...xa.*Qb.*[((1-((tethaab).^2)).*log(tauab))+ ((1-((tethaaa).^2)).*tauba.*log(tauba))];

Dab=exp(DDab)'


Este procedimiento se repite para las otras dos temperaturas cambiando el valor de la variableT y el valor de las difusividades diluidas (Dabo, Dbao).
Al ejecutar el anterior código, se obtuvo como resultado los valores de difusividad de la Tabla1 a diferente temperatura.

Tabla 1. DAB...