[论文解读] A Multilevel Correction Scheme for Nonsymmetric Eigenvalue Problems by Finite Element Methods

Hehu Xie, Zhimin Zhang|arXiv (Cornell University)|May 23, 2015

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

一句话总结

该论文提出了一种基于有限元法的非对称特征值问题多水平校正方案，可在无需假设特征值上升度为1的情况下，实现高效且精确的逼近。通过仅在最粗网格上求解小型特征值问题，并在更细的网格空间上求解源问题，该方法实现了最优收敛速率，并显著提升了计算效率，尤其在与多重网格或自适应加密结合时效果更佳。

ABSTRACT

A multilevel correction scheme is proposed to solve defective and nodefective of nonsymmetric partial differential operators by the finite element method. The method includes multi correction steps in a sequence of finite element spaces. In each correction step, we only need to solve two source problems on a finer finite element space and two eigenvalue problems on the coarsest finite element space. The accuracy of the eigenpair approximation is improved after each correction step. This correction scheme improves overall efficiency of the finite element method in solving nonsymmetric eigenvalue problems.

研究动机与目标

• 开发一种高效的有限元方法，用于求解流体动力学和稳定性分析中出现的非对称特征值问题。
• 将此前仅用于对称问题的多水平校正技术扩展至非对称算子，且无需假设特征值上升度为1。
• 通过减少在精细有限元空间中求解大规模特征值问题的成本，提升计算效率。
• 实现与多重网格法和自适应有限元法的无缝耦合，以实现最优收敛和并行化。

提出的方法

• 该方法采用一系列逐渐细化的嵌套有限元空间，从最粗的层次开始。
• 在每个校正步骤中，利用前一个特征函数作为载荷向量，在增强的有限元空间中求解源问题。
• 随后在包含已校正特征函数的增强空间中求解一个新的特征值问题，但仅使用最粗空间进行特征值计算。
• 通过一系列校正步骤，迭代地改进特征对的逼近，同时校正特征值和特征函数分量。
• 该方案通过在细网格上用更便宜的源问题替代完整的特征值问题，并在最粗空间上求解小型特征值问题，避免了在细网格上求解完整特征值问题。
• 该方法与边界值问题的多重网格求解器兼容，并与基于ZZ恢复的后验误差估计器结合使用自适应加密。

实验结果

研究问题

• RQ1 能否设计一种无需假设特征值上升度为1的非对称特征值问题多水平校正方案？
• RQ2 如何利用有限元方法降低求解大规模非对称特征值问题的计算成本？
• RQ3 当应用于具有奇异特征函数的问题（如L形区域）时，多水平校正方案能否保持最优收敛速率？
• RQ4 当微分算子中包含复数向量场时，该方法表现如何？
• RQ5 该方案能否与多重网格和自适应有限元技术有效结合以提升效率？

主要发现

• 即使在非凸区域（如L形区域）上存在奇异特征函数的问题中，多水平校正方案仍能实现特征值和特征函数逼近的最优收敛速率。
• 数值结果证实，该方法在拟均匀网格上对第一个特征值的收敛速率保持在$ h^{4/3} $的最优阶，与理论预期一致。
• 在自适应网格上，该方法实现了预期速率的最优收敛，初始网格对最终精度的影响极小。
• 该方法在微分算子中包含复数向量场时表现良好，数值结果表明前六个特征值均实现快速收敛。
• 在计算前六个特征值之和时，多空间方法的精度高于多重网格方法，表明其在高精度场景下具有优势。
• 该方案极为高效：仅在最粗空间上求解小型特征值问题，并在更细空间上求解源问题，即可达到在最细空间上求解完整特征值问题的精度。