ex8_12 - Ca=Cao/(1+tau*k1); kappa=UA/(vo*Cao)/Cp;

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

View Full Document Right Arrow Icon
% "ex8_12" % m-file to solve example 8.12 % % x(1)=T % xdot(1)=dT/dt % function xdot=ex8_12(t,x) f global tau DH1 DH2 kappa Tc Cp E1 E2 g Cp=200; Cao=0.3; To=283; tau=.01; DH1=-55000; DH2=-71500; vo=1000; E2=27000; E1=9900; UA=40000; Ta=330; k2=4.58*exp((E2/1.987)*(1/500-1/x(1))); k1=3.3*exp((E1/1.987)*(1/300-1/x(1)));
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: Ca=Cao/(1+tau*k1); kappa=UA/(vo*Cao)/Cp; G=-tau*k1/(1+k1*tau)*DH1-k1*tau*k2*tau*DH2/((1+tau*k1)*(1+tau*k2)); Tc=(To+kappa*Ta)/(1+kappa); Cb=tau*k1*Ca/(1+k2*tau); R=Cp*(1+kappa)*(x(1)-Tc); Cc=Cao-Ca-Cb; F=G-R; F xdot(1,:)=2;...
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