带上Robin边界条件的多维时间-空间分数阶扩散问题的数值求解方法

《Applied Numerical Mathematics》：Numerical solution of multidimensional time-space fractional diffusion problems with robin boundary conditions

　　

研究人员提出了一种二阶精度的数值格式，用于求解带上Robin边界条件的多维时间-空间分数阶扩散问题。该方法将矩阵转移技术（Matrix Transfer Technique, MTT）用于空间离散化，结合广义指数时间差分（Generalized Exponential Time Differencing, GETD）方法进行时间积分。所提出的框架扩展了现有方法，使其能够高效处理Robin边界条件、非光滑源项以及分数阶偏微分方程（Fractional Partial Differential Equations, FPDEs）的耦合系统。数值实验验证了该格式在一系列测试问题中的准确性、效率和稳健性。

## 研究背景与问题提出

时间-空间分数阶偏微分方程可用于刻画经典整数阶模型无法充分描述的异常扩散过程，这些方程能够有效捕捉众多物理、生物和工程系统中固有的非局部效应和记忆行为。现有研究中，文献[3]提出了一种求解含谱Laplacian算子的时空分数阶偏微分方程的数值格式，该方法采用截断谱展开结合广义指数时间差分（GETD）方法，但仅适用于Laplacian特征值和特征函数已知的情形，且涉及源项Fourier系数的计算。Ili?等人引入的矩阵转移技术（Matrix Transfer Technique, MTT）为谱算子的离散化提供了高效框架，该技术将非局部空间算子转化为矩阵函数，从而将问题转化为常微分方程组。后续Simpson等人（2008）证明了该逼近精度与底层Laplacian算子的离散化阶数直接相关。此后MTT被广泛应用于异常扩散、分数阶反应扩散系统及复合介质等众多领域。

为应对空间离散化导致的刚性问题，GETD格式提供了无需严格时间步长限制的显式时间积分方法，这些格式利用解的积分形式，精确积分齐次项，同时使用适当的插值或数值积分技术逼近非齐次项。通过全局Padé逼近可进一步提升计算效率，其中部分分式分解能够实现矩阵Mittag-Leffler函数（MLF）的高效求值。

在先前工作[12]中，研究人员针对具有光滑初始条件和源项的时空分数阶扩散模型提出了数值格式，这些模型由空间中的谱分数阶Laplacian算子和时间中的Caputo导数（Caputo derivative）控制。然而，实际应用频繁涉及Robin边界条件和非光滑源项，这些可能引入解的奇异性或不规则性；此外，许多实际问题还涉及多维区域或分数阶偏微分方程耦合系统带来的额外复杂性。尽管部分研究已处理了Robin边界条件的问题，但现有文献中尚未见将该分析扩展至多维或耦合分数阶偏微分方程的工作。为此，研究人员开发了针对Robin型多维时空分数阶扩散问题的二阶数值格式，以填补这一研究空白。

## 主要技术方法

研究所采用的核心技术方法包括以下四个方面：第一，矩阵转移技术（MTT）用于空间离散化，将分数阶谱Laplacian算子转化为离散Laplacian矩阵的分数幂；第二，广义指数时间差分（GETD）方法用于时间积分，采用积分形式的解并近似非齐次项；第三，利用全局Padé逼近结合部分分式分解实现矩阵Mittag-Leffler函数的高效计算；第四，针对Robin边界条件采用适当的离散化处理。研究样本来源于构造的代表性测试问题，涵盖不同边界条件类型、非光滑源项情形以及偏微分方程耦合系统。

## 研究结果

**Robin型时空分数阶边值问题**：研究人员考虑如下初边值问题，其中Ω?Rⁿ，n=1,2,3，为n个有界开区间的笛卡尔积；系数μ和ν为给定常数，满足ν≠0且μ/ν>0；n表示边界?Ω上的单位外法向量；参数κ>0为扩散系数，g为给定源项。该问题包含时间Caputo分数阶导数_c?_t^αu(x,t) = ?κ(?Δ)^β/2u(x,t) + g(x,t)，初始条件u(x,0)=u₀(x)，以及Robin边界条件μu(x,t)+ν?u(x,t)·n=0。研究人员指出Dirichlet和Neumann问题可作为特例处理。

**数值算法**：研究人员概述了求解该初边值问题的数值方法发展，首先采用MTT进行空间离散化，随后采用广义时间差分方法进行时间离散化。通过该两阶段离散过程，将连续问题转化为可数值求解的离散系统。

**算法实现**：研究人员分析了格式实现中的计算挑战，关键在于离散化矩阵A的分数幂计算，以及矩阵自变量Mittag-Leffler函数的高效求值，后者构成了另一主要复杂性来源。

**稳定性与收敛性分析**：研究人员分析了MTT-GETD格式求解问题时的稳定性和收敛速率。关于稳定性，在假设g(·,t)于[0,T]上连续、离散化矩阵A具有正特征值的条件下，对于具有次乘性的任意矩阵范数，均存在与时间步长无关的正常数M，使得数值解满足相应的稳定性估计。该稳定性结果为格式的可靠性提供了理论保障。

**应用与数值实验**：研究人员通过一组代表性测试问题验证算法的准确性和效率，测试涵盖不同类型边界条件（如Robin条件），并包含非光滑源项和偏微分方程组系统的实例。对于具有已知精确解的测试案例，通过公式e^(N):=max_1≤n≤N‖ū(t_n)?Uⁿ‖_∞计算最大绝对误差，其中ū(t_n)表示在时刻t_n按一致顺序排列的精确解向量。数值结果证实了格式的二阶精度、稳健稳定性，并展示了分级网格的有效性。

## 讨论与结论

研究人员在结论部分指出，该研究开发、分析并实现了一种求解Robin型多维时空分数阶扩散边值问题的二阶数值格式。所提出的方法将矩阵转移技术（MTT）用于空间离散化，结合广义指数时间差分（GETD）进行时间积分。数值实验确认了格式的二阶精度和稳健稳定性，并验证了分级网格（graded mesh）的有效性。该格式成功扩展至非光滑源项和耦合分数阶偏微分方程系统，为实际应用中的复杂问题提供了可靠的计算工具。