CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > ANSYS Meshing & Geometry

[ICEM] Negative volume elements

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 3 Post By bluebase

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 11, 2020, 04:52
Default Negative volume elements
  #1
New Member
 
Hazem
Join Date: Mar 2020
Posts: 11
Rep Power: 6
Hazem24 is on a distinguished road
Dear Members,

I'd been working on a wind turbine model (3 blades connected at a hub, all surrounded in an inner domain cylinder and an outer domain).

I thought I finally captured the geometry really well (pictures below) with blocks. Turns out there are some negative volume elements (negative determinants) in some areas of the blocks mainly in the connection regions between the inner and outer domain cylinders not at the turbine itself.
I tried PREMESH smoothing which took long hours and actually reduced them but still around 3000 negative determinants remained (out of 24mil) and that was enough to cause "divergence in the AMG solver" while initializing in fluent.

I then tried smoothing and repairing of MESH based on determinant criterion up to 0.2 and 0.0 but it also took almost a day for each and didn't remove any negative volumes.

I wonder if someone has a clue, many thanks in advance! If any other details or snapshot is needed please inform me.

Hazem.
happy massenger.JPG

happy massenger 2.JPG

quality 3x3 determinant (Mess2).jpg

quality angle (Mess2).jpg
Hazem24 is offline   Reply With Quote

Old   May 11, 2020, 06:56
Default
  #2
Senior Member
 
Sebastian Engel
Join Date: Jun 2011
Location: Germany
Posts: 566
Rep Power: 20
bluebase will become famous soon enough
Hi Hazem,

when you have negative elements still present, there is no point to run any solver with that mesh, almost any will fail.
The same is true for smoothing, whether the premesh or mesh smoother.

From the few image, your blocking looks quite nice and advanced. With such complexity of the blocking you might have to tweak some edges and vertices.

In regions where blocks stretch out and follow some curvature in neighbouring blocks the automatly created premesh can squeeze elements into each other till they become negative.

Using the qualitiy histogram is already a good head-start! A small tip is to change the plot boundaries to -1 and 1 and plot an even number of columns, so, you can draw the only the negative elements and find the source of the issue.

The next feature you likely need to resolve those areas is called "scan planes". Its a feature which you find in the context menu of the premesh object ( in the model tree)
It will allow you to draw layer of the premesh in specific directions.

In your case, i would draw a scan plan perpendicular to the "plane" of negative elements, and check which edge is deformed, or which vertex needs to be moved.

At the last step, i would play with "split edge" or "link edge" to prescribe a specifc curvature to edges to overcome the issues in the default premesh computation.

In conclusion, your blocking seems fixable to me. The solution is to find the entities which need individual manipulation.

Best,
Sebastian
redpsi, Sai Krishna and aero_head like this.
bluebase is offline   Reply With Quote

Old   May 11, 2020, 15:40
Default
  #3
Member
 
mCiFlDk's Avatar
 
mCiFlDk
Join Date: Feb 2020
Posts: 56
Rep Power: 6
mCiFlDk is on a distinguished road
Quote:
Originally Posted by Hazem24 View Post
I tried PREMESH smoothing which took long hours and actually reduced them but still around 3000 negative determinants remained (out of 24mil) and that was enough to cause "divergence in the AMG solver" while initializing in fluent.

I then tried smoothing and repairing of MESH based on determinant criterion up to 0.2 and 0.0 but it also took almost a day for each and didn't remove any negative volumes.

I wonder if someone has a clue, many thanks in advance! If any other details or snapshot is needed please inform me.

Hazem.
Hi Hazem.

I'm not a super expert in ICEM, that's why I can offer another idea: using another ANSYS pre-processor. What I use for complex 3D geometries is Fluent Meshing, it has a very good semi-manual way to fix cells with a wide variety of problems. I frequently obtain many errors in terms of free faces, disjointed cells, negative elements... And thanks to it I'm able to fix it 99% of time.

As a tip I wouldn't recommend to use global smoothers or similar tools, because they distort many cells that don't need to be fixed.

I hope it helped you a bit. Best of luck.
mCiFlDk is offline   Reply With Quote

Old   May 17, 2020, 06:10
Default
  #4
New Member
 
Hazem
Join Date: Mar 2020
Posts: 11
Rep Power: 6
Hazem24 is on a distinguished road
Dear Sebastian,

I'm very thankful for your answer !
I've also sent you a private message, could you please see that one too ?
Hazem24 is offline   Reply With Quote

Old   September 22, 2020, 14:43
Default Same negative volume elements
  #5
Member
 
PENANG
Join Date: Aug 2017
Location: Malaysia
Posts: 40
Rep Power: 8
musa19 is on a distinguished road
Hi members, I hope you are doing good. I am facing the same problem of negative elements. I tried several ways to fix but unfortunately the elements are not fixing. I request to please help me fix this problem. Thanks in advance. I am hereby attaching some pictures of the mesh. Thank you.
Attached Images
File Type: png pic1.PNG (8.0 KB, 32 views)
File Type: png pic2.PNG (7.7 KB, 27 views)
File Type: png pic3.PNG (16.5 KB, 22 views)
File Type: png pic4.PNG (20.0 KB, 21 views)
File Type: png pic6.PNG (88.4 KB, 31 views)
musa19 is offline   Reply With Quote

Old   September 22, 2020, 15:37
Default
  #6
Senior Member
 
Sebastian Engel
Join Date: Jun 2011
Location: Germany
Posts: 566
Rep Power: 20
bluebase will become famous soon enough
It seems the respective curve is not well defined.


You need to group all of these curves and associate them to the respective edge.
It seems likely to me, that the unwell defined edge has some sections which you either missed, or weren't obviously selectable. I guess, that's where the premesh's volume mesh detached and tried to get closer to the block's edge, which is straight.
bluebase is offline   Reply With Quote

Old   September 22, 2020, 15:52
Default
  #7
Member
 
PENANG
Join Date: Aug 2017
Location: Malaysia
Posts: 40
Rep Power: 8
musa19 is on a distinguished road
Hi Sebastian, Thank you for your quick reply. I did change the curves but it is messing up with the new curve association. I tried few things but I am not able to overcome this. If you can look into this case and update me about the issue. I will be thankful.
I can share the file with you and please let me know if it is ok for you. Thank you very much.
musa19 is offline   Reply With Quote

Old   September 22, 2020, 16:27
Default
  #8
Senior Member
 
Sebastian Engel
Join Date: Jun 2011
Location: Germany
Posts: 566
Rep Power: 20
bluebase will become famous soon enough
Hi Penang,


i am sorry, but i don't work on privately send material here. It contradicts the purpose of this forum in my opinion.


You may share what you can openly.


Best Sebastian




As a side note: I do offer my service as freelancer, in case you need it.
bluebase is offline   Reply With Quote

Old   September 23, 2020, 12:23
Default
  #9
Member
 
PENANG
Join Date: Aug 2017
Location: Malaysia
Posts: 40
Rep Power: 8
musa19 is on a distinguished road
Thank you, Sebastian, I understand your point. I have no issues attaching the file here. I am hereby attaching the link. Can you look into this and update me about the issue. Thank you.

https://drive.google.com/file/d/17iL...ew?usp=sharing
musa19 is offline   Reply With Quote

Old   November 15, 2020, 06:20
Default Dear Sebastian
  #10
New Member
 
Hazem
Join Date: Mar 2020
Posts: 11
Rep Power: 6
Hazem24 is on a distinguished road
I'm working on a multizone mesh for a 3 bladed wind turbine rotor on ICEM CFD, by doing 2D surface blocking then converting the block to 3D and creating an Ogrid.
This essentially requires creating topology first (which was tricky in the beginning to get all the curves red).


The 2D blocking went fine, then on converting to 3D the Ogrid failed!

This is possibly because a surface was meshing on the face of the airfoil at the blade tip (I've included pictures).
I started all over, this time trying to create surface using ICEM Geometry tools (surface from 2-4 curves). which resulted in 2 blue curves at the tip, trying to vary the tolerance results in other blue and yellow curves elsewhere!
I wonder if I can get any clues on how to create this surface and build topology.


The model is found in the following drive link
https://drive.google.com/drive/folde...i9?usp=sharing

Best Regards.
Attached Images
File Type: jpg 1.JPG (72.1 KB, 16 views)
File Type: jpg 2.JPG (106.9 KB, 14 views)
File Type: jpg 3.JPG (93.7 KB, 11 views)
File Type: jpg 4.JPG (39.8 KB, 9 views)
Hazem24 is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How to fix negative volume elements in icem cfd? shirish CFX 5 June 5, 2018 01:59
[blockMesh] blockMesh error - Negative Volume Block adoledin OpenFOAM Meshing & Mesh Conversion 2 June 22, 2016 10:44
[ICEM] Negative volume error in hybrid mesh siw ANSYS Meshing & Geometry 4 September 3, 2014 05:25
channelFoam for a 3D pipe AlmostSurelyRob OpenFOAM 3 June 24, 2011 13:06
[ICEM] Negative Volumes & 7 node HEXA elements mr_stoked ANSYS Meshing & Geometry 1 September 21, 2010 12:45


All times are GMT -4. The time now is 15:53.