DNS 1-3 Quantification of Resolution
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. 1 while the latter is reported in Fig. 2. 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. 2.
Figure 1: Stanford double diffuser, Alya DNS-250M DoF, relation between the mesh size and the Taylor microscale. |
Figure 2: 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. 3 and Fig. 4 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. 5 and Fig. 6 show a comparison of the contours of the span-wise velocity and span-wise rms fluctuations respectively.
Figure 3: Stanford double diffuser, Re=10000, validation. Average velocities, Cherry et al. (2008) and Ohlsson et al. (2010) vs data obtained with Alya. |
Figure 4: Stanford double diffuser, Re=10000, validation. Average streamwise velocity fluctuations, Cherry et al. (2008) and Ohlsson et al. (2010) vs data obtained with Alya. |
Fig. 3 and Fig. 4 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.
Figure 5: Stanford double diffuser, Re=10000, validation. Average streamwise velocity contours, Cherry et al. (2008) and Ohlsson et al. (2010) vs data obtained with Alya. |
Figure 6: Stanford double diffuser, Re=10000, validation. Average streamwise velocity fluctuation contours, Cherry et al. (2008) and Ohlsson et al. (2010) vs data obtained with Alya. |
Other turbulent statistics
Other relevant statistics from the database are presented in this section. In particular, the turbulent kinetic energy (TKE) is presented in Fig. 7 along the characteristic lines of the diffuser.
Figure 7: Stanford double diffuser, Re=10000. Turbulent kinetic energy (TKE) obtained with Alya on characteristic lines of the diffuser. |
Then, the characteristic length scales (i.e., the Taylor micro-scale <math>{\eta_T}<\math> and the Kolmogorov length-scale <math>{\eta_K}<\math>) are presented in Fig. 8.
Figure 8: Stanford double diffuser, Re=10000. Taylor microscale and Kolmogorov lengthscale obtained with Alya on characteristic lines of the diffuser. |
Reynolds stresses budget equation terms
Contributed by: Oriol Lehmkuhl, Arnau Miro — Barcelona Supercomputing Center (BSC)
© copyright ERCOFTAC 2024