[论文解读] New applications for the Boris Spectral Deferred Correction algorithm for plasma simulations

Kris Smedt, Daniel Ruprecht+3 more|arXiv (Cornell University)|Oct 15, 2021

Electromagnetic Simulation and Numerical Methods参考文献 46被引用 3

一句话总结

本文提出在粒子-网格（PIC）模拟中使用Boris谱延迟校正（Boris-SDC）算法作为高阶粒子推进器，相较于标准Boris方法，显著提升了精度并减少了数值漂移。该研究进一步提出了Boris-SDC的相对论扩展版本，在速度超过0.99c时仍保持高阶收敛性，并在无外力条件下表现出无漂移行为，尽管当前存在计算成本的权衡，但该方法在高精度等离子体模拟中展现出显著优势。

ABSTRACT

The paper investigates two new use cases for the Boris Spectral Deferred Corrections (Boris-SDC) time integrator for plasma simulations. First, we show that using Boris-SDC as a particle pusher in an electrostatic particle-in-cell (PIC) code can, at least in the linear regime, improve simulation accuracy compared with the standard second order Boris method. In some instances, the higher order of Boris-SDC even allows a much larger time step, leading to modest computational gains. Second, we propose a modification of Boris-SDC for the relativistic regime. Based on an implementation of Boris-SDC in the extsc{runko} PIC code, we demonstrate for a relativistic Penning trap that Boris-SDC retains its high order of convergence for velocities ranging from $0.5c$ to $>0.99c$. We also show that for the force-free case where acceleration from electric and magnetic field cancel, Boris-SDC produces less numerical drift than Boris.

研究动机与目标

• 评估Boris-SDC作为静电PIC模拟中高阶粒子推进器的性能，其中标准Boris方法为二阶精度。
• 解决PIC代码中二阶时间积分的局限性，该局限性尽管采用高阶空间方法，仍制约整体精度。
• 将Boris-SDC扩展至相对论区域，实现对相对论洛伦兹方程的高阶积分。
• 评估在无外力场景下的数值漂移，即电场与磁场力相互抵消的情况，这是长期模拟保真度的关键测试。
• 确定Boris-SDC的计算成本是否能通过其精度提升在真实等离子体构型中得到合理化。

提出的方法

• 将原本专用于单粒子动力学的Boris-SDC算法，适配为RUNKO PIC代码框架中的粒子推进器。
• 利用谱延迟校正（SDC）通过在配置点上进行多轮迭代，对低阶时间积分进行逐步校正，从而实现高阶精度。
• 通过重新表述速度更新步骤，实现Boris-SDC的相对论化改进，以在电磁场存在下正确处理相对论动量与洛伦兹力。
• 将该方法应用于两个测试案例：静电PIC模拟中的两束流不稳定性与兰道阻尼，以评估其精度与工作-精度效率。
• 在包含速度从0.5c到>0.99c的彭宁阱设置中测试相对论Boris-SDC，以验证其收敛阶数与稳定性。
• 在无外力场景下评估数值漂移，此时净加速度为零，比较Boris、Vay与Boris-SDC在不同节点数与迭代次数下的表现。

实验结果

研究问题

• RQ1 Boris-SDC能否在完整静电PIC代码中作为有效的高阶粒子推进器，相较于标准Boris方法显著提升精度？
• RQ2 在PIC模拟中使用Boris-SDC是否能带来可测量的计算效率增益，尤其是在结合高阶空间离散化时？
• RQ3 Boris-SDC能否成功推广至相对论区域，同时保持其高阶收敛性与稳定性？
• RQ4 在无外力条件下，Boris-SDC在数值漂移方面表现如何，此时长期能量与动量守恒至关重要？
• RQ5 为何Boris-SDC的某些配置表现出零或接近零的数值漂移？该现象是否具有理论可解释性？

主要发现

• Boris-SDC在非相对论与相对论区域均实现了高阶收敛，观测到的收敛阶数与理论预期一致，最高可达积分求积阶数。
• 在非相对论情况下，Boris-SDC在相同时间步长下相比标准Boris方法显著降低了位置与速度误差，但因主导的空间离散化误差，计算成本抵消了效率增益。
• 在相对论彭宁阱模拟中，Boris-SDC在0.5c至>0.99c的速度范围内均保持高阶收敛，展现出在极端条件下的鲁棒性。
• 在无外力情况下，使用M=3个节点与K=1次迭代的Boris-SDC未表现出任何数值漂移，而M=5个节点与K=1次迭代的配置仅表现出极小漂移——低于Boris与Vay积分器。
• 部分Boris-SDC配置完全无漂移，但作者目前尚无理论解释该现象。
• 当误差低于10−3时，相对论模拟中观察到计算增益，表明在高精度区域可能存在效率优势。