blob: 5682da97046014eb2c619230e7a4f25560fd3e07 (
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
28
29
30
31
32
33
34
|
//Example 5.7
clc;
clear;
close;
format('v',8);
//Given data :
Z1=0;//meter
Z2=5;//meter
Q=300*10^-3;//m/s
D1=0.3;//meter
D2=0.6;//meter
A1=%pi*D1^2/4;//m^2
A2=%pi*D2^2/4;//m^2
v1=Q/A1;//m/s
v2=Q/A2;//m/s
p1=100;//kN/m^2
p2=600;//kN/m^2
g=9.81;//constant
Vhead11=v1^2/2/g;//meter
Vhead22=v2^2/2/g;//meter
Phead11=p1/g;//meter
Phead22=p2/g;//meter
E1_11=Z1+Vhead11+Phead11;//meter
E2_22=Z2+Vhead22+Phead22;//meter
if E1_11>E2_22 then
disp("Energy at section 1-1("+string(E1_11)+" meter) is greater than energy at section 2-2("+string(E2_22)+" meter). Flow of water is from section 1-1 to 2-2.");
HeadLoss=E1_11-E2_22;//meter
disp(HeadLoss,"Head Loss in meter : ");
else
disp("Energy at section 2-2("+string(E2_22)+" meter) is greater than energy at section 1-1("+string(E1_11)+" meter). Flow of water is from section 2-2 to 1-1.");
HeadLoss=E2_22-E1_11;//meter
disp(HeadLoss,"Head Loss in meter : ");
end
|
