×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Log In

Come Join Us!

Are you an
Engineering professional?
Join Eng-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!
  • Students Click Here

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Students Click Here

Jobs

Curve fitting pump curves
10

Curve fitting pump curves

Curve fitting pump curves

(OP)
Dear all

Can anyone tell me which curve fit is the best fit to use when working with head versus flow rate data and head versus efficiency data on centrifugal slurry pumps with a low specific speed (100-700 or so)?

I am a Mechanical Engineer fresh out of university and am currently working for a slurry pump manufacturing company.

One of the projects I have been assigned is to help fix up pump selection software that the company currently employees in producing sales quotes. In its current form, the software tends to misquote on certain efficiency and %QBEP values when working with certain pump curves.

I believe that part of the problem lies with the way the software plots a curve through the pump input data it is given. Currently the software will accept 7 data points per speed setting (since these are slurry pumps, the impellers are not trimmed but rather belts are used to adjust the pump speeds). The data points inputted per speed setting consist of the following:

7 flow rate values
7 corresponding head values
7 corresponding efficiency values

Once these data points have been entered into the program it plots a 4th order polynomial through the data and generates curves of head versus flow rate and efficiency versus flow rate.

I have contested that plotting a 4th order polynomial through the data does not make the generated curves more accurate. I cannot remember from theory where I saw 4th order terms in any pumping related equations. I argue that a 2nd order curve is far more likely to represent the statistically most likely value of either head or efficiency that a customer may obtain for a given flow rate should they test the pump, even though the actual data points inputted into the program deviate from the curve more than they would for a 4th order curve.  

Thus I recommended a 2nd order curve fit for both head versus flow rate data and efficiency versus flow rate data and increasing the number of data points per curve from 7 to 10 so as to get more data values throughout the range of flow rates for a particular curve.

I have been challenged in this regard by the software developer, who is a senior mechanical engineer, and who tells me a 4th order curve will describe the data perfectly. He holds the position that a 4th order curve will fit irregularities in the actual performance of the pump that my university theory cannot account for and has dismissed my arguments on the grounds of inexperience in the field.

I have searched far and wide in literature for some guidance, including curve fitting techniques, correlation coefficients, residual plots etc but have not found anything worthwhile yet. I would appreciate the advice of anyone who has had experience with this type of problem.

Kind Regards,

Patrick               

RE: Curve fitting pump curves

The answer to which correlation fit is more accurate is that you could both be right, the physical curve will follow the laws of physics "Newtonian" so would tend to be a 2nd order polynomial since it would only have terms with flow squared in the physical world, however mathematically the more points you use for the fit the more accurately you fit a curve, hence the fourth order polynomial would more accurately fit the data - beware if you extrapolate past the first or last point, fourth porder will often curle away and give you less accuracy.

However I believe you are looking in the wrong area. A fact of life is that the efficiency of a working pump measured in the field will almost never match the one calculated at the bid phase for any of the following reasons and more:

1) Manufacturing tolerances
2) Installation and piping losses
3) Instrumentation calibration
4) Fluid physical properties (especially for a slurry)
5) location ambient conditions
6) Impeller or driver fouling or loss of efficiency

The goal of the bidding software would be to be as close as possible but state the tolerances to real world, since you dont know any of these factors before a pump is installed in the field generally some level of conservatism is needed based on experience.

RE: Curve fitting pump curves

I have done this with a little different approach.  We all know that 2nd order should work, but there are anomolies that keep this from being close enough.
First, select a pump.
When I did this it was with small multistage centrifugals.  We built 5 pumps, and tested 5 points on each one 10 times.
This will develop enough data to decribe each data point with a statistical distribution.
Now run your 2nd order fit and see what the mismatch is.  I usually ended up using a 3rd order fit.
Don't do anything that the statistics will not support.
4th order fits the data nicely, but it doesn't describe the real system.

= = = = = = = = = = = = = = = = = = = =
Rust never sleeps
Neither should your protection
http://www.trent-tube.com/contact/Tech_Assist.cfm

RE: Curve fitting pump curves

Whilst not any kind of expert in software modelling, I do however test real world pumps dynamically in my facility. This is done to compare refurbished pump performance  against manufacurers curves. It is absolutely correct that the more measurement points taken, the more accurate the result, particularly where small changes in head produce large changes in flow (or vice-versa)I do not use any less than 10 measurement points, and sometimes up to 15. I have found with the modest(free!)curve fitting software that I use, that quadratic fits have compared best with manufacturers published curves. Surely any software used to determine pump performance at various speeds/trims can only follow the affinity laws, which while mathematically hold true, are not always entirely correct in the real world as monaco8774 says, particularly as impeller vane angles stay the same whatever speed/trim is chosen. After all is said, a curve is only averaging out measurement points taken.
Regards
John   

RE: Curve fitting pump curves

Quote (monaco8774):

However I believe you are looking in the wrong area. A fact of life is that the efficiency of a working pump measured in the field will almost never match the one calculated at the bid phase for any of the following reasons and more:

I agree.

With respect to curve fitting, I tend to use a 4th order curve fit - as you say, it does a better job fitting the data that you have. As others have also said, the curve fitting is only as good as the data set. Maybe more data points would improve accuracy?

"Do not worry about your problems with mathematics, I assure you mine are far greater."   
Albert Einstein
Have you read FAQ731-376 to make the best use of Eng-Tips Forums?

RE: Curve fitting pump curves

I would suggest putting some of your data into Excel (or similar program) and letting Excel curve fit them.  Excel will curve fit for a variety of polynomal equations, give you the coefficients and the correlation coefficient.  Dead easy to change the polynominal equation and see how the curve changes to match the data.

Personally, I've had poor experience with 4th order polynomials on pump curves in that I see odd inflections at points not supported by the data.  Typically, I use a 2nd order polynomial or sometimes a 3rd order and get more than sufficient accuracy for my needs.  Now, I'm usually plotting data for 'common' centrifugal pumps, maybe that's the reason.

One caution with Excel, you can easily get round-off errors with the polynominal coefficients shown on the graph if you aren't careful.  I always use engineering notation for the coefficients, typically go with 3 significant factors and then calculate the predicted value(s) against the input value(s) to make sure I don't have any odd discrepancies.

RE: Curve fitting pump curves

(OP)


Thank you all for your suggestions. I am planning to test one of our newer pumps this week and get as much data as possible (running the pump up and down the selected test curves) and then use that to tell me what type of curve fit seems appropriate. Hopefully with sufficient test data I should be able to tell how the pump behaves at low flow rates and thus get a better idea of the curves to be used in general using free curve fitting software or a program such as MATLAB or Excel. I will also keep in mind the realities of the situation and try keep my analysis conservative. One of the irony’s of the slurry pumping business is that all out pump curves are generated in clear water, with adjustments made later on to accommodate for the slurry based on experience and a few correlation factors worked out empirically.

Anyway thank you again for your advice

Regards,
Patrick   

RE: Curve fitting pump curves

Another trick I've used is to just curve fit part of the curve if the operating range you are interested in is a relatively narrow part of the entire pump curve.  I used this on a cooling water pump (vertical) recently that I was not getting a good curve fit for trying to match the whole curve.  But, over the range I was interested it, a linear fit did the job.

Another trick is to model the curve in sections and then use Excel's "IF" statement to pick the right equation.  Lots of time this is overkill but it is an option.

RE: Curve fitting pump curves

A great way to do curve fitting using excel:

List your points (xi,yi).
Develop a functional form for your estimate of y, call it yhati = f(xi) with unknown parameters alpha, beta, gamma, delta etc.
For example yi = alpha + beta*xi + gamma * xi^2 + delta + xi^3 etc

Guess values for alpha and put them in a cell.

Type in your equation for yhati in terms of xi and the guessed paramters (guessed parameters as an absolute reference).  Copy it all the way down the column.

Compute residual Ri = yhati - yi
Compute square of residuals Ri^2
Sum the square of the residuals over all i.
Your objective is to find the values of unknown constants alpha, beta, gamma, delta etc which will minimize the sum of squares of the residuals.  Use the solver function to do this. The solver dialogue box has the general form:
Set target cell ____ to max/min/value by changing cells ____.
The target cell is your sum of squares of residuals...
You want to minimize (vs max or value) it...
By changing cells: alpha, beta, gamma, delta etc.

The advantage is you can pick whatever functional form you want. You don't have to choose between x^2 and x^4... you can have both as well as a cos function or anything else. It will pick coefficients to minimize the sum of squares of the residuals.

I think you will find a small investment in time in learning the solver tool (fairly simple) will pay off in many different types of problems, not just curve fitting.

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Curve fitting pump curves

(OP)
Thanks for giving me the link to www.systat.com. I am currently reviewing the software at the moment. In many ways it encompasses an automated version of the routine written out to me by electricpete. I agree with him in that the solver tool in excel is fairly powerful when utilised correctly. In general (outside the scope of my project) in terms of coefficients and powers for the equations used to describe data I quite like the fact that it has the ability to include non integer power values in the curve.   

This issue with regards to curve fitting has been quite educational for me in terms of trying to figure out what is accurate in terms of given data, what is realistic in terms of the system and what is acceptable by my company in terms of standards and time spent on the project.

Thank you all again for you help,

Regards,
Patrick


  

RE: Curve fitting pump curves

3
At the risk of saying some unpopular things, I would venture the following comments (please take these in a friendly spirit if you disagree, but do state your objections):

Fitting polynomials to data is one of the oldest procedures in the book.  You first need to decide just how accurate your data points are.  If the measurements are to be regarded as "perfect", i.e., you blithely ignore the fact that all measurements are subject to statistical error, the next question is:

How best do you find a curve of some kind that has the right mathematical properties, with respect to (a) passing through all the points, (b) having the proper first and second derivatives, and (c) a reasonable ability to extrapolate, however mildly, outside the range of the data.

In all these respects, polynomials (especially of a high order) are generally a very poor choice.

In my opinion, the proper answer is to use cubic splines to fit the data.  These have been around for decades, are extremely well researched, and are not difficult to understand or apply, although unfortunately many engineers have not had proper training in the use of this technique.  There are several choices in the available software to address item (c) above.  A cubic spline with the right "end conditions" specified will always do an excellent job and, in general, will out-perform a polynomial, no matter how high its order.

Another issue arises with respect to the handling of interpolation for compressor speeds different from the ones for which data have been provided.  Here, one is much better off using the "fan laws" to fit the original data in transformed coordinates (e.g., Hp/N^2 versus Qs/N for polytropic head, and PolEff versus  Qs/N for polytropic efficiency) to ensure that a rational interpolation is performed for any N.  Otherwise, you'd be condemned to use another arbitrary order polynomial to interpolate between the curves.  Of course, the fan laws can break down if the impeller tip velocity is too close to the sonic velocity, but this is not an issue generally with compressors used in the process industries (hydrogen compressors may be the main exception).

One of the nicest things about splines is that the vectors of spline interpolation coefficients can be saved, after the first time through, for future use.  This reduces the computational effort in subsequent interpolations.  The coefficients are revised only if the original data set changes.

A great reference for cubic splines (and a huge variety of other numerical problems that are relevant to engineers) is:

W. H. Press, et al: "Numerical Recipes in C", 2nd Edition (Cambridge Univ. Press, 1992).  There is a lot of published software in this book, beautifully explained, and also available at a nominal cost from the publisher.  Fortran 77, Fortran 90, and Pascal versions of the book and software are also available.

That said, it is sheer folly for any manufacturer to insist that his measurements are absolutely perfect.  Nothing in the universe is known perfectly, no matter how "accurate" your measurement techniques.  In a real-world system that is subject to dynamic disturbances, it is especially galling that anyone would adopt such a dogmatic position.  If a simple plot of the data shows abnormal bumps in the curvature, you can be pretty sure that the measurements are suspect.  Here, I would resort to use of the fan laws at least to eliminate the discrepancy.

RE: Curve fitting pump curves

Those are good points to broaden the discussion.

Each tool has it's place, and which tool is the best certainly depends on the data, the amount of noise/random error, and your knowledge of any underlying functional form.

Here is the book mentioned above available for free:

http://www.library.cornell.edu/nr/cbookfpdf.html

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Curve fitting pump curves

(OP)
I definitely agree with you UmeshMathur. There is a lot of valuable information for me contained within your post. I am trying to persuade my company to take into greater account error when testing pumps, developing pump databases and publishing pump curves (something that they currently do not do). I am hoping to develop a new set of standards with regards to the handling and processing of test data so that we can improve the utility of our performance curves and standardize them (in many ways they are not at the moment). I have a meeting with the software developer next week and I plan to look through the relevant sections of “Numerical recipes” before discussing the technicalities of the program and the way it handles data so that him and I may agree on methods to process test data in the program used for pump selections. I cannot agree more with the fact that no data points are ever known perfectly and this has been the basis of my argument to re-examine the way in which curves are applied through the data. Since I am relatively new to the field, it has been really interesting learning in depth, the methods and problems associated with curve fitting. For someone like me, who is inherently interested in the correlation between experimental and numerical values, this is a great project to develop my skills.     

I have noticed though that in general throughout out local slurry pump industry (South Africa and Africa) these issues have taken a back seat and have not been fully investigated. Even if we specify efficiency values with a relatively high confidence, in general our customers do not fully check these efficiency values (when the pump is fitted into the actual pipe network that it is to operate in with the slurry) or require us to pay penalties for discrepancies. There is also little information with regards to the change of efficiency values over time. I would have thought more information would have been supplied locally with regards to these factors and that independent certification tests would have been done on our pumps by our customers to second check our results. Currently we do have customer certification tests but they are conducted in our test facility (a clear water test facility) using our equipment and our data processing techniques. It may be that for the moment there is simply not enough of a requirement from the mining and process industries we sell our products to, to second guess our data.

RE: Curve fitting pump curves

As far as the cubic spline, what I understood from above is that it would be a good solution if we have high confidence the data is exact but not as good in the presence of noise.  For example imagine one bad data point outlieing from the smooth curve.  The spline will force the curve directly through that point.  The polynomial fit would generally draw a smooth curve that passes close to most of the points.


One thing to note about the excel/solver method I mentioned above. Sometimes it works great on the first try. Sometimes excel has problems selecting which terms to use first when you put a lot of terms in.  So you can do it in stages.  First fit yhat1 to the forme you think works best, maybe a*x^2 or a*x^n.  Then accept the solution for a and n and plot first set of residuals (y-yhat1).  Now look at that residual curve as your new objective function and guess the funcitonal form (maybe b*exp(-c*x).  User solver again to determine b and c to best match the first set of residuals (treating previoulsy-determined a and b as a constant).  Now you have a second set of residuals to keep repeating the process. Each step you are working with a smaller and smaller set of residuals so the longer you spend, the closer you will get to matching your data points.

One last thing to consider, if you have not just one set of test data, but multiple sets of test data, you can make far better guesses to eliminate random/non-repeatable variability (of course you can't get rid of bias errors that show up the same each time).

If you have any data sets to post, I would be glad to do a solution with excel and post a link to the file back here.

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Curve fitting pump curves

Patrick,

You may want to pay closer attention to how accurately the actual pump shaft speed is known during each test state.  I have found that seemingly minor variations in the shaft speed can have significant implications if they are not properly included in the mathematical adjustment of the data.  If the pump is driven by an induction motor, the variation in the actual shaft speed due to slip at different torque loadings can be troublesome.  If the pump is driven by an induction motor through v-belts or flat belts, then belt slip can introduce further speed variation from the nominal speed.

As others have stated above, I have found Excel to be very handy in resolving pump and fan problems.  Surprisingly often, I have found that adjusting for the slip of the driving induction motor was highly beneficial in the evaluation of apparent performance problems.

RE: Curve fitting pump curves

Patrick02:

My comments earlier mentioned compressors, but apply equally to centrifugal pumps.

SLURRY PUMPS
A major source of contention with slurries, of course, is the characterization of slurry properties.  With suspended solids, one must worry about accurate knowledge of the solids bulk density, percentage solids in the pumped fluid, adequate solids dispersion into the bulk, non-Newtonian behavior, and variability in solids content over time.  Such measurements are very hard to replicate in a test lab and all guarantees must reflect the lack of certainty in slurry characterization.

I would recommend the following as an excellent book for discussion of fluid flow problems for slurries:  R. P. King, "Introduction to Practical Fluid Flow" Chapters 3, 4, 5, and 6, (Butterworth-Heinemann, 2002).  King shows the very many intricacies in this kind of system that can waylay the casual user.

Back-calculation of pump or compressor efficiency in the field requires measurement of flows, temperatures, and pressures (all subject to previously mentioned measurement errors) and also user-specified values for fluid density and other transport properties.  I would say that field measurement of slurry flows is itself  a huge area of uncertainty.  For non-Newtonian fluids, several other parameter inputs are required that are themselves back-calculated from field data that are difficult to replicate in the laboratory.  The uncertainty in these measured and unmeasured inputs used for such calculations is often the source for a lot of wrangling between vendors and owners.  Consequently, I regard those slurry pump vendors who provide guarantees with monetary penalties attached as intrepid warriors.

CURVE FITTING
With respect to electricpete’s latest post, I agree that splines may not be appropriate if there is significant evidence of scatter in the underlying data.  However, most compressor vendor curves are smoothened before being provided to clients, so such discrepancies are not in evidence.

With pump curves that show much scattering in the data, the procedure I would recommend is first to reduce the data using the fan laws.  Then, select a simple, non-polynomial functional form that has the proper behavior with respect to curvature and boundary conditions.   Unfortunately, a simple parabola rarely suffices.  The final fitting of the function’s coefficients against the reduced data is performed a lot more easily than is possible in Excel using standardized statistical software, a lot of it available on the web free or for a nominal cost.  See, for example, http://www.nlreg.com/ for a very nice package that has a great non-linear solver, nice graphics, and an easy input language.

Whichever functional form you choose, it is essential to plot the predicted and actual values (Hp/N^2 or polytropic efficiency versus Qs/N).  This would show evidence of a systematic error in predictions owing to an erroneous choice of function.  Patrick02 may be pleased to hear that 4th order polynomials are assuredly poor performers in fitting such data in a huge proportion of cases.  This is because they exhibit up to two points of inflection which distorts the curve fit over much of the range of the data.

CHALLENGE PROBLEM
To make things interesting, and just to show you how hard the task of finding a proper fitting function can be, I extracted the following set of data for 17 points taken from a real compressor vendor’s curves:

Qs/N: 0.8,0.833333,0.866667,0.9,0.933333,0.966667,1,1.033333,1.066667,1.1,1.133333,1.166667,1.2,1.233333,1.266667,1.3,1.333333
1000*Hp/N^2:
0.297,0.296,0.294,0.293,0.292,0.29,0.288,0.286,0.283,0.28,0.276,0.272,0.267, 0.259,0.249,0.234,0.199
Pol. Eff.:
0.746,0.758,0.768,0.778,0.785,0.794,0.8,0.806,0.81,0.812,0.812,0.809,0.802, 0.794,0.777,0.75,0.676

The tasks are: (a) Fit a function to the Hp/N^2 v/s Qs/N data, and (b) Fit a function to the Pol. Eff. v/s Qs/N data.

This should be an interesting challenge problem for the “polynomial fitters” to try out.  The results of such efforts will, I assure you, be extremely revealing and should serve to reinforce my general point about use of polynomials for such work.

Needless to say, using a cubic spline to fit such data is a cinch and requires no guesswork at all about the type of function or polynomial order to use.  (As a rule, I specify that the slope at the end-points be equal to that of the neighboring points, so spline extrapolation is reasonable).

RE: Curve fitting pump curves

For your data with a very fine mesh and very smooth curve, spline is no doubt better. Since you have extracted the data from a curve, the best-fit has already been done for you in preparation of the curve. I would hope no noise would remain in the vendor's curve.

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Curve fitting pump curves

Actually there is a little bit of noise in the data. Perhaps based on the errors from reading off a curve.

I did a fit of 1000*hp/N^2  (Y) vs Qs/N  (X) here

http://154.6.66.55/electricpete/pumpcurvedata2.xls

The first estimate was Yhat1 which was a power law fit.

Residuals were R1.

I fit R1 using R1hat which was sum of two decaying exponentials

Result is the final estimate Yhat2 = Yhat1 + R1hat as follows:

Yhat2 =+$D$2+$D$3*(L9-$D$1)^$D$4 +$G$2*EXP(($G$1-L9)*$G$3)+$G$4*EXP(($G$1-L9)*$G$5)

where
column L is the X values
Constants as follows:
D1    0.499637258        G1    1.251207481
D2    0.297529817        G2    -0.039569578
D3    -0.176437834        G3    -17.73002606
D4    4.028911694        G4    0.053137589
            G5    -13.12118014

All of the calcs are done on the tabl leveled XYonly.

I had some difficulty fitting the right-most point until I increased it’s weighting.  (Objective function is the sum of all the residuals squared plus another 3 times that right-most point’s residual squared).

Results agree to within 1% for all the given data points.  (see cell H27 of tab XYonly)
You can see the overview of the results in the chart tab labeled “Y_andYhat2”

I think what you have is a smoother curve than if you used a spline for the same purpose.  Although I think the spline would probably be more accurate.    Also if there are known boundary conditions such as the curve crosses the axes at a right angle, it was indicated above that these could be incorporated into the spline method.  I don’t know of any easy way to incorporate that into the excel fit method I described. So I would think if you can incorporate those boundary conditions into the spline method, the spline method would be better for purposes of extrapolating beyond the data.

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Curve fitting pump curves

Actually if either/both the shutoff head or the 0-dp runout flow were known, these could easily be incoporated into the least square fit as another point....  using a similar weighting approach as mentioned above - heavier weighting of the points that are far away from the center cluster.

I don’t know enough about pumps to know if there are any other relevant boundary conditions.  (do the curves cross the axes at right anlges?)

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Curve fitting pump curves

Hi, electricpete:

First, many thanks for your intense effort.  I think your effort has demonstrated clearly that concocting a function of the type you used (power and exponential functions with 9 constants) would prove to be a painful undertaking for all but the most intrepid of practitioners.  Recall that this was necessary for a relatively innocuous data set.

As you suspected, I did read off the values manually from the manufacturer's graph, hence the slight jitter.  This is most often what users are reduced to doing as vendors invariably fail to provide either numerical values of the coordinates or an analytical equation for their curves.

The spline interpolation routines available in any decent software library can easily be required to provide any desired value for the terminal slope.  If you know what the stonewall flow is (recall that this is not really a measurable quantity), you should add it to the data set and then impose whatever terminal slope you desire.  Numerically, a very high slope would do the job, as a computer cannot "as yet" handle infinity!

For those curious about splines, the end conditions merely constrain the array of spline coefficient values so that the final curve behaves the way you would like it to, aesthetically speaking, if you were using a draftsman’s flexible strip to do the job.  

RE: Curve fitting pump curves

I would still say each method has it’s benefits.

Advantages of least squares fit:

- It gives a smoother curve in general

- If you do a real good job of least-square fitting a smooth curve to the noisy data, it will tend to be more accurate (not compared against the datapoints, but compared to the unknown “true” curve which presumably is a smooth curve through the "middle" of the data).

- It’s available in excel – a widely available program.  Spline is not (although I think if you select the “smoothed curve” option on a chart, excel draws a cubic spline through your data?  But the equation is not accessible?)

Advantage of cubic spline fit:

- If a terminal slope is known, that is an advantage to the spline method (can’t be put into the least squared method).

 - When I zoomed out I see the form of the my least-squares fit function extends fairly smoothly outwards to 0 DP units.  Unfortunately it only extends left to 0.5 in flow rate units due to my choice of fit functions.  I could have avoided this by eliminating the “zeroshifter” parameter in the first fit (yhat1).  Lesson learned by trial and error.

- Spline takes less time.  Once you decide cubic spline there are no more decisions to be made and the computer gives you the result.  The least-squared fit involves a little trial and error (although if you repeatedly fit pump curves I think it wouldn’t take long to know which terms to use).  It took me just a little over an hour to do this least squares fit (and I have never fit a pump curve before).  Then some additional time to write it up and post it on the board.

I have a question – what is innocous about this set of data compared to typical pump curve data? (bearing in mind the least-squares method is well equipped to deal with noisy data)

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Curve fitting pump curves

Also - 9 constants yes, but only 3 terms: a power-law and two exponentials

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Curve fitting pump curves

Hi, electricpete:

First, I believe that the curve fitting method anyone chooses is always subjective (driven by personal experience and belief).  My reason to enter this discussion was simply to mention that cubic splines are a much overlooked tool that can do the job very efficiently and without much controversy as to the quality of the fit.

I used the word innocuous to describe my data set simply because (a) despite not containing any weird kinks or discontinuities, it is incapable of being fit properly by any polynomial up to at least order 5, and (b) in all cases, the polynomial extrapolations are truly horrible.  This would be a nightmare for most people who don’t have much experience or knowledge beyond polynomial regression.

Your own very creditable effort demonstrated that concocting a well-behaved function took quite a bit of mathematical ingenuity, if not cunning.  I’m sorry to say that most engineers do not have the time or patience to do this, in my experience, even if they have the requisite skills.

RE: Curve fitting pump curves

Yes, that's a good point. I think if the data is smooth then spline will do the job very well and least-square curve fitting could waste a lot of time.  If done poorly, curve fitting gives worse results.  To do well may take a lot of time especially considering extrapolation requirements.

=====================================
Eng-tips forums: The best place on the web for engineering discussions.

RE: Curve fitting pump curves

You may be interested in the following curve-fitting program.  It has a table/spreadsheet interface, and data can easily be copied from excel.

http://curveexpert.webhop.biz/

RE: Curve fitting pump curves

(OP)
Hi again,

I have been monitoring the conversation between electricpete and UmeshMathur. It is interesting to note the outcome of trying to fit integer polynomial order fits to data versus spline fitting. I have also been evaluating the data I have been given and comparing the results of the various fits. I am still assessing which fit is "best" for my situation.

Essentially I have found through searching through our current database that I have 7 data points in general to describe a relationship between flow rate and head produced and flow rate and efficiency achieved. I initially hoped that I could increase the amount of data available to me and possibly get multiple sets of data for a single condition, thus increasing my confidence in the curve I eventually fit to the data. Since I am dealing with an entire database of over 200 pumps, retesting each pump is just not viable (My engineering manager would simply not allow for it). So for the immediate and practical future I have 7 data points to relate 2 variables, in general. These data points were extracted from the present pump selection software database. Due to a lack of documentation I also have no idea whether they represent test data or data that has been inputted from a previous published curve.

One of my fundamental problems is assessing the reliability of the data (but as the situation would have it – it is all that I have to work with).

I must ensure that the curves that are eventually fitted to the data represent the most likely relationship between the input and output variables in a typical situation for the pump. I have assumed points of inflection and turning points are not representative of the physical situation in terms of relating flow to both efficiency and pump head. Forgive the crude drawings but essentially the relationships are depicted below. Thus I do not expect any sudden turning points and so on within the data (as could happen when a 4th order polynomial tries to accommodate an outliner such as a datum with a high noise probability).  Sometime the 0 flow condition on Curve 1 does not correspond to the highest point on the curve (as depicted in this case by the 0 on Curve 1) and Curve 1 can “hump” at a certain flow rate, but the relationship should not produce any sudden points of inflection. The same holds true for Curve 2.  

      Ix
H    I       x
      I0           x
      I               x
      I                  x
      I                    x
      I                     x
      I_________________________                 (Curve 1)
            Flow

      I
Eff  I                x     
      I         x              x
      I      x                    x
      I                                                
      I   x                                                 
      I     
      Ix_____________________________    (Curve 2)
            Flow

I think one of my core problems is that I have little faith in the input data values. That does not mean to say I think they are non meaningful but rather I have no reliable assessment of just how meaningful they may or may not be. If I had reliable experimental data (including some error analysis) I would happily use the curve fitting techniques as well as the programs suggested to me by the respondents in this forum to obtain the “best” curve fit for my situation.

A good analogy would be to imagine that I was tasked to determine the relationship between Force and Deflection in a spring. For arguments sake I will assume that in reality the spring is linear under the range at which I am testing it, that is to say F=kX, where k is a constant – but at the start of investigation assume I do not know this fact and have been told that any relationship may apply. All I have to work with is a few data points. I have no idea of their origin or the test conditions (if they were derived experimentally as opposed to being read off a previously published graph). The relationship may be follows:



        I    
        I                                 x
        I
X      I                         x          
        I                  x
        I            x       
        I     x     
        I_________________________
                    F

I may know from university that the relationship should be linear over a certain range but I am told that in reality it is always slightly non linear. Since I find it hard to accept most things at face value without some proof I cannot accept this statement. If however I try to fit a curve (whether a polynomial, spline or other) to the data I find that, for arguments sake, a 3rd order polynomial fits the data very well (with corresponding R^2 vales and residual plots that look good). I find a linear fit does rather badly yet in my gut I simply feel it may more representative of the relationship.  In the end it may turn out that I was only partly correct and that a 2nd order curve could the most representative (albeit with an extremely small coefficient for the x^2 term). But from the given data I do not know. I could simply say that I will work with what I have been given and whatever curve appears to fit best, I will implement but I am apprehensive of causing trouble down the line, especially since my generated curves are to be used as masters that will be used as references in an extrapolation process that finds data between two curves.

Back to the pump curves, since I have little confidence in the data but testing the entire range of pumps is not feasible I have decided to try test one or two pumps thoroughly (up and down the curve etc) and obtain experimental data. I will apply what seems to be the most representative fit for the data, taking into account all the above comments and suggestions and then look at the resulting curve function. I hope to use that as a template for which to decide how the rest of the data should actually be fit and apply that fit type to the rest of the data set, regardless of what the initial data values say.

Of course in doing this I am assuming that the pumps behave hydraulically in the same manner (most of them are scaled versions of each other so it may be reasonable to assume this). I run the possible risk of forcing a fit to a data set that may in reality differ from the curve type I obtained from testing but for the time being I am willing to take that risk. If there are any major deviations for any of my pumps I will look at that pump individually and reassess the fit, possibly even testing that pump (if allowed).

Does that sound reasonable?

Regards,

Patrick02

P.S: All The curve fitting programs I have downloaded are obviously similar in function but the one suggested by spciesla is quite user friendly.

RE: Curve fitting pump curves

(OP)
Just a quick afterthought; I have data for multiple curves for the same pump as each curve as given in Curve 1 represents data at a particular speed setting and typically I have between 4 and 7 speed settings per pump. Thus I could convert all my equations to a non dimensional form and that would give me a lot more data on one curve. I could then try determine the best fit through that non-dimensional curve and use it as a template to generate my dimensional curves. This follows on from UmeshMathur’s comments on using the fan laws. The affinity laws are well known for certain parameters such as flow rate and pump head or pressure but I have never attempted to generate the non dimensional form of the relationship between flow and efficiency with respect to pump speed. It should be an interesting exercise to see if it is worthwhile or not. It’s worth a shot and I will know at the end whether it is useful to me or not. This does not alleviate my lack of faith in the data but it is another alternative to consider        

RE: Curve fitting pump curves

Patrick02:

A few suggestions:
(1) The general shapes of both your curves are quite reasonable.

(2) Since you have data for several RPMs, I think you really have no choice but to use the non-dimensional form of the curves.  By the way, for applying the similarity/affinity laws, the proper non-dimensional quantities when impeller diameter might also be a variable are: (a) Qs/N/D^3, and (b)g*Hp/N^2/D^2, where Qs=actual suction volumetric flow[m^3/s], N=speed[revolutions/s], D=diameter[m], g=gravitational constant[Newtons/kgf], and Hp=polytropic head[m*kgf/kg], in consistent metric units.  (I thought I would mention the diameter in case some of your test data was for geometrically similar but different size pumps.  Also, the quantities mentioned in previous posts are NOT dimensionless).

(2) Rest assured that the non-dimensional form for polytropic efficiency v/s Qs/N/D^3 is theoretically sound for pumps as well as compressors.

(3) Remember that these laws were developed for Newtonian fluids.  Their use for mildly non-Newtonian fluids is OK, provided the viscosity is not very large and you use the proper effective viscosity.  Another dimensionless quantity to check with varying fluid properties is (c) rho*N*D^2/mu, where rho=fluid density, and mu=dynamic viscosity.

(4) In case you see a systematic deviation in the non-dimensional curves for a particular RPM relative to the others, you might suspect an RPM measurement error (quite common) and consider applying an adjustment to the quoted RPM to make the values line up better with the rest of the curves.

(5) All in all, I feel that a spline is probably your best bet in the circumstances you have described with a limited amount of data.  Before doing any curve fits, however, do make sure that the curves for different RPMs were not generated from the curve for one RPM simply by using the affinity laws.  The same check should be applied for different pump diameters.

(6) Lastly, since you are dealing with slurries, there could easily be some uncertainty as to how representative the measured values really might be, especially in case slurry density changes as a result of a differing solids content.

Do advise us when you have finished your assignment, as the results will be interesting to all who have participated in this thread.

Have fun, and cheers.

RE: Curve fitting pump curves

U-M, great comments about fitting.  I will have to remember this when my next such task comes along.

Patrick, I thnk that early on in this discussion you hit on one of the root issues.  It isn't the type of fit that is used, it is the uncertanty in the data.  If you re-test the same pump over and over you will get different values (not to mention testing different pumps of the same model).  That is what you need to explore.  It is likely that some of your most basic measurements (shaft speed or pressure) are not nearly accuraate enough to support the conclusions that you are trying to make.

= = = = = = = = = = = = = = = = = = = =
Rust never sleeps
Neither should your protection
http://www.trent-tube.com/contact/Tech_Assist.cfm

RE: Curve fitting pump curves

(OP)
Hi to all those who helped me during my query on this thread.

My apologies with the delay on my reply.  I thought I would give you a concluding synopsis of what I had to do to fit these curves.

My immediate involvement in the project is now finished and it has been a learning curve for me.

Soon after the conclusion of this discussion I was informed by the marketing manager that the software selection program used by us only handles 4th order polynomial fits when generating data and nothing else. There was neither time nor money available to change this during the scope of the project. Why such a constraint was placed in the program in the first place is beyond me. Thus I was limited to using a 4th order fit to fit our data. The focus of the program as communicated to me and hence my assigned brief was to “fix” our database for the immediate future while a more comprehensive test program on the entire range of pumps was to be implemented. The input data source for the fixed database was our current published curve book. Essentially I had to ensure that what was in our selection program matched up to what was given in our curve book. My only guess was by inference this meant that what was in our curve book fairly represented the operational characteristics of the pump through customer witness tests (even though test data was not available to me). I devised a program to take digital images of our published curves, scale them and grab data points off them. Those points would be taken as what was seen in the curve book. My job was to fit 4th order polynomials through those points and generate an equation for the polynomial, which the software selection program could use. The whole process appeared a little bizarre to me but nevertheless it was my assigned brief. I could not successfully fit a 4th order curve through all my points, as discussed earlier, so I devised a program which would “add error” to certain points far to the left and right of the best efficiency point along a certain speed curve in order to increase the accuracy of the fit in the region close to the best efficiency point. After this was done for all speed curves on a particular pump curve I generated a so called “selectable region” where sales personal could safely make selections on a pump with errors between what they saw in a curve book and what the program told them of less than 1%. In regions outside of that the error increased dramatically. For example on some speed curves I was forced to add an error of up to 25% to the 0 flow condition (even though this is nonsensical) to make it fit well near the region of best efficiency. Thus if a salesman selected a pump at a flow condition near zero the program would tell him that he had a pump that was 35% efficient when in reality it was only 10% or less.

In this manner I managed to produce a short term solution to our problem. I did not, however, like the fundamental approach to this problem and method of “playing” with our data by adding and subtracting error in order to get a software program to work. Thus I documented where I added error and how I constructed the new curves so that sales personal were fully aware of the situation when making selections. I have since advocated a new testing program for our pumps with more documentation and control to avoid this type of situation. I may have fixed a problem but I did not consider it engineering. The one benefit to all of this is that we now have a lot more information on other aspects of our curves, which is documented electronically. This will help me when we eventually supplement this “fixed” database with up to date and controlled test data to increase our confidence in the data.

Thank you for all your help. This has been a learning experience for me.  
            

Kind regards,

Patrick


RE: Curve fitting pump curves

Sorry for the late entry into what is apparently a closing discussion, but I would like to know one thing.

Why is it important to know the efficiency of a pump that is flowing at near zero?  Is anybody really going to operate a pump for extended periods of time at a near zero flowrate.  Surely that the operator expects to pay a heavy penalty in efficiency when flows go lower than 50% BEP.  Should not the object be to develop highly accurate efficiency predictions at reasonable flowrates, or failing that, design a better pump station configuration?

   Going the Big Inch! worm
http://virtualpipeline.spaces.msn.com

RE: Curve fitting pump curves

(OP)
That is perfectly correct BigInch. There is no reason for any of our customers to wish to operate a pump far off its design point. Anything less than 60% QBEP or greater than 120% is considered a no go area. This served as the basis for my reasoning that I could decrease the accuracy in the regions < 60% or >120% QBEP in order to improve the accuracy of the fit though the region within 60% to 120% QBEP since no one would wish to operate outside the region anyway. The fundimental problem I had with this had more to do with the methodology behind the data in that I should not have been grabbing it from a curve book but rather be using reliable and controlled test data. We have reasonable efficiency predictions at reasonable flowrates but I was hoping that we will in time gather more test data on all our curves.

Regards,

Patrick  

RE: Curve fitting pump curves

Patric,
Thanks for that confirmation.  I couldn't imagine what all the fuss was about.  I thought my petroleum concentric experience had blinded me to some strange slurry application that required extended running off BEP.  IMO you made a good decision to concentrate your study towards the reasonable flow ranges.  I don't consider any efficiency data outside the margins you mention as significant for operating purposes, and start-up/shut-downs are such limited times, it makes no difference what's happening there.  If operating in those regions, there are problems that need looking after far worse than determining pump efficiency to the 4th order.

I always limit my pumps to running between 50% and your same 120% and will change a multiple pump configuration to accomodate operation within those ranges at all times.  If flow conditions change, I will try changing the process operation first, but if that's impossible I will try ripping out the existing pumps and get some that will operate within the proper margins.  The money saved in one year is considerable in relation to pump costs when pumps are operated at proper efficiencies and arn't being destroyed by poor operating practices.I

No matter what the published curves say or what performance is quoted, we will do shop flow tests to develop the as-built curves and then test in the field as proof of ultimate performance acceptability, so usually I get pretty good data, eventually.

For simulating pump performance, I will generally use 10 points, mainly because it makes an easy-read graph, but I will go to 16 points if extras are warranted for flowrates important to the objective of the simulation.

What I would like to be able to get to 4 signficant figures are the dry and wet moments of inertias.  Where's my wish list?

   Going the Big Inch! worm
http://virtualpipeline.spaces.msn.com

RE: Curve fitting pump curves

(OP)
BigInch,

What you are doing sounds quite awesome and even though I work for a pump manufacturer, were I to work for a process company I would also perform field tests of published curves simply because normally it is not that difficult to do and can save a lot of time and money should a problem be picked up. Gaining information on dry and wet moments of inertia seems interesting and would help were you able to conduct tests to get those figures. My particular field of expertise (even though I am fairly new to the game) is vibration analysis. What I want to study is the effect of sump level on the vibration response of our long vertical spindle pumps. Currently no one knows just how different slurry types and densities affect the critical speeds of these pumps. All quotes on pump critical speeds for the vertical spindle range are based on "dry" calculations. Due to the damping effect slurries have on the rotordynamics of the system they are bound to lower the critical speeds.

Just going back to what you mentioned on good operating practices and cost savings, you may be interested to check out an article called "How to Justify Machinery Improvements Using Reliability Engineering Principles." It is also about how best practices can lower the life cycle costs of the pump and the authors use reliability engineering simulations to test the theory. It is quite an interesting read. Go to :

http://www.barringer1.com/Papers.htm

and scroll down to article 13. and download it. Hope you enjoy the article

Regards,

Patrick  

RE: Curve fitting pump curves

Pat, I know the Barringer website and did "pick up" that doc. Its a good site.  I have pretty well cleaned it out.  I like simulations that incorporate Monte Carlo techniques and have gleened a lot of good info from them. (#15)

I'm doing some consulting off/on for a mideast state oil co where the pumps are generally high capacity, high pressure and quite high cost, so they get them with all the bells and whistles.  I've found that even the biggest pump mfgrs sometimes get it wrong and have seen a few cases where the as-built curves are a percent or two or three under the guoted curves ... over the entire Q range, and that was not for slurries.  I can only imagine what you're trying to deal with.  To expect them to be 4th order accurate way off BEP, just seemed to be too much to ask for in the first place, and in the second, not of much practical use at the end of the day.  Critical speeds of shafts, now that's useful.  Sounds like a high gray matter project.  I don't know much about slurries, but from what I have studied, they seem to suffer some radical behavior changes with what seems like tiny variations in densities and flows too.  At one point they're liquids, then shift to almost solids.... ugh!  There must be some really strange effective mass dampening going on around the shafts when the density, surface speed, and flowrates change.

Would you mind dropping me an e-mail (profile secton on my webspace page)?

   Going the Big Inch! worm
http://virtualpipeline.spaces.msn.com

Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Eng-Tips forums is a member-only feature.

Click Here to join Eng-Tips and talk with other members!


Resources