[论文解读] Explicit non-canonical symplectic algorithms for charged particle dynamics

He Yang, Yajuan Sun+3 more|arXiv (Cornell University)|Sep 25, 2015

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

一句话总结

本论文通过将哈密顿量分解为四个可解子系统，提出了一类显式、K-辛的数值方法，用于带电粒子动力学。每个子系统均保持非典型辛结构。所得到的方法具有二阶精度，长期稳定且能量守恒，相较于标准龙格-库塔方法，在回旋运动和导心漂移的长时间模拟中表现更优。

ABSTRACT

We study the non-canonical symplectic structure, or K-symplectic structure inherited by the charged particle dynamics. Based on the splitting technique, we construct non-canonical symplectic methods which is explicit and stable for the long-term simulation. The key point of splitting is to decompose the Hamiltonian as four parts, so that the resulting four subsystems have the same structure and can be solved exactly. This guarantees the K-symplectic preservation of the numerical methods constructed by composing the exact solutions of the subsystems. The error convergency of numerical solutions is analyzed by means of the Darboux transformation. The numerical experiment display the long-term stability and efficiency for these methods.

研究动机与目标

• 开发用于电磁场中带电粒子动力学的长期稳定数值方法，特别是针对托卡马克等多时间尺度问题。
• 解决传统典型辛方法的局限性，即对非典型系统而言隐式且计算成本高昂。
• 通过利用洛伦兹力系统内在的非典型辛（K-辛）结构，构建显式、结构保持的积分器。
• 通过哈密顿量分裂与精确解组合，保持K-辛结构，确保长期精度与稳定性。
• 通过回旋运动与导心漂移的长时间模拟，验证该方法相较于标准龙格-库塔格式的优越性。

提出的方法

• 利用反对称矩阵 $ K({f z}) $，将洛伦兹力系统重写为非典型哈密顿形式，揭示了在 $ ({f x}, {f v}) $ 坐标下的K-辛结构。
• 将哈密顿量划分为四部分：三个与速度相关的项 $ H_{v_i} $ 和一个势能项 $ H_{ ho} $，每一项对应一个具有相同K-辛结构的子系统。
• 通过利用磁场与电场结构导出的解析解，精确求解每个子系统，实现显式时间推进。
• 通过组合四个子系统的精确流构建完整数值方法，利用组合方法确保离散K-辛结构的保持。
• 通过达布变换进行误差分析，该变换将K-辛结构映射为典型结构，从而在二阶方法中以 $ h^2 $ 形式建立误差界。
• 在二维对称电磁场中开展数值实验，对比长期行为、能量误差与解误差，与四阶龙格-库塔方法进行比较。

实验结果

研究问题

• RQ1 能否构造出保持带电粒子动力学非典型辛（K-辛）结构的显式数值方法？
• RQ2 如何将哈密顿量分裂技术应用于洛伦兹力系统，将其分解为具有相同K-辛结构的可解子系统？
• RQ3 所得到的K-辛方法在长期稳定性与精度方面相较于标准龙格-库塔格式表现如何？
• RQ4 K-辛结构的保持在多大程度上限制了长时间积分下的能量与解误差？
• RQ5 在固定误差容限下，新方法的计算成本与经典方法相比如何？

主要发现

• 所构建的二阶K-辛方法在 $ 2 \times 10^4 $ 个时间步内保持能量误差在 $ O(h^2) $ 范围内，而四阶龙格-库塔方法则表现出显著的能量阻尼。
• 在190,000步后，K-辛方法仍能保持正确的螺旋轨道并实现稳定的导心漂移，而龙格-库塔方法因数值耗散导致回旋运动失真而失效。
• K-辛方法的解误差随时间缓慢增长，即使在长时间积分后，其误差仍低于龙格-库塔方法，尤其在较大时间步长下表现更优。
• 当解误差低于1%时，二阶K-辛方法的计算时间少于一阶变体，显示出更高的效率。
• 达布变换使理论误差分析成为可能，通过将K-辛方法与典型辛方法关联，解释了 $ O(h^2) $ 收敛率的合理性。
• 数值结果证实，由于其结构保持特性，K-辛方法在长时间模拟中表现更优，尤其在具有多时间尺度的系统中。