[论文解读] A Nonconforming Finite Element Method for Fourth Order Curl Equations in R^3

Bin Zheng, Qiya Hu+1 more|arXiv (Cornell University)|Jan 31, 2010

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

一句话总结

本文提出了一种用于求解三维磁流体动力学（MHD）模型中第四阶旋度方程的非协调有限元方法。该方法采用一种具有最小自由度的新单元，强制实现旋度的切向连续性，实现了最优收敛速率，并提供了显式基函数以支持实际实现。

ABSTRACT

In this paper we present a nonconforming finite element method for solving fourth order curl equations in three dimensions arising from magnetohydrodynamics models. We show that the method has an optimal error estimate for a model problem involving both curl^2 and curl^4 operators. The element has a very small number of degrees of freedom and it imposes the inter-element continuity along the tangential direction which is appropriate for the approximation of magnetic fields. We also provide explicit formulae of basis functions for this element.

研究动机与目标

• 开发一种稳定且高效的有限元方法，用于求解三维磁流体动力学（MHD）问题中的第四阶旋度算子。
• 解决在MHD模型中对 $(\nabla\times)^4$ 算子进行离散化的挑战，该问题由于自由度过高，使用标准协调元难以处理。
• 构建一种非协调有限元，确保在单元之间旋度的切向连续性，适用于磁场近似。
• 为新单元提供显式基函数，以支持实际实现。
• 在标准Sobolev正则性假设下，建立该方法的最优误差估计。

提出的方法

• 该方法基于分片二次多项式并引入与旋度相关的自由度，构建非协调有限元空间。
• 有限元空间强制实现旋度切向分量的单元间连续性，确保与磁场物理规律的一致性。
• 通过涉及 $ \nabla\times\mathbf{v}_h $ 的体积积分和面积分构造双线性形式 $ a_h(\cdot, \cdot) $，并引入参数 $ \alpha $ 和 $ \beta $ 的稳定项。
• 采用非协调Galerkin弱形式进行分析，通过分部积分和迹不等式推导出一致性和近似误差估计。
• 分析基于第二Strang引理，证明了旋度及其梯度在类似 $ H^1 $ 范数下的最优收敛性。
• 推导并提供了该单元的显式基函数，支持在有限元代码中直接实现。

实验结果

研究问题

• RQ1 能否构造一种非协调有限元方法，以最小自由度高效离散三维第四阶旋度算子？
• RQ2 在单元边界强制实现旋度的切向连续性，是否能为第四阶旋度问题带来最优收敛性？
• RQ3 如何利用非协调方法稳健地近似MHD模型中的 $(\nabla\times)^4$ 算子？
• RQ4 在标准正则性假设下，该方法可达到的最优误差估计是什么？
• RQ5 能否为新单元推导出显式基函数以支持实际实现？

主要发现

• 所提出的非协调有限元方法在能量范数下实现了最优收敛速率，误差被控制在 $ \lesssim h\|\mathbf{u}\|_{4,\Omega} $ 范围内，其中 $ \mathbf{u} \in (H^4(\Omega))^3 $。
• 该方法对解、其旋度以及旋度的梯度均表现出最优误差估计，确保高阶精度。
• 每个四面体单元仅具有30个自由度，相较于需要220个自由度的协调方法，显著降低了计算成本。
• 该方法通过在单元界面强制实现旋度的切向连续性，保持了磁场的物理结构。
• 推导出了基函数的显式公式，支持在有限元软件中直接实现。
• 一致性误差被控制在 $ \lesssim h(\|\nabla\times\Delta(\nabla\times\mathbf{u})\| + \|\nabla(\nabla\times\mathbf{u})\|_{1} + \|\nabla\times\mathbf{u}\| + \|\nabla\times\nabla\times\mathbf{u}\|) $ 范围内，证实了其最优逼近性质。