An efficient finite difference IFWENO-THINC hybrid scheme for capturing discontinuities

Abstract

The weighted essentially non-oscillatory (WENO) scheme is a famous shock-capturing scheme, but for high-precision numerical simulation, the classical WENO scheme is always considered to have excessive dissipation and insufficient resolution near discontinuities. Improving the accuracy order or reducing the dissipation of the WENO by optimizing the linear weights can effectively improve the resolution of this scheme in the smooth regions. However, under the same method, the improvement of shock-capturing abilities is generally limited and may lead to numerical oscillations. In order to simultaneously improve the resolution in discontinuous regions and ensure the performance of the smooth regions, this study uses the shock detector to combine the tangent of hyperbola for interface capturing (THINC) scheme, which is a better discontinuity-capturing scheme, with a given WENO discretization, such as the improved fast WENO (IFWENO) scheme that equips with high economy and does not loss accuracy at any critical point. The hybrid method in this study is an apriori algorithm, which is more efficient and universal than the method in (Takagi et al. J. Comput. Phys. 452(2022)110899). Additionally, the performance of the THINC scheme in different regions near the discontinuity is analyzed in detail through numerical experiments, and a shock detector that can fully exploit the advantages of the THINC scheme is provided. By using this shock detector, several one-dimensional and two-dimensional numerical cases are employed to verify that the hybrid scheme can greatly improve the performance in discontinuities without affecting the flow structure in the smooth regions, which proves the effectiveness of the proposed shock detector and the excellence of the hybrid IFWENO-THINC scheme.