Jump penalty stabilization techniques for under-resolved turbulence in discontinuous Galerkin schemes

Jump penalty stabilization techniques for under-resolved turbulence have been recently proposed for continuous and discontinuous high order Galerkin schemes [1–3]. The stabilization relies on the gradient or solution discontinuity at element interfaces to incorporate localised numerical diffusion in the numerical scheme. This diffusion acts as an implicit subgrid model and stabilizes under-resolved turbulent simulations. This paper investigates the effect of jump penalty stabilization methods (penalising gradient or solution) for stabilization and improvement of high-order discontinuous Galerkin schemes in turbulent regime. We analyze these schemes using an eigensolution analysis, a 1D non-linear Burgers equation (mimicking a turbulent cascade) and 3D turbulent Navier-Stokes simulations (Taylor-Green Vortex and Kelvin-Helmholtz Instability problems). We show that the two jump penalty stabilization techniques can stabilize under-resolved simulations thanks to the improved dispersion-dissipation characteristics (when compared to non-penalized schemes) and provide accurate results for turbulent flows. The numerical results indicate that the proposed jump penalty methods stabilize under-resolved simulations and improve the simulations, when compared to the original unpenalized scheme and to classic explicit subgrid models (Smagorinsky and Vreman).