blob: 7a9d32962fb68afd9fe5350992525c6e4063423a (
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
|
clc
//Initialization of variables
h=10 //m
g=9.81 //m/s^2
f1=0.019
f2=0.021
f3=0.020
z1=300 //m
z2=150 //m
z3=250 //m
d1=0.3 //m
d2=0.2 //m
d3=0.25 //m
//calculations
disp("part(a)")
Vbyg=h/(f1*z1/d1 +f2*z2/d2 *(d1/d2)^4 + f3*z3/d3 *(d1/d3)^4)
V1=sqrt(2*g*Vbyg)
Q=%pi/4 *d1^2 *V1
printf("\n Flow rate = %.3f m^3/s",Q)
disp('Part b')
Le2=z2*f2/f1 *(d1/d2)^5
Le3=z3*f3/f1 *(d1/d3)^5
Le1=z1
Le=Le1+Le2+Le3
V1byg=h*d1/Le/f1
V2=sqrt(2*g*V1byg)
Q1=%pi/4 *d1^2 *V2
printf("\n Flow rate = %.3f m^3/s",Q1)
|
