summaryrefslogtreecommitdiff
path: root/1962/CH12/EX12.2/example12_2.sce
blob: 6ed09195e2027eb02ec0f75d204d877b51fdfc0c (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

//example 12.2
//page 443
clc; funcprot(0);
//initialisation of variable
Q=0.4;
pi=3.14;
u2=31.4;//velocity
Gamma=9.81;//unit weight
g=9.81;
omega=2*pi*1500/60;//radial velocity
r2=0.2//m
b2=0.03//m
//part1
V2r=Q/2/pi/r2/b2;//radial velocity
V2t=u2-V2r*0.577;
v2=V2r/0.866;//speed
V2=sqroot(V2r^2+V2t^2);
disp(V2,"relative velocity (m/s):");
//part2
H=u2*V2t/g;
disp(H,"head of water (m)");
//part3
P=Gamma*Q*H;
disp(P,"power required(kW)");
clear