[논문 리뷰] Convergence proof for the multigrid method of the nonlocal model

Minghua Chen, Weihua Deng|arXiv (Cornell University)|2016. 05. 18.

Numerical methods in engineering참고 문헌 19인용 수 4

한 줄 요약

이 논문은 메쉬 크기 $ h \to 0 $ 일 때 비유계인 $ a_0/a_1 $ 를 갖는 페리다이내믹스(nonlocal models)에 적용된 이중그리드 방법(TGM)에 대한 엄밀한 수렴 증명을 제시한다. 대칭 토플리츠 강성행렬을 라플라시안 유사 행렬의 합으로 분해함으로써 저자들은 TGM의 균일 수렴성을 확립하고, 특정 조건 하에서 전반 다중그리드(FMGM) 및 V-사이클 스킴으로의 분석을 확장하며, 수치 실험을 통해 이차 수렴성과 $ \mathcal{O}(N\log N) $ 복잡도를 검증한다.

ABSTRACT

Recently, nonlocal models attract the wide interests of scientist. They mainly come from two applied scientific fields: peridyanmics and anomalous diffusion. Even though the matrices of the algebraic equation corresponding the nonlocal models are usually Toeplitz (denote a0 as the principal diagonal element, a1 as the trailing diagonal element, etc). There are still some differences for the models in these two fields. For the model of anomalous diffusion, a0/a1 is uniformly bounded; most of the time, a0/a1 of the model for peridyanmics is unbounded as the stepsize h tends to zero. Based on the uniform boundedness of a0/a1, the convergence of the two-grid method is well established [Chan, Chang, and Sun, SIAM J. Sci. Comput., 19 (1998), pp. 516--529; Pang and Sun, J. Comput. Phys., 231 (2012), pp. 693--703; Chen, Wang, Cheng, and Deng, BIT, 54 (2014), pp. 623--647]. This paper provides the detailed proof of the convergence of the two-grid method for the nonlocal model of peridynamics. Some special cases of the full multigrid and the V-cycle multigrid are also discussed. The numerical experiments are performed to verify the convergence.

연구 동기 및 목표

• 메쉬 크기 $ h \to 0 $ 일 때 $ a_0/a_1 $ 가 비유계인 페리다이내믹스 비국소 모델에 대해 이중그리드 방법(TGM)의 엄밀한 수렴 증명을 확립하는 것. 이는 비정상 확산 모델과는 다름.
• 기존 수렴 증명이 $ a_0/a_1 $ 의 유계성에 의존하는 한계를 극복하여, 페리다이내믹스 시스템에서는 실패하는 경우를 다루는 것.
• 특정 구조적 가정 하에 전반 다중그리드(FMGM) 및 V-사이클 다중그리드 방법의 수렴 분석을 확장하는 것.
• 이론적 결과를 수치 실험을 통해 검증하여 최적의 수렴 속도와 계산 효율성을 입증하는 것.

제안 방법

• 비국소 모델의 강성행렬을 라플라시안 유사 행렬의 합으로 표현함으로써 스펙트럼 및 조건수 성질 분석이 가능해진다.
• TGM의 수렴 증명은 이러한 라플라시안 유사 성분의 성질, 특히 특정 노름에서의 균일 유계성에 기반한다.
• 거푸질링 전략을 사용하여 거친 격자 연산자를 구성함으로써 수준 간 안정성과 근사 품질을 보장한다.
• 전반 다중그리드(FMGM) 및 V-사이클 다중그리드 스킴의 분석은 거친 격자 보정이 균일 수렴 조건을 만족한다고 가정하여 수행된다.
• 수치 실험은 $ \omega_{\text{pre}} = 1 $, $ \omega_{\text{post}} = 1/3 $ 인 V-사이클 알고리즘과 잔차 기반 정지 기준 $ \|r^{(i)}\| / \|r^{(0)}\| < 10^{-8} $ 을 사용한다.
• 계산 비용과 내구성 평가를 위해 갈레르킨과 메쉬 크기 두 배로 늘리는 두 가지 거친 격자 전략을 비교한다.

실험 결과

연구 질문

• RQ1 메쉬 크기 $ h \to 0 $ 일 때 $ a_0/a_1 $ 가 비유계인 비국소 페리다이내믹스 모델에 대해 이중그리드 방법이 균일하게 수렴할 수 있는가?
• RQ2 비유계인 $ a_0/a_1 $ 의 조건 하에서도 수렴 증명이 가능한 강성행렬의 구조적 성질은 무엇인가?
• RQ3 동일한 조건 하에서 전반 다중그리드 방법(FMGM)이 여전히 균일 수렴성을 유지하는가?
• RQ4 특정 클래스의 페리다이내믹스 문제에 대해 V-사이클 다중그리드 방법이 균일 수렴성을 달성할 수 있는가?
• RQ5 실제 적용 시 다른 거친 격자 전략이 수렴 속도와 계산 비용에 어떤 영향을 미치는가?

주요 결과

• 비유계인 $ a_0/a_1 $ 를 갖는 페리다이내믹스 비국소 모델에 대해 이중그리드 방법(TGM)이 균일하게 수렴함을 입증하였으며, 이는 라플라시안 유사 행렬로의 분해를 통해 이루어졌다.
• 수치 실험 결과, $ \delta = 1 $, $ \delta = \sqrt{h} $, $ \delta = 5h $, $ \delta = h $ 등 모든 $ \delta $ 값에 대해 $ l_\infty $ 노름에서 이차 수렴성을 확인하였으며, 수렴 속도는 약 2.0에 가까웠다.
• 다중그리드 솔버의 계산 비용은 거의 $ \mathcal{O}(N\log N) $ 의 비율로 증가하며, 각 격자 수준 간 평균 반복 횟수가 안정되어 있었다.
• 갈레르킨 전략과 메쉬 크기 두 배 전략이 거의 동일한 정확도와 수렴 행동을 보였으며, CPU 시간은 $ N $ 과 선형적으로 증가하였다.
• $ \delta = h $ 의 경우, TGM의 수준당 평균 반복 횟수는 18–23로 일정했고, FMGM은 47–56이었으며, 이는 내구성을 시사한다.
• 특수한 경우에서 거친 격자 보정이 요구되는 균일 수렴 조건을 만족함으로써 V-사이클 다중그리드의 수렴이 확인되었지만, 조건 (52)의 검증은 아직 미해결 상태이다.