[论文解读] Stochastic Variational Partitioned Runge-Kutta Integrators for Constrained Systems

Nawaf Bou‐Rabee, Houman Owhadi|ArXiv.org|Sep 14, 2007

Numerical methods for differential equations参考文献 19被引用 12

一句话总结

本文提出用于受白噪声影响的约束力学系统的随机变分分区龙格-库塔（VPRK）积分器，将变分积分器扩展至随机设置下的完整约束系统。主要贡献在于开发了一类一阶与二阶强收敛的RATTLE型积分器，能够保持辛结构与动量守恒，并通过后向误差分析及Vanden-Eijnden与Cicotti（2006）提出的方法，证明了其强收敛阶数。

ABSTRACT

Stochastic variational integrators for constrained, stochastic mechanical systems are developed in this paper. The main results of the paper are twofold: an equivalence is established between a stochastic Hamilton-Pontryagin (HP) principle in generalized coordinates and constrained coordinates via Lagrange multipliers, and variational partitioned Runge-Kutta (VPRK) integrators are extended to this class of systems. Among these integrators are first and second-order strongly convergent RATTLE-type integrators. We prove order of accuracy of the methods provided. The paper also reviews the deterministic treatment of VPRK integrators from the HP viewpoint.

研究动机与目标

• 为具有完整约束的随机力学系统开发结构保持型积分器，其中标准方法无法正确捕捉随机扰动。
• 通过在广义坐标与约束坐标下形式化随机哈密顿-庞特里亚金原理，将变分积分器理论扩展至随机设置。
• 构建一阶与二阶强收敛的RATTLE型积分器，在随机扰动下保持辛结构与动量守恒。
• 利用Vanden-Eijnden与Cicotti（2006）提出的方法，证明所提积分器的强收敛阶数，确保长时间模拟的可靠性。
• 实现对随机约束系统中能量行为、自相关性与经验分布等统计特性的精确计算。

提出的方法

• 在广义坐标下形式化随机哈密顿-庞特里亚金原理，并通过拉格朗日乘子证明其与约束形式的等价性。
• 通过离散变分原理推导随机变分积分器，确保辛结构保持性与守恒量的离散诺特定理。
• 通过将更新步骤分解为速度与位置步骤，并将结果投影至约束流形，构造RATTLE型积分器。
• 采用二阶精度的分区龙格-库塔格式构建离散变分积分器，对拉格朗日乘子与动量更新进行一致处理。
• 应用后向误差分析及Vanden-Eijnden与Cicotti（2006）提出的方法，证明强收敛阶数，包括至O(h²)的误差展开。
• 通过证明无约束更新步骤与投影步骤在均方意义下均保持辛形式，验证离散流映射的均方辛结构保持性。

实验结果

研究问题

• RQ1 能否将变分积分器扩展至具有完整约束的随机力学系统，同时保持辛结构与动量守恒？
• RQ2 随机VPRK积分器在约束系统中的强收敛阶数是多少？如何严格证明？
• RQ3 所提积分器在捕捉随机能量与动量扰动方面，与标准投影法或高阶方法相比有何表现？
• RQ4 SVI的结构保持特性（如正确能量行为与统计精度）在存在约束与噪声时是否仍能保持？

主要发现

• 成功构造了一类一阶强收敛、辛结构保持的RATTLE型积分器，用于约束随机力学系统，其强收敛阶数通过O(h²)误差展开得到严格证明。
• 该积分器的离散流映射为均方辛结构保持，确保在随机扰动下长期保持辛结构。
• 严格建立了随机哈密顿-庞特里亚金原理在广义坐标与约束坐标下的等价性，实现了变分原理的一致表述。
• 对随机扰动下的球面摆进行数值实验，验证了其强收敛阶数为1.0，弱收敛阶数为1.5，与理论预测一致。
• 积分器能正确捕捉能量与动量中的随机扰动，即使在大时间步长与非小噪声振幅下亦无虚假漂移。
• 该方法保持了随机效应的力学本质——即能量与动量的变化仅由噪声引起，而非数值误差所致。