CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   Open Source Meshers: Gmsh, Netgen, CGNS, ... (http://www.cfd-online.com/Forums/openfoam-meshing-open/)
-   -   Gmsh problem with cyclic bcbs (http://www.cfd-online.com/Forums/openfoam-meshing-open/61882-gmsh-problem-cyclic-bcbs.html)

sripplinger April 11, 2008 16:36

I'm trying to run a simple 3D
 
I'm trying to run a simple 3D rectangular domain. It has an inlet and outlet, walls at the top and on bottom, and periodic boundaries on left and right. The gmshToFoam utility seems to have worked fine. I manually edited the bc's in /constant/polymesh/boundary, as well as in /0/U and /0/p. When I try to run icoFoam I get the following message:

Reading field p



--> FOAM FATAL ERROR : face 0 and 475 areas do not match by 176.205% -- possible face ordering problem#0 Foam::error::printStack(Foam:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) in "/home/Scott/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/Scott/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam::cyclicFvPatch::makeWeights(Foam::Field<doubl e>&) const in "/home/Scott/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#3 Foam::surfaceInterpolation::makeWeights() const in "/home/Scott/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#4 Foam::surfaceInterpolation::weights() const in "/home/Scott/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::fvPatch::weights() const in "/home/Scott/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#6 Foam::coupledFvPatchField<double>::evaluate() in "/home/Scott/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#7 Foam::cyclicFvPatchField<double>::cyclicFvPatchFie ld(Foam::fvPatch const&, Foam::DimensionedField<double,> const&, Foam::dictionary const&) in "/home/Scott/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#8 Foam::fvPatchField<double>::adddictionaryConstruct orToTable<foam::cyclicfvpatchf ield<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<double,> const&, Foam::dictionary const&) in "/home/Scott/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#9 Foam::fvPatchField<double>::New(Foam::fvPatch const&, Foam::DimensionedField<double,> const&, Foam::dictionary const&) in "/home/Scott/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/icoFoam"
#10 Foam::GeometricField<double,>::GeometricBoundaryFi eld::GeometricBoundaryField(Fo am::fvBoundaryMesh const&, Foam::DimensionedField<double,> const&, Foam::dictionary const&) in "/home/Scott/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/icoFoam"
#11 Foam::GeometricField<double,>::readField(Foam::Ist ream&) in "/home/Scott/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/icoFoam"
#12 Foam::GeometricField<double,>::GeometricField(Foam ::IOobject const&, Foam::fvMesh const&) in "/home/Scott/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/icoFoam"
#13 main in "/home/Scott/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/icoFoam"
#14 __libc_start_main in "/lib64/libc.so.6"
#15 Foam::regIOobject::readIfModified() in "/home/Scott/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/icoFoam"


From function cyclicFvPatch::makeWeights(scalarField& w) const
in file fvMesh/fvPatches/constraint/cyclic/cyclicFvPatch.C at line 58.

FOAM aborting

Aborted


I also tried running couplePatches and got the following:

Mesh has coupled patches ...

Doing dummy mesh morph to correct face ordering ...
cyclicPolyPatch::order : Number of faces per zonehttp://www.cfd-online.com/OpenFOAM_D...lipart/sad.gif950 0)
cyclicPolyPatch::order : Writing half0 faces to OBJ file "rightSide_half0_faces.obj"
cyclicPolyPatch::order : Writing half1 faces to OBJ file "rightSide_half1_faces.obj"
cyclicPolyPatch::order : Writing half0 face centres to OBJ file "rightSide_half0.obj"
cyclicPolyPatch::order : Writing half1 face centres to OBJ file "rightSide_half1.obj"
--> FOAM Serious Error :
From function cyclicPolyPatch::order(const primitivePatch&, labelList&, labelList&) const
in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 596
patch:rightSide : Patch rightSide gets decomposed in two zones ofinequal size: 950 and 0
This means that the patch is either not two separate regions or one region where the angle between the different regions is not sufficiently sharp.
Please use topological matching or adapt the featureCos() setting
Continuing with incorrect face ordering from now on!
cyclicPolyPatch::order : Number of faces per zonehttp://www.cfd-online.com/OpenFOAM_D...lipart/sad.gif950 0)
cyclicPolyPatch::order : Writing half0 faces to OBJ file "leftSide_half0_faces.obj"
cyclicPolyPatch::order : Writing half1 faces to OBJ file "leftSide_half1_faces.obj"
cyclicPolyPatch::order : Writing half0 face centres to OBJ file "leftSide_half0.obj"
cyclicPolyPatch::order : Writing half1 face centres to OBJ file "leftSide_half1.obj"
--> FOAM Serious Error :
From function cyclicPolyPatch::order(const primitivePatch&, labelList&, labelList&) const
in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 596
patch:leftSide : Patch leftSide gets decomposed in two zones ofinequal size: 950 and 0
This means that the patch is either not two separate regions or one region where the angle between the different regions is not sufficiently sharp.
Please use topological matching or adapt the featureCos() setting
Continuing with incorrect face ordering from now on!
Mesh ordering ok. Nothing changed.
End


Has anybody been able to successfully use cyclic bc's on a gmsh converted mesh?

maysmech September 28, 2010 09:29

i have same problem. i used fluent3DMeshToFoam.


All times are GMT -4. The time now is 23:42.