# Rotational periodic Boundary condition

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

 May 26, 2010, 06:09 Rotational periodic Boundary condition #1 Member   achinta Join Date: May 2010 Location: Sydney Posts: 66 Rep Power: 9 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

 May 26, 2010, 06:16 #2 Member   achinta Join Date: May 2010 Location: Sydney Posts: 66 Rep Power: 9 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:rintStack(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:olyBoundaryMesh const&) in "/usr/logau/expsm/OpenFOAM/OpenFOAM- 1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #5 Foam:olyPatch::adddictionaryConstructorToTable::New( Foam::word const&, Foam::dictionary const&, int, Foam:olyBoundaryMesh const&) in "/usr/logau/expsm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #6 Foam:olyPatch::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:olyBoundaryMesh:olyBoundaryMesh(Foam::IOo bject const&, Foam:olyMe sh const&) in "/usr/logau/expsm/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpe nFOAM.so" #8 Foam:olyMesh:olyMesh(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

May 27, 2010, 04:04
suggestions
#3
Senior Member

Join Date: Mar 2009
Posts: 436
Rep Power: 15
Hello,
so your problems have been not fixed yet...
Quote:
 Originally Posted by achinta 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 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.
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,