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 CAb es la concentración del compuesto, z es el largo del lecho, DAB, es la dispersión, U es la velocidad superficial, rA es la velocidad de consumo de y rhob 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.
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 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.