Latest revision as of 14:32, 26 August 2020

The second step concerns the 3-D, non-reacting cold flow, used as validation by comparison with the published results of a pseudo-spectral solver ^[1]. The main quantities of interest for this comparison are:

1. Velocity profiles along centerlines of the domain at ${\displaystyle t=12.11\tau _{ref}}$, as illustrated in Fig. 5. 2. The evolution of kinetic energy (${\displaystyle k(t)={\frac {1}{2}}<u_{i}u_{i}>}$) and of its dissipation rate (${\displaystyle \epsilon (t)=dk/dt=2\nu <S_{ij}S_{ij}>}$) versus time, as shown in Fig. 6, where ${\displaystyle S_{ij}}$ is the symmetric strain-rate tensor,

${\displaystyle S_{ij}={\frac {1}{2}}({\frac {\partial u_{i}}{\partial x_{j}}}+{\frac {\delta u_{j}}{\delta x_{i}}})}$ ; (11)

Comparing the velocity fields at ${\displaystyle t=12.11\tau _{ref}}$ is done on purpose. As known from the TGV literature, this instant corresponds to a complex pseudo-turbulent field, before further turbulence decay due to dissipation (see also Figure 9 in ^[2]). Getting the correct velocity field in these conditions is challenging, since the obtained results are very sensitive with regard to the employed algorithms and discretization. Unlike the 2-D situation, no analytical solution is available there, and the results can only be compared to other numerical simulations. In the present case, the reference data are taken from a simulation relying on the pseudo-spectral code RLPK using ${\displaystyle 512^{3}}$ grid points ^[1].

Velocities along the centerlines of the domain for cold 3-D TGV case (Step 2) at ${\displaystyle t=12.11\tau _{ref}}$.

First, it appears on Fig. 5 that no differences can be identified visually from the velocity fields along the centerlines at ${\displaystyle t=12.11\tau _{ref}}$. Looking at Figs. 6 it can be observed that the three codes are able to reproduce the evolution of turbulence kinetic energy without any visible differences, whereas for the dissipation rate minute deviations appear at two instants (see enlargements in Fig. 6, right): (1) shortly after transition (${\displaystyle 11<t/\tau _{ref}<13.5}$) for YALES2, and (2) just before flow relaminarization (${\displaystyle 16.25<t/\tau _{ref}20}$) for both DINO and YALES2. The results of RLPK and of Nek5000 coincide visually at all times. These two time-instants are very sensitive moments at which the accuracy of the numerical methods and the resolution in time and space appear to play a major role. These small discrepancies are regarded as minor and considered acceptable with respect to the validation process of the codes. It must also be kept in mind that the data used as a reference have been obtained with a resolution of ${\displaystyle 512^{3}}$ grid points with a pseudo-spectral solver.

Here is the evolution with time of turbulent kinetic energy k (left, written KE in vertical axis) and of its dissipation rate ${\displaystyle \epsilon }$ (right), by comparison with published reference data from ^[1].

References

1. ↑ ^{1.0 1.1 1.2} W.M. van Rees, A. Leonard, D.I. Pullin,, P. Koumoutsakos, A comparison of vortex and pseudo-spectral methods for the simulation of periodic vortical flows at high Reynolds numbers., J. Comput. Phys. 230:2794--2805, 2011, Bibtex
   Author : W.M. van Rees, A. Leonard, D.I. Pullin,, P. Koumoutsakos
   Title : A comparison of vortex and pseudo-spectral methods for the simulation of periodic vortical flows at high Reynolds numbers.
   In : J. Comput. Phys. -
   Address :
   Date : 2011
   
2. ↑ A. Abdelsamie, G. Fru, F. Dietzsch, G. Janiga, D. Thévenin, Towards direct numerical simulations of low-Mach number turbulent reacting and two-phase flows using immersed boundaries, Comput. Fluids 131(5):123--141, 2016, Bibtex
   Author : A. Abdelsamie, G. Fru, F. Dietzsch, G. Janiga, D. Thévenin
   Title : Towards direct numerical simulations of low-Mach number turbulent reacting and two-phase flows using immersed boundaries
   In : Comput. Fluids -
   Address :
   Date : 2016