|September 19, 2008, 11:08||
Hi, I am trying to simulate
Join Date: Mar 2009
Location: Stuttgart, Germany
Posts: 77Rep Power: 8
I am trying to simulate the flow in a planar channel with streamwise periodic hills on the lower wall. The flow domain has planar, cyclic boundary conditions in streamwise and spanwise directions.
Using cyclic boundary conditions in streamwise direction works just fine, however, a face ordering problem exists when I use cyclic boundaries in spanwise direction.
face 0 area does not match neighbour 722 by 37.8947% -- possible face ordering problem.
patch:CYCLICSYMM my area:0.00280653 neighbour area:0.00411868 matching tolerance:0.001
Mesh face:41515 vertices:4((0 0.0526316 0) (0.0526316 0.0540166 0) (0.0526316 0 0) (0 0 0))
Neighbour face:42237 vertices:4((1 0 0) (1 0.0789474 0) (1.05263 0.0775623 0) (1.05263 0 0))
Rerun with cyclic debug flag set for more information.
From function cyclicPolyPatch::calcTransforms()
in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 150.
For the same case with the hills removed (simple planar channel flow with cyclic BC in streamwise and spanwise direction) this problem does not exist. The only difference between the two cases is that without the hills in place, the faces in spanwise directions are rectangular and with the hills in place, the faces are non-rectangular.
I also tried to export the grid with two separate boundary faces in spanwise direction and then tried to combine the faces into a cyclic patch using createPatch. Unfortunatley, this always resulted in a crash of createPatch.
Does anybody know why the face ordering problem suddenly arises when I am using non-rectangular patches as cyclic BC and more importantly, how can I solve this face ordering problem?
|September 26, 2008, 11:13||
I have the same problem, its a
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 550Rep Power: 18
I have the same problem, its always when I am trying to use a 3D mesh with a cyclic bc. Try renumberMesh and see what happens. Although I was not able to have any better luck when I did this. I did also try the createPatch utility and it hasn't worked in OF 1.5 for me yet, but did in 1.4.1-dev.
|August 17, 2011, 15:07||
Join Date: Feb 2011
Posts: 15Rep Power: 6
Check with the vertices values.. there is a difference in betwwen these points.
"4((0 0.0526316 0) (0.0526316 0.0540166 0) (0.0526316 0 0) (0 0 0))
4((1 0 0) (1 0.0789474 0) (1.05263 0.0775623 0) (1.05263 0 0)) ".
Hope it helps.
|Thread||Thread Starter||Forum||Replies||Last Post|
|Error in face ordering||anke||OpenFOAM Mesh Utilities||0||May 6, 2008 04:02|
|Cyclic Boundary conditions||johndeas||OpenFOAM Running, Solving & CFD||1||January 4, 2008 21:34|
|Cyclic boundary conditions||hadi||OpenFOAM Running, Solving & CFD||0||June 29, 2007 10:46|
|periodic and cyclic boundary conditions||Hadi||Main CFD Forum||2||June 29, 2007 07:19|
|cyclic boundary conditions..||Zerrin||CD-adapco||0||June 11, 2005 13:49|