[论文解读] On algebraic structures of numerical integration on vector spaces and manifolds

Alexander Selvikvåg Lundervold, Hans Munthe–Kaas|arXiv (Cornell University)|Jan 1, 2012

Numerical methods for differential equations被引用 12

一句话总结

本文提出了基于预李代数、后李代数、D-代数以及非交换霍普夫代数的代数与组合框架，用于分析向量空间与流形上的几何数值积分。它通过卷积代数建立了李-巴特勒级数、逆误差分析与代换律之间的联系，表明代数结构如何统一数值方法的阶理论与结构保持性质。

ABSTRACT

Abstract. Numerical analysis of time-integration algorithms has been applying ad-vanced algebraic techniques for more than fourty years. An explicit description of the group of characters in the Butcher–Connes–Kreimer Hopf algebra first appeared in Butcher’s work on composition of integration methods in 1972. In more recent years, the analysis of structure preserving algorithms, geometric integration techniques and in-tegration algorithms on manifolds have motivated the incorporation of other algebraic structures in numerical analysis. In this paper we will survey algebraic structures that have found applications within these areas. This includes pre-Lie structures for the ge-ometry of flat and torsion free connections appearing in the analysis of numerical flows on vector spaces. The much more recent post-Lie and D-algebras appear in the analysis of flows on manifolds with flat connections with constant torsion. Dynkin and Eule-rian idempotents appear in the analysis of non-autonomous flows and in backward error analysis. Non-commutative Bell polynomials and a non-commutative Faa ̀ di Bruno Hopf algebra are other examples of structures appearing naturally in the numerical analysis of integration on manifolds.

研究动机与目标

• 使用先进的代数与组合工具统一分析几何数值积分方法。
• 通过李-巴特勒级数与D-代数，将经典巴特勒级数理论从向量空间推广至流形。
• 阐明预李代数与后李代数等代数结构在建模具有挠率与曲率的流中的作用。
• 利用非交换霍普夫代数形式化数值积分器的逆误差分析与代换律。
• 通过代数恒等式与特征标，提供系统化框架以理解数值方法中的阶条件与结构保持性。

提出的方法

• 利用巴特勒级数与李-巴特勒级数（LB级数）作为以有根树为指标的正式展开，表示数值流。
• 应用预李代数来建模向量空间流中平坦无挠率的联络，捕捉复合规则的代数结构。
• 引入后李代数与D-代数以处理具有常数挠率的流形上的流，推广预李代数设定。
• 使用非交换法里-布龙尼霍普夫代数与非交换贝尔多项式来描述向量场的复合与代换。
• 利用霍普夫代数中的卷积积与特征标（例如巴特勒-康奈斯-克雷默）形式化逆误差分析与代换律。
• 通过修剪与分拆余乘法推导代换特征标的递归公式，实现修正向量场的计算。

实验结果

研究问题

• RQ1 预李代数与后李代数如何用于刻画具有或不具有挠率的流形上数值流的代数结构？
• RQ2 非交换霍普夫代数在李-巴特勒级数方法的逆误差分析中起什么作用？
• RQ3 狄尼金与欧拉型幂等元如何在非自治系统与结构保持积分器的分析中出现？
• RQ4 李-巴特勒级数的代换律背后的代数机制是什么？它与霍普夫代数中的卷积有何关联？
• RQ5 数值积分器的结构能否通过LB级数中换位子在给定阶数内消失来完全刻画？

主要发现

• 微分方程的精确解可表示为无树换位子的李-巴特勒级数，意味着结构保持方法必须在其阶数内换位子为零。
• 指数欧拉法对应平凡的LB级数 γEuler = ，其中代换特征标 αt∗( ) = 1/2，αt∗( ) = 1/4，等等。
• 李-隐式中点法导出系数的递归公式：α( ) = 1/2，α(τ) = 1/(2^j j!) α(τ1)⋯α(τj)，其中 τ = B+(τ1⋯τj)。
• LB级数的代换律由 BBf(α)(β) = Bf(α∗β) 给出，其中 α∗ 是通过对偶配对与修剪操作定义的D-代数同态。
• 代换特征标 αt∗ 满足涉及分拆余乘法与修剪的递归公式，从而可显式计算修正向量场。
• 该框架通过将阶条件与结构保持性嵌入霍普夫代数与操作代数的语言中，实现了统一处理。