[论文解读] On Triangle Inequality Based Approximation Error Estimation

A. K. Alekseev, А. Е. Бондарев+1 more|arXiv (Cornell University)|Aug 15, 2017

Computational Fluid Dynamics and Aerodynamics参考文献 19被引用 9

一句话总结

该论文提出了一种基于三角不等式的数值方法，用于估计偏微分方程数值解的离散化误差，通过不同精度阶数的解之间的距离来界定真实解的范围。该方法应用于二维欧拉方程的定常超音速流问题中，成功地在度量定义的区域内包围了精确解，其性能高度依赖于范数的选择。

ABSTRACT

The distance between the true and numerical solutions in some metric is considered as the discretization error magnitude. If error magnitude ranging is known, the triangle inequality enables the estimation of the vicinity of the approximate solution that contains the exact one (exact solution enclosure). The analysis of distances between the numerical solutions enables discretization error ranging, if solutions errors are significantly different. Numerical tests conducted using the steady supersonic flows, governed by the two dimensional Euler equations, demonstrate the properties of the exact solution enclosure. The set of solutions generated by solvers of different orders of approximation is used. The success of this approach depends on the choice of metric.

研究动机与目标

• 解决在无法获得精确解的情况下，量化偏微分方程数值解中离散化误差的挑战。
• 开发一种仅利用不同精度的数值解来估计真实解位置的方法。
• 在存在多个解的情况下，通过三角不等式提供严格的误差包围。
• 研究不同范数（度量）在构建紧密误差边界中的有效性。
• 在基准流体动力学问题中展示该方法的可行性与准确性。

提出的方法

• 利用三角不等式将数值解之间的距离与真实解到任一近似解的距离关联起来。
• 在选定的度量空间中，将离散化误差的大小定义为真实解与数值解之间的距离。
• 使用不同精度阶数的求解器生成多个数值解，以估计误差范围。
• 应用三角不等式，将精确解包围在由近似解定义的区域内。
• 选择适当的范数（度量）以最小化误差包围区域的大小，提升精度。
• 通过定常超音速流问题的数值实验验证该方法。

实验结果

研究问题

• RQ1 三角不等式能否有效用于仅基于不同精度的数值解来估计真实解的位置？
• RQ2 度量（范数）的选择如何影响误差包围的紧密性和可靠性？
• RQ3 该方法在实际流体动力学模拟（如由二维欧拉方程控制的模拟）中的表现如何？
• RQ4 不同精度阶数的解之间的差异如何影响误差范围的估计？
• RQ5 该方法能否提供一种可靠且计算高效的传统误差估计技术的替代方案？

主要发现

• 只要数值解足够精确，三角不等式即可用于构建严格包含精确解的解包围区域。
• 在由二维欧拉方程控制的定常超音速流问题的数值实验中，该方法成功地将真实解包围在内。
• 误差包围区域的大小对范数的选择极为敏感，某些度量可产生显著更紧的边界。
• 数值测试证实，当存在不同精度阶数的解且其误差显著不同时，该方法是有效的。
• 该方法提供了一套实用的后验误差估计框架，无需依赖精确解的知识。
• 该方法的成功依赖于存在足够宽范围的解误差，从而能够实现可靠的基于距离的误差估计。