R Graphical Manual

Browse All

Last data update: 2014.03.03

R: Performance of a centrifugal pump

C_fPump

R Documentation

Performance of a centrifugal pump

Description

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

Usage

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

Abella, M. A., Lorenzo, E. y Chenlo, F.: PV water pumping systems based on standard frequency converters. Progress in Photovoltaics: Research and Applications, 11(3):179â191, 2003, ISSN 1099-159X.
PerpiÃ±Ã¡n, O, EnergÃa Solar Fotovoltaica, 2015. (http://oscarperpinan.github.io/esf/)
PerpiÃ±Ã¡n, O. (2012), "solaR: Solar Radiation and Photovoltaic Systems with R", Journal of Statistical Software, 50(9), 1-32, http://www.jstatsoft.org/v50/i09/

Examples

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')