cde9_1 -...

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

View Full Document Right Arrow Icon
% "cde9_1" % m-file to solve CD example 9.1 % % x(1)=ca % x(2)=cb % x(3)=cc % x(4)=cm % x(5)=T % xdot(1)=d(ca)/dt, xdot(2)=d(cb)/dt, xdot(3)=d(cc)/dt, etc. % function xdot=cde9_1(t,x) f UA=70000; Ta=85; Fao=400; To=75; k=16.96e12*exp(-32400/1.987/(x(5)+460)); Fbo=5000; Fco=0; Fmo=20; ra=-k*x(1); rb=-k*x(1); rc=k*x(1); Thetacp=35+Fbo/Fao*18+Fmo/Fao*19.5; vol=Fao*58.08/(.859/454*2.54^3*12^3)+Fbo*18.02/(.9941/454*2.54^3*12^3)+Fco*76.11/
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: (1.036/454*2.54^3*12^3)+Fmo*32.04/(.7914/454*2.54^3*12^3); tau=.1; V=tau*vol; cao=Fao/vol; cbo=Fbo/vol; cmo=Fmo/vol; Na=x(1)*V; Nb=x(2)*V; Nc=x(3)*V; Nm=x(4)*V; NCp=Na*35+Nb*18+Nc*46+Nm*19.5; N xdot(1,:)=1/tau*(cao-x(1))+ra; xdot(2,:)=1/tau*(cbo-x(2))+rb; xdot(3,:)=1/tau*(0-x(3))+rc; xdot(4,:)=1/tau*(cmo-x(4)); xdot(5,:)=(UA*(Ta-x(5))-Fao*Thetacp*(x(5)-To)+(-36000)*ra*V)/NCp;...
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