CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

dangling edges in refined hexahedral grid

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

Like Tree1Likes
  • 1 Post By 31415

LinkBack Thread Tools Display Modes
Old   July 23, 2010, 08:09
Default dangling edges in refined hexahedral grid
New Member
Join Date: Jul 2010
Posts: 1
Rep Power: 0
31415 is on a distinguished road
dear foamers,
i would like to perform local refinement of hexahedral grid in octtree-like structure. my concern is that this approach generates dangling edges (or faces, in 3d) where a refined cell is adjacent to a coarser cell.

how do i handle the situation in the attached pic? to the right of the coarse face F0, two fine faces F1 and F2 have created an extra vertex V3.

if V3 is not a vertex of F0 this mesh is invalid, isn't it? do i have to replace the quadrilateral F0 with the pentagon V0-V1-V2-V3-V4? then the angle V4-V3-V2 is 180. does this impose an accuracy or mesh orthogonality problem?

in general, what is the best approach to locally refine hexahedral meshes? this seems to be standard in gerris (see below), but i have no idea how to work with dangling
edges in openfoam.

example from gerris:
Attached Images
File Type: png tmp.png (6.5 KB, 8 views)
mm.abdollahzadeh likes this.
31415 is offline   Reply With Quote

Old   July 23, 2010, 08:41
Senior Member
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 386
Rep Power: 15
deepsterblue will become famous soon enough
Well, this isn't exactly a 'dangling' entity. One of the benefits to having polyhedral mesh support is the fact that you can pull-off stunts like this quite easily. It looks like you've already grasped the general idea - the 'pentagon' V0-V1-V2-V3-V4 is exactly how it should be.

In general, finite-volume codes don't have the 180-degree restrictions for adjacent cell-faces. What's important is the mesh-orthogonality. In this case, there might be a slight deviation of the cell-to-cell vector away from the face centroid (like the F0-F2 vector for face V3-V4), but non-orthogonal corrections can account for that.

OpenFOAM already has this capability, using dynamicRefineFvMesh.
Sandeep Menon
University of Massachusetts Amherst
deepsterblue is offline   Reply With Quote

Old   July 24, 2010, 09:53
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,768
Rep Power: 21
hjasak will become famous soon enough

If you actually generated a mesh like this (with dangling edges) and it's causing you trouble, there is a way to fix it.

Ages ago I've written a utility to do just that, when we were moving FOAM onto full polyhedral support. It is called zipUpMesh and works without trouble.

By the bye, this is one of the most complex algorithm I've ever written - there is a special award if you can document how it does the job on the wiki (I'm not sure I would be able to write it down myself)


Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting:
hjasak is offline   Reply With Quote


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 for internal Flow vishwa OpenFOAM Native Meshers: snappyHexMesh and Others 23 August 6, 2014 03:50
external flow with snappyHexMesh chelvistero OpenFOAM 11 January 15, 2010 20:43
On grid independent solution for pulsatile flow David FLUENT 4 June 8, 2005 16:20
Combustion Convergence problems Art Stretton Phoenics 5 April 2, 2002 05:59
Grid Independent Solution Chuck Leakeas Main CFD Forum 2 May 26, 2000 11:18

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