Jump to content



Featured Articles

Check out the latest featured articles.

File Library

Check out the latest downloads available in the File Library.

New Article

Product Viscosity vs. Shear

Featured File

Vertical Tank Selection

New Blog Entry

Low Flow in Pipes- posted in Ankur's blog

Diseño De Un Reactor Tubular De Lecho Empacado

reactor lecho empacado packed bed

This topic has been archived. This means that you cannot reply to this topic.
No replies to this topic
Share this topic:
| More

#1 csandoval

csandoval

    Brand New Member

  • Members
  • 2 posts

Posted 28 October 2014 - 05:21 PM

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:

 

KKJkDWw.png

 

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.

 

% 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.

 






Similar Topics