2N存储Runge-Kutta方法:Butcher表的c反射对称性与分解
数值分析
2025-06-11 v1 数值分析
高能物理 - 格点
计算物理
摘要
Williamson型的低存储Runge-Kutta格式,即所谓的2N存储格式,作为近期工作的后续被进一步研究。研究发现,扩充的Butcher表分解为具有特殊性质的矩阵乘积。这些性质揭示出,全局精度阶数小于五的2N存储方法具有一种对称性,称为c反射对称性,即大多数方法成对存在。找到了关联成对Butcher表的一种变换,并证明了c反射方法满足与原方法相同的阶条件。呈现了验证解析结果的数值证据。由Carpenter和Kennedy首次探索的(5,4)方法的分支被数值构造。考察了四个新的(5,4)格式(其系数可用根式表示)和一个具有有理系数的格式以作说明。研究了八个新的(6,4)格式(部分可用有理数或根式表示)和一个(8,4)格式,以理解c反射对称性对阶段数更多的方法的实际意义。在四阶及以上的2N存储Runge-Kutta方法缺乏闭式解析解的情况下,一般对称性质以及此处给出的一些特定解析解,可能有助于2N存储格式的发展与优化。
引用
@article{arxiv.2506.08444,
title = {2N-storage Runge-Kutta methods: c-reflection symmetry and factorization of the Butcher tableau},
author = {Alexei Bazavov},
journal= {arXiv preprint arXiv:2506.08444},
year = {2025}
}
备注
41 pages, 7 figures