[论文解读] A second-order scheme with nonuniform time steps for a linear reaction-sudiffusion problem

Hong-lin Liao, William McLean+1 more|arXiv (Cornell University)|Jan 1, 2018

Fractional Differential Equations Solutions参考文献 16被引用 31

一句话总结

本文提出了一种误差卷积结构（ECS）分析方法，针对基于Alikhanov逼近的Caputo分数阶导数的二阶精度、非均匀时间步长格式。通过在拟均匀非均匀时间网格上采用分段二次和线性插值，作者建立了线性反应-亚扩散问题的L²-范数误差估计，其阶为O(τ² + h²)，证明了即使在初始奇异性存在的情况下仍具有稳定性和收敛性，并表明拟均匀网格能有效解析解的时间导数中的t^{α−1}奇异性。

ABSTRACT

Stability and convergence of a time-weighted discrete scheme with nonuniform time steps are established for linear reaction-subdiffusion equations. The Caupto derivative is approximated at an offset point by using linear and quadratic polynomial interpolation. Our analysis relies on two tools: a discrete fractional Grönwall inequality and the global consistency analysis. The new consistency analysis makes use of an interpolation error formula for quadratic polynomials, which leads to a convolution-type bound for the local truncation error. To exploit these two tools, some theoretical properties of the discrete kernels in the numerical Caputo formula are crucial and we investigate them intensively in the nonuniform setting. Taking the initial singularity of the solution into account, we obtain a sharp error estimate on nonuniform time meshes. The fully discrete scheme generates a second-order accurate solution on the graded mesh provided a proper grading parameter is employed. An example is presented to show the sharpness of our analysis.

研究动机与目标

• 为线性反应-亚扩散问题开发一种二阶时间离散化格式，以在非均匀时间网格上实现更高的精度。
• 解决解的时间导数中初始奇异性的问题，该问题限制了在均匀网格上的收敛性。
• 为具有可变时间步长的一般非均匀网格上的分数阶Crank–Nicolson方法建立严格的稳定性和收敛性理论。
• 证明拟均匀网格（如tk = T(k/N)^γ）即使在解在t=0处正则性较低时，也能恢复二阶收敛性。
• 将ECS分析框架扩展至更高阶逼近，特别是Alikhanov公式，以实现误差定位和全局一致性分析。

提出的方法

• 提出误差卷积结构（ECS）分析方法，以局部化并简化非均匀网格上Caputo导数离散卷积逼近的误差分析。
• 在子区间[tk−1, tk]（k = 1 到 n−1）使用分段二次插值，在[tn−1, tn−θ]使用线性插值，构造具有偏移量θ = α/2的Alikhanov逼近。
• 推导出适用于二次多项式的新型插值误差公式，其余项为积分形式，用于在ECS框架内控制局部截断误差。
• 建立全局一致性误差界，揭示长时间积分下误差的分布行为，表明当网格分级遵循误差等分布原则时，误差会变得有效局部化。
• 将ECS框架应用于证明分数阶Crank–Nicolson格式在局部时间步长比ρ ≤ 7/4的温和限制下的稳定性和收敛性。
• 在空间上采用Galerkin有限元方法，在时间上采用类似Crank–Nicolson的步长格式，实现时间方向二阶精度和空间方向二阶精度。

实验结果

研究问题

• RQ1 能否为线性反应-亚扩散问题构造一种二阶时间离散化格式，使其在具有可变时间步长的一般非均匀时间网格上保持稳定和收敛？
• RQ2 如何解决解中初始奇异性（∂tu ∼ t^{α−1}）的问题，以恢复时间方向的二阶收敛性？
• RQ3 网格分级（如tk = T(k/N)^γ）在实现非均匀时间步长格式的最优收敛性中起什么作用？
• RQ4 误差卷积结构（ECS）框架能否有效应用于Alikhanov公式等高阶逼近，以简化并实现误差分析的局部化？
• RQ5 非均匀Alikhanov格式的局部时间步长比ρ = max(τk/τk+1)的最优稳定性条件是什么？

主要发现

• ECS分析成功地将非均匀Alikhanov逼近的误差分析局部化，简化了对一般非均匀时间网格的处理。
• 为分数阶Crank–Nicolson格式导出了精确的L²-范数误差估计，阶为O(τ² + h²)，证实了时间方向的二阶收敛性。
• 在局部时间步长比满足ρ ≤ 7/4的条件下，证明了格式的稳定性和收敛性，允许自适应时间推进，且τ_next / τ_current ≥ 4/7。
• 当γ ≥ 1时，拟均匀网格能有效解析初始奇异性，即使在t=0处解不光滑，也能恢复二阶收敛性。
• 全局一致性误差界表明，当网格按误差等分布原则分级时，分数阶导数的非局部性会导致误差分布趋于局部化。
• 通过一个数值例子验证了理论误差界的有效性，表明其具有尖锐性。