Programa Para Determinar Temperatura de Burbuja

3
clear all %Calcular la Temperatura de Burbuja (Tb) para una mezcla de sustancias puras a una presión y fracciones molares conocidas %De la ecuación de Antoine A1=input('Introduzca el valor de la constante A para la sustancia 1') B1=input('Introduzca el valor de la constante B para la sustancia 1') C1=input('Introduzca el valor de la constante C para la sustancia 1') A2=input('Introduzca el valor de la constante A para la sustancia 2') B2=input('Introduzca el valor de la constante B para la sustancia 2') C2=input('Introduzca el valor de la constante C para la sustancia 2') %Presión del sistema y fracciones molares de los componentes P= input('Introduzca el valor de la presión del sistema en mm de Hg: '); X1=input('Introduzca el valor de la fracción molar del componente 1') X2=1-X1 %Determinación de la temperatura de ebullición normal Tem1=B1./(A1-log10(P))-C1; Tem2=B2./(A2-log10(P))-C2; Tem= [Tem1 Tem2] Min=min(Tem); Max=max(Tem); fprintf('La temperatura de tanteo debe estar entre T1=%g y T2=%g \ n',Min,Max); T= input('Temperatura inicial de tanteo en °C; ');

description

Programa de matlab para determinar la temperatura de burbuja

Transcript of Programa Para Determinar Temperatura de Burbuja

Page 1: Programa Para Determinar Temperatura de Burbuja

clear all

%Calcular la Temperatura de Burbuja (Tb) para una mezcla de sustancias puras a una presión y fracciones molares conocidas

%De la ecuación de Antoine

A1=input('Introduzca el valor de la constante A para la sustancia 1')

B1=input('Introduzca el valor de la constante B para la sustancia 1')

C1=input('Introduzca el valor de la constante C para la sustancia 1')

A2=input('Introduzca el valor de la constante A para la sustancia 2')

B2=input('Introduzca el valor de la constante B para la sustancia 2')

C2=input('Introduzca el valor de la constante C para la sustancia 2')

%Presión del sistema y fracciones molares de los componentes

P= input('Introduzca el valor de la presión del sistema en mm de Hg: ');

X1=input('Introduzca el valor de la fracción molar del componente 1')

X2=1-X1

%Determinación de la temperatura de ebullición normal

Tem1=B1./(A1-log10(P))-C1;

Tem2=B2./(A2-log10(P))-C2;

Tem= [Tem1 Tem2]

Min=min(Tem); Max=max(Tem);

fprintf('La temperatura de tanteo debe estar entre T1=%g y T2=%g \n',Min,Max);

T= input('Temperatura inicial de tanteo en °C; ');

% Presiones de vapor

for l=1:1000

P1o(l)=10.^(A1-(B1./(T+C1)));

P2o(l)= 10.^(A1-(B1./(T+C1)));

% Determinación de Y1 y Y2

Y1(l)=P1o(l)*X1/P; Y2(l)=P2o(l)*X2/P;

% Condición de Y1 + Y2 = 1

sumaY(l)= Y1(l)+Y2(l);

Page 2: Programa Para Determinar Temperatura de Burbuja

condicion(l)= abs(sumaY(l)-1);

if condicion(l)<= 0.01

Tb=T ; Iteracion=l;

break

else

T=T+0.1;

end

end

fprintf('A la presión de %g mm de Hg\n',P);

fprintf('y a las concentraciones de X1=%g y X2=%g\n',X1,X2);

fprintf('La temperatura de burbuja es %g °C\n',Tb);

fprintf('El número de iteraciones realizados para el cálculo es=%gn',Iteracion);