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

[blockMesh] pipe bend mesh imperfect

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By jochemvandenbosch

 
 
LinkBack Thread Tools Search this Thread Display Modes
Prev Previous Post   Next Post Next
Old   December 16, 2011, 10:13
Default pipe bend mesh imperfect
  #1
New Member
 
Jochem van den Bosch
Join Date: Oct 2011
Location: Rotterdam, The Netherlands
Posts: 20
Rep Power: 14
jochemvandenbosch is on a distinguished road
Dear fellow Foamers,

I'm still sort of a newbie in OpenFOAM and I'm trying to mesh a pipe bend in blockMesh. I've got no errors when I run blockMesh (and checkMesh) but for some reason the diameter of the bend seems smaller than the diameter of the rest of the pipe. I've checked every node and block like 5 times and everything seems to be in order. Does anybody see what mistake(s) I made? Thanks in advance! This is my blockMeshDict (in the attachment a picture of my mesh in paraFoam)

Quote:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.3 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;

root "";
case "";
instance "";
local "";

class dictionary;
object blockMeshDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 1;

vertices
(
// pipe horizontal
(0 0 0) // 0
(1 0 0) // 1
(1 0.12 0) // 2
(0 0.12 0) // 3
(0 0 0.12) // 4
(1 0 0.12) // 5
(1 0.12 0.12) // 6
(0 0.12 0.12) // 7

// bend region
(1.34 0 -0.34) // 8
(1.46 0 -0.34) // 9
(1.46 0.12 -0.34) // 10
(1.34 0.12 -0.34) // 11

// pipe vertical
(1.34 0 -1.34) // 12
(1.46 0 -1.34) // 13
(1.46 0.12 -1.34) // 14
(1.34 0.12 -1.34) // 15

// cylinder nodes
(0 -0.01071 -0.01071) // 16
(0 0.13071 -0.01071) // 17
(0 0.13071 0.13071) // 18
(0 -0.01071 0.13071) // 19

(1 -0.0107 -0.01071) // 20
(1 0.13071 -0.01071) // 21
(1 0.13071 0.13071) // 22
(1 -0.01071 0.13071) // 23

(1.32929 -0.01071 -0.34) // 24
(1.47071 -0.01071 -0.34) // 25
(1.47071 0.13071 -0.34) // 26
(1.32929 0.13071 -0.34) // 27

(1.32929 -0.01071 -1.34) // 28
(1.47071 -0.01071 -1.34) // 29
(1.47071 0.13071 -1.34) // 30
(1.32929 0.13071 -1.34) // 31
);

blocks
(
// inlet block
hex (0 1 2 3 4 5 6 7) (10 10 10) simpleGrading (1 1 1)
hex (4 5 6 7 19 23 22 18) (10 10 10) simpleGrading (1 1 1)
hex (3 2 21 17 7 6 22 18) (10 10 10) simpleGrading (1 1 1)
hex (16 20 21 17 0 1 2 3) (10 10 10) simpleGrading (1 1 1)
hex (16 20 1 0 19 23 5 4) (10 10 10) simpleGrading (1 1 1)

// bend block
hex (1 8 11 2 5 9 10 6) (10 10 10) simpleGrading (1 1 1)
hex (5 9 10 6 23 25 26 22) (10 10 10) simpleGrading (1 1 1)
hex (2 11 27 21 6 10 26 22) (10 10 10) simpleGrading (1 1 1)
hex (20 24 27 21 1 8 11 2) (10 10 10) simpleGrading (1 1 1)
hex (20 24 8 1 23 25 9 5) (10 10 10) simpleGrading (1 1 1)

// outlet block
hex (12 13 14 15 8 9 10 11) (10 10 10) simpleGrading (1 1 1)
hex (13 29 30 14 9 25 26 10) (10 10 10) simpleGrading (1 1 1)
hex (15 14 30 31 11 10 26 27) (10 10 10) simpleGrading (1 1 1)
hex (28 12 15 31 24 8 11 27) (10 10 10) simpleGrading (1 1 1)
hex (28 29 13 12 24 25 9 8) (10 10 10) simpleGrading (1 1 1)

);

edges
(
// cylinder 1
arc 16 17 (0 0.06 -0.04)
arc 17 18 (0 0.16 0.04)
arc 18 19 (0 0.06 0.16)
arc 19 16 (0 -0.04 0.06)
// cylinder 2
arc 20 21 (1 0.06 -0.04)
arc 21 22 (1 0.16 0.04)
arc 22 23 (1 0.06 0.16)
arc 23 20 (1 -0.04 0.06)
// cylinder 3
arc 24 25 (1.4 -0.04 -0.34)
arc 25 26 (1.5 0.06 -0.34)
arc 26 27 (1.4 0.16 -0.34)
arc 27 24 (1.3 0.06 -0.34)
// cylinder 4
arc 28 29 (1.4 -0.04 -1.34)
arc 29 30 (1.5 0.06 -1.34)
arc 30 31 (1.4 0.16 -1.34)
arc 31 28 (1.3 0.06 -1.34)
// bend
arc 1 8 (1.24041 0 -0.09958)
arc 5 9 (1.32527 0 -0.01473)
arc 2 11 (1.24041 0.12 -0.09958)
arc 6 10 (1.32527 0.12 -0.01473)

arc 20 24 (1.23284 -0.01071 -0.10716)
arc 21 27 (1.23284 0.13071 -0.10716)
arc 23 25 (1.33284 -0.01071 -0.00715)
arc 22 26 (1.33284 0.13071 -0.00715)
);


patches
(

patch inlet
(
(4 7 3 0)
(19 18 7 4)
(7 18 17 3)
(0 3 17 16)
(19 4 0 16)
)
patch outlet
(
(12 15 14 13)
(28 31 15 12)
(15 31 30 14)
(13 14 30 29)
(28 12 13 29)
)
wall walls
(
(19 23 22 18)
(17 18 22 21)
(16 17 21 20)
(23 19 16 20)
(20 21 27 24)
(21 22 26 27)
(23 25 26 22)
(23 20 24 25)
(24 27 31 28)
(27 26 30 31)
(26 25 29 30)
(25 24 28 29)
)
);

mergePatchPairs
(
);

// ************************************************** *********************** //
Attached Images
File Type: jpg pipe bend.jpg (24.2 KB, 282 views)
gchamma likes this.
jochemvandenbosch is offline   Reply With Quote

 


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
[mesh manipulation] extrudeMesh for pipe bend cfd4res OpenFOAM Meshing & Mesh Conversion 0 March 15, 2018 13:04
[GAMBIT] 3D Two pipe T-shape junction mesh Jing ANSYS Meshing & Geometry 8 April 26, 2013 12:27
[snappyHexMesh] snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Meshing & Mesh Conversion 2 March 27, 2011 21:11
[blockMesh] generate a graded pipe mesh. jenright OpenFOAM Meshing & Mesh Conversion 0 August 22, 2009 08:58
Icemcfd 11: Loss of mesh from surface mesh option? Joe CFX 2 March 26, 2007 18:10


All times are GMT -4. The time now is 17:06.