条件Vlasov–Ampère系统的渐近保持兼容离散化及其在等离子体物理中的应用

《Journal of Plasma Physics》：A novel conditional formulation of the Vlasov–Ampère equations: a conservative, positivity, asymptotic and Gauss law preserving scheme

　　

研究人员提出了一种针对条件Vlasov–Ampère系统的新型数值离散方法，该方法同时满足质量、动量和能量守恒，并在准中性极限下保持渐近一致性。通过将分布函数变换至无量纲速度变量w，并结合内部能量形式的流体矩方程，避免了总能量形式在非交错网格上的插值误差。研究采用隐式–显式(IMEX)时间积分策略，其中Vlasov方程采用显式二阶Runge–Kutta方法，矩与场方程采用隐式中点法，以确保离散守恒性。空间离散采用交错网格布局，标量量位于单元中心，矢量量位于单元界面，从而自然地保持离散高斯定律。数值通量采用单调保持的SMART重构及动能通量型格式，确保激波捕捉能力与单调性。理论证明表明，该格式在ε→0时收敛至准中性Ohm定律，并存在离散慢流形，从而实现从动力学尺度到流体尺度的光滑过渡。该工作为等离子体多尺度模拟提供了一种兼具守恒性、稳定性与渐近正确性的数值框架。

该研究发表于《Journal of Plasma Physics》，针对等离子体多尺度模拟中存在的数值刚性与守恒性难以兼顾的问题，提出了一种基于条件Vlasov–Ampère系统的渐近保持（Asymptotic-Preserving, AP）离散格式。传统全动力学Vlasov模拟在准中性极限下会因电子时间步长过小而面临严重的计算瓶颈，而现有简化模型往往无法同时保证质量、动量与能量的严格守恒，也难以在离散层面保持高斯定律。为此，研究人员将Vlasov方程通过无量纲速度变量w进行变换，得到分布函数??，使其前三个速度矩固定为不变量，从而避免直接从??计算流体矩，转而通过独立的矩方程演化密度、速度与温度。这一变换使得快尺度电子惯性项被显式分离，形成快慢变量系统，从而在ε→0时自然退化为准中性Ohm定律。

关键技术方法包括：（1）采用交错网格布局，标量量（密度、温度、分布函数）置于单元中心，矢量量（速度、电流、电场）置于单元界面，确保离散高斯定律自动满足；（2）引入隐式–显式（IMEX）时间积分，Vlasov方程采用显式二阶Runge–Kutta，矩与场方程采用隐式中点法，兼顾刚度处理与守恒性；（3）空间通量采用SMART单调重构与动能通量型格式，保证激波捕捉与单调性；（4）通过非Ampère电流方程替代电子动量方程，构造离散慢流形，实现准中性极限下的稳定推进。

研究结果如下：

4. 快慢变量形式

研究人员将电子动量方程替换为非Ampère电流方程，与Ampère定律共同构成快变量子系统，其余为慢变量。在ε→0极限下，快变量退化为代数关系，给出准中性条件与Ohm定律，并证明离散系统在单离子情形下存在唯一非退化不动点，从而形成离散慢流形。

5. 兼容渐近保持离散格式

提出的交错网格有限差分格式严格满足质量、动量与能量守恒，并通过数值通量的特殊构造保证激波稳定性。理论证明显示，离散高斯定律在所有时间层自动保持，且在ε→0时收敛至准中性Ohm定律。

5.1 高斯定律保持

通过电荷连续性方程与离散散度运算的结合，证明若初始场满足高斯定律，则所有后续时间层均满足，且在准中性极限下电荷密度趋于零。

5.2 质量守恒

由于连续性方程的保守形式与通量离散的一致性，总质量在时间推进中严格保持不变。

5.3 动量守恒

通过引入等效电子动量表达式，结合离散高斯定律与电场散度项，证明总动量（含离子与等效电子贡献）在时间上守恒。

5.4 能量守恒

将内部能量方程与动能方程结合，利用离散高斯定律消去电场功项，证明总能量（等离子体动能、热能、电场能）严格守恒。

5.5 离散慢流形的存在性

在线性重构与单离子情形下，证明离散系统在Δt→0与ε/Δt→0极限下存在光滑映射Φ_Δt,ε/Δt，其不动点对应连续慢流形，从而保证数值格式在准中性极限下的渐近正确性。

讨论部分指出，该格式避免了传统总能量形式在非交错网格上的插值不一致问题，同时通过内部能量方程与辅助动能方程的耦合，恢复了离散能量守恒。研究结论表明，所提出的AP格式能够在宽参数范围内稳定模拟从动力学到流体尺度的等离子体行为，并在准中性极限下保持物理一致性，为未来等离子体多尺度模拟提供了可靠数值工具。