"A Posteriori" Limited High Order and Robust Residual Distribution Schemes for Transient Simulations of Fluid Flows in Gas Dynamics

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) (Clain et al., Journal of Computational Physics 2011; Diot et al., Computer & Fluids 2012) approach, a candidate solution is computed at a next time level via a high-order accurate explicit scheme (Abgrall, Journal of Scientific Computing 2017; Abgrall et al., Computers & Mathematics with Applications 2018). 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 Clain et al. (Journal of Computational Physics 2011), Diot et al. (Computer & Fluids 2012) and of Vilar (Journal of Computational Physics 2019), 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.