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

[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 Search this Thread Display Modes
Old   November 5, 2013, 21:19
Default blockMesh gives "Inconsistent point locations" for axisymmetric problem
  #1
New Member
 
Geoffrey Irving
Join Date: Nov 2013
Posts: 3
Rep Power: 12
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 http://www.cfd-online.com/Forums/ope...anslation.html). 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, 19 views)
File Type: txt blockMeshDict.txt (757 Bytes, 8 views)
irving is offline   Reply With Quote

Old   November 5, 2013, 21:24
Default
  #2
New Member
 
Geoffrey Irving
Join Date: Nov 2013
Posts: 3
Rep Power: 12
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, 10:10
Default
  #3
New Member
 
Jason Lassaline
Join Date: Mar 2012
Posts: 3
Rep Power: 14
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 11:14. Reason: Correction to factor.
jvl001 is offline   Reply With Quote

Old   July 8, 2014, 01:11
Default
  #4
New Member
 
miladrakhsha
Join Date: Aug 2012
Posts: 29
Rep Power: 13
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, 03:28
Default Same trouble
  #5
Member
 
Gareth
Join Date: Jun 2010
Posts: 56
Rep Power: 15
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
http://www.cfd-online.com/Forums/ope...use-prism.html
*/
Attached Files
File Type: txt blocks.txt (21.3 KB, 3 views)

Last edited by bullmut; August 7, 2014 at 09:17. Reason: Used a work around to solve
bullmut is offline   Reply With Quote

Old   January 27, 2021, 12:23
Default
  #6
New Member
 
Sophie
Join Date: Jan 2021
Posts: 11
Rep Power: 5
So_LL is on a distinguished road
Many years have passed, I hope you are still there to answer

I would like to modify the source code of blockMeshMerge.C to correct an error similar to OP's reporter error (grading inconsistency), but the file is locked and can only be opened as readonly. Maybe this is a Linux thing, possibly a stupid question, but how can I modify the source code to test your solution?
So_LL is offline   Reply With Quote

Old   January 27, 2021, 13:27
Default Alas, 7 years is too much
  #7
New Member
 
Geoffrey Irving
Join Date: Nov 2013
Posts: 3
Rep Power: 12
irving is on a distinguished road
I've changed jobs 4 times since then, and also jumped from computational physics to AGI safety.
irving is offline   Reply With Quote

Old   January 27, 2021, 13:31
Default
  #8
New Member
 
Sophie
Join Date: Jan 2021
Posts: 11
Rep Power: 5
So_LL is on a distinguished road
Ho well ! Good for you !

I realized I just had to use sudo before opening the file with gedit ! Also found my issue had nothing to do with that, and just solved it 2 minutes ago.

Thank you for taking the time to reply though !

Have a good day !
So_LL is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem accessing data in point boundary condition Jibran OpenFOAM Programming & Development 0 February 28, 2018 13:12
[blockMesh] Problem using C in blockMesh DeAldeEiam OpenFOAM Meshing & Mesh Conversion 3 November 7, 2017 17:30
monitor point problem Omran CFX 12 April 3, 2017 10:35
[blockMesh] Inconsistent point locations between block pair 7 and 8 ordinary OpenFOAM Meshing & Mesh Conversion 1 November 27, 2016 12:07
[snappyHexMesh] snappyHexMesh and cyclic boundaries Ruli OpenFOAM Meshing & Mesh Conversion 2 December 9, 2013 07:51


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