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-8 Ex8.3.a Pg No.348
//Title:Apparent rate constant and consistency
//===========================================================================================================
clear
clc
//INPUT
LHSV_inv=[0.75 1.39];//Refer table 8.2 Test Results (Liquid Hourly Space Velocity)
X_S=[0.77 0.83];//Refer table 8.2 Percentage Sulphur removal
//CALCULATION
for i=1:2
kapp_rhob(i)=log((1/(1-X_S(i))))*(1/LHSV_inv(i));//Refer Equation 8.21
end
L=LHSV_inv(2)/LHSV_inv(1);
kapp_ratio=kapp_rhob(1)/kapp_rhob(2);
n=log10(kapp_ratio)/log10(L);
//OUTPUT
//Console Output
mprintf('\n\tThe Apparent rate constants are \n\t Run1: %0.2f hr-1 \n\t Run2: %0.2f hr-1 ',kapp_rhob(1),kapp_rhob(2))
mprintf('\n\tThe exponent value = %0.1f ',n);
mprintf('\n\tExponent value greater than 0.5 hence the Eq.(8.24) is not consistent');
mprintf('\n\tThe error may be due to error in assuming a first order reaction');
//File Output
fid = mopen('E:\Chapter8_Ex3_a_Output.txt', 'w');
mfprintf(fid,'\n\tThe Apparent rate constants are \n\t Run1: %0.2f hr-1 \n\t Run2: %0.2f hr-1 ',kapp_rhob(1),kapp_rhob(2))
mfprintf(fid,'\n\tThe exponent value %0.1f ',n);
mfprintf(fid,'\n\tExponent value greater than 0.5 hence the Eq.(8.24) is not consistent');
mfprintf(fid,'\n\tThe error maybe due to error in assuming a first order reaction');
mclose(fid);
//============================================END OF PROGRAM=================================================
|