blob: 9de3ab2dbc4a91684fd2d93bcf22abcfce997749 (
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
|
//exapple 4.2
clc; funcprot(0);
// Initialization of Variable
rd=[0 1 2.5 5 10 15 17.5]/100;//radial distance from pipe
dlv=[0 0.2 0.36 0.54 0.81 0.98 1]/100;//differnce in liquid levels
r=[.175 .165 .150 .125 .075 .025 0];//
g=9.81;
R=8.314;
rho=999;
temp=289;
P1=148*1000;
M=7.09/100;
pi=3.12
rhoCl2=P1*M/R/temp;//density of Cl2
nuCl2=1/rhoCl2;//specific volume of Cl2
function[y]=P2(x);
y=P1+x*(rho-rhoCl2)*g;
endfunction
for i=1:7
y=P2(dlv(i));
u(i)=sqrt(2*P1*nuCl2*log(y/P1));
a(i)=u(i)*r(i);
end
clf();
plot(r,a);
xtitle("","r (m)","u*r (m^2/s)");
s=0;
for i=1:6//itegration of the plotted graph
s=abs((r(i)-r(i+1))*.5*(a(i)+a(1+1)))+s;
end
s=s-0.01;
Q=2*pi*s;
disp(Q,"volumetric flow rate (m^3/s):");
disp(Q*rhoCl2,"mass flow rate of chlorine gas (kg/s)")
|
