CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Strange Result of ERCOFTAC Centrifugal Pump MRF case in OF2.2.x

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By linnemann

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 7, 2013, 21:03
Default Strange Result of ERCOFTAC Centrifugal Pump MRF case in OF2.2.x
  #1
Member
 
Masashi Ohbuchi
Join Date: Oct 2009
Posts: 74
Rep Power: 17
Ohbuchi is on a distinguished road
Hi,

Recently, I've tried ERCOFTAC Centrifugal Pump MRF case (ECPGgi2D) for OpenFOAM-2.2.x, and encountered strange results as attached images.
For comparison, I've also attached the results for OF2.1.x.
Results for OF2.1.x are in good agreement with original published results (OF-1.5-dev).

For OF2.2.x case, I used simpleFoam with following fvOptions.

-------fvOptions ----------------------------------------------------
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

MRF1
{
type MRFSource;
active true;
selectionMode cellZone;
cellZone rotor;
nonRotatingPatchs ( INLET OUTLET GGI_EXT GGI_INT BLADE_ROT);

MRFSourceCoeffs
{
origin (0 0 0);
axis (0 0 -1);
omega 209;
}
}
---------------------------------------------------------------------

If someone has the solution to correct this problem, please let me know.
Any comment would be appreciated.
Attached Images
File Type: png U_2.1.x.png (61.1 KB, 73 views)
File Type: png U_2.2.x.png (56.4 KB, 70 views)
File Type: png Wr_021x.png (7.7 KB, 58 views)
File Type: png Wr_022x.png (24.8 KB, 42 views)
Ohbuchi is offline   Reply With Quote

Old   August 12, 2013, 06:17
Default
  #2
Senior Member
 
linnemann's Avatar
 
Niels Nielsen
Join Date: Mar 2009
Location: NJ - Denmark
Posts: 556
Rep Power: 27
linnemann will become famous soon enough
Hi

I have no problems with commit 2.2.x-b6d83360ae0e

I have uploaded the working case here and removed the reliance on the old turbo stuff. Only 2.2 BC's.

LINK TO CASE
Attached Images
File Type: jpg 2013-08-12-13-23-000053.jpg (48.4 KB, 51 views)
Ohbuchi likes this.
__________________
Linnemann

PS. I do not do personal support, so please post in the forums.
linnemann is offline   Reply With Quote

Old   August 12, 2013, 20:48
Default
  #3
Member
 
Masashi Ohbuchi
Join Date: Oct 2009
Posts: 74
Rep Power: 17
Ohbuchi is on a distinguished road
Hi, linnemann.
I tried case file of you, I was able to get same result as you.
But for now, I could not find any important difference in your dictionary files with mine.
I'm going to continue to investigate in detail.

Best regards,

Masashi Ohbuchi
Ohbuchi is offline   Reply With Quote

Old   September 10, 2013, 05:18
Default
  #4
New Member
 
Brahim Aakti
Join Date: Mar 2009
Location: Lucerne, Switzerland
Posts: 23
Rep Power: 17
brahim is on a distinguished road
Hi,

here some comments on your fvOptions:
- there is a typo in the nonRotatingPatchs -> nonRotatingPatches
- I would not declare the patch BLADE_ROT as non-rotating

Here is the fvOptions I'm using:
Code:
MRF1
{
    type            MRFSource;
    active          true;
    selectionMode   cellZone;
    cellZone        rotor;

    MRFSourceCoeffs
    {
        nonRotatingPatches (INLET OUTLET GGI_EXT GGI_INT);
        origin      (0 0 0);
        axis        (0 0 -1);
        omega       209;
    }
}
Best Regards - Brahim.
brahim is offline   Reply With Quote

Old   December 6, 2013, 06:14
Default
  #5
Member
 
Join Date: Aug 2013
Posts: 50
Rep Power: 13
nash is on a distinguished road
hi everyone,

are the nonRotatingPatches needed to be included in the dict?
in my case (3D radial fan) the mrfZone is a cylinder and within the cylinder there are wall (from the fan blade)

so should i add the wall within the mrfZone in the nonRotatingPatches or just leave it without nonRotatingPatches?

thanks
regards,
nash
nash is offline   Reply With Quote

Old   December 6, 2013, 10:20
Default
  #6
Member
 
Join Date: Aug 2013
Posts: 50
Rep Power: 13
nash is on a distinguished road
i have tried without including the patch and with the patch included in fvOption.

the right one is not to include the patch within the mrfZone as nonRotatingPatches.
Hope this helps.

regards,
nash

Last edited by nash; December 16, 2013 at 10:59.
nash is offline   Reply With Quote

Old   December 13, 2013, 04:48
Default Help
  #7
New Member
 
Minh-VietNam
Join Date: Oct 2013
Posts: 7
Rep Power: 13
mingbn is on a distinguished road
Quote:
Originally Posted by linnemann View Post
Hi

I have no problems with commit 2.2.x-b6d83360ae0e

I have uploaded the working case here and removed the reliance on the old turbo stuff. Only 2.2 BC's.

LINK TO CASE
Hi
I'm a new openfoam user. I am interested in your case. (i'm studying about centrifugal pump). I tried your case but i can't get results the same your.
(i'm using ubuntu 13.04; openfoam 2.2.2)
Please show to me.
Thank you so much.
Attached Images
File Type: jpg Screenshot from 2013-12-13 15:44:06.jpg (46.9 KB, 33 views)
File Type: jpg Screenshot from 2013-12-13 15:45:15.jpg (45.8 KB, 27 views)
mingbn is offline   Reply With Quote

Old   January 15, 2014, 23:46
Default ercoftac Centrifugal Pump for OF2.2.x
  #8
Member
 
Masashi Ohbuchi
Join Date: Oct 2009
Posts: 74
Rep Power: 17
Ohbuchi is on a distinguished road
Thank you, Brahim.
And I appologize for my delayed response.
You were correct!.
My essential mistake was the setting of nonrotatingPaches.
I've attached case files adapted to OF2.2.x.
All validation cases; ECPGgi2D(MRF+Ggi), ECPStitchedMesh2D(MRF+StitchMesh)
and ECPMixerGgiFvMesh2D(Sliding+Ggi) will be working fine.
The Allrun script in the postProcessing directory also works.
In these archive, because of the limitation of maximum size, mesh files(rotor2D and stator2D) and measurement data are not included.

Best regards,
M.Ohbuchi
Attached Files
File Type: gz ercoftacCentrifugalPump1.tar.gz (27.0 KB, 40 views)
File Type: gz ercoftacCentrifugalPump2.tar.gz (30.9 KB, 39 views)
Ohbuchi is offline   Reply With Quote

Old   April 30, 2015, 06:03
Default
  #9
Senior Member
 
M. Montero
Join Date: Mar 2009
Location: Madrid
Posts: 155
Rep Power: 17
be_inspired is on a distinguished road
Hi,

Thank you for adapting the case to OF2.2.x. I am trying to use OF2.3.x for my case but with some problems.
I have some doubts:
It seems like now it is not needed to make faceZones. It is enough with cellZone definition performed inside topoSet. And now, inside OF code faceZone selection is done automatically.

Ohbuchi has defined BLADE_ROT within patches subDict. I thought that is not needed because by default all patches included in the MRF zone were rotated.

I am trying to simulate a rotor with a fixed part of the geometry but there are problem in the rotor surface. k, epsilon and so give very extrange values ( +12 or so at iteration 5).
For nonRotatingPatches I have added only AMI1 and AMI2 interfaces.

Mesh was generated for y+ around 100 so I am using nutkWallFunction.

I hope someone can help me.
be_inspired is offline   Reply With Quote

Old   May 3, 2015, 06:01
Default
  #10
Member
 
Masashi Ohbuchi
Join Date: Oct 2009
Posts: 74
Rep Power: 17
Ohbuchi is on a distinguished road
HI,

Please check and see attached tar archive. It'll work with OF2.3.
As you already mentioned above, cellZone is only needed for defininig MRF zone.
And faceZone was not used in my case files.
And I also think nutkWallFunction is more suitable for wall boundary for k.

Ohbuchi

Quote:
Originally Posted by be_inspired View Post
Hi,

Thank you for adapting the case to OF2.2.x. I am trying to use OF2.3.x for my case but with some problems.
I have some doubts:
It seems like now it is not needed to make faceZones. It is enough with cellZone definition performed inside topoSet. And now, inside OF code faceZone selection is done automatically.

Ohbuchi has defined BLADE_ROT within patches subDict. I thought that is not needed because by default all patches included in the MRF zone were rotated.

I am trying to simulate a rotor with a fixed part of the geometry but there are problem in the rotor surface. k, epsilon and so give very extrange values ( +12 or so at iteration 5).
For nonRotatingPatches I have added only AMI1 and AMI2 interfaces.

Mesh was generated for y+ around 100 so I am using nutkWallFunction.

I hope someone can help me.
Attached Files
File Type: gz simpleFoam.tar.gz (13.2 KB, 23 views)
Ohbuchi is offline   Reply With Quote

Old   May 10, 2015, 04:47
Default
  #11
Senior Member
 
Huang Xianbei
Join Date: Sep 2013
Location: Yangzhou,China
Posts: 302
Rep Power: 14
huangxianbei is on a distinguished road
Quote:
Originally Posted by linnemann View Post
Hi

I have no problems with commit 2.2.x-b6d83360ae0e

I have uploaded the working case here and removed the reliance on the old turbo stuff. Only 2.2 BC's.

LINK TO CASE
Dear Niels Nielsen:
Thank you for your sharing. But I can't open the link, could you please send a copy to my email znhuang@163.com?

Xianbei
huangxianbei is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
release of the ERCOFTAC centrifugal pump - Fourth OpenFOAM Workshop olivier OpenFOAM 8 October 29, 2018 08:49
Is Playstation 3 cluster suitable for CFD work hsieh OpenFOAM 9 August 16, 2015 15:53
strange behaviour of GGI in parallel on axis symmetrical case A.Devesa OpenFOAM Running, Solving & CFD 0 April 6, 2010 04:58
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 13:24
centrifugal pump john Main CFD Forum 0 April 22, 2004 07:13


All times are GMT -4. The time now is 08:09.