 |
 |
 |
|
|
|
[Sponsors] | ||||
[Other] Wedge patch '*' is not planar |
16Likes
 |
|
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 18, 2015, 13:41
|
Wedge patch '*' is not planar
|
#1 | |
|
New Member
Chris
Join Date: Jun 2015
Posts: 3
Rep Power: 11  |
Hi,
This is about a 2D axissymetric case. I was using a Matlab script to find the edges of the blocks I want to build. (you find the Matlab script in the attachment) checkMesh gives the output: Quote:
blockMeshDict: Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 1;
vertices
(
(6.497773e-02 -1.701502e-03 0)
(6.497773e-02 1.701502e-03 0)
(6.517266e-02 1.706606e-03 0)
(6.517266e-02 -1.706606e-03 0)
(6.497773e-02 -1.701502e-03 2.000000e-03)
(6.497773e-02 1.701502e-03 2.000000e-03)
(6.517266e-02 1.706606e-03 2.000000e-03)
(6.517266e-02 -1.706606e-03 2.000000e-03)
(6.497773e-02 -1.701502e-03 2.200000e-03)
(6.497773e-02 1.701502e-03 2.200000e-03)
(6.517266e-02 1.706606e-03 2.200000e-03)
(6.517266e-02 -1.706606e-03 2.200000e-03)
(6.537259e-02 -1.711842e-03 2.000000e-03)
(6.537259e-02 1.711842e-03 2.000000e-03)
(7.180039e-02 1.880159e-03 2.000000e-03)
(7.180039e-02 -1.880159e-03 2.000000e-03)
(6.537259e-02 -1.711842e-03 2.200000e-03)
(6.537259e-02 1.711842e-03 2.200000e-03)
(7.180039e-02 1.880159e-03 2.200000e-03)
(7.180039e-02 -1.880159e-03 2.200000e-03)
(6.497773e-02 -1.701502e-03 3.220000e-02)
(6.497773e-02 1.701502e-03 3.220000e-02)
(7.180039e-02 1.880159e-03 3.220000e-02)
(7.180039e-02 -1.880159e-03 3.220000e-02)
);
blocks
(
//Block 0
hex (0 3 2 1 4 7 6 5) (1 1 1) simpleGrading
(
1
1
1
)
//Block 1
hex (4 7 6 5 8 11 10 9) (1 1 1) simpleGrading
(
1
1
1
)
//Block 2
hex (7 12 13 6 11 16 17 10) (1 1 1) simpleGrading
(
1
1
1
)
//Block 3
hex (12 15 14 13 16 19 18 17) (1 1 1) simpleGrading
(
1
1
1
)
//Block 4
hex (8 19 18 9 20 23 22 21) (1 1 1) simpleGrading
(
1
1
1
)
);
edges();
boundary
(
inlet
{
type patch;
faces
(
(0 1 2 3) // von innen gesehen gegen den Uhrzeigersinn Aufgezaehlt
);
}
outlet
{
type patch;
faces
(
(20 23 22 21)
);
}
shaftInHole
{
type wall;
faces
(
(0 1 5 4)
);
}
shaftInFreestream
{
type wall;
faces
(
(4 5 9 8)
(8 9 21 20)
);
}
wallInHole
{
type wall;
faces
(
(3 7 6 2)
);
}
wallInFreestream
{
type wall;
faces
(
(14 15 19 18)
(18 19 23 22)
);
}
step
{
type wall;
faces
(
(6 7 12 13)
(14 13 12 15)
);
}
front
{
type wedge;
faces
(
(1 2 6 5)
(5 6 10 9)
(6 13 17 10)
(13 14 18 17)
(9 18 22 21)
);
}
back
{
type wedge;
faces
(
(0 4 7 3)
(4 8 11 7)
(7 11 16 12)
(12 16 19 15)
(8 20 23 19)
);
}
);
mergePatchPairs();
// ************************************************************************* //
Code:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.4.0-f0842aea0e77
Exec : blockMesh
Date : Jun 18 2015
Time : 18:36:28
Host : "ChrisHP"
PID : 11568
Case : /home/chris/OpenFOAM/nowak/Netzerzeugung/NeuesMesh
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Creating block mesh from
"/home/chris/OpenFOAM/nowak/Netzerzeugung/NeuesMesh/constant/polyMesh/blockMeshDict"
No non-linear edges defined
Creating topology blocks
Creating topology patches
Creating block mesh topology
--> FOAM Warning :
From function polyMesh::polyMesh(... construct from shapes...)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 903
Found 4 undefined faces in mesh; adding to default patch.
--> FOAM Warning :
From function wedgePolyPatch::calcGeometry(PstreamBuffers&)
in file meshes/polyMesh/polyPatches/constraint/wedge/wedgePolyPatch.C at line 72
Wedge patch 'front' is not planar.
At local face at (0.0650752 0.00170405 0.001) the normal (0.0261748 -0.999657 1.32563e-19) differs from the average normal (0.0261767 -0.999657 3.31431e-20) by 3.60957e-12
Either correct the patch or split it into planar parts
--> FOAM Warning :
From function wedgePolyPatch::calcGeometry(PstreamBuffers&)
in file meshes/polyMesh/polyPatches/constraint/wedge/wedgePolyPatch.C at line 72
Wedge patch 'front' is not planar.
At local face at (0.0650752 0.00170405 0.0021) the normal (0.0261748 -0.999657 1.32563e-18) differs from the average normal (0.0261767 -0.999657 3.31431e-20) by 3.60957e-12
Either correct the patch or split it into planar parts
... (Removed a few similiar Errormessages)
--> FOAM Warning :
From function wedgePolyPatch::calcGeometry(PstreamBuffers&)
in file meshes/polyMesh/polyPatches/constraint/wedge/wedgePolyPatch.C at line 72
Wedge patch 'back' is not planar.
At local face at (0.0685865 -0.001796 0.0021) the normal (0.0261768 0.999657 0) differs from the average normal (0.0261767 0.999657 -3.31431e-20) by 1.61774e-14
Either correct the patch or split it into planar parts
--> FOAM Warning :
From function wedgePolyPatch::calcGeometry(PstreamBuffers&)
in file meshes/polyMesh/polyPatches/constraint/wedge/wedgePolyPatch.C at line 72
Wedge patch 'back' is not planar.
At local face at (0.0683891 -0.00179083 0.0172) the normal (0.0261769 0.999657 0) differs from the average normal (0.0261767 0.999657 -3.31431e-20) by 2.83042e-14
Either correct the patch or split it into planar parts
Check topology
Basic statistics
Number of internal faces : 3
Number of boundary faces : 24
Number of defined boundary faces : 24
Number of undefined boundary faces : 0
Checking patch -> block consistency
Creating block offsets
Creating merge list .
Creating polyMesh from blockMesh
Creating patches
Creating cells
Creating points with scale 1
Block 0 cell size :
i : 0.000194997 .. 0.000194997
j : 0.003403 .. 0.003403
k : 0.002 .. 0.002
Block 1 cell size :
i : 0.000194997 .. 0.000194997
j : 0.003403 .. 0.003403
k : 0.0002 .. 0.0002
Block 2 cell size :
i : 0.000199999 .. 0.000199999
j : 0.00341321 .. 0.00341321
k : 0.0002 .. 0.0002
Block 3 cell size :
i : 0.00643 .. 0.00643
j : 0.00342368 .. 0.00342368
k : 0.0002 .. 0.0002
Block 4 cell size :
i : 0.006825 .. 0.006825
j : 0.003403 .. 0.003403
k : 0.03 .. 0.03
--> FOAM Warning :
From function wedgePolyPatch::calcGeometry(PstreamBuffers&)
in file meshes/polyMesh/polyPatches/constraint/wedge/wedgePolyPatch.C at line 72
Wedge patch 'front' is not planar.
At local face at (0.0650752 0.00170405 0.001) the normal (-0.0261748 0.999657 1.32563e-19) differs from the average normal (-0.0261767 0.999657 3.31431e-20) by 3.60957e-12
Either correct the patch or split it into planar parts
... (Removed a few similiar Errormessages)
--> FOAM Warning :
From function wedgePolyPatch::calcGeometry(PstreamBuffers&)
in file meshes/polyMesh/polyPatches/constraint/wedge/wedgePolyPatch.C at line 72
Wedge patch 'back' is not planar.
At local face at (0.0683891 -0.00179083 0.0172) the normal (-0.0261769 -0.999657 0) differs from the average normal (-0.0261767 -0.999657 -3.31431e-20) by 2.83042e-14
Either correct the patch or split it into planar parts
There are no merge patch pairs edges
Writing polyMesh
----------------
Mesh Information
----------------
boundingBox: (0.0649777 -0.00188016 0) (0.0718004 0.00188016 0.0322)
nPoints: 28
nCells: 5
nFaces: 27
nInternalFaces: 3
----------------
Patches
----------------
patch 0 (start: 3 size: 1) name: inlet
patch 1 (start: 4 size: 1) name: outlet
patch 2 (start: 5 size: 1) name: shaftInHole
patch 3 (start: 6 size: 2) name: shaftInFreestream
patch 4 (start: 8 size: 1) name: wallInHole
patch 5 (start: 9 size: 2) name: wallInFreestream
patch 6 (start: 11 size: 2) name: step
patch 7 (start: 13 size: 5) name: front
patch 8 (start: 18 size: 5) name: back
patch 9 (start: 23 size: 4) name: defaultFaces
End
Code:
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.4.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 2.4.0-f0842aea0e77 Exec : checkMesh Date : Jun 18 2015 Time : 18:36:40 Host : "ChrisHP" PID : 11574 Case : /home/chris/OpenFOAM/nowak/Netzerzeugung/NeuesMesh nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Allowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create polyMesh for time = 0 --> FOAM Warning : From function wedgePolyPatch::calcGeometry(PstreamBuffers&) in file meshes/polyMesh/polyPatches/constraint/wedge/wedgePolyPatch.C at line 72 Wedge patch 'front' is not planar. At local face at (0.0650752 0.00170405 0.001) the normal (-0.0261748 0.999657 1.32563e-19) differs from the average normal (-0.0261767 0.999657 3.31431e-20) by 3.60957e-12 Either correct the patch or split it into planar parts --> FOAM Warning : From function wedgePolyPatch::calcGeometry(PstreamBuffers&) in file meshes/polyMesh/polyPatches/constraint/wedge/wedgePolyPatch.C at line 72 Wedge patch 'front' is not planar. At local face at (0.0650752 0.00170405 0.0021) the normal (-0.0261748 0.999657 1.32563e-18) differs from the average normal (-0.0261767 0.999657 3.31431e-20) by 3.60957e-12 Either correct the patch or split it into planar parts --> FOAM Warning : From function wedgePolyPatch::calcGeometry(PstreamBuffers&) in file meshes/polyMesh/polyPatches/constraint/wedge/wedgePolyPatch.C at line 72 Wedge patch 'front' is not planar. At local face at (0.0652726 0.00170922 0.0021) the normal (-0.0261802 0.999657 -1.29248e-18) differs from the average normal (-0.0261767 0.999657 3.31431e-20) by 1.22804e-11 Either correct the patch or split it into planar parts ... (Removed a few similiar Errormessages) --> FOAM Warning : From function wedgePolyPatch::calcGeometry(PstreamBuffers&) in file meshes/polyMesh/polyPatches/constraint/wedge/wedgePolyPatch.C at line 72 Wedge patch 'back' is not planar. At local face at (0.0683891 -0.00179083 0.0172) the normal (-0.0261769 -0.999657 0) differs from the average normal (-0.0261767 -0.999657 -3.31431e-20) by 2.83042e-14 Either correct the patch or split it into planar parts Time = 0 Mesh stats points: 28 internal points: 0 faces: 27 internal faces: 3 cells: 5 faces per cell: 6 boundary patches: 10 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 5 prisms: 0 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 0 Checking topology... Boundary definition OK. ***Total number of faces on empty patches is not divisible by the number of cells in the mesh. Hence this mesh is not 1D or 2D. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. *Number of regions: 2 The mesh has multiple regions which are not connected by any face. <<Writing region information to "0/cellToRegion" <<Writing region 0 with 4 cells to cellSet region0 <<Writing region 1 with 1 cells to cellSet region1 Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology inlet 1 4 ok (non-closed singly connected) outlet 1 4 ok (non-closed singly connected) shaftInHole 1 4 ok (non-closed singly connected) shaftInFreestream 2 8 ok (non-closed singly connected) wallInHole 1 4 ok (non-closed singly connected) wallInFreestream 2 8 ok (non-closed singly connected) step 2 6 ok (non-closed singly connected) front 5 14 ok (non-closed singly connected) back 5 14 ok (non-closed singly connected) defaultFaces 4 12 ok (non-closed singly connected) Checking geometry... Overall domain bounding box (0.0649777 -0.00188016 0) (0.0718004 0.00188016 0.0322) Mesh (non-empty, non-wedge) directions (1 0 1) Mesh (non-empty) directions (1 1 0) Wedge front with angle 1.49998 degrees Wedge back with angle 1.49998 degrees All edges aligned with or perpendicular to non-empty directions. Boundary openness (2.05835e-17 0 0) OK. Max cell openness = 8.2747e-17 OK. Max aspect ratio = 30.7004 OK. Minimum face area = 3.89994e-08. Maximum face area = 0.00020475. Face area magnitudes OK. Min volume = 1.32868e-10. Max volume = 7.33094e-07. Total volume = 7.3931e-07. Cell volumes OK. Mesh non-orthogonality Max: 0 average: 0 Non-orthogonality check OK. Face pyramids OK. Max skewness = 0.0653982 OK. Coupled point location match (average 0) OK. Mesh OK. Time = 1 Mesh stats points: 28 internal points: 0 faces: 27 internal faces: 3 cells: 5 faces per cell: 6 boundary patches: 10 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 5 prisms: 0 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 0 Checking topology... Boundary definition OK. ***Total number of faces on empty patches is not divisible by the number of cells in the mesh. Hence this mesh is not 1D or 2D. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. *Number of regions: 2 The mesh has multiple regions which are not connected by any face. <<Writing region information to "1/cellToRegion" <<Writing region 0 with 1 cells to cellSet region0 <<Writing region 1 with 4 cells to cellSet region1 Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology inlet 1 4 ok (non-closed singly connected) outlet 1 4 ok (non-closed singly connected) shaftInHole 1 4 ok (non-closed singly connected) shaftInFreestream 2 8 ok (non-closed singly connected) wallInHole 1 4 ok (non-closed singly connected) wallInFreestream 2 8 ok (non-closed singly connected) step 2 6 ok (non-closed singly connected) front 5 14 ok (non-closed singly connected) back 5 14 ok (non-closed singly connected) defaultFaces 4 12 ok (non-closed singly connected) Checking geometry... Overall domain bounding box (0.0649777 -0.00188016 0) (0.0718004 0.00188016 0.0322) Mesh (non-empty, non-wedge) directions (1 0 1) Mesh (non-empty) directions (1 1 0) Wedge front with angle 1.49998 degrees Wedge back with angle 1.49998 degrees All edges aligned with or perpendicular to non-empty directions. Boundary openness (2.05835e-17 0 0) OK. Max cell openness = 8.2747e-17 OK. Max aspect ratio = 30.7004 OK. Minimum face area = 3.89994e-08. Maximum face area = 0.00020475. Face area magnitudes OK. Min volume = 1.32868e-10. Max volume = 7.33094e-07. Total volume = 7.3931e-07. Cell volumes OK. Mesh non-orthogonality Max: 0 average: 0 Non-orthogonality check OK. Face pyramids OK. Max skewness = 0.0653982 OK. Coupled point location match (average 0) OK. Mesh OK. End |
||
|
|
||
|
June 20, 2015, 04:35
|
|
#2 |
|
New Member
Chris
Join Date: Jun 2015
Posts: 3
Rep Power: 11 |
Ok, Problem was solved by setting the angle to 0 and all y values to -1 or +1. I also changed the type of the front and back surfaces to patch. So I got a 2D plane mesh. In a new case dir I executed extrudeMesh with the wedge options using the front and back patches as input. It worked =)
|
|
|
|
|
|
|
October 6, 2015, 11:26
|
makeAxialMesh
|
#3 |
|
New Member
Antony Premkumar
Join Date: Jul 2015
Posts: 19
Rep Power: 11 |
I couldn't really understood how you solved it.
Could you please explain me what angle did you change to 0. I am also facing the same problem. Plus my check mesh is showing mesh not OK with seven errors |
|
|
|
|
|
|
November 24, 2015, 19:32
|
|
#4 |
|
Member
Alex
Join Date: Jun 2011
Posts: 33
Rep Power: 15 |
Antony: I had some similar issues recently with axisymmetric meshes. I had success switching the writeFormat in controlDict from ascii to binary.
I guess there is some rounding introduced when you convert back and forth from human readable format. It would be nice if there was some place to adjust the precision requirement for wedge planar-ness. |
|
|
|
|
|
|
May 21, 2016, 08:50
|
Ignore and proceed to computation
|
#5 |
|
New Member
Smrutimaya Pattanaik
Join Date: May 2016
Location: Chennai, India
Posts: 11
Rep Power: 0 |
Those error messages are not very serious, ignore and proceed to computation.
|
|
|
|
|
|
|
February 1, 2018, 16:51
|
|
#6 |
|
Member
Foad
Join Date: Aug 2017
Posts: 58
Rep Power: 9 |
I had a same issue with an axisymmetric geometry/mesh I created with Gmsh. I reran the gmshToFoam and edited the constant/polyMesh/boundary file manually, changing the wedges type and physicalType to wedge and the noslip walls to wall. it solved the problem and ran smoothly.
|
|
|
|
|
|
|
January 12, 2021, 06:55
|
|
#7 |
|
New Member
Naveen kumar
Join Date: Feb 2019
Posts: 1
Rep Power: 0 |
I want to do 2d axisymmertic meshing in OpenFoam. But it is showing this error "Either correct the patch or split it into planar parts". Please have a look on it.
|
|
|
|
|
|
|
September 18, 2024, 06:52
|
|
#8 |
|
New Member
Join Date: Nov 2021
Posts: 1
Rep Power: 0 |
You can try increasing the writePrecision in controlDict (writePrecision 10; ). It worked for me.
|
|
|
|
|
|
|
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| y+ and u+ values with low-Re RANS turbulence models: utility + testcase | florian_krause | OpenFOAM | 114 | August 23, 2023 06:37 |
| steadyUniversalMRFFoam Tutorial fails in MixingPlane | HenrikJohansson | OpenFOAM Bugs | 0 | February 14, 2019 05:48 |
| Near wall treatment in k-omega SST | Arnoldinho | OpenFOAM Running, Solving & CFD | 38 | March 8, 2017 14:48 |
| [mesh manipulation] Importing Multiple Meshes | thomasnwalshiii | OpenFOAM Meshing & Mesh Conversion | 18 | December 19, 2015 19:57 |
| [Gmsh] Import gmsh msh to Foam | adorean | OpenFOAM Meshing & Mesh Conversion | 24 | April 27, 2005 09:19 |
Linear Mode
