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

cyclicGgi discontinuity problems

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 19, 2009, 04:55
Default cyclicGgi discontinuity problems
  #1
New Member
 
Martin Becker
Join Date: Sep 2009
Posts: 5
Rep Power: 7
mbecker is on a distinguished road
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
Attached Images
File Type: jpg cyclicGgi_360deg.jpg (90.8 KB, 94 views)
File Type: jpg cyclicGgi_channel.jpg (47.9 KB, 84 views)
File Type: jpg cyclicGgi_channel_2.jpg (61.4 KB, 78 views)
File Type: jpg cyclicGgi_channel_problem.jpg (91.4 KB, 73 views)
Attached Files
File Type: txt simpleSRFFoam_out.txt (35.9 KB, 20 views)
mbecker is offline   Reply With Quote

Old   October 19, 2009, 15:17
Default
  #2
Senior Member
 
Martin Beaudoin
Join Date: Mar 2009
Posts: 330
Rep Power: 13
mbeaudoin will become famous soon enough
Hello Martin,

One little clarification is necessary here:

Quote:
Originally Posted by mbecker View Post
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
mbeaudoin is offline   Reply With Quote

Old   March 12, 2010, 13:01
Default
  #3
Member
 
Nick Gardiner
Join Date: Apr 2009
Location: Chichester, UK
Posts: 93
Rep Power: 8
NickG is on a distinguished road
Hi Martin

Is this the case with turbDyMFoam as I have a similar problem where the fluid does not cross the interface

Regards

Nick
NickG is offline   Reply With Quote

Old   January 27, 2011, 00:14
Default
  #4
Member
 
Kurne
Join Date: Aug 2010
Location: Pune, INDIA
Posts: 88
Rep Power: 6
kurne is on a distinguished road
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
kurne is offline   Reply With Quote

Old   February 25, 2011, 13:18
Default
  #5
Member
 
Claudio
Join Date: Mar 2010
Posts: 57
Rep Power: 7
claco is on a distinguished road
Quote:
Originally Posted by mbecker View Post
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

Calculating processor boundary addressing

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?

Thank You in advance.


Claudio


P.S.: Does the cyclicGGi b.c.,( nowaday, in Your opinion) work properly or not?
claco is offline   Reply With Quote

Old   February 27, 2011, 12:36
Default
  #6
New Member
 
Martin Becker
Join Date: Sep 2009
Posts: 5
Rep Power: 7
mbecker is on a distinguished road
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
mbecker is offline   Reply With Quote

Reply

Tags
cyclicggi

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
Needed Benchmark Problems for FSI Mechstud Main CFD Forum 4 July 26, 2011 12:13
Some problems with Star CD Micha CD-adapco 0 August 6, 2003 13:55
Airpak: Problems with heat flow contour conditions Ivan de Miguel Urain FLUENT 1 November 16, 2001 10:36
Stepwise Discontinuity problem Alan Kang Main CFD Forum 1 February 23, 2000 00:46
Inverse problems Aleksey Alekseev Main CFD Forum 0 May 12, 1999 15:38


All times are GMT -4. The time now is 14:12.