[论文解读] Entropy Decay Rates for Conservative Spectral Schemes Modeling Fokker-Planck-Landau Type Flows in the Mean Field Limit

Clark A. Pennie, Irene M. Gamba|arXiv (Cornell University)|Oct 7, 2019

Gas Dynamics and Kinetic Theory参考文献 21被引用 4

一句话总结

本文提出了一种保守的混合谱-不连续伽辽金格式，用于模拟描述库仑、麦克斯韦和硬球粒子相互作用的福克-普朗克-兰道方程中的熵衰减。结果表明，库仑情况表现出退化的谱，其衰减速率接近理论上的三分之二幂律，而麦克斯韦和硬球情况则显示出谱隙，证实了分析预测，并验证了库仑情况估计的紧致性。

ABSTRACT

The focus of this work is to create benchmark simulations of decay rates to statistical equilibrium in transport plasma models for Coulomb particle interactions given by a coupled Vlasov-Poisson Fokker-Planck-Landau equation, as well as with Maxwell type and hard sphere interactions. The qualitative decay to the equilibrium Maxwell-Boltzmann distribution through relative entropy is studied in detail for all three types of particle interactions by means of a conservative hybrid spectral and discontinuous Galerkin scheme adapted from previous work. More precisely, the Coulomb case shows that there is a degenerate spectrum, with a decay rate close to the law of two thirds predicted by upper estimates in a work of Strain and Guo in 2006, while the Maxwell type and hard sphere examples both exhibit a spectral gap as predicted by Desvillettes and Villani in 2000. Such decay rate behavior indicates that the analytical estimates for the Coulomb case is sharp while, still to this date, there is no analytical proof of sharp degenerate spectral behaviour for the Fokker-Planck-Landau operator. Simulations are presented, both for the space-homogeneous case of just particle potential interactions and the space-inhomogeneous case for the mean field coupling through the Poisson equation for total charges in periodic domains. New explicit derivations of spectral collisional weights are presented in the case of Maxwell type and hard sphere interactions and the stability of all three scenarios, including Coulomb interactions, is investigated.

研究动机与目标

• 开发一种保守的数值格式，用于模拟平均场极限下维拉索-泊松福克-普朗克-兰道方程的长时间动力学。
• 研究在三种粒子相互作用类型（库仑、麦克斯韦和硬球）下，相对熵向麦克斯韦-玻尔兹曼平衡态的衰减速率。
• 验证分析熵衰减速率的紧致性，特别是库仑情况下退化谱的紧致性。
• 通过最小化过程强制保持碰撞不变量，确保长时间模拟中的数值精度和稳定性。
• 将先前针对空间均匀情况的结果扩展至具有周期性域中泊松耦合的空间非均匀情况。

提出的方法

• 采用一种保守的混合谱方法与不连续伽辽金方法，结合谱方法计算碰撞算子，DG方法处理对流项。
• 谱方法在 O(N³log N) 次操作内计算碰撞算子的傅里叶变换，通过 O(N³) 的校正步骤强制保持守恒。
• 基于 Gamba 等人（2017）对玻尔兹曼方程的质量浓度结果，施加速度空间截断。
• 给出了麦克斯韦型和硬球相互作用的谱碰撞权重的显式推导。
• 该格式保持了总质量、动量和能量守恒，确保长时间稳定性及正确的平衡态动力学。
• 在空间均匀和空间非均匀设置下均进行了数值模拟，包括与泊松方程耦合以实现自洽电场。

实验结果

研究问题

• RQ1 该保守谱格式能否准确捕捉福克-普朗克-兰道方程在库仑相互作用下的熵向平衡态的衰减？
• RQ2 在平均场极限下，库仑、麦克斯韦型和硬球相互作用的熵衰减速率如何比较？
• RQ3 库仑情况下预测的三分之二幂律衰减速率是否在数值上得到验证，是否表明分析估计的紧致性？
• RQ4 麦克斯韦型和硬球相互作用是否如 Desvillettes 和 Villani（2000）所预测的那样在熵衰减中表现出谱隙？
• RQ5 该格式能否在长时间模拟中保持不变量守恒，确保正确趋近于平衡态？

主要发现

• 库仑情况表现出退化的谱，熵衰减速率接近理论上的三分之二幂律，支持了 Strain 和 Guo（2006）所得上界估计的紧致性。
• 对于麦克斯韦型和硬球相互作用，模拟结果证实了熵衰减中存在谱隙，与 Desvillettes 和 Villani（2000）的预测一致。
• 保守谱格式在长时间模拟中保持了质量、动量和能量守恒，从而实现了对麦克斯韦-玻尔兹曼平衡态的准确逼近。
• 该方法在空间均匀和空间非均匀情况下均实现了长达 400 个时间单位的稳定且精确的动力学，包括泊松耦合。
• 显式推导并数值验证了麦克斯韦和硬球相互作用的谱权重，增强了格式的适用性。
• 数值结果表明，尽管在库仑情况下数值上已确认了其紧致性，但福克-普朗克-兰道算子退化谱行为的分析证明仍是开放问题。