ex4_8 - V=pi*(z*R^2-1/3*(z-L)^3-1/3*L^3); G=m/Ac;

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
% "ex4_8" % m-file to solve example 4.8 % % x(1)=X % x(2)=y % xdot(1)=dX/dz, xdot(2)=dy/dz % function xdot=ex4_8(z,x) f Fa0=440; P0=2000; Ca0=0.32; R=30; phi=0.4; kprime=0.02; L=27; rhocat=2.6; m=44; Ca=Ca0*(1-x(1))*x(2)/(1+x(1)); Ac=pi*(R^2-(z-L)^2);
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: V=pi*(z*R^2-1/3*(z-L)^3-1/3*L^3); G=m/Ac; ra=-kprime*Ca*rhocat*(1-phi); beta=(98.87*G+25630*G^2)*0.01; W=rhocat*(1-phi)*V; W xdot(1,:)=-ra*Ac/Fa0; xdot(2,:)=-beta/P0/x(2)*(1+x(1));...
View Full Document

This note was uploaded on 09/02/2011 for the course CHEM 102 taught by Professor Kmn during the Spring '11 term at Babson College.

Ask a homework question - tutors are online