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/CH11 | |
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/CH11')
-rwxr-xr-x | 249/CH11/EX11.1/11_01.sce | 26 | ||||
-rwxr-xr-x | 249/CH11/EX11.2/11_02.sce | 33 | ||||
-rwxr-xr-x | 249/CH11/EX11.3/11_03.sce | 23 | ||||
-rwxr-xr-x | 249/CH11/EX11.4/11_04.sce | 17 | ||||
-rwxr-xr-x | 249/CH11/EX11.5/11_05.sce | 11 |
5 files changed, 110 insertions, 0 deletions
diff --git a/249/CH11/EX11.1/11_01.sce b/249/CH11/EX11.1/11_01.sce new file mode 100755 index 000000000..24554079f --- /dev/null +++ b/249/CH11/EX11.1/11_01.sce @@ -0,0 +1,26 @@ +clear
+clc
+//Given Time(min) and Tracer Output Concentration(g/litre)
+T=[0;5;10;15;20;25;30;35];
+Cpulse=[0;3;5;5;4;2;1;0];
+dt=5;
+//Mean Residence time(t)
+sum1=0;
+sum2=0;
+Area=0;//Initialization
+for i=1:8
+sum1=sum1+T(i)*Cpulse(i)*dt;
+sum2=sum2+Cpulse(i)*dt;
+//Area Under Concentration-Time Curve
+Area=Area+Cpulse(i)*dt;
+end
+t=sum1/sum2;
+printf("\n The mean residence time(min) is %f \n",t)
+for j=1:8
+ E(j)=Cpulse(j)/Area;
+end
+plot(T,E)
+xlabel('time(min)')
+ylabel('E(min^-1)')
+title('Exit age distribution E vs time')
+
diff --git a/249/CH11/EX11.2/11_02.sce b/249/CH11/EX11.2/11_02.sce new file mode 100755 index 000000000..37859180a --- /dev/null +++ b/249/CH11/EX11.2/11_02.sce @@ -0,0 +1,33 @@ +clear
+clc
+M=150;//Molecular mass(gm)
+v=5;//litre/sec
+v=5*60;//litre/min
+V=860;//litres
+//From Material Balance
+Area1=M/v;//gm.min/litre
+//From the tracer curve
+A1=0.375;
+Area2=A1*(1+1/4+1/16+1/64+1/256+1/1024+1/4096);//Taking Significant Areas
+printf("\n From material balance Area(gm.min/litre) is %f",Area1)
+printf("\n From Tracer Curve Area(gm.min/litre) is %f",Area2)
+printf("\n Part a")
+printf("\n As the two areas are equal,this is a properly done experiment \n")
+//For the liquid,calculating t
+sum1=0;
+for i=1:10
+ sum1=sum1+2*i*A1/(4^(i-1));
+ t=sum1/Area1;
+end
+//liquid volume in vessel
+Vl=t*v;
+//Fraction of liquid
+f=Vl/V;
+printf("\n Part b")
+printf("\n Fraction of liquid is %f",f)
+//E=Cpulse/M/v
+printf("\n Part c")
+printf("\n The E curve is 1.5C")
+printf("\n Part d")
+printf("\n The vessel has a strong recirculation of liquid,probably induced by the rising bubbles")
+
diff --git a/249/CH11/EX11.3/11_03.sce b/249/CH11/EX11.3/11_03.sce new file mode 100755 index 000000000..d611752b8 --- /dev/null +++ b/249/CH11/EX11.3/11_03.sce @@ -0,0 +1,23 @@ +clear
+clc
+//From the given graph
+Cin(1)=0;Cin(2)=8;Cin(3)=4;Cin(4)=6;Cin(5)=0;
+E(5)=0;E(6)=0.05;E(7)=0.5;E(8)=0.35;E(9)=0.1;E(10)=0;
+for t=8:14
+sum1=0;
+for p=5:t-1
+ if p>10 |(t-p)>5
+ h=2;
+ else
+sum1=sum1+Cin(t-p)*E(p);
+Cout(t)=sum1;
+end
+end
+end
+t=[1:1:14];
+Cout=Cout';
+t
+Cout
+plot(t,Cout)
+xlabel('t');ylabel('Cout')
+
diff --git a/249/CH11/EX11.4/11_04.sce b/249/CH11/EX11.4/11_04.sce new file mode 100755 index 000000000..05c117fc1 --- /dev/null +++ b/249/CH11/EX11.4/11_04.sce @@ -0,0 +1,17 @@ +clear
+clc
+k=0.307;
+//Given mean residence time(min)
+t=15;
+//For plug flow with negligible density
+fr_unconverted=exp(-k*t);
+printf("\n The fraction of reactant unconverted in a plug flow reactor is %f",fr_unconverted)
+//For the real reactor
+T=[5;10;15;20;25;30];//given time
+E=[0.03;0.05;0.05;0.04;0.02;0.01];//given
+dt=5;
+sum1=0;
+for i=1:6
+ sum1=sum1+exp(-k*T(i))*E(i)*dt;
+end
+printf("\n The fraction of reactant unconverted in a real reactor is %f",sum1)
\ No newline at end of file diff --git a/249/CH11/EX11.5/11_05.sce b/249/CH11/EX11.5/11_05.sce new file mode 100755 index 000000000..cf613db22 --- /dev/null +++ b/249/CH11/EX11.5/11_05.sce @@ -0,0 +1,11 @@ +clear
+clc
+k=0.5;//litre/mol.min
+CAo=2;//mol/litre
+//From the batch eqn
+//CA/CAo=1/(1+kCAo*t)
+to=1;t1=3;
+E=0.5;
+//Using eqn 13
+XA_avg=1-(E*integrate('1/(1+t)','t',to,t1));
+printf("\n Average concentration of A remaining in the droplet is %f",XA_avg)
|