[논문 리뷰] The weak Galerkin method for eigenvalue problems

Hehu Xie, Qilong Zhai외 1명|arXiv (Cornell University)|2015. 08. 21.

Advanced Numerical Methods in Computational Mathematics참고 문헌 1인용 수 6

한 줄 요약

이 논문은 약한 함수와 약한 도함수를 활용하여 고차수 수렴을 달성하는 약한 갈레르킨(WG) 유한요소 방법을 제안한다. 이 방법은 라플라스 고유값 문제를 해결하는 데 사용되며, 고유값에 대한 하한 근사값을 제공하고 수렴 속도가 2를 초과하여 고전적인 비일치 방법이 제한하는 2차 수렴보다 크게 향상된다.

ABSTRACT

This article is devoted to computing the eigenvalue of the Laplace eigenvalue problem by the weak Galerkin (WG) finite element method with emphasis on obtaining lower bounds. The WG method is on the use of weak functions and their weak derivatives defined as distributions. Weak functions and weak derivatives can be approximated by polynomials with various degrees. Different combination of polynomial spaces leads to different WG finite element methods, which makes WG methods highly flexible and efficient in practical computation. We establish the optimal-order error estimates for the WG finite element approximation for the eigenvalue problem. Comparing with the classical nonconforming finite element method which can just provide lower bound approximation by linear elements with only the second order convergence, the WG methods can naturally provide lower bound approximation with a high order convergence (larger than $2$). Some numerical results are also presented to demonstrate the efficiency of our theoretical results.

연구 동기 및 목표

• 라플라스 연산자의 고유값에 대해 신뢰할 수 있는 하한 근사값을 제공하는 약한 갈레르킨 유한요소 방법을 개발하는 것.
• 고전적인 비일치 방법이 하한 근사값에 대해 오직 2차 수렴에만 도달할 수 있다는 한계를 극복하는 것.
• 라플라스 고유값 문제에 적용된 WG 방법에 대해 최적의 오차 추정치를 수립하는 것.
• 다각형 및 L자형 영역에서의 수치 실험을 통해 이 방법의 고차수 수렴성과 강건성을 입증하는 것.
• 다각형 메쉬에서 불연속 다항식 근사를 사용하는 탄력적이고 고차수의 수치 프레임워크를 고유값 문제에 제공하는 것.

제안 방법

• 이 방법은 유한요소 분할 위에서 분포로 정의된 약한 함수와 약한 도함수를 사용하여 다각형 요소에서 불연속 근사를 가능하게 한다.
• 다양한 차수의 다항식 공간을 사용하여 약한 함수와 약한 도함수를 근사함으로써 높은 유연성과 고차수 정확도를 달성한다.
• 약한 연속성을 강제하고 이산 형식의 안정성을 확보하기 위해 안정화 항을 도입한다.
• 이산 고유값 문제를 약한 기울기와 약한 함수의 내적을 포함하는 변분 형식에서 유도한다.
• 다각형 및 L자형 영역에서 동일한 딜레르흐 경계 조건을 가진 표준 고유값 문제에 이 방법을 적용한다.
• 수렴성과 하한 성질을 평가하기 위해 다항식 차수(k=1,2)와 메쉬 크기(h=1/4에서 1/128)를 다양하게 설정한 수치적 방법을 구현한다.

실험 결과

연구 질문

• RQ1 약한 갈레르킨 방법은 고차수 수렴을 보이며 라플라스 고유값에 대한 하한 근사값을 제공할 수 있는가?
• RQ2 WG 방법의 수렴 속도는 하한 정확도 측면에서 고전적인 비일치 유한요소 방법과 비교해 어떻게 다른가?
• RQ3 라플라스 고유값 문제의 맥락에서 고유쌍에 대한 WG 근사의 최적 오차 추정치는 무엇인가?
• RQ4 L자형 영역과 같이 볼록이 아닌 영역에서도 WG 방법은 고유값에 대해 하한 성질을 유지하는가?
• RQ5 WG 방법은 하한 성질을 유지하면서도 2차 수렴을 초월하는 수렴 속도를 달성할 수 있는가?

주요 결과

• 약한 갈레르킨 방법은 고유값 근사에 대해 최적의 오차 추정치를 달성하며, 고차수 다항식 근사에 대해 수렴 속도가 2를 초과한다.
• k=1(선형 요소)인 경우, 이 방법은 고전적인 비일치 방법과 일치하는 2차 수렴을 보이며 하한 근사값을 제공한다.
• k=2(이차 요소)인 경우, 이 방법은 수렴 속도가 2를 초과하여 선형 요소보다 더 높은 정확도를 달성한다.
• L자형 영역에서의 수치 결과는 메쉬 크기 h가 감소함에 따라 이산 고유값 λj,h가 단조적으로 증가함을 보여주며, 하한 성질을 확인한다.
• 계산된 고유값은 정확한 고유값으로 아래에서 수렴하며, k=2와 ε=0일 때 h→0으로 갈수록 λ1,h는 약 9.6392에 수렴한다.
• 이 방법은 다양한 메쉬 세분화와 영역 기하학, 특히 특이 영역인 L자형 영역에서도 강건한 성능을 유지한다.