[论文解读] The Lower Bounds for Eigenvalues of Elliptic Operators --By Nonconforming Finite Element Methods

Jun Hu, Yunqing Huang+1 more|arXiv (Cornell University)|Dec 6, 2011

Advanced Numerical Methods in Computational Mathematics参考文献 53被引用 14

一句话总结

本文提出了一种系统方法，利用非协调有限元法为对称椭圆算子的特征值提供下界。关键洞见在于，当非协调有限元空间的局部逼近性质优于其整体连续性时，所得离散特征值被保证为下界，并为此识别并验证了一个充分条件，适用于已知及新构造的单元。

ABSTRACT

The aim of the paper is to introduce a new systematic method that can produce lower bounds for eigenvalues. The main idea is to use nonconforming finite element methods. The general conclusion herein is that if local approximation properties of nonconforming finite element spaces $V_h$ are better than global continuity properties of $V_h$, corresponding methods will produce lower bounds for eigenvalues. More precisely, under three conditions on continuity and approximation properties of nonconforming finite element spaces we first show abstract error estimates of approximate eigenvalues and eigenfunctions. Subsequently, we propose one more condition and prove that it is sufficient to guarantee nonconforming finite element methods to produce lower bounds for eigenvalues of symmetric elliptic operators. As one application, we show that this condition hold for most nonconforming elements in literature. As another important application, this condition provides a guidance to modify known nonconforming elements in literature and to propose new nonconforming elements. In fact, we enrich locally the Crouzeix-Raviart element such that the new element satisfies the condition; we propose a new nonconforming element for second order elliptic operators and prove that it will yield lower bounds for eigenvalues. Finally, we prove the saturation condition for most nonconforming elements.

研究动机与目标

• 开发一种系统化方法，用于计算对称椭圆算子特征值的下界。
• 识别非协调有限元法产生下界的条件，克服协调方法仅能提供上界这一局限。
• 为观察到的数值现象提供理论基础：某些非协调单元产生特征值下界，而另一些则不产生。
• 通过确保局部逼近优于整体连续性，指导新非协调单元的构造，以保证下界。
• 证明大多数非协调单元的饱和条件，支持所提方法的可靠性。

提出的方法

• 基于非协调有限元空间的逼近性和连续性性质，建立离散特征值与特征函数的抽象误差估计。
• 引入一个充分条件，涉及非协调有限元空间中局部逼近与整体连续性性质之间的平衡，以保证特征值下界。
• 应用投影平均插值算子，将非协调空间与协调有限元空间关联，通过泡函数技术实现误差分析。
• 利用涉及单元界面跳跃和边界项的插值误差估计，以界定非协调方法中的一致性误差。
• 验证已知非协调单元（如Crouzeix-Raviart、Morley、Wilson及增强旋转Q1单元）的充分条件。
• 通过增强Crouzeix-Raviart单元构造一个新的非协调单元，使其满足充分条件，从而保证特征值的下界。

实验结果

研究问题

• RQ1 在何种条件下，非协调有限元法能为对称椭圆算子的特征值提供下界？
• RQ2 为何某些非协调单元产生下界，而另一些则产生上界，尽管其逼近性质相似？
• RQ3 能否建立一个统一的理论框架，以预测并保证使用非协调方法的特征值下界？
• RQ4 如何修改现有非协调单元或设计新单元，以满足下界条件？
• RQ5 局部逼近与整体连续性在决定特征值下界方向中的作用是什么？

主要发现

• 大多数标准非协调单元（包括Crouzeix-Raviart、Morley、Wilson及增强旋转Q1单元）均满足特征值下界所需的充分条件。
• 通过增强Crouzeix-Raviart单元构造的新非协调单元满足充分条件，因而能保证特征值的下界。
• 证明了大多数非协调单元的特征值误差饱和条件成立，支持该方法的可靠性。
• 该方法提供了一套系统化的设计原则：当非协调单元的局部逼近优于整体连续性时，将产生特征值下界。
• 理论分析证实，某些非协调方法产生下界的观察到的数值行为并非偶然，而是具有数学基础。
• 投影平均插值算子通过连接非协调与协调有限元空间，实现了有效的误差估计，促进了关键估计的证明。