"LESSON 8 EXAMPLE -  PART (b)"

"Constants and properties"
g = 9.81
rho = density (water, T=20, P = 101.325)
mu = viscosity(water,T=20, P = 101.325)

"Given info"
z_1 = 400
z_2 = 500
L = 1200
D = 0.05252
epsilon = 0.00026

"Volume flow rate"
Vol_dot = V_1*pi*D^2/4
Vol_dot = 0.005


"M.E.E."
P_1/(rho) + V_1^2/2 + g*z_1 = 0 + 0 + g*z_2 +e_loss
H_1 = p_1/(rho*g)

"Loss term"
e_loss = f*(L/D)*V_1^2/2 + 2*K_45_elbow*V_1^2/2 + 4*K_90_elbow*V_1^2/2 + K_valve*V_1^2/2 + K_exit*V_1^2/2

"Major loss"
Re = (rho*V_1*D)/mu
f = ((-2*log10((epsilon/(D*3.7065))-5.0452/Re*log10(((epsilon/D)^1.1098)/2.8257+5.8506/Re^0.8981))))^(-2)
f_T = ((-2*log10((epsilon/(D*3.7065))-5.0452/(1*10^(9))*log10(((epsilon/D)^1.1098)/2.8257+5.8506/1*10^(9)^0.8981))))^(-2) "Note: This is the fully turbulent friction for this pipe's relative roughness"


"Minor losses"
K_45_elbow = 16*f_T
K_90_elbow = 30*f_T
K_valve = 340*f_T
K_exit = 1