summaryrefslogtreecommitdiff
path: root/1268/CH8/EX8.2/8_2.sce
blob: 607dcee41f27d977a05722834e415c6732e3ba39 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
clc;
disp("Example 8.2")
U=1.5 // in m/s
d=0.025 // in m
density= 1000 // in kg/m^3
mew=0.001 // in kg/ms
Re=d*density*U/mew
f=0.079/(Re^0.25)
l=25 // length of the pipe in m
delP=2*f*density*U*U*l/d
h1=delP/(density*9.81)
disp(h1,"Head loss is ")
h2=15
h=h1+h2
flow=%pi*d*d*U*density/4
power=flow*h*9.81
disp(power,"Theoretical Power required is ")