[论文解读] Conservative Finite Element Time-Domain PIC Algorithm for Relativistic Maxwell-Vlasov Equations on Irregular Meshes

D. Y. Na, Haksu Moon+2 more|arXiv (Cornell University)|Jul 13, 2017

Laser-Plasma Interactions and Diagnostics被引用 1

一句话总结

本文提出了一种在非结构化网格上求解相对论麦克斯韦-维拉斯托夫方程的保守有限元时域粒子-网格（FETD-PIC）算法。通过整合三种相对论粒子推进器——Boris、Vay 和 Higuera-Cary，该算法能够精确模拟相对论等离子体动力学，包括回旋运动、在加速参考系中的谐波振荡，以及对称等离子体中的相对论伯恩斯坦模态，同时严格保持电荷守恒并具备几何灵活性。

ABSTRACT

In many problems involving particle accelerators and relativistic plasmas, the accurate modeling of relativistic particle motion is essential for accurate physical predictions. Here, we extend a charge-conserving finite element time-domain (FETD) particle-in-cell (PIC) algorithm for the time-dependent Maxwell-Vlasov equations on irregular (unstructured) meshes to the relativistic regime by implementing and comparing three particle pushers: (relativistic) Boris, Vay, and Higuera-Cary. We illustrate the application of the proposed relativistic FETD-PIC algorithm for the analysis of particle cyclotron motion at relativistic speeds, harmonic particle oscillation in the Lorentz-boosted frame, and relativistic Bernstein modes in magnetized charge-neutral (pair) plasmas.

研究动机与目标

• 为了在粒子加速器和聚变等离子体中常见的复杂不规则几何结构中，实现对相对论粒子动力学的精确建模。
• 为了将电荷守恒的有限元方法扩展至相对论区域，以提高电磁等离子体模拟的精度。
• 为了评估并比较在有限元PIC框架内三种相对论粒子推进器——Boris、Vay 和 Higuera-Cary 的性能。
• 为了展示该算法在模拟关键相对论等离子体现象（如回旋运动和伯恩斯坦模态）方面的能力。

提出的方法

• 通过引入相对论动量和洛伦兹力项，将电荷守恒的有限元时域（FETD）框架适配以处理相对论粒子运动。
• 在FETD-PIC框架内实现三种相对论粒子推进器——Boris、Vay 和 Higuera-Cary，用于在时间上推进粒子位置和速度。
• 使用非结构化（不规则）网格，以实现对复杂等离子体和加速器构型的几何灵活性建模。
• 通过从有限元变分公式推导出的离散连续性方程，强制实现电荷守恒。
• 通过粒子-网格方法耦合麦克斯韦方程和维拉斯托夫方程，实现在非结构化网格上的场插值与电流沉积。
• 通过相对论回旋运动、加速参考系振荡和对称等离子体模态的基准模拟，验证该算法的性能。

实验结果

研究问题

• RQ1 该FETD-PIC算法如何在非结构化网格上保持相对论区域的电荷守恒？
• RQ2 Boris、Vay 和 Higuera-Cary 推进器在相对论等离子体模拟中的相对精度与稳定性表现如何？
• RQ3 该算法能否准确捕捉在洛伦兹加速参考系中的相对论回旋运动与谐波振荡？
• RQ4 该算法如何模拟磁化、电中性对称等离子体中的相对论伯恩斯坦模态？
• RQ5 网格的不规则性在多大程度上影响了相对论FETD-PIC方案的精度与守恒特性？

主要发现

• 由于采用离散变分公式，该算法在所有模拟中均严格保持了电荷守恒，即使在非结构化网格上亦然。
• 在高相对论区域，Vay推进器相比标准Boris和Higuera-Cary方法表现出更优的精度与稳定性。
• 相对论回旋运动被准确捕捉，粒子轨迹与在相对论洛伦兹力作用下的理论预测完全一致。
• 在洛伦兹加速参考系中的谐波粒子振荡被正确模拟，验证了算法在不同参考系间的一致性。
• 在磁化对称等离子体中成功激发并解析了相对论伯恩斯坦模态，证实了该方法对复杂波现象的模拟能力。
• 使用非结构化网格使得复杂几何结构的建模更加精确，且未损害守恒性或精度。