Difference between revisions of "Analysis of Step 2"

From CFD Benchmark
Jump to: navigation, search
Line 16: Line 16:
 
pseudo-spectral code RLPK using 5123 grid points <ref name="vanrees2011"/>.
 
pseudo-spectral code RLPK using 5123 grid points <ref name="vanrees2011"/>.
  
[[File:vx_3d_incomp.pdf]]
+
[[File:vx_3d_incomp.png]]
[[File:vy_3d_incomp.pdf]]
+
[[File:vy_3d_incomp.png]]
  
 
Velocities along the centerlines of the domain for cold 3-D TGV case ('''Step 2''') at <math>t = 12.11 \tau_{ref}</math>.
 
Velocities along the centerlines of the domain for cold 3-D TGV case ('''Step 2''') at <math>t = 12.11 \tau_{ref}</math>.
Line 32: Line 32:
 
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 5123 grid points with a pseudo-spectral solver'''.
 
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 5123 grid points with a pseudo-spectral solver'''.
  
[[File:ke_3d_incomp.pdf]]
+
[[File:ke_3d_incomp.png]]
[[File:eps_3d_incomp.pdf]]
+
[[File:eps_3d_incomp.png]]
  
 
Here is the evolution with time of turbulent kinetic energy k (left, written KE in vertical axis) and of its dissipation rate � (right), by comparison with published reference data from <ref name="vanrees2011"/>.
 
Here is the evolution with time of turbulent kinetic energy k (left, written KE in vertical axis) and of its dissipation rate � (right), by comparison with published reference data from <ref name="vanrees2011"/>.

Revision as of 20:15, 25 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 , as illustrated in Fig. 5. 2. The evolution of kinetic energy () and of its dissipation rate () versus time, as shown in Fig. 6, where is the symmetric strain-rate tensor,

 ; (11)

Comparing the velocity fields at 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 5123 grid points [1].

Vx 3d incomp.png Vy 3d incomp.png

Velocities along the centerlines of the domain for cold 3-D TGV case (Step 2) at .

First, it appears on Fig. 5 that no differences can be identified visually from the velocity fields along the centerlines at . 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 () for YALES2, and (2) just before flow relaminarization () 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 5123 grid points with a pseudo-spectral solver.

Ke 3d incomp.png Eps 3d incomp.png

Here is the evolution with time of turbulent kinetic energy k (left, written KE in vertical axis) and of its dissipation rate � (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