Prediction of Conversion with an Isothermal Axially-Dispersed Plug Flow Reactor: first_order.m - File Exchange

No BSD License

Highlights from
Prediction of Conversion with an Isothermal Axially-Dispersed Plug Flow Reactor

f=shoot1(N)
f=shoot2(N)
ydot=goveqn1(t,y)
ydot=goveqn2(t,y)
first_order.m
second_order.m
View all files

from Prediction of Conversion with an Isothermal Axially-Dispersed Plug Flow Reactor by Housam Binous
Prediction of the conversion of an isothermal axially dispersed plug flow reactor.

first_order.m

global Pe Da N


Di=0.005;
k=0.01;
L=1;
u=0.01;

Pe=L*u/Di;
Da=k*L/u;

Pe 
Da

N=fsolve('shoot1',[0.1 -0.1],optimset('TolFun',1e-3,'Display','iter'));

[t y]=ode45('goveqn1',[0 1],N);

1-y(end,1)
plot(t,y(:,1),'g')

Highlights from Prediction of Conversion with an Isothermal Axially-Dispersed Plug Flow Reactor

Highlights from
Prediction of Conversion with an Isothermal Axially-Dispersed Plug Flow Reactor