summaryrefslogtreecommitdiff
path: root/1268/CH8/EX8.3/8_3.sce
blob: 48987eb25a6bde4b34a50c8f34797c1c91ee5a1b (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
clc;
disp("Example 8.3")
d=0.025 // diameter in m
l=120 // length in m
density= 1000
Q=2.525e-3  // volumetric flow rate in m^3/s
U=4*Q/(%pi*d*d)
Re=density*U*d/mew
f=0.0014+(0.125/(Re^0.32))
delP=2*f*l*U*U*density/d
disp(delP,"Pressure head is ")
Power=delP*Q
disp(Power,"Power required to overcome friction is ")