[论文解读] Continuous-stage Runge-Kutta methods based on weighted orthogonal polynomials

Wensheng Tang|arXiv (Cornell University)|May 25, 2018

Numerical methods for differential equations参考文献 43被引用 9

一句话总结

本文提出了一种利用加权正交多项式（如雅可比多项式）构建连续阶 Runge-Kutta（csRK）方法的新框架，以避免求解复杂的非线性阶条件。通过利用成熟的加权插值型求积理论，该方法实现了几何积分器（如辛和对称方法）的系统构造，同时将 csRK 方法的应用范围从勒让德多项式扩展至更广泛的正交多项式。

ABSTRACT

We develop continuous-stage Runge-Kutta methods based on weighted orthogonal polynomials in this paper. There are two main highlighted merits for developing such methods: Firstly, we do not need to study the tedious solution of multi-variable nonlinear algebraic equations associated with order conditions; Secondly, the well-known weighted interpolatory quadrature theory appeared in every numerical analysis textbook can be directly and conveniently used. By introducing weight function, various orthogonal polynomials can be used in the construction of Runge-Kutta-type methods. It turns out that new families of Runge-Kutta-type methods with special properties (e.g., symplectic, symmetric etc.) can be constructed in batches, and hopefully it may produce new applications in numerical ordinary differential equations.

研究动机与目标

• 克服以往 csRK 方法仅依赖勒让德多项式的局限性。
• 提出一种使用任意加权正交多项式构造 csRK 方法的一般性框架。
• 通过避免求解多变量非线性代数系统，简化阶条件的推导过程。
• 通过基于求积的设计方法，系统化构造保持结构的积分器（如辛、对称方法）。
• 将正交多项式技术的应用范围扩展至数值 ODE 中的连续阶方法。

提出的方法

• 引入一种使用权函数定义正交多项式基以构造 Butcher 系数的 csRK 形式。
• 采用加权插值型求积规则（数值分析中的标准方法）直接构造系数 $A_{\tau,\sigma}$、$B_\tau$ 和 $C_\tau$。
• 将 Butcher 核表示为正交多项式展开，例如 $A_{\tau,\sigma} = \frac{1}{2} + \sum_{i+j \text{ odd}} \alpha_{(i,j)} L_i(\tau)L_j(\sigma)$，其中 $\alpha_{(i,j)}$ 为实系数。
• 利用加权多项式（如雅可比、勒让德）的正交性和求积性质，确保方法的稳定性和精度。
• 应用高斯-克里斯特夫求积规则理论，推导出具有最优收敛特性的高阶方法。
• 证明该方法不依赖于特定正交多项式族的选择，从而实现广泛推广。

实验结果

研究问题

• RQ1 能否使用加权正交多项式构造 csRK 方法，而无需求解复杂的非线性阶条件？
• RQ2 如何利用现有的求积理论设计具有所需几何性质的 csRK 方法？
• RQ3 能否在此框架下系统地生成超越勒让德多项式的辛和对称 csRK 方法？
• RQ4 权函数在推广 csRK 方法类别中起什么作用？
• RQ5 是否可能使用这种基于多项式的方法构造高阶几何积分器？

主要发现

• 通过使用正交多项式展开，该方法消除了与阶条件相关的多变量非线性代数系统求解需求。
• 加权插值型求积理论使得 csRK 系数能够被直接且系统地构造，显著简化了设计过程。
• 可使用雅可比多项式及其他加权正交多项式构造出新的辛和对称 csRK 方法族。
• 该框架具有通用性，不限于勒让德多项式，可广泛应用于不同多项式族。
• 该方法为未来开发其他几何积分器（如能量保持和共轭辛方法）提供了支持。
• 该方法揭示了谱方法与 csRK 方法之间的深层联系，提示其具有高阶精度和超收敛的潜力。