blob: d57275982e4cc7238742992da41c9c2fd321e29a (
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
|
//Example 20.2
clc
L=100;//ft
D=1;//ft
D1=D/12;//inches
D2=D1*2.42;//centimetres
rho=62.4;//lb/ft^3
mu=1.5;//cp
Cv=4;
pv=100;//psi
G=1;
q=Cv*sqrt(pv/G);//maximum flow
disp('gpm',q,'q=')
printf("Let us start flow from q=30 gpm\n")
q=30;//gpm
q1=q/(60*7.48);//ft^3/sec
q2=q1*60*60;//ft^3/hr
Re=4*q2*rho/(%pi*mu*D2)//Reynolds number
//For this value of Reynolds number and for smooth pipe fanning friction factor is 0.005
f=0.005;//fanning friction factor
gc=32.2;
p=32*f*L*rho*q1^2/(144*%pi^2*gc*D1^5);//psi
P=pv-p
qmax=Cv*sqrt(P/G);
disp('gpm',qmax,'qmax=')
x=q/qmax//lift
|