A dynamic multiscale viscosity algorithm for shock capturing in Runge Kutta Discontinuous Galerkin methods

Abstract

In order to improve the performance of higher-order Discontinuous Galerkin finite element solvers, a shock capturing procedure has been developed for hyperbolic equations. The Dynamic Multiscale Viscosity method, originally presented by Oberai and Wanderer [8, 9] in a Fourier Galerkin context, is adapted to the Discontinuous Galerkin discretization. The notions of diffusive model term, artificial viscosities, and the Germano identity are introduced. A general technique for the evaluation of the multiscale model term's parameters is then presented. This technique is used to perform efficient shock capturing on an one-dimensional stationary Burgers' equation with 1-parameter and 2-parameter model terms. Corresponding numerical results are shown.