An Adaptive Flux Reconstruction Scheme for Robust Shock Capturing

Abstract: In the case of hyperbolic conservation laws, high-order methods, such as the classical DG method, experience the phenomenon of unwanted high-frequency oscillations in the vicinity of a shock. Shock-capturing methods such as artificial dissipation, solution, flux, or TVD limiting are generally used to eliminate non-physical oscillations and provide bounds on physical quantities. For entropy-stable schemes, the additional objective would be to retain provable entropy dissipation guarantees of the underlying scheme, i.e. subcell limiting or entropy filtering [1, 2, 3, 4]. The nonlinearly-stable flux reconstruction (NSFR) semi-discretization given in Eq. 7 with a suitable flux reconstruction scheme has been demonstrated to mitigate spurious oscillations in the presence of shock discontinuities and at CFL values substantially larger than the DG variant of the NSFR scheme whilst retaining the property of entropy stability [5]. NSFR schemes achieve this by introducing an alternative lifting operator for surface numerical flux penalization, albeit at the expense of accuracy. In this technical note, we present an adaptive approach to the choice of the lifting operator employed, which maintains higher accuracy and allows for larger CFL values while retaining the underlying provable attributes of the scheme. While it cannot eliminate oscillations such as the aforementioned shock-capturing methods, together with a positivity-preserving limiter, the scheme provides for solutions that are essentially oscillation-free.