fPump: Performance of a centrifugal pump

Description Usage Arguments Value Author(s) References See Also Examples

View source: R/fPump.R

Description

Compute the performance of the different parts of a centrifugal pump fed by a frequency converter following the affinity laws.

Usage

1
fPump(pump, H)

Arguments

pump

list containing the parameters of the pump to be simulated. It may be a row of pumpCoef.

H

Total manometric head (m).

Value

lim

Range of values of electrical power input

fQ

Function constructed with splinefun relating flow and electrical power

fPb

Function constructed with splinefun relating pump shaft power and electrical power of the motor

fPh

Function constructed with splinefun relating hydraulical power and electrical power of the motor

fFreq

Function constructed with splinefun relating frequency and electrical power of the motor

Author(s)

Oscar Perpiñán Lamigueiro.

References

See Also

NmgPVPS, prodPVPS, pumpCoef, splinefun.

Examples

 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
library(latticeExtra)

data(pumpCoef)
CoefSP8A44 <- subset(pumpCoef, Qn == 8 & stages == 44)

fSP8A44 <- fPump(pump = CoefSP8A44,H = 40)
SP8A44 = with(fSP8A44,{
                Pac = seq(lim[1],lim[2],by = 100)
                Pb = fPb(Pac)
                etam = Pb/Pac
                Ph = fPh(Pac)
                etab = Ph/Pb
                f = fFreq(Pac)
                Q = fQ(Pac)
                result = data.frame(Q,Pac,Pb,Ph,etam,etab,f)})

#Efficiency of the motor, pump and the motor-pump
SP8A44$etamb = with(SP8A44,etab*etam)
lab = c(expression(eta[motor]), expression(eta[pump]), expression(eta[mp]))
p <- xyplot(etam + etab + etamb ~ Pac,data = SP8A44,type = 'l', ylab = 'Efficiency')
p+glayer(panel.text(x[1], y[1], lab[group.number], pos = 3))

#Mechanical, hydraulic and electrical power
lab = c(expression(P[pump]), expression(P[hyd]))
p <- xyplot(Pb + Ph ~ Pac,data = SP8A44,type = 'l', ylab = 'Power (W)', xlab = 'AC Power (W)')
p+glayer(panel.text(x[length(x)], y[length(x)], lab[group.number], pos = 3))

#Flow and electrical power
xyplot(Q ~ Pac,data = SP8A44,type = 'l')

solaR documentation built on Oct. 19, 2021, 9:06 a.m.