ex10_7 - Pb=Pao*x(1)/(1+eps*x(1)); vo=Uo*3.1416*D*D/4;...

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

View Full Document Right Arrow Icon
% "ex10_7" % m-file to solve example 10.7 % % x(1)=X % xdot(1)=d(X)/dz % function xdot=ex10_7(z,x) f global A g Ka=0.05; Kb=.15; Pao=12; eps=1; A=7.6; R=0.082; T=400+273; rho=80; kprime=0.0014; D=1.5; Uo=2.5; Kc=0.1; U=Uo*(1+eps*x(1)); Pa=Pao*(1-x(1))/(1+eps*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: Pb=Pao*x(1)/(1+eps*x(1)); vo=Uo*3.1416*D*D/4; Cao=Pao/R/T; Kca=Ka*R*T; Pc=Pb; a=1/(1+A*(z/U)^0.5); raprime=a*(-kprime*Pa/(1+Ka*Pa+Kb*Pb+Kc*Pc)); ra=rho*raprime; r xdot(1,:)=-ra/U/Cao;...
View Full Document

Ask a homework question - tutors are online