Development of nonlinear weighted compact schemes with increasingly higher order accuracy

  • Authors:

  • Shuhai Zhang;Shufen Jiang;Chi-Wang Shu

  • Affiliations:

  • China Aerodynamics Research and Development Center, Mianyang, Sichuan 621000, China;China Aerodynamics Research and Development Center, Mianyang, Sichuan 621000, China;Division of Applied Mathematics, Brown University, Providence, RI 02912, USA

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2008

Quantified Score

Hi-index 31.48

Visualization

Abstract

In this paper, we design a class of high order accurate nonlinear weighted compact schemes that are higher order extensions of the nonlinear weighted compact schemes proposed by Deng and Zhang [X. Deng, H. Zhang, Developing high-order weighted compact nonlinear schemes, J. Comput. Phys. 165 (2000) 22-44] and the weighted essentially non-oscillatory schemes of Jiang and Shu [G.-S. Jiang, C.-W. Shu, Efficient implementation of weighted ENO schemes, J. Comput. Phys. 126 (1996) 202-228] and Balsara and Shu [D.S. Balsara, C.-W. Shu, Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy, J. Comput. Phys. 160 (2000) 405-452]. These nonlinear weighted compact schemes are proposed based on the cell-centered compact scheme of Lele [S.K. Lele, Compact finite difference schemes with spectral-like resolution, J. Comput. Phys. 103 (1992) 16-42]. Instead of performing the nonlinear interpolation on the conservative variables as in Deng and Zhang (2000), we propose to directly interpolate the flux on its stencil. Using the Lax-Friedrichs flux splitting and characteristic-wise projection, the resulted interpolation formulae are similar to those of the regular WENO schemes. Hence, the detailed analysis and even many pieces of the code can be directly copied from those of the regular WENO schemes. Through systematic test and comparison with the regular WENO schemes, we observe that the nonlinear weighted compact schemes have the same ability to capture strong discontinuities, while the resolution of short waves is improved and numerical dissipation is reduced.