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

Cyclic GGi

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 13, 2011, 05:57
Default Cyclic GGi
  #1
Member
 
Kurne
Join Date: Aug 2010
Location: Pune, INDIA
Posts: 88
Rep Power: 6
kurne is on a distinguished road
Dear All

My case is of Turbo machinery and i am using OpenFOAM 1.5-Dev.I have successfully transfer the mesh from the Ansys ICEM to Open Foam and now i want to use cyclicggi to two patches.Can any one tell me what necessary changes to make and where.I am in great need of it.

Please Help.

Thanks In Advance.
__________________
Simulation Is Determination Of Imagination Towards Approximation ®


Best Regards

Mubeen K Kurne
kurne is offline   Reply With Quote

Old   April 15, 2011, 06:06
Default
  #2
Member
 
Kurne
Join Date: Aug 2010
Location: Pune, INDIA
Posts: 88
Rep Power: 6
kurne is on a distinguished road
Dear All

Please tell me. I am in need of it.
__________________
Simulation Is Determination Of Imagination Towards Approximation ®


Best Regards

Mubeen K Kurne
kurne is offline   Reply With Quote

Old   April 15, 2011, 07:15
Default
  #3
New Member
 
amgode's Avatar
 
Amol
Join Date: Jul 2009
Location: Pune, INDIA
Posts: 23
Blog Entries: 2
Rep Power: 8
amgode is on a distinguished road
Check here for a test case with regards to GGI (Thanks to Olivier Petit, Maryse Page, Håkan Nilsson and Martin Beaudoin)

for cyclicGGI the procedure for faceSet and faceZone creation is same as GGI
The definition in constant/boundary file is something like:


cylic_left
{
type cyclicGgi;
nFaces 100;
startFace 1211;
shadowPatch cyclic_right;
zone cyclic_left_zone;
bridgeOverlap off;
rotationAxis (0 1 0);
rotationAngle 90;
separationOffset (0 0 0);
}

cylic_right
{
type cyclicGgi;
nFaces 100;
startFace 1311;
shadowPatch cyclic_left;
zone cyclic_right_zone;
bridgeOverlap off;
rotationAxis (0 1 0);
rotationAngle 90;
separationOffset (0 0 0);
}

Where the entries are self explanatory..........

Also,

The boundary type in the 0/ directory for all fields is cyclicGGI.


Hope that helps !!


Cheers !

Amol
__________________
Amol Gode
amgode is offline   Reply With Quote

Old   April 18, 2011, 06:19
Default
  #4
Member
 
Kurne
Join Date: Aug 2010
Location: Pune, INDIA
Posts: 88
Rep Power: 6
kurne is on a distinguished road
Dear Amol Gode

Thanks a lot for guiding me.I have done all the necessary changes in constant/polymesh/boundary file.My case is of turbo machinery and i am using simpleSRFFoam solver.

I have applied the correct boundary conditions also but still i am getting the error which is as follows.

CCreate time

Create mesh for time = 0

Reading field p

Reading field Urel

Reading/calculating face flux field phi

Face zone name PERI1_zone not found. Please check your GGI interface definition.

From function label ggiPolyPatch::zoneIndex() const
in file meshes/polyMesh/polyPatches/constraint/ggi/ggiPolyPatch.C at line 404.

FOAM aborting

Aborted


The error is showing that PERI1 zone not found but the PERI1 patch is there.

Will anybody help me.
Thanks In Advance.
__________________
Simulation Is Determination Of Imagination Towards Approximation ®


Best Regards

Mubeen K Kurne

Last edited by kurne; April 18, 2011 at 06:35.
kurne is offline   Reply With Quote

Old   April 18, 2011, 09:12
Default
  #5
New Member
 
amgode's Avatar
 
Amol
Join Date: Jul 2009
Location: Pune, INDIA
Posts: 23
Blog Entries: 2
Rep Power: 8
amgode is on a distinguished road
Quote:
Originally Posted by kurne View Post
Face zone name PERI1_zone not found. Please check your GGI interface definition.
There does not seem to be a face zone with name 'PERI1_zone".
Can be checked in constant/polyMesh/faceZones

Also the face zone name 'PERI1_zone' seems to present in the definition of GGI / cyclicGGI interfaces, though actually not available as a face zone.



Amol
__________________
Amol Gode
amgode is offline   Reply With Quote

Old   April 19, 2011, 00:52
Default
  #6
Member
 
Kurne
Join Date: Aug 2010
Location: Pune, INDIA
Posts: 88
Rep Power: 6
kurne is on a distinguished road
Dear Amol

I have check the constant/polyMesh/faceZones directory and in this directory only numbers are present means there is no names of zones as i have transfer the mesh from the ICEM to OpenFOAM.In boundary file i have edited the patches to be made cyclicGgi and which is given below

PERI1
{
type cyclicGgi;
nFaces 2517;
startFace 758483;
shadowPatch periodic_sh;
zone PERI1;
bridgeOverlap off;
rotationAxis (0 0 1);
rotationAngle -72;
separationOffset (0 0 0);
}
periodic_sh
{
type cyclicGgi;
nFaces 2517;
startFace 761000;
shadowPatch PERI1;
zone periodic_sh;
bridgeOverlap off;
rotationAxis (0 0 1);
rotationAngle 72;
separationOffset (0 0 0);
}

I have correctly edited the cyclicGgi interfaces and you can check once again.What is the problem i am unable to recognize.Please help to sort out this problem.Can anybody help me.

Thanks in Advance.
__________________
Simulation Is Determination Of Imagination Towards Approximation ®


Best Regards

Mubeen K Kurne

Last edited by kurne; April 19, 2011 at 07:33.
kurne is offline   Reply With Quote

Old   April 19, 2011, 06:11
Default
  #7
New Member
 
amgode's Avatar
 
Amol
Join Date: Jul 2009
Location: Pune, INDIA
Posts: 23
Blog Entries: 2
Rep Power: 8
amgode is on a distinguished road
Quote:
Originally Posted by kurne View Post

PERI1
{
type cyclicGgi;
nFaces 2517;
startFace 758483;
shadowPatch periodic_sh;
zone PERI1_ZONE;
bridgeOverlap off;
rotationAxis (0 0 1);
rotationAngle -72;
separationOffset (0 0 0);
}
periodic_sh
{
type cyclicGgi;
nFaces 2517;
startFace 761000;
shadowPatch PERI1;
zone periodic_sh_ZONE;
bridgeOverlap off;
rotationAxis (0 0 1);
rotationAngle 72;
separationOffset (0 0 0);
}

The faceZones file should look something like:



FoamFile
{
version 2.0;
format ascii;
class regIOobject;
location "constant/polyMesh";
object faceZones;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

2
(
periodic_sh_ZONE
{
type faceZone;
faceLabels List<label>
36
(
...
//faces in periodic_sh_ZONE
....
)
;
flipMap List<bool> 36{0};
}
PERI1_ZONE
{
type faceZone;
faceLabels List<label>
36
(
....
//faces in PERI1_ZONE
.....
)
;
flipMap List<bool> 36{0};
}
)




If not then please check the steps taken for creation of faceSets and faceZones.

Check the tutorial on OpenFoamWiki referred to above.


Amol
__________________
Amol Gode
amgode is offline   Reply With Quote

Old   April 19, 2011, 07:32
Default
  #8
Member
 
Kurne
Join Date: Aug 2010
Location: Pune, INDIA
Posts: 88
Rep Power: 6
kurne is on a distinguished road
Dear Amol

Thank you very much brother.I created the faceZone from the tutorial file.But i am facing the same problem as you face and mention here GGI implementation in MRFSimpleFoam

When i run the solver i get the error

Create time
Create mesh for time = 0
Reading field p
Reading field Urel
Reading/calculating face flux field phi

Problem with patch-to zone addressing: some patch faces not found in interpolation zone
From function void ggiPolyPatch::calcZoneAddressing() const
in file meshes/polyMesh/polyPatches/constraint/ggi/ggiPolyPatch.C at line 77.

FOAM aborting
Aborted

Hope this would be the final error and as you know the solution of it also please help brother.
Thank You very Much Once Again.
__________________
Simulation Is Determination Of Imagination Towards Approximation ®


Best Regards

Mubeen K Kurne
kurne is offline   Reply With Quote

Old   April 20, 2011, 03:18
Default
  #9
New Member
 
amgode's Avatar
 
Amol
Join Date: Jul 2009
Location: Pune, INDIA
Posts: 23
Blog Entries: 2
Rep Power: 8
amgode is on a distinguished road
You could find the reply in the same thread you quoted
__________________
Amol Gode
amgode is offline   Reply With Quote

Old   April 20, 2011, 03:24
Default
  #10
Member
 
Kurne
Join Date: Aug 2010
Location: Pune, INDIA
Posts: 88
Rep Power: 6
kurne is on a distinguished road
Dear Amol

I have check all the things mention in that thread but still getting the error.Hope you help me now.
I get the faceZones file as given below

)
;
flipMap List<bool> 731419{0};
}
periodic_sh_zone
{
type faceZone;
faceLabels 0();
flipMap 0();
}
PERI1_zone
{
type faceZone;
faceLabels 0();
flipMap 0();
}
)

And my boundary file is as follows

PERI1
{
type cyclicGgi;
nFaces 2517;
startFace 758483;
shadowPatch periodic_sh;
zone PERI1_zone;
bridgeOverlap off;
rotationAxis (0 0 1);
rotationAngle -72;
separationOffset (0 0 0);
}
periodic_sh
{
type cyclicGgi;
nFaces 2517;
startFace 761000;
shadowPatch PERI1;
zone periodic_sh_zone;
bridgeOverlap off;
rotationAxis (0 0 1);
rotationAngle 72;
separationOffset (0 0 0);
}


Thanks In Advance.
__________________
Simulation Is Determination Of Imagination Towards Approximation ®


Best Regards

Mubeen K Kurne

Last edited by kurne; April 20, 2011 at 07:17.
kurne is offline   Reply With Quote

Old   April 20, 2011, 05:12
Default
  #11
New Member
 
Wei Zhao
Join Date: Mar 2010
Posts: 28
Rep Power: 7
zhaowei is on a distinguished road
c this,

http://www.tfd.chalmers.se/~hani/kur..._kurs_2010.tgz

i think u can get what u want
zhaowei is offline   Reply With Quote

Old   April 20, 2011, 05:35
Default
  #12
Member
 
Kurne
Join Date: Aug 2010
Location: Pune, INDIA
Posts: 88
Rep Power: 6
kurne is on a distinguished road
Dear Wei Zhao

Thanks a lot for it and it is new case file but the problem solution is not in it.Anyway thanks a lot for new case files.
__________________
Simulation Is Determination Of Imagination Towards Approximation ®


Best Regards

Mubeen K Kurne
kurne is offline   Reply With Quote

Old   April 21, 2011, 05:19
Default
  #13
Member
 
Kurne
Join Date: Aug 2010
Location: Pune, INDIA
Posts: 88
Rep Power: 6
kurne is on a distinguished road
Dear Amol

Thank you very much brother my case got work.

Thanks A Lot Once Again.
__________________
Simulation Is Determination Of Imagination Towards Approximation ®


Best Regards

Mubeen K Kurne

Last edited by kurne; April 21, 2011 at 06:15.
kurne is offline   Reply With Quote

Old   October 23, 2013, 02:07
Unhappy turboPassageRotating2d
  #14
New Member
 
Minh-VietNam
Join Date: Oct 2013
Posts: 7
Rep Power: 3
mingbn is on a distinguished road
Hi. Everybody.
I'm a new openFoam user. so i can not launch this case. I don't know how to fix the error. Can you help please.?
This is case (case of an master thesis which downloaded from web http://www.tfd.chalmers.se/~hani/kurser/OS_CFD/)
thanks a lot!
======
This is error:
--> FOAM FATAL IO ERROR:

Cannot find 'value' entry on patch rotor_cyclic_upper of field p in file "/home/minh/Desktop/turboPassageRotating2D/0/p"
which is required to set the values of the generic patch field.
(Actual type cyclicGgi)

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


file: /home/minh/Desktop/turboPassageRotating2D/0/p.boundaryField.rotor_cyclic_upper from line 43 to line 43.

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

FOAM exiting

+++
Please show to me, thankS!
Attached Files
File Type: gz turboPassageRotating2D.tar.gz (57.2 KB, 8 views)

Last edited by mingbn; November 1, 2013 at 02:40.
mingbn 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
Ggi FabOr OpenFOAM 17 May 9, 2013 10:19
Cyclic vs ggi vs directMapped Patches jens_klostermann OpenFOAM Running, Solving & CFD 33 May 3, 2013 02:45
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
Pressure instability with rhoSimpleFoam daniel_mills OpenFOAM Running, Solving & CFD 44 February 17, 2011 18:08


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