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

[ICEM] Eliminating Negative Volume Elements

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 14, 2010, 18:15
Default Eliminating Negative Volume Elements
  #1
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 10
Josh is on a distinguished road
Hello all -

Despite that I have used minimal smoothing in ICEM (I just smoothed for orthogonality, being that it is a hexahedral grid), my grid appears to have negative volume elements, as shown below.



Apparently, as shown on the left in blue writing and in the histogram, there are almost 63 negative volume elements (8 + 55), with the smallest being -18.7586 mm^3. I need to get rid of these, but am not sure how.

Any ideas?
Josh is offline   Reply With Quote

Old   April 14, 2010, 18:54
Default
  #2
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 10
Josh is on a distinguished road
The ICEM Help file mentions that it does not consider "degenerate hex elements, and will assign a negative value for them."

Could this be the problem? If so, what is a degenerate hex element?
Josh is offline   Reply With Quote

Old   April 15, 2010, 06:28
Default
  #3
New Member
 
Mark Bell
Join Date: Apr 2010
Location: Holywood
Posts: 12
Rep Power: 8
mbell10 is on a distinguished road
A hexa cell should have six surfaces, ie 8 nodes. A degenerate will have less, for instance, 7 nodes, or even 6 nodes (giving you a prism)
mbell10 is offline   Reply With Quote

Old   April 15, 2010, 14:52
Default
  #4
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 10
Josh is on a distinguished road
Thanks, Mark. I couldn't find an answer to that in the Help file or online. Upon closer inspection, I can see that there are some degenerate elements (e.g. the one shown below), as well as some negative volume elements.

http://img515.imageshack.us/i/t106am...mpt4degen.jpg/

I'm thinking the problem arises from the interaction of my expansive boundary layer on the blade surface and the larger structured H-grid preceding the blade, shown below (the second picture is zoomed-in to the stagnation point).

http://img138.imageshack.us/i/t106ameshcattempt4.jpg/

http://img59.imageshack.us/i/t106ame...mpt4stagn.jpg/

How could I make this interaction feasible? Maybe I should move the blocking vertexes slightly ahead of the stagnation point? Any additional input would be much appreciated.

Last edited by wyldckat; September 3, 2015 at 17:33. Reason: disabled embedded images
Josh is offline   Reply With Quote

Old   April 16, 2010, 05:42
Default
  #5
New Member
 
Mark Bell
Join Date: Apr 2010
Location: Holywood
Posts: 12
Rep Power: 8
mbell10 is on a distinguished road
Hi Josh,

Off the top of my head, I would suggest making a split in the blocks near to the surface of your blade that runs through the inlet and outlet sections as well. This would allow you to manipulate the shape of the block around the leading edge of the blade.

I have attached an image which shows what to do. Hope it helps!
Attached Images
File Type: png fix1.PNG (10.0 KB, 204 views)
mbell10 is offline   Reply With Quote

Old   April 16, 2010, 12:19
Default Suggestions...
  #6
Retired from CFD Online
 
PSYMN's Avatar
 
Simon Pereira
Join Date: Mar 2009
Location: Ann Arbor, MI
Posts: 2,660
Blog Entries: 1
Rep Power: 38
PSYMN has a spectacular aura aboutPSYMN has a spectacular aura about
Excellent suggestion by Mark,

Zooming in on the problem you can see that the mesh is tripping.
Josh_1.jpg
You could do a block split as mark suggested or simply split the edge like this...
Josh_2.jpg
Of course it doesn't help that your periodic boundaries meet the blade below the stagnation point...

you could also try smoothing. This is with the orthogonality smoother, the Multiblock smoother would be much better, but I am out of time for today...
Josh_3.jpg

Probably your best bet would be a different topology that included ogrids around the airfoil blade (Cgrids since you have only half a blade on each side)
A6614_Channel_OGrid_Smoothed_Perioidic.jpg

I will try to get back to you with more later...
PSYMN is offline   Reply With Quote

Old   April 16, 2010, 14:37
Default
  #7
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 10
Josh is on a distinguished road
Hi Mark -

Thanks for the reply.

I assume that to perform that split in your picture, I should use the Split Block > Split Method > Curve Parameter method to allow the block to "follow" my geometry. Below is my blocking without the split.



The picture below shows the blocking with the suggested split (assuming I did it correctly).



Finally, here is the "finished" product blocking scheme.



With this scheme, my pre-mesh looks like this:


Last edited by Josh; April 16, 2010 at 14:58. Reason: Updating
Josh is offline   Reply With Quote

Old   April 16, 2010, 16:19
Default
  #8
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 10
Josh is on a distinguished road
Hi Simon -

Using your method, I split the edge with Control Points located near the vertex. I obtained the following:



One thing I'm having problems with is using the Link Edge tool to make both block edges in the z-direction equal. In the above picture, I did both edges manually, and the overlay is pretty good, but I'd prefer it to be exact. The Link Edge tool, however, does not seem to overlay them properly.

The mesh looks like this:



I think this is the right idea. Obviously, I still have to refine the sizing and smooth with orthogonality.
Josh is offline   Reply With Quote

Old   April 16, 2010, 16:56
Default
  #9
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 10
Josh is on a distinguished road
To follow-up, before I've done any smoothing...

With the Split Edge method, I've obtained negative volume elements along the airfoil suction surface. In the picture below, you can see the mesh "colliding" with the airfoil surface, forming negative elements.

Josh is offline   Reply With Quote

Old   April 17, 2010, 13:31
Default OGrids would avoid all this hassle.
  #10
Retired from CFD Online
 
PSYMN's Avatar
 
Simon Pereira
Join Date: Mar 2009
Location: Ann Arbor, MI
Posts: 2,660
Blog Entries: 1
Rep Power: 38
PSYMN has a spectacular aura aboutPSYMN has a spectacular aura about
Very strange... I will try to take another look some time this weekend...

I can probably block it out with Ogrids in a few minutes and save you some hassle.

Simon
PSYMN is offline   Reply With Quote

Old   April 18, 2010, 22:52
Default Something anyway...
  #11
Retired from CFD Online
 
PSYMN's Avatar
 
Simon Pereira
Join Date: Mar 2009
Location: Ann Arbor, MI
Posts: 2,660
Blog Entries: 1
Rep Power: 38
PSYMN has a spectacular aura aboutPSYMN has a spectacular aura about
Hey Josh... I spend a few minutes with it and created a 2D blocking (less work). You can always extrude this for CFX, or simply write it out as a 2D Fluent file and then read it into CFX (CFX will turn it into a 2.5D model).

I thought I recorded a replay script for you... but then I realized I had forgotten to hit record... Oh well.

Anyway... here are some pics.

First, here is the OGrid Setup (Blocks and Faces to select) for 3D...

Josh_5_OgridSetup.jpg

Then here it is for 2D... (The Yellow blocks are VORFN blocks activated when I selected the edges)

Josh_6_OgridSetup2D.jpg

Then I just used edge splits to get roughly what I wanted, and used edge params to get the distribution that I wanted... You can see that in the blocking file (sent privately).

Then I smoothed it with orthogonality and Laplace (and fudged a few nodes with "move nodes")

here are the final shots... If these won't run we can relax the periodic sides (so they can move but still be periodic) and get a better mesh.

Josh_7_TGridSkew1.jpg

Josh_7_TGridSkew2.jpg

Simon
PSYMN is offline   Reply With Quote

Old   April 19, 2010, 03:50
Default
  #12
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 10
Josh is on a distinguished road
Hi Simon -

You've been an incredible help. Thank you.

One thing, though. You mentioned you sent me (privately) the blocking file, but I never received it. Did you use my Gmail account?
Josh is offline   Reply With Quote

Old   April 20, 2010, 00:45
Default
  #13
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 10
Josh is on a distinguished road
Hi all -

I'm thinking of improving the quality of the post-blade mesh, as a quality check showed problems with the aspect ratio, skew, and overall quality in the angled region just after the blades. You can see the skewing quality in Simon's two bottom pictures in his last post. As shown, the quality could be improved.

I'm thinking there are two ways to go about this:

1) Adjust the blocking (edge distributions, vertex locations) in that region. This seems like the better option, though I'd have to re-mesh and, thus, eliminate all of the fine smoothing and node adjustments Simon performed.

2) Smoothen for quality. I know this isn't a great idea, especially for hexa meshes, and could cause overlap and degenerate/negative elements.

What do you think?

Also, just to follow-up: Simon's mesh ran perfectly in CFX.

Last edited by Josh; April 20, 2010 at 00:46. Reason: added material
Josh is offline   Reply With Quote

Old   April 20, 2010, 13:03
Default Go for it.
  #14
Retired from CFD Online
 
PSYMN's Avatar
 
Simon Pereira
Join Date: Mar 2009
Location: Ann Arbor, MI
Posts: 2,660
Blog Entries: 1
Rep Power: 38
PSYMN has a spectacular aura aboutPSYMN has a spectacular aura about
Sure, go for it, but keep in mind that you need a periodic model... This makes it difficult because you will have to deal with sloping channel behind the airfoil.

Sorry, I haven't had time to do much more than type a quick message. My machine is cranking right now and I just can't run any more tasks...

For smoothing, I just used orthogonal for a bunch of iterations and a few iterations of laplace, then back to orthogonal for 10 more. Then I used move nodes to straighten up some of the nodes ahead of the lower airfoil... Perhaps it would have been a bit better if I had relaxed the orthogonality requirement on the periodic walls...
PSYMN is offline   Reply With Quote

Old   April 20, 2010, 13:04
Default Shifted Periodic
  #15
Retired from CFD Online
 
PSYMN's Avatar
 
Simon Pereira
Join Date: Mar 2009
Location: Ann Arbor, MI
Posts: 2,660
Blog Entries: 1
Rep Power: 38
PSYMN has a spectacular aura aboutPSYMN has a spectacular aura about
You could try a "shifted periodic" blocking...
PSYMN is offline   Reply With Quote

Old   April 20, 2010, 17:14
Default
  #16
Senior Member
 
Joshua Counsil
Join Date: Jul 2009
Location: Halifax, Nova Scotia, Canada
Posts: 366
Rep Power: 10
Josh is on a distinguished road
Thanks, Simon. I'll see how I'm doing for time.
Josh is offline   Reply With Quote

Old   July 9, 2014, 13:01
Default
  #17
Senior Member
 
Join Date: Mar 2010
Location: Cape Town, SA
Posts: 156
Rep Power: 9
Jonathan is on a distinguished road
Hi Josh,

I hope you are well. I cam across your post just the other day when i was looking for blocking strategies for a turbine blade i am modelling. I usually, use the blade-in-the-middle approach with an o-grid, but for specific reasons i need to model it with the blade surfaces at the outer edge, like you did here.

however, having read your posts, i was wondering how / whether, you managed to get a fully periodic (cyclic) mesh out using this approach i.e. did you manage to run your simulations ok using this mesh.

i know it is only 2D, but i assume you needed to set periodic vertices and the periodic rotation angle in the global mesh setup? Would you be willing / able to give me an idea of how you did it?

many thanks and best regards
cheers
jonathan

Quote:
Originally Posted by Josh View Post
Hi Simon -

You've been an incredible help. Thank you.

One thing, though. You mentioned you sent me (privately) the blocking file, but I never received it. Did you use my Gmail account?
Jonathan 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
BlockMesh FOAM warning gaottino OpenFOAM Native Meshers: blockMesh 7 July 19, 2010 14:11
On the damBreak4phaseFine cases paean OpenFOAM Running, Solving & CFD 0 November 14, 2008 22:14
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 12:55
Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Native Meshers: blockMesh 10 April 2, 2007 14:00
CFX4.3 -build analysis form Chie Min CFX 5 July 12, 2001 23:19


All times are GMT -4. The time now is 00:32.