blob: b938d35021652384bc5b6c7ee55d5ee3727afdc7 (
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
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
|
clear
clc
CAo=[2;5;6;6;11;14;16;24];//mmol/m3
CA=[0.5;3;1;2;6;10;8;4];//mmol/m3
t=[30;1;50;8;4;20;20;4];//min
vo=0.1;//m3/min
for i=1:8
inv_rA(i)=t(i)/(CAo(i)-CA(i));
end
//Sorting CA and accordingly changing -1/rA for plotting graph between CA and -1/rA
for i=1:8
for j=i:8
if CA(i)>CA(j)
temp=CA(i);
CA(i)=CA(j);
CA(j)=temp;
temp1=inv_rA(i);
inv_rA(i)=inv_rA(j);
inv_rA(j)=temp1;
end
end
end
plot(CA,inv_rA)
xlabel('CA(mmol/m3)');ylabel('-1/rA(m3.min/m mol)');
disp('From the graph,we can see that we should use plug flow with recycle')
//From fig
CAin=6.6;//mmol/m3
R=(10-6.6)/(6.6-1);
//V=t*vo=area*vo
V=(10-1)*1.2*vo;
vr=vo*R;
printf("\n Part a")
printf("\n The vol of reactor is %f",V)
printf("m3 \n The recycle flow rate is %f",vr)
printf("m3/min")
//Part b,from fig
t=(10-1)*10;
t1=(10-2.6)*0.8;
t2=(2.6-1)*10;
//For 1 tank
V=t*vo;
//For 2 tank
V1=t1*vo;
V2=t2*vo;
Vt=V1+V2;
printf("\n Part b")
printf("\n For 1 tank volume is %f",V)
printf("m3 \n For 2 tank the volume is %f",Vt)
printf("m3")
printf("\n Part c")
disp('We should use mixed flow followed by plug flow')
//For MFR
tm=(10-4)*0.2;
Vm=tm*vo;
//For PFR
tp=5.8;//by graphical integration
Vp=tp*vo;
Vtotal=Vp+Vm;
printf("\n For MFR volume(m3) is %f",Vm)
printf("\n For PFR volume(m3) is %f",Vp)
printf("\n Total volume is %f",Vtotal)
printf("m3")
|
