[论文解读] Conservative discontinuous Galerkin/Hermite Spectral Method for the Vlasov-Poisson System

Francis Filbet, Tao Xiong|arXiv (Cornell University)|Apr 6, 2020

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

一句话总结

本论文提出了一种用于Vlasov-Poisson系统的保守间断伽辽金/埃尔米特谱方法，通过在速度空间使用埃尔米特多项式、在相空间使用间断伽辽金格式离散化，实现了高阶精度，同时严格保证质量与总能量守恒。数值结果验证了高阶收敛性以及在等离子体不稳定性基准测试中关键守恒量的机器精度守恒。

ABSTRACT

We propose a class of conservative discontinuous Galerkin methods for the Vlasov-Poisson system written as a hyperbolic system using Hermite polynomials in the velocity variable. These schemes are designed to be systematically as accurate as one wants with provable conservation of mass and possibly total energy. Such properties in general are hard to achieve within other numerical method frameworks for simulating the Vlasov-Poisson system. The proposed scheme employs discontinuous Galerkin discretization for both the Vlasov and the Poisson equations, resulting in a consistent description of the distribution function and electric field. Numerical simulations are performed to verify the order of accuracy and conservation properties.

研究动机与目标

• 开发一种用于Vlasov-Poisson系统的高阶数值格式，严格保证质量与总能量守恒。
• 克服粒子-网格（PIC）方法的局限性，后者存在随$1/\sqrt{N}$增长的噪声，以及光谱方法在非周期条件下出现的吉布斯振荡问题。
• 结合间断伽辽金方法的优势——局部守恒性与高阶精度——与埃尔米特谱方法，实现高效的速度空间表示。
• 通过一种滤波技术抑制数值振荡，同时不破坏守恒性质，从而在强非线性区域保持数值稳定性。
• 实现对复杂动力学现象（如兰道阻尼、两束流不稳定性和拖尾不稳定性）的精确模拟。

提出的方法

• 将Vlasov-Poisson系统通过在速度变量中使用埃尔米特多项式作为基函数，重述为双曲系统，从而在速度空间实现谱精度。
• 对Vlasov方程和Poisson方程均应用间断伽辽金（DG）方法，确保在物理空间与速度空间中的局部守恒性与高阶精度。
• 构建半离散DG格式，使其在离散层面上保持质量与总能量守恒，通过弱形式与数值通量实现守恒性的理论证明。
• 应用一种滤波技术——具体为Hou-Li滤波器——以抑制分布函数中的虚假振荡，同时不破坏守恒性质。
• 时间积分采用强稳定保形Runge-Kutta方法，以保持空间离散化的守恒性与精度。
• 该方法在1D相空间中实现，采用$N_x \times N_H$个网格点，其中$N_H$对应埃尔米特模态数，并与半拉格朗日有限差分方法（SLFDM）进行对比。

实验结果

研究问题

• RQ1 间断伽辽金/埃尔米特谱方法是否能在Vlasov-Poisson系统中实现高阶精度，同时保证质量与总能量守恒？
• RQ2 所提出的格式在模拟强非线性等离子体不稳定性（如拖尾不稳定性与两束流不稳定性）时表现如何？
• RQ3 该滤波技术在多大程度上能有效抑制虚假振荡，同时不损害守恒性质？
• RQ4 在相同分辨率下，该格式的精度与守恒性与成熟的半拉格朗日有限差分方法（SLFDM）相比如何？
• RQ5 该方法是否能在长时间模拟中保持质量、动量与能量的机器精度守恒？

主要发现

• 所提出的DG/埃尔米特谱方法实现了高阶精度，在所有测试基准中收敛率均与理论方案阶数一致。
• 在所有模拟中，质量、动量与总能量均被守恒至机器精度，包括在强非线性区域（如$\alpha = 0.04$的拖尾不稳定性）中。
• 电场的$L^2$与$L^\infty$范数演化与SLFDM参考解一致，表现出相似的定性与定量行为。
• 总电场能量随时间保持良好，无非物理解增长，证实了该格式的长期稳定性。
• 分布函数的表面图显示，DGFSM能准确捕捉相空间结构，即使在长时间（$T=200$）后也未出现数值发散或过度扩散现象。
• 尽管使用的速度模态数较少（$N_H = 512$），低于SLFDM参考解的$N_v = 320$，DGFSM在相空间结构的分辨率上仍达到相当或更优水平，展现出高效性。