We revisit the problem of spurious modes that are sometimes encountered in partial differential equations discretizations. It is generally suspected that one of the causes for spurious modes is due to how boundary conditions are treated, and we use this as the starting point of our investigations. By regarding boundary conditions as algebraic constraints on a differential equation, we point out that any differential equation with homogeneous boundary conditions also admits a typically infinite number of hidden or implicit boundary conditions. In most discretization schemes, these additional implicit boundary conditions are violated, and we argue that this is what leads to the emergence of spurious modes. These observations motivate two definitions of the quality of computed eigenvalues based on violations of derivatives of boundary conditions on the one hand, and on the Grassmann distance between subspaces associated with computed eigenspaces on the other. Both of these tests are based on a standardized treatment of boundary conditions and do not require a priori knowledge of eigenvalue locations. The effectiveness of these tests is demonstrated on several examples known to have spurious modes. In addition, these quality tests show that in most problems, about half the computed spectrum of a differential operator is of low quality. The tests also specifically identify the low accuracy modes, which can then be projected out as a type of model reduction scheme.
翻译:我们重新审视了偏微分方程离散化中有时遇到的虚假模态问题。通常认为,虚假模态的成因之一在于边界条件的处理方式,本文以此作为研究的出发点。通过将边界条件视为微分方程的代数约束,我们指出,任何具有齐次边界条件的微分方程通常还包含无限多个隐藏或隐式边界条件。在大多数离散化方案中,这些额外的隐式边界条件被违反,我们认为这正是导致虚假模态出现的原因。这些观察结果促使我们基于两方面定义了计算特征值的质量:一方面是基于边界条件导数违反程度,另一方面是基于与计算特征空间相关的子空间之间的Grassmann距离。这两种检验方法均基于标准化的边界条件处理,无需预先知道特征值的位置。我们在多个已知存在虚假模态的示例上证明了这些检验方法的有效性。此外,这些质量检验显示,在大多数问题中,微分算子计算谱中约一半为低质量谱。检验还能具体识别出低精度的模态,这些模态可作为模型降阶方案的一种类型被投影剔除。