DNS 1-3 Quantification of Resolution: Difference between revisions
Line 31: | Line 31: | ||
[[lib:DNS_1-3_statistical_#figure11|Fig. 11]] and [[lib:DNS_1-3_statistical_#figure12|Fig. 12]] show a fair agreement between both numerical simulations, proving that the approach used is an optimal strategy. Moreover, the points where the DNS performed with Alya deviates from the results of [[lib:DNS_1-3_description#1|Ohlsson ''et al.'' (2010)]], this is in favour to move closer to the experimental data of [[lib:DNS_1-3_description#5|Cherry ''et al.'' (2008)]]. Thus, a good agreement is also obtained with the experiment. The inclusion of a roughness element on the long inlet duct, to trigger the turbulence, is shown to have a negligible effect on the quality of the results. | [[lib:DNS_1-3_statistical_#figure11|Fig. 11]] and [[lib:DNS_1-3_statistical_#figure12|Fig. 12]] show a fair agreement between both numerical simulations, proving that the approach used is an optimal strategy. Moreover, the points where the DNS performed with Alya deviates from the results of [[lib:DNS_1-3_description#1|Ohlsson ''et al.'' (2010)]], this is in favour to move closer to the experimental data of [[lib:DNS_1-3_description#5|Cherry ''et al.'' (2008)]]. Thus, a good agreement is also obtained with the experiment. The inclusion of a roughness element on the long inlet duct, to trigger the turbulence, is shown to have a negligible effect on the quality of the results. | ||
===Other turbulent statistics=== | ===Other turbulent statistics=== |
Revision as of 12:58, 28 January 2022
Quantification of resolution
Mesh resolution
The mesh resolution is quantified by obtaining a relation between the mesh characteristic length () and the characteristics lengths, i.e., the Taylor microscale () and Kolmogorov length scale (). The former relation is shown in Fig. 9 while the latter is reported in Fig. 10. As it can be seen, both relations indicate that the resolution achieved by the present grid is at DNS level. In particular, it is commonly accepted that DNS is achieved when , as shown in Fig. 10.
Figure 9: Stanford double diffuser, Alya DNS-250M DoF, relation between the mesh size and the Taylor microscale. |
Figure 10: Stanford double diffuser, Alya DNS-250M DoF, relation between the mesh size and the Kolmogorov length scale. |
Solution verification
This section aims to provide an insight on the dataset freely available in this wiki. Simulation results are compared with the DNS of Ohlsson et al. (2010) and the experimental data of Cherry et al. (2008). The grid resolution has been proven to be fine enough to be at DNS level. An animation of the fluctuations of this flow can be seen in this link.
Velocity profiles and RMS quantities
The 10000 Reynolds number diffuser has been simulated and compared with the results with previous DNS data from Ohlsson et al. (2010) and experimental data from Cherry et al. (2008). Fig. 11 and Fig. 12 show the averaged velocities and the RMS quantities along the characteristic lines of the diffuser for the experimental data and both DNS simulations. Finally, Fig. 13 and Fig. 14 show a comparison of the contours of the span-wise velocity and span-wise rms fluctuations respectively.
Fig. 11 and Fig. 12 show a fair agreement between both numerical simulations, proving that the approach used is an optimal strategy. Moreover, the points where the DNS performed with Alya deviates from the results of Ohlsson et al. (2010), this is in favour to move closer to the experimental data of Cherry et al. (2008). Thus, a good agreement is also obtained with the experiment. The inclusion of a roughness element on the long inlet duct, to trigger the turbulence, is shown to have a negligible effect on the quality of the results.
Other turbulent statistics
Other relevant statistics from the database are presented in this section. In particular, the turbulent kinetic energy (TKE) and the characteristic length scales (i.e., the Taylor micro-scale and the Kolmogorov length-scale ) are presented in Fig. 15 along some characteristic contours of the diffuser.
Turbulent kinetic energy budget equation terms
The turbulent kinetic energy budget equation terms are presented in this section. As it can be seen in Fig. 16, the magnitudes of the residuals for the TKE budget equations is very small, of the order of , which is deemed sufficient for the present case.
Figure 16: Stanford double diffuser, Re=10000. Turbulent kinetic energy (TKE) budget balance obtained with Alya on characteristic lines of the diffuser. |
Contributed by: Oriol Lehmkuhl, Arnau Miro — Barcelona Supercomputing Center (BSC)
© copyright ERCOFTAC 2024