[论文解读] Multi-stage waveform Relaxation and Multisplitting Methods for Differential Algebraic Systems

Jürgen Geiser|arXiv (Cornell University)|Jan 4, 2016

Numerical methods for differential equations参考文献 6被引用 4

一句话总结

本文提出了一种用于求解微分代数方程（DAE）的多阶段波形松弛（MSWR）与多分裂方法，结合了分层矩阵分裂与并行迭代求解器。通过将系统分解为外层与内层迭代——最终简化为三阶段求解器——并利用单位分解矩阵实现多分裂，该方法在保持与单阶段方法相当精度的同时显著加速了收敛，通过更小的矩阵求逆和并行化实现了显著的计算节省。

ABSTRACT

We are motivated to solve differential algebraic equations with new multi-stage and multisplitting methods. The multi-stage strategy of the waveform relaxation (WR) methods are given with outer and inner iterations. While the outer iterations decouple the initial value problem of differential algebraic equations (DAEs) in the form of $A \frac{d y(t)}{dt} + B y(t) = f(t)$ to $M_A \frac{d y^{k+1}(t)}{dt} + M_1 y^{k+1}(t) = N_1 y^k(t) + N_A \frac{d y^{k}(t)} + f(t)$, where $A = M_A - N_A$, $B = M_1 - N_1$. The inner iterations decouple further $M_1 = M_2 - N_2$ and $M_2 = M_3 - N_3$ with additional iterative processes, such that we result to invert simpler matrices and accelerate the solver process. The multisplitting method use additional a decomposition of the outer iterative process with parallel algorithms, based on the partition of unity, such that we could improve the solver method. We discuss the different algorithms and present a first experiment based on a DAE system.

研究动机与目标

• 解决传统迭代求解器在求解大规模微分代数方程（DAE）时计算效率低下的问题。
• 通过引入多级分裂策略，降低波形松弛（WR）方法中矩阵求逆的计算成本。
• 通过基于单位分解的多分裂技术，增强并行化能力。
• 研究在不同分裂层次与多分裂方案下，计算效率与精度之间的权衡。
• 在半离散化的PDAE系统上，验证所提方法的可行性和性能。

提出的方法

• 采用三阶段波形松弛（MSWR）策略：外层迭代通过 $ A = M_A - N_A $，$ B = M_1 - N_1 $ 分解系统，随后进行两轮内层迭代，形式为 $ M_1 = M_2 - N_2 $，$ M_2 = M_3 - N_3 $，从而在每一级简化矩阵求逆。
• 在内层阶段使用迭代求解器：$ \frac{dy^{k+1}}{dt} + M_1 y^{k+1} = N_1 y^k + f(t) $，以及后续的 $ z^{\nu+1} $，$ \tilde{z}^{\mu+1} $ 系统，避免对大规模矩阵进行直接求逆。
• 通过单位分解矩阵 $ E_p $ 实现多分裂，其中 $ \sum_p E_p = I $，支持 $ L $ 个子系统的并行求解，并采用加权更新 $ y^{k+1} = \sum_p E_p y^{p,k+1} $。
• 在多分裂中引入重叠域分解，重叠度为 $ o=1 $（最小重叠）和 $ o=m/2-1 $（最大重叠），以研究其对收敛性和精度的影响。
• 采用雅可比与高斯-赛德尔型多分裂变体，对比其收敛行为与并行效率。
• 将方法应用于半离散化的PDAE系统，采用时间窗积分，参数为 $ h=0.1 $，$ J=20 $，并基于误差范数 $ \leq 10^{-3} $ 或固定 $ K=20 $ 次外层迭代作为停止准则。

实验结果

研究问题

• RQ1 与单阶段WR方法相比，多阶段波形松弛是否能在保持精度的同时降低计算成本？
• RQ2 引入内层迭代（如两级或三级分裂）如何影响DAE求解器的收敛性与效率？
• RQ3 在DAE求解器中，采用重叠子域的多分裂在多大程度上提升了并行化与收敛性？
• RQ4 不同类型的多分裂（雅可比与高斯-赛德尔）对方法的精度与收敛速度有何影响？
• RQ5 所提出的多阶段与多分裂方案是否能在显著降低计算量的同时，达到与单阶段方法相当的精度？

主要发现

• 三阶段MSWR方法在保持与单阶段方法相同精度的同时，通过分层矩阵分裂显著降低了计算成本。
• 两阶段与三阶段方法在 $ L_2 $ 与 $ L_\infty $ 范数下均与单阶段方法达到相同的误差水平，如图1所示。
• 采用高斯-赛德尔型更新的多分裂方法在重叠配置下精度高于雅可比型更新。
• 最大重叠（$ o = m/2 - 1 $）配置相比最小重叠（$ o = 1 $）表现出更优的收敛性，但最佳性能介于两者之间。
• 多分裂方法在保持与单阶段及多阶段WR方法相同精度的同时，实现了在多处理器上的完全并行化。
• 计算节省显著：通过在每一级将大矩阵分解为更小、更易处理的块，大幅减少了对大矩阵求逆的需求。