We propose, analyze mathematically, and study numerically a novel approach for the finite element approximation of the spectrum of second-order elliptic operators. The main idea is to reduce the stiffness of the problem by subtracting a least-squares penalty on the gradient jumps across the mesh interfaces from the standard stiffness bilinear form. This penalty bilinear form is similar to the known technique used to stabilize finite element approximations in various contexts. The penalty term is designed to dampen the high frequencies in the spectrum and so it is weighted here by a negative coefficient. The resulting approximation technique is called softFEM since it reduces the stiffness of the problem. The two key advantages of softFEM over the standard Galerkin FEM are to improve the approximation of the eigenvalues in the upper part of the discrete spectrum and to reduce the condition number of the stiffness matrix. We derive a sharp upper bound on the softness parameter weighting the stabilization bilinear form so as to maintain coercivity for the softFEM bilinear form. Then we prove that softFEM delivers the same optimal convergence rates as the standard Galerkin FEM approximation for the eigenvalues and the eigenvectors. We next compare the discrete eigenvalues obtained when using Galerkin FEM and softFEM. Finally, a detailed analysis of linear softFEM for the 1D Laplace eigenvalue problem delivers a sensible choice for the softness parameter. With this choice, the stiffness reduction ratio scales linearly with the polynomial degree. Various numerical experiments illustrate the benefits of using softFEM over Galerkin FEM.
翻译:我们提议, 从数学角度分析, 并用数字学来研究一种新颖的方法, 用于第二阶椭圆操作员频谱的有限元素近似值。 主要的想法是减少问题的僵硬度, 从标准硬度双线形界面中减去对网状界面斜度跳跃的最小平方的罚款。 这种惩罚双线形式类似于用来稳定不同背景中有限元素近似值的已知技术。 惩罚用词旨在抑制频谱中的高频度, 因此它在此被负系数加权。 由此产生的近似技术被称为软FEM, 因为它可以降低问题的僵硬度。 软FEM 相对于标准 Galerkin FEM 界面, 软性FEMEM 的两大关键优势在于改善离心率的近似度, 并且用软软软FEM 度的软度递增软度递增。 当我们用软性FEM 软度的软度递增时, 我们用软性FEM 软度的软度递减率 和直线形的直径直径直径直径直径直径直径直径直径直径直径直径直径直径直径直径直径直径的直径直径直径直径直径,,, 用软直径直度直度直度直向右的FEM, 的直径直径直径直径直度直度直度直径直径直径直径直径直径直度直度直度直度直度直径直径直度直度直度直, 。