|
[Sponsors] |
[Other] Non-Orthogonality calculation in high aspect elements |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
Senior Member
Vangelis Skaperdas
Join Date: Mar 2009
Location: Thessaloniki, Greece
Posts: 287
Rep Power: 21 ![]() |
Hello everyone,
I am in the process of creating high aspect ratio layer elements for OpenFOAM and have been surprised by the sensitivity of the calculation of non-orthogonality in cases where the aspect ratio of the elements is high. Consider the following example: I have two perfect hexa layer elements of width 1000 and height 1, so an aspect ratio of 1000. All angles are 90deg. In this case I have 0 non-orthogonality. However if I begin to move one of the four common nodes along the height direction, then non-orthogonality increases at an amazing rate. So for a movement of 1.5% of the height, non-orthogonality reaches 60 and for a movement of about 10% orthogonality reaches 85 This is an idealized case where all edge lengths are the same, angles are orthgonal etc. In a real geometry where layers are extruded it is practically impossible to satisfy non-orthogonality criteria due to this extreme sensitivity. Should it be this way? Thanks Vangelis Last edited by vangelis; November 7, 2017 at 03:00. |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 403
Rep Power: 25 ![]() |
It might be worthwhile to plot the centroids of cells on either side of the face in question (i.e., the values computed by OpenFOAM in primitiveMesh), to see if they actually make sense.
|
|
![]() |
![]() |
![]() |
![]() |
#3 |
Senior Member
Vangelis Skaperdas
Join Date: Mar 2009
Location: Thessaloniki, Greece
Posts: 287
Rep Power: 21 ![]() |
Hi Sandeep
Thanks for your reply. Indeed the calculated centroids move and the vector that connects them forms an angle with the normal of the common facets of the two hexas. I was wondering however if the calculation of centroids is good for such cases. Why should a movement in the normal direction, affect sideways the centroids of the two elements? I mean these elements are near perfect to the eye after the nodel movement. In the end it is impossible to create mesh with high aspect ratio that can satisfy orthogonality, in real cases. |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 403
Rep Power: 25 ![]() |
Are you using the OpenCFD version or OpenFOAM-extend? Might be useful to compare the two as well.
|
|
![]() |
![]() |
![]() |
![]() |
#5 |
Senior Member
Vangelis Skaperdas
Join Date: Mar 2009
Location: Thessaloniki, Greece
Posts: 287
Rep Power: 21 ![]() |
I am using OpenCFD version but I believe it is the same for all.
|
|
![]() |
![]() |
![]() |
![]() |
#6 |
Senior Member
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 403
Rep Power: 25 ![]() |
I had changed the centroid calculation a bit on the extend side, so I'm curious to see if there's a difference, which is why I mentioned it. If you look at primitiveMeshCellCentresAndVols.C, you'll notice there's a difference between versions.
|
|
![]() |
![]() |
![]() |
![]() |
#7 |
Senior Member
Vangelis Skaperdas
Join Date: Mar 2009
Location: Thessaloniki, Greece
Posts: 287
Rep Power: 21 ![]() |
Thank you Sandeep we will look it up.
So in your custom version do you see the same sensitivity and tendency for the orthogonality to shoot up when making small such nodal movements? |
|
![]() |
![]() |
![]() |
![]() |
#8 |
Senior Member
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 403
Rep Power: 25 ![]() |
I've never tried a test like this to be honest, which is why I'm curious about your results. It looks like your face is essentially unchanged, so its normal / centroid would be unchanged as well. Which leaves the adjacent cell centroid as the only variable.
|
|
![]() |
![]() |
![]() |
![]() |
#9 |
Senior Member
M. Montero
Join Date: Mar 2009
Location: Madrid
Posts: 153
Rep Power: 17 ![]() |
Any news about this topic? I am facing a similar problem where orthogonality reaches 89.5 but due to high aspect ratio cells with y+ less than 1.
|
|
![]() |
![]() |
![]() |
Tags |
aspect ratio, non-orthogonal |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Ansys SIG$ILL error | loth | ANSYS | 3 | December 24, 2015 05:31 |
aspect ratio of boundary layer elements in COMSOL | TryMyBest | Mesh Generation & Pre-Processing | 0 | March 11, 2013 06:53 |
High aspect ratio in viscous sublayer? | Marius | Main CFD Forum | 0 | December 22, 2005 08:15 |
high aspect ratio elements | Olivier | CFX | 5 | April 3, 2003 06:55 |
CFX4.3 -build analysis form | Chie Min | CFX | 5 | July 12, 2001 23:19 |