PNNL

Abstract

Hierarchical reconstruction (HR) has been successfully applied to prevent oscillations in solutions computed by finite volume, discontinuous Galerkin, spectral volume schemes when solving hyperbolic conservation laws. In this paper, we demonstrate that HR can also be combined with spectral/hp element methods for solving hyperbolic conservation laws. We show that HR preserves the order of accuracy of spectral/hp element methods for smooth solutions and generate essentially non-oscillatory solution profiles for shock wave problems.