Shock Capturing via Limiting for High-Order Methods Including Discontinuous Galerkin

As is well-known, popular high-order methods such as discontinuous Galerkin (DG) and flux reconstruction (FR) tend to generate oscillations near shocks and discontinuities, which can lead to negative pressure or density, ultimately causing code breakdown.

For standard second-order methods, one approach to mitigating oscillations is to impose constraints on the calculated slopes, ensuring they do not become excessively steep. This limiting process is rooted in the idea of preserving monotonicity introduced by Van Leer (1974): when the data are monotone, limiting the slopes result in a monotone piecewise linear solution. However, a drawback of this approach is the loss of accuracy near extrema, where the non-monotone solution is somewhat flattened and achieves only first-order accuracy.