[论文解读] A posteriori error estimates for space-time IgA approximations to parabolic initial boundary value problems

Ulrich Langer, Svetlana Matculevich+1 more|arXiv (Cornell University)|Dec 28, 2016

Advanced Numerical Methods in Computational Mathematics参考文献 41被引用 5

一句话总结

本论文利用功能后验估计技术，为抛物型初边值问题的空间-时间等几何分析（IgA）逼近方法推导了后验误差估计。这些估计不依赖于网格的常数，为能量误差范数提供了可靠且高效的界，并通过自由参数实现灵活性，可针对不同误差范数进行优化。

ABSTRACT

This work is concerned with a posteriori error estimates of the functional type for approximations constructed by space-time IgA scheme presented in paper by Langer, Neumueller, and Moore (2016). We consider approxima- tions in the corresponding IgA spaces based on elliptic and bounded bilinear form (associated with the spatial part). It is proved that the approximations satisfy classic a priori error estimates. Also, we deduce a posteriori error estimates for a stabilized weak formulation of the considered parabolic initial boundary value problem (I-BVP). They are derived by a general functional method and do not contain mesh dependent constants. The estimates are valid for a wide class of approximations. In particular, they imply estimates for the discrete norm of IgA approximations. Moreover, we introduce different forms of a posteriori error estimates (error majorants) and establish equivalence of majorants and energy error norm. This property justifies efficiency and reliability of a posteriori error estimates. Another important property of the estimates is their flexibility with respect to a certain amount of free parameters. Using these parameters we can obtain estimates for different error norms and minimize the respective majorant in order to find the best possible bound of the error.

研究动机与目标

• 推导抛物型初边值问题空间-时间IgA逼近方法的后验误差估计。
• 确保估计不依赖于网格相关常数，并适用于能量类中的任意逼近。
• 建立误差主要估计量与能量误差范数之间的等价性，以保证可靠性和效率。
• 提供一个具有自由参数的灵活框架，以最小化误差界并针对特定误差范数。
• 将功能后验估计方法的适用性扩展至空间-时间IgA格式，以支持自适应计算与并行计算。

提出的方法

• 利用泛函分析与对偶论证推导后验误差主要估计量，避免使用依赖于网格的常数。
• 采用基于B样条、NURBS或局部样条的高光滑度空间-时间IgA逼近方法，以实现高精度。
• 构建两种形式的误差主要估计量 $ \overline{\rm M}^{\rm I} $ 和 $ \overline{\rm M}^{\rm II} $，其中包含可调自由参数 $ \gamma, \alpha_i, \rho_i, \beta_i, \epsilon $。
• 在适当的函数空间上，建立能量误差范数与主要估计量下确界的等价性。
• 利用变分恒等式与柯西-施瓦茨型估计，将误差以残差和相容性项的形式进行界定。
• 证明主要估计量是可计算的，并可通过在参数上最小化来获得最紧致的误差界。

实验结果

研究问题

• RQ1 能否为抛物问题的空间-时间IgA逼近方法推导出可靠且高效的后验误差估计？
• RQ2 如何使估计独立于网格相关常数，同时仍适用于IgA逼近？
• RQ3 所提出的误差主要估计量与真实能量误差范数之间存在何种关系？
• RQ4 主要估计量中的自由参数如何影响误差界的质量及其最小化？
• RQ5 该框架是否可通过参数调节用于推导不同误差范数的界？

主要发现

• 所提出的后验误差主要估计量与能量误差范数等价，从而确保了可靠性和效率。
• 等价常数 $ C^{ m I}_{ m eq} $ 和 $ C^{ m II}_{ m eq} $ 可显式计算，且仅依赖于自由参数。
• 该估计适用于能量类中的任意逼近，无需满足符合性或依赖于网格的假设。
• 主要估计量 $ \overline{\rm M}^{\rm I} $ 和 $ \overline{\rm M}^{\rm II} $ 可通过对参数进行最小化，以获得能量误差的最紧致上界。
• 该方法通过提供可计算且可调参数的误差指示器，支持自适应加密与并行计算。
• 通过调节参数 $ \gamma, \rho_1, \rho_2, \epsilon $ 等，该框架可支持不同误差范数，实现针对性的误差估计。