cyclicGgi discontinuity problems

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

October 19, 2009, 04:55
cyclicGgi discontinuity problems
#1
New Member

Martin Becker
Join Date: Sep 2009
Posts: 5
Rep Power: 10
I am using OpenFoam to simulate the waterflow in an axial turbine runner. I constructed a mesh for one blade channel and tried to simulate it using the cyclicGgi interface (OpenFoam 1.5-dev rev. 1238 as suggested by the SIG Turbomachinery).

Unfortunately, I see errors on the cyclic interface. Pressure as well as velocity is discontinous on the interface, that leads to unphysical results with bad convergence.

I'm using the simpleSRFFoam solver, but I also saw the same problems in my wicked gate simulation using simpleFoam.

I attached some pictures and the output of the simulation:
- cyclicGgi_channel.jpg: the four blades and the computational domain. cyclicGgi interfaces are coloured green and red

- cyclicGgi_channel_2.jpg: one side of the interface is rotated onto the other -> surfaces are congruent

- cyclicGgi_channel_problem.jpg: flow calculated with cyclicGgi interface, see the discontinuities at the interface (pressure on cylindrical cut)

- cyclicGgi_360deg.jpg: same case, same boundary conditions, but calculated with 360° mesh (without cyclicGgi) -> very good results (here: pressure on cylindrical cut)

simpleSRFFoam_out.txt: simulation output

My Ggi weighting factors look fine:
Initializing the GGI interpolator between master/shadow patches:
perio1/perio2
Evaluation of GGI weighting factors:
Largest slave weighting factor correction : 3.33067e-16 average:
6.04589e-17
Largest master weighting factor correction: 1.26787e-13 average:
4.43631e-16

Here is my cyclicGgi patch definition in (constant/polymesh/boundary):
perio1
{
type cyclicGgi;
nFaces 4600;
startFace 614376;
zone perio1_zone;
bridgeOverlap off;
rotationAxis (0 0 1);
rotationAngle 90;
separationOffset (0 0 0);
}
perio2
{
type cyclicGgi;
nFaces 4600;
startFace 618976;
zone perio2_zone;
bridgeOverlap off;
rotationAxis (0 0 1);
rotationAngle 90;
separationOffset (0 0 0);
}

What is the current status of cyclicGgi?
Are the problems I described common ("known bug")?
What about enableGgiNonOrthogonalCorrection? Could that help me?

Any help is greatly appreciated!

Martin
Attached Images
 cyclicGgi_360deg.jpg (90.8 KB, 99 views) cyclicGgi_channel.jpg (47.9 KB, 86 views) cyclicGgi_channel_2.jpg (61.4 KB, 83 views) cyclicGgi_channel_problem.jpg (91.4 KB, 76 views)
Attached Files
 simpleSRFFoam_out.txt (35.9 KB, 29 views)

October 19, 2009, 15:17
#2
Senior Member

Martin Beaudoin
Join Date: Mar 2009
Posts: 332
Rep Power: 15
Hello Martin,

One little clarification is necessary here:

Quote:
 Originally Posted by mbecker I am using OpenFoam to simulate the waterflow in an axial turbine runner. I constructed a mesh for one blade channel and tried to simulate it using the cyclicGgi interface (OpenFoam 1.5-dev rev. 1238 as suggested by the SIG Turbomachinery).
The recommendation to stick to revision 1238 is uniquely related to the ERCOFTAC Centrifugal Pump test case published on the Wiki, and the usage of MRFSimpleFoam.

This has nothing to do with the GGI, and you should upgrade your 1.5-dev installation to the latest 1.5-dev release available in order to get the best GGI implementation available.

The problem lies with MRFSimpleFoam that basically stopped working for the ECP test case when we upgraded to a version newer that rev. 1238.

My group is working on a different MFR implementation. See here for a little sneak peek preview:
http://openfoamwiki.net/index.php/Si...leTurboMFRFoam

We hope to release simpleTurboMFRFoam in a near future.

Regards,

Martin

 March 12, 2010, 13:01 #3 Member   Nick Gardiner Join Date: Apr 2009 Location: Chichester, UK Posts: 93 Rep Power: 10 Hi Martin Is this the case with turbDyMFoam as I have a similar problem where the fluid does not cross the interface Regards Nick

 January 27, 2011, 00:14 #4 Member   Kurne Join Date: Aug 2010 Location: Pune, INDIA Posts: 88 Rep Power: 9 Dear All Will anyone please tell me whether the GGI is available with the OpenFOAM 1.7.1 or not.Also the simpleTurboMRFFoam will be Run with the OpenFOAM 1.7.1 or not. __________________ Simulation Is Determination Of Imagination Towards Approximation ® Best Regards Mubeen K Kurne

February 25, 2011, 13:18
#5
Member

Claudio
Join Date: Mar 2010
Posts: 57
Rep Power: 9
Quote:
 Originally Posted by mbecker I am using OpenFoam to simulate the waterflow in an axial turbine runner. I constructed a mesh for one blade channel and tried to simulate it using the cyclicGgi interface (OpenFoam 1.5-dev rev. 1238 as suggested by the SIG Turbomachinery). Unfortunately, I see errors on the cyclic interface. Pressure as well as velocity is discontinous on the interface, that leads to unphysical results with bad convergence. I'm using the simpleSRFFoam solver, but I also saw the same problems in my wicked gate simulation using simpleFoam. I attached some pictures and the output of the simulation: - cyclicGgi_channel.jpg: the four blades and the computational domain. cyclicGgi interfaces are coloured green and red - cyclicGgi_channel_2.jpg: one side of the interface is rotated onto the other -> surfaces are congruent - cyclicGgi_channel_problem.jpg: flow calculated with cyclicGgi interface, see the discontinuities at the interface (pressure on cylindrical cut) - cyclicGgi_360deg.jpg: same case, same boundary conditions, but calculated with 360° mesh (without cyclicGgi) -> very good results (here: pressure on cylindrical cut) simpleSRFFoam_out.txt: simulation output My Ggi weighting factors look fine: Initializing the GGI interpolator between master/shadow patches: perio1/perio2 Evaluation of GGI weighting factors: Largest slave weighting factor correction : 3.33067e-16 average: 6.04589e-17 Largest master weighting factor correction: 1.26787e-13 average: 4.43631e-16 Here is my cyclicGgi patch definition in (constant/polymesh/boundary): perio1 { type cyclicGgi; nFaces 4600; startFace 614376; shadowPatch perio2; zone perio1_zone; bridgeOverlap off; rotationAxis (0 0 1); rotationAngle 90; separationOffset (0 0 0); } perio2 { type cyclicGgi; nFaces 4600; startFace 618976; shadowPatch perio1; zone perio2_zone; bridgeOverlap off; rotationAxis (0 0 1); rotationAngle 90; separationOffset (0 0 0); } What is the current status of cyclicGgi? Are the problems I described common ("known bug")? What about enableGgiNonOrthogonalCorrection? Could that help me? Any help is greatly appreciated! Martin

Dear Martin,

I am investigating a case similar to Yours, that is a blade channel of a radial pump with seven blades.
So I have to use cyclicGgi patches (together with some ggi between the
rotor and the stator part).
I use the 1.5-dev version.

The problem is that when I decompose my case, the following error
occurs.

/*---------------------------------------------------------------------------*\
| ========= |
|
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
| \\ / O peration | Version: 1.5-dev
|
| \\ / A nd | Revision: 1861
|
| \\/ M anipulation | Web: http://www.OpenFOAM.org
|
\*---------------------------------------------------------------------------*/
Exec : decomposePar
Date : Feb 25 2011
Time : 16:20:37
Host : claudio-laptop
PID : 6163
Case : /root/OpenFOAM/root-1.5-dev/run/pompa
nProcs : 1

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
* * //
Create time

Time = 0
Create mesh for region region0

Calculating distribution of cells
Selecting decompositionMethod simple

Finished decomposition in 0.02 s

Calculating original mesh data

Distributing cells to processors

Distributing faces to processors
Preserving global face zone Per2Zone
Preserving global face zone Per1Zone
Preserving global face zone OutInter2Zone
Preserving global face zone OutInter1Zone
Preserving global face zone InInter2Zone
Preserving global face zone InInter1Zone

Distributing points to processors

Constructing processor meshes
Segmentation fault

I run the setSet -batch setBatch and setsToZones -noFlipMap
applications before decomposing the domain, and they work properly.

The case runs in serial mode, so there aren't macro problems as regards the settings.

Could You kindly help me to understand the reason for this problem?

Claudio

P.S.: Does the cyclicGGi b.c.,( nowaday, in Your opinion) work properly or not?

 February 27, 2011, 12:36 #6 New Member   Martin Becker Join Date: Sep 2009 Posts: 5 Rep Power: 10 Dear Claudio, as far as I can see from your text, your procedure looks OK. You can be sure that cyclicGgi works as it should. Your segmentation fault can have various reasons. Did you try to debug it? Maybe it is useful to post more details of your case (boundary file etc.). Did you try another version of OpenFoam as well (e.g. 1.6-ext)? Martin

 Tags cyclicggi

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Mechstud Main CFD Forum 4 July 26, 2011 12:13 Micha Siemens 0 August 6, 2003 13:55 Ivan de Miguel Urain FLUENT 1 November 16, 2001 10:36 Alan Kang Main CFD Forum 1 February 23, 2000 00:46 Aleksey Alekseev Main CFD Forum 0 May 12, 1999 15:38

All times are GMT -4. The time now is 07:37.