ex9_2 - Qr=UA*(x(1)-298); Theata=Nbo/Nao;

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

View Full Document Right Arrow Icon
% "ex9_2" % m-file to solve example 9.2 % % x(1)=T % x(2)=X % xdot(1)=d(T)/dt, xdot(2)=d(X)/dt, etc % function xdot=ex9_2(t,x) global TT g NCp=2504; V=3.265+1.854; Nao=9.0448; UA=35.83; dH=-590000; then=0; Nbo=33; TT=x(1)-273; k=.00017*exp(11273/(1.987)*(1/461-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: Qr=UA*(x(1)-298); Theata=Nbo/Nao; ra=-k*(Nao^2)*(1-x(2))*(Theata-2*x(2))/V^2; rate=-ra; Qg=ra*V*(dH); Q xdot(1,:)=(Qg-Qr)/NCp; xdot(2,:)=(-ra)*V/Nao; x...
View Full Document

Ask a homework question - tutors are online