|
[Sponsors] |
March 15, 2012, 05:05 |
3D axisymmetric meshing
|
#1 |
Member
Aru
Join Date: Feb 2012
Location: Chennai, India
Posts: 40
Rep Power: 14 |
hi Foamers
I am a beginner in OpenFOAM and am trying to create a mesh for solving roundjet problem. I intend to make a unit cube with a 1 cm diameter inlet and 1 cm diameter outlet (trial purpose). I created the blockMeshdict as follows: convertToMeters 1; vertices ( (0 0 0) //0 // (1 0 0) //1 // (1 1 0) //2 // (0 1 0) //3 // (0 0 1) //4 // (1 0 1) //5 // (1 1 1) //6 // (0 1 1) //7 // (0 0.5 0.505) //8// (0 0.505 0.5) // 9// (0 0.5 0.455) //10// (0 0.455 0.5) // 11 // (1 0.5 0.505) //12// (1 0.505 0.5) // 13// (1 0.5 0.455) //14// (1 0.455 0.5) // 15 // (0 0.5 0.5) //16// (1 0.5 0.5) // 17 // ); blocks ( hex (0 4 11 10 1 5 15 14) (10 10 10) simpleGrading (1 1 1) hex (4 7 8 11 5 6 12 15) (10 10 10) simpleGrading (1 1 1) hex (7 3 9 8 6 2 13 12) (10 10 10) simpleGrading (1 1 1) hex (3 9 10 0 2 14 15 1) (10 10 10) simpleGrading (1 1 1) hex (8 12 13 9 16 17 17 16) (50 2 2) simpleGrading (1 1 1) hex (9 13 14 10 16 17 17 16) (50 2 2) simpleGrading (1 1 1) hex (10 14 15 11 16 17 17 16) (50 2 2) simpleGrading (1 1 1) hex (11 15 12 8 16 17 17 16) (50 2 2) simpleGrading (1 1 1) ); edges ( arc 8 9 (0 0.50353553 0.50353553) arc 9 10 (0 0.50353553 0.49646447) arc 10 11 (0 0.49646447 0.49646447) arc 11 8 (0 0.49646447 0.50353553) arc 12 13 (1 0.50353553 0.50353553) arc 13 14 (1 0.50353553 0.49646447) arc 14 15 (1 0.49646447 0.49646447) arc 15 12 (1 0.49646447 0.50353553) ); boundary ( /*movingWall { type wall; faces ( (3 7 6 2) ); }*/ fixedWalls { type wall; faces ( (0 1 2 3) (4 5 6 7) (7 6 2 3) (4 5 1 0) (3 9 10 0) (0 10 11 4) (4 11 8 7) (7 8 9 3) (2 13 14 1) (1 14 15 5) (5 15 12 6) (6 12 13 2) ); } inlet { type patch; faces ( //(0 3 7 4)// (8 9 10 11) ); } oulet { type patch; faces ( //(2 1 5 6)// (12 13 14 15) ); } ); mergePatchPairs ( ); When i executed blockMesh, i encountered a few errors. The first one says, Creating block mesh topology --> FOAM Warning : From function cellModel::mag(const labelList&, const pointField&) in file meshes/meshShapes/cellModel/cellModel.C at line 128 zero or negative pyramid volume: -0.0452479 for face 0 --> FOAM Warning : From function cellModel::mag(const labelList&, const pointField&) in file meshes/meshShapes/cellModel/cellModel.C at line 128 zero or negative pyramid volume: -0.0377479 for face 1 --> FOAM Warning : From function cellModel::mag(const labelList&, const pointField&) in file meshes/meshShapes/cellModel/cellModel.C at line 128 zero or negative pyramid volume: -0.0795833 for face 2 --> FOAM Warning : From function cellModel::mag(const labelList&, const pointField&) in file meshes/meshShapes/cellModel/cellModel.C at line 128 zero or negative pyramid volume: -0.0034125 for face 3 --> FOAM Warning : From function cellModel::mag(const labelList&, const pointField&) in file meshes/meshShapes/cellModel/cellModel.C at line 128 zero or negative pyramid volume: -0.0414979 for face 4 --> FOAM Warning : From function cellModel::mag(const labelList&, const pointField&) in file meshes/meshShapes/cellModel/cellModel.C at line 128 zero or negative pyramid volume: -0.0414979 for face 5 --> FOAM Warning : From function blockMesh::createTopology(IOdictionary&) in file blockMesh/blockMeshTopology.C at line 255 negative volume block : 0, probably defined inside-out I encountered the same error for all the 8 blocks defined. How can i solve this error? What is the reason for this? Does this say the volume is negative? I am defenitely not sure whether the blockMeshDict i made is right. I apologise for such a long description.My intention was to make everything clear at the very beginning itself. |
|
March 25, 2012, 04:05 |
|
#2 |
Member
Vishal Achasrya
Join Date: Nov 2011
Posts: 38
Rep Power: 15 |
look at the documentation webpage about blockmesh, in that you will see that edges have a direction sense. If you create a block, that implies a certain direction sense for the edges it contains. Now when you create an edge using arc, the direction sense needs to be maintained.. which means that
arc 0 1 () and arc 1 0 () are not the same... so based on the way you have set up your hex blocks, figure out the direction of the arcs and appropriately use the right vertex order. I had the same problem and i fixed it for my case. Let me know if this helped. |
|
March 25, 2012, 16:43 |
|
#3 |
Disabled
Join Date: Mar 2011
Posts: 174
Rep Power: 15 |
It also helps a lot to try creating one block at a time.
Comment all the blocks except the first one, run blockMesh, note the warning, open the geometry in paraFoam regardless of the warning and see what's wrong. If something is inside out, it is coloured in grey instead of white in paraFoam, as far as I remember. Once you are finished, proceed with the rest. |
|
Tags |
blockmesh, meshing 3d, openfoam 2.1 |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[ANSYS Meshing] Icem meshing queries | Mfaizan | ANSYS Meshing & Geometry | 21 | July 31, 2015 00:10 |
[ANSYS Meshing] Axisymmetric Meshing | archer545 | ANSYS Meshing & Geometry | 0 | March 3, 2012 08:36 |
[ANSYS Meshing] Meshing strategy for External Flows | Hybrid | ANSYS Meshing & Geometry | 0 | January 24, 2012 15:27 |
Best Meshing scheme for Cylinder | Nutrex | Main CFD Forum | 4 | July 29, 2008 12:03 |
Volume Meshing & Face Meshing? singularity of grid | ken | FLUENT | 0 | September 4, 2003 12:08 |