“A posteriori” limited high order and robust schemes for transient simulations of fluid flows in gas dynamics

Abstract

In this paper, we propose a novel approximation strategy for time-dependent hyperbolic systems of conservation laws for the Euler system of gas dynamics that aims to represent the dynamics of strong interacting discontinuities. The goal of our method is to allow an approximation with a high-order of accuracy in smooth regions of the flow, while ensuring robustness and a non-oscillatory behaviour in the regions of steep gradients, in particular across shocks. Following the Multidimensional Optimal Order Detection (MOOD) ([15,17]) approach, a candidate solution is computed at a next time level via a high-order accurate explicit scheme ([3,5]). A so-called detector determines if the candidate solution reveals any spurious oscillations or numerical issue and, if so, only the troubled cells are locally recomputed via a more dissipative scheme. This allows to design a family of “a posteriori” limited, robust and positivity preserving, as well as high accurate, non-oscillatory and effective scheme. Among the detecting criteria of the novel MOOD strategy, two different approaches from literature, based on the work of [15,17] and on [35], are investigated. Numerical examples in 1D and 2D, on structured and unstructured meshes, are proposed to assess the effective order of accuracy for smooth flows, the non-oscillatory behaviour on shocked flows, the robustness and positivity preservation on more extreme flows.