[论文解读] Error Estimates for Approximations of Distributed Order Time Fractional Diffusion with Nonsmooth Data

Bangti Jin, Raytcho Lazarov+2 more|arXiv (Cornell University)|Apr 7, 2015

Fractional Differential Equations Solutions参考文献 19被引用 4

一句话总结

本文针对具有非光滑初值的分布式阶时间分数阶扩散方程，进行了严格的数值分析，提出了一种空间半隐式伽辽金有限元方法及两种全隐式格式——基于拉普拉斯变换求积和基于向后欧拉时间离散化的卷积求积。在 $L^2(\nabla)$ 和 $H^1$ 范数下建立了最优误差估计，其中基于拉普拉斯变换的格式呈现指数收敛，而卷积求积格式呈现一阶收敛，且通过大量数值实验得到验证。

ABSTRACT

In this work, we consider the numerical solution of an initial boundary value problem for the distributed order time fractional diffusion equation. The model arises in the mathematical modeling of ultra-slow diffusion processes observed in some physical problems, whose solution decays only logarithmically as the time $t$ tends to infinity. We develop a space semidiscrete scheme based on the standard Galerkin finite element method, and establish error estimates optimal with respect to data regularity in $L^2(D)$ and $H^1(D)$ norms for both smooth and nonsmooth initial data. Further, we propose two fully discrete schemes, based on the Laplace transform and convolution quadrature generated by the backward Euler method, respectively, and provide optimal convergence rates in the $L^2(D)$ norm, which exhibits exponential convergence and first-order convergence in time, respectively. Extensive numerical experiments are provided to verify the error estimates for both smooth and nonsmooth initial data, and to examine the asymptotic behavior of the solution.

研究动机与目标

• 开发并分析适用于具有非光滑初值的分布式阶时间分数阶扩散方程的数值格式，此类初值在超慢扩散过程中常见。
• 针对光滑与非光滑初值，在 $L^2(\nabla)$ 和 $H^1$ 范数下建立最优误差估计，弥补文献中的空白。
• 提出并分析两种全隐式格式：一种基于拉普拉斯变换求积，实现指数收敛；另一种基于卷积求积与向后欧拉时间离散化，实现一阶收敛。
• 通过大量数值实验验证理论误差估计，涵盖光滑与非光滑初值。
• 研究解的渐近行为，确认理论误差界（包括对数项）的紧致性。

提出的方法

• 采用标准 $H^1_0$-协调有限元，构建空间半隐式伽辽金有限元方法，误差分析依赖于福岛与铃木提出的一种算子技巧。
• 提出两种全隐式格式：一种基于变形双曲路径上的逆拉普拉斯变换求积；另一种基于向后欧拉方法生成的卷积求积。
• 基于拉普拉斯变换的格式利用路径积分与有理逼近，实现时间方向的指数收敛。
• 卷积求积格式利用向后欧拉方法生成稳定的一阶收敛时间离散化。
• 误差估计基于非光滑初值下解的正则性结果，其依赖于终时间 $T$，且在 $L^2$ 初值情形下包含对数项。
• 在 $\mathbb{R}^d$（$d=1,2,3$）的凸多边形区域上进行数值实验，采用不同初值与权函数 $\mu(\alpha)$，包括 $\chi_{[1/2,1]}(\alpha)$。

实验结果

研究问题

• RQ1 空间半隐式伽辽金有限元方法在具有非光滑初值的分布式阶时间分数阶扩散问题中，其最优收敛率是多少？
• RQ2 基于拉普拉斯变换求积的全隐式格式能否实现该问题的指数收敛？其与卷积求积格式相比表现如何？
• RQ3 误差估计如何依赖于初值的正则性，特别是当 $v \in L^2(\Omega)$ 时？
• RQ4 误差估计中的对数因子 $\ell_2(t) = \log(1/t)$ 是否紧致？其是否在数值实验中显现？
• RQ5 不同初值类型（光滑与奇异）下收敛行为如何变化？这些格式能否高效处理大时间区间？

主要发现

• 空间半隐式伽辽金有限元方法在 $L^2(\Omega)$ 和 $H^1(\Omega)$ 范数下均实现最优误差估计，收敛率取决于初值正则性。
• 基于拉普拉斯变换的全隐式格式在时间方向呈现指数收敛，因此在计算大量或大时间点的解时极为高效。
• 基于向后欧拉时间离散化的卷积求积格式实现时间方向的一阶收敛，所有测试案例的误差率均通过数值验证。
• 对于 $L^2(\Omega)$ 初值，误差估计包含对数因子 $\ell_2(t) = \log(1/t)$，数值实验确认其紧致性，显示 $\|U_h^1 - u(\tau)\|_{L^2} \leq c\tau \log(1/\tau)$。
• 数值结果表明，向后欧拉格式的误差主要受 $\tau / t$ 主导，且在小 $t$ 时收敛率显著下降，与理论预测一致。
• 大量数值实验验证了光滑与非光滑初值下的理论误差估计，并展示了与超慢扩散模型一致的解的渐近对数衰减行为。