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

---
 249/CH11/EX11.1/11_01.sce | 26 ++++++++++++++++++++++++++
 249/CH11/EX11.2/11_02.sce | 33 +++++++++++++++++++++++++++++++++
 249/CH11/EX11.3/11_03.sce | 23 +++++++++++++++++++++++
 249/CH11/EX11.4/11_04.sce | 17 +++++++++++++++++
 249/CH11/EX11.5/11_05.sce | 11 +++++++++++
 5 files changed, 110 insertions(+)
 create mode 100755 249/CH11/EX11.1/11_01.sce
 create mode 100755 249/CH11/EX11.2/11_02.sce
 create mode 100755 249/CH11/EX11.3/11_03.sce
 create mode 100755 249/CH11/EX11.4/11_04.sce
 create mode 100755 249/CH11/EX11.5/11_05.sce

(limited to '249/CH11')

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)
-- 
cgit