The stability priority of spatial-temporal coupled compact element methods over decoupled compact element methods

TL;DR

This paper compares the stability of spatial-temporal coupled compact element methods with decoupled methods in CFD, showing coupled methods offer better stability and flow information utilization.

Contribution

It provides the first Fourier stability analysis comparing RK and LW methods, highlighting the advantages of spatial-temporal coupled schemes in CFD.

Findings

Coupled methods demonstrate superior stability over decoupled methods.

Decoupled methods require smaller time steps due to flow information loss.

Coupled methods effectively utilize initial-value information for better stability.

Abstract

With the increasing industrial demands, two families of high-order numerical schemes are widely used within the computational fluid dynamics community. One is the method of line, which relies on Runge-Kutta (RK) time-stepping applied to a semi-discrete, spatio-temporally decoupled formulation. The other is the family of Lax-Wendroff (LW) type method, which are inherently spatial-temporal coupled and are constructed within a multi-stage multi-derivative (MSMD) framework. This paper, for the first time, conducted a comparative Fourier stability analysis of RK and LW method to distinguish the dispersion and dissipation effects of numerical schemes respectively. Through rigorous theoretical derivation and consistent numerical validation, we draw the following conclusions: While explicit RK line methods are straightforward like Discontinuous Galerkin (DG) method and flux reconstruction (FR) method, they employ from a decoupling of spatial and temporal accuracy, thus discarding flow field evolution information and requiring small time steps. In contrast, spatial-temporal coupled compact methods, such as the gas-kinetic scheme (GKS) and the generalized Riemann problem (GRP) solver, utilize initial-value information from space far more effectively for time evolution. Even with just one additional order of spatial-temporal coupled information, they show better stability compared to RK methods. This provides new insights for CFD algorithm design, emphasizing the need for consistency between the dependence in the physical domain and that in numerical domain.