Mesh study external aerodynamics problem

Hi @math_marques96

The problem is with the inflation layer. These bad elements are present where the prism layer cells end. I recommend you to improve the inflation layers.

Thanks
Ani

4 Likes

Aha, if it is ‘bad’ Snappy prism cells causing these bumps then this may be a case where you could try INCREASING certain default quality checks that Snappy uses when deciding which of the perfect prism cells to keep, otherwise is there an actual ‘bump’ in the geometry there?

1 Like

No sir, I verified that there is no imperfection in the geometry. The problem is with the inflation layer. I made a new mesh with better inflation layer settings and it solved the problem.

OLD MESH


NEW MESH

OLD MESH TRAILING EDGE


NEW MESH TRAILING EDGE

Thanks
Ani

4 Likes

Great! Looks good but I see still some pretty nasty gradients from prism cells to the volume cells which I have been led to believe causes issues with results accuracy (maybe another reason to use ‘Relative layering’???).

EDIT: Would be interesting to see Y+ at surface and Eddy Viscosity Ratio plots there…

and

I guess I should have deduced that from your last image :wink: of post #6

The yplus value on the surface seems to be good :slightly_smiling_face:.


The simulation got converged in 425 iterations

1 Like

For me I would have to say the jury is still out on that Y+ mapping, is this for Full resolution Y+ values or Wall functions at 30<Y+<300

Looks like a lot of Y+ could be in bad territory between using either Full res or Wall functions…

I think that Bad area for us using κ-ω SST is 10<Y+<30 according to this (but I just recently learned this and have not run any validity test on this range)

I hope you or someone else tries to make those Eddy Viscosity Ratio plots, saving me the time to figure them out :wink:

Yup @Lwhitson2 is absolutely correct about the yplus value for k-omega SST model. The shear stress transport (SST) k-ω model is a two-equation model that combines the advantages of the k-ε and k-ω models. This model uses a blending function to gradually transition from standard k-ω in the inner boundary layer to standard k-ε in the outer boundary layer. The standard k-ω model performs better in the near-wall region while the standard k-ε performs better in the free-stream region.

Also, there are not much cells with yplus between 10 &30.

These Wall treatment method are used by some other CFD software. Just wanted to share this with you guyzz.

7 Likes

@DaleKramer you’re never going to get a perfect Y+ value everywhere on your model. The key is to get Y+ correct in the most important parts of the flow. Overall I think the Y+ looks really good on this simulation after @anirudh2821998 made his updates to the mesh.

4 Likes

Yes but I do not think I could easily see that from Ani’s first Y+ mapping presentation. The second mapping of 10<Y+<30 does confirm this for me.

Very good Ani !!!

EDIT: And just to confirm, this will be a Full res Y+ simulation, correct?

also, I like to show mappings like that with a Grey area for Y+<10 , I know we can do that I just don’t remember off hand how I did it, or Y+<1 like this:

1 Like

Sir, Is this what you are suggesting??

Not really, I added an example while you were posting, it still shows rainbow 1 to 10

Well, there is a white to brown color scheme but I can’t see any raibow type color scheme. Did you use some other post processing tool?

I don’t believe so but it was with old workbench, maybe can’t do this anymore…

I did add the ‘Gray is y+<1’ text in an image processor tho.

Hello Ani,

For this inflation layers parameters you used relative layering?
And this layer calculation is for a full resolution right?
Did you just changed the inflation layers parameters?

Would it be good to alter the surface refinement level, region refinement level or the bounding box resolution in Ani´s mesh to have a smoothier transition from prism cells to region cells?

No

Yes

I just modified the size of the bounding box(made it a bit larger) and changed only the inflation layer parameter.

It is a must to simulate your problem with atleast 2 different quality mesh. This will provide you with an idea of the deviation in results with mesh fineness and will also serve the purpose for mesh independence study.

1 Like

To add to this, I would suggest running a full grid convergence study using at least three grids and verifying the order of convergence for your model. This would satisfy the scientific/research level requirements and is great experience to have in validating your model. Here is a link to a great article by NASA that will help you in the calculations. Near the bottom is an example of the grid refinement study.

Examining Spatial (Grid) Convergence

4 Likes

I have some questions about the layer refinement. In this mesh you did, what value of Y+ were you aiming? And what value of expansion rate you used?

@math_marques96 it depends on what turbulence model you plan on using. I would start with the \kappa \text - \epsilon model and target a y^+ of about 50 with a growth factor of around 1.2 - 1.3 and a last layer thickness somewhere around 83\% of your standard surface mesh size. After running a few cases with \kappa \text - \epsilon you will know whether or not separation is an issue. If so, I would switch to the \kappa \text - \omega\,SST model and then target 1 \leq y^+ \leq 5 in the areas where the separation occurs.

The first model will get you about 90\% of the way there and the second will get you that last 10\%. It all depends on your time requirements and what you are looking to get out of the simulation.

3 Likes

Hello,
I made some changes in the layers refinement and acquired a better mesh. I am in doubt about a layer parameter.
I found this forum Boundary layer inflation thickness - #4 by AnetaNovak and it was not clear for me if what parameter is the “final layer thickness”. On this forum there is also a “first layer thickness” parameter that not appear in my simulation.
I uploaded this draw of a boundary layer, to show my question, the parameter “final layer thickness” is the A or B? Is it possible to choose between using the “final layer thickness” or “first layer thickness”, I did not understand if it is the same also.