From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 800/CH8/EX8.10/8_10.sce | 40 ++++++++++++++++++++++++++++++++++++++++ 1 file changed, 40 insertions(+) create mode 100755 800/CH8/EX8.10/8_10.sce (limited to '800/CH8/EX8.10/8_10.sce') diff --git a/800/CH8/EX8.10/8_10.sce b/800/CH8/EX8.10/8_10.sce new file mode 100755 index 000000000..2be7eaad9 --- /dev/null +++ b/800/CH8/EX8.10/8_10.sce @@ -0,0 +1,40 @@ +//clear// +clc +clear +//eY(2)ec("8.6data.sci"); +W = 0:1:28.58; +W0=0; +function w=f(W,Y) + w =zeros(3,1); + + +fao=.188 +visc=.090 +Ta=1264.67 +deltah=-42471-1.563*(Y(3)-1260)+.00136*(Y(3)**2-1260**2)-(2.459*10**(-7))*(Y(3)**3-1260**3); +summ= 57.23+.014 * Y(3)-1.94 *10**(-6.)*Y(3)**2 +dcp=-1.5625+2.72*10**(-3)*Y(3)-7.38*10**(-7)*Y(3)**2 +k=360D*exp(-176008/Y(3)-(110.1*log(Y(3)))+912.8) +thetaso=0; +Po=2 +Pao=.22 +thetao=.91 +eps=-.055 +R=1.987; +Kp=exp(42311/R/Y(3)-11.24); +if(Y(2)< =.05) + + ra=(-k*(.848-.012/(Kp**2))); +else + ra=(-k*(1-Y(2))/(thetaso+Y(2)))**.5*(Y(1)/Po*Pao*((thetao-.5*Y(2))/((1+eps*Y(2)))-((thetaso+Y(2))/(1-Y(2)))**2/(Kp**2))); +end + +w(1)=(-1.12*10**(-8)*(1-.055*Y(2))*Y(3))*(5500*visc+2288)/Y(1) ; +w(2)=-(ra)/fao ; +w(3)=(5.11*(Ta-Y(3))+(-ra)*(-deltah) )/(fao*(summ+Y(2)*dcp)) +endfunction + +X=ode([2;0;1400],W0,W,f); + +plot2d(W,X(1,:)); +plot2d(W,X(3,:)); -- cgit