[论文解读] Finite difference method for a fractional porous medium equation

Félix del Teso|arXiv (Cornell University)|Jan 18, 2013

Fractional Differential Equations Solutions参考文献 7被引用 5

一句话总结

本文提出了一种有限差分法，用于求解带有 1/2 阶分数阶拉普拉斯算子的分数阶多孔介质方程，利用 Caffarelli-Silvestre 扩展将非局部问题转化为高一维中的局部问题。该方法证明了解的存在性、唯一性及对真实解的收敛性，数值实验验证了二阶收敛率，并准确捕捉了重尾解。

ABSTRACT

We formulate a numerical method to solve the porous medium type equation with fractional diffusion \[\frac{\partial u}{\partial t}+(-Δ)^{1/2} (u^m)=0.\] The problem is posed in $x\in \mathbb{R}^N$, $m\geq 1$ and with nonnegative initial data. The fractional Laplacian is implemented via the so-called Caffarelli-Silvestre extension. We prove existence and uniqueness of the solution of this method and also the convergence to the theoretical solution of the equation. We run numerical experiments on typical initial data.

研究动机与目标

• 开发一种针对具有 (-Δ)^1/2 扩散的分数阶多孔介质方程的鲁棒数值格式。
• 通过使用 Caffarelli-Silvestre 扩展将分数阶扩散中的非局部性问题转化为高一维中的局部问题，以应对该挑战。
• 利用最大值原理，证明有限差分解的存在性、唯一性及收敛性。
• 通过针对不同初始数据（包括狄拉克δ函数和紧支集函数）的数值实验验证该方法。
• 将方法扩展至一般单调非线性项 φ(u)，超越 u^m 的情形，包括对数扩散作为极限情况。

提出的方法

• 通过 Caffarelli-Silvestre 扩展实现分数阶拉普拉斯算子，将非局部问题转化为 R^N × R^+ 中的局部边值问题。
• 在有界区域 Ω = (-X,X)^N × (0,Y) 内数值求解扩展问题，侧边界 Γ_h 上施加齐次狄利克雷条件。
• 对扩展方程应用有限差分格式，采用半隐式时间离散化及空间上的中心差分。
• 时间演化通过单调格式处理，以保持解的正性并满足离散最大值原理。
• 通过引入 φ 及其逆函数，将时间导数条件重新表述，从而将方法推广至一般 φ(u)。
• 通过离散最大值原理与稳定性估计证明收敛性，误差分析基于精细网格与 m=1 时的精确解。

实验结果

研究问题

• RQ1 能否为具有 (-Δ)^1/2 扩散的分数阶多孔介质方程构造一种有限差分法，以确保收敛性与稳定性？
• RQ2 Caffarelli-Silvestre 扩展如何实现对非局部问题的实际数值求解？
• RQ3 对于不同的 m 值与初始数据，所提出的有限差分格式的收敛速率如何？
• RQ4 该方法能否扩展至一般非线性项 φ(u)，包括对数扩散作为极限情形？
• RQ5 该方法在多大程度上准确捕捉了分数阶扩散下多孔介质方程所特有的重尾解行为？

主要发现

• 有限差分法通过离散最大值原理，证明了离散解的存在性与唯一性。
• 当网格尺寸 Δx 与时间步长 Δt 趋近于零时，该方法收敛于分数阶多孔介质方程的真实解。
• 数值实验表明收敛率为二阶：当 m=1 时，Δx 从 1 降至 0.25（在 200 个节点网格上），误差从 0.011 降至 0.00043。
• 当 m=2 时，误差从 0.0277 降至 0.00085，类似网格细化下仍确认收敛，尽管衰减速率较慢。
• 该方法准确捕捉了解的重尾行为，尤其在 m 较小时表现优异，并成功再现了 m=1 和 m=10 时的 Barenblatt 型解。
• 该格式成功扩展至 φ(u) = log(u+1)，证明了其在一般单调非线性项下的鲁棒性。