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

---
 1040/CH3/EX3.3/Ex3_3.sce | 66 ++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 66 insertions(+)
 create mode 100644 1040/CH3/EX3.3/Ex3_3.sce

(limited to '1040/CH3/EX3.3/Ex3_3.sce')

diff --git a/1040/CH3/EX3.3/Ex3_3.sce b/1040/CH3/EX3.3/Ex3_3.sce
new file mode 100644
index 000000000..5c6d5ae52
--- /dev/null
+++ b/1040/CH3/EX3.3/Ex3_3.sce
@@ -0,0 +1,66 @@
+//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc.,USA,pp 436.
+//Chapter-3 Ex3.3 Pg No. 97
+//Title:Effect of temperature on yield
+//================================================================================================================
+clear
+clc
+//INPUT
+C_A0=1;//Initial concentration of A
+C_B0=5;//Initial concentration of B
+E1=15;//Activation energy for first reaction(kcal)
+E2=20;//Activation energy for second reaction(kcal)
+X_A=0.88;// Total conversion of reactant A
+Y=0.81;//Yield for the reaction to produce C
+R=1.987;//Gas Constant(cal/K^-1 mol^-1)
+T_0=350;//Temperature (K)
+
+//CALCULATION
+//Assuming first order by taking concentration of B constant since B is in Excess
+C_A= C_A0*(1-X_A);//Unreacted amount of A
+C_B=C_B0-Y;//Unreacted amount of B
+k1_plus_k2_t=(X_A/(1-X_A));
+S=Y/X_A;//At 350K
+k1_by_k2=11.57;
+k1_plus_k2_by_k2=k1_by_k2+1;//Refer Ex3.3 for the coded equations
+k2_t=k1_plus_k2_t/k1_plus_k2_by_k2;
+k1_t=k1_plus_k2_t-k2_t;
+T=345;
+for i=1:7 
+T=T+5;
+Temp(i)=T;
+k1_dash_t(i)=k1_t*exp(((E1*1000/R)*((1/T_0)-(1/T))));//Arrhenius law
+k2_dash_t(i)=k2_t*exp(((E2*1000/R)*((1/T_0)-(1/T))));//Arrhenius law
+k1_plus_k2_t_new(i)=k1_dash_t(i)+k2_dash_t(i);
+X_A_new(i)=k1_plus_k2_t_new(i)/(1+k1_plus_k2_t_new(i));
+S_new(i)=((k1_dash_t(i)/k2_dash_t(i))/(1+(k1_dash_t(i)/k2_dash_t(i))));
+Y_new(i)=S_new(i)*X_A_new(i);
+end
+
+//OUTPUT
+mprintf('=======================================');
+mprintf('\n\t T \t  X_A \t  S \t  Y');
+mprintf('\n\t K \t (-) \t  (-) \t (-)');
+mprintf('\n======================================');
+for i=1:7
+    mprintf('\n\t %d \t %0.3f \t %0.3f \t %0.3f',Temp(i),X_A_new(i),S_new(i),Y_new(i));
+end
+  maximum=max(Y_new);
+  mprintf('\n\t\nThe maximum value of yield is %f ',maximum);
+  mprintf('\n\t\nHigh yield is obtained between 365K to 375K');
+
+//FILE OUTPUT
+fid=mopen('.\Chapter3-Ex3-Output.txt','w');
+mfprintf(fid,'=======================================');
+mfprintf(fid,'\n\t T \t  X_A \t  S \t  Y');
+mfprintf(fid,'\n\t K \t (-) \t  (-) \t (-)');
+mfprintf(fid,'\n======================================');
+for i=1:7
+    mfprintf(fid,'\n\t %d \t %0.3f \t %0.3f \t %0.3f',Temp(i),X_A_new(i),S_new(i),Y_new(i));
+end
+  maximum=max(Y_new);
+  mfprintf(fid,'\n\t\nThe maximum value of yield is %f ',maximum);
+  mfprintf(fid,'\n\t\nHigh yield is obtained between 365K to 375K');
+  mclose(fid);
+//======================================================END OF PROGRAM===================================================
+//Disclaimer:Refer Ex3.3 in the textbook The Arrhenius law equation has a typo error. Exponential term missing in the textbook 
+
-- 
cgit