[论文解读] Multi-term fractional linear equations modeling oxygen subdiffusion through capillaries

Vittorino Pata, Sergii V. Siryk+1 more|arXiv (Cornell University)|Oct 10, 2022

Fractional Differential Equations Solutions被引用 4

一句话总结

本文提出了一种新颖的分析与数值框架，用于建模毛细血管中氧气运输的多术语分数阶亚扩散方程，采用Caputo导数和时变系数。通过引入一种专门的正则化技术，该框架在一般条件下建立了全局经典可解性，克服了对非负卷积核的依赖，且通过高精度有限差分格式与数值示例验证了该方法的有效性。

ABSTRACT

For $0<ν_2<ν_1\leq 1$, we analyze a linear integro-differential equation on the space-time cylinder $Ω imes(0,T)$ in the unknown $u=u(x,t)$ $$\mathbf{D}_{t}^{ν_1}(\varrho_{1}u)-\mathbf{D}_{t}^{ν_2}(\varrho_2 u)-\mathcal{L}_{1}u-\mathcal{K}*\mathcal{L}_{2}u =f$$ where $\mathbf{D}_{t}^{ν_i}$ are the Caputo fractional derivatives, $\varrho_i=\varrho_i(x,t)$ with $\varrho_1\geq μ_0>0$, $\mathcal{L}_{i}$ are uniform elliptic operators with time-dependent smooth coefficients, $\mathcal{K}$ is a summable convolution kernel, and $f$ is an external force. Particular cases of this equation are the recently proposed advanced models of oxygen transport through capillaries. Under suitable conditions on the given data, the global classical solvability of the associated initial-boundary value problems is addressed. To this end, a special technique is needed, adapting the concept of a regularizer from the theory of parabolic equations. This allows us to remove the usual assumption about the nonnegativity of the kernel representing fractional derivatives. The problem is also investigated from the numerical point of view.

研究动机与目标

• 解决具有变系数的多术语时间分数阶亚扩散方程初边值问题的适定性问题。
• 克服分数阶导数模型中传统限制——即要求卷积核非负。
• 开发一种高精度与高稳定性的数值格式，用于求解此类方程。
• 为所提模型在真实生物情境下的经典可解性提供理论与计算证据。
• 为未来关于反问题与模型非线性扩展的研究奠定基础。

提出的方法

• 作者采用源自抛物型PDE理论的广义正则化技术，以处理多术语分数阶方程的积分微分结构。
• 模型使用阶数为 $\nu_1$ 和 $\nu_2$ 的Caputo分数阶导数，满足 $0 < \nu_2 < \nu_1 \leq 1$，从而在氧气亚扩散中体现记忆效应。
• 方程包含时变扩散系数 $\varrho_1, \varrho_2$，一致椭圆微分算子 $\mathcal{L}_1, \mathcal{L}_2$，以及可 summable 的卷积核 $\mathcal{K}$，以表征记忆效应。
• 采用Crank-Nicolson型时间离散与空间中心差分相结合的有限差分格式，并利用卷积四分法处理分数阶导数。
• 通过已知解析形式的制造解验证该格式，确保其一致性和收敛性。
• 推导出理论先验估计，以在系数光滑性与椭圆性条件下证明全局经典可解性。

实验结果

研究问题

• RQ1 在何种条件下，多术语分数阶亚扩散方程的初边值问题具有经典可解性？
• RQ2 是否可在不假设分数阶导数项中卷积核非负的条件下，建立经典可解性？
• RQ3 所提出的有限差分格式能以多高精度逼近多术语分数阶亚扩散方程的解？
• RQ4 对于不同分数阶阶数 $\nu_1$ 与 $\nu_2$，该数值格式的收敛行为如何？
• RQ5 该模型能否推广至反问题或具有退化系数的非线性变体？

主要发现

• 所提出的正则化技术成功在不依赖卷积核非负性条件下，建立了初边值问题的全局经典可解性。
• 有限差分格式实现了高精度，其离散 $L^2$-范数误差 $\gimel$ 在 $\nu_1 = 0.9$ 时低至 $1.8727 \times 10^{-4}$。
• 当 $\nu_1 = 0.5$ 时，数值解与制造的解析解高度吻合，$\gimel = 3.1859 \times 10^{-4}$，验证了格式的可靠性。
• 该方法在不同分数阶阶数下均表现稳健，误差从 $\nu_1 = 0.1$ 时的 $6.4793 \times 10^{-4}$ 降低至 $\nu_1 = 0.9$ 时的 $1.8727 \times 10^{-4}$，表明高阶数下收敛性更优。
• 数值结果证实了理论收敛速率，表明该格式具有良好的稳定性和效率。
• 该模型成功捕捉了氧气运输中复杂的亚扩散动力学，且在二维时空域中实现了全系统的精确求解。