[论文解读] A self adjusting multirate algorithm based on the TR-BDF2 method

Luca Bonaventura, Francesco Casella+2 more|arXiv (Cornell University)|Jan 27, 2018

Numerical methods for differential equations参考文献 13被引用 3

一句话总结

本文提出一种基于TR-BDF2求解器的自适应多速率时间积分方法，用于刚性常微分方程（ODEs）和双曲型偏微分方程（PDEs）。通过基于误差估计器对各组件动态调整时间步长，并利用TR-BDF2的嵌入式误差控制与三次样条插值，该方法在保持刚性与双曲型问题的精度与稳定性的同时，显著提升了计算效率——尤其适用于具有局部动力学（如激波）的问题。

ABSTRACT

We propose a self adjusting multirate method based on the TR-BDF2 solver. The potential advantages of using TR-BDF2 as the key component of a multirate framework are highlighted. A linear stability analysis of the resulting approach is presented and the stability features of the resulting algorithm are analysed. The analysis framework is completely general and allows to study along the same lines the stability of self adjusting multirate methods based on a generic one step solver. A number of numerical experiments demonstrate the efficiency and accuracy of the resulting approach also the time discretization of hyperbolic partial differential equations.

研究动机与目标

• 解决具有广泛时间尺度差异的系统中单速率时间积分的计算低效问题。
• 克服现有多速率方法在刚性与双曲型问题中的稳定性与效率限制。
• 利用TR-BDF2方法的优良特性——L-稳定性、嵌入式误差估计与连续输出——实现鲁棒的多速率积分。
• 通过自适应时间步长细化实现对活跃（快速）与潜在（慢速）组件的自动识别。
• 在基准问题上验证该方法的有效性，包括具有激波与膨胀波的非线性双曲型PDE。

提出的方法

• 采用自适应多速率策略：先以TR-BDF2进行全局时间步长的初步计算，仅在局部误差容限超出时对特定组件进行细化。
• 选用TR-BDF2作为基础单速率积分器，因其具备二阶精度、L-稳定性及嵌入式三阶误差估计器。
• 基于局部误差估计设计划分准则，将组件分类为活跃（需小步长）或潜在（可使用大步长）。
• 采用三次样条插值实现连续输出，以在时间步长细化过程中保持精度，尤其对刚性与瞬态现象至关重要。
• 实施递归时间步长策略，通过仅在必要区域进行局部细化，避免全局刚性约束。
• 采用适用于任意单步法的线性稳定性分析框架，并在此处应用于多速率环境下的TR-BDF2。

实验结果

研究问题

• RQ1 基于TR-BDF2的自适应多速率方法是否能在刚性与双曲型问题上相比单速率求解器实现更高的稳定性与效率？
• RQ2 在不同耦合强度与刚性条件下，特别是针对具有纯虚特征值的问题时，该多速率TR-BDF2方法的稳定性表现如何？
• RQ3 与线性插值相比，使用三次样条插值在多速率格式中在提升精度与减少步长拒绝方面有多大改善？
• RQ4 该方法能否在最小用户干预下自动检测并解析双曲型PDE中的局部现象（如激波与膨胀波）？
• RQ5 在应用该方法于基准非线性双曲型PDE（如Burgers方程）时，其计算节省与精度之间的权衡如何？

主要发现

• 该自适应多速率TR-BDF2方法通过减少潜在组件的时间步长数量，在具有局部动力学（如激波）的问题中实现了显著的计算效率提升。
• 数值实验表明，三次样条插值对性能至关重要：使用线性插值导致步长拒绝过多，且效率无法超越单速率TR-BDF2。
• 采用三次样条插值的多速率TR-BDF2、线性插值版本及单速率TR-BDF2的相对误差相近，表明各方法在精度上保持一致。
• 与高精度参考解（使用MATLAB的ode45求得）相比，所有变体的误差均降低数个数量级，证实空间与时间离散化误差得到良好控制。
• 即使在具有纯虚特征值的问题中，该方法也表现出良好的稳定性特性——表明其在双曲型PDE与结构动力学中具有潜在应用价值，尽管尚未提供完整的理论稳定性证明。
• 在膨胀波情形中，活跃组件的比例更高，计算节省效果弱于激波情形，反映出解的活动范围更广。