[论文解读] Homogenization of the Vlasov Equation and of the Vlasov - Poisson System with a Strong External Magnetic Field

Emmanuel Frénod, Éric Sonnendrücker|arXiv (Cornell University)|Nov 18, 2010

Advanced Mathematical Modeling in Engineering参考文献 12被引用 67

一句话总结

该论文利用双尺度收敛方法，为强外部磁场下的Vlasov方程和Vlasov-Poisson方程建立了严格的均质化框架。证明了粒子分布收敛到一个描述导心运动的极限方程，其中沿磁感线的快速回旋运动被平均化，而自洽电场仍与平均粒子运动保持一致。

ABSTRACT

Motivated by the difficulty arising in the numerical simulation of the movement of charged particles in presence of a large external magnetic field, which adds an additional time scale and thus imposes to use a much smaller time step, we perform in this paper a homogenization of the Vlasov equation and of the Vlasov-Poisson system which yield approximate equations describing the mean behavior of the particles. The convergence proof is based on the two scale convergence tools introduced by N'Guetseng and Allaire. We also consider the case where, in addition to the magnetic field, a large external electric field orthogonal to the magnetic field and of the same magnitude is applied.

研究动机与目标

• 为粒子模拟中强外部磁场引入的刚性时间尺度带来的数值挑战提供解决方案。
• 推导一个简化且渐近一致的系统，以捕捉带电粒子的平均行为，同时滤除快速回旋运动。
• 严格证明Vlasov-Poisson系统中的粒子相互作用可完全用导心运动表示，无需额外校正项。
• 将经典的导心近似从单粒子动力学推广到Vlasov-Poisson系统的动力学非线性情形。

提出的方法

• 应用N’Guetseng和Allaire提出的双尺度收敛技术，分析当磁场强度ε⁻¹趋于无穷大时Vlasov方程的极限。
• 考虑具有强恒定外部磁场B^ε = B + M/ε的Vlasov方程，其中M为单位向量。
• 利用L∞(0,T;L²(Ω))中的弱-*收敛和双尺度收敛，识别出极限分布F(t,x,v,τ)，其中τ代表快速回旋相位。
• 通过在快速时间尺度τ上平均，推导出有效方程，使得速度v被其平行分量v∥所取代，动力学由导心运动控制。
• 通过Aubin-Lions引理，利用泊松方程的正则性以及电荷和电流密度的有界性，证明电场E^ε在L∞(0,T;L²_loc(R³_x))中强收敛。
• 证明电荷密度和电流密度也实现双尺度收敛，且电荷密度的极限与快速时间尺度τ无关，从而确保与平均场近似的自洽性。

实验结果

研究问题

• RQ1 能否在Vlasov-Poisson系统中，通过均质化方法严格平均掉强磁场下带电粒子的快速回旋运动？
• RQ2 Vlasov-Poisson系统中粒子间的非线性耦合在极限下是否保持导心近似，而无需额外项？
• RQ3 当磁场强度趋于无穷大时，自洽电场的行为如何？
• RQ4 双尺度收敛技术能否用于推导Vlasov-Poisson系统的自洽磁流体力学型方程？
• RQ5 为确保电场在极限下强收敛，电场和粒子分布需要具备何种正则性和收敛性质？

主要发现

• 强磁场下Vlasov方程的解f^ε在L∞(0,T;L²(Ω))中弱-*收敛于一个极限f，该极限满足仅包含速度平行分量v∥的简化方程。
• 极限分布f的初始条件通过沿磁场方向快速旋转的初始数据f₀的平均给出，即f(0) = (1/(2π)) ∫₀²π f₀(x, u(v,τ)) dτ。
• 由于在W¹,⁷/⁵中有界且时间导数属于L∞(0,T;L⁷/⁶)，通过Aubin-Lions引理的紧嵌入，电场E^ε在L∞(0,T;L²_loc(R³_x))中强收敛。
• 电荷密度ρ^ε双尺度收敛于一个与快速时间尺度τ无关的极限ρ̄，意味着平均密度中无振荡校正项。
• 电流密度J^ε双尺度收敛于一个与τ无关的极限J̄，确保了平均场动力学的一致性。
• 极限f的最终有效方程为导心方程：∂f/∂t + v∥·∇ₓf + (E∥ + v×B∥)·∇ᵥf = 0，该方程描述了沿磁感线的平均运动，并已平均化回旋效应。