[论文解读] A unified formulation of splitting-based implicit time integration schemes

S. González‐Pinto, D. Hernández‐Abreu+4 more|arXiv (Cornell University)|Mar 1, 2021

Numerical methods for differential equations参考文献 71被引用 7

一句话总结

本文提出了一种基于隐式-隐式广义结构加法龙格-库塔（IMIM-GARK）方法的统一框架，用于分裂型隐式时间积分格式。该框架将经典方法如交替方向隐式法（ADI）、算子分裂法和分数步长格式统一起来，推导了其阶次条件，并开发了具有优化稳定性与误差特性的新型高阶（3阶与4阶）IMIM-GARK格式，通过抛物型PDE的数值实验验证了其最优收敛性与计算效率。

ABSTRACT

Splitting-based time integration approaches such as fractional steps, alternating direction implicit, operator splitting, and locally one-dimensional methods partition the system of interest into components and solve individual components implicitly in a cost-effective way. This work proposes a unified formulation of splitting time integration schemes in the framework of general-structure additive Runge-Kutta (GARK) methods. Specifically, we develop implicit-implicit (IMIM) GARK schemes, provide the order conditions and stability analysis for this class, and explain their application to partitioned systems of ordinary differential equations. We show that classical splitting methods belong to the IMIM GARK family, and therefore can be studied in this unified framework. New IMIM-GARK splitting methods are developed and tested using parabolic systems.

研究动机与目标

• 将多种基于分裂的隐式时间积分方法——如分数步长法、交替方向隐式法（ADI）、算子分裂法和LOD法——统一于单一数学框架之下。
• 为适用于分区ODE系统的一类隐式-隐式GARK（IMIM-GARK）格式建立系统化的阶次条件理论。
• 构建新型高阶（3阶与4阶）IMIM-GARK格式，其稳定性与误差常数经过优化，以提升精度与效率。
• 在具有精确解的抛物型PDE上验证所提方法，展示其最优收敛性与计算性能。

提出的方法

• 在广义结构加法龙格-库塔（GARK）框架内构建基于分裂的时间积分格式，实现各分量的独立隐式积分。
• 将IMIM-GARK格式定义为其中但彻表可通过置换化为下三角形式的GARK格式，从而实现顺序、解耦的隐式阶段求解。
• 推导了IMIM-GARK格式至四阶的阶次条件，支持高阶格式的系统化构造。
• 将该框架应用于经典格式：ADI、算子分裂法与分数步长龙格-库塔法，表明它们均为IMIM-GARK的特例。
• 基于推导的阶次条件，构造了新型3阶与4阶IMIM-GARK格式，其稳定性与误差常数已优化。
• 采用二阶中心有限差分法进行空间半离散化，并以ℓ2范数测量时间误差与精确解的差异。

实验结果

研究问题

• RQ1 能否将经典的基于分裂的隐式时间积分格式（如ADI、算子分裂法与分数步长法）统一于单一数学框架之下？
• RQ2 隐式-隐式GARK（IMIM-GARK）格式的阶次条件是什么？如何利用这些条件构造高阶格式？
• RQ3 如何在基于分裂的IMIM-GARK格式中优化抛物型PDE的稳定性与误差特性？
• RQ4 在刚性不断增强的情况下，新IMIM-GARK格式在实际中在多大程度上实现了其名义阶次的精度？
• RQ5 在基准抛物型问题上，新IMIM-GARK格式与现有GLM-ADI和IMEX-RK方法相比，在效率与精度方面表现如何？

主要发现

• 所提出的IMIM-GARK框架成功地将ADI、算子分裂法与分数步长法等经典分裂格式统一于单一数学表达形式之下。
• 构建了新型3阶与4阶IMIM-GARK格式，其稳定性与误差常数已优化，在光滑问题上实现了预期的经典阶次收敛性。
• 在二维与三维抛物型测试问题中，新IMIM-GARK格式在粗网格上实现了名义阶次收敛，但在更细网格上因刚性与时间依赖边界条件出现观测到的收敛阶下降。
• 并行ADIGARK变体实现了经典阶次收敛，并展现出良好的计算效率，其中4阶格式在误差-计算时间权衡上优于通用的IMEX-RK 4格式。
• 在效率比较中，ADIGARK 3方法的性能接近IMEX-RK 4方法，而ADIGARK 4方法的效率几乎与最优的GLM-ADI 4方法相当。
• 数值结果证实，IMIM-GARK框架能够支持高阶、稳定且高效的分区系统时间积分格式设计，尤其适用于刚性与多维PDE问题。