blob: a1cb511b2a8369f88cccd0ec9a79ade1d714c12d (
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
|
clc;
warning("off");
printf("\n\n example10.15 - pg466");
// given
Q=400; //[gpm]
p=1.1*62.4; //[lbm/ft^3]
mu=2*(6.72*10^-4); //[lb/ft*sec]
e=1.5*10^4;
// 4 inch schedule pipe
d=0.3355;
S=(%pi*(d^2))/4;
U4=Q/39.6;
ebyd=e/d;
w=3671/60;
pm=13.45*62.4;
g=32.1;
gc=32.174;
deltaz=2.5;
deltap=(g/gc)*(pm-p)*(deltaz);
betaa=((1)/(1+[(2*p*gc)*(deltap)]*(((0.61*S)/w)^2)))^(1/4);
d2=betaa*d;
Nre2=(4*w)/(%pi*d2*mu);
a=(1/30)*4.026;
b=(1/4)*(2.013-1.21);
c=(1/8)*(2.42);
if a<b then
if a<c then
opt=a;
else
opt=c;
end
else
if b<c then
opt=b;
else
opt=c;
end
end
printf("\n\n The pertinent orifice details are \n orifice diameter = %f in \n corner taps, square edge\n orifice plate not over %f in thick",d2*12,opt);
|
