summaryrefslogtreecommitdiff
path: root/1040/CH8/EX8.3.a/Chapter8_Ex3_a.sce
blob: d3e500ede4bbb87513c49700fe9d06f17127952d (plain)
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=================================================