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

Ggi

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 4, 2010, 08:55
Default Ggi
  #1
Member
 
Join Date: May 2010
Posts: 33
Rep Power: 16
FabOr is on a distinguished road
Hi Foamers,

i´ve got some strange results using ggi-patches => see attached picture and plot-over-line, y-axis.

I use ggi-patches for minX/maxX and minY/maxY similar to a simple construction. This construction looks like a cuboid.
So the values of p,U,... of the patch minX (= master) should be transferred to maxX (=slave), same for minY/maxY. Correct me if i am wrong.
But in the plot you can see that the ggi-patch doesnt work like i want: the values of U and p on both sides dont coincide.

For better understanding and looking for mistakes here is my boundary-file:

Code:
    minX
    {
        type            ggi;
        nFaces          6549;
        startFace       182291;
    shadowPatch    maxX;
    rotationAxis    (0 0 1);
    rotationAngle    0;
    zone        minX_Zone;
    bridgeOverlap    on;//on
    separationOffset    (0 0 0);
    }
    maxX
    {
        type            ggi;
        nFaces          6460;
        startFace       188840;
    shadowPatch    minX;
    rotationAxis    (0 0 1);
    rotationAngle    0;
    zone        maxX_Zone;
    bridgeOverlap    on;//on
    separationOffset    (0 0 0);
    }
    minY
    {
        type            ggi;
        nFaces          600;
        startFace       195300;
    shadowPatch    maxY;
    rotationAxis    (0 0 1);
    rotationAngle    0;
    zone        minY_Zone;
    bridgeOverlap    on;//on
    separationOffset    (0 0 0);
    }
    maxY
    {
        type            ggi;
        nFaces          570;
        startFace       195900;
    shadowPatch    minY;
    rotationAxis    (0 0 1);
    rotationAngle    0;
    zone        maxY_Zone;
    bridgeOverlap    on;//on
    separationOffset    (0 0 0);
    }
The patches in the U,p,..-files are described as type ggi.

Maybe this kind of mistake can be related to this error-warning running simpleFoam:

Code:
Reading/calculating face flux field phi

Initializing the GGI interpolator between master/shadow patches: minX/maxX
Evaluation of GGI weighting factors:

    From function void GGIInterpolation<MasterPatch, SlavePatch>::rescaleWeightingFactors() const
    in file /home/fabian/OpenFOAM/OpenFOAM-1.5-dev/src/OpenFOAM/lnInclude/GGIInterpolationWeights.C at line 534
    Uncovered faces found.  On master: 6549 on slave: 6460
  Largest slave weighting factor correction : 0 average: 0
  Largest master weighting factor correction: 0 average: 0

Initializing the GGI interpolator between master/shadow patches: minY/maxY
Evaluation of GGI weighting factors:

    From function void GGIInterpolation<MasterPatch, SlavePatch>::rescaleWeightingFactors() const
    in file /home/fabian/OpenFOAM/OpenFOAM-1.5-dev/src/OpenFOAM/lnInclude/GGIInterpolationWeights.C at line 534
    Uncovered faces found.  On master: 600 on slave: 570
  Largest slave weighting factor correction : 0 average: 0
  Largest master weighting factor correction: 0 average: 0

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kEpsilon
Any help and guidance is much appreciated! If you need more files/information to help me to solve my problem feel free to ask!

regards
Attached Images
File Type: jpg screen1.jpg (19.0 KB, 72 views)
Attached Files
File Type: gz screen-1.png.tar.gz (97.5 KB, 18 views)
FabOr is offline   Reply With Quote

Old   November 4, 2010, 09:27
Default
  #2
Senior Member
 
Chris Sideroff
Join Date: Mar 2009
Location: Ottawa, ON, CAN
Posts: 434
Rep Power: 22
cnsidero is on a distinguished road
Ya, in the simpleFoam output, it's telling you their are uncovered faces (actually all of them), meaning the ggi wasn't able to match up the patches.

One question, did you create the face zones (minX_zone, maxX_zone, minY_zone, maxY_zone)?
cnsidero is offline   Reply With Quote

Old   November 4, 2010, 09:38
Default
  #3
Member
 
Join Date: May 2010
Posts: 33
Rep Power: 16
FabOr is on a distinguished road
Hi,

thx for your answer cnsidero!

i´ve created the zones by: 1) faceSet for each patch minX/maxX/...
2) setsToZones -noFlipMap

best regards
FabOr is offline   Reply With Quote

Old   November 4, 2010, 09:49
Default
  #4
Senior Member
 
Chris Sideroff
Join Date: Mar 2009
Location: Ottawa, ON, CAN
Posts: 434
Rep Power: 22
cnsidero is on a distinguished road
Ok, I was just checking.

I too have had this happen to me before but I was using cyclicGgi and didn't have the rotation angles correct. How about removing the unnecessary stuff, e.g.:

minX
{
type
ggi;
shadowPatch maxX;
bridgeOverlap false;
zone minX_zone;
}

I'm not sure ... just grasping.
cnsidero is offline   Reply With Quote

Old   November 4, 2010, 10:09
Default
  #5
Member
 
Join Date: May 2010
Posts: 33
Rep Power: 16
FabOr is on a distinguished road
Hi,

i´ve deleted the ones i dont need => nothing changed
and when i set bridgeOverlap to off => case doesnt start because of uncovered faces.

regards
FabOr is offline   Reply With Quote

Old   November 6, 2010, 12:46
Default
  #6
Senior Member
 
Martin Beaudoin
Join Date: Mar 2009
Posts: 332
Rep Power: 22
mbeaudoin will become famous soon enough
Looks like you are using GGI patches separated by some translation offset.

You need to set the separationOffset to something else than (0 0 0).

That value would be the translation offset necessary for the 2 GGI patches to overlap each other.

Martin

Quote:
Originally Posted by FabOr View Post
Hi Foamers,

i´ve got some strange results using ggi-patches => see attached picture and plot-over-line, y-axis.

I use ggi-patches for minX/maxX and minY/maxY similar to a simple construction. This construction looks like a cuboid.
So the values of p,U,... of the patch minX (= master) should be transferred to maxX (=slave), same for minY/maxY. Correct me if i am wrong.
But in the plot you can see that the ggi-patch doesnt work like i want: the values of U and p on both sides dont coincide.

For better understanding and looking for mistakes here is my boundary-file:

Code:
    minX
    {
        type            ggi;
        nFaces          6549;
        startFace       182291;
    shadowPatch    maxX;
    rotationAxis    (0 0 1);
    rotationAngle    0;
    zone        minX_Zone;
    bridgeOverlap    on;//on
    separationOffset    (0 0 0);
    }
    maxX
    {
        type            ggi;
        nFaces          6460;
        startFace       188840;
    shadowPatch    minX;
    rotationAxis    (0 0 1);
    rotationAngle    0;
    zone        maxX_Zone;
    bridgeOverlap    on;//on
    separationOffset    (0 0 0);
    }
    minY
    {
        type            ggi;
        nFaces          600;
        startFace       195300;
    shadowPatch    maxY;
    rotationAxis    (0 0 1);
    rotationAngle    0;
    zone        minY_Zone;
    bridgeOverlap    on;//on
    separationOffset    (0 0 0);
    }
    maxY
    {
        type            ggi;
        nFaces          570;
        startFace       195900;
    shadowPatch    minY;
    rotationAxis    (0 0 1);
    rotationAngle    0;
    zone        maxY_Zone;
    bridgeOverlap    on;//on
    separationOffset    (0 0 0);
    }
The patches in the U,p,..-files are described as type ggi.

Maybe this kind of mistake can be related to this error-warning running simpleFoam:

Code:
Reading/calculating face flux field phi

Initializing the GGI interpolator between master/shadow patches: minX/maxX
Evaluation of GGI weighting factors:

    From function void GGIInterpolation<MasterPatch, SlavePatch>::rescaleWeightingFactors() const
    in file /home/fabian/OpenFOAM/OpenFOAM-1.5-dev/src/OpenFOAM/lnInclude/GGIInterpolationWeights.C at line 534
    Uncovered faces found.  On master: 6549 on slave: 6460
  Largest slave weighting factor correction : 0 average: 0
  Largest master weighting factor correction: 0 average: 0

Initializing the GGI interpolator between master/shadow patches: minY/maxY
Evaluation of GGI weighting factors:

    From function void GGIInterpolation<MasterPatch, SlavePatch>::rescaleWeightingFactors() const
    in file /home/fabian/OpenFOAM/OpenFOAM-1.5-dev/src/OpenFOAM/lnInclude/GGIInterpolationWeights.C at line 534
    Uncovered faces found.  On master: 600 on slave: 570
  Largest slave weighting factor correction : 0 average: 0
  Largest master weighting factor correction: 0 average: 0

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kEpsilon
Any help and guidance is much appreciated! If you need more files/information to help me to solve my problem feel free to ask!

regards
mbeaudoin is offline   Reply With Quote

Old   November 8, 2010, 10:47
Default
  #7
Member
 
Join Date: May 2010
Posts: 33
Rep Power: 16
FabOr is on a distinguished road
hi,

thanks for your answer, martin.

my boundary-file:

Code:
minX
    {
        type            ggi;
        nFaces          5511;
        startFace       279799;
        shadowPatch    maxX;
        zone        minX_Zone;
        bridgeOverlap    on;
        separationOffset    (60e-06 0 0);
    }

    maxX
    {
        type            ggi;
        nFaces          5514;
        startFace       285310;
        shadowPatch    minX;
        zone        maxX_Zone;
        bridgeOverlap    on;
        separationOffset    (60e-06 0 0);
    }

    minY
    {
        type            ggi;
        nFaces          372;
        startFace       290824;
        shadowPatch    maxY;
        zone        minY_Zone;
        bridgeOverlap    on;
        separationOffset    (600e-06 0 0);
    }

    maxY
    {
        type            ggi;
        nFaces          372;
        startFace       291196;
        shadowPatch    minY;
        zone        maxY_Zone;
        bridgeOverlap    on;
       separationOffset    (600e-06 0 0);
i´ve also tried to set the separationOffset to:

minX: (60e-06 0 0)
maxX: (-60e-06 0 0)
minY: (0 600e-06 0)
maxY: (0 -600e-06 0)

but the problem remains. i dont get the same values of p,U on opposite sites. the uncovered faces message when starting simpleFoam also remains.

best regards
FabOr is offline   Reply With Quote

Old   November 9, 2010, 10:43
Default
  #8
Member
 
Join Date: May 2010
Posts: 33
Rep Power: 16
FabOr is on a distinguished road
Hi,

i solved the problem.

best regards
FabOr is offline   Reply With Quote

Old   November 9, 2010, 15:10
Default
  #9
Senior Member
 
Chris Sideroff
Join Date: Mar 2009
Location: Ottawa, ON, CAN
Posts: 434
Rep Power: 22
cnsidero is on a distinguished road
Can you share with us what your solution was?
cnsidero is offline   Reply With Quote

Old   November 10, 2010, 04:26
Default
  #10
Member
 
Join Date: May 2010
Posts: 33
Rep Power: 16
FabOr is on a distinguished road
Of course

i switched the patch type from ggi to cyclicGgi. then i´ve done the same procedure like Thomas in his test-case. some uncovered faces still remain, but that refers to my .stl-file.
as you can see in the attached plot cyclicGgi works fine!

look here for further guidance:

http://www.cfd-online.com/Forums/ope...cyclicggi.html

best regards
Attached Images
File Type: jpg plot_cyclicGgi.jpg (17.6 KB, 50 views)
FabOr is offline   Reply With Quote

Old   January 19, 2012, 07:50
Default
  #11
Member
 
Aqua
Join Date: Oct 2011
Posts: 96
Rep Power: 15
aqua is on a distinguished road
Quote:
Originally Posted by FabOr View Post
Of course

i switched the patch type from ggi to cyclicGgi. then i´ve done the same procedure like Thomas in his test-case. some uncovered faces still remain, but that refers to my .stl-file.
as you can see in the attached plot cyclicGgi works fine!

look here for further guidance:

http://www.cfd-online.com/Forums/ope...cyclicggi.html

best regards
Hello, I have similar problem with you, so could you please help me on this? my error is:
Uncovered faces found. On master: 4800 on slave: 4800


--> FOAM FATAL ERROR:
Found uncovered faces for GGI interface interfacei/interfaceo
This is an unrecoverable error. Aborting.

From function void overlapGgiPolyPatch::calcPatchToPatch() const
in file meshes/polyMesh/polyPatches/constraint/overlapGgi/overlapGgiPolyPatchGeometry.C at line 186.

FOAM aborting


my geometry is as seen in the attachment, the two cubes in side stand for two cars, moving towards each other. so i only set the interface between the two parts as overlapGgi, others as patch, the icube.stl and ocube.stl as wall.

I am not sure I should change the interface type to cyclicGgi, or remain it as overlapGgi..

Thank you so much!

Aqua
Attached Images
File Type: jpg two blocks with two cubes.jpg (20.0 KB, 40 views)
aqua is offline   Reply With Quote

Old   January 19, 2012, 10:04
Default
  #12
Member
 
Timo K.
Join Date: Feb 2010
Location: University of Stuttgart
Posts: 66
Rep Power: 16
timo_IHS is on a distinguished road
Hallo,

I'd take "only" ggi.
Did you create faceZones of your patches?

Best
timo_IHS is offline   Reply With Quote

Old   January 19, 2012, 12:16
Default
  #13
Member
 
Aqua
Join Date: Oct 2011
Posts: 96
Rep Power: 15
aqua is on a distinguished road
Quote:
Originally Posted by timo_IHS View Post
Hallo,

I'd take "only" ggi.
Did you create faceZones of your patches?

Best
Hey,
Thank you so much for your suggestion! I changed it to ggi, when perform "decomposePar", it worked for a while then there is the error:
"Attempt to cast type genericPatch to type lduInterface
From function refCast<To>(From&) in file /bb/civ/liuyu/OpenFOAM/OpenFOAM-1.6-ext/src/OpenFOAM/lnInclude/typeInfo.H at line 115.

FOAM aborting"

yes, I had faceZones of the interfaces which are ggi type.

Could you please give more advice about the new error?
Thank you so much!
Aqua
aqua is offline   Reply With Quote

Old   January 19, 2012, 12:41
Default
  #14
Member
 
Timo K.
Join Date: Feb 2010
Location: University of Stuttgart
Posts: 66
Rep Power: 16
timo_IHS is on a distinguished road
Do you have
globalFaceZones (ggiPatchA_zone ggiPatchB_zone)
in your decomposeParDict?

and do you use GAMG? at a first starting point I wouldn't, cause there might be other problems.
timo_IHS is offline   Reply With Quote

Old   January 19, 2012, 13:53
Default
  #15
Member
 
Aqua
Join Date: Oct 2011
Posts: 96
Rep Power: 15
aqua is on a distinguished road
Quote:
Originally Posted by timo_IHS View Post
Do you have
globalFaceZones (ggiPatchA_zone ggiPatchB_zone)
in your decomposeParDict?

and do you use GAMG? at a first starting point I wouldn't, cause there might be other problems.
I changed something, now when run "decomposePar", there was no "foam aborting", but really "End".

Then I perform "icoDyMFoam", it's now simulating.... but i can see there is still some information like "uncovered faces found"... but it's still simulating without aborting..

Do you think it's right?

Thank you!

Aqua
aqua is offline   Reply With Quote

Old   January 20, 2012, 04:53
Default
  #16
Member
 
Timo K.
Join Date: Feb 2010
Location: University of Stuttgart
Posts: 66
Rep Power: 16
timo_IHS is on a distinguished road
Could you precise that?
timo_IHS is offline   Reply With Quote

Old   January 20, 2012, 06:24
Default
  #17
Member
 
Aqua
Join Date: Oct 2011
Posts: 96
Rep Power: 15
aqua is on a distinguished road
Quote:
Originally Posted by timo_IHS View Post
Could you precise that?
Good morning!
I want to tell that I use the tutorial case turbopassageRotating.
after some modification for my own case, I run setSet -batch setBatch; regionCellSets; setsToZones -noFliiMap.... Then I run decomposePar, OF didn't abort, but give some information like "unconvered faces found"(althought there was information like this, please notice that OF didn't abort). So after that,


I run icoDyMFoam, OF didn't abort, and now I got two files: 0.01 0.02, as in attachment.
Then I run paraFoam,Error accured like this:

Cannot find 'value' entry on patch interfacei of field p in file "/bb/civ/li
which is required to set the values of the generic patch field.

could you please help on this?
Attached Images
File Type: jpg 12.jpg (6.4 KB, 29 views)
aqua is offline   Reply With Quote

Old   May 9, 2013, 11:19
Default
  #18
Member
 
Dogan
Join Date: Nov 2012
Location: Bochum/Germany
Posts: 42
Rep Power: 14
dogan is on a distinguished road
Hi Fabor,

I am also getting error messages with my case, and it is somehow similar with yours. i have 2 cyclicGgi patches which contain uncovered faces, and there supposed to be a relative movement between those patches. my simulation starts, but crashes with a floating point exception error message after a while. (i can see the mesh is rotating like 1 degree). i am not sure what exactly seperationOffset is, but how should i calculate this value?

and who do you mean when you say you have done the same precedure like Thomas? can you please describe the steps one more time?
thanks and regards
Dogan

Quote:
Originally Posted by FabOr View Post
Of course

i switched the patch type from ggi to cyclicGgi. then i´ve done the same procedure like Thomas in his test-case. some uncovered faces still remain, but that refers to my .stl-file.
as you can see in the attached plot cyclicGgi works fine!

look here for further guidance:

http://www.cfd-online.com/Forums/ope...cyclicggi.html

best regards
dogan 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
Floating point exception error Alan OpenFOAM Running, Solving & CFD 11 July 1, 2021 22:51
Difference between ggi and overlapGgi? GGI Tips and Tricks? philippose OpenFOAM Running, Solving & CFD 7 January 16, 2013 10:40
GGI in OpenFOAM-1.5-dev philippose OpenFOAM Running, Solving & CFD 14 November 13, 2011 15:55
Problem using GGI besto OpenFOAM 13 October 30, 2010 08:34
GGI bug?? rieuk OpenFOAM Running, Solving & CFD 2 April 25, 2010 11:47


All times are GMT -4. The time now is 00:28.