CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[Other] Wedge patch '*' is not planar

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree7Likes
  • 6 Post By liquidspoon
  • 1 Post By foadsf

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 18, 2015, 12:41
Default Wedge patch '*' is not planar
  #1
New Member
 
Chris
Join Date: Jun 2015
Posts: 3
Rep Power: 7
LilumDaru is on a distinguished road
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:
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
I can't get rid of these error messages. I already tried to change the writePrecision and the writeFormat without any success. Do these errors matter? Can you please help me to eliminate these errors?


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();

// ************************************************************************* //
blockMesh:
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
checkMesh:
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
Attached Files
File Type: zip matlab_scripts.zip (1.3 KB, 48 views)
LilumDaru is offline   Reply With Quote

Old   June 20, 2015, 03:35
Default
  #2
New Member
 
Chris
Join Date: Jun 2015
Posts: 3
Rep Power: 7
LilumDaru is on a distinguished road
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 =)
LilumDaru is offline   Reply With Quote

Old   October 6, 2015, 10:26
Default makeAxialMesh
  #3
New Member
 
Antony Ashley Raja B
Join Date: Jul 2015
Posts: 16
Rep Power: 7
Antony_Bosco is on a distinguished road
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
Antony_Bosco is offline   Reply With Quote

Old   November 24, 2015, 18:32
Default
  #4
Member
 
Alex
Join Date: Jun 2011
Posts: 32
Rep Power: 11
liquidspoon is on a distinguished road
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.
Artur, Icemaan, Danillo and 3 others like this.
liquidspoon is offline   Reply With Quote

Old   May 21, 2016, 07:50
Smile Ignore and proceed to computation
  #5
New Member
 
Smrutimaya Pattanaik
Join Date: May 2016
Location: Chennai, India
Posts: 11
Rep Power: 0
smrutimaya is on a distinguished road
Those error messages are not very serious, ignore and proceed to computation.
smrutimaya is offline   Reply With Quote

Old   February 1, 2018, 15:51
Default
  #6
Member
 
Foad
Join Date: Aug 2017
Posts: 40
Rep Power: 5
foadsf is on a distinguished road
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.
Sun likes this.
foadsf is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
y+ and u+ values with low-Re RANS turbulence models: utility + testcase florian_krause OpenFOAM 113 June 25, 2019 00:58
steadyUniversalMRFFoam Tutorial fails in MixingPlane HenrikJohansson OpenFOAM Bugs 0 February 14, 2019 04:48
Near wall treatment in k-omega SST Arnoldinho OpenFOAM Running, Solving & CFD 38 March 8, 2017 13:48
[mesh manipulation] Importing Multiple Meshes thomasnwalshiii OpenFOAM Meshing & Mesh Conversion 18 December 19, 2015 18:57
[Gmsh] Import gmsh msh to Foam adorean OpenFOAM Meshing & Mesh Conversion 24 April 27, 2005 08:19


All times are GMT -4. The time now is 04:37.