[论文解读] Spectral Analysis of Continuous FEM for Hyperbolic PDEs: Influence of Approximation, Stabilization, and Time-Stepping

Sixtine Michel, Davide Torlo+2 more|arXiv (Cornell University)|Jan 1, 2021

Advanced Numerical Methods in Computational Mathematics参考文献 49被引用 11

一句话总结

本文对一维双曲型PDE的连续有限元方法进行了全面的谱分析，评估了近似空间（伯恩斯坦、等距、高斯-勒让德-高斯勒布托）与稳定化技术（SUPG、CIP、LPS）以及时间积分格式（RK、SSPRK、DeC）之间的相互作用。主要贡献在于提出了一种完全离散的傅里叶分析，识别出最优的CFL数与稳定化参数，结论为：采用SSPRK时间积分与CIP或LPS稳定化的分段积分单元（cubature elements）在精度、稳定性与计算效率之间实现了最佳平衡。

ABSTRACT

We study continuous finite element dicretizations for one dimensional hyperbolic partial differential equations. The main contribution of the paper is to provide a fully discrete spectral analysis, which is used to suggest optimal values of the CFL number and of the stabilization parameters involved in different types of stabilization operators. In particular, we analyze the streamline-upwind Petrov–Galerkin stabilization technique, the continuous interior penalty (CIP) stabilization method and the orthogonal subscale stabilization (OSS). Three different choices for the continuous finite element space are compared: Bernstein polynomials, Lagrangian polynomials on equispaced nodes, and Lagrangian polynomials on Gauss-Lobatto cubature nodes. For the last choice, we only consider inexact quadrature based on the formulas corresponding to the degrees of freedom of the element, which allows to obtain a fully diagonal mass matrix. We also compare different time stepping strategies, namely Runge–Kutta (RK), strong stability preserving RK (SSPRK) and deferred correction time integration methods. The latter allows to alleviate the computational cost as the mass matrix inversion is replaced by the high order correction iterations. To understand the effects of these choices, both time-continuous and fully discrete Fourier analysis are performed. These allow to compare all the different combinations in terms of accuracy and stability, as well as to provide suggestions for optimal values discretization parameters involved. The results are thoroughly verified numerically both on linear and non-linear problems, and error-CPU time curves are provided. Our final conclusions suggest that cubature elements combined with SSPRK and CIP or OSS stabilization are the most promising combinations.

研究动机与目标

• 分析连续伽辽金FEM在一维双曲守恒律中，不同离散化选择下的稳定性和精度。
• 通过谱分析，确定CFL数与稳定化参数（SUPG、CIP、LPS）的最优值。
• 比较不同有限元空间（伯恩斯坦多项式、等距拉格朗日、高斯-勒让德-高斯勒布托分段积分单元）的性能，尤其关注质量矩阵结构与积分方法的影响。
• 评估时间积分格式：龙格-库塔（RK）、强稳定保模RK（SSPRK）与延迟校正（DeC），重点关注质量矩阵处理方式与计算成本。
• 在线性与非线性问题（包括浅水方程）上数值验证理论结果。

提出的方法

• 采用时间连续与完全离散的傅里叶分析（冯·诺依曼方法），评估所有方案组合下的相位误差与耗散误差。
• 分析三种稳定化方法：SUPG（迎风型Petrov-Galerkin）、CIP（连续内部惩罚）、LPS（局部投影稳定化），其对质量矩阵与一致性的影响各不相同。
• 比较三种有限元空间：伯恩斯坦多项式、等距节点上的拉格朗日多项式、高斯-勒让德-高斯勒布托节点（分段积分）上的拉格朗日多项式，后者通过非精确积分实现对角质量矩阵。
• 评估三种时间积分方法：标准RK、SSPRK（用于稳定性）、DeC（通过迭代校正避免质量矩阵求逆）。
• 采用谱分析框架，最小化相对误差 ηu，以确定每种方案的最优稳定化与CFL参数。
• 在对流方程、Burgers方程与浅水系统上进行数值验证，绘制误差-计算时间曲线以比较性能。

实验结果

研究问题

• RQ1 在连续FEM求解一维双曲型PDE时，CFL数与稳定化参数（SUPG、CIP、LPS）的最优值是什么，可确保稳定性和精度？
• RQ2 不同有限元空间（尤其是采用非精确积分的分段积分单元）如何影响方案的谱性质与计算成本？
• RQ3 时间积分方法的选择（RK、SSPRK、DeC）如何影响稳定性、精度与效率，特别是与质量矩阵求逆的关系？
• RQ4 完全离散格式的谱性质（相位与耗散误差）在多大程度上与线性与非线性问题上的实际数值收敛率相关？
• RQ5 在双曲问题中，哪种单元类型、稳定化方法与时间积分器的组合能在精度、稳定性与计算成本之间实现最佳权衡？

主要发现

• 分段积分单元（高斯-勒让德-高斯勒布托节点配合非精确积分）可生成完全对角的质量矩阵，实现无需矩阵求逆的高效求解，在谱精度与稳定性方面优于其他单元类型。
• 分段积分单元结合SSPRK时间积分与CIP或LPS稳定化，实现了精度、稳定性与计算效率的最佳平衡，该结论得到误差-计算时间曲线的验证。
• DeC时间积分通过用高阶校正迭代替代质量矩阵求逆，降低了计算成本，但引入了更高的相位与耗散误差，尤其在伯恩斯坦多项式与高阶单元中更为显著。
• 对于所有方案，最优CFL数约为SSPRK的2.0，高阶DeC的1.5，该值通过谱分析中最小化相对误差 ηu 确定。
• 谱分析表明，分段积分单元的相位误差低于等效质量矩阵方法（如DeC），且稳定性优于完整质量矩阵方法，尤其在高频模态中表现更优。
• 浅水方程的数值结果验证了理论预测：采用LPS与CIP稳定化的分段积分单元收敛率与预期阶数一致，而伯恩斯坦多项式在P2与P3阶配合DeC时表现出次优收敛。