[论文解读] Asymptotically stable particle-in-cell methods for the Vlasov-Poisson system with a strong external magnetic field

Francis Filbet, L. Miguel Rodrigues|arXiv (Cornell University)|Nov 23, 2015

Gas Dynamics and Kinetic Theory参考文献 43被引用 63

一句话总结

本文提出了一种用于强外部磁场下Vlasov-Poisson系统的渐近保持半隐式粒子-网格（PIC）格式，可在无严格时间步长限制的情况下实现稳定且精确的模拟。该方法在Larmor半径较小时渐近捕捉导向中心动力学，一致地离散了涡度形式的不可压缩Euler方程。

ABSTRACT

This paper deals with the numerical resolution of the Vlasov-Poissonsystem with a strong external magnetic field by Particle-In-Cell(PIC) methods. In this regime, classical PIC methods are subject tostability constraints on the time and space steps related to the smallLarmor radius and plasma frequency. Here, we propose anasymptotic-preserving PIC scheme which is not subjected to theselimitations. Our approach is based on first and higher order semi-implicit numericalschemes already validated on dissipative systems. Additionally, when the magnitude of the external magneticfield becomes large, this method provides a consistent PICdiscretization of the guiding-center equation, that is, incompressibleEuler equation in vorticity form. We propose several numerical experiments which provide a solid validation of the method and its underlying concepts.

研究动机与目标

• 解决经典PIC方法在模拟强外部磁场下电浆时出现的不稳定性和严格的时间步长限制问题。
• 开发一种数值格式，即使在磁场强度增加（对应Larmor半径减小）时仍保持稳定与精确。
• 确保该方法在Larmor半径趋于零的极限下渐近恢复导向中心近似。
• 在强磁场区域提供不可压缩Euler方程在涡度形式下的粒子基一致离散化。
• 通过单粒子运动和双极子不稳定性数值模拟验证该格式，展示其长期稳定性和精确性。

提出的方法

• 为强磁场下的Vlasov-Poisson系统提出一阶、二阶和三阶半隐式Runge-Kutta时间格式。
• 引入一种时间分裂方法，将快速回旋运动与慢速导向中心动力学解耦，从而允许更大的时间步长。
• 采用粒子基方法离散分布函数，粒子轨迹通过Vlasov方程的特征曲线推进。
• 采用有限差分法求解导向中心模型中的电势Poisson方程，以供对比。
• 确保该格式在ε → 0极限下保持初始的精度阶数，其中ε与Larmor频率的倒数成正比。
• 形式化证明该半隐式格式在ε → 0极限下与导向中心模型保持稳定性和一致性。

实验结果

研究问题

• RQ1 半隐式时间格式能否稳定模拟强外部磁场下的Vlasov-Poisson系统？
• RQ2 随着磁场强度增加，所提格式是否渐近收敛至导向中心模型？
• RQ3 该方法能否在无时间步长限制的情况下准确捕捉长时动态，如双极子不稳定性？
• RQ4 在ε较小时，能量守恒与电场范数如何演化？该格式是否保持这些性质？
• RQ5 在模拟不可压缩Euler方程的涡度形式时，粒子-网格方法在多大程度上保持精度与稳定性？

主要发现

• 一阶半隐式格式在光滑电磁场下无条件稳定，并与导向中心模型渐近一致。
• 当ε = 10⁻²时，长时间模拟下相对能量变化保持在10⁻³以下，表明能量守恒性能极佳。
• Vlasov-Poisson系统与导向中心模型中电场范数‖E(t)‖∞的时间演化高度一致，证实了格式的一致性。
• 在t = 120时，Vlasov-Poisson系统与导向中心模型的密度分布高度吻合，即使采用大时间步长和高粒子数亦然。
• 粒子-网格方法成功捕捉了双极子不稳定性的发展过程，即使在ε → 0区域也未出现过度数值振荡。
• 三阶半隐式格式保持了高精度与稳定性，在长时间模拟中未观察到振荡或发散。