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
49
50
51
52
53
54
55
56
57
58
59
60
61
62
|
//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc.,USA,pp 436.
//Chapter-3 Ex3.2 Pg No. 96
//Title:Residence time and heat generation for four STR's in series
//========================================================================================
clear
clc
// COMMON INPUT
X_A=0.95;//Given conversion
t_batch=6;//Batch time to reach the desired conversion
N=4//No.of reactors in series
X_final=X_A;
//CALCULATION (Ex3.2.a)
k=log((1/(1-X_A)))/t_batch;//Refer equation 3.29 Pg No. 90
t_1=((1/(1-X_A))^(1/N)-1)/k;//Refer equation 3.40 Pg No. 94
t_Tot=N*t_1;
//OUTPUT (Ex3.2.a)
mprintf('\n OUTPUT Ex3.2.a');
mprintf('\n==================================================================');
mprintf('\nThe total residence time of the four reactors in series= %f hr',t_Tot);
//=======================================================================================
//Title:Heat generation in CSTR in Series
//=============================================================================================================
//CALCULATION (Ex3.2.b)
t_1=((1/(1-X_final))^(1/N)-1)/k;//Refer equation 3.40 Pg No. 94
for i=1:N
X(i)=1-(1/(1+k*t_1)^(i));
end
delQ_by_Q(1)=(X(1))/X_final; // Ratio of heat generated in 1st reactor
for i=1:N-1
delQ_by_Q(i+1)=(X(i+1)-X(i))/X_final; // Ratio of heat generated in 2nd, 3rd and 4th reactors
end
//OUTPUT (Ex3.2.b)
mprintf('\n========================================================================================\n')
mprintf('\n OUTPUT Ex3.2.b');
mprintf('\n==================================================================');
mprintf('\nReactor vessel \t Conversion \t Fraction of total heat released \n')
mprintf('\n==================================================================')
for i=1:N
mprintf('\n %d \t \t %0.3f \t \t \t %0.3f \n',i,X(i),delQ_by_Q(i))
end
//FILE OUTPUT
fid=mopen('.\Chapter3-Ex2-Output.txt','w');
mfprintf(fid,'\n OUTPUT Ex3.2.a');
mfprintf(fid,'\n==================================================================');
mfprintf(fid,'\nThe total residence time of the four reactors in series= %f hr',t_Tot);
mfprintf(fid,'\n==================================================================')
mfprintf(fid,'\nReactor vessel \t Conversion \t Fraction of total heat released \n')
mfprintf(fid,'\n==================================================================')
for i=1:N
mfprintf(fid,'\n %d \t \t %0.3f \t \t \t %0.3f \n',i,X(i),delQ_by_Q(i))
end
mclose(fid);
//=============================================================END OF PROGRAM================================
|
