[论文解读] Convergence proof for the multigrid method of the nonlocal model

Minghua Chen, Weihua Deng|arXiv (Cornell University)|May 18, 2016

Numerical methods in engineering参考文献 19被引用 4

一句话总结

本文为非局部模型在双网格法（TGM）中的收敛性提供了严格的证明，该模型源自微分动力学，其中当网格尺寸 $ h \to 0 $ 时，比值 $ a_0/a_1 $ 无界。通过将对称托普利茨刚度矩阵分解为若干拉普拉斯型矩阵之和，作者建立了双网格法的统一收敛性，并在特定条件下将分析扩展至全多重网格（FMGM）和V型循环格式。数值实验验证了其二阶精度和 $ \mathcal{O}(N\log N) $ 的计算复杂度。

ABSTRACT

Recently, nonlocal models attract the wide interests of scientist. They mainly come from two applied scientific fields: peridyanmics and anomalous diffusion. Even though the matrices of the algebraic equation corresponding the nonlocal models are usually Toeplitz (denote a0 as the principal diagonal element, a1 as the trailing diagonal element, etc). There are still some differences for the models in these two fields. For the model of anomalous diffusion, a0/a1 is uniformly bounded; most of the time, a0/a1 of the model for peridyanmics is unbounded as the stepsize h tends to zero. Based on the uniform boundedness of a0/a1, the convergence of the two-grid method is well established [Chan, Chang, and Sun, SIAM J. Sci. Comput., 19 (1998), pp. 516--529; Pang and Sun, J. Comput. Phys., 231 (2012), pp. 693--703; Chen, Wang, Cheng, and Deng, BIT, 54 (2014), pp. 623--647]. This paper provides the detailed proof of the convergence of the two-grid method for the nonlocal model of peridynamics. Some special cases of the full multigrid and the V-cycle multigrid are also discussed. The numerical experiments are performed to verify the convergence.

研究动机与目标

• 为非局部微分动力学模型中的双网格法（TGM）建立严格的收敛性证明，其中当 $ h \to 0 $ 时 $ a_0/a_1 $ 无界，这与异常扩散模型的情形不同。
• 克服现有收敛性证明依赖于 $ a_0/a_1 $ 有界的局限性，该假设在微分动力学系统中不成立。
• 在特定结构假设下，将收敛性分析扩展至全多重网格（FMGM）和V型循环多重网格方法。
• 通过数值实验验证理论结果，展示最优收敛速率与计算效率。

提出的方法

• 将非局部模型的刚度矩阵表示为若干拉普拉斯型矩阵之和，从而能够分析其谱性质与条件数。
• 双网格法的收敛性证明基于这些拉普拉斯型分量的分解及其性质，特别是其在某些范数下的统一有界性。
• 采用伽辽金插值策略构造粗网格算子，确保各层次上的稳定性与逼近质量。
• 在粗网格校正满足统一收敛条件的假设下，分析全多重网格（FMGM）与V型循环多重网格方案。
• 数值实验采用V型循环算法，平滑参数为 $ \omega_{\text{pre}} = 1 $，$ \omega_{\text{post}} = 1/3 $，并采用基于残差的停止准则 $ \|r^{(i)}\| / \|r^{(0)}\| < 10^{-8} $。
• 对比了两种粗化策略——伽辽金法与网格尺寸加倍法，以评估其鲁棒性与计算成本。

实验结果

研究问题

• RQ1 当 $ h \to 0 $ 时 $ a_0/a_1 $ 无界，是否可证明双网格法在非局部微分动力学模型中实现统一收敛？
• RQ2 刚度矩阵的何种结构特性使得在 $ a_0/a_1 $ 无界时仍能实现收敛性证明？
• RQ3 全多重网格方法在相同条件下是否仍保持统一收敛性？
• RQ4 V型循环多重网格格式是否可对微分动力学问题的特定类别实现统一收敛？
• RQ5 不同粗化策略在实际中如何影响收敛速率与计算成本？

主要发现

• 尽管 $ a_0/a_1 $ 无界，双网格法（TGM）在微分动力学非局部模型中仍实现统一收敛，其证明基于分解为拉普拉斯型矩阵。
• 数值实验确认，对所有测试的 $ \delta $ 值（包括 $ \delta = 1 $，$ \delta = \sqrt{h} $，$ \delta = 5h $，$ \delta = h $）在 $ l_\infty $ 范数下均呈现二阶精度，收敛率接近 2.0。
• 多重网格求解器的计算成本近似为 $ \mathcal{O}(N\log N) $，各层级的平均迭代次数保持稳定。
• 伽辽金法与网格尺寸加倍粗化策略在精度与收敛行为上几乎一致，CPU 时间随 $ N $ 线性增长。
• 当 $ \delta = h $ 时，TGM 每层的平均迭代次数稳定在 18–23 之间，FMGM 为 47–56 之间，表明其鲁棒性。
• 在特定情况下，V型循环多重网格的收敛性已得到验证，前提是粗网格校正满足所需的统一收敛条件，但条件 (52) 的验证仍为开放问题。