DNS 1-5 Quantification of Resolution: Difference between revisions
No edit summary |
|||
(27 intermediate revisions by 3 users not shown) | |||
Line 1: | Line 1: | ||
=HiFi-TURB-DLR rounded step= | |||
{{ | {{DNSHeader | ||
|area=1 | |area=1 | ||
|number=5 | |number=5 | ||
Line 11: | Line 11: | ||
==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>). | ||
Here, the mesh characteristic length takes into account | Here, the mesh characteristic length takes into account the degree of the DG polynomial approximation. In particular, it is defined as the cubic root of the ratio between the mesh element volume <math>\left(V\right)</math> and the number of DoFs <math>\left(N_{DoF}\right)</math> within the mesh element per equation | ||
<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 [[ | The comparison with respect to the Taylor microscale is shown in [[DNS_1-5_quantification_#figure6|Fig. 6]]. The maximum ratio <math>{\Delta}/{\eta_{T}}</math> within the outer layer of the boundary layer is approximately 0.4. | ||
In [[ | This outcome suggests that the current space resolution is sufficient to capture turbulence scales in the intertial range. | ||
In [[DNS_1-5_quantification_#figure7|Fig. 7]] is reported the comparison with respect to the Kolmogorov length scale. | |||
It is commonly accepted that DNS requirements are achieved when <math>{\Delta}/{\eta_{K}}\leq 5</math>. | |||
Current simulation shows above the flat plate upstream the rounded step a ratio <math>{\Delta}/{\eta_{K}}</math> below 5.5, while above the rounded step a ratio lower than 7.5. | |||
DNS requirements are thus not fulfilled in this last region. This is the reason why the present study is referred to as under-resolved DNS (uDNS). | |||
For future highly resolved simulations of this test case a mesh refinement is advised above and downstream the rounded step. | |||
<div id=" | <div id="figure6"></div> | ||
{|align="center" width=1200 | {|align="center" width=1200 | ||
|[[Image:DNS1-5_rounded_step_Taylor_scale.png|1200px]] | |[[Image:DNS1-5_rounded_step_Taylor_scale.png|1200px]] | ||
|- | |- | ||
|align="center"|'''Figure | |align="center"|'''Figure 6:''' HiFi-TURB-DLR rounded step, Re=78,490. Relation between the mesh size and the Taylor microscale at midspan using MIGALE with DG P3 (~300 million DoF/eqn). | ||
|} | |} | ||
<div id=" | <div id="figure7"></div> | ||
{|align="center" width=1200 | {|align="center" width=1200 | ||
|[[Image:DNS1-5_rounded_step_Kolmogorov_scale.png|1200px]] | |[[Image:DNS1-5_rounded_step_Kolmogorov_scale.png|1200px]] | ||
|- | |- | ||
|align="center"|'''Figure | |align="center"|'''Figure 7:''' HiFi-TURB-DLR rounded step, Re=78,490. 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 | The average wall resolution in streamwise (<math>{\Delta x^{+}}</math>), spanwise (<math>{\Delta z^{+}}</math>) and wall-normal (<math>{\Delta y^{+}}</math>) directions at different streamwise locations is reported in [[DNS_1-5_quantification#table2|Tab. 2]]. | ||
<div id="table2"></div> | |||
{|align="center" border="1" cellpadding="10" | {|align="center" border="1" cellpadding="10" | ||
|- | |- | ||
| | |align="center"|'''<math>{x/H}</math>'''||align="center"|'''<math>{\Delta x^{+}}</math>'''||align="center"|'''<math>{\Delta z^{+}}</math>||align="center"|'''<math>{\Delta y_{1}^{+}}</math>''' | ||
|| '''<math>{ | |||
| | |||
|| '''<math>{1 | |||
|- | |- | ||
| | |align="center"|'''<math>{-3.5}</math>'''||align="center"|'''<math>{17.8}</math>'''||align="center"|'''<math>{18.6}</math>'''||align="center"|'''<math>{1.30}</math>''' | ||
|| '''<math>{ | |||
| | |||
|| '''<math>{ | |||
|- | |- | ||
| | |align="center"|'''<math>{3}</math>'''||align="center"|'''<math>{6.53}</math>'''||align="center"|'''<math>{8.19}</math>'''||align="center"|'''<math>{0.55}</math>''' | ||
|| '''<math>{ | |||
| | |||
|| '''<math>{0. | |||
|- | |- | ||
| | |align="center"|'''<math>{4}</math>'''||align="center"|'''<math>{3.90}</math>'''||align="center"|'''<math>{5.71}</math>'''||align="center"|'''<math>{0.40}</math>''' | ||
|| '''<math>{ | |||
| | |||
|| '''<math>{0. | |||
|- | |- | ||
| | |align="center"|'''<math>{5}</math>'''||align="center"|'''<math>{9.22}</math>'''||align="center"|'''<math>{9.28}</math>'''||align="center"|'''<math>{0.66}</math>''' | ||
|| '''<math>{12.8}</math>''' | |- | ||
| | |align="center"|'''<math>{6}</math>'''||align="center"|'''<math>{12.8}</math>'''||align="center"|'''<math>{13.0}</math>'''||align="center"|'''<math>{0.93}</math>''' | ||
|| '''<math>{0.93}</math>''' | |||
|} | |} | ||
<center>'''Table 2:''' Wall space resolution at different streamwise locations</center> | |||
<br/> | <br/> | ||
<br/> | <br/> | ||
Line 72: | Line 62: | ||
As first step, an assessment of code MIGALE in closing the budgets is performed. | As first step, an assessment of code MIGALE in closing the budgets is performed. | ||
[[ | [[DNS 1-5 quantification#figure8|Fig. 8]] reports the budget of streamwise Reynolds-stress and TKE equations in a channel flow at <math>{Re_\tau = 180}</math> using a DG polynomial approximation of degree 5 on a mesh of <math>91\times 43\times 48</math> hexahedral elements (10.5 million DoF/eqn.). | ||
Domain dimension and reference results are given by DNS of [[DNS_1-5 quantification#1|Moser ''et al.'' (1999)]]. | |||
As the maximum value of the residual is <math>0.22\%</math> and <math>0.64\%</math> of the production peak for the Reynolds-stress xx and the TKE budgets, respectively, the results show that MIGALE code can close well the budgets when sufficient spatial and time resolution is considered. | |||
<div id="figure8"></div> | |||
{|align="center" width=500 | |||
|[[Image:Channel UniBG MIGALE Re xx budget.png|500px]]||[[Image:Channel UniBG MIGALE TKE budget.png|500px]] | |||
|- | |||
!align="center" colspan="2"|'''Figure 8:''' Channel flow at <math>Re_{\tau}=180</math> Reynolds-stress xx and TKE budgets: dissipation <math>\varepsilon</math>, production <math>P</math>, turbulent diffusion <math>D^{1}</math>, pressure diffusion <math>D^{2}</math>, viscous diffusion <math>D^{3}</math> and pressure strain <math>\Phi</math>. Solid lines from the DG P5 computation and symbols from [[DNS_1-5 quantification#1|Moser ''et al.'' (1999)]]. | |||
|} | |||
Considering now the budgets for the HiFi-TURB-DLR rounded step simulation, the conclusion obtained comparing the mesh characteristic length and the Kolmogorov length scale (see [[DNS_1-5_quantification_#figure7|Fig. 7]]) still holds. | |||
Indeed, on [[DNS_1-5_quantification_#figure9|Fig. 9]] and [[DNS_1-5_quantification_#figure10|Fig. 10]] are reported the budgets for two different steamwise locations, i.e., <math>x/H=-3</math> (near the checkpoint) and <math>x/H=4</math> (above the rounded step, in the region of incipient separation), respectively. As can be seen, budget closure has not been achieved, especially above the rounded step, where the residual is larger than the production term. The oscillating behaviour of the profiles of the different terms can be ascribed to a lack of spatial resolution. | |||
This result suggests that for future campaigns it is recommended to increase further the grid density. | |||
<div id=" | <div id="figure9"></div> | ||
{|align="center" width= | {|align="center" width=500 | ||
|[[Image: | |[[Image:DNS1-5 rounded step Re xx budget -300.png|500px]]||[[Image:DNS1-5 rounded step TKE budget -300.png|500px]] | ||
|[[Image: | |||
|- | |- | ||
!align="center" colspan="2"|'''Figure 9:''' HiFi-TURB-DLR rounded step, Reynolds-stress xx and TKE budgets at location <math>x/H=-3</math>: convection <math>C</math>, production <math>P</math>, turbulent diffusion <math>D^{1}</math>, pressure diffusion <math>D^{2}</math>, viscous diffusion <math>D^{3}</math>, pressure strain <math>\Phi</math>, pressure dilatation <math>\Phi'</math>, dissipation <math>\varepsilon</math> and density fluctuation <math>K</math>. | |||
|} | |} | ||
<div id="figure10"></div> | |||
{|align="center" width=500 | |||
|[[Image:DNS1-5 rounded step Re xx budget 400.png|500px]]||[[Image:DNS1-5 rounded step TKE budget 400.png|500px]] | |||
|- | |||
!align="center" colspan="2"|'''Figure 10:''' HiFi-TURB-DLR rounded step, Reynolds-stress xx and TKE budgets at location <math>x/H=4</math>: convection <math>C</math>, production <math>P</math>, turbulent diffusion <math>D^{1}</math>, pressure diffusion <math>D^{2}</math>, viscous diffusion <math>D^{3}</math>, pressure strain <math>\Phi</math>, pressure dilatation <math>\Phi'</math>, dissipation <math>\varepsilon</math> and density fluctuation <math>K</math>. | |||
|} | |||
<br/> | |||
==References== | |||
#<div id="1">'''Moser, R. D., Kim, J., Mansour, N. N. (1999)''': Direct numerical simulation of turbulent channel flow up to Re_tau 590. ''Physics of Fluids'', Vol. 11(4), pp.943-945.</div> | |||
---- | ---- | ||
{{ACContribs | {{ACContribs | ||
Line 100: | Line 103: | ||
| organisation=Università degli studi di Bergamo (UniBG) | | organisation=Università degli studi di Bergamo (UniBG) | ||
}} | }} | ||
{{ | {{DNSHeader | ||
|area=1 | |area=1 | ||
|number=5 | |number=5 |
Latest revision as of 16:12, 17 February 2023
HiFi-TURB-DLR rounded step
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 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. 6. The maximum ratio within the outer layer of the boundary layer is approximately 0.4. This outcome suggests that the current space resolution is sufficient to capture turbulence scales in the intertial range. In Fig. 7 is reported the comparison with respect to the Kolmogorov length scale. It is commonly accepted that DNS requirements are achieved when . Current simulation shows above the flat plate upstream the rounded step a ratio below 5.5, while above the rounded step a ratio lower than 7.5. DNS requirements are thus not fulfilled in this last region. This is the reason why the present study is referred to as under-resolved DNS (uDNS). For future highly resolved simulations of this test case a mesh refinement is advised above and downstream the rounded step.
Figure 6: HiFi-TURB-DLR rounded step, Re=78,490. Relation between the mesh size and the Taylor microscale at midspan using MIGALE with DG P3 (~300 million DoF/eqn). |
Figure 7: HiFi-TURB-DLR rounded step, Re=78,490. 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 streamwise (), spanwise () and wall-normal () directions at different streamwise locations is reported in Tab. 2.
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. 8 reports the budget of streamwise 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.). Domain dimension and reference results are given by DNS of Moser et al. (1999). As the maximum value of the residual is and of the production peak for the Reynolds-stress xx and the TKE budgets, respectively, the results show that MIGALE code can close well the budgets when sufficient spatial and time resolution is considered.
Figure 8: Channel flow at Reynolds-stress xx and TKE budgets: dissipation , production , turbulent diffusion , pressure diffusion , viscous diffusion and pressure strain . Solid lines from the DG P5 computation and symbols from Moser et al. (1999). |
---|
Considering now the budgets for the HiFi-TURB-DLR rounded step simulation, the conclusion obtained comparing the mesh characteristic length and the Kolmogorov length scale (see Fig. 7) still holds. Indeed, on Fig. 9 and Fig. 10 are reported the budgets for two different steamwise locations, i.e., (near the checkpoint) and (above the rounded step, in the region of incipient separation), respectively. As can be seen, budget closure has not been achieved, especially above the rounded step, where the residual is larger than the production term. The oscillating behaviour of the profiles of the different terms can be ascribed to a lack of spatial resolution. This result suggests that for future campaigns it is recommended to increase further the grid density.
References
- Moser, R. D., Kim, J., Mansour, N. N. (1999): Direct numerical simulation of turbulent channel flow up to Re_tau 590. Physics of Fluids, Vol. 11(4), pp.943-945.
Contributed by: Francesco Bassi, Alessandro Colombo, Francesco Carlo Massa — Università degli studi di Bergamo (UniBG)
© copyright ERCOFTAC 2024