[论文解读] Contour Dynamics for One-Dimensional Vlasov-Poisson Plasma with the Periodic Boundary

Hiroki Sato, T. Watanabe+1 more|arXiv (Cornell University)|Feb 1, 2021

Fluid Dynamics and Turbulent Flows参考文献 13被引用 3

一句话总结

本文提出了一种将轮廓动力学（CD）方法新颖地应用于具有周期性边界条件的一维Vlasov-Poisson离子体的方法，利用周期性格林函数处理跨越边界的轮廓，避免了切割或重新分配节点。该方法成功模拟了线性朗道阻尼及朗缪尔波引起的非线性电子捕获现象，能量和粒子数守恒误差分别低于2.5×10⁻⁵和10⁻⁴，表明其在无网格限制条件下对长时间动力学模拟具有强健性。

ABSTRACT

We revisit the contour dynamics (CD) simulation method which is applicable to large deformation of distribution function in the Vlasov-Poisson plasma with the periodic boundary, where contours of distribution function are traced without using spatial grids. Novelty of this study lies in application of CD to the one-dimensional Vlasov-Poisson plasma with the periodic boundary condition. A major difficulty in application of the periodic boundary is how to deal with contours when they cross the boundaries. It has been overcome by virtue of a periodic Green's function, which effectively introduces the periodic boundary condition without cutting nor reallocating the contours. The simulation results are confirmed by comparing with an analytical solution for the piece-wise constant distribution function in the linear regime and a linear analysis of the Landau damping. Also, particle trapping by Langmuir wave is successfully reproduced in the nonlinear regime.

研究动机与目标

• 将轮廓动力学（CD）方法扩展至具有周期性边界条件的一维Vlasov-Poisson离子体系统，其中轮廓可穿越系统边界。
• 克服轮廓穿越周期性边界时无法切割或重新分配节点的挑战。
• 在线性区域中，通过与解析解对比验证CD方法的准确性，并与朗道阻尼的基准结果进行比较。
• 展示该方法在模拟非线性现象（如朗缪尔波引起的电子捕获）方面的能力。
• 确保在长时间模拟中保持能量与粒子数的高守恒性，且无需节点重分布或轮廓手术操作。

提出的方法

• 该方法采用定义在无限相空间上的周期性格林函数，以在不改变轮廓拓扑结构的前提下施加周期性边界条件。
• 分布函数的轮廓由离散节点表示，其运动通过格林函数的线积分所导出的哈密顿方程计算。
• 电场和电势通过满足∇²G(x;ξ) = 1/L − δ(x−ξ)的周期性格林函数计算，从而在周期性区域内一致求解泊松方程。
• 通过使用由格林函数和f的类似涡量跳跃所导出的速度场，推进节点位置来计算分布函数的时间演化。
• 未使用轮廓手术或节点重分布，从而在整个模拟过程中保持轮廓表示的完整性。
• 该方法利用了李维维尔定理和不可压缩性，确保相空间中体积与粒子数的守恒。

实验结果

研究问题

• RQ1 如何将轮廓动力学方法适配至具有周期性边界条件的一维Vlasov-Poisson系统，特别是当轮廓穿越域边界时？
• RQ2 周期性格林函数是否能有效实现周期性边界条件，而无需切割或重新分配轮廓？
• RQ3 CD方法是否能准确再现具有分段常数初始f的一维Vlasov-Poisson系统中的线性朗道阻尼？
• RQ4 该方法是否能在无数值伪影的情况下捕捉朗缪尔波引起的非线性电子捕获动力学？
• RQ5 在该CD实现中，能量与粒子数在长时间模拟中的守恒程度如何？

主要发现

• CD方法成功模拟了线性朗道阻尼，与解析解高度一致，验证了其在线性区域的准确性。
• 该方法在非线性区域中准确再现了朗缪尔波引起的电子捕获现象，证实了其对复杂相空间动力学的模拟能力。
• 总能量守恒误差在整个模拟过程（τ = 30）中保持在2.5×10⁻⁵以下，表明具有高度数值稳定性。
• 粒子（或面积）守恒误差随时间保持在10⁻⁴以下，且仅缓慢增长，表明即使在轮廓发生形变时仍具有强健性。
• 使用周期性格林函数消除了轮廓切割或节点重分配的需要，实现了高效且拓扑结构保持的模拟。
• 该方法保持了高分辨率并避免了网格限制，适用于具有强烈相空间形变的系统。