[논문 리뷰] Fully Computable Error Bounds for Eigenvalue Problem

Hehu Xie, Meiling Yue외 1명|arXiv (Cornell University)|2016. 01. 07.

Advanced Numerical Methods in Computational Mathematics인용 수 7

한 줄 요약

이 논문은 일반 메esh에서 conforming 방법을 사용한 고유값 문제의 유한요소 근사에 대해 완전히 계산 가능한 오차 한계를 제시한다. 보완 에너지 방법과 후행 오차 추정 기법을 조합함으로써, 고유함수 오차에 대한 보장된 상한과 첫 번째 고유값에 대한 보장된 하한을 확보하였으며, 준균일 및 정규 메쉬에서 점점 정확해지고 강건함을 입증하였다.

ABSTRACT

This paper is concerned with the computable error estimates for the eigenvalue problem which is solved by the general conforming finite element methods on the general meshes. Based on the computable error estimate, we can give an asymptotically lower bound of the general eigenvalues. Furthermore, we also give a guaranteed upper bound of the error estimates for the first eigenfunction approximation and a guaranteed lower bound of the first eigenvalue based on computable error estimator. Some numerical examples are presented to validate the theoretical results deduced in this paper.

연구 동기 및 목표

• 일반 메쉬에서 conforming 유한요소 방법을 사용한 고유쌍 근사에 대해 계산 가능한 후행 오차 추정을 개발한다.
• 첫 번째 고유함수 근사의 오차에 대한 보장된 상한을 제공한다.
• 계산 가능한 추정기법을 사용하여 첫 번째 고유값에 대한 보장된 하한을 유도한다.
• 일반 메쉬에서 일반 고유값에 대해 오차 추정의 점점 정확해짐과 하한 추정의 점점 정확해짐을 보장한다.
• 선형 및 이차 유한요소를 사용하여 적응형 및 정규 메쉬에서 방법을 수치적으로 검증한다.

제안 방법

• 고유함수 오차에 대한 계산 가능한 오차 추정기법을 구성하기 위해 보완 에너지 방법을 사용한다.
• 에너지 노름 오차 $\|u - u_h\|_a$에 대한 상한을 도출하기 위해 후행 오차 추정 기법을 적용한다.
• 오차 추정기법을 계산하기 위해 유한요소 공간 $\mathbf{W}_h^0$ 및 $\mathbf{W}_h^1$에서 이중 문제를 해결한다.
• 계산된 오차 추정기법과 $\lambda_2$에 대한 하한을 사용하여 첫 번째 고유값 $\lambda_1$에 대한 보장된 하한을 도출한다.
• 국소 요소 기반 기여로부터 전체 후행 오차 추정기 $\eta_{\rm ad}(\lambda_h, u_h)$를 구성한다.
• 선형 및 이차 요소를 사용하여 적응형 및 균일 메쉬에서 L자형 도메인을 대상으로 방법을 검증한다.

실험 결과

연구 질문

• RQ1 일반 메쉬에서 conforming 유한요소를 사용한 고유쌍 근사에 대해 완전히 계산 가능한 후행 오차 한계를 구성할 수 있는가?
• RQ2 이 방법은 고유함수 오차에 대한 보장된 상한과 첫 번째 고유값에 대한 보장된 하한을 제공할 수 있는가?
• RQ3 비균일 또는 정규 메쉬에서 오차 추정은 점점 정확해지고 강건한가?
• RQ4 이중 유한요소 공간의 선택($\mathbf{W}_h^0$ 대비 $\mathbf{W}_h^1$)은 오차 추정기와 하한의 정확성에 어떤 영향을 미치는가?
• RQ5 이 방법은 밀도함수 이론에서 발생하는 비선형 고유값 문제로 확장될 수 있는가?

주요 결과

• 제안된 계산 가능한 오차 추정기 $\eta_{\rm ad}(\lambda_h, u_h)$는 일반 메쉬에서 고유쌍 근사에 대해 효율적이고 점점 정확하다.
• 선형 및 이차 유한요소에 대해 진짜 에너지 오차 $\|u - u_h\|_a$를 신뢰성 있게 봉인하는 상한 $\eta_h^U(\lambda_h, u_h, \mathbf{y}_h^*)$가 성립한다.
• 계산된 하한 $\lambda_h^L$은 $\mathbf{W}_h^0$ 또는 $\mathbf{W}_h^1$에서 이중 문제를 풀었는지 여부와 관계없이 첫 번째 고유값 $\lambda_1$에 대한 보장된 하한이다.
• L자형 도메인에서의 수치 결과는 후행 오차 추정기의 효율성과 선형 및 이차 요소에 대한 보장된 하한의 신뢰성을 확인한다.
• 이 방법은 적응형 및 정규 메쉬에서 높은 정확도를 유지하며, 그 일반성과 실용성을 입증한다.
• 이 방법은 스테클로프, 스토크스 문제 및 밀도함수 이론에서 발생하는 비선형 고유값 문제와 같은 다른 고유값 문제로도 확장 가능하다.