CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Rotational periodic Boundary condition (http://www.cfd-online.com/Forums/openfoam/76458-rotational-periodic-boundary-condition.html)

achinta May 26, 2010 06:09

Rotational periodic Boundary condition
 
hi,
I need to use rotational periodic boundary condition for my model. I used 'cyclic' boundary condition for 2 planes and OF gave the following error.
----------------------------------
Create time


Create mesh for time = 0






face 0 area does not match neighbour 1021 by 123.634% -- possible face ordering problem.
patch:SYM1 my area:3.13674e-06 neighbour area:1.32932e-05 matching tolerance:0.001
Mesh face:2900402 vertices:3((0.00766235 0.0434244 0) (0.00210264 0.0462494 0) (0.00291487 0.0469651 0))
Neighbour face:2901423 vertices:3((0.00390903 0.00637109 0) (0.00794004 0 0) (0.00376707 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 180.


FOAM exiting
---------------------------------------------------



How can I improve my mesh to use cyclic boundary condition?

Are there any other methods to define rotational periodic boundary condition?

Thank You.

Regards,
Achinta

achinta May 26, 2010 06:16

Hi,
I tried to use 'wedge' boundary condition on 2 planes to define rotational periodic condition for my model. OF gave the following error.
------------------------
Create time


Create mesh for time = 0


#0 Foam::error::printStack(Foam::Ostream&) in "/usr/logau/expsm/OpenFOAM/OpenFO
AM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/usr/logau/expsm/OpenFOAM/OpenFOAM-1.6/
lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::wedgePolyPatch::initTransforms() in "/usr/logau/expsm/OpenFOAM/OpenFOA
M-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 Foam::wedgePolyPatch::wedgePolyPatch(Foam::word const&, Foam::dictionary con st&, int, Foam::polyBoundaryMesh const&) in "/usr/logau/expsm/OpenFOAM/OpenFOAM- 1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#5 Foam::polyPatch::adddictionaryConstructorToTable<F oam::wedgePolyPatch>::New( Foam::word const&, Foam::dictionary const&, int, Foam::polyBoundaryMesh const&) in "/usr/logau/expsm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#6 Foam::polyPatch::New(Foam::word const&, Foam::dictionary const&, int, Foam:: polyBoundaryMesh const&) in "/usr/logau/expsm/OpenFOAM/OpenFOAM-1.6/lib/linux64G ccDPOpt/libOpenFOAM.so"
#7 Foam::polyBoundaryMesh::polyBoundaryMesh(Foam::IOo bject const&, Foam::polyMe sh const&) in "/usr/logau/expsm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpe nFOAM.so"
#8 Foam::polyMesh::polyMesh(Foam::IOobject const&) in "/usr/logau/expsm/OpenFOA M/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#9 Foam::fvMesh::fvMesh(Foam::IOobject const&) in "/usr/logau/expsm/OpenFOAM/Op enFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so"
#10 main in "/usr/logau/expsm/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64Gcc DPOpt/simpleFoam"
#11 __libc_start_main in "/lib64/libc.so.6"
#12 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start .S:116
Floating decimal point exception
---------------------------------------------



What could be the reason? I am new to OF and I am not able to understand what the error says.



Thank You



Regards,
Achinta

maddalena May 27, 2010 04:04

suggestions
 
Hello,
so your problems have been not fixed yet...
Quote:

Originally Posted by achinta (Post 260314)
I tried to use 'wedge' boundary condition on 2 planes to define rotational periodic condition for my model

This will not work for sure if your model is 3D! wedge requires that the mesh is only one cell thick!

Quote:

Originally Posted by achinta (Post 260314)
hi,
I need to use rotational periodic boundary condition for my model. I used 'cyclic' boundary condition for 2 planes and OF gave the following error.
----------------------------------
Create time


Create mesh for time = 0


face 0 area does not match neighbour 1021 by 123.634% -- possible face ordering problem.
patch:SYM1 my area:3.13674e-06 neighbour area:1.32932e-05 matching tolerance:0.001
Mesh face:2900402 vertices:3((0.00766235 0.0434244 0) (0.00210264 0.0462494 0) (0.00291487 0.0469651 0))
Neighbour face:2901423 vertices:3((0.00390903 0.00637109 0) (0.00794004 0 0) (0.00376707 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 180.

Please check that:
  1. your plane mesh must be the equal; use a constraint on your meshing software to keep the plane identical.
  2. check that your cyclic domain is formed only by two coupled planes, and not by multiple planes. Every cyclic plane must be connected to its counterpart! i.e. if you have 2 couples of 2 cyclic planes, define them in polyMesh/boundary separately.
If you have already done that, try to increase the matching tolerance in createPatchDict. Try to increase it little by little, as soon as the two planes match.
Hope this help,
cheers,

maddalena


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