diff options
Diffstat (limited to '1040/CH3/EX3.3')
| -rw-r--r-- | 1040/CH3/EX3.3/Chapter3_Ex3_Output.txt | 15 | ||||
| -rw-r--r-- | 1040/CH3/EX3.3/Ex3_3.sce | 66 |
2 files changed, 81 insertions, 0 deletions
diff --git a/1040/CH3/EX3.3/Chapter3_Ex3_Output.txt b/1040/CH3/EX3.3/Chapter3_Ex3_Output.txt new file mode 100644 index 000000000..66f0546e5 --- /dev/null +++ b/1040/CH3/EX3.3/Chapter3_Ex3_Output.txt @@ -0,0 +1,15 @@ +=======================================
+ T X_A S Y
+ K (-) (-) (-)
+======================================
+ 350 0.880 0.920 0.810
+ 355 0.909 0.913 0.830
+ 360 0.931 0.905 0.843
+ 365 0.948 0.896 0.849
+ 370 0.961 0.887 0.852
+ 375 0.970 0.878 0.851
+ 380 0.977 0.868 0.848
+
+The maximum value of yield is 0.852042
+
+High yield is obtained between 365K to 375K
\ No newline at end of file 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
+
|