[论文解读] A Hybrid Nodal-Staggered Pseudo-Spectral Electromagnetic Particle-In-Cell Method with Finite-Order Centering

Edoardo Zoni, Remi Lehé+6 more|arXiv (Cornell University)|Jun 23, 2021

Laser-Plasma Interactions and Diagnostics参考文献 61被引用 4

一句话总结

本论文提出了一种新颖的混合伪谱电磁PIC方法，通过有限阶中心化技术将节点网格与交错网格相结合，以提升数值稳定性和精度。该方法在交错Yee网格上求解麦克斯韦方程组，同时通过高阶插值在节点网格上沉积电流并收集力，从而减少数值切伦科夫不稳定性（NCI），并在包括激光等离子体加速和正负电子对产生在内的多种等离子体物理场景中，优于纯节点或纯交错方案。

ABSTRACT

Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon sources driven by ultra-intense lasers. In a PIC code, Maxwell's equations are solved on a grid with a numerical method of choice. This article focuses on pseudo-spectral analytical time-domain (PSATD) algorithms and presents a novel hybrid PSATD PIC scheme that combines the respective advantages of standard nodal and staggered methods. The novelty of the hybrid scheme consists in using finite-order centering of grid quantities between nodal and staggered grids, in order to combine the solution of Maxwell's equations on a staggered grid with the deposition of charges and currents and the gathering of electromagnetic forces on a nodal grid. The correctness and performance of the novel hybrid scheme are assessed by means of numerical tests that employ different classes of PSATD equations in a variety of physical scenarios, ranging from the modeling of electron-positron pair creation in vacuum to the simulation of laser-driven and particle beam-driven plasma wakefield acceleration. It is shown that the novel hybrid scheme offers significant numerical and computational advantages, compared to purely nodal or staggered methods, for all the test cases presented.

研究动机与目标

• 解决电磁PIC模拟中的数值切伦科夫不稳定性（NCI），特别是在相对论性和低密度等离子体中的问题。
• 克服纯节点或纯交错PIC方案的局限性，后者因网格交错导致高色散或不稳定性。
• 开发一种混合方法，利用交错网格在场求解中的稳定性，以及节点网格在粒子推进器中力收集的高精度。
• 引入基于Fornberg系数的有限阶插值，实现节点网格与交错网格之间场和电流数据的中心化。
• 在多种物理条件下展示改进的性能，包括激光驱动和束流驱动的等离子体尾场加速，以及电子-正电子对产生。

提出的方法

• 采用混合PIC循环，其中电荷和电流沉积在节点网格上，而麦克斯韦方程组在交错Yee网格上求解。
• 应用基于Fornberg系数的有限阶插值，以高精度在节点网格与交错网格之间转移场和电流数据。
• 将三类PSATD方案——标准、伽利略变换和平均伽利略变换——适配至交错网格，以提升时域精度。
• 在傅里叶空间中对麦克斯韦方程组进行解析积分，以最小化数值色散并保持稳定性。
• 使用Fornberg插值系数定义一个转移算子，调制网格量的傅里叶变换，有效控制插值系数ζF偏离1的程度。
• 通过更高阶插值最小化0 < ζF < 1的区域，从而确保动量守恒并减少NCI。

实验结果

研究问题

• RQ1 结合节点与交错网格的混合PSATD-PIC方法是否能比纯节点或纯交错方案更有效地减少数值切伦科夫不稳定性？
• RQ2 在节点与交错网格之间进行有限阶中心化，如何影响相对论性等离子体模拟中电磁场解的稳定性和精度？
• RQ3 在不同物理条件下，PSATD变体（标准、伽利略变换、平均伽利略变换）的选择在混合方案中的性能影响有多大？
• RQ4 插值系数ζF在决定NCI抑制中的作用是什么？其偏离1的程度如何依赖于有限阶中心化的阶数？
• RQ5 该混合方案在极端条件下（如电子-正电子对产生和高强度激光-等离子体相互作用）是否仍能保持精度与稳定性？

主要发现

• 与纯节点或纯交错PIC方法相比，该混合方案在低密度和相对论性区域显著减少了数值切伦科夫不稳定性（NCI）。
• 更高阶的有限阶中心化（最高至2m = 8）减少了插值系数ζF偏离1的程度，从而最小化了不稳定性易发区域（0 < ζF < 1）。
• 该方法在多种物理场景中保持了高精度与稳定性，包括激光驱动和粒子束驱动的等离子体尾场加速。
• 数值测试表明，该混合方案在场精度和长期模拟稳定性方面优于标准节点和交错PSATD方法。
• 基于Fornberg的插值方法实现了节点与交错网格之间高阶数据转移，有效平衡了精度与计算效率。
• 色散分析证实，当ζF接近1时，NCI得到抑制，而更高阶中心化在混合配置中更有效地实现了这一点。