基于次序关系的Caputo分数阶传播子近似及其相关数值方法
Subordination based approximation of Caputo fractional propagator and related numerical methods
摘要 Abstract
本文提出了一种针对Caputo分数阶传播子$S_\alpha(t)$以及与时间无关的扇形算子系数$A$的Cauchy问题的温和解的指数收敛数值方法,其中时间变量的Caputo分数阶导数的阶数为$\alpha \in (0,2)$。所提出的数值方法通过推广先前针对$S_\alpha(t)$构建的次序原则近似方法而得来。该技术能够消除误差估计主项对$\alpha$的依赖性,同时保留原近似方法的其他计算相关特性:支持多级并行性、处理初始数据时对空间光滑性的最低要求,以及对于所有$t \in [0, T]$的稳定指数收敛性。最终,次序原则的应用显著改善了方法的收敛行为,特别是在$\alpha < 0.5$的小值情况下,并为进一步提高数据重用效率提供了更多可能性。为了验证理论结果,我们考虑了所开发方法在直接问题(解的逼近)和反问题(分数阶识别)中的应用。
In this work, we propose an exponentially convergent numerical method for the Caputo fractional propagator $S_\alpha(t)$ and the associated mild solution of the Cauchy problem with time-independent sectorial operator coefficient $A$ and Caputo fractional derivative of order $\alpha \in (0,2)$ in time. The proposed methods are constructed by generalizing the earlier developed approximation of $S_\alpha(t)$ with help of the subordination principle. Such technique permits us to eliminate the dependence of the main part of error estimate on $\alpha$, while preserving other computationally relevant properties of the original approximation: native support for multilevel parallelism, the ability to handle initial data with minimal spatial smoothness, and stable exponential convergence for all $t \in [0, T]$. Ultimately, the use of subordination leads to a significant improvement of the method's convergence behavior, particularly for small $\alpha < 0.5$, and opens up further opportunities for efficient data reuse. To validate theoretical results, we consider applications of the developed methods to the direct problem of solution approximation, as well as to the inverse problem of fractional order identification.