# ex9_4 - cbo=Fbo/vo; cmo=Fmo/vo; Na=x(1)*V; Nb=x(2)*V;...

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

% "ex9_4" % m-file to solve example 9.4 % % 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=ex9_4(t,x) f UA=16000; Ta1=60; Fao=80; To=75; k=16.96e12*exp(-32400/1.987/(x(5)+460)); Fbo=1000; Fmo=100; ra=-k*x(1); rb=-k*x(1); rc=k*x(1); Thetacp=35+Fbo/Fao*18+Fmo/Fao*19.5; vo=Fao/0.923+Fbo/3.45+Fmo/1.54; V=500/7.484; tau=V/vo; cao=Fao/vo;
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: cbo=Fbo/vo; cmo=Fmo/vo; Na=x(1)*V; Nb=x(2)*V; Nc=x(3)*V; Nm=x(4)*V; mc=1000; NCp=Na*35+Nb*18+Nc*46+Nm*19.5; Ta2=x(5)-((x(5)-Ta1)*exp(-UA/(18*mc))); Q=mc*18*(Ta1-Ta2); Q 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,:)=(Q-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