[논문 리뷰] The $p$- and $hp$-versions of the virtual element method for elliptic eigenvalue problems

Certik, O., Francesca Gardini외 3명|arXiv (Cornell University)|2018. 12. 21.

Advanced Numerical Methods in Computational Mathematics참고 문헌 48인용 수 9

한 줄 요약

이 논문은 다각형 메esh에서 변수 확산계수와 잠재항을 가진 타원형 고유값 문제에 대해 가상요소법(Virtual Element Method, VEM)의 $p$- 및 $hp$-버전을 개발하고 분석한다. 분석 결과, 해가 해석적일 경우 $p$-버전에서 지수 수렴을 입증하였고, 유한한 정규성 조건을 만족하는 해에 대해 $hp$-초기화에서 자유도의 세제곱근에 대한 지수 수렴을 보였다. 이는 체스보드 형태의 확산계수를 포함한 도전적인 문제들에 대한 광범위한 수치 실험을 통해 검증되었다.

ABSTRACT

We discuss the $p$- and the $hp$-versions of the virtual element method for the approximation of eigenpairs of elliptic operators with a potential term on polygonal meshes. An application of this model is provided by the Schrödinger equation with a pseudo-potential term. We present in details the analysis of the p-version of the method, proving exponential convergence in the case of analytic eigenfunctions. The theoretical results are supplied with a wide set of experiments. We also show numerically that, in the case of eigenfunctions with finite Sobolev regularity, an exponential approximation of the eigenvalues in terms of the cubic root of the number of degrees of freedom can be obtained by employing $hp$-refinements. Importantly, the geometric flexibility of polygonal meshes is exploited in the construction of the $hp$-spaces.

연구 동기 및 목표

• 변수 확산계수와 잠재항을 가진 타원형 고유값 문제에 대해 $p$- 및 $hp$-버전의 가상요소법을 확장한다.
• 특히 $p$-최적 근사 추정과 향상된 가상요소 공간에서의 안정화 항을 중심으로 $p$-버전의 VEM을 분석한다.
• 해석적 고유함수에 대해 $p$-버전이 지수 수렴을 보이며, 유한한 소볼레프 정규성을 가진 해에 대해 $hp$-버전이 자유도의 세제곱근에 대해 지수 수렴을 보임을 입증한다.
• 특이점과 불연속 계수를 포함한 문제들에 대해 이론적 결과를 수치 실험을 통해 검증한다. 특히 체스보드 형태의 확산 텐서를 포함한다.
• 다각형 메쉬의 기하학적 유연성을 활용하여 고유값 문제에서 $hp$-초기화 전략을 효과적으로 구현한다.

제안 방법

• VEM의 $p$-버전은 메쉬를 고정한 채로 국소 다항식 차수 $p$를 증가시키며, 국소 PDE에 의해 정의된 암묵적인 고차원 함수를 포함하는 향상된 가상요소 공간을 사용한다.
• 안정화 항은 다항식 공간 위로의 사영자와 안정성과 계산 가능성을 보장하는 이차형식을 통해 구성되며, 가상 형상 함수의 암묵적 성격에도 불구하고 이를 보장한다.
• 분석은 향상된 가상요소 공간에서의 $p$-최적 근사 추정에 기반하며, 안정화 항을 신중히 다루어 최적 수렴률을 유지한다.
• $hp$-버전의 경우, 국소 메쉬 초기화를 적응적으로 적용하며, 특히 정규성이 낮은 영역을 향해 집중적으로 수행한다. 이는 다각형 요소의 유연성을 활용한 것이다.
• 수치 실험은 내부 모멘트에 대해 $L^2$-정규직교 다항식 기저를 사용하며, 균일하고 계단형 메쉬에서 $h$-, $p$-, $hp$-초기화를 비교한다.
• 오차는 이산 고유값과 기준 고유값 간의 상대적 차이로 측정되며, 자유도의 세제곱근에 따른 수렴성을 분석한다.

실험 결과

연구 질문

• RQ1 가상요소법의 $p$-버전이 변수 확산계수와 잠재항을 가진 고유값 문제에서 지수 수렴을 달성할 수 있는가?
• RQ2 향상된 가상요소 공간에서의 $p$-최적 근사 추정은 어떻게 행동하는가? 안정화 항은 수렴 분석에서 어떤 역할을 하는가?
• RQ3 소볼레프 정규성이 유한한 고유값 문제에 대해 $hp$-버전의 VEM이 자유도의 세제곱근에 대해 지수 수렴을 보일 수 있는가?
• RQ4 불연속적이거나 극적으로 변하는 계수를 가진 문제, 예를 들어 체스보드 형태의 확산 텐서에 대해 이 방법은 어떻게 성능을 발휘하는가?
• RQ5 다각형 메쉬는 고유값 문제에서 효과적인 $hp$-초기화 전략을 지원할 수 있으며, 기존의 $h$-초기화와 비교해 어떤가?

주요 결과

• 해석적 고유함수에 대해 $p$-버전의 VEM은 철저한 $p$-최적 근사 추정과 안정화 분석을 통해 다항식 차수 $p$에 대해 지수 수렴을 달성한다.
• 소볼레프 정규성이 유한한 해에 대해 $hp$-버전의 VEM은 자유도의 세제곱근에 대해 지수 수렴을 보이며, 이는 이론적 기대와 일치한다.
• 수치 실험을 통해 균일 및 계단형 메쉬에서 $hp$-버전이 지수 수렴을 확인하였으며, 불연속적인 체스보드 형태의 확산 계수를 가진 문제에서도 동일하게 성능을 발휘한다.
• 체스보드 모델의 $ u$-파rameter가 매우 큰 경우(예: $ u = 10^8$)에도 불구하고 메서드는 계수 변화에 민감하지 않게 유지되며, 강건성을 보인다.
• 두 번째 및 세 번째 고유값의 수렴 성능 저하는 이전 연구에서 기준 고유값의 정확도 부족에서 기인하며, 방법론적 결함 때문이 아니다.
• 다각형 메쉬의 기하학적 유연성 덕분에 $hp$-초기화가 효과적으로 수행되며, 특히 정규성이 낮은 영역으로 자유도를 집중함으로써 수렴 속도를 향상시킨다.