Vlasov–Poisson 方程的低秩投影分裂积分器
数值分析
2018-06-12 v2 计算物理
摘要
等离子体物理中遇到的许多问题需要用动力学方程描述,这些方程定义在高达六维的相空间中。对该相空间的直接离散化(常称为欧拉方法)具有诸多优点,但从计算角度看极其昂贵。本文中,我们针对 Vlasov–Poisson 方程提出一种动态低秩近似,并采用特定的分裂方法进行时间积分。该近似通过将动力学约束到低秩函数流形(经由切空间投影)并将此投影拆分为构成它的子投影而导出。这将六维(或四维)Vlasov–Poisson 方程的一个时间步简化为在该时间步内求解两组三维(或二维)平流方程,一组在位置变量中、一组在速度变量中,每组平流方程的大小等于所选秩。通过分层动态低秩近似,Vlasov–Poisson 方程的一个时间步可进一步简化为六(或四)组一维平流方程,每组平流方程的大小仍等于秩。所得平流方程组随后可用半拉格朗日或谱方法等标准技术求解。在二维和四维下对线性朗道阻尼、双流不稳定性及等离子体回波问题的数值模拟凸显了该数值方法的良好特性,并表明所提算法能够大幅降低所需的计算代价。
引用
@article{arxiv.1801.01103,
title = {A low-rank projector-splitting integrator for the Vlasov--Poisson equation},
author = {Lukas Einkemmer and Christian Lubich},
journal= {arXiv preprint arXiv:1801.01103},
year = {2018}
}