[论文解读] High-order generalized-$\alpha$ methods

Quanling Deng, Pouria Behnoudfar+1 more|arXiv (Cornell University)|Feb 14, 2019

Numerical methods for differential equations参考文献 8被引用 4

一句话总结

该论文提出了一种三阶精度的广义-$\alpha$时间积分格式，保持无条件稳定性和用户可控的高频耗散特性，通过在离散格式中采用更高阶的泰勒展开，将二阶广义-$\alpha$方法推广至三阶精度。该方法在保持原始格式单步、无条件稳定结构的同时实现了三阶精度，通过参数$\alpha_m$和$\alpha_f$控制数值耗散。

ABSTRACT

The generalized-$\\alpha$ method encompasses a wide range of time integrators. The method possesses high-frequency dissipation while minimizing unwanted low-frequency dissipation and the numerical dissipation can be controlled by the user. The method is unconditionally stable and is of second-order accuracy in time. We extend the second-order generalized-$\\alpha$ method to third-order in time while the numerical dissipation can be controlled in a similar fashion. We establish that the third-order method is unconditionally stable. We discuss a possible path to the generalization to higher order schemes. All these high-order schemes can be easily implemented into programs that already contain the second-order generalized-$\\alpha$ method.

研究动机与目标

• 将二阶广义-$\alpha$方法推广至时间方向的三阶精度，同时保持其关键特性。
• 在新型高阶格式中保持无条件稳定性。
• 保留广义-$\alpha$方法的标志性特征——对高频数值耗散的用户可控性。
• 实现与现有二阶广义-$\alpha$代码库的无缝集成。
• 为推广至更高阶格式提供理论框架。

提出的方法

• 通过将$U_{n}^{\alpha_f}$和$V_{n}^{\alpha_m}$视为更高阶泰勒展开而非线性子步骤，重新表述广义-$\alpha$方法。
• 引入额外参数$\alpha_f$和$\alpha_m$，以表示$U_{n+\alpha_f}$和$V_{n+\alpha_m}$泰勒展开中的高阶项，从而实现三阶精度。
• 利用泰勒级数分析推导时间推进格式，以确定自由参数$\alpha_m$和$\alpha_f$的最优取值以实现最佳精度。
• 构造放大矩阵并分析其特征值，以确定无条件稳定区域。
• 通过将高频特征值$\rho_\infty$设置为$[0,1]$区间内的用户指定值来控制数值耗散。
• 采用符号分析与数值分析相结合的方法，验证无条件稳定性，并确定稳定性和期望耗散特性的参数区域。

实验结果

研究问题

• RQ1 广义-$\alpha$方法能否在不牺牲无条件稳定性的前提下推广至时间方向的三阶精度？
• RQ2 在更高阶时间积分格式中，如何保持对高频数值耗散的用户可控性？
• RQ3 哪些参数选择$\alpha_m$和$\alpha_f$能够实现三阶精度与无条件稳定性？
• RQ4 所得到的格式能否通过最小改动集成到现有的二阶广义-$\alpha$代码中？
• RQ5 三阶广义-$\alpha$方法的稳定区域是什么？与二阶版本相比有何差异？

主要发现

• 所提出的三阶广义-$\alpha$方法通过对方程解和速度项进行更高阶泰勒展开，实现了时间方向的正式三阶精度。
• 当$\alpha_m \geq \frac{7}{12}$且$\frac{1}{2} \leq \alpha_f \leq \alpha_m - \frac{1}{12}$时，该方法对所有参数组合均保持无条件稳定，其稳定区域通过符号与数值方法共同验证。
• 高频下的数值耗散由参数$\rho_\infty$控制，其中$\alpha_m = \frac{13 + 20\rho_\infty - 5\rho_\infty^2}{12(\rho_\infty + 1)^2}$且$\alpha_f = \frac{1 + 3\rho_\infty}{2(\rho_\infty + 1)^2}$，适用于$0 \leq \rho_\infty \leq 1$。
• 放大矩阵的特征值模长在无条件稳定区域内始终不超过1，证实了方法的稳定性。
• 数值结果表明，三阶格式的无条件稳定区域大于二阶版本，如数值图所示。
• 该方法仍为单步隐式格式，可仅通过极少代码修改便轻松集成至现有的二阶广义-$\alpha$求解器中。