High-Resolution Weighted Essentially Non-Oscillatory Compact Least-Squares Schemes with Implicit Time Integration for Compressible Navier-Stokes Equations on Curvilinear Grids
This work offers an incremental improvement for computational fluid dynamics practitioners seeking more accurate and efficient shock-capturing schemes on complex grids.
The paper introduces improved weighted essentially non-oscillatory compact least-squares schemes for compressible Navier-Stokes equations on curvilinear grids, addressing high-frequency oscillations and efficiency. The proposed method demonstrates robust shock-capturing and high resolution in smooth regions across 1D, 2D, and 3D inviscid and viscous flows.
This paper presents a family of high-resolution weighted essentially non-oscillatory compact least-squares schemes with implicit time integration for the compressible Navier-Stokes equations on curvilinear grids. Compared with the original compact least-squares schemes, the proposed method introduces two main improvements. First, instead of enforcing the accuracy constraints over the entire computational domain, including discontinuous regions, it constructs the reconstruction matrix only along smooth reconstruction lines using an accuracy-preserving weighting strategy. This treatment effectively suppresses the persistent high-frequency oscillations observed in the original compact least-squares schemes and yields sharper profiles near discontinuities. Second, the method simplifies the shock-capturing procedure and improves efficiency by using a single set of polynomials, whereas the original compact least-squares schemes require both unlimited and limited polynomials. Combined with spectral optimization, the proposed method exhibits more favorable spectral properties than conventional weighted essentially non-oscillatory schemes. The smoothness indicators and penalty matrices are constructed through an efficient iterative procedure with modest additional cost. Numerical results for inviscid and viscous one-, two-, and three-dimensional flows demonstrate that the proposed method provides robust shock-capturing capability while maintaining high resolution in smooth regions and across contact discontinuities.