[论文解读] On the velocity space discretization for the Vlasov-Poisson system: comparison between Hermite spectral and Particle-in-Cell methods. Part 2: fully-implicit scheme

Enrico Camporeale, Gian Luca Delzanno+2 more|arXiv (Cornell University)|Dec 17, 2013

Magnetic confinement fusion research参考文献 68被引用 9

一句话总结

本文提出了一种用于Vlasov-Poisson系统的全隐式、Hermite谱方法，通过正交Hermite基底对速度空间进行离散化，并通过傅里叶分解对配置空间进行离散化，实现了电荷、动量和能量的精确守恒。该方法在计算成本仅为粒子-网格（PIC）方法一小部分的情况下，实现了比PIC高达10⁷倍的精度，显著优于高分辨率PIC模拟在标准等离子体问题中的表现。

ABSTRACT

We describe a spectral method for the numerical solution of the Vlasov-Poisson system where the velocity space is decomposed by means of an Hermite basis, and the configuration space is discretized via a Fourier decomposition. The novelty of our approach is an implicit time discretization that allows exact conservation of charge, momentum and energy. The computational efficiency and the cost-effectiveness of this method are compared to the fully-implicit PIC method recently introduced by Markidis and Lapenta (2011) and Chen et al. (2011). The following examples are discussed: Langmuir wave, Landau damping, ion-acoustic wave, two-stream instability. The Fourier-Hermite spectral method can achieve solutions that are several orders of magnitude more accurate at a fraction of the cost with respect to PIC. This paper concludes the study presented in Camporeale et al. (2013) where the same method has been described for a semi-implicit time discretization, and was compared against an explicit PIC.

研究动机与目标

• 开发一种用于傅里叶-Hermite谱方法的全隐式时间离散化，确保在Vlasov-Poisson模拟中电荷、动量和能量的精确守恒。
• 将全隐式Hermite谱方法的计算效率和精度与近期开发的全隐式PIC方法进行比较。
• 评估谱速度空间离散化是否能够优于PIC方法，以解析弱非线性动力学现象，如回声和高次谐波生成。
• 在标准等离子体基准测试中，评估Hermite谱方法相对于PIC方法的成本效益，包括朗缪尔波、兰道阻尼、离子声波和两束流不稳定性。

提出的方法

• 通过正交Hermite基底对速度空间进行离散化，使分布函数可表示为Hermite多项式的级数展开。
• 通过傅里叶谱分解对配置空间进行离散化，实现空间表示的高精度。
• 采用Crank-Nicolson时间离散化，确保二阶精度和无条件稳定性。
• 通过全隐式公式同时求解耦合方程组，实现电荷、动量和能量的精确守恒。
• 采用无雅可比牛顿-Krylov（JFNK）求解器处理隐式时间推进产生的非线性系统，实现高精度收敛。
• 通过使用谱表示而非粒子采样，避免了PIC方法固有的粒子噪声。

实验结果

研究问题

• RQ1 全隐式Hermite谱方法是否能在Vlasov-Poisson模拟中实现电荷、动量和能量的精确守恒？
• RQ2 在标准等离子体基准测试中，全隐式Hermite谱方法的精度和计算成本与全隐式PIC方法相比如何？
• RQ3 Hermite谱方法是否能够解析标准PIC模拟即使在高粒子数下也难以捕捉的弱非线性动力学效应（如高次谐波生成和等离子体回声）？
• RQ4 在无碰撞等离子体的长期模拟中，以单位CPU时间内的误差衡量，Hermite谱方法与PIC方法相比的相对有效性如何？

主要发现

• 在离子声波基准测试中，全隐式Hermite谱方法成功捕捉了主导模式的高阶谐波，而即使使用51,200个粒子/单元的PIC模拟也未能解析这些谐波。
• 在约550秒的CPU时间下，Hermite谱方法的解比表现最佳的PIC解精确约10⁷倍，展现出显著的成本效益优势。
• Hermite谱方法实现了电荷、动量和能量的精确守恒，而PIC方法——即使为全隐式——也未能守恒动量，且需要持续监控。
• 该方法准确解析了两束流不稳定性线性增长速率，并以高保真度捕捉了非线性演化过程，在精度和计算效率上均优于PIC方法。
• 该方法在所有测试案例中均表现出一致的性能优势：朗缪尔波、兰道阻尼、离子声波和两束流不稳定性，证实了其鲁棒性和通用性。
• 结果表明，Hermite谱方法比PIC方法高效数个数量级，即使在PIC使用高粒子数时亦如此，且能够解析标准PIC因噪声限制而无法捕捉的现象，如回声和波-波耦合。