Parallel implicit DNS of temporally-evolving turbulent shear layer instability

In this study, temporally-evolving incompressible and compressible Turbulent Shear Layer (TSL) instability problem is solved using an all-speed (allMach), implicit, non-dissipative and kinetic energy conserving algorithm. An in-house, fully parallel, finite-volume Direct Numerical Simulation (DNS) solver was developed using PETSc. Convergence characteristics at low-Mach numbers were also improved using a relaxation procedure. We aim here to assess the performance and behavior of the present algorithm for complex flows which contain multiscale physics and gradually evolve into turbulence. The results show that the algorithm is able to produce correct physical mechanisms and capture the evolution of the turbulent fluctuations for both incompressible and compressible cases. It is observed that the non-dissipative and kinetic energy conserving properties make the algorithm powerful and applicable to challenging problems. For higher Mach numbers, a shock-capturing or a dissipative mechanism is required for robustness.