CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

release of the ERCOFTAC centrifugal pump - Fourth OpenFOAM Workshop

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By olivier

LinkBack Thread Tools Display Modes
Old   May 13, 2009, 13:00
Default release of the ERCOFTAC centrifugal pump - Fourth OpenFOAM Workshop
olivier Petit
Join Date: Mar 2009
Location: Göteborg, Sweden
Posts: 67
Rep Power: 10
olivier is on a distinguished road

The ERCOFTAC centrifugal pump Case Study has recently been updated.
The centrifugal pump with a vaned diffuser problem is a testcase of the ERCOFTAC Turbomachinery Special Interest Group.
It will be presented at the Fourth OpenFOAM Workshop in Montreal (1-4 June 2009). Prof. Marina Ubaldi, from the Università di Genova, allowed us to use the measurements data for distribution and set-up of the analysis within OpenFOAM.

You can find on the wiki the following informations:

* How to check out the files from the SourceForge svn.
* How to generate and simulate the cases.
* Different informations on some useful utilities (mergeMesh, stitchMesh, MRFSimpleFoam, GGi interface, and even more ...).
* Automatic post-processing and validation using sample and gnuplot.
* Automatic plotting of residuals using foamLog and gnuplot.

The case study can be found in the svn, in the following branch:

You can find the ERCOFTAC Centrifugal Pump Case Study at:

Please help us improve this Case Study by giving feedback or by contributing!

Best regards,
Olivier, Maryse, Håkan and Martin.

zhujj09 likes this.
olivier is offline   Reply With Quote

Old   May 28, 2009, 04:26
Default ERCOFTAC centrifugal pump test case - problem
New Member
Ivana B.
Join Date: Mar 2009
Posts: 4
Rep Power: 10
ivana is on a distinguished road
Hello everyone,

I have been running ERCOFTAC centrifugal pump test case (running with no changes and with Allrun) and I face the following problem:
In the Time = 10 after Solving Uy I get “Floating point exception”. Actually I have this problem since I have made update to revision 1266. I have returned back the revision 1241 and it is working fine. I compiled the both versions in the same way and I have not noticed any problems. I have tried on different hardware and OS platforms (but the same compiler). I always get the same error.

Anyone else experience the same, or I should further search the problem by us?
Thanks for your reply.

ivana is offline   Reply With Quote

Old   May 28, 2009, 08:31
olivier Petit
Join Date: Mar 2009
Location: Göteborg, Sweden
Posts: 67
Rep Power: 10
olivier is on a distinguished road
Hello Ivana,

There is indeed a problem with the latest version of the svn, hopefully it will be solved soon.
What I would recommend is to stay for now at the working version (svn 1240) until this is solved. However, it is possible to make the test cases converge with the latest update of the dev version, but you have to play with some parameters a bit: the reason why the simulations are stopping is that k and epsilon are bounding above the acceptable limit, so the simulation stopps. A bounding k and epsilon is usually quite common at the begining of a simulation, as a lot happens at that time, but it should stay under a certain level. If it does not, as it is now, you can choose to increase your under-relaxation, so that you take only a small amount of the solution into account.

So here, if you want it to work, you can put in system/fvSolution a under-relaxing parameter of 0.5 for k and epsilon, and from what I have seen, it should work. However, the convergence is not great.

Ultimately, I recommend to stay at a working svn version.

olivier is offline   Reply With Quote

Old   July 13, 2009, 05:11
David Hora
Join Date: Mar 2009
Location: Zürich, Switzerland
Posts: 63
Rep Power: 10
david is on a distinguished road
Hi all

Thanks a lot for this interesting test case. At the moment I'm testing the GGI with turbDyMFoam on a similar, but 3D case. This 2D case will be interesting for me as I can get results in a much shorter time.

I have two questions concerning the ECPGgi2D case:

1) checkMesh reports the following error:

***Number of edges not aligned with or perpendicular to non-empty directions: 24545
<<Writing 49090 points on non-aligned edges to set nonAlignedEdges

The simulation seems to work without problems but could it be that the non-aligned edges cause difficulties in 2D cases? Or is this error negligible and not critical?

2) I've seen that a new solver called simpleTurboMFRFoam was developed. What was the reason for this? Will it have any advantages over MRFSimpleFoam?

Best regards,
david is offline   Reply With Quote

Old   June 25, 2013, 03:54
Join Date: Mar 2009
Posts: 85
Rep Power: 10
husker is on a distinguished road

Thanks for building the case study and contributing to it.

Although I'm not allowed to employ OpenFOAM in my office, I'm deeply interested in a centrifugal pump case study.

Could anyone provide geometry and performance data for which to be used in commercial CFD codes such as FLUENT or Star-CCM. Following to my assessment, I will be glad to take your attention and discuss the results.

husker is offline   Reply With Quote

Old   June 26, 2014, 11:10
Default Reference pressure - CHALMERS study
Lisandro Maders
Join Date: Feb 2013
Posts: 83
Rep Power: 6
Lisandro Maders is on a distinguished road

I am currently performing the validation case of ERCOFTAC centrifugal pump with vane diffusers. Looking at a CHALMERS pdf, I found the way they calculated the Cp coefficient was using the standard equation:


The p0 is the static pressure at the suction pipe. They made an assumption of such value by trying to obtain similar levels of Cp as the experimental results of Ubaldi had. The value they found is 700 Pascal.

What I am wondering is if this value is in absolute or Gauge pressure. Also, if they are using an incompressible solver in OpenFOAM, shouldn't they be using the pressure in m^2/s^2 (pressure/roh) ?? Ok, maybe they multiplied the pressure value taken from the simulation by the density and then figured out the Cp coefficient.
Anyway, the absolute/gauge pressure issue keeps in my mind.

Really thankful for any help.

Regards, Lisandro Maders
Lisandro Maders is offline   Reply With Quote

Old   October 6, 2015, 04:48
Default How is ensemble average achieved?
Senior Member
Huang Xianbei
Join Date: Sep 2013
Location: Yangzhou,China
Posts: 287
Rep Power: 7
huangxianbei is on a distinguished road
Hi, olivier:
I read your paper "Numerical Investigations of Unsteady Flow in a Centrifugal
Pump with a Vaned Diffuser" and found that the calculation results are ensemble-averaged as I read from the paper "Fig-ures 3 and 4 show the ensemble-averaged experimental and
calculated radial and tangential velocities in the radial gap at ...". So did you use the same procedure in the experiment (phase-locked sampling and ensemble-average)?

huangxianbei is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Fourth OpenFOAM workshop in Montreal page OpenFOAM 4 June 18, 2009 14:54
Turbomachinery at OpenFOAM Workshop Milan 2008 hani OpenFOAM 0 March 21, 2008 03:33
Second OpenFOAM Workshop in Zagreb Croatia 79Jun2007 hjasak OpenFOAM 5 June 10, 2007 12:33
OpenFOAM Training and Workshop Zagreb 2628Jan2006 hjasak OpenFOAM 1 February 2, 2006 22:07
OpenFOAM Training and Workshop Hrvoje Jasak Main CFD Forum 0 October 7, 2005 07:14

All times are GMT -4. The time now is 03:06.