blob: e834f47b4a3d3aadcc5d63c963289d250a52dd47 (
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
|
clc
u_A=1.35; // m/s
d_A=0.225; // m
d_B=0.150; // m
d_C=0.150; // m
d=5.6; //m
friction=2.5; // kW
power_req=12.7; // kW
rho=1000; // kg/m^3
rho_m=13560; // kg/m^3
g=9.81; // m/s^2
pC=35000; // Pa
pA=rho_m*g*(-d_B);
Area_A=%pi*d_A^2/4;
Area_B=%pi*d_B^2/4;
Area_C=%pi*d_C^2/4;
u_B=u_A*(Area_A/Area_B);
u_C=u_A*(Area_A/Area_C);
// Energy_added_by_pump/time = (Mass/time)*((pC-pA)/rho+(u_C^2-u_A^2)/2+g*(zC-zA))
Energy_added = Area_A*u_A*(pC-pA+rho/2*(u_C^2-u_A^2)+rho*g*d)/1000+friction;
Efficiency=Energy_added/power_req*100;
disp("Overall efficiency of the pump =")
disp(Efficiency)
disp("%")
|
