An Adaptive Flux Reconstruction Scheme for Robust Shock Capturing

TL;DR

This paper introduces an adaptive flux reconstruction scheme for hyperbolic conservation laws that improves shock capturing by balancing accuracy, stability, and larger CFL values, reducing non-physical oscillations.

Contribution

It presents an adaptive approach to the lifting operator in NSFR schemes, enhancing accuracy and stability while maintaining entropy stability and shock-capturing capabilities.

Findings

Achieves larger CFL values than traditional NSFR schemes.

Maintains provable entropy stability with adaptive lifting.

Provides oscillation-free solutions with positivity-preserving limiter.

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.