CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Open Source Meshers: Gmsh, Netgen, CGNS, ...

[wing mesh with gmsh]

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 26, 2013, 04:16
Default [wing mesh with gmsh]
  #1
New Member
 
David Doose
Join Date: Jan 2013
Posts: 5
Rep Power: 4
huitetquatre is on a distinguished road
Hi,

My main objective is to "create" a wing. I've already analysed different Naca profiles without any problems. The second step is to analyse a complete wing.
I'm actually generating points with a specific pg, then the .geo file is generated, then the mesh is created by gmsh and I finally apply the gmshToFoam command.

The generated mesh seams to be good but the checkMesh failed.

A simple geo sample:
Code:
xmin =  0;
xmax =  1;
ymin = -1;
ymax =  1;
zmin = -1;
zmax =  2;

box_lc = 0.5;
lc = 0.01;

Point(1) = { xmin, ymin, zmin, box_lc};
Point(2) = { xmin, ymax, zmin, box_lc};
Point(3) = { xmin, ymax, zmax, box_lc};
Point(4) = { xmin, ymin, zmax, box_lc};

Line(1) = {1, 2};
Line(2) = {2, 3};
Line(3) = {3, 4};
Line(4) = {4, 1};

Line Loop(1) = {1, 2, 3, 4};

Plane Surface(1) = {1};

tmp[] = Extrude {xmax-xmin, 0, 0} {
  Surface{1};
};
Delete { Volume{tmp[1]}; }
Delete { Surface{1}; }

// ----------
l = 0.2;
h = 0.2;
w = 0.5;

Point(101) = {xmin, -l, 0, lc};
Point(102) = {xmin,  l, 0, lc};
Point(103) = {xmin,  0, h, lc};

Point(201) = {xmin+w, -l, 0, lc};
Point(202) = {xmin+w,  l, 0, lc};
Point(203) = {xmin+w,  0, h, lc};

Line(101) = {101, 102};
Line(102) = {102, 103};
Line(103) = {103, 101};

Line(201) = {201, 202};
Line(202) = {202, 203};
Line(203) = {203, 201};

Line(301) = {101, 201};
Line(302) = {102, 202};
Line(303) = {103, 203};

Line Loop(2) = {101, 102, 103};
Line Loop(3) = {201, 202, 203};

Line Loop(11) = {301, 201, -302, -101};
Line Loop(12) = {302, 202, -303, -102};
Line Loop(13) = {303, 203, -301, -103};


Plane Surface(100) = {1, 2}; // Left
Plane Surface(101) = {3};    // Wing right

Plane Surface(201) = {11};
Plane Surface(202) = {12};
Plane Surface(203) = {13};

// Right = 26
// Front = 25
// Top   = 21
// Back  = 17
// Bot   = 13

Surface Loop(1000) = {26, 25, 21, 17, 13, 100, 201, 202, 203, 101};
Volume (1000) = {1000};

Physical Surface("front")   = {25};
Physical Surface("back")    = {17};
Physical Surface("top")     = {21};
Physical Surface("bottom")  = {13};
Physical Surface("left")    = {100};
Physical Surface("right")   = {26};

Physical Surface("foil")    = {201, 202, 203, 101};

Physical Volume("internal") = {1000};
checkMesh -allTopology -allGeometry :
Code:
Create polyMesh for time = 0

Time = 0

Mesh stats
    points:           8990
    faces:            72001
    internal faces:   58059
    cells:            32515
    boundary patches: 7
    point zones:      0
    face zones:       0
    cell zones:       1

Overall number of cells of each type:
    hexahedra:     0
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    32515
    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.
    Topological cell zip-up check OK.
    Face-face connectivity OK.
  <<Writing 103 cells with with two non-boundary faces to set twoInternalFacesCells
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
    Patch               Faces    Points   Surface topology                   Bounding box
    bottom              8        9        ok (non-closed singly connected)   (0 -1 -1) (1 1 -1)
    back                16       15       ok (non-closed singly connected)   (0 1 -1) (1 1 2)
    top                 6        8        ok (non-closed singly connected)   (0 -1 2) (1 1 2)
    front               16       15       ok (non-closed singly connected)   (0 -1 -1) (1 -1 2)
    right               44       31       ok (non-closed singly connected)   (1 -1 -1) (1 1 2)
    left                607      359      ok (non-closed singly connected)   (0 -1 -1) (0 1 2)
    foil                13245    6671     ok (non-closed singly connected)   (0 -0.2 0) (0.5 0.2 0.2)

Checking geometry...
    Overall domain bounding box (0 -1 -1) (1 1 2)
    Mesh (non-empty, non-wedge) directions (0 1 1)
    Mesh (non-empty) directions (0 1 1)
 ***Number of edges not aligned with or perpendicular to non-empty directions: 45693
  <<Writing 8990 points on non-aligned edges to set nonAlignedEdges
    Boundary openness (-2.48715e-18 1.85141e-18 6.14172e-19) OK.
    Max cell openness = 1.42973e-16 OK.
    Max aspect ratio = 31.7252 OK.
    Minumum face area = 1.67583e-05. Maximum face area = 0.388435.  Face area magnitudes OK.
    Min volume = 8.84634e-10. Max volume = 0.043476.  Total volume = 5.98.  Cell volumes OK.
    Mesh non-orthogonality Max: 86.2608 average: 29.571
   *Number of severely non-orthogonal faces: 314.
    Non-orthogonality check OK.
  <<Writing 314 non-orthogonal faces to set nonOrthoFaces
    Face pyramids OK.
 ***Max skewness = 6.06065, 96 highly skew faces detected which may impair the quality of the results
  <<Writing 96 skew faces to set skewFaces
    Coupled point location match (average 0) OK.
    Face tets OK.
    Min/max edge length = 0.00553551 1.20185 OK.
    All angles in faces OK.
    All face flatness OK.
    Cell determinant (wellposedness) : minimum: 3.80836e-06 average: 0.494671
 ***Cells with small determinant found, number of cells: 133
  <<Writing 133 under-determined cells to set underdeterminedCells
    Concave cell check OK.

Failed 3 mesh checks.

End

What's wrong ? How can I solve this problem ?
huitetquatre is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[ICEM] surface mesh merging problem everest ANSYS Meshing & Geometry 39 June 5, 2013 19:02
Inner geometry gets lost exporting mesh from ICEM CFD to CFX-Pre powpow CFX 3 December 20, 2012 10:14
[ICEM] Problem making structural mesh on a surface froztbear ANSYS Meshing & Geometry 1 November 10, 2011 09:52
snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Native Meshers: snappyHexMesh and Others 2 March 27, 2011 21:11
basic of mesh refinement arya CFX 4 June 19, 2007 12:21


All times are GMT -4. The time now is 21:35.