[논문 리뷰] The Lower Bounds for Eigenvalues of Elliptic Operators --By Nonconforming Finite Element Methods

Jun Hu, Yunqing Huang외 1명|arXiv (Cornell University)|2011. 12. 06.

Advanced Numerical Methods in Computational Mathematics참고 문헌 53인용 수 14

한 줄 요약

이 논문은 대칭 타원형 연산자의 고유값에 대한 하한을 도출하기 위해 비일치 유한요소법(nonconforming finite element methods)을 사용하는 체계적인 방법을 제시한다. 핵심 통찰은 비일치 유한요소 공간의 국소 근사 성질이 전역 연속성 성질을 초월할 경우, 이에 대응하는 이산 고유값이 하한으로 보장된다는 것이다. 이는 충분한 조건이 특정되어 있으며, 기존 및 새로 구성된 요소들에 대해 검증되었다.

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.

연구 동기 및 목표

• 대칭 타원형 연산자의 고유값에 대한 하한을 체계적으로 계산하는 방법을 개발하기 위해.
• 일치 유한요소법이 오직 상한만을 제공하는 데에 한계가 있음에도 불구하고, 비일치 유한요소법이 하한을 도출할 수 있는 조건을 규명하기 위해.
• 일부 비일치 요소는 하한을 생성하지만 다른 요소는 그렇지 않은 관측된 수치 현상의 이면에 이론적 기반을 제공하기 위해.
• 국소 근사 성질이 전역 연속성 성질을 초월하도록 보장함으로써, 하한을 보장하는 새로운 비일치 요소를 설계할 수 있도록 이끌기 위해.
• 대부분의 비일치 요소에 대해 포화 조건(saturation condition)을 증명하여 제안된 방법의 신뢰성을 뒷받침하기 위해.

제안 방법

• 비일치 유한요소 공간의 근사 및 연속성 성질에 기반하여 이산 고유값과 고유함수에 대한 추상적 오차 추정을 수립하기 위해.
• 비일치 유한요소 공간의 국소 근사 성질과 전역 연속성 성질 간의 균형을 포함하는 충분한 조건을 도입하여 고유값 하한을 보장하기 위해.
• 비일치 공간을 일치 유한요소 공간과 연결하기 위해 투영 평균 보간 연산자를 도입하여 버블 함수 기법을 통한 오차 분석을 가능하게 하기 위해.
• 요소 간의 점프와 경계 항을 포함하는 보간 오차 추정을 사용하여 비일치 방법의 일致성 오차를 제한하기 위해.
• Crouzeix-Raviart, Morley, Wilson, 그리고 풍부화된 기울인 Q1 요소와 같은 잘 알려진 비일치 요소들에 대해 충분한 조건을 검증하기 위해.
• Crouzeix-Raviart 요소를 풍부화시켜 충분한 조건를 만족시키는 새로운 비일치 요소를 구성함으로써 고유값에 대한 하한을 보장하기 위해.

실험 결과

연구 질문

• RQ1 비일치 유한요소법이 대칭 타원형 연산자의 고유값에 대해 하한을 생성하는 조건는 무엇인가?
• RQ2 유사한 근사 성질을 가진 비일치 요소들 사이에서 일부는 하한을 생성하고 다른 일부는 상한을 생성하는 이유는 무엇인가?
• RQ3 비일치 방법을 사용하여 하한 고유값을 예측하고 보장할 수 있는 통합된 이론적 프레임워크를 수립할 수 있는가?
• RQ4 기존 비일치 요소들을 어떻게 수정하거나 새로운 요소를 설계하여 하한 조건을 만족시킬 수 있는가?
• RQ5 고유값 하한의 방향을 결정하는 데 있어 국소 근사 성질과 전역 연속성 성질 간의 역할은 무엇인가?

주요 결과

• 대부분의 표준 비일치 요소, 즉 Crouzeix-Raviart, Morley, Wilson, 그리고 풍부화된 기울인 Q1 요소들은 고유값 하한을 위한 충분한 조건를 만족한다.
• Crouzeix-Raviart 요소를 풍부화시켜 구성한 새로운 비일치 요소는 충분한 조건를 만족하므로 고유값에 대한 하한을 보장한다.
• 대부분의 비일치 요소에 대해 고유값 오차의 포화 조건이 성립함을 증명하여 제안된 방법의 신뢰성을 뒷받침한다.
• 이 방법은 체계적인 설계 원칙을 제공한다: 전역 연속성 대비 국소 근사 성질이 뛰어난 비일치 요소는 고유값 하한을 도출한다.
• 이론적 분석은 일부 비일치 방법이 하한을 생성하는 관측된 수치적 행동이 우연이 아니며 수학적으로 기반을 지닌다는 것을 확인한다.
• 투영 평균 보간 연산자는 비일치 및 일치 유한요소 공간을 연결함으로써 효과적인 오차 추정을 가능하게 하며, 핵심 추정의 증명을 촉진한다.