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

[blockMesh] blockMesh error for 3d pipe

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By hxaxtma

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 2, 2015, 18:41
Default blockMesh error for 3d pipe
  #1
Member
 
nsreddy
Join Date: Sep 2010
Posts: 40
Rep Power: 15
nsreddysrsit is on a distinguished road
hi,
I am creating a blockMesh for pipe, here i am attaching geometry as well , I am getting error.
here is the Mesh file.
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.01;

vertices
(
(0 0 0) //0
(0.15 0.15 0)
(-0.15 0.15 0)
(-0.15 -0.15 0)
(0.15 -0.15 0)
(0.3 0.3 0)//5
(-0.3 0.3 0)
(-0.3 -0.3 0)
(0.3 -0.3 0)
(0.3535 0.3535 0)//9
(-0.3535 0.3535 0)
(-0.3535 -0.3535 0)
(0.3535 -0.3535 0)
(0.15 0.15 30)//13
(-0.15 0.15 30)
(-0.15 -0.15 30)
(0.15 -0.15 30)
(0.3 0.3 30) //17
(-0.3 0.3 30)
(-0.3 -0.3 30)
(0.3 -0.3 30)
(0.3535 0.3535 0) //21
(-0.3535 0.3535 0)
(-0.3535 -0.3535 0)
(0.3535 -0.3535 0)
(0 0 30) //25

);

blocks
(
hex (1 2 3 4 13 14 15 16) (10 10 150) simpleGrading (1 1 1)
hex (1 5 6 2 13 17 18 14) (10 10 150) simpleGrading (1 1 1)
hex (2 6 7 3 14 18 19 15) (10 10 150) simpleGrading (1 1 1)
hex (3 7 8 4 15 19 20 16) (10 10 150) simpleGrading (1 1 1)
hex (4 8 5 1 16 20 17 13) (10 10 150) simpleGrading (1 1 1)
hex (5 9 10 6 17 21 22 18) (10 10 150) simpleGrading (1 1 1)
hex (6 10 11 7 18 22 23 19) (10 10 150) simpleGrading (1 1 1)
hex (7 11 12 8 19 23 24 20) (10 10 150) simpleGrading (1 1 1)
hex (12 9 5 8 24 21 17 20) (10 10 150) simpleGrading (1 1 1)


);

edges
(
arc 5 6 (0 0.42426 0)
arc 6 7 (-0.42426 0 0)
arc 7 8 (0 -0.42426 0)
arc 8 5 (0.42426 0 0)
arc 9 10 (0 0.5 0)
arc 10 11 (-0.5 0 0)
arc 11 12 ( 0 -0.5 0)
arc 12 9 (0.5 0 0)
arc 17 18 (0 0.42426 30)
arc 18 19 (-0.42426 0 30)
arc 19 20 (0 -0.42426 30)
arc 20 17 (0.42426 0 30)
arc 21 22 (0 0.5 30)
arc 22 23 (-0.5 0 30)
arc 23 24 (0 -0.5 30)
arc 24 21 (0.5 0 30)

);
boundary
(
inlet
{
type inlet;
faces
(
(13 16 15 14)
(13 14 18 17)
(14 15 19 18)
(15 16 20 19)
(16 13 17 20)
(17 18 22 21)
(18 19 23 22)
(19 20 24 23)
(20 17 21 24)

);
}
outlet
{
type outlet;
faces
(
(1 5 6 2)
(2 6 7 3)
(3 7 8 4)
(4 8 5 1)
(1 2 3 4)
(5 9 10 6)
(6 10 11 7)
(7 11 12 8)
(8 12 9 5)


);
}
fixwdWalls
{
type fixedWalls;
faces
(
(21 22 10 9)
(22 23 11 10)
(23 24 12 11)
(24 21 9 12)
);
}
);

mergePatchPairs
(
);

// ************************************************** *********************** //
Regards
Attached Images
File Type: png Screenshot from 2015-02-03 05:00:23.png (70.0 KB, 120 views)
nsreddysrsit is offline   Reply With Quote

Old   February 3, 2015, 04:38
Default
  #2
Senior Member
 
Join Date: Jan 2014
Posts: 179
Rep Power: 12
hxaxtma is on a distinguished road
You had some mistakes in your blockMeshDict

- Forgotten curles )
and there are no patch types like inlet and outlet or fxdwalls.mist commonly used are patch, wall, emtpy, symmetryPlane, e.g.....

Here is the corrected Dictionary.

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 0.1;

vertices
(
    (0 0 0) //0
    (0.15 0.15 0)
    (-0.15 0.15 0)
    (-0.15 -0.15 0)
    (0.15 -0.15 0)
    (0.3 0.3 0)//5
    (-0.3 0.3 0)
    (-0.3 -0.3 0)
    (0.3 -0.3 0)
    (0.3535 0.3535 0)//9
    (-0.3535 0.3535 0)
    (-0.3535 -0.3535 0)
    (0.3535 -0.3535 0)
    (0.15 0.15 30)//13
    (-0.15 0.15 30)
    (-0.15 -0.15 30)
    (0.15 -0.15 30)
    (0.3 0.3 30) //17
    (-0.3 0.3 30)
    (-0.3 -0.3 30)
    (0.3 -0.3 30)
    (0.3535 0.3535 0) //21
    (-0.3535 0.3535 0)
    (-0.3535 -0.3535 0)
    (0.3535 -0.3535 0)
    (0 0 30) //25
);

blocks
(
    hex (1 2 3 4 13 14 15 16) (10 10 150) simpleGrading (1 1 1)
    hex (1 5 6 2 13 17 18 14) (10 10 150) simpleGrading (1 1 1)
    hex (2 6 7 3 14 18 19 15) (10 10 150) simpleGrading (1 1 1)
    hex (3 7 8 4 15 19 20 16) (10 10 150) simpleGrading (1 1 1)
    hex (4 8 5 1 16 20 17 13) (10 10 150) simpleGrading (1 1 1)
    hex (5 9 10 6 17 21 22 18) (10 10 150) simpleGrading (1 1 1)
    hex (6 10 11 7 18 22 23 19) (10 10 150) simpleGrading (1 1 1)
    hex (7 11 12 8 19 23 24 20) (10 10 150) simpleGrading (1 1 1)
    hex (12 9 5 8 24 21 17 20) (10 10 150) simpleGrading (1 1 1)
);

edges
(
    arc 5 6 (0 0.42426 0)
    arc 6 7 (-0.42426 0 0)
    arc 7 8 (0 -0.42426 0)
    arc 8 5 (0.42426 0 0)
    arc 9 10 (0 0.5 0)
    arc 10 11 (-0.5 0 0)
    arc 11 12 ( 0 -0.5 0)
    arc 12 9 (0.5 0 0)
    arc 17 18 (0 0.42426 30)
    arc 18 19 (-0.42426 0 30)
    arc 19 20 (0 -0.42426 30)
    arc 20 17 (0.42426 0 30)
    arc 21 22 (0 0.5 30)
    arc 22 23 (-0.5 0 30)
    arc 23 24 (0 -0.5 30)
    arc 24 21 (0.5 0 30)
);

boundary
(
    inlet
    {
        type patch;
        faces
        (
        (13 16 15 14)
        (13 14 18 17)
        (14 15 19 18)
        (15 16 20 19)
        (16 13 17 20)
        (17 18 22 21)
        (18 19 23 22)
        (19 20 24 23)
        (20 17 21 24)    
        );
    }

    fixedWalls
    {
        type wall;
        faces
        (
        (21 22 10 9)
        (22 23 11 10)
        (23 24 12 11)
        (24 21 9 12)
        );
    }

    outlet
    {
        type patch;
        faces
        (
        (1 5 6 2)
        (2 6 7 3)
        (3 7 8 4)
        (4 8 5 1)
        (1 2 3 4)
        (5 9 10 6)
        (6 10 11 7)
        (7 11 12 8)
        (8 12 9 5)
        );
    
    }
);

mergePatchPairs
(
);
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
hxaxtma is offline   Reply With Quote

Old   February 3, 2015, 04:43
Default
  #3
Senior Member
 
Join Date: Jan 2014
Posts: 179
Rep Power: 12
hxaxtma is on a distinguished road
And just a hint: I use following tools while creating blockMeshDicts
-pyFoamDisplayBlockMesh.py
-paraFoam -blocks
hxaxtma is offline   Reply With Quote

Old   February 3, 2015, 06:50
Default hi
  #4
Member
 
nsreddy
Join Date: Sep 2010
Posts: 40
Rep Power: 15
nsreddysrsit is on a distinguished road
thanks for the reply,
After checkMesh I am getting following error.

Checking geometry...
Overall domain bounding box (-1.43321 -1.43321 0) (1.43321 1.43321 3)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-8.33742e-18 -1.47677e-16 -3.01599e-18) OK.
***High aspect ratio cells found, Max aspect ratio: 4896.76, number of cells 2704
<<Writing 2704 cells with high aspect ratio to set highAspectRatioCells
Minimum face area = 7.56604e-07. Maximum face area = 0.247873. Face area magnitudes OK.
Min volume = 2.46125e-08. Max volume = 1.77893e-05. Total volume = 0.2571. Cell volumes OK.
Mesh non-orthogonality Max: 180 average: 52.447
*Number of severely non-orthogonal (> 70 degrees) faces: 147240.
***Number of non-orthogonality errors: 5716.
<<Writing 152956 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 11432 faces are incorrectly oriented.
<<Writing 5716 faces with incorrect orientation to set wrongOrientedFaces
***Max skewness = 9.38368, 36 highly skew faces detected which may impair the quality of the results
<<Writing 36 skew faces to set skewFaces
Coupled point location match (average 0) OK.

Failed 4 mesh checks.

Best Regards
nsreddysrsit is offline   Reply With Quote

Old   February 3, 2015, 13:49
Default
  #5
Senior Member
 
Join Date: Jan 2014
Posts: 179
Rep Power: 12
hxaxtma is on a distinguished road
Well you defined your geometry wrong...

this is what I got, looks like a space fighter from star wars.

I am missing time at the moment. Maybe during the week I can fix it.
Attached Images
File Type: jpg mesh.jpg (14.2 KB, 79 views)
hxaxtma is offline   Reply With Quote

Old   February 3, 2015, 13:53
Default
  #6
Senior Member
 
Join Date: Jan 2014
Posts: 179
Rep Power: 12
hxaxtma is on a distinguished road
Quick Google search,

example Code for "good" pipe mesh. Just adapt it to your data:

Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.7.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 1;

vertices        
(
//internal square
    (0 -1 -1) //0
    (1 -1 -1) //1
    (1 -1  1) //2
    (0 -1  1) //3
    (0  1 -1) //4
    (1  1 -1) //5
    (1  1  1) //6
    (0  1  1) //7
//Going clockwise with sections
//North section
    (0  2  2) //8
    (0 -2  2) //9
    (1  2  2) //10
    (1 -2  2) //11
//East section
    (0  2 -2) //12
    (1  2 -2) //13
//South section
    (0 -2 -2) //14
    (1 -2 -2) //15
);

edges           
(
//North section
   arc  8  9 (0  0     2.82)
   arc 10 11 (1  0     2.82)
//East section
   arc  8 12 (0  2.82  0)
   arc 10 13 (1  2.82  0)
//South section
   arc 12 14 (0  0 -2.82 )
   arc 13 15 (1  0 -2.82 )
//West section
   arc 14  9 (0  -2.82 0)
   arc 15 11 (1  -2.82 0)
);

blocks
(
//Centre
    hex (4 5 6 7 0 1 2 3) (32 48 48) simpleGrading (1 1 1)
//East section
    hex (8 10 11 9 7 6 2 3) (32 48 92) simpleGrading (1 1 0.25)
//East section
    hex (12 13 10 8 4 5 6 7) (32 48 92) simpleGrading (1 1 0.25)
//South section
    hex (14 15 13 12 0 1 5 4) (32 48 92) simpleGrading (1 1 0.25)
//West section
    hex (9 11 15 14 3 2 1 0) (32 48 92) simpleGrading (1 1 0.25)
);

patches         
(
    wall pipeWall
    (
        (9 11 10 8)
        (8 10 13 12)
        (12 14 15 13)
        (14 9 11 15)
    )
    cyclic cNorth
    (
        (9  3 7  8)
        (11 2 6 10)
    )
    cyclic cEast
    (
        (7  4 12  8)
        (6  5 13 10)
    )
    cyclic cSouth
    (
        (0 14 12 4)
        (1 5 13 15)
    )
    cyclic cWest
    (
        (0 14 9 3)
        (1 15 11 2)
    )
    cyclic cCentre
    (
        (0 4 7 3)
        (1 2 6 5)
    )
);

mergePatchPairs
(
);

// ************************************************************************* //
Attached Images
File Type: jpg mesh2.jpg (67.3 KB, 129 views)
alethea24 likes this.
hxaxtma is offline   Reply With Quote

Old   February 5, 2015, 10:12
Default hi
  #7
Member
 
nsreddy
Join Date: Sep 2010
Posts: 40
Rep Power: 15
nsreddysrsit is on a distinguished road
thank you,
i will try
nsreddysrsit is offline   Reply With Quote

Reply


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
[blockMesh] Quarter pipe with blockMesh kaber OpenFOAM Meshing & Mesh Conversion 1 January 10, 2018 21:51
need help about double pipe heat exchanger with chtMultiRegionSimpleFoam wuyangzhen OpenFOAM 10 December 12, 2017 00:19
[blockMesh] blockMesh pipe problem xiyuqiu OpenFOAM Meshing & Mesh Conversion 6 July 22, 2014 23:50
Double Walled Pipe Boundary dahvqaz FLUENT 2 December 5, 2012 10:14
[GAMBIT] Meshing a pipe vedravi ANSYS Meshing & Geometry 1 March 25, 2010 13:19


All times are GMT -4. The time now is 02:16.