No replies to this topic

[csandoval]

Hola, saludos a todos. Soy estudiante de Ingeniería Civil Química y tengo una consulta sobre el diseño de un reactor tubular de lecho empacado para catálisis heterogénea en estado estacionario en Matlab. Estoy utilizando el siguiente modelo:

 

 

Donde C_Ab es la concentración del compuesto, z es el largo del lecho, D_AB, es la dispersión, U es la velocidad superficial, r_A es la velocidad de consumo de y rho_b es la densidad del lecho catalítico.

 

El asunto es que e revisado y aún no encuentro en qué me e equivocado y bastante tiempo e perdido en esto. Sería de mucha ayuda si alguien pudiera revisar lo que tengo y dar algún consejo.

 

% script

clear all; clc
%butanotiol+heptano
%T=25ºC


Deff=0.046E-8        %m2/s
ro_b=167.9400;      %kg/m3
L=4.27;                   %m
U=0.0024/0.327;     %m/s
k=4.68E-2                %m6/(mol*s*kg_cat)
Ca0=6.05;                %mol/m3


%ODE
z=0:0.000001:L;
z1=z';
y0=[Ca0; 0]
ode_react = @(z,y)react(z,y,U,Deff,ro_b,k);
[z,y]=ode15s(ode_react,z1,y0)
plot(z,y(:,1))
ylabel('Concentración Ca [mol/m3]')
xlabel('z   [m]')

 

function

 

function dy=react(z,y,U,Deff,ro_b,k)
dy = zeros(2,1);
dy(1)=y(2);
dy(2)=(U./Deff).*y(2)-k.*ro_b.*(y(1).^2)./Deff;
end

Muchas gracias de antemano.