[论文解读] Metriplectic particle-in-cell integrators for the Landau collision operator

Eero Hirvijoki, M. Kraus+1 more|arXiv (Cornell University)|Feb 14, 2018

Gas Dynamics and Kinetic Theory被引用 6

一句话总结

本文提出了一种用于兰道碰撞算子的新型度量-哈密顿粒子-网格积分器，通过宏观粒子权重的有限维度度量-哈密顿系统，保持守恒律与熵产生。通过将宏观粒子离散化与离散梯度相结合，该方法在保持密度、动量与能量代数守恒的同时，确保熵产生为半正定，且与现有的Vlasov-Maxwell泊松积分器完全兼容。

ABSTRACT

In this paper, we present a new framework for addressing the nonlinear Landau collision operator in terms of particle-in-cell methods. We employ the underlying metriplectic structure of the collision operator and, using a macro particle discretization for the distribution function, we transform the infinite-dimensional system into a finite-dimensional time-continuous metriplectic system for advancing the macro particle weights. Temporal discretization is accomplished using the concept of discrete gradients. The conservation of density, momentum, and energy, as well as the positive semi-definite production of entropy in both the time-continuous and the fully discrete system is demonstrated algebraically. The new algorithm is fully compatible with the existing particle-in-cell Poisson integrators for the Vlasov-Maxwell system.

研究动机与目标

• 为粒子-网格方法中兰道算子缺乏结构保持型碰撞算法的问题提供解决方案。
• 解决蒙特卡罗碰撞方法中使用可变粒子权重时导致守恒律不一致的问题。
• 建立一个用于通过兰道碰撞算子演化宏观粒子权重的连续时间、有限维度度量-哈密顿系统。
• 在连续与离散时间下，确保密度、动量与能量的代数守恒，以及熵产生的半正定性。
• 实现与现有结构保持型Vlasov-Maxwell粒子-网格积分器的完全兼容。

提出的方法

• 利用兰道碰撞算子的度量-哈密顿结构，推导出宏观粒子权重的有限维连续时间系统。
• 通过空间与速度空间的形状函数对分布函数进行粒子离散化，实现连续形式与离散形式之间的泛函导数等价。
• 采用离散梯度进行时间离散化，以在完全离散系统中保持度量-哈密顿结构。
• 通过伪逆算子推导离散泛函导数，实现从连续泛函导数到粒子权重导数的精确映射。
• 确保有限维熵保持为离散泊松括号的卡西米尔不变量，从而在离散设置中保持H定理。
• 通过共享的度量-哈密顿结构，将碰撞系统无缝集成至现有Vlasov-Maxwell系统的泊松积分器中。

实验结果

研究问题

• RQ1 能否为兰道碰撞算子构建一个结构保持型粒子-网格积分器，以实现密度、动量与能量的守恒？
• RQ2 如何在宏观粒子权重的有限维系统中对兰道算子的度量-哈密顿结构进行离散化？
• RQ3 能否使用离散梯度在时间离散系统中保持度量-哈密顿结构，以确保熵产生非负？
• RQ4 离散系统中的有限维熵是否仍为卡西米尔不变量，从而确保离散H定理成立？
• RQ5 所得到的碰撞积分器能否与现有结构保持型Vlasov-Maxwell粒子-网格积分器实现完全兼容？

主要发现

• 所提方法在连续时间与完全离散系统中均实现密度、动量与能量的代数守恒。
• 熵产生在连续与离散时间下均被证明为半正定，满足离散H定理。
• 有限维熵泛函在离散泊松括号中仍为卡西米尔不变量，确保长期稳定与物理一致性。
• 通过伪逆算子推导出离散泛函导数，实现从连续泛函导数到粒子权重导数的精确映射。
• 该方法与现有结构保持型Vlasov-Maxwell粒子-网格积分器完全兼容，可实现离子体动力学的统一结构保持模拟。
• 该框架为三步粒子-网格碰撞算法提供了数学上严谨的基础，解决了插值误差与守恒律违反问题。