summaryrefslogtreecommitdiff
path: root/650/CH8/EX8.9/9.sce
blob: 330887acdebc68d7e240e73aebff7358ab1b36e1 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
clc
d2=0.06; // m
d1=0.12; // m
k=0.44;
f=0.05;
L1=500; // m
g=9.81; // m/s^2

// v1=d2^2/d1^2*v2

// H_f=4*f*L1/16/d*v2^2/2/g
// H_c=k*v2^2/2/g
// H_f=4*f*L2/d*v2^2/2/g
// H_exit=v2^2/2/g

v2=sqrt(30*2*g/173.4);

Q=%pi*d2^2/4*v2;
disp("The rate of flow =")
disp(Q)
disp("m^3/s")