[论文解读] A Semi-Lagrangian Spectral Method for the Vlasov-Poisson System based on Fourier, Legendre and Hermite Polynomials

Lorella Fatone, Daniele Funaro+1 more|arXiv (Cornell University)|Jul 5, 2018

Gas Dynamics and Kinetic Theory参考文献 74被引用 8

一句话总结

本文通过在速度空间中引入勒让德多项式和埃尔米特多项式，并结合二阶BDF时间推进格式，将半拉格朗日谱方法扩展至Vlasov-Poisson系统。该方法实现了高精度，并表现出优异的质量、动量和能量守恒性能；数值结果表明，在埃尔米特情形下，最优性能在很大程度上依赖于高斯权重函数中缩放参数的调优。

ABSTRACT

In this work, we apply a semi-Lagrangian spectral method for the Vlasov-Poisson system, previously designed for periodic Fourier discretizations, by implementing Legendre polynomials and Hermite functions in the approximation of the distribution function with respect to the velocity variable. We discuss second-order accurate-in-time schemes, obtained by coupling spectral techniques in the space-velocity domain with a BDF time-stepping scheme. The resulting method possesses good conservation properties, which have been assessed by a series of numerical tests conducted on the standard two-stream instability benchmark problem. In the Hermite case, we also investigate the numerical behavior in dependence of a scaling parameter in the Gaussian weight. Confirming previous results from the literature, our experiments for different representative values of this parameter, indicate that a proper choice may significantly impact on accuracy, thus suggesting that suitable strategies should be developed to automatically update the parameter during the time-advancing procedure.

研究动机与目标

• 将半拉格朗日谱方法从傅里叶-傅里叶离散化扩展至速度空间中包含勒让德和埃尔米特谱近似的Vlasov-Poisson系统。
• 通过将谱空间-速度离散化与BDF时间推进方法结合，发展二阶时间精度的格式。
• 利用两流体不稳定性基准测试，评估新离散化下的守恒性质，特别是质量、动量和能量的守恒性能。
• 研究傅里叶-埃尔米特方法对高斯权重函数中缩放参数α的敏感性及其对数值精度的影响。
• 为未来发展动态调整时间积分过程中α的自适应策略奠定基础。

提出的方法

• 该方法在空间使用傅里叶多项式，在速度空间使用勒让德或埃尔米特函数进行谱近似，从而实现高阶空间-速度分辨率。
• 采用二阶BDF时间推进格式以实现时间推进，确保数值稳定性和精度。
• 半拉格朗日方法通过时间反向追踪特征线计算解，利用任意阶数的泰勒展开进行插值。
• 空间和速度方向的谱导数以谱收敛速率高效计算，保持精度。
• 埃尔米特基函数加权于exp(−α²v²)，通过改变参数α以研究其对收敛性和稳定性的影响。
• 通过从展开系数导出的离散不变量，实现质量、动量和能量守恒定律的强制满足。

实验结果

研究问题

• RQ1 在速度空间中采用勒让德和埃尔米特谱近似，如何影响半拉格朗日方法在Vlasov-Poisson系统中的精度和守恒性质？
• RQ2 埃尔米特基函数权重函数exp(−α²v²)中的缩放参数α对数值解精度有何影响？
• RQ3 二阶BDF时间推进与谱空间-速度离散化结合，能否保持质量、动量和能量的良好守恒性？
• RQ4 在能量守恒和解质量方面，傅里叶-埃尔米特方法相较于傅里叶-傅里叶方法表现如何？
• RQ5 是否需要发展自适应策略，在时间积分过程中动态调整α以维持最优精度？

主要发现

• 傅里叶-勒让德和傅里叶-埃尔米特方法在所有数值测试中均表现出优异的质量和动量守恒性能，未观察到显著漂移。
• 能量守恒略有削弱，表现为依赖于时间步长的微小增长率，但BDF格式相比一阶方法提升了稳定性。
• 对于傅里叶-埃尔米特方法，α的选择显著影响精度：当α=1时达到最优性能，而偏离（如α=0.9或1.2）则导致收敛性下降。
• 在相同分辨率（N=M=2⁴）下，α=1的傅里叶-埃尔米特方法在能量守恒和解质量方面优于傅里叶-傅里叶方法。
• 数值结果证实，合理选择α可显著提升收敛性，表明自适应参数控制对实现鲁棒性能至关重要。
• 该方法在空间和速度方向均保持谱精度，半拉格朗日方法通过时间反向特征追踪实现稳定的时推进。