[论文解读] On Runge-Kutta discontinuous Galerkin schemes for Vlasov-Poisson systems

Yingda Cheng, Irene M. Gamba+1 more|arXiv (Cornell University)|Sep 28, 2012

Gas Dynamics and Kinetic Theory参考文献 49被引用 3

一句话总结

本文提出了一种用于求解一维无碰撞等离子体模型的Vlasov-Poisson系统的高阶Runge-Kutta不连续伽辽金（RKDG）格式。该方法保证了守恒性质，通过限制器维持分布函数的正性，并通过傅里叶分析和多项式逼近在基准问题（如Landau阻尼和两束不稳定性）上实现高精度解。

ABSTRACT

In this paper we consider Runge-Kutta discontinuous Galerkin (RKDG) schemes for Vlasov-Poisson systems that model collisionless plasmas. One-dimensional systems are emphasized. The RKDG method, originally devised to solve conservation laws, is seen to have excellent conservation properties, be readily designed for arbitrary order of accuracy, and capable of being used with a positivity-preserving limiter that guarantees positivity of the distribution functions. The RKDG solver for the Vlasov equation is the main focus, while the electric field is obtained through the classical representation by Green's function for the Poisson equation. A rigorous study of recurrence of the DG methods is presented by Fourier analysis, and the impact of different polynomial spaces and the positivity-preserving limiters on the quality of the solutions is ascertained. Several benchmark test problems, such as Landau damping, two-stream instability and the KEEN (Kinetic Electrostatic Electron Nonlinear) wave, are given.

研究动机与目标

• 开发一种用于一维无碰撞等离子体动力学中Vlasov-Poisson系统的高阶、守恒数值格式。
• 通过正性保持限制器确保在粒子动力学模拟中分布函数的正性。
• 通过傅里叶分析研究DG格式的周期性与稳定性特性。
• 在标准基准问题（如Landau阻尼和两束不稳定性）上验证该格式。
• 研究多项式阶数与限制器选择对解质量与精度的影响。

提出的方法

• 将Runge-Kutta不连续伽辽金（RKDG）方法应用于Vlasov方程，实现高阶精度与局部守恒。
• 通过Poisson方程的Green函数表示法精确重构势能，计算电场。
• 将正性保持限制器集成到RKDG框架中，防止分布函数出现负值。
• 采用傅里叶分析研究DG离散化的周期性与谱特性。
• 采用不同阶次的多项式空间，评估在求解动力学方程时的收敛性与精度。
• 在基准问题上开展数值实验，评估稳定性、精度与守恒性质。

实验结果

研究问题

• RQ1 不同多项式空间如何影响RKDG格式在Vlasov-Poisson系统上的精度与稳定性？
• RQ2 正性保持限制器在长时间模拟中在多大程度上保持了物理一致性？
• RQ3 DG离散化中的周期性起什么作用，它如何影响解的质量？
• RQ4 RKDG格式在多大程度上能准确捕捉Landau阻尼与两束不稳定性等关键等离子体现象？
• RQ5 高阶空间与时间离散化对解的精度与守恒性质有何影响？

主要发现

• RKDG格式在空间与时间上均表现出高阶精度，其收敛率与光滑解的理论预期一致。
• 正性保持限制器在所有测试基准问题中成功维持了非负的分布函数。
• 傅里叶分析表明，RKDG方法具有有利的周期性特性，有效抑制了虚假振荡。
• 该方法以极小的数值耗散准确捕捉了Landau阻尼，保持了电场预期的指数衰减特性。
• 在两束不稳定性问题中，该格式能够以高分辨率解析不稳定性的发展与非线性饱和过程。
• KEEN波解被准确模拟，证实了该格式处理复杂非线性动力学现象的能力。