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

Cyclic vs ggi vs directMapped Patches

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 4, 2009, 16:05
Default
  #21
Senior Member
 
Martin Beaudoin
Join Date: Mar 2009
Posts: 332
Rep Power: 14
mbeaudoin will become famous soon enough
Hello,
  • Check the sign of your offset values.
  • Make sure you are running with the latest version available from openfoam-extend. (At least with the svn revision 1505).

Martin

Quote:
Originally Posted by mauscher View Post
Hello,

I'm just writing my thesis and I have a problem with the cyclic ggi boundary. First I want described what I've done and then what the problem is.

1. createPatch -> sidewallright and sidewallleft cyclicGgi
from the boundary-file ->

sidewallleft
{
type cyclicGgi;
nFaces 10176;
startFace 3583008;
shadowPatch sidewallright;
zone zonesidewallleft;
bridgeOverlap true;
rotationAxis (0 0 1);
rotationAngle 0;
separationOffset (0 0 0.075);
}
sidewallright
{
type cyclicGgi;
nFaces 10176;
startFace 3593184;
shadowPatch sidewallleft;
zone zonesidewallright;
bridgeOverlap true;
rotationAxis (0 0 1);
rotationAngle 0;
separationOffset (0 0 -0.075);
}

2. faceSet from Patches sidewallright and sidewallleft

3. setsToZones -noFlipMap

4. decomposePar (2 processors)

5. mpirun ....

6. From the log-file->

Evaluation of GGI weighting factors:

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

The calculation works, but one half of the sidewallpatch is calculated wrong. So where is the problem?


\\mauscher
mbeaudoin is offline   Reply With Quote

Old   December 8, 2009, 06:27
Default
  #22
Member
 
Join Date: Sep 2009
Posts: 45
Rep Power: 9
AirS is on a distinguished road
Hello,

I'm struggling with directMapped. Below are my 0/U and changeDictionary files:
*************************************
INLET
{
type directMappedFixedValue;
value uniform (1 0 0);
setAverage true;
average (1 0 0);
}
*************************************
*************************************
dictionaryReplacement
{
boundary
{
INLET
{
type directMappedPatch;
offset (0.76 0 0);
sampleRegion region0;
sampleMode nearestCell;
samplePatch none;
}
}
}
*************************************
Then when I run pisoFoam, I've got this message error:

*************************************
Time = 0.0005

Courant Number mean: 0 max: 0.0202144



gradientInternalCoeffs cannot be called for a genericFvPatchField (actual type directMappedFixedValue)
on patch INLET of field U in file "/nbu/eee/soulierr/OpenFOAM/fineTrainLESCyclicDirectMapped/0/U"
You are probably trying to solve for a field with a generic boundary condition.

From function genericFvPatchField<Type>::gradientInternalCoeffs( ) const
in file fields/fvPatchFields/basic/generic/genericFvPatchField.C at line 782.

FOAM exiting
**************************************

I do not know how to sort this issue out ...I had a look at the source code mentionned above but It didn't help me. I hope someone will be able to help me dealing with.
Thank you in advance.
AirS is offline   Reply With Quote

Old   December 8, 2009, 07:21
Default
  #23
Member
 
olivier Petit
Join Date: Mar 2009
Location: Göteborg, Sweden
Posts: 67
Rep Power: 9
olivier is on a distinguished road
Hello Mauscher,

From the message that you get from the GGI weighting factors, it looks as you have 5088 faces uncovered on your GGI. If you look at the number of cells that you have in your interface, which seems to be 10176, it seems that half your GGI interface is not covered.
I believe that is why you have a problem with your GGI interface. Have look at it, something is not right there.


Olivier

Quote:
Originally Posted by mauscher View Post
Hello,

I'm just writing my thesis and I have a problem with the cyclic ggi boundary. First I want described what I've done and then what the problem is.

1. createPatch -> sidewallright and sidewallleft cyclicGgi
from the boundary-file ->

sidewallleft
{
type cyclicGgi;
nFaces 10176;
startFace 3583008;
shadowPatch sidewallright;
zone zonesidewallleft;
bridgeOverlap true;
rotationAxis (0 0 1);
rotationAngle 0;
separationOffset (0 0 0.075);
}
sidewallright
{
type cyclicGgi;
nFaces 10176;
startFace 3593184;
shadowPatch sidewallleft;
zone zonesidewallright;
bridgeOverlap true;
rotationAxis (0 0 1);
rotationAngle 0;
separationOffset (0 0 -0.075);
}

2. faceSet from Patches sidewallright and sidewallleft

3. setsToZones -noFlipMap

4. decomposePar (2 processors)

5. mpirun ....

6. From the log-file->

Evaluation of GGI weighting factors:

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

The calculation works, but one half of the sidewallpatch is calculated wrong. So where is the problem?


\\mauscher
olivier is offline   Reply With Quote

Old   December 8, 2009, 18:43
Default
  #24
New Member
 
Join Date: Apr 2009
Location: Finspång, Sweden
Posts: 4
Rep Power: 9
mauscher is on a distinguished road
Hi,

thanks for your help oliver, but I already noticed that. The problem is, if I run the simulation on just one processor all faces are uncovered. If I decompose the case for two processors, as it is mentioned here, one half of the face is uncovered. Thats quite strang. Right now, I try to run the extionsions from sf and hope that the problem will be solved.

\\mauscher
mauscher is offline   Reply With Quote

Old   January 4, 2010, 20:03
Default
  #25
New Member
 
Kyungrock Kim
Join Date: Jun 2009
Posts: 20
Rep Power: 9
Kr_kim is on a distinguished road
Hi Airs

You need to set type, 'directMapped' instead of directMappedFixedValue.
"directMappedFixedValue" is not correct type following error message.

Hope It be helpful for you

kkr
Quote:
Originally Posted by AirS View Post
Hello,

I'm struggling with directMapped. Below are my 0/U and changeDictionary files:
*************************************
INLET
{
type directMappedFixedValue;
value uniform (1 0 0);
setAverage true;
average (1 0 0);
}
*************************************
*************************************
dictionaryReplacement
{
boundary
{
INLET
{
type directMappedPatch;
offset (0.76 0 0);
sampleRegion region0;
sampleMode nearestCell;
samplePatch none;
}
}
}
*************************************
Then when I run pisoFoam, I've got this message error:

*************************************
Time = 0.0005

Courant Number mean: 0 max: 0.0202144



gradientInternalCoeffs cannot be called for a genericFvPatchField (actual type directMappedFixedValue)
on patch INLET of field U in file "/nbu/eee/soulierr/OpenFOAM/fineTrainLESCyclicDirectMapped/0/U"
You are probably trying to solve for a field with a generic boundary condition.

From function genericFvPatchField<Type>::gradientInternalCoeffs( ) const
in file fields/fvPatchFields/basic/generic/genericFvPatchField.C at line 782.

FOAM exiting
**************************************

I do not know how to sort this issue out ...I had a look at the source code mentionned above but It didn't help me. I hope someone will be able to help me dealing with.
Thank you in advance.
Kr_kim is offline   Reply With Quote

Old   January 5, 2010, 07:14
Default
  #26
Member
 
Join Date: Sep 2009
Posts: 45
Rep Power: 9
AirS is on a distinguished road
Thanks Kr_kim,
I changed it like you said. But I had the previous error which deals with an ordering problem .
Thanks though.
AirS is offline   Reply With Quote

Old   January 23, 2012, 07:13
Default GGi
  #27
Member
 
Aqua
Join Date: Oct 2011
Posts: 96
Rep Power: 7
aqua is on a distinguished road
Quote:
Originally Posted by mbeaudoin View Post
Hello Jens,

Please provide your complete test case.

Martin
Hello Martin,
I am trying to perform a simulation about two cars passing by each other. Please find attached picture about the mesh.
so there are two cubes standing for two cars. iblock and oblock, containing the two cars(cubes), will move towards each other. So, interfaces between iblock and oblock is defined as ggi. I used the tutorial turboPassengeRotating. the boundary file is like this:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

14
(
imaxz
{
type patch;
nFaces 400;
startFace 1464469;
}
iminz
{
type patch;
nFaces 400;
startFace 1464869;
}
iminy
{
type patch;
nFaces 400;
startFace 1465269;
}
imaxy
{
type patch;
nFaces 400;
startFace 1465669;
}
interfacei
{
type ggi;
nFaces 400;
startFace 1466069;
rotationAxis (0 0 1);
bridgeOverlap true;
// nCopies 12;
shadowPatch interfaceo;
zone interfacei_faces;
}
imaxx
{
type patch;
nFaces 400;
startFace 1466469;
}
icube_ICUBE
{
type wall;
nFaces 72018;
startFace 1466869;
}
ominy
{
type patch;
nFaces 703;
startFace 1538887;
}
omaxy
{
type patch;
nFaces 400;
startFace 1539590;
}
ominz
{
type patch;
nFaces 400;
startFace 1539990;
}
omaxz
{
type patch;
nFaces 400;
startFace 1540390;
}
interfaceo
{
type ggi;
nFaces 400;
startFace 1540790;
rotationAxis (0 0 1);
bridgeOverlap true;
// nCopies 12;
shadowPatch interfacei;
zone interfaceo_faces;
}
omaxx
{
type patch;
nFaces 400;
startFace 1541190;
}
ocube_ocube
{
type wall;
nFaces 58255;
startFace 1541590;
}
)

two blocks with two cubes1.JPG
I run
" setSet -batch setBatch
regionCellSets
setsToZones -noFlipMap" there is no error.
then i run decomposePar, OF didn't abort, but there is some information like “uncovered faces found”.
then i run icoDyMFoam, on Friday, until today(Monday), obviously there are some time steps file created (please see another picture in attachment, and OF is still running).
but when i run paraFoam, there is the error:
Cannot find 'value' entry on patch interfacei of field p in file "/bb/civ/liuyu/OpenFOAM/liuyu-1.6-ext/run/Ggicase1/0/p"
which is required to set the values of the generic patch field.
(Actual type ggi)

Please add the 'value' entry to the write function of the user-defined boundary-condition


file: /bb/civ/liuyu/OpenFOAM/liuyu-1.6-ext/run/Ggicase1/0/p::boundaryField::interfacei from line 66 to line 66.

From function genericFvPatchField<Type>::genericFvPatchField(con st fvPatch&, const Field<Type>&, const dictionary&)
in file genericFvPatchField/genericFvPatchField.C at line 71.

FOAM exiting

could you please, help me on this? Thank you so much!!!!!
Attached Images
File Type: jpg 13.jpg (20.9 KB, 15 views)
aqua is offline   Reply With Quote

Old   January 23, 2012, 08:15
Default
  #28
Member
 
olivier Petit
Join Date: Mar 2009
Location: Göteborg, Sweden
Posts: 67
Rep Power: 9
olivier is on a distinguished road
Hi aqua,

What version of OpenFOAM are you using when you open paraFoam? To me it seems that you are using a none-extended version, which is why it is unable to recognize the patch name.
olivier is offline   Reply With Quote

Old   January 23, 2012, 10:31
Default
  #29
Member
 
Aqua
Join Date: Oct 2011
Posts: 96
Rep Power: 7
aqua is on a distinguished road
Quote:
Originally Posted by olivier View Post
Hi aqua,

What version of OpenFOAM are you using when you open paraFoam? To me it seems that you are using a none-extended version, which is why it is unable to recognize the patch name.
hi, Olivier,
Thank you so much for your reply!
I do use OF1.6ext...

Please help me out..

Aqua
aqua is offline   Reply With Quote

Old   January 30, 2012, 07:34
Default
  #30
Member
 
Aqua
Join Date: Oct 2011
Posts: 96
Rep Power: 7
aqua is on a distinguished road
Quote:
Originally Posted by jens_klostermann View Post
Hi Francesco,

yes!
I just read something that only works for planar faces and I figured out that it is possible to give OpenFoam-1.5 the choice whether the cyclic is rotational, translational or unknown in the boundary file. Now I have to figure out to use it the right way!

Jens
Hello Jens,
Did you figure out how to tell cyclic to rotate or translate? Could you please share the method?
Thank you so much!
aqua is offline   Reply With Quote

Old   April 21, 2013, 10:26
Default
  #31
New Member
 
Esmaeil Ahmadiyan
Join Date: Sep 2011
Posts: 26
Rep Power: 7
e.ahmadiyan is on a distinguished road
Hi dear friends
I am using of openFoam2.0.1
My geometry and createPatchDict file are as below:

createPatchDict.zip

When running the case this error reported:

This probably means your geometry is not consistent with the specified separation and might lead to problems.
Continuing with specified separation vector (0.0002 0 0.0002)
patchutlet1 neighbour:inlet1
--> FOAM Warning :
From function cyclicPolyPatch::calcTransforms()
in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 368



In your mind are there any incorrect setting in createPatchDict file proportional to geometry?
What are separationVector and rotationAxis and rotationCentre ? please demonstrate for me. In fact I do not know how I should adjust this parameters in createPatchDict file proportional to my geometry.

All my case is as below
cyclic.zip

Thanks
__________________
Best regards

e.ahmadiyan
e.ahmadiyan is offline   Reply With Quote

Old   May 3, 2013, 02:22
Default
  #32
New Member
 
Esmaeil Ahmadiyan
Join Date: Sep 2011
Posts: 26
Rep Power: 7
e.ahmadiyan is on a distinguished road
Hi dear foamers
I am working with OF2.0.1 .My geometry and createPatchDict file are as below:

createPatchDict.zip

When running the case this error reported:

This probably means your geometry is not consistent with the specified separation and might lead to problems.
Continuing with specified separation vector (0.0002 0 0.0002)
patchutlet1 neighbour:inlet1
--> FOAM Warning :
From function cyclicPolyPatch::calcTransforms()
in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 368



In your mind are there any incorrect setting in createPatchDict file proportional to geometry?
What are separationVector and rotationAxis and rotationCentre ? please demonstrate for me. In fact I do not know how I should adjust this parameters in createPatchDict file proportional to my geometry.

All my case is as below
cyclic.zip

Please help me
Thanks
__________________
Best regards

e.ahmadiyan
e.ahmadiyan is offline   Reply With Quote

Old   May 3, 2013, 02:25
Default
  #33
New Member
 
Esmaeil Ahmadiyan
Join Date: Sep 2011
Posts: 26
Rep Power: 7
e.ahmadiyan is on a distinguished road
Hi dear foamers


Please help me
Thanks
__________________
Best regards

e.ahmadiyan
e.ahmadiyan is offline   Reply With Quote

Old   May 3, 2013, 02:45
Default
  #34
New Member
 
Esmaeil Ahmadiyan
Join Date: Sep 2011
Posts: 26
Rep Power: 7
e.ahmadiyan is on a distinguished road
11111111111111111111111111111111111111111111111111 1111111111111
__________________
Best regards

e.ahmadiyan
e.ahmadiyan is offline   Reply With Quote

Reply

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
MRFSimpleFoam amp cyclic patches david OpenFOAM Running, Solving & CFD 36 October 21, 2008 21:55
How to define these patches zhoubinwx Open Source Meshers: Gmsh, Netgen, CGNS, ... 6 September 15, 2008 08:52
Application of cyclic patches klaus OpenFOAM Meshing Format & General Technical 9 January 18, 2008 05:28
Cyclic boundary conditions for non conformed Patches turnow OpenFOAM Running, Solving & CFD 1 October 19, 2007 01:17
rediation patches novice CD-adapco 2 February 25, 2004 04:53


All times are GMT -4. The time now is 16:42.