|
[Sponsors] |
April 11, 2008, 17:36 |
Gmsh problem with cyclic bcbs
|
#1 |
Member
Scott Ripplinger
Join Date: Mar 2009
Posts: 30
Rep Power: 17 |
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:stream&) 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 zone950 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 zone950 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? |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
LES in OF 3.0.1, no Output for nuSgs | Mirage | OpenFOAM Programming & Development | 1 | October 4, 2016 19:00 |
[mesh manipulation] Problem with using createPatchDict to define cyclic boundaries | kaszt | OpenFOAM Meshing & Mesh Conversion | 0 | April 1, 2016 22:18 |
Possible createPatch/createBaffles bug? | simpomann | OpenFOAM Bugs | 2 | July 15, 2014 08:07 |
[Gmsh] gmsh or gmshToFoam problem? | Madeleine P. Vincent | OpenFOAM Meshing & Mesh Conversion | 0 | June 19, 2012 09:52 |
[Gmsh] Import problem | ARC | OpenFOAM Meshing & Mesh Conversion | 0 | February 27, 2010 11:56 |