[论文解读] Energy and implicit discretization of the Fokker-Planck and Keller-Segel type equations

Luís Almeida, Federica Bubba+2 more|arXiv (Cornell University)|Mar 28, 2018

Mathematical Biology Tumor Growth参考文献 17被引用 5

一句话总结

本文提出了两种针对Fokker-Planck方程和Keller-Segel方程的新型隐式有限体积格式，可在离散层面保持能量耗散、质量守恒、解的正性以及稳态结构。通过结合梯度流方法与Scharfetter-Gummel离散化策略并引入迎风格式，该方法确保了无时间步长限制的适定性与稳定性，从而能够准确模拟模式形成过程，并有效区分物理不稳定性与数值伪影。

ABSTRACT

The parabolic-elliptic Keller-Segel equation with sensitivity saturation, because of its pattern formation ability, is a challenge for numerical simulations. We provide two finite-volume schemes whose goals are to preserve, at the discrete level, the fundamental properties of the solutions, namely energy dissipation, steady states, positivity and conservation of total mass. These requirements happen to be critical when it comes to distinguishing between discrete steady states, Turing unstable transient states, numerical artifacts or approximate steady states as obtained by a simple upwind approach. These schemes are obtained either by following closely the gradient flow structure or by a proper exponential rewriting inspired by the Scharfetter-Gummel discretization. An interesting feature is that upwind is also necessary for all the expected properties to be preserved at the semi-discrete level. These schemes are extended to the fully discrete level and this leads us to tune precisely the terms according to explicit or implicit discretizations. Using some appropriate monotony properties (reminiscent of the maximum principle), we prove well-posedness for the scheme as well as all the other requirements. Numerical implementations and simulations illustrate the respective advantages of the three methods we compare.

研究动机与目标

• 开发能够保持Fokker-Planck方程和Keller-Segel方程基本物理性质（能量耗散、质量守恒、正性及稳态结构）的数值格式。
• 解决由于强非线性项与对流项导致的模式形成模拟中的数值挑战，特别是区分物理Turing不稳定性与数值伪影。
• 构建隐式时间格式，避免CFL条件限制，同时保持所有关键离散性质。
• 比较基于梯度流与基于Scharfetter-Gummel的格式在准确捕捉复杂模式方面的性能。
• 利用次解-超解理论证明格式的适定性与基于单调性的收敛性。

提出的方法

• 该格式基于两种不同的对称化策略：梯度流结构与漂移-扩散方程的Scharfetter-Gummel方法。
• 在两种格式中均引入迎风格式，以确保在半离散层面保持正性、质量守恒与能量耗散。
• 通过隐式时间推进构建全离散格式，确保无时间步长限制且无条件稳定。
• 利用类似于最大值原理的单调性性质，通过次解-超解理论证明解的存在性与唯一性。
• 通过构造一个与连续能量不等式相对应的李雅普诺夫泛函，严格证明了离散能量耗散。
• 通过数值实验验证了格式性能，结果表明其在解析稳态解与瞬态模式方面优于标准迎风格式，具有更高的精度。

实验结果

研究问题

• RQ1 隐式有限体积格式能否在Fokker-Planck方程和Keller-Segel方程中同时保持能量耗散、质量守恒、正性及稳态结构？
• RQ2 基于梯度流与基于Scharfetter-Gummel的离散化在处理强非线性问题时，对解的定性行为保持效果如何比较？
• RQ3 迎风格式在半离散层面能多大程度上确保关键物理性质的保持？
• RQ4 隐式时间离散化能否在保持所有离散守恒与稳定性性质的同时消除CFL条件限制？
• RQ5 在这些系统中，如何区分数值伪影（如虚假振荡或错误稳态）与真实的模式形成？

主要发现

• 所提出的格式在半离散与全离散层面均保持能量耗散，能量泛函随时间单调递减。
• 两种格式均严格保持解的正性并守恒总质量，即使在非光滑或尖峰状初始数据下亦成立。
• 格式无条件稳定且适定，解的存在性与唯一性通过次解-超解理论得到证明。
• 基于Scharfetter-Gummel的格式为梯度流方法提供了稳健替代方案，尤其在一般敏感函数φ(u)下表现优异。
• 数值模拟表明，该格式能准确捕捉Turing型不稳定性与类峰状稳态解，避免了标准迎风格式中常见的虚假振荡。
• 隐式时间离散化允许采用大时间步长而不损失稳定性或精度，使该方法适用于长期模式形成模拟。