[论文解读] Robust A Posteriori Error Estimates for Stabilized Finite Element Methods

Lutz Tobiska, R. Verfürth|arXiv (Cornell University)|Feb 24, 2014

Advanced Numerical Methods in Computational Mathematics参考文献 17被引用 3

一句话总结

本文为对对流-扩散方程应用的一类广泛稳定的有限元方法，建立了鲁棒的后验误差估计。结果表明，Streamline 污染、局部投影、子网格尺度和连续内罚方法均产生相同的误差指标，且常数与粘性、对流或反应项无关，确保了所有方案在统一可靠性与效率方面表现优异。

ABSTRACT

There is a wide range of stabilized finite element methods for stationary and non-stationary convection-diffusion equations such as streamline diffusion methods, local projection schemes, subgrid-scale techniques, and continuous interior penalty methods to name only a few. We show that all these schemes give rise to the same robust a posteriori error estimates, i.e. the multiplicative constants in the upper and lower bounds for the error are independent of the size of the convection or reaction relative to the diffusion. Thus, the same error indicator can be used modulo higher order terms caused by data approximation.

研究动机与目标

• 统一对对流-扩散问题中多种稳定有限元方法的后验误差估计。
• 证明不同稳定化方案导出的误差指标可产生等价的鲁棒误差界。
• 确保误差估计在对流主导区域中对粘性 ε 和反应参数 β 保持统一性。
• 推导出与对流或反应项相对大小无关的残差可计算上界与下界。
• 通过基于残差与一致性误差控制的统一框架，将分析扩展至定常与非定常问题。

提出的方法

• 采用 [39] 的通用框架，建立误差与残差之间鲁棒等价性的通用理论。
• 通过引理 2.3–2.6，为每种稳定化方案推导出一致性误差的显式、可计算上界。
• 采用针对双线性形式对称部分量身定制的能量范数与对偶范数，以确保鲁棒性。
• 采用基于残差的误差估计器，对不同方案引入局部贡献项 η、θ 与 Θ，包括连续内罚方法中的参数 σ_cip。
• 针对非定常问题，引入基于时间平均数据与离散时间步长的时间离散误差估计器。
• 在空间与时间上建立全局上界与局部下界，且常数与 ε 和 β 无关。

实验结果

研究问题

• RQ1 能否为多种稳定有限元方法在对流-扩散问题中应用单一鲁棒的后验误差估计器？
• RQ2 不同稳定化方案导出的误差估计是否在对流或反应项相对大小变化时仍保持统一的可靠性与效率？
• RQ3 每种稳定化方法引入的一致性误差能否以网格参数为变量，实现统一且显式的有界？
• RQ4 相同的误差指标是否适用于定常与非定常对流-扩散方程，并具备相同的鲁棒性特征？
• RQ5 能否使基于残差的误差估计器在所有稳定化方案中对粘性 ε 与反应参数 β 均保持鲁棒？

主要发现

• 所考虑的所有稳定有限元方法——包括 Streamline 污染、局部投影、子网格尺度与连续内罚方法——均产生相同的鲁棒后验误差估计。
• 误差的上界与下界对 ε 与 β 保持统一，即界中常数与对流或反应强度无关。
• 对于定常情形，残差范数被误差估计项 η、θ 与 Θ 上下界控制，常数 c♭ 与 c♭ 与 ε 和 β 无关。
• 对于非定常情形，全局误差被时间步长上估计量之和上界控制，常数 c* 与 c* 与 T、ε 和 β 无关。
• 时间离散误差估计器包含初始数据、时间演化与数据逼近的贡献，所有项均被统一有界。
• 连续内罚方法中参数 σ_cip 取值为 1，其余方法取值为 0，反映了各方案特异性的一致性误差贡献。