diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /249/CH21 | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '249/CH21')
| -rwxr-xr-x | 249/CH21/EX21.1/21_01.sce | 39 | ||||
| -rwxr-xr-x | 249/CH21/EX21.2/21_02.sce | 39 |
2 files changed, 78 insertions, 0 deletions
diff --git a/249/CH21/EX21.1/21_01.sce b/249/CH21/EX21.1/21_01.sce new file mode 100755 index 000000000..71fe097be --- /dev/null +++ b/249/CH21/EX21.1/21_01.sce @@ -0,0 +1,39 @@ +clear
+clc
+t=[0;2;4;6];
+XA=[0.75;0.64;0.52;0.39];
+t1=4000;//kg.s/m3
+density_s=1500;//kg/m3
+De=5*10^-10;
+d=2.4*10^-3;
+//Assuming -rA=kCA*a,-da/dt=kd*a
+//For this rate a plot of ln(CAo/CA-1)vs t should give a straight line
+for i=1:4
+ y(i)=log((1/(1-XA(i)))-1);
+end
+plot(y,t)
+xlabel('t')
+ylabel('ln(CAo/CA-1)')
+//Guessing No Intrusion of Diffusional Resistance
+//ln(CAo/CA-1)=ln(k*t1)-kd*t
+coeff =regress(t,y);
+kd=coeff(2);
+k=exp(coeff(1))/t1;
+L=d/6;
+Mt=L*sqrt(k*density_s/De);
+//Assuming Runs were made in regime of strong resistance to pore diffusion
+k1=((exp(coeff(1)))^2)*(L^2)*density_s/(t1*t1*De);
+kd1=-2*coeff(2);
+Mt=L*sqrt(k1*density_s/De);
+printf("\n Rate equation(mol/kg.s) in diffusion free regime with deactivation is %f ",k1)
+printf("CA*a with \n -da/dt(hr-1) is %f",kd1)
+printf("a")
+//In strong pore diffusion
+k2=k1*sqrt(De/(k1*density_s));
+printf("\n Rate equation(mol/kg.s) in strong pore diffusion resistance regime with deactivation is %f ",k2)
+printf("CA*a^0.5/L with \n -da/dt(hr-1) is %f",kd1)
+printf("a")
+
+
+
+
diff --git a/249/CH21/EX21.2/21_02.sce b/249/CH21/EX21.2/21_02.sce new file mode 100755 index 000000000..74b961c15 --- /dev/null +++ b/249/CH21/EX21.2/21_02.sce @@ -0,0 +1,39 @@ +clear
+clc
+PAo=3;//atm
+R=82.06*10^-6;//m3.atm/mol.k
+T=730;//k
+W=1000;//kg
+FAo=5000;//mol/hr
+CAo=PAo/(R*T);
+tau=W*CAo/FAo;
+i=0;
+for t=0:5:120
+ i=i+1;
+ //Part a
+a(i)=1-(8.3125*10^-3)*t;
+XA(i)=(tau^2)*a(i)/(1+(tau^2)*a(i));
+//Part b
+a1(i)=exp(-0.05*t);
+XA1(i)=(tau^2)*a1(i)/(1+(tau^2)*a1(i));
+//Part c
+a2(i)=1/(1+3.325*t);
+XA2(i)=(tau^2)*a2(i)/(1+(tau^2)*a2(i));
+//Part d
+a3(i)=1/(sqrt(1+1333*t));
+XA3(i)=(tau^2)*a3(i)/(1+(tau^2)*a3(i));
+end
+t=[0:5:120];
+plot(t,XA,t,XA1,t,XA2,t,XA3)
+xlabel('Time(days)')
+ylabel('XA')
+legend('Zero Order','1st Order','2nd Order','3rd Order');
+XA_avg=(1/120)*integrate('(100*(1-(8.3125*10^-3)*t))/(1+100*(1-(8.3125*10^-3)*t))','t',0,120);
+XA1_avg=(1/120)*integrate('(100*exp(-0.05*t))/(1+100*exp(-0.05*t))','t',0,120);
+XA2_avg=(1/120)*integrate('(100*(1/(1+3.325*t)))/(1+100*(1/(1+3.325*t)))','t',0,120);
+XA3_avg=(1/120)*integrate('(100*1/(sqrt(1+1333*t)))/(1+100*(1/sqrt(1+1333*t)))','t',0,120);
+printf("\n for d=0,the mean conversion is % f",XA_avg)
+printf("\n for d=1,the mean conversion is % f",XA1_avg)
+printf("\n for d=2,the mean conversion is % f",XA2_avg)
+printf("\n for d=3,the mean conversion is % f",XA3_avg)
+
|