CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Native Meshers: blockMesh

blockMesh gives "Inconsistent point locations" for axisymmetric problem

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

Like Tree1Likes
  • 1 Post By jvl001

Reply
 
LinkBack Thread Tools Display Modes
Old   November 5, 2013, 20:19
Default blockMesh gives "Inconsistent point locations" for axisymmetric problem
  #1
New Member
 
Geoffrey Irving
Join Date: Nov 2013
Posts: 2
Rep Power: 0
irving is on a distinguished road
blockMesh is bailing when I give it two ungraded prisms. The output from blockMesh and the blockMeshDict are attached (I had to add a .txt extension to the latter to get it to upload).

The unusual bit of blockMesh output is me adding print statements to blockMesh to check whether it was a precision issue (based on Inconsistent Grading Caused by Simple Translation). It doesn't appear to be precision related. It's also clearly not grading, since I'm asking for one cell in all dimensions.

Thanks in advance for any suggestions!
Attached Files
File Type: txt blockMesh.txt (2.8 KB, 4 views)
File Type: txt blockMeshDict.txt (757 Bytes, 2 views)
irving is offline   Reply With Quote

Old   November 5, 2013, 20:24
Default
  #2
New Member
 
Geoffrey Irving
Join Date: Nov 2013
Posts: 2
Rep Power: 0
irving is on a distinguished road
Here are excepts of both files for convenience, since looking at attachments seems bothersome.

Code:
# blockMeshDict exerpt
vertices (
  (0 0 0)
  (0 1.375 -0.0119991)
  (0 1.375 0.0119991)
  (3 0 0)
  (3 1.375 -0.0119991)
  (3 1.375 0.0119991)
  (6 0 0)
  (6 1.375 -0.0119991)
  (6 1.375 0.0119991)
);
blocks (
  prism (0 1 2 3 4 5) (1 1 1) simpleGrading (1 1 1)
  prism (3 4 5 6 7 8) (1 1 1) simpleGrading (1 1 1)
);
Code:
# blockMesh output
--> FOAM FATAL ERROR:
Inconsistent point locations between block pair 0 and 1, point 0
    probably due to inconsistent grading.

    From function blockMesh::calcMergeInfo()
    in file blockMesh/blockMeshMerge.C at line 313.
irving is offline   Reply With Quote

Old   November 11, 2013, 09:10
Default
  #3
New Member
 
Jason Lassaline
Join Date: Mar 2012
Posts: 3
Rep Power: 4
jvl001 is on a distinguished road
I've looked into this problem as I was experiencing similar issues. The problem appears to be with the merging of points in
src/mesh/blockMesh/blockMesh/blockMeshMerge.C. I found that not all points along the sharp edge of the wedge were being merged correctly, perhaps due to a truncation issue.

I increased mergeSqrDist at line 121 in Foam::blockMesh::calcMergeInfo by a factor of 100 to guarantee all points along this edge were merged correctly
const scalar mergeSqrDist = magSqr(10*SMALL*bb.span());
and now all wedge-shaped blocks are handled correctly. A factor of SMALL seems to be just a bit too small.
miladrakhsha likes this.

Last edited by jvl001; November 11, 2013 at 10:14. Reason: Correction to factor.
jvl001 is offline   Reply With Quote

Old   July 8, 2014, 00:11
Default
  #4
New Member
 
miladrakhsha
Join Date: Aug 2012
Posts: 26
Rep Power: 4
miladrakhsha is on a distinguished road
Quote:
Originally Posted by jvl001 View Post
I've looked into this problem as I was experiencing similar issues. The problem appears to be with the merging of points in
src/mesh/blockMesh/blockMesh/blockMeshMerge.C. I found that not all points along the sharp edge of the wedge were being merged correctly, perhaps due to a truncation issue.

I increased mergeSqrDist at line 121 in Foam::blockMesh::calcMergeInfo by a factor of 100 to guarantee all points along this edge were merged correctly
const scalar mergeSqrDist = magSqr(10*SMALL*bb.span());
and now all wedge-shaped blocks are handled correctly. A factor of SMALL seems to be just a bit too small.
Dear jvl001
Thank you so much for your invaluable comment. I literally wasted 2 days on this issue
Best
miladrakhsha is offline   Reply With Quote

Old   August 7, 2014, 02:28
Default Same trouble
  #5
New Member
 
Gareth
Join Date: Jun 2010
Posts: 12
Rep Power: 6
bullmut is on a distinguished road
Hi guys

I am having a similar issue

"Creating merge list

--> 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
"

I tried, blindly i might added, adjusting the blockMEshMerge.C file as indicated above but it didnt change the results.

This in an exert from my blockMeshDict:

blocks
(
hex (0 6 60 1 30 36 64 31) blade (10 10 1) simpleGrading (1 1 1)
prism (6 7 60 36 37 64) blade (10 10 1) simpleGrading (1 1 1)
hex (7 8 61 60 37 38 65 64) blade (10 10 1) simpleGrading (1 1 1)
hex (8 9 62 61 38 39 66 65) blade (10 10 1) simpleGrading (1 1 1)
prism (9 10 62 39 40 66) blade (10 10 1) simpleGrading (1 1 1)
hex (62 10 2 1 66 40 32 31) blade (10 10 1) simpleGrading (1 1 1)
prism (63 62 1 67 66 31) blade (10 10 1) simpleGrading (1 1 1)
prism (60 63 1 64 67 31) blade (10 10 1) simpleGrading (1 1 1)
hex (2 10 11 3 32 40 41 33) blade (10 10 1) simpleGrading (1 1 1)
hex (3 11 12 4 33 41 42 34) blade (10 10 1) simpleGrading (1 1 1)
hex (4 12 13 5 34 42 43 35) blade (10 10 1) simpleGrading (1 1 1)
hex (5 13 6 0 35 43 36 30) blade (10 10 1) simpleGrading (1 1 1)
hex (14 22 23 15 44 52 53 45) (10 10 1) simpleGrading (1 1 1)
hex (15 23 24 16 45 53 54 46) (10 10 1) simpleGrading (1 1 1)
hex (16 24 25 17 46 54 55 47) (10 10 1) simpleGrading (1 1 1)
hex (17 25 26 18 47 55 56 48) (10 10 1) simpleGrading (1 1 1)
hex (18 26 27 19 48 56 57 49) (10 10 1) simpleGrading (1 1 1)
hex (19 27 28 20 49 57 58 50) (10 10 1) simpleGrading (1 1 1)
hex (20 28 29 21 50 58 59 51) (10 10 1) simpleGrading (1 1 1)
hex (21 29 22 14 51 59 52 44) (10 10 1) simpleGrading (1 1 1)
);

If anyone can help me out...
I also attached the full file to this reply.
/* No real solution but since i asked somethign similar elsewhere i, thought i should link it
blockMesh-how to use prism
*/
Attached Files
File Type: txt blocks.txt (21.3 KB, 1 views)

Last edited by bullmut; August 7, 2014 at 08:17. Reason: Used a work around to solve
bullmut 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
determining displacement for added points CFDnewbie147 OpenFOAM Native Meshers: snappyHexMesh and Others 1 October 22, 2013 09:53
Problem with blockMesh and my shape TneurolF OpenFOAM Native Meshers: blockMesh 4 June 25, 2013 13:52
a reconstructPar issue immortality OpenFOAM Post-Processing 8 June 16, 2013 11:25
Run blockMesh problem fpetit OpenFOAM Installation 9 January 10, 2008 05:40
CFX4.3 -build analysis form Chie Min CFX 5 July 12, 2001 23:19


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