|
[Sponsors] |
July 23, 2010, 08:09 |
dangling edges in refined hexahedral grid
|
#1 |
New Member
Lars
Join Date: Jul 2010
Posts: 1
Rep Power: 0 |
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: |
|
July 23, 2010, 08:41 |
|
#2 |
Senior Member
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 403
Rep Power: 25 |
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. |
|
July 24, 2010, 09:53 |
|
#3 |
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,905
Rep Power: 33 |
Heya,
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) Enjoy, Hrv
__________________
Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
On grid independent solution for pulsatile flow | David | FLUENT | 5 | March 25, 2022 03:33 |
[snappyHexMesh] SnappyHexMesh for internal Flow | vishwa | OpenFOAM Meshing & Mesh Conversion | 24 | June 27, 2016 08:54 |
[snappyHexMesh] external flow with snappyHexMesh | chelvistero | OpenFOAM Meshing & Mesh Conversion | 11 | January 15, 2010 19:43 |
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 |