initial commit / add all books

7 files changed, 120 insertions, 0 deletions

diff --git a/249/CH9/EX9.1/9_01.sce b/249/CH9/EX9.1/9_01.sce
new file mode 100755
index 000000000..82eb1849a
--- /dev/null
+++ b/249/CH9/EX9.1/9_01.sce
@@ -0,0 +1,15 @@
+clear

+clc

+//Cp values(J/mol.k) given

+CpA=35;CpB=45;CpR=70;

+T1=25;T2=1025;

+Hr=-50000;

+//Enthalpy balance for 1mol A,1 mol B,2 mol R

+nA=1;nB=1;nR=2;

+dH=nA*CpA*(T1-T2)+nB*CpB*(T1-T2)+(Hr)+nR*CpR*(T2-T1);

+printf("\n dH(J) at temperature 1025C is %f \n",dH)

+if dH>0 then

+    printf("Reaction is Exothermic")

+else

+     printf("Reaction is Endothermic")

+end

diff --git a/249/CH9/EX9.2/9_02.sce b/249/CH9/EX9.2/9_02.sce
new file mode 100755
index 000000000..2b290f19e
--- /dev/null
+++ b/249/CH9/EX9.2/9_02.sce
@@ -0,0 +1,22 @@
+clear

+clc

+//Standard heat of reaction(J/mol) and Gibbs free energy(J/mol) 

+Ho=-75300;Go=-14130;

+R=8.3214;T1=298;

+//With all specific heais alike,dCp=0

+Hr=-Ho;

+K298=exp(-Go/(R*T1));

+//Taking different values of T

+T1=[2;15;25;35;45;55;65;75;85;95];//degree celcius

+T=[278;288;298;308;318;328;338;348;358;368];//kelvin

+for i=1:10

+K=K298*exp((Hr/R)*((1/T(i))-(1/298)));

+XAe(i)=K/(K+1);

+end

+plot(T1,XAe)

+xlabel('Temperature(C)')

+ylabel('XAe')

+disp(" From the graph we see temp must stay below 78 C if conversion of 75% or above is expected")

+

+

+

diff --git a/249/CH9/EX9.3/9_03.sce b/249/CH9/EX9.3/9_03.sce
new file mode 100755
index 000000000..e83b694ad
--- /dev/null
+++ b/249/CH9/EX9.3/9_03.sce
@@ -0,0 +1,24 @@
+clear

+clc

+//At 338 k

+XA=0.581;

+t=1;//min

+//From ex 9.2 at 65 degree celcius

+XAe=0.89;

+//For a batch reactor,k1t/XAe=-ln(1-XA/Xae)

+k1_338=-(XAe/t)*log(1-(XA/XAe));

+//At 25 degree celcius

+XAe1=0.993;

+T1=338;T2=298;

+R=8.314;

+//At 298 k

+XA1=0.6;

+t1=10;//min

+k1_298=-(XAe1/t1)*log(1-(XA1/XAe1));

+E1=(R*log(k1_338/k1_298))*(T1*T2)/(T1-T2)

+ko=k1_338/(exp(-E1/(R*T1)))

+//k1=ko*exp(-E1/RT)

+//k2=k1/k

+printf("\n The rate constants are k=exp[(75300/RT)-24.7] min-1")

+printf("\n k1=exp[17.2-(48900/RT)] min-1")

+printf("\n k2=exp[41.9-(123800/RT)] min-1 ")
\ No newline at end of file
diff --git a/249/CH9/EX9.4/9_04.sce b/249/CH9/EX9.4/9_04.sce
new file mode 100755
index 000000000..ea2565f34
--- /dev/null
+++ b/249/CH9/EX9.4/9_04.sce
@@ -0,0 +1,14 @@
+clear

+clc

+CAo=4;//mol/litre

+FAo=1000;//mol/min

+//Drawing locus of max rates on conversion-temp graph

+//tgen drawing optimum path for this system

+//integrating graphicaaly,we ger

+A=0.405;//litre/mol.min

+t=CAo*A;

+V=FAo*A;

+printf("\n Part a")

+printf("\n The space time needed is %f",t)

+printf(" min \n The Volume needed is %f",V)

+printf(" litres")
\ No newline at end of file
diff --git a/249/CH9/EX9.5/9_05.sce b/249/CH9/EX9.5/9_05.sce
new file mode 100755
index 000000000..3ee87e373
--- /dev/null
+++ b/249/CH9/EX9.5/9_05.sce
@@ -0,0 +1,28 @@
+clear

+clc

+//Concentration of A(mol/litre)

+CAo=4;

+//Flow rate of A(mol/min)

+FAo=1000;

+XA=0.8;

+Cp=250;//cal/molA.K

+Hr=18000;//cal/molA

+//Using Xa vs T chart of fig 9.3 at 80% conversion

+//Reaction Rate has the value 0.4 mol/min.litre

+rA=0.4;

+//From the performance eqn for mixed flow,Volume(l) is

+V=FAo*XA/rA;

+printf("\n Part a")

+printf("\n The size of reactor(litres) needed is %f",V)

+slope=Cp/Hr;

+//Using graph

+Qab1=Cp*20;//cal/molA

+Qab=Qab1*1000;//cal/min

+Qab=Qab*0.000070;//KW

+printf("\n Part b")

+printf("\n Heat Duty(KW) of precooler is %f",Qab)

+Qce1=Cp*37;//cal/molA fed

+Qce=Qce1*1000;//cal/min

+Qce=Qce*0.000070;//KW

+printf("\n Heat Duty(KW) of postcooler is %f",Qce)

+

diff --git a/249/CH9/EX9.6/9_06.sce b/249/CH9/EX9.6/9_06.sce
new file mode 100755
index 000000000..654c349e6
--- /dev/null
+++ b/249/CH9/EX9.6/9_06.sce
@@ -0,0 +1,9 @@
+clear

+clc

+FAo=1000;//mol/min

+//Drawing trial operating lines with a slope of 1/72 and for each evaluating integral dXA/-rA

+//From graph 

+Area=1.72;

+V=FAo*Area;

+printf("\n The volume of adiabatic plug flow reactor is %f",V)

+printf("litres")
\ No newline at end of file
diff --git a/249/CH9/EX9.7/9_07.sce b/249/CH9/EX9.7/9_07.sce
new file mode 100755
index 000000000..221354808
--- /dev/null
+++ b/249/CH9/EX9.7/9_07.sce
@@ -0,0 +1,8 @@
+clear

+clc

+//Using ex 9.6 and finding optimum recycle area

+FAo=1000;//mol/min

+Area=(0.8-0)*1.5;

+V=FAo*Area;

+printf("\n The volume required is %f",V)

+printf("litre")
\ No newline at end of file