[논문 리뷰] An efficient multigrid solver for isogeometric analysis

Álvaro Pé de la Riva, Carmen Rodrigo외 1명|arXiv (Cornell University)|2018. 06. 15.

Advanced Numerical Analysis Techniques참고 문헌 22인용 수 5

한 줄 요약

이 논문은 스퍼링 차수 $p$에 따라 변하지 않는 수렴성을 확보하기 위해 스퍼링 차수 $p$에 맞추어 블록 크기를 조정하는 겹치는 승법 색소프트너를 사용하는 이소지오메트릭 해석(IGA)을 위한 강건한 기하학적 멀티그리드 솔버를 제안한다. 이 방법은 V-사이클에서 단 한 번의 프리-스무딩 단계만으로도 스퍼링 차수 $p$와 메쉬 크기 $h$에 관계없이 수렴성을 확보하며, 국소 푸리에 분석과 복잡한 기하구조에서의 수치 실험을 통해 검증되었다.

ABSTRACT

The design of fast solvers for isogeometric analysis is receiving a lot of attention due to the challenge that offers to find an algorithm with a robust convergence with respect to the spline degree. Here, we analyze the application of geometric multigrid methods to this type of discretizations, and we propose a multigrid approach based on overlapping multiplicative Schwarz methods as smoothers. The size of the blocks considered within these relaxation procedures is adapted to the spline degree. A simple multigrid V-cycle with only one step of pre-smoothing results to be a very efficient algorithm, whose convergence is independent on the spline degree and the spatial discretization parameter. Local Fourier analysis is shown to be very useful for the understanding of the problems encountered in the design of a robust multigrid method for IGA, and it is performed to support the good convergence properties of the proposed solver. In fact, an analysis for any spline degree and an arbitrary size of the blocks within the Schwarz smoother is presented for the one-dimensional case. The efficiency of the solver is also demonstrated through several numerical experiments, including a two-dimensional problem on a nontrivial computational domain.

연구 동기 및 목표

• 스퍼링 차수 $p$가 증가함에 따라 악화되는 이소지오메트릭 해석에서 멀티그리드 수렴성의 문제를 해결한다.
• 고차원 B-스플라인 및 NURBS 이산화에 대해 효율적이고 강건한 기하학적 멀티그리드 방법을 개발한다.
• 스퍼링 차수 $p$ 증가에 따라 수렴성이 떨어지는 표준 스무딩 기법(예: 가우스-세이델)의 한계를 극복한다.
• 스퍼링 차수에 따라 블록 크기를 조정하는 스무딩 기법을 설계하여 스퍼링 차수에 적응하도록 한다.
• 복잡한 기하구조(예: 1/4 원환체)에서의 강건성을 입증하여 실용적 적용 가능성을 검증한다.

제안 방법

• 단 한 번의 프리-스무딩 단계만을 포함하는 V-사이클 아키텍처를 가진 기하학적 멀티그리드 프레임워크를 사용한다.
• 스퍼링 차수 $p$에 따라 동적으로 블록 크기를 조정하는 겹치는 승법 색소프트너를 사용한다.
• 수렴 속도 향상을 위해 승법 색소프트너의 색상화된 버전을 구현한다.
• 국소 푸리에 분석(LFA)을 사용하여 1차원에서 임의의 $p$와 블록 크기에 대해 수렴 행동을 분석하고 검증한다.
• B-스플라인과 NURBS를 사용한 정확한 기하구조 표현을 통해 다항식의 연속성 $C^{p-1}$까지 유지하는 멀티그리드 솔버를 구축한다.
• 선형 시스템을 해결하기 위해 갈레르킨 이소지오메트릭 이산화에서 유도된 강성 행렬을 사용한다.

실험 결과

연구 질문

• RQ1 이소지오메트릭 해석에서 스퍼링 차수 $p$에 대해 강건한 기하학적 멀티그리드 방법을 개발할 수 있는가?
• RQ2 색소프트너의 블록 크기 선택이 고차원 IGA 이산화에서 수렴성에 어떤 영향을 미치는가?
• RQ3 적응형 색소프트너를 사용할 경우 단 한 번의 프리-스무딩 단계만으로도 강건한 수렴성이 확보되는가?
• RQ4 국소 푸리에 분석은 1차원에서 임의의 $p$와 블록 크기에 대해 멀티그리드 솔버의 수렴 행동을 예측할 수 있는가?
• RQ5 제안된 솔버는 NURBS 기하구조를 가진 비트리비얼한 계산 영역(예: 1/4 원환체)에서도 효과적으로 작동하는가?

주요 결과

• 제안된 멀티그리드 솔버는 스퍼링 차수 $p$와 메쉬 크기 $h$에 관계없이 수렴성을 확보하며, 잔차를 $10^{-8}$로 감소시키기 위해 단지 네 또는 다섯 번의 $V(1,0)$-사이클이 필요하다.
• 1/4 원환체 영역에서 테스트된 모든 스퍼링 차수 $p=2$에서 $p=8$까지의 경우, 반복 횟수는 네 번 이내로 일정하게 유지되어 강건성을 입증한다.
• 스퍼링 차수에 적응하는 블록 크기를 가진 색상화된 승법 색소프트너를 사용할 경우, 표준 스무딩 기법보다 수렴성이 크게 향상된다.
• 국소 푸리에 분석은 1차원에서 임의의 스퍼링 차수 $p$와 임의의 블록 크기에 대해 방법의 강건성을 확인한다.
• 2차원 문제에 대한 수치 실험, 특히 비트리비얼한 1/4 원환체 영역에서의 실험 결과, 다양한 $p$와 $h$ 값에서도 일관된 성능을 보였다.
• 계산 시간은 효율적으로 스케일링되며, 고차수 $p$와 미세한 메쉬에서도 다소의 증가만을 보여, 솔버의 효율성을 확인한다.