[论文解读] Chebyshev symplectic methods based on continuous-stage Runge-Kutta methods

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

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

一句话总结

本文提出了一类基于第一类和第二类切比雪夫正交多项式的新型辛龙格-库塔方法，通过一种新的加权连续阶段龙格-库塔方法框架构建。所提方法实现了高阶精度并精确保持辛结构，在哈密顿系统数值实验中表现出优越的长期稳定性和能量守恒性，尽管初始误差增长略高，但在定性行为上优于经典的高斯-勒让德方法。

ABSTRACT

We develop Chebyshev symplectic methods based on Chebyshev orthogonal polynomials of the first and second kind separately in this paper. Such type of symplectic methods can be conveniently constructed with the newly-built theory of weighted continuous-stage Runge-Kutta methods. A few numerical experiments are well performed to verify the efficiency of our new methods.

研究动机与目标

• 开发适用于哈密顿系统的新辛积分器，以在长时间区间内保持几何结构。
• 将连续阶段龙格-库塔方法理论扩展至使用第一类和第二类切比雪夫多项式构造辛方法。
• 提供一种替代经典W-变换技术的方法，该技术仅限于勒让德多项式，而本方法可利用任意加权正交多项式进行构造。
• 通过扰动开普勒问题的数值实验，验证所提方法的效率和几何精度。

提出的方法

• 利用新发展的加权连续阶段龙格-库塔（csRK）方法理论，构建辛积分器。
• 使用第一类和第二类切比雪夫多项式作为基函数，定义csRK方法的系数。
• 应用简化条件（B̌(ξ), Č(η), Ď(ζ)）以确保高阶精度和辛性。
• 采用基于切比雪夫节点的插值型求积规则，确保离散流中精确辛性。
• 基于正交多项式根，推导出三阶段四阶和五阶段六阶的Butcher型表。
• 通过满足阶次和辛性积分条件的函数Aτ,σ、Bτ和Cτ，实现连续阶段形式的算法实现。

实验结果

研究问题

• RQ1 能否在连续阶段龙格-库塔框架内，系统性地利用切比雪夫多项式构造辛积分器？
• RQ2 基于切比雪夫的辛方法在能量守恒和误差增长方面与经典高斯-勒让德方法相比表现如何？
• RQ3 该框架能否超越勒让德多项式，例如推广至切比雪夫多项式，以生成新的辛方法？
• RQ4 当应用于刚性或扰动哈密顿系统时，所得辛方法的阶次和稳定性行为如何？

主要发现

• 所提出的四阶和六阶切比雪夫辛方法在长时间积分区间内表现出线性误差增长和有界能量误差，证实了其几何精度。
• 数值结果表明，切比雪夫I型和II型四阶方法与四阶高斯-勒让德方法性能相当，初始解误差略高，但能量守恒性相似。
• 六阶切比雪夫辛方法（I型和II型）在解误差和能量守恒性方面与六阶高斯-勒让德方法几乎完全一致。
• 这些方法为精确辛性，经验证在离散流中保持了辛结构，而谱配点方法仅能近似实现。
• 该构造框架具有通用性，可扩展至切比雪夫以外的其他加权正交多项式，如雅可比或勒让德多项式。