[论文解读] Propagation of internal errors in explicit Runge--Kutta methods and internal stability of SSP and extrapolation methods

David I. Ketcheson, L Lajos+1 more|King Abdullah University of Science and Technology Repository (King Abdullah University of Science and Technology)|Sep 5, 2013

Numerical methods for differential equations参考文献 27被引用 5

一句话总结

本文分析了显式Runge--Kutta（ERK）方法中内部误差（如舍入误差和求解器不准确性）的传播，表明实现形式对误差放大具有关键影响。文章引入了最大内部放大因子作为关键度量，并推导了SSP和外推法的边界，揭示了高阶数方法可能遭受严重的舍入误差放大，从而破坏高精度目标。

ABSTRACT

In practical computation with Runge--Kutta methods, the stage equations are not satisfied exactly, due to roundoff errors, algebraic solver errors, and so forth. We show by example that propagation of such errors within a single step can have catastrophic effects for otherwise practical and well-known methods. We perform a general analysis of internal error propagation, emphasizing that it depends significantly on how the method is implemented. We show that for a fixed method, essentially any set of internal stability polynomials can be obtained by modifying the implementation details. We provide bounds on the internal error amplification constants for some classes of methods with many stages, including strong stability preserving methods and extrapolation methods. These results are used to prove error bounds in the presence of roundoff or other internal errors.

研究动机与目标

• 研究显式Runge--Kutta方法单步内内部误差（如舍入误差和代数求解器不准确性）的放大机制。
• 证明内部误差传播在很大程度上取决于实现形式，而不仅仅是方法系数，从而挑战了代数等价方法在扰动下行为一致的假设。
• 建立一个通用框架，用于计算内部稳定多项式并量化高阶数方法的误差放大。
• 为强稳定保持（SSP）和外推法提供内部误差放大的边界，尤其关注高阶、高精度计算场景。
• 确立最大内部放大因子作为实用度量，以评估内部稳定性是否可能在高精度计算中损害数值精度。

提出的方法

• 基于ERK方法的实现形式，推导内部稳定多项式的通用表达式，区分Butcher形式与Shu--Osher形式。
• 引入最大内部放大因子作为衡量单步内内部误差可被放大的程度的指标，类似于条件数。
• 分析扰动方程（5）中的误差传播，该方程模拟了实际实现中由舍入和迭代求解器残差引起的误差。
• 将该框架应用于特定类别：强稳定保持（SSP）方法和外推法（EM），推导解析和数值边界。
• 利用复分析和稳定性区域上的优化，计算或估计不同阶数和阶段数下的最大放大因子。
• 使用符号计算（如Mathematica）推导特定情况下的精确代数次数和有理边界，特别是针对p=2,4,6,8阶的EM方法。

实验结果

研究问题

• RQ1 实现形式的选择（如Butcher与Shu--Osher形式）如何影响显式Runge--Kutta方法中内部误差的传播？
• RQ2 强稳定保持（SSP）和外推法的最大内部放大因子是多少？其随阶数和阶段数如何变化？
• RQ3 内部误差放大能否独立于截断误差进行边界控制？在何种条件下其会在高精度计算中占主导地位？
• RQ4 在不改变方法阶数或稳定性特性的情况下，通过修改实现细节在多大程度上可改善内部稳定性？
• RQ5 对于高阶外推法，内部误差放大因子在复平面稳定性区域（特别是左半平面）内如何表现？

主要发现

• 在高阶段数的显式Runge--Kutta方法中，即使方法本身条件良好且稳定，内部误差放大也可能造成灾难性后果。
• 二阶显式中点外推法的最大内部放大因子约为2.198，随阶数迅速增长：p=8时约为88.755，p=20时约为1.11×10⁶。
• 对于外推法，左半平面稳定性区域内的最大放大因子渐近增长为~1.1524·1.509^p / p，表明其随阶数呈指数增长。
• 当限制在z=0时，p=20阶EM方法的放大因子达到约199.9，表明即使在原点，高阶格式的放大效应也可能非常显著。
• 通过改变实现细节，即使对于代数等价的方法，内部稳定多项式（以及由此产生的误差放大）也可被改变，证明实现选择至关重要。
• 本研究揭示，忽略舍入误差的传统误差边界在高精度设置下不再有效，此时内部误差放大可能占主导地位并破坏计算精度。