[论文解读] High order and energy preserving discontinuous Galerkin methods for the Vlasov-Poisson system

Blanca Ayuso de Dios, Soheil Hajian|arXiv (Cornell University)|Sep 18, 2012

Gas Dynamics and Kinetic Theory参考文献 43被引用 11

一句话总结

本文提出了一类针对一维Vlasov-Poisson系统的高阶不连续伽辽金（DG）方法的改进家族，该方法在保持高阶精度的同时，精确保持总质量和能量。该方法在相空间中采用不连续多项式逼近，通过弱形式施加边界条件，并采用一种新型通量构造方法以确保能量守恒，其有效性通过离子体物理模拟中的数值基准测试得到验证。

ABSTRACT

We present a computational study for a family of discontinuous Galerkin methods for the one dimensional Vlasov-Poisson system that has been recently introduced. We introduce a slight modification of the methods to allow for feasible computations while preserving the properties of the original methods. We study numerically the verification of the theoretical and convergence analysis, discussing also the conservation properties of the schemes. The methods are validated through their application to some of the benchmarks in the simulation of plasma physics.

研究动机与目标

• 开发高阶不连续伽辽金格式，以在Vlasov-Poisson系统中精确保持关键物理不变量——特别是质量和能量。
• 通过引入一种确定性的、高阶的欧拉型方法（基于DG方法），解决粒子-网格方法中的数值噪声问题。
• 对现有DG格式进行修改，以在保持能量守恒和最优收敛率的同时实现可行的计算效率。
• 在标准离子体物理基准测试上验证该方法，展示其在长时间模拟中的鲁棒性与精度。
• 探索使用伯恩斯坦多项式以高效检测电场的符号变化与极值，从而增强稳定性与正性控制。

提出的方法

• 该方法采用高阶多项式基函数（拉格朗日与伯恩斯坦）在位置和速度变量上进行不连续伽辽金空间离散化。
• 采用半离散DG格式，通过弱形式施加边界条件，构造的数值通量可保证质量和能量守恒。
• 利用伯恩斯坦多项式重构电场，从而通过系数分析实现高效可靠的正性与极值检测。
• 该格式在先前工作的基础上进行了改进，以提升计算效率，特别是通过预计算范德蒙德矩阵的逆矩阵以实现基变换。
• 时间积分采用标准显式时间推进格式，利用块对角质量矩阵结构以提升效率。
• 该方法在分布函数与电场上均保证了最优收敛率，理论分析与数值验证均支持此结论。

实验结果

研究问题

• RQ1 能否设计出在Vlasov-Poisson系统中保持能量守恒的同时维持最优收敛率的高阶不连续伽辽金方法？
• RQ2 如何在不牺牲守恒性质的前提下降低能量守恒DG格式的计算成本？
• RQ3 伯恩斯坦多项式表示在多大程度上提升了DG格式中电场评估的效率与鲁棒性？
• RQ4 改进后的DG格式在标准离子体物理基准测试中的表现如何，其精度与长时间稳定性如何？
• RQ5 基函数选择（拉格朗日与伯恩斯坦）对电场重构条件数与正性保持有何影响？

主要发现

• 所提出的DG方法在$ L^2 $-范数下对分布函数与电场均实现了$ k+1 $阶的最优收敛率。
• 由于采用局部不连续伽辽金格式并弱形式施加边界条件，总质量被精确守恒。
• 经理论证明与数值验证，改进后的格式在长时间模拟中仍能保持系统的总能量守恒。
• 使用伯恩斯坦多项式可实现高效可靠的电场符号变化与极值检测，从而提升数值鲁棒性。
• 通过预计算基变换的范德蒙德矩阵逆矩阵，显著降低了计算成本，使得时间积分过程中可快速执行系数检查。
• 数值基准测试证实了该方法在模拟兰道阻尼与两束不稳定性等关键离子体现象时的精度与稳定性。