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
37
38
39
40
41
42
43
44
45
46
47
48
|
//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc.,USA,pp 436.
//Chapter-9 Ex9.1 Pg No.376
//Title:Model II- Volumetric Mass Transfer Coefficient (K)
//============================================================================================================
clear
clc
//INPUT
u0=[ 0.1 0.3 0.5 0.75 0.95 1.15];//Fluid Velocities (m/sec)
X=[0.923 0.872 0.846 0.775 0.728 0.664];//Conversion
h_by_h0=[1.26 1.44 1.66 2.0 2.3 2.7];//Height of bed under fluidized condition by height of packed bed
Epsilon_m=0.456;//Fraction of voids in packed bed
h0=0.75;//Height of packed bed (m)
k_r=4.45 ;//Reaction rate constant(sec-1)
W=5;//Weight of the bed (kg)
//CALCULATION
n=length(X);
for i=1:n
K0_L_by_u0(i)=log(1/(1-X(i)));//Refer equation 9.21 Pg No.371
L(i)=h_by_h0(i)*h0;
one_minus_epsilon(i)=(1-Epsilon_m)/h_by_h0(i);
k_rhob(i)=k_r*one_minus_epsilon(i);
K0(i)=K0_L_by_u0(i)*u0(i)/L(i);
K(i)=1/((K0(i).^(-1))-(1/k_rhob(i))); //Refer equation 9.19 Pg No.371
end
//OUTPUT
mprintf('\nThe values of K for given velocities')
mprintf('\n u (m/sec) \t K (sec-1) ');
mprintf('\n===================================================================');
for i=1:n
mprintf('\n %.3g \t \t %0.3f',u0(i),K(i));
end
//FILE OUTPUT
fid= mopen('.\Chapter9-Ex1-Output.txt','w');
mfprintf(fid,'\nThe values of K for given velocities')
mfprintf(fid,'\n u (m/sec) \t K (sec-1) ');
mfprintf(fid,'\n===================================================================');
for i=1:n
mfprintf(fid,'\n %.3g \t \t %0.3f',u0(i),K(i));
end
//==============================================END OF PROGRAM=============================================
|
