[论文解读] Splitting methods for Fourier spectral discretizations of the strongly magnetized Vlasov-Poisson and the Vlasov-Maxwell system

J. Ameres|arXiv (Cornell University)|Jul 11, 2019

Gas Dynamics and Kinetic Theory参考文献 40被引用 5

一句话总结

本文提出了一种用于强磁场下Vlasov-Poisson与Vlasov-Maxwell系统傅里叶谱离散化的新型分裂方法，实现了四维与三维高精度模拟。在Kelvin–Helmholtz与Weibel流体不稳定性等挑战性测试案例中，展示了稳定性和几何守恒性，即使在低分辨率下使用Strang分裂与哈密顿积分器，能量误差也低于10⁻¹²。

ABSTRACT

Fourier spectral discretizations belong to the most straightforward methods for solving the unmagnetized Vlasov--Poisson system in low dimensions. In this article, this highly accurate approach is extended two the four-dimensional magnetized Vlasov--Poisson system with new splitting methods suited for strong magnetic fields. Consequently, a comparison to the asymptotic fluid model is provided at the example of a turbulent Kelvin--Helmholtz instability. For the three dimensional electromagnetic Vlasov--Maxwell system different novel charge conserving implementations of a Hamiltonian splitting are discussed and simulation results of the Weibel streaming instability are presented.

研究动机与目标

• 解决强磁场下高维Vlasov系统中傅里叶谱离散化缺乏几何分裂方法的问题。
• 为强磁场下四维磁化Vlasov-Poisson系统开发高效、结构保持的时间积分器。
• 将哈密顿分裂技术扩展至三维Vlasov-Maxwell系统，实现电荷与能量守恒。
• 在基准等离子体不稳定性（包括Kelvin–Helmholtz与Weibel流体不稳定性）上验证方法，涵盖低分辨率与高分辨率条件。
• 证明所提格式在粗网格与强磁场条件下仍能保持长期稳定性和几何精度。

提出的方法

• 应用哈密顿分裂将位置空间与速度空间中的对流解耦，通过傅里叶变换实现精确求解。
• 采用傅里叶谱方法对空间与速度空间进行离散化，利用FFTW实现高效计算。
• 实现指数分裂格式（如Scovel方法与指数Boris方法），以在傅里叶空间中精确积分常系数对流项。
• 通过保持底层哈密顿结构的几何积分器，确保Vlasov-Maxwell系统中的电荷守恒。
• 采用Strang分裂进行时间积分，实现二阶精度，同时保持能量与动量守恒特性。
• 利用贝塞尔函数与高斯求积法处理Ampère方程中的磁旋平均，实现Scovel方法中精确的时间积分。

实验结果

研究问题

• RQ1 能否通过新型分裂格式有效将傅里叶谱方法扩展至四维强磁场磁化Vlasov-Poisson系统？
• RQ2 所提出的分裂方法在强磁场下模拟湍流不稳定性（如Kelvin–Helmholtz）时表现如何？
• RQ3 哈密顿分裂能否适配至三维Vlasov-Maxwell系统，以确保电荷与能量守恒？
• RQ4 分辨率与时间步长对Weibel型不稳定性长期模拟中能量与动量误差的影响如何？
• RQ5 相较于标准伪谱方法，新格式在多大程度上保持了几何结构与稳定性？

主要发现

• 所提出的Scovel分裂方法在广泛磁场强度范围内保持了高精度与稳定性，即使在无渐近极限条件下亦然。
• 在Kelvin–Helmholtz不稳定性中，强磁场下的动力学模型与流体模型高度一致，验证了渐近区的合理性。
• 采用哈密顿分裂的Vlasov-Maxwell求解器在高分辨率（128×128×128网格点）下能量误差低于1.71×10⁻¹²，与高阶GEMPIC方法相当。
• 即使在低分辨率（32×32×32）下，求解器在长时间模拟中仍保持稳定，能量误差低于10⁻¹²，展现出强鲁棒性。
• Weibel流体不稳定性模拟显示了动能与磁场的稳定增长，对称与非对称情形分别呈现出不同的纤维化结构。
• 贝塞尔函数与高斯求积法的使用使Scovel方法中Ampère方程的精确时间积分成为可能，是实现长期守恒的关键技术突破。