blob: 95d7ea47ed3cd2915e7b5c319777a4feea528459 (
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
|
clear
clc
//Initial Concentration(mol/litre)eactant in combined feed
CAo=10;CBo=10;
//For 90% Conversion
XA=0.9;
CAf=CAo*(1-XA);
//Instantaneous fractional yield of desired compound is
//Q(R/A)=dCR/(dCR+dCS)=CA/(CA+CB^1.5)
//But CA=CB so Q(R/A)=CA/(CA+CA^1.5)
//For Plug Flow
//Overall Fractional Yield(Qp)is
CA=CAf;
Qp=(-1/(CAo-CAf))*integrate('1/(1+CA^0.5)','CA',CAo,CAf);
CRf=9*Qp;
printf("\n Part a")
printf("\n For Plug Flow")
printf("\n Concentration of R(mol/litre) in the product stream is %f \n",CRf)
//Mixed Flow
//Overall Fractional Yield(Qm)is
Qm=CA/(CA+CA^1.5);
CRf=9*Qm;
printf("\n Part b")
printf("\n For Mixed Flow")
printf("\n Concentration of R(mol/litre) in the product stream is %f \n",CRf)
//Plug flow A,Mixed flow B
CAo=19;CB=1;
Q=-1/(CAo-CAf)*integrate('CA/(CA+CB^1.5)','CA',CAo,CAf);
CRf=9*Q;
printf("\n Part c")
printf("\n For Plug flow A,Mixed flow B")
printf("\n Concentration of R(mol/litre) in the product stream is %f \n",CRf)
disp('The result for plug flow varies as there seems to be typographical error in integration done in book')
|
