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

Félix del Teso, Juan Luís Vázquez|arXiv (Cornell University)|Jul 9, 2013

Differential Equations and Numerical Methods参考文献 14被引用 21

一句话总结

本文提出了一种用于广义分数阶多孔介质方程（$ m \geq 1 $ 且 $ \sigma \in (0,2) $）的有限差分格式，利用 Caffarelli-Silvestre 扩展方法处理非局部分数阶拉普拉斯算子。该方法实现了 $ O(\Delta t + \Delta x^{2-\sigma} + \frac{\Delta x^a}{\Delta t}) $ 的收敛阶，其中 $ a $ 依赖于 $ \sigma $，并严格证明了数值解的存在性、唯一性及对理论解的收敛性。

ABSTRACT

We formulate a numerical method to solve the porous medium type equation with fractional diffusion \[ \frac{\partial u}{\partial t}+(-Δ)^{σ/2} (u^m)=0 \] posed for $x\in \mathbb{R}^N$, $t>0$, with $m\geq 1$, $σ\in (0,2)$, and nonnegative initial data $u(x,0)$. We prove existence and uniqueness of the solution of the numerical method and also the convergence to the theoretical solution of the equation with an order depending on $σ$. We also propose a two points approximation to a $σ$-derivative with order $O(h^{2-σ})$.

研究动机与目标

• 为具有非局部扩散的广义分数阶多孔介质方程开发一种鲁棒的数值格式。
• 通过扩展方法解决 $ \sigma \neq 1 $ 时分数阶拉普拉斯算子的实现挑战。
• 证明数值解对理论解的存在性、唯一性及收敛性。
• 建立依赖于分数阶参数 $ \sigma $ 的收敛速率。
• 确保在有界区域上的数值解随区域扩大而收敛至 $ \mathbb{R}^N $ 上的解。

提出的方法

• 通过 Caffarelli-Silvestre 扩展将分数阶扩散项 $ (-\Delta)^{\sigma/2}(u^m) $ 重表述，将非局部问题转化为高一维空间中的局部 PDE。
• 在有界计算区域 $ \Omega = [-X,X] \times [0,X] $ 上，通过有限差分格式求解扩展后的方程，并施加适当的边界条件。
• 提出一种二点有限差分逼近，阶数为 $ O(h^{2-\sigma}) $，用于处理扩展形式中出现的加权导数 $ \partial / \partial y^\sigma $。
• 证明该数值格式具有稳定性与收敛性，局部截断误差被控制在 $ O(\Delta t + \Delta x^{2-\sigma}) $ 以内。
• 在 $ \Omega \to \mathbb{R}^N $ 且 $ \Delta x \to 0 $ 的条件下，建立了有界区域上数值解对全空间解的收敛性，边界误差通过 Barenblatt 解估计进行控制。
• 在假设 $ \Delta t \leq C(m,f) \Delta x^\sigma $ 的条件下，证明了理论收敛性，确保了格式的稳定性和相容性。

实验结果

研究问题

• RQ1 能否构造一种有限差分方法，以准确且稳定的方式逼近广义 $ \sigma \in (0,2) $ 和 $ m \geq 1 $ 条件下分数阶多孔介质方程的解？
• RQ2 如何在数值格式中高效且准确地逼近非局部分数阶拉普拉斯算子？
• RQ3 所提出的数值格式的收敛速率是多少，其如何依赖于分数阶参数 $ \sigma $？
• RQ4 在有界区域上的数值解如何收敛至整个空间 $ \mathbb{R}^N $ 上的解？
• RQ5 为确保数值解收敛至理论解，对区域大小和时间步长有何要求？

主要发现

• 所提出的有限差分格式被证明具有数值解的存在性与唯一性。
• 该方法实现了对理论解的收敛，收敛阶为 $ O(\Delta t + \Delta x^{2-\sigma} + \frac{\Delta x^a}{\Delta t}) $，其中 $ a $ 依赖于 $ \sigma $：当 $ \sigma \in (0,1] $ 时 $ a=2 $，当 $ \sigma \in (1,3/2] $ 时 $ a=3 $，当 $ \sigma \in (3/2,2) $ 时 $ a=4 $。
• 对 $ \sigma $-导数的二点逼近实现了 $ O(h^{2-\sigma}) $ 阶的局部截断误差，该阶数对给定的离散化是最优的。
• 在满足 $ |X| \geq L \cdot \Delta x^{-a/(N+\sigma)} $ 的条件下，证明了有界区域上数值解对 $ \mathbb{R}^N $ 上解的收敛性，确保边界误差以 $ \Delta x^a $ 速度衰减。
• 在稳定性条件 $ \Delta t \leq C(m,f) \Delta x^\sigma $ 下，收敛速率达到最优，该条件保证了格式的一致性与稳定性。
• 理论收敛性通过数值实验得到验证，总误差被控制在时间、空间与边界截断误差之和范围内，各项误差均通过 Barenblatt 解估计实现控制。