search

ucns3d-team / UCNS3D

Unstructured Compressible Navier Stokes 3D code (UCNS3D)
https://ucns3d.com
GNU General Public License v3.0
247 stars 90 forks source link

Info on capability of iLES in the code #51

Closed SRkumar97 closed 1 year ago

SRkumar97 commented 1 year ago

Hello!

I went through the 2022 whitepaper doc provided in the Git, and came across its features - RANS, iLES and DDES

In testcases (both TESTS and TESTS2), the iLES capability is demonstrated for the 3D Taylor Green Vortex case. In the whitepaper an additional SD7003 airfoil case iLES at Mach 0.2 is demonstrated.

I just wish to know if the iLES module is capable of supersonic turbulent boundary layers? Towards Mach 2 and above. Also, is there a way to introduce synthetic turbulence signal at the inlet for iLES?

SRkumar97 commented 1 year ago

I have a minor doubt wrt the parameter "Initial conditions profile (4 DEFAULT)" in the UCNS3D.DAT file. I presume that the parameter is explained in the inputs.md file with just a different name "CODE_PROFILE" https://github.com/ucns3d-team/UCNS3D/blob/master/docs/inputs.md My doubt here is that value of 95 is specified in the TGV iLES case - however the description list for CODE_PROFILE doesn't include this value. It would be of great help, if someone could explain the relevance of this value = 95 for iLES.

TakisCFD commented 1 year ago

the 95 corresponds to the initial condition setup for the case you want to run. (have a look at other initial conditions codes in profile.f90). The code profile corresponds on how you want the code to perform (in terms of some additional settings found in parameters.f90) unless you want to stick with the default behaviour of the code (option 0).

SRkumar97 commented 1 year ago

In the UCNS3D whitepaper it is mentioned that - The physical domain $\Omega$ consists of any combination of conforming tetrahedral, hexahedral, prism or pyramids in 3D, and quadrilateral or triangular in 2D.

My doubt here is whether the code supports purely structured meshes - quadrilateral in 2D, hexahedral in 3D? I see that the TGV testcase has got combination of hexahedral and tetrahedral cells; can the entire domain be purely hexahedral?

SRkumar97 commented 1 year ago

In the UCNS3D whitepaper it is mentioned that - The physical domain Ω consists of any combination of conforming tetrahedral, hexahedral, prism or pyramids in 3D, and quadrilateral or triangular in 2D.

My doubt here is whether the code supports purely structured meshes - quadrilateral in 2D, hexahedral in 3D? I see that the TGV testcase has got combination of hexahedral and tetrahedral cells; can the entire domain be purely hexahedral?

Tested out pure hexahedral grids. The code runs smoothly without any trouble.