Hybrid Compact-WENO Finite Difference Scheme with Conjugate Fourier Shock Detection Algorithm for Hyperbolic Conservation Laws.
Wai-Sun Don,Zhen Gao,Peng Li,Xiao Wen
DOI: https://doi.org/10.1137/15m1021520
IF: 2.968
2016-01-01
SIAM Journal on Scientific Computing
Abstract:For discontinuous solutions of hyperbolic conservation laws, a hybrid scheme, based on the high order nonlinear characteristicwise weighted essentially nonoscillatory (WENO) conservative finite difference scheme and the high resolution spectral-like linear compact finite difference (Compact) scheme, is developed for capturing shocks and strong gradients accurately and resolving smooth scale structures efficiently. The key issue in any hybrid scheme is the design of an accurate, robust, and efficient high order shock detection algorithm that is capable of determining the smoothness of the solution at any given grid point. The conjugate Fourier (cF) partial sum and its derivative are investigated for its applicability as a shock detector due to its unique property, namely, the cF partial sum converges to the location and strength of an isolated jump. For a nonperiodic problem, the data are first evenly extended before the derivative of the cF partial sum and its mean are computed. The mean allows one to partition the domain into subdomains containing strong gradients or smooth solutions. The locations of shocks are then accurately identified and flagged for special treatment using the WENO scheme. The matrix-multiply, even-odd decomposition, and cosine/sine fast transform algorithms of the cF analysis are derived, and their advantages and disadvantages in their implementations, usage, and technical issues are discussed in detail. The cF shock detector and its iterative version for detecting jumps of large difference in scales are presented. The Hybrid-cF scheme is applied to one-dimensional shock-density wave interaction problem, two-dimensional (2D) Riemann initial value problems, and 2D Mach 10 double Mach reflection problems. The preliminary results are in good agreement with those obtained in the literature. They demonstrate the spatial and temporal adaptivity of the Hybrid-cF scheme for problems containing strong shocks, multiple developing shocklets, and high frequency waves. A speedup of the CPU times by a factor up to 2-3 is obtained showing the potential efficiency of the Hybrid-cF scheme over the pure WENO-Z scheme.