Spectral Difference method with a posteriori limiting: Application to the Euler equations in one and two space dimensions

TL;DR

This paper introduces a high-order spectral difference numerical scheme with an a-posteriori limiting strategy that combines accuracy in smooth regions with robust shock capturing for Euler equations in 1D and 2D.

Contribution

The paper presents a novel combination of spectral difference and a-posteriori limiting using a MUSCL-Hancock fallback, enhancing accuracy and robustness in solving Euler equations.

Findings

Achieves very high order accuracy in smooth flow regions.

Effectively captures shocks without spurious oscillations.

Demonstrates robustness on 1D and 2D Euler test problems.

Abstract

We present a new numerical scheme which combines the Spectral Difference (SD) method up to arbitrary high order with \emph{a-posteriori} limiting using the classical MUSCL-Hancock scheme as fallback scheme. It delivers very accurate solutions in smooth regions of the flow, while capturing sharp discontinuities without spurious oscillations. We exploit the strict equivalence between the SD scheme and a Finite-Volume (FV) scheme based on the SD control volumes to enable a straightforward limiting strategy. At the end of each stage of our high-order time-integration ADER scheme, we check if the high-order solution is admissible under a number of numerical and physical criteria. If not, we replace the high-order fluxes of the troubled cells by fluxes from our robust second-order MUSCL fallback scheme. We apply our method to a suite of test problems for the 1D and 2D Euler equations. We demonstrate that this combination of SD and ADER provides a virtually arbitrary high order of accuracy, while at the same time preserving good sub-element shock capturing capabilities.