[논문 리뷰] MultiGrid Preconditioners for Mixed Finite Element Methods of Vector Laplacian

Long Chen, Yongke Wu외 2명|arXiv (Cornell University)|2016. 01. 15.

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

한 줄 요약

이 논문은 점근적 가우스-세이델 스무딩을 사용하는 가변 V-사이클 다중격자 방법을 통해 벡터 라플라스 연산자 및 마ク스웰 방정식의 혼합 유한요소 이산화에 대한 블록 대각 및 블록 상삼각 조건자(preconditioner)를 제안한다. 이 방법은 스무딩을 통해 슈어 여백을 근사함으로써 메esh 독립적 수렴성을 달성하며, 기하학적 스무딩 기법의 복잡한 국소 해법과 기하학적 정보 의존도를 줄여, 부분적인 정규성 조건 하에서도 강력한 성능을 보이는 블랙박스, 대수적 대체 방법을 제공한다.

ABSTRACT

Due to the indefiniteness and poor spectral properties, the discretized linear algebraic system of the vector Laplacian by mixed finite element methods is hard to solve. A block diagonal preconditioner has been developed and shown to be an effective preconditioner by Arnold, Falk, and Winther [Acta Numerica, 15:1--155, 2006]. The purpose of this paper is to propose alternative and effective block diagonal and block triangular preconditioners for solving this saddle point system. A variable V-cycle multigrid method with the standard point-wise Gauss-Seidel smoother is proved to be a good preconditioner for a discrete vector Laplacian operator. This multigrid solver will be further used to build preconditioners for the saddle point systems of the vector Laplacian and the Maxwell equations with divergent free constraint. The major benefit of our approach is that the point-wise Gauss-Seidel smoother is more algebraic and can be easily implemented as a black-box smoother.

연구 동기 및 목표

• 벡터 라플라스 연산자의 혼합 유한요소 이산화에서 발생하는 사다리꼴 행렬계(system)에 대해 효과적이고 구현 가능한 조건자를 개발한다.
• 기하학적 스무딩 기법이 복잡한 국소 해법과 기하학적 정보를 필요로 하는 한계를 극복하기 위해 더 대수적인 스무딩 기법을 도입한다.
• 표준 점근적 가우스-세이델 이완을 사용하여 시스템 행렬의 슈어 여백에 대해 최적의 수렴 속도를 보장하는 다중격자 프레임워크를 수립한다.
• 완전한 정규성 조건이 성립하지 않을 경우에도 강건한 성능을 유지하면서, 수렴성 보장이 가능한 조건자를 마크스웰 방정식에 확장한다.
• 수치적 실험을 통해 제안된 조건자가 메esh 독립적 반복 횟수와 높은 CPU 성능을 달성함을 보여준다.

제안 방법

• 벡터 라플라스 연산자의 이산화를 나타내는 슈어 여백 $A = B^T M_v^{-1} B + C^T M_f C$ 에 대해 가변 V-사이클 다중격자 방법을 적용한다.
• 점근적 가우스-세이델 스무딩을 블랙박스 스무딩으로 사용하여 전용 기하학적 또는 패치 기반 스무딩 기법이 필요 없도록 한다.
• 실용성과 희소성 확보를 위해 질량 행렬의 역행렬 $M_v^{-1}$ 을 희소이고 뭉친 질량 행렬 $\tilde{M}_v^{-1}$ 으로 근사한다.
• 다중격자 근사 $\tilde{A}^{-1}$ 를 기반으로 블록 대각 및 블록 상삼각 조건자를 구성하며, $(2,2)$ 블록은 $\tilde{M}_v^{-1}$ 또는 스칼라 라플라스 연산자의 V-사이클로 근사한다.
• 브람블-파스카크-수 프레임워크를 기반으로 이론적 수렴성을 입증하며, 정규성 및 근사 가정에 의존하며, 대칭 정부호 행렬성과 가우스-세이델 이론에 의해 스무딩 성질이 보장된다.
• 수치 실험은 iFEM 패키지를 사용하여 2D 및 3D 영역, 특히 전체 정규성이 성립하지 않는 특이 영역에서 성능을 검증한다.

실험 결과

연구 질문

• RQ1 표준 점근적 가우스-세이델 스무딩기가 벡터 라플라스 연산자의 사다리꼴 행렬계의 슈어 여백에 대해 다중격자 방법에서 효과적인 블랙박스 스무딩으로 사용될 수 있는가?
• RQ2 전체 타원형 정규성 가정이 위반될 경우, 제안된 다중격자 조건자가 메esh 독립적 수렴성을 유지하는가?
• RQ3 편미분 방정식의 발산 제약 조건이 있는 마크스웰 방정식에서, 블록 대각 및 블록 상삼각 조건자 간의 반복 횟수와 계산 비용은 어떻게 비교되는가?
• RQ4 슈어 여백에 대한 다중격자 솔버가 기하학적 정보나 국소 패치 해법이 필요 없이 강력한 대수적 조건자로 사용될 수 있는가?
• RQ5 완전한 선형 다항식 공간이 이산 공간에 포함되지 않을 경우, 유한요소 근사 성질(예: 최저차수 모서리 요소)이 수렴 행동에 미치는 영향은 무엇인가?

주요 결과

• 블록 상삼각 조건자가 블록 대각 조건자보다 반복 횟수에서 뛰어난 성능을 보이며, 특히 3차원 및 볼록하지 않은 영역에서 두드러진다.
• 대각 조건자의 반복 횟수는 정규성 조건이 약화될수록 서서히 증가하여 타원형 정규성 가정에 민감함을 시사한다.
• 삼각 조건자는 전체 정규성이 실패하는 L자형 및 균열 영역에서도 강건성을 유지하며, 내부 정규성과 추가 스무딩이 경계 특이성을 보완함을 시사한다.
• 마크스웰 방정식의 경우, 두 조건자가 수렴을 달성하며 반복 횟수는 서서히 증가한다. 삼각 조건자는 $(2,2)$ 블록 해법이 더 비용이 많이 들지만 CPU 시간은 유사하다.
• 단일 V-사이클($m_J=2$ (2D), $m_J=3$ (3D))이 효과적인 조건자를 제공하여 다중격자 근사의 효율성을 입증한다.
• 유한요소 공간이 전체 선형 다항식을 포함하지 않을 경우에도 방법이 효과를 유지함으로써, 이론적 가정을 초월한 강건성을 입증한다.