针对正常或强磁场中带电粒子动力学的滤波二步变分积分器
A filtered two-step variational integrator for charged-particle dynamics in a normal or strong magnetic field
摘要 Abstract
本文提出了一种新的针对正常或强磁场中略非均匀磁场下带电粒子动力学问题的滤波二步变分积分器,其中磁场强度由无量纲参数$\epsilon$反比例表示。对于正常磁场情形($\epsilon \approx 1$),该方法获得了二阶误差界以及长时间的能量和动量守恒性,并通过后向误差分析完成了长时间分析的证明。对于强磁场情形($0<\epsilon \ll 1$),本文探讨了滤波变分积分器在大步长$h^2 \geq \epsilon$和小步长$h \sim \epsilon$下的行为。误差界分析基于精确解与数值解的调制傅里叶展开比较。研究表明,对于大步长该积分器在位置和并行速度上的精度达到$\mathcal{O}(h^2)$,而对于小步长的精度为$\mathcal{O}(\epsilon)$。此外,本文还通过发展所提方案的调制傅里叶展开,得到了强磁场条件下的长时间能量和磁矩守恒性。所有理论结果均通过两个数值实验进行了验证。
This article is concerned with a new filtered two-step variational integrator for solving the charged-particle dynamics in a mildly non-homogeneous normal or strong magnetic field with a dimensionless parameter $\epsilon$ inversely proportional to the strength of the magnetic field. In the case of a normal magnetic field ($\epsilon \approx 1$), second-order error bounds and long time energy and momentum conservations are obtained. Moreover, the proof of the long-term analysis is accomplished by the backward error analysis. For the strong magnetic field ($0<\epsilon \ll1$), this paper clarifies the behaviour of the filtered variational integrator for both a large stepsize $h^2 \geq \epsilon$ and a smaller stepsize $ h \sim \epsilon$. The approach to analysing the error bounds for these two stepsizes is based on comparing the modulated Fourier expansions of the exact and the numerical solutions. It is shown that the proposed integrator achieves a second-order accuracy $\mathcal{O}(h^2)$ in the position and in the parallel velocity for a large step size and an $\mathcal{O}(\epsilon)$ accuracy for a smaller stepsize. This paper also yields the long time energy and magnetic moment conservations for the strong magnetic field by developing the modulated Fourier expansion of the proposed scheme. All the theoretical results of the error behaviour and long-term conservations are numerically demonstrated by two numerical experiments.