Gmsh problem with cyclic bcbs
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? |
i have same problem. i used fluent3DMeshToFoam.
|
All times are GMT -4. The time now is 16:17. |