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

SnappyHexMesh hangs during refinement

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 8, 2009, 15:51
Default SnappyHexMesh hangs during refinement
  #1
New Member
 
Simon Rees
Join Date: Mar 2009
Posts: 12
Rep Power: 8
sjrees is on a distinguished road
I am trying to mesh a geometry consisting of a single background mesh box created with blockMesh and a single STL surface that will become the bottom of the box. This is a rough surface with a number of building outlines protruding into the flow.

I can get the mesh to refine to the first level but snappyHexMesh always hangs when trying to refine to the second level. I can an error:

>Surface refinement iteration 2
>------------------------------
>
>[5] meshRefinement::checkData() : Found 0 duplicate pairs of faces.
>[7] meshRefinement::checkData() : Found 0 duplicate pairs of faces.
>Marked for refinement due to surface intersection : 9 cells.

The application is running flat out continually (100% cpu shown using top) from this point. Only a modest amount of memory is being used so I know it is not swapping at this point.

Any ideas please?
sjrees is offline   Reply With Quote

Old   September 11, 2010, 10:43
Default
  #2
Senior Member
 
Ziad Boutanios
Join Date: Mar 2009
Location: Montréal, Canada
Posts: 113
Rep Power: 8
ziad is on a distinguished road
Same problem here.

Meshing worked fine before then I shifted the outlet plane further downstream (flow around a 3D rectangular building structure) and sHM just hangs after the 4th mesh refinement iteration 2. This is the last "mesh refinement iteration 2" since I am running 4 threads in parallel. The other 3 complete no problem. Shifting the outlet domain further downstream merely adds a rectangular chunk of structured hex on the blockMesh side. sHM doesn't even do anything down there.

Tried refining/coarsening around the point where it was hanging but no luck.

Output from debug level 1 is attached and shows where sHM is hanging.
Attached Files
File Type: gz output.log.gz (5.3 KB, 2 views)
ziad is offline   Reply With Quote

Old   September 11, 2010, 12:29
Default obviously a snappyHexMesh bug
  #3
Senior Member
 
Ziad Boutanios
Join Date: Mar 2009
Location: Montréal, Canada
Posts: 113
Rep Power: 8
ziad is on a distinguished road
As posted previously, refining or coarsening the area of the problem did not help. Also tried playing with the mesh quality parameters such as allowing more or less concavity, using more smoothing iterations, etc. No luck there either.

Here's what worked: pushing the outlet plane of the domain a few more units downstream! Weird but consistent with the source of the problem. sHM started hanging when I initially pushed the outlet plane downstream.

Definitely a (rather exotic) bug!
ziad is offline   Reply With Quote

Old   September 11, 2010, 17:05
Default
  #4
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,253
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings Ziad,

Could you report your findings about this bug here: http://www.openfoam.com/bugs/

By the way, are you using the latest OpenFOAM 1.7.x version? If not, perhaps this bug might already have been fixed.

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   September 11, 2010, 21:59
Default
  #5
Senior Member
 
Ziad Boutanios
Join Date: Mar 2009
Location: Montréal, Canada
Posts: 113
Rep Power: 8
ziad is on a distinguished road
Hi Bruno,

I'll report the bug.

Meanwhile I am using 1.7.0. Not planning to upgrade just for this though.

Cheers,
Ziad
ziad 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
SnappyHexMesh - How to set-up to avoid skewfaces? pcaron OpenFOAM Native Meshers: snappyHexMesh and Others 17 September 13, 2014 15:58
SnappyHexMesh with local refinement of ONE STLfile wolle1982 OpenFOAM Mesh Utilities 16 August 23, 2013 12:12
snappyHexMesh "killed" Bob OpenFOAM 4 December 1, 2011 17:48
snappyHexMesh - Floating point error derjames OpenFOAM Mesh Utilities 4 September 19, 2009 14:14
SnappyHexMesh hangs wile mesh refinement bastil OpenFOAM Bugs 3 November 4, 2008 07:08


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