Abstract:The shock instability problem commonly arises in flow simulations involving strong shocks, particularly when employing high-order schemes, limiting their applications in hypersonic flow simulations. This study focuses on exploring the numerical characteristics and underlying mechanisms of shock instabilities in fifth-order finite-volume WENO schemes. To this end, for the first time, we have established the matrix stability analysis method for the fifth-order scheme. By predicting the evolution of perturbation errors in the exponential growth stage, this method provides quantitative insights into the behavior of shock-capturing and helps elucidate the mechanisms that cause shock instabilities. Results reveal that even dissipative solvers also suffer from shock instabilities when the spatial accuracy is increased to fifth-order. Further investigation indicates that this is due to the excessively high spatial accuracy of the WENO scheme near the numerical shock structure. Moreover, the shock instability problem of fifth-order schemes is demonstrated to be a multidimensional coupling problem. To stably capture strong shocks, it is crucial to have sufficient dissipation on transverse faces and ensure at least two points within the numerical shock structure in the direction perpendicular to the shock. The source location of instability is also clarified by the matrix stability analysis method, revealing that the instability arises from the numerical shock structure. Additionally, stability analysis demonstrates that local characteristic decomposition helps mitigate shock instabilities in high-order schemes, although the instability still persists. These conclusions pave the way for a better understanding of the shock instability in fifth-order schemes and provide guidance for the development of more reliable high-order shock-capturing methods for compressible flows with high Mach numbers.

A robust shock-capturing scheme based on rotated Riemann solvers

Shock-Stable Roe Scheme Combining Entropy Fix and Rotated Riemann Solver

Robust and accurate Roe-type Riemann solver with compact stencil: Rotated-RoeM scheme

Shock Capturing Scheme Based on Two- Shock Approximation to Riemann Problem

Mechanism-Derived Shock Instability Elimination for Riemann-Solver-Based Shock-Capturing Scheme

An Improved Meshfree Scheme Based on Radial Basis Functions for Solving Incompressible Navier–Stokes Equations

A Second-Order Rotational Upwind Transport Scheme for Solving Multi-Dimensional Compressible Euler Equations

The Role of Momentum Interpolation Mechanism of the Roe Scheme in the Shock Instability

A robustness-enhanced method for Riemann solver

Discontinuity-resolving shock-capturing schemes on unstructured grids

Role of the Momentum Interpolation Mechanism of the Roe Scheme in Shock Instability

General Procedure for Riemann Solver to Eliminate Carbuncle and Shock Instability

A shock-stable numerical scheme accurate for contact discontinuities: Applications to 3D compressible flows

Mechanism Study of Shock Instability in Riemann-Solver-Based Shock-Capturing Scheme

A Second-order Hybrid Rotated Riemann Solver for Multi-dimensional Euler Equations on Unstructured Meshes

Hybrid Compact-WENO Finite Difference Scheme with Radial Basis Function Based Shock Detection Method for Hyperbolic Conservation Laws.

A rotated lattice Boltzmann flux solver with improved stability for the simulation of compressible flows with intense shock waves at high Mach number

Shock-capturing particle hydrodynamics with reproducing kernels

Numerical stability analysis of shock-capturing methods for strong shocks II: high-order finite-volume schemes

Modifying an Unstructured Roe Solver for Large Eddy Simulation

A Robust WENO Type Finite Volume Solver for Steady Euler Equations on Unstructured Grids