Okay there is now a completely clean simulation called “For Dale” in the list there. Have at it, big guy.
I’m cooking on these mesh reducing techniques. Check this out:
Surface refinement 7, relative layer sizing, Final layer thickness=0.005, ER=2
Surface refinement 4-7, absolute sizing, Final layer thickness=0.005, ER=2. Note BL uniformity through three levels of surface refinement.
Great!
I am glad that absolute may work here, I wish I could remember what project I was on when I chose to start using the default relative layering all the time, oh well.
ER 2 is likely pushing it but I hope we have that new metric of Eddy Viscosity Ratio (EVR from now on ) plots to confirm…
Maybe I should wait on that big sim until you get sorted out, looks like you are close to trying to get a whole plane BL refinement with the fruits of your cat and mousing
I’m running a sim on that reduced mesh now and it’s showing a similar instability to the 25MM cell mesh, which is both a good and a bad sign. Good that it’s acting like the much desnser mesh, bad that it’s not even remotely converging.
How are you going to plot Eddy Viscosity Ratio? Paraview?
Likely Paraview, but I am a total novice at that, which is why I have been hinting about EVR plots since @LWhitson2 (hint hint) pointed us to them, hoping someone else would pioneer them here at SimScale…
I still think we should do the big sim on your best full aircraft full res mesh, you just say the word… I will look now at sim you are running or last sim if it is still in Job Status…
EDIT: Wow I thought it had blown up, looks like it recovered nicely and will eventually let us have a look at EVR plots somehow. I hope there wasn’t some result we needed to add to figure EVR with Paraview…
EDIT: Interesting all looks pretty stable at 400s EXCEPT for pressure moments in pitch and roll ??
But its looks like they will stabilize by the end…
EDIT: I’d say that rudder is still pretty effective there with left yaw moment of ~2800 Nm and there is a left roll moment of ~1900 Nm with pitch moment about ~0
EDIT: IT FINISHED YEAH, SUCCESS !!!
Can’t wait for your next post
Darn, I just took a look for the first time and I do not think we will have any EVR results to plot. I was hoping that we had the numerator (nut) and denominator (nu) already available but I guess everyone needs to vote on this Turbulent Viscosity Ratio Calculations since we only have nut in the results… .
@jousefm how do we get this feature quickly, I think it will substantially reduce all the speculation about whether we have adequately resolved our boundary layers.
@DaleKramer the nu value is taken directly from the fluid properties. It is the Kinematic Viscosity. You’ve got everything you need to make a nut/nu field using Paraview. I would love to see custom calculations in the post processing available in SimScale though. Also, it would be great if they would add a Line Interval Convolution (LIC) option to the slices.
Ahhhh! thanks, makes sense, should have realized that, dah, you know forest for the trees etc…
I see you have a new image
Yes somewhat disappointing. I have new results from my meshing experiments though and am rebuilding an even thinner mesh on the Thrush with 12 layers. It may fail though, as it ends up being a LOT of cells, but the results are very good. The smallest cells below are refinement level 8 and you can see the entire BL has been generated. The BL inflation algorithm really seems to like small cells whenever there is any curvature to the surface. Having a variable surface refinement and absolute BL sizing lets you do this without having an insanely large mesh.
Interestingly, I have found that the BL cells are best retained when you use an FLT ratio somewhat lower… if the small cells are too small it pukes too. I did a couple trials varying the max refinement level with absolute BL sizing enabled. I found good results if I went to one surface refinement level LOWER than your typical ratio, which IIRC was .30-.40. Here I’m .22 with surface refinement from 4-8.
Sorry, its late and I am a little confused. The absolute BL testing ended up being not as good so you went back to using lower FLT ratio ???
If so, I think that 0.3-0.4 FLT ratio range was simply my own derived range. It is good that you may be optimizing that… Although it is hard to stay too tight to it because FLT ratio changes so much from level to level. Intuitively a lower FLT ratio (larger size gradient) would seem to be expect worse sim results but perhaps a better layering (that is a way out there thought from a tired mind) ??
I am glad you are becoming such an expert at this Isn’t this frustrating but fun
Yeah maybe a bit confusing. I found that while the layering algo liked smaller cells, that getting too small tended to disrupt their formation almost as much as being too coarse. The variable surface refinement lets it add detail where necessary, but I found that limiting the max surface refinement level to the next step below your 0.3-0.4 range gave the best results. In the above case we’re kind of in the middle, but your range would have recommended surface level 9 so I went with level 8.
I’m just about complete on the best mesh yet at about 15MM cells. Interestingly, I’m finding the more quality checks I just turn off the better it is and it still doesn’t fail the SimScale quality checks. Just now I took face twist from 0.0005 to -1 (turning it off) and my retained cells jumped by 7%.
Very interesting, I did not much research on that, I just warned people that they were on their own with quality after you start changing the defaults
As I understand it 114 illegal faces is getting up where we should be concerned. I just try for 0 illegals (like Trump )
I have never figured out this newish and sometimes missing ‘Mesh Quality Metrics’ at the end…
Mesh quality metrics:
volumeRatio
min: 1
max: 2543.5148081883954
average: 1.9350681767938023
standard deviation 2.330299409787477
median: 1.1802326178923026
0-th percentile: 1
20-th percentile: 1.0036717141371905
40-th percentile: 1.0476639184279137
60-th percentile: 1.4999964640160461
80-th percentile: 1.74773126868174
100-th percentile: 2543.5148081883954
EDIT: At least I eventually just found this with a Google openfoam “volumeratio” search:
‘- VolumeRatio is a new metric, which computes the relative volume change between two adjacent cells’
AND further searching says it is from SnappyHexMesh.
aha…
Which means that overall the mesh transition is good(looking at the 80th percentile) but at some places there is a big transition between big elements and small elements.
To shed some light on the percentile issue…
The volumeRatio measures ratio of volumes between each neighboring cells. (So for each mesh face, it measures the volume ratio between the two mesh cells that are adjacent to the face.)
So then each mesh has a lot of volume ratios. For each neighboring cells we get one ratio.
Since we have a set of volume ratios, we can compute some statistics. The min, max, average, standard deviation and median are the general statistics.
The percentiles are another way of representing the distribution of the volume ratios, also see:
Examples:
0th percentile: 1.00000018 - it means that 0% of all volume ratios is lower than 1.00000018 - so this is the same as the minimum - because no value is lower than minimum
80th percentile: 1.068 - it means that 80% of all volume ratios is lower than 1.068
Hope that makes sense.
Jousef
Thanks Jousef.
That statement relates the most to the issue I have with volumeRatio.
I am not sure of how each cell face has only two mesh cells adjacent to it, or of the definition of CFD adjacent cells, or of the definition or CFD neighboring cells
Maybe I have to get my head around the definition of adjacent cells in CFD terms.
In CFD terms, I believe that our cells have faces, edges and corner points that define them in our mind (3 items).
So, the combination of items that cells can share would be, pick 1, 2 or 3 from a group of 3 (7 permutations I think).
What permutation would CFD consider adjacent cells to be?
My simple brain relates the concept only to the most common use of adjacent, as in adjacent rooms. Which is the simple permutation of faces only, described as follows :
In adjacent room terms, I think that they are most commonly thought of as rooms that have a common wall (this common wall could have different dimensions for each room).
Although rooms that are on different floors would not commonly be considered adjacent, likely they would be if they had an access door in the ceiling or floor.
I also do not think that rooms that only share a common edge or point of a wall, would be considered adjacent, since even though there is some point/points in common, there would be no way to have a door in a common face large enough to have that secret rendezvous (without creating some sort of private additional room/hallway type structure ).
Wow, so much gymnastics to explain my issue…
If I get all the above straight then maybe I can ‘see’ volumeRatio.
After all I think we need to ‘see’ things in our mind before we can describe them, use the concept wisely or even truly understand them…
I am pretty much a robot right now with my SimScaling but I am trying…
My interpretation of this for hexahedron meshes below. By this definition, each 3D cell will have six volume ratios.
@jousefm, I’m confused why volume ratio is a measure of quality. How does a high VR affect the solver? It is not surprising that there are some very high volume ratios in this mesh because undoubtedly there are some areas with my very tiny boundary layer inflated on level 4 surface refinement. Recall that I’m using a variable surface refinement which has resulted in a much more efficient mesh (~50% smaller than using only the highest level of refinement) while still maintaining good BL coverage.
Because I like pictures, here’s an example. Is this a bad thing?
The bigger question that keeps coming up for me is what fundamentally makes a good mesh? Lack of certain types of cells? Volume ratios? I keep eliminating quality checks and the BL gets better and it still passes. As in the last case, increasing numbers of “illegal” faces exist but for instance on the last mesh that’s 114 faces out of millions.
Are we sure that our meshes only create cells of 6 faces? (I thought I saw some that might be 4… maybe even 5)