ex9_3 - cp=75240; cpa=170700; T0=300; dh=-7.9076e7;...

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

View Full Document Right Arrow Icon
% "ex9_3" % m-file to solve example 9.3 % % x(1)=ca; x(2)=cb; x(3)=cc=cd; x(4)=Nw; x(5)=T % xdot(1)=d(ca)/dt, etc. % function xdot=ex9_3(t,x) f Kc=10^(3885.44/x(5)); k=0.39175*exp(5472.7*((1/273)-(1/x(5)))); ra=-k*((x(1)*x(2))-((x(3)*x(3))/Kc)); v0=0.004; Vi=0.2; V=Vi+v0*t; rb=ra; rc=-ra; cb0=1; cw0=55; UA=3000; Ta=290; Fb0=cb0*v0;
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: cp=75240; cpa=170700; T0=300; dh=-7.9076e7; Na=x(1)*V; Nb=x(2)*V; Nc=x(3)*V; Nd=x(3)*V; thetaw=55; t xdot(1,:)=ra-(v0*x(1)/V); xdot(2,:)=rb+(v0*(cb0-x(2))/V); xdot(3,:)=rc-(x(3)*v0)/V; xdot(4,:)=v0*cw0; xdot(5,:)=(UA*(Ta-x(5))-Fb0*cp*(1+thetaw)*(x(5)-T0)+ra*V*dh)/(cp*(Nb+Nc+Nd+x(4)) +cpa*Na); +...
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