[论文解读] On converting any one-step method to a variational integrator of the same order

George W. Patrick, Raymond J. Spiteri+2 more|ArXiv.org|Mar 5, 2009

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

一句话总结

本文提出了一种通用算法，可将任意单步数值积分器转换为约束力学系统的同阶变分积分器。通过利用离散变分原理，并使用标准积分器的输出构建辛层，该方法在保持原方法精度阶数的同时，确保了辛性与动量守恒等结构保持性质。

ABSTRACT

In the formalism of constrained mechanics, such as that which underlies the SHAKE and RATTLE methods of molecular dynamics, we present an algorithm to convert any one-step integration method to a variational integrator of the same order. The one-step method is arbitrary, and the conversion can be automated, resulting in a powerful and flexible approach to the generation of novel variational integrators with arbitrary order.

研究动机与目标

• 开发一种系统化方法，将任意单步数值积分器转换为约束力学系统的同阶变分积分器。
• 在所得积分器中保持辛性与动量守恒等关键几何性质。
• 实现从现有标准积分器自动生成新型高阶变分积分器。
• 确保转换过程保持与原始单步方法相同的精度阶数。
• 提供适用于具有完整约束的系统（包括分子动力学与刚体动力学）的框架。

提出的方法

• 该方法使用标准单步积分器作为“标准层”，生成曲线段并近似作用量值，用于离散拉格朗日量的构建。
• 通过定义偏置参数 α⁻ 和 α⁺ 构建辛层，以控制离散时间演化并确保与变分原理的一致性。
• 利用离散拉格朗日量及其导数的近似值推导离散欧拉-拉格朗日方程，其中关键近似基于标准层的输出。
• 通过矩阵形式如 [I + D²(θᵀg), Dgᵀ; Dg, 0] 解决涉及约束力与乘子的线性方程组，以保证一致性。
• 使用近似式 (4.8) 迭代求解非线性系统，其计算上等价于求解完整线性系统。
• 通过乘子 ν₁, ν₂ 跟踪约束力，其非唯一，但其物理含义（如约束力）在乘积 ν₂ᵀDg(q) 中得以保留。

实验结果

研究问题

• RQ1 能否系统性地将任意单步数值积分器转换为同阶变分积分器？
• RQ2 在标准积分器转换过程中，如何保持辛性与动量守恒？
• RQ3 偏置参数 α⁻ 与 α⁺ 在构建一致离散变分积分器中的作用是什么？
• RQ4 当原始积分器本身不具备辛性时，如何强制实施离散变分原理？
• RQ5 在什么条件下可确保所得积分器保持与原始单步方法相同的精度阶数？

主要发现

• 所提出的算法成功地将任意单步积分器转换为同阶变分积分器，同时保持辛性与动量守恒。
• 该方法完全自动化，适用于任意单步方法，可实现新型高阶变分积分器的生成。
• 离散拉格朗日量基于标准层的输出构建，确保其与经典作用量在所需阶数内一致。
• 离散变分方程的求解依赖于通过 (4.8) 近似求解线性系统，其计算上等价于精确系统。
• 约束力由非唯一的乘子编码，但其物理意义在如 ν₂ᵀDg(q) 的乘积中得以保留，该乘积代表实际的约束力。
• 迭代求解过程收敛至正确的状态变量 q₂ 与 v₂，收敛性通过物理量而非乘子本身进行监控。