[论文解读] Structure preserving schemes for nonlinear Fokker-Planck equations with anisotropic diffusion

Nadia Loy, Mattia Zanella|arXiv (Cornell University)|May 8, 2019

Advanced Thermodynamics and Statistical Mechanics被引用 1

一句话总结

本文提出了一类高阶、结构保持的数值格式，用于具有各向异性扩散的非线性福克-普朗克方程，能够保持解的非负性、熵耗散性以及精确的稳态近似。该格式在瞬态阶段至少达到二阶精度，在长时间阶段实现高阶精度，从而能够稳健模拟集体行为与生命科学模型中的异常扩散。

ABSTRACT

In this work we propose novel numerical schemes for nonlinear Fokker-Planck-type equations with anisotropic diffusion matrix that preserve fundamental structural properties like non negativity of the solution, entropy dissipation and which guarantees an arbitrarily accurate approximation of the steady state of the problem. All the methods presented are at least second order accurate in the transient regimes and high order for large times. Applications of the schemes to models for collective phenomena and life sciences are considered, in these examples anomalous diffusion is often observed and must be taken into account in realistic models.

研究动机与目标

• 开发能够保持关键结构性质——非负性、熵耗散性以及精确稳态近似——的数值格式，用于具有各向异性扩散的非线性福克-普朗克方程。
• 解决在集体现象与生命科学中建模异常扩散的挑战，其中标准格式可能无法保持物理一致性。
• 确保在瞬态阶段至少达到二阶精度，并在长时间阶段实现高阶精度，从而提升长时间模拟的保真度。

提出的方法

• 采用半隐式时间离散化结合有限体积空间离散化，以保持非负性与熵结构。
• 对解变量施加非线性变换，以确保非负性，并在各向异性扩散存在时稳定格式。
• 采用高阶重构技术计算数值通量，实现在时间上二阶精度、在长时间阶段空间上高阶精度。
• 通过从连续方程结构推导出的离散熵不等式，强制实现熵稳定性。
• 该格式在构造上可收敛至正确稳态，且精度可任意提高，即使在复杂的各向异性扩散张量下亦成立。
• 通过在集体动力学与生物系统模型中的应用，验证了该方法的有效性。

实验结果

研究问题

• RQ1 对于具有各向异性扩散的非线性福克-普朗克方程，如何设计数值格式在保持高阶精度的同时确保解的非负性？
• RQ2 哪些离散的结构保持技术可确保数值格式中熵的耗散，从而与连续问题保持一致？
• RQ3 在存在复杂各向异性扩散的情况下，如何实现对稳态的任意精确近似？
• RQ4 该格式在瞬态与长时间阶段中，能在多大程度上保持高阶精度？
• RQ5 所提出的格式能否有效模拟集体行为与生命科学模型中的异常扩散？

主要发现

• 所提出的格式在任意各向异性扩散矩阵下均能保持解的非负性，确保物理一致性。
• 该格式满足离散熵耗散律，复现了连续系统中的热力学行为。
• 在瞬态阶段实现二阶精度，长时间阶段观察到高阶精度，显著提升了长期模拟的可靠性。
• 由于结构设计，格式可收敛至正确稳态，且精度可任意提高，与时间步长无关。
• 在集体运动与生命科学模型中的应用表明，该方法在捕捉异常扩散现象方面具有优异的鲁棒性与精度。