[논문 리뷰] An auxiliary space multigrid preconditioner for the weak Galerkin method

Long Chen, Junping Wang외 2명|arXiv (Cornell University)|2014. 10. 04.

Advanced Numerical Methods in Computational Mathematics참고 문헌 12인용 수 5

한 줄 요약

이 논문은 2D/3D 단체 메esh에서 두 번째 차수의 확산 방정식에 적용된 약한 갈레르킨 유한요소법에 대해 보조 공간 다중격자 조건부 행렬을 제안한다. H¹-연속인 일차 다항형 유한요소 공간을 보조 거친 공간으로 사용함으로써, 메쉬 크기와 무관하게 조건수의 유계성을 확보하였으며, 수치 결과는 균일하고 적응형 메쉬에서의 강건한 성능을 확인한다.

ABSTRACT

In this paper, the authors constructed an auxiliary space multigrid preconditioner for the weak Galerkin finite element method for second-order diffusion equations, discretized on simplicial 2D or 3D meshes. The idea of the auxiliary space multigrid preconditioner is to use an auxiliary space as a "coarse" space in the multigrid algorithm, where the discrete problem in the auxiliary space can be easily solved by an existing solver. In this construction, the authors conveniently use the $H^1$ conforming piecewise linear finite element space as an auxiliary space. The main technical difficulty is to build the connection between the weak Galerkin discrete space and the $H^1$ conforming piecewise linear finite element space. The authors successfully constructed such an auxiliary space multigrid preconditioner for the weak Galerkin method, as well as a reduced system of the weak Galerkin method involving only the degrees of freedom on edges/faces. The preconditioned systems are proved to have condition numbers independent of the mesh size. Numerical experiments are conducted to support the theoretical results.

연구 동기 및 목표

• 두 번째 차수의 확산 방정식에 대한 약한 갈레르킨 유한요소 이산화에서 발생하는 선형 시스템에 대해 효율적인 조건부 행렬을 개발하기 위해.
• 약한 갈레르킨 방법에서 발생하는 큰 희소 선형 시스템을 최적의 확장성으로 해결하는 데 도전하기 위해.
• 메쉬 크기와 무관한 조건수를 갖는 강건한 조건부 행렬을 구축하여 반복 해법의 빠른 수렴을 보장하기 위해.
• 모서리/면의 자유도만 유지하는 축소된 시스템 설정을 도입하고 철저히 분석하여 계산 효율성을 향상시키기 위해.
• 기존의 조건부 행렬 공간 다중격자 프레임워크를 약한 갈레르킨 맥락으로 확장하여, 연속 유한요소 공간에서의 기존 해법을 활용하기 위해.

제안 방법

• 메쉬의 거친 문제를 효율적으로 해결할 수 있도록 H¹-연속인 일차 다항형 유한요소 공간을 보조 거친 공간으로 사용한다.
• 약한 갈레르킨 이산 공간과 연속 H¹ 공간 사이에 안정적인 전이 연산자를 설정하여 세밀한 문제와 거친 문제를 연결한다.
• 세밀한 수준에서는 약한 갈레르킨 이산화를, 거친 수준에서는 동일한 메쉬에서 연속 P1 요소를 사용하는 V-사이클 다중격자 프레임워크를 적용한다.
• 보조 공간 조건부 행렬을 사용하여 공액 기울기(PCM) 해법을 적용하여 유도된 선형 시스템을 해결한다.
• 내부 요소 자유도를 제거하여 축소된 시스템을 유도하고, 모서리/면 자유도만 유지함으로써 효율성을 향상시킨다.
• 조건부 행렬 시스템의 조건수가 메쉬 크기와 무관하게 유계임을 증명한다. 이는 메쉬 크기와 무관한 수렴 보장에 기여한다.

실험 결과

연구 질문

• RQ1 약한 갈레르킨 방법에 대해 메쉬 크기에 영향을 받지 않는 수렴성을 보장하는 보조 공간 다중격자 조건부 행렬을 구성할 수 있는가?
• RQ2 효과적인 조건부 행렬을 위해 비연속적인 약한 갈레르킨 공간과 연속 H¹ 유한요소 공간 간의 연결을 어떻게 설정할 수 있는가?
• RQ3 모서리/면 자유도만 포함하는 축소된 시스템 설정이 최적의 수렴성과 효율성을 유지하는가?
• RQ4 불연속 계수를 가진 문제나 적응형 메쉬에서 이 조건부 행렬의 강건성은 어떠한가?
• RQ5 조건부 행렬의 성능은 메쉬 크기와 무관한가? 그리고 다양한 문제 유형과 기하학적 형태에서도 효과를 유지하는가?

주요 결과

• 조건부 행렬 시스템의 조건수가 메쉬 크기와 무관하게 유계이므로, PCG 해법의 메쉬 크기와 무관한 수렴이 보장된다.
• 원래의 약한 갈레르킨 시스템에서, 예제 1~3의 모든 메쉬 세분화에서 PCG 반복 횟수가 일정하게 유지됨(예: 13~14단계), 이는 강건성을 나타낸다.
• 축소된 시스템 설정은 시스템 크기를 약 2/3 감소시키고, 원래 시스템 대비 해법 시간을 약 절반으로 줄였다.
• 계수의 변화가 극심한 예제 5(ε = 10⁻⁴에서 10⁴)에서 축소된 시스템의 반복 횟수는 항상 13으로 일정하게 유지되어 계수의 급격한 변화에 대한 강건성을 입증한다.
• 적응형 메쉬에서의 수치 실험(예제 3) 결과, 반복 횟수와 성능이 유사하게 유지되어, 준균일 메쉬를 초월한 강건성을 보여준다.
• 특히 축소된 시스템의 새로운 처리 방식 덕분에 일반성과 효율성 면에서 기존 접근법을 초월한다.