[论文解读] A uniform inf--sup condition with applications to preconditioning

Kent‐André Mardal, Joachim Schöberl+1 more|arXiv (Cornell University)|Jan 6, 2012

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

一句话总结

该论文在不依赖区域凸性的前提下，为参数依赖的斯托克斯问题建立了统一的 inf–sup 条件，从而能够构造出在不同扰动参数和网格加密下均保持鲁棒性的统一预条件子。其关键贡献是提出了一种针对二维泰勒-胡德元的新型 $L^2$-有界、与散度算子可交换的投影算子，该算子通过与最低阶内德éléc 边元之间的联系推导得出。

ABSTRACT

A uniform inf-sup condition related to a parameter dependent Stokes problem is established. Such conditions are intimately connected to the construction of uniform preconditioners for the problem, i.e., preconditioners which behave uniformly well with respect to variations in the model parameter as well as the discretization parameter. For the present model, similar results have been derived before, but only by utilizing extra regularity ensured by convexity of the domain. The purpose of this paper is to remove this artificial assumption. As a byproduct of our analysis, in the two dimensional case we also construct a new projection operator for the Taylor-Hood element which is uniformly bounded in $L^2$ and commutes with the divergence operator. This construction is based on a tight connection between a subspace of the Taylor-Hood velocity space and the lowest order Nedelec edge element.

研究动机与目标

• 消除以往为实现参数依赖斯托克斯问题统一预条件子所必需的非物理凸性假设。
• 为奇异扰动问题 (1.1) 建立一个与扰动参数 $\epsilon$ 和网格尺寸 $h$ 无关的统一 inf–sup 条件。
• 为二维泰勒-胡德有限元方法构造一个在 $L^2$ 中一致有界、与散度算子可交换的投影算子。
• 通过数值证据支持理论解释，说明预条件子在非凸区域中的鲁棒性。
• 通过利用 Bogovskiï¿ 算子和形状正则三角剖分，将统一预条件子技术的适用范围扩展至一般的利普希茨区域。

提出的方法

• 利用 Bogovskiï¿ 算子作为散度算子的右逆，在连续情形下建立连续性与稳定性。
• 通过揭示其速度空间与最低阶内德éléc 边元之间的紧密联系，为泰勒-胡德元构造新的投影算子。
• 在宏观单元上定义一个对称正定矩阵 $M$，通过精心选择参数 $\gamma$，确保其特征值不随网格尺寸和单元形状变化。
• 利用形状正则性对参数 $\gamma$ 从上下两方面进行有界控制，确保三角剖分中所有单元的一致性。
• 应用 Gershgorin 圆定理证明矩阵 $M$ 的特征值下界，从而保证双线性形式的强制性。
• 将宏观单元上的局部投影与全局构造的与散度可交换、$L^2$-有界的映射 $\tilde{\Phi}_h$ 相结合，实现统一的 inf–sup 稳定性。

实验结果

研究问题

• RQ1 是否可以在不假设区域凸性的前提下，为参数依赖的斯托克斯问题建立统一的 inf–sup 条件？
• RQ2 是否可能为二维泰勒-胡德元构造一个在 $L^2$ 中一致有界、与散度算子可交换的投影算子？
• RQ3 如何利用泰勒-胡德速度空间与最低阶内德éléc 边元之间的联系来构建此类投影算子？
• RQ4 Bogovskiï¿ 算子在确保连续问题统一稳定性方面发挥什么作用？
• RQ5 所得到的预条件子是否能在非凸区域中保持对 $\epsilon$ 和网格加密的统一条件数？

主要发现

• 在一般利普希茨区域上，为参数依赖的斯托克斯问题 (1.1) 建立了统一的 inf–sup 条件，无需依赖凸性假设。
• 统一的 inf–sup 条件确保了解算子在 $\epsilon$ 和 $h$ 中的一致有界性，从而支持统一预条件子的构造。
• 为泰勒-胡德元构造了一种在 $L^2$ 中一致有界、且在二维中与散度算子可交换的新投影算子。
• 该构造依赖于泰勒-胡德速度空间与最低阶内德éléc 边元之间的精确关系，从而实现一个对称正定矩阵 $M$，其特征值下界为 $c_1|T|$。
• 证明了矩阵 $M$ 是对称正定的，且其特征值下界为 $c_1|T|$，其中 $c_1 = \alpha \min(\gamma, \gamma^{-1})/2$，且 $\gamma$ 的取值独立于 $h$ 和单元形状。
• 全局投影映射 $\tilde{\Phi}_h$ 实现了统一的 $L^2$ 稳定性，满足 $\sup_{z \in Z_h^0} \frac{\langle v, z \rangle}{\|z\|_{L^2}} \geq c_0 \|v\|_{L^2}$，从而验证了统一的 inf–sup 条件。