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

[Other] Graded mesh of hole in square plate

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 7, 2015, 20:50
Default Graded mesh of hole in square plate
  #1
Member
 
Alex
Join Date: Jan 2014
Posts: 54
Rep Power: 12
H0T_S0UP is on a distinguished road
I have finished the tutorial for stress analysis of a square plate with a hole through the center. I am trying to beyond the tutorial and would like to run a case where the mesh is graded near the hole boundary. I have made a few attempted editing the blockMesh but get fatal errors such as:

Code:
--> FOAM FATAL ERROR: 
Inconsistent point locations between block pair 0 and 1
    probably due to inconsistent grading.

    From function blockMesh::calcMergeInfo()
    in file blockMesh/blockMeshMerge.C at line 294.

FOAM exiting
Quite frankly, my gut feeling tells me that an alternative method from editing the blockmeshdict file will have to be implemented. At the very least the mesh wont be uniform and has to be graded in the normal direction to the hole boundary. I am not aware of a way to do this using edgegrading or simplegrading alone. Any advice is appreciated, thanks.
H0T_S0UP is offline   Reply With Quote

Old   February 8, 2015, 03:46
Default
  #2
Member
 
Peter
Join Date: Feb 2015
Location: New York
Posts: 73
Rep Power: 11
opedrofunk is on a distinguished road
Not sure exactly what you're going for, but you can definitely control the mesh density/gradation at the hole with a few simple edits. You can actually do quite a bit with blockMesh :-)

For example (the original):

Code:
blocks
(
    hex (5 4 9 10 16 15 20 21) (10 10 1) simpleGrading (1 1 1)
    hex (0 1 4 5 11 12 15 16) (10 10 1) simpleGrading (1 1 1)
    hex (1 2 3 4 12 13 14 15) (20 10 1) simpleGrading (1 1 1)
    hex (4 3 6 7 15 14 17 18) (20 20 1) simpleGrading (1 1 1)
    hex (9 4 7 8 20 15 18 19) (10 20 1) simpleGrading (1 1 1)
);
Yields:




With changing the density and gradation:

Code:
blocks
(
    hex (5 4 9 10 16 15 20 21) (20 20 1) simpleGrading (3 1 1)
    hex (0 1 4 5 11 12 15 16) (20 20 1) simpleGrading (3 1 1)
    hex (1 2 3 4 12 13 14 15) (20 20 1) simpleGrading (3 1 1)
    hex (4 3 6 7 15 14 17 18) (20 20 1) simpleGrading (3 1 1)
    hex (9 4 7 8 20 15 18 19) (20 20 1) simpleGrading (1 1 1)
);
Yields:



Cheers,
Peter
opedrofunk is offline   Reply With Quote

Old   February 8, 2015, 16:32
Default
  #3
Member
 
Alex
Join Date: Jan 2014
Posts: 54
Rep Power: 12
H0T_S0UP is on a distinguished road
Thanks! It worked. So, I thought when defining the expansion ratio that it was in absolute x,y,z coordinates. I am guessing since you only changed the first values adjacent to the hole they are in n-t coordinates?

Also, how come you had to grade blocks 4 and 2 (the blocks above and to the right of the hole), instead of just blocks 0 and 1, which comprise the hole?
H0T_S0UP is offline   Reply With Quote

Old   February 18, 2015, 11:05
Default
  #4
New Member
 
Karen Fitzgerald
Join Date: Nov 2014
Location: Dublin, Ireland
Posts: 2
Rep Power: 0
kazzyfitz is on a distinguished road
"Also, how come you had to grade blocks 4 and 2 (the blocks above and to the right of the hole), instead of just blocks 0 and 1, which comprise the hole?"

If I'm interpreting your question properly, adjoining block edges mush have the same number of cells and the same grading. Hope that helps.
kazzyfitz 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
[snappyHexMesh] Creating multiple multiple cell zones with snappyHexMesh - a newbie in deep water! divergence OpenFOAM Meshing & Mesh Conversion 0 January 23, 2019 04:17
[ICEM] Mesh is leaking (hole) at trailing edge of wind turbine blade Shahzad123 ANSYS Meshing & Geometry 5 September 13, 2018 05:52
[blockMesh] Is it possible to have a mesh to represent a domain with a hole kmurphy OpenFOAM Meshing & Mesh Conversion 3 October 13, 2014 17:23
[ICEM] ICEM CFD Cooling Hole Hybrid Mesh RealENG22 ANSYS Meshing & Geometry 1 November 22, 2013 22:27
[snappyHexMesh] snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Meshing & Mesh Conversion 2 March 27, 2011 21:11


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