CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Bugs

Floating point error in BlockMesh for really simple mesh

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 8, 2009, 05:47
Default Floating point error in BlockMesh for really simple mesh
  #1
New Member
 
Coen Wit
Join Date: Jul 2009
Posts: 5
Rep Power: 8
CoenW is on a distinguished road
As a first trial with OpenFOAM from the ground up I am creating a really simple simulation: 2D flow across a cylinder. So far I've only edited the blockMeshDict to create this geometry. I've tried two blocking strategies on this problem (the first one is a little weird):


Instead of meshing this geometry correctly blockMesh gives the following error (I've left out the header) :
Quote:
Creating curved edges

Creating blocks
Floating point exception
and crashes back to the terminal.

To determine what the cause of this problem is I've meshed the lower half of the problem, which works without any problem. The same goes when I mesh the complete problem, but without the blocks and faces on the upper half of the cylinder (see attached figures).


With the other (normal) blocking strategy, I can't even get half the mesh to work correctly.
Changing the order of the points in the blocking hasn't helped me so far and I'm lost for a solution of this problem. I think it has something to do with the use of arcs, but I'm not sure.
Any solutions or ideas would be greatly appreciated. I'd like to get the whole case to run in OpenFOAM, instead of having to resort to a commercial program for meshing.

I've attached the following files:
blockMeshDict.fail <-- the complete mesh, but with the top blocks near the cylinder commented out
blockMeshDicthalf.work <-- the lower half of the problem, which meshes correctly
blockMeshDict <-- dictionary for a more correct blocking strategy, which also gives the floating point error
Attached Files
File Type: txt blockMeshDict.txt (4.3 KB, 21 views)
File Type: txt blockMeshDicthalf.work.txt (5.6 KB, 8 views)
File Type: txt blockMeshDict.fail.txt (5.9 KB, 5 views)
CoenW is offline   Reply With Quote

Old   July 9, 2009, 04:47
Default
  #2
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
I tried blockMeshDict.txt in 1.5.x and that gives

face 0 in patch 0 does not have neighbour cell face: 4(12 7 27 32).

You can visualise the blocks and points of your blockMesh with the -blockTopology option to dump out an .obj file which you can postprocess with e.g. javaview or with Paraview (using objToVTK).

The blockMeshDict.fail seems to work ok.
mattijs is offline   Reply With Quote

Old   July 9, 2009, 06:56
Default
  #3
New Member
 
Coen Wit
Join Date: Jul 2009
Posts: 5
Rep Power: 8
CoenW is on a distinguished road
Quote:
Originally Posted by mattijs View Post
I tried blockMeshDict.txt in 1.5.x and that gives

face 0 in patch 0 does not have neighbour cell face: 4(12 7 27 32).

You can visualise the blocks and points of your blockMesh with the -blockTopology option to dump out an .obj file which you can postprocess with e.g. javaview or with Paraview (using objToVTK).

The blockMeshDict.fail seems to work ok.
Thanks for having a look at it. Apparently I didn't mention it clearly, but I commented out the sections of the files that cause the floating point error. So that's why it reports the neighbour cell face problem: the face is commented out.

the .fail file creates the mesh depicted in the second image. Again: the offending faces (the missing ones in the middle) have already been commented out. I've attached the uncommented file below. Could you check if it does the same thing on your system?

I didn't know about the -blockTopology option, where does it dump the files? I ran the option while I also got the floating point error and I can't find any .obj files.
Attached Files
File Type: txt blockMeshDict.txt (4.3 KB, 13 views)

Last edited by CoenW; July 9, 2009 at 07:00. Reason: attach file
CoenW is offline   Reply With Quote

Old   July 22, 2009, 06:33
Default
  #4
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
There was a problem with having zero points along the edge (your blocks are 1x1x1). It works with more cells along the edge. I pushed a fix to 1.5.x so it works with 1 cell as well.

Thanks,

Mattijs
mattijs 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
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 12:55
Very simple moving mesh Pete FLUENT 4 February 10, 2006 01:12
Floating point error Chico Mbanu FLUENT 1 July 15, 2004 08:10
Floating point error Chico Mbanu FLUENT 0 July 14, 2004 13:56
unstructured vs. structured grids Frank Muldoon Main CFD Forum 1 January 5, 1999 11:09


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