Sign-preserving of principal eigenfunctions in P1 finite element approximation of eigenvalue problems of second-order elliptic operators

Abstract

This paper is concerned with the P1 finite element approximation of the eigenvalue problem of second-order elliptic operators subject to the Dirichlet boundary condition. The focus is on the preservation of basic properties of the principal eigenvalue and eigenfunctions of continuous problems. It is shown that when the stiffness matrix is an irreducible $M$-matrix, the algebraic eigenvalue problem maintains those properties such as the smallest eigenvalue being real and simple and the corresponding eigenfunctions being either positive or negative inside the physical domain. Mesh conditions leading to such a stiffness matrix are also studied. A sufficient condition is that the mesh is simplicial, acute when measured in the metric specified by the inverse of the diffusion matrix, and interiorly connected. The acute requirement can be replaced by the Delaunay condition in two dimensions. Numerical results are presented to verify the theoretical findings.

22 pages

References in corpus (6)

• Preserving Monotonicity in Anisotropic Diffusion
• Tiling space and slabs with acute tetrahedra
• Computing the lower and upper bounds of Laplace eigenvalue problem: by combining conforming and nonconforming finite element methods
• Maximum principle in linear finite element approximations of anisotropic diffusion-convection-reaction problems
• Conditioning of Finite Element Equations with Arbitrary Anisotropic Meshes
• The Lower Bounds for Eigenvalues of Elliptic Operators --By Nonconforming Finite Element Methods

Cited by in corpus (4)

• Conditioning of Finite Volume Element Method for Diffusion Problems with General Simplicial Meshes
• A study on anisotropic mesh adaptation for finite element approximation of eigenvalue problems with anisotropic diffusion operators
• On enforcing maximum principles and achieving element-wise species balance for advection-diffusion-reaction equations under the finite element method
• On mesh restrictions to satisfy comparison principles, maximum principles, and the non-negative constraint: Recent developments and new results