[论文解读] Multigrid Methods for Hellan-Herrmann-Johnson Mixed Method of Kirchhoff Plate Bending Problems

Long Chen, Jun Hu+1 more|arXiv (Cornell University)|Aug 18, 2015

Advanced Numerical Methods in Computational Mathematics参考文献 45被引用 8

一句话总结

本文为基尔霍夫板弯曲问题的赫兰-赫尔曼-约翰逊（HHJ）混合有限元离散化提出了一种V型多重网格法，采用基于核空间稳定多水平分解的乘法Schwarz平滑子。关键贡献在于：仅通过一次平滑步骤即可实现均匀收敛，且无需假设完整的椭圆正则性，该成果通过精确序列和强化柯西-施瓦茨不等式得以实现。

ABSTRACT

A V-cycle multigrid method for the Hellan-Herrmann-Johnson (HHJ) discretization of the Kirchhoff plate bending problems is developed in this paper. It is shown that the contraction number of the V-cycle multigrid HHJ mixed method is bounded away from one uniformly with respect to the mesh size. The uniform convergence is achieved for the V-cycle multigrid method with only one smoothing step and without full elliptic regularity. The key is a stable decomposition of the kernel space which is derived from an exact sequence of the HHJ mixed method, and the strengthened Cauchy Schwarz inequality. Some numerical experiments are provided to confirm the proposed V-cycle multigrid method. The exact sequences of the HHJ mixed method and the corresponding commutative diagram is of some interest independent of the current context.

研究动机与目标

• 为基尔霍夫板弯曲问题HHJ混合方法离散化所产生的鞍点系统开发一种高效的多重网格求解器。
• 消除对完整椭圆正则性的依赖，使方法在非凸区域也能收敛。
• 实现仅需一次平滑步骤的V型多重网格法的均匀收敛，避免通常需要大量平滑步骤的限制。
• 建立与HHJ格式中散度-散度算子相关的核空间的稳定多水平分解。

提出的方法

• 设计一种基于散度算子核空间多水平分解的乘法Schwarz平滑子的V型多重网格算法。
• 利用HHJ混合有限元方法在连续和离散层面的精确序列结构。
• 通过离散精确序列和赫兰-赫尔曼-约翰逊混合格式的性质，构建核空间的稳定分解。
• 证明强化柯西-施瓦茨不等式，以支持稳定性与收敛性分析。
• 由于HHJ方法中有限元空间具有嵌套性，采用标准插值 prolongation 算子。
• 收敛性分析基于文献[21]的框架，经适配后应用于HHJ混合方法并仅使用最少的平滑步骤。

实验结果

研究问题

• RQ1 V型多重网格法能否在不假设完整椭圆正则性的前提下，对基尔霍夫板弯曲问题的HHJ混合方法实现均匀收敛？
• RQ2 是否可能仅通过一次平滑步骤即实现均匀收敛，而非需要多次或大量平滑步骤？
• RQ3 如何为HHJ混合方法中的散度-散度算子核空间构造稳定的多水平分解？
• RQ4 精确序列与强化柯西-施瓦茨不等式在确保多重网格法均匀收敛中起到何种作用？
• RQ5 所提出的方法在非凸区域（如L形区域）是否仍能保持最优收敛速率？

主要发现

• V型多重网格法的收缩数关于网格尺寸均匀远离1，确保了均匀收敛。
• 该方法仅通过一次预平滑和一次后平滑步骤即可实现均匀收敛，与通常需要大量平滑步骤的方法相比，显著降低了计算成本。
• 在正方形和L形区域上的数值实验表明，随着网格细化，迭代次数几乎保持不变，验证了均匀收敛性。
• 对于单位正方形且ν=0.3的情况，‖u−uₕ‖₀ 和 |u−uₕ|₁ 的收敛率分别为 O(h²) 和 O(h)，证实了最优收敛阶。
• 表1中的迭代次数在网格细化过程中保持稳定：正方形区域第6层为23步，L形区域(1,1)平滑下为20步。
• 即使在缺乏完整正则性的情况下，方法仍表现良好，如f=1且ν=0的L形区域测试中，迭代次数保持低且稳定。