DNS 1-5 Quantification of Resolution: Difference between revisions

From KBwiki
Jump to navigation Jump to search
Line 6: Line 6:
__TOC__
__TOC__
= Quantification of resolution =
= Quantification of resolution =
This section provides details of the solution accuracy obtained by tackling the rounded step DNS with MIGALE.
After providing details of the mesh resolution in comparison with spatial turbulent scales, a discussion on the closure of the Reynolds stress equations budgets is given.
==Mesh resolution==
==Mesh resolution==
The mesh resolution is quantified by comparing the mesh characteristic length (<math>{\Delta}</math>) with the characteristic lengths of the turbulence, i.e., the Taylor microscale (<math>{\eta_T}</math>) and Kolmogorov length scale (<math>{\eta_K}</math>).
The mesh resolution is quantified by comparing the mesh characteristic length (<math>{\Delta}</math>) with the characteristic lengths of the turbulence, i.e., the Taylor microscale (<math>{\eta_T}</math>) and Kolmogorov length scale (<math>{\eta_K}</math>).
Line 12: Line 15:
<math>\Delta=\sqrt[3]{\dfrac{V}{N_{DoF}}}</math>
<math>\Delta=\sqrt[3]{\dfrac{V}{N_{DoF}}}</math>


The comparison with respect to the Taylor microscale is shown in [[lib:DNS1-5_quantification_Taylor_#figure3|Fig. 3]]. The maximum ratio <math>{\Delta}/{\eta_{T}}</math> within the outer layer of the boundary layer is approximately 0.4.
The comparison with respect to the Taylor microscale is shown in [[lib:DNS_1-5_quantification_#figure3|Fig. 3]]. The maximum ratio <math>{\Delta}/{\eta_{T}}</math> within the outer layer of the boundary layer is approximately 0.4.
In [[lib:DNS1-5_quantification_Kolmogorov_#figure4|Fig. 4]] is reported the comparison with respect to the Kolmogorov length scale. Above the flat plate upstream the rounded step the ratio <math>{\Delta}/{\eta_{K}}</math> is below 5.5, while above the rounded step is lower than 7.5.
In [[lib:DNS_1-5_quantification_#figure4|Fig. 4]] is reported the comparison with respect to the Kolmogorov length scale. Above the flat plate upstream the rounded step the ratio <math>{\Delta}/{\eta_{K}}</math> is below 5.5, while above the rounded step is lower than 7.5.


<div id="figure3"></div>
<div id="figure3"></div>
Line 29: Line 32:
|}
|}


The average wall resolution at the checkpoint location <math>\left(x_{ref}/H=-3.5\right)</math> is <math>x^{+}=18</math> in stream direction, <math>y_{1}^{+}=1</math> in normal direction, and <math>z^{+}=19.5</math> in span direction.
The average wall resolution in stream (<math>{x^{+}}</math>), span (<math>{z^{+}}</math>) and wall normal (<math>{y^{+}}</math>) directions at different streamwise locations is reported in the following.
 
{|align="center" border="1" cellpadding="10"
|  '''<math>{x/H}</math>'''
|| '''<math>{x^{+}}</math>'''
|  '''<math>{z^{+}}</math>'''
|| '''<math>{y_{1}^{+}}</math>'''
|-
|  '''<math>{-3.5}</math>'''
|| '''<math>{17.8}</math>'''
|  '''<math>{18.6}</math>'''
|| '''<math>{1.30}</math>'''
|-
|  '''<math>{3}</math>'''
|| '''<math>{6.53}</math>'''
|  '''<math>{8.19}</math>'''
|| '''<math>{0.55}</math>'''
|-
|  '''<math>{4}</math>'''
|| '''<math>{3.90}</math>'''
|  '''<math>{5.71}</math>'''
|| '''<math>{0.40}</math>'''
|-
|  '''<math>{5}</math>'''
|| '''<math>{9.22}</math>'''
|  '''<math>{9.28}</math>'''
|| '''<math>{0.66}</math>'''
|-
|  '''<math>{6}</math>'''
|| '''<math>{12.8}</math>'''
|  '''<math>{13.0}</math>'''
|| '''<math>{0.93}</math>'''
|}
<br/>
<br/>


==Solution verification==
==Solution verification==

Revision as of 11:06, 2 December 2022


Front Page

Description

Computational Details

Quantification of Resolution

Statistical Data

Instantaneous Data

Storage Format

Quantification of resolution

This section provides details of the solution accuracy obtained by tackling the rounded step DNS with MIGALE. After providing details of the mesh resolution in comparison with spatial turbulent scales, a discussion on the closure of the Reynolds stress equations budgets is given.

Mesh resolution

The mesh resolution is quantified by comparing the mesh characteristic length () with the characteristic lengths of the turbulence, i.e., the Taylor microscale () and Kolmogorov length scale (). Here, the mesh characteristic length takes into account of the degree of the DG polynomial approximation. In particular, it is defined as the cubic root of the ratio between the mesh element volume and the number of DoFs within the mesh element per equation

The comparison with respect to the Taylor microscale is shown in Fig. 3. The maximum ratio within the outer layer of the boundary layer is approximately 0.4. In Fig. 4 is reported the comparison with respect to the Kolmogorov length scale. Above the flat plate upstream the rounded step the ratio is below 5.5, while above the rounded step is lower than 7.5.

DNS1-5 rounded step Taylor scale.png
Figure 3: Rounded step case, Re=78490. Relation between the mesh size and the Taylor microscale at midspan using MIGALE with DG P3 (~300 million DoF/eqn).
DNS1-5 rounded step Kolmogorov scale.png
Figure 4: Rounded step case, Re=78490. Relation between the mesh size and the Kolmogorov length scale at midspan using MIGALE with DG P3 (~300 million DoF/eqn).

The average wall resolution in stream (), span () and wall normal () directions at different streamwise locations is reported in the following.



Solution verification

One way to verify that the DNS are properly resolved is to examine the budget of the Reynolds-stress equations and the turbulent kinetic energy (TKE) equation.

As first step, an assessment of code MIGALE in closing the budgets is performed. Fig. 6 reports the budget of Reynolds-stress and TKE equations in a channel flow at using a DG polynomial approximation of degree 5 on a mesh of hexahedral elements (10.5 million DoF/eqn.). Reference results are the DNS data of Moser-Kim-Mansour


To this aim, the figures below show the turbulent kinetic energy equation budgets and the Reynolds stress streamwise direction budget in a channel flow at . The results are compared with the DNS data of Hoyas and Jimenez (2008)

(in circles) and Moser et al. (1999) (in triangles) that can be found freely available online. As it can be seen, the residuals close well within an order of magnitude with respect to the production. This proves that the methodology used in this entry is sound and that Alya can close the budgets if given the necessary grid resolution and integration time.


One way to verify that the DNS are properly resolved is to examine the residuals of the Reynolds- stress budget equations. These residuals are among the statistical volume data to be provided as described in Statistical Data section.



Contributed by: Francesco Bassi, Alessandro Colombo, Francesco Carlo Massa — Università degli studi di Bergamo (UniBG)

Front Page

Description

Computational Details

Quantification of Resolution

Statistical Data

Instantaneous Data

Storage Format


© copyright ERCOFTAC 2024