[olle]

Dear Foamers,

As a newbie, I've been working with OpenFoam 2.0.x a couple of months and thought it was enough until I realized I needed cyclic boundary conditions on the faces surrounding my domain in spanwise and normal directions. I'm working with incompressible fluid passing a cylinder in free stream.

Have created my hex mesh in ICEM, exported it as a Fluent 6 mesh, imported and exported it again in Gambit until I finally converted it with fluent3DMeshToFoam. Tried to use cyclic boundaries first which were available in 2.0.x. Got problems with faces that were not matching, so I started to look at the cyclicGgi boundary condition which are implemented in the 1.5-dev and 1.6-ext versions.

It looked promising at first when I tried with 1.5-dev. Coupled the boundaries that were set up as cyclicGgi in my constant/polyMesh/boundary file (see below):

front
{
type cyclicGgi;
nFaces 52299;
startFace 6038409;
shadowPatch new_back;
zone front_zone;
bridgeOverlap false;
rotationAxis (0 1 0);
rotationAngle 0;
separationOffset (0 0 0.0628317);
}
new_back
{
type cyclicGgi;
nFaces 52299;
startFace 6090708;
shadowPatch front;
zone new_back_zone;
bridgeOverlap false;
rotationAxis (0 1 0);
rotationAngle 0;
separationOffset (0 0 0.0628317);
}
bottom
{
type cyclicGgi;
shadowPatch top;
zone bottom_zone;
bridgeOverlap false;
rotationAxis (1 0 0);
rotationAngle 0;
separationOffset (0 -0.360 0);
nFaces 11388;
startFace 6143007;
}
top
{
type cyclicGgi;
shadowPatch bottom;
zone top_zone;
bridgeOverlap false;
rotationAxis (1 0 0);
rotationAngle 0;
separationOffset (0 -0.360 0);
nFaces 11388;
startFace 6154395;
}

All preprocessing steps went fine, i.e. I was able to execute the setSet, setsToZones and decomposePar commands. However, I ran into trouble since I wanted to use pimpleFoam which is not implemented in 1.5-dev. So, after installation of 1.6-ext and some minor changes to system/fvSolution to make it compatible with 1.6-ext I started with the preprocessing steps again. Run into problems already in the setSet step:

Create time

Create polyMesh for time = 0

Initializing the GGI interpolator between master/shadow patches: front/new_back

From function void GGIInterpolation<MasterPatch, SlavePatch>::rescaleWeightingFactors() const
in file lnInclude/GGIInterpolationWeights.C at line 544
Uncovered faces found. On master: 52299 on slave: 52299


--> FOAM FATAL ERROR:
Found uncovered faces for GGI interface front/new_back while the bridgeOverlap option is not set in the boundary file.
This is an unrecoverable error. Aborting.

From function void ggiPolyPatch::calcPatchToPatch() const
in file meshes/polyMesh/polyPatches/constraint/ggi/ggiPolyPatch.C at line 221.

FOAM aborting

I guess my main questions are:

1. Should I setup my boundaries in constant/polyMesh/boundary differently in 1.5-dev and 1.6-ext?

2. If not, are there other differences between the two OF versions in how they handle the cyclicGgi boundary condition?

Best regards
Olle

[declanj]

Hi All,

I am having the same problems. A simple test case with 2 cyclicGgi patches works perfectly in 1.5-dev but returns uncovered faces when using 1.6-ext.

Has anyone found a solution to this problem?

[yvyan]

Same problem!

Did somebody solve?