Agit |
November 23, 2018 08:45 |
1 Attachment(s)
Hi guys,
I know, that this thread is pretty old, but I'm having the same issues.
I am using OpenFOAM 6 and snappyHexMesh.
Here is how I get the same matching Error:
1. blockMesh
2. decomposePar
3. runParallel snappyHexMesh -overwrite
SnappyHexMesh is running until it is snapping and displacing the mesh to fit to my surfaceFeatureExtract-Object (or in some cases until layerAddition).
Here is the last bit of my log.snappyHexMesh file:
Code:
Morph iteration 3
-----------------
Calculating patchDisplacement as distance to nearest surface point ...
Wanted displacement : average:0.000594714 min:5.62277e-07 max:0.00214618
Calculated surface displacement in = 0.11 s
Detecting near surfaces ...
Overriding nearest with intersection of close gaps at 0 out of 27372 points.
Overriding displacement on features :
implicit features : false
explicit features : true
multi-patch features : false
Morph iteration 3
-----------------
Calculating patchDisplacement as distance to nearest surface point ...
Wanted displacement : average:0.000594714 min:5.62277e-07 max:0.00214618
Calculated surface displacement in = 0.11 s
Detecting near surfaces ...
Overriding nearest with intersection of close gaps at 0 out of 27372 points.
Overriding displacement on features :
implicit features : false
explicit features : true
multi-patch features : false
Detected 0 baffle edges out of 54250 edges.
--> FOAM Warning :
From function Foam::treeBoundBox::treeBoundBox(const Foam::UList<Foam::Vector<double> >&)
in file meshes/treeBoundBox/treeBoundBox.C at line 136
cannot find bounding box for zero-sized pointField, returning zero
Initially selected 0 points out of 27372 for reverse attraction.
Selected 0 points out of 27372 for reverse attraction.
Stringing feature edges : changed 0 points
Attraction:
linear : max:(0.00183956 -0.00110546 -7.51662e-06) avg:(4.50171e-05 -2.00879e-05 -2.80543e-09)
feature : max:(0 0 0) avg:(0 0 0)
Feature analysis : total master points:27048 attraction to :
feature point : 0
feature edge : 0
nearest surface : 0
rest : 27048
Smoothing displacement ...
Iteration 0
Iteration 10
Iteration 20
Displacement smoothed in = 6.6 s
Moving mesh ...
Iteration 0
Moving mesh using displacement scaling : min:1 max:1
Correcting 2-D mesh motion--> FOAM Warning :
From function void Foam::motionSmootherAlgo::modifyMotionPoints(Foam::pointField&) const
in file motionSmoother/motionSmootherAlgo.C at line 657
2D mesh-motion probably not correct in parallel
...done
[1] processorPolyPatch::calcGeometry : Writing my 14323 faces to OBJ file "/media/agit/Volume/OpenFOAM/_baseCaseTest01/1_cases/0.507/5E-02/processor1/procBoundary1to2_faces.obj"
[2] processorPolyPatch::calcGeometry : Writing my 14323 faces to OBJ file "/media/agit/Volume/OpenFOAM/_baseCaseTest01/1_cases/0.507/5E-02/processor2/procBoundary2to1_faces.obj"
[2] processorPolyPatch::calcGeometry : Dumping cell centre lines between corresponding face centres to OBJ file"/media/agit/Volume/OpenFOAM/_baseCaseTest01/1_cases/0.507/5E-02/processor2/procBoundary2to1_faceCentresConnections.obj"
[1] processorPolyPatch::calcGeometry : Dumping cell centre lines between corresponding face centres to OBJ file"/media/agit/Volume/OpenFOAM/_baseCaseTest01/1_cases/0.507/5E-02/processor1/procBoundary1to2_faceCentresConnections.obj"
[2]
[2]
[2] --> FOAM FATAL ERROR:
[2] face 5738 area does not match neighbour by 0.0267405% -- possible face ordering problem.
patch:procBoundary2to1 my area:3.99941e-06 neighbour area:3.99835e-06 matching tolerance:1.03479e-09
Mesh face:1177078 vertices:4((-0.0929725 0.00598575 0.04625) (-0.0993708 0.00603835 0.04625) (-0.0993676 0.00603726 0.045625) (-0.0929683 0.00598432 0.045625))
If you are certain your matching is correct you can increase the 'matchTolerance' setting in the patch dictionary in the boundary file.
Rerun with processor debug flag set for more information.
[2]
[2] From function virtual void Foam::processorPolyPatch::calcGeometry(Foam::PstreamBuffers&)
[2] in file meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C at line 284.
[2]
FOAM parallel run exiting
[2]
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 2 in communicator MPI COMMUNICATOR 3 SPLIT FROM 0
with errorcode 1.
NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
[1]
[1]
[1] --> FOAM FATAL ERROR:
[1] face 5738 area does not match neighbour by 0.0267405% -- possible face ordering problem.
patch:procBoundary1to2 my area:3.99835e-06 neighbour area:3.99941e-06 matching tolerance:1.03424e-09
Mesh face:1187038 vertices:4((-0.0929725 0.00598575 0.04625) (-0.0929683 0.00598432 0.045625) (-0.0993664 0.00603686 0.045625) (-0.0993686 0.00603759 0.04625))
If you are certain your matching is correct you can increase the 'matchTolerance' setting in the patch dictionary in the boundary file.
Rerun with processor debug flag set for more information.
[1]
[1] From function virtual void Foam::processorPolyPatch::calcGeometry(Foam::PstreamBuffers&)
[1] in file meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C at line 284.
[1]
FOAM parallel run exiting
[1]
Detected 0 baffle edges out of 54250 edges.
--> FOAM Warning :
From function Foam::treeBoundBox::treeBoundBox(const Foam::UList<Foam::Vector<double> >&)
in file meshes/treeBoundBox/treeBoundBox.C at line 136
cannot find bounding box for zero-sized pointField, returning zero
Initially selected 0 points out of 27372 for reverse attraction.
Selected 0 points out of 27372 for reverse attraction.
Stringing feature edges : changed 0 points
Attraction:
linear : max:(0.00183956 -0.00110546 -7.51662e-06) avg:(4.50171e-05 -2.00879e-05 -2.80543e-09)
feature : max:(0 0 0) avg:(0 0 0)
Feature analysis : total master points:27048 attraction to :
feature point : 0
feature edge : 0
nearest surface : 0
rest : 27048
Smoothing displacement ...
Iteration 0
Iteration 10
Iteration 20
Displacement smoothed in = 6.6 s
Moving mesh ...
Iteration 0
Moving mesh using displacement scaling : min:1 max:1
Correcting 2-D mesh motion--> FOAM Warning :
From function void Foam::motionSmootherAlgo::modifyMotionPoints(Foam::pointField&) const
in file motionSmoother/motionSmootherAlgo.C at line 657
2D mesh-motion probably not correct in parallel
...done
[1] processorPolyPatch::calcGeometry : Writing my 14323 faces to OBJ file "/media/agit/Volume/OpenFOAM/_baseCaseTest01/1_cases/0.507/5E-02/processor1/procBoundary1to2_faces.obj"
[2] processorPolyPatch::calcGeometry : Writing my 14323 faces to OBJ file "/media/agit/Volume/OpenFOAM/_baseCaseTest01/1_cases/0.507/5E-02/processor2/procBoundary2to1_faces.obj"
[2] processorPolyPatch::calcGeometry : Dumping cell centre lines between corresponding face centres to OBJ file"/media/agit/Volume/OpenFOAM/_baseCaseTest01/1_cases/0.507/5E-02/processor2/procBoundary2to1_faceCentresConnections.obj"
[1] processorPolyPatch::calcGeometry : Dumping cell centre lines between corresponding face centres to OBJ file"/media/agit/Volume/OpenFOAM/_baseCaseTest01/1_cases/0.507/5E-02/processor1/procBoundary1to2_faceCentresConnections.obj"
[2]
[2]
[2] --> FOAM FATAL ERROR:
[2] face 5738 area does not match neighbour by 0.0267405% -- possible face ordering problem.
patch:procBoundary2to1 my area:3.99941e-06 neighbour area:3.99835e-06 matching tolerance:1.03479e-09
Mesh face:1177078 vertices:4((-0.0929725 0.00598575 0.04625) (-0.0993708 0.00603835 0.04625) (-0.0993676 0.00603726 0.045625) (-0.0929683 0.00598432 0.045625))
If you are certain your matching is correct you can increase the 'matchTolerance' setting in the patch dictionary in the boundary file.
Rerun with processor debug flag set for more information.
[2]
[2] From function virtual void Foam::processorPolyPatch::calcGeometry(Foam::PstreamBuffers&)
[2] in file meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C at line 284.
[2]
FOAM parallel run exiting
[2]
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 2 in communicator MPI COMMUNICATOR 3 SPLIT FROM 0
with errorcode 1.
NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
[1]
[1]
[1] --> FOAM FATAL ERROR:
[1] face 5738 area does not match neighbour by 0.0267405% -- possible face ordering problem.
patch:procBoundary1to2 my area:3.99835e-06 neighbour area:3.99941e-06 matching tolerance:1.03424e-09
Mesh face:1187038 vertices:4((-0.0929725 0.00598575 0.04625) (-0.0929683 0.00598432 0.045625) (-0.0993664 0.00603686 0.045625) (-0.0993686 0.00603759 0.04625))
If you are certain your matching is correct you can increase the 'matchTolerance' setting in the patch dictionary in the boundary file.
Rerun with processor debug flag set for more information.
[1]
[1] From function virtual void Foam::processorPolyPatch::calcGeometry(Foam::PstreamBuffers&)
[1] in file meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C at line 284.
[1]
FOAM parallel run exiting
[1]
I have tried to use different type of decomposition (simple, scotch), but still the same problem. I have also tried to reduce the mergeTolerance to 1e-04 instead of using the recommendet value of 1e-06.
Would someone be so kind and enlighten me, please? Since I want to make a parameter study, I need to run the meshing in parallel, otherwise it would take too long for me...
I have compressed my blockMeshDict, decomposeParDict and my snappyHexMeshDict to system.zip:
|