[論文レビュー] Positivity-preserving and energy-dissipative finite difference schemes for the Fokker-Planck and Keller-Segel equations

Jingwei Hu, Xiangxiong Zhang|arXiv (Cornell University)|Mar 31, 2021

Mathematical Biology Tumor Growth参考文献 22被引用数 4

ひとこと要約

本稿では、Fokker-PlanckおよびKeller-Segel方程式に対する、空間方向で2次および4次精度の半陽的有限差分スキームを提示する。これらのスキームは正値性を保ち、エネルギーを散逸させる特性を持つ。4次精度スキームは、時間刻みとメッシュの制限が緩い条件下でも同様の性質を達成でき、正値性とエネルギー散逸の両方を同時に保つ最初の高次空間離散化手法である。

ABSTRACT

In this work, we introduce semi-implicit or implicit finite difference schemes for the continuity equation with a gradient flow structure. Examples of such equations include the linear Fokker-Planck equation and the Keller-Segel equations. The two proposed schemes are first order accurate in time, explicitly solvable, and second order and fourth order accurate in space, which are obtained via finite difference implementation of the classical continuous finite element method. The fully discrete schemes are proven positivity-preserving and energy-dissipative: the second order scheme can achieve so unconditionally; the fourth order scheme only requires a mild time step and mesh size constraint. Furthermore, the fourth order scheme is the first high order spatial discretization that can achieve both positivity and energy decay properties, which is suitable for long time simulation and to obtain accurate steady state solutions.

研究の動機と目的

• Fokker-PlanckやKeller-Segelのような勾配流れ方程式に対して、長期安定性を有する構造保存型数値スキームを開発すること。
• 離散解が非負であり、エネルギーが単調に減少することを保証し、連続的な変分構造を再現すること。
• 正値性とエネルギー散逸を維持する高次空間離散化（2次および4次）を構築すること。時間刻みに厳しい制約が不要である。
• 化学走性および拡散モデルの正確な長期シミュレーションと信頼性の高い定常状態の計算を可能にすること。

提案手法

• 安定性を確保するために、時間方向に1次精度の半陽的スキームを用い、系を完全に陰的または半陽的に解く。
• 均等メッシュ上での線形および2次多項式を用いた連続有限要素法の有限差分近似により、2次および4次空間精度を達成する。
• スキームは、連続方程式の等価な係数付き拡散形 ∂ₜρ = ∇·(ℳ∇(ρ/ℳ)) から導出される。ここで ℳ = exp(−(𝒱 + 𝒲∗ρ)) である。
• 正値性は、離散スキームにおける単調性条件によって強制され、離散最大原理の議論により証明される。
• エネルギー散逸は、離散自由エネルギー関数が時間とともに減少することを示すことで確立される。これは連続的なエネルギー減少を模倣する。
• 4次精度スキームでは、単調性を保つために時間刻みに下限（Δtの下限）が必要であり、上限制約は不要である。

実験結果

リサーチクエスチョン

• RQ1 Fokker-PlanckおよびKeller-Segel方程式に対して、高次有限差分スキームが正値性とエネルギー散逸の両方を保存できるか。
• RQ2 高次スキームの単調性と正値性を保証するためには、時間および空間離散化にどのような制約が必要か。
• RQ3 提案されたスキームは、長期シミュレーションおよび爆発付近の状態でどのように性能を発揮するか。
• RQ4 メッシュ制約を緩和しても、4次精度スキームが正値性とエネルギー散逸を維持できるか。

主な発見

• 2次精度スキームは、時間刻みに制限がなく、正値性を保ち、エネルギーを散逸させる。したがって、制限なしに安定性が保証される。
• 4次精度スキームは、時間刻みが下限に制限されるという緩い制約のもとで、正値性とエネルギー散逸を達成する。
• 数値結果から、両スキームとも、臨界以下初期質量の場合、T ≈ 13.52 で ∥ρⁿ⁺¹ − ρⁿ∥∞ ≤ 10⁻⁸ に収束することが示された。
• 臨界を超える初期質量の場合、4次精度スキームは2次精度スキームよりも、特にT = 0.8における爆発ダイナミクスをより正確に捉えている。
• エネルギーの時間発展プロットから、両スキームとも爆発後も単調な減少を示しており、エネルギー散逸や正値性の違反は認められない。
• 4次精度スキームは、正値性とエネルギー散逸の両方を同時に保存する最初の高次空間離散化手法であり、正確な長期シミュレーションを可能にする。