非均匀介质中时间分数阶反常输运的快速C-Bézier(C-Bézier)有限元法(FEM)的最优误差估计

《Axioms》：Optimal Error Estimates of a Fast C-Bézier Finite Element Method for Time-Fractional Anomalous Transport in Heterogeneous Media

　　

摘要：时间分数阶扩散方程（Time-Fractional Diffusion Equations, TFDEs）是模拟非均匀介质中非经典反常输运的重要模型，但高精度长时间模拟面临两大瓶颈：分数阶导数非局部性导致的O(N2)时间计算复杂度，以及基于多项式的有限元法（Finite Element Method, FEM）在刻画振荡羽流或奇异源时的空间截断误差。本文提出一种将C-Bézier FEM用于空间离散、快速L1格式用于时间离散的框架。通过将C-Bézier基函数的形状参数μ与网格尺寸h耦合（μ=πh），该格式可精确再现对应频率的三角函数型解剖面；当解的空间部分属于C-Bézier空间时，该方法消除空间截断误差并使相关误差常数趋于近零。采用指数和（Sum-of-Exponentials, SOE）逼近将时间方向计算复杂度由O(N2)降至O(N)、存储量降至O(1)，可实现可扩展的三维模拟。研究人员证明了格式具有最优收敛阶O(τ2?α+hk+1)，数值实验验证了该收敛速率。对于能被基函数匹配的解剖面，该方法较Lagrange FEM误差显著减小；对于不属于C-Bézier空间的普通解，两方法具有相同阶数与相近误差幅值，故精度提升仅限于基函数可再现的场。研究人员进一步考察了低正则性情形，包括间断界面与狄拉克δ源注入。

论文解读：非均匀介质中时间分数阶反常输运的快速C-Bézier有限元法的最优误差估计

---

一、研究背景与意义

时间分数阶扩散方程（Time-Fractional Diffusion Equation, TFDE）被广泛用于描述多孔介质、裂隙岩体等非均匀介质中的反常亚扩散输运过程，其时间分数阶导数能刻画记忆效应与长时关联。然而，开展高保真长时间数值模拟存在两个主要障碍：（1）Caputo分数阶时间导数具有全域非局部性，经典L1离散在N个时间层上需O(N²)运算量与O(N)存储，难以推广至三维大规模计算；（2）传统基于多项式形函数的Lagrange有限元法（FEM）对含有振荡羽流或奇异源的TFDE解逼近能力不足，空间截断误差受限于多项式逼近能力。若待求场具三角函数或指数—三角函数特征（如某些本征模态或周期强迫下的稳态分量），多项式基无法精确匹配，引入不必要的离散误差。因此，设计一种可降低时间复杂度且能在特定函数类消除空间离散误差的数值格式具有重要理论与应用价值。该文发表于《Axioms》。

研究人员提出将C-Bézier曲线基函数引入有限元空间离散，并与快速L1时间离散及指数和（Sum-of-Exponentials, SOE）加速相结合，构建时间—空间耦合的数值格式；通过理论分析给出最优误差估计，并通过数值算例验证收敛率及在低正则性情形下的表现。

---

二、主要关键技术方法

研究人员采用如下关键方法开展研究：（1）构造以形状参数μ控制的C-Bézier基函数作为一维/二维有限元试探—检验函数，令μ=πh使之匹配解的三角函数分量；（2）对Caputo时间分数阶导数用L1公式离散，并引入SOE逼近历史核函数将卷积递归压缩为O(N)复杂度与O(1)存储；（3）在C-Bézier有限元空间建立全离散变分格式，利用能量范数技巧与Gronwall不等式推导L²与H¹范数下的最优误差估计O(τ^2?α+h^k+1)（α为分数阶阶数，τ为时间步长，h为网格尺寸，k为C-Bézier空间插值阶数）；（4）设计光滑三角函数解、一般多项式解、带间断系数界面及含Dirac-δ点源四类数值算例，对比C-Bézier FEM与Lagrange FEM的误差与收敛阶。

---

三、研究结果

（1）格式构造与C-Bézier基函数匹配性质

通过将C-Bézier基的形状参数取为μ=πh，证明当真解空间部分属C-Bézier函数空间时可被基函数精确表示，从而消除空间离散引入的截断误差，空间误差常数趋于近似零，此性质为常规多项式Lagrange元所不具备。

（2）快速L1–SOE时间离散与复杂度分析

应用SOE将Caputo导数的核( t?α )展开为若干指数项之和，把历史项卷积转化为递推更新，理论证明时间方向计算复杂度由经典O(N²)降为O(N)，存储需求由O(N)降为O(1)，适用于三维延展模拟。

（3）最优误差估计证明

在齐次Dirichlet边界条件下，基于变分形式、椭圆投影及分数阶离散余项估计，严格导出全离散解在L²(Ω)与H¹₀(Ω)范数下满足误差界O(τ^2?α+h^k+1)，即时间方向达2?α阶、空间方向达拟最优k+1阶收敛，且该阶数与网格无关。

（4）数值实验验证

—光滑三角函数解算例：C-Bézier FEM误差远小于同阶Lagrange FEM，且随h减小误差急剧下降接近机器精度，验证"基函数匹配→误差消除"结论；

—一般光滑解算例：两者收敛阶一致、误差幅值相当，表明优势仅在解可被C-Bézier空间表示时显现；

—间断系数/界面算例与Dirac-δ点源算例：格式保持稳定，收敛阶略有