[논문 리뷰] A multigrid algorithm for the $p$-version of the Virtual Element Method

Paola F. Antonietti, Lorenzo Mascotto외 1명|arXiv (Cornell University)|2017. 03. 07.

Advanced Numerical Methods in Computational Mathematics참고 문헌 37인용 수 7

한 줄 요약

이 논문은 2차원 푸아송 문제에 적용된 $p$-버전 가상요소법(Virtual Element Method, VEM)에 대해 다중격자 알고리즘을 제시한다. W-사이클 구조를 사용하며, 굴곡된 수준에서는 다항식 차수를 감소시킨다. 보조 VEM 공간을 도입하여 계산 가능한 보간 연산자를 구성함으로써, 방법은 균일 수렴성을 확보하고 PCG에 대해 균일한 조건수를 제공한다. 수치 결과는 다양한 다각형 메쉬와 다항식 차수에서 안정적인 성능을 확인한다.

ABSTRACT

We present a multigrid algorithm for the solution of the linear systems of equations stemming from the $p-$version of the Virtual Element discretization of a two-dimensional Poisson problem. The sequence of coarse spaces are constructed decreasing progressively the polynomial approximation degree of the Virtual Element space, as in standard $p$-multigrid schemes. The construction of the interspace operators relies on auxiliary Virtual Element spaces, where it is possible to compute higher order polynomial projectors. We prove that the multigrid scheme is uniformly convergent, provided the number of smoothing steps is chosen sufficiently large. We also demonstrate that the resulting scheme provides a uniform preconditioner with respect to the number of degrees of freedom that can be employed to accelerate the convergence of classical Krylov-based iterative schemes. Numerical experiments validate the theoretical results.

연구 동기 및 목표

• 2차원 푸아송 문제의 $p$-버전 VEM 이산화에서 발생하는 선형계를 효율적으로 풀 수 있는 반복적 해법을 개발한다.
• 일반적인 다각형 메쉬에서 표준 VEM 자유도가 보간을 위해 충분하지 않은 상황에서 일관된 굴곡된 수준의 해법을 구성하는 데 도전한다.
• 충분한 스무딩 단계가 확보된 조건 하에서 $p$-다중격자 W-사이클 알고리즘의 균일 수렴성을 증명한다.
• 다중격자 방법이 프리컨dition드 코시-그래디언트(PCG) 방법에 대해 균일한 조건수를 제공함을 보여준다.
• 구조적 및 비구조적 다각형 메쉬에서의 수치 실험을 통해 이론적 결과를 검증한다.

제안 방법

• W-사이클 $p$-다중격자 알고리즘을 사용하며, 굴곡된 수준은 VEM 공간의 다항식 차수를 점진적으로 감소시켜 구성한다. 최소 선형 요소까지 내려간다.
• 표준 자유도만으로는 보간 연산자를 구성할 수 없기 때문에, 계산 가능한 보간 연산자를 가능하게 하기 위해 보조 VEM 공간을 도입한다.
• 인터스페이스 연산자는 보조 공간에 존재하는 고차 다항식 프로젝터를 사용하여 구성되며, 수준 간 일관성을 확보한다.
• 이론적 분석에서는 균일 수렴성을 증명하기 위해 리처드슨 스무딩을 사용하지만, 수치 실험에서는 수렴 속도 향상을 위해 대칭화된 가우스-세이델 스무딩을 사용한다.
• 다중격자 방법은 PCG 방법의 조건수로 테스트되었으며, 잔차 감소율이 $10^{-6}$ 이하일 때 수렴 기준으로 삼는다.
• 정규 다각형 메쉬, 보로노이-로이드 메쉬, 정사각형 메쉬에서 수치 실험을 수행하여 다양한 메쉬 유형에서의 강건성 평가를 수행한다.

실험 결과

연구 질문

• RQ1 일반적인 다각형 메쉬에서 $p$-버전 VEM에 대해 $p$-다중격자 알고리즘을 구성할 수 있는가?
• RQ2 표준 VEM 자유도가 보간을 위해 충분하지 않을 경우, 어떻게 일관된 보간 연산자를 설계할 수 있는가?
• RQ3 제안된 다중격자 방법은 자유도 수에 관계없이 균일하게 수렴하는가?
• RQ4 다중격자 방법이 PCG와 같은 코시-그래디언트 계열의 Krylov 해법에 대해 균일한 조건수를 제공할 수 있는가?
• RQ5 스무딩 방법의 선택(Richardson vs. Gauss–Seidel)이 알고리즘의 수렴 행동에 어떤 영향을 미치는가?

주요 결과

• 충분한 스무딩 단계가 확보된 조건 하에서, 고전적 다중격자 프레임워크를 사용해 $p$-다중격자 W-사이클 알고리즘이 균일하게 수렴함을 증명하였다.
• 수렴을 위해 필요한 최소 후방 스무딩 단계 수는 다항식 차수 $p$에 따라 달라지며, 다양한 메쉬 유형에서 $p=2$부터 $p=6$까지 2에서 4 사이의 값을 가진다.
• 대칭화된 가우스-세이델 스무딩을 사용할 경우 리처드슨 스무딩보다 수렴 인자가 크게 감소하며, 보로노이-로이드 메쉬에서 $p=5$ 및 $m=14$일 때 수렴 인자가 0.35 이하로 떨어진다.
• 다중격자 조건수 방법은 PCG 반복 횟수가 $p$ 증가함에 따라 거의 일정하게 유지됨을 보여주며, 특히 가우스-세이델 스무딩을 사용할 경우 균일한 조건수를 제공함을 확인하였다.
• 부분 코시-데컴포지션 조건수 방법은 제안된 다중격자 조건수 방법과 달리 균일 수렴을 제공하지 못한다.
• 이 방법은 정규형(육각형) 및 비정규형(보로노이-로이드) 다각형 메쉬에서 모두 강건한 성능을 유지하며, 일반성과 확장 가능성을 입증한다.