7 files changed, 175 insertions, 0 deletions
diff --git a/249/CH18/EX18.1/18_01.sce b/249/CH18/EX18.1/18_01.sce
new file mode 100755
index 000000000..6e102ae11
--- /dev/null
+++ b/249/CH18/EX18.1/18_01.sce
@@ -0,0 +1,53 @@
+clear
+clc
+dp=2.4*(10^-3);L=dp/6;
+//Effective mass conductivity(m3/hr.mcat)
+De=5*10^-5;
+//Effective thermal conductivity(KJ/hr.mcat.K)
+Keff=1.6;
+//For the gas film surrounding the pellet
+h=160;//heat transfer coefficient(KJ/hr.m2cat.K)
+kg=300;//mass transfer coefficient(m3/hr.m2cat)
+//For the reaction
+Hr=-160;//KJ/molA
+CAg=20;//mol/m3
+rA_obs=10^5;//mol/hr.m3cat
+kobs=rA_obs/CAg;
+Vp=3.14*(dp^3)/6;
+S=3.14*(dp^2);
+//Observed rate/rate if film resistance controls
+ratio=kobs*Vp/(kg*S);
+printf("\n Part a")
+if ratio<0.01 
+    printf("\n Resistance to mass transport to film should not influence rate of reaction")
+else
+      printf("\n Resistance to mass transport to film should influence rate of reaction")
+end
+printf("\n Part b")
+
+Mw=rA_obs*(L^2)/(De*CAg);
+printf("\n Mw= %f",Mw)
+if Mw>4
+  printf("\n Pore diffusion is influencing and hence strong pore diffusion")
+else
+      printf("\n Pore diffusion is  not influencing and hence weak pore diffuusion")
+end
+//Temp variation within pellet
+dt_max_pellet=De*(CAg-0)*(-Hr)/Keff;
+//Temp variation Across the gas film
+dt_max_film=L*rA_obs*(-Hr)/h;
+printf("\n Part c")
+printf("\n dTmax,pellet is %f",dt_max_pellet)
+printf(" degree C \n dTmax,film is %f",dt_max_film)
+printf(" degree C")
+if dt_max_pellet<1
+    printf("\n Pellet is close to uniform in temperature")
+else
+    printf("\n There is a variation in temp within pellet")
+    end
+    if dt_max_film<1
+    printf("\n Film is close to uniform in temperature")
+else
+    printf("\n There is a variation in temp within Film")
+
+end
+\ No newline at end of file
diff --git a/249/CH18/EX18.2/18_02.sce b/249/CH18/EX18.2/18_02.sce
new file mode 100755
index 000000000..2725ac3f5
--- /dev/null
+++ b/249/CH18/EX18.2/18_02.sce
@@ -0,0 +1,30 @@
+clear
+clc
+//Pressure(atm)
+PAo=3.2;
+R=0.082;//litre.atm/mol.k
+T=390;//k
+v=20;//litre/hr
+W=0.01;///kg
+CA_in=[0.1;0.08;0.06;0.04];
+CA_out=[0.084;0.07;0.055;0.038];
+CAo=PAo/(R*T);
+FAo=CAo*v;
+eA=3;
+for i=1:4
+XA_in(i)=(1-CA_in(i)/CAo)/(1+eA*CA_in(i)/CAo);
+XA_out(i)=(1-CA_out(i)/CAo)/(1+eA*CA_out(i)/CAo);
+dXA(i)=XA_out(i)-XA_in(i);
+rA(i)=dXA(i)/(W/FAo);
+CA_avg(i)=(CA_in(i)+CA_out(i))/2;
+end
+plot(CA_avg,rA)
+xlabel('CA(mol/litre)')
+ylabel('-rA(mol/hr.kg)')
+coeff1=regress(CA_avg,rA)
+k=coeff1(2)
+printf("\n The rate of reaction(mol/hr.kg) is %f",k)
+printf("CA")
+disp('The answer slightly differs from those given in book as regress fn is used for calculating slope and intercept')
+
+
diff --git a/249/CH18/EX18.3/18_03.sce b/249/CH18/EX18.3/18_03.sce
new file mode 100755
index 000000000..a281ee41e
--- /dev/null
+++ b/249/CH18/EX18.3/18_03.sce
@@ -0,0 +1,33 @@
+clear
+clc
+CAo=0.1;//mol/litre
+FAo=2;//mol/hr
+eA=3;
+CA=[0.074;0.06;0.044;0.029];//mol/litre
+W=[0.02;0.04;0.08;0.16];//kg
+//Gussing 1st order,plug flow rxn
+//(1+eA)*log(1/(1-XA))-eA*XA=k*(CAo*W/FAo)
+for i=1:4
+XA(i)=(CAo-CA(i))/(CAo+eA*CA(i));
+y(i)=(1+eA)*log(1/(1-XA(i)))-eA*XA(i);
+x(i)=CAo*W(i)/FAo;
+W_by_FAo(i)=W(i)/FAo;
+CAout_by_CAo(i)=CA(i)/CAo;
+XA1(i)=(1-CAout_by_CAo(i))/(1+eA*CAout_by_CAo(i));
+end
+plot(x,y)
+coeff3=regress(x,y);
+xlabel('CAoW/FAo'),ylabel('4ln(1/1-XA)-3XA')
+k=coeff3(2);
+printf("\n Part a,using integral method of analysis")
+printf("\n The rate of reaction(mol/litre) is %f",k)
+printf("CA")
+//Part b
+//plotting W_by_FAo vs XA1,the calculating rA=dXA/d(W/FAo) for last 3 points,we get
+rA=[5.62;4.13;2.715];
+coeff2=regress(CA(2:4),rA);
+printf("\n Part b,using differential method of analysis")
+printf("\n The rate of reaction(mol/litre) is %f",coeff2(2))
+printf("CA")
+
+
diff --git a/249/CH18/EX18.4/18_04.sce b/249/CH18/EX18.4/18_04.sce
new file mode 100755
index 000000000..30b7da7f6
--- /dev/null
+++ b/249/CH18/EX18.4/18_04.sce
@@ -0,0 +1,8 @@
+clear
+clc
+XA=0.35;
+FAo=2000;//mol/hr
+eA=3;k=96;
+CAo=0.1;
+W=((1+eA)*log(1/(1-XA))-eA*XA)*(FAo/(k*CAo));
+printf("\n The amount of catalyst(kg) needed in a packed bed reactor is %f",W)
diff --git a/249/CH18/EX18.5/18_05.sce b/249/CH18/EX18.5/18_05.sce
new file mode 100755
index 000000000..deab73fbb
--- /dev/null
+++ b/249/CH18/EX18.5/18_05.sce
@@ -0,0 +1,37 @@
+clear
+clc
+CAo=0.1;
+eA=3;
+rA=[3.4;5.4;7.6;9.1];
+CA=[0.039;0.0575;0.075;0.092];
+XA=zeros(4,1);
+inv_rA=zeros(4,1);
+for i=1:4
+XA(i)=(1-CA(i)/CAo)/(1+eA*CA(i)/CAo);
+inv_rA(i)=1/rA(i);
+end
+//W=FAo*integral(dXA/-rA) from 0 to 0.35
+//Using Trapezoidal rule to find area,XA must be in increasing order
+//Sorting XA and accordingly inv_rA
+for i=1:4
+    small=XA(i);
+    for j=i:4
+        next=XA(j);
+        if small>next
+            temp=XA(i);
+            XA(i)=XA(j);
+            XA(j)=temp;
+            temp1=inv_rA(i);
+            inv_rA(i)=inv_rA(j);
+            inv_rA(j)=temp1;
+        end
+    end
+end
+plot(XA,inv_rA)
+xlabel('XA');ylabel('-1/rA');
+//extending points to include XA=0.35
+XA(5)=0.35;inv_rA(5)=0.34;
+Area=inttrap(XA,inv_rA);
+W=Area*2000;
+printf("Amount of catalyst needed(kg) is %f",W)
+disp('The answer slightly differs from those given in book as trapezoidal rule is used for calculating area')
diff --git a/249/CH18/EX18.6/18_06.sce b/249/CH18/EX18.6/18_06.sce
new file mode 100755
index 000000000..ce199c8ac
--- /dev/null
+++ b/249/CH18/EX18.6/18_06.sce
@@ -0,0 +1,11 @@
+clear
+clc
+XA=0.35;
+FAo=2000;//mol/hr
+CAo=0.1//mol/litre
+eA=3;k=96;
+CA=CAo*((1-XA)/(1+eA*XA))
+rA=k*CA;
+//For mixed flow
+W=FAo*XA/rA;
+printf("The amount of catalyst needed(kg) is %f",W)
+\ No newline at end of file
diff --git a/249/CH18/EX18.7/18_07.sce b/249/CH18/EX18.7/18_07.sce
new file mode 100755
index 000000000..9a9922636
--- /dev/null
+++ b/249/CH18/EX18.7/18_07.sce
@@ -0,0 +1,3 @@
+clear
+clc
+printf("Its a theorotical qn")
+\ No newline at end of file