[论文解读] Overcoming timestep limitations in boosted-frame Particle-In-Cell simulations of plasma-based acceleration

Olga V. Shapoval, Remi Lehé+4 more|arXiv (Cornell University)|Apr 28, 2021

Laser-Plasma Interactions and Diagnostics参考文献 29被引用 7

一句话总结

本文提出平均伽利略 PSATD 算法，以克服等离子体加速器的增强帧粒子-网格（PIC）模拟中的时间步长限制。通过在傅里叶空间中对电磁场进行时间上的解析平均，该方法可在不引发数值不稳定性的情况下实现大时间步长（c∆t ≫ ∆x, ∆y），从而充分发挥增强帧技术在高分辨率横向模拟中的计算加速潜力。

ABSTRACT

Explicit electromagnetic Particle-In-Cell (PIC) codes are typically limited by the Courant- Friedrichs-Lewy (CFL) condition, which implies that the timestep multiplied by the speed of light must be smaller than the smallest cell size. In the case of boosted-frame PIC simulations of plasma-based acceleration, this limitation can be a major hinderance as the cells are often very elongated along the longitudinal direction and the timestep is thus limited by the small, transverse cell size. This entails many small-timestep PIC iterations, and can limit the potential speed-up of the boosted-frame technique. Here, by using a CFL-free analytical spectral solver, and by mitigating additional numerical instabilities that arise at large timestep, we show that it is possible to overcome traditional limitations on the timestep and thereby realize the full potential of the boosted-frame technique over a much wider range of parameters.

研究动机与目标

• 解决由于横向网格尺寸限制导致的增强帧 PIC 模拟中时间步长受限的问题，尽管纵向分辨率已放宽。
• 解决在基于 PSATD 的增强帧模拟中，大时间步长下数值切伦科夫不稳定性（NCI）重现的问题。
• 在 PSATD 模拟中实现大时间步长（c∆t ≫ ∆x, ∆y）的同时，保持模拟的精度与稳定性。
• 将增强帧技术的计算加速潜力扩展到更多需要高横向分辨率的等离子体加速器构型中。

提出的方法

• 提出平均伽利略 PSATD 算法，该算法在傅里叶空间中对时间区间 [(n−1/2)∆t, (n+1/2)∆t] 内的电磁场进行解析平均。
• 将 PSATD 框架适配为利用恒定电流下麦克斯韦方程的解析解，计算时间平均场 ⟨ˆE⟩n 和 ⟨ˆB⟩n。
• 推导出包含复系数（Ψ1, Ψ2, A1, A2, Cρ）的新场更新方程，以同时考虑伽利略漂移与谱精度。
• 修改标准 PIC 循环，使用时间平均场进行粒子推进，确保场收集与动量更新的一致性。
• 在傅里叶空间中采用 p 阶离散化，以在不同网格分辨率下保持高精度与稳定性。
• 在开源 WarpX 代码中实现该算法，支持 CPU 与 GPU 平台。

实验结果

研究问题

• RQ1 是否可以在增强帧 PSATD 模拟中安全地使用大时间步长（c∆t ≫ ∆x, ∆y）而不引发数值不稳定性？
• RQ2 为何标准伽利略 PSATD 算法在 c∆t > ∆x, ∆y 时会变得不稳定，尽管其为 CFL 自由？
• RQ3 是否可通过在傅里叶空间中进行时间平均来修改场收集过程，从而缓解该不稳定性？
• RQ4 新算法在真实等离子体加速器场景中在多大程度上保持物理完整性与计算加速潜力？
• RQ5 该方法是否可推广至其他需要高横向分辨率但无需高时间分辨率的 PIC 模拟？

主要发现

• 平均伽利略 PSATD 算法成功抑制了大时间步长（c∆t ≫ ∆x, ∆y）下数值切伦科夫不稳定性（NCI）的重现，实现了稳定模拟。
• 在 c∆t = ∆z（远大于横向网格尺寸）的模拟中，NCI 模式无增长现象，而标准伽利略 PSATD 算法则出现不稳定性。
• 在二维均匀电浆模拟中，新算法在 c∆t 达到 ∆z 时仍保持稳定，而标准方法在 c∆t > ∆x 时即变得不稳定。
• 在二维激光尾场加速（LWFA）模拟中，平均 PSATD 方法在大时间步长下仍能正确再现尾场结构与电子束动力学，与参考模拟结果一致。
• 在三维电浆尾场加速（PWFA）模拟中，该方法实现了 c∆t ≫ ∆x, ∆y 的稳定高速模拟，显著提升了计算效率。
• 该算法已在 WarpX 代码中实现并验证，证明其与自适应网格加密及多平台（CPU/GPU）高性能计算环境兼容。