absorber

9/05/11 12:20 PM

C:\Users\Martha\Documents\MATLAB\absorber.m

function dy=absorber(z,y)
dy=zeros(6,1);

%Crea la matriz de la variacion de y (variables)
%con z (altura de la torre)
%Davalores iniciales de cero a la matriz de 6 columnas en
%y (variables a calcular)y una en z (altura de la torre)

%Datos de variables conocidas:
XAGB=.17;
%Concentración de A en el gas que entra(Bottom)
XALB=.026; %Concentración de A en el agua que sale (Bottom) calculado del
%50% de saturación
XAGT=0.01; %Concentración de A en el gas que sale (Top)
NGB=26.6;
%Flux de gas que entra (Bottom)NLB=-164;
%Flux de líquido que sale (Bottom). El signo negativo indica
%que va en la dirección contraria e z positivo (baja)
DHV=55.6;
%Entalpía de vaporización del agua
DHS=1.3E4; %Entalpía dedisolución del amoníaco en el agua
CPG=7.0;
%Cp del gas
CPL=18;
%Cp del líquido
PTOT=1.0;
%Presión total del sistema
TB=35;
%Temperatura en el fondo de la torre (Bottom)
KOGAI=19;
%KOG*aI%Parámetros del equilibrio de amoníaco:
A1=1.6E9; A2=6276;
%Parámetros de la ecuación de Antoine para calcular la presión de vapor del
%agua:
B1=8.10765; B2=1750.286; B3=235;
NXAGB=4.522;
%Fluxde A en el gas de entrada
NXALB=-4.28;
%Flux de A en el líquido a la salida
NXBGB=.146;
%Flux de B en el gas de entrada (10% de humedad relativa)
%Definiciones:
%A partir de las 6 ecuacionesdiferenciales de la tabla 13-2:
%y(1)=NG*XAG
%y(2)=NL*XAL
%y(3)=T
El resto del factor que acompaña a T se pasa a dividir
%y(4)=NG*XGB
%y(5)=NG
%y(6)=y(1)/y(5) Esta se define convenientemente parahacer explícito el
%valor de XAG que es el que nos interesa.
%Fórmulas:
XAG=y(1)/y(5);
y(6)=XAG;
NL=-y(5)+NGB+NLB;

MA=A1*exp(-A2/(y(3)+273));
XAL=y(2)/NL;
XAGS=MA*XAL;NAI=KOGAI*(XAG-XAGS);
ALPHA=B1-B2/(B3+y(3));
PBST=10^ALPHA;

%Hace explícito XAG
%Hace explicito XAG
%Última ecuación de la tabla 13-2, no requiere
%ser definida como Y(i) porque no es una
%ecuación diferencial...