[论文解读] Amiable mixed schemes for fourth order curl equations

Shuo Zhang|arXiv (Cornell University)|Jul 17, 2016

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

一句话总结

本文通过引入辅助变量，将四阶旋度方程转化为适定的混合弱形式，提出了一类稳定且低阶的混合有限元格式。该方法利用 $H^1$、$H({\rm curl})$ 和 $L^2$ 有限元空间，确保了最优收敛率，其正则性结果在凸多面体区域上得到验证，收敛性通过标准变分分析得以证明。

ABSTRACT

In this paper, amiable mixed schemes are presented for two variants of fourth order curl equations. Specifically, mixed formulations for the problems are constructed, which are well-posed in Babuska-Brezzi's sense and admit stable discretizations by finite element spaces of low smoothness and of low degree. The regularities of the mixed formulations and thus equivalently the primal problems are established, and some finite elements examples are given which can exploit the regularity of the solutions to an optimal extent.

研究动机与目标

• 为两类四阶旋度方程变体构造在最小光滑性要求下稳定的、适定的混合弱形式。
• 在凸多面体区域上建立混合弱形式的正则性结果，并证明其与原问题的等价性。
• 设计能充分利用解的正则性、采用低次、低光滑性有限元空间的有限元离散格式。
• 通过确保嵌套网格上代数与拓扑嵌套关系，实现高效可扩展的求解策略。
• 为未来工作在参数鲁棒离散化与特征值问题方面提供基础。

提出的方法

• 引入辅助变量，将四阶旋度问题转化为一组二阶方程，形成混合弱形式。
• 利用 $L^2$、$H({\rm curl})$ 和 $H^1$ 类型有限元空间，构造在 Babuška-Brezzi 意义下适定的混合变分格式。
• 通过证明在凸多面体上向量势及其旋度具有 $H^2$ 正则性，建立混合格式的正则性。
• 使用标准有限元（如 Nédélec 元和分片多项式）实现稳定且收敛的离散化。
• 应用对偶论证与先验误差估计，推导能量范数与 $L^2$ 范数下的收敛速率。
• 确保在嵌套网格上有限元空间具有代数与拓扑嵌套性，以兼容多水平方法。

实验结果

研究问题

• RQ1 能否使用低光滑性、低次有限元构造四阶旋度方程的稳定混合格式？
• RQ2 混合格式解的正则性如何？其是否确认了在凸多面体区域上原问题解的 $H^2$ 正则性？
• RQ3 混合格式能否以一种支持最优收敛率且可借助现有预conditioner 高效求解的方式离散化？
• RQ4 在变分框架中，如何正确强制解的旋度的边界条件？
• RQ5 该格式未来能否扩展至参数鲁棒或特征值问题？

主要发现

• 对于两种问题变体（A）和（B），混合格式在 Babuška-Brezzi 条件下均适定且稳定，使用 $L^2$、$H({\rm curl})$ 和 $H^1$ 空间。
• 解 $\undertilde{u}$ 及其旋度 $\nabla\times\undertilde{u}$ 在凸多面体区域上具有 $H^2$ 正则性，验证了先前研究中的假设。
• 实现了最优收敛率：$\|\undertilde{\varphi} - \undertilde{\varphi}_h\|_{0,\Omega} \leq Ch^2\|\undertilde{f}\|_{0,\Omega}$，其中 $\undertilde{\varphi}$ 为混合变量。
• 对于问题（B），在 $H({\rm curl})$-范数下的收敛率为 $O(h)$，若 $\undertilde{f} \in \undertilde{H}^1(\Omega)$，则可提升至 $O(h^2)$。
• 所用有限元空间（如 $N^2_{h0}$ 和 $\mathcal{L}^2_{h0}$）在嵌套网格上支持代数与拓扑嵌套，从而支持高效多级求解器。
• 该格式可使用标准有限元软件包实现，并与现有最优预conditioner 兼容。