A new high-order RKDG method based on the TENO-THINC scheme for shock-capturing

In recent years, Runge-Kutta Discontinuous Galerkin (RKDG) methods have gained substantial attention in solving hyperbolic conservation laws, attributed to their high-order accuracy and adaptability to unstructured meshes. However, standard RKDG methods cannot capture discontinuities without oscillation unless they are supplemented with troubled cell indicators and limiters. Existing indicators, such as the total variation bounded (TVB) minmod indicator and the KXRCF indicator, typically depend on a critical parameter tied to the equation's solution, necessitating tuning for different cases. In terms of limiters, the popular limiters even fail to guarantee the high-order property in the smooth regions. The advent of Weighted Essentially Non-Oscillatory (WENO) schemes prompts the implementation of WENO limiters for RKDG methods, preserving high-order properties. However, WENO-family schemes exhibit significant numerical dissipation, potentially smearing small-scale flow structures. In this work, the Targeted Essentially Non-Oscillatory (TENO) indicator [1] is utilized, which leverages the nonlinear weighting strategy of the TENO scheme to separate high-wavenumber physical fluctuations and genuine discontinuities from smooth regions with a unified set of parameters. For troubled cells, a novel limiter is proposed for structured meshes, which combines a TENO scheme for resolving high-wavenumber physical fluctuations and a novel non-polynomial Tangent of Hyperbola for the INterface Capturing (THINC) scheme for resolving genuine discontinuities with extremely low numerical dissipation. Furthermore, the shifting between the TENO and THINC schemes is based on a new boundary variation diminishing (BVD) strategy, which only relies on compact neighborhoods and is significantly simpler than its predecessors. Meanwhile, a new strategy is proposed to ensure the consistency of the new limiter applied for 1D and 2D cases. A set of 1D and 2D benchmark cases including strong shockwaves and a broad range of flow length scales is simulated with uniform meshes for 1D cases and structured quadrilateral meshes for 2D cases to demonstrate the performance of the new numerical scheme. The indicator does not activate any limiters in the accuracy test cases to ensure the high-order property of the whole numerical scheme. Other cases with discontinuities demonstrate the low-dissipation properties of the new limiter in comparison to the simple WENO limiter.