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). • Stabilisation of DG for under-resolved turbulence through penalising gradient or solution jumps at element interfaces. • Implementation for 1D advection-diffusion and 3D Navier-Stokes equations. • Eigensolution analysis for dispersion-dissipation behaviour for various polynomial orders and Riemann fluxes. • Validation on Burgers and Taylor-Green vortex problems and comparison with classic explicit subgrid models. [ABSTRACT FROM AUTHOR]