initial commit / add all books

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