1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
|
//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc.,USA,pp 436.
//Chapter-3 Ex3.6 Pg No. 114
//Title: Optimum reaction temperature
//===========================================================================================================
clear
clc
//INPUT
del_H=-20*10^3;//Heat of reaction(cal)
T_eq=[500 700];//Equivalent temperatures (K)
R=1.987;//Gas Constant (cal/mol K)
E2_by_E1=2;//Ratio of activation energy
//CALCULATION
T_opt(1)=T_eq(1)/(1+(log(E2_by_E1)*(R/(-del_H)))*T_eq(1));//Refer equation 3.63 Pg No. 113
T_opt(2)=T_eq(2)/(1+(log(E2_by_E1)*(R/(-del_H)))*T_eq(2));
delta_T(1)=T_eq(1)-T_opt(1);
delta_T(2)=T_eq(2)-T_opt(2);
//OUTPUT
mprintf('\n \t \t Temperature_1\t Temperature_2 ');
mprintf('\n \t \t==================================');
mprintf('\n(T_eq - T_opt)(K): \t%0.0f \t\t%0.0f',delta_T(1),delta_T(2));
mprintf('\n T_opt(K):\t \t%0.0f\t\t%0.0f', T_opt(1),T_opt(2));
fid= mopen('.\Chapter3-Ex6-Output.txt','w');
mfprintf(fid,'\n \t \t Temperature_1\t Temperature_2 ');
mfprintf(fid,'\n \t \t==================================');
mfprintf(fid,'\n(T_eq - T_opt)(K): \t%0.0f \t\t%0.0f',delta_T(1),delta_T(2));
mfprintf(fid,'\n T_opt(K):\t \t%0.0f\t\t%0.0f', T_opt(1),T_opt(2));
mclose(fid);
//=========================================================END OF PROGRAM=====================================
//Disclaimer:There is an arithmetic error in the optimum temperatures obtained in the textbook.
// Based on the values (T_eq - T_opt)1=17 and (T_eq - T_opt)2=32 the optimum temperatures obtained are
// T_opt1=483 K and T_opt2=668 K respectively.
|
