Improvement of Stability and Accuracy for Weighted Essentially Nonoscillatory Scheme

Abstract

This paper studies the weights stability and accuracy of the implicit fifth-order weighted essentially nonoscillatory finite difference scheme. It is observed that the weights of the Jiang-Shu weighted essentially nonoscillatory scheme oscillate even for smooth flows. An increased e value of 10 -2 is suggested for the weighted essentially nonoscillatory smoothness factors, which removes the weights oscillation and significantly improves the accuracy of the weights and solution convergence. With the improved e value, the weights achieve the optimum value with minimum numerical dissipation in smooth regions and maintain the sensitivity to capture nonoscillatory shock profiles for the transonic flows. The theoretical justification of this treatment is given in the paper. The wall surface boundary condition uses a half-point mesh so that the conservative differencing can be enforced. A third-order accurate finite difference scheme is given to treat wall boundary conditions. The implicit time-marching method with unfactored Gauss-Seidel line relaxation is used with the high-order schemes to achieve a high convergence rate. Several transonic cases are calculated to demonstrate the robustness, efficiency, and accuracy of the methodology.