EQ_pumpstation_p(): power and total of a pumpstation

About:

Calculate power and total of a pumpstation.

©WETpython, 2015

\begin{equation} P = \frac{\gamma \times Q \times H_t}{\eta} \end{equation} \begin{equation} H_t = Z_{r2} - Z_{r1} + f_s \times ls + f_c \times lc + \sum \delta hs + \sum \delta hc \end{equation}

 Where,

$P$, pump power, W

$\gamma$, fluid specific weight, $\frac{N}{m^3}$;

$Q$, pumped flow, $\frac{m^3}{s}$

$h_t$, pump total head, $m$;

$Z_{r1}$, uptake reservoir water level, $m$

$Z_{r2}$, downstream reservoir water level, $m$

$f_s \times ls$, head loss in suction pipe, $m$

$f_c \times ls$, head loss in pressure pipe, $m$

$\sum \delta hs$, sum of minor losses in suction pipe,$m$

$\sum \delta hc$, sum of minor losses in compression pipe, $m$

Module: Equiment

Function: EQ_pumpstation_p(zr1,zr2,q,ypump,t,s,ds,ns,sumks,ls,dc,nc,sumkc,lc)

Definition of variables:

zr1,upstream reservoir water level, m;
zr2,donstream reservoir water level, m;
q , pumped flow, m3/s;
ypump, pump yeld;
t , water temperature, oC;
s, water salinity, g/kg;
ds, suction pipe diameter, m;
ns, suction pipe Manning coefficient, m;
sumks, sum of the minor losses K factor in suction pipe;
ls, suction pipe length, m;
dc, compression pipe diamter, m;
nc, compression pipe Manning coefficient, m;
sumkc, sum of the minor losses K factor in compression pipe;
lc, compression pipe length, m.

Sample code:

from Equipment import *

zr1 = 20.0 #upstream reservoir water level, m
zr2 = 60.0 #donstream reservoir water level, m
q = 0.1 #pumped flow, m3/s
ypump = 0.75 # pump yeld
t = 25 # temperature, oC
s = 2.5 #water salinity, g/kg
ds = 0.3 #suction pipe diameter, m
ns = 0.013 # suction pipe Manning coefficient, m
sumks = 5.0 # sum of the minor losses K factor in suction pipe
ls = 10.0 # suction pipe length, m
dc = 0.4 #compression pipe diamter, m
nc = 0.013 # compression pipe Manning coefficient, m
sumkc = 25.0 # sum of the minor losses K factor in compression pipe
lc = 2000.0 # compression pipe length, m

(power, ht) = EQ_pumpstation_p(zr1,zr2,q,ypump,t,s,ds,ns,sumks,ls,dc,nc,sumkc,lc)

print ("Pump power:", power,"W",", total head:", ht, "m")

Result:

Pump power: 60131.5187194 W , total head: 46.0357086847 m