变系数时间-空间分数阶扩散方程的高效预处理子

《Mathematics》：An Efficient Preconditioner for Time–Space Fractional Diffusion Equations with Variable Coefficients

　　

本文针对变系数时间-空间分数阶扩散方程提出了一类新型高效预处理子。Caputo时间分数阶导数采用L2-1σ格式离散，Riesz空间分数阶导数则通过移位Grünwald-Letnikov公式近似，由此得到一个适合并行求解的全耦合（all-at-once）系统。由于所得离散系统的系数矩阵??不再具有严格的Toeplitz结构，求解该系统将带来巨大的存储与计算开销。为解决这一困难，研究人员构造了一类Strang-τ结构预处理子：首先利用扩散系数的算术平均近似变系数对角矩阵，进而分别用Strang循环矩阵和τ矩阵逼近时间和空间的Toeplitz矩阵。该预处理子可通过快速傅里叶变换（FFT）、离散正弦变换（DST）以及Sherman-Morrison秩一更新公式高效实现。进一步地，研究人员推导了谱分布估计，表明预处理后的系数矩阵???-1??的特征值位于以1为中心的圆盘内，其半径由时间低秩校正误差、循环逼近误差、空间τ逼近误差以及变系数均值逼近误差共同控制。数值实验表明，所提出的预处理子能够显著加速GMRES(m)迭代的收敛速度。

论文解读：变系数时间-空间分数阶扩散方程的高效Strang-τ结构预处理子

一、研究背景与问题提出

分数阶扩散方程因其能更精确地刻画反常扩散现象而广泛应用于物理、生物及工程领域。然而，该类方程的数值求解面临两大挑战：一是分数阶导数的非局部性导致离散后形成的代数系统系数矩阵稠密，存储与计算复杂度极高；二是当扩散系数随空间位置变化时，系数矩阵丧失严格的Toeplitz结构，传统的基于快速傅里叶变换（FFT）的预处理技术难以直接适用。现有预处理方法多针对常系数情形设计，对变系数问题效果有限。因此，开发一种既能保持结构化矩阵运算效率又能有效处理变系数影响的预处理子具有重要的理论意义与应用价值。

二、研究内容与结论

研究人员针对变系数时间-空间分数阶扩散方程，提出了一类新颖的Strang-τ结构预处理子。该预处理子通过三步近似策略构建：首先用扩散系数的算术平均值替代变系数对角矩阵，其次用Strang循环矩阵近似时间方向的Toeplitz矩阵，最后用τ矩阵近似空间方向的Toeplitz矩阵。研究人员证明，经过该预处理子作用后的系数矩阵的特征值全部集中在以1为中心的一个圆盘内，且圆盘半径由四类误差——时间低秩校正误差、循环逼近误差、空间τ逼近误差以及变系数均值逼近误差——联合控制。数值结果表明，该预处理子能显著降低GMRES(m)迭代次数，大幅提升求解效率。

三、关键技术方法

研究人员采用的主要关键技术包括：（1）Caputo时间分数阶导数的L2-1_σ离散格式；（2）Riesz空间分数阶导数的移位Grünwald-Letnikov近似；（3）Strang循环矩阵与τ矩阵构造方法，用于逼近非Toeplitz结构中的Toeplitz块；（4）基于FFT、离散正弦变换（DST）及Sherman-Morrison秩一更新公式的快速算法，以实现预处理子的高效应用。

四、研究结果

4.1 离散系统的建立与结构分析

研究人员利用L2-1_σ格式和移位Grünwald-Letnikov公式对原方程进行全离散，得到全耦合代数系统。分析表明，尽管变系数的引入使系数矩阵??失去严格Toeplitz结构，但其仍可分解为常数系数Toeplitz部分与变系数修正项之和，为后续预处理子的构造提供了基础。

4.2 Strang-τ预处理子的构造

针对变系数带来的非Toeplitz特性，研究人员提出用扩散系数的算术平均替代变系数对角矩阵，从而将原矩阵近似为一个块Toeplitz矩阵。在此基础上，分别用Strang循环矩阵和τ矩阵近似时间和空间方向的Toeplitz矩阵，最终形成Strang-τ结构预处理子???。该预处理子可通过FFT和DST快速求逆。

4.3 谱性质分析

研究人员推导了预处理后矩阵???^-1??的谱分布估计，证明其特征值均落在以1为中心的圆盘内，且圆盘半径受控于四项误差之和。该结果从理论上保证了预处理子的有效性，即随着网格细化，特征值聚类程度更高，GMRES收敛速度更快。

4.4 数值实验验证

研究人员通过多组数值算例验证了所提方法的性能。结果表明，与无预处理的情形相比，采用Strang-τ预处理子后，GMRES(m)迭代步数显著减少，且计算时间大幅降低。在不同分数阶参数和变系数函数下，该方法均表现出良好的鲁棒性和加速效果。

五、讨论与结论

在讨论部分，研究人员指出，所提出的Strang-τ预处理子框架具有良好的扩展性，可推广至其他类型的分数阶偏微分方程及高维问题。此外，预处理子的构造不依赖于特定边界条件，具有一定的通用性。研究结论如下：本文针对变系数时间-空间分数阶扩散方程，提出了一类高效的Strang-τ结构预处理子。通过算术平均近似变系数、Strang循环矩阵近似时间Toeplitz矩阵、τ矩阵近似空间Toeplitz矩阵，并结合FFT、DST及Sherman-Morrison公式实现了预处理子的快速应用。理论分析和数值实验共同证实，该预处理子能使预处理后矩阵的特征值高度聚集，从而显著加速GMRES(m)迭代收敛，为变系数分数阶扩散方程的高效数值求解提供了一种有力工具。