“…Originally developed by Huynh for 1D advection problems in [3], the FR method is a framework which allows to develop new high-order schemes while being simple and computationally efficient, particularly on graphical processor units (GPU) [4,5].Despite those advantages, the FR method suffers from the same pacing items as the other high-order methods such as the slow convergence to steady state and the lack of robust shock capturing methods [6]. When trying to approximate a discontinuity (such as a jump in a physical property due to a shock) with a polynomial representation, the Gibbs phenomenon causes a decrease in accuracy and leads to the appearance of spurious oscillations which can trigger instabilities in the numerical computations [7]. To solve this issue, different limiting and artificial viscosity (AV) strategies have been developed as summarized in [8].…”