[论文解读] The stability priority of spatial-temporal coupled compact element methods over decoupled compact element methods

Qihui Gao, Xing Ji+4 more|arXiv (Cornell University)|Jan 24, 2026

Advanced Numerical Methods in Computational Mathematics被引用 0

一句话总结

论文比较时空耦合（Lax-Wendroff 型）与解耦（Runge-Kutta 型）紧凑单元方法的傅里叶稳定性，结果显示耦合方法通常具有更好的稳定性。它分析耗散和色散以指导CFD方案设计。

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.

研究动机与目标

• 通过考察时空耦合与解耦方法的稳定性来激发高阶CFD方案的研究兴趣。
• 制定两矩傅里叶稳定性分析，以量化 RK 与 LW 方案之间的耗散和色散差异。
• 表明时空耦合 enhances 稳定性并为CFD求解器设计提供指导。

提出的方法

• 在 GRP/GKS 与 DG/FR 框架内表述代表性 RK 型与 LW 型紧凑单元求解器。
• 使用模态傅里叶模态进行两矩傅里叶稳定性分析，以研究耗散和稳定性特性。
• 推导并比较 RK2 与 S1O2（单阶段、二阶）方案的两矩更新关系与放大因子。
• 在完全离散化的有限体积上下文中，利用相同的 CFL 设置评估稳定性。
• 通过复平面中的特征值分布讨论数值耗散的谱特性。

实验结果

研究问题

• RQ1 时空耦合的 LW 型与解耦的 RK 型紧凑单元方法在色散与耗散特性方面有何差异？
• RQ2 时空耦合是否提高稳定裕度并允许高阶CFD方案使用更大的时间步长？
• RQ3 在相似离散化下，CGKS/GRP/GKS 与 RK-DG/FR 方案的谱稳定特征（特征值）是什么？
• RQ4 在物理域依赖性与数值域传播之间，是否出现实现一致性的设计原则？
• RQ5 两矩（解与导数）形式是否比传统的一矩 RK 基方法在稳定性上更具优势？

主要发现

• 显式 RK 线性方法将时空解耦，放弃了一部分流动演变信息，强制较小的稳定时间步长。
• 时空耦合紧凑方法（如 GKS、GRP）更有效地利用初始值信息进行时间演化。
• 即使多获得一个数量级的时空耦合信息，耦合方法在可比设置下也显示出比 RK 方法更好的稳定性。
• 两矩方案的傅里叶分析表明耦合方法提供更好的稳定性特性和耗散性质。
• 稳定性洞见强调在CFD方法设计中物理域依赖性与数值域传播的一致性的重要性。