[论文解读] Inexact spectral deferred corrections using single-cycle multigrid

Robert Speck, Daniel Ruprecht+3 more|arXiv (Cornell University)|Jan 30, 2014

Advanced Numerical Methods in Computational Mathematics参考文献 11被引用 3

一句话总结

本文提出了一种不精确谱延迟校正（SDC）方法，利用单周期多重网格法高效求解隐式配点法中产生的大型线性系统。通过将SDC的迭代校正框架与多重网格预条件技术相结合，该方法在降低计算成本的同时实现了高阶精度，显著提升了刚性初值问题的收敛速度，优于标准SDC方法。

ABSTRACT

Implicit integration methods based on collocation are attractive for a num-ber of reasons, e.g. their ideal (for Gauss-Legendre nodes) or near ideal(Gauss-Radau or Gauss-Lobatto nodes) order and stability properties. How-ever, straightforward application of a collocation formula with M nodes toan initial value problem with dimension d requires the solution of one largeMd×Md system of nonlinear equations.Spectral deferred correction (SDC) methods, introduced by Dutt et al.[2000],arean attractiveapproachforiterativelycomputing the solution to thecollocation problem using a low-order method (like implicit or IMEX Euler)as a building block. Instead of solving one huge system of size Md × Md,SDC iteratively solves M smaller d×d systems to approximate the solutionof the full system (see also the discussion in Huang et al. [2006]). It has beenshown e.g. by Xia et al. [2007] that each iteration/sweep of SDC raises theorder by one, so that SDC with k iterations and a first-order base method isof order k, up to the order of the underlying collocation formula. Therefore,to achieve formal order p, SDC requires p/2 nodes and p iterations and thusp

研究动机与目标

• 为解决刚性初值问题中隐式配点法所产生大型耦合系统带来的高计算成本问题。
• 通过减少收敛所需的迭代次数，提升谱延迟校正（SDC）方法的效率。
• 将多重网格技术整合到SDC中，通过有效预条件化每次SDC扫描中的线性系统，加速收敛。
• 在降低多维问题中求解大型系统计算负担的同时，保持高阶精度。
• 证明不精确SDC结合单周期多重网格可在最小额外开销下实现最优收敛速率。

提出的方法

• 以谱延迟校正（SDC）作为基础框架，通过在M个配点上使用低阶基方法（如隐式欧拉法）迭代校正解。
• 在每次SDC扫描中，应用单周期多重网格法求解所生成的线性系统，以多重网格预条件方法替代直接求解器或迭代求解器。
• 通过构建粗网格层次结构，加速线性系统求解的收敛速度，减少每次SDC扫描所需的迭代次数。
• 通过容忍多重网格解中的可控误差，引入不精确求解，实现精度与计算成本之间的平衡。
• 通过确保各次扫描中校正步骤的一致性，保持底层配点公式的精度阶数。
• 结合SDC的迭代特性与多重网格的快速收敛性，实现计算工作量显著降低的高阶时间积分。

实验结果

研究问题

• RQ1 单周期多重网格能否有效预条件SDC中的线性系统，在不损失精度的前提下加速收敛？
• RQ2 与标准SDC相比，结合多重网格的不精确SDC方法在收敛速度和计算成本方面表现如何？
• RQ3 在多重网格循环中容忍不精确求解对SDC方法整体精度阶数和稳定性有何影响？
• RQ4 SDC与多重网格的结合能否在显著减少所需迭代次数的同时保持高阶精度？
• RQ5 所提出方法在问题维度和配点数量增加时是否具备良好的可扩展性？

主要发现

• 与精确SDC相比，不精确SDC结合单周期多重网格在每次迭代的计算成本显著降低的同时，实现了相近的收敛速率。
• 采用多重网格预条件可减少每次SDC扫描所需的迭代次数，尤其在条件数较高的刚性问题中效果显著。
• 即使在线性系统求解过程中采用不精确解，该方法仍能保持底层配点公式的形式精度阶数。
• 单周期多重网格方法在精度与效率之间提供了有利的权衡，使得收敛速度优于标准SDC，且不损失精度阶数。
• 每次SDC迭代的计算工作量按多重网格收敛速率成比例减少，从而在大规模问题中实现显著节省。
• 数值实验表明，不精确SDC结合多重网格在各类刚性常微分方程及不断增大的问题规模下，均保持了稳定性和准确性。