[论文解读] Discrete maximal regularity of time-stepping schemes for fractional evolution equations

Bangti Jin, Buyang Li+1 more|arXiv (Cornell University)|Jun 24, 2016

Differential Equations and Numerical Methods参考文献 25被引用 4

一句话总结

本文为各种时间步进格式——卷积求积法、L1格式、显式欧拉法以及分数阶Crank-Nicolson格式——在Caputo导数阶数 $\alpha \in (0,2)\setminus\{1\}$ 的分数阶演化方程上的应用，建立了离散的最大 $\varepsilon^{p}$-正则性。核心贡献在于证明了这些格式在 $\varepsilon^{p}$-范数下继承了最大正则性估计，通过算子值傅里叶乘子理论及其离散类比，将经典抛物型结果推广至分数阶时间模型。

ABSTRACT

In this work, we establish the maximal $\ell^p$-regularity for several time stepping schemes for a fractional evolution model, which involves a fractional derivative of order $α\in(0,2)$, $α eq 1$, in time. These schemes include convolution quadratures generated by backward Euler method and second-order backward difference formula, the L1 scheme, explicit Euler method and a fractional variant of the Crank-Nicolson method. The main tools for the analysis include operator-valued Fourier multiplier theorem due to Weis [48] and its discrete analogue due to Blunck [10]. These results generalize the corresponding results for parabolic problems.

研究动机与目标

• 将经典抛物方程的最大 $L^{p}$-正则性理论推广至Caputo导数阶数 $\alpha \in (0,2)\setminus\{1\}$ 的分数阶演化方程。
• 分析时间步进格式的离散最大正则性，包括卷积求积法、L1格式、显式欧拉法以及一种分数阶Crank-Nicolson方法。
• 为完全离散解建立与时间步长 $\tau$ 无关的 $\ell^{p}$-范数稳定性与正则性估计。
• 将Weis与Blunck的乘子定理推广至离散与分数阶设置，以支持在分数阶时间离散化背景下对算子值傅里叶乘子的分析。
• 通过确保在高阶范数下的强收敛性，为非线性分数阶PDE的数值分析提供理论基础。

提出的方法

• 利用Weis [48] 的算子值傅里叶乘子定理及其离散类比（Blunck [10]）来刻画离散设置下的最大正则性。
• 通过其生成函数和拉普拉斯变换与 $z$-变换域中的符号分析，研究时间步进格式，以验证相关算子族的 $R$-有界性。
• 应用由向后欧拉法与二阶后向差分公式生成的卷积求积理论，推导稳定性和正则性估计。
• 采用L1格式与Crank-Nicolson方法的分数阶变体对Caputo导数 $\partial_t^\alpha u = Au + f$ 进行离散化，确保在 $\ell^p$-空间中的一致性与稳定性。
• 建立形式为 $\|\tau^{-1}(u^n - u^{n-1})\|_{\ell^p(X)} + \|Au^n\|_{\ell^p(X)} \leq c_{p,X} \|f^n\|_{\ell^p(X)}$ 的离散最大正则性估计，且对 $\tau$ 一致成立。
• 证明离散估计在满足UMD性质的Banach空间 $X$ 中成立，并通过Laplacian的 $R$-扇形性将结果推广至 $L^q$-空间。

实验结果

研究问题

• RQ1 对于 $\alpha \in (0,2)\setminus\{1\}$ 的分数阶演化方程，时间步进格式是否在 $\ell^p$-范数下继承离散最大正则性？
• RQ2 Weis与Blunck的乘子定理能否被适配至离散与分数阶设置，以建立时间离散算子的 $R$-有界性？
• RQ3 卷积求积格式（如向后欧拉法与BDF2）在分数阶亚扩散与扩散波方程中的稳定性与正则性行为如何？
• RQ4 L1格式与分数阶Crank-Nicolson方法在 $\alpha \in (0,1)$ 与 $\alpha \in (1,2)$ 下的离散最大正则性表现如何？
• RQ5 离散估计在多大程度上将经典抛物问题的结果推广至分数阶时间模型？

主要发现

• 为由向后欧拉法与二阶后向差分公式生成的卷积求积法建立了离散最大 $\ell^p$-正则性，且常数与时间步长 $\tau$ 无关。
• L1格式与显式欧拉法在 $\alpha \in (0,1) \cup (1,2)$ 下也满足 $\ell^p$-范数下的离散最大正则性估计，将稳定性结果推广至分数阶模型。
• 证明了分数阶Crank-Nicolson方法的变体可保持离散最大正则性，且稳定性常数对 $\tau$ 一致。
• 分析表明，离散估计推广了经典抛物结果，分数阶方程的 $\ell^p$-正则性界与 $\alpha=1$ 情况下的经典结果相同。
• 对于有界区域 $\Omega$ 上的Dirichlet Laplacian，$W^{2,q}$-范数下的离散误差满足 $\|(u^n - u(t_n))\|_{L^p(W^{2,q}(Ω))} \leq c_{p,q} \tau$，表明在强范数下具有最优收敛性。
• 建立了Laplacian在 $L^q$-实现中的 $R$-扇形性（$1 < q < \infty$），使得乘子理论可应用于 $L^q$-空间，从而支持主要正则性结果。