Personal tools
You are here: Home Forum Use cycle boundary condition for OpenFoam

cycle boundary condition for OpenFoam

Up to Use

cycle boundary condition for OpenFoam

Posted by Makoto Yamamoto at August 03. 2012

Hello, 

I have a simple box for a flow simulation with inlet and outlet normal to x-axis. normal to z-axis there are walls. And the sidewalls normal to the y-axis are supposed to have later in OpenFoam cycle boundary conditions. 

OpenFoam needs coincident nodes on both sidewalls for the cycle boundary condition. How can I achieve that in Salome?

Thank you very much!!

/Yato

 

Re: cycle boundary condition for OpenFoam

Posted by Makoto Yamamoto at August 06. 2012

Sorry. I mean of course cyclic boundary conditions. 

I still would be glad for some help on this issue....

Re: cycle boundary condition for OpenFoam

Posted by Chrys Lides at August 06. 2012

Hello,

If you need for each node of a face to be in relation to a node on the opposite face, you have to use "Projection 2D" algorithm: http://docs.salome-platform.org/salome_6_5_0/gui/SMESH/projection_algos_page.html

Define a submesh on the opposite face, with the base face as source face. You can leave the other parameters blank.

If the nodes of the faces wires are not in vis-à-vis, you should use "Projection 1D-2D" instead.

Chrys

Re: cycle boundary condition for OpenFoam

Posted by Makoto Yamamoto at August 07. 2012
Hello Chrys,
 
thanks a lot for your help. 
 
I have tried to follow your instruction, but I got this message in OpenFoam:
 
--> FOAM FATAL ERROR: 
For patch side1 there are 244 face centres, for the neighbour patch side2 there are 236
 
What I have done:
 
1) create a block in Salome and create groups of faces. walls, inlet, outlet, side1 and side2. side1 and side2 are the opposite faces which should be physically related to each other. 
 
2) create a mesh (Mesh_1). Geometry is the block and I use Netgen 1D-2D-3D algorithm. 
 
3) create a submesh: geometry: side2; mesh: Mesh_1; algorithm: Projection 2D. in hypothesis construction I picked for face: side1 and leave the other blank. 
 
4) compute the mesh and assign the patches for inlet....
 
5) export to unv and import to openfoam...
     
 
Can you tell me what I am doing wrong?
 
Thanks a lot in advance. 
 
Regards,
yato
 

Re: cycle boundary condition for OpenFoam

Posted by Chrys Lides at August 07. 2012

Projection algorithm seems not to work with NETGEN 1D-2D-3D.

 

Please define the algorithm for each dimension separately:

* 3D: NETGEN

* 2D: NETGEN

* 1D: wire discretisation, local length or number of segments.

 

Then set the projection 2D submesh as you did.

 

See attached script for a working example.

 

You can check the number of faces by calling Mesh information on the two groups of faces.

 

Chrys

Attachments

Re: cycle boundary condition for OpenFoam

Posted by Makoto Yamamoto at August 09. 2012

Hello Chrys,

 

thank you very much for your fast response. 

 

After using cyclic boundary conditions on the box. I tried to set them also on my actual model. There I have the problem, that the faces are displaced and I get following error message in openFoam while running checkMesh

--> FOAM FATAL ERROR: 

face 0 area does not match neighbour by 23.7135% -- possible face ordering problem.

patch:side1 my area:4.45648e-06 neighbour area:3.51171e-06 matching tolerance:0.0001

Mesh face:360585 fc:(0.0885241 0.114046 0.001)

Neighbour fc:(0.041739 0.0930121 0.00372698)

 

I have no idea how to solve this issue. Do you have any suggestions?

 

yato

Re: cycle boundary condition for OpenFoam

Posted by Chrys Lides at August 10. 2012

You can try to use the OpenFOAM utility createPatch to order the faces as expected by OpenFOAM. A user had a similar issue on http://www.openfoam.org/mantisbt/view.php?id=154

Re: cycle boundary condition for OpenFoam

Posted by Makoto Yamamoto at August 22. 2012

Ok I managed it by using cyclicAMI BCs. They are availabe since OpenFoam v2.1.0 and you don't need exact matching points anymore.

If anyone has a similar problem, this might be helpful: 

http://www.cfd-online.com/Forums/openfoam-pre-processing/105531-setting-side-wedge-boundary-condition.html#post378087

Thank you for your help. 

Powered by Ploneboard
Document Actions