|
[Sponsors] |
[General] Repair holes in surface (from haning nodes) |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 21, 2014, 12:00 |
Repair holes in surface (from haning nodes)
|
#1 |
New Member
Adam Læssøe
Join Date: Nov 2013
Posts: 29
Rep Power: 12 |
I have hexa-mesh produced with ICEM. The mesh is block-refined in ICEM, so there are hanging nodes in places. It seems Paraview cannot handle hanging nodes, but is there some filter that will either
a) do as the answer in the link suggests - break hanging nodes into tetras, and interpolate values or b) simply fill the holes? The data there isn't important, but it doesn't look very nice. |
|
February 21, 2014, 16:50 |
|
#2 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Greetings Adam,
Can you provide a simple example file? It would make it easier to diagnose and test how this can be achieved. Best regards, Bruno |
|
February 21, 2014, 17:38 |
|
#3 |
New Member
Adam Læssøe
Join Date: Nov 2013
Posts: 29
Rep Power: 12 |
Sure - this is the smallest working example(18MB) I could get from my data.
To reproduce: Import with EnSight CellDataToPointData Contour[phase_2_vof,0.5] |
|
February 22, 2014, 14:13 |
|
#4 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Hi Adam,
I noticed that you've also asked about this at the mailing list: http://www.paraview.org/pipermail/pa...ry/030632.html I believe the trick is to use the filter "Tetrahedralize" Keep in mind that this will substantially increase the computational requirements of the resulting mesh (680 kcell went up to 3.9 Mcell). Attached is shown the results of my investigation into this:
edit: Sorry, I didn't properly read your question today You had already tried this and it was too heavy, as you indicated. OK, try using the filter "Extract cells by region", which can allow you to get only a certain section of the mesh. Then apply the "Tetrahedralize" filter to it. You can use multiple "Extract cells by region" filters with the "box" sub-filter (instead of "plane") and then use the "Group Datasets" filter or one of the "Append" filters to group all extractions together. Best regards, Bruno Last edited by wyldckat; February 22, 2014 at 14:19. Reason: see "edit:" |
|
February 22, 2014, 14:35 |
|
#5 |
New Member
Adam Læssøe
Join Date: Nov 2013
Posts: 29
Rep Power: 12 |
Hi Bruno
Thanks a lot for looking into it! Yeah, I wasn't sure how much attention was given here so after a day of waiting, I figured I'd try the mailing list As you say, the Tetrahedralize increases computation a lot. For animating transient data, that might well be a deal breaker. On top of that, the output seems more faceted. Could there be a way to only tetrahedralize the cells with hanging nodes? For that to work I suppose something like the following would be necessary (with [ ] meaning taken input) : Full data set (1); |-Extract cells w/ hanging nodes (2); | |-Tetrahedralize (3); | |-'Mask' using data set [(2)] (4); |-Merge with dataset [(3)] (5); But I have no idea how to accomplish (2) - let alone (4) and (5). /adam |
|
February 22, 2014, 14:46 |
|
#6 |
New Member
Adam Læssøe
Join Date: Nov 2013
Posts: 29
Rep Power: 12 |
Oh, and I tried VeronoiMergeConnected, thinking it sounded maybe applicable - but ParaView crashes when I apply it.
|
|
February 22, 2014, 15:23 |
|
#7 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
The steps list you wrote gave me an idea and attached is the resulting pipeline. I'm not sure how much lighter it truly is and the results seem a bit... strange.
The pipeline is essentially:
edit: never mind, the "Pass Arrays" filter seems to be broken... at least in another test I did. edit 2: The third image shows another possible pipeline, which may or may not be a little bit lighter. Last edited by wyldckat; February 22, 2014 at 15:29. Reason: see "edit:" and "edit 2" |
|
February 22, 2014, 16:47 |
|
#8 |
New Member
Adam Læssøe
Join Date: Nov 2013
Posts: 29
Rep Power: 12 |
That's _exactly_ what I had in mind.
It does look a little weird still, but I suspect no more so than what all the data would have looked like if tetrahedralizing everything. As for memory, the Information tab reports the full Tetrahedralize as 290MB, while the MergeBlocks is reported as 140 - but then the mesh quality reported 120MB - which sounds weird, so I'm wondering how accurate that panel really is. For some reason, I couldn't select AppendDatasets - 'twas grayed out. |
|
March 1, 2014, 05:03 |
|
#9 | |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Hi Adam,
Quote:
Best regards, Bruno |
||
January 10, 2019, 04:02 |
|
#10 | |
New Member
Join Date: Feb 2013
Posts: 21
Rep Power: 13 |
Dear Bruno,
I tried to follow your instructions to create a decent looking visualisation of my domain, which contains hanging nodes (see attachment 1). Although the procedure you mentioned helps to improve the iso-surface quality, it is still far from decent representation and even produces false representation (see attachments 2 where the proposed fix leads to showing a continuous surface, when it shouldn't according to screenshot 3). Do you know of any workaround, since the solution was proposed in 2014? Quote:
Last edited by Shahr; January 10, 2019 at 04:10. Reason: changing one of the screenshots |
||
January 12, 2019, 08:56 |
|
#11 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Quick answer @Shahr: I'm a bit confused on what I should be seeing in the first attached image, because all I can understand is that this is a surface mesh of the actual meshed domain, it's not an iso-surface. Therefore, there should not be any hanging nodes on that mesh... unless I'm not understanding it correctly.
Beyond that, I believe that the more recent versions of ParaView have improved their capabilities for iso-surface generation, so try using the more advanced options that are revealed when you turn on the gear icon. |
|
January 12, 2019, 09:22 |
|
#12 |
New Member
Join Date: Feb 2013
Posts: 21
Rep Power: 13 |
Dear Bruno,
thanks for your reply. You are right. The first snapshot is just the domain in opaque solid colour with "surface with edges" option selected. The inner conical structure visible in that image is the irregular shapes that are visible even in the mesh visualisation without any field selected. I attached another shot of another region of my domain which has local refinement, that leads to hanging nodes, without showing cell edges. So the problem is not only specific to iso-surface generation, but in general for the mesh representation. Of course this issue leads to messed up representation of any iso-surface/iso-volume that intersects those blocks. In order to reproduce this problem, just locally refine/coarsen any block in ICEM (see attachment 2) and import the mesh in OF using the generated fluent mesh. |
|
January 13, 2019, 13:33 |
|
#13 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Quick answer: OK, so the mesh was converted to OpenFOAM format... in that case, please try to use the vtk polyhedron representation, instead of the default option, as instructed here: https://openfoamwiki.net/index.php/F...is_in_ParaView - because it looks like that the mesh you are seeing has the more complex cells decomposed into tetrahedral cells... and the problem with that is that it may interpolate values incorrectly into the decomposed mesh...
|
|
January 16, 2019, 03:22 |
|
#14 |
New Member
Join Date: Feb 2013
Posts: 21
Rep Power: 13 |
Thank you Mr. Santos very much. Your suggested solution fixed our problem.
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[snappyHexMesh] Surface triangulation using snappyHexMesh | shaileshbg | OpenFOAM Meshing & Mesh Conversion | 4 | October 17, 2019 04:42 |
How To Create Fowcs William Hawkings Surface in Fluent? | marinkobezu | FLUENT | 7 | October 25, 2018 10:51 |
Pierced Faces, Meshing, Repair Surface | O.Lilienthal | STAR-CCM+ | 1 | July 7, 2013 04:37 |
To make sure nodes on surface | Twiti | CFX | 0 | December 23, 2004 20:06 |
CFX4.3 -build analysis form | Chie Min | CFX | 5 | July 12, 2001 23:19 |