[论文解读] Avoiding order reduction when integrating diffusion-reaction boundary value problems with exponential splitting methods

I. Alonso-Mallo, B. Cano+1 more|arXiv (Cornell University)|May 4, 2017

Numerical methods for differential equations参考文献 9被引用 3

一句话总结

本论文提出了一种新颖技术，用于在求解具有非齐次 Dirichlet、Neumann 和 Robin 边界条件的扩散-反应两点边值问题时，防止指数分裂方法中的阶数降低。通过使用精确边界校正的指数积分器和基于 Krylov 的矩阵函数，对时间积分过程进行修改，该方法在一般空间离散化和光滑性假设下，恢复了完整的收敛阶——对于 Lie-Trotter 分裂法为一阶，对于 Strang 分裂法为二阶。

ABSTRACT

In this paper, we suggest a technique to avoid order reduction in time when integrating reaction-diffusion boundary value problems under non-homogeneous boundary conditions with exponential splitting methods. More precisely, we consider Lie-Trotter and Strang splitting methods and Dirichlet, Neumann and Robin boundary conditions. Beginning from an abstract framework in Banach spaces, a thorough error analysis after full discretization is performed and some numerical results are shown which corroborate the theoretical results.

研究动机与目标

• 解决在求解具有非齐次边界条件的非线性扩散-反应问题时，指数分裂方法中长期存在的阶数降低问题。
• 将先前针对具有非齐次边界条件的线性问题的研究结果，扩展至非线性反应-扩散系统。
• 在一般空间离散化和最大范数下，提供严格的误差分析，以确保方法的广泛适用性。
• 通过推导 Lie-Trotter 和 Strang 分裂法在边界校正下的精确公式，确保方法的实际可实施性。
• 证明在完成全空间-时间离散化后，局部误差和全局误差的收敛阶均得以保持。

提出的方法

• 在抽象 Banach 空间框架下建立问题模型，引入算子 A（微分算子）和 ∂（边界算子），以适用于一般 PDE。
• 引入一种校正机制，利用预解算子 K(z) 处理指数积分器中的非齐次边界数据。
• 将 Lie-Trotter 和 Strang 分裂法应用于抽象问题 u' = Au + f(t,u)，其中边界条件为 ∂u = g(t)，并推导出修正的时间推进公式。
• 使用 Krylov 类方法高效且稳定地计算矩阵函数 e^{kA} 和 φ 函数 φ₁(kA)、φ₂(kA)。
• 为两种分裂方法推导出精确的更新公式：Lie-Trotter 方法为公式 (75)–(76)，Strang 方法为公式 (95)、(97)、(99)，均包含边界校正。
• 将该方法应用于有限差分空间离散化，利用块对角结构实现计算效率。

实验结果

研究问题

• RQ1 对于具有非齐次 Dirichlet、Neumann 和 Robin 边界条件的非线性扩散-反应问题，是否可以避免指数分裂方法中的阶数降低？
• RQ2 在边界条件非齐次时，对标准 Lie-Trotter 和 Strang 分裂法进行何种修改，才能在时间方向保持完整的收敛阶？
• RQ3 在完成全空间-时间离散化后，误差行为如何？在何种条件下全局误差的收敛阶得以保持？
• RQ4 所提出的算法是否能够通过 Krylov 子空间技术实现高效的矩阵函数计算，从而实现实际应用？
• RQ5 该理论误差分析是否在一般空间离散化下成立，包括有限差分法和配点法？

主要发现

• 所提出的方法成功避免了在具有非齐次边界条件的非线性扩散-反应问题中，Lie-Trotter 和 Strang 指数分裂方法的阶数降低。
• 对于 Lie-Trotter 分裂法，全局误差以一阶收敛，局部误差也以一阶收敛，该结果通过 k = 5×10⁻³、2.5×10⁻³、1.25×10⁻³ 的数值实验得到验证。
• 对于 Strang 分裂法，当满足公式 (91) 中的空间离散化界时，全局误差以二阶收敛，表 8 中的收敛阶分别为 1.9560 和 1.8946，表明其收敛阶接近理论值。
• Strang 分裂法的局部误差表现出二阶收敛性（收敛阶为 1.8169 和 1.8307），表明该方法保持了理论收敛阶。
• 由于空间算子具有块三对角结构，该方法计算高效，可通过 N×N 块实现快速矩阵-向量乘法。
• 表 7 和表 8 中的数值结果证实，对于一维和二维问题，采用二阶有限差分法时，理论收敛阶得以实现。