[论文解读] $H^1$-norm stability and convergence of an L2-type method on nonuniform meshes for subdiffusion equation

Chaoyu Quan, Xu Wu|arXiv (Cornell University)|May 12, 2022

Differential Equations and Numerical Methods被引用 4

一句话总结

本论文在一般非均匀网格上建立了 L2 类方法求解亚扩散方程的 H¹-范数稳定性与收敛性。在时间步长比约束较宽松（0.457 ≤ ρₖ ≤ 3.562）的条件下，证明了关键双线性形式的半正定性，从而实现了长时间 H¹-稳定性，并在分级参数满足 r > 5/α − 1 的改进型分级网格上，获得了 H¹-范数下 (5−α)/2 阶的最优收敛率。

ABSTRACT

This work establishes $H^1$-norm stability and convergence for an L2 method on general nonuniform meshes when applied to the subdiffusion equation. Under mild constraints on the time step ratio $ρ_k$, such as $0.4573328\leq ρ_k\leq 3.5615528$ for $k\geq 2$, the positive semidefiniteness of a crucial bilinear form associated with the L2 fractional-derivative operator is proved. This result enables us to derive long time $H^1$-stability of L2 schemes. These positive semidefiniteness and $H^1$-stability properties hold for standard graded meshes with grading parameter $15/α-1$. To the best of our knowledge, this study is the first work on $H^1$-norm stability and convergence of L2 methods on general nonuniform meshes for the subdiffusion equation.

研究动机与目标

• 在一般非均匀时间网格上，建立 L2 类格式求解亚扩散方程的 H¹-范数稳定性与收敛性。
• 解决 L2 方法在非均匀网格上缺乏 H¹-稳定性分析的问题，特别是针对低正则性解的情形。
• 在较宽松的时间步长比约束下，证明与 L2 分数阶导数算子相关的双线性形式的半正定性。
• 将稳定性与收敛性结果推广至具有特定分级参数的分级网格与改进型分级网格。
• 提供严格的 H¹-范数误差分析，并通过数值实验验证收敛速率。

提出的方法

• 在时间步长比满足 $\rho_k \in [0.4573328, 3.5615528]$（$k \geq 2$）的条件下，证明双线性形式 $\mathcal{B}_n(v,w) = \sum_{k=1}^n \langle L_k^\alpha v, \delta_k w \rangle$ 的半正定性。
• 采用基于离散分数阶 Grönwall 不等式与能量方法的新型分析框架，推导长时间 H¹-稳定性。
• 将该分析方法应用于标准分级网格（分级参数 $r \leq 3.2016538$），建立 H¹-稳定性。
• 推导一般非均匀网格下的 H¹-范数误差估计，并证明当 $r > 5/\alpha - 1$ 时，改进型分级网格下收敛阶为 $(5 - \alpha)/2$。
• 在空间方向采用谱配置法，使用切比雪夫–高斯–勒让德配点以实现高阶空间精度。
• 通过在改进型分级网格上进行数值实验，验证理论结果，测量最大 H¹-误差，参数包括 $\alpha$、$r$ 和 $N$。

实验结果

研究问题

• RQ1 在一般非均匀网格上，当时间步长比约束较宽松时，与 L2 分数阶导数算子相关的双线性形式是否为半正定？
• RQ2 能否在非均匀网格上为亚扩散方程的隐式 L2 格式建立长时间 H¹-稳定性？
• RQ3 在改进型分级网格上，L2 方法在 H¹-范数下的最优收敛率是多少？其与分级参数 $r$ 的关系如何？
• RQ4 在数值测试中，双线性形式的半正定性结果是否仍保持紧致性？
• RQ5 理论 H¹-范数收敛率 $(5 - \alpha)/2$ 是否能在实际中观测到，且与观测到的收敛阶是否一致？

主要发现

• 在较宽松约束 $0.4573328 \leq \rho_k \leq 3.5615528$（$k \geq 2$）下，双线性形式 $\mathcal{B}_n$ 被证明为半正定，从而支持 H¹-稳定性分析。
• 在一般非均匀网格上，对具有齐次 Dirichlet 边界条件的亚扩散方程，L2 格式的长时间 H¹-稳定性得以建立。
• 对于标准分级网格，当分级参数满足 $1 < r \leq 3.2016538$ 时，H¹-稳定性成立。
• 当 $r > 5/\alpha - 1$ 时，改进型分级网格下 H¹-范数收敛阶为 $(5 - \alpha)/2$。
• 数值结果表明，H¹-范数下的实际收敛阶约为 $\min\{r\alpha, 3 - \alpha\}$，与理论预期一致。
• 数值验证了最小特征值下界 $\lambda_{\min} \geq C N^{\alpha}$ 的最优性，且最大特征值满足 $\lambda_{\max} \geq C N^{r\alpha}$。